Search NASA⌕ Search

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.

At least 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Hourglass control in staggered-grid hydrodynamics using virtual element stabilization techniques

Numerical simulations using the staggered-grid hydrodynamics (SGH) discretization suffer from hourglass instabilities. In this work, we develop a stabilization method to suppress the hourglass instabilities using techniques from the virtual element method (VEM). The stiffness matrix of the VEM consists of two terms: the consistency matrix which is rank deficient and the stability matrix. Here, we first show that in two dimensions and on general polygons, the stiffness matrix of the SGH is identical to the consistency matrix of the linear VEM for both the diffusion equation and the linear elasticity equation. These analyses explain the origin of the hourglass instabilities of the SGH discretization method, and establish a theoretical foundation for our proposed stabilization method by augmenting the stiffness matrix of the SGH discretization using the VEM stability matrix. Then, we present numerical examples using Lagrangian SGH simulations. The numerical experiments demonstrate that the proposed VEM stabilization method is effective at eliminating hourglass modes in the SGH discretization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A lumped particle direct simulation Monte-Carlo method combined with the collisional-radiative model for simulations of non-equilibrium laser-induced plasma plumes

Collisional plasma plumes induced by laser irradiation of material targets exhibit large variations in local density as well as ionization and excitation states, making purely hydrodynamic or kinetic simulations inaccurate or infeasible. To address this challenge and capture non-equilibrium effects in laser-induced plasma plumes at arbitrary degrees of ionization, we develop a hybrid computational approach that combines the kinetic direct simulation Monte Carlo (DSMC) method with a collisional-radiative model (CRM). This ℓDSMC-CRM approach utilizes a lumped particle method to represent minor fractions of excited ions in particle-based simulations and a special coarse-graining technique for atomic spectra and photoionization rates, ensuring numerical convergence at reduced computational cost. The hybrid approach is applied to simulate spatially homogeneous relaxation as well as one- and two-dimensional expansions of plasma plumes induced by irradiation of a copper target by a nanosecond laser pulse in a vacuum or background gas. The comparison with an equilibrium model, where local Saha-Boltzmann equilibrium is enforced, shows that the non-equilibrium effects play a dominant role. The equilibrium model can fail to predict the flow structure and strongly underestimate the degree of absorption of laser radiation by the plume. The ℓDSMC-CRM approach is validated against experimental data demonstrating reasonable agreement with the experimental electron density and temperature, while the equilibrium model is found to dramatically underestimate electron density and temperature. The flexibility of the ℓDSMC-CRM approach allows for its seamless integration into existing DSMC frameworks, making it a valuable tool for high-fidelity plasma modeling in laser-material interactions, laser-based manufacturing, and beyond.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1995

A modified thermal conductivity for low density plasma magnetic flux tubes

In response to inconsistencies which have arisen in results from a hydrodynamic model in simulation of high ion temperature (1-2 eV) observed in low density, outer plasmasphere flux tubes, we postulate a reduced thermal conductivity coefficient in which only particles in the loss cone of the quasi-collisionless plasma contribute to the thermal conduction. Other particles are assumed to magnetically mirror before they reach the topside ionosphere and therefore not to remove thermal energy from the plasmasphere. This concept is used to formulate a mathematically simple, but physically limiting model for a modified thermal conductivity coefficient. When this modified coefficient is employed in the hydrodynamic model in a case study, the inconsistencies between simulation results and observations are largely resolved. The high simulated ion temperatures are achieved with significantly lower ion temperatures in the topside ionosphere. We suggest that this mechanism may be operative under the limited low density, refilling conditions in which high ion temperatures are observed.

Comfort, R. H.↗
NASA NTRS1994

Powder agglomeration in a microgravity environment

This is the final report for NASA Grant NAG3-755 entitled 'Powder Agglomeration in a Microgravity Environment.' The research program included both two types of numerical models and two types of experiments. The numerical modeling included the use of Monte Carlo type simulations of agglomerate growth including hydrodynamic screening and molecular dynamics type simulations of the rearrangement of particles within an agglomerate under a gravitational field. Experiments included direct observation of the agglomeration of submicron alumina and indirect observation, using small angle light scattering, of the agglomeration of colloidal silica and aluminum monohydroxide. In the former class of experiments, the powders were constrained to move on a two-dimensional surface oriented to minimize the effect of gravity. In the latter, some experiments involved mixture of suspensions containing particles of opposite charge which resulted in agglomeration on a very short time scale relative to settling under gravity.

Cawley, James D.↗
DOE OSTI2026

Signatures of QCD conductivities in heavy-ion collisions

Dissipative processes are pivotal for understanding the hydrodynamic evolution of hot and dense quantum chromodynamics (QCD) matter created in relativistic nuclear collisions. The interplay of multiple conserved charges—net baryon, strangeness, and electric charge—is of particular interest. Here, we simulate the longitudinal hydrodynamic evolution with the three diffusion currents in a hydrodynamic model with a lattice-QCD-based equation of state, NEOS -4 D , and estimate rapidity distributions including diffusive corrections to the phase-space distribution in the presence of multiple charges, which ensure charge conservation at particlization. We determine the response of particle yields at midrapidity to changes in the diagonal and off-diagonal conductivities. Inversely, we find that most components of the conductivity matrix can be constrained experimentally using identified particle multiplicities at different collision energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2010

Vortex Shedding Inside a Baffled Air Duct

Common in the operation of both segmented and un-segmented large solid rocket motors is the occurrence of vortex shedding within the motor chamber. A portion of the energy within a shed vortex is converted to acoustic energy, potentially driving the longitudinal acoustic modes of the motor in a quasi-discrete fashion. This vortex shedding-acoustic mode excitation event occurs for every Reusable Solid Rocket Motor (RSRM) operation, giving rise to subsequent axial thrust oscillations. In order to better understand this vortex shedding/acoustic mode excitation phenomena, unsteady CFD simulations were run for both a test geometry and the full scale RSRM geometry. This paper covers the results from the subscale geometry runs, which were based on work focusing on the RSRM hydrodynamics. Unsteady CFD simulation parameters, including boundary conditions and post-processing returns, are reviewed. The results were further post-processed to identify active acoustic modes and vortex shedding characteristics. Probable locations for acoustic energy generation, and subsequent acoustic mode excitation, are discussed.

Davis, Philip↗
DOE OSTI2025

Multiphase computational fluid dynamics modeling of reacting flows in absorption columns for carbon capture

First-principles derived computational fluid dynamics (CFD) simulations have been proposed as a fundamental tool for investigating solvent-based CO 2 absorption in packed columns due to their ability to accurately represent the underlying nonlinear, multiscale dynamics. Numerous studies have previously utilized such CFD simulations to investigate hydrodynamics of columns with structured and random packings by assessing the key hydrodynamic metrics such as the interfacial and wetted areas. While mapping such metrics for different conditions is essential to the optimization of absorption columns, it is not sufficient, as the CO 2 capture rate depends also on the coupled, nonlinear dynamics from the underlying chemical reaction kinetics, thermodynamics, and heat-transfer rates. In this work, we present detailed CFD simulation results augmented by incorporating the effects of interfacial physical mass transfer of CO 2 , heat release from chemical reaction kinetics, and thermophysical property variations from resulting temperature gradients. We demonstrate the applicability of the proposed approach in numerically assessing the performance of packed columns by evaluating key hydrodynamic quantities, CO 2 absorption rates, and temperature rise in a reference column with packings that are structurally similar to the Sulzer Mellapak™ 250.Y packing, for different solvent inflow velocities and temperatures. Predictions from simulation results are found to be consistent with the trends in experimental observations from the literature, suggesting that the predictive capabilities of the simulation framework can be leveraged to guide the future development of absorber-column designs and optimized process flowsheets.

Absorption columns↗
NASA NTRSDate not supplied

Nanoflare Heating of an X-Ray Bright Point

Nanoflares are thought to be one of the prime candidates that can keep the solar corona to its multimillion kelvin temperature. Individual nanoflares are difficult to detect with the present generation instruments, however their presence can be inferred by comparing the nanoflare heated simulated plasma emissions with the observed emission. Here, we present a simulation of emission from an X-ray Bright Point (XBP) that was observed by the Marshall Grazing Incidence X-ray Spectrometer (MaGIXS), along with concurrent observations from SDO/AIA and Hinode/XRT. We use EBTEL hydrodynamic code to simulate the XBP loops. Length and magnetic field strength of these loops are derived from the potential field extrapolation of the observed photospheric magnetogram by HMI/SDO. Each loop is assumed to be heated by random nanoflares, whose magnitude and frequency are determined by the looplength and magnetic field strength. The simulated outputs are used to predict the intensity of spectrally pure map of Fe-18, Fe-17, Ne-9 ,O-8, O-9, Ne-9 etc, which are then compared with the derived intensity from MaGIXS observation. Further we have predicted the intensity map as observed by AIA and XRT and compared them with the observation. We also estimated the temperature distribution of the XBP from the simulation and found a good agreement with the derived distribution from MaGIXS observation.

coronal heating↗
NASA NTRS2006

Unsteady Flow Simulation of High-speed Turbopumps

Computation of high-speed hydrodynamics requires high-fidelity simulation to resolve flow features involving transient flow, cavitation, tip vortex and multiple scales of unsteady fluctuations. One example of this type in aerospace is related to liquid-fueled rocket turbopump. Rocket turbopumps operate under severe conditions at very high rotational speeds typically at thousands of rpm. For example, the Shuttle orbiter low-pressure-fuel-turbopump creates transient flow features associated with reverse flows, tip clearance effects, secondary flows, vortex shedding, junction flows, and cavitation effects. Flow unsteadiness originating from the orbiter Low-Pressure-Fuel-Turbopump (LPFTP) inducer is one of the major contributors to the high frequency cyclic loading that results in high cycle fatigue damage to the flow liners just upstream of the LPFTP. The reverse flow generated at the tip of the inducer blades travels upstream and interacts with the bellows cavity. Simulation procedure for this type high-speed hydrodynamic problems requires a method for quantifying multi-scale and multi-phase flow as well as an efficient high-end computing strategy. The current paper presents a high-fidelity computational procedure for unsteady hydrodynamic problems using a high-speed liquid-fueled rocket turbopump.

Kiris, Cetin C.↗
NASA NTRS2018

Investigations of Pole Erosion Mechanisms in the 12.5 kW HERMeS Hall Thruster Using Numerical Simulations and Ion Velocity Measurements

When the Hall Effect Rocket with Magnetic Shielding (HERMeS) was tested at the NASA GRC’s VF-5 the erosion rate on the pole covers was found to be at least a factor of two higher at 300 V and 20.8 A than at other operating conditions in the 400 V to 600 V range. Simulations using our hydrodynamics code Hall2De do not predict increased erosion rates at 300 V but accurately compute erosion rates similar to the measurements at the other operating conditions. We investigate the source of the discrepancy between measurements and simulations at 300 V using a combination of numerical simulations and experimental measurements of the ion velocity fields in the acceleration region of the thruster. By examining previous simulations that predicted the measured erosion rates for other thrusters and operating conditions, we determine that sputtering of the pole surfaces by high energy ions is the most likely mechanism behind the erosion rates at 300 V. High energy ions sputter the pole surfaces when the acceleration region of the thruster is downstream of the pole surface plane. When the latter occurs, the curvature of the plasma potential contours at the edges of the channel accelerates a small fraction of the high energy ions radially. We find that neither our Hall2De simulations nor the experimental measurements produce ions of sufficiently high energy and flux to the inner pole cover to explain the measured erosion there. Only at the outer edge of the cover is where we find the ions needed to yield simulation results that are comparable to the erosion measurements. We argue that one possible mechanism for the higher erosion is high energy ions that graze the channel corner of the pole and become trapped in a sheath that develops downstream. Within the sheath, we find that the electric field is large enough to turn the ions towards the surface of the pole. A simplified calculation shows that the erosion rates produced by sputtering of trapped ions are similar to the measurements. We also propose that the presence of local oscillations and high-energy cathode ions may be sources of the enhanced erosion We conclude this article by proposing a combination of analysis, simulation, and experimental measurements that can be used to address the validity of these hypotheses.

Mikellides, Ioannis G.↗
NASA NTRS1991

Cometary activity, discrete outgassing areas, and dust-jet formation

Conceptual models for various types of features observed in cometary comae (jets, spirals, halos, fans, etc.), their computer simulation, and the hydrodynamic models for jet formation are critically reviewed, and evidence for anisotropic, strongly collimated flows of ejecta emanating from discrete active regions (vents) on the rotating cometary nuclei is presented. Techniques employed to generate synthetic comet images that simulate the features observed are described, and their relevance to the primary objects of coma-morphology studies is discussed. Modeling of temporal variations in the water emission from discrete active regions suggests that production curves asymmetric with respect to perihelion should be commonplace. Critical comparisons with the activity profiles of Enke's comet and with light curves of disappearing comets and comets that undergo outbursts are presented. Recent developments in the understanding of the processes that cause the nongravitational perturbations of cometary motions are reviewed, and the observed discontinuities are identified with the birth of new sources and/or deactivation of old vents.

Sekanina, Z.↗
NASA NTRS1992

Simulating hypervelocity impact effects on structures using the smoothed particle hydrodynamics code MAGI

Analysis of interaction occurring between space debris and orbiting structures is of great interest to the planning and survivability of space assets. Computer simulation of the impact events using hydrodynamic codes can provide some understanding of the processes but the problems involved with this fundamental approach are formidable. First, any realistic simulation is necessarily three-dimensional, e.g., the impact and breakup of a satellite. Second, the thickness of important components such as satellite skins or bumper shields are small with respect to the dimension of the structure as a whole, presenting severe zoning problems for codes. Thirdly, the debris cloud produced by the primary impact will yield many secondary impacts which will contribute to the damage and possible breakup of the structure. The problem was approached by choosing a relatively new computational technique that has virtues peculiar to space impacts. The method is called Smoothed Particle Hydrodynamics.

Libersky, Larry↗
DOE OSTI2025

Exploring Flood Predictability in Taiwan through Coupled Atmospheric–Hydrological and High-Performance Hydrodynamic Models

Effective flood simulation capabilities can tremendously support early warning and disaster prevention. To examine the applicability of a fully physics-based and high-performance flood simulation and forecasting modeling framework for a flood-prone region in Taiwan, we conduct a numerical experiment that couples the Weather Research and Forecasting (WRF) Model, WRF-Hydrological modeling system (WRF-Hydro), and the Two-Dimensional Runoff Inundation Toolkit for Operational Needs (TRITON) to perform integrated rainfall, streamflow, and flood simulations. Furthermore, we first use the coupled WRF and WRF-Hydro (WWH) to predict rainfall and streamflow and then drive TRITON with the predicted streamflow hydrographs to simulate flood depth and inundation area. With the refined spatial resolution and parameterization, this framework can better predict rainfall with reasonable spatial patterns. Although WWH could overestimate the amount of rainfall in some areas, the uncertain rainfall–streamflow predictions produce reasonable flood maps able to pinpoint regions at risk of flooding. In terms of model efficiency, the graphics processing unit–based computation can yield a speed-up factor as high as ∼13 compared to the central processing unit–based computation, promoting the efficacy of the coupled modeling framework in practical real-time flood forecasting.

Coupled models↗
NASA NTRS2019

SI IV Resonance Line Emission during Solar Flares: Non-LTE, Non-equilibrium, Radiation Transfer Simulations

The Interface Region Imaging Spectrograph routinely observes the Si IV resonance lines. When analyzing quiescent observations of these lines, it has typically been assumed that they form under optically thin conditions. This is likely valid for the quiescent Sun, but this assumption has also been applied to the more extreme flaring scenario. We used 36 electron-beam-driven radiation hydrodynamic solar flare simulations, computed using the RADYN code, to probe the validity of this assumption. Using these simulated atmospheres, we solved the radiation transfer equations to obtain the non-LTE, nonequilibrium populations, line profiles, and opacities for a model silicon atom, including charge exchange processes. This was achieved using the "minority species" version of RADYN. The inclusion of charge exchange resulted in a substantial fraction of Si IV at cooler temperatures than those predicted by ionization equilibrium. All simulations with an injected energy flux F> 5 x 10(exp 10) erg cm(exp -2) s(exp -1) resulted in optical depth effects on the Si IV emission, with differences in both intensity and line shape compared to the optically thin calculation. Weaker flares (down to F ~ 10(exp 9) erg cm(exp -2) s(exp -1)) also resulted in Si IV emission forming under optically thick conditions, depending on the other beam parameters. When opacity was significant, the atmospheres generally had column masses in excess of 5 x 10-6 g cm-2 over the temperature range 40-100 kK, and the Si IV formation temperatures were between 30 and 60 kK. We urge caution when analyzing Si IV flare observations, or when computing synthetic emission without performing a full radiation transfer calculation.

Kerr, Graham S.↗
NASA NTRS2008

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 1960% 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. Relocated test equipment was dated and in need of upgrade. 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. Future ITF improvements will be focused on continued instrumentation and performance enhancements. These enhancements will allow further, more in-depth, characterization of rain drop demise characterization and evaluation of ice crystal impact. Performance enhancements also include increasing the upper velocity limit of the current environmental guns to allow direct environmental simulation for missile components. 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 grins, 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.

Evans, Steve↗
NASA NTRS2008

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. Relocated test equipment was dated and in need of upgrade. 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. Future ITF improvements will be focused on continued instrumentation and performance enhancements. These enhancements will allow further, more in-depth, characterization of rain drop demise characterization and evaluation of ice crystal impact. Performance enhancements also include increasing the upper velocity limit of the current environmental guns to allow direct environmental simulation for missile components. 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.

Evans, Steve↗
NASA NTRS2001

Jet Stability and the Generation of Superluminal and Stationary Components

We present a numerical simulation of the response of an expanding relativistic jet to the ejection of a superluminal component. The simulation has been performed with a relativistic time-dependent hydrodynamical code from which simulated radio maps are computed by integrating the transfer equations for synchrotron radiation. The interaction of the superluminal component with the underlying jet results in the formation of multiple conical shocks behind the main perturbation. These trailing components can be easily distinguished because they appear to be released from the primary superluminal component instead of being ejected from the core. Their oblique nature should also result in distinct polarization properties. Those appearing closer to the core show small apparent motions and a very slow secular decrease in brightness and could be identified as stationary components. Those appearing farther downstream are weaker and can reach superluminal apparent motions. The existence of these trailing components indicates that not all observed components necessarily represent major perturbations at the jet inlet; rather, multiple emission components can be generated by a single disturbance in the jet. While the superluminal component associated with the primary perturbation exhibits a rather stable pattern speed, trailing components have velocities that increase with distance from the core but move at less than the jet speed. The trailing components exhibit motion and structure consistent with the triggering of pinch modes by the superluminal component. The increase in velocity of the trailing components is an indirect consequence of the acceleration of the expanding fluid, which is assumed to be relativistically hot; if observed, such accelerations would therefore favor an electron-positron (as opposed to proton rest mass) dominated jet.

Agudo, Ivan↗