Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Code of Accounts”

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

FSISPH: An SPH formulation for impacts between dissimilar materials

Here, we present an SPH formulation with several new features designed to better model the fully-compressible interaction of dissimilar materials. We developed the new method to simulate the atmospheric entry and break-up of small celestial bodies in planetary atmospheres. The formulation uses a unity-based, density-energy discretization of the hydrodynamic conservation laws with linear-corrected kernel gradients. To account for variations in compressibility, we use an HLLC approximate Riemann solver to adjust the velocity gradient at material interfaces. To handle large transverse velocity discontinuities, we introduce a simple slip interface model that limits the artificial viscosity at material interfaces. Diffusion is optionally applied through the velocity gradient and this allows the density and specific thermal energy to evolve in a manner more consistent with the first law of thermodynamics in comparison to other more direct diffusion schemes. We also introduce a material-local second-order artificial conduction scheme used to smooth the specific thermal energy field. Material damage fits neatly under this framework by treating the damage front as a material interface. The method has been implemented as a solver, FSISPH, within the code, Spheral++, and is publicly available on github. We test our new solver on a number of classic shock, mixing, and multi-material problem. The components we outline can significantly improve accuracy of SPH for problems with sharp contact discontinuities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Toroidal modeling of runaway electron loss due to 3D fields in ITER

Mitigation of runaway electrons (REs) by three-dimensional (3D) magnetic field perturbations is numerically investigated for the ITER 15 MA baseline D–T scenario, utilizing the MARS-F code (Liu et al Phys. Plasmas 7 3681) with a drift orbit test particle tracing module. Considered are two types of 3D fields: the n = 3 (n is the toroidal mode number) resonant magnetic perturbation (RMP) utilized for the purpose of controlling the edge localized modes in ITER, and perturbations generated by the n = 1 magneto-hydrodynamic (MHD) instabilities in a post-disruption plasma. The RMP field, applied to a pre-disruption plasma, is found to be moderately effective in mitigating the RE seeds in ITER when vacuum field model is assumed. Up to ~40% loss fraction is possible at 90 kA-turn coil current. The mitigation efficiency is however substantially reduced, down to less than 5%, when the plasma response is taken into account. This is due to strong screening of the resonant magnetic field components by the plasma response resulting in much less field line stochasticity. On the other hand, the MARS-F modeling, based on the DINA-simulated post-disruption equilibria, shows that the n = 1 resistive kink instabilities develop in these plasmas, as the edge safety factor qa evolves and drops below integer numbers. RE mitigation by these MHD instabilities is sensitive to the eigenmode structure. The best mitigation is achieved as q a drops below 3, when a global kink instability occurs that encompasses both internal and external components. This global instability is found to be capable of mitigating over 80% MeV-level passing RE orbits at a field perturbation |δB|/B 0 that is comparable to that observed in DIII-D experiments, and full mitigation if the perturbation amplitude is doubled. The 'wetted' area on the ITER limiting surface, due to MHD instability induced RE loss, generally increases with the perturbation amplitude (together with increasing loss fraction). At the highest perturbation level assumed in this study, the wetted area reaches ~60% of the total limiting surface area. Lastly, the lost RE orbits mainly strike the outer divertor region of the limiting surface, with some fraction also hitting a wide area along the inboard side of the surface.

3D fields↗
DOE OSTI2025

Slicing Solutions for Wire Arc Additive Manufacturing

Both commercial and research applications of wire arc additive manufacturing (WAAM) have seen considerable growth in the additive manufacturing of metallic components. However, there remains a clear lack of a unified paradigm for toolpath generation when slicing parts for WAAM deposition. Existing toolpath generation options typically lack the appropriate features to account for all complexities of the WAAM process. This manuscript explores the key slicing challenges specific to toolpaths for WAAM geometry and pairs each consideration with multiple solutions to mitigate most negative effects on completed components. These challenges must be addressed to minimize voids, prevent bead collapse, and ensure deposited components accurately approximate the desired geometry. Slicing considerations are grouped into four general categories: geometric, process, thermal, and productivity. Geometric considerations are addressed with overhang compensation, corner-sharpening, and toolpath-smoothing features. Process considerations are addressed with start point configuration and controls for the bead lengths and end points. Thermal and productivity considerations are addressed with island optimization, multi-material printing, and connected insets. Finally, tools for the post-processing of generated G-code are explored. Overall, these solutions represent a critical set of slicing features used to improve generated toolpaths and the quality of the components deposited with those toolpaths.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Porous Flow Modeling of Axial Gas Redistribution in Fragmented LWR Fuel Rods using MOOSE

Understanding how gas axially redistributes within fragmented fuel pellets is crucial for predicting the behavior of Light Water Reactor (LWR) fuel rods, particularly during transient and accidental scenarios. The time scale of this phenomenon plays a fundamental role in determining the progression and hazard of a Loss Of Coolant Accident (LOCA), especially when high burn-up fuel in a severe state of fragmentation is involved. Here, this study presents a Computational Fluid Dynamics (CFD) model developed within the Multiphysics Object-Oriented Simulation Environment (MOOSE) to predict the time-scale of plenum depressurization in Light-Water Reactor (LWR) fuel rods driven by axial gas transport through fragmented pellets. The model examines the effects of incorporating non-linearities in the friction term by comparing the results with experimental data. The experiment employed surrogate fuel rods containing pellets subjected to mechanical and/or thermal loadings to simulate various severity of cracking, and aimed at studying the influence of fuel conditions on axial gas redistribution. The results of this analysis indicate that under certain flow regime conditions - determined by the value of an equivalent Reynolds number - accounting for the non-linear friction term in Navier-Stokes equations guarantees better predictions for the time-scale of plenum depressurization. Also, the model enabled the simulation of the pressure decay by assigning distinct permeability values to each pellet instead of a single uniform value. Multiple simulations were run across all possible pellet position combinations, having each pellet assigned with values of permeability extracted from the experimental data. This allows to quantify the impact of the considering various non-uniform distributions of permeability on the dynamics of axial gas redistribution. The present work findings enhance the understanding of axial gas transport, and provide valuable insights for the integration of a model for predicting the axial gas redistribution during a LOCA scenario into the BISON fuel performance code.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Resiliency of Degraded Built Infrastructure

Infrastructure resiliency depends on the ability of infrastructure systems to withstand, adapt, and recover from chronic and extreme stresses. In this white paper, we address the resiliency of infrastructure assets and discuss improving infrastructure stability through development of our understanding of cement and concrete degradation. The resiliency of infrastructure during extreme events relies on the condition, adaptability, and recoverability of built infrastructure (roads, bridges, dams), which serves as the backbone of existing infrastructure systems. Much of the built infrastructure in the US has consistently been rated D+ by the American Society of Civil Engineers (ASCE). Aged infrastructure introduces risk to the system, since unreliable infrastructure increases the likelihood of failures under chronic and extreme stress and are particularly concerning when extreme events occur. To understand and account for this added risk from poor infrastructure quality, more research is needed on (i) how the changing environment alters the aging of new and existing built infrastructure and (ii) how degradation causes unique failure mechanisms. The aging of built infrastructure is based on degradation of the structural materials, such as concrete and steel supports, which causes failure. Current work in cement/concrete degradation is based on (i) the development of high strength and degradation resistance concrete mixtures, (ii) methods of assessing the age and reliability of existing structures, and (3) modeling of structural stability and the microstructural evolution of concrete/cement from degradation mechanisms (sulfide attack, carbonation, decalcification). Sandia National Laboratories (SNL) has made several investments in studying the durability and degradation of cement based materials, including using SNL-developed codes and methodologies (peridynamics, PFLOTRAN) to focus on chemo-mechanical fracture of cement for energy applications. Additionally, a recent collaboration with the University of Colorado Boulder has included fracture of concrete gravity dams, scaling the existing work to applications in full sized infrastructure problems. Ultimately, SNL has the experience in degradation of cementitious materials to extend the current research portfolio and answer concerns about the resilience of aging built infrastructure.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Analytic model of dislocation density evolution in fcc polycrystals accounting for dislocation generation, storage, and dynamic recovery mechanisms

Here, an analytic model of the evolution of dislocation density in fcc polycrystals is described. The evolution equations approximately account for most known dislocation storage, dynamic recovery, and dislocation generation mechanisms in fcc polycrystals. Specifically, the model incorporates network (forest) and grain boundary storage, mobile-network and mobile–mobile annihilation, screw–screw annihilation via athermal and thermal single cross-slip, generation by double cross-slip (Koehler mechanism, including dipole formation), Frank-Read sources, grain boundary nucleation, and mobile–immobile dislocation nucleation due to shock loading. Single cross-slip is assumed to proceed through the Friedel–Escaig (FE) mechanism; the corresponding activation energy is calculated using a modified FE model. The activation energy for double cross-slip is calculated for the first time by extending the FE model. The exact evolution equations are integro-differential equations, and as such are difficult to implement in a code; hence, the evolution equations are simplified by making several approximations. Preliminary results on copper are presented, including comparisons to experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Inertial Confinement Fusion Design Search Using Bayesian Optimization

Inertial confinement fusion (ICF) experiments rely on complex multi-physics simulation codes such as the Lawrence Livermore National Laboratory-developed HYDRA to guide design work. However, these simulations have several dozen tunable parameters and can be computationally expensive. This makes searching the parameter space challenging and time-consuming. Recently developed automated tools utilize Bayesian optimization to search these high-dimensional parameter spaces for optimal designs. The optimization tools run 2D integrated simulations in HYDRA to converge on a design that produces specified scalar or vector outputs. In this paper, we apply the Bayesian optimization tools to two common tuning scenarios. First, we tune simulation inputs to match measurements of a well-characterized experiment at the National Ignition Facility. This type of tuning is commonly performed to compensate for the use of simplified simulation settings (e.g. reduced resolution) or to account for missing physics in the simulations. Second, we search for an ICF simulation design that has a particular radiation drive profile. These optimizations replicate the kinds of tuning researchers routinely perform, but do so with significantly reduced manual effort. This approach demonstrates a powerful and efficient pathway toward autonomous, high-fidelity design optimization for future ICF experiments.

Bayesian optimization↗
DOE OSTI2025

A novel approach to increase accuracy in remotely sensed evapotranspiration through basin water balance and flux tower constraints

Remote sensing-derived evapotranspiration (RSET) products capture the spatiotemporal variations of evapotranspiration (ET) from field to basin scales with unprecedented details. However, their accuracy varies across RSET estimation methods and diverse hydroclimate regions. While ET modeling efforts to account for biophysical processes and controlling parameters have made good progress in recent years, a parallel approach of integrating in-situ ET with RSET could reduce biases in RSET products. Basin water balance ET (WBET) and flux tower ET are widely applied to evaluate RSET accuracy, yet such ET measurements are rarely used for RSET bias corrections, especially for large area applications. To address this issue, we propose a novel approach: the water balance equivalence (WABE) method, which generates spatially continuous WBET for correcting biases in RSET products. The WABE method computes synthetic WBET by integrating observed WBET and flux tower-derived FLUXCOM ET, which fills the spatial gaps of observed WBET and generates a spatially continuous WBET dataset. Synthetic WBET (2002–2015 annual average) of eight-digit hydrologic unit code (HUC8) basins across the conterminous United States (CONUS), constituting 44 % (887 out of 2035 basins) of CONUS basins, was determined within 2.0 % (RMSE = 12 %) of observed WBET at CONUS and between 1–12 % (RMSE = 3–33 %) across 18 regions in CONUS. With WABE-based bias corrections, the overall annual bias of RSET decreased from 10 % (RMSE = 34 %) to 6 % (RMSE = 26 %) across 37 flux tower sites. The WABE method offers a new approach for RSET accuracy improvement and shows great promise for large area implementations with a potential to yield substantial benefits for building accurate basin water budgets and water management decisions.

Khand, Kul↗
DOE OSTI2024

Simulations of nozzle gas flow and gas-puff Z-pinch implosions on the Weizmann Z-pinch

We present simulations of an oxygen gas puff Z-pinch on a University scale generator at the Weizmann Institute of Science. The work accounts for the detailed geometry of the nozzle, the initial neutral gas density distribution, and the subsequent implosion. The modeling results show significant improvement with data for the current at the time of stagnation in comparison with a previous effort [Rosenzweig et al., Phys. Plasmas 27, 022705 (2020)]. As a first step, we performed simulations of the flow of neutral diatomic oxygen from a plenum through a nozzle within a recessed cathode, across a gap, and into the anode with a recessed grounded honeycomb. These simulations show an agreement with the measured initial gas density profiles within the region not blocked by the recesses and accessible to visible measurements. The computed neutral gas flow profile serves as the initial condition for a radiation magnetohydrodynamic simulation of the implosion using the MACH2-TCRE code. By considering the specific details of the nozzle and chamber geometry, we find agreement with the measured current profile, including the inductive notch. The simulations predict that the plasma undergoes a strong pinch within the hidden anode recess. The simulations also predict the strongest radiation pulse occurs within the anode recess and at the time of the observed inductive notch.

Physics↗
DOE OSTI2025

Identification of functional non-coding variants associated with orofacial cleft

Oral facial cleft (OFC) comprises cleft lip with or without cleft palate (CL/P) or cleft palate only. Genome wide association studies (GWAS) of isolated OFC have identified common single nucleotide polymorphisms (SNPs) in many genomic loci where the presumed effector gene (for example, IRF6 in the 1q32 locus) is expressed in embryonic oral epithelium. To identify candidates for functional SNPs at eight such loci we conduct a massively parallel reporter assay in a fetal oral epithelial cell line, revealing SNPs with allele-specific effects on enhancer activity. We filter these SNPs against chromatin-mark evidence of enhancers and test a subset in traditional reporter assays, which support the candidacy of SNPs at loci containing FOXE1, IRF6, MAFB, TFAP2A, and TP63. For two SNPs near IRF6 and one near FOXE1, we engineer the genome of induced pluripotent stem cells, differentiate the cells into embryonic oral epithelium, and discover allele-specific effects on the levels of effector gene expression, and, in two cases, the binding affinity of transcription factors FOXE1 or ETS2. Conditional analyses of GWAS data suggest the two functional SNPs near IRF6 account for the majority of risk for CL/P at this locus. This study connects genetic variation associated with OFC to mechanisms of pathogenesis.

Kumari, Priyanka↗
DOE OSTI2025

2D end-to-end modelling of kilonovae from binary neutron star merger remnants

We investigate the kilonova emission resulting from outflows produced in a 3D general-relativistic magnetohydrodynamic (GRMHD) simulation of a hypermassive neutron star (HMNS) remnant. We map the outflows into the flash hydrodynamics code to model their expansion in axisymmetry, and study the effects of employing different r-process heating rates. Except for the highest heating rate prescription, we find no significant differences with respect to overall ejecta dynamics and morphology compared to the simulation without heating. Once homologous expansion is attained, typically after ~2s for these ejecta, we map the outflows to the sedona radiative transfer code and compute the spectral evolution of the kilonova and broad-band light curves in various Legacy Survey of Space and Time (LSST) bands. The kilonova properties depend on the remnant lifetime, with peak luminosities and peak time-scales increasing for longer lived remnants that produce more massive ejecta. For all models, there is a strong dependence of both the bolometric and broad-band light curves on the viewing angle. While the short-lived (12 ms) remnant produces higher luminosities when viewed from angles closer to the pole, longer lived remnants (240 ms and 2.5 s) are more luminous when viewed from angles closer to the equator. Our results highlight the importance of self-consistent, long-term modelling of merger ejecta, and taking viewing-angle dependence into account when interpreting observed kilonova light curves. We find that magnetized outflows from an HMNS – if it survives long enough – could explain blue kilonovae, such as the blue emission seen in AT2017gfo.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

MARS-F/K modeling of plasma response and fast ion losses due to RMP in KSTAR

Abstract The toroidal single-fluid magnetohydrodynamic (MHD) code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) and the MHD-kinetic hybrid code MARS-K (Liu et al 2008 Phys. Plasmas 15 112503) are utilized to study the plasma response to the n = 1 ( n is the toroidal mode number) resonant magnetic perturbation (RMP), applied to suppress the type-I edge localized mode (ELM) in a KSTAR discharge. Both the resistive-rotating and ideal-static plasma models identify strong screening of the resonant radial field harmonics of the applied RMP due to the plasma response, and predict a strong edge-peeling response of the plasma which is consistent with the optimal ELM control coil current configuration adopted in experiment. The RMP-induced radial displacement of the plasma, computed by the resistive-rotating plasma model, agrees reasonably well with that reconstructed from the measured data in the plasma core. Taking into account the drift kinetic response of fast ions, MARS-K hybrid modeling also finds quantitative agreement of the plasma core fluid pressure perturbation with experiment. Based on the MARS-F computed plasma response, a guiding-center orbit-tracing simulation finds about 0.3% of fast ion losses due the n = 1 RMP in the KSTAR ELM control experiment considered. Most losses are associated with counter-current fast ions located near the plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Modeling Radionuclide Inventories in MSR Off-Gas Systems with Radiochemical Transport Analysis

Molten salt reactors (MSR) contain unique characteristics that may require enhancements to modeling tools to accurately predict phenomena. One characteristic that may be advantageous to leverage during normal operation is on-line processing of the circulating fuel salt, such as an off-gas system (OGS) to remove volatile fission products. Therefore, new modeling tools must be developed to integrate spatial resolution and chemistry effects into fuel depletion tools to be able to account for these non-core sources of radioactivity. Such types of radiochemical transport analysis tools were used to estimate the removal rates for 12 elements within a flow model of the Molten Salt Reactor Experiment (MSRE) by optimizing against legacy experimental data of the gas-borne (GB) percentages of 12 nuclides. The removal rates were used in a depletion model to calculate the FP inventory that enters the MSRE OGS. Calculations are in good agreement with the empirical GB percentages reported for the 12 nuclides, which validates the approach and verifies each tool’s treatment of the radiochemical flow effects. The OGS inventory is discussed in terms of the largest nuclide contributors to activity, dose consequence, decay heat, and elemental composition. Finally, insights from the study allow recommendations to be made for future code development activities.

Shahbazi, Shayan↗
DOE OSTI2019

Solar soft cost ontology: a review of solar soft costs

The continued deployment of renewable energy is critical to the energy transition necessary for climate change mitigation. Since cost is one of the primary drivers of solar adoption, we must continue to reduce the cost of solar to hasten the transition to clean energy. Currently, solar soft costs account for 52%–70% of the cost of an installed residential solar photovoltaic (PV) system in the U.S. These costs are persistent and an increasing share of the costs, since hardware costs continue to decline. The consensus on a definition of soft costs is 'non-hardware costs', but defining something by what it is not does not tell us what it is, which makes it difficult to study and understand. To this end, we developed the solar soft cost ontology (SSCO) to provide a comprehensive view of the complex and variable knowledge domain of solar soft costs. In this work, we identify the main categories of solar soft costs, how they relate to each other, how that knowledge is created, shared, and acquired, and the relevant actors in the solar soft cost knowledge ecosystem. Developed by coding a corpus of nearly 130 academic articles, the resulting ontology in web ontology language (OWL) includes 136 classes and 87 properties. The SSCO can facilitate and enhance knowledge transfer of research findings and best practices within and between stakeholder communities, including researchers, practitioners, and policy makers, that seek to address solar soft costs. We conclude with a discussion of potential applications and use cases for the SSCO among this diverse stakeholder community.

14 SOLAR ENERGY↗
DOE OSTI2020

Wave optics simulations for new soft x-ray beamlines at NSLS-II

In this work, we performed fully- and partially-coherent synchrotron emission and propagation simulations with the "Synchrotron Radiation Workshop" computer code to analyze the performance of two soft X-ray beamlines under development at the National Synchrotron Light Source II: Soft X-ray Nanoprobe (SXN), and Angle-Resolved Photoemission Spectroscopy (ARPES) and Resonant Inelastic X-ray Scattering (RIXS) Imaging (ARI). The SXN beamline intends to provide high flux and high spatial resolution coherent soft X-ray imaging capabilities using both zone plate and lensless coherent imaging techniques. The ARI beamline aims to perform high flux ARPES and RIXS experiments with a focal spot size at the sample approaching 100 nm using highly-demagnifying mirrors in Kirkpatrick-Baez geometry. To accurately calculate the resolution and the degree of X-ray coherence provided by the two state-of-the-art beamlines, partial coherence effects are required to be taken into account in wave optics simulations for these two beamlines. In this talk, beamline performance parameters such as spot size, degree of coherence, flux, and energy resolution at the sample are presented. The effects of mirror surface slope errors on beamline performance were studied and some suggestions for further optimization are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

2-D soil zymography: Accounting for the spatial variation of pH

Soil zymography is commonly used to quantify spatial distribution of hydrolytic enzyme activities on soil and plant root surfaces. It is recommended to adjust pH in zymography substrates and calibration solutions with respect to soil/root pH. However, pH values may vary greatly within a few mm of plant rhizosphere, potentially altering the distribution of pH in zymography membranes. Despite the fact that the effect of pH on the calibration of zymography membranes is generally known, its potential impact on zymography results is unaccounted for in processing zymography images and calculations of enzyme activity. In this study we assessed the effect of pH variations on the persistency of the methylumbelliferone (MUF) calibration. The studied pH values ranged from 4.5 to 7.5. The MUF calibration curves greatly deviated from that at a reference pH of 6.5, with a marked nonlinear increase of deviation with greater membrane brightness. We suggest that the problem can be partially alleviated by reducing the membrane incubation time. However, such deviations suggest the need for a more comprehensive resolution via mapping pH and using pH-specific calibrations to process zymography images. Here, we developed a MATLAB code to implement a pixel-based correction of enzyme activity for pH in processing time-lapse zymography images.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Ab Initio Direct Dynamics

The reactivity and dynamics of molecular systems can be explored computationally by classical trajectory calculations. The traditional approach involves fitting a functional form of a potential energy surface (PES) to the energies from a large number of electronic structure calculations and then integrating numerous trajectories on this fitted PES to model the molecular dynamics. The ever-decreasing cost of computing and continuing advances in computational chemistry software have made it possible to use electronic structure calculations directly in molecular dynamics simulations without first having to construct a fitted PES. In this “on-the-fly” approach, every time the energy and its derivatives are needed for the integration of the equations of motion, they are obtained directly from quantum chemical calculations. This approach started to become practical in the mid-1990s as a result of increased availability of inexpensive computer resources and improved computational chemistry software. The application of direct dynamics calculations has grown rapidly over the last 25 years and would require a lengthy review article. The present Account is limited to some of our contributions to methods development and various applications. To improve the efficiency of direct dynamics calculations, we developed a Hessian-based predictor-corrector algorithm for integrating classical trajectories. Hessian updating made this even more efficient. Furthermore, this approach was also used to improve algorithms for following the steepest descent reaction paths. For larger molecular systems, we developed an extended Lagrangian approach in which the electronic structure is propagated along with the molecular structure. Strong field chemistry is a rapidly growing area, and to improve the accuracy of molecular dynamics in intense laser fields, we included the time-varying electric field in a novel predictor-corrector trajectory integration algorithm. Since intense laser fields can excite and ionize molecules, we extended our studies to include electron dynamics. Specifically, we developed code for time-dependent configuration interaction electron dynamics to simulate strong field ionization by intense laser pulses. Our initial application of ab initio direct dynamics in 1994 was to CH 2 O → H 2 + CO; the calculated vibrational distributions in the products were in very good agreement with experiment. In the intervening years, we have used direct dynamics to explore energy partitioning in various dissociation reactions, unimolecular dissociations yielding three fragments, reactions with branching after the transition state, nonstatistical dynamics of chemically activated molecules, dynamics of molecular fragmentation by intense infrared laser pulses, selective activation of specific dissociation channels by aligned intense infrared laser fields, angular dependence of strong field ionization, and simulation of sequential double ionization.

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

Enhancement of ECCD by the current condensation effect for stabilizing large magnetic islands caused by neoclassical tearing modes in tokamak plasmas

The radio frequency current condensation effect reported in Reiman and Fisch (2018 Phys. Rev. Lett. 121 225001) is modeled in the nonlinear resistive magnetohydrodynamic code. A series of numerical investigations have been performed to investigate the enhancement of electron cyclotron current drive (ECCD) by the current condensation effect during the control of neoclassical tearing mode (NTM) in tokamak plasmas. In the numerical model, both the parallel transport and the perpendicular transport of electron temperature are considered. The EC driven current and driven perturbed electron temperature can nonlinearly evolve within the given magnetic configuration and eventually reach saturation states. The input power threshold of ECCD and the fold bifurcation phenomenon are numerically verified via nonlinear simulations. The numerical results show good agreements with the analytical results. Moreover, spatial distributions of EC current for the two solutions at different condensed level are displayed. The control effectiveness of ECCD for large NTM islands has been evaluated while considering the current condensation effect. While taking into account current condensation effect, for a sufficiently large input power, a larger island can be more effectively stabilized than a smaller one, which suggests a reassessment of the previous idea that the ECCD should always be turned on as early as possible. The potential physics mechanism behind the ECCD control have all been discussed in detail. Furthermore, the condensation effect is found to have favorable effects on the radial misalignment of ECCD. In the consideration of the situation for extremely localized control needs, a highly peaked heating profile is adopted to verify that the fold bifurcation phenomenon still exists and the current condensation effect can still take effect in this extreme condition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗