Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Temperature model”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Thermally driven flow in a rotating spherical shell Axisymmetic states

A spherical analogue of the rotating annulus experiments modeling atmospheric motion, in which a liquid is contained between two rigid, corotating and concentric hemispheres upon both of which thermal gradients are imposed, is presently studied by means of numerical models. Temperatures are lower on the inner than on the outer sphere, and decrease towards the pole. Using Navier-Stokes equations which assume symmetry about the polar axis, finite difference numerical models yield steady-state solutions to the equations. Hydrostatic and nonhydrostatic solutions are compared for cylindrical and spherical cases, and it is found in the case of the spherical shell that the differences between hydrostatic and nonhydrostatic solutions are small and largely confined to the regions near the pole and equator. It is suggested that nonhydrostatic effects on the axisymmetric state will not affect the flow's baroclinic stability.

Miller, T. L.↗
DOE OSTI2013

Phase field simulation of kinetic superheating and melting of aluminum nanolayer irradiated by pico- and femtosecond laser

Two melting mechanisms are reproduced and quantified for superheating and melting of Al nanolayer irradiated by pico- and femtosecond laser using the advanced phase-field approach coupled with mechanics and a two-temperature model. At heating rates Q≤79.04 K/ps induced by picosecond laser, two-sided barrierless surface melting forms two solid-melt interfaces, which meet near the center of a sample. The temperature for surface melting is a linear function, and for complete melting it is a cubic function, of logQ. At Q≥300 K/ps induced by femtosecond laser, barrierless and homogeneous melting (without nucleation) at the sample center occurs faster than due to interface propagation. Good agreement with experimental melting time was achieved in a range of 0.95≤Q≤1290 K/ps without fitting of material parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Role of electronic thermal transport in amorphous metal recrystallization: A molecular dynamics study

Recrystallization of glasses is important in a wide range of applications including electronics and reactive materials. Molecular dynamics (MD) has been used to provide an atomic picture of this process, but prior work has neglected the thermal transport role of electrons, the dominant thermal carrier in metallic systems. We characterize the role of electronic thermal conductivity on the velocity of recrystallization in Ni using MD coupled to a continuum description of electronic thermal transport via a two-temperature model. Our simulations show that for strong enough coupling between electrons and ions, the increased thermal conductivity removes the heat from the exothermic recrystallization process more efficiently, leading to a lower effective temperature at the recrystallization front and, consequently, lower propagation velocity. We characterize how electron-phonon coupling strength and system size affect front propagation velocity. Interestingly, we find that initial recrystallization velocity increases with decreasing system size due to higher overall temperatures. Overall, we show that a more accurate description of thermal transport due to the incorporation of electrons results in better agreement with experiments.

McClure, Zachary D. (ORCID:0000000183300217)↗
NASA NTRS1998

Flow Characteristics of UHV Nitrogen Arcjet

Flow characteristics in the nozzle of a nitrogen arcjet have been simulated by a model. Electrons, ions, and neutral atoms and molecules are represented as separate species. Dissociation, ionization, and recombination are treated as separate reactions. Thermal non-equilibrium is represented by a two-temperature model. Energy input to plasma from the arc is modeled as a source determined by local current and electrical conductivity. Momentum and energy losses and recombinations at the nozzle wall are included as source/sink determined by fluxes to the wall. The results show essentially frozen gas-phase reactions and thermal nonequilibrium due to Joule heating of electrons and in the expansion part of the nozzle. Wall interactions have strong effects on the results, indicating that they play important roles in the flow due to very high surface area relative to the volume of the plasma. For example, plasma velocity decreases from 6000 m per s without wall interactions to 3000 m per s with modest amount of wall interactions. These results also show reasonable agreement with optical emission measurements, which confirms that the arcjet plasma is far from LTE. The spectra suggest nitrogen dissociation levels of 0.3% - 9%, depending on nitrogen flow rate and arc plasma current, which also determine the relative amounts of excited atomic and molecular nitrogen. Langmuir probe studies of the source show that electron and ion fluxes increase with increasing power, and that the ion energy distribution shifts to lower energies. Typical ion fluxes were on the order of 4E-9 A/cm**2 with a maximum ion kinetic energy of 3.5eV. The median electron energy was 1eV, with a maximum of less than 4eV.

Chang, C. H.↗
DOE OSTI2016

VERA-CS Modeling and Simulation of PWR Main Steam Line Break Core Response to DNB

The Virtual Environment for Reactor Applications core simulator (VERA-CS) being developed by the Consortium for the Advanced Simulation of Light Water Reactors (CASL) includes coupled neutronics, thermal-hydraulics, and fuel temperature components with an isotopic depletion capability. The neutronics capability employed is based on MPACT, a three-dimensional (3-D) whole core transport code. The thermal-hydraulics and fuel temperature models are provided by the COBRA-TF (CTF) subchannel code. As part of the CASL development program, the VERA-CS (MPACT/CTF) code system was applied to model and simulate reactor core response with respect to departure from nucleate boiling ratio (DNBR) at the limiting time step of a postulated pressurized water reactor (PWR) main steamline break (MSLB) event initiated at the hot zero power (HZP), either with offsite power available and the reactor coolant pumps in operation (high-flow case) or without offsite power where the reactor core is cooled through natural circulation (low-flow case). The VERA-CS simulation was based on core boundary conditions from the RETRAN-02 system transient calculations and STAR-CCM+ computational fluid dynamics (CFD) core inlet distribution calculations. The evaluation indicated that the VERA-CS code system is capable of modeling and simulating quasi-steady state reactor core response under the steamline break (SLB) accident condition, the results are insensitive to uncertainties in the inlet flow distributions from the CFD simulations, and the high-flow case is more DNB limiting than the low-flow case.

CTF↗
DOE OSTI2014

Modelling morphology evolution during solidification of IPP in processing conditions

During polymer processing, crystallization takes place during or soon after flow. In most of cases, the flow field dramatically influences both the crystallization kinetics and the crystal morphology. On their turn, crystallinity and morphology affect product properties. Consequently, in the last decade, researchers tried to identify the main parameters determining crystallinity and morphology evolution during solidification In processing conditions. In this work, we present an approach to model flow-induced crystallization with the aim of predicting the morphology after processing. The approach is based on: interpretation of the FIC as the effect of molecular stretch on the thermodynamic crystallization temperature; modeling the molecular stretch evolution by means of a model simple and easy to be implemented in polymer processing simulation codes; identification of the effect of flow on nucleation density and spherulites growth rate by means of simple experiments; determination of the condition under which fibers form instead of spherulites. Model predictions reproduce most of the features of final morphology observed in the samples after solidification.

36 MATERIALS SCIENCE↗
NASA NTRS1975

The NATA code: Theory and analysis, volume 1

A computer program for calculating quasi-one-dimensional gas flow in axisymmetric and two-dimensional nozzles and rectangular channels is presented. Flow is assumed to start from a state of thermochemical equilibrium at a high temperature in an upstream reservoir. The program provides solutions based on frozen chemistry, chemical equilibrium, and nonequilibrium flow with finite reaction rates. Electronic nonequilibrium effects can be included using a two-temperature model. An approximate laminar boundary layer calculation is given for the shear and heat flux on the nozzle wall. Boundary layer displacement effects on the inviscid flow are considered also. Chemical equilibrium and transport property calculations are provided by subroutines. The code contains precoded thermochemical, chemical kinetic, and transport cross section data for high-temperature air, CO2-N2-Ar mixtures, helium, and argon. It provides calculations of the stagnation conditions on axisymmetric or two-dimensional models, and of the conditions on the flat surface of a blunt wedge. The primary purpose of the code is to describe the flow conditions and test conditions in electric arc heated wind tunnels.

Bade, W. L.↗
DOE OSTI2022

Evaluation of a collisional radiative model for electron temperature determination in hydrogen plasma

A collisional-radiative (CR) model that extracts the electron temperature, T e , of hydrogen plasmas from Balmer-line-ratio measurements is examined for the plasma electron density, n e , and T e ranges of 10 10 –10 15 cm –3 and 5–500 eV, respectively. The CR code, developed and implemented in Python, has a forward component that computes the densities of excited states up to n = 15 as functions of T e , n e , and the molecular-to-atomic neutral ratio r(H 2 /H). The backward component provides ne and r(H 2 /H) as functions of the Balmer ratios to predict the T e . The model assumes Maxwellian electrons. Furthermore, the density profiles of the electrons and of the molecular and atomic hydrogen neutrals are shown to be of great importance, as is the accuracy of the line-ratio measurement method.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Uncertainty in Reconstructing Paleo‐Elevation of the Antarctic Ice Sheet From Temperature‐Sensitive Ice Core Records

Paleotemperature reconstructions from ice cores are mixed signals of changes in climate and ice‐surface elevation. A common, temperature‐based paleoaltimetry method suggests these signals can be disentangled by comparing two proxy locations with similar climates. The difference between the records is assumed to be due to elevation, which is estimated by scaling the temperature difference by a lapse rate. We investigate the uncertainty associated with this approach using a case study of the Antarctic Ice Sheet during the Last Glacial Maximum. From an ensemble of climate simulations, we extract modeled temperatures at locations of real ice cores. We find uncertainty on the order of hundreds of meters that results from spatial heterogeneity in non‐adiabatic temperature change, which itself stems in part from elevation‐induced atmospheric circulation change. Our findings suggest that caution is needed when interpreting temperature‐based paleoaltimetry results for ice sheets.

58 GEOSCIENCES↗
DOE OSTI2022

Equilibrium reconstruction of DIII-D plasmas using predictive modeling of the pressure profile

New workflows have been developed for predictive modeling of magnetohydrodynamic (MHD) equilibrium in tokamak plasmas. The goal of this work is to predict the MHD equilibrium in tokamak discharges without having measurements of the kinetic profiles. The workflows include a cold start tool, which constructs all the profiles and power flows needed by transport codes; a Grad–Shafranov equilibrium solver; and various codes for the sources and sinks. For validation purposes, a database of DIII-D tokamak discharges has been constructed that is comprised of scans in the plasma current, toroidal magnetic field, and triangularity. Initial efforts focused on developing a workflow utilizing an empirically derived pressure model tuned to DIII-D discharges with monotonic safety factor profiles. This workflow shows good agreement with experimental kinetic equilibrium calculations, but is limited in that it is a single fluid (equal ion and electron temperatures) model and lacks H-mode pedestal predictions. The best agreement with the H-mode database is obtained using a theory-based workflow utilizing pressure profile predictions from a coupled TGLF turbulent transport and EPED pedestal models together with external magnetics and Motional Stark Effect (MSE) data to construct the equilibrium. Here, we obtain an average root mean square error of 5.1% in the safety factor profile when comparing the predicted and experimental kinetic equilibrium. We also find good agreement with the plasma stored energy, internal inductance, and pressure profiles. Including MSE data in the theory-based workflow results in noticeably improved agreement with the q-profiles in high triangularity discharges in comparison with the results obtained with magnetic data only. The predictive equilibrium workflow is expected to have wide applications in experimental planning, between-shot analysis, and reactor studies.

Kinsey, J. E. (ORCID:0000000193347473)↗
DOE OSTI2022

The Variability of the Black Hole Image in M87 at the Dynamical Timescale

The black hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5–61 days) is comparable to the 6 day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expected structural changes of the images but are free of station-based atmospheric and instrumental errors. We explored the day-to-day variability in closure-phase measurements on all six linearly independent nontrivial baseline triangles that can be formed from the 2017 observations. We showed that three triangles exhibit very low day-to-day variability, with a dispersion of ∼3°–5°. The only triangles that exhibit substantially higher variability (∼90°–180°) are the ones with baselines that cross the visibility amplitude minima on the u–v plane, as expected from theoretical modeling. We used two sets of general relativistic magnetohydrodynamic simulations to explore the dependence of the predicted variability on various black hole and accretion-flow parameters. We found that changing the magnetic field configuration, electron temperature model, or black hole spin has a marginal effect on the model consistency with the observed level of variability. On the other hand, the most discriminating image characteristic of models is the fractional width of the bright ring of emission. Models that best reproduce the observed small level of variability are characterized by thin ring-like images with structures dominated by gravitational lensing effects and thus least affected by turbulence in the accreting plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Spread of Ocean Heat Uptake Efficiency Traced to Ocean Salinity

Abstract Ocean heat uptake (OHU) plays an important role in determining the rate of surface warming under CO 2 forcing. We examined the transient response to CO 2 forcing in a set of 30 climate models from Coupled Model Intercomparison Project Phase 6 (CMIP6) and found that the intermodel spread in OHU efficiency largely depends on the base‐state ocean stratification, especially in the Southern Ocean. This dependence is primarily traced to ocean salinity rather than ocean temperature. Models with weaker ocean stratification primarily due to higher upper‐ocean salinity tend to sequester heat into the deeper ocean, leading to a lower rate of surface warming, a primary reason for higher OHU efficiency; the stratification impact on total OHU is secondary. We applied the relationship between salinity and OHU efficiency for an emergent constraint on OHU efficiency, suggesting an OHU efficiency higher than the CMIP6 multimodel average and arguing against models with extremely low efficiency.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Ab initio evidence for nonthermal characteristics in ultrafast laser melting

Laser melting of semiconductors has been observed for almost 40 years; surprisingly, it is not well understood where most theoretical simulations show a laser-induced thermal process. Ab initio nonadiabatic simulations based on real-time time-dependent density functional theory reveal intrinsic nonthermal melting of silicon, at a temperature far below the thermal melting temperature of 1680 K. Both excitation threshold and time evolution of diffraction intensity agree well with experiment. Nonthermal melting is attributed to excitation-induced drastic changes in bonding electron density, and the subsequent decrease in the melting barrier, rather than lattice heating as previously assumed in the two-temperature models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2004

Analysis of 1970-1995 Trends in Tropospheric Ozone at Northern Hemisphere Midlatitudes with the GEOS-CHEM Model

I ] The causes of trends in tropospheric ozone at Northern Hemisphere midlatitudes from 1970 to 1995 are investigated with the GEOS-CHEM model, a global three-dimensional model of the troposphere driven by assimilated meteorological observations from the Goddard Earth Observing System (GEOS). This model is used to investigate the sensitivity of tropospheric ozone with respect to (1) changes in the anthropogenic emission of nitrogen oxides and nonmethane hydrocarbons, (2) increases in methane concentrations, (3) variations in the stratospheric source of ozone, (4) changes in solar radiation resulting from stratospheric ozone depletion, and ( 5 ) increases in tropospheric temperatures. Model results indicate that local increases in NO, emissions have caused most of the increases seen in lower tropospheric ozone over Europe and Japan. Increases in methane are responsible for roughly one fifth of the anthropogenically induced increase in tropospheric ozone at northern midlatitudes. However, changes in ozone precursors do not adequately explain either the spatial differences in observed ozone trends across midlatitudes or the observed decreases in ozone over Canada throughout the troposphere. We argue that ozone depletion in the lowermost stratosphere is likely to have reduced the stratospheric source by as much as 30% from the early 1970s to the mid 1990s. Model simulations that account for such a reduction along with reported changes in anthropogenic emissions show steep declines of ozone in the upper troposphere and variable increases in the lower troposphere that are more consistent with observations. Differential temperature trends in summer between North America and Europe may account for at least some of the remaining spatial variation in tropospheric ozone trends. Increases in ultraviolet (UV) radiation due to stratospheric ozone depletion do not appear to significantly reduce tropospheric ozone, except at midlatitudes in the Southern Hemisphere following the breakup of the ozone hole.

Fusco, Andrew C.↗
NASA NTRS2013

Analysis of Present Day and Future OH and Methane Lifetime in the ACCMIP Simulations

Results from simulations performed for the Atmospheric Chemistry and Climate Modeling Intercomparison Project (ACCMIP) are analysed to examine how OH and methane lifetime may change from present day to the future, under different climate and emissions scenarios. Present day (2000) mean tropospheric chemical lifetime derived from the ACCMIP multi-model mean is 9.8+/-1.6 yr (9.3+/-0.9 yr when only including selected models), lower than a recent observationally-based estimate, but with a similar range to previous multi-model estimates. Future model projections are based on the four Representative Concentration Pathways (RCPs), and the results also exhibit a large range. Decreases in global methane lifetime of 4.5 +/- 9.1% are simulated for the scenario with lowest radiative forcing by 2100 (RCP 2.6), while increases of 8.5+/-10.4% are simulated for the scenario with highest radiative forcing (RCP 8.5). In this scenario, the key driver of the evolution of OH and methane lifetime is methane itself, since its concentration more than doubles by 2100 and it consumes much of the OH that exists in the troposphere. Stratospheric ozone recovery, which drives tropospheric OH decreases through photolysis modifications, also plays a partial role. In the other scenarios, where methane changes are less drastic, the interplay between various competing drivers leads to smaller and more diverse OH and methane lifetime responses, which are difficult to attribute. For all scenarios, regional OH changes are even more variable, with the most robust feature being the large decreases over the remote oceans in RCP8.5. Through a regression analysis, we suggest that differences in emissions of non-methane volatile organic compounds and in the simulation of photolysis rates may be the main factors causing the differences in simulated present day OH and methane lifetime. Diversity in predicted changes between present day and future OH was found to be associated more strongly with differences in modelled temperature and stratospheric ozone changes. Finally, through perturbation experiments we calculated an OH feedback factor (F) of 1.24 from present day conditions (1.50 from 2100 RCP8.5 conditions) and a climate feedback on methane lifetime of 0.33+-0.13 yr/K, on average. Models that did not include interactive stratospheric ozone effects on photolysis showed a stronger sensitivity to climate, as they did not account for negative effects of climate-driven stratospheric ozone recovery on tropospheric OH, which would have partly offset the overall OH/methane lifetime response to climate change.

atmospheric composition↗
NASA NTRS1990

Filament winding cylinders. III - Selection of the process variables

By using the Lee-Springer filament winding model temperatures, degrees of cure, viscosities, stresses, strains, fiber tensions, fiber motions, and void diameters were calculated in graphite-epoxy composite cylinders during the winding and subsequent curing. The results demonstrate the type of information which can be generated by the model. It is shown, in reference to these results, how the model, and the corresponding WINDTHICK code, can be used to select the appropriate process variables.

Lee, Soo-Yong↗
DOE OSTI2018

Nuclear level density, thermalization, chaos, and collectivity

The knowledge of the level density is necessary for understanding nuclear reactions involving excited nuclear states. In particular, it is an important element in description of astrophysical processes and in technological applications. This review article explains main ideas of physics forming the level density in complex nuclei that grows very fast due to combinatorial complexity of total excitation energy shared by many constituents. This can be translated into a language of statistical physics by the Darwin–Fowler method. We briefly go through the historical development from the nuclear Fermi-gas model to the self-consistent mean field including the pairing effects. At the next step we introduce the ideas of thermalization in a closed mesoscopic system and quantum chaos with very complicated eigenfunctions. This is supported by the experience of the shell model in a limited orbital space that either provides an exact solution or uses the Monte Carlo approach. The statistical method of moments allows one to avoid the exact diagonalization keeping intact the quality of the results. We discuss the popular “constant temperature model” that describes well available data and the shell-model results; it is shown that its success cannot be explained by the phase transition from superfluid to a normal phase. The interpretation is suggested, supported by the numerical studies, in terms of dynamical chaotization including the collective enhancement of the level density. The role of incoherent collision-like interactions is stressed as a necessary element of the thermalization process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Fit-free optical determination of electronic thermalization time in nematic iron-based superconductors

We present a nematic response function model (NRFM) for fit-free direct extraction of the characteristic time of ultrafast electronic thermalization in iron-based superconductors, materials with electronic nematicity. By combining the NRFM for polarization-dependent pump–probe measurements of electronic nematic response with the two-temperature model (TTM) for sub-picosecond quasiparticle relaxation, we quantify the electronic thermalization timescales and their anisotropy. The nematic response function is modeled as the difference in normalized reflectivity signals, revealing a pronounced sub-picosecond extremum in signal evolution that directly yields the characteristic electronic thermalization time. This method demonstrates that the NRFM is consistent with TTM fits of transient optical response, yielding electronic thermalization time constants on the order of 110–230 fs for the FeSe 1−x Te x and Ba(Fe 0.92 Co 0.08 ) 2 As 2 thin films. Here, the proposed approach can be applied to any material that exhibits electronic nematicity, providing a powerful tool for direct mapping of the relaxation time in nematic materials, avoiding complex experimental data-fitting procedures.

Bartenev, Alexander [University of Puerto Rico, Ma↗