Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal equilibrium”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Predictive Modeling of Carbon Ablators Using Micro and Macro-Scale Modeling

Efforts to build a Predictive Material Modeling (PMM) framework from the micro-scale to the macro-scale are presented in this abstract. To reduce the need for extensive testing, accelerate the design cycle process, and reduce uncertainty margins applied to final designs, NASA is developing simulation and modeling tools that enable characterization of material properties and response to high-enthalpy environments. The Porous Microstructure Analysis (PuMA) code has been developed for computing macroscale (volume averaged) properties of porous materials using microscale images from micro-computed tomography (micro-CT). Microscale modeling requires a realistic representation of a material microstructure; these are obtained either synthetically during the design of the material or through X-ray micro-CT. Volume averaged properties are then used to inform macroscale material response models, such as those implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO) software, also actively developed by NASA. The computational model in PATO is a generic heat and mass transfer model for porous reactive materials containing several solid phases and a single gas phase. The detailed chemical interactions occurring between the solid phases and the gas phase are modeled at the pore scale assuming local thermal equilibrium. These tools were developed to efficiently interface with other pre-existing codes such as SPARTA (direct simulation Monte Carlo), DPLR (hypersonic CFD), NEQAIR (radiative transport) and DAKOTA (uncertainty quantification and optimization). Detailed flight data (Mars Science Laboratory [MSL] Entry Descent and Landing Instrument [MEDLI]) is critical for validating these computational tools for NASA applications. Examples of modeling ablative material response using these codes will be presented including 3D simulations of the full-scale heatshield of the MSL capsule. The simulations demonstrate the ability of the modern material response code, PATO, to handle the material response of geometrically complex and large domains, through the use of massively parallel computations.

Thermal Protection Systems↗
NASA NTRSDate not supplied

Predictive Modeling of Carbon Ablators

Efforts to build a Predictive Material Modeling (PMM) framework from the micro-scale to the macro-scale are presented in this abstract. To reduce the need for extensive testing, accelerate the design cycle process, and reduce uncertainty margins applied to final designs, NASA is developing simulation and modeling tools that enable characterization of material properties and response to high-enthalpy environments. The Porous Microstructure Analysis (PuMA) code has been developed for computing macroscale (volume averaged) properties of porous materials using microscale images from micro-computed tomography (micro-CT). Microscale modeling requires a realistic representation of a material microstructure; these are obtained either synthetically during the design of the material or through X-ray micro-CT. Volume averaged properties are then used to inform macroscale material response models, such as those implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO) software, also actively developed by NASA. The computational model in PATO is a generic heat and mass transfer model for porous reactive materials containing several solid phases and a single gas phase. The detailed chemical interactions occurring between the solid phases and the gas phase are modeled at the pore scale assuming local thermal equilibrium. These tools were developed to efficiently interface with other pre-existing codes such as SPARTA (direct simulation Monte Carlo), DPLR (hypersonic CFD), NEQAIR (radiative transport) and DAKOTA (uncertainty quantification and optimization). Detailed flight data (Mars Science Laboratory [MSL] Entry Descent and Landing Instrument [MEDLI]) is critical for validating these computational tools for NASA applications. Examples of modeling ablative material response using these codes will be presented including 3D simulations of the full-scale heatshield of the MSL capsule. The simulations demonstrate the ability of the modern material response code, PATO, to handle the material response of geometrically complex and large domains, through the use of massively parallel computations.

Thermal Protection Systems↗
NASA NTRS1999

Abundant Cool Magnesium-Rich Pyroxene Crystals in Comet Hale-Bopp

Modeling of the observed dust emission from Comet Hale-Bopp over a large range of heliocentric distances (2.8 AU - 0.93 AU -1.7 AU) led to the discovery of Mg-rich pyroxene crystals in the coma These pyroxene crystals are apparent in the 10 micron spectrum only when the comet is close to perihelion (r(sub h) = 1.2 AU) because they are cooler than the other silicate minerals. The pyroxene crystals are cooler than the other species because of their high Mg-content. They do not absorb as efficiently as the other silicate minerals. Given the same high Mg content of Mg/(Mg + Fe) = 0.9, radiative equilibrium computations show that pyroxene crystals are expected to be 150 K cooler than olivine crystals. The pyroxene crystals are also about 10x more abundant than the other silicate mineral species. Their high Mg content and relatively large abundance are in agreement with the preponderance of pyroxene interplanetary dust particles (IDPs) and the recent reanalysis of the PUMA-I flyby of Comet Halley. Before Hale-Bopp, only olivine crystals were detected spectroscopically in comets, probably because the pyroxene crystals are less optically active, hence significantly cooler and harder to detect in contrast to the warmer silicate species. Determining the relative abundances of silicate minerals depends on their Mg contents. If the pyroxene crystals in Comet Hale-Bopp are solar nebula condensates, then they probably had to form during the early FU Orionis epoch when the inner disk was hot enough and be transported out to the region of formation of icy planetesimals without being reheated. Reheating events appear to reincorporate Fe back into the crystals or form Fe-rich rims, which are not consistent with the high-Mg-content crystals. The condensation of Mg-rich pyroxene crystals is expected from solar nebula thermal equilibrium computations. However, their subsequent transport to the outer solar nebula unaltered has yet to be theoretically demonstrated. The discovery of Mg-rich crystals in Comet Hale-Bopp and in AGB stars opens the possibility that these crystals are relic interstellar grains. One-third of IDPs have been shown to have significant deuterium enrichments, thus indicating that they contain presolar material. By spectroscopic analogy to IDPs the Mg-rich pyroxene crystals in Comet Hale-Bopp may be presolar grains. If so, then the comet contains largely ISM silicates. ISM grains may have been the dominant source of dust in the outer early solar nebula.

Wooden, D. H.↗
DOE OSTI2022

Numerical Investigation of the Spark Discharge Process in a Crossflow

The present study numerically investigates the spark discharge process under crossflow conditions using a thermal equilibrium plasma solver that fully couples the electromagnetic physics and fluid dynamics in a computational framework. Numerical results are validated by the comparison with experimental data. Here, the spark discharge experiment is performed in a constant volume vessel using an inductive coil ignition system for automotive applications, and the evolution of the spark channel is measured using high-speed imaging. The crossflow in the gap between the spark-plug electrodes is generated by a rotating fan with two different fan speeds, and the flow velocity across the gap is characterized by particle image velocimetry (PIV) measurement. A computational fluid dynamics (CFD) solver is employed to simulate the crossflow and provide the flow field variables (velocity, pressure, temperature) to the plasma solver. The crossflow velocity predicted in the flow simulation agrees well with the PIV data in that the non-uniform velocity profiles at monitoring points are reproduced by the CFD code. With the crossflow initialization in the plasma solver, the simulated spark discharge process from the breakdown to spark discharge matches the experimental data, including the voltage and circuit waveforms and the high-speed images of the spark channel evolution. The stretch of spark channel captured by plasma simulations agrees with the measured data. The plasma simulation reveals that the mean temperature of the spark channel is maintained at 5000 K during the discharge phase, and the temperature varies along the spark channel so that the highest value is obtained at the spark root on the center electrode. Overall, the results presented in this paper are meant to provide valuable information about the properties of the plasma generated by the spark discharge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Ultrafast Carrier and Lattice Dynamics in Plasmonic Nanocrystalline Copper Sulfide Films

Excited carrier dynamics in plasmonic nanostructures determines many important optical properties such as nonlinear optical response and photocatalytic activity. Here it is shown that mesoscopic plasmonic covellite nanocrystals with low free-carrier concentration exhibit a much faster carrier relaxation than in traditional plasmonic materials. A nonequilibrium hot-carrier population thermalizes within first 20 fs after photoexcitation. A decreased thermalization time in nanocrystals compared to a bulk covellite is consistent with the reduced Coulomb screening in ultrathin films. The subsequent relaxation of thermalized, equilibrium electron gas is faster than in traditional plasmonic metals due to the lower carrier concentration and agrees well with that in a bulk covellite showing no evidence of quantum confinement or hot-hole trapping at the surface states. The excitation of coherent optical phonon modes in a covellite is also demonstrated, revealing coherent lattice dynamics in plasmonic materials, which until now was mainly limited to dielectrics, semiconductors, and semimetals. These findings show advantages of this new mesoscopic plasmonic material for active control of optical processes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Islands in the fluid: islands are common in cosmology

We discuss the possibility of entanglement islands in cosmological spacetimes with a general perfect fluid with an equation of state w. We find that flat universes with time-symmetric slices where the Hubble parameter vanishes always have islands on that slice. We then move away from such slices, considering still universes with a general perfect fluid. Under the local thermal equilibrium assumption, the comoving entropy density s c is constant. As a result, the conditions for an island become an inequality between the energy density (or Hubble parameter) and the temperature at some time of normalization. The consequences are that islands can exist for practically all fluids that are not radiation, i.e. w ≠ 1/3. We also discuss the ramifications of our results for universes with spatial curvature. Finally, we show that islands occur in the Simple Harmonic Universe model which has no classical singularity at the background level, in contrast to all previous examples where islands occurred only in space-times with singularities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The price of curiosity: information recovery in de Sitter space

Recent works have revealed that quantum extremal islands can contribute to the fine-grained entropy of black hole radiation reproducing the unitary Page curve. In this paper, we use these results to assess if an observer in de Sitter space can decode information hidden behind their cosmological horizon. By computing the fine-grained entropy of the Gibbons-Hawking radiation in a region where gravity is weak we find that this is possible, but the observer’s curiosity comes at a price. At the same time the island appears, which happens much earlier than the Page time, a singularity forms which the observer will eventually hit. We arrive at this conclusion by studying de Sitter space in Jackiw-Teitelboim gravity. We emphasize the role of the observer collecting radiation, breaking the thermal equilibrium studied so far in the literature. By analytically solving for the backreacted geometry we show how an island appears in this out-of-equilibrium state.

2D gravity↗
DOE OSTI2020

Sterile neutrino dark matter in left-right theories

SU(2) L × SU(2) R gauge symmetry requires three right-handed neutrinos (N i ), one of which, N 1 , can be sufficiently stable to be dark matter. In the early universe, W R exchange with the Standard Model thermal bath keeps the right-handed neutrinos in thermal equilibrium at high temperatures. N 1 can make up all of dark matter if they freeze-out while relativistic and are mildly diluted by subsequent decays of a long-lived and heavier right-handed neutrino, N 2 . We systematically study this parameter space, constraining the symmetry breaking scale of SU(2) R and the mass of N 1 to a triangle in the (v R , M 1 ) plane, with v R = (10 6 - 3 × 10 12 ) GeV and M 1 = (2 keV–1 MeV). Much of this triangle can be probed by signals of warm dark matter, especially if leptogenesis from N 2 decay yields the observed baryon asymmetry. The minimal value of v R is increased to 10 8 GeV for doublet breaking of SU(2) R , and further to 10 9 GeV if leptogenesis occurs via N 2 decay, while the upper bound on M 1 is reduced to 100 keV. In addition, there is a component of hot N 1 dark matter resulting from the late decay of N 2 → N 1 ℓ+ℓ- that can be probed by future cosmic microwave background observations. Interestingly, the range of v R allows both precision gauge coupling unification and the Higgs Parity understanding of the vanishing of the Standard Model Higgs quartic at scale vR. Finally, we study freeze-in production of N 1 dark matter via the WR interaction, which allows a much wider range of (v R , M 1 ).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Black hole induced false vacuum decay: the role of greybody factors

We study false vacuum decay catalyzed by black holes. We consider a toy two-dimensional model of a scalar field with an unstable potential in the background of a dilaton black hole. A realistic black hole in four dimensions possesses the potential barrier for linear field perturbations. We model this barrier — the greybody factor — for spherically-symmetric perturbations in the toy model by adding a coupling between the scalar field and dilaton. We compute analytically the decay rate for the black hole in thermal equilibrium (Hartle-Hawking state) and for the radiating black hole in empty space (Unruh state). Our results show that, contrary to the Hartle-Hawking vacuum, the decay probability of the Unruh vacuum remains exponentially suppressed at all black hole temperatures. We argue that this result holds also in four dimensions.

2D Gravity↗
DOE OSTI2025

The discriminant power of bubble wall velocities: gravitational waves and electroweak baryogenesis

A precise determination of the bubble wall velocity v$_{w}$ is crucial for making accurate predictions of the baryon asymmetry and gravitational wave (GW) signals in models of electroweak baryogenesis (EWBG). Working in the local thermal equilibrium approximation, we exploit entropy conservation to present efficient algorithms for computing v$_{w}$, significantly streamlining the calculation. We then explore the parameter dependencies of v$_{w}$, focusing on two sample models capable of enabling a strong first-order electroweak phase transition: a ℤ$_{2}$-symmetric singlet extension of the SM, and a model for baryogenesis with CP violation in the dark sector. We study correlations among v$_{w}$ and the two common measures of phase transition strength, α$_{n}$ and v$_{n}$/T$_{n}$. Interestingly, we find a relatively model-insensitive relationship between v$_{n}$/T$_{n}$ and α$_{n}$. We also observe an upper bound on α$_{n}$ for the deflagration/hybrid wall profiles naturally compatible with EWBG, the exact value for which varies between models, significantly impacting the strength of the GW signals. In summary, our work provides a framework for exploring the feasibility of EWBG models in light of future GW signals.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Submerging islands through thermalization

We illustrate scenarios in which Hawking radiation collected in finite regions of a reservoir provides temporary access to the interior of black holes through transient entanglement “islands.” Whether these islands appear and the amount of time for which they dominate — sometimes giving way to a thermalization transition — is controlled by the amount of radiation we probe. In the first scenario, two reservoirs are coupled to an eternal black hole. The second scenario involves two holographic quantum gravitating systems at different temperatures interacting through a Rindler-like reservoir, which acts as a heat engine maintaining thermal equilibrium. The latter situation, which has an intricate phase structure, describes two eternal black holes radiating into each other through a shared reservoir.

2D gravity↗
DOE OSTI2020

Freeze-in dark matter from secret neutrino interactions

We investigate a simplified freeze-in dark-matter model in which the dark matter only interacts with the standard-model neutrinos via a light scalar. The extremely small coupling for the freeze-in mechanism is naturally realized in several neutrino-portal scenarios with the secret neutrino interactions. We study possible evolution history of the hidden sector: the dark sector would undergo pure freeze-in production if the interactions between the dark-sector particles are negligible, while thermal equilibrium within the dark sector could occur if the reannihilation of the dark matter and the scalar mediator is rapid enough. We investigate the relic abundance in the freeze-in and dark freeze-out regimes, calculate evolution of the dark temperature, and study its phenomenological aspects on BBN and CMB constraints, the indirect-detection signature, as well as the potential to solve the small scale structure problem.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Effect of laser melt schedule on the microstructure of additively manufactured IN718 Superalloy

Laser powder bed fusion (L-PBF) has enabled the fabrication of geometrically complex metallic structures and components that are challenging to producing using conventional manufacturing approaches. The site-specific and far from equilibrium thermal conditions of L-PBF offer the potential to facilitate multi-length scale design of structure and properties across the atomic-through macro-levels. However, L-PBF systems face scalability challenges due to throughput constraints. Laser rotary powder bed fusion (L-RPBF) systems are being investigated as a solution to enhance the deposition rates compared to conventional L-PBF. Rotary systems also offer additional flexibility for controlling the time structure of melting through laser interleaving on alternating layers. Here, in this study, IN718 test samples were printed using single-laser or interleaved dual-laser configuration in a L-RPBF system to investigates the effect of process settings and melt-interleaving on as-fabricated microstructure. The microstructural evolution, such as grain size and crystallographic texture, was assessed by determining variations in the melt-pool shapes. Laser interleaving leads to a reduction in average grain size compared to single laser by ∼ 40 % at high power (400 W) and by ∼36 % at medium power (370 W). Results presented here identify key challenge for obtaining uniform microstructures and barriers for the broader adoption of high-deposition rate L-RPBF.

Dual-laser↗
DOE OSTI2024

Catalytic production of δ-valerolactone (DVL) from biobased 2-hydroxytetrahydropyran (HTHP) – Combined experimental and modeling study

δ-Valerolactone (DVL) is a five-carbon (C5) cyclic ester that can undergo ring-opening polymerization to yield high-performance, biocompatible polyesters. But current market prices of C5 chemicals like DVL are very high due to poor availability of C5 feedstock in petroleum. Herein, we demonstrate a novel route to DVL synthesis via dehydrogenation of biomass-derived 2-hydroxytetrahydropyran (HTHP) over Cu/SiO 2 without the use of toxic reagents. Since HTHP exists in thermal equilibrium with 3,4-dihydropyran (DHP) via dehydration, and with 2,2’-oxybis(tetrahydropyran) and 5-(tetrahydropyran-2-yloxy)pentanal via acetalization, we have also determined the thermochemistry (ΔH rxn and ΔG rxn ) of each competing reaction using density functional theory (DFT) calculations at the M06-2X/cc-pVTZ level. Finally, by developing a kinetic model of all 8 reactions involved, we have achieved 84% selectivity to DVL at 150°C in a packed bed reactor for over 72 hours of time on stream.

42 ENGINEERING↗
DOE OSTI2022

Silicon alloying enhances fast heating rate combustion of aluminum particles

Alloys and composites of aluminum (Al) have shown promise in regulating and enhancing particle combustion for energy generation applications. Recent work on aluminum-silicon (Al-Si) spherical alloy particles has shown improved combustion at low heating rates through enhanced diffusion accompanying a lower melting temperature. This study extends reactivity analysis to higher heating rates, comparing oxidation of Al-Si with Al. Flame speeds of Al-Si powder mixed with molybdenum trioxide (MoO 3 ) powder (Al-Si+MoO 3 ) exhibited a faster transition to steady propagation relative to Al+MoO 3 . Bomb calorimetry experiments revealed up to 5.8% greater early temperature rise for Al-Si powder. Rapid steady propagation for Al-Si mixture and faster temperature rise for Al-Si particles were attributed to accelerated kinetics evidenced in thermal equilibrium analysis of the mixtures using a differential scanning calorimeter (DSC). Larger DSC exotherms in the early stages of oxidation (i.e., 480–720 °C) correlated with early heat release and promoted steady flame propagation for Al-Si+MoO 3 compared with Al+MoO 3 . Furthermore, in multiple heating rate DSC studies, Al-Si+MoO 3 ignited at a lower heating rate (i.e., 15 ⁰C/min) than Al+MoO 3 (i.e., 20 ⁰C/min). Both mixtures ignited at temperatures less than 625 °C which is above the melting temperature of Al-Si (574 °C) and below Al (659 °C). Thus, alloy particle fuels with enhanced diffusion-controlled kinetics promote steady flame propagation and show promise for energy generation applications. This is especially promising for technologies driven by enhanced condensed phase combustion such as some propellant additives or applications dependent on reliable burn rates such as pyrotechnic time delay formulations and primers.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A fully implicit, scalable, conservative nonlinear relativistic Fokker–Planck 0D-2P solver for runaway electrons

Upon application of a sufficiently strong electric field, electrons break away from thermal equilibrium and approach relativistic speeds. These highly energetic ‘runaway’ electrons (~ MeV) play a significant role in tokamak disruption physics, and therefore their accurate understanding is essential to develop reliable mitigation strategies. As such, we have developed a fully implicit solver for the 0D-2P (i.e., including two momenta coordinates) relativistic nonlinear Fokker–Planck equation (rFP). As in earlier implicit rFP studies (NORSE, CQL3D), electron–ion interactions are modeled using the Lorentz operator, and synchrotron damping using the Abraham–Lorentz–Dirac reaction term. However, our implementation improves on these earlier studies by (1) ensuring exact conservation properties for electron collisions, (2) strictly preserving positivity, and (3) being scalable algorithmically and in parallel. Key to our proposed approach is an efficient multigrid preconditioner for the linearized rFP equation, a multigrid elliptic solver for the Braams–Karney potentials, and a novel adaptive technique to determine the associated boundary values. We verify the accuracy and efficiency of the proposed scheme with numerical results ranging from small electric-field electrical conductivity measurements to the accurate reproduction of runaway tail dynamics when strong electric fields are applied.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Thermodynamic modeling of MDEA(aq)-NH 3 (aq)-K 2 CO 3 (aq)-CO 2 (aq) using the Extended UNIQUAC model

Here, the development of advanced processes for carbon dioxide capture requires reliable thermodynamic models to ensure accurate process design. This work presents an optimized parameter set of the Extended UNIQUAC model for advanced mixed solvent solutions containing CO 2 -MDEA-NH 3 -K 2 CO 3 in water. The model was parameterized using vapor–liquid equilibrium, pure component saturation pressure, excess enthalpy, heat of absorption, molar heat capacity, apparent molar heat capacity, and solid–liquid equilibrium. These thermal and equilibrium properties cover temperatures between 253.6 K to 600 K and pressures up to 75.6 bar. The current model represents an extension of previous versions that were suitable for modeling aqueous NH 3 and K 2 CO 3 solutions. This parametrization accounts for molecular interactions between MDEA and these two components, which allows for the design of advanced capture processes. This paper demonstrates that the new model parameters are valid within the range of operational conditions typically found in CO 2 capture processes.

42 ENGINEERING↗
DOE OSTI2022

The bands method for tabulating NLTE material properties

The use of Non-Local Thermal Equilibrium (NLTE) material properties within radiation-hydrodynamic simulations remains a major challenge. The plasma radiative and EOS properties in NLTE can depend on the detailed local radiation field and electron energy distribution. Fully characterizing each of these quantities may require 10’s to 100’s of parameters, making tabulation effectively impossible and requiring expensive calculations for each set of conditions within the simulation. In this work we present a new method that characterizes the local radiation field with a limited number of parameters. Furthermore, this permits pre-calculation and tabulation of the material properties, allowing codes to perform NLTE radiation-hydrodynamic simulations which would otherwise be computationally prohibitive.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗