Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrodynamics methods”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Simulation of coupled multiphase flow and geomechanics in porous media with embedded discrete fractures

In fractured natural formations, the equations governing fluid flow and geomechanics are strongly coupled. Hydrodynamical properties depend on the mechanical configuration, and they are therefore difficult to accurately resolve using uncoupled methods. In recent years, significant research has focused on discretization strategies for these coupled systems, particularly in the presence of complicated fracture network geometries. In this work, we explore a finite-volume discretization for the multiphase flow equations coupled with a finite-element scheme for the mechanical equations. Fractures are treated as lower dimensional surfaces embedded in a background grid. Interactions are captured using the embedded discrete fracture model (EDFM) and the embedded finite element method (EFEM) for the flow and the mechanics, respectively. This nonconforming approach significantly alleviates meshing challenges. EDFM considers fractures as lower dimension finite volumes that exchange fluxes with the rock matrix cells. The EFEM method provides, instead, a local enrichment of the finite-element space inside each matrix cell cut by a fracture element. Both the use of piecewise constant and piecewise linear enrichments are investigated. They are also compared to an extended finite element approach. One key advantage of EFEM is the element-based nature of the enrichment, which reduces the geometric complexity of the implementation and leads to linear systems with advantageous properties. Synthetic numerical tests are presented to study the convergence and accuracy of the proposed method. It is also applied to a realistic scenario, involving a heterogeneous reservoir with a complex fracture distribution, to demonstrate its relevance for field applications.

58 GEOSCIENCES↗
DOE OSTI2025

River Dissolved Oxygen Prediction Using Machine Learning Models and Wireless Sensor Measurements

Simultaneous flooding&heat and droughts&heat events can potentially destabilize hydro-meteorological conditions to deteriorate the water quality of Neches River. Machine learning (ML) models utilizing wireless sensor measurements have been applied to predict water quality and optimize various water management strategies. This study aims to develop ML models to predict dissolved oxygen (DO) prediction under various hydro-meteorological conditions and enhance water management decision-making. Wireless sensor measurements of DO, water temperature, sample depth, conductivity, turbidity, and pH, along with discharge from the United States Geological Survey stations, are collected for model inputs at the Pine Island Bayou C749 station (PIB-C749) and Neches River Saltwater Barrier (SWB). Multilayer perceptron neural networks, recurrent neural networks, long short-term memory (LSTM), and bidirectional LSTM (BiLSTM) with and without attention mechanism (AT) are tested to determine the best model, which is applied the rolling forecast method to predict 14-day DO. Traditional and recurrent transfer learning (TL and RTL) methods are adopted to overcome insufficient data at the SWB. The input feature importance analysis using the integrated gradients (IG) algorithm is applied to determine dominant inputs. The results show LSTM-based models are capable handling long sequential data. AT-BiLSTM and RTL-LSTM demonstrate the best performance at the PIB-C749 (RMSE=0.054) and the SWB (RMSE=0.028), respectively. TL and RTL methods significantly improve model performance at the SWB. DO, temperature, and pH show higher importance, consistent with hydrodynamics and water chemistry. Both best models are applied to predict 14-day DO and demonstrate reasonable performance for decision-making. Hydro-meteorological conditions of 2017 flood and 2012 drought events are simulated and reveal that possible hypoxia occurs after flooding due to increasing temperature and turbidity, and DO concentration decreases significantly under heat and drought conditions. In conclusion, LSTM-based models utilizing wireless sensor data can be a timely and effective approach to make appropriate decisions on water resource management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Disorder, interactions, and their interplay in novel narrow-gap Dirac materials and Weyl semimetals

Progress of the modern day condensed matter physics is to a large extent driven by the synthesis of new materials, advances in their experimental characterization and theoretical description. Recent discoveries of novel gapless Weyl semimetals, such as NaBi,CdAs, and BiTe-based films, in which magnetic dopants essentially suppress the gap, have added to the family of graphene and topological insulators actively investigated over the past decade. With the field of novel semimetals rapidly maturing, its focus necessarily shifts from demonstrations of the feasibility of such materials to their quantitative characterization. While the transport and optical properties of graphene and topological insulators are well captured within the picture of free non-interacting electrons, gapless 3D Weyl semimetals and narrow-gap 2D semiconductors with Dirac spectrum are known to be extremely susceptible to disorder and electron-electron interactions. This susceptibility obscures the manifestations of nontrivial band structure -- like quantum anomalous Hall effect -- of the new topological materials. Among particular projects to be addressed are: 1) optical conductivity of 3D gapless Dirac fermions in the presence of smooth disorder, 2) interplay of disorder and Coulomb interactions in the spectral properties of such fermions, 3) formation and structure of the impurity band with Coulomb supercritical clusters, 4) Coulomb interaction-driven renormalization of the electron spectrum and of the transport response in the presence of strong magnetic field, 5) instanton approach to the disorder-induced fluctuation states in zero-gap 3D materials, and 6) the role of disorder in quantum anomalous Hall effect. The proposal relies upon the investigators' previous broad expertise in interacting and disordered electron systems. The methods to be employed include perturbative diagrammatic technique, non-perturbative instanton and self-consistent approximations, hydrodynamics of electron liquid. Both analytical as well as numerical approaches are to be employed. The anticipated broader outcome of the proposal includes gaining an in-depth understanding of the interplay of the disorder and interactions under the conditions when this interplay has the most dramatic impact on observables. Traditionally, interaction effects are among the most challenging and interesting problems of condensed matter physics. Similarly, disordered systems typically present very difficult but extremely rich problems in the description of various materials. Importantly, understanding the spectral and transport properties of such materials not only presents the fundamental objective, but is also of particular interest for many applications, such as computation, memory, optics, plasmonics. In particular realization of the quantum anomalous Hall effect may lead to the development of low-power-consumption electronics. Indeed, a major constraint for practical use of the quantum Hall effect is limited by the requirement of the quantizing magnetic field. At the same time, the quantum anomalous Hall effect samples exhibit non-dissipative edge quantum transport in a zero magnetic field.

36 MATERIALS SCIENCE↗
DOE OSTI2025

3D Convective Urca Process in a Simmering White Dwarf

Abstract A proposed setting for thermonuclear (Type Ia) supernovae is a white dwarf that has gained mass from a companion to the point of carbon ignition in the core. In the early stages of carbon burning, called the simmering phase, energy released by the reactions in the core drive the formation and growth of a core convection zone. One aspect of this phase is the convective Urca process, a linking of weak nuclear reactions to convection, which may alter the composition and structure of the white dwarf. The convective Urca process is not well understood and requires 3D fluid simulations to properly model the turbulent convection, an inherently 3D process. Because the neutron excess of the fluid both sets and is set by the extent of the convection zone, the realistic steady state can only be determined in simulations with real 3D mixing processes. Additionally, the convection is relatively slow (Mach number less than 0.005) and thus a low Mach number method is needed to model the flow over many convective turnovers. Using the MAESTROeX low Mach number hydrodynamic software, we present the first full-star 3D simulations of the A = 23 convective Urca process, spanning hundreds of convective turnover times. Our findings on the extent of mixing across the Urca shell, the characteristic velocities of the flow, the energy-loss rates due to neutrino emission, and the structure of the convective boundary can be used to inform 1D stellar models that track the longer-timescale evolution.

Boyd, Brendan (ORCID:0000000254199751)↗
NASA NTRS1976

Interplanetary disturbances in the solar wind produced by density, temperature, or velocity pulses at 0.08 AU

Time-dependent solutions of a one-fluid model of the interplanetary medium are investigated. This set of unsteady hydrodynamic equations has been written in conserved form in order to apply the Lax-Wendroff (1960) method for the solution of this problem. The initial disturbance is specified by a pulse at 0.08 AU. Physically, this pulse can be interpreted as having been caused by a solar flare, surge, or any other solar disturbance. The equilibrium condition is determined to be the steady solution of the governing equations and represents the quiet solar wind. Results are presented in terms of density, temperature, and velocity profiles of the interplanetary gas flow at heliocentric distances up to 6 AU at several times. Also, the trajectories of disturbances for various initial pulses are shown. Some June 1972 interplanetary observational data are compared with these theoretical calculations. On the basis of these results, the effects of solar disturbances on the interplanetary environment (such as the generation of large nonlinear wave trains in the shock wakes) can be inferred.

Wu, S. T.↗
NASA NTRS1985

Formation of the prelunar accretion disk

Considerations related to the angular momentum of the earth-moon system have led Ward and Cameron (1978) to the suggestion of a collisional origin of the moon. The required projectile would be about the mass of Mars (i.e., about 0.1 earth masses). The possibility has been considered that both the proto-earth and the projectile would be remnants of the formation of giant gaseous protoplanets. Processes of lunar formation studied take into account the viscous dissipation of an accretion disk. The present paper is concerned with the formation of the prelunar accretion disk. The hydrodynamics of a compressible viscous gas are simulated with the aid of the particle-in-cell method. The obtained results imply that the projectile is in an elliptical earth-crossing orbit.

Cameron, A. G. W.↗
NASA NTRS1990

The vibrating ribbon problem revisited

A revised formal solution of the vibrating ribbon problem of hydrodynamic stability is presented. The initial formulation of Gaster (1965) is modified by application of the Briggs method and a careful treatment of the complex double Fourier transform inversions. Expressions are obtained in a natural way for the discrete spectrum as well as for the four branches of the continuous spectra. These correspond to discrete and branch-cut singularities in the complex wavenumber plane. The solutions from the continuous spectra decay both upsteam and downstream of the ribbon, with the decay in the upstream direction being much more rapid than that in the downstream direction. Comments and clarification of related prior work are made.

Ashpis, David E.↗
DOE OSTI2024

WEC-Sim Modeling of Laminar Scientific's Patented Seesaw Wave Energy Converter: Preprint

Laminar Scientific's patented seesaw wave energy converter was modeled in WEC-Sim to predict performance. The device operates by utilizing ocean surface waves to rotate a truss in pitch about a pivot. The pivot is located at the top of two pylons, which are embedded in the seafloor. The seesaw has a float on either end, and the buoyancy forces from each float cause the system to rise or fall with passing waves. Device performance relies upon seesaw length and ocean wavelength creating an antiphase effect. The seesaw truss has an adjustable length intended to achieve this effect. The operation method enforces a narrow band of wavelengths which induce the largest rotational motion from the device. The hydrodynamic analysis of the device was performed using Capytaine, and the results confirmed that the device operates best in a narrow frequency band. Four float-to-float spacing cases and three pylon radii were examined. The hydrodynamic results indicate a match between the model and the physical expectations for the device, and that varying the pylon radii by 0.1-m increments for three instances creates minimal changes in hydrodynamic properties. Power matrices for three float spacing cases of the device were simulated with Joint North Sea Wave Project spectra waves and optimal power take-off damping in WEC-Sim. The maximum average power production for the 15-m spacing case was 14.1 kW with a 5.0-s peak wave period and 4-m significant wave heights. Plots of capture-width ratios indicated that the device performance was linear and confirmed that the device is optimal in a narrow frequency band. The maximum percentage of the available wave power produced by the 15-m device was approximately 16%. Simulations of the device in regular waves were used to produce plots of average power compared to a ratio of float spacing to wavelength. These plots indicate that the power production is maximized at a ratio of 0.5, and further confirm that the device has a narrow frequency response. The device was simulated at an example field location, where the device produced an annual average power rating of 1.6 kW given an average omnidirectional wave climate of 10.3 kW m-1 and an optimal, linearized power take-off model. While the maximum predicted device performance is reliant upon a narrow band of wave frequencies, the conducted analysis provides an opportunity to improve device design prior to prototyping and testing. Modifying the design to respond to a broader frequency range would improve device performance.

numerical modeling↗
DOE OSTI2021

Sondheimer oscillations as a probe of non-ohmic flow in WP 2 crystals

As conductors in electronic applications shrink, microscopic conduction processes lead to strong deviations from Ohm’s law. Depending on the length scales of momentum conserving (l MC ) and relaxing (l MR ) electron scattering, and the device size (d), current flows may shift from ohmic to ballistic to hydrodynamic regimes. So far, an in situ methodology to obtain these parameters within a micro/nanodevice is critically lacking. In this context, we exploit Sondheimer oscillations, semi-classical magnetoresistance oscillations due to helical electronic motion, as a method to obtain l MR even when l MR $\gg$ d. We extract l MR from the Sondheimer amplitude in WP 2 , at temperatures up to T ~ 40 K, a range most relevant for hydrodynamic transport phenomena. Our data on μm-sized devices are in excellent agreement with experimental reports of the bulk l MR and confirm that WP 2 can be microfabricated without degradation. These results conclusively establish Sondheimer oscillations as a quantitative probe of l MR in micro-devices.

36 MATERIALS SCIENCE↗
NASA NTRS2022

Finite Difference Methods for Turbulence Simulations

The optimal finite difference discretization used in simulations of turbulent flows is influenced by both, the type of the scale resolving simulation (DNS or LES), as well as the flow-physics (hydrodynamic instabilities, shocks, acoustics, etc.) one expects to resolve. Insight into dispersion and dissipation error requirements for some common scale-resolving simulation scenarios help to highlight the issues faced in selecting a scheme.

Finite Difference Methods↗
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↗
DOE OSTI2025

Gas and dust dynamics in γ Cephei-type disks

Context. Giant planets are observed orbiting the primary stars of close binary systems. Such planets may have formed in compact circumprimary disks, which once surrounded these stars, under conditions much different than those encountered around single stars. Aims. In order to quantify the effects of the strong gravitational perturbations exerted on circumprimary disk material, the three-dimensional (3D) dynamics of gas and dust in orbit around the primary star of a compact and eccentric binary system was modeled by applying the stellar and orbital parameters of γ Cephei, a well-known system that can be representative of a class of close binaries. Methods. Circumprimary gas was approximated as an Eulerian viscous and compressible fluid and modeled by means of 3D hydrodynamical simulations, assuming locally isothermal conditions in the medium around the primary star. Dust grains were modeled as Lagrangean particles, subjected to gravity and aerodynamic drag forces. Models that include a giant planet were also considered. Results. Models indicate that spiral density waves excited around pericenter passage propagate toward the inner boundary of the disk, through at least a few pressure scale-heights from the mid-plane, inducing radial and vertical mixing in the gas. However, perturbations imparted to gas, both in terms of eccentricity and precession, are far weaker than previously estimated by two-dimensional (2D) simulations. Models predict small eccentricities, ≲0.03, and slow retrograde precession. The addition of a giant planet does not change the low eccentricity state of the disk. The parameters applied to the disk would lead to the formation of a massive planet, many times the mass of Jupiter, in agreement with some observations. Micron to mm-size dust grains are well coupled to the gas, resulting in similar dynamics and statistically similar distributions of orbital elements. The planet only affects the dust distributions locally. In agreement with outcomes of recent 2D models, the lifetime of an isolated circumprimary disk would be brief, ~10 5 years, because of its compact nature, requiring a long-term external supply of mass to allow for the in situ formation of a giant planet.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Cross-beam energy transfer saturation by ion trapping-induced detuning

The performance of direct-drive inertial confinement fusion implosions relies critically on the coupling of laser energy to the target plasma. Cross beam energy transfer (CBET), the resonant exchange of energy between intersecting laser beams mediated by ponderomotively driven ion-acoustic waves (IAW), inhibits this coupling by scattering light into unwanted directions. Furthermore, the variety of beam intersection angles and varying plasma conditions in an implosion results in IAWs with a range of phase velocities. Here we show that CBET saturates through a resonance detuning that depends on the IAW phase velocity and that results from trapping-induced modifications to the ion distribution functions. For smaller phase velocities, the modifications to the distribution functions can rapidly thermalize in the presence of mid-Z ions, leading to a blueshift in the resonant frequency. For larger phase velocities, the modifications can persist, leading to a redshift in the resonant frequency. Ultimately, these results may reveal pathways towards CBET mitigation and inform reduced models for radiation hydrodynamics codes to improve their predictive capability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Enhanced electron acceleration by high-intensity lasers in extended (confined) preplasma in cone targets

Here we report on experimental results from a high-intensity laser interaction with cone targets that increase the number (×3) and temperature (×3) of the measured hot electrons over a traditional planar target. This increase is caused by a substantial increase in the plasma density within the cone target geometry, which was induced by 17 ± 9 mJ prepulse that arrived 1.5 ns prior to the main high intensity (>10 19 W/cm 2 ). Three-dimensional hydrodynamic simulations are conducted using hydra which show that the cone targets create substantially longer and denser plasma than planar targets due to the geometric confinement of the expanding plasma. The density within the cone is a several hundred-micron plasma “shelf” with a density of approximately 10 20 n e /cc. The HYDRA simulated plasma densities are used as the initial conditions for two-dimensional particle-in-cell simulations using EPOCH. These simulations show that the main acceleration mechanism is direct-laser-acceleration, with close agreement between experimentally measured and simulated electron temperatures. Further analysis is conducted to investigate the acceleration of the electrons within the long plasma generated within a compound parabolic concentrator by the prepulse.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Laser wavelength dependence of particle acceleration mechanisms in high intensity laser–solid density plasma interactions

We investigate the generation of relativistic electrons and the subsequent ion acceleration due to target-normal sheath acceleration when ultra-intense (⁠ I > 10 18 W/cm 2 ⁠) short pulse (⁠ τ L < 10ps⁠) lasers are incident onto solid density targets as laser wavelength is varied. Scaling laws for the hot electron temperature, T hot ⁠, and the maximum ion energy, E max ⁠, are recast as a function of laser wavelength. These predictions are compared to results from particle-in-cell computer simulations in a variety of geometries, including cases where realistic plasma density profiles as determined by a radiation hydrodynamics code are used. It is found that the wavelength dependence observed in simulation is less pronounced than what is predicted from the well-established scaling laws. An assessment of how switching to longer laser wavelengths, specifically 2 μm Tm:YLF technology, would impact current high energy density science applications and diagnostics is made.

Electromagnetism↗
NASA NTRS1998

Excitation of Continuous and Discrete Modes in Incompressible Boundary Layers

This report documents the full details of the condensed journal article by Ashpis & Reshotko (JFM, 1990) entitled "The Vibrating Ribbon Problem Revisited." A revised formal solution of the vibrating ribbon problem of hydrodynamic stability is presented. The initial formulation of Gaster (JFM, 1965) is modified by application of the Briggs method and a careful treatment of the complex double Fourier transform inversions. Expressions are obtained in a natural way for the discrete spectrum as well as for the four branches of the continuous spectra. These correspond to discrete and branch-cut singularities in the complex wave-number plane. The solutions from the continuous spectra decay both upstream and downstream of the ribbon, with the decay in the upstream direction being much more rapid than that in the downstream direction. Comments and clarification of related prior work are made.

Ashpis, David E.↗
NASA NTRS2003

Minimum impulse thruster valve design and development

The design and development of a minimum impulse thruster valve was conducted, by Moog, under contract by NASA's Jet Propulsion Laboratory, California Institute of Technology, for deep space propulsion systems. The effort was focused on applying known solenoid design techniques scaled to provide a 1 -millisecond response capability for monopropellant, hydrazine ACS thruster applications. The valve has an extended operating temperature range of 20(deg)F to +350(deg)F with a total mass of less than 25 grams and nominal power draw of 7 watts. The design solution resulted in providing a solenoid valve that is one-tenth the scale of the standard product line. The valve has the capability of providing a mass flow rate of 0.0009 pounds per second hydrazine. The design life of 1,000,000 cycles was demonstrated both dry and wet. Not all design factors scaled as expected and proved to be the focus of the final development effort. These included the surface interactions, hydrodynamics and driver electronics. The resulting solution applied matured design approaches to minimize the program risk with innovative methods to address the impacts of scale.

impulse thruster valve↗