Stochastic hydrodynamics of sediment transport
Three dimensional unsteady and steady sediment distribution obtained by Monte Carlo simulation method, using stochastic hydrodynamics model for transport
Chen, K.-C.↗
SEARCH · Search NASA
Results for “hydrodynamic simulation”
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.
Analytic Modeling of the Hydrodynamic, Thermal, and Structural Behavior of Foil Thrust Bearings
A simulation and modeling effort is conducted on gas foil thrust bearings. A foil bearing is a self acting hydrodynamic device capable of separating stationary and rotating components of rotating machinery by a film of air or other gaseous lubricant. Although simple in appearance these bearings have proven to be complicated devices in analysis. They are sensitive to fluid structure interaction, use a compressible gas as a lubricant, may not be in the fully continuum range of fluid mechanics, and operate in the range where viscous heat generation is significant. These factors provide a challenge to the simulation and modeling task. The Reynolds equation with the addition of Knudsen number effects due to thin film thicknesses is used to simulate the hydrodynamics. The energy equation is manipulated to simulate the temperature field of the lubricant film and combined with the ideal gas relationship, provides density field input to the Reynolds equation. Heat transfer between the lubricant and the surroundings is also modeled. The structural deformations of the bearing are modeled with a single partial differential equation. The equation models the top foil as a thin, bending dominated membrane whose deflections are governed by the biharmonic equation. A linear superposition of hydrodynamic load and compliant foundation reaction is included. The stiffness of the compliant foundation is modeled as a distributed stiffness that supports the top foil. The system of governing equations is solved numerically by a computer program written in the Mathematica computing environment. Representative calculations and comparisons with experimental results are included for a generation I gas foil thrust bearing.
Bruckner, Robert J.↗
A New Cell-Centered Implicit Numerical Scheme for Ions in the 2-D Axisymmetric Code Hall2de
We present a new algorithm in the Hall2De code to simulate the ion hydrodynamics in the acceleration channel and near plume regions of Hall-effect thrusters. This implementation constitutes an upgrade of the capabilities built in the Hall2De code. The equations of mass conservation and momentum for unmagnetized ions are solved using a conservative, finite-volume, cell-centered scheme on a magnetic-field-aligned grid. Major computational savings are achieved by making use of an implicit predictor/multi-corrector algorithm for time evolution. Inaccuracies in the prediction of the motion of low-energy ions in the near plume in hydrodynamics approaches are addressed by implementing a multi-fluid algorithm that tracks ions of different energies separately. A wide range of comparisons with measurements are performed to validate the new ion algorithms. Several numerical experiments with the location and value of the anomalous collision frequency are also presented. Differences in the plasma properties in the near-plume between the single fluid and multi-fluid approaches are discussed. We complete our validation by comparing predicted erosion rates at the channel walls of the thruster with measurements. Erosion rates predicted by the plasma properties obtained from simulations replicate accurately measured rates of erosion within the uncertainty range of the sputtering models employed.
electric propulsion↗
Tethers as Debris: Simulating Impacts of Tether Fragments on Shuttle Tiles
The SPHC hydrodynamic code was used to simulate impacts of Kevlar and aluminum projectiles on a model of the LI-900 type insulating tiles used on Space Shuffle Orbiters The intent was to examine likely damage that such tiles might experience if impacted by orbital debris consisting of tether fragments. Projectile speeds ranged from 300 meters per second to 10 kilometers per second. Damage is characterized by penetration depth, tile surface-hole diameter, tile body-cavity diameter, coating fracture diameter, tether and cavity wall material phases, and deformation of the aluminum backwall.
Evans, Steven W.↗
Using Models for How Energetic Electrons Heat the Atmosphere During Flares
Using models for how energetic electrons heat the atmosphere during flares, we simulate the radiative-hydrodynamic response of the lower solar atmosphere to flare heating. The simulations account for much of the non-LTE, optically thick radiative transfer that occurs in the chromosphere. Our models predict an increase in white light continuum during the flare on the order of 20%, but this is highly sensitive to the electron beam flux used in the simulation. We find that a majority of the white light continuum originates in the chromosphere as a result of Balmer and Paschen recombinations, but a significant portion also forms in the photosphere which has been heated by radiative backwarming.
Allred, Joel↗
Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows
High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10 −3 to 10 −1 , which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.
79 ASTRONOMY AND ASTROPHYSICS↗
Simulation of Binary-single Interactions in AGN Disk. I. Gas-enhanced Binary Orbital Hardening
Stellar-mass binary black hole (BBH) mergers within the accretion disks of active galactic nuclei may contribute to gravitational wave (GW) events detected by ground-based GW detectors. In particular, the interaction between a BBH and a single stellar-mass black hole (sBH), known as the binary-single interaction (BSI) process, can potentially lead to GW events with detectable nonzero eccentricity. Previous studies of the BSI process, which neglected the effects of gas, showed that BSIs contribute non-negligibly to GW events in a coplanar disk environment. In this work, we conduct a series of two-dimensional hydrodynamical and N-body simulations to explore the BSI in a gas environment by coupling REBOUND with Athena++. We perform 360 simulation runs, spanning parameters in disk surface density Σ 0 and impact parameter b. We find that the gas-induced energy dissipation within the three-body system becomes significant if the encounter velocity between the sBHs is sufficiently large (≫c s ). Our simulation results indicate that approximately half of the end states of the BSI are changed by gas. Furthermore, at higher gas density, the number of close encounters during the BSI process will increase, and the end-state BBHs tend to be more compact. Consequently, the presence of gas may shorten the GW merger timescale for end-state BBHs and increase the three-body merger rate.
79 ASTRONOMY AND ASTROPHYSICS↗
The Effects of Material-Filled Voids on Detonation Wave Shape in Rubberized RDX Explosives
The sensitivity of explosives is affected by inhomogeneities within the material. This is evident in the increased shock sensitivity of explosives with slightly lower densities resulting from an increased number of hotspots. The influence of hotspots on explosive initiation has been well studied; however, few studies have been conducted on the effect of intermediate-sized voids (0.1–10 mm) on a propagating detonation wave. Cylindrical voids filled with air have been studied for diameters ranging from 0.3 mm to 0.8 mm for both 1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX) and 1,3,5-trinitro-1,3,5-triazinane (RDX)-based rubberized explosives. Continuing the investigation into single cylindrical voids, this study examined the effects of 0.5 mm diameter voids filled with different inert cylindrical metals on the detonation wave shape for an RDX-based rubberized explosive. The metals selected for experiments were 1066 aluminum, brass, copper, and tungsten. The propagation of the detonation wave was captured using a digital streak camera. Experimental results showed that the extent of detonation wave shaping was closely tied to the density differential between the bulk explosive and metal insert. Forty-four different filler materials, including non-metals, were simulated using a hydrodynamic code to further analyze material inclusion effects. The main factors hypothesized to be of interest were bulk sound speed, shock impedance, and filler material density. We found that the local detonation delay could be correlated fairly well to a ratio of bulk sound speed and density. Understanding the influence of material inclusions on detonation performance and wave shape allows for tailoring of detonations.
45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
Development of a coupled experimental–computational approach for engineering optimization of spout-fluidized bed particle coating systems
The design of spout-fluidized bed (SFB) coating systems for nuclear particle fuels typically relies on trial-and-error processes, comprising iterative and time-consuming coating deposition experiments and post-deposition characterization. At an engineering scale, this approach to guided SFB system design is inefficient, highlighting the need for streamlined experimental methodologies which can correlate fluidization conditions to downstream coating outcomes. In this study, we combine time-resolved particle image velocimetry (PIV) with CFD–DEM simulations to benchmark hydrodynamic behavior in a 3D spout-fluidized bed. By exploiting easily accessible optical measurements of particle motion at the bed wall and within the spouting region, we obtain quantitative velocity fields that can be directly compared with model predictions of the occluded bed region, without resorting to complex imaging and characterization techniques such as X-ray or magnetic resonance tomography. Experimental benchmarking reveals strong agreement between CFD–DEM and PIV in the spout and annulus regions, while discrepancies near the wall highlight areas for future model development. Here, the proposed integrated experimental–numerical framework will enable a direct connection between measured variables and numerically predicted fluidization performance of dense, surrogate nuclear particle fuel feedstock such that experimental SFB component design can be rapidly evaluated, informing design decisions for nozzle geometry and operating conditions. Future work will extend this framework by correlating quantified fluidization metrics across nozzle geometries and operating conditions with the resulting coating morphology, microstructure, and uniformity. Establishing these correlations will enable predictive links between hydrodynamic performance and coating quality, providing a rational, scalable basis for optimizing SFB design prior to coating deposition.
CFD/DEM↗
Violent mergers can explain the inflated state of some of the fastest stars in the Galaxy
A significant number of hypervelocity stars with velocities between 1500 − 2500 km s −1 have recently been observed. The only plausible explanation so far is that they were produced through thermonuclear supernovae in white dwarf binaries. Since these stars are thought to be surviving donors of Type Ia supernovae, a surprising finding was that these stars are inflated, with radii an order of magnitude higher than expected for Roche-lobe-filling donors. Recent attempts at explaining them have combined 3D hydrodynamical supernova explosion simulations with 1D stellar modelling to explain the impact of supernova shocks on runaway white dwarfs. However, only the hottest and most compact of those runaway stars can so far marginally be reproduced by detailed models of runaways from supernova explosions. In this and a companion paper, we introduce a new AREPO simulation of two massive CO white dwarfs that explode via a violent merger. During the merger, the primary white dwarf ignites when the secondary is on its last orbit and plunging towards the primary. In the corresponding aftermath, the core of the secondary white dwarf of 0.16 M ⊙ remains bound, moving at a velocity of ∼2800 km s −1 . We mapped this object into MESA and show that this runaway star can explain the observations of two hypervelocity stars that were dubbed D6-1 and D6-3 based on their original discovery motivated by the D6 scenario, though the violent merger scenario presented here is somewhat distinct from the D6 scenario.
Astronomy and AstroPhysics↗
athenapk
The LANL fork of the AthenaPK project, which contains curvlinear coordinates The Code is a set of modifications to AthenaPK, an existing open source 3D hydrodynamics and magnetohydrodynamics simulation software for astrophysics, which enable simulation with uniform curvilinear geometries. Specifically it extends the existing Cartesian coordinates to uniform spherical coordinate and uniform cylindrical coordinate systems
Miller, Jonah [@LANL]↗
Theory and modeling of atmospheric turbulence, part 2
Two dimensional geostrophic turbulence driven by a random force is investigated. Based on the Liouville equation, which simulates the primitive hydrodynamical equations, a group-kinetic theory of turbulence is developed and the kinetic equation of the scaled singlet distribution is derived. The kinetic equation is transformed into an equation of spectral balance in the equilibrium and non-equilibrium states. Comparison is made between the propagators and the Green's functions in the case of the non-asymptotic quasi-linear equation to prove the equivalence of both kinds of approximations used to describe perturbed trajectories of plasma turbulence. The microdynamical state of fluid turbulence is described by a hydrodynamical system and transformed into a master equation analogous to the Vlasov equation for plasma turbulence. The spectral balance for the velocity fluctuations of individual components shows that the scaled pressure strain correlation and the cascade transfer are two transport functions that play the most important roles.
Chen, C. M.↗
Liquid oxygen sloshing in Space Shuttle External Tank
This paper describes a numerical simulation of the hydrodynamics within the liquid oxygen tank of the Space Shuttle External Tank during liftoff. Before liftoff, the tank is filled with liquid oxygen (LOX) to approximately 97 percent with the other 3 percent containing gaseous oxygen (GOX) and helium. During liftoff, LOX is drained from the bottom of the tank, and GOX is pumped into the tank's ullage volume. There is a delay of several seconds before the GOX reaches the tank which causes the ullage pressure to decrease for several seconds after liftoff; this pressure 'slump' is a common phenomenon in rocket propulsion. When four slosh baffles were removed from the tank, the ullage gas pressure dropped more rapidly than in all previous flights. The purpose of this analysis was to determine whether the removal of the baffles could have caused the increased pressure 'slump' by changing the LOX surface dynamics. The results show that the LOX surface undergoes very high vertical accelerations (up to 5 g) and, therefore, splashing almost certainly occurs. The number of baffles does not affect the surface if the structural motion is assumed; but, the number of baffles may affect the structural motion of the tank.
Kannapel, M. D.↗
A prediction of the gamma-ray flux from Nova Herculis 1991
UV observations made with the IUE satellite of Nova Her 1991 are presented and discussed. A distance of 3.4 kpc and reddening of 0.6 E(B-V) to Nova Her 1991 are determined by comparing its UV and optical outburst behavior with those of the extragalactic ONeMg nova LMC 1990 No. 1. Her 1991 was also detected by Rosat five days into its outburst, and it is argued that a small fraction of the X-rays may be the product of Comptonized 1.2 and 0.5 MeV emission from the decay of Na-22. If the mass of ejected shell is indeed as high as was estimated by Woodward et al. (1992), namely 0.0001 solar mass, then the present results predicted that Na-22 gamma-ray emission from the ejecta of Nova Her 1991 would still be detectable by GRO in the spring of 1992. The detection level was expected to be about 5 sigma in a two-week observation. Such a detection would strongly guide and constrain the hydrodynamic and nucleosynthesis simulations of nova outbursts on massive white dwarfs.
Starrfield, Sumner↗
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.↗
Testing the pressure-confined Ly alpha cloud model
The Ly-alpha absorption line forest, seen in quasar spectra, is generally interpreted as being due to cosmologically distributed 'clouds' of primordial gas. Analyses of the observations reveal that the number distribution can be described by power laws: dN/dz is proportional to (1 + z)(sup gamma) and dN/dN(sub HI) is proportional to N(sub HI)(sup -beyda), where N(sub HI) is the HI column density. The typical values for power law indices range between 2 is approximately less than gamma is approximately less than 2.6 and 1.7 is approximately less than gamma is approximately less than 1.9. One model postulates that the Ly-alpha clouds are optically thin entities, photoionized by the background UV flux, J(sub nu) is proportional to (1 + z)(sup j), and confined by an adiabatically evolving intercloud medium (ICM): P(z) is proportional to (1 + z)(sup 5). Analytic studies of this model suggest that the ensuing Ly-alpha line statistics can account for the observations (in particular, the dN/dz and the dN/dN(sub HI) distributions) if the cloud mass spectrum is a power law dN/dN is proportional to M(sup -delta), delta is approximately 1.9, and j is approximately 4. One of the simplifying assumptions incorporated into these studies is the existence of a large mass range for the clouds at all epochs, the validity of which is questionable. The pressure-confined model is investigated using a 1-D spherically symmetric hydrodynamical code to simulate cloud evolution over the epoch 1.8 less than z less than 6. This enables us to relax many of the assumptions incorporated in the analytic studies.
Babul, Arif↗
LDEF data: Comparisons with existing models
The relationship between the observed cratering impact damage on the Long Duration Exposure Facility (LDEF) versus the existing models for both the natural environment of micrometeoroids and the man-made debris was investigated. Experimental data was provided by several LDEF Principal Investigators, Meteoroid and Debris Special Investigation Group (M&D SIG) members, and by the Kennedy Space Center Analysis Team (KSC A-Team) members. These data were collected from various aluminum materials around the LDEF satellite. A PC (personal computer) computer program, SPENV, was written which incorporates the existing models of the Low Earth Orbit (LEO) environment. This program calculates the expected number of impacts per unit area as functions of altitude, orbital inclination, time in orbit, and direction of the spacecraft surface relative to the velocity vector, for both micrometeoroids and man-made debris. Since both particle models are couched in terms of impact fluxes versus impactor particle size, and much of the LDEF data is in the form of crater production rates, scaling laws have been used to relate the two. Also many hydrodynamic impact computer simulations were conducted, using CTH, of various impact events, that identified certain modes of response, including simple metallic target cratering, perforations and delamination effects of coatings.
Coombs, Cassandra R.↗
Marshall Space Flight Center's Impact Testing Facility Capabilities
Marshall Space Flight Center s (MSFC) Impact Testing Facility (ITF) serves as an important installation for space and missile related materials science research. The ITF was established and began its research in spacecraft debris shielding in the early 1960s, then played a major role in the International Space Station debris shield development. As NASA became more interested in launch debris and in-flight impact concerns, the ITF grew to include research in a variety of impact genres. Collaborative partnerships with the DoD led to a wider range of impact capabilities being relocated to MSFC as a result of the closure of Particle Impact Facilities in Santa Barbara, California. The Particle Impact Facility had a 30 year history in providing evaluations of aerospace materials and components during flights through rain, ice, and solid particle environments at subsonic through hypersonic velocities. The facility s unique capabilities were deemed a "National Asset" by the DoD. The ITF now has capabilities including environmental, ballistic, and hypervelocity impact testing utilizing an array of air, powder, and two-stage light gas guns to accommodate a variety of projectile and target types and sizes. Numerous upgrades including new instrumentation, triggering circuitry, high speed photography, and optimized sabot designs have been implemented. Other recent research has included rain drop demise characterization tests to obtain data for inclusion in on-going model development. The current and proposed ITF capabilities range from rain to micrometeoroids allowing the widest test parameter range possible for materials investigations in support of space, atmospheric, and ground environments. These test capabilities including hydrometeor, single/multi-particle, ballistic gas guns, exploding wire gun, and light gas guns combined with Smooth Particle Hydrodynamics Code (SPHC) simulations represent the widest range of impact test capabilities in the country.
Finchum, Andy↗