Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fluid wave”

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 505 records · Page 28

Symmetry Determining Equations of the Rankine-Hugoniot Equations for Variable Velocity Shock Waves

The “constant velocity piston” problem (Fig. 1), also known as the “piston problem,” is a standard model for a one dimensional, in our case linear, symmetric shock wave moving through an inviscid, perfect gas. The model can be divided into two regions - a perturbed section on the left and an unperturbed section on the right - by a moving shock wave moving left to right. Both the perturbed and unperturbed sections, i.e. the shocked and unshocked regions, respectively, obey the Eulerian conservation equations; however, at the exact location of the shock, there is a mathematical discontinuity not satisfied by the Euler equations. To ensure continuity and conservation of certain quantities when crossing between the unshocked and shocked regions, we evoke a series of equations derived from the Eulerian conservation equations, called the Rankine-Hugoniot equations, or “jump” equations as it is often referred to in the literature on the topic. The classical constant-velocity piston problem assumes the piston features a constant driving velocity (among many other willing suspensions of belief required in the pursuit of a first principles equation model); consequent to this assumption is a constant-velocity shock and a constant-velocity shocked flow state. However, using Lie Group Theory (LGT), also known as symmetry analysis, we can attempt to reinterpret the model with a shock wave of variable velocity in time and space. An extension of the model in this way opens up the possibility for obtaining new analytical solutions to the piston problem for certain shock velocity models. In this report, we use LGT to derive the symmetry determining equations (SDEs), whose solutions are Lie groups, which permit analytical solutions. In the future, we can then use the SDEs to define constraint equations on the shock velocity model and what the successive solutions to the Euler equations might be based off such constraints. This report is structured as follows: Section 2 provides a brief derivation of the Rankine-Hugoniot (“jump”) equations; Section 3 gives an overview of Lie group theory; Section 4 derives the SDEs of the jump equations; Section 5 derives the Euler conservation equations for fluids; and Section 6 presents concluding remarks and opportunities for future studies.

42 ENGINEERING↗

Fluid dynamic mathematical aspects of supernova remnants

Supernovae—explosions of stars—are a central problem in astrophysics since they contain information on the entire process of stellar evolution and nucleosynthesis. Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities, developing during the supernova blast, lead to intense interfacial RT/RM mixing of the star's materials and couple astrophysical to atomic scales. This work analyzes some fluid dynamic mathematical aspects of the titanic task of supernova's blast. We handle mathematical challenges of RT/RM dynamics in supernova relevant conditions by directly linking the conservation laws governing RT/RM dynamics to symmetry-based momentum model, by exactly deriving the model parameters in the scale-dependent and scale-invariant regimes, and by exploring the special self-similar class for RT/RM interfacial mixing with variable accelerations. Here we reveal that RT/RM dynamics is strongly influenced by deterministic (the initial and the flow) conditions in the scale-dependent linear and nonlinear regimes and in the self-similar mixing regime. The theory outcomes are consistent with the observations of supernova remnants, explain the results of the scaled laboratory experiments in high energy density plasmas, and yield the design of future experiments for the accurate quantification of RT/RM dynamics in supernova relevant conditions. We find that from fluid dynamic mathematical perspectives, supernovae can be regarded as an astrophysical initial value problem. Along with the guidance of what explodes at microscopic scales, supernova remnants encapsulate information on the explosion hydrodynamics and the associated deterministic conditions at macroscopic scales. We urge such effects be considered in interpretations of the observational data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Self-consistent interaction of linear gravitational and electromagnetic waves in non-magnetized plasma

Here, this paper explores the hybridization of linear metric perturbations with linear electromagnetic (EM) perturbations in non-magnetized plasma for a general background metric. The local wave properties are derived from first principles for inhomogeneous plasma, without assuming any symmetries of the background metric. First, we derive the effective (“oscillation-center”) Hamiltonian that governs the average dynamics of plasma particles in a prescribed quasimonochromatic wave that involves metric perturbations and EM fields simultaneously. Then, using this Hamiltonian, we derive the backreaction of plasma particles on the wave itself and obtain gauge-invariant equations that describe the resulting self-consistent gravito-electromagnetic (GEM) waves in a plasma. The transverse tensor modes of gravitational waves are found to have no interaction with the plasma and the EM modes in the geometrical-optics limit. However, for $longitudinal$ GEM modes with large values of the refraction index, the interplay between gravitational and EM interactions in plasma can have a strong effect. In particular, the dispersion relation of the Jeans mode is significantly affected by electrostatic interactions. As a spin-off, our calculation also provides an alternative resolution of the so-called Jeans swindle.

79 ASTRONOMY AND ASTROPHYSICS↗

Multidimensional Numerical Modeling of Combustion Dynamics in a Non-Premixed Rotating Detonation Engine With Adaptive Mesh Refinement

In the present work, a novel computational fluid dynamics (CFD) methodology was developed to simulate full-scale non-premixed rotating detonation engines (RDEs). A unique feature of the modeling approach was the incorporation of adaptive mesh refinement (AMR) to achieve a good trade-off between model accuracy and computational expense. Here, unsteady Reynolds-averaged Navier–Stokes (RANS) simulations were performed for an Air Force Research Laboratory (AFRL) non-premixed RDE configuration with hydrogen as fuel and air as the oxidizer. The finite-rate chemistry model, along with a ten-species detailed kinetic mechanism, was employed to describe the H 2 -Air combustion chemistry. Three distinct operating conditions were simulated, corresponding to the same global equivalence ratio of unity but different fuel/air mass flowrates. For all conditions, the capability of the model to capture essential detonation wave dynamics was assessed. An exhaustive verification and validation study was performed against experimental data in terms of a number of waves, wave frequency, wave height, reactant fill height, oblique shock angle, axial pressure distribution in the channel, and fuel/air plenum pressure. The CFD model was demonstrated to accurately predict the sensitivity of these wave characteristics to the operating conditions, both qualitatively and quantitatively. A comprehensive heat release analysis was also conducted to quantify detonative versus deflagrative burning for the three simulated cases. The present CFD model offers a potential capability to perform rapid design space exploration and/or performance optimization studies for realistic full-scale RDE configurations.

42 ENGINEERING↗

Deuterium and tritium anomalous transport in the generalized Hasegawa–Wakatani resistive drift wave turbulence model with finite ion Larmor radius

Anomalous transport of mixed deuterium–tritium plasma in the edge of magnetic fusion reactors is investigated using numerical solutions of resistive drift wave turbulence model equations, including finite Larmor radius effects, that are derived within the generalized Hasegawa–Wakatani framework. The anomalous cross field diffusivities of deuterium and tritium are compared in turbulence regimes with different values of the electron adiabaticity parameter controlling the existence of zonal flow. The dependence of the tritium-to-deuterium diffusivity ratio on the deuterium and tritium densities and the logarithmic density gradients is analyzed, and a scaling relation is obtained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Regolith, rock and fluid distributions at the Upper Colorado River Basin via a multicomponent seismic imaging approach

In October 2018, we conducted an exploratory seismic survey of the East River Watershed Function Science Focus Area, near Crested Butte, CO. Our survey goals were to map seismic bedrock, regolith and sediment properties and relate them to geomorphic and hydrologic processes. We acquired ~12 km of seismic data along roads with a 72-channel, 1.25m spaced streamer, an average shot spacing of 2.5m and an average acquisition rate of 3 km/day. We also obtained ~3km of planted geophone data with geophone and shot spacing between 2.5 and 5 meters and an acquisition rate of 0.5 km/day. Vp results show a sharp transition between regolith and bedrock along hillslopes suggesting a similarly abrupt change in porosity. Vp anisotropy in the Mancos Shale is observable in the lower reaches of the watershed, but is absent in regions near crystalline intrusions where contact metamorphism has likely modified pore structure within the shale. Alluvium thickness varies along the river axis from 0 to > 20 meters reflecting variations in depositional processes as function of stream gradient. We show that when strong discontinuities at the regolith base are present, inhomogeneous waves that travel along subsurface interfaces are generated and can be measured on vertical component geophones. We present a summary manuscript in the report appendix.

58 GEOSCIENCES↗

Design of High-deflection Foils for Marine Hydrokinetic Applications (Final Report)

This project used model-scale tank testing and fluid-structure-interaction (FSI) simulations to investigate the behavior of foils with large deflections and the effect of these deflections on crossflow turbine performance with the goal of determining the maximum allowable deflections consonant with efficiency and a robust, durable structure. A validated modeling and simulation approach was used in the design of an Ocean Renewable Power Company (ORPC) full-scale turbine. The methodology is applicable to the design of other marine hydrokinetic (MHK) devices, and to hydrokinetic turbines which experience significant deflections during operation.

16 TIDAL AND WAVE POWER↗

Relating Hydro–Mechanical and Elastodynamic Properties of Dynamically Stressed Tensile–Fractured Rock in Relation to Applied Normal Stress, Fracture Aperture, and Contact Area

We exploit nonlinear elastodynamic properties of fractured rock to probe the micro-scale mechanics of fractures and understand the relation between fluid transport and fracture aperture under dynamic stressing. Experiments were conducted on rough, tensile-fractured Westerly granite subject to triaxial stresses. We measure fracture permeability for steady-state fluid flow with deionized water. Pore pressure oscillations are applied at amplitudes ranging from 0.2 to 1 MPa at 1 Hz frequency. During dynamic stressing we transmit ultrasonic signals through the fracture using an array of piezoelectric transducers (PZTs) to monitor evolution of interface properties. We examine the influence of fracture aperture and contact area by conducting measurements at effective normal stresses of 10–20 MPa. Additionally, the evolution of contact area with stress is characterized using pressure sensitive film. These experiments are conducted separately with the same fracture and map contact area at stresses from 9 to 21 MPa. The measurements are a proxy for “true” contact area for the fracture surface and we relate them to elastic properties using the calculated PZT sensor footprints via numerical modeling of Fresnel zones. We compare the elastodynamic response of the fracture using the stress-induced changes in ultrasonic wave velocities for transmitter-receiver pairs to image spatial variations in contact properties. We show that nonlinear elasticity and permeability enhancement decrease with increasing normal stress. Additionally, post-oscillation wave velocity and permeability exhibit quick recoveries toward pre-oscillation values. Estimates of fracture contact area (global and local) demonstrate that the elastodynamic and permeability responses are dominated by fracture topology.

58 GEOSCIENCES↗

Autonomous sensor suite for evaluating fish-turbine interactions and environmental impacts in marine renewable energy and hydropower

Marine renewable energy (MRE) harnesses ocean-based resources such as waves, tides, currents, and thermal or salinity gradients for sustainable power generation. It has the potential to complement existing renewable resources, support remote communities, and contribute to decarbonization efforts. However, understanding the hydrodynamic forces created by MRE devices and their impacts on marine life is critical for responsible deployment. Here, to address these concerns, advanced sensor devices, including the Marine Sensor Fish (MSF), Sensor Fish Mini (SF Mini), and Flexible Sensor Fish (FSF), were developed to measure interactions between aquatic organisms and MRE systems. This paper details the design, manufacturing, calibration, and field deployment of these sensor suites, highlighting their ability to capture key physical stressors such as shear forces, pressure changes, and collision impacts. The MSF successfully evaluated turbine interactions at a tidal turbine in the Salish Sea, capturing data on turbulence, collision impact, and pressure gradients. The SF Mini validated hydrodynamic conditions in scaled hydraulic models, supporting computational fluid dynamics simulations. The FSF, with its flexible silicone body, measured species-specific impacts in turbulent environments. This research demonstrates the potential of Sensor Fish technology to advance sustainable marine energy systems by reducing biological impacts and informing environmentally sustainable designs.

Ecological impacts↗

Techniques for studying materials under extreme states of high energy density compression

The properties of materials under extreme conditions of pressure and density are of key interest to a number of fields, including planetary geophysics, materials science, and inertial confinement fusion. In geophysics, the equations of state of planetary materials, such as hydrogen and iron, under ultrahigh pressure and density provide a better understanding of their formation and interior structure [Celliers et al., “Insulator-metal transition in dense fluid deuterium,” Science 361, 677–682 (2018) and Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)]. The processes of interest in these fields occur under conditions of high pressure (100 GPa–100 TPa), high temperature (>3000 K), and sometimes at high strain rates (>103 s−1) depending on the process. With the advent of high energy density (HED) facilities, such as the National Ignition Facility (NIF), Linear Coherent Light Source, Omega Laser Facility, and Z, these conditions are reachable and numerous experimental platforms have been developed. To measure compression under ultrahigh pressure, stepped targets are ramp-compressed and the sound velocity, measured by the velocity interferometer system for any reflector diagnostic technique, from which the stress-density of relevant materials is deduced at pulsed power [M. D. Knudson and M. P. Desjarlais, “High-precision shock wave measurements of deuterium: Evaluation of exchange-correlation functionals at the molecular-to-atomic transition,” Phys. Rev. Lett. 118, 035501 (2017)] and laser [Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)] facilities. To measure strength under high pressure and strain rates, experimenters measure the growth of Rayleigh–Taylor instabilities using face-on radiography [Park et al., “Grain-size-independent plastic flow at ultrahigh pressures and strain rates,” Phys. Rev. Lett. 114, 065502 (2015)]. The crystal structure of materials under high compression is measured by dynamic x-ray diffraction [Rygg et al., “X-ray diffraction at the national ignition facility,” Rev. Sci. Instrum. 91, 043902 (2020) and McBride et al., “Phase transition lowering in dynamically compressed silicon,” Nat. Phys. 15, 89–94 (2019)]. Medium range material temperatures (a few thousand degrees) can be measured by extended x-ray absorption fine structure techniques, Yaakobi et al., “Extended x-ray absorption fine structure measurements of laser-shocked V and Ti and crystal phase transformation in Ti,” Phys. Rev. Lett. 92, 095504 (2004) and Ping et al., “Solid iron compressed up to 560 GPa,” Phys. Rev. Lett. 111, 065501 (2013), whereas more extreme temperatures are measured using x-ray Thomson scattering or pyrometry. This manuscript will review the scientific motivations, experimental techniques, and the regimes that can be probed for the study of materials under extreme HED conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Virtual Flow Solver - Geophysics: A 3D Incompressible Navier-Stokes Solver

Virtual Flow Solver - Geophysics (VFS-Geophysics) is a three-dimensional (3D) incompressible Navier-Stokes solver based on the Curvilinear Immersed Boundary (CURVIB) method. The CURVIB is a sharp interface type of immersed boundary (IB) method that enables the simulation of fluid flows in the presence of geometrically complex moving bodies. The CURVIB method can be applied to wind/MHK turbine simulations and energy applications. VFS-Geophysics is the result of many years of research work by several graduate students, post-docs, and research associates that have been involved in the Computational Hydrodynamics and Biofluids Laboratory directed by Professor Fotis Sotiropoulos. The preparation of the present manual has been supported by the U.S. Department of Energy (DE-EE 0005482).

16 TIDAL AND WAVE POWER↗

Zonal Flow and Self-regulating Mechanism in a Hydrodynamic Disk

This study addresses key aspects of momentum transport in hydrodynamic disks, which is critical for understanding zonal flow generation and turbulence in compressible hydrodynamic disks. We find that nonlinear momentum/density transport leads to the formation of zonal flows from the Rossby wave instability in disks. We analytically derive the generation and location of zonal flows and describe a modified Taylor identity applicable to compressible disk flows. We further present a self-regulation model, revealing a dynamic interplay between zonal flow and fluctuations driven by Rossby wave instability that regulates the nonlinear saturation state. This theoretical framework contributes insights into the dynamics of disks such as protoplanetary disks, shedding light on the intricate processes governing momentum/density transport and the emergence of zonal flows in the saturation of protoplanetary disks.

Astrophysical fluid dynamics↗

The kinetic analog of the pressure–strain interaction

Energy transport in weakly collisional plasma systems is often studied with fluid models and diagnostics. However, the applicability of fluid models is limited when collisions are weak or absent, and using a fluid approach can obscure kinetic processes that provide key insights into the physics of energy transport. Kinetic diagnostics retain all of the information in 3D-3V phase space and thereby reach beyond the insights of fluid models to elucidate the mechanisms responsible for collisionless energy transport. In this work, we derive the Kinetic Pressure–Strain (KPS): a kinetic analog of the pressure–strain interaction, which is the channel between flow energy density and internal energy density in fluid models. Through two case studies of electron Landau damping, we demonstrate that the KPS diagnostic can elucidate kinetic mechanisms that are responsible for energy transport in this channel, just as the related field–particle correlation is known to identify kinetic mechanisms of transport between electromagnetic field energy density and kinetic energy density in particle flows. In addition, we show that resonant electrons play a major role in transferring energy between fluid flows and internal energy during the process of Landau damping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The effect of blue and infrared laser melting frequency on oxide morphology in 304L

Continuous Wave Laser Beam (LB) melting offers control over the localised heating and cooling of melt pools for welding, brazing, additive manufacturing and solidification. Research into laserliquid metal interactions has primarily focused on the heat and mass transport under large thermal gradients imposed by the localised melting conditions. Recent commercial availability of blue (450 nm) wavelength lasers has enabled the printing of material with poor IR (1064) nm attenuation like Cu. However, little research has been conducted into varying input laser frequency and scan rate at fixed absorption to understand the effects of laser light on surface oxide formation and the underlying liquid metal thermodynamic response. Here, this article conducts laser melting of 304L under Blue (450 nm) and Infrared (IR, 1064 nm) laser frequencies and examines their impact on laser oxide thickness, chemistry and coloration. We find that laser frequencies induce changes in the oxide layer thickness and chemistry that cannot be explained using conventional thermal absorption shifts and fluid dynamics. We suggest that light coupling may have thermodynamic implications for the chemical potential of the liquid metal, which may drive the observed phase behaviour. Identifying mechanisms for the tuning of the chemical potential via laser frequency control could enable the ability to control solidification microstructure and atomic ordering in laser welding and additive manufacturing processes.

304L↗

Simulation and flow physics of a shocked and reshocked high-energy-density mixing layer

This work describes a computational investigation of multimode instability growth and multimaterial mixing induced by multiple shock waves in a high-energy-density (HED) environment, where pressures exceed 1 Mbar. The simulations are based on a series of experiments performed at the National Ignition Facility (NIF) and designed as an HED analogue of non-HED shock-tube studies of the Richtmyer–Meshkov instability and turbulent mixing. A three-dimensional computational modelling framework is presented. It treats many complications absent from canonical non-HED shock-tube flows, including distinct ion and free-electron internal energies, non-ideal equations of state, radiation transport and plasma-state mass diffusivities, viscosities and thermal conductivities. The simulations are tuned to the available NIF data, and traditional statistical quantities of turbulence are analysed. Integrated measures of turbulent kinetic energy and enstrophy both increase by over an order of magnitude due to reshock. Large contributions to enstrophy production during reshock are seen from both the baroclinic source and enstrophy–dilatation terms, highlighting the significance of fluid compressibility in the HED regime. Dimensional analysis reveals that Reynolds numbers and diffusive Péclet numbers in the HED flow are similar to those in a canonical non-HED analogue, but conductive Péclet numbers are much smaller in the HED flow due to efficient thermal conduction by free electrons. It is shown that the mechanism of electron thermal conduction significantly softens local spanwise gradients of both temperature and density, which causes a minor but non-negligible decrease in enstrophy production and small-scale mixing relative to a flow without this mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Semi-implicit method to solve compressible multiphase fluid flows without acoustic time step restrictions

A new 3D multiphase numerical capability is presented here for simulating multiphase flow regimes at all Mach numbers (M). The new method is a semi-implicit extension of the finite volume discrete equation method (DEM) of Chinnaya et al. (2004), which originally used explicit time-stepping. The capability is also developed to work with another extension of the DEM to moving grids for arbitrary Lagrangian-Eularian (ALE) methods detailed in Dunn (2011). Rather than solving all phase equations simultaneously, the DEM reduces the equations to a system of single-phase Riemann solves, where each phase has its own velocity and thermodynamic state. Exchanges between the phases are modeled through source terms accounting for the phase interactions. Furthermore, since the original multiphase scheme uses an explicit time-advancement scheme, it has time step restrictions dictated by the speed of sound, which limits the model's ability to simulate weakly compressible flows which typically need to be integrated for longer time periods. Here, we extend the current multiphase formulation by implementing a pressure-correcting step to enable implicit calculations and remove acoustic time step limitations. The new semi-implicit algorithm allows use of relatively large time steps compared to an explicit method. Validation and benefits of the new implicit time-step method are illustrated using several examples including weakly compressible flows and strong shock waves. The scheme presented here is general and may be used for a variety of applications which require capabilities for handling multiphase flow at a wide range of Mach numbers. However, the goal of this research is to simulate all stages of high energy explosions, including the shock formation (high Mach numbers) and evolution of the buoyant cloud (low Mach numbers).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Closure models for the feedback of energetic particles on plasma turbulence

Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.

79 ASTRONOMY AND ASTROPHYSICS↗

Flow Channels and the Generation of Alfvénic Turbulence Along Storm‐Time Inner Magnetospheric Field‐Lines

Abstract A feature of Earth's storm‐time magnetosphere outside the plasmapause is the occurrence of broad‐spectrum Alfvénic fluctuations. In this letter observations from the Van Allen Probes are compared with 3‐D fluid‐kinetic simulations of an evolving convective flow channel to investigate the mechanisms generating the observed spectrum. It is shown how narrow channels of fast convection are unstable to the Kelvin‐Helmholtz instability which on closed field‐lines initiates a cascade to small scales. Sustained driving of the flow combined with reflection from the topside ionosphere leads to the generation of an intensified spectrum of electromagnetic structures having similar spectral and morphological characteristics to those observed. This process couples enhanced magnetospheric convection to kinetic scale electromagnetic fluctuations that drive particle transport, scattering and energization through the outer radiation belt and ring current during geomagnetic storms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗