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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

Experimental Results and Interfacial Lift-off Model Predictions of Critical Heat Flux for Flow Boiling with Subcooled Inlet Conditions – In Preparation for Experiments Onboard the International Space Station

This study investigates critical heat flux (CHF) for subcooled flow boiling of n-Perfluorohexane based on results of pre-launch Earth-gravity Mission Sequence Tests (MSTs) of the Flow Boiling and Condensation Experiment (FBCE), which was launched to the International Space Station (ISS) in August 2021. CHF measurements were made in a rectangular channel having a 2.5 mm by 5 mm cross-section and a 114.6-mm long heated segment. Both single-sided and double-sided heating were tested in vertical upflow in Earth gravity for a variety of inlet conditions. The inlet subcooling was varied in the range of 0.4 – 32.0°C and encompassed both near-saturated and highly subcooled conditions. Experimental trends and high-speed video records were investigated to better understand the mechanism of CHF. Overall trends show CHF increases as flow rate and/or inlet subcooling are increased. Flow features from the events around CHF justify the applicability of the Interfacial Lift-off Model and the determination of limiting criteria for its application. The present experimental data are combined with prior databases for various flow orientations with respect to Earth gravity and microgravity data collected on parabolic flights. Predictions are made using the Interfacial Lift-off Model for this consolidated subcooled-inlet FBCE-CHF database. A heat utility ratio was included in the model to capture the effects of subcooling and corresponding thermodynamic non-equilibrium. An overall mean absolute error of 19.04% indicates good predictive capability of the model for both heating configurations, different gravity environments, and a wide range of inlet subcooling.

flow boiling↗
NASA NTRS1998

WARP3D-Release 10.8: Dynamic Nonlinear Analysis of Solids using a Preconditioned Conjugate Gradient Software Architecture

This report describes theoretical background material and commands necessary to use the WARP3D finite element code. WARP3D is under continuing development as a research code for the solution of very large-scale, 3-D solid models subjected to static and dynamic loads. Specific features in the code oriented toward the investigation of ductile fracture in metals include a robust finite strain formulation, a general J-integral computation facility (with inertia, face loading), an element extinction facility to model crack growth, nonlinear material models including viscoplastic effects, and the Gurson-Tver-gaard dilatant plasticity model for void growth. The nonlinear, dynamic equilibrium equations are solved using an incremental-iterative, implicit formulation with full Newton iterations to eliminate residual nodal forces. The history integration of the nonlinear equations of motion is accomplished with Newmarks Beta method. A central feature of WARP3D involves the use of a linear-preconditioned conjugate gradient (LPCG) solver implemented in an element-by-element format to replace a conventional direct linear equation solver. This software architecture dramatically reduces both the memory requirements and CPU time for very large, nonlinear solid models since formation of the assembled (dynamic) stiffness matrix is avoided. Analyses thus exhibit the numerical stability for large time (load) steps provided by the implicit formulation coupled with the low memory requirements characteristic of an explicit code. In addition to the much lower memory requirements of the LPCG solver, the CPU time required for solution of the linear equations during each Newton iteration is generally one-half or less of the CPU time required for a traditional direct solver. All other computational aspects of the code (element stiffnesses, element strains, stress updating, element internal forces) are implemented in the element-by- element, blocked architecture. This greatly improves vectorization of the code on uni-processor hardware and enables straightforward parallel-vector processing of element blocks on multi-processor hardware.

Koppenhoefer, Kyle C.↗
NASA NTRS1988

Modelling the slow climatic attractor

Attempts to model the evolution of the slow-response parts of the climatic system such as the global ice mass are reviewed, and the fundamental laws governing the complete complex multidomain slow system, called the slow climatic attractor, are analyzed. The uses of paleoclimatic records are discussed, and a low-order model for the late Quaternary variations that constitutes a credible global representation of the slow climatic attractor is presented. It is shown that although according to the model the climatic system is centered on an unstable equilibrium, the points of the limit cycle attractor are quite stable and robust. It is also noted that one major deficiency of the model is its inability to portray the geographic spatial variability of climatic mean states. However, at least with regard to atmospheric CO2, estimates of the geographic fields of all significant atmospheric and surface boundary layer climate variables are deducible.

Saltzman, Barry↗
NASA NTRS2017

3D Material Response Analysis of PICA Pyrolysis Experiments

Primarily interested in improving ablation modeling for use in inverse reconstruction of flight environments on ablative heat shields. Ablation model is essentially a component of the heat flux sensor, so model uncertainties lead to measurement uncertainties. Non-equilibrium processes have been known to be significant in low density ablators for a long time, but increased accuracy requirements of the reconstruction process necessitates incorporating this physical effect. Attempting to develop a pyrolysis model for implementation in material response based on the PICA data produced by Bessire and Minton. Pyrolysis gas species molar yields as a function of temperature and heating rate. Several problems encountered while trying to fit Arrhenius models to the data led to further investigation of the experimental setup.

Oliver, Brandon A.↗
NASA NTRS1992

A quasi-static magnetospheric convection model in two-dimensions

A self-consistent 2D model of subsonic plasma sheet convection is presented. Specifically, time sequences of static equilibrium solutions for the 2D magnetospheric magnetic field are constructed consistent with adiabatic convection. This model self-consistently includes a dipole field and a reasonable accounting for the effects of inner magnetospheric shielding. Starting from a relaxed magnetospheric equilibrium, the earthward convection of plasma sheet flux tubes results in the stretching of inner plasma sheet field lines, the development of a local minimum in the equatorial magnetic field Be in the near-earth plasma sheet, and an increasing lobe magnetic field. This evolution in time occurs independent of the specific magnetopause or far-tail boundary conditions, provided the plasma sheet is not at marginal interchange stability. This behavior results solely from the convection of flux tubes of increasing entropy into the near-earth plasma sheet. These results are discussed in the general context of earthward convection in earth's plasma sheet and in the specific context of magnetospheric substorms.

Erickson, G. M.↗
NASA NTRS1981

The Voyager 1 EUV spectrum of the Io plasma torus

The Voyager 1 EUV spectrum of the hot Io plasma torus obtained before encounter in 1979 at a range of 4,800,000 km is analyzed in terms of a collisional model with improved collision strengths. Model calculations are based entirely on steady-state equilibrium between electron excitation-deactivation and radiative decay in the subspecies structure, and plasma temperature and species number densities are derived by fitting model spectra to the observations. Lines of S II, S III, S IV, O II, and O III are found to be dominant in the plasma torus, at torus half width densities of 44, 160, 220, 50, and 340/cu cm, respectively, and lines of S V and K III may also be present. An effective electron temperature of 80,000 K is estimated for the central dense region near 6 Jupiter radii. Model calculations are in substantial agreement with other spectroscopic observations of the hot torus; however, the observations exhibit strong morphological differences with recent ground-based observations of S II emission, and calculated ion densities differ from in situ Voyager 1 plasma science measurements.

Shemansky, D. E.↗
NASA NTRSDate not supplied

HCN emissions from the explosive volcanic eruption of Mt. Pinatubo, Philippines, in June 1991

In June 1991, Mt. Pinatubo explosively erupted magma and overlying rock, with ejecta reaching stratospheric altitudes. The sulphate aerosol burden resulting from the 1991 eruption was still measurable in 1992 and its complete removal took 5-6 years. Here we present new analyses of measurements by the ATMOS instrument on the Atlas-1 Space Shuttle mission and the MkIV balloon interferometer in 1992 that show about 40-50 ppt of unexpected excess HCN in the middle stratosphere. HCN has no stratospheric sources and a lifetime of several years. Multi-year runs with the GEOS-Chem Chemistry-Transport model have been performed in order to derive the most plausible HCN injection amounts. The overall structure of the modeled HCN agrees reasonably well with observations by MkIV and ATMOS. The current best estimate is an emission of about 6 kt of HCN into the stratosphere at an altitude of about 23 km. This compares to a total gas emission in the order of 70-700 Mt, based on an emission of 14 Mt SO2 and an SO2/total gas ratio of 2% - 20%. We will discuss possible sources and formation mechanisms of the HCN. These include abiotic formation during or prior to the eruption and may involve HCN accumulated in the subaerial volcanic-hydrothermal system of Mt. Pinatubo. We present results from thermochemical equilibrium calculations as well as photochemical plume modeling in order to provide plausibility constraints on the source of the HCN.

Armin Kleinboehl↗
NASA NTRS1995

Fully-Explicit and Self-Consistent Algebraic Reynolds Stress Models

A fully-explicit, self-consistent algebraic expression for the Reynolds stress, which is the exact solution to the Reynolds stress transport equation in the 'weak equilibrium' limit for two-dimensional mean flows for all linear and some quasi-linear pressure-strain models, is derived. Current explicit algebraic Reynolds stress models derived by employing the 'weak equilibrium' assumption treat the production-to-dissipation (P/epsilon) ratio implicitly, resulting in an effective viscosity that can be singular away from the equilibrium limit. In the present paper, the set of simultaneous algebraic Reynolds stress equations are solved in the full non-linear form and the eddy viscosity is found to be non-singular. Preliminary tests indicate that the model performs adequately, even for three dimensional mean flow cases. Due to the explicit and non-singular nature of the effective viscosity, this model should mitigate many of the difficulties encountered in computing complex turbulent flows with the algebraic Reynolds stress models.

Girimaji, Sharath S.↗
NASA NTRS1982

Simple phenomenological modeling of transition-region capacitance of forward-biased p-n junction diodes and transistor diodes

The derivation of a simple expression for the capacitance C(V) associated with the transition region of a p-n junction under a forward bias is derived by phenomenological reasoning. The treatment of C(V) is based on the conventional Shockley equations, and simpler expressions for C(V) result that are in general accord with the previous analytical and numerical results. C(V) consists of two components resulting from changes in majority carrier concentration and from free hole and electron accumulation in the space-charge region. The space-charge region is conceived as the intrinsic region of an n-i-p structure for a space-charge region markedly wider than the extrinsic Debye lengths at its edges. This region is excited in the sense that the forward bias creates hole and electron densities orders of magnitude larger than those in equilibrium. The recent Shirts-Gordon (1979) modeling of the space-charge region using a dielectric response function is contrasted with the more conventional Schottky-Shockley modeling.

Lindholm, F. A.↗
NASA NTRS1995

Assimilation of Altimeter Data into a Quasigeostrophic Model of the Gulf Stream System: Dynamical Considerations - Part 1

The dynamical consequences of constraining a numerical model with sea surface height data have been investigated. The model used for this study is a quasigeostrophic model of the Gulf Stream region. The data that have been assimilated are maps of sea surface height obtained as the superposition of sea surface height variability deduced from the Geosat altimeter measurements and a mean field constructed from historical hydrographic data. The method used for assimilating the data is the nudging technique. Nudging has been implemented in such a way as to achieve a high degree of convergence of the surface model fields toward the observations. The assimilation of the surface data is thus equivalent to the prescription of a surface pressure boundary condition. The authors analyzed the mechanisms of the model adjustment and the characteristics of the resultant equilibrium state when the surface data are assimilated. Since the surface data are the superposition of a mean component and an eddy component, in order to understand the relative role of these two components in determining the characteristics of the final equilibrium state, two different experiments have been considered: in the first experiment only the climatological mean field is assimilated, while in the second experiment the total surface streamfunction field (mean plus eddies) has been used. It is shown that the model behavior in the presence of the surface data constraint can be conveniently described in terms of baroclinic Fofonoff modes. The prescribed mean component of the surface data acts as a 'surface topography' in this problem. Its presence determines a distortion of the geostrophic contours in the subsurface layers, thus constraining the mean circulation in those layers. The intensity of the mean flow is determined by the inflow/outflow conditions at the open boundaries, as well as by eddy forcing and dissipation.

Capotondi, Antonietta↗
NASA NTRS1975

Calculation of three-dimensional compressible laminar and turbulent boundary layers. An implicit finite-difference procedure for solving the three-dimensional compressible laminar, transitional, and turbulent boundary-layer equations

An implicit finite-difference procedure is presented for solving the compressible three-dimensional boundary-layer equations. The method is second-order accurate, unconditionally stable (conditional stability for reverse cross flow), and efficient from the viewpoint of computer storage and processing time. The Reynolds stress terms are modeled by (1) a single-layer mixing length model and (2) a two-layer eddy viscosity model. These models, although simple in concept, accurately predicted the equilibrium turbulent flow for the conditions considered. Numerical results are compared with experimental wall and profile data for a cone at an angle of attack larger than the cone semiapex angle. These comparisons clearly indicate that the numerical procedure and turbulence models accurately predict the experimental data with as few as 21 nodal points in the plane normal to the wall boundary.

Harris, J. E.↗
NASA NTRS2022

Constraining Clouds and Convective Parameterizations in A Climate Model Using Paleoclimate Data

Cloud and convective parameterizations strongly influence uncertainties in equilibrium climate sensitivity. We provide a proof-of-concept study to constrain these parameterizations in a perturbed parameter ensemble of the atmosphere-only version of the Goddard Institute for Space Studies Model E2.1 simulations by evaluating model biases in the present-day runs using multiple satellite climatologies and by comparing simulated δ 18 O of precipitation (δ 18 O p ), known to be sensitive to parameterization schemes, with a global database of speleothem δ 18 O records covering the Last Glacial Maximum (LGM), mid-Holocene (MH) and pre-industrial (PI) periods. Relative to modern interannual variability, paleoclimate simulations show greater sensitivity to parameter changes, allowing for an evaluation of model uncertainties over a broader range of climate forcing and the identification of parts of the world that are parameter sensitive. Certain simulations reproduced absolute δ 18 O p values across all time periods, along with LGM and MH δ 18 O p anomalies relative to the PI, better than the default parameterization. No single set of parameterizations worked well in all climate states, likely due to the non-stationarity of cloud feedbacks under varying boundary conditions. Future work that involves varying multiple parameter sets simultaneously with coupled ocean feedbacks will likely provide improved constraints on cloud and convective parameterizations.

cloud and convective parameterization↗
NASA NTRS2015

Spin-Up and Tuning of the Global Carbon Cycle Model Inside the GISS ModelE2 GCM

Planetary carbon cycle involves multiple phenomena, acting at variety of temporal and spacial scales. The typical times range from minutes for leaf stomata physiology to centuries for passive soil carbon pools and deep ocean layers. So, finding a satisfactory equilibrium state becomes a challenging and computationally expensive task. Here we present the spin-up processes for different configurations of the GISS Carbon Cycle model from the model forced with MODIS observed Leaf Area Index (LAI) and prescribed ocean to the prognostic LAI and to the model fully coupled to the dynamic ocean and ocean biology. We investigate the time it takes the model to reach the equilibrium and discuss the ways to speed up this process. NASA Goddard Institute for Space Studies General Circulation Model (GISS ModelE2) is currently equipped with all major algorithms necessary for the simulation of the Global Carbon Cycle. The terrestrial part is presented by Ent Terrestrial Biosphere Model (Ent TBM), which includes leaf biophysics, prognostic phenology and soil biogeochemistry module (based on Carnegie-Ames-Stanford model). The ocean part is based on the NASA Ocean Biogeochemistry Model (NOBM). The transport of atmospheric CO2 is performed by the atmospheric part of ModelE2, which employs quadratic upstream algorithm for this purpose.

Climate model↗
NASA NTRS1975

Finite-rate chemistry effects upon convective and radiative heating of an atmospheric entry vehicle

A mathematical model of the aerothermochemical environment along the stagnation line of a planetary return spacecraft using an ablative thermal protection system was developed and solved for conditions typical of atmospheric entry from planetary missions. The model, implemented as a FORTRAN 4 computer program, was designed to predict viscous, reactive and radiative coupled shock layer structure and the resulting body heating rates. The analysis includes flow field coupling with the ablator surface, binary diffusion, coupled line and continuum radiative and equilibrium or finite rate chemistry effects. The gas model used includes thermodynamic, transport, kinetic and radiative properties of air and ablation product species, including 19 chemical species and 16 chemical reactions. Specifically, the impact of nonequilibrium chemistry effects upon stagnation line shock layer structure and body heating rates was investigated.

Guillermo, P.↗
NASA NTRS1975

Three-dimensional dynamics of scientific balloon systems in response to sudden gust loadings

A mathematical model was developed of the three-dimensional dynamics of a high-altitude scientific research balloon system perturbed from its equilibrium configuration by an arbitrary gust loading. The platform is modelled as a system of four coupled pendula, and the equations of motion were developed in the Lagrangian formalism assuming a small-angle approximation. Three-dimensional pendulation, torsion, and precessional motion due to Coriolis forces are considered. Aerodynamic and viscous damping effects on the pendulatory and torsional motions are included. A general model of the gust field incident upon the balloon system was developed. The digital computer simulation program is described, and a guide to its use is given.

Dorsey, D. R., Jr.↗
NASA NTRS1980

Volatile element chemistry in the solar nebula - Na, K, F, Cl, Br, and P

The results of the most extensive set to date of thermodynamic calculations on the equilibrium chemistry of several hundred compounds of the elements Na, K, F, Cl, Br, and P in a solar composition system are reported. Two extreme models of accretion are investigated. In one extreme complete chemical equilibrium between condensates and gases is maintained because the time scale for accretion is long compared to the time scale for cooling or dissipation of the nebula. Condensates formed in this homogeneous accretion model include several phases such as whitlockite, alkali feldspars, and apatite minerals which are found in chondrites. In the other extreme complete isolation of newly formed condensates from prior condensates and gases occurs due to a time scale for accretion that is short relative to the time required for nebular cooling or dissipation. The condensates produced in this heterogeneous accretion model include alkali sulfides, ammonium halides, and ammonium phosphates. None of these phases are found in chondrites. Available observations of the Na, K, F, Cl, Br, and P elemental abundances in the terrestrial planets are found to be compatible with the predictions of the homogeneous accretion model.

Fegley, B., Jr.↗
NASA NTRS1978

Self-acting shaft seals

Self-acting seals are described in detail. The mathematical models for obtaining a seal force balance and the equilibrium operating film thickness are outlined. Particular attention is given to primary ring response (seal vibration) to rotating seat face runout. This response analysis reveals three different vibration models with secondary seal friction being an important parameter. Leakage flow inlet pressure drop and affects of axisymmetric sealing face deformations are discussed. Experimental data on self-acting face seals operating under simulated gas turbine conditions are given. Also a spiral groove seal design operated to 244 m/sec (800 ft/sec) is described.

Ludwig, L. P.↗
NASA NTRS2008

Models of Coronal Loops and Those Annoying Observations

It is often said that coronal loops are the fundamental building blocks of the magnetically-closed corona. This is certainly true, especially when one realizes that the diffuse component of the corona can be thought of as a collection of indistinguishable loops. During the Skylab era, it was believed that loops are in states of quasi static equilibrium. However, more recent observations combined with numerical models have revealed that this interpretation is not correct, at least for many loops. The concept of loops as bundles of impulsively-heated strands is very appealing and solves many of the problems faced by static models, but it too faces observational challenges. Yet another idea involves the fascinating phenomenon of thermal nonequilibrium. I will review the various attempts to model coronal loops and discuss how they agree and disagree with observations, paying particular attention to the latest results from Hinode.

Klimchuk, James A.↗