Search NASA⌕ Search

SEARCH · Search NASA

Results for “Turbulent Heating”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A Hybrid Bulk Algorithm to Predict Turbulent Fluxes over Dry and Wet Bare Soils

Measurements made in the Columbia River basin (Oregon) in an area of irregular terrain during the second Wind Forecast Improvement Project (WFIP2) field campaign are used to develop an optimized hybrid bulk algorithm to predict the surface turbulent fluxes from readily measured or modeled quantities over dry and wet bare or lightly vegetated soil surfaces. The hybrid (synthetic) algorithm combines (i) an aerodynamic method for turbulent flow, which is based on the transfer coefficients (drag coefficient and Stanton number), roughness lengths, and Monin–Obukhov similarity; and (ii) a modified Priestley–Taylor (P-T) algorithm with physically based ecophysiological constraints, which is essentially based on the surface energy budget (SEB) equation. Soil heat flux in the latter case was estimated from measurements of soil temperature and soil moisture. In the framework of the hybrid algorithm, bulk estimates of the momentum flux and the sensible heat flux are derived from a traditional aerodynamic approach, whereas the latent heat flux (or moisture flux) is evaluated from a modified P-T model. Direct measurements of the surface fluxes (turbulent and radiative) and other ancillary atmospheric/soil parameters made during WFIP2 for different soil conditions (dry and wet) are used to optimize and tune the hybrid bulk algorithm. The bulk flux estimates are validated against the measured eddy-covariance fluxes. We also discuss the SEB closure over dry and wet surfaces at various time scales based on the modeled and measured fluxes. Although this bulk flux algorithm is optimized for the data collected during the WFIP2, a hybrid approach can be used for similar flux-tower sites and field campaigns.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Evolution of intermittent filaments in the scrape-off layer of NSTX

Filamentary structures naturally arise from background turbulence in the scrape-off layer (SOL) of plasmas, leading to significant particle and heat transport that can degrade overall plasma confinement. This enhanced transport can contribute to unacceptably high heat loads on plasma-facing components. As such, understanding the physics of SOL plasma filaments is critical for predicting and mitigating their effects in future fusion devices. On the National Spherical Torus Experiment (NSTX), plasma filaments—commonly referred to as “blobs”—were investigated using the gas-puff imaging (GPI) diagnostic in the edge and SOL regions. The analysis involved identifying, segmenting, and tracking the characteristic contours of the blobs in each frame of the GPI video sequences. Their evolution was characterized through shape descriptors, velocity, and angular velocity derived from their contour coordinates. The results indicate that as the blob area increases, their shapes become more concave and less circular, suggesting reduced structural stability in larger blobs. This result aligns with previous theoretical results where it was shown that larger blobs are more susceptible to instabilities [Krasheninnikov et al., J. Plasma Phys. 74, 679–717 (2008) and D'Ippolito et al., Phys. Plasmas 18, 060501 (2011)]. A positive correlation was observed between radial velocity and radial position, suggesting radially outward acceleration of the filaments, potentially driven by decreasing viscous drag toward the far SOL. Interestingly, blobs in background SOL turbulence exhibited minimal spinning in contrast to filaments originating from edge localized modes, which show substantial rotation during their paths [Lampert et al., Phys. Plasmas 29, 102502 (2022)]. Statistical analysis of the solidity and total curvature shape descriptors, along with their temporal evolution, revealed relatively broad, near-Gaussian distributions. This suggests that blob morphology is strongly influenced by stochastic turbulent processes in the surrounding plasma environment. Blob parameters were also compared with bulk plasma and radial profile measurements. Notable trends were found between blob rotation and poloidal velocity with collisionality and line-integrated density. These findings contribute to a deeper understanding of blob dynamics and provide valuable insights for refining SOL turbulence models.

Covariance and correlation↗
DOE OSTI2022

Reconnection-driven energy cascade in magnetohydrodynamic turbulence

Magnetohydrodynamic turbulence regulates the transfer of energy from large to small scales in many astrophysical systems, including the solar atmosphere. We perform three-dimensional magnetohydrodynamic simulations with unprecedentedly large magnetic Reynolds number to reveal how rapid reconnection of magnetic field lines changes the classical paradigm of the turbulent energy cascade. By breaking elongated current sheets into chains of small magnetic flux ropes (or plasmoids), magnetic reconnection leads to a previously undiscovered range of energy cascade, where the rate of energy transfer is controlled by the growth rate of the plasmoids. As a consequence, the turbulent energy spectra steepen and attain a spectral index of -2.2 that is accompanied by changes in the anisotropy of turbulence eddies. The omnipresence of plasmoids and their consequences on, for example, solar coronal heating, can be further explored with current and future spacecraft and telescopes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Intracore Natural Circulation Study in the High Temperature Test Facility

The development of the Modular High-Temperature Gas-Cooled Reactor is a significant milestone in advanced nuclear reactor technology. One of the concerns for the reactor’s safe operation is the effects of a loss-of-flow accident (LOFA) where the coolant circulators are tripped, and forced coolant flow through the core is lost. Depending on the steam generator placement, loop or intracore natural circulation develops to help transfer heat from the core to the reactor cavity, cooling system. This paper investigates the fundamental physical phenomena associated with intracore coolant natural circulation flow in a one-sixth Computational Fluid Dynamics (CFD) model of the Oregon State University High Temperature Test Facility (OSU HTTF) following a loss-of-flow accident transient. This study employs conjugate heat transfer and steady-state flow along with an SST k-ω turbulence model to characterize the phenomenon of core channel-to-channel natural convection. Previous studies have revealed the importance of complex flow distribution in the inlet and outlet plenums with the potential to generate hot coolant jets. For this reason, complete upper and lower plenum volumes are included in the analyzed computational domain. CFD results also include parametric studies performed for a mesh sensitivity analysis, generated using the STAR-CCM+ software. The resulting channel axial velocities and flow directions support the test facility scaling analysis and similarity group distortions calculation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Computational and Experimental Investigation of Discrete “W” Turbulators for sCO2 Turbine Internal Cooling

Direct-fired supercritical CO2 (sCO2 ) cycles are of particular interest for power generation due to the potential for higher thermal efficiencies with inherent CO2 capture. Direct fired sCO2 turbines must withstand high temperatures and the density of the working fluid is much higher than conventional Brayton cycles. Based on conventional Brayton cycle applications, internal serpentine cooling channels with ribbed turbulator strips are a common approach to achieve the required cooling, while minimizing coolant utilization. Most of the previous studies of coolant channel heat transfer enhancements have focused on serpentine channels using air as the cooling media at Reynold’s numbers typical of Brayton cycles. This paper investigates the Nusselt number augmentation of coolant channels with no features, ribbed turbulators, and discrete “W” turbulators with surface roughness effects included and CO2 as the coolant. The numerical results are compared to experimental results using CO2 as the cooling media at Reynold’s numbers ranging from 100,000 to 300,000. The results show the Nusselt number augmentation plateaus with increasing Reynolds number, in contrast to the decay observed in prior publications with air as the coolant media. Additionally, the CFD results indicate the discrete “W” geometry provides a 10% increase in Nusselt number augmentation relative to the angled ribbed turbulators.

computational fluid dynamics (CFD)↗
DOE OSTI2022

Ice-shelf ocean boundary layer dynamics from large-eddy simulations

Abstract. Small-scale turbulent flow below ice shelves is regionally isolated and difficult to measure and simulate. Yet these small-scale processes, which regulate heat and salt transfer between the ocean and ice shelves, can affect sea-level rise by altering the ability of Antarctic ice shelves to “buttress” ice flux to the ocean. In this study, we improve our understanding of turbulence below ice shelves by means of large-eddy simulations at sub-meter resolution, capturing boundary layer mixing at scales intermediate between laboratory experiments or direct numerical simulations and regional or global ocean circulation models. Our simulations feature the development of an ice-shelf ocean boundary layer through dynamic ice melting in a regime with low thermal driving, low ice-shelf basal slope, and strong shear driven by the geostrophic flow. We present a preliminary assessment of existing ice-shelf basal melt parameterizations adopted in single component or coupled ice-sheet and ocean models on the basis of a small parameter study. While the parameterized linear relationship between ice-shelf melt rate and far-field ocean temperature appears to be robust, we point out a little-considered relationship between ice-shelf basal slope and melting worthy of further study.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Anisotropic Heating and Parallel Heat Flux in Electron-only Magnetic Reconnection with Intense Guide Fields

Electron-only reconnection (E-REC) is a process recently observed in the Earth’s magnetosheath, where magnetic reconnection occurs at electron kinetic scales, and ions do not couple to the reconnection process. Electron-only reconnection is likely to have a significant impact on the energy conversion and dissipation of turbulence cascades at kinetic scales in some settings. This paper investigates E-REC under different intensities of strong guide fields (the ratio between the guide field and the in-plane asymptotic field strength is 5, 10 and 20, respectively) via two-dimensional fully kinetic particle-in-cell simulations, focusing on electron heating. The simulations are initialized with a force-free current sheet equilibrium under various intensities of strong guide fields. Similarly to previous experimental studies, electron temperature anisotropy along separatrices is observed, which is found to be mainly caused by the variations of parallel temperature. Both regions of anisotropy and parallel temperature increase/decrease along separatrices become thinner with increasing guide fields. Besides, we find a transition from a quadrupolar to a hexapolar (six-polar) to an octopolar (eight-polar) structure in temperature anisotropy and parallel temperature as the guide field intensifies. Non-Maxwellian electron velocity distribution functions (EVDFs) at different locations in the three simulations are observed. Our results show that parallel electron velocity varies notably with different guide field intensities and finite parallel electron heat flux density is observed. The three simulations exhibit features of the Chew–Goldberger–Low theory, with the level of consistency increasing as the guide field strength increases. This explains the electron parallel temperature variations and the shape of the EVDFs observed along the separatrices. This work may provide insights into the understanding of electron heating and parallel heat flux density in E-REC observed in the turbulent magnetosheath.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Ion-temperature- and density-gradient-driven instabilities and turbulence in Wendelstein 7-X close to the stability threshold

Electrostatic gyrokinetic instabilities and turbulence in the Wendelstein 7-X stellarator are studied. Particular attention is paid to the ion-temperature-gradient (ITG) instability and its character close to marginal stability [Floquet-type turbulence (Zocco et al., Phys. Rev. E, vol. 106, 2022, p. L013202) with no electron temperature gradient]. The flux tube version of the δƒ code stella (Barnes et al., J. Comput. Phys., vol. 391, 2019, pp. 365–380) is used to run linear and nonlinear gyrokinetic simulations with kinetic electrons. The nature of the dominant instability depends on the wavelength perpendicular to the magnetic field, and the results are conveniently displayed in stability diagrams that take this dependence into account. This approach highlights the presence of universal instabilities, which are less unstable but have longer wavelengths than other modes. A quasi-linear estimate of the heat flux suggests they are relevant for transport. Close to the stability threshold, the linear eigenmodes and turbulence form highly extended structures along the computational domain if the magnetic shear is small. Numerical experiments and diagnostics are undertaken to assess the resulting radial localisation of the turbulence, which affects the interaction of the latter with zonal flows. Increasing the amplitude of the magnetic shear (e.g.through current drive) has a stabilising effect on the turbulence and, thus, reduces the nonlinear energy transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Numerical Simulations of Orographic Convection across Multiple Gray Zones

Idealized simulations are used to determine the sensitivity of moist orographic convection to horizontal grid spacing Δ h . In simulated mechanically (MECH) and thermally (THERM) forced convection over an isolated ridge, Δ h is varied systematically over both the deep-convection (Δ h ~ 10 –1 km) and turbulence (Δ h ~ 1 km–100 m) gray zones. To aid physical interpretation, a new parcel-based bulk entrainment/detrainment diagnosis for horizontally heterogeneous flows is developed. Within the deep-convection gray zone, the Δ h sensitivity is dominated by differences in parameterized versus explicit convection; the former initiates convection too far upstream of the ridge (MECH) and too early in the diurnal heating cycle (THERM). These errors stem in part from a large underprediction of parameterized entrainment and detrainment. Within the turbulence gray zone, sensitivities to Δh arise from the representation of both subcloud- and cloud-layer turbulence. As Δ h is decreased, MECH exhibits stronger cloud-layer entrainment to enhance the convective mass flux M co , while THERM exhibits stronger detrainment to suppress Mco and delay convection initiation. The latter is reinforced by increased subcloud turbulence at smaller Δ h , which leads to drying and diffusion of the central updraft responsible for initiating moist convection. Numerical convergence to a robust solution occurs only in THERM, which develops a fully turbulent flow with a resolved inertial subrange (for Δ h ≤ 250 m). In MECH, by contrast, turbulent transition occurs within the orographic cloud, the details of which depend on both physical location and Δ h .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Experimental Study of Alfvén Wave Reflection from an Alfvén-speed Gradient Relevant to the Solar Coronal Holes

Abstract We report the first experimental detection of a reflected Alfvén wave from an Alfvén-speed gradient under conditions similar to those in coronal holes. The experiments were conducted in the Large Plasma Device at the University of California, Los Angeles. We present the experimentally measured dependence of the coefficient of reflection versus the wave inhomogeneity parameter, i.e., the ratio of the wavelength of the incident wave to the length scale of the gradient. Two-fluid simulations using the Gkeyll code qualitatively agree with and support the experimental findings. Our experimental results support models of wave heating that rely on wave reflection at low heights from a smooth Alfvén-speed gradient to drive turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Assessment of turbulent combustion models for simulating pre-chamber ignition in a natural gas engine

Lean combustion in an internal combustion engine is a promising strategy to increase thermal efficiency by leveraging a more favorable specific heat ratio of the fresh mixture and simultaneously suppressing the heat losses to the cylinder wall. However, unstable ignition events and slow flame propagation at fuel-lean condition lead to high cycle-to-cycle variability and hence limit the high-efficiency engine operating range. Pre-chamber ignition is considered an effective concept to extend the lean operating limit, by providing spatially distributed ignition with multiple turbulent flame-jets and enabling faster combustion rate compared to the conventional spark ignition approach. From a numerical modeling perspective, to date, still the science base and available simulation tools are inadequate for understanding and predicting the combustion processes in pre-chamber ignited engines. In this paper, conceptually different RANS combustion models widely adopted in the engine modeling community were used to simulate the ignition and combustion processes in a medium-duty natural gas engine with a pre-chamber spark-ignition system. A flamelet-based turbulent combustion model, i.e., G-equation, and a multi-zone well-stirred reactor model were employed for the multi-dimensional study. Simulation results were compared with experimental data in terms of in-cylinder pressure and heat release rate. Finally, the analysis of the performance of the two models is carried out to highlight the strengths and limitations of the two formulations respectively.

combustion model evaluation↗
DOE OSTI2020

On the Parameterization of Convective Downdrafts for Marine Stratocumulus Clouds

The role of nonlocal transport on the development and maintenance of marine stratocumulus (Sc) clouds in coarse-resolution models is investigated, with a special emphasis on the downdraft contribution. A new parameterization of cloud-top-triggered downdrafts is proposed and validated against large-eddy simulation (LES) for two Sc cases. The applied nonlocal mass-flux scheme is part of the stochastic multiplume eddy-diffusivity/mass-flux (EDMF) framework decomposing the turbulent transport into local and nonlocal contributions. The complementary local turbulent transport is represented with the Mellor–Yamada–Nakanishi–Niino (MYNN) scheme. This EDMF version has been implemented in the Weather Research and Forecasting (WRF) single-column model (SCM) and tested for three model versions: without mass flux, with updrafts only, and with both updrafts and downdrafts. In the LES, the downdraft and updraft contributions to the total heat and moisture transport are comparable and significant. The WRF SCM results show a good agreement between the parameterized downdraft turbulent transport and LES. While including updrafts greatly improves the modeling of Sc clouds over the simulation without mass flux, the addition of downdrafts is less significant, although it helps improve the moisture profile in the planetary boundary layer.

Wu, Elynn↗
DOE OSTI2023

Near-wall model for compressible turbulent boundary layers based on an inverse velocity transformation

In this work, a near-wall model, which couples the inverse of a recently developed compressible velocity transformation (Griffin et al., Proc. Natl Acad. Sci., vol. 118, 2021, p. 34) and an algebraic temperature–velocity relation, is developed for high-speed turbulent boundary layers. As input, the model requires the mean flow state at one wall-normal height in the inner layer of the boundary layer and at the boundary-layer edge. As output, the model can predict mean temperature and velocity profiles across the entire inner layer, as well as the wall shear stress and heat flux. The model is tested in an a priori sense using a wide database of direct numerical simulation high-Mach-number turbulent channel flows, pipe flows and boundary layers (48 cases, with edge Mach numbers in the range 0.77–11, and semi-local friction Reynolds numbers in the range 170–5700). The present model is significantly more accurate than the classical ordinary differential equation (ODE) model for all cases tested. The model is deployed as a wall model for large-eddy simulations in channel flows with bulk Mach numbers in the range 0.7–4 and friction Reynolds numbers in the range 320–1800. When compared to the classical framework, in the a posteriori sense, the present method greatly improves the predicted heat flux, wall stress, and temperature and velocity profiles, especially in cases with strong heat transfer. In addition, the present model solves one ODE instead of two, and has a computational cost and implementation complexity similar to that of the commonly used ODE model.

42 ENGINEERING↗
DOE OSTI2020

Numerical Investigation of Fuel Property Effects on Mixed-Mode Combustion in a Spark-Ignition Engine

In this research, lean mixed-mode combustion is numerically investigated using computational fluid dynamics (CFD) in a spark-ignition engine. A new E30 fuel surrogate is developed using a neural network model with matched octane numbers. A skeletal mechanism is also developed by automated mechanism reduction and by incorporating a NO x submechanism. A hybrid approach that couples the G-equation model and the well-stirred reactor model is employed for turbulent combustion modeling. The developed CFD model is shown to well predict pressure and apparent heat release rate (AHRR) traces compared with experiment. Two types of combustion cycles (deflagration-only and mixed-mode cycles) are observed. The mixed-mode cycles feature early flame propagation and subsequent end-gas auto-ignition, leading to two distinctive AHRR peaks. The validated CFD model is then employed to investigate the effects of NO x chemistry. The NO x chemistry is found to promote auto-ignition through the residual gas, while the deflagration phase remains largely unaffected. Sensitivity analysis is finally performed to understand effects of fuel properties, including heat of vaporization (HoV) and laminar flame speed (S L ). An increased HoV tends to suppress auto-ignition through charge cooling, while the impact of HoV on flame propagation is insignificant. In contrast, an increased S L is found to significantly promote both flame propagation and end-gas auto-ignition. The promoting effect of S L on auto-ignition is not a direct chemical effect; it is rather caused by an advancement of the combustion phasing, which increases compression heating of the end-gas.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Investigation of Gas Turbine Internal Cooling Using Supercritical CO 2 —Effect of Surface Roughness and Channel Aspect Ratio

In this paper, an experimental and numerical investigation of internal cooling channels with rib turbulators is presented with sCO 2 as the working fluid at process conditions (pressure-20.7 MPa and temperature up to 150 °C). The effect of channel aspect ratio up to 2:1 on thermal-hydraulic performance is explored in additively manufactured rectangular channels and square channels, both with and without 60 deg ribs on the top and bottom sides. The Wilson-plot method is employed to experimentally measure channel-averaged Nusselt number over a Reynolds number range up to 370,000. The friction factor is calculated from pressure drop and mass flow rate and additionally, the overall thermal performance factor (TPF) is reported. A companion computational fluid dynamics (CFD) simulation is performed for the rib turbulated cooling configurations reported in the experiments using the Reynolds average Navier–Stokes-based turbulence model. The objective of the numerical study is to gain insight into the local heat transfer augmentation in the ribbed channels as a result of varying the aspect ratio, channel configuration (square versus rectangular), operating conditions (Reynolds number) and the surface roughness, an inherent outcome of the additive manufacturing process. Surface roughness is simulated using sand grain roughness height (K S ) calculated from the experimental data, and a comparison is presented with the corresponding channel configuration with varying surface roughness heights starting from smooth surfaces (K S = 0). Experimental results indicate that the heat transfer augmentation is negligible in the rectangular channels with ribs on the long side compared to the square channel. However, it is enhanced by 60% in comparison to placing ribs on the shorter side. The TPF remains constant at around 1 for the entire range of Reynolds numbers consistent with prior work at the National Energy Technology Laboratory (NETL). Here, the simulation results highlight that increased surface roughness can have a favorable considerable influence on Nusselt number and overall thermal performance enhancement.

42 ENGINEERING↗
DOE OSTI2022

Examination of stiff ion temperature gradient mode physics in simulations of DIII-D H-mode transport

A systematic evaluation of gyrokinetic and gyrofluid model predictions of ion temperature gradient (ITG) stability and transport using parameters from DIII-D high confinement mode (H-mode) plasmas has been performed. The nonlinear CGYRO code is used to make the gyrokinetic predictions, and the quasilinear TGLF model for the corresponding gyrofluid predictions. The assessments are made at three radii (normalized toroidal flux ρtor = 0.4, 0.55, and 0.7) in three different plasma scenarios with varying levels of neutral beam heating and torque. For each of the nine cases (3 radii × 3 scenarios) considered, ITG turbulence is found to be the dominant long-wavelength instability and transport mechanism. The inclusions of both transverse magnetic fluctuations and dynamic fast beam ions are stabilizing for all cases considered, with strongest effects seen at ρor = 0.4 where the fast ion population and normalized plasma pressure β = 2μ0nT/B2 are highest. The further inclusion of parallel magnetic fluctuations does not have a meaningful impact on the ITG turbulence in these scenarios, but does destabilize (in combination with fast ions) new high-frequency instabilities at ρtor = 0.4 in the high power scenarios. In each case the linear and nonlinear ITG critical gradients are predicted to be lower than the measured ITG scale lengths and their associated uncertainties. Inclusion of equilibrium flow shear in the transport predictions generally leads to an upshift in effective critical gradient rather than a qualitative change in the predicted stiffness, with stronger responses typically seen in the gyrokinetic predictions than in the gyrofluid results. However, in most cases these upshifted gradients still remain below the measured values and their uncertainties. Although the predicted critical gradients are below the measured gradients, both models predicted flux-matching gradients consistent with measured values in six of the nine cases considered, with no clear systematic over- or underprediction. Thus, while the experimental ion temperature profiles do not appear to be closely pinned to the ITG critical gradient, both gyrokinetic and gyrofluid models are able to accurately match the measured gradients reasonably well in most cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Capturing Secondary Kinetic Instabilities in Three‐Dimensional Dayside Reconnection Using an Improved Gradient‐Based Closure

Magnetic reconnection is a highly dynamic process that excites a wide variety of kinetic waves and instabilities. Transverse current sheet instabilities such as the lower-hybrid drift and secondary drift-kink instabilities in particular have been shown by kinetic simulations to modify the reconnection and introduce significant turbulence and mixing to the reconnection layer. Past studies using the ten-moment fluid model to capture important kinetic physics such as the electron inertia and full representation of the pressure tensor proved advantageous to a two-fluid representation of reconnection, but the model struggled when using a local relaxation closure for the heat flux to replicate the current sheet instabilities and subsequent mixing seen in kinetic simulations. This work uses the Gkeyll software framework to perform simulations of asymmetric reconnection based on the 16 October 2015 MMS crossing of a diffusion region, the Burch event. An improved gradient-based heat flux closure is implemented, showing significant improvement in secondary kinetic instabilities that grow in the current sheet. These instabilities generate turbulence which leads to growth of secondary magnetic islands and flux ropes.

Bradshaw, K. [Princeton University, NJ (United Sta↗
DOE OSTI2023

Use of intrinsic hysteresis for the active control of internal transport barriers in magnetically confined fusion plasmas

In magnetically confined fusion devices, control of internal transport barriers (ITBs) is important both to enhance and suppress the turbulent transport to improve confinement control. Barrier control should allow for the improvement of confinement to aid in achieving the needed fusion criteria while also permitting the degradation of confinement to control profiles and clean the device by moving out the impurities accumulated near the core. In this work, we present a novel control scenario that takes advantage of the hysteresis intrinsic to transport barriers to easily cycle through enhanced and degraded confinement regimes. The control scenario is illustrated using a five-field simplified transport model for an ITB using typical parameters of a neutral beam injection-heated DIII-D tokamak discharge. Pellets and ion cyclotron resonance frequency power are used as control knobs for this active control scenario. These knobs adequately modify at will the local gradients and, therefore, the growth rates and shearing rates, allowing for easy and efficient control of the barrier by taking advantage of the barrier hysteresis. The result is a control cycle that could be operated with a relatively small amount of power in high performance regimes which, nowadays, typically require large power to control. It may also have advantages to avoid, or at least ameliorate, the appearance of magnetohydrodynamic instabilities in the barrier region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗