Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics simulations”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

Anomalous Ion Heating and Pole Erosion in the 12.5-kW Hall Effect Rocket with Magnetic Shielding (HERMeS)

Wear measurements of the inner pole cover in HERMeS have shown as much as ten times higher erosion rates than past numerical simulations predicted. Also, measurements of ion velocity distribution functions (IVDF) suggest that ion heating occurs near the chamfered regions of the acceleration channel and around the inner pole surface. In our companion paper we present theoretical evidence that the plasma conditions in this region can excite instabilities in the lower hybrid frequency. Here, we provide further evidence of their presence through comparisons between numerical simulations and measurements. After including in the 2-D axisymmetric hydrodynamic-PIC code Hall2De an idealized model of the ion temperature based on a scaling with the lower hybrid frequency, we find the agreement between the computed and measured IVDF is significantly improved. In addition, the computed erosion rates at the inner pole surface are now in close agreement with wear test measurements.

Frieman, Jason↗
NASA NTRS2017

Early Hydrodynamic Escape Limits Rocky Planets to Less Than or Equal to 1.6 Earth Radii

In the past decade thousands of exoplanet candidates and hundreds of confirmed exoplanets have been found. For sub-Neptune-sized planets, those less than approx. 10 Earth masses, we can separate planets into two broad categories: predominantly rocky planets, and gaseous planets with thick volatile sheaths. Observations and subsequent analysis of these planets show that rocky planets are only found with radii less than approx. 1.6 Earth radii. No rocky planet has yet been found that violates this limit. We propose that hydrodynamic escape of hydrogen rich protoatmospheres, accreted by forming planets, explains the limit in rocky planet size. Following the hydrodynamic escape model employed by Luger et al. (2015), we modelled the XUV driven escape from young planets (less than approx.100 Myr in age) around a Sun-like star. With a simple, first-order model we found that the rocky planet radii limit occurs consistently at approx. 1.6 Earth radii across a wide range of plausible parameter spaces. Our model shows that hydrodynamic escape can explain the observed cutoff between rocky and gaseous planets. Fig. 1 shows the results of our model for rocky planets between 0.5 and 10 Earth masses that accrete 3 wt. % H2/He during formation. The simulation was run for 100 Myr, after that time the XUV flux drops off exponentially and hydrodynamic escape drops with it. A cutoff between rocky planets and gaseous ones is clearly seen at approx. 1.5-1.6 Earth radii. We are only interested in the upper size limit for rocky planets. As such, we assumed pure hydrogen atmospheres and the highest possible isothermal atmospheric temperatures, which will produce an upper limit on the hydrodynamic loss rate. Previous work shows that a reasonable approximation for an upper temperature limit in a hydrogen rich protoatmosphere is 2000-3000 K, consistent with our assumptions. From these results, we propose that the observed dichotomy between mini-Neptunes and rocky worlds is simply explained by an early episode of thermally-driven hydrodynamic escape when host stars have saturated XUV fluxes.

Lehmer, O. R.↗
NASA NTRS1991

A simulation of the collapse and fragmentation of cooling molecular clouds

The application of the Smoothed Particle Hydrodynamics method to the fragmentation of rotating cloud and disk systems is described, allowing for molecular cooling due to H2 and CO. A novel approach to solving Poisson's equation for disklike structures which exploits the multigrid algorithm is also described. Numerical studies are presented which investigate the evolution of both rotating clouds and Maclaurin disks, in each case with both an isothermal equation of state and with molecular cooling. The results establish the influence of molecular cooling on the fragmentation of molecular clouds. The isothermal sequences, if they fragment at all, do so into far fewer lumps than the cooling sequences. This is not due to a cooling instability as such, but rather to the reduced thermal support. One of the sequences shows a remarkable similarity to the W49A star-forming region.

Monaghan, Joe J.↗
NASA NTRS2013

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↗
NASA NTRS2004

Tethers as Debris: Simulating Impacts of Kevlar Tethers on Shuttle Tiles

In a previous paper I examined the effects of impacts of polymer tethers on aluminum plates using the SPHC hydrodynamic code. In this paper I apply tether models to a new target - models of Space Shuttle tiles developed during the STS 107 accident investigation. In this three-dimensional simulation, a short tether fragment strikes a single tile supported on an aluminum backing plate. A tile of the LI-900 material is modeled. Penetration and damage to the tile and the backwall are characterized for three normal impact velocities. The tether is modeled as a bundle of eight 1-mm strands, with the bundle having dimensions 2-mm x 4-mm x 20-cm. The bulk material properties used are those of Kevlar(TradeMark) 49, for which a Mie-Gruneisen multiphase equation of state (eos) is used. In addition, the strength model is applied in a linear sense, such that tensile loads along the strand length are supported, but there is no strength in the lateral directions. Tile models include the various layers making up the tile structure. The outermost layer is a relatively dense borosilicate glass, known as RCG, 0.5-mm thick. The RCG layer is present on the top and four sides of the tile. Below this coating is the bulk of the tile, 1.8- in thick, made of LI-900, a product consisting of rigidized fiberous silica with a density of 9 lWft3. Below the main insulating layer is a bottom layer of the same material that has been treated to increase its density by approximately 69% to improve its strength. This densified layer is bonded to a Strain Isolation Pad (SIP), modeled as a refractory felt fabric. The SIP is bonded to an aluminum 2024 wall 0.1-in thick. The tile and backwall materials use a Me-Gruneisen multiphase eos, with the exception of the SIP felt, which uses a fabric equation of state. Fabrics must be crushed to the full bulk material density before bulk material properties and a Mie-Gruneisen eos are applied. Tether fragment impact speeds of 3,7, and 10 km/s are simulated, with impact velocities normal to the tile face. Damage results are presented in tabular format.

Evans, Steven W.↗
NASA NTRS1985

Shock acceleration of diffuse ions at the earth's bow shock Acceleration efficiency and A/Z enhancement

Observations of particle spectra, intensity, and enhancement of alpha particles over protons at diffuse ion events at the quasi-parallel earth bow shock are compared to a Monte Carlo simulation of diffusive shock acceleration. The simulation includes the back reaction of accelerated particles on the shock structure, particle escape at an upstream free escape boundary, and a low energy per nucleon threshold for thermal leakage of downstream, shock-heated particles into the upstream region. The simulation assumes that the same scattering operator that gives rise to shock acceleration can also describe a viscous shock governed by hydrodynamic turbulence. This implies that accelerated ions can be drawn directly from the thermal solar wind with no separate superthermal seed population. Good agreement between the simulation and observations made during nearly radial magnetic field configurations lends support to thermal leakage of downstream, shock-heated ions as the mode of injection for diffusion ion events.

Ellison, D. C.↗
NASA NTRS1981

Numerical simulation of the magnetospheric gate model for X-ray bursters

A Lagrangian, fully implicit, one dimensional hydrodynamic computer code was used to investigate the evolution of a gas cloud impacting the surface of a 20 km, 1 Msub solar neutron star. This gas is initially at rest with respect to the surface of the neutron star, extends to 185 km above the surface, and is optically thick. The infall results in a burst which lasts about 0.1 seconds and reached a peak luminosity and effective temperature of 240,000 Lsub solar and 9 million; respectively. The burst was followed by a phase of oscillations with a period 0.2 seconds.

Starrfield, S. G.↗
NASA NTRS2008

An Axisymmetric, Hydrodynamical Model for the Torus Wind in Active Galactic Nuclei

We report on time-dependent axisymmetric simulations of an X-ray-excited flow from a parsec-scale, rotating, cold torus around an active galactic nucleus. Our simulations account for radiative heating and cooling and radiation pressure force. The simulations follow the development of a broad biconical outflow induced mainly by X-ray heating. We compute synthetic spectra predicted by our simulations. The wind characteristics and the spectra support the hypothesis that a rotationally supported torus can serve as the source of a wind which is responsible for the warm absorber gas observed in the X-ray spectra of many Seyfert galaxies.

Dorodnitsyn, A.↗
NASA NTRS1991

Preferential concentration of particles by turbulence

Direct numerical simulation of isotropic turbulence was used to investigate the effect of turbulence on the concentration fields of heavy particles. The hydrodynamic field was computed using 64-cubed points, and a statistically stationary flow was obtained by forcing the low-wave-number components of the velocity field. The particles used in the simulations were time advanced according to Stokes drag law and were also assumed to be much more dense than the fluid. Properties of the particle cloud were obtained by following the trajectories of 1,000,000 particles through the simulated flow fields. Three values of the ratio of the particle time constant to large-scale turbulence time scale were used in the simulations: 0.075, 0.15, and 0.52. The simulations show that the particles collect preferentially in regions of low vorticity and high strain rate. This preferential collection was most pronounced for particles with the instantaneous number density as much as 25 times the mean value. The fact that dense particles collect in regions of low vorticity and high strain in turn implies that turbulence may actually inhibit rather than enhance mixing of particles.

Squires, Kyle D.↗
NASA NTRS2000

Theory of the lattice Boltzmann Method: Dispersion, Dissipation, Isotropy, Galilean Invariance, and Stability

The generalized hydrodynamics (the wave vector dependence of the transport coefficients) of a generalized lattice Boltzmann equation (LBE) is studied in detail. The generalized lattice Boltzmann equation is constructed in moment space rather than in discrete velocity space. The generalized hydrodynamics of the model is obtained by solving the dispersion equation of the linearized LBE either analytically by using perturbation technique or numerically. The proposed LBE model has a maximum number of adjustable parameters for the given set of discrete velocities. Generalized hydrodynamics characterizes dispersion, dissipation (hyper-viscosities), anisotropy, and lack of Galilean invariance of the model, and can be applied to select the values of the adjustable parameters which optimize the properties of the model. The proposed generalized hydrodynamic analysis also provides some insights into stability and proper initial conditions for LBE simulations. The stability properties of some 2D LBE models are analyzed and compared with each other in the parameter space of the mean streaming velocity and the viscous relaxation time. The procedure described in this work can be applied to analyze other LBE models. As examples, LBE models with various interpolation schemes are analyzed. Numerical results on shear flow with an initially discontinuous velocity profile (shock) with or without a constant streaming velocity are shown to demonstrate the dispersion effects in the LBE model; the results compare favorably with our theoretical analysis. We also show that whereas linear analysis of the LBE evolution operator is equivalent to Chapman-Enskog analysis in the long wave-length limit (wave vector k = 0), it can also provide results for large values of k. Such results are important for the stability and other hydrodynamic properties of the LBE method and cannot be obtained through Chapman-Enskog analysis.

Lallemand, Pierre↗
NASA NTRS2004

High-Order Simulation of Non-Linear Oscillations and Shocks in the Solar Atmosphere

The solar atmosphere presents a rich source of highly non-linear magneto-hydrodynamic phenomena: strong gradients and forcing terms result in both large shocks and oscillations. The additional requirements of energy balance and initialization in hydrostatic equilibrium compound the challenge of this problem. A wealth of observational data allows us to check the results of our simulations. The problem of simulating the solar atmosphere provides, in addition to an interesting system in its own right, a challenging testbed for high-order shock-capturing methods. We discuss the challenge of simulating solar atmospheric phenomena, concentrating on various high-order central methods ranging from second to fourth order. Our method is based on the central-upwind scheme of Kurganov, Noelle and Petrova, which we extend to high order via various interpolants. We investigate various initial data for our simulations, corresponding to observed conditions in different regions of the solar surface: the normal quiet sun and sunspots. When non-oscillatory using second- and third-order methods, we are able to reproduce non-trivial observational results. In particular we find a correlation between initial data and both the shock speeds and particle oscillation spectra that match observations in the corresponding regions. When using fourth-order WENO interpolants, we find that while the individual shock profiles at any given time appear non-oscillatory, spurious oscillations appear in the fields after long time integrations. The issue of initialization in hydrostatic equilibrium raises difficult issues. Careful treatment of the gravitational source term can reduce violations of hydrostatic equilibrium, but difficulties remain, primarily due to discontinuities in the piecewise-polynomial reconstructions. In the case of the Euler equations in gravity, the use of high-order methods reduces the violation of hydrostatic equilibrium to a sufficiently low level for the above described results to be obtained. Including various energy flux terms in the equations, however, leads to severe instability when hydrostatic balance is violated. We discuss various strategies to address this problem, including adaptive grid methods.

Bryson, S.↗
NASA NTRS2009

Simulation of Relativistic Shocks and Associated Radiation from Turbulent Magnetic Fields

Plasma instabilities (e.g., Buneman, Weibel and other two-stream instabilities) excited in collisionless shocks are responsible for particle (electron, positron, and ion) acceleration. Using a new 3-D relativistic particle-in-cell code, we have investigated the particle acceleration and shock structure associated with an unmagnetized relativistic electron-positron jet propagating into an unmagnetized electron-positron plasma. The simulation has been performed using a long simulation system in order to study the nonlinear stages of the Weibel instability, the particle acceleration mechanism, and the shock structure. Cold jet electrons are thermalized and slowed while the ambient electrons are swept up to create a partially developed hydrodynamic (HD) like shock structure. In the leading shock, electron density increases by a factor of <_ 3.5 in the simulation frame. Strong electromagnetic fields are generated in the trailing shock and provide an emission site. We discuss the possible implication of our simulation results within the AGN and GRB context. We have calculated the time evolution of the spectrum from two electrons propagating in a uniform parallel magnetic field to verify the technique. The same technique will be used to calculate radiation from accelerated electrons (positrons) in turbulent magnetic fields generated by Weibel instability.

Nishikawa, K.↗
NASA NTRS1990

Hollow conical jet models for SS 433 - A paradigm lost?

A precessing jet such as that in SS 433 may be approximated as an axisymmetric flow, if the precession time is short by comparison to the propagation time. A series of simulations has been conducted for precessing jets using an axisymmetric finite-difference hydrodynamics code. Examinations are made of hollow cylindrical jets, which lack the complication of a growing interior volume, and conical jets, which model the behavior of a precessing jet propagating on the surface of its precession cone.

Kochanek, Christopher S.↗
NASA NTRS1999

Bubble Dynamics on a Heated Surface

In this work, we study steady and oscillatory thermocapillary and natural convective flows generated by a bubble on a heated solid surface. The interaction between gas and vapor bubbles with the surrounding fluid is of interest for both space and ground-based processing. A combined numerical-experimental approach is adopted here. The temperature field is visualized using Mach-Zehnder and/or Wollaston Prism Interferometry and the flow field is observed by a laser sheet flow visualization technique. A finite element numerical model is developed which solves the transient two-dimensional continuity, momentum, and energy equations and includes the effects of temperature-dependent surface tension and bubble surface deformation. Below the critical Marangoni number, the steady state low-g and 1-g temperature and velocity fields predicted by the finite element model are in excellent agreement with both the visualization experiments in our laboratory and recently published experimental results in the literature. Above the critical Marangoni number, the model predicts an oscillatory flow which is also closely confirmed by experiments. It is shown that the dynamics of the oscillatory flow are directly controlled by the thermal and hydrodynamic interactions brought about by combined natural and thermocapillary convection. Therefore, as numerical simulations show, there are considerable differences between the 1-g and low-g temperature and flow fields at both low and high Marangoni numbers. This has serious implications for both materials processing and fluid management in space.

Kassemi, M.↗
NASA NTRS2004

Multi-Element Unstructured Analyses of Complex Valve Systems

The safe and reliable operation of high pressure test stands for rocket engine and component testing places an increased emphasis on the performance of control valves and flow metering devices. In this paper, we will present a series of high fidelity computational analyses of systems ranging from cryogenic control valves and pressure regulator systems to cavitating venturis that are used to support rocket engine and component testing at NASA Stennis Space Center. A generalized multi-element framework with sub-models for grid adaption, grid movement and multi-phase flow dynamics has been used to carry out the simulations. Such a framework provides the flexibility of resolving the structural and functional complexities that are typically associated with valve-based high pressure feed systems and have been difficult to deal with traditional CFD methods. Our simulations revealed a rich variety of flow phenomena such as secondary flow patterns, hydrodynamic instabilities, fluctuating vapor pockets etc. In the paper, we will discuss performance losses related to cryogenic control valves, and provide insight into the physics of the dominant multi-phase fluid transport phenomena that are responsible for the choking like behavior in cryogenic control elements. Additionally, we will provide detailed analyses of the modal instability that is observed in the operation of the dome pressure regulator valve. Such instabilities are usually not localized and manifest themselves as a system wide phenomena leading to an undesirable chatter at high flow conditions.

Sulyma, Peter↗
NASA NTRS2000

Airbreathing Propulsion System Analysis Using Multithreaded Parallel Processing

In this paper, parallel processing is used to analyze the mixing, and combustion behavior of hypersonic flow. Preliminary work for a sonic transverse hydrogen jet injected from a slot into a Mach 4 airstream in a two-dimensional duct combustor has been completed [Moon and Chung, 1996]. Our aim is to extend this work to three-dimensional domain using multithreaded domain decomposition parallel processing based on the flowfield-dependent variation theory. Numerical simulations of chemically reacting flows are difficult because of the strong interactions between the turbulent hydrodynamic and chemical processes. The algorithm must provide an accurate representation of the flowfield, since unphysical flowfield calculations will lead to the faulty loss or creation of species mass fraction, or even premature ignition, which in turn alters the flowfield information. Another difficulty arises from the disparity in time scales between the flowfield and chemical reactions, which may require the use of finite rate chemistry. The situations are more complex when there is a disparity in length scales involved in turbulence. In order to cope with these complicated physical phenomena, it is our plan to utilize the flowfield-dependent variation theory mentioned above, facilitated by large eddy simulation. Undoubtedly, the proposed computation requires the most sophisticated computational strategies. The multithreaded domain decomposition parallel processing will be necessary in order to reduce both computational time and storage. Without special treatments involved in computer engineering, our attempt to analyze the airbreathing combustion appears to be difficult, if not impossible.

Schunk, Richard Gregory↗
NASA NTRS2003

Water Impact Test and Simulation of a Composite Energy Absorbing Fuselage Section

In March 2002, a 25-ft/s vertical drop test of a composite fuselage section was conducted onto water. The purpose of the test was to obtain experimental data characterizing the structural response of the fuselage section during water impact for comparison with two previous drop tests that were performed onto a rigid surface and soft soil. For the drop test, the fuselage section was configured with ten 100-lb. lead masses, five per side, that were attached to seat rails mounted to the floor. The fuselage section was raised to a height of 10-ft. and dropped vertically into a 15-ft. diameter pool filled to a depth of 3.5-ft. with water. Approximately 70 channels of data were collected during the drop test at a 10-kHz sampling rate. The test data were used to validate crash simulations of the water impact that were developed using the nonlinear, explicit transient dynamic codes, MSC.Dytran and LS-DYNA. The fuselage structure was modeled using shell and solid elements with a Lagrangian mesh, and the water was modeled with both Eulerian and Lagrangian techniques. The fluid-structure interactions were executed using the fast general coupling in MSC.Dytran and the Arbitrary Lagrange-Euler (ALE) coupling in LS-DYNA. Additionally, the smooth particle hydrodynamics (SPH) meshless Lagrangian technique was used in LS-DYNA to represent the fluid. The simulation results were correlated with the test data to validate the modeling approach. Additional simulation studies were performed to determine how changes in mesh density, mesh uniformity, fluid viscosity, and failure strain influence the test-analysis correlation.

Fasanella, Edwin L.↗