SEARCH · Engineering Papers

Results for “Flow Conditions”

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 721 records · Page 40

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Pore architecture controls on mineral reactivity

Mineral dissolution rates measured in natural environments are much slower than those measured in laboratory settings. This project tested the hypothesis that the way fluid flows through rocks in natural systems creates areas where mineral dissolution is fast and areas where mineral dissolution is slow. This hypothesis was tested with a combination of laboratory experiments and numerical simulation. We demonstrated a separation of fluid flow pathways and rates of mineral dissolution in laboratory experiments for the first time using an experimental approach where we created synthetic rocks that have different ratios of connected and dead-end pathways for fluid flow. In laboratory experiments with higher proportions of dead-end pathways, the mineral dissolution rates were slower. We found that where fluid flows through connected pathways the continuous refreshing of fluid at the mineral surface creates conditions where dissolution is fast. In contrast, where fluid either flows slowly through poorly connected pathways or is stagnant in dead-end pathways, mineral dissolution is slow. The results from this project suggest that the overall slowing of rates of mineral dissolution is important when the proportion of dead-end pathways is greater than ~40%. This project informs our understanding of the way that fluids react with rocks in carbon dioxide sequestration and enhanced geothermal projects where fluids are purposefully injected into rocks for energy applications.

58 GEOSCIENCES↗

DOE OSTI2026

Time-resolved spray characterization via unified optical flow and binarization technique

This work leverages an unsupervised machine learning and advanced image processing techniques to characterize the breakup of fuel sprays in a small-scale combustor under reacting conditions, providing valuable insights into near-nozzle flow phenomenology. The proposed methodology integrates an improved optical flow model on a convolutional neural network to extract flow vectors with a binarization technique to assess droplets’ size and shape across the region of interest. The velocimetry approach demonstrates superior performance compared to a state-of-the-art optical flow model when applied to high-speed X-ray phase contrast spray images, achieving more accurate and reliable flow predictions. Moreover, breakup processes are quantified by breakup length and sphericity in accordance with velocity estimations, allowing a more complete characterization of the flow. This study establishes a robust methodology for analyzing spray morphology and primary breakup in compact combustors, contributing valuable means of understanding and optimizing fuel spray behavior in advanced combustion systems.

42 ENGINEERING↗

DOE OSTI2022

Mass Transport in Membrane Systems: Flow Regime Identification by Fourier Analysis

The numerical calculation of local mass distributions in membrane systems by computational fluid dynamics (CFD) offers indispensable benefits. However, the concept to calculate such distributions in response to separate variations of operation conditions (OCs) makes it difficult to address overall, flow-physics-related questions, which require the consideration of the collective interaction of OCs. It is shown that such understanding-related relationships can be obtained by the analytical solution of the advection–diffusion equation considered. A Fourier series model (FSM) is presented, which provides exact solutions of an advection–diffusion equation for a wide range of OCs. On this basis, a new zeroth-order model is developed, which is very simple and as accurate as the complete FSM for all conditions of practical relevance. Advection-dominated blocked and diffusion-dominated unblocked flow regimes are identified (depending on a Péclet number which compares the flow geometry with a length scale imposed by the flow), which implies relevant requirements for the use of lab results for pilot- and full-scale applications. Analyses reveal the equivalence of variations of OCs, which offers a variety of options to accomplish desired flow regime changes.

Heinz, Stefan (ORCID:0000000248712416)↗

DOE OSTI2021

Hydraulic Response to Thermal Stimulation Efforts at Raft River Based on Stepped Rate Injection Testing

The injection well stimulation project at the Raft River geothermal field tests the effect of long-term cold water injection and high pressure injection on well injectivity, improvements to which could reduce operating costs. The primary data for analysis and interpretation of the injection test are step-rate flow tests run before each new phase of the injection. These tests were analyzed using a combination of standard pump-test analytical solution methods and methods developed expressly for the observed conditions. The stepped rate injection tests, combined with long-term flow and pressure response data suggest that the well is located within a fractured formation of low transmissivity but high storativity. These calculated parameters appeared to increase with pressure during the first injection test and the higher values were reproduced during the second stepped rate test. Calculated transmissivity and storativity are on the order of 4E-5 m cm 2 and 1E-4 m Pa -1 , respectively. The apparent pressure dependence of fitted hydraulic parameters may reflect near-well fracture compliance that increased the effective radius of the wellbore during the first test. While the type curve fit analysis also suggests that the reservoir behaves as a uniformly fractured reservoir with a radial flow regime, the hydraulic parameters indicate that condition may exist only a very limited distance (<10 m) from the well. Longer-term pressure response suggests that flow in the system effectively reaches steady state in a period of less than a day, which may reflect pressure stabilization resulting from pressure-dependent permeability or a region of much higher permeability located with a few meters of the well. The transmissivity estimates obtained from this analysis, converted to approximate fracture density and aperture, provide useful constraints on the distance to which the thermal front may migrate from the well during the cold water injection phase of the stimulation project. We estimate that the cooling front will migrate less than a tenth of a kilometer over an approximately one-year injection period. Here, the effects of that cooling, however, may be substantial, because increases in permeability have maximum effect nearest the well.

cold water injection↗

DOE OSTI2021

Particle charge distributions in the effluent of a flow-through atmospheric pressure low temperature plasma

Atmospheric pressure low-temperature plasmas are often utilized to perform particle synthesis, treatment, and removal. It is well-known that dust particles are highly negatively charged in these plasmas; however, little is known about dust particle charging behavior as particles leave the plasma volume and pass through the spatial afterglow region. In this work, monodisperse particles of various sizes and work functions were introduced into an atmospheric pressure radiofrequency capacitively coupled flow-through plasma. Dust particle electrical mobility distributions downstream of the flow-through plasma were measured utilizing a differential mobility analyzer in conjunction with a condensation particle counter at various gas flow velocities. Charge distributions were determined from the measured electrical mobility distributions. Experiments confirm that particles become less negatively charged, and even net-positively charged after leaving the plasma volume, with a distribution that follows a shifted Boltzmann charge distribution. Additionally, particle charge in the effluent of the flow-through plasma is negligibly dependent on work function but highly size and flow velocity dependent. Larger particles were shown to have a higher magnitude of charge under all studied conditions; however, particle polarity was switchable by varying gas flow velocity. The charging dynamics were simulated utilizing a constant number Monte Carlo model that accounts for electron temperature decay and the transition from ambipolar to free diffusion of electrons and ions in the spatial afterglow. Simulation results also suggest that, at the same flow velocity, larger particles obtain a greater magnitude of charge, negative or positive. Furthermore, the decrease in electron mobility and the difference between ion and electron convective loss rates create an ion-rich region in the plasma effluent that promotes ion–particle collisions and drives particle charge removal and even reversal of polarity. Larger particles more favorably collide with energetic species in these environments, which results in higher charge states.

atmospheric pressure low temperature plasma↗

DOE OSTI2024

Complementary Subsurface Characterization Methods to Develop a Geologic Model for the EGS Collab Experiment, Sanford Underground Research Facility

The EGS (Enhanced Geothermal Systems) Collab project was performed within the Sanford Underground Research Facility (SURF) with a goal of understanding processes and evaluation of models related to hydraulic stimulation of rock at depth. The present work deals with the development of Testbed 2 where experiments were conducted at a depth of 1.25 km and were located within a well-characterized testbed in a metamorphic, amphibolite host rock. A total of eleven boreholes varying in length between 10.6 m and 81.2 m were continuously cored to develop the testbed. In addition to the continuous coring of the amphibolite host rock, geophysical instrumentation supporting electrical resistivity tomography (ERT), microearthquake (MEQ) detection, and optical fiber providing distributed temperature sensing (DTS) and distributed acoustic sensing (DAS) were installed in the monitoring boreholes, all of which produced a comprehensive complementary suite of characterization and monitoring technologies. Pre-stimulation characterization of the groundwater conditions identified only two hydraulically significant fractures, neither of which transects the central portion of the testbed. Flow and pressure monitoring indicated low preexisting pore pressure conditions likely affected by the mine openings.

15 GEOTHERMAL ENERGY↗

DOE OSTI2021

Influence of Agricultural Managed Aquifer Recharge (AgMAR) and Stratigraphic Heterogeneities on Nitrate Reduction in the Deep Subsurface

Agricultural managed aquifer recharge (AgMAR) is a strategy whereby surface water is used to intentionally flood croplands to recharge underlying aquifers. However, nitrate (NO 3 - ) contamination in agriculturally intensive regions poses a threat to groundwater resources under AgMAR. We use a reactive transport model to understand the effects of AgMAR management strategies (i.e., by varying the frequency, duration between flooding events, and amount of water) on NO 3 - leaching to groundwater under different stratigraphic configurations and antecedent moisture conditions. We examine the potential of denitrification and nitrogen retention in deep vadose zone sediments (~15 m) using variable AgMAR application rates on two-dimensional representations of differently textured soils, soils with discontinuous bands/channels, and with preferential flow paths characteristic of agricultural fields. Simulations indicate finer textured sediments, alone or embedded within/adjacent to high flow regions, are important reducing zones providing conditions needed for denitrification. Simulation results suggest that applying water all-at-once rather than in increments transports higher concentrations of NO 3 - deeper into the profile, which may exacerbate groundwater quality. This transport into deeper depths can be aggravated by wetter antecedent soil moisture conditions. However, applying water all-at-once also increases denitrification within the vadose zone by promoting anoxic conditions. We conclude that AgMAR can be designed to enhance denitrification in the subsurface and reduce NO 3 - leaching to groundwater, while specifically accounting for lithologic heterogeneity, antecedent soil moisture conditions, and depth to the water table. Our findings are potentially relevant to other systems that experience flooding inundation such as floodplains and dedicated recharge basins.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Flow and arrest in stressed granular materials

Flowing granular materials often abruptly arrest if not driven by sufficient applied stresses. Such abrupt cessation of motion can be economically expensive in industrial materials handling and processing, and is significantly consequential in intermittent geophysical phenomena such as landslides and earthquakes. Using discrete element simulations, we calculate states of steady flow and arrest for granular materials under the conditions of constant applied pressure and shear stress, which are also most relevant in practice. Here the material can dilate or compact, and flow or arrest, in response to the applied stress. Our simulations highlight that under external stress, the intrinsic response of granular materials is characterized by uniquely-defined steady states of flow or arrest, which are highly sensitive to interparticle friction. While the flowing states can be equivalently characterized by volume fraction, coordination number or internal stress ratio, to characterize the states of shear arrest, one needs to also consider the structural anisotropy in the contact network. Here, we highlight the role of dilation in the flow-arrest transition, and discuss our findings in the context of rheological transitions in granular materials.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Hydrodynamic characterization of the redox chemistry of crown-encapsulated uranyl complexes

The redox properties of actinide-containing species strongly influence their reactivity, speciation, and interfacial behavior, but the experimental quantification of the electrochemical characteristics of molecular actinide complexes in nonaqueous media has not received the attention it deserves. Here, results from hydrodynamic methods and electrochemical simulations of U(VI)/U(V) redox are reported, including quantification of heterogeneous electron-transfer kinetics and estimation of chemical reversibility of U(VI)/U(V) interconversion at electrodes in acetonitrile-based electrolyte. The complexes investigated are recently reported U(VI) and U(V) complexes in which the uranyl ion (UO 2 n+ ) is encapsulated in a macrocyclic 18-crown-6-like moiety templated by a Pt(II) center. These complexes feature the most positive value U VI /U V reduction potential yet reported and are thus particularly relevant to study of facile U(V) generation from U(VI) precursors as well as uranium electroanalysis. Rotating disk electrode (RDE) studies have been used to quantify the diffusion coefficients of the U(VI) and U(V) complexes, and standard heterogeneous electron transfer rate constants ( k 0 ) for the redox have been determined using a conventional Koutecký-Levich analysis. Rotating ring-disk electrode (RRDE) studies have been used to directly interrogate the chemical reversibility of U(VI)-U(V) interconversion, confirming that reduction of the U(VI) complex at an Au disk is associated with formation of the U(V) analogue that can be readily re-oxidized at a Pt ring under hydrodynamic (rotating) conditions. Because measurements of the type reported here are generally associated with current flows that are larger than those found in corresponding quiescent (unstirred) conditions, our findings suggest that hydrodynamic methods could be advantageous for design of electroanalytical approaches to detection of actinide species and study of their redox properties.

actinides↗

DOE OSTI2022

The influence of intake flow and coolant temperature on gasoline spray morphology during early-injection DISI engine operation

Multi-hole gasoline injectors operating at conditions spanning throttled early-intake stroke operation produce spray plumes that either remained separated or merge and collapse due to flash boiling. Flash boiling occurs due to the sudden expansion of gas bubbles in the liquid fuel at high fuel temperature and low ambient pressure. This study records high-speed images of spray-morphology changes due to in-cylinder flow, thereby revealing operating conditions that do and do not affect the self-induced morphology observed in quiescent vessels. Specifically, in a central-injection, four-valve, high-tumble engine, where the thermodynamic state and in-cylinder cross flow are dynamic. Motivated by cold start and hot restart operation, the fuel pressure, coolant temperature, in-cylinder air pressure, and engine rpm were systematically varied over relevant operating conditions, which bracketed the range from non- to flash-boiling sprays. The results reveal the operating conditions at which the in-cylinder cross flow disrupts the spray morphology as well as the extent of the disruption. At 650 rpm, the spray morphology was similar to that observed in quiescent vessels at nominally equivalent fuel temperature and in-cylinder pressure, indicating that the spray’s self-induced entrainment flow dominated the in-cylinder flow. However, for fuel temperature and ambient pressure near the transition between non- and flash-boiling, the intake cross flow at higher engine speed (1950 rpm) significantly disrupted the spray morphology. The high cross-flow velocity appears to induce plume merging and collapse, whereas none was evident at low rpm (650 rpm). This study led to the postulate that the spray merging and collapse are governed by the rate of atomization near the nozzle exit, presumed to be controlled by either or both aerodynamic atomization and flash-boiling intensity. It would then follow that spray modeling in CFD requires atomization models that blend the effects of both physical processes.

42 ENGINEERING↗

DOE OSTI2022

Design and execution of a verification, validation, and uncertainty quantification plan for a numerical model of left ventricular flow after LVAD implantation

Left ventricular assist devices (LVADs) are implantable pumps that act as a life support therapy for patients with severe heart failure. Despite improving the survival rate, LVAD therapy can carry major complications. Particularly, the flow distortion introduced by the LVAD in the left ventricle (LV) may induce thrombus formation. While previous works have used numerical models to study the impact of multiple variables in the intra-LV stagnation regions, a comprehensive validation analysis has never been executed. The main goal of this work is to present a model of the LV-LVAD system and to design and follow a verification, validation and uncertainty quantification (VVUQ) plan based on the ASME V&V40 and V&V20 standards to ensure credible predictions. The experiment used to validate the simulation is the SDSU cardiac simulator, a bench mock-up of the cardiovascular system that allows mimicking multiple operation conditions for the heart-LVAD system. The numerical model is based on Alya, the BSC’s in-house platform for numerical modelling. Alya solves the Navier-Stokes equation with an Arbitrary Lagrangian-Eulerian (ALE) formulation in a deformable ventricle and includes pressure-driven valves, a 0D Windkessel model for the arterial output and a LVAD boundary condition modeled through a dynamic pressure-flow performance curve. The designed VVUQ plan involves: (a) a risk analysis and the associated credibility goals; (b) a verification stage to ensure correctness in the numerical solution procedure; (c) a sensitivity analysis to quantify the impact of the inputs on the four quantities of interest (QoIs) (average aortic root flow $Q^{avg}_{Ao}$, maximum aortic root flow $Q^{max}_{Ao}$, average LVAD flow $Q^{avg}_{VAD}$, and maximum LVAD flow $Q^{max}_{VAD}$; (d) an uncertainty quantification using six validation experiments that include extreme operating conditions. Numerical code verification tests ensured correctness of the solution procedure and numerical calculation verification showed a grid convergence index (GCI)95% <3.3%. The total Sobol indices obtained during the sensitivity analysis demonstrated that the ejection fraction, the heart rate, and the pump performance curve coefficients are the most impactful inputs for the analysed QoIs. The Minkowski norm is used as validation metric for the uncertainty quantification. It shows that the midpoint cases have more accurate results when compared to the extreme cases. The total computational cost of the simulations was above 100 [core-years] executed in around three weeks time span in Marenostrum IV supercomputer. This work details a novel numerical model for the LV-LVAD system, that is supported by the design and execution of a VVUQ plan created following recognised international standards. We present a methodology demonstrating that stringent VVUQ according to ASME standards is feasible but computationally expensive.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Riverscape heterogeneity in estimated Chinook Salmon emergence phenology and implications for size and growth

Abstract Many salmonid‐bearing rivers exhibit thermal and hydrologic heterogeneity at multiple spatial and temporal scales, but how this translates into spatiotemporal patterns of fry emergence is poorly understood. Understanding this variability is important because emergence timing determines the biophysical conditions fish first experience (e.g., temperature, flow, and food supply), thereby influencing growth opportunities and survival during this critical life stage. We predicted spring Chinook Salmon ( Oncorhynchus tshawytscha ) emergence phenology across four northeastern Oregon subbasins over 5–9 years using empirical spawning and temperature data. We then related interannual emergence timing estimates to juvenile salmon size and growth rates at consistent sampling locations. There were clear longitudinal patterns of predicted emergence timing in each subbasin: The shape of these patterns was consistent among years, but not among subbasins. In two subbasins, emergence occurred progressively later with distance upstream, whereas in the other two subbasins emergence was earliest at upstream sites. Within each year, median emergence dates among sites within each subbasin ranged between 44 and 58 days. This spatial variation was comparable to interannual variation, with median emergence dates for a given location in each subbasin ranging between 47 and 74 days among years. Contrary to our expectations, juvenile salmon were not larger in years with earlier emergence, owing to slower estimated spring and summer growth rates compared to years with later emergence. Despite large interannual variation in estimated emergence dates, these results suggest that other factors (e.g., stream flow, temperature, and density‐dependence) were more important than growth duration in determining juvenile salmon growth rates and size among years. We demonstrated considerable spatial and interannual variation in emergence phenology within these subbasins. Understanding how this variation translates to spatiotemporal patterns of juvenile salmon habitat use, growth, and survival has important implications for guiding restoration efforts and understanding how climate change may impact these populations.

Kaylor, Matthew J.↗

DOE OSTI2023

Dynamic modeling and simulation of pressure swing adsorption processes using toPSAil

Pressure swing adsorption (PSA) has attracted significant recent interest for chemical process intensification due to its potential for high energy efficiency and amenability to small, modular designs. However, the lack of simulation tools that are readily available, transparent, and trusted, has been identified as a serious impediment to widespread adoption of PSA, as well as to further research on PSA modeling, numerical solution, optimization, and control. This paper presents a complete framework for dynamic modeling and simulation of PSA processes and its implementation in an open-source simulator called toPSAil. Further, the presentation is tutorial and includes many modeling and implementation details often overlooked in existing literature. Novel methods are presented for handling flow reversals and implementing various pressure–flow relationships, along with controlled boundary conditions. Finally, the code contains several innovations designed to improve efficiency and reduce the extensive trial-and-error tuning often required to produce a working PSA cycle.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Chromium enhances the mechanical performance of 3d transition metal high entropy alloys

A comprehensive examination of the compositional effects on the deformation behavior of high entropy alloys (HEAs) was conducted through nanoindentation, indentation creep, and stress relaxation experiments at ambient temperature. The alloys investigated included NiCoFe, NiCoCr, NiCoCrFe, and NiCoCrFeMn, all of which underwent identical thermomechanical processing. From the experimental results, NiCoCr exhibited the highest maximum shear stress for dislocation nucleation (13 GPa) and nanoindentation hardness (3.8 GPa) among the alloys tested. It also showed the lowest stress sensitivity (0.02) and activation volume (3–4 b 3 ) under creep and stress relaxation conditions, indicative of its superior plastic flow properties. Density functional theory (DFT) calculations further revealed an uneven charge density and larger bonding-length variation in NiCoCr due to the presence of Cr, which increased lattice distortion, impeding dislocation movement. Furthermore, in this study, post-deformation microscopy using transmission electron microscopy (TEM) revealed a high density of stacking faults and twins in NiCoCr that effectively enhanced its creep strength. The results highlight the significance of specific elemental effects, particularly from Cr in this study, over configurational entropy effect (i.e., compositional complexity) in governing the deformation microstructure and mechanical properties of HEAs. Finally, these insights will be instrumental for the design and optimization of advanced alloys for load-bearing applications.

36 - MATERIALS SCIENCE↗

DOE OSTI2022

Characterizing the initial conditions of heavy-ion collisions at the LHC with mean transverse momentum and anisotropic flow correlations

Correlations between mean transverse momentum [pT] and anisotropic flow coefficients v 2 or v 3 are measured as a function of centrality in Pb–Pb and Xe–Xe collisions at $\sqrt{s{NN}}$ = 5.02 TeV and 5.44 TeV, respectively, with ALICE. In addition, the recently proposed higher-order correlation between [p T ], v 2 , and v 3 is measured for the first time, which shows an anticorrelation for the presented centrality ranges. These measurements are compared with hydrodynamic calculations using IP-Glasma and T R ENTo initialstate shapes, the former based on the Color Glass Condensate effective theory with gluon saturation, and the latter a parameterized model with nucleons as the relevant degrees of freedom. The data are better described by the IP-Glasma rather than the T R ENTo based calculations. In particular, Trajectum and JETSCAPE predictions, both based on the T R ENTo initial state model but with different parameter settings, fail to describe the measurements. As the correlations between [p T ] and v n are mainly driven by the correlations of the size and the shape of the system in the initial state, these new studies pave a novel way to characterize the initial state and help pin down the uncertainty of the extracted properties of the quark–gluon plasma recreated in relativistic heavy-ion collisions

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Reactive Modified Epoxy Resin and Its Miscible Blends Based on Recycled Oligomers from Solvolysis

Chemical depolymerization of fully cured epoxy resin with 20% reactive modifier was successfully performed via a solvent-assisted solvolysis process into low molecular weight recyclable oligomers (RO) at 240 °C in a pressure vessel at 650 psi for 4 h. The thermoset epoxy resin was depolymerized into transparent brown viscous fluid with a higher viscosity than the uncured epoxy resin with approximately 93% yield. Different concentrations of the RO were homogeneously mixed with the pure epoxy resin, and their curing kinetics, viscosity, FTIR, mechanical properties, DMA, and cross-link density were investigated. The curing kinetics of the pure reactive modified epoxy resin (baseline) and its mixtures with RO of different concentrations were investigated under both isothermal and nonisothermal conditions using small amplitude oscillatory shear flow. The elastic and viscous moduli (G′ and G″), complex viscosity (η*), and tan δ values were evaluated at different curing times and temperatures. The G′, G″, and η* increased dramatically, while tan δ decreased strongly by several orders of magnitude at the gel point. The zero-shear viscosity (η 0 ) was determined from the angular frequency dependent on η* based on the Cross model for different blend compositions in the liquid state before curing. The composition dependence of η 0 showed a positive deviation from the linear mixing rule and was well described by the Lecyar model. Here, the apparent activation energy of curing (E a ) was also evaluated according to the Arrhenius equation and was found to be 46 ± 2 kJ/mol regardless of the different contents of RO. For all blends up to 40 wt % RO, only one tan δ peak systematically shifting to lower temperatures with increasing content of RO was observed in the DMA measurements, indicating that the epoxy resin and the RO are miscible with up to 40 wt % RO.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Identification of a non-axisymmetric mode in laboratory experiments searching for standard magnetorotational instability

Abstract The standard magnetorotational instability (SMRI) is a promising mechanism for turbulence and rapid accretion in astrophysical disks. It is a magnetohydrodynamic (MHD) instability that destabilizes otherwise hydrodynamically stable disk flow. Due to its microscopic nature at astronomical distances and stringent requirements in laboratory experiments, SMRI has remained unconfirmed since its proposal, despite its astrophysical importance. Here we report a nonaxisymmetric MHD instability in a modified Taylor-Couette experiment. To search for SMRI, a uniform magnetic field is imposed along the rotation axis of a swirling liquid-metal flow. The instability initially grows exponentially, becoming prominent only for sufficient flow shear and moderate magnetic field. These conditions for instability are qualitatively consistent with SMRI, but at magnetic Reynolds numbers below the predictions of linear analyses with periodic axial boundaries. Three-dimensional numerical simulations, however, reproduce the observed instability, indicating that it grows linearly from the primary axisymmetric flow modified by the applied magnetic field.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Simulations of nozzle gas flow and gas-puff Z-pinch implosions on the Weizmann Z-pinch

We present simulations of an oxygen gas puff Z-pinch on a University scale generator at the Weizmann Institute of Science. The work accounts for the detailed geometry of the nozzle, the initial neutral gas density distribution, and the subsequent implosion. The modeling results show significant improvement with data for the current at the time of stagnation in comparison with a previous effort [Rosenzweig et al., Phys. Plasmas 27, 022705 (2020)]. As a first step, we performed simulations of the flow of neutral diatomic oxygen from a plenum through a nozzle within a recessed cathode, across a gap, and into the anode with a recessed grounded honeycomb. These simulations show an agreement with the measured initial gas density profiles within the region not blocked by the recesses and accessible to visible measurements. The computed neutral gas flow profile serves as the initial condition for a radiation magnetohydrodynamic simulation of the implosion using the MACH2-TCRE code. By considering the specific details of the nozzle and chamber geometry, we find agreement with the measured current profile, including the inductive notch. The simulations predict that the plasma undergoes a strong pinch within the hidden anode recess. The simulations also predict the strongest radiation pulse occurs within the anode recess and at the time of the observed inductive notch.

Physics↗