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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Force-free magnetic fields - Generating functions and footpoint displacements

This paper presents analytic and numerical calculations that explore equilibrium sequences of bipolar force-free magnetic fields in relation to displacments of their magnetic footpoints. It is shown that the appearance of magnetic islands - sometimes interpreted as marking the loss of equilibrium in models of the solar atmosphere - is likely associated only with physically unrealistic footpoint displacements such as infinite separation or 'tearing' of the model photosphere. The work suggests that the loss of equilibrium in bipolar configurations, sometimes proposed as a mechanism for eruptive solar events, probably requires either fully three-dimensional field configurations or nonzero plasma pressure. The results apply only to fields that are strictly bipolar, and do not rule out equilibrium loss in more complex structures such as quadrupolar fields.

Wolfson, Richard↗
NASA NTRS1986

Magnetospheric equilibrium configurations and slow adiabatic convection

This review paper demonstrates how the magnetohydrostatic equilibrium (MHE) theory can be used to describe the large-scale magnetic field configuration of the magnetosphere and its time evolution under the influence of magnetospheric convection. The equilibrium problem is reviewed, and levels of B-field modelling are examined for vacuum models, quasi-static equilibrium models, and MHD models. Results from two-dimensional MHE theory as they apply to the Grad-Shafranov equation, linear equilibria, the asymptotic theory, magnetospheric convection and the substorm mechanism, and plasma anisotropies are addressed. Results from three-dimensional MHE theory are considered as they apply to an intermediate analytical magnetospheric model, magnetotail configurations, and magnetopause boundary conditions and the influence of the IMF.

Voigt, Gerd-Hannes↗
NASA NTRS2022

Atomic-Scale Structure and Non-Stoichiometry of Meteoritic Hibonite: A Transmission Electron Microscope Study

Hibonite (CaAl_12O_19) is a common refractory mineral in Ca-Al-rich inclusions (CAIs) in primitive meteorites. Transmission electron microscope (TEM) studies have identified enigmatic planar defects in different occurrences of hibonite in the Allende meteorite that give rise to strong streaking along c* in electron diffraction patterns. Atomic resolution high-angle annular dark-field (HAADF) imaging and energy-dispersive X-ray (EDX) analyses were used to determine the nature and origin of these planar features. HAADF images of hibonite grains reveal lamellar intergrowths of common 1.6 nm spacing, and less commonly 2.0 and 2.5 nm spacings, interspersed in stoichiometric hibonite showing 1.1 nm (002) spacing. Stoichiometric hibonite consists of alternating Ca-containing (“R”) and spinel-structured (“S”) blocks stacked in a sequence RS. In contrast, the 1.6 nm layers result from a doubled S block such that the stacking sequence is RSS, while in the widest defect observed, the stacking sequence is RSSSS. These intergrowths are epitaxial and have coherent, low-strain boundaries with the host hibonite Meteoritic hibonite shows common Ti and Mg substitution for Al in its structure. Atomic-resolution EDX maps of hibonite grains in the Allende CAI confirm the preferred site occupancy of Mg on tetragonal M3 sites in S blocks and of Ti on trigonal bipyramidal M2 and octahedral M4 sites in R blocks. Mg is highly concentrated, but Ti is absent in the planar defects where wider S blocks show Al-rich compositions compared to stoichiometric MgAl_2O_4 spinel. Therefore, Mg likely played the major role in the formation and metastability of planar defects in hibonite. Electron energy loss spectroscopy data from the Ti L_2,3 edge show the presence of mixed Ti oxidation states with ~ 15–20% of Ti as Ti^(3+) in hibonite, suggesting a direct substitution of Ti^(3+) ↔ Al^(3+) in hibonite. The remaining ~ 80–85% of Ti is present as Ti^(4+) and corresponding EDX analyses are consistent with the well-known coupled substitution 2Al^(3+) ↔ Ti^(4+) + Mg^(2+) being the major mechanism for Ti and Mg substitution in hibonite. The formation of planar defects in hibonite occurred during high-temperature nebular condensation or melting/crystallization processes. The occurrence of non-stoichiometric hibonite in the Allende CAI deviates from the mineral formation sequence predicted from equilibrium condensation models. Overall, our atomic resolution TEM observations signify non-equilibrium, kinetic-controlled crystal growth during the high-temperature formation of refractory solids in the early solar nebula.

hibonite↗
NASA NTRS1995

Far-infrared lines from H II regions: Abundance variations in the galaxy

Far-infrared lines of (N III) (57 microns), (O III) (52, 88 microns), (Ne III) (36 microns), and (S III) (19, 33 microns) have been measured in the H II regions G1.13 - 0.11, W31B, G23.95 + 0.15, G25.38 - 0.18, G29.96 - 0.02, W43, W51e, S156, S158, NGC 3576, NGC 3603, and G298.22-0.34. These observations were made with the facility Cryogenic Grating Spectrometer on the Kuiper Airborne Observatory to examine variations in abundances throughout the Galaxy. Previously published observations of G0.095 + 0.012, G333.60 - 0.21, G45.13 + 0.14A, K3-50, and M17 are also discussed. The giant H II region 30 Doradus in the Large Magellanic Cloud (LMC) was observed for comparison. Fluxes for (Ne II) (12.8 microns), (S IV) (10.5 microns), and the radio free-free continuum were collected from the literature for those sources. Electron densities were estimated from FIR line-pair ratios, and ionic abundances were estimated from the FIR line and radio fluxes. The excitation was estimated from the O(2+)/S(2+) ratio. Corrections for unseen ionization stages were calculated with the use of constnat-density H II region models. The validity and range of applicability of such semiempirical ionization correction schemes are discussed. The abundances with respect to hydrogen exhibit gradients with R(sub G) comparable to those previously measured for our Galaxy and for other galaxies. The overall gradients are d (log N/H)/dR = -0.10 +/- 0.02 dex/kpc, d(log Ne/H)/dR = -0.08 +/- 0.02 dex/kpc and d(log S/H)/dR = 0.07 +/- 0.02 dex/kpc. Compared to the Orion Nebula, the intermediate R(sub G) H II regions with 6 is less than R(sub G) is less than 11 kpc have similar or lower S/H and N/O ratios. The N/O ratios in the inner Galaxy are more than twice those observed in the Orion Nebula and intermediate R(sub G) H II regions. In fact, all the abundance ratios are as well or better fitted by a step fit with two levels than by a linear gradient. As has been noted in previous studies, the N/O ratio estimated from infrared observations of the doubly ionized N and O lines in H II regions is larger than the ratio estimated from optical observations of the singly ionized N and O lines. The Ne(2+)/O(2+) ratio is observed to be essentially constant over a wide range of excitation. This contradicts predictions of model H II regions calculated with the use of Local Thermodynamic Equilibrium (LTE) model stellar atmospheres. We conclude that these stellar atmospheres significantly underestimate the actual emergent fluxes for energies greater than 41 eV.

Simpson, Janet P.↗
NASA NTRS1989

Recent advances in Runge-Kutta schemes for solving 3-D Navier-Stokes equations

A thin-layer Navier-Stokes has been developed for solving high Reynolds number, turbulent flows past aircraft components under transonic flow conditions. The computer code has been validated through data comparisons for flow past isolated wings, wing-body configurations, prolate spheroids and wings mounted inside wind-tunnels. The basic code employs an explicit Runge-Kutta time-stepping scheme to obtain steady state solution to the unsteady governing equations. Significant gain in the efficiency of the code has been obtained by implementing a multigrid acceleration technique to achieve steady-state solutions. The improved efficiency of the code has made it feasible to conduct grid-refinement and turbulence model studies in a reasonable amount of computer time. The non-equilibrium turbulence model of Johnson and King has been extended to three-dimensional flows and excellent agreement with pressure data has been obtained for transonic separated flow over a transport type of wing.

Vatsa, Veer N.↗
NASA NTRS1993

Towards a generalized computational fluid dynamics technique for all Mach numbers

Currently there exists no single unified approach for efficiently and accurately solving computational fluid dynamics (CFD) problems across the Mach number regime, from truly low speed incompressible flows to hypersonic speeds. There are several CFD codes that have evolved into sophisticated prediction tools with a wide variety of features including multiblock capabilities, generalized chemistry and thermodynamics models among other features. However, as these codes evolve, the demand placed on the end user also increases simply because of the myriad of features that are incorporated into these codes. In order for a user to be able to solve a wide range of problems, several codes may be needed requiring the user to be familiar with the intricacies of each code and their rather complicated input files. Moreover, the cost of training users and maintaining several codes becomes prohibitive. The objective of the current work is to extend the compressible, characteristic-based, thermochemical nonequilibrium Navier-Stokes code GASP to very low speed flows and simultaneously improve convergence at all speeds. Before this work began, the practical speed range of GASP was Mach numbers on the order of 0.1 and higher. In addition, a number of new techniques have been developed for more accurate physical and numerical modeling. The primary focus has been on the development of optimal preconditioning techniques for the Euler and the Navier-Stokes equations with general finite-rate chemistry models and both equilibrium and nonequilibrium thermodynamics models. We began with the work of Van Leer, Lee, and Roe for inviscid, one-dimensional perfect gases and extended their approach to include three-dimensional reacting flows. The basic steps required to accomplish this task were a transformation to stream-aligned coordinates, the formulation of the preconditioning matrix, incorporation into both explicit and implicit temporal integration schemes, and modification of the numerical flux formulae. In addition, we improved the convergence rate of the implicit time integration schemes in GASP through the use of inner iteration strategies and the use of the GMRES (General Minimized Resisual) which belongs to the class of algorithms referred to as Krylov subspace iteration. Finally, we significantly improved the practical utility of GASP through the addition of mesh sequencing, a technique in which computations begin on a coarse grid and get interpolated onto successively finer grids. The fluid dynamic problems of interest to the propulsion community involve complex flow physics spanning different velocity regimes and possibly involving chemical reactions. This class of problems results in widely disparate time scales causing numerical stiffness. Even in the absence of chemical reactions, eigenvalue stiffness manifests itself at transonic and very low speed flows which can be quantified by the large condition number of the system and evidenced by slow convergence rates. This results in the need for thorough numerical analysis and subsequent implementation of sophisticated numerical techniques for these difficult yet practical problems. As a result of this work, we have been able to extend the range of applicability of compressible codes to very low speed inviscid flows (M = .001) and reacting flows.

Walters, R. W.↗
NASA NTRS2003

Ising model of cardiac thin filament activation with nearest-neighbor cooperative interactions

We have developed a model of cardiac thin filament activation using an Ising model approach from equilibrium statistical physics. This model explicitly represents nearest-neighbor interactions between 26 troponin/tropomyosin units along a one-dimensional array that represents the cardiac thin filament. With transition rates chosen to match experimental data, the results show that the resulting force-pCa (F-pCa) relations are similar to Hill functions with asymmetries, as seen in experimental data. Specifically, Hill plots showing (log(F/(1-F)) vs. log [Ca]) reveal a steeper slope below the half activation point (Ca(50)) compared with above. Parameter variation studies show interplay of parameters that affect the apparent cooperativity and asymmetry in the F-pCa relations. The model also predicts that Ca binding is uncooperative for low [Ca], becomes steeper near Ca(50), and becomes uncooperative again at higher [Ca]. The steepness near Ca(50) mirrors the steep F-pCa as a result of thermodynamic considerations. The model also predicts that the correlation between troponin/tropomyosin units along the one-dimensional array quickly decays at high and low [Ca], but near Ca(50), high correlation occurs across the whole array. This work provides a simple model that can account for the steepness and shape of F-pCa relations that other models fail to reproduce.

Evaluation Studies↗
NASA NTRS1985

Predictions of transonic separated flow with an eddy-viscosity/Reynolds-shear-stress closure model

In previous work, a new turbulence-closure model was specifically developed for two-dimensional turbulent boundary layers subjected to strong, adverse-pressure gradients and the attendant separation. This closure model was shown to perform well, but the inverse, boundary-layer calculations used in evaluating the model were not truly predictive since the mass-flux parameter distribution obtained from the experimental data had to be specified. The objective of this study was to evaluate this closure model for a series of strong, transonic, inviscid-viscous interactions with varying degrees of separation using a fully predictive calculation method based on the Reynolds-averaged, Navier-Stokes equations. Calculated results are compared with experimental results and are shown to be in excellent agreement, even for the interactions with massive separation. This closure model has the very favorable property of requiring little more computational effort than equilibrium, algebraic models.

Johnson, D. A.↗
NASA NTRS1989

The effect of electron temperature and impact ionization on Martian return AOTV flowfields

Various electron impact ionization models in conjunction with a quasi-equilibrium electron temperature model have been investigated and applied to the stagnation region of a hypothetical 2.3 m nose radius Martian return AOTV. For the conditions considered, U = 12 km/sec at 80 km, both multi-temperature inviscid and viscous results indicate that a two-step ionization impact model predicts ionization distances in agreement with experimental data, that nonequilibrium chemistry and radiation effects are important throughout the stagnation zone, and that the quasi-equilibrium electron temperature model is reasonable. Also, using a non-grey emission-absorption radiation step model, it is shown that nonequilibrium causes a reduction in radiative heating from that predicted for equilibrium conditions and that compared to an adiabatic wall a cool wall (1650 deg K) results in a 28 to 45 percent reduction in radiative heating due to absorption near the wall.

Carlson, Leland A.↗
NASA NTRS1985

Separated transonic airfoil flow calculations with a nonequilibrium turbulence model

Navier-Stokes transonic airfoil calculations based on a recently developed nonequilibrium, turbulence closure model are presented for a supercritical airfoil section at transonic cruise conditions and for a conventional airfoil section at shock-induced stall conditions. Comparisons with experimental data are presented which show that this nonequilibrium closure model performs significantly better than the popular Baldwin-Lomax and Cebeci-Smith equilibrium algebraic models when there is boundary-layer separation that results from the inviscid-viscous interactions.

King, L. S.↗
NASA NTRS1990

Numerical solution of 3D Navier-Stokes equations with upwind implicit schemes

An upwind MUSCL type implicit scheme for the three-dimensional Navier-Stokes equations is presented. Comparison between different approximate Riemann solvers (Roe and Osher) are performed and the influence of the reconstructions schemes on the accuracy of the solution as well as on the convergence of the method is studied. A new limiter is introduced in order to remove the problems usually associated with non-linear upwind schemes. The implementation of a diagonal upwind implicit operator for the three-dimensional Navier-Stokes equations is also discussed. Finally the turbulence modeling is assessed. Good prediction of separated flows are demonstrated if a non-equilibrium turbulence model is used.

Marx, Yves P.↗
NASA NTRS1992

Model free simulations of a high speed reacting mixing layer

The effects of compressibility, chemical reaction exothermicity and non-equilibrium chemical modeling in a combusting plane mixing layer were investigated by means of two-dimensional model free numerical simulations. It was shown that increased compressibility generally had a stabilizing effect, resulting in reduced mixing and chemical reaction conversion rate. The appearance of 'eddy shocklets' in the flow was observed at high convective Mach numbers. Reaction exothermicity was found to enhance mixing at the initial stages of the layer's growth, but had a stabilizing effect at later times. Calculations were performed for a constant rate chemical rate kinetics model and an Arrhenius type kinetics prototype. The Arrhenius model was found to cause a greater temperature increase due to reaction than the constant kinetics model. This had the same stabilizing effect as increasing the exothermicity of the reaction. Localized flame quenching was also observed when the Zeldovich number was relatively large.

Steinberger, Craig J.↗
NASA NTRS1991

Continuous system modeling

A comprehensive and systematic introduction is presented for the concepts associated with 'modeling', involving the transition from a physical system down to an abstract description of that system in the form of a set of differential and/or difference equations, and basing its treatment of modeling on the mathematics of dynamical systems. Attention is given to the principles of passive electrical circuit modeling, planar mechanical systems modeling, hierarchical modular modeling of continuous systems, and bond-graph modeling. Also discussed are modeling in equilibrium thermodynamics, population dynamics, and system dynamics, inductive reasoning, artificial neural networks, and automated model synthesis.

Cellier, Francois E.↗
NASA NTRS1991

Numerical models of giant planet formation with rotation

A 1D quasi-spherical approximation is presently used under the assumption of hydrostatic equilibrium to model giant-planet formation processes that encompass the accretion and transport of angular momentum for both radiative and convective zones. These calculations have been conducted up to masses comparable to those of Saturn. Angular momentum transport is strong, leading to strong concentration of angular momentum in the outer convective zones which develop envelope masses greater than 30 earth masses. When the models are allowed to contract near the calculation's end, the ratio of centrifugal force to gravity at the outer radius rises sufficiently to validate the quasi-spherical approximation.

Korycansky, D. G.↗
NASA NTRS1990

Ballooning instability of the earth's plasma sheet region in the presence of parallel flow

Stability of the plasma sheet and plasma sheet boundary layer against the ballooning mode instability is investigated. The equilibrium state of a two-dimensional plasma sheet configuration with parallel sheared flow is modeled. This equilibrium is shown to be ballooning unstable when delta-W is not positive definite, where delta-W is the potential energy. The eigenmode structure of the ballooning mode is found by imposing the boundary conditions that the waves are totally reflected from the ionosphere, and that no waves are coming in from infinity. The eigenmode structure of the unstable balloning modes is highly oscillatory, extending beyond about 100 R(E). The ballooning modes are thus a possible candidate for explaining the MHD waves and other dynamical events observed in the magnetotail by ISEE 3 and other spacecraft.

Lakhina, G. S.↗
NASA NTRS1993

Solar flare model atmospheres

Solar flare model atmospheres computed under the assumption of energetic equilibrium in the chromosphere are presented. The models use a static, one-dimensional plane parallel geometry and are designed within a physically self-consistent coronal loop. Assumed flare heating mechanisms include collisions from a flux of non-thermal electrons and x-ray heating of the chromosphere by the corona. The heating by energetic electrons accounts explicitly for variations of the ionized fraction with depth in the atmosphere. X-ray heating of the chromosphere by the corona incorporates a flare loop geometry by approximating distant portions of the loop with a series of point sources, while treating the loop leg closest to the chromospheric footpoint in the plane-parallel approximation. Coronal flare heating leads to increased heat conduction, chromospheric evaporation and subsequent changes in coronal pressure; these effects are included self-consistently in the models. Cooling in the chromosphere is computed in detail for the important optically thick HI, CaII and MgII transitions using the non-LTE prescription in the program MULTI. Hydrogen ionization rates from x-ray photo-ionization and collisional ionization by non-thermal electrons are included explicitly in the rate equations. The models are computed in the 'impulsive' and 'equilibrium' limits, and in a set of intermediate 'evolving' states. The impulsive atmospheres have the density distribution frozen in pre-flare configuration, while the equilibrium models assume the entire atmosphere is in hydrostatic and energetic equilibrium. The evolving atmospheres represent intermediate stages where hydrostatic equilibrium has been established in the chromosphere and corona, but the corona is not yet in energetic equilibrium with the flare heating source. Thus, for example, chromospheric evaporation is still in the process of occurring.

Hawley, Suzanne L.↗