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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,495 records · Page 83

All sources NASA NTRS DOE OSTI

NASA NTRS1993

The role of large-scale eddies in the climate equilibrium. Part 2: Variable static stability

Lorenz's two-level model on a sphere is used to investigate how the results of Part 1 are modified when the interaction of the vertical eddy heat flux and static stability is included. In general, the climate state does not depend very much on whether or not this interaction is included, because the poleward eddy heat transport dominates the eddy forcing of mean temperature and wind fields. However, the climatic sensitivity is significantly affected. Compared to two-level model results with fixed static stability, the poleward eddy heat flux is less sensitive to the meridional temperature gradient and the gradient is more sensitive to the forcing. For example, the logarithmic derivative of the eddy flux with respect to the gradient has a slope that is reduced from approximately 15 on a beta-plane with fixed static stability and approximately 6 on a sphere with fixed static stability, to approximately 3 to 4 in the present model. This last result is more in line with analyses from observations. The present model also has a stronger baroclinic adjustment than that in Part 1, more like that in two-level beta-plane models with fixed static stability, that is, the midlatitude isentropic slope is very insensitive to the forcing, the diabatic heating, and the friction, unless the forcing is very weak.

Zhou, Shuntai↗

NASA NTRS1995

Non-LTE spectral analysis and model constraints on SN 1993J

We present non-Local Thermodynamic Equilibrium (LTE) synthetic spectra for a time series of observations of SN 1993J obtained on 1993 March 30-31, April 7, April 13-15, and June 13 UT. The spectra are dominated by hydrogen Balmer lines; neutral helium lines, which have been nonthermally excited; and Fe II features. The density profile evolves from an extremely steep 'brick wall' structure with an equivalent power-law index of about 50 on March 30 to a more typical SN II profile with a power law index of about 10. The early spectra are well fitted by a solar composition of metals, although an enhanced abundance of helium is required in order to fit the neutral helium lines. By June 13, the photosphere has receded deep into the helium layer, although there appears to be a layer of hydrogen at higher velocity. The distance is estimated for each epoch. While consistent results are found for spectra obtained in the month of April, the spread in distances from March to June is quite large. Our value for April is mu = 28.0 +/- 0.3 mag, consistent with the recent Cepheid distance to the host galaxy M81. We also compare our results to other implementations of the expanding photosphere method.

Baron, E.↗

NASA NTRS1985

Evidence of solar wind energy deposition into the ionosphere of Mars

Suprathermal electron fluxes measured in the ionosphere of Mars by the retarding potential analyzer (RPA) on Viking lander 1 are presented and compared with the photoelectron flux that is produced by the absorption of the solar EUV. The calculation of the equilibrium photoelectron population on Mars is based on the multistream electron transport theory and a model neutral atmosphere and ionosphere that was actually observed by Viking lander 1. From the theoretical equilibrium photoelectron population, the expected RPA volt-ampere characteristic curves are computed and compared with those recorded by the instrument. The theoretical and the observed RPA currents below about 170 km are in agreement, confirming that the solar EUV is the main source of suprathermal electrons at these altitudes. Above about 170 km an additional source of suprathermal electrons is required to explain the observations.

Mantas, G. P.↗

NASA NTRS1974

Reacting viscous-shock-layer solutions with multicomponent diffusion and mass injection

Numerical solutions of the viscous-shock-layer equation where the chemistry is treated as being either frozen, equilibrium, or nonequilibrium are presented. Also the effects of the diffusion model, surface catalysis, and mass injection on surface transport and flow parameters are considered. The flow is treated as a mixture of five inert and thermally perfect species. The viscous-shock-layer equations are solved by using an implicit-difference scheme. All calculations are for hyperboloids with included angles of 20 and 45. The flight conditions are those for various altitudes and velocities in the earth's atmosphere. Data are presented to show the effects of the chemical models; diffusion models; surface catalysis; and mass injection of air on heat transfer; skin friction; shock standoff distance; wall pressure distribution; and tangential velocity, temperature, and species profiles. The results show that an equilibrium analysis can substantially overpredict the heat-transfer rates for flow conditions experienced by earth-orbital entry vehicles. Moreover, at such conditions surface catalysis significantly influences heat-transfer and flow-field properties. If a binary rather than a multicomponent diffusion model is assumed, negligible errors in most flow properties result. Quantitative results are presented that show the effect of mass injection on flow properties within and downstream of the injection region.

Moss, J. N.↗

NASA NTRS1992

Improved accuracy for finite element structural analysis via a new integrated force method

A comparative study was carried out to determine the accuracy of finite element analyses based on the stiffness method, a mixed method, and the new integrated force and dual integrated force methods. The numerical results were obtained with the following software: MSC/NASTRAN and ASKA for the stiffness method; an MHOST implementation method for the mixed method; and GIFT for the integrated force methods. The results indicate that on an overall basis, the stiffness and mixed methods present some limitations. The stiffness method generally requires a large number of elements in the model to achieve acceptable accuracy. The MHOST method tends to achieve a higher degree of accuracy for course models than does the stiffness method implemented by MSC/NASTRAN and ASKA. The two integrated force methods, which bestow simultaneous emphasis on stress equilibrium and strain compatibility, yield accurate solutions with fewer elements in a model. The full potential of these new integrated force methods remains largely unexploited, and they hold the promise of spawning new finite element structural analysis tools.

Patnaik, Surya N.↗

NASA NTRS1992

Improved accuracy for finite element structural analysis via an integrated force method

Patnaik, S. N.↗

DOE OSTI2025

A quasi-linear model of electromagnetic turbulent transport and its application to flux-driven transport predictions for STEP

A quasi-linear reduced transport model is developed from a database of high-β electromagnetic nonlinear gyrokinetic simulations performed with spherical tokamak for energy production (STEP) relevant parameters. The quasi-linear model is fully electromagnetic and accounts for the effect of equilibrium flow shear using a novel approach. Its flux predictions are shown to agree quantitatively with predictions from local nonlinear gyrokinetic simulations across a broad range of STEP-relevant local equilibria. This reduced transport model is implemented in the T3D transport solver that is used to perform the first flux-driven simulations for STEP to account for transport from hybrid kinetic ballooning mode turbulence, which dominates over a wide region of the core plasma. Nonlinear gyrokinetic simulations of the final transport steady state from T3D return turbulent fluxes that are consistent with the reduced model, indicating that the quasi-linear model may also be appropriate for describing the transport steady state. Within the assumption considered here, our simulations support the existence of a transport steady state in STEP with a fusion power comparable to that in the burning flat top of the conceptual design, but do not demonstrate how this state can be accessed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Model based co-simulation platform for integrated building system control and design optimization

Both steady-state and dynamic simulations have been widely used by HVAC&R industry to support product/equipment development for decades. Steady-state simulation focuses on the system mass, energy and momentum balance of an equilibrium state. It is based on high-fidelity components models, and thus is suitable for system and component design optimization. Dynamic simulation studies the system transient response and is generally used for controls development and verification. It usually does not require rigorous component models of high accuracy because 1) the commonly used PID control is feedback control whose control performance evaluation doesn’t require high fidelity system/plant model; 2) high-fidelity dynamic model significantly increases the number of equations and variables and creates tremendous challenge for math solver. For supervisory control, transactive control or optimization of an integrated building system, the HVAC&R equipment is often one of the sub-components to be controlled. High-fidelity equipment models are required for accurately evaluating control strategies. In addition, building equipment manufacturers have developed a lot of high-fidelity steady-state equipment/component models per their expertise. Thus, a platform that can integrate OEM high-fidelity steady-state model with dynamic building simulation and/or electric power system & grid simulation to support the development and verification of supervisory control for integrated building systems is necessary. In this study, ORNL’s heat pump design tool (HPDM) is utilized to develop a co-simulation platform for supervisory control and optimization in integrated building systems. It is based on a model that integrates high-fidelity steady-state simulation equipment models with dynamic building simulation. A practical case of using the proposed co-simulation platform to develop and evaluate the supervisory control and optimization is presented and discussed.

Sun, Jian↗

DOE OSTI2026

Avoidance of disruptions on KSTAR due to vertical displacement events via novel real-time stability assessment

Disruption avoidance via the DECAF approach has been achieved on KSTAR using a novel real-time vertical stability assessment and a multiactuator feedback control strategy. The development of disruption avoidance strategies with reactor-relevant reliability is an urgent activity, enabling future fusion power plants. The stability metric employed is based on a new formulation of a vertical force gradient balance metric evaluated across the poloidal cross section of the plasma, with parameters tuned using historical data. Evaluation of this metric on a validation set of 400 recent KSTAR shots indicates >82% of Vertical displacement events can be avoided via feedback control. Essential to its calculation is the two-dimensional toroidal current density distribution in the plasma. Measurement of this profile faster than fully-converged equilibrium reconstructions can deliver is found to improve forecaster performance and is achieved with a surrogate model that takes as input magnetic diagnostic measurements and outputs the current profile on a basis comprising the top principal components of historical current profiles (from past equilibrium reconstructions). This method solves the non-uniqueness problem typically faced when reconstructing current profiles directly from diagnostics, while improving computational time and accuracy. On average, profiles produced by this model reach coefficients of determination of >0.99 with respect to those from equilibrium reconstructions. The avoidance actuators employed include poloidal field coils and an electron cyclotron current drive system. The multiactuator approach, as shown in this first demonstration, allows disruption avoidance while minimizing impact to operational performance. This ability, along with its flexibility and speed, makes this new approach an attractive option for avoiding these types of disruptions in reactors.

Tobin, Matthew [Columbia Univ., New York, NY (Unit↗

NASA NTRS1986

Influence of transient baroclinic eddies on planetary-scale waves

Low-order models with a mean zonal flow forced by zonally symmetric heating and one planetary-scale wave forced by topography have multiple flow equilibria. Two stable equilibria are characterized by strong zonal flow with low wave amplitude (high index circulation) and weak zonal flow with a high wave amplitude (low index circulation) fixed in phase with the topography. These two states presumably represent normal zonal circulation and a blocking configuration, respectively. When a shorter, baroclinically unstable wave is dded to this low-order model, the planetary-scale wave no longer stays in stable equilibrium states. Instead, the long wave remains in weather regimes or preferred regions in phase space. These regimes are, in general, different from the equilibria of the model with the planetary-scale wave only. Thus, the short unstable wave adds some randomness to planetary-scale circulation and changes its position with respect to the planetary-scale topography. The Reinhold and Pierrehumbert model is further explored here by adding a long-wave in the thermal forcing and exploring wider parameter space. When the symmetric thermal forcing is weak so that the short wave is stable, the planetary scale wave is fixed in amplitude and phase by the asymmetric forcing. However, when the symmetric forcing is increased to an unstable level, the planetary-scale wave becomes less organized. With asymmetric thermal forcing only (i.e., no topography) the long-wave is randomly distributed through phase space. Thermal forcing seems to be less effective in organizing weather regimes than topographic forcing.

Branscome, L. E.↗

NASA NTRS2008

Climate Response of Direct Radiative Forcing of Anthropogenic Black Carbon

The equilibrium climate effect of direct radiative forcing of anthropogenic black carbon (BC) is examined by 100-year simulations in the Goddard Institute for Space Studies General Circulation Model II-prime coupled to a mixed-layer ocean model. Anthropogenic BC is predicted to raise globally and annually averaged equilibrium surface air temperature by 0.20 K if BC is assumed to be externally mixed. The predicted increase is significantly greater in the Northern Hemisphere (0.29 K) than in the Southern Hemisphere (0.11 K). If BC is assumed to be internally mixed with the present day level of sulfate aerosol, the predicted annual mean surface temperature increase rises to 0.37 K globally, 0.54 K for the Northern Hemisphere, and 0.20 K for the Southern Hemisphere. The climate sensitivity of BC direct radiative forcing is calculated to be 0.6 K W (sup -1) square meters, which is about 70% of that of CO2, independent of the assumption of BC mixing state. The largest surface temperature response occurs over the northern high latitudes during winter and early spring. In the tropics and midlatitudes, the largest temperature increase is predicted to occur in the upper troposphere. Direct radiative forcing of anthropogenic BC is also predicted to lead to a change of precipitation patterns in the tropics; precipitation is predicted to increase between 0 and 20 N and decrease between 0 and 20 S, shifting the intertropical convergence zone northward. If BC is assumed to be internally mixed with sulfate instead of externally mixed, the change in precipitation pattern is enhanced. The change in precipitation pattern is not predicted to alter the global burden of BC significantly because the change occurs predominantly in regions removed from BC sources.

Chung, Serena H.↗

NASA NTRS1971

Laminar or turbulent boundary-layer flows of perfect gases or reacting gas mixtures in chemical equilibrium

Turbulent boundary layer flows of non-reacting gases are predicted for both interal (nozzle) and external flows. Effects of favorable pressure gradients on two eddy viscosity models were studied in rocket and hypervelocity wind tunnel flows. Nozzle flows of equilibrium air with stagnation temperatures up to 10,000 K were computed. Predictions of equilibrium nitrogen flows through hypervelocity nozzles were compared with experimental data. A slender spherically blunted cone was studied at 70,000 ft altitude and 19,000 ft/sec. in the earth's atmosphere. Comparisons with available experimental data showed good agreement. A computer program was developed and fully documented during this investigation for use by interested individuals.

Anderson, E. C.↗

DOE OSTI2023

Modeling Henry's law and phase separations of water–NaCl–organic mixtures with solvation and ion-pairing

Empirical measurements of solution vapor pressure of ternary acetonitrile (MeCN) H 2 O–NaCl–MeCN mixtures were recorded, with NaCl concentrations ranging from zero to the saturation limit, and MeCN concentrations ranging from zero to an absolute mole fraction of 0.64. After accounting for speciation, the variability of the Henry's law coefficient at vapor–liquid equilibrium (VLE) of MeCN ternary mixtures decreased from 107% to 5.1%. Solute speciation was modeled using a mass action solution model that incorporates solute solvation and ion-pairing phenomena. Two empirically determined equilibrium constants corresponding to solute dissociation and ion pairing were utilized for each solute. When speciation effects were considered, the solid–liquid equilibrium of H 2 O–NaCl–MeCN mixtures appear to be governed by a simple saturation equilibrium constant that is consistent with the binary H 2 O–NaCl saturation coefficient. Further, our results indicate that the precipitation of NaCl in the MeCN ternary mixtures was not governed by changes in the dielectric constant. Our model indicates that the compositions of the salt-induced liquid–liquid equilibrium (LLE) boundary of the H 2 O–NaCl–MeCN mixture correspond to the binary plateau activity of MeCN, a range of concentrations over which the activity remains largely invariant in the binary water–MeCN system. Broader comparisons with other ternary miscible organic solvent (MOS) mixtures suggest that salt-induced liquid–liquid equilibrium exists if: (1) the solution displays a positive deviation from the ideal limits governed by Raoult's law; and (2) the minimum of the mixing free energy profile for the binary water-MOS system is organic-rich. Finally, this work is one of the first applications of speciation-based solution models to a ternary system, and the first that includes an organic solute.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2003

Internal Flow Thermal/Fluid Modeling of STS-107 Port Wing in Support of the Columbia Accident Investigation Board

As part of the aero-thermodynamics team supporting the Columbia Accident Investigation Board (CAB), the Marshall Space Flight Center was asked to perform engineering analyses of internal flows in the port wing. The aero-thermodynamics team was split into internal flow and external flow teams with the support being divided between shorter timeframe engineering methods and more complex computational fluid dynamics. In order to gain a rough order of magnitude type of knowledge of the internal flow in the port wing for various breach locations and sizes (as theorized by the CAB to have caused the Columbia re-entry failure), a bulk venting model was required to input boundary flow rates and pressures to the computational fluid dynamics (CFD) analyses. This paper summarizes the modeling that was done by MSFC in Thermal Desktop. A venting model of the entire Orbiter was constructed in FloCAD based on Rockwell International s flight substantiation analyses and the STS-107 reentry trajectory. Chemical equilibrium air thermodynamic properties were generated for SINDA/FLUINT s fluid property routines from a code provided by Langley Research Center. In parallel, a simplified thermal mathematical model of the port wing, including the Thermal Protection System (TPS), was based on more detailed Shuttle re-entry modeling previously done by the Dryden Flight Research Center. Once the venting model was coupled with the thermal model of the wing structure with chemical equilibrium air properties, various breach scenarios were assessed in support of the aero-thermodynamics team. The construction of the coupled model and results are presented herein.

Sharp, John R.↗

DOE OSTI2022

Influence of triangularity on the plasma response to resonant magnetic perturbations

The effect of the axisymmetric plasma shape on the non-axisymmetric plasma response to resonant magnetic perturbations is investigated in experiment and modeling for the DIII-D, EAST and ASDEX Upgrade tokamaks. Systematically modeling the effect of the triangularity, whilst keeping other equilibrium quantities largely unchanged, reveals that the plasma response is strongly suppressed at high triangularity compared to that at low triangularity. This is validated through targeted comparison with experiments at DIII-D, EAST, and ASDEX Upgrade. DIII-D and EAST magnetic measurements are used to validate simulations, while at ASDEX Upgrade the plasma edge displacement is measured. Both experiments and modeling find a reduced magnetic plasma response on the high-field side at high triangularity across devices. Multi-modal analysis of the simulation results extracts the mode structure and applied spectrum dependence of each mode. The amplitude of the dominant mode reveals similar trends with the edge resonance and radial displacement near the X-point, which suggests that the multi-mode plasma response provides another way to understand the edge localized mode (ELM) control physics. The plasma response is strongly reduced at high triangularity compared to that at low triangularity, which implies different ELM control effects as shaping is varied. Furthermore, these findings indicate that the plasma shape should be taken into consideration when designing an RMP-ELM control strategy in experiment, and that predictive plasma response calculations can be used to maximize access to RMP-ELM control in future devices by maximizing the coupling between coils and the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1973

Interaction of pollution abatement with world dynamics

The world dynamics model of Jay W. Forrester was modified to account for pollution abatement. In the modified model, it is assumed that as pollution increases, efforts are made to control pollution. There is a competition between food supply, material standard of living, and pollution abatement for capital, and time is required for diversion of capital toward pollution abatement. Inclusion of pollution abatement in the model drastically alters the response of the world system for the case in which depletion of natural resources is not considered. Instead of undergoing a pollution catastrophe, all system levels move more or less smoothly toward an equilibrium. A FORTRAN program listing of the modified world dynamics model is included.

Smith, G. L.↗

NASA NTRS1975

Experimental aerodynamic characteristics of vehicles traveling in tubes

A simplified theoretical model for a vehicle traveling through an unvented tube under equilibrium incompressible conditions was used to guide the test program, reduce the data, and determine the self-consistency of the results. The results were then used to establish values for the arbitrary coefficients in the theoretical model. Substantial progress was made in understanding the aerodynamic characteristics of vehicles traveling in tubes as exemplified by the good agreement of the theoretical model predictions with the experimental data throughout the Reynolds number range (three orders of magnitude, up to that for an actual full-scale system) and the many geometric variables tested.

Kurtz, D. W.↗

NASA NTRS1990

Liquid-Vapor Equilibrium of Multicomponent Cryogenic Systems

Liquid-vapor and solid-vapor equilibria at low to moderate pressures and low temperatures are important in many solar system environments, including the surface and clouds of Titan, the clouds of Uranus and Neptune, and the surfaces of Mars and Triton. The familiar cases of ideal behavior are limiting cases of a general thermodynamic representation for the vapor pressure of each component in a homogeneous multicomponent system. The fundamental connections of laboratory measurements to thermodynamic models are through the Gibbs-Duhem relation and the Gibbs-Helmholtz relation. Using laboratory measurements of the total pressure, temperature, and compositions of the liquid and vapor phases at equilibrium, the values of these parameters can be determined. The resulting model for vapor-liquid equilibrium can then conveniently and accurately be used to calculate pressures, compositions, condensation altitudes, and their dependencies on changing climatic conditions. A specific system being investigated is CH4-C2H6-N2, at conditions relevant to Titan's surface and atmosphere. Discussed are: the modeling of existing data on CH4-N2, with applications to the composition of Titan's condensate clouds; some new measurements on the CH4-C2H6 binary, using a high-precision static/volumetric system, and on the C2H6-N2 binary, using the volumetric system and a sensitive cryogenic flow calorimeter; and describe a new cryogenic phase-equilibrium vessel with which we are beginning a detailed, systematic study of the three constituent binaries and the ternary CH4-C2H6-N2 system at temperatures ranging from 80 to 105 K and pressures from 0.1 to 7 bar.

Thompson, W. Reid↗