Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Supercooling”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Molecular Dynamics Simulation and Theoretical Analysis of Structural Relaxation, Bond Exchange Dynamics, and Glass Transition in Vitrimers

Vitrimers are a class of polymer networks featuring dynamic covalent cross-links that can undergo associative bond exchange. There has been recent interest in these materials due to their promise as recyclable thermosets or self-healing polymers because of the ability of vitrimer networks to rearrange at the molecular level and undergo macroscopic flow. However, the practical use of these materials often occurs in the supercooled regime or glassy state, where the implications of dynamic bonds are complicated by the interplay between slow activated segmental dynamics, cross-link (i.e., bond-exchange) kinetics, and ultimately material properties. Here, in this paper, we combine coarse-grained molecular dynamics simulation and microscopic statistical mechanical theory to understand how cross-linking kinetics affect material dynamics and how this couples to segmental relaxation of the polymeric network strands across a spectrum of length and time scales, especially in the supercooled regime. We characterize the Kuhn segmental alpha relaxation time and bond exchange time for vitrimer systems across various cross-link densities, temperatures, and bond exchange rates. Simulation and theory both exhibit a bending-up behavior for bond exchange time upon cooling, suggesting a coupling between bond exchange dynamics and segmental relaxation that intensifies with faster bond exchange kinetics. We also found bond exchange dynamics have an impact on Kuhn segment alpha relaxation time, which is most significant at higher cross-link densities. Both these effects are most prominent when the bond exchange time is similar to the Kuhn segment alpha relaxation time, and the resulting coupling of these two relaxation processes is tied to both the probability of a free end to find a bonded pair and the time scale of the constraints imposed by the dynamic cross-links. This relationship is reflected by a cross-link dependence of a theoretical parameter which represents the quantitative degree of coupling between bond exchange and segmental dynamics. Overall, the combination of simulation and theory clarifies the intricate interaction between bond kinetics and segmental relaxation and demonstrates the ability to provide molecular-level insights into vitrimer dynamics over a wide temperature range.

dynamic relaxation↗
DOE OSTI2025

Evaluating Cloud Properties at Scott Base: Comparing Ceilometer Observations With ERA5, JRA55, and MERRA2 Reanalyses Using an Instrument Simulator

This study compares CL51 ceilometer observations made at Scott Base, Antarctica, with statistics from the ERA5, JRA55, and MERRA2 reanalyses. To enhance the comparison we use a lidar instrument simulator to derive cloud statistics from the reanalyses which account for instrumental factors. The cloud occurrence in the three reanalyses is slightly overestimated above 3 km, but displays a larger underestimation below 3 km relative to observations. Unlike previous studies, we see no relationship between relative humidity and cloud occurrence biases, suggesting that the cloud biases do not result from the representation of moisture. We also show that the seasonal variation of cloud occurrence and cloud fraction, defined as the vertically integrated cloud occurrence, are small in both the observations and the reanalyses. We also examine the quality of the cloud representation for a set of weather states derived from ERA5 surface winds. The variability associated with grouping cloud occurrence based on weather state is much larger than the seasonal variation, highlighting weather state is a strong control of cloud occurrence. All the reanalyses continue to display underestimates below 3 km and overestimates above 3 km for each weather state. But the variability in ERA5 statistics matches the changes in the observations better than the other reanalyses. We also use a machine learning scheme to estimate the quantity of supercooled liquid water cloud from the ceilometer observations. Ceilometer low-level supercooled liquid water cloud occurrences are considerably larger than values derived from the reanalyses, further highlighting the poor representation of low-level clouds in the reanalyses.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Exposing and Reducing Biases of Simulating Mixed-Phase Clouds in the Convection-Permitting E3SM Atmosphere Model: Lessons From an Arctic Cold-Air Outbreak

Mixed-phase clouds modulate the water and energy cycles of high-latitude regions, yet their liquid-ice phase partitioning has long been poorly simulated in climate models. Here, simulations of Arctic mixed-phase clouds by the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) are assessed against large-eddy simulations, satellite data, and ground-based observations during the Cold-Air Outbreaks in the Marine Boundary Layer Experiment field campaign. SCREAM simulates nearly completely frozen clouds, which is attributed largely to the unreasonably strong Wegener–Bergeron–Findeisen (WBF) process that converts liquid to ice excessively and partly to the early over-abundant ice production at cold temperatures from a temperature-deterministic deposition ice nucleation scheme. Assuming no subgrid variation for the WBF process in the original formulation particularly conflicts with the instantaneous saturation adjustment assumption in the condensation scheme that assumes subgrid variability, leading to exaggerated WBF process rates. A proposed simple physically-based improvement on the treatment of subgrid cloud overlap substantially increases supercooled liquid water content and notably improves cloud-top phase partitioning, aligning better with observations. Improvement of supercooled liquid water content also converges with increasing horizontal resolution. The deposition ice nucleation scheme is found responsible for a falsely-produced ice cloud aloft that is not observed, biasing the simulated cloud radiative effects and top-of-atmosphere radiative fluxes. This study identifies key deficiencies in cloud parameterizations that continue to challenge convection-permitting models.

Geosciences↗
DOE OSTI2025

Primordial black holes from first-order phase transition in the singlet-extended SM

Supercooled first-order phase transition (FOPT) can lead to the formation of primordial black holes (PBHs). This scenario imposes stringent requirements on the profile of the effective potential. In this work, we use the singlet extended Standard Model (xSM) as a benchmark model to investigate this possibility at the electroweak scale. The PBHs formed during a supercooled FOPT have a narrow mass distribution around the mass of Earth. This distribution is closely tied to the temperature at which the PBHs form, corresponding to the FOPT at the electroweak scale. This scenario can be probed with microlensing experiments, space-based gravitational wave detectors, and collider experiments. Remarkably, the future space-based gravitational wave detector LISA will hold the potential to either confirm this scenario that leads to PBH formation across the observable Universe in the xSM or completely rule it out. Interestingly, our findings suggest that PBHs within the xSM framework may align with observations of the six ultrashort timescale events reported by the OGLE microlensing experiment. Published by the American Physical Society 2025

Gonçalves, Dorival (ORCID:0000000342234238)↗
DOE OSTI2025

Crystal nucleation rates in one-component Yukawa systems

Nucleation in the supercooled Yukawa system is relevant for addressing current challenges in understanding a range of crystallizing systems including white dwarf (WD) stars. We use both brute force and seeded molecular dynamics simulations to study homogeneous nucleation of crystals from supercooled Yukawa liquids. With our improved approach to seeded simulations, we obtain quantitative predictions of the crystal nucleation rate and cluster size distributions as a function of temperature and screening length. These quantitative results show trends towards fast nucleation with short-ranged potentials. They also indicate that for temperatures T > 0.9⁢T m , where T m is the melt temperature, classical homogeneous nucleation is too slow to initiate crystallization but transient clusters of ~100 particles should be common. As a result, we apply these general results to a typical WD model and obtain a delay of ~0.6 Gyr in the onset of crystallization that may be observable.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Glaciation of liquid clouds, snowfall, and reduced cloud cover at industrial aerosol hot spots

The ability of anthropogenic aerosols to freeze supercooled cloud droplets remains debated. Here, in this work, we present observational evidence for the glaciation of supercooled liquid-water clouds at industrial aerosol hot spots at temperatures between -10° and -24°C. Compared with the nearby liquid-water clouds, shortwave reflectance was reduced by 14% and longwave radiance was increased by 4% in the glaciation-affected regions. There was an 8% reduction in cloud cover and an 18% reduction in cloud optical thickness. Additionally, daily glaciation-induced snowfall accumulations reached 15 millimeters. Glaciation events downwind of industrial aerosol hot spots indicate that anthropogenic aerosols likely serve as ice-nucleating particles. However, rare glaciation events downwind of nuclear power plants indicate that factors other than aerosol emissions may also play a role in the observed glaciation events.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

High Energy Density Physics of Inertial Confinement Fusion Ablator Materials

The historic December 5, 2022 experiment at Lawrence Livermore National Lab’s (LLNL) National Ignition Facility (NIF) reached fusion energy ignition for the first time. This is the most important scientific breakthrough of the 21st century paves the way to future clean inertial fusion energy (IFE). The diamond (high density carbon (HDC)) ablator material used in this experiment displays detrimental effects due to the development of hydrodynamic instabilities at the diamond/fuel interface under shock compression. New alternatives to diamond ablators are required to step up the energy yield in ICF experiments. The unique combination of mechanical strength (approaching that of diamond), the ability to accommodate high-Z dopants (in contrast to diamond), and the tunability of the properties (through synthesis material with varying sp 3 content) make amorphous carbon (a-C) a promising material for next-generation IFE ablative capsules. However, despite its critical importance to the IFE program, the behavior of a-C carbon at extreme temperatures and pressures remains largely unexplored. The primary goals of this project were to perform groundbreaking dynamic compression experiments and predictive simulations to uncover the fundamental high-energy-density physics of amorphous carbon. Our goals were (1) to uncover the metastability range of amorphous carbon and probe phase transitions to diamond or metastable supercooled liquid carbon; (2) to acquire high-quality equation of state (EOS) data and develop an experimentally validated EOS from machine-learning MD simulations of the complex states of carbon; and (3) to uncover the complex behavior of carbon liquid in both thermodynamically stable and metastable supercooled states by accessing large areas of carbon phase diagram with amorphous samples with variable sp 3 content. Our proposed experimental program included measurements of equation of state and diffraction measurements using the Omega EP laser at the Laboratory of Laser Energetics at the University of Rochester. The theoretical/simulation program involved the development of machine-learning models of the complex response of amorphous carbon under dynamic compression by performing molecular dynamics simulations at experimental time and length scales using leadership class DOE supercomputers. Simulations guided experiments to observe predicted phenomena and acquire critical experimental data in specific pressure-temperature domains to validate theoretical models. This research delivered fundamental properties of novel amorphous carbon IFE ablator material, including phase diagram and EOS. These results will aid in IFE target design and implosion experiments. A unique combination of predictive simulations and dynamic and static experiments provided a highly inspirational intellectual environment for graduate students and postdocs involved in this project.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Thermodynamics of Liquid Uranium from Atomistic and Ab Initio Modeling

We present thermodynamic properties for liquid uranium obtained from classical molecular dynamics (MD) simulations and the first-principles theory. The coexisting phases method incorporated within MD modeling defines the melting temperature of uranium in good agreement with the experiment. The calculated melting enthalpy is in agreement with the experimental range. Classical MD simulations show that ionic contribution to the total specific heat of uranium does not depend on temperature. The density of states at the Fermi level, which is a crucial parameter in the determination of the electronic contribution to the total specific heat of liquid uranium, is calculated by ab initio all electron density functional theory (DFT) formalism applied to the atomic configurations generated by classical MD. The calculated specific heat of liquid uranium is compared with the previously calculated specific heat of solid γ-uranium at high temperatures. The liquid uranium cannot be supercooled below T sc ≈ 800 K or approximately about 645 K below the calculated melting point, although, the self-diffusion coefficient approaches zero at T D ≈ 700 K. Uranium metal can be supercooled about 1.5 times more than it can be overheated. The features of the temperature hysteresis are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The impact of cloud microphysics and ice nucleation on Southern Ocean clouds assessed with single-column modeling and instrument simulators

Abstract. Supercooled liquid clouds are common at higher latitudes (especially over the Southern Ocean) and are critical for constraining climate projections. We take advantage of the Macquarie Island Cloud and Radiation Experiment (MICRE) to perform an analysis of observed and simulated cloud processes over the Southern Ocean in a region and season dominated by supercooled liquid clouds. Using a single-column version of the European Centre for Medium-Range Weather Forecasts (ECMWF) Integrated Forecast System (IFS), we compare two different cloud microphysical schemes to ground-based observations of cloud, precipitation, and radiation over a 2.5-month period (1 January–17 March 2017). Both schemes are able to reproduce aspects of the cloud and radiation observations during MICRE to within the uncertainty of the data when the thermodynamic profile is prescribed with relaxation. There are differences in water mass and representation of reflectivity between the schemes. A sensitivity study of the cloud microphysics schemes, one a bulk one-moment scheme and the other a two-moment scheme with prediction of mass and number, indicates that several key processes create differences between the schemes. Surface radiative fluxes and total water path are highly sensitive to the formation and fall speed of precipitation. The prediction of hydrometeor number with the two-moment scheme yields a better comparison with observed reflectivity and radiative fluxes, despite predicting higher liquid water contents than observed. With the two-moment scheme, we are also able to test the sensitivity of the results to the input of liquid cloud condensation nuclei (CCN) and ice nuclei (IN). The cloud properties and resulting radiative effects are found to be sensitive to the CCN and IN concentrations. More CCN and IN increase liquid and ice water paths, respectively. Thus, both the dynamic environment and aerosols, integrated through the cloud microphysics, are important for properly representing Southern Ocean cloud radiative effects.

58 GEOSCIENCES↗
NASA NTRS1971

Fog dispersal technology.

The state-of-the-art in fog dispersal technology is briefly discussed. Fog is categorized as supercooled fog, occurring in air temperatures below freezing, and warm fog, occurring at above-freezing temperatures. Operational techniques are available to disperse supercooled fog in the airport area. It is much more difficult to cope with warm fog. Various known concepts to disperse warm fog are evaluated as to their operational merits. The most effective concept for immediate use involves heating the air to cause fog evaporation. Use of helicopter downwash has some application, possibly complementing the promising concept of seeding with sized hygroscopic particles. These latter two concepts appear to have future application, pending further research. The concept using polyelectrolytes is of uncertain value, lacking both a scientific explanation and a substantive evaluation of reported operational successes.

Mcgowan, W. A.↗
NASA NTRS1976

An experimental investigation of olivine morphology

Results are reported for a morphological study of olivine and an experimental investigation performed to determine the degrees of supercooling and the cooling rates necessary to crystallize particular morphologies. Ten arbitrary categories of three-dimensional olivine crystal shape are identified: polyhedral, granular, hopper, chain, lattice, plate, branching, radiate, feather, and swallow-tail. The morphological study establishes that equant and tabular crystals are the common shapes of olivine, nonequant crystals are elongate parallel to the a or c axis, and skeletal crystals result when a particular form is missing or only partially developed. In the experiment, olivine crystals were grown by melting rock samples above their liquidus temperatures before initiating crystallization. The results show that olivine morphology changes systematically as a function of the degree of melt supercooling, the melt cooling rate, and the normative olivine and water contents of the melt. It is also found that each shape has a specific range of temperature stability which is essentially independent of melt composition.

Donaldson, C. H.↗
NASA NTRS1975

Experimental modeling of the cooling history of Apollo 12 olivine basalts

An analog of the Apollo 12 olivine vitrophyres has been crystallized in a 1-atm gas-mixing furnace at cooling rates ranging between 1250 and 0.7 C/hr and isothermally at degrees of supercooling ranging from 10 C to 325 C. Mineral chemistry, crystal shapes, grain sizes and textures are systematically related to cooling rate and degree of supercooling. At linear cooling rates not exceeding 40 C/hr the texture is porphyritic - large olivine crystals are set in a groundmass of finer-grained pyroxene, plagioclase, and ilmenite; a later generation of olivine never crystallizes. There are three shapes of olivine crystals in the Apollo 12 olivine vitrophyres: glomerocrysts of subhedral crystals, large subequant skeletons, and highly elongate skeletons. These result from three generations of nucleation and a two-stage cooling history - a slow preeruption stage and a rapid continuously increasing posteruption stage. It seems likely that the Apollo 12 olivine basalt magmas were erupted with olivine crystals in suspension. The nucleation temperature of olivine in cooling-rate experiments is dependent on the experimental technique, and hence results of cooling experiments should be applied with caution.

Donaldson, C. H.↗
NASA NTRS1979

Convection and diffusion effects during dendritic solidification

A report is presented of the first quantitative measurements of dendritic growth at supercooling levels where convection instead of diffusion is the controlling heat transfer mechanism. Precautions similar to that used in an investigation conducted by Glicksman et al. (1976) were taken to insure 'free' dendritic growth conditions. Dendritic growth velocity was measured as a function of growth orientation at seventeen supercoolings which ranged from 0.043 C to 2 C. Selected but representative measurements of velocity versus orientation angle are shown in a graph. The relative growth velocity of a downward growing dendrite is found to be greater than that of a diffusion-limited dendrite. This result is consistent with that expected from the enhanced heat transfer arising from natural convection.

Glicksman, M. E.↗
NASA NTRS1978

Green glass vitrophyre 78526 - An impact of very low-Ti mare basalt composition

Rake sample 78526 is an 8.77 g rock consisting primarily of vitrophyric pale green glass with subordinate mineral and lithic relics. Petrographic and compositional evidence leads to the following conclusions: (1) the bulk composition represents that of a mixture formed by impact melting of at least two different textural and compositional varieties of VLT mare basalt that are now present in the rock as lithic relics and a poorly defined low-Ti mare basalt component observed in thin section only in the form of isolated mineral relics; (2) the admixed VLT mare basalts had REE abundances lower than those found in other mare basalts (but probably higher than emerald green glass) and REE patterns showing significant enrichment of the heavy relative to light REE's, suggesting that they were derived by comparatively high degrees of partial melting of a clinopyroxene-rich source region; and (3) the impact melt supercooled to produce the vitrophyre, with rather sharply contrasting textural domains present in the vitrophyre resulting from differences in nucleation kinetics and degrees of supercooling in various portions of the sample.

Warner, R. D.↗
NASA NTRS1980

Study of LPE methods for growth of InGaAsP/InP CW lasers

Two methods for liquid phase growth of InGaAsP/InP lasers were studied. Single phase growth, based on saturated melts and 5 C supercooling, was compared to two phase growth excess InP and 20 C nominal supercooling. Substrates cut on the (100) plane were used, and morphology in both cases was excellent and comparable to that obtainable in AlGaAs materials. A high degree of reproducibility was obtained in the materials grown by the two phased method, which is therefore presently preferred for the preparation of laser material. A refractive index step of 0.28 and an index n = 3.46 were obtained for In.81Ga.19As,5P5 lasing at 1.3 microns. Oxide-stripe lasers with typical room temperature cw threshold currents of 180 mA were obtained and some of them showed single mode behavior without lateral cavity modifications. COntinuous operation of 800 h at room temperature was obtained without noticeable degradation.

Ladany, I.↗
NASA NTRS1981

Homogeneous condensation - Freezing nucleation rate measurements for small water droplets in an expansion cloud chamber

Experimental data on ice nucleation, presented in an earlier paper, are analyzed to yield information about the homogeneous nucleation rate of ice from supercooled liquid and the heights of energy barriers to that nucleation. The experiment consisted of using an expansion cloud chamber to nucleate from the vapor a cloud of supercooled pure water drops and the observation of the fraction of drops which subsequently froze. The analysis employed standard classical homogeneous nucleation theory. The data are used to extract the first experimental measurement (albeit indirect) of the activation energy for the transfer of a water molecule across the liquid-ice interface at temperatures near -40 C. The results provide further evidence that the local liquid structure becomes more icelike as the temperature is lowered.

Hagen, D. E.↗
NASA NTRS1982

The interaction of radio frequency electromagnetic fields with atmospheric water droplets and applications to aircraft ice prevention

The feasibility of computerized simulation of the physics of advanced microwave anti-icing systems, which preheat impinging supercooled water droplets prior to impact, was investigated. Theoretical and experimental work performed to create a physically realistic simulation is described. The behavior of the absorption cross section for melting ice particles was measured by a resonant cavity technique and found to agree with theoretical predictions. Values of the dielectric parameters of supercooled water were measured by a similar technique at lambda = 2.82 cm down to -17 C. The hydrodynamic behavior of accelerated water droplets was studied photograhically in a wind tunnel. Droplets were found to initially deform as oblate spheroids and to eventually become unstable and break up in Bessel function modes for large values of acceleration or droplet size. This confirms the theory as to the maximum stable droplet size in the atmosphere. A computer code which predicts droplet trajectories in an arbitrary flow field was written and confirmed experimentally. The results were consolidated into a simulation to study the heating by electromagnetic fields of droplets impinging onto an object such as an airfoil. It was determined that there is sufficient time to heat droplets prior to impact for typical parameter values. Design curves for such a system are presented.

Hansman, R. J., Jr.↗
NASA NTRS1979

Influence of convection on free growth of dendrite crystals from solution

The free growth of dendrites in a uniformly supercooled solution was examined using cine photography with a Schlieren optical system. Crystals were grown in the bulk of the solution from a centrally located capillary tube, nucleated at the interface with a liquid nitrogen cooled wire. Crystals propagated along the tube, the slower growing orientations eliminated, and emerged at the tip, usually growing parallel to the tube direction. For both sodium sulfate decahydrate from its solution and ice from sodium chloride solution, growth rate and fineness of dendrites increased with supercooling. In sodium sulfate, upward convection of the less dense depleted solution occurs; downward convection was observed for the rejected, more concentrated sodium chloride solution. In both cases, there was a spatial and temporal delay in the release of the convective plume from the moving dendrite tip. The role of this convection on the growth characteristics and the production of secondary crystals is examined. A proposed low-g experiment to examine differences in growth rate, crystal texture, and secondary nucleation in a reduced convective regime where molecular diffusion is the dominant transfer process is discussed.

Hallett, J.↗