Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “etc)”

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 667 records · Page 37

Solidification and grain formation in alloys: a 2D application of the grand-potential-based phase-field approach

Solidification is a significant step in the forming of crystalline structures during various manufacturing and material processing techniques. Solidification characteristics and the microstructures formed during the process dictate the properties and performance of the materials. Hence, understanding how the process conditions relate to various microstructure formations is paramount. In this work, a grand-potential-based multi-phase, multi-component, multi-order-parameter phase-field model is used to demonstrate the solidification of alloys in 2D. This model has several key advantages over other multi-phase models such as it decouples the bulk energy from the interfacial energy, removes the constraints for the phase concentration variable, and prevents spurious third-phase formation at the two phase interfaces. Here, the model is implemented in a finite-element-based phase-field modeling code. The role of various modeling parameters in governing the solidification rate and the shape of the solidified structure is evaluated. It is demonstrated that the process conditions such as temperature gradient, thermal diffusion, cooling rate, etc, influence the solidification characteristics by altering the level of undercooling. Furthermore, the capability of the model to capture directional solidification and polycrystalline structure formation exhibiting various grain shapes is illustrated. In both these cases, the process conditions have been related to the growth rate and associated shape of the dendritic structure. As a result, this work serves as a stepping stone towards resolving the larger problem of understanding the process–structure–property–performance correlation in solidified materials.

36 MATERIALS SCIENCE↗

Electric-field-induced crossover of polarization reversal mechanisms in Al 1–x Sc x N ferroelectrics

Scandium-doped aluminum nitride, Al 1–x Sc x N, represents a new class of displacive ferroelectric materials with high polarization and sharp hysteresis along with high-temperature resilience, facile synthesizability and compatibility with standard CMOS fabrication techniques. The fundamental physics behind the transformation of unswitchable piezoelectric AlN into switchable Al–Sc–N ferroelectrics depends upon important atomic properties such as local structure, dopant distributions and the presence of competing mechanism of polarization switching in the presence of an applied electric-field that have not been understood. We computationally synthesize Al 1–x Sc x N to quantify the inhomogeneity of Sc distribution and phase segregation, and characterize its crystal and electronic structure as a function of Sc-doping. Nudged elastic band calculations of the potential energy surface and quantum molecular dynamics simulations of direct electric-field-driven ferroelectric switching reveal a crossover between two polarization reversal mechanisms—inhomogeneous nucleation-and-growth mechanism originating near Sc-rich regions in the limit of low applied fields and nucleation-limited-switching in the high-field regime. Finally, understanding polarization reversal pathways for these two mechanisms as well as the role of local Sc concentration on activation barriers provides design rules to identify other combinations of dopant elements, such as Zr, Mg etc. to synthesize superior AlN-based ferroelectric materials.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Kinetic simulations of collision-less plasmas in open magnetic geometries

Laboratory plasmas in open magnetic geometries can be found in many different applications such as (a) scrape-of-layer (SOL) and divertor regions in toroidal confinement fusion devices, (b) linear divertor simulators, (c) plasma-based thrusters and (d) magnetic mirrors etc. A common feature of these plasma systems is the need to resolve, in addition to velocity space, at least one physical dimension (e.g. along flux lines) to capture the relevant physics. In general, this requires a kinetic treatment. Fully kinetic particle-in-cell (PIC) simulations can be applied but at the expense of large computational effort. A common way to resolve this is to use a hybrid approach: kinetic ions and fluid electrons. In the present work, the development of a hybrid PIC computational tool suitable for open magnetic geometries is described which includes (a) the effect of non-uniform magnetic fields, (b) finite fully-absorbing boundaries for the particles and (c) volumetric particle sources. Analytical expressions for the momentum transport in the paraxial limit are presented with their underlying assumptions and are used to validate the results from the PIC simulations. A general method is described to construct discrete particle distribution functions in a state of mirror-equilibrium. This method is used to obtain the initial state for the PIC simulation. Collisionless simulations in a mirror geometry are performed. The results show that the effect of magnetic compression is correctly described and momentum is conserved. The self-consistent electric field is calculated and is shown to modify the ion velocity distribution function in a manner consistent with analytic theory. Based on this analysis, the ion distribution function is understood in terms of a loss-cone distribution and an isotropic Maxwell-Boltzmann distribution driven by a volumetric plasma source. Finally, the inclusion of a Monte Carlo based Fokker-Planck collision operator is discussed in the context of future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Developing solid-surface plasma facing components for pilot plants and reactors with replenishable wall claddings and continuous surface conditioning. Part B: required research in present tokamaks

Abstract The companion part A paper (Stangeby et al 2022) reports a number of independent estimates indicating that high-duty-cycle DT tokamaks starting with pilot plants will likely experience rates of net erosion and deposition of solid PFC, plasma facing component, material in the range of 10 3 to 10 4 kg yr −1 , regardless of the material used. The subsequent redeposition of such large quantities of material has the potential for major interference with tokamak operation. Similar levels and issues will be involved if ∼continuous low-Z powder dropping is used for surface conditioning of DT tokamaks, independent of the material used for the PFC armor. In Stangeby et al (2022) (part A) it is proposed that for high-duty-cycle DT tokamaks, non-metallic low-Z refractory materials such as ceramics (graphite, SiC, etc) used as in situ replenishable, relatively thin—of order mm—claddings on a substrate which is resistant to neutron damage could provide a potential solution for protecting the main walls, while reducing the risk of degrading the confined plasma. Assessment of whether such an approach is viable will require information, much of which is not available today. Section 6 of part A identifies a partial list of major physics questions that will need to be answered in order to make an informed assessment. This part B report describes R&D needed to be done in present tokamaks in order to answer many of these questions. Most of the required R&D is to establish better understanding of low-Z slag generation and to identify means to safely manage it. Powder droppers provide a unique opportunity to carry out controlled studies on the management of low-Z slag in current tokamaks, independent of whether their protection tiles use low-Z or high-Z material.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Developing solid-surface plasma facing components for pilot plants and reactors with replenishable wall claddings and continuous surface conditioning. Part A: concepts and questions

Abstract It is estimated that pilot plants and reactors may experience rates of net erosion and deposition of solid plasma facing component (PFC) material of 10 3 –10 5 kg yr −1 . Even if the net erosion (wear) problem can be solved, the redeposition of so much material has the potential for major interference with operation, including disruptions due to so-called ‘unidentified flying objects (UFOs)’ and unsafe dust levels. The potential implications appear to be no less serious than for plasma contact with the divertor target: a dust explosion or a major UFO-disruption could be as damaging for an actively-cooled deuterium-tritium (DT) tokamak as target failure. It will therefore be necessary to manage material deposits to prevent their fouling operation. This situation appears to require a fundamental paradigm shift with regard to meeting the challenge of taming the plasma–material interface: it appears that any acceptable solid PFC material will in effect be flow-through , like liquid–metal PFCs, although at far lower mass flow rates. Solid PFC material will have to be treated as a consumable , like brake pads in cars. ITER will use high-Z (tungsten) armor on the divertor targets and low-Z (beryllium) on the main walls. The ARIES-AT reactor design calls for a similar arrangement, but with SiC cladding on the main walls. Non-metallic low-Z refractory materials such as ceramics (graphite, SiC, etc) used as in situ replenishable, relatively thin—of order mm—claddings on a substrate which is resistant to neutron damage could provide a potential solution for the main walls, while reducing the risk of degrading the confined plasma. Separately, wall conditioning has proven essential for achieving high performance. For DT devices, however, standard methods appear to be unworkable, but recently powder droppers injecting low-Z material ∼continuously into discharges have been quite effective and may be usable in DT devices as well. The resulting massive generation of low-Z debris, however, has the same potential to seriously disrupt operation as noted above. Powder droppers provide a unique opportunity to carry out controlled studies on the management of low-Z slag in all current tokamaks, independent of whether their protection tiles use low-Z or high-Z material.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The new ITER baseline, research plan and open R&D issues

A new baseline (NB) has been proposed by the ITER Project to ensure a robust achievement of the Projects’ goals, in view of past challenges including delays incurred due to the Covid-19 pandemic, technical challenges in completing first-of-a-kind components and in nuclear licensing. The NB includes modifications to the configuration of the ITER device and its ancillaries (e.g. change from beryllium to tungsten as first wall material, modification of the heating and current drive mix, etc.) as well as additional testing of components (e.g. toroidal field coils) or phased installation (start with inertially cooled first wall before later installation of the final actively water-cooled components) to minimise operational risks. In the NB, the ITER research plan (IRP) will be divided into three main phases: (a) start of research operation, with 40 MW of ECH and 10 MW of ICH, which will focus on the demonstration of 15 MA operation in L-mode, commissioning of all required systems, including disruption mitigation, and the demonstration of H-mode plasma operation in deuterium; (b) DT-1, with 60–67 MW of ECH, 33 MW of neutral beam injection (NBI) and 10–20 MW of ICH, which will demonstrate robust operation in high confinement H-mode plasmas in DT up to Q ⩾ 10 and for burn durations of 300–500 s within an accumulated neutron fluence of ∼1% of the ITER machine’s lifetime total, and; (c) DT-2, with up to 67 MW of ECH, up to 49.5 MW of NBI and up to 20 MW of ICH, with the ITER tokamak and ancillaries in their final configuration to demonstrate routine operation in DT plasmas at high Q and the Q ⩾ 5 long-pulse and steady-state scenarios to the final neutron fluence and to perform R&D on nuclear fusion reactor issues. The logic, physics basis, modelling and experimental evaluations carried out to support the NB and the associated IRP are described. These include the impact of the tungsten wall on plasma scenarios and associated risk mitigation measures, as well as the optimisation of the tokamak components and ancillaries to minimise Project risks. Open R&D issues related to these evaluations and mitigation measures are also described together with experimental, modelling and validation activities required to address them.

ITER↗

The principle of minimal power

This article is devoted to the memory of Yuri P Raizer, who passed away in 2021. He left a noticeable trace in gas discharge physics. The principle of minimal power (the state that requires minimal power is most probable) is thoroughly used in his books. Although the fundamental laws of physics do not imply this ad hoc principle, a detailed analysis of underlying phenomena can often reveal why nature prefers this path. Raizer illustrated this principle for plasma stratification, formation of electrode spots, discharge constriction, the shape of an arc channel, etc. We argue that the nonlinearity of equations describing gas discharges can often justify the realization of a plasma state maintained at minimal electric power. This nonlinearity appears because small groups of energetic electrons often control the ionization processes. The number of these electrons depends strongly on the ratio of the electric field to gas density, E/N. Under certain conditions, the ionization rate can also depend nonlinearly on electron density due to stepwise ionization and Coulomb collisions. Here, we use the principle of minimal power to illustrate some of Raizer's contributions to gas discharge physics from a single point of view. We demonstrate that nonlinearity of ionization processes in gas discharges can substantiate this principle for plasma stratification. However, striations of s, p, and r types in neon could exist with minimal or no ionization enhancement. This reminds us of Raizer's warning that applying the minimal power principle could lead to erroneous predictions, and a proper theory is required in each case to justify its use. 'The phenomenon of striations satisfies the principle of minimal power' – Yuri Raizer

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electron cyclotron drift instability and anomalous transport: two-fluid moment theory and modeling

In the presence of a strong electric field perpendicular to the magnetic field, the electron cross-field (E × B) flow relative to the unmagnetized ions can cause the so-called electron cyclotron drift instability (ECDI) due to resonances of the ion acoustic mode and the electron cyclotron harmonics. This occurs in, for example, collisionless shock ramps in space, and in E × B discharge devices such as Hall thrusters. A prominent feature of ECDI is its capability to induce an electron flow parallel to the background >E field at a speed greatly exceeding predictions by classical collision theory. Such anomalous transport is important due to its role in particle thermalization at space shocks, and in causing plasma flows towards the walls of E × B devices, leading to unfavorable erosion and performance degradation, etc. The development of ECDI and anomalous transport is often considered requiring a fully kinetic treatment. In this work, however, we demonstrate that a reduced variant of this instability, and more importantly, the associated anomalous transport, can be treated self-consistently in a collisionless two-fluid framework without any adjustable collision parameter. By treating both electron and ion species on an equal footing, the free energy due to the inter-species velocity shear allows the growth of an anomalous electron flow parallel to the background E field. We will first present linear analyses of the instability in the two-fluid five- and ten-moment models, and compare them against the fully-kinetic theory. At low temperatures, the two-fluid models predict the fastest-growing mode in good agreement with the kinetic result. Also, by including more (> = 10) moments, secondary (and possibly higher) unstable branches can be recovered. The dependence of the instability on ion-to-electron mass ratio, plasma temperature, and background B field strength is also thoroughly explored. We then carry out direct numerical simulations of the cross-field setup using the five-moment model. The development of the instability, as well as the anomalous transport, is confirmed and in excellent agreement with theoretical predictions. The force balance properties are also studied using the five-moment simulation data. Here this work casts new insights into the nature of ECDI and the associated anomalous transport and demonstrates the potential of the two-fluid moment model in efficient modeling of E × B plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Decoupling plasma, catalyst, and gaseous mechanisms for non-oxidative methane conversion

Direct non-oxidative methane (CH 4 ) conversion to value-added hydrogen (H 2 ) and C 2 products remains hindered by fundamental catalytic scaling constraints and rapid surface deactivation at elevated temperatures. Plasma-enabled catalysis offers a promising route to overcome the thermodynamic and kinetic barriers of direct non-oxidative methane upgrading at mild conditions, yet control over C–C product selectivity and catalyst stability remains elusive. Here, we establish a unified mechanistic framework including Langmuir–Hinshelwood (L–H) and Langmuir–Rideal (L–R) mechanisms that disentangles the roles of plasma excitation (including vibrationally activated species and radicals), surface temperature (T sur ), and catalyst binding energy in steering CH 4 conversion to H 2 and C 2 hydrocarbons. Through a combination of density functional theory (DFT) informed microkinetic modeling, in situ and ex situ surface characterization, and product quantification under dielectric barrier discharge conditions, we show that vibrationally excited CH 4 lowers activation barriers selectively for dissociative chemisorption, enabling surface activation across a wide range of transition metal catalysts at low thermal energy input. We find that once CH 4 is dissociatively chemisorbed, the branching between C 2 H 2 , C 2 H 4 , and C 2 H 6 is governed by surface properties (carbon binding energy, T sur , etc), regardless of plasma excitation. The DFT informed microkinetic model decouples the effects of molecular activation from surface properties and indentifies operating windows that maximize target yields while suppressing carbon accumulation and subsequent catalytic inactivation. Experiments on polycrystalline Cu/Al 2 O 3 , Ni/Al 2 O 3 , and Pt/Al 2 O 3 validate these predictions, revealing catalyst-dependent branching toward ethane or ethylene and distinct deactivation profiles. We unify these trends into a generalized three-dimensional plasma-thermal-catalytic design space, from which reduced descriptors such as T vib /T sur in the limit of vibrationally excited L–H pathways emerge as predictive metrics. These results enable rational tuning of methane conversion pathways and unlock selective C 2 formation using earth-abundant metals under mild plasma conditions.

catalyst inactivation↗

Self-assembly of wood-based shape memory composites triggered by solar-thermal energy

Transporting and assembling large, complex structures poses significant challenges due to their size, geometry, and cost. Additionally, the installation sites are often inaccessible or hazardous for humans, necessitating self-assembling capabilities in these structures. To mitigate these challenges, we propose using 3D printing materials with shape memory effect (SME) for both transport and construction. This approach involves developing 3D modular components into flat sheets for easier transportation, and then self-assembling into 3D structures on-site using solar energy. To gain a deeper understanding of the factors influencing material memory performance, we have chosen a composite PLA/WF, which is polylactic acid (PLA) with 20 wt% wood flour (WF) for this purpose, leveraging its high tensile modulus at 0.966 GPa, low cost, and sustainability. Printed shapes with this material can maintain a recovery ratio over 90% after 3 cycles. While traditional composites fillers (e.g. glass or carbon fiber) are added to enhance mechanical and thermal properties, the addition of bio-based fillers like WF accomplish similar goals without compromising sustainability. We conducted multiple experiments to demonstrate how environmental conditions (i.e. temperature) maximize the material’s SME. Although still at an early stage, this study provides initial insights into bridging the gap between the small-scale nature of shape memory polymers (SMPs) and their potential for large-scale additive manufacturing, addressing a critical need for efficient and sustainable construction. In the long term, we hope our study contributes to the design vision of utilizing SMPs for transportation, assembly, and deployment of complex structures, providing a new pathway for sustainable construction and transportation of large-scale structures to hard-to-access locations such as disaster-affected areas and remote deserts, etc.

4D printing↗

Polarization angle requirements for CMB B-mode experiments. Application to the LiteBIRD satellite

A methodology to provide the polarization angle requirements for different sets of detectors, at a given frequency of a CMB polarization experiment, is presented. The uncertainties in the polarization angle of each detector set are related to a given bias on the tensor-to-scalar ratio r parameter. The approach is grounded in using a linear combination of the detector sets to obtain the CMB polarization signal. In addition, assuming that the uncertainties on the polarization angle are in the small angle limit (lower than a few degrees), it is possible to derive analytic expressions to establish the requirements. The methodology also accounts for possible correlations among detectors, that may originate from the optics, wafers, etc. The approach is applied to the LiteBIRD space mission. We show that, for the most restrictive case (i.e., full correlation of the polarization angle systematics among detector sets), the requirements on the polarization angle uncertainties are of around 1 arcmin at the most sensitive frequency bands (i.e., ≈ 150 GHz) and of few tens of arcmin at the lowest (i.e., ≈ 40 GHz) and highest (i.e., ≈ 400 GHz) observational bands. Conversely, for the least restrictive case (i.e., no correlation of the polarization angle systematics among detector sets), the requirements are ≈ 5 times less restrictive than for the previous scenario. At the global and the telescope levels, polarization angle knowledge of a few arcmins is sufficient for correlated global systematic errors and can be relaxed by a factor of two for fully uncorrelated errors in detector polarization angle. The reported uncertainty levels are needed in order to have the bias on r due to systematics below the limit established by the LiteBIRD collaboration.

79 ASTRONOMY AND ASTROPHYSICS↗

Did we hear the sound of the Universe boiling? Analysis using the full fluid velocity profiles and NANOGrav 15-year data

Abstract In this paper, we analyse sound waves arising from a cosmic phase transition where the full velocity profile is taken into account as an explanation for the gravitational wave spectrum observed by multiple pulsar timing array groups. Unlike the broken power law used in the literature, in this scenario the power law after the peak depends on the macroscopic properties of the phase transition, allowing for a better fit with pulsar timing array (PTA) data. We compare the best fit with that obtained using the usual broken power law and, unsurprisingly, find a better fit with the gravitational wave (GW) spectrum that utilizes the full velocity profile. Even more importantly, the thermal parameters that produce the best fit are quite different. We then discuss models that can produce the best-fit point and complementary probes using CMB experiments and searches for light particles in DUNE, IceCUBE-Gen2, neutrinoless doubleβ-decay, and forward physics facilities (FPF) at the LHC like FASERν, etc.

Astronomy & Astrophysics↗

Constraining explosive nucleosynthesis by indirect reaction methods at storage rings using unstable beams in batch mode

Nuclear reaction studies on unstable isotopes can strongly help in improving our understanding of nucleosynthesis in stars. Indirect approaches to determining astrophysical reaction rates are increasingly common-place and undergoing continuous refinement. Of particular interest is the use of such indirect techniques at storage rings, which, among other aspects, allow to recycle rare unstable beams. Here, we propose to investigate the reaction rates of astrophysical interest using indirect methods (surrogate, Trojan horse, etc.) in reverse kinematics at the IMP-CAS storage ring. Long lived radioactive ion beams, produced remotely, can be accelerated, and made to interact with light targets. The proposed reactions are 85 Kr(p, p’γ), 85 Kr(d, pγ), constraining the neutron flux in an s-process branching point, 79 Se(p, p’γ), 79 Se(d, pγ), constraining the temperature in s-process nucleosyntheses, and 59 Fe(d, pγ), constraining core collapse supernovae.

Angelis, G. de [National Inst. of Nuclear Physics ↗

LAMOST Time-Domain survey: first results of four K2 plates

From Oct. 2019 to Apr. 2020, LAMOST performed a time-domain (TD) spectroscopic survey of four K2 plates with both low- and medium-resolution observations. The low-resolution spectroscopic survey acquired 282 exposures ( ≈ 46.6 h) over 25 nights, yielding a total of about 767 000 spectra, and the medium-resolution survey took 177 exposures ( ≈ 49.1 h) over 27 nights, collecting about 478 000 spectra. More than 70%/50% of low-resolution/medium-resolution spectra have signal-to-noise ratio higher than 10. We determine stellar parameters (e.g., T eff , log g, [Fe/H]) and radial velocity (RV) with different methods, including LASP, DD-Payne and SLAM. In general, these parameter estimations from different methods show good agreement, and the stellar parameter values are consistent with those of APOGEE. We use the Gaia DR2 RV values to calculate a median RV zero point (RVZP) for each spectrograph exposure by exposure, and the RVZP-corrected RVs agree well with the APOGEE data. The stellar evolutionary and spectroscopic masses are estimated based on the stellar parameters, multi-band magnitudes, distances and extinction values. Finally, we construct a binary catalog including about 2700 candidates by analyzing their light curves, fitting the RV data, calculating the binarity parameters from medium-resolution spectra and cross-matching the spatially resolved binary catalog from Gaia EDR3. In conclusion, the LAMOST TD survey is expected to represent a breakthrough in various scientific topics, such as binary systems, stellar activity, stellar pulsation, etc.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantification of the effect of uncertainty on impurity migration in PISCES-A simulated with GITR

A Bayesian inference strategy has been used to estimate uncertain inputs to global impurity transport code (GITR) modeling predictions of tungsten erosion and migration in the linear plasma device, PISCES-A. This allows quantification of GITR output uncertainty based on the uncertainties in measured PISCES-A plasma electron density and temperature profiles (n e , T e ) used as inputs to GITR. The technique has been applied for comparison to dedicated experiments performed for high (4 × 10 22 m –2 s –1 ) and low (5 × 10 21 m –2 s –1 ) flux 250 eV He–plasma exposed tungsten (W) targets designed to assess the net and gross erosion of tungsten, and corresponding W impurity transport. The W target design and orientation, impurity collector, and diagnostics, have been designed to eliminate complexities associated with tokamak divertor plasma exposures (inclined target, mixed plasma species, re-erosion, etc) to benchmark results against the trace impurity transport model simulated by GITR. The simulated results of the erosion, migration, and re-deposition of W during the experiment from the GITR code coupled to materials response models are presented. Specifically, the modeled and experimental W I emission spectroscopy data for a 429.4 nm line and net erosion through the target and collector mass difference measurements are compared. Furthermore, the methodology provides predictions of observable quantities of interest with quantified uncertainty, allowing estimation of moments, together with the sensitivities to plasma temperature and density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of zonal flows on ion temperature gradient instability in the scrape-off layer of a field-reversed configuration

Gyrokinetic simulations of long wavelength ion temperature gradient (ITG) turbulence in the scrape-off layer (SOL) of a field-reversed configuration (FRC) find that zonal flows are nonlinearly generated and are the dominant mechanism for the nonlinear saturation of the ITG instability. After the ITG saturation, zonal flows remain undamped and gradually suppress the turbulent transport to a very low level. In the simulations with collisions, collisional damping gradually reduces zonal flow amplitude to a lower level, which allows finite ITG turbulence intensity and ion heat transport in the SOL. The steady state turbulence intensity and ion heat transport are found to be proportional to the collision frequency. Furthermore, this favorable scaling suggests that minimizing collisions (e.g. increasing temperature, reducing impurity content, etc) and preserving toroidal symmetry could improve plasma confinement in the FRC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A reduced model for the ITER divertor based on SOLPS solutions for ITER Q = 10 baseline conditions: A. identifying options for the control parameters*

Edge codes such as SOLPS coupled to neutral codes such as EIRENE have become so comprehensive and sophisticated that they now constitute, in effect, 'code-experiments' that, as for actual experiments, can benefit from interpretation using simple models and conceptual frameworks, i.e. reduced models. The first task is the identification of options for the reduced model control parameters that are best suited for control of the action of the divertor, i.e. for control of target power loading and sputter-erosion, primarily. A strong correlation between the electron temperature at the divertor target, T e,t , and the neutral deuterium D 2 density at the target, n D2,t , flux-tube resolved, has recently been reported for a number of code studies including SOLPS-4.3 modeling of a set of ~50 ITER baseline cases: Q DT = 10, q 95 = 3, P SOL = 100 MW, metallic walls, and Ne seeding. This part A of the present study reports new results for largely the same ITER cases, confirming the strong correlations reported earlier between local values of T e,t , and (i) n D2,t , and (ii) normalized volumetric losses of power and pressure in the divertor. Strong correlations have now also been found, and are reported here for the first time, between T e,t and all of the divertor target quantities of practical interest. A physical explanation for this surprising result has not been fully identified; nevertheless it has encouraging implications for reduced modeling of the ITER divertor. For such ITER conditions, (i) the global Ne injection rate, Inj Ne [Ne/s], and (ii) the electron temperature at the location on the target where the peak power deposition occurs, T e,t @q ⊥,pk (eV), are found to be promising reduced model control parameters. In the companion report, part B, a reduced model for the ITER divertor is developed and described in detail, based on reversed-direction 2 point modelling, Rev2PM. The input to the reduced model is a value of the variable pair (T et @q ⟂,pk ,Inj Ne ) and the output are values of the various target as well as divertor-entrance quantities of practical interest, e.g. q ⊥,pk , n e,Xpt (the electron density at the poloidal location of the X-point), etc. In part B the reduced model is quantitatively characterized using one half of the code cases; it is then used to successfully predict (replicate) the code values of e.g. n e,Xpt for the other half of the cases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A reduced model for the ITER divertor based on SOLPS solutions for ITER Q = 10 baseline conditions: B. A reduced model based on reversed-direction two point modeling

Edge codes such as SOLPS coupled to neutral codes such as EIRENE have become so comprehensive and sophisticated that they now constitute, in effect, 'code-experiments' that, as for actual experiments, can benefit from interpretation using simple models and conceptual frameworks, i.e. reduced models. The first task is the identification of options for the reduced model control parameters that are best suited for control of the action of the divertor, i.e. for control of target power loading and sputter–erosion, primarily. A strong correlation between the electron temperature at the divertor target, T e,t , and the neutral deuterium D 2 density at the target, n D2,t , flux-tube resolved, has recently been reported for a number of code studies including SOLPS-4.3 modeling of a set of ~50 ITER baseline cases: Q DT = 10, q 95 = 3, P SOL = 100 MW, metallic walls, and Ne seeding (Pitts et al 2019). Part A of the present study reports new results for largely the same ITER cases, confirming the strong correlation reported earlier between local values of T e,t , and (i) n D2,t , and (ii) normalized volumetric losses of power and pressure in the divertor. Strong correlations have now also been found, and are reported here for the first time, between T e,t and all of the divertor target quantities of practical interest. A physical explanation for this surprising result has not as yet been fully identified; nevertheless it has encouraging implications for reduced modeling of the ITER divertor. For such ITER conditions, (i) the global Ne injection rate, Inj Ne (Ne s –1 ), and (ii) the electron temperature at the location on the target where the peak power deposition occurs, T e,t @q ⊥,pk (eV), are found to be promising reduced model control parameters. In this part B, a reduced model for the ITER divertor is developed and described in detail, based on reversed-direction two point modeling, Rev2PM. The input to the reduced model is a value of the variable pair (T e,t @q ⊥,pk , Inj Ne )for a chosen case and the output are values of the various target as well as divertor-entrance quantities of practical interest, e.g. q ⊥,pk , the electron density at the X-point, n e,Xpt , etc. The reduced model was quantitatively characterized using one half of the code cases; it was then used to successfully predict (replicate) the code values of e.g. n e,Xpt for the other half.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗