SEARCH · Engineering Papers

Results for “fluid instabilities”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Contribution to: UNM-UNR-VT Project Research Performance Progress Report 'Seeding and Evolution of Magnetohydrodynamic Instabilities of a Metal Surface Driven by Intense Current' Year 4 (no-cost extension for FY 2022): Oct 2021 - Sept 2022

The multiphysics code Ares has been used to perform resistive MHD fluid simulations of the METI-II experiments (described in Section 1.1.2). The simulations include both 1D Lagrangian simulations and axisymmetric 2D ALE simulations. The results of these simulations are summarized here and a full manuscript is nearing completion for submission.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Increasing CO 2 Capture Rate in Liquid-Solvent Direct-Air Carbon Capture via Additive Manufacturing

Carbon capture is essential to meeting climate change mitigation goals. One approach currently being commercialized utilizes liquid-based solvents to capture CO 2 directly from the atmosphere but is limited by slow absorption of CO 2 into the liquid. Improved air/solvent liquid mixing increases CO 2 absorption rate, and this increased CO 2 absorption efficiency allows for smaller carbon capture systems with lower capital costs and better economic viability. In this project, we study the use of passive micromixers fabricated by metal additive manufacturing. The micromixer’s small-scale surface geometric features perturb and mix the liquid film to enhance mass transfer and CO 2 absorption. In this project, we evaluated this hypothesis through computational and experimental studies. Computational investigations focused on developing capabilities to simulate thin film (~ 100μm) fluid flow on rough surfaces. Such thin films are in a surface-tension dominated regime and simulations in this regime are prone to instabilities. Improvements to the Nalu code completed in this project resulted in a 10x timestep stability improvement for these problems.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Planetary waves can activate resonant drag instabilities in 3D dusty gaseous discs

Resonant drag instabilities (RDIs) in protoplanetary discs are driven by the aerodynamic back-reaction of dust on gas and occur when the relative dust-gas motion resonate with a wave mode intrinsic to the gas fluid. Axisymmetric models indicate that the RDI generates filamentary perturbations, leading to grain clumping and planetesimal formation. Motivated by these findings, we investigate the dust–gas interaction in a non-axisymmetric inviscid protoplanetary disc with an embedded low-mass planet ($M_{\mathrm{p}}\in [0.3, 3] \, \rm M_{\oplus }$, here $\, \rm M_{\oplus }$ is the Earth mass). We conduct global 3D high-resolution two-fluid simulations, with the dust being parametrized by the Stokes number $\mathrm{St}\in [0.01,0.5]$. We find that planetary waves (PWs; also known as Rossby waves), which propagate along the downstream separatrices of the horseshoe region, resonate with the streaming motion and trigger the RDI. The consequent development of a global-scale filamentary dust distribution does not sensitively depend on the Stokes number, nor does it depend on the fast dust settling that takes place in an inviscid disc. The rapid onset of this instability, which is comparable to the dynamical orbital time-scale, suppresses the formation of asymmetric structures in the dust in the vicinity of the planet (such as dust voids and filaments). Additionally, we find that the dust feedback enables buoyancy resonances in an otherwise non-buoyant (globally isothermal) disc. Therefore, our results provide the first numerical evidence of RDIs generation driven by planetary waves.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Theoretical study of various nonlinear phenomena in plasma systems and scaling of magneto-inertial-fusion targets

Plasma physics is an exciting field of study with a wide variety of nonlinear processes that come into play. Examples of such processes include the interaction of small-scale turbulence with large-scale plasma structures and the nonlinear saturation of plasma instabilities, for example those of magneto-hydrodynamical nature. During this Truman LDRD project, I studied a collection of nonlinear problems that are of interest to the field of plasma physics. This LDRD report summarizes four main research accomplishments. First, a new statistical model for describing inhomogeneous drift-wave turbulence inter- acting with zonal flows was developed. This new model includes the effects of nonlinear wave-wave collisions, which are expected to change the spectrum of the underlying DW turbulence and therefore the generation of zonal flows. Second, a new mathematical formalism was proposed to systematically apply the non- linear WKB approximation to general field theories, including those often used in fluid dynamics. This formalism represents an interesting tool for studying physical systems that show an explicit scale separation. Third, a weakly nonlinear model was developed to describe the magneto-Rayleigh-Taylor instability. This instability is of paramount importance to understand as it can reduce the performance of magnetic-inertial-fusion (MIF) platforms. The developed models captures the effects of harmonic generation and saturation of the linear growth of the instability. Finally, a framework was proposed for scaling magneto-inertial fusion (MIF) targets to larger pulsed-power drivers. From this framework, a set of scaling rules were derived that conserve the physical regimes of MIF systems when scaling up in peak current. By doing so, deleterious nonlinear processes that affect MIF performance may be kept at bay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Frictional Behavior Downdip Along the Subduction Megathrust: Insights From Laboratory Experiments on Exhumed Samples at In Situ Conditions

Subduction megathrusts exhibit a spectrum of slip modes, including catastrophic earthquakes. Although the mechanical and frictional properties of materials sampled from subduction zones have been studied extensively, few datasets have been collected for compositions and at pressure and temperature conditions representative of those in situ. The Nankai subduction zone in southwest Japan is a well-studied margin, and abundant data provide an opportunity to advance our understanding of fault and earthquake physics. Here, we use samples exhumed in the Shimanto and Sanbagawa Belts on Shikoku Island of southwest Japan that represent analogs for materials along the present-day megathrust at depths of ~5–>25 km, and we shear these at their peak in situ pressure-temperature (P-T) conditions. Effective normal stresses range from 28 to 192 MPa, and temperatures from 105°C to 470°C. We used pore fluid pressures of 45–240 MPa, corresponding to fluid overpressure ratios λ of 0.65 and 0.90. Slip velocities of 0.1–100 μm/s were used, in order to focus on the nucleation of instability and earthquakes. We found predominantly velocity-strengthening (inherently stable) behavior under all conditions for λ = 0.65. For λ = 0.90, velocity-weakening behavior was observed at 350°C, with velocity-strengthening behavior at lower and higher temperatures. The rate/state frictional stability parameter (a-b) increases with slip velocity at temperatures up to ~200°C and remains constant or decreases with slip velocity at higher temperatures. Altogether, our results demonstrate the potentially important roles of both temperature and slip velocity in controlling the distribution of stress and frictional rheology along subduction thrusts.

58 GEOSCIENCES↗

DOE OSTI2020

Kinetic simulations of anomalous resistivity in high-temperature current carrying plasmas

Anomalous resistivity caused by lower hybrid drift instabilities (LHDIs) has been extensively studied in the literature, and has been invoked to explain the rates of magnetic diffusion and plasma profile evolution in the low-density plasma periphery of theta and screw pinch plasmas. In this article, we present a fully kinetic simulation study of LHDIs with a realistic mass ratio using the particle-in-cell code Chicago over a large range of drift speeds, 0.5 < v$^{2}_{di}$/v$^{2}_{ti}$ < 14.5, plasma β, β ≥ 2×10 –5 , and temperature ratios, 0.1 < T e /T i < 10. The resistivity quantified from the simulations is compared with the analytic estimates from the literature and a generalized resistivity expression is presented, which can be applied over the entire range of plasma environments that were simulated in this study. This expression uses the amplitude of the LHDI fluctuations resulting from instability saturation by electron resonance broadening rather than plateau formation, current relaxation, or ion trapping. The generalized resistivity expression is appropriate for full-scale fluid simulations that encounter a large range of plasma conditions, and has better agreement with the resistivity quantified from kinetic simulations than previous expressions in the literature, particularly at low plasma β. Anomalous resistivity has also been considered as an explanation of the rapid magnetic field dissipation during magnetic reconnection events at large plasma β. At a large β, the magnetic field fluctuations of LHDIs can dominate the electrostatic fluctuations and this article investigates enhanced collisionality due to the nonlinear magnetic force, δJ×δB/c.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Simulating Fracture Sealing by Granular LCM Particles in Geothermal Drilling

Lost circulation occurs when the returned fluid is less than what is pumped into the well due to loss of fluid to pores or fractures. A lost-circulation event is a common occurrence in a geothermal well. Typical geothermal reservoirs are often under-pressured and have larger fracture apertures. A severe lost-circulation event is costly and may lead to stuck pipe, well instability, and well abandonment. One typical treatment is adding lost-circulation materials (LCMs) to seal fractures. Conventional LCMs fail to properly seal fractures because their mechanical limit is exceeded at elevated temperatures. In this paper, parametric studies in numerical simulations are conducted to better understand different thermal effects on the sealing mechanisms of LCMs. The computational fluid dynamics (CFDs) and the discrete element method (DEM) are coupled to accurately capture the true physics of sealing by granular materials. Due to computational limits, the traditional Eulerian–Eulerian approach treats solid particles as a group of continuum matter. With the advance of modern computational power, particle bridging is achievable with DEM to track individual particles by modeling their interactive forces between each other. Particle–fluid interactions can be modeled by coupling CFD algorithms. Fracture sealing capability is investigated by studying the effect of four individual properties including fluid viscosity, particle size, friction coefficient, and Young’s modulus. It is found that thermally degraded properties lead to inefficient fracture sealing.

15 GEOTHERMAL ENERGY↗

DOE OSTI2020

A novel, magnetically driven convergent Richtmyer–Meshkov platform

Here, we introduce a novel experimental platform for the study of the Richtmyer–Meshkov instability in a cylindrically converging geometry using a magnetically driven cylindrical piston. Magnetically driven solid liner implosions are used to launch a shock into a liquid deuterium working fluid and, ultimately, into an on-axis rod with a pre-imposed perturbation. The shock front trajectory is tracked through the working fluid and up to the point of impacting the rod through the use of on axis photonic Doppler velocimetry. This configuration allows for precise characterization of the shock state as it impacts the perturbed rod interface. Monochromatic x-ray radiography is used to measure the post-shock interface evolution and rod density profile. The ALEGRA MHD model is used to simulate the dynamics of the experiment in one dimension. We show that late in time the perturbation growth becomes non-linear as evidenced by the observation of high-order harmonics, up to n = 5. Two dimensional simulations performed using a combination of the GORGON MHD code and the xRAGE radiation hydrodynamics code suggest that the late time non-linear growth is modified by convergence effects as the bubbles and spikes experience differences in the pressure of the background flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Townsend to glow discharge transition for a nanosecond pulse plasma in helium: space charge formation and resulting electric field dynamics

Stark polarization spectroscopy is used to investigate the temporal evolution of the electric field distribution in the cathode region of a nanosecond pulsed discharge in helium at 120 Torr. The measurements are performed on the He I transition at 492.19 nm, during the early stages of the discharge formation. The experimental results are compared with the predictions of a 1D fluid model. Time-resolved ICCD images show that the discharge develops as a diffuse, cathode-directed ionization wave with a Townsend-like feature before transitioning into a glow-like structure. Near anode instabilities characterized by filament formation were observed near the high voltage electrode. Within 30 ns, a reduction of the sheath thickness to about 250 μm is observed, coinciding with a gradual increase of the discharge current and proportional increase in electric field at the cathode. Additionally, the cathode electric field corresponding to this sheath with a thickness of 250 μm is about 40 kV cm -1 . A subsequent steep increase of the discharge current leads to a further reduction of the sheath width. The electric field evolution as obtained by the fluid model is in excellent agreement with the measurements and shows that an enhanced ionization near the cathode is causing the space charge formation responsible for the increase in electric field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Numerical modeling of pedestal stability and broadband turbulence of wide-pedestal QH-mode plasmas on DIII-D

The wide-pedestal quiescent high confinement mode discovered on DIII-D in recent years is a stationary and quiescent H-mode (QH-mode) with the pedestal width exceeding EPED prediction by at least 25%. Its characteristics, such as low rotation, high energy confinement and edge localized mode-free operation, make it an attractive operation mode for future reactors. Linear and nonlinear simulations using BOUT++ reduced two fluid MHD models and awere carried out to investigate the bursty broadband turbulence often observed in the edge of wide-pedestal QH-mode plasmas. Two kinds of MHD-scale instabilities in different spatial locations within the pedestal were found in the simulations: one mild peeling–ballooning (PB) mode γ PB < 0.04ω A ) located near the minimum in E r</:sub> well propagating in ion diamagnetic drift direction; and one drift-Alfvén wave locates at smaller radius compared to E r</:sub> well propagating in the electron diamagnetic drift direction and unstable only when the parallel electron dynamics is included in the simulation. The coupling between drift wave and shear Alfvén wave provides a possible cause of the experimentally observed local profile flattening in the upper-pedestal. The rotation direction, mode location, as well as the wavenumber of these two modes from BOUT++ simulations agree reasonably well with the experimental measurements, while the lack of quantitative agreement is likely due to the lack of trapped electron physics in current fluid model. This work presents improved physics understanding of the pedestal stability and turbulence dynamics for wide-pedestal QH-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

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↗

DOE OSTI2021

Magnetic-Confinement Fusion—Plasma Theory: Tokamak Magnetohydrodynamic Equilibrium and Stability

Magnetohydrodynamics (MHD) provides a useful model to describe the crucial plasma macroscopic equilibrium and stability behaviors in toroidal tokamak devices by considering the plasma as a conducting fluid interacting with a surrounding confining electromagnetic field. MHD is the most basic plasma model, incorporating most large-scale phenomena, including plasma equilibrium and all major instabilities. MHD equations are obtained by taking moments of the Boltzmann equations for different plasma species. They provide a set of comprehensive physics constrains to compute and optimize the equilibrium plasma shape and pressure and current profiles that are critical to its stability and performance. In the ideal case, the equations have special properties that lead to efficient numerical calculation schemes, the most important of which is the ideal MHD energy principle for linear stability against small departures from equilibrium. In a tokamak plasma, equilibrium pressure is mostly destabilizing for MHD modes, whereas equilibrium current is also often a major driving force. Plasma resistivity creates new freedom for a MHD instability to grow, but there are also cases where the plasma resistivity plays a stabilizing role. Equilibrium toroidal flow and/or flow shear can affect MHD instabilities. Principal MHD instabilities include the internal kink mode, sawtooth, fishbone, external kink, resistive wall mode, resistive interchange, tearing and neoclassical tearing modes (NTMs), locked modes, toroidal Alfven eigenmodes (TAEs), and edge localized modes (ELMs). Fast-growing MHD instabilities can lead to an abrupt plasma disruption and termination that can potentially damage the device plasma facing components (PFCs) and in-vessel structures. Furthermore, an important MHD application is to develop robust techniques to mitigate and control MHD instabilities.

Dispersive pellet injection↗

DOE OSTI2020

Effect of the numerical discretization scheme in Shock-Driven turbulent mixing simulations

In this work, we evaluate the effects of distinct numerical strategies in simulations of shock-driven turbulent mixing. An air-SF 6 -air gas curtain subjected to Mach 1.2 shock-waves is computed with Implicit Large-Eddy Simulation based on three numerical schemes: directional-split with Harten-Lax-van Leer (HLL) solver; directional-unsplit using HLL-Contact (HLLC) solver; and directional-unsplit with HLLC solver and a Low Mach number Correction (LMC). The results illustrate the importance of the numerical strategy to the accuracy of the predictions. Whereas both split and unsplit schemes result in a similar spatial development of the initial shock-driven instability, only the unsplit schemes can predict the turbulent mixing transition after reshock observed by the laboratory experiments. Such feature increases the mixing rate of the two fluids, this being particularly pronounced when the LMC is active due to i) the reduced flow Mach number, and ii) the larger effective Reynolds number. Since the selected mixing problem is driven by the deposition of vorticity at the fluids’ interface, the resultant flow physics is analyzed by investigating the contribution of distinct inviscid mechanisms to the production of vorticity: baroclinicity, stretching, and dilatation. As expected, it is observed that the production of vorticity is initially dominated by the baroclinicity mechanism. Yet, the relevance of the remainder mechanisms is enhanced after the reshock and may even surpass the baroclinicity term.

42 ENGINEERING↗

DOE OSTI2020

A machine learning study on spinodal clumping in heavy ion collisions

Possible observables of baryon number clustering due to the instabilities occurring at a first order QCD phase transition are discussed. The dynamical formation of baryon clusters at a QCD phase transition can be described by numerical fluid dynamics, augmented with a gradient term and an equation of state with a mechanically unstable region. It is shown that the dynamical description of this phase transition, in nuclear collisions, will lead to the formation of dense baryon clusters at the phase boundary. State-of-the-art machine learning methods find that the coordinate space clumping leaves characteristic imprints on the spatial net density distribution in almost every event. On the other hand the momentum distributions do not show any clear event-by-event features. Lastly, it is shown that the 'third order' cumulant, the skewness, shows a peak at the beam energy where the system, created in the heavy ion collision, reaches the deconfinement phase transition.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Resistive hose modes in tokamak runaway electron beams

Beams of energetic runaway electrons are generated during disruptions in tokamaks, and fluid models are used to study their effects on macroscale dynamics. Linear computations of a massless, runaway electron beam coupled to MHD plasma show that resistive hose instabilities grow faster than tearing modes at large resistivity. Eigenvalue results with reduced models of the resistive hose instability are compared with results from the full MHD and beam system, showing that the resistive hose decouples from any plasma response. An estimate of plasma temperature at which growth of the resistive hose dominates tearing for post-disruption DIII-D plasma parameters is in a physically relevant regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Coupling of Laminar-Turbulent Transition with RANS Computational Fluid Dynamics

This project combines several new concepts to create a boundary layer transition prediction capability that is suitable for analyzing modern hypersonic flight vehicles. The first new concept is the use of ''optimization'' methods to detect the hydrodynamic instabilities that cause boundary layer transition; the use of this method removes the need for many limiting assumptions of other methods and enables quantification of the interactions between boundary layer instabilities and the flow field imperfections that generate them. The second new concept is the execution of transition analysis within a conventional hypersonics CFD code, using the same mesh and numerical schemes for the transition analysis and the laminar flow simulation. This feature enables rapid execution of transition analysis with less user oversight required and no interpolation steps needed.

42 ENGINEERING↗

DOE OSTI2022

Modeling of thermal-ion-driven internal kink in DIII-D high- T i plasmas

Toroidal modeling utilizing the non-perturbative magnetohydrodynamic-kinetic hybrid code MARS-K (Liu Y.Q. et al 2008 Phys. Plasmas 15 112503), to a DIII-D high-performance discharge with the core thermal ion temperature of T i ~ 13 keV, finds multiple branches of unstable internal kink (IK) modes triggered by the bounce resonance of trapped thermal ions. The range of the computed mode frequency aligns well with the bounce frequency of trapped thermal ions, confirming the nature of resonance drive for the instability. The computed mode frequency is close to that observed in experiments. Plasma toroidal flow, in particular the fast flow ($ > $10% of Alfvén speed) as in the considered DIII-D discharge, strongly stabilizes the fluid IK branch but not all of the thermal-ion-driven internal kink (TIK) branches, leaving one dominant unstable TIK branch under the experimental conditions. Plasma resistivity also significantly affects stability of the TIKs (but not much on the IK which is ideally strongly unstable) in combination with the drift kinetic drive by thermal ions. High resistivity can drive additional unstable TIK branches. Scanning the on-axis safety factor over a large range confirms robustness of the TIK triggering by thermal ions. The computed TIKs contain both the IK and the tearing components as the dominant contributions to the eigenmode structure. Furthermore, these findings, together with the experimental results reported in reference (Du X.D. et al 2021 Phys. Rev. Lett. 127 025001), strongly suggest the possibility of thermal-ion driven instabilities in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Turbulence theories and statistical closure approaches

When discussing research in physics and in science more generally, it is common to ascribe equal importance to the three components of the scientific trinity: theoretical, experimental, and computational studies. This review will explore the future of modern turbulence theory by tracing its history, which began in earnest with Kolmogorov’s 1941 analysis of turbulence cascade and inertial range [A.N. Kolmogorov, Dokl. Akad. Nauk SSSR, 30, 299, (1941); 32, 19, (1941)]. The 80th Anniversary of Kolmogorov’s landmark study is a welcome opportunity to survey the achievements and evaluate the future of the theoretical approach of turbulence research. Over the years, turbulence theories have been critically important in laying the foundation of our understanding of the nature of turbulent flows. In particular, the Direct Interaction Approximation (DIA) [R.H. Kraichnan, J. Fluid Mech., 5, 497 (1959)] and its subsequent development, known as the statistical closure approach, can be identified as perhaps the most profound single advancement. The remarkable success of the statistical closure has furnished a platform to study such essential concepts as the energy transfer process and interacting scales, and the roles of the straining and sweeping motions. More recently, the quasi-Lagrangian formulation of V. L’vov & I. Procaccia and Kraichnan’s solvable passive scalar model provided powerful ways to explore another fundamental aspect of turbulent flows, the phenomena of intermittency, and the associated anomalous scaling exponents. In the meantime, the theory of fluid equilibria has been developed to describe the large-scale structures that can emerge from turbulent cascades of two-dimensional and geophysical flows at a later time. And yet, despite all these successes, analytical treatments suffer from mathematical complexities. As a result, the utility of theoretical approaches has been limited to relatively idealized flows. On the other hand, in recent decades, computational abilities and experimental facilities have reached an unprecedented scale. Looking beyond the horizon, the imminent deployment of exascale supercomputers will generate complete datasets of the entire flow field of key benchmark flows, allowing researchers to extract additional measurements concerning fully developed, complex turbulent flow fields far beyond those available from the statistical closure theories. Some other developments that could potentially influence the future course of turbulence theories include the advancement of machine learning, artificial intelligence, and data science; likely disruptions arising from the advent of quantum computation; and the increasingly prominent role of turbulence research in providing more accurate climate scientific data. Finally, turbulence theorists can leverage these developments by asking the right questions and developing advanced, sophisticated frameworks that will be able to predict and correlate vast amounts of data from the other two components of the trinity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗