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.

101 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effects of collisional ion orbit loss on neoclassical tokamak radial electric fields

Ion orbit loss is considered important for generating the radially inward electric field Er in a tokamak edge plasma. In particular, this effect is emphasized in diverted tokamaks with a magnetic X point. In neoclassical equilibria, Coulomb collisions can scatter ions onto loss orbits and generate a radially outward current, which in steady state is balanced by the radially inward current from viscosity. To quantitatively measure this loss-orbit current in an edge pedestal, an ion-orbit-flux diagnostic has been implemented in the axisymmetric version of the gyrokinetic particle-in-cell code XGC. As the first application of this diagnostic, a neoclassical DIII-D H-mode plasma is studied using gyrokinetic ions and adiabatic electrons. The validity of the diagnostic is demonstrated by studying the collisional relaxation of Er in the core. After this demonstration, the loss-orbit current is numerically measured in the edge pedestal in quasisteady state. In this plasma, it is found that the radial electric force on ions from Er approximately balances the ion radial pressure gradient in the edge pedestal, with the radial force from the plasma flow term being a minor component. The effect of orbit loss on Er is found to be only mild.

gyrokinetic simulations↗
DOE OSTI2024

Liquid lithium divertor analysis using coupled plasma material interaction model

A liquid lithium divertor can improve performance of future fusion devices by creating efficient power exhaust and improving the energy confinement via pumping of the hydrogen isotopes. In addition, significantly higher heat fluxes can be handled if controlled vapor shielding is used to redistribute the divertor heat flux over a wider area. Design and optimization of such a system calls for an analysis model which includes a strong two-way coupling between the plasma and divertor material. The incoming plasma heat and particle flux will affect the divertor surface temperature, which is a defining factor of the lithium evaporative and sputtered flux going into the plasma. Results of the coupled model based on the plasma edge code SOLPS-ITER and the computational fluid dynamics (CFD) code ANSYS-CFX will be presented for different configurations. An analytical slab flow model is used as a heat transfer boundary condition for SOLPS, defining particle flux from the wall via calculation of the surface temperature. At the final step, results of the SOLPS analysis are verified using a 3D CFD magnetohydrodynamics (MHD) analysis which uses heat and particle flux from SOLPS as a boundary condition. In addition to plasma heat flux, both analytical and CFD temperature models include several plasma material interaction effects, such as lithium evaporation, condensation and sputtering based on deuterium target flux. New adatom sputtering model based on the available experimental data is presented. Analytical model is expanded to include free surface axisymmetric configurations. Results of parametric studies of the divertor configurations with different lithium inlet temperature and velocity will be presented leading to the optimal design resulting in the lowest possible lithium contamination in the core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Effects of collisional ion orbit loss on neoclassical tokamak radial electric fields

We report that ion orbit loss is considered important for generating the radially inward electric field E r in a tokamak edge plasma. In particular, this effect is emphasized in diverted tokamaks with a magnetic X point. In neoclassical equilibria, Coulomb collisions can scatter ions onto loss orbits and generate a radially outward current, which in steady state is balanced by the radially inward current from viscosity. To quantitatively measure this loss-orbit current in an edge pedestal, an ion-orbit-flux diagnostic has been implemented in the axisymmetric version of the gyrokinetic particle-in-cell code XGC. As the first application of this diagnostic, a neoclassical DIII-D H-mode plasma is studied using gyrokinetic ions and adiabatic electrons. The validity of the diagnostic is demonstrated by studying the collisional relaxation of E r in the core. After this demonstration, the loss-orbit current is numerically measured in the edge pedestal in quasisteady state. In this plasma, it is found that the radial electric force on ions from E r approximately balances the ion radial pressure gradient in the edge pedestal, with the radial force from the plasma flow term being a minor component. The effect of orbit loss on E r is found to be only mild.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

RANS Simulation of Variable Density Turbulent Round Jets with Coflow using xRAGE Hydrodynamic Code and BHR Turbulence Models

This work for the fiscal year 2023 (FY23) is a continuation of previous efforts to evaluate the BHR turbulence models for their ability to accurately simulate variable density turbulent round jets with coflow. As before, RANS simulations are carried out using the xRAGE hydrodynamic code. The following are some of the previous findings. Israel showed that i) the symmetry boundary conditions for the BHR models in axisymmetric simulations were in error, ii) three grids of different resolutions did not lead to converging solutions, and iii) BHR 3.1 simulation exhibited instabilities and did not reach a steady state. Saenz and Rauenzahn derived and implemented into xRAGE the correct BHR boundary conditions at the symmetry axis. Cline conducted sensitivity studies with various parameters including the BHR models (versions 2, 3.1 and 4), gravity, material pressure, specific heat, initial turbulent kinetic energy and initial turbulent length scale, and found that the largest factor impacting on simulation results was the BHR model version, followed by the initial turbulent length scale. In addition, freeze boundary conditions at the exit and the side wall of the computational domain were used to remove anomalous flow behavior. Cline adjusted the inlet jet velocity, initial turbulent kinetic energy and initial turbulent length scale to obtain the best reasonable match with the experimental data of Charonko and Prestridge. It was found that the BHR 2 and 3.1 models performed in a similar manner, but the BHR 2 model produced much lower levels of density-specific-volume covariance and turbulent kinetic energy. The main focus for the FY23 is to study the effects of computational parameters related to the boundary conditions, mesh, domain size and timestep size. The reasoning behind this is that, unless simulation results can be shown to be reasonably independent from the aforementioned computational parameters, it would be difficult to attribute any discrepancies between simulation and experimental results to turbulence models. This important aspect has largely been overlooked in the previous years, and therefore will be studied comprehensively here. Additionally, effects of varying the initial turbulent length scale will be examined because it was previously identified as a major factor affecting the flow fields.

42 ENGINEERING↗
DOE OSTI2024

Inverse aspect-ratio expanded tokamak equilibria

Following Greene et al. [Phys. Fluids 14, 671 (1971)] and Connor et al. [Phys. Plasmas 31, 577 (1988); Plasma Phys. Control. Fusion 34, 161 (1992); and Nucl. Fusion 33, 1533 (1993)], the Grad-Shafranov equation for an axisymmetric tokamak plasma equilibrium is solved via an expansion in the, supposedly small, inverse aspect-ratio of the plasma, ϵ. The displacements of equilibrium magnetic flux-surfaces due to plasma shaping are assumed to be $\mathcal{O}$(ϵ) smaller than the minor radii of the surfaces, but no other restriction is placed on the nature of the shaping. The solution of the Grad-Shafranov equation is matched to a vacuum solution that extends to infinity, and consists of an expansion in toroidal functions. The external poloidal magnetic field generated by a finite set of discrete external poloidal magnetic field-coils is calculated, and incorporated into the toroidal function expansion. In this manner, the shape of a large aspect-ratio tokamak plasma is directly related to the currents flowing in the external poloidal field-coils. Finally, a pedestal in the plasma pressure, and the associated spike in the bootstrap current, are incorporated into the model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The influence of boundary and edge-plasma modeling in computations of axisymmetric vertical displacement

A number of previously published studies apply extended-MHD computations to model macroscopic dynamics of tokamak vertical displacement events (VDEs). The characteristic plasma-wall contact and resistive-wall diffusion imply sensitivity to boundary conditions in a general sense, but which conditions have significant influence depends on how the plasma is modeled. This work examines these dependencies by comparing results of axisymmetric extended-MHD computations with different sets of boundary conditions, plasma thermal-conduction models, and wall resistivity values. The geometry and plasma parameters of the computations represent a modest-sized tokamak. The forced-displacement scenario and computational setup are similar to those used in previous 3D computations. The results show that for a given wall resistivity, the VDE time-scale is most sensitive to variations in the boundary and thermal-conduction parameters that affect electron thermal transport. The electrical conductivity depends on electron temperature, and the dependence on thermal transport stems from its influence on the electrical circuit that includes the open-field halo current. Lastly, conditions that lead to hotter, broader halo regions slow the evolution. Significant sensitivity to the boundary condition on plasma flow-velocity exists when electron thermal conduction is restricted and electron energy loss is convective, which is expected for conditions at the entrance of the magnetic presheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Effect of triangularity on ion-temperature-gradient-driven turbulence

The linear and nonlinear properties of ion-temperature-gradient-driven turbulence with adiabatic electrons are modeled for axisymmetric configurations for a broad range of triangularities δ, both negative and positive. Peak linear growth rates decrease with negative δ but increase and shift toward a finite radial wavenumber kx with positive δ. The growth-rate spectrum broadens as a function of kx with negative δ and significantly narrows with positive δ. Here, the effect of triangularity on linear instability properties can be explained through its impact on magnetic polarization and curvature. Nonlinear heat flux is weakly dependent on triangularity for |δ|≤0.5, decreasing significantly with extreme δ, regardless of sign. Zonal modes play an important role in nonlinear saturation in the configurations studied, and artificially suppressing zonal modes increased nonlinear heat flux by a factor of about four for negative δ, increasing with positive δ by almost a factor of 20. Proxies for zonal-flow damping and drive suggest that zonal flows are enhanced with increasing positive δ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Non-thermal electron acceleration from magnetically driven reconnection in a laboratory plasma

Magnetic reconnection rapidly converts magnetic energy into some combination of plasma flow energy, thermal energy and non-thermal energetic particles. Various reconnection acceleration mechanisms have been theoretically proposed and numerically studied in different collisionless and low-β environments, where β refers to the plasma-to-magnetic pressure ratio. These mechanisms include Fermi acceleration, betatron acceleration, parallel electric field acceleration along magnetic fields and direct acceleration by the reconnection electric field. However, none of them have been experimentally confirmed, as the direct observation of non-thermal particle acceleration in laboratory experiments has been difficult due to short Debye lengths for in situ measurements and short mean free paths for ex situ measurements. Here we report the direct measurement of accelerated non-thermal electrons from magnetically driven reconnection at low β in experiments using a laser-powered capacitor coil platform. We use kilojoule lasers to drive parallel currents to reconnect megagauss-level magnetic fields in a quasi-axisymmetric geometry. The angular dependence of the measured electron energy spectrum and the resulting accelerated energies, supported by particle-in-cell simulations, indicate that the mechanism of direct electric field acceleration by the out-of-plane reconnection electric field is at work. Scaled energies using this mechanism show direct relevance to astrophysical observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Optimum Injector Parameters for Thermoacoustic Stability in a Multi-Nozzle Can Combustion System

Abstract High-frequency transverse instabilities are an important concern in can combustor configurations. In these configurations which are typically operated with multiple injectors around a central injector, each injector is subjected to different parts of the acoustic mode shape and thus respond differently for the same instability mode. Recent work by the author has modeled the response of premixed flames to excitation by natural high-frequency transverse modes in a can combustor both in the center and outer nozzles. The stability of these acoustically non-compact flames was assessed using the Rayleigh criterion (Rayleigh Integral denoted as RI) and not the overall unsteady heat release as is the case for compact flames. Several key control parameters were studied, namely — flame angle, swirling strength, nozzle location. For non-axisymmetric modes such as the commonly occurring 1-T mode, both radial and azimuthal offsets of the nozzle location affected stability. The framework was applied to an optimization study to identify the optimal combination of parameters that minimizes RI for the different nozzles in the multi-nozzle system. In this study, a N-around-1 configuration was studied, and the results indicated that the different nozzles needed to be operated at different flame angles and swirl numbers to result in an overall minimum RI. However, the specific response of the different injectors was not considered. The helical mode distribution at each injector varies as we azimuthally go around the combustor’s injector distribution and thus the most amplified mode and the resulting flame response would be different. To minimize RI, it is important to determine the injector configurations that result in a hydrodynamic profile that minimizes the individual RI for each nozzle. The resulting relationship between the injector’s flow and local hydrodynamics can then be used in a hydrodynamics study of an individual injector so that the most optimal injector is chosen depending on its location in the combustor dump plane.

Acharya, Vishal↗
DOE OSTI2020

Validation of two-layer model for underexpanded hydrogen jets

Previous studies have shown that the two-layer model more accurately predicts hydrogen dispersion than the conventional notional nozzle models without significantly increasing the computational expense. However, the model was only validated for predicting the concentration distribution and has not been adequately validated for predicting the velocity distributions. In the present study, particle imaging velocimetry (PIV) was used to measure the velocity field of an underexpanded hydrogen jet released at 10 bar from a 1.5 mm diameter orifice. The two-layer model was the used to calculate the inlet conditions for a two-dimensional axisymmetric CFD model to simulate the hydrogen jet downstream of the Mach disk. The predicted velocity spreading and centerline decay rates agreed well with the PIV measurements. The predicted concentration distribution was consistent with data from previous planar Rayleigh scattering measurements used to verify the concentration distribution predictions in an earlier study. The jet spreading was also simulated using several widely used notional nozzle models combined with the integral plume model for comparison. These results show that the velocity and concentration distributions are both better predicted by the two-layer model than the notional nozzle models to complement previous studies verifying only the predicted concentration profiles. Thus, this study shows that the two-layer model can accurately predict the jet velocity distributions as well as the concentration distributions as verified earlier. Though more validation studies are needed to improve confidence in the model and increase the range of validity, the present work indicates that the two-layer model is a promising tool for fast, accurate predictions of the flow fields of underexpanded hydrogen jets.

08 HYDROGEN↗
DOE OSTI2024

Model-based optimization strategies for direct hydrogenation of carbon dioxide to dimethyl ether

Here, this study discusses model-based optimization strategies for CO 2 hydrogenation to dimethyl ether (DME) over aCuZnZr (CZZ) and ferrite (FER) mixed catalyst system in a packed-bed reactor configuration. A two-dimensional axisymmetric, nonisothermal packed-bed reactor model was developed using COMSOL Multiphysics 6.2 software. The model solves two-dimensional (radial and axial) heat and mass transport equations in the packed-bed and integrates intraparticle diffusion and heat transfer in a 1D approach. This powerful feature differs from a traditional porous media approach and takes into account any heat and mass transfer limitations that may exist. Analysis shows that the heat transfer limitations are negligible, but strong internal mass transfer limitations were observed at 10 ≤ WHSV ≤ 90 h -1 on the FER catalyst and at 240 °C. The optimum catalyst composition (i.e., mixing ratio) varies depending on the operating regime. The FER catalyst weight in the mixture can be as low as 5 wt.%, but the ideal composition de-pends on the internal mass transfer limitation and its relationship with the operating regime (i.e., weight hourly space velocity, temperature). A catalyst composition of 80 wt.% CZZ and 20 wt.% FER was suggested; this composition can provide high CO 2 conversion and DME production rates at a wide range of temperatures and flow rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗