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 289 records · Page 16

Apparently superluminal superfluids

We consider the superfluid phase of a specific renormalizable relativistic quantum field theory. We prove that, within the regime of validity of perturbation theory and of the superfluid effective theory, there are consistent and regular vortex solutions where the superfluid’s velocity field as traditionally defined smoothly interpolates between zero and arbitrarily large superluminal values. We show that this solution is free of instabilities and of superluminal excitations. We show that, in contrast, a generic vortex solution for an ordinary fluid does develop an instability if the velocity field becomes superluminal. All this questions the characterization of a superfluid velocity field as the actual velocity of “something”.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Mean-field dynamo as a quantum-like modulational instability

Presented here is a novel formulation of the mean-field dynamo as a modulational instability of magnetohydrodynamic (MHD) turbulence. This formulation, termed mean-field wave kinetics (MFWK), is based on the Weyl symbol calculus and allows describing the interaction between the mean fields (magnetic field and fluid velocity) and turbulence without requiring scale separation that is commonly assumed in the literature. The turbulence is described by the Wigner–Moyal equation for the spectrum of the two-point correlation matrix (Wigner matrix) of magnetic-field and velocity fluctuations and depicts the turbulence as an effective plasma of quantum-like particles that interact via the mean fields. Eddy–eddy interactions, which serve as ‘collisions’ in this effective plasma, are modelled within the standard minimal tau approximation to aid comparison with existing theories. Using MFWK, the non-local electromotive force is calculated for generic turbulence from first principles, modulo the limitations of MFWK. This result is then used to study, both analytically and numerically, the modulational modes of MHD turbulence, which appear as linear instabilities of the said effective quantum-like plasma of fluctuations. The standard α 2 -dynamo and other known results are reproduced as special cases. A new dynamo effect is predicted that is driven by correlations between the turbulent flow velocity and the turbulent current.

astrophysical plasmas↗

Development of five-moment two-fluid modeling for Z-pinch physics

The Z-pinch m = 0 instability as well as its stabilization by radially sheared axial flow is studied using the nonlinear ideal five-moment twofluid (5M2F) model with an extension of that model to include Braginskii heat and momentum transport. Using the ideal 5M2F model, linear growth rate results are compared with prior work using MHD and Hall MHD. At small normalized wavenumber, 1 < k a < 4 , where a is the effective pinch radius, 5M2F results agree with Hall MHD within 20% in scenarios without radially sheared axial flow. With the sheared flow and focusing on ka = 10/3, agreement with Hall MHD is excellent. In the limit of small ion inertial length, results also match with MHD. A comparison with PIC modeling of shear-free m = 0 stability focuses on a plasma scenario based on recent experimental results. In a scan of mode wavenumber, ideal 5M2F results are qualitatively similar to PIC: the growth rate rises to a peak at a moderate wavenumber and declines at a large wavenumber in contrast to MHD results, which show the saturation of the growth rate with the increasing wavenumber rather than a decline. The peak normalized 5M2F growth rate is &#x3B3; &#x3C4; A = 1.5 , where sA is the Alfven transit time across the pinch. The peak occurs at normalized wavenumber ka = 10. For comparison, PIC results have a peak growth of &#x3B3; &#x3C4; A = 0.77 at ka = 5. Including Braginskiibased closure of the 5M2F model does not qualitatively change the ideal results in this particular case. Nonlinear saturation is studied using the 5M2F model with the dissipative Braginskii-based closure in cases with pinch-edge sheared-flow speed equal to half the Alfven speed. Nonlinear mixing due to the sheared flow yields a quasi-steady state after modest losses of pinch ion inventory and pinch thermal energy, approximately 30% and 10%, respectively. 5M2F modeling captures the essential physics of m = 0 instability and offers a computationally tractable route to high-fidelity modeling of 3D Z-pinch behavior, including m = 1 instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The dynamics, mixing, and thermonuclear burn of compressed foams with varied gas fills

Inertial confinement fusion (ICF) implosions involve highly coupled physics and complex hydrodynamics that are challenging to model computationally. Due to the sensitivity of such implosions to small features, detailed simulations require accurate accounting of the geometry and dimensionality of the initial conditions, including capsule defects and engineering features such as fill tubes used to insert gas into the capsule, yet this is computationally prohibitive. It is therefore difficult to evaluate whether discrepancies between the simulation and experiment arise from inadequate fidelity to the capsule geometry and drive conditions, uncertainties in physical data used by simulations, or inadequate physics. We present results from detailed high-resolution three-dimensional simulations of ICF implosions performed as part of the MARBLE campaign on the National Ignition Facility [Albright et al., Phys. Plasmas 29, 022702 (2022)]. These experiments are foam-filled separated-reactant experiments, where deuterons reside in the foam and tritons reside in the capsule gas fill and deuterium–tritium (DT) fusion reactions only occur in the presence of mixing between these materials. Material mixing in these experiments is primarily seeded by shock interaction with the complex geometry of the foam and gas fill, which induces the Richtmyer–Meshkov instability. We compare results for experiments with two different gas fills (ArT and HT), which lead to significant differences in the hydrodynamic and thermodynamic developments of the materials in the implosion. Our simulation results show generally good agreement with experiments and demonstrate a substantial impact of hydrodynamic flows on measured ion temperatures. The results suggest that viscosity, which was not included in our simulations, is the most important unmodeled physics and qualitatively explains the few discrepancies between the simulation and experiment. The results also suggest that the hydrodynamic treatment of shocks is inadequate to predict the heating and yield produced during shock flash, when the shock converges at the center of the implosion. Alternatively, underestimation of the level of radiative preheat from the shock front could explain many of the differences between the experiment and simulation. Nevertheless, simulations are able to reproduce many experimental observables within the level of experimental reproducibility, including most yields, time-resolved X-ray self-emission images, and an increase in burn-weighted ion temperature and neutron down-scattered ratio in the line of sight that includes a jet seeded by the glue spot that joins capsule hemispheres.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

HERO WEC V1.0 - WEC-Sim Model

This zip file contains the files that are needed to simulate NREL's HERO WEC (hydraulic and electric reverse osmosis wave energy converter). This requires the user to have already installed WEC-Sim. In addition to the standard toolboxes that are required to run WEC-Sim the user will also need the Simscape Fluids and Simscape Driveline packages. In the zip file you will find the following: - HEROV1_HPTO.slx: Simulink-based WEC Sim model of the first gen (V1.0) Hydraulic PTO (power take-off) that was designed for the HERO WEC - wecSimInputFile.m: Input file needed to run the model - userDefinedFunctionsMCR.m: MCR (multi condition run) script that is needed if a use wants to simulate multiple wave conditions. - geometry (folder): Includes the geometry file that is needed for visualization - hydroData (folder): Includes the required WAMIT data to run WEC-Sim

16 TIDAL AND WAVE POWER↗

Validation of numerical models for seismic fluid-structure-interaction analysis of nuclear, safety-related equipment

Seismic design, qualification, and risk assessment of nuclear safety-related vessels filled with liquid will have to consider the interaction between the vessel (tank), the contained liquid, and submerged components, if any. Seismic fluid-structure-interaction (FSI) analysis of nuclear vessels will rely on numerical models, which are required to be verified and validated. This study validates previously-verified numerical models using test data generated from earthquake-simulator experiments involving a cylindrical tank. Two solvers in LS-DYNA for FSI simulations are used for the numerical analysis: Arbitrary-Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD). Numerical and test results are compared for responses critical to seismic design of advanced reactor vessels: hydrodynamic pressures on the tank wall, reactions at the support, and wave heights of the contained liquid. Analysis is performed for one-, two-, and three-directional seismic inputs with a range of intensity, and rocking motions. The accuracy of the numerical results and different methods for outputting wave heights are discussed. Recommendations for validation of seismic FSI numerical models of advanced reactors are provided. As a result, the validation exercise presented in this study is broadly applicable to cylindrical tanks, regardless of industry sector.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

A Plasma Modeling Hierarchy and Verification Approach

This report reviews a hierarchy of formal mathematical models for describing plasma phenomena. Starting with the Boltzmann equation, a sequence of approximations and modeling assumptions can be made that progressively reduce to the equations for magnetohydrodynamics. Understanding the assumptions behind each of these models and their mathematical form is essential to appropriate use of each level of the hierarchy. A sequence of moment models of the Boltzmann equation are presented, then focused into a generalized three-fluid model for neutral species, electrons, and ions. This model is then further reduced to a two-fluid model, for which Braginskii described a useful closure. Further reduction of the two-fluid model yields a Generalized Ohm's Law model, which provides a connection to magnetohydrodynamic approaches. A verification approach based on linear plasma waves is presented alongside the model hierarchy, which is intended as an initial and necessary but not sufficient step for verification of plasma models within this hierarchy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Identification of Kelvin-Helmholtz generated vortices in magnetised fluids

The Kelvin-Helmholtz Instability (KHI), arising from velocity shear across the magnetopause, plays a significant role in the viscous-like transfer of mass, momentum, and energy from the shocked solar wind into the magnetosphere. While the KHI leads to growth of surface waves and vortices, suitable detection methods for these applicable to magnetohydrodynamics (MHD) are currently lacking. A novel method is derived based on the well-established λ-family of hydrodynamic vortex identification techniques, which define a vortex as a local minimum in an adapted pressure field. The J × B Lorentz force is incorporated into this method by using an effective total pressure in MHD, including both magnetic pressure and a pressure-like part of the magnetic tension derived from a Helmholtz decomposition. The λ MHD method is shown to comprise of four physical effects: vortical momentum, density gradients, fluid compressibility, and the rotational part of the magnetic tension. A local three-dimensional MHD simulation representative of near-flank magnetopause conditions (plasma β’s 0.5 – 5 and convective Mach numbers M f ∼ 0.4) under northward interplanetary magnetic field (IMF) is used to validate λ MHD . Analysis shows it correlates well with hydrodynamic vortex definitions, though the level of correlation decreases with vortex evolution. Overall, vortical momentum dominates λ MHD at all times. During the linear growth phase, density gradients act to oppose vortex formation. By the highly nonlinear stage, the formation of small-scale structures leads to a rising importance of the magnetic tension. Compressibility was found to be insignificant throughout. Finally, a demonstration of this method adapted to tetrahedral spacecraft observations is performed.

79 ASTRONOMY AND ASTROPHYSICS↗

Rayleigh–Taylor and Richtmyer–Meshkov instabilities: A journey through scales

Hydrodynamic instabilities such as Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities usually appear in conjunction with the Kelvin–Helmholtz (KH) instability and are found in many natural phenomena and engineering applications. They frequently result in turbulent mixing, which has a major impact on the overall flow development and other effective material properties. This can either be a desired outcome, an unwelcome side effect, or just an unavoidable consequence, but must in all cases be characterized in any model. The RT instability occurs at an interface between different fluids, when the light fluid is accelerated into the heavy. The RM instability may be considered a special case of the RT instability, when the acceleration provided is impulsive in nature such as that resulting from a shock wave. Here in this pedagogical review, we provide an extensive survey of the applications and examples where such instabilities play a central role. First, fundamental aspects of the instabilities are reviewed including the underlying flow physics at different stages of development, followed by an overview of analytical models describing the linear, nonlinear and fully turbulent stages. RT and RM instabilities pose special challenges to numerical modeling, due to the requirement that the sharp interface separating the fluids be captured with fidelity. These challenges are discussed at length here, followed by a summary of the significant progress in recent years in addressing them. Examples of the pivotal roles played by the instabilities in applications are given in the context of solar prominences, ionospheric flows in space, supernovae, inertial fusion and pulsed-power experiments, pulsed detonation engines and Scramjets. Progress in our understanding of special cases of RT/RM instabilities is reviewed, including the effects of material strength, chemical reactions, magnetic fields, as well as the roles the instabilities play in ejecta formation and transport, and explosively expanding flows. The article is addressed to a broad audience, but with particular attention to graduate students and researchers who are interested in the state-of-the-art in our understanding of the instabilities and the unique issues they present in the applications in which they are prominent.

42 ENGINEERING↗

Drift kinetic effects on plasma response to resonant magnetic perturbation for EU DEMO design

Abstract A systematic investigation of the plasma response to an externally applied resonant magnetic perturbation (RMP) field, for the purpose of controlling edge localized modes, is carried out for an EU DEMO reference plasma. Particular emphasis is placed on the role of kinetic effects associated with both thermal particles and fusion-born alphas. The single fluid, resistive model predicts a large peak amplification of the n = 1 ( n is the toroidal mode number) plasma response to the target equilibrium, which is found to be close to the Troyon no-wall limit. A more advanced response model, including kinetic resonances between the RMP perturbation and drift motions of thermal and energetic particles, on the other hand, finds a strong suppression of the n = 1 field amplification. A major role is played by the precessional drift resonance of fusion-born alphas. A strong parallel sound wave damping model is found to well reproduce the full kinetic response results for the DEMO plasma, in terms of both the resonant field response amplitude and the plasma displacement. Finally, both fluid and kinetic models produce similar responses for the n = 2 and 3 RMP fields for the considered DEMO plasma, whilst kinetic effects again become important for the n = 4 RMP due to proximity of the reference plasma to the no-wall limit for the n = 4 ideal kink instability.

Physics↗

Wave steepening and shock formation in ultracold neutral plasmas

We present observations of wave steepening and signatures of shock formation during expansion of ultracold neutral plasmas formed with an initial density distribution that is centrally peaked and decays exponentially with distance. The plasma acceleration and velocity decrease at large distance from the plasma center, leading to central ions overtaking ions in the outer regions and the development of a steepening front that is narrow compared to the size of the plasma. The density and velocity change dramatically across the front, and significant heating of the ions is observed in the region of steepest gradients. For a reasonable estimate of electron temperature, the relative velocity of ions on either side of the front modestly exceeds the local sound speed (Mach number M≳1). This indicates that by sculpting steep density gradients, it is possible to create the conditions for shock formation, or very close to it, opening a new avenue of research for ultracold neutral plasmas.

Computational fluid dynamics↗

Feasibility of an experiment on clumping induced by the Crow instability along a shocked cylinder

The growth of three-dimensional perturbations subject to the Crow instability along a vortex dipole resulting from the passage of a shock wave through a heavy gaseous cylinder is examined numerically. A linear stability analysis is performed based on geometric parameters extracted from two-dimensional simulations to determine the range of unstable wavenumbers, which is found to extend from 0.0 to 1.3 when normalized by the core separation distance. The analysis is then verified by comparison to three-dimensional simulations, which clearly show the development of the instability and the pinch-off of the vortex dipole into isolated vortex rings, which manifest as clumps of the original cylinder material. A scaling law is developed to determine the relevant spatiotemporal scales of the instability development, which is then used to assess the feasibility of a high-energy-density experiment visualizing clump formation. Specifically, a shocked cylinder with an initial diameter of 100 μm consisting of a perturbation of approximate wavelength and amplitude of 600 and 10 μm, respectively, is expected to form clumps resulting from the Crow instability approximately 40 ns after it is shocked, with dynamics which can be readily visualized on the Omega EP laser facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

HERO WEC V1.0 - WEC-Sim Model (July 2024)

**This submission supersedes submission MHKDR-483** This submission file contains the files that are needed to simulate NREL's HERO WEC (hydraulic and electric reverse osmosis wave energy converter). This requires the user to have already installed WEC-Sim. In addition to the standard toolboxes that are required to run WEC-Sim the user will also need the Simscape Fluids and Simscape Driveline packages. The zip file (HERO_V1_WECSim_2024.zip) contains the following: - HERO_HPTO_2024.slx: Simulink-based WEC Sim model of the first gen (V1.0) Hydraulic PTO (power take-off) that was designed for the HERO WEC. This model has been updated since submission #483 based on in-laboratory experimental results. - wecSimInputFile.m: Input file needed to run the model - userDefinedFunctionsMCR.m: MCR (multi condition run) script that is needed if a use wants to simulate multiple wave conditions. - geometry (folder): Includes the geometry file that is needed for visualization - hydroData (folder): Includes the required WAMIT data to run WEC-Sim -HydVisualization.mlx: Visualization script to plot simulation results (not needed to run)

16 TIDAL AND WAVE POWER↗

Impact of Numerical Hydrodynamics in Turbulent Mixing Transition Simulations

Underresolved simulations are unavoidable in high Reynolds (Re) and Mach (Ma) number turbulent flow applications at scale. Implicit large-Eddy simulation (ILES) often becomes the effective strategy to capture the dominating effects of convectively driven flow instabilities. We evaluate the impact of three distinct numerical strategies in simulations of transition and turbulence decay with ILES: the Harten–Lax–van Leer (HLL) Riemann solver applying Strang splitting and a Lagrange-plus-Remap formalism to solve the directional sweep—denoted split; the Harten–Lax–Van Leer-Contact (HLLC) Riemann solver using a directionally unsplit strategy and parabolic reconstruction—denoted unsplit; and the HLLC Riemann solver using unsplit and a low-Ma correction (LMC)—denoted unsplit*. Three case studies are considered: (1) a shock tube problem prototyping shock-driven turbulent mixing, (2) the Taylor–Green Vortex (TGV) prototyping transition to turbulence, and, (3) an homogeneous isotropic turbulence (HIT) case, focusing on the impact of discretization on transition and decay from fixed well-characterized initial conditions. Significantly more accurate predictions are provided by the unsplit schemes, in particular, when augmented with the LMC. For given resolution, only the unsplit schemes predict the turbulent mixing transition after reshock observed in the shock tube experiments. Relevant comparisons of ILES based on Euler and Navier–Stokes equations addressing potential occurrence of low-Re regimes in the applications are presented. Unsplit schemes are instrumental in allowing to capture the spatial development of the TGV flow and its validation at prescribed Re with significantly less resolution. HIT analysis confirms higher simulated turbulence Re and increased small-scale content associated with the unsplit discretizations.

74 ATOMIC AND MOLECULAR PHYSICS↗

Simulating the filamentation of smoothed laser beams with three-dimensional nonlinear dynamics

In a plasma, the ponderomotive force of an inhomogeneous electromagnetic field expels plasma from regions of high intensity. When a laser propagates through a plasma, this force creates density wells that subsequently modify the index of refraction. The beam refracts and focuses into these wells and may filament. In extreme cases, the laser beam will spray due to increasing angular divergence of the beam. The threshold for ponderomotive self-focusing is well established for isolated laser hotspots or speckles. Here, we define a practical threshold for characterizing the filamentation of thousands of speckles that are found in the focal plane of high-power laser beams spatially smoothed with random phase plates as used at high energy and power laser facilities studying inertial confinement fusion. This threshold is tested against three-dimensional simulations of speckled laser light propagating through plasma. Four metrics are applied to assess filamentation: the fraction of power above five times the average intensity, an effective f-number, the mean-squared perpendicular wavenumber, and the fraction of rarefied density with deviation from the initial density exceeding | δn/ n| = 0.1. The speckled beams studied are generated by random phase plates, both with and without additional polarization smoothing, in a parameter regime of relevance to indirect drive experiments. While filamentation has been discussed extensively in the literature, we believe this to be the first published simulation study with three-dimensional nonlinear hydrodynamics that addresses the onset threshold of ponderomotive filamentation and establishes the lengths and time scales necessary to reach a statistical steady state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fabrication of freeform potassium dihydrogen phosphate crystals by belt-on-wheel polishing for spatially tailored polarization control in high-power lasers

Large-aperture (>30 cm) optical components that provide spatially tailored control of the properties of a laser beam, such as phase, polarization, and amplitude, have been the focus of much research and development in recent years. Such optics can improve energy throughput and target implosion efficiency in high-power laser systems conducting inertial confinement fusion research. Here, a method is demonstrated for fabricating freeform surface topographies into the birefringent crystalline material, potassium dihydrogen phosphate (KDP). Belt-on-wheel polishing with an oil-based fluid, developed to mitigate the deliquescent properties of KDP, is used to deterministically polish a freeform topography into KDP, enabling the fabrication of a wave plate with spatially arbitrary retardance. Testing of the laser-induced damage threshold of belt-on-wheel polished KDP was conducted at a wavelength of 351 nm and a pulse duration of 1 ns. The results revealed that the polished surfaces are highly resistant to laser-induced damage, making them suitable for large-aperture, high-power laser applications.

Urban, Nathaniel D. [Univ. of Rochester, NY (Unite↗

Numerical Investigation of a Small Low-Flow Marine Hydrokinetic (MHK) Turbine for Small Autonomous Unmanned Mobile Recharge Stations (CRADA Final Report)

This project will improve the hydrodynamic performance of the undershot waterwheel turbine with a flow concentrator designed by the Participant, using a computational fluid dynamics study. The influence of different result parameters including the number of blades, turbine profile, submerged depth, tip speed ratio, concentrator geometry and wave + current conditions will be assessed to maximize power performance and minimize structural loads.

16 TIDAL AND WAVE POWER↗

Erosion-Corrosion Detection in Carbon Steel Pipe Loops using SRNLs Thickness and Mass Loss Measurement Coupons - 20464

Structural health of waste transfer infrastructure is of utmost importance to DOE and communities surrounding their active sites. A Fitness-for-Service (FFS) program for the Waste Transfer System has been implemented by the US Department of Energy (DOE) for the purposes of evaluating and quantifying the degradation of critical components. Currently, the method of approach is to section and evaluate pipelines that have been removed at end service. In order to augment data collection methods applied to making lifespan estimates of waste transfer system components, Florida International University (FIU) has been investigating new-technology based strategies for quantifying the rate of degradation. Present research considers the option of a coupon system installed on a section of pipe at points of interest that will provide easy access to high resolution measurements of erosion/corrosion rates during operation. Researchers at Savannah River National Laboratory (SRNL) have developed this novel solution for in situ erosion and corrosion measurements in many applicable locations such as tanks, pipes, ducts, and injection points. The coupons are approximately 12.7 mm (1/2 inch) in diameter and 25.4 mm (1 inch) in length and can be placed in the pipe sections using Swagelok compression fittings. The coupon design allows for monitoring of the total degradation resulting from corrosion and erosion. High-resolution Ultrasonic Thickness (UT) testing allows in-situ monitoring of both degradation mechanisms. When removed from the system mass loss resulting from corrosion and erosion may also be evaluated. Experimental testing at FIU implemented the SRNL coupon system in a 6.10 meter (20 foot) pipe loop system comprised of 50.8 mm (2 inch) and 76.2 mm (3 inch) pipe sections in order to replicate expected use cases at the sites. Operational conditions were emulated by pumping abrasive materials through the pipeline at full capacity and collecting measurements from the SRNL coupons. Initial testing was conducted by circulating water followed by a mixture of sand and water for uniform wear rates. Thickness changes and mass loss were measured, and the surface exposed to the process fluid was observed for visual evidence of characteristic wear patterns resulting from an erosive environment. The test matrix consisted of several hours of operation and using fine, medium and coarse grain sands. Test results were analyzed to check for all three parameters (mass loss/gain, thickness change and surface patters) which could provide estimates of the remaining useful life of the pipe sections. Thus, the research presents the findings of this evaluation of the SRNL coupons as a means for conducting high accuracy, cost effective in-situ analysis of pipeline and waste transfer systems for erosion and corrosion susceptibility. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗