Search NASA⌕ Search

SEARCH · Search NASA

Results for “AXISYMMETRIC FLOW”

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 19 records

One-dimensional temperature measurement of supersonic jet flow by resonantly ionized photoemission thermometry of molecular nitrogen

As the field of fluid dynamics progresses, the demand for sophisticated diagnostic methods to accurately assess flow conditions rises. In this work, resonantly ionized photoemission thermometry (RIPT) has been used to directly target and ionize diatomic nitrogen (N 2 ) to measure one-dimensional (1D) temperature profiles in a supersonic jet flow. This technique can be considered non-intrusive as the premise uses resonantly enhanced multiphoton ionization (REMPI) to target molecular nitrogen. This resonance excites N 2 into absorption bands of the P, Q, and R rotational branches of N 2 ( b 1 Π u ). The ideal (3 + 1) REMPI scheme excites from the ground state and ionizes N 2 ( b 1 Π u ← X 1 Σ g + ) where de-excitation results in photoemission from the first negative band of ionized N 2 + ( B 2 Σ u + → X 2 Σ g + ) as nitrogen returns to the ground state. The resulting emission can be observed using an intensified camera, thus permitting inference of the rotational temperature of ground-state molecular nitrogen. A linearly regressive Boltzmann distribution is applied based on previous calibration data for this technique to quantify the temperature along the ionized line. This work applies this technique to a pure N 2 supersonic jet in cross-flow and counter-flow orientations to demonstrate N 2 RIPT’s applications in a supersonic flow. Temperature variations are observed at different locations downstream of the exit in cross-flow, and axisymmetric in counter-flow, to generate profiles characterizing the flow dynamics. Due to the collisional effects resulting from the number density of N 2 at higher pressures, a (3 + 2) REMPI scheme is observed throughout this text.

Clark, Aleksander (ORCID:0009000855600042)↗

Global Non-axisymmetric Hall Instabilities in a Rotating Plasma

Non-axisymmetric, flow-driven instabilities in the incompressible Hall-MHD model are studied in a differentially rotating cylindrical plasma. It is found that, in the Hall-MHD regime, both whistler waves and ion-cyclotron waves can extract energy from the flow shear, resulting in two distinct branches of global instability. The non-axisymmetric whistler modes grow significantly faster than non-axisymmetric, ideal MHD modes. A discussion of the global whistler instability mechanism is presented in the large-ion-skin-depth, “electron-MHD” limit. When the magnetic field is azimuthal, a subset of the whistler modes having zero axial wave number are uncovered to be destabilized by the “corotation amplifier” mechanism. It is observed that the effect of the Hall term on the non-axisymmetric modes can be appreciable when d i is on the order of a few percent of the width of the cylindrical annulus. Distinct global modes emerge in the strong Hall-MHD regime at significantly stronger magnetic fields than those required for unstable global MHD modes, as the Hall effect weakens the stabilizing “field-line bending” by decoupling ion motion from the magnetic field. These global non-axisymmetric modes may play an important role in weakly ionized accretion disks.

Alfven waves↗

Boundary-layer receptivity to oblique freestream vorticity waves for a high-enthalpy hypersonic flow

The receptivity of a Mach 15 straight-cone boundary layer to oblique freestream vorticity waves is investigated using direct numerical simulation (DNS) alongside linear stability theory and the linear parabolized stability equations. A thermochemical nonequilibrium gas model is used. Oblique freestream vorticity waves at frequencies of 400, 800, and 1200 kHz are considered, with incident angles ranging from 0° to 29.4° at 400 kHz, 0° to 15.7° at 800 kHz, and 0° to 20.6° at 1200 kHz. The 400 kHz case is of primary interest due to the strong second-mode amplification at this frequency. Although the underlying base flow is axisymmetric, the oblique vorticity waves lead to a fully three-dimensional boundary-layer disturbance whose characteristics vary depending on the azimuthal ray relative to the freestream wave. Moving downstream within the second-mode instability region, some clear trends emerge in terms of the boundary-layer disturbance amplitudes; that is, disturbance amplitudes are highest at the leeward ray (relative to the freestream wave), but weakest about halfway between the windward and leeward rays. Increasing the incident angle causes the amplitudes to increase on the leeward ray and decrease on the windward ray. Moreover, the boundary-layer disturbance throughout contains a wide spectrum of azimuthal wavenumbers in which the disturbance energy falls off at higher wavenumbers. Increasing the incident angle causes the azimuthal spectrum of the boundary-layer disturbance to broaden overall. Qualitatively similar results are found for the two higher frequencies leading up to the peak-amplitude locations corresponding to the second-mode instability.

42 ENGINEERING↗

Analysis of the density field around a supersonic conical projectile using quantitative schlieren

Quantitative schlieren imaging is a flow measurement technique that is capable of measuring density fields throughout refractive flowfields. The technique was applied here to measure the density field surrounding supersonic conical projectiles in free flight. Shock waves attached to a supersonic conical projectile offer a simple geometry with the well-established Taylor–Maccoll analytical flow solution to which these experimental measurements were compared. The schlieren images recorded a projection of the index of refraction field surrounding the 10° half-angle cones which was converted to density first through an Abel inversion and then the Gladstone–Dale law. Three Abel inversion methods—two-point, three-point, and arbitrary ray axisymmetric projection (ARAP)—were applied to deconvolute the three-dimensional flow within the constrained axisymmetric flow field. The resulting reconstructed density profiles were compared to the Taylor–Maccoll solution, parameterized by cone geometry and Mach number. The experimental density fields demonstrated strong agreement with the theoretical profiles. Experimental consistency was confirmed across various projectile speeds, demonstrating quantitative schlieren’s capability to accurately reconstruct the density of the flow field, even within the resolution constraints imposed by high-speed imaging. In conclusion, an assessment of experimental uncertainties in the density reconstruction was performed.

Abel Inversion↗

Observation of Nonaxisymmetric Standard Magnetorotational Instability Induced by a Free-Shear Layer

The standard magnetorotational instability (SMRI) with a magnetic field component parallel to the rotation axis is widely believed to be responsible for the fast accretion in astronomical disks. In conventional base flows with a Keplerian profile or an ideal Couette profile, most studies focus on axisymmetric SMRI, since excitation of nonaxisymmetric SMRI in such flows requires a magnetic Reynolds number (Rm) more than an order of magnitude larger. Here, we report that, in a magnetized Taylor-Couette flow, nonaxisymmetric SMRI with an azimuthal mode number m = 1 can be triggered by a free-shear layer in the base flow at Rm ≳1, the same threshold as for axisymmetric SMRI. Global linear analysis reveals that the free-shear layer reduces the required Rm, possibly by introducing an extremum in the vorticity of the base flow. Nonlinear simulations validate the results from linear analysis and confirm that a novel instability recently discovered experimentally [Wang et al., Nat. Commun. 13, 4679 (2022)] is the nonaxisymmetric m =1 SMRI. Further, our finding has astronomical implications as free-shear layers are ubiquitous in celestial systems, such as the disk-star boundary layer, the solar tachocline, and the edge of planet-opened gaps in protoplanetary disks.

79 ASTRONOMY AND ASTROPHYSICS↗

Generalized Integral Boundary Layer Equations in Streamline Curvilinear Coordinates

The MEIT equations are derived in orthogonal streamline curvilinear coordinates defined by the inviscid velocity at the wall and the surface-normal direction. The resulting formulations is (1) is not restricted to flat-plate or axisymmetric geometries, (2) accommodates vortical inviscid flows (e.g., induced by curved shocks) and wall mass blowing, and (3) allows for an arbitrary definition of the edge state. The streamline ODEs are also extended to a surface PDE framework.

97 MATHEMATICS AND COMPUTING↗

Simplifying the creation of thermal decomposition mechanisms for designing mixed fibre composites

Among the most challenging aspects of simulating thermal decomposition of fibre reinforced polymers is the determination of appropriate reaction parameters. Thermogravimetric analysis is typically used to generate decomposition data, to which the reaction parameters are then fit. When designing mixed fibre materials (e.g. combinations of glass and carbon fibres), the number of TGA experiments needed to explore the entire design space may be intractable. Here, we demonstrate the creation of two candidate proxy mixed fibre mechanisms and compare them to a mechanism created by fitting parameters from TGA on the mixed fibre composite. These mechanisms are then demonstrated in a 2D axisymmetric numerical decomposition, heat transfer, and porous flow model. We find a maximum 11% uncertainty in mass and 4% in temperature difference when using a proxy mechanism.

Scott, Sarah N. [Sandia National Lab. (SNL-CA), Li↗

CFD simulations of Molten Salt Fast Reactor core cavity flows

Computational Fluid Dynamics (CFD) has become increasingly important in the research and development of advanced nuclear reactors. Here, in the current study, extensive CFD simulations were conducted for the coolant flow in Molten Salt Fast Reactor (MSFR) core models using the state-of-the-art spectral element flow solver Nek5000 and multiscale coarse-mesh thermal-hydraulic software Pronghorn. The underlying motivation is to seek an in-depth understanding of how the internal velocity distribution can be influenced by the MSFR core cavity shape, the Reynolds number, turbulence modeling options and the inlet boundary conditions. The CFD techniques involved in this investigation range from coarse-mesh CFD, RANS modeling, to the high-fidelity LES calculations. Specifically, a series of RANS simulations were performed for the 2-D axisymmetric core model and 3-D wedge domains to study the flow distribution inside the MSFR core. It is observed that a proper representation of the MSFR inlet channel duct is important for the prediction of internal flow distribution. It is also showcased here how researchers can leverage the Nek5000 CFD results to calibrate more efficient coarse-mesh CFD tools, like Pronghorn, for the actual MSFR design needs. Moreover, this paper highlights a 3-D LES model for an entire MSFR core using the spectral element method and demonstrates the feasibility of this modeling approach. The readiness and potential limitations of the RANS approach are examined with respect to the high-fidelity LES simulations. The present investigation lays a solid foundation as we are leveraging the high-fidelity CFD capabilities to inform MSFR design efforts.

97 MATHEMATICS AND COMPUTING↗

Effect of parallel flow on resonant layer responses in high beta plasmas

Abstract Resonant layers in a tokamak respond to non-axisymmetric magnetic perturbations by amplifying the mode amplitude and balancing the plasma rotation through magnetic reconnection and force balance, respectively. This resonant response can be characterized by local layer parameters and especially by a single quantity in the linear regime, the so-called inner-layer Δ. The computation of Δ under two-fluid drift-MHD formalism has been progressed by reducing the order of the system in the phase space, where the shielding current is approximated as being only carried by electrons, a posteriori . In this study, we relax the approximation and compute Δ accounted for by the parallel flow associated with the ion shielding current. The posteriori is numerically verified in great agreement with the original SLAYER developed in a previous paper (J.-K. Park 2022 Phys. Plasmas 29 072506). Extending the resonant layer response theory to high β plasmas, our research findings answer two important questions: how the parallel flow influences the resonant layer response and why the parallel flow effect appears in high β plasmas. The complicated plasma compression in high β regime allows the parallel flow response to give rise to the ion shielding current, which not only shifts the zero-crossing condition of the ExB flow but also enhances the field penetration threshold. Technically, the Riccati matrix transformation method is adapted to handle the numerical stiffness due to the increased order of the system. The high fidelity of this numerical method makes use of further extension of the model to higher-order systems to take other physical phenomena into account. This work is envisaged to predict the resonant layer response under high β fusion reactor conditions.

Lee, Yeongsun (ORCID:000000034474416X)↗

Particle Acceleration in Collisionless Magnetically Arrested Disks

We present the first collisionless realization of two-dimensional axisymmetric black hole accretion consistent with a persistent magnetically arrested disk state. The accretion flow, consisting of an ion-electron disk plasma combined with magnetospheric pair creation effects, is simulated using first-principles general-relativistic particle-in-cell methods. The simulation is evolved over significant dynamical timescales during which a quasisteady accretion state is reached with several magnetic flux eruption cycles. We include a realistic treatment of inverse Compton scattering and pair production, which allows for studying the interaction between the collisionless accretion flow and pair-loaded jet. Our findings indicate that magnetic flux eruptions associated with equatorial magnetic reconnection within the black hole magnetosphere and the formation of spark gaps are locations of maximal particle acceleration. Flux eruptions, starting near the central black hole, can trigger Kelvin-Helmholtz-like vortices at the jet-disk interface that facilitate efficient mixing between disk and jet plasma in this region. Transient periods of increased pair production following magnetic flux eruptions and reconnection events are responsible for most of the highly accelerated particles.

Accretion disk & black-hole plasma↗

Simulations of vertical displacement oscillatory modes and global Alfvén Eigenmodes in JET geometry

Vertical Displacement Oscillatory Modes (VDOM), with frequency in the Alfvén range, are natural modes of oscillation of magnetically confined laboratory plasmas with elongated cross-section. These axisymmetric modes arise from the interaction between the plasma current, which is in equilibrium with currents flowing in external coils, and perturbed currents induced on a nearby conducting wall. The restoring force exerted by these perturbed currents on the vertical motion of the plasma column leads to its oscillatory behavior. An analytic model for VDOM was proposed based on an idealized 'straight tokamak' equilibrium with uniform equilibrium current density. This article introduces the first numerical simulations of VDOM in a realistic JET tokamak configuration, using the extended-MHD code NIMROD and drawing comparisons with Global Alfvén Eigenmodes (GAE). The results show qualitative agreement with analytic predictions regarding mode frequency and radial structure, supporting the identification of VDOM as a fundamental oscillation mode in tokamak plasmas. VDOM and GAE are modeled in a representative JET discharge, where axisymmetric perturbations with toroidal mode number n = 0 driven unstable by fast ions were observed. The two modes are examined separately using a forced oscillator within the NIMROD code, which enables a comparison of their characteristics and helps identify the experimentally observed mode possibly as a GAE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Validation of predictive performance models for supersonic gas-jet nozzles at the Laboratory for Laser Energetics

We present results characterizing the neutral-density distributions produced by the supersonic nozzles used in experiments on the OMEGA-60 and OMEGA-EP laser systems at the University of Rochester’s Laboratory for Laser Energetics (LLE). Axisymmetric Fluent® simulations using LLE nozzle specifications capture the viscous effects, gas expansion, and shock waves that complicate flow predictions for offsets above the nozzle exit. Here, these simulations show good agreement with neutral-density measurements obtained using a four-wave shearing interferometer. An analytical form is given for the plateau length. Fits to simulation data for boundary layer thickness, mean plateau density, and density ramps are given as functions of nozzle offset and nozzle backing pressure for a number of nozzles and gases.

47 OTHER INSTRUMENTATION↗

Investigating nonlinear magnetohydrodynamics in an optimized, reactor-scale quasi-axisymmetric stellarator

We use the new simulation capabilities of the extended-magnetohydrodynamic (MHD) code, M3D-C1, to investigate the nonlinear MHD properties of a reactor-scale quasisymmetric stellarator equilibrium. Our model captures the self-consistent evolution of the magnetic field, temperature, density, and flow profiles without imposing restrictions on the structure of the first. We include the effects of resistivity using a realistic temperature-dependent Spitzer model, along with a model for heat transport that captures the key physical characteristic, namely, strongly anisotropic diffusion in directions perpendicular and parallel to the magnetic field. We consider a quasi-axisymmetric, finite-pressure equilibrium that was optimized for self-consistent bootstrap current, quasi-symmetry, and energetic particle confinement. Our assessment finds that the equilibrium is highly unstable to interchange-like pressure-driven instabilities near the plasma edge. The initially unstable modes rapidly destabilize other modes in the direction of the N-fold rotational symmetry (toroidal, in this case). For this equilibrium, N = 2, meaning destabilization of a large number of even-numbered toroidal Fourier modes. Thus, field-periodicity is likely to be an important factor in the nonlinear MHD stability characteristics of optimized stellarators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Application of the portable diagnostic package to the Wisconsin high-temperature-superconducting axisymmetric mirror (WHAM)

We present an application of the Portable Diagnostic Package (PDP) on the Wisconsin High-temperature-superconducting Axisymmetric Mirror (WHAM), which integrates an optical emission spectroscopy (OES) system and an active Thomson scattering (TS) system. The OES system facilitates a comprehensive impurity line survey and enables flow measurements through the Doppler effect observed on impurity lines. Plasma rotation profiles were successfully derived from doubly charged carbon lines. The TS system enabled the first measurements of the electron temperature in commissioning plasmas on WHAM. Notably, the PDP was installed, commissioned, and used to obtain OES and TS data within ∼6 months from project start, enabled by its designed portability and standardized interfaces. These results demonstrate the PDP’s potential to accelerate diagnostic readiness and advance experimental plasma studies.

Doppler effect↗

Structure Formation Through Magnetohydrodynamical Instabilities in Primordial Disks

The shear flow instabilities under the presence of magnetic fields in the primordial disk can greatly facilitate the formation of density structures that serve as seeds prior to the onset of the gravitational Jeans instability. We evaluate the effects of the Parker, magnetorotational and kinematic dynamo instabilities by comparing the properties of these instabilities. We calculate the mass spectra of coagulated density structures by the above mechanism in the radial direction for an axisymmetric magnetohydrodynamic (MHD) torus equilibrium and power density profile models. Our local three-dimensional MHD simulation indicates that the coupling of the Parker and magnetorotational instability creates spiral arms and gas blobs in an accretion disk, reinforcing the theory and model. Such a mechanism for the early structure formation may be tested in a laboratory. The recent progress in experiments involving shear flows in rotating tokamak, field reversed configuration (FRC) and laser plasmas may become a key element to advance in nonlinear studies.

Noguchi, Koichi↗

Zonal magnetic fields regulate nonlinear edge-localized-mode dynamics via self-consistent force balance

Edge-localized modes (ELMs) eject intense bursts of heat and particles that threaten plasma-facing components in fusion reactors. Nonlinear full-torus BOUT++ simulations show that turbulence-driven zonal magnetic fields (ZMFs) play an essential role in nonlinear ELM evolution by maintaining self-consistent force balance. Zonal flows mitigate the initial crash through shear but do not prevent continued radial transport. When ZMFs are self-consistently included, turbulence-driven zonal currents modify the parallel current distribution and magnetic tension and are associated with a reduction of the axisymmetric (𝑛 = 0) perturbed radial force imbalance. This coincides with a transition from convective, bursty propagation to more localized, diffusive transport. Similar behavior is observed across the regimes considered, including both resistive-ballooning and peeling-ballooning cases. Finally, associated signatures, including radial electric field shear and parallel current redistribution, provide experimentally accessible diagnostics for present devices and ITER-relevant conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enabling "Fuel" Switching through On-Demand Wind Control

This work demonstrates that classical shear-flow stability theory can be successfully applied to modify wind turbine wakes and also explains the success of several emerging, empirically-arrived control methods (i.e., dynamic induction and helix control). Linear stability theory predictions are extended to include the effects of non-axisymmetric inflow profiles, such as wind shear, which is shown to not strongly affect the primary forcing frequency. The predictions, as well as idealized large-eddy simulations using actuator-line representation of the turbine blades, agree that the n = 0 and ±1 modes have faster initial growth rates than higher-order modes, suggesting the lower-order modes are more appropriate for wake control. Exciting the lower-order modes with periodic pitching of the blades produces higher entrainment into the wake and consequently faster wake recovery.

17 WIND ENERGY↗

Toward Understanding the Differences between Mesoscale and Large-Eddy Simulations of Tropical Cyclones

In this work, we investigate the ability of mesoscale and large-eddy simulation (LES) model configurations to predict the mean wind speed profile within the boundary layer of tropical cyclones (TCs). To this end, we perform idealized simulations of five hypothetical intense storms ranging from categories 1 to 5 on the Saffir–Simpson scale and extract time-averaged quantities near the eyewall region. We compare the model-generated data against mean wind speed profiles compiled from dropsondes launched from reconnaissance aircraft operating in the North Atlantic basin. Our analysis shows that mesoscale- and LES-generated mean wind fields display important differences in the boundary layer, including the magnitude of shear as well as the height where their low-level wind speed maxima are located. In addition, a comparison between the two model configurations with the dropsonde data shows that both modeling approaches are unable to capture the typical structure of mean winds in the lower part of the TC boundary layer (10–500 m), calling into question the use of simulations of near-axisymmetric storms for investigating the wind structure of past events. To better understand these differences, we conduct a momentum-budget analysis and show that modeled turbulent fluxes are underestimated in the mesoscale boundary layer parameterization compared to the LES model. Based on the analysis of the horizontal turbulent fluxes and their potential impact on mean flow quantities, a TC-specific boundary layer parameterization may be needed.

17 WIND ENERGY↗