Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Power Iteration”

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 217 records · Page 12

Thermal feedback coupling in a transient Monte Carlo high-order/low-order scheme

Transient simulations of nuclear systems face the computational challenge of resolving both space and time during reactivity changes. A common strategy for tackling this issue is to split the neutron flux into shape and amplitude functions. This split can be solved with high- order/low-order methods. While this multi-fidelity approach has traditionally been reserved for deterministic methods, it is also possible to implement in Monte Carlo as an efficient alternative to Dynamic Monte Carlo. This work implements the frequency transform method with thermal feedback in high-order/low-order Monte Carlo by blending static coupling methods such as single-batch Monte Carlo, with a simple thermal-fluids calculation. While previous work focused solely on prescribed transients, the addition of time-dependent thermal-fluids allows transients to be self-propagating. Tests were run in a fluids-initiated transient to showcase the basic functionalities of this methodology. Preliminary results behave as expected, paving the way for studying more sophisticated thermal-fluids coupling methods. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Numerical assessment of the new V-shape small-angle slot divertor on DIII-D

The small-angle slot (SAS) divertor of the DIII-D tokamak, and its upcoming upgrade, the V-shape small-angle slot (SAS-V) divertor, are numerically investigated using the SOLPS-ITER code package, including the effect of particle drifts, for a range of plasma density, heating power, strike point position in the slot, and for both magnetic field directions. The simulations show that the electron temperature near the strike point is reduced in SAS-V compared to SAS, for both magnetic field directions, such that SAS-V achieves divertor detachment at a lower value of the outboard mid-plane separatrix electron density. The detachment threshold is lower because the V-shape focuses recycling neutrals on the V-end, densifying and cooling the plasma in the slot. At sufficiently high density, the V-shape also reduces the radial gradient of the temperature profile at the target, which in turns reduces the radial electric field and the E x B drift velocities, further densifying and cooling the plasma in the slot and leading to detachment. The V-shape effect, however, is reduced for higher heating power. With more heating power, the detachment density increases, reducing the ionization mean free path of recycled neutrals, which therefore become less sensitive to target shape changes. Furthermore, this suggests that in a fusion reactor, where the heating power is high, optimization of the divertor target shape needs to be combined with other strategies to lower the detachment density, such as in-slot injection of low-Z impurities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DART-PFLOTRAN: An ensemble-based data assimilation system for estimating subsurface flow and transport model parameters

Ensemble-based Data Assimilation (EDA), based on the Monte Carlo approach, has been effectively applied to estimate model parameters through inverse modeling in subsurface flow and transport problems. However, implementation of EDA approach involves a complicated workflow that include setting up and executing ensemble forward model simulations, processing observations and model simulation results for parameter updates, and repeat for sequential or iterative EDA. To facilitate the management of such workflow and lower the barriers for adopting EDA-based parameter estimation in subsurface science, we develop a generic software frame-work linking the Data Assimilation Research Testbed (DART) with a massively parallel subsurface FLOw and TRANsport code PFLOTRAN. The new DART-PFLOTRAN leverages both the core data assimilation engines in DART and the computational power afforded by PFLOTRAN. In addition to the standard smoother and filtering options, DART-PFLOTRAN enables an iterative EDA workflow based on the Ensemble Smoother for Multiple Data Assimilation method (ES-MDA) to improve estimation accuracy for nonlinear forward problems. Here, we verify the implementation of ES-MDA in DART-PFLOTRAN using two synthetic cases designed to estimate static permeability and dynamic exchange fluxes across the riverbed, respectively, from continuous temperature measurements made across a depth profile. One-dimensional hydro-thermal simulations are performed in both cases to relate temperature responses with the parameters of interest. In the case of estimating dynamic parameters, we demonstrate the flexibility of DART-PFLOTRAN in automating sequential ES-MDA workflow, which will significantly reduce the time researchers spend on managing complex workflows in similar applications. Both studies yield accurate estimations of the parameters compared to their synthetic truth, while ES-MDA leads to more accurate estimation when a high level of nonlinearity exist between observed responses and unknown parameters. With a code base in Python and Fortran, DART-PFLOTRAN paves the way for applications in large-scale subsurface inverse modeling by automating the complex workflow of sequential ES-MDA that can be executed on various computing platforms.

97 MATHEMATICS AND COMPUTING↗

Numerical Modeling and Optimization of the iProTech Pitching Inertial Pump (PIP) Wave Energy Converter (WEC) (Cooperative Research and Development Final Report, CRADA Number: CRD-22-22968)

This work generated a first-of-its-kind automated workflow to couple time-domain simulations of wave energy converters written in one software language with a set of design generation and evaluation scripts written in another software language. This automated workflow used an existing optimization package to analyze the sensitivity of different design parameters on the power output of a specific WEC, iProTech’s Pitching Inertial Pump (PIP). Geometric, inertial, and power take-off variables were all varied and optimized to find values that produced the highest amount of power generated over varying wave conditions. The findings on these parameter sensitivity studies are used to inform future design iterations of the PIP WEC. Including more design variables in the optimizations will only increase computational run time and further software development is needed to analyze a larger optimization.

16 TIDAL AND WAVE POWER↗

A Unified Design Theory for Multi-Port Polyphase Transformers Enabling Scalable Power-Multiplexed EV Fleet Charging Systems

This paper presents a unified analytical design theory for multi-port polyphase transformers, targeting scalable and isolated high-power Electric Vehicle (EV) fleet charging systems with power multiplexing capability. As fleet electrification accelerates, conventional one-to-one charger architectures face significant challenges in infrastructure cost, peak power demand, and low utilization of installed power electronics. Power-multiplexed charging architectures, which dynamically distribute power from a shared pool of converter modules across multiple vehicles, have emerged as a promising solution. However, such architectures require scalable, isolated multi-port power interfaces capable of routing energy among multiple inputs and outputs, whose design remains complex and dependent on iterative modeling. To address this gap, the proposed theory provides closed-form expressions for self-inductance, leakage inductance, and mutual coupling terms for arbitrary multi-phase, multi-port transformer structures. The formulation enables direct synthesis of isolated multi-input and multi-output resonant converter systems without reliance on geometry-specific finite-element analysis or extensive parameter extraction. This capability is particularly critical for power-multiplexed systems, where modular converter structures must interface with multiple vehicles while maintaining galvanic isolation and flexible power allocation. The effectiveness of the proposed framework is demonstrated through the design of a 360 kW multi-phase system operating over a 700–900 VDC input and 400–1250 VDC output range. PLECS simulation results confirm accurate prediction of system behavior and validate the applicability of the approach to multi-port, power-multiplexed charging scenarios. The proposed method significantly reduces design complexity while enabling scalable, cost-effective, and fully utilized EV fleet charging infrastructure.

Asa, Erdem [ORNL] (ORCID:0000000190884812)↗

Recent progress in the design of the K-DEMO divertor

The preliminary conceptual design of the Korean fusion demonstration reactor (K-DEMO) with a major radius of 6.8 m and the fusion power of 2200 MW has been studied since 2012. The overall configuration of the K-DEMO divetor system based on the ITER-like water-cooled tungsten technology is a double-null type symmetric divertor subdivided into 32 toroidal modules for the vertical maintenance. A detached divertor scenario with impurity seeding was considered as the primary approach for the power exhaust to reduce the peak heat flux lower than the engineering limit of 10 MW/m 2 . The power exhaust performance at the scrap off layer was estimated by using UEDGE-2D code, a two-dimensional fluid transport code for collisional edge plasma and neutral species like N, Ne, and Ar. Particle and heat flux on inboard and outboard divertor targets were calculated for the detached cased depending on parameters such as the impurity seeding rate, pumping rate, and the pedestal density. On the other hand, a magnetic solution like X-divertor, snowflake divertor, and super X-divertor to expand the plasma wet area was considered for K-DEMO since the detached divertor increasing a radiation fraction by impurity seeding might be able to be unstable. However, the extremely high current of poloidal coils was required more than the engineering limit, 20 MA, to form magnetic field lines for the alternative divertors. Based on the physical calculation of the edge plasma, engineering analyses were carried out to find out the thermal and structural reliability. The thermo-hydraulic analysis confirmed thermal stability, whether all comprising materials are operating within their allowable temperature windows when the case of the peak heat flux is set to 10 MW/m 2 on the outboard divertor target. The response surface optimization method derived two optimal design candidates employing two kinds of heat sink materials, respectively: the reduced activation ferritic martensitic (RAFM) steel and CuCrZr alloy. Additionally, the drawbacks and merits of the two materials were definite. The optimal design with applying RAFM steel was vulnerable to withstand thermal and mechanical loads since low thermal conductivity caused too thin thickness of the heat sink. On the other hand, the CuCrZr alloy has critical drawbacks in terms of activation and radioactive waste despite its high thermal conductivity. Meanwhile, preliminary electromagnetic (EM) analysis was carried out to estimate the EM loads caused by the abnormal behaviors of plasma since EM loads are one of the most critical external loads for designing a DEMO divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ion-Pair Proton Exchange Membrane Fuel Cells for Heavy-duty Transportation

The purpose of this project is to build a fuel cell prototype that demonstrates high temperature proton exchange membrane fuel cell technology (HT PEMFC) commercial readiness. It advances from a technology baseline established in a prior CRADA in which these same Parties collaborated during the years 2020-2024. These efforts are aligned with the intentions of the DOE Hydrogen Fuel Cell Technology Office (HFTO) L’Innovator Pilot Program. Participant observes that its US-based competitive advantage may be improved by further optimizing Laboratory inventions. This CRADA is intended to further improve National Laboratory compositions and to develop commercial HT PEMFC membrane electrode assemblies (MEAs) that perform over a range of operating temperatures and relative humidity, that generate power from both low grade and pure hydrogen, and that are resistant to air impurities. Participant’s goal is to advance the commercial competitiveness of National Laboratory technologies in fuel cell products for heavy-duty truck integrators, for stationary power generation, marine applications, and for emerging interest from commercial aviation. The Parties will collaborate on producing a HT PEMFC MEA that integrates LANL ion pair technology and Brookhaven catalyst technology. The Laboratories plan to research and develop compositions for testing and evaluation. The participant will manufacture multiple MEAs based on the Laboratory compositions and will enable selected customers to validate performance characteristics. The Laboratories will convey to Participant the methods of making the National Laboratory compositions such that Participant may reproduce and then scale Laboratory methods into commercial production. Throughout the project, the Laboratories and Participant will iteratively test and optimize the compositions of materials of the MEAs, intending to optimize characteristics such as power density and lifecycle durability. Participant’s intended customers include fuel cells for large trucks, aviation, power back up, marine, and portable power.

99 GENERAL AND MISCELLANEOUS↗

Optimal Frequency-domain Analysis for Spacecraft Time Series: Introducing the Missing-data Multitaper Power Spectrum Estimator

While the Lomb–Scargle periodogram is foundational to astronomy, it has a significant shortcoming: the variance in the estimated power spectrum does not decrease as more data are acquired. Statisticians have a 60 yr history of developing variance-suppressing power spectrum estimators, but most are not used in astronomy because they are formulated for time series with uniform observing cadence and without seasonal or daily gaps. Here we demonstrate how to apply the missing-data multitaper power spectrum estimator to spacecraft data with uniform time intervals between observations but missing data during thruster fires or momentum dumps. The F-test for harmonic components may be applied to multitaper power spectrum estimates to identify statistically significant oscillations that would not rise above a white noise–based false alarm probability. Multitapering improves the dynamic range of the power spectrum estimate and suppresses spectral window artifacts. We show that the multitaper–F-test combination applied to Kepler observations of KIC 6102338 detects differential rotation without requiring iterative sinusoid fitting and subtraction. Significant signals reside at harmonics of both fundamental rotation frequencies and suggest an antisolar rotation profile. Next we use the missing-data multitaper power spectrum estimator to identify the oscillation modes responsible for the complex "scallop-shell" shape of the K2 light curve of EPIC 203354381. We argue that multitaper power spectrum estimators should be used for all time series with regular observing cadence.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimizing Power Generation of a Bottomraised Oscillating Surge Wave Energy Converter Using a Theoretical Model

Preliminary sizing of an oscillating surge wave energy converter (OSWEC) is an iterative process that relies on knowledge of the relevant hydrodynamic coefficients for a given geometry. Often, the complex definition of the device geometry requires coefficients to be obtained through experiments or numerical boundary element solvers such as WAMIT. These techniques demand significant user and computational effort, therefore inhibiting the fine scale parametrization of object dimensions. In this study, a theoretical model, originally presented in Michele et al. (2016), is developed and demonstrated to efficiently optimize the power production for an OSWEC device (subjected to certain environmental conditions) with variations in device widths, heights, and distances from the seabed. Assuming negligible device thickness, the OSWEC motions are modelled as a bottom-raised two-dimensional flap in regular waves using potential flow theory formulated in elliptical coordinates. The solutions to this diffraction-radiation problem are obtained using Mathieu functions with appropriate boundary conditions. The resulting potentials are then used to derive frequency-dependent expressions for the added mass and radiation damping coefficients, along with wave excitation magnitude in surge, pitch, and coupled surge-pitch motions. Good agreement in hydrodynamic coefficient curves is shown between the theoretical model and the numerical results obtained from the boundary element-based program WAMIT. The theoretical model is then employed to maximize the time-averaged output power while maintaining or reducing the hinge reaction force, with variations in device dimensions, wave frequency and amplitude.

ACE↗

Spectrally accelerated edge and scrape-off layer gyrokinetic turbulence simulations

This paper presents the first gyrokinetic (GK) simulations of edge and scrape-off layer (SOL) turbulence accelerated by a velocity-space spectral approach in the full-f GK code GENE-X. Building upon the original grid velocity-space discretization, we derive and implement a new spectral formulation and verify the numerical implementation using the method of manufactured solution. We conduct a series of spectral turbulence simulations focusing on the TCV-X21 reference case (Oliveira et al., 2022 [26]) and compare these results with previously validated grid simulations (Ulbl et al., 2023 [25]). The spectral approach reproduces the outboard midplane (OMP) profiles (density, temperature, and radial electric field), dominated by trapped electron mode (TEM) turbulence, with excellent agreement and significantly lower velocity-space resolution. As a consequence, the spectral approach reduces the computational cost (CPUh) by at least an order of magnitude, of approximately 50 for the TCV-X21 case. This enables high-fidelity GK simulations to be performed within a few days on modern CPU-based supercomputers for medium-sized devices and establishes GENE-X as a powerful tool for studying edge and SOL turbulence, moving towards reactor-relevant devices like ITER.

Gyrokinetic↗

Effect of magnetic perturbations on turbulence-flow dynamics at the L-H transition on DIII-D

Detailed 2D turbulence measurements from the DIII-D tokamak provide an explanation for how resonant magnetic perturbations (RMPs) raise the L-H power threshold P LH [P. Gohil et al. , Nucl. Fusion 51 , 103020 (2011)] in ITER-relevant, low rotation, ITER-similar-shape plasmas with favorable ion ∇ B direction. RMPs simultaneously raise the turbulence decorrelation rate Δ ω D and reduce the flow shear rate ω shear in the stationary L-mode state preceding the L-H transition, thereby disrupting the turbulence shear suppression mechanism. RMPs also reduce the Reynolds stress drive for poloidal flow, contributing to the reduction of ω shear . On the ~100 μ s timescale of the L-H transition, RMPs reduce Reynolds-stress-driven energy transfer from turbulence to flows by an order of magnitude, challenging the energy depletion theory for the L-H trigger mechanism. In contrast, non-resonant magnetic perturbations, which do not significantly affect P LH , do not affect Δ ω D and only slightly reduce ω shear and Reynolds-stress-driven energy transfer.<

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Novel SOLPS-ITER simulations of X-point target and snowflake divertors

Abstract The design and understanding of alternative divertor configurations may be crucial for achieving acceptable steady-state heat and particle material loads for magnetic confinement fusion reactors. Multiple X-point alternative divertor geometries such as snowflakes and X-point targets have great potential in reducing power loads, but have not yet been simulated widely in codes with kinetic neutrals. This paper discusses recent changes made to the SOLPS-ITER code to allow for the simulation of X-point target and low-field side snowflake divertor geometries. Snowflake simulations using this method are presented, in addition to the first SOLPS-ITER simulation of the X-point target. Analysis of these results show reasonable consistency with the simple modelling and theoretical predictions, supporting the validity of the methodology implemented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Overset meshes for incompressible flows: On preserving accuracy of underlying discretizations

This study on overset meshes for incompressible-flow simulations is motivated by accurate prediction of wind farm aerodynamics involving large motions and deformations of components with complex geometry. Specifically, using first-order hyperbolic and elliptic equation proxies for the incompressible Navier-Stokes (NS) equations, we investigate the influence of information exchange between overset meshes on numerical performance where the underlying discretization is second-order accurate. The first aspect of information exchange surrounds interpolation of solution where we examine Lagrange and point-cloud-based interpolation for creating constraint equations between overset meshes. To maintain overall second-order accuracy, higher-order interpolation is required for elliptic problems, but linear interpolation is sufficient for hyperbolic problems in first-order form. Higher-order point-cloud-based interpolation provides a pathway to maintaining accuracy in unstructured meshes, but at higher complexity. The second aspect of information exchange focuses on comparing the approaches of overset single system (OSS) and overset Additive Schwarz (OAS) for coupling the linear systems of the overlapping meshes. While the former involves a single linear system, in the latter the discrete linear systems are solved separately, and solving the global system is accomplished through outer iterations and sequential information exchange in a Jacobi fashion. For the test cases studied, accuracy for hyperbolic problems is maintained by performing two outer iterations, whereas many outer iterations are required for elliptic systems. The order-of-accuracy studies explored here are critical for verifying the overset-mesh coupling algorithms used in engineering simulations. Accuracy of these simulations themselves is, however, quantified using engineering quantities of interest such as drag, power, etc. Consequently, we conclude with numerical experiments using NS equations for incompressible flows where we show that linear interpolation and few outer iterations are sufficient for achieving asymptotic convergence of engineering quantities of interest.

97 MATHEMATICS AND COMPUTING↗

Analysis of Needlet Internal Linear Combination performance on B -mode data from sub-orbital experiments

The observation of primordial B modes in cosmic microwave background (CMB) polarisation data represents the main scientific goal of most of the future CMB experiments. This signal is predicted to be much lower than polarised Galactic emission (foregrounds) in any region of the sky, pointing to the need for effective component separation methods. Aims. Among all the techniques, the blind Needlet Internal Linear Combination (NILC) is of great relevance given our current limited knowledge of the B-mode foregrounds. In this work, we explore the possibility of employing NILC for the analysis of B modes reconstructed from partial-sky data, specifically addressing the complications that such an application yields such as E–B leakage, needlet filtering, and beam convolution. We consider two complementary simulated datasets of future experiments: the balloon-borne Short Wavelength Instrument for the Polarisation Explorer (SWIPE) of the Large Scale Polarisation Explorer, which targets the observation of both reionisation and recombination peaks of the primordial CMB B-mode angular power spectrum, and the ground-based Small Aperture Telescope of Simons Observatory, which, instead, is designed to observe only the recombination bump at ℓ ~ 80. We assessed the performance of the following two alternative techniques to correct for the CMB E–B leakage: the recycling technique and the Zhao-Baskaran method. We find that both techniques reduce the E–B leakage residuals at a negligible level given the sensitivity of the considered experiments, except for the recycling method in the SWIPE footprint at ℓ < 20. Thus, we implemented two extensions of the pipeline, the iterative B decomposition and the diffusive inpainting, which enabled us to recover the input CMB B-mode power for ℓ ≥ 5. For the considered experiments, we demonstrate that needlet filtering and beam convolution do not affect the CMB B-mode reconstruction. Finally, with an appropriate masking strategy, we find that NILC foregrounds subtraction allows one to achieve sensitivities on the tensor-to-scalar ratio in agreement with the targets of the considered CMB experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

Subtask 1.5 – CO2 Injection Monitoring with an Optimized Scalable, Automated, Semipermanent Seismic Array

The scalable, automated, semipermanent seismic array (SASSA) method is a flexible and relatively cost-effective surface geophysical method for regular time-lapse monitoring of the movement of injected carbon dioxide (CO2) in a reservoir for CO2 enhanced oil recovery (EOR) or geologic CO2 storage operations. It has the advantages of a low-environmental-footprint while monitoring regions of a reservoir from the surface without the need for a regular grid distribution of receivers. Automated data collection is possible. As only time-lapse amplitude changes at the reservoir level due to CO2 movement within the reservoir are monitored, the turnaround time to deliver results from the SASSA method can be short, without the need for long, time-consuming data-processing workflows. As data is collected and processed, incremental information can be provided to the field operator. The Energy & Environmental Research Center (EERC) conducted a SASSA field test from September 2018 to November 2020 in a portion of the Bell Creek Field in Montana, which implemented new CO2 EOR field activities during the study period. Lessons learned from a proof-of-concept study were incorporated to improve the data quality of the SASSA method and demonstrate the viability of the technology. The EERC implemented several enhancements to improve data quality, including 1) an iterative survey design, which allowed placing the receivers in strategic locations where the movement of the CO2 in the reservoir could be tracked with minimum interference by the cultural noise in the study area; 2) the use of powerful seismic sources in the form of surface orbital vibrators, and 3) data acquisition during optimal periods. History-matched reservoir simulation was performed to predict gas saturation and pressure response induced by CO2 injection in the study area. The results were compared with the SASSA-measured responses to CO2 injection as a partial validation technique. A match between the two methods was observed for most of the SASSA points predicted to have intersected a CO2 saturation change. The validated results provide confidence that the SASSA method can be used independently as a CO2 saturation monitoring technique. As data are collected and processed, incremental information can be provided to the field operator. The critical components of the SASSA workflow for a successful application of the method are the following: Iterative survey design with information about CO2 injection activities from the oilfield operator. A detailed CO2 injection plan is the key driver to select the strategic monitoring location of the SASSA sensors. After this information is incorporated in the initial distribution of sources and receivers in the study area, high-resolution satellite images are used to identify ground locations not affected by cultural noise sources, such as power lines, pipelines/flow lines, or roadways. In the next iteration of the survey design, a scouting trip to the study area is needed to understand more details of the noise sources identified in the previous step and the intensity of the field activities that can also generate noise during the monitoring. Integrating the information from the scouting trip into the survey design to select the optimum source and receiver locations is the final step. Noise attenuation. The variety of noise types during seismic monitoring of an oil field is enormous. Tailored noise characterization and processing at a node-by-node level can enhance the performance and sensitivity of the SASSA technique. Future advancements that could improve the efficiency and application of the SASSA technology include: Gaining a better understanding of the noise field produced by the seismic source to aid the choice of receiver location. Surface noise from the source can overwhelm the small signal changes due to CO2 that the SASSA method measures. Improved data-processing workflow to automatically analyze and adapt to dynamic noise conditions associated with industrial settings. This subtask was funded through the EERC–DOE Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE- FE0024233.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Long-pulse high-performance H-mode plasmas achieved on EAST

A record duration of a 310 s H-mode plasma (H 98y2 ~ 1.3, n e /n GW ~ 0.7, f BS > 50%) has been recently achieved on experimental advanced superconducting tokamak (EAST) with metal walls, exploiting the device's improved long-pulse capabilities. The experiment demonstrates good control of tungsten concentration, core/edge MHD stability, and particle and heat exhaust with an ITER-like tungsten divertor and zero injected torque, establishing a milestone on the path to steady-state long-pulse high-performance scenarios in support of ITER and CFETR. Important synergistic effects are leveraged toward this result, which relies purely on radio frequency (RF) powers for heating and current drive (H&CD). On-axis electron cyclotron heating enhances the H&CD efficiency from lower hybrid wave injection, increasing confinement quality and enabling fully non-inductive operation at high density (n e /n GW ~ 70%) and high poloidal beta (β P ~ 2.5). A small-amplitude grassy edge localized mode regime facilitates the RF power coupling to the H-mode edge and reduces divertor sputtering/erosion. The high energy confinement quality (H 98y2 ~ 1.3) is achieved with the experimental and simulated results pointing to the strong effect of Shafranov shift on turbulence. Transport analysis suggests that trapped electron modes dominate in the core region during the record discharge. The detailed physics processes (RF synergy, core-edge integration, confinement properties, etc.) of the steady-state operation will be illustrated in the content. In the future, EAST will aim at accessing more relevant dimensionless parameters to develop long-pulse high-performance plasma toward ITER and CFETR steady-state advanced operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Towards Realistic and High Fidelity Models for Nuclear Reactor Power Synthesis Simulation with Self-Powered Neutron Detectors

As presented in this report, a weighting function–based inferencing method is being applied to synthesize the power distribution in next-generation and university research reactors based on simulated self power neutron detector (SPND) responses. The overall goal is to assess the impacts of sensor uncertainty and true power distribution perturbations on the error in the synthesized power distribution. Regarding sensor uncertainty, the NuScale Small Modular Reactor (SMR) and the Westinghouse AP1000 serve as testbeds for analyzing the impact of varying the sensor uncertainty, as well as varying the number of sensors per sensor string in the reactor core. The reactor models are informed by Monte Carlo N-Particle (MCNP) neutron flux tallies. For the NuScale SMR and Westinghouse AP1000, the SPND response functions (i.e., the response of the SPNDs to individual segments of fuel) were determined homogeneously. Regarding an analysis of power distribution perturbation detection, the Texas A&M Testing, Research, Isotopes, General Atomics Reactor (TAMU TRIGA) reactor was used as a demonstration case with one particular arrangement of SPNDs; the response functions for this reactor model were determined heterogeneously, making this a uniquely high-fidelity demonstration of perturbation detection. Finally, SPND models generated in the Geometry and Tracking 4 (Geant4) code have been generated and tested for comparison with traditionally implemented analytical SPND models, with the intent for Geant4 integration with the full methodological framework. SPND current outputs as a function of distance from some fuel assembly segment in the NuScale SMR are compared with the analytically determined currents. Results from the sensor uncertainty simulations for the NuScale SMR and AP1000 indicate that the average error in the inferred power distribution on the fuel assembly segment level is reasonably low, being slightly less than the random uncertainty applied to all respective SPNDs in both cores. For example, if all SPNDs in the core have a random uncertainty of 5%, then the corresponding fuel assembly segment level error (i.e. difference between the true and inferred local power) is ~2–3%. However, the maximum error in the inferred power distribution on the fuel assembly segment level can be considerably high (>15%) when SPND random uncertainties start to exceed ~3%. In general, the average and maximum errors in the inferred power distribution were slightly higher in the AP1000 as opposed to the NuScale SMR for the sensor string configurations considered herein. Another result determined from analysis of the sensor uncertainty simulations was that increasing the number of SPNDs per string does not clearly reduce inferred power distribution error and can in fact make the error large in some cases; however, this assessment may skewed due to imposed iteration limits. Results from the perturbation detection demonstration using the high-fidelity TAMU TRIGA model indicate that, given the arrangement of 17 SPND strings and 4 SPNDs per string considered herein, there is a clear, provable ability to infer a localized Gaussian-type peak perturbation in the 3D power distribution. Such a perturbation was detected with an average fuel assembly segment level error of 0.19%, and the general visualization of the detected perturbation clearly indicates that the magnitude and shape were appropriately resolved. Finally, the electrical current output generated by the Geant4 modeled SPND indicates significant magnitude differences than the analytically modeled SPND, demonstrating the need for accurate SPND models which account for finite sensor geometry effects to inform the power synthesis work described herein.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Alcator C-Mod: Research in support of ITER and steps beyond

This paper presents an overview of recent highlights from research on Alcator C-Mod. Significant progress has been made across all research areas over the last two years, with particular emphasis on divertor physics and power handling, plasma-material-interaction studies, ELM-suppressed pedestal dynamics, core transport and turbulence, and RF heating and current drive utilizing Ion Cyclotron and Lower Hybrid tools. Specific results of particular relevance to ITER include: inner wall SOL transport studies that have led, together with results from other experiments, to the change of the detailed shape of the inner wall in ITER; runaway electron studies showing that the critical electric field required for runaway generation is much higher than predicted from collisional theory; core tungsten impurity transport studies reveal that tungsten accumulation is naturally avoided in typical C-Mod conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗