Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “field optimization”

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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Characterizing the Impact of Heliostat Size, Focus Method, and Optical Error on Concentrating Solar Tower Systems

Concentrating Solar Power (CSP) tower systems use heliostat fields to direct solar energy to a central receiver. The efficiency of the heliostat field is crucial for achieving thermal generation targets and is also the largest cost component, making up 30-40% of total project costs. Developers face a tradeoff: reducing heliostat costs may decrease optical performance and increase land use, while improving performance could raise manufacturing costs (through higher precision). The optimal design depends on land, labor, and manufacturing costs. Additionally, inherent performance limitations exist - optical errors result in larger solar images as heliostats move farther from the receiver, increasing spillage and attenuation losses. Ultimately, there exists a point where receiver design power can no longer be achieved at design conditions for a given solar field area. This point depends on key heliostat field design parameters including field layout, heliostats size, focus method, and heliostat optical error. In this work, we explore and characterize the performance trade-offs between heliostat size, focus method, and optical error using Monte Carlo ray-tracing simulations using NREL's high performance computing resources. The goal of this fundamental investigation is to understand practical heliostat design performance impacts that can be used to produce more cost-effective heliostat fields while still achieving plant thermal generation targets.

14 SOLAR ENERGY↗
DOE OSTI2025

Dopant Optimization of Donors in Semiconductor Opening Switches to Eliminate Prepulse

Semiconductor opening switches are solid-state devices capable of delivering nanosecond, hundreds of kilovolts pulses by interrupting kiloamps of current. The interruption of the current occurs in a moderately doped p-region when a high electric field region (HFR) is formed. The HFR occurs because the reverse pumping current cannot be supported by the saturation velocity and majority carrier concentration of the doping level. However, the donor profile also significantly affects the pulse performance. A secondary prepulse occurs if a secondary HFR is formed at the interface of the background n-doping and N+ doping (X n ) . By moving the location of X n deeper into the diode, the effect of the prepulse is reduced. This article investigates the effect of the donor doping profile on the performance metrics of semiconductor opening switches through technology computer-aided design (TCAD) simulations and experimental results. Through a SILVACO TCAD optimization, we designed a P + /p/n - base/n/N + where the intersection of the moderate p-region and intrinsic n-base region (X p ) is at 160 μm and X n is at 220 μm. This profile is fabricated via silicon epitaxy. Experimentally, it is shown that a deep X n (220 μm) compared with a shallow X n (300 μm) reduces the rise time by >5× . In addition, the magnitude of current density during interruption affects the prepulse foot and pulse shape. At lower current densities without the graded donor profile, high peak voltages are not achieved. Comparing the experimental results to the TCAD simulations shows that the model is predictive under high-current densities in the semiconductor opening switch (SOS) regime.

nanosecond pulse power↗
DOE OSTI2025

The Effect of Hydrothermal Alteration and Microcracks on Hydraulic Properties and Poroelastic Deformation: A Case Study of the Blue Mountain Geothermal Field

Abstract Geothermal energy plays a vital role in decarbonizing electricity and heat supply. Effective utilization of geothermal resources hinges on identifying or generating permeable reservoir zones and understanding how effective pressure variations affect fluid circulation and reservoir properties by poroelastic deformation. Hydrothermal alteration can modify the petrophysical properties of geothermal reservoir rocks, which may increase or decrease its productivity. Understanding these alteration effects is essential to predict and optimize long‐term sustainable geothermal operations. Here, we investigate the impact of hydrothermal alteration on poroelastic and hydraulic properties of diverse lithologies in a series of deformation tests performed at several confining (0–80 MPa) and pore pressure (10–30 MPa) levels. Experimental results of hydrothermally altered dikes and phyllites obtained from the Blue Mountain geothermal field (Nevada, USA) are compared to thermally cracked La Peyratte granite (France) and correlated with petrophysical properties, mineral composition, and microstructures. Argillic alteration of dikes increases porosity and storage capacity but lowers thermal conductivity and increases pore compressibility. Conversely, silicate precipitation in phyllites increases stiffness and thermal conductivity but also reduces porosity and permeability. Experimentally determined effective pressure coefficients range from 0.1 to 0.9, differ for permeability and volumetric strain and decrease with increasing effective pressure. The presence of compliant microcracks and crack‐like pores significantly increases the stress sensitivity of La Peyratte granite and silicified phyllites. This study demonstrates how thermal and chemical alteration impacts poromechanical and petrophysical characteristics of geothermal targets, which ultimately govern reservoir stability and subsidence, induced seismicity as well as fluid and heat extraction efficiency during geothermal operations.

Schuster, Valerian [Helmholtz Centre Potsdam GFZ G↗
DOE OSTI2024

Epitaxial Growth of Large-Scale 2D CrTe 2 Films on Amorphous Silicon Wafers With Low Thermal Budget

2D van der Waals (vdW) magnets open landmark horizons in the development of innovative spintronic device architectures. However, their fabrication with large scale poses challenges due to high synthesis temperatures (>500 °C) and difficulties in integrating them with standard complementary metal-oxide semiconductor (CMOS) technology on amorphous substrates such as silicon oxide (SiO 2 ) and silicon nitride (SiN x ). Here, a seeded growth technique for crystallizing CrTe 2 films on amorphous SiN x /Si and SiO 2 /Si substrates with a low thermal budget is presented. This fabrication process optimizes large-scale, granular atomic layers on amorphous substrates, yielding a substantial coercivity of 11.5 kilo-oersted, attributed to weak intergranular exchange coupling. Field-driven Néel-type stripe domain dynamics explain the amplified coercivity. Moreover, the granular CrTe 2 devices on Si wafers display significantly enhanced magnetoresistance, more than doubling that of single-crystalline counterparts. Current-assisted magnetization switching, enabled by a substantial spin–orbit torque with a large spin Hall angle (85) and spin Hall conductivity (1.02 × 10 7 ℏ/2e Ω -1 m -1 ), is also demonstrated. These observations underscore the proficiency in manipulating crystallinity within integrated 2D magnetic films on Si wafers, paving the way for large-scale batch manufacturing of practical magnetoelectronic and spintronic devices, heralding a new era of technological innovation.

2D magnetic thin films↗
DOE OSTI2025

A digital twin platform for building performance monitoring and optimization: Performance simulation and case studies

Advancements in sensor technology, data analytics, affordable compute, and communication infrastructure have paved the way for Digital Twin technology in optimizing building operations and controls. This study presents the development of an open and interoperable web-based Digital Twin platform for integrating diverse data streams and facilitating effective user interactions. The platform utilizes modern technologies for the web framework and time-series data management, ensuring scalability and responsiveness. The backend supports seamless integration of diverse data sources and emulators, incorporating data from building sensors and meters, external weather Application Programming Interfaces, and advanced EnergyPlus simulation models of the building and its energy systems including the Distributed Energy Resources that are formulated in Functional Mockup Units. A simulation case study was conducted with FlexLab, a test facility on Lawrence Berkeley National Laboratory campus. The case study includes normal operations, Distributed Energy Resource integration, and power outage scenarios, to illustrate the Digital Twin’s ability to provide critical insights into energy performance and thermal resilience. The results demonstrated the platform’s potential as a decision-support tool for optimizing building energy performance and enhancing resilience against extreme weather events. Future work will focus on deploying the Digital Twin platform to a real building for field validation, extending its capabilities to cover more scenarios such as bidirectional Electric Vehicle interactions, and enhancing user engagement.

EnergyPlus↗
DOE OSTI2025

Vector-level feedforward control of LPBF melt pool area using a physics-based thermal model

Laser powder bed fusion (LPBF) is an additive manufacturing technique that has gained popularity thanks to its ability to produce geometrically complex, fully dense metal parts. However, these parts are prone to internal defects and geometric inaccuracies, stemming in part from variations in the melt pool. Here, this paper proposes a novel vector-level feedforward control framework for regulating melt pool area in LPBF. By decoupling part-scale thermal behavior from small-scale melt pool physics, the controller provides a scale-agnostic prediction of melt pool area and efficient optimization over it. This is done by operating on two coupled lightweight models: a finite-difference thermal model that efficiently captures vector-level temperature fields and a reduced-order, analytical melt pool model. Each model is calibrated separately with minimal single-track and 2D experiments, and the framework is validated on a complex 3D geometry in both Inconel 718 and 316L stainless steel. Results showed that feedforward vector-level laser power scheduling reduced geometric inaccuracy in key dimensions by 62%, overall porosity by 16.5%, and photodiode root-mean-squared deviation by 38.5% on average. Overall, this modular, data-efficient approach demonstrates that proactively compensating for known thermal effects can significantly improve part quality while remaining computationally efficient and readily extensible to other materials and machines.

Additive manufacturing↗
DOE OSTI2025

Characterization of dangling bond defects at the crystalline Si/SiO x interface in a polycrystalline Si passivating contact solar cell at room temperature with electrically detected magnetic resonance spectroscopy

Monocrystalline silicon solar cells can achieve photoconversion efficiencies exceeding 26%; however, performance-limiting defects that trap carriers continue to be a challenge. In this work, we have characterized Si solar cells with tunneling SiO x /polycrystalline-Si (poly-Si) passivating contacts (TOPCon) on As-doped Czochralski Si wafers with electrically detected magnetic resonance (EDMR) spectroscopy. We fabricated 2 × 20 mm 2 TOPCon-like mini solar cells with edge passivation alongside larger 4 cm 2 sister cells and obtained similar device characteristics. We performed EDMR spectroscopy at 300 K on two minicells with different degrees of surface passivation based on the recombination parameter, J o , values of 40 and 310 fA/cm2. We optimized the resolution and the signal-to-noise ratio of the EDMR response of the minicells by varying the forward bias voltage and the magnetic field modulation amplitude. We detect two distinct signals with EDMR spectroscopy, an axial-like signal at g = 2.009, 2.0087, and 2.0015, and an isotropic signal at g = 2.0024, which we attribute to Si dangling bonds (P b0 and P b centers) and boron–oxygen related defects, respectively, at or near the c-Si/SiO x interface. The EDMR signals were lower for the cell with a lower value of J o , while the ratio of the two defect populations was very similar. The EDMR signal increases with forward bias but drops to zero at bias voltages >0.5 V, consistent with interface defects within or near the boron-doped emitter depletion region. Our study demonstrates a method to fabricate minicells that can be characterized with EDMR spectroscopy to detect industrially relevant defects in TOPCon cells.

14 SOLAR ENERGY↗
DOE OSTI2025

SymProp: Scaling Sparse Symmetric Tucker Decomposition via Symmetry Propagation

Sparse symmetric tensors are an important class of tensors, and their decompositions serve as powerful tools for revealing low-rank structures. This paper introduces SymProp, a novel approach for scaling sparse symmetric Tucker decomposition by propagating symmetry through intermediate computations. SymProp optimizes two key computational kernels: Sparse Symmetric Tensor Times Same Matrix chain (S3 TTMc) for Higher-Order Orthogonal Iteration (HOOI) and Sparse Symmetric Tensor Times Same Matrix chain Times Core (S3 TTMcTC) for Higher-Order QR Iteration (HOQRI). Our method employs a metaprogramming-based index iteration approach to efficiently handle the upper triangular parts of intermediate dense symmetric tensors. SymProp achieves up to 50.9× speedup over SPLATT and up to 360.8× over Compressed Sparse Symmetric (CSS) format on the S3 TTMc operation. Moreover, our S3 TTMc and S3 TTMcTC implementations support tensor orders four levels higher than state-of-the-art methods. Our HOQRI demonstrates superior scalability and up to a 33.6× speedup over optimized HOOI. By enabling more scalable Tucker decompositions for higher orders, decomposition ranks, and dimension sizes, SymProp opens new possibilities for analyzing complex hypergraph structures in fields such as network science, data mining, and machine learning.

Li, Zecheng [North Carolina State University]↗
DOE OSTI2023

Improving Signal-to-Noise Ratio (SNR) for Readout Signals Using Adaptive Filters on Reconfigurable Controls Hardware

This study investigates the optimization of Signal-to-Noise Ratio (SNR) in superconducting quantum computing readout signals through adaptive filtering. Quantum computing technology has the potential to revolutionize various fields by delivering exponential speedup in solving certain computational problems. However, the technology's practical implementation is hindered by the difficulty of extracting clean, reliable signals during the readout phase, with various sources of noise presenting a significant barrier to clean signals. This noise, often present in readout profiles due to imperfect isolation, degrades the system's overall SNR, thus impeding the ability to extract the quantum state accurately. The research leverages the power of adaptive filtering to improve the SNR of quantum computing readout signals. Specifically, an adaptive filter is implemented in a PYNQ overlay on an FPGA, and eventually will be connected to a quantum computing system. The system models the noise with a Least Mean Squares (LMS) adaptive filter, and then subtracts the estimated noise from the received signal to improve the SNR. A Direct Memory Access (DMA) channel is used to handle the signal processing, delivering efficient, high-speed data transfer between the PYNQ system and the hardware. The study explores the benefits of this adaptive filtering technique, potentially providing a significant contribution to practical and fast quantum computing.

Johnson, Hans↗
DOE OSTI2024

Projective Integral Updates for High-Dimensional Variational Inference

Variational inference is an approximation framework for Bayesian inference that seeks to improve quantified uncertainty in predictions by optimizing a simplified distribution over parameters to stand in for the full posterior. Capturing model variations that remain consistent with training data enables more robust predictions by reducing parameter sensitivity. This work introduces a fixed-point optimization for variational inference that is applicable when every feasible log density can be expressed as a linear combination of functions from a given basis. In such cases, the optimizer becomes a fixed-point of projective integral updates. When the basis spans univariate quadratics in each parameter, the feasible distributions are Gaussian mean-fields and the projective integral updates yield quasi-Newton variational Bayes (QNVB). Other bases and updates are also possible. Since these updates require high-dimensional integration, this work begins by proposing an efficient quasirandom sequence of quadratures for mean-field distributions. Each iterate of the sequence contains two evaluation points that combine to correctly integrate all univariate quadratic functions and, if the mean-field factors are symmetric, all univariate cubics. More importantly, averaging results over short subsequences achieves periodic exactness on a much larger space of multivariate polynomials of quadratic total degree. The corresponding variational updates require four loss evaluations with standard (not second-order) backpropagation to eliminate error terms from over half of all multivariate quadratic basis functions. Furthermore, this integration technique is motivated by first proposing stochastic blocked mean-field quadratures, which may be useful in other contexts. A PyTorch implementation of QNVB allows for better control over model uncertainty during training than competing methods. Experiments demonstrate superior generalizability for multiple learning problems and architectures.

Gaussian mean-field↗
DOE OSTI2026

PyJMAK: An Open-Source Python Toolkit for Modeling Solid-State Metallurgical Phase Transformations

Accurate prediction of metallurgical phase transformations is an essential basis for autonomous optimization and rapid part qualification. Several methods can be used to estimate the evolution of phase fractions such as JMAK kinetics-based models, phase-field models, thermodynamic models, and data-driven machine learning models. Thermodynamic and phase-field-based methodologies solve multiphysics equations requiring numerous calibration parameters and significant computational resources. As a result, the computation domain is limited to a point or on order of micron-meters. The data-driven models rely on large datasets from experiments and simulations. While the JMAK model only provides information about phase fraction evolution, it can predict this evolution in near real-time using thermal history and thermodynamic data without restriction on the domain. JMAK models have been popularly used by researchers to model phase transformations occuring during additive manufacturing or over arbitrary temperature profiles. Commercial proprietary software such as Abaqus and Ansys or closed-source in-house implementations offer the ability to model JMAK based kinetics to predict phase transformation. However, these software packages are not open-source or freely available for use and development in conjunction with manufacturing machines, sensors, and machine learning algorithms. In addition, the use of the model is restricted by a license token. In contrast, given temperature profiles at multiple points in the domain, this Python-based PyJMAK model can compute phase evolution in parallel due to its stand-alone modular, voxel-based structure, and it can be executed on high-performance computing resources without any license restrictions.

Prabhune, Bhagya [Oak Ridge National Laboratory (O↗
DOE OSTI2020

Learning to Solve Large-Scale Security-Constrained Unit Commitment Problems

Security-constrained unit commitment (SCUC) is a fundamental problem in power systems and electricity markets. In practical settings, SCUC is repeatedly solved via mixed-integer linear programming (MIP), sometimes multiple times per day, with only minor changes in input data. In this work, we propose a number of machine learning techniques to effectively extract information from previously solved instances in order to significantly improve the computational performance of MIP solvers when solving similar instances in the future. Based on statistical data, we predict redundant constraints in the formulation, good initial feasible solutions, and affine subspaces where the optimal solution is likely to lie, leading to a significant reduction in problem size. Computational results on a diverse set of realistic and large-scale instances show that using the proposed techniques, SCUC can be solved on average 4.3 times faster with optimality guarantees and 10.2 times faster without optimality guarantees, with no observed reduction in solution quality. Out-of-distribution experiments provide evidence that the method is somewhat robust against data-set shift. Summary of Contribution. The paper describes a novel computational method, based on a combination of mixed-integer linear programming (MILP) and machine learning (ML), to solve a challenging and fundamental optimization problem in the energy sector. The method advances the state-of-the-art, not only for this particular problem, but also, more generally, in solving discrete optimization problems via ML. We expect that the techniques presented can be readily used by practitioners in the energy sector and adapted, by researchers in other fields, to other challenging operations research problems that are solved routinely.

Machine Learning↗
DOE OSTI2024

Final cooling with thick wedges for a muon collider

In the final cooling stages for a muon collider, the transverse emittances are reduced while the longitudinal emittance is allowed to increase. In previous studies, Final 4-D cooling used absorbers within very high field solenoids to cool low-momentum muons. Simulations of the systems did not reach the desired cooling design goals. In this study, we develop and optimize a different conceptual design for the final 4D cooling channel, which is based on using dense wedge absorbers. We used G4Beamline to simulate the channel and Python to generate and analyze particle distributions. We optimized the design parameters of the cooling channel and produced conceptual designs (corresponding to possible starting points for the input beam) which achieve transverse cooling in both x and y by a factor of ~3.5. These channels achieve a lower transverse and longitudinal emittance than the best design previously published.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

INVESTIGATION OF THE MULTILAYER SHIELDING EFFECT THROUGH NbTiN-AlN COATED BULK NIOBIUM

Superconducting radio frequency (SRF) cavities, composed of bulk niobium (Nb), are extensively utilized in particle accelerators. To optimize their performance and extend their operational range, thin films of superconducting materials are applied to the Nb surfaces. This study focuses on investigating the magnetic field penetration characteristics of NbTiN-AlN coating on bulk niobium. The experimental setup involves measuring the magnetic field penetration across flat superconducting samples at cryogenic temperatures. One key parameter of interest is the field just before the first flux penetration, Bp which enables the observation of the multilayer shielding effect on bulk niobium under high magnetic fields. The study reveals a significant enhancement in Bp of the NbTiN-AlN coated Nb samples compared to bare Nb samples. This property holds critical importance for achieving high field performance in SRF cavities. Furthermore, the study explores the relationship between coating thickness and Bp corroborating the theoretical predictions.

Senevirathne, I. H.↗
DOE OSTI2020

LDMS-GPU: Lightweight Distributed Metric Service (LDMS) for NVIDIA GPGPUs

GPUs are now a fundamental accelerator for many high-performance computing applications. They are viewed by many as a technology facilitator for the surge in fields like machine learning and Convolutional Neural Networks. To deliver the best performance on a GPU, we need to create monitoring tools to ensure that we optimize the code to get the most performance and efficiency out of a GPU. Since NVIDIA GPUs are currently the most commonly implemented in HPC applications and systems, NVIDIA tools are the solution for performance monitoring. The Light-Weight Distributed Metric System (LDMS) at Sandia is an infrastructure widely adopted for large-scale systems and application monitoring. Sandia has developed CPU application monitoring capability within LDMS. Therefore, we chose to develop a GPU monitoring capability within the same framework. In this report, we discuss the current limitations in the NVIDIA monitoring tools, how we overcame such limitations, and present an overview of the tool we built to monitor GPU performance in LDMS and its capabilities. Also, we discuss our current validation results. Most of the performance counter results are the same in both vendor tools and our tool when using LDMS to collect these results. Furthermore, our tool provides these statistics during the entire runtime of the tool as a time series and not just aggregate statistics at the end of the application run. This allows the user to see the progress of the behavior of the applications during their lifetime.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Improving Signal-to-Noise Ratio (SNR) for Readout Signals Using Adaptive Filters on Reconfigurable Controls Hardware

This study investigates the optimization of Signal-to-Noise Ratio (SNR) in superconducting quantum computing readout signals through adaptive filtering. Quantum computing technology has the potential to revolutionize various fields by delivering exponential speedup in solving certain computational problems. However, the technology's practical implementation is hindered by the difficulty of extracting clean, reliable signals during the readout phase, with various sources of noise presenting a significant barrier to clean signals. This noise, often present in readout profiles due to imperfect isolation, degrades the system's overall SNR, thus impeding the ability to extract the quantum state accurately. The research leverages the power of adaptive filtering to improve the SNR of quantum computing readout signals. Specifically, an adaptive filter is implemented in a PYNQ overlay on an FPGA, and eventually will be connected to a quantum computing system. The system models the n oise with a Least Mean Squares (LMS) adaptive filter, and then subtracts the estimated noise from the received signal to improve the SNR. A Direct Memory Access (DMA) channel is used to handle the signal processing, delivering efficient, high-speed data transfer between the PYNQ system and the hardware. The study explores the benefits of this adaptive filtering technique, potentially providing a significant contribution to practical and fast quantum computing.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Final Cooling with Thick Wedges for a Muon Collider

A muon collider is a particle accelerator that collides muons rather than protons or electrons. An important component of such a collider is cooling, in which the emittance of the muon beam is reduced to make it suitable for use in a collider. In the last stage of this process, called final 4D cooling, emittance in the transverse axes is reduced while emittance in the longitudinal axis is allowed to grow. In previous studies, Final 4-D cooling used absorbers within very high field solenoids to cool low-momentum muons. Simulations of the systems did not reach the desired cooling design goals. In this study, we develop and optimize a different conceptual design for the final 4D cooling channel, which is based on single thick wedges. We used G4Beamline to simulate the channel and Python to generate and analyze particle distributions. We characterized the effects of various parameters on the performance of this system, including the initial emittances and Twiss parameters, the length and angle of the wedges, the length of the drift channel, and the length, frequency, gradient, and phase of the RF cavities. We optimized the design parameters of the cooling channel and produced two conceptual designs (corresponding to two possible starting points for the input beam) which achieve transverse cooling in both x and y by a factor of ≈ 3.5. These channels achieve a lower transverse and longitudinal emittance than the best design previously published.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

FES 2023 Theory & Simulation Performance Target

The Report for the 2023 Theory & Simulation Performance Target (TSPT) is presented herein. The 1st, 2nd and 3rd Quarter Reports are given in Sec. I, Sec. II and Sec. III. These are unchanged from those previously submitted, except that in this final document the equations, figures and citations have been automatically renumbered. The 4th Quarter Report is given in Sec. IV. The following material presents verification calculations of several widely used three-dimensional (3D) magnetohydrodynamic (MHD) equilibrium codes in strongly shaped configurations; for both fixed and free-boundary equilibria; for the so-called vacuum case and for the case of high plasma pressure; and, for magnetic fields that have nested magnetic surfaces and for fields that have a magnetic island. The configurations studied includes the recently discovered “precise-QA” case [Phys. Rev. Lett., 128, 035001 (2002)] and other highly optimized configurations. Because of the strong (extreme) shaping, the latter of these pose challenges to the 3D MHD codes. A new method is introduced for constructing coordinates in arbitrarily shaped toroidal domains, and this can improve the robustness of 3D equilibrium calculations. Beyond the equilibrium calculation, this report demonstrates the capability of performing non-linear extended-MHD modelling of optimized stellarator equilibria. This capability allows the assessment not only of linear stability thresholds but also of the nonlinear evolution of unstable equilibria to new configurations that might be far from the original configuration, with significantly different transport properties and magnetic topology. They also enable the self-consistent determination of pressure profiles in non-integrable magnetic topology through the use of highly anisotropic thermal-transport equations. The calculations presented here represent the first demonstration of extended-MHD calculations in strongly-shaped, optimized quasi-symmetric configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗