Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lagrangian hydrodynamics”

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.

33 records · Page 2

Thermonuclear runaways in thick hydrogen rich envelopes of neutron stars

A Lagrangian, fully implicit, one dimensional hydrodynamic computer code was used to evolve thermonuclear runaways in the accreted hydrogen rich envelopes of 1.0 Msub solar neutron stars with radii of 10 km and 20 km. Simulations produce outbursts which last from about 750 seconds to about one week. Peak effective temeratures and luninosities were 26 million K and 80 thousand Lsub solar for the 10 km study and 5.3 millison and 600 Lsub solar for the 20 km study. Hydrodynamic expansion on the 10 km neutron star produced a precursor lasting about one ten thousandth seconds.

Starrfield, S. G.↗

Thermonuclear processes and accretion onto neutron star envelopes - X-ray burst and transient sources

A Lagrangian, fully implicit, one-dimensional hydrodynamic computer code is used to investigate the evolution of thermonuclear runaways in the thick accreted hydrogen-rich envelopes of 1.0-solar-mass neutron stars having radii of 10 km and 20 km. The simulations produce outbursts that range in time scale from about 2000 seconds to more than a day. For the 10-km study, the peak effective temperature is 3.3 x 10 to the 7th K, and the peak luminosity is 2 x 10 to the 5th solar luminosities. The 20-km neutron star produces a peak effective temperature and luminosity of 5.3 x 10 to the 6th K and 5.9 x 10 to the 2nd solar luminosities. Also investigated are the effects of changes in the rates of the O-14(alpha, proton) and O-15(alpha, gamma) reactions on the evolution. Hydrodynamic expansion on the 10-km neutron star produces a precursor lasting about 10 to the -6th sec. The evolution of a gas cloud impacting the surface of a 20-km, 1-solar-mass neutron star is studied in an attempt to simulate the magnetospheric gate model of the X-ray burst sources.

Starrfield, S.↗

Thermonuclear runaways in thick hydrogen rich envelopes of neutron stars

A Lagrangian, fully implicit, one-dimensional hydrodynamic computer code is used to evolve thermonuclear runaways in the accreted hydrogen-rich envelopes of 1.0-solar-mass neutron stars with radii of 10 km and 20 km. The simulations produce outbursts lasting from approximately 750 seconds to approximately one week. The peak effective temperatures and luminosities are 2.6 x 10 to the 7th K and 8 x 10 to the 4th solar luminosities for the 10 km study and 5.3 x 10 to the 6th K and 600 solar luminosities for the 20 km study. It is found that hydrodynamic expansion on the 10 km neutron star produced a precursor lasting approximately 0.0001 second. The study assumes that the bursters and transient X-ray sources occur as a result of mass transfer from a secondary onto a neutron star in a fashion analogous to the nova phenomena. The peak temperatures and luminosities are found to be inversely proportional to the radius of the neutron stars and the calculations here, together with those in the literature, indicate that the actual radii of most neutron stars must be closer to 10 km than 20 km.

Starrfield, S.↗

Comparison of ALE and SPH Simulations of Vertical Drop Tests of a Composite Fuselage Section into Water

Simulation of multi-terrain impact has been identified as an important research area for improved prediction of rotorcraft crashworthiness within the NASA Subsonic Rotary Wing Aeronautics Program on Rotorcraft Crashworthiness. As part of this effort, two vertical drop tests were conducted of a 5-ft-diameter composite fuselage section into water. For the first test, the fuselage section was impacted in a baseline configuration without energy absorbers. For the second test, the fuselage section was retrofitted with a composite honeycomb energy absorber. Both tests were conducted at a nominal velocity of 25-ft/s. A detailed finite element model was developed to represent each test article and water impact was simulated using both Arbitrary Lagrangian Eulerian (ALE) and Smooth Particle Hydrodynamics (SPH) approaches in LS-DYNA, a nonlinear, explicit transient dynamic finite element code. Analytical predictions were correlated with experimental data for both test configurations. In addition, studies were performed to evaluate the influence of mesh density on test-analysis correlation.

Jackson, Karen E.↗

Dynamics of the baryonic component in hierarchical clustering universes

I present self-consistent 3-D simulations of the formation of virialized systems containing both gas and dark matter in a flat universe. A fully Lagrangian code based on the Smoothed Particle Hydrodynamics technique and a tree data structure has been used to evolve regions of comoving radius 2-3 Mpc. Tidal effects are included by coarse-sampling the density of the outer regions up to a radius approx. 20 Mpc. Initial conditions are set at high redshift (z greater than 7) using a standard Cold Dark Matter perturbation spectrum and a baryon mass fraction of 10 percent (omega(sub b) = 0.1). Simulations in which the gas evolves either adiabatically or radiates energy at a rate determined locally by its cooling function were performed. This allows us to investigate with the same set of simulations the importance of radiative losses in the formation of galaxies and the equilibrium structure of virialized systems where cooling is very inefficient. In the absence of radiative losses, the simulations can be rescaled to the density and radius typical of galaxy clusters. A summary of the main results is presented.

Navarro, Julio↗

Simulation of Hypervelocity Impact on Aluminum-Nextel-Kevlar Orbital Debris Shields

An improved hybrid particle-finite element method has been developed for hypervelocity impact simulation. The method combines the general contact-impact capabilities of particle codes with the true Lagrangian kinematics of large strain finite element formulations. Unlike some alternative schemes which couple Lagrangian finite element models with smooth particle hydrodynamics, the present formulation makes no use of slidelines or penalty forces. The method has been implemented in a parallel, three dimensional computer code. Simulations of three dimensional orbital debris impact problems using this parallel hybrid particle-finite element code, show good agreement with experiment and good speedup in parallel computation. The simulations included single and multi-plate shields as well as aluminum and composite shielding materials. at an impact velocity of eleven kilometers per second.

Fahrenthold, Eric P.↗

Sensitivity of Single-Scattering Properties to Precipitation Particle Meltwater Geometry for GPM Passive Microwave and Radar Remote Sensing Applications

The objective of the current work is to describe the microwave single-scattering properties of partially melted ice-phase precipitation particles using physically-based computational methods, taking into full account the varied geometries of the initially “dry” ice particles and the distributions of liquid water that develop during the melting process. Ultimately, the bulk properties of ensembles of these partially-melted particles will be included in “scattering tables” to support radar and combined radar-radiometer remote sensing of precipitation. A meshless Lagrangian melting procedure (snowMELT) based on smoothed-particle hydrodynamics is applied to both spherical and finely-structured ice particle models to describe the full evolution of the particles from dry ice particles to liquid drops. The melting of spherical ice particles using snowMELT is compared to an alternative continuum physics model to validate melt times and internal thermodynamics. The discrete dipole approximation is then utilized to calculate the single-scattering properties of different mixed-phase particles throughout the melting process. The sensitivities of particle single-scattering properties to thermodynamic and hydrodynamic assumptions in snowMELT are explored, and the implications for combined radar-radiometer precipitation remote sensing from GPM are discussed.

William S Olson↗

Variations on a spherically symmetrical accretion flow

Spherically symmetrical accretion onto a compact gravitating object is considered from a purely hydrodynamical viewpoint while disregarding heating and cooling processes. The study derives the radial perturbation equation, expresses it in a Eulerian form which is simpler than the Lagrangian form considered, for example, by Balazs (1972), and discusses the relation between the two formulations. This perturbation equation is used for the stability analysis of standing waves as well as for the study of short wavelength traveling waves which are a generalization of Parker's (1966) traveling wave solutions in the sense that they are not restricted to isothermal flows. These traveling waves can cause fluctuations in the mass accretion rate to grow large, but will not drive the unperturbed solution away from its assumed stationary profile. The waves conserve the total energy contained in the perturbation.

Petterson, J. A.↗

On the stream-accretion disk interaction - Response to increased mass transfer rate

The time-dependent interaction between the stream of mass from the inner Lagrangian point and the accretion disk, resulting from an increasing mass transfer rate is calculated. The calculation is fully three-dimensional, using a pseudoparticle description of the hydrodynamics. It is demonstrated that the results of such calculations, when combined with specific observations, have the potential of both determining essential parameters, such as the viscosity parameter alpha, and can distinguish between different models of dwarf nova eruptions.

Dgani, Ruth↗

An implicit method for two-dimensional hydrodynamics

An implicit method for compressible multidimensional flows is presented. The method, which is strongly oriented toward astrophysical applications, enables one to simulate very subsonic flows by removing the Courant condition upon time steps. It consists of an implicit purely Lagrangian step, followed by an explicit and second-order accurate (at least in one dimension) remapping step, which is optional. When the remapping step is performed the time step is limited by the 'particle crossing time' and otherwise it is limited only by accuracy considerations. The suggested method, which results from a compromise between accuracy and efficiency, is very efficient relative to other methods. It enables the computation of many multidimensional problems in stellar evolution, such as those governed by very subsonic flows, which were not calculable with existing explicit methods.

Livne, Eli↗

PowderSim: Lagrangian Discrete and Mesh-Free Continuum Simulation Code for Cohesive Soils

PowderSim is a calculation tool that combines a discrete-element method (DEM) module, including calibrated interparticle-interaction relationships, with a mesh-free, continuum, SPH (smoothed-particle hydrodynamics) based module that utilizes enhanced, calibrated, constitutive models capable of mimicking both large deformations and the flow behavior of regolith simulants and lunar regolith under conditions anticipated during in situ resource utilization (ISRU) operations. The major innovation introduced in PowderSim is to use a mesh-free method (SPH-based) with a calibrated and slightly modified critical-state soil mechanics constitutive model to extend the ability of the simulation tool to also address full-scale engineering systems in the continuum sense. The PowderSim software maintains the ability to address particle-scale problems, like size segregation, in selected regions with a traditional DEM module, which has improved contact physics and electrostatic interaction models.

Johnson, Scott↗

Computing supernova collapse to neutron stars and black holes

We present a new numerical code for spherical hydrodynamics in general relativity. The code can handle gravitational collapse to a neutron star or to a black hole without the appearance of singularities. Moreover, the variables and equations in the code are very similar to those appearing in traditional Lagrangian supernova codes. Any such existing code can thus be easily adapted to treat collapse where the final fate is uncertain and may be either a neutron star or a black hole. The code is based on the formulation of Hernandez & Misner, in which retarded time is used as coordinate. This prevents the computational grid from penetrating inside any black hole that may form. We present the equations and a complete finite difference scheme for the adiabatic evolution of a fluid that obeys a gamma-law equation of state. We summarize the results of several testbed calculations performed to check our code. We also give the transformation of the analytic Oppenheimer-Snyder solution for homogeneous dust collapse to our coordinate system.

Baumgarte, Thomas W.↗

Numerical Investigation of Two-Phase Flows With Charged Droplets in Electrostatic Field

A numerical method to solve two-phase turbulent flows with charged droplets in an electrostatic field is presented. The ensemble-averaged Navier-Stokes equations and the electrostatic potential equation are solved using a finite volume method. The transitional turbulence field is described using multiple-time-scale turbulence equations. The equations of motion of droplets are solved using a Lagrangian particle tracking scheme, and the inter-phase momentum exchange is described by the Particle-In-Cell scheme. The electrostatic force caused by an applied electrical potential is calculated using the electrostatic field obtained by solving a Laplacian equation and the force exerted by charged droplets is calculated using the Coulombic force equation. The method is applied to solve electro-hydrodynamic sprays. The calculated droplet velocity distributions for droplet dispersions occurring in a stagnant surrounding are in good agreement with the measured data. For droplet dispersions occurring in a two-phase flow, the droplet trajectories are influenced by aerodynamic forces, the Coulombic force, and the applied electrostatic potential field.

Kim, Sang-Wook↗

Modeling of Water Injection into a Vacuum

A loosely coupled two-phase vacuum water plume model has been developed. This model consists of a source flow model to describe the expansion of water vapor, and the Lagrangian equations of motion for particle trajectories. Gas/Particle interaction is modeled through the drag force induced by the relative velocities. Particles are assumed traveling along streamlines. The equations of motion are integrated to obtain particle velocity along the streamline. This model has been used to predict the mass flux in a 5 meter radius hemispherical domain resulting from the burst of a water jet of 1.5 mm in diameter, mass flow rate of 24.2 g/s, and stagnation pressure of 21.0 psia, which is the nominal Orbiter water dump condition. The result is compared with an empirical water plume model deduced from a video image of the STS-29 water dump. To further improve the model, work has begun to numerically simulate the bubble formation and bursting present in a liquid stream injected into a vacuum. The technique of smoothed particle hydrodynamics was used to formulate this simulation. A status and results of the on-going effort are presented and compared to results from the literature.

Alred, John W.↗

Gaseous Vortices in Barred Spiral Galaxies

During the course of examining many two-dimensional, as well as a smaller sample of three-dimensional, models of gas flows in barred spiral galaxies, we have been impressed by the ubiquitous presence fo vortex pairs, oriented roughly perpendicular to their bars, with one vortex on each side. The vortices are obvious only when viewed in the bar frame, and the centers of their velocity fields usually are near Lagrangian points L(sub 4,5). In all models that we have studied, the vortices form on essentially the same time scale as that for the development of gaseous spiral arms, typically two bar rotations. Usually the corotation radius, r(sub c), lies slightly beyond the end of the bar. Depending upon the mass distributions of the various components, gas spirals either into, or out of, the vortices: In the former case, the vortices become regions of high density, whereas the opposite is true if the gas spirals out of a vortex. The models described in this paper have low-density vortices, as do most of the models we have studied. Moreover, usually the vortex centers lie approximately within +/- 15 deg of L(sub 4,5). In the stellar dynamic limit, when pressure and viscous forces are absent, short-period orbits exist, centered on L(sub 4,5). These orbits need not cross and therefore their morphology is that of gas streamlines, that is, vortices. We believe that the gas vortices in our models are hydrodynamic analogues of closed, short-period, libration orbits centered on L(sub 4,5).

England, Martin N.↗