Hybrid rocket motors with paraffin-based fuels are of interest due to higher regression rates compared to other polymers. During paraffin combustion, a liquid layer forms on the fuel surface that, together with shearing forces from the oxidizer flow, results in the formation of instabilities at the fuel-oxidizer interface. These instabilities lead to the formation and entrainment of heterogeneous sized liquid droplets into the main flow and the combusting droplets result in higher motor output. The atomization process begins with droplet formation and ends with droplet pinch-off. The goal of this paper is to conduct an uncertainty quantification (UQ) analysis of the pinch-off process characterized by a pinch-off volume ($V_{po}$) and time ($t_{po}$). We study these
quantities of interest (QoIs) in the context of a slab burner setup. We have developed a computationally expensive mathematical model that describes droplet formation under external forcing and trained an inexpensive Gaussian Process surrogate of the model to facilitate UQ. We use the pinch-off surrogate to forward propagate uncertainty of the model inputs to the QoIs and conduct two studies: one with gravity present and one without gravity effects. After forward-propagating the uncertainty of the inputs using the surrogate, we concluded that both QoIs have right-skewed distributions, corresponding to larger probability densities towards smaller pinch-off volumes and times. Specifically, for the pinch-off times, the resulting distributions reflect the effect of gravity acting against droplet formation, resulting in longer pinch-off times compared to the case where there is no gravity.
more »
« less
Simulating Paraffin Wax Droplets Using Mixed Finite Element Method
Paraffin wax is a prominent solid fuel for hybrid rockets. The atomization process of the paraffin wax fuel into he hybrid rocket combustion involves the droplets pinching off from the fuel surface. Therefore, droplet formation and pinch-off dynam- ics is analyzed using a one-dimensional axisymmetric approximation to understand droplet size distribution and pinch-off time. A mixed finite element formulation is used to solve the numerical problem. The computational algorithm uses adaptive mesh refinement to capture singularity and runs self-consistently to calculate droplet elongation. The code is verified using the Method of Manufactured Solution (MMS) and validated against laboratory experiments. Moreover, paraffin wax simulations are explored for varying inlet radius and it is found that the droplet size increases very slightly with the increasing inlet radius. Also, the pinch-off time increases up to a point where it starts to decrease as we increase the inlet radius. This behavior leads to a conjecture for the theoretical maximum radius that the droplet approaches as the inlet radius increases, which is a motivation for the future work.
more »
« less
- Award ID(s):
- 1931524
- NSF-PAR ID:
- 10401275
- Editor(s):
- Brehm, Christop; Pandya, Shishir
- Date Published:
- Journal Name:
- ICCFD 11 PROCEEDINGS
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Droplet formation happens in finite time due to the surface tension force. The linear stability analysis is useful to estimate the size of a droplet but fails to approximate the shape of the droplet. This is due to a highly nonlinear flow description near the point where the first pinch-off happens. A one-dimensional axisymmetric mathematical model was first developed by Eggers and Dupont [“Drop formation in a one-dimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] using asymptotic analysis. This asymptotic approach to the Navier–Stokes equations leads to a universal scaling explaining the self-similar nature of the solution. Numerical models for the one-dimensional model were developed using the finite difference [Eggers and Dupont, “Drop formation in a one-dimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] and finite element method [Ambravaneswaran et al., “Drop formation from a capillary tube: Comparison of one-dimensional and two-dimensional analyses and occurrence of satellite drops,” Phys. Fluids 14, 2606–2621 (2002)]. The focus of this study is to provide a robust computational model for one-dimensional axisymmetric droplet formation using the Portable, Extensible Toolkit for Scientific Computation. The code is verified using the Method of Manufactured Solutions and validated using previous experimental studies done by Zhang and Basaran [“An experimental study of dynamics of drop formation,” Phys. Fluids 7, 1184–1203 (1995)]. The present model is used for simulating pendant drops of water, glycerol, and paraffin wax, with an aspiration of extending the application to simulate more complex pinch-off phenomena.more » « less
-
We experimentally investigate the interfacial instability that emerges when a water droplet is deposited on a bath of glycerol-water solution. Despite the absence of surface tension to stabilize short-wavelength undulations, we observe finite-size fingers that grow and pinch off from the drop. We show that the fingering patterns formed in the experiments resultes from a balance between the outward buoyancy effect and inward Marangoni flow. This induced Marangoni flow inhibits small perturbations and acts as an effective surface tension on the miscible interface of the spreading drop. To characterize the final size and shape of the drop, we perform systematic experiments by varying the drop volume and the glycerol-water volume fraction. In addition, we have developed scaling arguments for the drop’s final radius using key physical forces, and show that the final wavelength is inversely proportional to the Bond number.more » « less
-
null (Ed.)When a droplet is generated, the ligament connecting the drop to the nozzle thins down and eventually pinches off. Adding solid particles to the liquid phase leads to a more complex dynamic, notably by increasing the shear viscosity. Moreover, it introduces an additional length scale to the system, the diameter of the particles, which eventually becomes comparable to the diameter of the ligament. In this paper, we experimentally investigate the thinning and pinch-off of drops of suspensions with two different sizes of particles. We characterize the thinning for different particle size ratios and different proportions of small particles. Long before the pinch-off, the thinning rate is that of an equivalent liquid whose viscosity is that of the suspension. Later, when the ligament thickness approaches the size of the large particles, the thinning accelerates and leads to an early pinch-off. We explain how the bidisperse particle size distribution lowers the viscosity by making the packing more efficient, which speeds up the thinning. This result can be used to predict the dynamics of droplet formation with bidisperse suspensions.more » « less
-
ASME (Ed.)Research was conducted to determine combustion characteristics such as: ignition delay (ID), combustion delay (CD), combustion phasing (CA 50), combustion duration, derived cetane number (DCN) and ringing intensity (RI) of F24, for its compatibility in Common Rail Direct Injection (CRDI) compression ignition (CI) engine. The first part of this study is investigating the performance of Jet-A, F24, and ultra-low sulfur diesel #2 (ULSD) using a constant volume combustion chamber (CVCC) followed by experiments in a fired CRDI research engine. Investigations of the spray atomization and droplet size distribution of the neat fuels were conducted with a Malvern Mie scattering He-Ne laser. It was found that the average Sauter Mean Diameter (SMD) for Jet-A and F24 are similar, with both fuels SMD droplet range between 25–29 micrometers. Meanwhile, ULSD was found to have a larger SMD particle size in the range of 34–40 micrometers. It was observed during the study, utilizing the CVCC, that the ID and CD for neat ULSD and Jet-A are nearly identical while the combustion of F24 is delayed. F24 was found to have longer durations of both ID and CD by approx. 0.5 ms. This results in a lower DCN for the fuel of 43.5, whereas ULSD and Jet-A have DCNs of 45 and 47 respectively. The peak AHRR for ULSD and Jet-A are nearly identical, whereas F24 has a peak magnitude of approx. 20% lower than ULSD and Jet-A. It was found that both aviation fuels had significantly fewer ringing events occurring after peak high temperature heat release (HTHR), a trend also observed in the CRDI research engine. Neat F24, Jet-A and ULSD were researched in the experimental engine at the same thermodynamic parameters: 5 bar indicated mean effective pressure (IMEP), 50°C (supercharged and EGR) inlet air temperature, 1500 RPM, start of injection (SOI) 16°BTDC, and 800 bar of fuel rail injection pressure as the baseline parameters in order to observe their ignition behavior, low temperature heat release, combustion phasing, and combustion duration. It was found that the ignition delay of F24 and Jet-A was greater than ULSD, approx. 5% for both aviation fuels. This ignition delay also affected the combustion phasing, or CA 50, of the aviation fuels. The CA 50 of the aviation fuels was delayed by approx. 2% compared to ULSD. Jet-A had a nearly identical combustion duration compared to ULSD, however F24 had an extended combustion duration which was approx. 3% longer than that of ULSD and Jet-A. It was discovered with the accumulations of these delays in ID, CD, CA50, that the RI of the aviation fuels were reduced. F24 was discovered to have more delays, and the RI correlates with these results having a 70% reduction in RI compared to ULSD.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Conference Paper:
- The DOI is not currently available.
