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,549 records · Page 86

Computer Simulation of Motion of Heliospheric Termination Shock

A time-dependent computer model is used to study a non-equilibrium structure of the heliosphere, resulting from the interaction between the solar wind and the interstellar plasma. An interaction of the heliospheric termination shock with various disturbances such as density or velocity jumps or interplanetary shocks in upstream solar wind as well as pressure jumps in local interstellar medium is investigated. Consequences of latitudinal variation in dynamic pressure are also discussed. The initial boundary conditions of unsteady calculations are given by the solar wind parameters as functions of time and/or latitude on an inner boundary, and/or by the local interstellar pressure as a function of time (roughly simulating the effect of the local interstellar medium) on an outer boundary.

Ratkiewicz, R. E.↗

Non-Invasive Investigation of Bone Adaptation in Humans to Cumulative Daily Mechanical Loading

The goal of our research is to better understand the functional relationship between cumulative daily skeletal loading generated by daily activity and the regulation of bone density and bone structure. We have proposed the calcaneus and tibia as useful model bone sites loaded by internal forces in equilibrium with the ground reaction force during gait. The daily history of the ground reaction force is a good relative measure of daily lower limb and calcaneal loading that can be compared to bone density and structure of the calcaneus and cross-sectional geometry of the tibia and fibula. Over the past several years, we have developed image-processing technologies to improve our ability to measure bone density and structure in the calcaneus and lower leg non-invasively with computed tomography and bone densitometry, or DXA. The objective of our current research effort is to determine the accuracy and precision of our CT and DXA image processing methods.

Whalen, Robert↗

The Long-Term Performance of Small-Cell Batteries Without Cell-Balancing Electronics

Tests approx.8 yrs ago showed Sony HC do not imbalance. AEA developed a theory (ESPC 2002): a) Self-discharge (SD) decreases with state-of-charge (SOC); b) Cells diverge to a state of dynamic equilibrium; c) Equilibrium spread depends on cell SD uniformity. Balancing model verified against test data. Short-term measures of SD difficult in Sony cells and very small values, depends on technique. Long-term evidence supports lower SD at low SD. Battery testing best proof of performance, typically mission specific tests.

Pearson, C.↗

Climate Response to the Increase in Tropospheric Ozone since Preindustrial Times: A Comparison between Ozone and Equivalent CO2 Forcings

We examine the characteristics of the climate response to anthropogenic changes in tropospheric ozone. Using a general circulation model, we have carried out a pair of equilibrium climate simulations with realistic present-day and preindustrial ozone distributions. We find that the instantaneous radiative forcing of 0.49 W m(sup -2) due to the increase in tropospheric ozone since preindustrial times results in an increase in global mean surface temperature of 0.28 C. The increase is nearly 0.4 C in the Northern Hemisphere and about 0.2 C in the Southern Hemisphere. The largest increases (greater than 0.8 C) are downwind of Europe and Asia and over the North American interior in summer. In the lower stratosphere, global mean temperatures decrease by about 0.2 C due to the diminished upward flux of radiation at 9.6 micrometers. The largest stratospheric cooling, up to 1.0 C, occurs over high northern latitudes in winter, with possibly important implications for the formation of polar stratospheric clouds. To identify the characteristics of climate forcing unique to tropospheric ozone, we have conducted two additional climate equilibrium simulations: one in which preindustrial tropospheric ozone concentrations were increased everywhere by 18 ppb, producing the same global radiative forcing as present-day ozone but without the heterogeneity; and one in which CO2 was decreased by 25 ppm relative to present day, with ozone at present-day values, to again produce the same global radiative forcing but with the spectral signature of CO2 rather than ozone. In the first simulation (uniform increase of ozone), the global mean surface temperature increases by 0.25 C, with an interhemispheric difference of only 0.03 C, as compared with nearly 0.2 C for the heterogeneous ozone increase. In the second simulation (equivalent CO2), the global mean surface temperature increases by 0.36 C, 30% higher than the increase from tropospheric ozone. The stronger surface warming from CO2 is in part because CO2 forcing (obscured by water vapor) is shifted relatively poleward where the positive ice-albedo feedback amplifies the climate response and in part because the magnitude of the CO2 forcing in the mid-troposphere is double that of ozone. However, we find that CO2 is far less effective than tropospheric ozone in driving lower stratospheric cooling at high northern latitudes in winter.

Mickley L. J.↗

One Dimensional Analysis Model of a Condensing Spray Chamber Including Rocket Exhaust Using SINDA/FLUINT and CEA

Modeling droplet condensation via CFD codes can be very tedious, time consuming, and inaccurate. CFD codes may be tedious and time consuming in terms of using Lagrangian particle tracking approaches or particle sizing bins. Also since many codes ignore conduction through the droplet and or the degradating effect of heat and mass transfer if noncondensible species are present, the solutions may be inaccurate. The modeling of a condensing spray chamber where the significant size of the water droplets and the time and distance these droplets take to fall, can make the effect of droplet conduction a physical factor that needs to be considered in the model. Furthermore the presence of even a relatively small amount of noncondensible has been shown to reduce the amount of condensation [Ref 1]. It is desirable then to create a modeling tool that addresses these issues. The path taken to create such a tool is illustrated. The application of this tool and subsequent results are based on the spray chamber in the Spacecraft Propulsion Research Facility (B2) located at NASA's Plum Brook Station that tested an RL-10 engine. The platform upon which the condensation physics is modeled is SINDAFLUINT. The use of SINDAFLUINT enables the ability to model various aspects of the entire testing facility, including the rocket exhaust duct flow and heat transfer to the exhaust duct wall. The ejector pumping system of the spray chamber is also easily implemented via SINDAFLUINT. The goal is to create a transient one dimensional flow and heat transfer model beginning at the rocket, continuing through the condensing spray chamber, and finally ending with the ejector pumping system. However the model of the condensing spray chamber may be run independently of the rocket and ejector systems detail, with only appropriate mass flow boundary conditions placed at the entrance and exit of the condensing spray chamber model. The model of the condensing spray chamber takes into account droplet conduction as well as the degrading effect of mass and heat transfer due to the presence of noncondensibles. The one dimension model of the condensing spray chamber makes no presupposition on the pressure profile within the chamber, allowing the implemented droplet physics of heat and mass transfer coupled to the SINDAFLUINT solver to determine a transient pressure profile of the condensing spray chamber. Model results compare well to the RL-10 engine pressure test data.

Condensing Spray Model↗

Spillage of Lunar Polar Crater Volatiles onto Adjacent Terrains: The Case for Dynamic Processes

We present an investigation of the release and transport of lunar polar crater volatiles onto topside regions surrounding the cold traps. The volatiles are liberated via surface energization processes associated with the harsh space environment, including solar wind plasma sputtering and impact vaporization. We find that some fraction of these volatiles can migrate from crater floors onto topside regions (those regions directly adjacent to and above the polar crater floors), and that these surrounding terrains should contain a sampling of the material originating within the crater itself. It is concluded that the nature of the volatile content on crater floors can be obtained by sampling the surface volatiles that have migrated or "spilled out" onto the adjacent terrain. This "spillage" effect could make human or robotic prospecting for crater resources significantly easier, since an assessment may not require direct entry into the very harsh polar crater environment. We also suggest that there are dynamic processes actively operating on the crater floors, and we estimate their source rates assuming dynamic equilibrium of the observed water frost and our modeled loss rates.

Vaporization↗

Turbulence Simulations of Transonic Flows over an NACA-0012 Airfoil

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES),and wall-modeled large-eddy simulation (WMLES)are employed to simulate transonic flow over an NACA-0012 airfoil at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75,and the Reynolds number based on the chord length is 10million. These conditions are the same as the wind-tunnel experimental conditions of McDevitt and Okuno (1985).1The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model.The unsteady RANS simulations,only with the compressibility correction,predict the pre- and post-buffet characteristics,which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment.The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment

Transonic Buffet↗

Experiments and numerical modeling of surrogate droplet combustion for a 93 octane gasoline calibration fuel

Gasoline is the world’s most prolific liquid transportation fuel for ground transportation vehicles. Yet, a fundamental understanding of its burning characteristics is still lacking. Complicating this matter is the highly multicomponent nature of gasoline (comprised of hundreds of chemical species) and the seasonable variability of its composition due to variations in the crude oil feedstocks used in refining. The need to bring order to this complex problem has led to development of calibration fuels (Tier II in presented study) which on average will provide repeatable metrics that replicate performance of refined gasoline for various engine designs. This study presents results of experiments using the configuration of an isolated droplet to study the combustion dynamics of calibration fuel and three different surrogates (S4, S5, S6) for the calibration fuel. Droplet burning is under microgravity conditions, where gas transport is determined entirely by fuel evaporation in order to promote spherical symmetry in the burning process. The initial droplet diameters ranged between 0.55 mm and 0.60 mm. The droplets were spark-ignited and their burning histories were recorded by digital video cameras to obtain quantitative data. For each fuel, the evolutions of droplet, flame, and soot-shell diameters were measured to calculate the burn rates, flame-droplet diameter ratios (FSR) and soot-droplet diameter ratios. Additionally, a detailed numerical modeling is presented. This simulation includes multicomponent phase equilibrium, unsteady gas and liquid transport, radiative heat transfer, internal liquid diffusion and detailed combustion chemistry of different fuels during their burning process. Experiment results show that the calibration fuel and surrogates have a close burn rates and FSR, and S5 presents a relatively similar results with calibration fuel. Predictions from the numerical simulation using the surrogate are compared with experiment measurements, which showed a good agreement with the experiment results, including the droplet burning rate and FSR.

Guo, Songtao↗

Modelling the pole tide and its effect on the earth's rotation

The pole tide is the response of the ocean to incremental centrifugal forces associated with the Chandler wobble. The tide has a potentially important effect on the period and damping of the wobble, but it is at present not well constrained by observations. Here, both analytical and numerical models for the pole tide are constructed. The analytical models consider the tide first in a global ocean and then in an enclosed basin on a beta-plane. The results are found to approach equilibrium linearly with decreasing frequency and inversely with increasing basin depth. The numerical models solve Laplace's tidal equations over the world's oceans using realistic continental boundaries and bottom topography. The results indicate that the effects of non-equilibrium portion of the deep ocean tide on the Chandler wobble period and damping are negligible.

Carton, J. A.↗

Analysis of simultaneous Skylab and ground based flare observations

HeI and HeII resonance line data from Skylab were reduced, analyzed and compared with HeI D3 line intensities taken simultaneously from the Lockheed Rye Canyon Solar Observatory. Computer codes were developed for the calculation of total He line intensities and line profiles from model flare regions. These codes incorporate simultaneous solution of the line and continuum transport equations as needed together with the statistical equilibrium equations for a 30 level HeI, HeII, HeIII system. The energy level model consists of all terms through principal quantum number four. Interpretation of the observed data in terms of these parametric solutions and with simultaneous solution of the transport equations are discussed.

Kulander, J. L.↗

Lattice Boltzmann model for compressible fluids

A lattice Boltzmann model is derived which simulates compressible fluids. By choosing the parameters of the equilibrium distribution appropriately, the sound speed (which may be set arbitrarily low), bulk viscosity, and kinematic viscosity can be selected. This model simulates compressible flows and can include shocks. With a proper rescaling and zero-sound speed, this model simulates Burgers's equation. The viscosity determined by a Chapman-Enskog expansion compares well with that measured form simulations. The exact solutions of Burgers's equation on the unit circle are compared to solutions of lattice Boltzmann model finding reasonable agreement.

Alexander, F. J.↗

SN 1993J: A Type IIb supernova

The evolution of the bright Type II supernova discovered last year in M81, SN 1993J, is consistent with that expected for the explosion of a star which on the main sequence had a mass of 13-16 Solar Mass but which, owing to mass exchange with a binary companion (a intially approximately 3-5 AU, depending upon the actual presupernova radius and the masses of the two stars) lost almost all of its hydrogen-rich envelope during late helium burning. At the time of explosion, the helium core mass was 4.0 +/- 0.5 Solar Mass and the hydrogen envelope, 0.20 +/- 0.05 Solar Mass. The envelope was helium and nitrogen-rich (carbon-deficient) and the radius of the star, 4 +/- 1 x 10(exp 13) cm. The luminosity of the presupernova star was 3 + 1 x 10(exp 38) ergs/s, with the companion star contributing an additional approximately 10(exp 38) ergs/s. The star may have been a pulsating variable at the time of the explosion. For an explosion energy near 10(exp 51) ergs (KE at infinity) and an assumed distance of 3.3 Mpc, a mass of Ni-56 in the range 0.07 +/- 0.01 Solar Mass was produced and ejected. This presciption gives a light curve which compares favorably with the bolomatric observations. Color photometry is more restrictive and requires a model in which the hydrogen-envelope mass is low and the mixing of hydrogen inward has been small, but in which appreciable Ni-56 has been mixed outward into the helium and heavy-element core. It is possible to obtain good agreement with B and V light curves during the first 50 days, but later photometry, especially in bands other than B and V, will require a non-local thermo-dynamic equilibrium (LTE) spectral calculation for comparison. Based upon our model, we predict a flux of approximately 10(exp -5)(3.3 Mpc/D)(exp 2) photons/sq cm/s in the 847 keV line of CO-56 at peak during 1993 August. It may be easier to detect the Computonized continuum which peaks at a few times 10(exp -4) photons /s/sq cm/MeV at 40 keV a few months after the explosion (though neither of these signals were, or should have been, detected by the Compton Gamma-Ray observatory (CGRO). The presupernova star was filing its Roche lobe at the time of the explosion and thus its envelope was highly deformed (about 3:2). The companion star is presently embedded in the supernova, but should becopme visable at age 3 yr (perhaps earlier in the ultraviolet) when the supernova has faded below 10(exp 38) ergs/s. Indeed, if 'kicks' have not played an important role, it is still bound to the neutron star.

Woosley, S. E.↗

Thermal-capillary model for Czochralski growth of semiconductor materials

The success of efficiently calculating the temperature field, crystal radius, melt mensicus, and melt/solid interface in the Czochralski crystal growth system by full finite-element solution of the government thermal-capillary model is demonstrated. The model predicts realistic response to changes in pull rate, melt volume, and the thermal field. The experimentally observed phenomena of interface flipping, bumping, and the difficulty maintaining steady-state growth as the melt depth decreases are explained by model results. These calculations will form the basis for the first quantitative picture of Cz crystal growth. The accurate depiction of the melt meniscus is important in calculating the crystal radius and solidification interface. The sensitivity of the results to the equilibrium growth angle place doubt on less sophisticated attempts to model the process without inclusion of the meniscus. Quantitative comparison with experiments should be possible once more representation of the radiation and view factors in the thermal system and the crucible are included. Extensions of the model in these directions are underway.

Derby, J. J.↗

Small Signal Modelling and Control of the Hydrogen Mixer for Facility E1

We have undertaken the theoretical modelling of an existing liquid hydrogen (LH2) and gas hydrogen (GH2) mixer subsystem of the E1 Ground Test Facility at NASA John C. Stennis Space Center. The E1 test facility carries out comprehensive ground-based testing and certification of various liquid rocket engines and their components. The mixer described in this work is responsible for combining high pressure LH2 and GH2 to produce a hydrogen flow that meets certain thermodynamic properties before it is fed into a test article. The desired properties are maintained by precise control of the mixture of LH2 and GH2 flows. The mixer is modelled as a general multi-flow lumped volume for single constituent fluids using density and internal energy as states. The set of nonlinear differential equations is linearized about an equilibrium point and the resulting two-state, 3-input linear model is analyzed as a possible candidate for control design.

Barbieri, Enrique↗

Small Signal Modeling and Control of the Hydrogen Mixer for Facility E1

We have undertaken the theoretical modelling of an existing liquid hydrogen (LH2) and gas hydrogen (GH2) mixer subsystem of the E1 Ground Test Facility at NASA John C. Stennis Space Center. The E1 test facility carries out comprehensive ground-based testing and certification of various liquid rocket engines and their components. The mixer described in this work is responsible for combining high pressure before it is fed into a test article. The desired properties are maintained by precise control of the mixture of LH2 and GH2 flows. The mixer is modelled as a general multi-flow lumped volume for single constituent fluids using density and internal energy as states. The set of nonlinear differential equations is linearized about an equilibrium point and the resulting two-state, 3-input linear model is analyzed as a possible candidate for control design.

Barbieri, Enrique↗

Global X-Ray Properties of the Vela and Puppis A Supernova Remnants

The Vela and Puppis A supernova remnants(SNRs)comprise a large emission region of∼8°diameter in the soft X-ray sky. The HaloSat CubeSat mission provides the first soft X-ray (0.4–7 keV) observation of the entire Vela SNR and Puppis A SNR region with a single pointing and moderate spectral resolution. HaloSat observations of the Vela SNR are best fit with a two-temperature thermal plasma model consisting of a cooler component with=-+kT0.1910.010.01keV in collisional ionization equilibrium and a hotter component with=-+kT1.0620.270.45 keV in nonequilibrium ionization. Observations of the Puppis A SNR are best fit with a single-component plane-parallel shocked plasma model with=-+kT0.860.050.06 keV in nonequilibrium ionization. For the first time, we find the total X-ray luminosities of both components of the Vela SNR spectrum in the 0.5–7 keV energy band to be= ́-+L4.410X1.41.434 erg s−1 for the cooler component and= ́-+L4.110X1.51.834 erg s−1for the hotter component. We find the total X-ray luminosities of the Vela and Puppis A SNRs to be= ́L8.4 10X34 erg s−1 and= ́-+L6.710X0.91.136 erg s−1

E. M. Silich↗

A new ply model for interlaminar stress analysis

An accurate estimate of interlaminar stresses is crucial to understanding, as well as predicting, many delamination-related failures in composite materials. A new model for ply-level sublaminate analysis is presented and applied. The homogeneous plate theory developed earlier by the authors (Valisetty and Rehfield, 1983) is further refined, and the equations are reduced appropriately for the classical finite-width free-edge laminate elasticity problem and a related delamination crack growth problem. It is applied to the laminate on a ply-by-ply basis. This theory incorporates all the essential physical effects and appears to be an adequate model for predicting the behavior of individual layers in equilibrium. On the basis of the number of equations and boundary conditions required for the implementation of layer equilibrium, this theory also appears to be the simplest of its kind presented so far. The stress induced in the free-edge region of a (0,90,90,0) laminate in uniform extension and the energy release rates for the delamination between the -30 deg and 90 deg plies of a (+, -30,+, -30, 90,90)s laminate are computed using the new analysis. The results are in excellent agreement with the existing numerical solutions. The new ply behavioral model appears to be very promising; it yields stresses and displacements that are statically and kinematically compatible at interlaminar surfaces.

Rao Valisetty, R.↗

Enhancing Grid Resilience with HIVE: Decentralized V2G Coordination for Black Starts

This paper proposes the HIVE (Harmonized Integration of Vehicle Energy for Grid Support) model, a novel game-theoretic framework for decentralized coordination of electrified vehicles to enable black start and load restoration during grid outages. In the absence of a central controller, HIVE employs a cooperative game to model vehicle interactions, allowing autonomous decision-making while admitting to a Nash equilibrium for grid restoration. The framework addresses the heterogeneity of vehicles and their operational constraints, selecting a lead vehicle for grid-forming and coordinating grid-following vehicles to support prioritized loads. Applied to a hospital blackout scenario, HIVE demonstrates robust performance in forming an islanded microgrid and sustaining critical loads under varying vehicle availability, state of charge, and power constraints. Simulation results highlight the model’s effectiveness in ensuring decentralized coordination of energy allocation and prioritizing loads, offering a scalable solution for resilient grid operations.

32 - ENERGY CONSERVATION, CONSUMPTION, AND UTILIZA↗