Search NASA⌕ Search

SEARCH · Search NASA

Results for “procedures”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Report on the AAPM grand challenge on deep generative modeling for learning medical image statistics

Abstract Background The findings of the 2023 AAPM Grand Challenge on Deep Generative Modeling for Learning Medical Image Statistics are reported in this Special Report. Purpose The goal of this challenge was to promote the development of deep generative models for medical imaging and to emphasize the need for their domain‐relevant assessments via the analysis of relevant image statistics. Methods As part of this Grand Challenge, a common training dataset and an evaluation procedure was developed for benchmarking deep generative models for medical image synthesis. To create the training dataset, an established 3D virtual breast phantom was adapted. The resulting dataset comprised about 108 000 images of size 512 512. For the evaluation of submissions to the Challenge, an ensemble of 10 000 DGM‐generated images from each submission was employed. The evaluation procedure consisted of two stages. In the first stage, a preliminary check for memorization and image quality (via the Fréchet Inception Distance [FID]) was performed. Submissions that passed the first stage were then evaluated for the reproducibility of image statistics corresponding to several feature families including texture, morphology, image moments, fractal statistics, and skeleton statistics. A summary measure in this feature space was employed to rank the submissions. Additional analyses of submissions was performed to assess DGM performance specific to individual feature families, the four classes in the training data, and also to identify various artifacts. Results Fifty‐eight submissions from 12 unique users were received for this Challenge. Out of these 12 submissions, 9 submissions passed the first stage of evaluation and were eligible for ranking. The top‐ranked submission employed a conditional latent diffusion model, whereas the joint runners‐up employed a generative adversarial network, followed by another network for image superresolution. In general, we observed that the overall ranking of the top 9 submissions according to our evaluation method (i) did not match the FID‐based ranking, and (ii) differed with respect to individual feature families. Another important finding from our additional analyses was that different DGMs demonstrated similar kinds of artifacts. Conclusions This Grand Challenge highlighted the need for domain‐specific evaluation to further DGM design as well as deployment. It also demonstrated that the specification of a DGM may differ depending on its intended use.

Radiology, Nuclear Medicine & Medical Imaging↗
DOE OSTI2024

A Corrected Score Function Framework for Modelling Circadian Gene Expression

Many biological processes display oscillatory behaviour based on an approximately 24 h internal timing system specific to each individual. One process of particular interest is gene expression, for which several circadian transcriptomic studies have identified associations between gene expression during a 24 h period and an individual's health. A challenge with analysing data from these studies is that each individual's internal timing system is offset relative to the 24 h day-night cycle, where day–night cycle time is recorded for each collected sample. Laboratory procedures can accurately determine each individual's offset and determine the internal time of sample collection. However, these laboratory procedures are labour-intensive and expensive. Here, in this paper, we propose a corrected score function framework to obtain a regression model of gene expression given internal time when the offset of each individual is too burdensome to determine. A feature of this framework is that it does not require the probability distribution generating offsets to be symmetric with a mean of zero. Simulation studies validate the use of this corrected score function framework for cosinor regression, which is prevalent in circadian transcriptomic studies. Illustrations with data from three circadian transcriptomic studies further demonstrate that the proposed framework consistently mitigates bias relative to using a score function that does not account for this offset.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Thermal partition function of $$ {J}_3{\overline{J}}_3 $$ deformed AdS3

Abstract We derive a compact formula for the one-loop, bosonic string partition function of Euclideanized$$ {J}_3{\overline{J}}_3 $$ J 3 J ¯ 3 deformedAdS 3 with periodic Euclidean time as an integral transform of the partition function of the undeformed EuclideanizedAdS 3 . Such a deformation is interpretable as an irrelevant “single-trace$$ T\overline{T} $$ T T ¯ deformation” of the boundary. We will do this by first establishing a formal procedure to compute a worldsheet torus zero point function for an exactly marginal$$ J\overline{J} $$ J J ¯ deformation of a sigma model with U(1) L × U(1) R global symmetry. We then describe how this procedure is implemented on SL(2,R) sigma model and its Euclidean continuation. Finally, we describe the embedding of the deformed SL(2,R) torus amplitude into critical string theory and interpret the result as the leading perturbative contribution to the thermal partition function of the deformed theory.

Physics↗
DOE OSTI2024

Quantum entanglement and Bell inequality violation in semi-leptonic top decays

Abstract Quantum entanglement is a fundamental property of quantum mechanics. Recently, studies have explored entanglement in the$$ t\overline{t} $$ t t ¯ system at the Large Hadron Collider (LHC) when both the top quark and anti-top quark decay leptonically. Entanglement is detected via correlations between the polarizations of the top and anti-top and these polarizations are measured through the angles of the decay products of the top and anti-top. In this work, we propose searching for evidence of quantum entanglement in the semi-leptonic decay channel where the final state includes one lepton, one neutrino, twob-flavor tagged jets, and two light jets from theWdecay. We find that this channel is both easier to reconstruct and has a larger effective quantity of data than the fully leptonic channel. As a result, the semi-leptonic channel is 60% more sensitive to quantum entanglement and a factor of 3 more sensitive to Bell inequality violation, compared to the leptonic channel. In 139 fb −1 (3 ab −1 ) of data at the LHC (HL-LHC), it should be feasible to measure entanglement at a precision of ≲ 3% (0.7%). Detecting Bell inequality violation, on the other hand, is more challenging. With 300 fb −1 (3 ab −1 ) of integrated luminosity at the LHC Run-3 (HL-LHC), we expect a sensitivity of 1.3σ(4.1σ). In our study, we utilize a realistic parametric fitting procedure to optimally recover the true angular distributions from detector effects. Compared to unfolding this procedure yields more stable results.

Physics↗
DOE OSTI2025

Orientation reversal and the Chern-Simons natural boundary

We show that the fundamental property of preservation of relations, underlying resurgent analysis, provides a new perspective on crossing a natural boundary, an important general problem in theoretical and mathematical physics. This reveals a deeper rigidity aspect of resurgence in a quantum field theory path integral. The physical context here is the non-perturbative completion of complex Chern-Simons theory that associates to a 3-manifold a collection of q-series invariants labeled by Spinc structures, for which crossing the natural boundary corresponds to orientation reversal of the 3-manifold. Our new resurgent perspective leads to a practical numerical algorithm that generates q-series which are dual to unary q-series composed of false theta functions. Until recently, these duals were only known in a limited number of cases, essentially based on Ramanujan’s mock theta functions, and the common belief was that the duals might not even exist in the general case. Resurgence analysis identifies as primary objects Mordell integrals: up to changes of variables, they are Laplace transforms of resurgent functions. Their unique Borel summed transseries decomposition on either side of the Stokes line is simply the unique decomposition into real and imaginary parts. In turn, the latter are combinations of unary q-series in terms of q and its modular counterpart $\overset{\sim }{q}$ , and are resurgent by construction. The Mordell integral is analytic across the natural boundary of the q and $\overset{\sim }{q}$ series, and uniqueness of a similar decomposition which preserves algebraic relations on the other side of the boundary defines the unique boundary crossing of the q series. We demonstrate that this continuation can be efficiently implemented numerically. In the cases where unique mock modular identities are known, they are found by this numerical procedure, but the procedure can go well beyond the known list of identities. A particularly interesting feature of the resurgent approach is that it reveals new aspects, and is very different from other known approaches based on indefinite theta series, Appell-Lerch sums, and representation theory of logarithmic vertex operator algebras.

Chern-Simons theories↗
DOE OSTI2025

10-th order of accuracy for numerical solution of 3-D elasticity equations for heterogeneous materials on unfitted Cartesian meshes

We have developed the Optimal Local Truncation Error Method (OLTEM) with 10-th order of accuracy on unfitted Cartesian meshes for a system of 3-D elasticity equations with smooth irregular interfaces. 5 x 5 x 5 = 125-point stencils (similar to those for quadratic finite elements) for elastic heterogeneous materials are used for OLTEM. There are no unknowns at the interface points between different materials; the structure of the global discrete equations is the same for homogeneous and heterogeneous materials. The calculation of unknown stencil coefficients is based on the minimization of the local truncation error of the stencil equations and yields the optimal 10-th order of accuracy for OLTEM on unfitted Cartesian meshes, i.e., the increase by 7 orders in accuracy compared to quadratic finite elements on conformal meshes. A new post-processing procedure provides the 9-th order of accuracy for stresses in the 3-D case. Similar to basic computations it uses OLTEM with the 125-point stencils, the interface conditions and the elasticity equations. It was shown that the use of the elasticity equations for post-processing improves the accuracy of 0.1% stresses by 6 orders compared to post-processing without the use of PDEs. At an accuracy of for stresses, OLTEM with the new post-processing procedure reduces the number of degrees of freedom by 360 - 8000 times compared to quadratic finite elements with similar stencils. OLTEM with the 125-point stencils yields even more accurate results than high-order finite elements with much wider stencils. OLTEM provides accurate numerical results for compressible and nearly incompressible materials.

elasticity equations↗
DOE OSTI2024

Simple separation of technetium from molybdenum for tracer isotope production

A procedure for the separation of technetium isotopes from bulk molybdenum was developed in nitric acid media. After irradiation of nat Mo foils, and dissolution in H 2 O 2 , the solution is acidified to 2 M HNO 3 , and anion exchange resin is used to separate technetium from molybdenum. The procedure is simple, requiring only basic laboratory equipment and, with a two-column separation, the technetium yield is high (~ 85%) with extremely high purity (< 0.1 ppm nat Mo). In conclusion, this is ideally suited for laboratory production of technetium tracer isotopes, particularly 95m Tc.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Capacitated p -hub approach for park-and-ride facility location problem under nested logit demand function: polyhedral approaches

By generalizing the unconstrained p-hub approach for the park-and-ride (P&R) facility location problem under the multinomial logit demand function, the capacitated p-hub approach for the problem under the nested logit demand function captures a broader range of real-world cases. To solve this problem optimally, we introduce a mixed-integer linear program and accelerate its solution by enhancing the branch-and-cut procedure. To address the problem at a large scale, we introduce two other polyhedral approaches: variable neighborhood search (VNS) and adaptive randomized rounding (ARR). Downtown areas in Seoul have a high modal share of public transportation and congested road traffic, yet P&R has not been widely implemented. Therefore, we apply the ARR procedure to solve a real-world problem using traffic and geographic data from the Seoul metropolitan area. ARR performs better than VNS and addresses real-world cases. The solutions obtained by ARR present a phased expansion plan that encourages policymakers to start installing a small number of P&Rs immediately.

Capacitated p-hub approach↗
DOE OSTI2024

Simultaneous Optimization of Crystal Plasticity Hardening Parameters

Crystal plasticity models relate macroscopic deformation behavior to the evolution of slip systems strength, but their parameterization is often non-unique, with multiple parameter sets being able to describe the same macroscopic behavior. To address this issue, the present work adopts a Bayesian optimization framework for the parameterization of face-centered cubic plasticity models while simultaneously considering multiple experimental datasets from the literature. For single crystal Cu, parameter optimization was guided by the tensile stress–strain curves along several crystallographic orientations, with an adequate fit being found for five orientations at once. While additional parameters allowed for the consideration of more physical mechanisms, like different slip system interaction strengths or misorientations inherent to the experimental data, the extra dimensionality was found to limit the efficiency of the global minimization procedure. For polycrystalline Ni, multiple grain sizes were considered together in a representative polycrystalline model, with the optimization able to reconcile the model with the data for three grain sizes at once. As meaningful interpretation of parameters relies on the uniqueness of their values, incorporating multiple datasets into this discerning parameterization procedure enables more robust prediction and application of crystal plasticity models.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Reducing systematic bias in machine learning applications to J/ψ signal extraction in high-energy nuclear physics

Machine learning techniques are increasingly used in high-energy nuclear physics because they can exploit multivariate correlations more efficiently than conventional cut-based analyses. A central challenge is the construction of training samples that faithfully reproduce the detector response observed in data. Signal samples are usually derived from detector simulations; therefore, mismatches between simulation and data can degrade classifier performance and introduce systematic biases. This work presents two practical correction procedures, namely cumulative distribution function (CDF) mapping and a shift-and-scale transformation, to align simulated signal features with those measured in data. Their performance is demonstrated with $J$/$\psi$ yield measurements in $\sqrt{s_{nn}}$ = 200 GeV Ru+Ru and Zr+Zr collisions recorded by STAR. A set of self-consistency tests shows that these procedures substantially suppress the systematic bias associated with data-simulation discrepancies in machine-learning-based signal extraction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

An iterative dynamic chemical stiffness removal method for reacting flow simulations

Abstract An iterative dynamic chemical stiffness removal method (IDCSR) based on quasi-steady-state approximation (QSSA) is proposed. The IDCSR method is built on a previously developed non-iterative method which has proved to work well for small timestep sizes. A novel iterative procedure is designed in IDCSR to enable explicit time integration of stiff chemistry at relatively large timestep sizes relevant to practical reacting flow simulations. The effectiveness of the iterative procedure is first demonstrated with a toy problem and homogeneous auto-ignition with fixed integration step sizes, showing that larger timestep sizes can be allowed for explicit time integration using IDCSR compared with the previous non-iterative method. IDCSR is then compared with existing explicit chemistry solvers for simulations of homogeneous auto-ignition and shows similar or lower computational cost but significantly higher accuracy across a wide range of timestep sizes. IDCSR is further combined with an automatic adaptive time-stepping scheme for simulations of 0-D homogeneous auto-ignition and a 2-D laminar lifted n -dodecane jet flame. For the 0-D auto-ignition simulations, IDCSR is shown to reduce both the error (by 43%–90%) and computational cost (by 6–15 times) compared with existing explicit solvers, while achieving speed-up factors of up to 400 compared with VODE for a wide range of timestep sizes and reaction mechanisms. For the 2-D jet flame simulations, speed-up factors of 15 and 31 for chemistry integration, and 5 and 9 for overall simulation, are achieved by IDCSR compared with CVODE with and without analytic Jacobian, respectively.

Xu, Chao (ORCID:0000000153074159)↗
DOE OSTI2025

Role of Mg:Al ratio and Pt promotion on mitigating the deactivation of Ni-based methane steam reforming catalysts

Ni-based catalysts for CH 4 steam reforming require fine-tuning to withstand deactivation under dynamic operation conditions relevant to emerging H 2 -driven technologies. This study investigates the impact of the Mg:Al ratio and Pt presence in catalyst composition on the activity and stability of industrially relevant Ni-based mono- and bimetallic catalysts during simulated daily start-up and shut-down cycles in various gas atmospheres. The evolution of the catalyst structure during extensive testing procedures was investigated in detail by complementary electron microscopy, in situ/operando XAS and XRD. The results obtained revealed that catalyst deactivation is promoted at high Mg:Al ratios and especially affects the monometallic catalysts. By converting CH 4 at lower temperatures, Pt regulates the extent of Ni oxidation and its incorporation into MgO lattice. Further prevention of catalyst deactivation was achieved by optimizing the reactor shut-down procedure to minimize the simultaneous exposure to high temperatures and H 2 O vapors. By flushing the reactor with N 2 only around the reaction extinction temperature, a fraction of Ni species is maintained in metallic state, which is beneficial for the long-term activity and reaction operation economy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A new perspective on density and strength loss profiles at the surface of thermally oxidized nuclear graphite

Oxidation of graphite components could influence their designed life in a high-temperature nuclear reactor. The oxidized regions could potentially lower the allowed stress capacity. The American Society of Mechanical Engineers rules for the design and construction of graphite-moderated reactors recommend that subsurface regions that might become excessively damaged by oxidation during reactor operation be identified and excluded from geometry and stress calculations. Identification of oxidation-affected regions is possible, in principle, through complex modeling exercises of reactor behavior during hypothetical accident scenarios coupled with graphite oxidation models, but this procedure may not have the precision needed for informed decisions. Here, this paper proposes an alternate method, based on interpretation of a series of well-designed oxidation experiments, which could augment the designer's tools. The procedure is illustrated by data on oxidation by air of several graphite grades (NBG-18, PCEA, IG-110, R4-650) that are corroborated with independent literature information, when available. The Wichner model for graphite oxidation used for this analysis provides conservative results that could be quickly implemented in the design process.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Large eddy simulation of lean blow-off in swirl-stabilized flame with the subgrid dissipation concept

The lean blow-off mechanism of the premixed swirl flame is numerically investigated by large eddy simulation (LES) with the subgrid dissipation concept (SDC) combustion model. Three simulated cases cover stable, near blow-off, and transient conditions. Compared with the experiment, the LES-SDC approach captures the flow and combustion features for stable and near blow-off conditions. More importantly, the predictions of the blowoff procedure and duration agree satisfactorily with the experiment, indicating that the LES-SDC approach is a promising tool for predicting strong, unsteady turbulent combustion processes. Further, the numerical results are used to investigate the blow-off mechanism. Two stages in the blow-off procedure are specified. The first is the necking and extinction of the downstream flame surface, and the second is the shrinking of the upstream flame surface. The blow-off mechanism is well explained by the theory of stretched flame extinction. At the end of the recirculation zone, the large negative radial velocity pushes the flame to the central line. The combustion process here can be abstracted as the stretched counter-flame of the reactant-reactant configuration. The excessive flame stretch dominates the flame extinction and triggers the blow-off event. The upstream flame resists the intense stretch with the help of hot product recirculation, and the combustion here can be idealized as the counter-flame of the reactant-product configuration. The alignment of the temperature gradient and flow velocity, together with the excessive stretch, clearly indicates the tendency of flame local extinction. A Damköhler number-based blow-off criterion is raised from the mean flow strain rate and laminar flame bulk extinction strain rate.

42 ENGINEERING↗
DOE OSTI2024

Interpretable, extensible linear and symbolic regression models for charge density prediction using a hierarchy of many-body correlation descriptors

Here, density functional theory (DFT) is routinely used to make electronic structure predictions for high-throughput screening of materials and molecules for technologically relevant areas, like the identification of better catalysts, electronic materials, and drug discovery. However, the DFT formalism is limited by (a) its poor (quadratic-to-quartic) scaling, and (b) the need to perform repeated eigenvalue computations of the electronic Hamiltonian as part of its self-consistent field (SCF) iteration procedure to obtain the converged ground state electron density, ρ (r). Approaches that directly predict ρ (r) of a structure with high accuracy can accelerate conventional SCF calculations and can also be used in linearly scaling methods such as orbital-free DFT. To this end, we present a procedure to predict the ground state electron density of molecular and periodic three-dimensional systems directly from the atomic structure with a particular emphasis on physical interpretability. In our framework, ρ (r) is modeled using many-body correlation descriptors that accurately capture the effects of local atomic arrangements in the neighborhood of a grid point. Our use of a linear regression scheme to fit to charge density data enables transparent analysis of the relative contributions of various types of local atomic correlations. By systematically including increasingly complex correlations, our model is shown to accurately predict ρ (r) for a variety of chemically and electronically diverse systems — amorphous Ge, Al(001) slab, crystalline Ga 2 O 3 , molecular benzene, and polyethylene. We then demonstrate a symbolic regression-based protocol to construct easily computable, interpretable features from lower-order correlations that significantly improves our electron density predictions with effectively no increase in the computational cost.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Poromechanical cohesive interface element with combined Mode I-II cohesive zone elastoplasticity for simulating fracture in fluid-saturated porous media

A combined Mode I-II cohesive zone (CZ) elasto-plastic constitutive model, and a two-dimensional (2D) cohesive interface element (CIE) are formulated and implemented at small strain within an ABAQUS User Element (UEL) for simulating 2D crack nucleation and propagation in fluid-saturated porous media. Here, the CZ model mitigates problems of convergence for the global Newton-Raphson solver within ABAQUS, which when combined with a viscous stabilization procedure allows for simulation of post-peak response under load control for coupled poromechanical finite element analysis, such as concrete gravity dam stability analysis. Verification examples are presented, along with a more complex ambient limestone-concrete wedge fracture experiment, water-pressurized concrete wedge experiment, and concrete gravity dam stability analyses. A calibration procedure for estimating the CZ parameters is demonstrated with the limestone-concrete wedge fracture process. For the water-pressurized concrete wedge fracture experiment it is shown that the inherent time-dependence of the poromechanical CIE analysis provides a good match with experimental force versus displacement results at various crack mouth opening rates, yet misses the pore water pressure evolution ahead of the crack tip propagation. This is likely a result of the concrete being partially-saturated in the experiment, whereas the finite element analysis assumes fully water saturated concrete. For the concrete gravity dam analysis, it is shown that base crack opening and associated water uplift pressure leads to a reduced Factor of Safety, which is confirmed by separate analytical calculations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Component-based SHGC determination of BIPV glazing for product comparison

Building-integrated photovoltaic (BIPV) systems are intrinsically designed to generate electricity and to provide at least one building-related function. When BIPV modules act as glazing products in windows, skylights or curtain walls, their ability to control the transmission of solar energy into the building must be characterised by a Solar Heat Gain Coefficient (SHGC) or g value (also known as Total Solar Energy Transmittance – TSET – or “solar factor”). For the comparison of BIPV glazing products consisting of one PV laminate and possibly further, conventional glazing layers separated by gas-filled cavities, the procedures documented in international standards for architectural glazing (e.g. ISO 9050 and EN 410) form a suitable starting point. Easily implemented modifications to these procedures are proposed to take both optical inhomogeneity (if relevant) and extraction of electricity from BIPV glazing units into account. Geometrically complex glazing and shading devices, and light-scattering glazing layers, are outside the scope of the proposed methodology; SHGC determination for obliquely incident solar radiation is also excluded. For these cases, the experimental calorimetric approach documented in [ISO 19467:2017; ISO 19467-2:2021] is recommended. The paper also presents results and conclusions from an implementation exercise and sensitivity study carried out by participants of the IEA-PVPS Task 15 on BIPV. The cell coverage ratio in the PV laminate, the thermal resistance offered by the glazing configuration, the choice of boundary conditions and the effect of extracting electricity were all identified as parameters which significantly affect the SHGC value determined for a given type of BIPV glazing. A practicable approach to accommodate the great variety of dimensions typical for BIPV glazing is also proposed. These findings should pave the way for modifying the existing component-based standards for architectural glazing to take the specific features of BIPV glazing into account.

14 SOLAR ENERGY↗
DOE OSTI2024

Models implemented in the methodological approach to design the initial STEP first wall contour

The official Spherical Tokamak for Energy Production mission aims to demonstrate the ability to generate net electricity from fusion with the STEP Prototype Power plant. One of the key technological and engineering challenges in fusion power plants is managing the loads on the first wall within acceptable limits. Therefore, the conceptual design development of the STEP Prototype Power plant needs to be based on load estimates derived using legitimate plasma physics assumptions through dynamic and flexible tools. The current design foresees the STEP main chamber first wall to withstand steady-state heat loads of up to ~1 MW/m 2 , excluding critical regions expected to receive higher heat loads such as the baffle regions approaching the divertors. These critical areas will require ad hoc assessments and will be designed with the presence of limiters. This article focuses on the models and methodology adopted for designing the 2-D poloidal contour of the STEP first wall, based on the anticipated charged particle and radiation heat loads during normal operation. Firstly, the models adopted for calculating the charged particle and radiation heat loads are introduced. The first model is validated through benchmarking against the particle tracing code SMARDDA, while the second model is verified by comparing it with data from the MAST-U experiment. Secondly, the model used to design the 2-D first wall contour according to the heat loads is explained. We acknowledge that this preliminary design stage assumes certain simplifications, notably an axisymmetric geometry, for computational efficiency and clarity in presentation. It is understood that subsequent design phases will address the complexities of real-world engineering, including non-axisymmetric effects, transient plasma scenarios, and the impact of disruptions on the first wall design. Finally, an automatic procedure based on these models is presented for defining the 2-D poloidal contour of the STEP first wall to minimize heat loads, taking into account the need to radiate most of the alpha-particle and auxiliary heating power. Here, by providing an overview of the models, methodology, and an automatic procedure, this paper contributes to the design process of the STEP first wall, addressing the engineering challenges associated with fusion power plant development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗