Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “droplets”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Machine learning accelerated turbulence modeling of transient flashing jets

We report modeling the sudden depressurization of superheated liquids through nozzles is a challenge because the pressure drop causes rapid flash boiling of the liquid. The resulting jet usually demonstrates a wide range of structures, including ligaments and droplets, due to both mechanical and thermodynamic effects. As the simulation comprises increasingly numerous phenomena, the computational cost begins to increase. One way to moderate the additional cost is to use machine learning surrogacy for specific elements of the calculation. This study presents a machine learning-assisted computational fluid dynamics approach for simulating the atomization of flashing liquids accounting for distinct stages, from primary atomization to secondary breakup to small droplets using the Σ - Y model coupled with the homogeneous relaxation model. Notably, the models for thermodynamic non-equilibrium (HRM) and Σ - Y are coupled, for the first time, with a deep neural network that simulates the turbulence quantities, which are then used in the prediction of superheated liquid jet atomization. The data-driven component of this method is used for turbulence modeling, avoiding the solution of the two-equation turbulence model typically used for Reynolds-averaged Navier-Stokes simulations for these problems. Both the accuracy and speed of the hybrid approach are evaluated, demonstrating adequate accuracy and at least 25% faster computational fluid dynamics simulations than the traditional approach. This acceleration suggests that perhaps additional components of the calculation could be replaced for even further benefit. Published under an exclusive license by AIP Publishing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Measuring the density, viscosity, and surface tension of molten titanates using electrostatic levitation in microgravity

Rare earth and barium titanates are useful as ferroelectric, dielectric, and optical materials. Measurements of their thermophysical properties in the liquid state can help guide melt processing technologies for their manufacture and advance understanding of fragile liquids' behavior and glass formation. Here, we report the density, thermal expansion, viscosity, and surface tension of molten BaTi 2 O 5 , BaTi 4 O 9 , and 83TiO 2 -17RE 2 O 3 (RE = La or Nd). Measurements were made using electrostatic levitation and droplet oscillation techniques in microgravity, which provide access to quiescent liquid droplets and deep supercooling of 510–815 K below the equilibrium melting points. Densities were measured over 900–2400 K. Viscosities were similar for all four compositions, increasing from ~10 mPa s near 2100 K to ~30 mPa s near 1750 K. Surface tensions were 450–490 dyn cm -1 for the rare earth titanates and 383–395 dyn cm -1 for the barium titanates; surface tensions of all compositions had small or negligible temperature dependence over 1700–2200 K. In conclusion, for solids recovered after melt quenching, x-ray microtomography revealed the fracture mechanics in crystalline products and minimal internal porosity in glass products, likely arising from entrapped gas bubbles. Internal microstructures were generally similar for products processed either in microgravity or in a terrestrial aerodynamic levitator.

42 ENGINEERING↗
DOE OSTI2025

Universality in diffusion-controlled nucleation and growth

Nucleation and growth are studied in a system that undergoes diffusion-controlled condensation under gradual changes in parameters, such as cooling. Here, it is demonstrated that when the Gibbs–Thompson effect becomes negligible, the system falls into a universal regime. i.e., the final droplet size distribution remains invariant under certain rescaling of system parameters. An approximate yet very accurate analytic form is obtained for the final droplet size distribution in this regime.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Photo-printing of faceted DNA patchy particles

Patchy particles with shape complementarity can serve as building blocks for assembling colloidal superstructures. Alternatively, encoding information on patches using DNA can direct assembly into a variety of crystalline or other preprogrammed structures. Here, we present a tool where DNA is used both to engineer shape and to encode information on colloidal particles. Two reactive oil emulsions with different but complementary DNA (cDNA) brushes are assembled into CsCl-like crystalline lattices. The DNA brushes are recruited to and ultimately localized at the junctions between neighboring droplets, which gives rise to DNA-encoded faceted patches. The emulsions are then solidified by ultraviolet (UV) polymerization, producing faceted patchy particles. The facet size and DNA distribution are determined by the balance between the DNA binding energy and the elastic deformation energy of droplets. This method leads to a variety of new patchy particles with directional interactions in scalable quantities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Enhancement and maximum in the isobaric specific-heat capacity measurements of deeply supercooled water using ultrafast calorimetry

Knowledge of the temperature dependence of the isobaric specific heat (C p ) upon deep supercooling can give insights regarding the anomalous properties of water. If a maximum in C p exists at a specific temperature, as in the isothermal compressibility, it would further validate the liquid–liquid critical point model that can explain the anomalous increase in thermodynamic response functions. The challenge is that the relevant temperature range falls in the region where ice crystallization becomes rapid, which has previously excluded experiments. Here, we have utilized a methodology of ultrafast calorimetry by determining the temperature jump from femtosecond X-ray pulses after heating with an infrared laser pulse and with a sufficiently long time delay between the pulses to allow measurements at constant pressure. Evaporative cooling of ~15-µm diameter droplets in vacuum enabled us to reach a temperature down to ~228 K with a small fraction of the droplets remaining unfrozen. We observed a sharp increase in C p , from 88 J/mol/K at 244 K to about 218 J/mol/K at 229 K where a maximum is seen. The C p maximum is at a similar temperature as the maxima of the isothermal compressibility and correlation length. From the C p measurement, we estimated the excess entropy and self-diffusion coefficient of water and these properties decrease rapidly below 235 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Asymmetric fluctuations and self-folding of active interfaces

We study the structure and dynamics of the interface separating a passive fluid from a microtubule-based active fluid. Turbulent-like active flows power giant interfacial fluctuations, which exhibit pronounced asymmetry between regions of positive and negative curvature. Experiments, numerical simulations, and theoretical arguments reveal how the interface breaks up the spatial symmetry of the fundamental bend instability to generate local vortical flows that lead to asymmetric interface fluctuations. The magnitude of interface deformations increases with activity: In the high activity limit, the interface self-folds invaginating passive droplets and generating a foam-like phase, where active fluid is perforated with passive droplets. These results demonstrate how active stresses control the structure, dynamics, and break-up of soft, deformable, and reconfigurable liquid–liquid interfaces.

active fluid↗
DOE OSTI2020

Annular Flow Simulation Supported by Iterative In-Memory Mesh Adaptation

Various flow regimes exist in a boiling water reactor (BWR) as the steam quality increases in the uprising coolant flow, from bubbly flow, slug/churn flow, to annular flow. The annular flow is characterized by the presence of a fast-moving gas core and the surrounding liquid film flowing on the conduit wall. Additionally, entrained droplets can be observed in the gas core with ingested bubbles in the liquid film. The dynamics occurring on the wavy interface between the liquid film and gas core plays a crucial role in affecting the heat transfer rate and pressure drop within the BWR core. However, a fundamental understanding of annular flow is still lacking, partly due to the difficulty in obtaining detailed local data in annular flow experiments. In the current study, a novel simulation framework is developed for the annular flow by coupling a computational fluid dynamics flow solver with state-of-the-art meshing software. The gas-liquid interface is tracked with the level set method. Based on the computed flow solutions, the computational mesh is dynamically adapted in memory to meet the local mesh resolution requirement. This iterative simulation-adaptation framework can ensure the fine mesh resolution across the interface, which not only helps mitigate the mass conservation degradation known to level set methods but also improves the representation of dramatic interface topological changes such as wave breaking and droplet entrainment. The present investigation will shed light onto the complex interfacial processes involved in annular flow and generate much needed simulation data for annular flow modeling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Liquid-to-gas transfer of sodium in a liquid cathode glow discharge

Abstract Plasma-liquid interactions have been extensively studied with a focus on the transport of reactive species from the plasma to the liquid phase and their induced liquid phase chemistry and resulting applications. While solute transfer from the liquid to the gas phase in plasmas has been widely used in analytical chemistry, the underlying processes remain relatively unexplored. We report spatially and temporally resolved absolute density measurements of sodium in a plasma with a NaCl solution cathode using two-photon absorption laser induced fluorescence (TaLIF). The observed non-linear increase in sodium density with solution conductivity is shown to correlate with droplet generation as visualized by Mie scattering. The findings are explained by droplet generation by electrospray induced by Taylor cone formation as underpinning mechanism for the introduction of sodium in the plasma. An analytical sheath model combined with a scaling law shows an increase in electric field force with solution conductivity that is consistent with the observed non-linear increase in sodium density in the plasma with solution conductivity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Isolating the effect of biomass burning aerosol emissions on 20th century hydroclimate in South America and Southeast Asia

Abstract Biomass burning is a significant source of aerosol emissions in some regions and has a considerable impact on regional climate. Earth system model simulations indicate that increased biomass burning aerosol emissions contributed to statistically significant decreases in tropical precipitation over the 20th century. In this study, we use the Community Earth System Model version 1 Large Ensemble (CESM1-LENS) experiment to evaluate the mechanisms by which biomass burning aerosol contributed to decreased tropical precipitation, with a focus on South America and Southeast Asia. We analyze the all-but-one forcing simulations in which biomass burning aerosol emissions are held constant while other forcings (e.g., greenhouse gas concentrations) vary throughout the 20th century. This allows us to isolate the influence of biomass burning aerosol on processes that contribute to decreasing precipitation, including cloud microphysics, the radiative effects of absorbing aerosol particles, and alterations in regional circulation. We also show that the 20th century reduction in precipitation identified in the CESM1-LENS historical and biomass burning experiments is consistent across Coupled Model Intercomparison Project Phase 5 models with interactive aerosol schemes and the CESM2 single-forcing experiment. Our results demonstrate that higher concentrations of biomass burning aerosol increases the quantity of cloud condensation nuclei and cloud droplets, limiting cloud droplet size and precipitation formation. Additionally, absorbing aerosols (e.g., black carbon) contribute to a warmer cloud layer, which promotes cloud evaporation, increases atmospheric stability, and alters regional circulation patterns. Corresponding convectively coupled circulation responses, particularly over the tropical Andes, contribute to further reducing the flow of moisture and moisture convergence over tropical land. These results elucidate the processes that affect the water cycle in regions prone to biomass burning and inform our understanding of how future changes in aerosol emissions may impact tropical precipitation over the 21st century.

Environmental Sciences & Ecology↗
DOE OSTI2021

Vortex rings from high energy central 𝑝 + 𝐴 collisions

Relativistic p + A collisions may produce droplets of quark gluon plasma (QGP) that quickly develop a toroidal vortex structure similar to that of an expanding smoke ring. We present viscous relativistic hydrodynamic calculations of ultra-central p + A collisions and develop an experimental observable to probe the structure, correlating the polarization and momentum of hyperons emitted from the collision. This effect is robust against changes in the definition of vorticity used to calculate the polarization. Experiments at RHIC and LHC may test the existence and strength of the vortex toroids, bringing new evidence to bear on the question of collectivity in the smallest QGP droplets.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Hypertriton puzzle in relativistic heavy-ion collisions

The yields of hadrons and light nuclei in relativistic collisions of heavy-nuclei at a center of mass energy of $\sqrt{s_{NN}}$ = 2.6 TeV can be described remarkably well by a thermal distribution of an ideal gas of hadrons and light nuclei interacting only via the decay of resonances. Given the particularly small binding energy of hypertritons relative to the temperature describing the yields (about 156 MeV), one might naturally expect hypertrions to dissociate in medium, making the agreement of hypertriton yields with thermal predictions highly puzzling. The puzzle is compounded by the fact that small binding energy is associated with the large size of the hypertriton. This size is on a similar scale to the overall size of the fireball and much larger than the length scale over which temperatures in the fireball vary over phenomenologically relevant amounts. Here, this paper quantifies the tension this effect causes and shows that it is sufficiently large to render the thermal model inconsistent: its natural assumptions are in conflict with its outputs. The possibility that hypertritons are formed at freeze out as compact objects, quark droplets, that subsequently evolve into hypertritons is considered as a way to resolve the puzzle. It is noted that beyond making the assumption that compact quark droplets form, additional detailed dynamical assumptions which have not been justified are needed to make the thermal model work. The issue of why, despite these issues, the hypertriton is well described by a simple statistical description at freeze out is unresolved. Resolving the hypertriton puzzle is important as it may clarify whether the phenomenological success of the simple thermal model for yields accurately reflects the simple picture of the underlying physics on which it is based.

hydrodynamic models↗
DOE OSTI2020

Ring patterns generated by an expanding colloidal meniscus

The drop-and-dry is a common technique allowing for creation of periodic nanoparticle (NP) structures for sensing, photonics, catalysis, etc. However, the reproducibility and scalability of this approach for fabrication of NP-based structures faces serious challenges due to the complexity of the simple, at first glance, evaporation process. In this work we study the effect of the spatial confinement on the NP self-assembly under slow solvent evaporation, when the air-liquid-substrate contact line (CL) expands from the center towards the walls of a cylindrical cell, forming a toroid. Using in situ video monitoring of the stick-slip CL motion, we find regular hydrodynamic perturbations in the meniscus, and reveal fine details of the formation of quasiperiodic rings of close packed NP layers. We report that drying of the toroidal NP droplet has a number of important differences from drying of the classical hemispherical colloidal drops. In toroidal drops we observe linear-in-time average meniscus motion, in contrast to the hemispherical drops where the meniscus moves as a square root of time. While both droplet geometries produce NP ring patterns, the ring width for the toroidal drop decreases with increasing ring radius, while it decreases with decreasing the radius of the hemispherical drop. We suggest that free ligands are the main cause of the Marangoni instabilities driving the periodic vorticity in the meniscus. Finally, we show that the usually ignored contact line tension may yield a considerable contribution to the CL pinning causing the CL slip-stick motion and the ring formation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Strong Constraints on Jet Quenching in Centrality-Dependent p + Pb Collisions at 5.02 TeV from ATLAS

Jet quenching is the process of color-charged partons losing energy via interactions with quark-gluon plasma droplets created in heavy-ion collisions. The collective expansion of such droplets is well described by viscous hydrodynamics. Similar evidence of collectivity is consistently observed in smaller collision systems, including p p and p + Pb collisions. In contrast, while jet quenching is observed in Pb + Pb collisions, no evidence has been found in these small systems to date, raising fundamental questions about the nature of the system created in these collisions. The ATLAS experiment at the Large Hadron Collider has measured the yield of charged hadrons correlated with reconstructed jets in 0.36 nb -1 of p + Pb and 3.6 pb -1 of pp collisions at 5.02 TeV. The yields of charged hadrons with ${p}^{ch}_{T}$ > 0.5 GeV near and opposite in azimuth to jets with ${p}^{jet}_{T}$ > 30 or 60 GeV, and the ratios of these yields between p + Pb and pp collisions, I pPb , are reported. The collision centrality of p + Pb events is categorized by the energy deposited by forward neutrons from the struck nucleus. The I pPb values are consistent with unity within a few percent for hadrons with ${p}^{ch}_{T}$ > 4 GeV at all centralities. These data provide new, strong constraints that preclude almost any parton energy loss in central p + Pb collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Self-Cleaning Hybrid Hydrophobic–Hydrophilic Surfaces: Durability and Effect of Artificial Soilant Particle Type

Dew accelerates the soiling rates and increases the dust adhesion. To use dew for self-cleaning, a fluorinated ethylene propylene coating was applied to suppress the reactions between dust and glass as well as to facilitate the dew to condense as mobile droplets. An array of rectangular hydrophilic channels in the coating increases the condensation rates and droplet slide-off diameters. Furthermore, the durability of the coating was evaluated by artificial UV weathering. Four different types of soilants were used in the artificial soiling tests to assess the effect of soilant type and surface properties on the soiling rates and self-cleaning efficacy under the simulated dew conditions. Soil deposition and self-cleaning mechanisms are being reported.

42 ENGINEERING↗
DOE OSTI2022

Numerical Modeling of Atmospheric Rime Ice Accretion on an Airfoil Using an Eulerian Approach

This paper presents an approach to numerically simulate the inherently unsteady rime ice accretion problem on a two-dimensional airfoil and elucidate the associated variations under different icing conditions. The airflow field and the water impingement on the airfoil are obtained based on an Eulerian two-phase model. A dynamic mesh strategy is employed to unsteadily account for the changes in the ice profile and its impact on the air and droplet flow by continuously reconstructing the computational grid at each time-step through smoothing and layering mechanisms. All main icing modules including the airflow field, droplet trajectory, icing thickness profile, and mesh management are fully coupled within the same computational framework without resorting to any external tools. Classical icing theory is employed to model the rime ice roughness, and it is assumed that the ice accretes in a direction normal to the airfoil surface. The governing Reynolds-averaged Navier–Stokes (RANS) conservation equations along with the energy and continuity equations are solved to produce the velocity and temperature fields. A convective film heat transfer coefficient is computed based on the surface heat flux and a recovery temperature which takes into account the dissipative heat release in the boundary layer in the vicinity of the airfoil surface. With the implemented strategy and calculating the convective heat transfer coefficient, the water film thickness is also calculated along with the ice shape. Furthermore, the model is validated by comparing the local collection efficiency distribution and ice shape with experimental data, and the results show that the implemented approach provides acceptable predictions of ice accretion profiles and rates.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Large Eddy Simulation of Gasoline Sprays in a Lagrangian–Eulerian Framework Using the High-Order Spectral Element Method

Predicting the spray evolution using simulations requires accurate modeling of the turbulent gas-phase flow field. Here, in this study, the high-order spectral-element method (SEM), implemented in the code Nek5000, was used to provide highly resolved solutions to the turbulent flow field. Spray modeling capabilities were implemented into the Nek5000 code. The spray is modeled in a Lagrangian–Eulerian (LE) framework, where the liquid is represented by discrete parcels of droplets. The method for coupling liquid and gas in the context of SEM is described, which allows for very fine meshes to be used without affecting the stability of the solution. Large-eddy simulations (LES) of the eight-hole ECN Spray G gasoline injector were conducted. Numerical results are compared against experimental data for liquid penetration, droplet size and gas velocity. The morphology of the multiplume spray is compared against experimental data. The effect of different spray injection inputs is analyzed. It was found that using a plume direction of 33 deg and an injection cone angle of 30 deg produced the best results overall. This work shows the applicability of SEM for spray modeling applications, where use of a high-order flow solver can help us understand the multiplume spray aerodynamics and how it leads to plume collapse under certain conditions. Results also highlight the need for tuning spray input parameters in the LE framework, even when high-fidelity gas flow solutions are possible.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Evaluation of Spray and Combustion Models for Simulating Dilute Combustion in a Direct-Injection Spark-Ignition Engine

Dilute combustion in spark-ignition engines has the potential to improve thermal efficiency by mitigating knock and by reducing throttling and wall heat losses. However, ignition and combustion processes can become unstable for dilute operation due to a lowered laminar flame speed, resulting in excessive cycle-to-cycle variability (CCV) of the combustion process. To compensate for the slower combustion in less reactive mixtures, a modified intake port geometry can be employed to generate a strong tumble flow in the cylinder and elevate turbulence levels around the spark plug, thereby promoting a faster transition to turbulent deflagration. Consequently, optimizing combustion chamber geometry and operating strategy is crucial to maximizing the benefits of using dilute combustion with enhanced in-cylinder turbulence across a wide range of operating conditions. Computational fluid dynamics (CFD) simulations can be utilized for virtual engine optimization tasks, but this would require the models to be truly predictive regarding the impact of changes to the engine design and operational parameters.In this study, multicycle large-eddy simulations (LES) are performed for a direct-injection spark-ignition engine to investigate the model performance in predicting engine combustion characteristics with respect to changes in the intake configuration. A tumble plate that blocks the lower part of the intake port inlet is used to vary the tumble. A set of CFD models that have been recently developed are employed, which takes into account the drag of nonspherical droplets, flash-boiling behavior of liquid sprays, spray-wall interaction, surrogate formulation of a research-grade E10 gasoline, and fast chemical kinetic solvers. Simulation results are compared to experimental engine data in terms of cylinder pressure, apparent heat release rate, mass fraction burned timing, and flame images. It is found that LES employing the state-of-the-art CFD models are capable of properly predicting the spray processes and reproducing the measured mean cylinder pressure for the case with the tumble plate. On the other hand, the LES over-predicts the combustion rate during the early combustion stage and under-estimates the CCV, and these discrepancies become larger when the tumble plate is removed.

computational fluid dynamics simulation↗
DOE OSTI2024

Mesofluidic Oil-Water Separation

Separating flowing multiphase solutions in process at industrial scale is challenging. A rapid continuous approach is desired. These mesofluidic oil-water separation tests were motivated by reports by Dijkshoorn et al [1] who reported separating droplets and particles using displacement with regularly spaced channels using a cross-flow microsieve.. Droplets, including deformable particles, larger than the channel migrate in one direction and thereby deplete in the other direction. Test flowrates were motivated by Burns et al. [4] who showed that separation could occur under laminar conditions, and Pease et al [8, 9] who extended the range to turbulent flow. Testing was conducted using two two-component oil-water mixtures at room temperature. Analysis showed separation of both mineral oils from their oil water mixtures at flowrates from 1-2 L/min and concentrations from 2-10 wt% oil in water.

cloud height, effective cleaning radius, effective↗