Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equilibrium modeling”

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 1,387 records · Page 77

Dissecting Anvil Cloud Response to Sea Surface Warming

Abstract We derive an anvil cloud diagnostic from the continuity equation of cloud ice and apply it to the output of convection‐permitting Energy Exascale Earth System Model (E3SM) simulations run in radiative‐convective equilibrium mode. This diagnostic shows that anvil cloud fraction can be reliably diagnosed as a product of cloud detrainment and lifetime. Detrainment is found to be approximated well by a product of clear sky convergence and cloud ice mixing ratio, while cloud lifetime is dominated by sedimentation. Taken together, this diagnostic expresses anvil cloud fraction as a function of five physically measurable quantities. Of these, clear‐sky convergence changes drive the anvil cloud reduction with warming while an increase in cloud ice mixing ratio buffers the decrease. Accordingly, this study provides a theoretical foundation upon which the Stability‐Iris hypothesis can be tested.

54 ENVIRONMENTAL SCIENCES↗

Phase equilibrium of liquid water and hexagonal ice from enhanced sampling molecular dynamics simulations

We study the phase equilibrium between liquid water and ice Ih modeled by the TIP4P/Ice interatomic potential using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice Ih-liquid free energy differences from simulations that visit reversibly both phases. The reversible interconversion is achieved by introducing a static bias potential as a function of an order parameter. The order parameter was tailored to crystallize the hexagonal diamond structure of oxygen in ice Ih. Furthermore, we analyze the effect of the system size on the ice Ih-liquid free energy differences, and we obtain a melting temperature of 270 K in the thermodynamic limit. This result is in agreement with estimates from thermodynamic integration (272 K) and coexistence simulations (270 K). Since the order parameter does not include information about the coordinates of the protons, the spontaneously formed solid configurations contain proton disorder as expected for ice Ih.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Invariant regimes of Spencer scaling law for magnetic compression of rotating FRC plasma

Abstract The scaling laws for the magnetic compression of a toroidally rotating field reversed configuration (FRC) have been investigated in this work. The magnetohydrodynamics (MHD) simulations of the magnetic compression on rotating FRCs employing the NIMROD code (Sovinec et al 2004 J. Comput. Phys. 195 355), are compared with the Spencer’s one-dimensional (1D) theory (Spencer et al 1983 Phys. Fluids 26 1564) for a wide range of initial flow speeds and profiles. The toroidal flow can influence the scalings directly through the alteration of the compressional work as also evidenced in the 1D adiabatic model, and indirectly by reshaping the initial equilibrium. However, in comparison to the static initial FRC equilibrium cases, the pressure and the radius scalings remain invariant for the magnetic compression ratio B w 2 / B w 1 up to 6 in presence of the initial equilibrium flow, suggesting a broader applicable regime of the Spencer scaling law for FRC magnetic compression. The invariant scaling has been proven a natural consequence of the conservation of angular momentum of both fluid and magnetic field during the dynamic compression process.

Ma, Yiming↗

Wealth inequality in the prehispanic northern US Southwest: from Malthus to Tyche

Persistent differences in wealth and power among prehispanic Pueblo societies are visible from the late AD 800s through the late 1200s, after which large portions of the northern US Southwest were depopulated. In this paper we measure these differences in wealth using Gini coefficients based on house size, and show that high Ginis (large wealth differences) are positively related to persistence in settlements and inversely related to an annual measure of the size of the unoccupied dry-farming niche. We argue that wealth inequality in this record is due first to processes inherent in village life which have internally different distributions of the most productive maize fields, exacerbated by the dynamics of systems of balanced reciprocity; and second to decreasing ability to escape village life owing to shrinking availability of unoccupied places within the maize dry-farming niche as villages get enmeshed in regional systems of tribute or taxation. We embed this analytical reconstruction in the model of an ‘Abrupt imposition of Malthusian equilibrium in a natural-fertility, agrarian society’ proposed by Puleston et al. (Puleston C, Tuljapurkar S, Winterhalder B. 2014 PLoS ONE 9, e87541 (doi:10.1371/journal.pone.0087541)), but show that the transition to Malthusian dynamics in this area is not abrupt but extends over centuries. This article is part of the theme issue ‘Evolutionary ecology of inequality’.

99 GENERAL AND MISCELLANEOUS↗

Stretched exponential magnetic relaxation dynamics in artificial square ice revealed through x-ray photon correlation spectroscopy

X-ray photon correlation scattering measurements are undertaken on a thermally active artificial spin ice based on the square lattice, referred to as artificial square ice, to probe the fluctuation timescales as a function of temperature as the system passes through the paramagnetic-antiferromagnetic phase transition, which belongs to the two-dimensional Ising universality class. In the paramagnetic regime, a single exponential timescale is seen, whereas at and below the critical temperature, a stretched exponential decorrelation is observed, with the stretching exponent decreasing from unity down to below one-half as the temperature reduces. This trend is confirmed by kinetic Monte Carlo simulations of a simplified point-dipolar square ice system, and is in agreement with past theoretical work on the kinetic Ising model where stretched exponential relaxation due to equilibrium domain wall dynamics below the critical temperature was found.

36 MATERIALS SCIENCE↗

Dynamics of the O ( 4 ) critical point in QCD: Critical pions and diffusion in model G

We present a detailed study of the finite momentum dynamics of the O ( 4 ) critical point of QCD, which lies in the dynamic universality class of “model G.” The critical scaling of the model is analyzed in multiple dynamical channels. For instance, the finite momentum analysis allows us to precisely extract the pion dispersion curve below the critical point. The pion velocity is in striking agreement with the predictions relation and static universality. The pion damping rate and velocity are both consistent with the dynamical critical exponent ζ = 3 / 2 of model G. Similarly, although the critical amplitude for the diffusion coefficient of the conserved O ( 4 ) charges is small, it is clearly visible both in the restored phase and with finite explicit symmetry breaking, and its dynamical scaling is again consistent with ζ = 3 / 2 . We determine a new set of universal dynamical critical amplitude ratios relating the diffusion coefficient to a suitably defined order parameter relaxation time. We also show that in a finite volume simulation, the chiral condensate diffuses on the coset manifold in a manner consistent with dynamical scaling, and with a diffusion coefficient that is determined by the transport coefficients of hydrodynamic pions. Finally, the amplitude ratios (together with other nonuniversal amplitudes also reported here) compile all relevant information for further studies of model G both in and out of equilibrium. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Command of active and responsive elastomers by topological defects and patterns

The project resulted in the development of stimuli-responsive elastomer coatings formed by photopolymerized liquid crystal molecules. The molecular orientation of the liquid crystal elastomers is coupled to rubber-like elasticity. The project developed an approach to produce complex patterns of molecular orientation in elastomers by employing photopatterning technique based on plasmonic metamasks that convert unpolarized incoming light into a transmitted light beam with spatially-varying linear polarization. The polarization-modulated light beam aligns the film of azodye molecules, which serves as the substrate for liquid crystal elastomer coating. The alignment of the azodye molecules follows the polarization pattern with a phase shift of the alignment direction by 90 degrees. The azodye substrate imposes orientation of molecules in the adjacent liquid crystal layer which is preserved after the monomers are polymerized into a liquid crystal elastomer. By using differently aligned substrates one could produce various patterns of molecular orientation in the elastomer, which dictate the mechanical properties of the elastomers and their response to external stimuli such as temperature, humidity, and ultraviolet irradiation. In particular, the predesigned molecular orientation produces deterministically defined topography response of the liquid crystal coatings, in which their free surface changes from flat to either locally elevated or locally compressed, depending on whether the pre-inscribed molecular orientation is predesigned with a splay or bend deformation. The research established that the mechanism of this deterministic relation between the dynamic topography and pre-inscribed molecular orientation is rooted in the change of the tensorial order parameter that characterizes the degree of orientational order, and in the contraction/expansion of polymer molecules in response to the order parameter changes. In particular, lowering the order causes contraction of polymer globules. The work demonstrated that the dynamic topography of elastomer coatings with three-dimensional pattern of orientational order includes up and down motion of the material (along the normal to the coating) and also shifts in lateral directions. We demonstrate that the deterministic relationship between the orientational patterns and variations of coatings’ profile is caused by forces that are mathematically equivalent to active forces in the system of “pullers” or “pushers” in active matter, demonstrating universality of the out-of-equilibrium behavior. The project resulted in theoretical models that predictively describe the dynamic response of coatings. The topography of liquid crystal elastomers is sensitive to stimuli such as temperature, ultraviolet irradiation, and humidity. The applicability of the patterned approach to produce dynamic coatings was expanded from nematic elastomers to their smectic counterparts. The project combined a battery of imaging techniques to unveil the dynamic properties of the coatings, ranging from the state-of-the-art dynamic holographic microscopy to fluorescent confocal polarizing microscopy and polscope microscopy. The research uncovered potential applications of dynamic elastomer coatings, such as control placement of colloidal particles, engineering of living tissues with alignment patterns of cells that follow the orientational order of elastomer coatings. The approach developed in the project can be used to design programmable dynamic coatings with functionalities that mimic biological tissues such as skin.

36 MATERIALS SCIENCE↗

Ablation and radiation coupled viscous hypersonic shock layers, volume 1

The results for a stagnation-line analysis of the radiative heating of a phenolic-nylon ablator are presented. The analysis includes flow field coupling with the ablator surface, equilibrium chemistry, a step-function diffusion model and a coupled line and continuum radiation calculation. This report serves as the documentation, i e. users manual and operating instructions for the computer programs listed in the report.

Engel, C. D.↗

Global exchange of water vapor on Mars

Observations of the global distribution and seasonal variation of the Martian atmospheric water vapor have been made from the Viking orbiters for a continuous period covering more than half of a Martian year. The seasonal dependence of the latitude distribution of the column abundance of vapor is consistent with a model in which the vapor is in equilibrium with the regolith at polar and middle latitudes. The results indicate the existence of a permanent reservoir of water ice buried at a depth of 10 cm to 1 meter at all latitudes poleward of 40 deg. The behavior of the vapor in the summer hemisphere suggests an annual net transport of the vapor phase from the southern to the northern hemisphere, with deposition of ice of thickness on the order of a few milligrams per sq cm at northern polar latitudes. The hemispheric asymmetry is the result of the propagation of global dust storms originating in the south.

Farmer, C. B.↗

Self-acting shaft seals

Report reviews operating principles and design of self-acting seals. Influences of adverse operating conditions are considered also. Elements of analysis used in seal performance predictions are described and evaluated. Mathematical models for obtaining seal force balance and equilibrium film thickness are outlined. Self-acting seals are nonrubbing, have lower leakage rates than labyrinth seals, and are well suited for advanced aircraft engines.

Ludwig, L. P.↗

A probabilistic approach to radiative energy loss calculations for optically thick atmospheres - Hydrogen lines and continua

An approximate probabilistic radiative transfer equation and the statistical equilibrium equations are simultaneously solved for a model hydrogen atom consisting of three bound levels and ionization continuum. The transfer equation for L-alpha, L-beta, H-alpha, and the Lyman continuum is explicitly solved assuming complete redistribution. The accuracy of this approach is tested by comparing source functions and radiative loss rates to values obtained with a method that solves the exact transfer equation. Two recent model solar-flare chromospheres are used for this test. It is shown that for the test atmospheres the probabilistic method gives values of the radiative loss rate that are characteristically good to a factor of 2. The advantage of this probabilistic approach is that it retains a description of the dominant physical processes of radiative transfer in the complete redistribution case, yet it achieves a major reduction in computational requirements.

Canfield, R. C.↗

The aluminum I autoionization doublet in the quiet solar spectrum

Observations are presented of the Al I autoionization doublet 1932 A and 1936 A in the quiet solar spectrum, obtained from the NRL slit spectrograph aboard Skylab and from the University of Hawaii Echelle Rocket Spectrograph. The observed profiles are compared with theoretical spectra computed for the Harvard Smithsonian Reference Atmosphere and the Vernazza, Avrett and Loeser (1976) solar models. It is found that nonlocal thermodynamic equilibrium effects are important in the line-formation problem and the synthetic spectra are in good agreeement with the data.

Heasley, J. N.↗

The bulk-Moon MgO/FeO ratio: A highlands perspective

Compositional data for nonmare (highlands) samples suggest that the Moon's mg ratio (MgO/FeO) is higher than general estimates. Geochemically representative highlands soils have mg ratios of 0.66 (Apollo 16), 0.69 (Luna 20) and 0.73 (ALHA81005). These soils are mixtures of unrelated pristine nonmare rocks, of which there are at least three groups: Mg-rich rocks, ferroan anorthosites, and KREEP. Other than Mg-rich rocks, virtually all pristine rocks have mg 0.65. Thus, assuming the mixing process that sampled Mg-rich materials was random, the average mg of Mg-rich parent magmas was probably at least 0.70. More direct evidence can be derived from the Mg-rich rocks themselves. Nine of them have bulk-rock mg 0.87. Two (15445 A and 67435 PST) contain Fo(92) olivine. Production of melts that crystallized Fo(92) olivine implies that the mg ratios of source regions in lunar mantle were commensurably high. A quantification of this constraint is developed assuming that the parent melts formed by equilibrium (batch) partial melting. Implications of the model are discussed.

Warren, P. H.↗

On the Role Played by Lines in Radiatively Driven Stellar Winds Depending on the Position of the Stars in the HR Diagram

The radiative force due to transfer in ultraviolet lines is always an important mechanism in hot star wind dynamics. However, it is not clear when it is the dominant mechanism and which are the noise parameters. To investigate the efficiency of purely radiative momentum/energy transfer in hot star winds and in various regions of the HR diagram, the Leroy and Lafon model was improved and put to its limits; correlations between the mass loss rate, the luminosity and other parameters and the theoretical and the observational results, looking for observed stars violating the model were compared. It is concluded that in widespread region of the HR diagram, line driven models are consistent with observations, the radiative equilibrium physics is relevant throughout the expanding atmospheres and the mass loss rate is quasilinearly correlated with the luminosity.

Migozzi, M. C.↗

A theoretical investigation of laser-sustained plasma thruster

A numerical code has been successfully developed for the investigation of thruster performance using a laser-sustained hydrogen plasma as the propellant. The plasma was sustained using a 10.6-micron CO2 laser beam focused at different positions within the thruster. The physical model assumed that plasma is in thermodynamical equilibrium (LTE), and geometric ray tracing was adopted to describe the laser beam. The steady-state, axisymmetric, Navier-Stokes equations coupled with the laser power absorption process have been solved numerically. A pressure based Navier-Stokes numerical solver using body-fitted coordinates was used to calculate the laser-supported rocket flow which includes both subsonic and supersonic flow regions. From the limited parametric study, which did not try to optimize the rocket performance, it was found that better performance was obtained when the laser beam was focused closer to the rocket throat.

Jeng, San-Mou↗

Solar transition region and coronal response to heating rate perturbations

Observations of Doppler shifts in UV emission lines formed in the solar transition region show continual plasma downflows and impulsive plasma upflows. Using numerical simulations, the authors examine the conjecture that areas of downflowing plasma are the base regions of coronal loops in which the heating is gradually decreasing and that areas of upflowing plasma are the base regions of coronal loops in which the heating rate is gradually increasing. Beginning with a coronal loop in equilibrium, the heating rate is reduced on time scales of 100, 1000, and 2000 s to 10 percent and 1 percent of the initial value, and the loop is allowed to evolve to a new equilibrium. The heating rate for the cooled models is then increased back to the initial value on the same time scales. While significant mass motions do develop in the simulations, both the emission measure and the velocity at 100,000 K do not show the characteristics present in UV observations.

Mariska, John T.↗

Plasma diffusion at the magnetopause - The case of lower hybrid drift waves

The diffusion expected from the quasi-linear theory of the lower hybrid drift instability at the earth's magnetopause is recalculated. The resulting diffusion coefficient is marginally large enough to explain the thickness of the boundary layer under quiet conditions, based on observational upper limits for the wave intensities. Thus, one possible model for the boundary layer could involve equilibrium between the diffusion arising from lower hybrid waves and various loss processes.

Treumann, R. A.↗

Estimation of the initial equilibrium constants in the formation of tetragonal lysozyme nuclei

Results are presented from a study of the equilibria, kinetic rates, and the aggregation pathway which leads from a lysozyme monomer crystal to a tetragonal crystal, using dialyzed and recrystallized commercial hen eggwhite lysozyme. Relative light scattering intensity measurements were used to estimate the initial equilibrium constants for undersaturated lysozyme solutions in the tetragonal regime. The K1 value was estimated to be (1-3) x 10 exp 4 L/mol. Estimates of subsequent equilibrium constants depend on the crystal aggregation model chosen or determined. Experimental data suggest that tetragonal lysozyme crystal grows by addition of aggregates preformed in the bulk solution, rather than by monomer addition.

Pusey, Marc L.↗