SEARCH · Search NASA

Results for “Hydrodynamic codes”

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 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS2022

Flows in Enthalpy-based Thermal Evolution of Loops

Plasma-filled loop structures are common in the solar corona. Because detailed modeling of the dynamical evolution of these structures is computationally costly, an efficient method for computing approximate but quick physics-based solutions is to rely on space-integrated 0D simulations. The enthalpy-based thermal evolution of loops (EBTEL) framework is a commonly used method to study the exchange of mass and energy between the corona and transition region. EBTEL solves for density, temperature, and pressure, averaged over the coronal part of the loop, velocity at coronal base, and the instantaneous differential emission measure distribution in the transition region. The current single-fluid version of the code, EBTEL2, assumes that at all stages the flows are subsonic. However, sometimes the solutions show the presence of supersonic flows during the impulsive phase of heat input. It is thus necessary to account for this effect. Here, we upgrade EBTEL2 to EBTEL3 by including the kinetic energy term in the Navier–Stokes equation. We compare the solutions from EBTEL3 with those obtained using EBTEL2, as well as the state-of-the-art field-aligned hydrodynamics code HYDRAD. We find that the match in pressure between EBTEL3 and HYDRAD is better than that between EBTEL2 and HYDRAD. Additionally, the velocities predicted by EBTEL3 are in close agreement with those obtained with HYDRAD when the flows are subsonic. However, EBTEL3 solutions deviate substantially from HYDRAD’s when the latter predicts supersonic flows. Using the mismatches in the solution, we propose a criterion to determine the conditions under which EBTEL can be used to study flows in the system.

Solar coronal heating↗

NASA NTRS2003

Foam on Tile Impact Modeling for the Space Shuttle Program

Following the breakup of the Space Shuttle Columbia during reentry a NASA-wide investigation team was formed to examine the probable damage inflicted on Orbiter Thermal Protection System (TPS) elements by impact of External Tank insulating foam projectiles. Our team was to apply rigorous, physics-based analysis techniques to help determine parameters of interest for an experimental test program, utilize validated codes to investigate the full range of impact scenarios, and use analysis derived models to predict aero-thermal-structural responses to entry conditions. We were to operate on a non-interference basis with the j Team, and were to supply significant findings to that team and to the Orbiter Vehicle Engineering Working Group, being responsive to any solicitations for support from these entities. The authors formed a working sub-group within the larger team to apply the Smooth Particle Hydrodynamics code SPHC to the damage estimation problem. Numerical models of the LI-900 TPS tiles and of the BX-250 foam were constructed and used as inputs into the code. Material properties needed to properly model the tiles and foam were obtained from other working sub-groups who performed tests on these items for this purpose. Two- and three- dimensional models of the tiles were constructed, including the glass outer layer, the densified lower layer of LI-900 insulation, the Nomex felt Strain Isolation Pad (SIP) mounting layer, and the underlying aluminum 2024 vehicle skin. A model for the BX-250 foam including porous compression, elastic rebound, and surface erosion was developed. Code results for the tile damage and foam behavior were extensively validated through comparison with the Southwest Research Institute (SwRI) foam-on-tile impact experiments carried out in 1999. These tests involved small projectiles striking individual tiles and small tile arrays. Following code and model validation we simulated impacts of larger ET foam projectiles on the TPS tile systems used on the wings of the orbiter. Tiles used on the Wing Acreage, the Main Landing Gear Door, and the Carrier Panels near the front edge of the wing were modeled. Foam impacts shot for the CAB investigation were modeled, as well as impacts at larger angles, including rapid rotation of the projectile, and with varying foam properties. General results suggest that foam impacts on tiles at about 500 mph could cause appreciable damage if the impact angle is greater than about 20 degrees. Some variations of the foam properties, such as increased brittleness or increased density could increase damage in some cases. Rapid (17 rps) rotation failed to increase the damage for the two cases considered. This does not rule out other cases in which the rotational energy might lead to an increase in tile damage, but suggests that in most cases rotation will not be an important factor. Similar models will be applied for other impacting materials, other velocities, and other geometries as part of the Return to Flight process.

Stellingwerf, R. F.↗

NASA NTRS1975

Hydrodynamic effects in the atmosphere of variable stars

Numerical models of variable stars are established, using a nonlinear radiative transfer coupled hydrodynamics code. The variable Eddington method of radiative transfer is used. Comparisons are for models of W Virginis, beta Doradus, and eta Aquilae. From these models it appears that shocks are formed in the atmospheres of classical Cepheids as well as W Virginis stars. In classical Cepheids, with periods from 7 to 10 days, the bumps occurring in the light and velocity curves appear as the result of a compression wave that reflects from the star's center. At the head of the outward going compression wave, shocks form in the atmosphere. Comparisons between the hydrodynamic motions in W Virginis and classical Cepheids are made. The strong shocks in W Virginis do not penetrate into the interior as do the compression waves formed in classical Cepheids. The shocks formed in W Virginis stars cause emission lines, while in classical Cepheids the shocks are weaker.

Davis, C. G., Jr.↗

NASA NTRS1980

Curvilinear grids for sinuous river channels

In order to effectively analyze the flow in sinuous river channels, a curvilinear grid system was developed for use in the appropriate hydrodynamic code. The CENTERLINE program was designed to generate a two dimensional grid for this purpose. The Cartesian coordinates of a series of points along the boundaries of the sinuous channel represent the primary input to CENTERLINE. The program calculates the location of the river centerline, the distance downstream along the centerline, and both radius of curvature and channel width as a function of such distance downstream. These parameters form the basis for the generation of the curvilinear grid. Based on input values for longitudinal and lateral grid spacing, the corresponding grid system is generated and a file is created containing the appropriate parameters for use in the associated explicit finite difference hydrodynamic programs. Because of the option for a nonuniform grid, grid spacing can be concentrated in areas containing the largest flow gradients.

Tatom, F. B.↗

NASA NTRS2023

The F-CHROMA Grid of 1D RADYN Flare Models

Context. Solar flares are the result of the sudden release of magnetic energy in the corona. Much of this energy goes into accelerating charged particles to high velocity. These particles travel along the magnetic field and the energy is dissipated when the density gets high enough, primarily in the solar chromosphere. Modelling this region is difficult because the radiation energy balance is dominated by strong, optically thick spectral lines. Aims. Our aim is to provide the community with realistic simulations of a flaring loop with an emphasis on the detailed treatment of the chromospheric energy balance. This will enable a detailed comparison of existing and upcoming observations with synthetic observables from the simulations, thereby elucidating the complex interactions in a flaring chromosphere. Methods. We used the 1D radiation hydrodynamics code RADYN to perform simulations of the effect of a beam of electrons injected at the apex of a solar coronal loop. A grid of models was produced, varying the total energy input, the steepness, and low-energy cutoff of the beam energy spectrum. Results. The full simulation results for a grid of models are made available online. Some general properties of the simulations are discussed.

hydrodynamics↗

NASA NTRS1992

Effects of convection electric field on upwelling and escape of ionospheric O(+)

A Monte Carlo code is used to explore the full effects of the convection electric field on distributions of upflowing O(+) ions from the cusp/cleft ionosphere. Trajectories of individual ions/neutrals are computed as they undergo multiple charge-exchange collisions. In the ion state, the trajectories are computed in realistic models of the magnetic field and the convection, corotation, and ambipolar electric fields. The effects of ion-ion collisions are included, and the trajectories are computed with and without simultaneous stochastic heating perpendicular to the magnetic field by a realistic model of broadband, low frequency waves. In the neutral state, ballistic trajectories in the gravitational field are computed. The initial conditions of the ions, in addition to ambipolar electric field and the number densities and temperatures of O(+), H(+), and electrons as a function of height in the cusp/cleft region were obtained from the results of Gombosi and Killeen (1987), who used a hydrodynamic code to simulate the time-dependent frictional-heating effects in a magnetic tube during its motion though the convection throat. The distribution of the ion fluxes as a function of height are constructed from the case histories.

Cladis, J. B.↗

NASA NTRS1974

Carbon detonation supernovae

Arguments aimed at linking main sequence stars in the mass range 4-8 solar masses with type II supernovae and pulsars are presented. In this regard, the problematic evolution of the growing carbon cores before ignition is found to be best rationalized in terms of convective Urca losses. Some calculations on the parameters of carbon detonation models are listed, and the use of a standard numerical hydrodynamic code to deal with these models is discussed.

Buchler, J.-R.↗

NASA NTRS1974

Measurements and calculations of detonation propulsion performance in helium

Previously obtained results are extended to helium ambient gas to help in the assessment of performance in the atmospheres of the major planets. Measurements confirm benefits derived form detonating propellant over conventional chemical rocket propulsion. Benefits however, are reduced at high pressures in low molecular weight gas. Numerical calculations by means of a monodimensional hydrodynamic code follow the different trends obtained experimentally for high and low molecular weight gas and also offer new insights on the time behavior of the process.

Kim, K.↗

NASA NTRS1975

Hydrodynamic models of a Cepheid atmosphere. I - Deep envelope models

The implicit hydrodynamic code of Kutter and Sparks has been modified to include radiative transfer effects. This modified code has been used to compute deep envelope models of a classical Cepheid with a period of 12 days. It is shown that in this particular model the hydrogen ionization region plays only a small role in producing the observed phase lag between the light and velocity curves. The cause of the bumps on the model's light curve is examined, and a mechanism is presented to explain those Cepheids with two secondary features on their light curves. This mechanism is shown to be consistent with the Hertzsprung sequence only if the evolutionary mass-luminosity law is used.

Karp, A. H.↗

NASA NTRS1975

Hydrodynamic models of a Cepheid atmosphere

Instead of computing a large number of coarsely zoned hydrodynamic models covering the entire atmospheric instability strip, the author computed a single model as well as computer limitations allow. The implicit hydrodynamic code of Kutter and Sparks was modified to include radiative transfer effects in optically thin zones.

Karp, A. H.↗

NASA NTRS1976

Numerical simulation of tornado wind loading on structures

A numerical simulation of a tornado interacting with a building was undertaken in order to compare the pressures due to a rotational unsteady wind with that due to steady straight winds used in design of nuclear facilities. The numerical simulations were performed on a two-dimensional compressible hydrodynamics code. Calculated pressure profiles for a typical building were then subjected to a tornado wind field and the results were compared with current quasisteady design calculations. The analysis indicates that current design practices are conservative.

Maiden, D. E.↗

NASA NTRS1978

Time-dependent spherically symmetric accretion onto compact X-ray sources

Analytical arguments and a numerical hydrodynamic code are used to investigate spherically symmetric accretion onto a compact object, in an attempt to provide some insight into gas flows heated by an outgoing X-ray flux. It is shown that preheating of spherically symmetric accretion flows by energetic radiation from an X-ray source results in time-dependent behavior for a much wider range of source parameters than was determined previously and that there are two distinct types of instability. The results are compared with observations of X-ray bursters and transients as well as with theories on quasars and active galactic nuclei that involve quasi-spherically symmetric accretion onto massive black holes. Models based on spherically symmetric accretion are found to be inconsistent with observations of bursters and transients.

Cowie, L. L.↗

NASA NTRS1979

Collisional and dynamical processes in moon and planet formation

The collisional and dynamical processes in moon and planet formation are discussed. A hydrodynamic code of collision calculations, the orbital element changes due to gravitational scattering, a validation of the mass shifting algorithm, a theory of rotations, and the origin of asteroids are studied. A numerical model of planet growth is discussed and a methodology to evaluate the rate at which megaregolith increases its depth as a function of total accumulate number of impacts on an initially smooth, coherent surface is described.

Source record↗

NASA NTRS1979

The hydrodynamics of off-center explosions

The behavior of off-center supernova explosions is investigated using a two-dimensional hydrodynamic code. An important application of these calculations is the possible formation of high-velocity pulsars. The dependence of the final velocity of the collapsed remnant on the location and energy of the explosion is computed. The largest remnant velocities result from explosions located at a mass fraction of 0.5. An explosion energy 50% greater than the binding energy of the star ejects 0.51 solar masses, producing a 1.4 solar mass remnant with a velocity of 400 km/s. However, this energy must be generated in a very small region of the star in order to create the required asymmetry in the explosion. Because of this, a specific energy of about 10 to the 20th ergs/g is needed. Nuclear reactions can produce no more than about 5 x 10 to the 17th erg/g, and it is unclear how the energy produced in gravitational collapse models can be sufficiently localized. Unless a supernova mechanism can be found which can produce enough energy in a small region of the star, off-center explosions do not provide a satisfactory explanation for high-velocity pulsars.

Fryxell, B. A.↗

NASA NTRS1980

Protostellar formation in rotating interstellar clouds. II - Axially symmetric collapse

A two spatial dimension gravitational hydrodynamics code has been used to calculate the initial isothermal dynamic collapse phase of axially symmetric, rotating interstellar clouds. The Eulerian code has been constructed so as to conserve angular momentum both locally (approximately) and globally (exactly). An axially symmetric, rotating cloud collapses to form a rotating, near-equilibrium ring. The rings so formed are compared with those found previously by Black and Bodenheimer, and Bodenheimer and Tscharnuter, and found to agree in structure fairly well. Numerical tests with the code, as well as an analytic calculation of the collapse of a pressureless, rotating, axially symmetric cloud in a fixed gravitational potential, support the assertion that the observed ring formation is physically realistic.

Boss, A. P.↗

NASA NTRS1980

Protostellar formation in rotating interstellar clouds. III - Nonaxisymmetric collapse

The paper discusses a full three spatial-dimension gravitational hydrodynamic code used to follow the collapse of isothermal rotating clouds subjected to various nonaxially symmetric perturbations (NAP). An initially axially symmetric cloud collapsed to form a ring which then fragmented into a binary protostellar system; a low thermal energy cloud with a large bar-shaped NAP collapsed and fragmented into a binary, and higher thermal energy clouds damp out such NAPs while higher rotational energy clouds produce binaries with wider separations. The three-dimensional calculations indicate that isothermal interstellar clouds may fragment into protostellar objects while still in the isothermal regime. Interstellar clouds and their fragments may pass through collapse phases with fragmentation and reduction of spin angular momentum terminating in the formation of pre-main-sequence stars with the observed pre-main-sequence rotation rates.

Boss, A. P.↗

NASA NTRS1980

Time-dependent, optically thick accretion onto a black hole

A fully relativistic hydrodynamics code which incorporates diffusive radiation transport is used to study time-dependent, spherically symmetric, optically thick accretion onto a black hole. It is found that matter free-falls into the hole regardless of whether the diffusion time scale is longer or shorter than the dynamical time. Nonadiabatic heating due to magnetic field reconnection is included. The internal energy thus generated affects the flow in a purely relativistic way, again ensuring free-fall collapse of the inflowing matter. Any matter enveloping a black hole will thus be swallowed on a dynamical time scale with relatively small net release of energy. The inclusion of angular momentum will not necessarily affect this conclusion.

Gilden, D. L.↗

NASA NTRS1980

A numerical study of the NOL large-scale test

The large-scale gap test was modeled numerically using the reactive hydrodynamics code 2DE with forest fire burn rates. The model gives good agreement between calculated and experimental critical gap values. The calibration of peak pressure in the gap versus gap length is obtained from these calculations, and is in good agreement with the published experiment curve. The two dimensional nature of the gap test is evident from the curvature of the shock wave in the gap and test sample, and from the observed distance of run to detonation in the test samples.

Bowan, A. L.↗