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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
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 OSTI2018

Solution Self-Assemblies of Sequence-Defined Ionic Peptoid Block Copolymers

A series of amphiphilic ionic peptoid block copolymers where the total number (1 or 3) and position of ionic monomers along the polymer chain are precisely controlled have been synthesized by the submonomer method. Upon dissolution in water at pH = 9, the amphiphilic peptoids self-assemble into small spherical micelles having hydrodynamic radius in ~5–10 nm range and critical micellar concentration (CMC) in the 0.034–0.094 mg/mL range. Small-angle neutron scattering (SANS) analysis of the micellar solutions revealed unprecedented dependence of the micellar structure on the number and position of ionic monomers along the chain. It was found that the micellar aggregation number (N agg ) and the micellar radius (R m ) both increase as the ionic monomer is positioned progressively away from the junction of the hydrophilic and hydrophobic segments along the polymer chain. By defining an ionic monomer position number ( n ) as the number of monomers between the junction and the ionic monomer, N agg exhibited a power law dependence on n with an exponent of ~1/3 and ~3/10 for the respective singly and triply charged series. On the other hand, R m exhibited a weaker dependence on the ionic monomer position by a power law relationship with an exponent of ~1/10 and ~1/20 for the respective singly and triply charged series. Moreover, R m was found to scale with N agg in a power-law relationship with an exponent of 0.32 for the singly charged series, consistent with a weakly charged ionic star-like polymer model in the unscreened regime. This study demonstrated a unique method to precisely tailor the structure of small spherical micelles based on ionic block copolymers through controlling the sequence and position of the ionic monomer.

36 MATERIALS SCIENCE↗
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↗
DOE OSTI2016

Hydrodynamic trapping for rapid assembly and in situ electrical characterization of droplet interface bilayer arrays

The droplet interface bilayer (DIB) is a modular technique for assembling planar lipid membranes between water droplets in oil. The DIB method thus provides a unique capability for developing digital, droplet-based membrane platforms for rapid membrane characterization, drug screening and ion channel recordings. This paper demonstrates a new, low-volume microfluidic system that automates droplet generation, sorting, and sequential trapping in designated locations to enable the rapid assembly of arrays of DIBs. The channel layout of the device is guided by an equivalent circuit model, which predicts that a serial arrangement of hydrodynamic DIB traps enables sequential droplet placement and minimizes the hydrodynamic pressure developed across filled traps to prevent squeeze-through of trapped droplets. Furthermore, the incorporation of thin-film electrodes fabricated via evaporation metal deposition onto the glass substrate beneath the channels allows for the first time in situ, simultaneous electrical interrogation of multiple DIBs within a sealed device. Combining electrical measurements with imaging enables measurements of membrane capacitance and resistance and bilayer area, and our data show that DIBs formed in different trap locations within the device exhibit similar sizes and transport properties. Simultaneous, single channel recordings of ion channel gating in multiple membranes are obtained when alamethicin peptides are incorporated into the captured droplets, qualifying the thin-film electrodes as a means for measuring stimuli-responsive functions of membrane-bound biomolecules. Furthermore, this novel microfluidic-electrophysiology platform provides a reproducible, high throughput method for performing electrical measurements to study transmembrane proteins and biomembranes in low-volume, droplet-based membranes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Imaging the initial condition of heavy-ion collisions and nuclear structure across the nuclide chart

High-energy nuclear collisions encompass three key stages: the structure of the colliding nuclei, informed by low-energy nuclear physics, the initial condition , leading to the formation of quark–gluon plasma (QGP), and the hydrodynamic expansion and hadronization of the QGP, leading to final-state hadron distributions that are observed experimentally. Recent advances in both experimental and theoretical methods have ushered in a precision era of heavy-ion collisions, enabling an increasingly accurate understanding of these stages. However, most approaches involve simultaneously determining both QGP properties and initial conditions from a single collision system, creating complexity due to the coupled contributions of these stages to the final-state observables. To avoid this, we propose leveraging established knowledge of low-energy nuclear structures and hydrodynamic observables to independently constrain the QGP’s initial condition. By conducting comparative studies of collisions involving isobar-like nuclei—species with similar mass numbers but different ground-state geometries—we can disentangle the initial condition’s impacts from the QGP properties. This approach not only refines our understanding of the initial stages of the collisions but also turns high-energy nuclear experiments into a precision tool for imaging nuclear structures, offering insights that complement traditional low-energy approaches. Opportunities for carrying out such comparative experiments at the Large Hadron Collider and other facilities could significantly advance both high-energy and low-energy nuclear physics. Additionally, this approach has implications for the future electron-ion collider. While the possibilities are extensive, we focus on selected proposals that could benefit both the high-energy and low-energy nuclear physics communities. Originally prepared as input for the long-range plan of U.S. nuclear physics, this white paper reflects the status as of September 2022, with a brief update on developments since then.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Three-dimensional Boltzmann-Hydro Code for Core-collapse in Massive Stars. II. The Implementation of Moving-mesh for Neutron Star Kicks

We present a newly developed moving-mesh technique for the multi-dimensional Boltzmann-Hydro code for the simulation of core-collapse supernovae (CCSNe). What makes this technique different from others is the fact that it treats not only hydrodynamics but also neutrino transfer in the language of the 3 + 1 formalism of general relativity (GR), making use of the shift vector to specify the time evolution of the coordinate system. This means that the transport part of our code is essentially general relativistic, although in this paper it is applied only to the moving curvilinear coordinates in the flat Minknowski spacetime, since the gravity part is still Newtonian. The numerical aspect of the implementation is also described in detail. Employing the axisymmetric two-dimensional version of the code, we conduct two test computations: oscillations and runaways of proto-neutron star (PNS). We show that our new method works fine, tracking the motions of PNS correctly. We believe that this is a major advancement toward the realistic simulation of CCSNe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

VALIDATION AND DETERMINATION OF SIGNIFICANT SIMULATION PARAMETERS USING THE SMOOTHED PARTICLE HYDRODYNAMIC CODE NEUTRINO

In recent years there has been an increase in flooding concern to nuclear power plants; thus, tools to reduce the cost of mitigation efforts are valuable to industry. Simulating flooding events is achieved through different methods and at varied levels of detail. However, these simulations must be able to match application conditions, which requires code validation as well as knowledge of which parameters affect the accuracy of the simulation. To validate and determine the significant parameters for these tools, a sloshing experiment was conducted at the George Washington University and then simulated with the smoothed particle hydrodynamic code Neutrino. The simulation setup was constructed to match the experiment as closely as possible, and the end-wall pressure results were analysed. The 90% and 50% end-wall pressure bounds were compared along with the pressure impulse. The simulation bounds fell mostly within the experimental bounds, with exceptions in insignificant low-pressure areas, such as at the beginning (~40% difference) and the end of the cycle (~40% difference). Another exception was seen at a point in the peak pressure (~75% difference). The simulated pressure impulse matched the experiments well and was within 10%. Additionally, seven parameters were qualitatively assessed and then investigated to determine their significance. The parameters’ significance was determined by sampling five values across the parameters’ range using the Risk Analysis Virtual Environment, then comparing the average pressure for each run. All seven parameters were identified as having some significance. However, the particle size, interaction-radius to particle-size ratio, and fluid settling seem to have a greater effect due to the larger fluctuation of the results. Now that the significant parameters are identified, further research is being done to quantify the significance of each parameter.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

First-principles investigations on ionization and thermal conductivity of polystyrene for inertial confinement fusion applications

Using quantum molecular-dynamics (QMD) methods based on the density functional theory, we have performed first-principles investigations of the ionization and thermal conductivity of polystyrene (CH) over a wide range of plasma conditions (ρ = 0.5 to 100 g/cm{sup 3} and T = 15 625 to 500 000 K). The ionization data from orbital-free molecular-dynamics calculations have been fitted with a “Saha-type” model as a function of the CH plasma density and temperature, which gives an increasing ionization as the CH density increases even at low temperatures (T < 50 eV). The orbital-free molecular dynamics method is only used to gauge the average ionization behavior of CH under the average-atom model in conjunction with the pressure-matching mixing rule. The thermal conductivities (κ{sub QMD}) of CH, derived directly from the Kohn–Sham molecular-dynamics calculations, are then analytically fitted with a generalized Coulomb logarithm [(lnΛ){sub QMD}] over a wide range of plasma conditions. When compared with the traditional ionization and thermal conductivity models used in radiation–hydrodynamics codes for inertial confinement fusion simulations, the QMD results show a large difference in the low-temperature regime in which strong coupling and electron degeneracy play an essential role in determining plasma properties. Hydrodynamic simulations of cryogenic deuterium–tritium targets with CH ablators on OMEGA and the National Ignition Facility using the QMD-derived ionization and thermal conductivity of CH have predicted ∼20% variation in target performance in terms of hot-spot pressure and neutron yield (gain) with respect to traditional model simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗