Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ITER design”

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 649 records · Page 36

KSSOLV 2.0: An efficient MATLAB toolbox for solving the Kohn-Sham equations with plane-wave basis set

KSSOLV (Kohn-Sham Solver) is a MATLAB toolbox for performing Kohn-Sham density functional theory (DFT) calculations with a plane-wave basis set. KSSOLV 2.0 preserves the design features of the original KSSOLV software to allow users and developers to easily set up a problem and perform ground-state calculations as well as to prototype and test new algorithms. Furthermore, it includes new functionalities such as new iterative diagonalization algorithms, k-point sampling for electron band structures, geometry optimization and advanced algorithms for performing DFT calculations with local, semi-local, and hybrid exchange-correlation functionals. It can be used to study the electronic structures of both molecules and solids. So, we describe these new capabilities in this work through a few use cases. We also demonstrate the numerical accuracy and computational efficiency of KSSOLV on a variety of examples.

97 MATHEMATICS AND COMPUTING↗

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↗

Testing of Fusion Machine Components Against Cryogenic Pellet Impacts

Shattered pellet injection (SPI), currently the most effective method of disruption mitigation, is currently implemented on tokamaks worldwide for experimental purposes. Cryogenic pellets are formed and fired into an angled surface before entering the plasma. The impact with the angled surface causes the pellets to fragment into a cloud of particles with the purpose of increasing the surface area for ablation. As pellets traverse guide tubes, depending on design, there is a chance of an off-normal pellet impact. Pellet impacts are also, depending on design, possible in the plasma chamber if not fully ablated and assimilated, or if the fragment plume is not directed in the proper direction. This paper outlines a series of pellet impact tests on various tiles and components relevant to the ITER, JET, and ASDEX Upgrade SPI systems. In conclusion, testing was done to assess the potential damage from pellet and fragment plume impacts through high-speed imaging and the visual inspection of components.

Cryogenic pellets↗

Enhancing Gaussian Process Surrogates for Optimization and Posterior Approximation via Random Exploration

This paper proposes novel noise-free Bayesian optimization strategies that rely on a random exploration step to enhance the accuracy of Gaussian process surrogate models. The new algorithms retain the ease of implementation of the classical GP-UCB algorithm, but the additional random exploration step accelerates their convergence, nearly achieving the optimal convergence rate. Furthermore, to facilitate Bayesian inference with intractable likelihoods, we propose to utilize optimization iterates for maximum a posteriori estimation to build a Gaussian process surrogate model for the unnormalized log-posterior density. We provide bounds for the Hellinger distance between the true and the approximate posterior distributions in terms of the number of design points. We demonstrate the effectiveness of our Bayesian optimization algorithms in nonconvex benchmark objective functions, in a machine learning hyperparameter tuning problem, and in a black-box engineering design problem. The effectiveness of our posterior approximation approach is demonstrated in two Bayesian inference problems for parameters of dynamical systems.

Bayesian inference↗

Motion tomography via occupation kernels

The goal of motion tomography is to recover a description of a vector flow field using measurements along the trajectory of a sensing unit. In this paper, we develop a predictor corrector algorithm designed to recover vector flow fields from trajectory data with the use of occupation kernels developed by Rosenfeld et al. [9,10]. Specifically, we use the occupation kernels as an adaptive basis; that is, the trajectories defining our occupation kernels are iteratively updated to improve the estimation in the next stage. Initial estimates are established, then under mild assumptions, such as relatively straight trajectories, convergence is proven using the Contraction Mapping Theorem. We then compare the developed method with the established method by Chang et al. [5] by defining a set of error metrics. Here, we found that for simulated data, where a ground truth is available, our method offers a marked improvement over [5]. For a real-world example, where ground truth is not available, our results are similar results to the established method.

97 MATHEMATICS AND COMPUTING↗

The Relationship Between Dose Rate and Decay Heat for Spent Nuclear Fuel Casks

Decay heat and dose rate are two important limits used for determining the allowable contents of spent fuel (SNF) in dry storage systems and transportation packages. While the decay heat limit is used to maintain fuel cladding integrity and ensure retrievability, dose rates are used to demonstrate compliance with regulatory requirements on radiation protection. Because both dose rate and decay heat result from decay of radioisotopes in SNF, this study is an attempt to examine the relationship between dose rate and decay heat for a given cask design. Dose rates were evaluated for 198 cask configurations, that include various SNF system designs (e.g., storage, transfer, transport), SNF characteristics (e.g., fuel types, burnup, cooling time), and loading maps (e.g., uniform loading, zone loading), while a constant decay heat was maintained. The decay heat was calculated using US Nuclear Regulatory Commission (NRC) Regulatory Guide (RG) 3.54, Revision 2, and verified using ORIGEN sequence within SCALE code system. The ORIGEN outputs were used as source terms in the dose analysis using 198 different configurations. A computer script was developed to calculate the cooling time necessary to achieve a given decay heat for a given enrichment, assembly average burnup, assembly mass, and in-core history using a rootfinder algorithm. Initially the computer script was developed to provide cooling time and burnup calculations directly to the analysis of dose rates and decay heats, so a comparison between Used Nuclear Fuel-Storage Transportation and Disposal Analysis Resource Data System (UNF-ST&DARDS) results and RG3.54 data was made; results are included in the appendix to this document. However, an iterative approach was used to compute cooling time, and its accuracy did not depend on the results of the RG3.54r2 algorithm, although the algorithm was still used. The results of the evaluation presented herein clearly demonstrate that a given decay heat does not correspond to a unique dose rate for a variety of cask and package designs. There is no clear pattern to develop a correlation between decay heat and the source terms. Depending on burnup, enrichment, cask type, and loading pattern, dose rates varied for the exact same decay heat—in some cases by 400% for a given cask. For cases in which decay heat was held constant through selection of the appropriate cooling time, dose rates would decrease with increasing burnup, and in other cases, dose rates would increase. The large variation in dose rates for a constant decay heat indicates that casks loaded based on decay heat—that is allowing any burnup, cooling time, and enrichment combinations that yield the qualified decay heat limit(s) —cannot ensure that an Independent Spent Fuel Storage Installation or a spent fuel transportation package will meet the regulatory limits set forth by the respective regulations, i.e., 10 CFR 72 or 10 CFR 71.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Hybrid real- and reciprocal-space full-field imaging with coherent illumination

In this work, we present a novel diffractive imaging method that harnesses a low-resolution real-space image to guide the phase retrieval. A computational algorithm is developed to utilise such prior knowledge as a real-space constraint in the iterative phase retrieval procedure. Numerical simulations and proof-of-concept experiments are carried out, demonstrating our method's capability of reconstructing high-resolution details that are otherwise inaccessible with traditional phasing algorithms. With the present method, we formulate a conceptual design for the coherent imaging experiments at a next-generation x-ray light source.

42 ENGINEERING↗

Risk-Informed Applications and Science – Past, Present, & Future

This is a presentation to be made at the University of Tokyo. The Idaho National Laboratory is creating the next-generation reliability- and risk-assessment methods and tools that support risk-informed decision-making by combining physics-based models with probabilistic quantification approaches. Integrating these two worlds of physics and probability using a simulation framework leads us to predictions based upon an approach called “computational risk assessment” which will serve as the technical basis for the future of reliability and risk approaches. The driving factors for this new approach includes: temporal (timing issues), spatial (location issues), mechanistic (physics issues), and topological (complexity issues). By combining phenomenology directly with stochastic quantification, we are able to perform advanced uncertainty analysis directly on both parameters and models. In addition, new tools allow for the creation of reduced order models that mimic high-fidelity engineering models while still permitting the realization of thousands of simulation iterations in short run-times on current workstation-class computers. While these advanced methods and tools can provide increased realism in our engineering safety and risk approaches, their greater benefit is to provide a risk-informed engineering framework for design and operation. This talk will briefly review the history of risk assessment and risk-informed applications, describe some of the current research and development found at the Idaho National Laboratory, and discuss potential future applications and approaches for advanced methods and tools.

97 MATHEMATICS AND COMPUTING↗

Flexible and Modular Simultaneous Modeling of Flow and Reactive Transport in Rivers and Hyporheic Zones

Investigations of coupled multiphysics processes in rivers and hyporheic zones have extensively used numerical models. Most existing models use a sequential, one-way coupling between the surface and subsurface domains. Such one-way coupling potentially introduces error. To overcome this, a fully coupled model, hyporheicFoam, was developed using the open-source computational platform OpenFOAM. It captures the coupled flow and multicomponent reactive transport processes within both surface and subsurface domains and across their interface. The coupling between two domains is implemented by mapping conservative flux boundary conditions at the interface through an iterative algorithm. Reactive transport is enabled by specifying a reaction network. To start, we have implemented reaction kinetics following the double Monod-type model with inhibition. The model capability is illustrated through modeling of both conservative and reactive hyporheic flow and transport through dune bedforms. With the novel coupled model, it is now possible to quantify reactions wherein the reactants and products are constantly exchanging between domains and have feedbacks. hyporheicFoam can simulate large, three-dimensional cases owing to the computational flexibility and power offered by the code structure and parallel design of OpenFOAM.

58 GEOSCIENCES↗

Advances in the long-pulse steady-state high beta H-mode scenario with active controls of divertor heat and particle fluxes in EAST

Since the last IAEA-Fusion Energy Conference, the Experimental Advanced Superconducting Tokamak (EAST) research program has been, in support of ITER and CFETR, focused on development in terms of the long-pulse steady-state (fully noninductive) high beta H-mode scenario with active controls of the stationary and transient divertor heat and particle fluxes. The operational domain of the steady-state H-mode plasma scenario has been significantly extended with ITER-like tungsten mono-block divertor, plasma control and heating schemes. EAST has achieved several important milestones in the development of high β p H-mode scenario and its key physics and technologies. A 60 s-scale long-pulse steady-state high β p H-mode discharge with the major normalized plasma parameters similar to the designed performance of the CFETR 1GW fusion power operation scenario has been successfully established and sustained by pure RF heating and current drive. Several feedback control schemes have been developed for a sustained detachment with good core confinement. This includes control of the total radiation power, target electron temperature, and particle flux measured using divertor Langmuir probes or a combination of the control of target electron temperature and AXUV radiation near the X point. The detachment feedback control schemes have been integrated with small-ELM regimes and high β p scenario via neon seeding, enabling a core and edge compatible integrated high-beta scenario applicable to long-pulse operations. ELM suppression has been achieved using various methods, including resonant magnetic perturbations and impurity seeding. Full suppression of ELMs by using n=4 RMPs has been demonstrated for ITER for the first time in low input torque plasmas in EAST. Finally, EAST has been operated with helium to support the ITER research requirements for the first time. For a long-pulse, high bootstrap current fraction operation, a new lower tungsten divertor with active water-cooling has been installed, along with improvements in the heating and current drive capability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Schwarz alternating method for transient solid dynamics

Abstract In our earlier work, we formulated the Schwarz alternating method as a means for concurrent multiscale coupling in finite deformation solid mechanics for quasi‐static problems. Herein, we advance this method for the study of transient dynamic multiscale solid mechanics problems where information is exchanged back and forth between small and large scales. The extension to dynamics relies on the notion of a global time stepper. Within each global time step, the subdomains are coupled by the standard Schwarz iterative process. Remarkably, each subdomain can use its own time step or even its own time integrator to advance its solution in time, provided that they synchronize at each global time step. We study the performance of the Schwarz method on several examples designed for this purpose. Our numerical experiments demonstrate that the method is capable of coupling regions with different mesh resolutions, different element types, and different time integration schemes (e.g., implicit and explicit), all without introducing any artifacts that afflict other coupling methods for transient dynamics. Finally, we apply the dynamic Schwarz alternating method to the simulation of a bolted joint subjected to dynamic loading, as a demonstration of the performance of the method in a realistic scenario.

Mota, Alejandro↗

Rapid 3D nanoscale coherent imaging via physics-aware deep learning

Phase retrieval, the problem of recovering lost phase information from measured intensity alone, is an inverse problem that is widely faced in various imaging modalities ranging from astronomy to nanoscale imaging. The current process of phase recovery is iterative in nature. As a result, the image formation is time consuming and computationally expensive, precluding real-time imaging. Here, we use 3D nanoscale X-ray imaging as a representative example to develop a deep learning model to address this phase retrieval problem. We introduce 3D-CDI-NN, a deep convolutional neural network and differential programing framework trained to predict 3D structure and strain, solely from input 3D X-ray coherent scattering data. Our networks are designed to be “physics-aware” in multiple aspects; in that the physics of the X-ray scattering process is explicitly enforced in the training of the network, and the training data are drawn from atomistic simulations that are representative of the physics of the material. We further refine the neural network prediction through a physics-based optimization procedure to enable maximum accuracy at lowest computational cost. 3D-CDI-NN can invert a 3D coherent diffraction pattern to real-space structure and strain hundreds of times faster than traditional iterative phase retrieval methods. Our integrated machine learning and differential programing solution to the phase retrieval problem is broadly applicable across inverse problems in other application areas.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Coupled Charge and Radiation Transport Processes in Thermophotovoltaic and Thermoradiative Cells

Accurate modeling of charge transport and both thermal and luminescent radiation is crucial to the understanding and design of radiative thermal energy converters. Charge-carrier dynamics in semiconductors are well-described by the Poisson-drift-diffusion equations, and thermal radiation in emitter-absorber structures can be computed using multilayer fluctuational electrodynamics. These two types of energy flows interact through radiation absorption or luminescence and charge-carrier generation or recombination. However, past research has typically only assumed limited interaction, with thermal radiation absorption as an input for charge-carrier models to predict device performance. To examine this assumption, we develop a fully coupled iterative model of charge and radiation transport in semiconductor devices, and we use our model to analyze near-field and far-field GaSb thermophotovoltaic and thermoradiative systems. By comparing our results with past methods that do not consider cross-influences between charge and radiation transport, we find that a fully coupled approach is necessary to accurately model photon recycling and near-field enhancement of external luminescence. As these effects can substantially alter device performance, our modeling approach can aid in the design of efficient thermophotovoltaic and thermoradiative systems.

30 DIRECT ENERGY CONVERSION↗

Evaluation Toolkit for Technical Assistance Programs

Across the United States, every home, business, industry, and government depends on abundant, reliable, and affordable electricity. Energy stakeholders have many options to improve how energy is generated, distributed, and used, but choosing the right path can be complex and challenging. To support more informed decision-making about local electricity systems, the U.S. Department of Energy (DOE) and its national laboratories provide customized technical assistance (TA) through several different programs. These TA programs are designed to support a wide range of stakeholders with a variety of needs. Technical assistance may include brief consultations with subject-matter experts, in-depth technical modeling and analysis, stakeholder engagement, and peer-to-peer exchange. Delivering effective TA is an iterative process that requires evaluation and adaptation. A comprehensive and robust evaluation framework provides the structure needed to ensure that TA programs and practitioners remain effective, responsive to industry trends, and aligned with local needs. This toolkit provides a framework—including a clear, structured process—that TA program staff can adapt to their program's goals to evaluate success.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Innovative laser-based methods for monitoring fuel retention in ITER

This paper addresses the challenge of tritium inventory management in ITER and future fusion reactors, highlighting the importance of accurate tritium measurement and its spatial distribution within the vacuum vessel. Given ITER’s operational constraints, especially the limit on tritium retention, precise measurement is essential for both safety and regulatory compliance. To tackle these questions, the paper presents the T-monitor diagnostic system developed by Forschungszentrum Jülich, which uses Laser-Induced Desorption (LID) in combination with Diagnostic Residual Gas Analysis (DRGA) to measure hydrogen isotope concentrations on the surface of divertor tiles. The system integrates a high-power laser, advanced optical components, and a Fast Scanning Mirror Unit (FSMU) for accurate laser spot positioning with rapid response. Designed to measure in situ tritium retention, the diagnostic provides high-resolution spatial mapping, vital for evaluating detritiation strategies. The laser heating process increases the divertor surface temperature to 1600 K within the laser spot, promoting hydrogen isotope desorption. Accurate measurements require the precise control of laser parameters, including pulse duration and spot size, with a target relative accuracy of 20%. The optical design includes both in-vessel and ex-vessel components, such as durable high-reflectivity mirrors made of gold and copper, selected not only for their infrared performance but also for their transmission of visible wavelengths for observation purposes. To protect optical components from contamination, a pneumatic shutter is used.

Hydrogen isotopes↗

On the Sampling-Based Computation of Nash Equilibria Under Uncertainty via the Nikaido–Isoda Function

We consider the computation of an equilibrium of a stochastic Nash equilibrium problem, where the player objectives are assumed to be L 0 -Lipschitz continuous and convex, given rival decisions with convex and closed player-specific feasibility sets. To address this problem, we consider minimizing a suitably defined value function defined using the Nikaido–Isoda function. Such an avenue does not necessitate either monotonicity properties of the concatenated gradient map or potentiality requirements on the game but does require a suitable regularity requirement under which a stationary point is a Nash equilibrium. We design and analyze a sampling-enabled projected-gradient-response method, reliant on inexact resolution of a player-level best-response subproblem. Here, by deriving suitable Lipschitzian guarantees on the value function, we derive both asymptotic guarantees for the sequence of generated iterates as well as rate and complexity guarantees for computing a stationary point by appropriate choices of the sampling rate and inexactness sequence.

Nikaido-Isoda function↗

Deep Reinforcement Learning Based Volt-VAR Optimization in Smart Distribution Systems

This paper develops a model-free volt-VAR optimization (VVO) algorithm via multi-agent deep reinforcement learning (DRL) in unbalanced distribution systems. This method is novel since we cast the VVO problem in distribution networks to an intelligent deep Q-network (DQN) framework, which avoids solving a specific optimization model directly when facing time-varying operating conditions in the systems. We consider statuses/ratios of switchable capacitors, voltage regulators, and smart inverters installed at distributed generators as the action variables of the agents. A delicately designed reward function guides these agents to interact with the distribution system, in the direction of reinforcing voltage regulation and power loss reduction simultaneously. The forward-backward sweep method for radial three-phase distribution systems provides accurate power flow results within a few iterations to the DRL environment. The proposed method realizes the dual goals for VVO. We test this algorithm on the unbalanced IEEE 13-bus and 123-bus systems. Numerical simulations validate the excellent performance of this method in voltage regulation and power loss reduction.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Coherent diffraction imaging in the undergraduate laboratory

We present an undergraduate optics instructional laboratory designed to teach skills relevant to a broad range of modern scientific and technical careers. In this laboratory project, students image a custom aperture using coherent diffraction imaging, while learning principles and skills related to digital image processing and computational imaging, including multidimensional Fourier analysis, iterative phase retrieval, noise reduction, finite dynamic range, and sampling considerations. After briefly reviewing these imaging principles, we describe the required experimental materials and setup for this project. Our experimental apparatus is both inexpensive and portable, and a software application we developed for interactive data analysis is freely available.

Porter, J. Nicholas↗