Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “REYNOLDS NUMBER”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Shape effects on the local dynamics of suspensions of spheroidal particles

The effect of shape on the dynamics of suspensions of non-spherical heavy particles is examined by fully resolved numerical simulations of oblate and prolate spheroids, as well as spheres, for a density ratio of ten, volume fractions ranging from 0.5% to 5%, and Reynolds numbers between 20 and 30. The dynamics is determined both by the interactions of the particles with the fluid as well as by collisions, with the number and importance of collisions increasing with volume fractions. A single isolated oblate or prolate spheroid falling under gravity generally falls broadside on, for the governing parameters examined here, and at low-volume fractions, the majority of particles in a suspension fall that way. At higher-volume fractions, the orientation is more random. The slip velocity decreases as the volume fraction increases for all shapes, as expected, but the effect of the shape is much less than seen for a single particle. Furthermore, this seems to be due to two effects. For all volume fractions, the most deformed particles cluster more than spheres and less deformed particles, which increases their slip velocity. As the concentration increases, the increased particle interactions also causes more particles to fall short side-on, which reduces the frontal area and the resulting drag, increasing the slip velocity. This second effect is, of course, absent for spherical particles.

42 ENGINEERING↗
DOE OSTI2025

FreeMHD: Validation and verification of the open-source, multi-domain, multi-phase solver for electrically conductive flows

The extreme heat fluxes in the divertor region of tokamaks may require an alternative to solid plasma-facing components, for the extraction of heat and the protection of the surrounding walls. Flowing liquid metals are proposed as an alternative, but raise additional challenges that require investigation and numerical simulations. Free surface designs are desirable for plasma-facing components, but steady flow profiles and surface stability must be ensured to limit undesirable interactions with the plasma. Previous studies have mainly used steady-state, 2D, or simplified models for internal flows and have not been able to adequately model free-surface liquid metal (LM) experiments. Therefore, FreeMHD has been recently developed as an open-source magnetohydrodynamics (MHD) solver for free-surface electrically conductive flows subject to a strong external magnetic field. The FreeMHD solver computes incompressible free-surface flows with multi-region coupling for the investigation of MHD phenomena involving fluid and solid domains. The model utilizes the finite-volume OpenFOAM framework under the low magnetic Reynolds number approximation. FreeMHD is validated using analytical solutions for the velocity profiles of closed channel flows with various Hartmann numbers and wall conductance ratios. Next, experimental measurements are then used to verify FreeMHD, through a series of cases involving dam breaking, 3D magnetic fields, and free-surface LM flows. These results demonstrate that FreeMHD is a reliable tool for the design of LM systems under free surface conditions at the reactor scale. Furthermore, it is flexible, computationally inexpensive, and can be used to solve fully 3D transient MHD flows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

On the accuracy of compressibility transformations

This study highlights the importance of satisfying the eddy viscosity equivalence below the logarithmic layer, to deriving accurate compressibility transformations. First, we analyze the ability of known transformations to satisfy the eddy viscosity equivalence and show that the accuracy of these transformations is strongly dependent on this ability. Second, in a step-by-step manner, we devise new transformations that satisfy this hypothesis. An approach based on curve fitting of the incompressible Direct Numerical Simulation data for eddy viscosity profiles below the logarithmic layer provides an extremely accurate transformation, which motivates self-contained methods, making use of mixing length formulas in the inner region. It is shown that the accuracy of existing transformations can be significantly improved by applying these ideas, below the logarithmic layer. Motivated by the effectiveness of the formulations derived from eddy viscosity equivalence, we introduce a new integral transformation based on Reynolds number equivalence between compressible and incompressible flows. This approach is based on defining a new compressible velocity scale, which affects the accuracy of transformations. Several choices for the velocity scale are tested, and in each attempt, it is shown that the eddy viscosity equivalence plays a very important role for the accuracy of compressibility transformations.

42 ENGINEERING↗
DOE OSTI2024

An enrichment wall modeling framework for spectral element methods

In the present work, a first-of-its-kind enrichment wall-model is developed within the spectral element method (SEM) framework for large-eddy simulations (LES) of wall-bounded turbulent flows. The method augments the polynomial solution in the wall-adjacent elements with an analytical law-of-the-wall enrichment function representing the mean velocity near the wall. In the solution representation, this enrichment function captures the large gradients in the boundary layer, which allows the polynomial modes to represent the turbulent fluctuations. The enriched solution is able to resolve the shear stress at the wall without any modification to the no-slip wall boundary conditions, which allows for greater accuracy in the near-wall region compared to traditional methods. The enrichment wall modeling approach is implemented in a high-order SEM computational fluid dynamics solver, Nek5000, and its performance is assessed in turbulent channel flow wall-modeled LES for a range of Reynolds numbers. It is demonstrated that the enrichment wall-model improves solution accuracy on under-resolved near-wall grids as compared to traditional shear stress wall-models.

42 ENGINEERING↗
DOE OSTI2023

Lagrangian large eddy simulations via physics-informed machine learning

High-Reynolds number homogeneous isotropic turbulence (HIT) is fully described within the Navier–Stokes (NS) equations, which are notoriously difficult to solve numerically. Engineers, interested primarily in describing turbulence at a reduced range of resolved scales, have designed heuristics, known as large eddy simulation (LES). LES is described in terms of the temporally evolving Eulerian velocity field defined over a spatial grid with the mean-spacing correspondent to the resolved scale. This classic Eulerian LES depends on assumptions about effects of subgrid scales on the resolved scales. Here, we take an alternative approach and design LES heuristics stated in terms of Lagrangian particles moving with the flow. Our Lagrangian LES, thus L-LES, is described by equations generalizing the weakly compressible smoothed particle hydrodynamics formulation with extended parametric and functional freedom, which is then resolved via Machine Learning training on Lagrangian data from direct numerical simulations of the NS equations. The L-LES model includes physics-informed parameterization and functional form, by combining physics-based parameters and physics-inspired Neural Networks to describe the evolution of turbulence within the resolved range of scales. The subgrid-scale contributions are modeled separately with physical constraints to account for the effects from unresolved scales. We build the resulting model under the differentiable programming framework to facilitate efficient training. We experiment with loss functions of different types, including physics-informed ones accounting for statistics of Lagrangian particles. We show that our L-LES model is capable of reproducing Eulerian and unique Lagrangian turbulence structures and statistics over a range of turbulent Mach numbers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Inertia-induced mixing and reaction maximization in laminar porous media flows

Solute transport and biogeochemical reactions in porous and fractured media flows are controlled by mixing, as are subsurface engineering operations such as contaminant remediation, geothermal energy production, and carbon sequestration. Porous media flows are generally regarded as slow, so the effects of fluid inertia on mixing and reaction are typically ignored. Here, we demonstrate through microfluidic experiments and numerical simulations of mixing-induced reaction that inertial recirculating flows readily emerge in laminar porous media flows and dramatically alter mixing and reaction dynamics. An optimal Reynolds number that maximizes the reaction rate is observed for individual pore throats of different sizes. This reaction maximization is attributed to the effects of recirculation flows on reactant availability, mixing, and reaction completion, which depend on the topology of recirculation relative to the boundary of the reactants or mixing interface. Recirculation enhances mixing and reactant availability, but a further increase in flow velocity reduces the residence time in recirculation, leading to a decrease in reaction rate. The reaction maximization is also confirmed in a flow channel with grain inclusions and randomized porous media. Interestingly, the domain-wide reaction rate shows a dramatic increase with increasing Re in the randomized porous media case. This is because fluid inertia induces complex three-dimensional flows in randomized porous media, which significantly increases transverse spreading and mixing. This study shows how inertial flows control reaction dynamics at the pore scale and beyond, thus having major implications for a wide range of environmental systems.

Chen, Michael A. (ORCID:0000000294614710)↗
DOE OSTI2021

Pressure drop in a prototypical 3D magnetohydrodynamic flow across contraction of a fusion blanket manifold

Predicting 3D magnetohydrodynamic (MHD) pressure losses associated with liquid metal flows in complex geometry ducts is required to design breeding blankets for fusion reactors. Of such components, manifolds exhibit major pressure losses. Recent correlations [T. Rhodeset al. Magnetohydrodynamic pressure drop and flow balancing of liquid metal flow in a prototypic fusion blanket manifold. Phys. Fluids. 2018; 30: 057101.] demonstrate promise for predicting pressure drops in electrically insulated inlet manifolds. In the present work, the extend to which these correlations can be applied to outlet manifolds, which feature sudden contractions, is investigated in both viscous-electromagnetic (VE) and inertial-electromagnetic (IE) regimes for Reynolds numbers 50 < Re < 1500, Hartmann numbers 2500 < Ha < 5475, and expansion/contraction ratio 10. Numerical computations have shown that the MHD pressure drops in the eutectic lead-lithium (PbLi) flows in contractions are almost identical to those in the flows featuring expansion with slightly higher magnitudes in the contraction cases. The small discripancy in the MHD pressure drop between contractions and expansions (<8%) suggests that the earlier obtained correlations for the 3D MHD pressure drop in a duct flow with a sudden expansion can also be applied to flows with a sudden contraction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Toward Exascale: Overview of Large Eddy Simulations and Direct Numerical Simulations of Nuclear Reactor Flows with the Spectral Element Method in Nek5000

At the beginning of the last decade, Petascale supercomputers (i.e., computers capable of more than 1 petaFLOP) emerged. Now, at the dawn of exascale supercomputing, we provide a review of recent landmark simulations of portions of reactor components with turbulence-resolving techniques that this computational power has made possible. In fact, these simulations have provided invaluable insight into flow dynamics, which is difficult or often impossible to obtain with experiments alone. We focus on simulations performed with the spectral element method, as this method has emerged as a powerful tool to deliver massively parallel calculations at high fidelity by using large eddy simulation or direct numerical simulation. We also limit this paper to constant-property incompressible flow of a Newtonian fluid in the absence of other body or external forces, although the method is by no means limited to this class of flows. We briefly review the fundamentals of the method and the reasons it is compelling for the simulation of nuclear engineering flows. We review in detail a series of Petascale simulations, including the simulations of helical coil steam generators, fuel assemblies, and pebble beds. Even with Petascale computing, however, limitations for nuclear modeling and simulation tools remain. In particular, the size and scope of turbulence-resolving simulations are still limited by computing power and resolution requirements, which scale with the Reynolds number. In the final part of this paper, we discuss the future of the field, including recent advancements in emerging architectures such as GPUbased supercomputers, which are expected to power the next generation of high-performance computers.

computational fluid dynamics↗
DOE OSTI2025

High-Fidelity CFD Simulation of Mixed Convection and Forced Convection in a Pebble Bed Test Reactor Core

The Hermes low-power [35-MW(thermal)] reactor will be built and operated by Kairos Power LLC (KP) to demonstrate its fluoride salt-cooled high-temperature reactor (FHR) technology. In the KP FHR, the reactor core is composed of randomly packed pebbles with TRISO fuel particles inside with FLiBe flow upward through the core acting as a coolant. Previous numerical and experimental studies have been limited to either a small-size bed or to a lack of detailed measurements for heat transfer. Here, to address the lack of high-fidelity heat transfer data in a real-size FHR core, in this study, we simulated a pebble bed core with 34 374 pebbles randomly packed, similar to the Hermes reactor's size. The core radius was 14 times that of the pebble diameter, while the core height was 45 times. In this work, we were particularly interested in a mixed convection regime, where buoyancy is important. Therefore, we performed several large-eddy simulations at different Reynolds numbers (160 to 1000) with gravitational force included. The spectral element computational fluid dynamics code NekRS with graphics processing unit acceleration was used for this study. The low-Mach number approximation was applied to address property changes in the FLiBe and to account for buoyancy. A pure hexahedral mesh with 60 million elements was generated by the Voronoi cell method. At the polynomial order of 5, the total degrees of freedom was 7.5 billion. The developed case in this work is the first of its kind in terms of size and complexity. The local numerical data across the domain were obtained and compared with empirical correlations. After examining the data, we found the following conclusions. For pressure drop, the Reger correlation predicted less than a 5% error. On the other hand, for heat transfer, the Wakao correlation outperformed the others. Based on our findings, we recommend the use of the Wakao correlation for the Nusselt number calculation, and for pressure drop, the KTA (Kerntechnischer Ausschuss) correclation, among the available experimental correlations. In conclusion, the Reger direct numerical simulation-driven correlation for pressure drops should also be considered, given its best agreement with our calculations.

Mixed Convection↗
DOE OSTI2025

Dynamic Behavior of Oval-Twisted Helical Tube Heat Exchanger: Numerical Study with RELAP5-3D

Convective heat transfer characteristics and theoretical thermal stress behaviors are numerically calculated using RELAP5-3D for the helical-coiled once-through steam generator (H-OTSG) and the novel heat exchanger design known as the oval-twisted helically coiled heat exchanger (OTHCHX) under (1) fluctuating wall temperature conditions, (2) square-wave pulsating flow conditions, and (3) the combined effects of fluctuating wall temperature and square-wave pulsating flow conditions. Heat transfer coefficient models for the H-OTSG and OTHCHX were developed based on existing data and implemented into RELAP5-3D, successfully capturing the N⁢uavg behavior within 8% to 10% of the reported data. Under fluctuating wall temperature conditions, the OTHCHX displayed higher N⁢u avg behavior than the H-OTSG. As 𝑓 increased, the $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ decreased. The $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ was higher for the OTHCHX than for the H-OTSG under fluctuating wall temperature conditions. Under pulsating flow conditions, the H-OTSG and OTHCHX displayed much higher 𝑁⁢𝑢 𝑎𝑣𝑔 than under constant flow conditions. The H-OTSG displayed a higher $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ over the OTHCHX. Under combined fluctuating wall temperature and pulsating flow conditions, the augmented heat transfer behavior from the pulsating flow was counteracted by the wall temperature fluctuations, producing slightly higher 𝑁⁢𝑢 𝑎𝑣𝑔 over constant wall temperature, constant flow conditions, but much lower than only constant pulsating flow under constant wall temperature conditions. The effects of simultaneous wall temperature fluctuations and square-wave pulsating flow caused higher $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ than that of only wall temperature fluctuations or pulsating flow. As the Reynolds number (Re) increased, $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ increased. However, when 𝑓=𝑓$_{\dot{m}}$, the $𝜎^{𝑚⁢𝑎⁢𝑥}_{𝑡⁢ℎ}$ showed decreasing values as Re increased. In conclusion, the results indicate that thermal-fluid resonance can help mitigate thermal stresses.

Thermal stress↗
DOE OSTI2026

Modeling Flow in a Microreactor Core: From High-Fidelity CFD to Subchannel Analysis

The accurate prediction of pressure drop in tightly packed, low–Reynolds number (Re) bare rod bundles is essential for the thermal-hydraulic design of the Microreactor Applications Research Validation and EvaLuation (MARVEL) reactor and other microreactor concepts. However, existing friction factor correlations, particularly the upgraded Cheng-Todreas (UCTD) correlation, have limited validation for the small pitch-to-diameter ratios (P/Ds) and transitional flow conditions characteristic of these systems. In this work, we perform high-fidelity large-eddy simulations (LESs) of both an infinite bare rod bundle and a finite bare 37-pin scalloped rod bundle across the range 1000≤ 𝑅𝑒 ≤5500. The simulations reveal strong gap vortex–driven transitional behavior and indicate that the UCTD may underpredict the friction factor by up to 28% at 𝑃/𝐷 = 1.05. Using the LES-calculated pressure drops, we formulate a new friction factor correlation that follows the Cheng-Todreas functional form but is calibrated for low-Re and tightly packed geometries representative of a MARVEL-like reactor. The correlation is implemented in the MOOSE (Multiphysics Object-Oriented Simulation Environment) subchannel module and compared against both the LES and UCTD predictions. Across all subchannel types, the proposed model reduces the streamwise velocity differences from as high as 44% (UCTD) to below 9%, and decreases the pressure gradient differences from 13% to 25% (UCTD) to 0.7% to 7% relative to the LES results. These results suggest that the new correlation has the potential to improve the pressure drop and flow field predictions for such geometries, highlighting the importance of high-fidelity simulations in supporting microreactor thermal-hydraulic model development and motivating future pressure drop experiments for compact rod bundles to further validate these findings.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2022

Pressure Drop in Seven-Pin Wire-Wrapped Rod Bundle for the Sodium Cartridge Loop in Versatile Test Reactor

This work studies the hydrodynamics of the seven pin wire-wrapped rod bundle in the sodium cartridge loop for the Versatile Test Reactor (VTR) through scaled water experiments and computational fluid dynamics (CFD) simulations. Here, the scaling analysis is first performed to demonstrate the hydrodynamic similarity between water and sodium flows at the same Reynolds number Re. A separate-effects test facility is designed and constructed based on the scaling analysis. Detailed experimental data on the pressure drop covering a wide range of Re values (1165 to 27 689) are obtained, which are used to evaluate existing correlations for friction factor and to benchmark CFD simulations. The experimentally determined friction factors f Exp agree well with the Upgraded Cheng and Todreas Detailed Correlation and Pacio-Chen-Todreas Detailed Model within ±7% but are significantly underpredicted by Rehme's correlation by 25%. Various CFD near-wall treatment methods are tested using ANSYS Fluent and evaluated by experimental data. It is found that when the recommended wall y + values are met, most of the near-wall treatment methods can give accurate friction factor predictions. The resolved near-wall method (y + ~ 1) with the Shear Stress Transport k - ω turbulence model and the scalable wall functions (y + > 12) with the realizable k - ε turbulence model can predict f Exp within ±5% The standard wall functions (y + > 30) and nonequilibrium wall functions (y + > 30) with the realizable k - ε model can predict f Exp within ±10%.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Experimental and Numerical Evaluation of the Thermal-Fluid Characteristics of the Helium-Cooled Modular Divertor with Multiple Jets Using a Reversed Heat Flux Approach

A single-finger unit of the Helium-Cooled Modular Divertor with Multiple Jets (HEMJ) with a plasma-facing surface (PFS) area of about 2 cm 2 has been studied in a helium (He) loop at He mass flow rates m˙ H ≤ 8 g/s and nearly prototypical conditions. Based on previous studies of the single finger of the HEMJ, our Georgia Institute of Technology group is planning to experimentally study larger divertors. Given that the HEMJ test section was heated with an induction heater and that it is impractical to scale this up to divertors with larger PFS areas, a reversed heat flux approach is being considered to measure heat transfer coefficients (HTCs). In this approach, the direction of the heat flux is reversed with water cooling and high-temperature He heating of the outer shell attached to the PFS. This work presents an initial experimental and numerical evaluation of this approach for a single HEMJ finger. Experiments with brass and copper-chromium-zirconium outer shells were conducted at dimensionless He mass flow rates or Reynolds numbers Re = 1 × 10 4 to 4.7 × 10 4 , an inlet pressure of 10 MPa, temperatures as great as 673 K, and maximum heat flux of 8.4 MW/m 2 . The experiments verify that the He-side HTCs are independent of the direction of the heat flux. Here, the results agree well with previous Nusselt number correlation and pressure loss coefficients for the HEMJ obtained using the normal heating approach.

42 ENGINEERING↗
DOE OSTI2023

Bayesian Parameter Estimation of the k-ω Shear Stress Transport Model for Accurate Simulations of Impinging-Jet Heat Transfer

Given the lack of fusion-relevant component test facilities, current estimates of the thermo-fluid performance of plasma-facing components are based for the most part on numerical simulations. A major source of uncertainty in these simulations is the semiempirical turbulence (closure) models for the Reynolds stresses appearing in the governing Reynolds-averaged Navier-Stokes equations, which involve a set of constants that depend upon the flow. The objective of this study is to evaluate Bayesian parameter estimation of turbulence closure constants in ANSYS Fluent to model heat transfer in impinging jets. The Bayesian statistical calibration produces a probability distribution for these constants from experimental data; the maximum a posteriori estimates are then taken to be the calibrated constants, or parameters. The turbulence model constants are calibrated using an experimental study of a submerged jet of air impinging on a flat heated surface at Reynolds numbers Re = O(10 4 ) and impingement distance in jet diameters H/d = 2. Numerical predictions using the calibrated model parameters are then compared with those generated using the default constants. Predictions obtained with model parameters calibrated on datasets of two different sizes are compared to evaluate the effect of the number of calibration samples. Lastly, the extrapolative ability of the calibrated model is examined by predictions at a Re beyond the calibration values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Simulations of Heat Transfer Using Tight-Fitting Twisted Tape Inserts for First Wall Cooling in Molten Salt Breeder Blankets

One of the major components in fusion energy systems is the fusion blanket, which has a vacuum vessel to contain the plasma. As part of the fusion blanket/vacuum vessel, the first wall and plasma-facing components require sufficient cooling to prevent material degradation during operation from the superheated plasma. Most fusion blanket concepts involve first wall and divertor coolant channels with heat transfer enhancements (HTEs) that are intended to withstand the incident high heat fluxes of 1 to 5 MW/m 2 . Twisted tape inserts are a proposed HTE that have been investigated previously for first wall cooling and monoblock divertor cooling channels and in other nonfusion heat transfer components. By inserting twisted tapes into straight pipes, the amount of turbulence in the system can be increased at lower Reynolds numbers by swirling the flow. This results in better heat transfer characteristics with marginal increases in frictional pressure losses. In particular, simulations of high-Prandtl-number fluids such as the proposed molten salt FLiBe in twisted tapes, which is prototypic to liquid immersion blankets, have not been previously explored. Here, in this study, we simulate various Prandtl numbers in pipes with twisted tape inserts using large eddy simulations to determine the effects of increasing Prandtl numbers on heat transfer performance. The quantities of particular interest are the Nusselt number and the friction factor, which were recovered using data reduction techniques to determine impacts on heat transfer and pressure losses. This work serves as a starting point for determining the feasibility of twisted tape inserts for liquid immersion blanket concepts.

LES↗
DOE OSTI2020

Field-emission plasma enhancement of H₂–O₂ micro-combustion

In this paper, a detailed three-dimensional computational model is developed to assess the potential of micron-scale field emission dielectric barrier discharge (FE-DBD) plasma actuators in non-premixed microburners with channel heights in the range of L z = 0.25–0.75 mm. Results for H 2 –O 2 combustion with Reynolds number Re = 100 are studied for different configurations of plasma actuator arrays. Without plasma actuation, results agree with corresponding experiments, although in this case the burner typically suffers from incomplete combustion due to thermal quenching, radical quenching, and slow diffusive mixing. For the range of conditions simulated, plasma actuation increases combustion completeness from 25%–70% to 40%–80% by enhancing mixing, generating radicals, and Joule heating. Ionic wind resulting from the FE-DBD plasma disrupts the otherwise simple laminar flow in the burner, accelerating the growth of a mixing layer and enhancing heat transfer to channel walls. The larger, lower temperature flames reduce thermal gradients (and thereby thermal stresses) in the walls. Radical generation due to electron-impact dissociation extends the reaction zone in the diffusion flame and increases combustion completeness. In all cases, the additional heat release with plasma actuation (100–190 W) is significantly larger than the power needed to sustain the plasma (≲0.1 W), suggesting that plasma actuation may provide a significant advantage in energy generation using realistic portable micro-combustion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Novel plasma actuator for mitigation of dynamic stall

A novel plasma actuator, the Linear Counter-flow using a Point Embedded Electrode (LCPEE), is developed for the prevention of dynamic stall for a sinusoidal pitching movement between α = 4° and α = 18°. The LCPEE is implemented on a NACA0012 airfoil and tested at a Reynolds Number of Re c = 2 × 105 and reduced frequency of k = π/16. Prior investigations using a standard linear actuator showed that the exposed electrode introduced perturbations passively which delayed dynamic stall when the actuator was off. For the LCPEE actuator, when turned off, there is least passive delay. When the LCPEE is turned on at St f = 50, the dynamic stall is prevented for the sinusoidal pitching motion of the airfoil. The LCPEE actuator is also tested for the same sinusoidal motion between α = 6° and α = 20°. Four cases are considered for the higher α range of motion: actuator off, actuator on at St f = 50 with a sinusoidal input waveform, actuator on at St f = 50 with a triangular input waveform, and actuator on at Stf = 100 with a sinusoidal input waveform. Specifically for the last case, experiment shows no flow reversal demonstrating the efficacy of LCPEE in controlling the dynamic stall. The effects of LCPEE on the flow energy distribution have also been studied by using proper orthogonal decomposition (POD) method.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Lagrangian particle model for 3D simulation of pellets and SPI fragments in tokamaks

A 3D numerical model for the ablation of pellets and shattered pellet injection (SPI) fragments in tokamaks in the plasma disruption mitigation and fueling parameter space has been developed based on the Lagrangian particle code [R. Samulyak, X. Wang, H.-S. Chen, Lagrangian Particle Method for Compressible Fluid Dynamics, J. Comput. Phys., 362 (2018), 1-19]. The pellet code implements the low magnetic Reynolds number MHD equations, kinetic models for the electronic heating, a pellet surface ablation model, an equation of state that supports multiple ionization states, radiation, and a model for grad-B drift of the ablated material across the magnetic field. The Lagrangian particle algorithm is highly adaptive, capable of simulating a large number of fragments in 3D while eliminating numerical difficulties of dealing with the tokamak background plasma. The code has achieved good agreement with theory for spherically symmetric ablation flows. Axisymmetric simulations of neon and deuterium pellets in magnetic fields ranging from 1 to 6 Tesla have been compared with previous simulations using the FronTier code, and very good agreement has also been obtained. Furthermore, the main physics contribution of the paper is a detailed study of the influence of 3D effects, in particular grad-B drift, on pellet ablation rates and properties of ablation clouds. Smaller reductions of ablation rates in magnetic fields compared to axially symmetric simulations have been demonstrated because the ablated material is not confined to narrowing channels in the presence of grad-B drift. Contribution of various factors in the grad-B drift model has also been quantified.

97 MATHEMATICS AND COMPUTING↗