Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geometry 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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A GPU-based Approach for Path Planning Optimization via Travel Length Reduction

Typically, before constructing an object with an additive manufacturing system, the 3D object must be sent through a process called slicing. Slicing converts a 3D object commonly in the form of an STL file into a set of layers by horizontally intersecting a plane with the object at various heights. At each height, called a layer, multiple 2D polygons can be generated. Each polygon represents a boundary for solid geometry and is called an island. Each island is then comprised of multiple path types in an attempt to optimally fill the polygon. To move between each island and each islands’ paths, travels are inserted. Travels are simply motion by the system to move from one area of construction to another. Travels do not contribute to the construction of the object, and so, are considered wasted motion. In large-scale additive manufacturing, objects can be quite large and the distance between islands can be large as well. As a result, these travels can waste a significant amount of time. Ideally, travels would be as short as possible, however, computing global minimal travel paths is computationally expensive. To combat this problem, researchers at Oak Ridge National Lab developed a GPU-based approach to travel insertion based on a unique factoradic representation. This representation was then utilized by the GPU to solve the Traveling Salesman Problem (TSP). This algorithm was able to compute global minimal travel paths quickly resulting in faster object construction. A general investigation was also carried out to determine when a GPU vs CPU implementation would be beneficial.

Borish, Michael↗
DOE OSTI2023

The multichannel i -propyl + O2 reaction system: A model of secondary alkyl radical oxidation

The i-propyl + O2 reaction mechanism has been investigated by definitive quantum chemical methods to establish this system as a benchmark for the combustion of secondary alkyl radicals. Focal point analyses extrapolating to the ab initio limit were performed based on explicit computations with electron correlation treatments through coupled cluster single, double, triple, and quadruple excitations and basis sets up to cc-pV5Z. The rigorous coupled cluster single, double, and triple excitations/cc-pVTZ level of theory was used to fully optimize all reaction species and transition states, thus, removing some substantial flaws in reference geometries existing in the literature. The vital i-propylperoxy radical (MIN1) and its concerted elimination transition state (TS1) were found 34.8 and 4.4 kcal mol−1 below the reactants, respectively. Two β-hydrogen transfer transition states (TS2, TS2′) lie above the reactants by (1.4, 2.5) kcal mol−1 and display large Born–Oppenheimer diagonal corrections indicative of nearby surface crossings. An α-hydrogen transfer transition state (TS5) is discovered 5.7 kcal mol−1 above the reactants that bifurcates into equivalent α-peroxy radical hanging wells (MIN3) prior to a highly exothermic dissociation into acetone + OH. The reverse TS5 → MIN1 intrinsic reaction path also displays fascinating features, including another bifurcation and a conical intersection of potential energy surfaces. An exhaustive conformational search of two hydroperoxypropyl (QOOH) intermediates (MIN2 and MIN3) of the i-propyl + O2 system located nine rotamers within 0.9 kcal mol−1 of the corresponding lowest-energy minima.

Chemistry↗
DOE OSTI2025

On measuring stimulated photon–photon scattering using multiple ultraintense lasers

Stimulated photon–photon scattering is a predicted consequence of quantum electrodynamics that has yet to be measured directly. Measuring the cross section for stimulated photon–photon scattering is the aim of a flagship experiment for NSF OPAL, a proposed laser user facility with two, 25-PW beamlines. We present optimized experimental designs for achieving this challenging and canonical measurement. A family of experimental geometries is identified that satisfies the momentum- and energy-matching conditions for two selected laser frequency options. Numerical models predict a maximum signal exceeding 1000 scattered photons per shot at the experimental conditions envisaged at NSF OPAL. Experimental requirements on collision geometry, polarization, cotiming and copointing, background suppression, and diagnostic technologies are investigated numerically. These results confirm that a beam cotiming shorter than the pulse duration and control of the copointing on a scale smaller than the shortest laser wavelength are needed to robustly scatter photons on a per-shot basis. Finally, we assess the bounds that a successful execution of this experiment may place on the mass scale of Born–Infeld nonlinear electrodynamics beyond the standard model of physics.

Beyond the Standard Model↗
DOE OSTI2023

Multi-Probe Campaign FY24 Omega EP Shot Request [Slides]

The ultimate goal of the Multi-Probe experimental campaign is to assess the feasibility of using multiple probe species (x-ray, proton, neutron, and electron) in future LANL radiographic facilities. To use Omega EP most efficiently, two different experiments are proposed for each of the two requested days in FY24. In this way, shots may be interleaved on the two short pulse beams and the shot rate doubled. In FY24, we propose to continue to develop the neutron source by first optimizing the “pitcher” source of deuterons/protons and then adding the “catcher” neutron converter on day 2. The second series of experiments continues examining the origin and nature of the x-ray background found in our last two experiments. This background is not well understood or characterized by the community. The first day will measure the electrons emanating from the proton and x-ray targets while also optimizing the background shielding first used in FY23. The second day will test a geometry that enables simultaneous radiographs using both short-pulse beams.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Critical problems of energy frontier Muon Colliders: optics, magnets and radiation

This White Paper brings together our previous studies on a Muon Collider (MC) and presents a design concept of the 6 TeV MC optics, the superconducting (SC) magnets, and a preliminary analysis of the protection system to reduce magnet radiation loads as well as particle backgrounds in the detector. The SC magnets and detector protection considerations impose strict limitations on the lattice choice, hence the design of the collider optics, magnets and Machine Detector Interface (MDI) are closely intertwined. As a first approximation we use the Interaction Region (IR) design with beta-star=3 mm, whereas for the arcs we re-scale the arc cell design of the 3 TeV MC. Traditional cos-theta coil geometry and Nb3Sn superconductor were used to provide field maps for the analysis and optimization of the arc lattice and IR design, as well as for studies of beam dynamics and magnet protection against radiation. The stress management in the coil will be needed to avoid large degradation or even damage of the brittle SC coils. In the assumed IR designs, the dipoles close to the Interaction Point (IP) and tungsten masks in each IR (to protect magnets) help reducing background particle fluxes in the detector by a substantial factor. The tungsten nozzles in the 6 to 600 cm region from the IP, assisted by the detector solenoid field, trap most of the decay electrons created close to the IP as well as most of the incoherent electron-positron pairs generated in the IP. With sophisticated tungsten, iron, concrete and borated polyethylene shielding in the MDI region, the total reduction of background loads by more than three orders of magnitude can be achieved.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Realization of a high luminosity muon beam at the Jefferson Lab

The proton is the major component of the visible matter of our universe, so being able to determine its fundamental properties, such as its charge radius, is of vital importance. In recent years, several experiments have been performed for the accurate measurement of the proton radius with different experimental techniques. The main and most common are based on two different techniques: ep scattering and Lamb shift spectroscopy of regular hydrogen atoms. The latter can also be performed on muonic hydrogen atoms with a significant improvement in accuracy. Recent measurements that use this method have highlighted a discrepancy between the results obtained with the diffusion of electrons and those obtained with the spectroscopic technique leading to the so-called ?proton radius puzzle?. The explanation of why this difference between the measured values exists is still a matter of discussion. The last attempt to address the puzzle was the PRad Experiment. It was performed with a unique experimental setup that gave the scientists an excellent control on the systematic errors that are usually present in electron-proton scattering experiments. This has significantly increased the precision of the measurement, bringing it in the same area as the spectroscopic ones done on muonic hydrogen atoms, which are considered the most precise. However, to solve the puzzle once and for all, even more precise electron scattering measurements are needed. Alternatively, new methods can be used, such as eg. muon-proton scattering at low transferred momenta in order to lower the systematic errors that comes with the measurement. This thesis tries to address this problem by studying the set up of an intense muon beam at the Jefferson Laboratory (JLab) that can be used to carry out a high-precision m ? p scattering experiment. Using Montecarlo simulations I studied the muon production resulting from the interaction of the primary electron beam with the beam-dump at JLab Hall-A. The simulations were performed using the FLUKA package via the FLAIR graphical interface. In this study, I modified an existing geometry (built for the ?Beam Dump eXperiment?) and performed several studies in order to optimize the muon beam parameters. In order to estimate the beam intensity, I studied the muon generation and transportation through the concrete bunker surrounding the Hall-A beam-dump. Detailed studies on muon attenuation and dispersion at different depth allowed me to track the beam profile in space and momentum as a function of the distance from the beam-dump. Last, but not least, I studied a possible tagging system for the muon beam. The tagging system will allow us to know with good precision (O(1=1000)) the energy of each muon of the beam. Based on the results obtained, I created an event generator that incorporates the detailed description of the muon beam (energy, position, emittance) useful to define and optimize the experimental setup for a future measurement of the proton radius at the JLab.

Fulci, Antonino↗
DOE OSTI2020

CONFIDANTE Demonstration Prototype Report

This report summarizes design optimization and performance evaluation studies for a new prototype CONFIDANTE (CONfirmation using a Fast-neutron Imaging Detector with Anti-image Null-positive Time Encoding) warhead confirmation system. It was found that a spherical mask geometry and a 2” diameter cylindrical central detector is expected to best discriminate between the parametrically varied source distributions that were evaluated. The optimized design as fabricated and its performance was evaluated in a series of laboratory measurements. The performance was in good agreement with the design studies, with demonstrated discrimination between objects with differences in scale on the order of 5 centimeters or better at 1 meter stand-off

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

A Block-Structured Adaptive Mesh Framework to Solve Radiation Transfer Equation in Irregular Embedded Geometries

Radiation transport arises in various scientific, industrial, and medical fields, and understanding its effect in applications is needed to make accurate predictions, safety assessments and performance optimizations. Solving the Radiation Transport Equation (RTE) is challenging due to its integro-differential nature, which involves both differential and integral terms. The differential term describes the change in radiation intensity due to absorption and emission, while the integral term accounts for scattering. The accurate modeling of radiation is further complicated in many applications due to the complex, irregular geometries. Various methods exist for solving the RTE, including the zonal, Monte Carlo, spherical harmonics, discrete ordinates, and finite volume methods. Traditional mesh-based approaches, which rely on structured or unstructured meshes, struggle with irregular geometries due to: a) the difficulty of conforming structured grids to irregular domains, b) challenges in enforcing boundary conditions correctly, and c) the additional computational cost of unstructured mesh methods. This work presents a second-order accurate method for solving the RTE in irregular geometries. The radiation intensity is discretized using the finite-volume method in both spatial and angular directions on regular Cartesian grid blocks. Leveraging the block-structured adaptive mesh refinement (AMR) framework provided by AMReX, our method refines the grid locally to reduce spatial discretization error, ensuring a converged numerical solution while minimizing computational costs elsewhere. A two-stage deferred correction approach is employed: First, a first-order discretization on grid blocks is solved using an algebraic multigrid method in HYPRE. Second, a correction term is applied explicitly to achieve second-order accuracy. The correction term is calculated by approximating the radiation flux on cell faces using a Total Variation Diminishing (TVD) scheme. This approach ensures quick convergence of the multigrid method while preserving higher-order accuracy of the numerical solution. Irregular geometries are resolved as embedded boundaries (EB), resulting in both cut cells and regular cells. In cut cells, we modify the fluxes using face fractions and incorporate additional contributions from EB boundary conditions. To ensure higher-order convergence near the EB interface, the correction term is modified by interpolating the radiation intensity to fictitious ghost points. The implementation takes advantage of modern supercomputers by leveraging AMReX’sMPI/X parallelization strategy where X can be MPI or a GPU accelerator including CUDA, HIP and DPC++. We validate our solver using classical test cases, both with and without EB, demonstrating accuracy and efficiency. Additionally, we analyze the impact of adaptive mesh refinement on solution accuracy and computational cost, highlighting the advantages of our approach for high-resolution radiation transport simulations.

computational fluid dynamics (CFD)↗
DOE OSTI2024

3D high-fidelity automated neutronics guided optimization of fusion blanket designs

The compact Fusion Pilot Plant (FPP) is defined in the recent National Academies of Sciences, Engineering, and Medicine report as the next step of fusion energy demonstration with a $50$ MWe peak net electricity production, $Q_e$ greater than $1$, and at least $3$ hours of continuous operation. This fusion pilot plant will be a test bed enabling materials, designs, and fuel management assessment, and it will represent an engineering challenge because of its high-fusion power and compact design targets. Previous reactor data is limited to experiments operating in different design space ranges. Therefore, design iterations and assessments should rely on high-fidelity first-principle theoretical and computational models. The high-fidelity integrated modeling of the plasma is a fundamental part of fusion energy research. However, the whole device modeling is often neglected, utilizing low-fidelity, system-level analysis. Recently, the need for high-fidelity multi-physics modeling was recognized, resulting in a selection of integrated tools. Further, autonomous design optimization requires a streamlined framework that perturbs the design point, reruns the analysis, and examines the outputs. However, high-fidelity analysis requires complex geometry specification that is difficult to perturb. This work presents the parametric CAD generation tool TRACER and a new neutronic workflow. TRACER allows the perturbation of the geometry representation, creating geometry files ready for further analysis. The streamlined neutronic workflow allows efficient and accurate calculations. The two new tools coupled together were used to perform a 3D high-fidelity multi-objective, multi-input optimization of an "ARC Class" compact tokamak design. The workflow was driven by an optimization driver for full automation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Proximal Galerkin: A Structure-Preserving Finite Element Method for Pointwise Bound Constraints

The proximal Galerkin finite element method is a high-order, low iteration complexity, nonlinear numerical method that preserves the geometric and algebraic structure of pointwise bound constraints in infinite-dimensional function spaces. This paper introduces the proximal Galerkin method and applies it to solve free boundary problems, enforce discrete maximum principles, and develop a scalable, mesh-independent algorithm for optimal design with pointwise bound constraints. This paper also introduces the latent variable proximal point (LVPP) algorithm, from which the proximal Galerkin method derives. When analyzing the classical obstacle problem, we discover that the underlying variational inequality can be replaced by a sequence of second-order partial differential equations (PDEs) that are readily discretized and solved with, e.g., the proximal Galerkin method. Throughout this work, we arrive at several contributions that may be of independent interest. These include (1) a semilinear PDE we refer to as the entropic Poisson equation; (2) an algebraic/geometric connection between high-order positivity-preserving discretizations and certain infinite-dimensional Lie groups; and (3) a gradient-based, bound-preserving algorithm for two-field, density-based topology optimization. The complete proximal Galerkin methodology combines ideas from nonlinear programming, functional analysis, tropical algebra, and differential geometry and can potentially lead to new synergies among these areas as well as within variational and numerical analysis. Open-source implementations of our methods accompany this work to facilitate reproduction and broader adoption.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

De novo design of pH-responsive self-assembling helical protein filaments

Abstract Biological evolution has led to precise and dynamic nanostructures that reconfigure in response to pH and other environmental conditions. However, designing micrometre-scale protein nanostructures that are environmentally responsive remains a challenge. Here we describe the de novo design of pH-responsive protein filaments built from subunits containing six or nine buried histidine residues that assemble into micrometre-scale, well-ordered fibres at neutral pH. The cryogenic electron microscopy structure of an optimized design is nearly identical to the computational design model for both the subunit internal geometry and the subunit packing into the fibre. Electron, fluorescent and atomic force microscopy characterization reveal a sharp and reversible transition from assembled to disassembled fibres over 0.3 pH units, and rapid fibre disassembly in less than 1 s following a drop in pH. The midpoint of the transition can be tuned by modulating buried histidine-containing hydrogen bond networks. Computational protein design thus provides a route to creating unbound nanomaterials that rapidly respond to small pH changes.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Role of initial conditions in plasma-current coupling of gas-puff Z-pinches

Azimuthal magnetic field measurements obtained during the implosion phase of an oxygen gas-puff Z-pinch on a 500 kA peak current and 180 ns rise time linear transformer driver are presented. While a fraction of the driver current was measured within the imploding plasma, key initial conditions were found to significantly impact the delivery of current to the plasma load. The electrode geometry was modified to assist the initial dielectric breakdown and resulted in improved shot reproducibility. Optimization of the gas injection plenum pressure and timing resulted in an increase in the current coupling parameter, defined as the ratio of the measured value of $B$ $θ$ to the expected value, from 50% to 75%. Furthermore, the degree of radial expansion of the gas puff in the load region, which is suspected to lead to the observed current loss during the implosion, was reduced by shortening the valve opening duration. Additionally, a pre-embedded axial magnetic field of up to 0.2 T was found to have no significant impact on the plasma-current coupling of the oxygen implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Multiphysics Running-In Simulations for Pebble-Bed Reactors with Griffin

Griffin, a Multiphysics Object-Oriented Simulation Environment (MOOSE)–based application targeting transient modeling of advanced reactors, has been used recently to model pebble-bed reactors (PBRs). The modeling effort has focused thus far on equilibrium core calculations. A new capability to simulate the running-in phase of PBR operation has been added to Griffin. This work demonstrates the new capability with a coupled multiphysics running-in simulation. Griffin computes power densities in the core at each time step of the running-in simulation and passes these to Pronghorn, which models fluid flow and heat transfer to calculate pebble surface temperatures. These surface temperatures are used along with the power densities in a heat conduction model to compute average fuel and moderator temperatures, which are passed back to Griffin and accounted for with temperature-dependent cross sections. This work also describes a novel methodology for determining appropriate pebble feed rates and control rod positioning during the running-in simulation. Furthermore, the RZ-geometry model used in this work requires minimal computational resources and can be used for optimization and uncertainty studies in future works.

Griffin↗
DOE OSTI2023

Monolithically Integrated $\varepsilon$-Ge/In x Ga 1-x As Quantum Well Laser Design: Experimental and Theoretical Investigation

Here, we have analyzed the electrical and optical phenomenon occurring in a $\varepsilon$-Ge/In x Ga 1-x -Ge/In x Ga 1-x As quantum well (QW) laser through self-consistent physical solvers calibrated using in-house experimental results. A separate confinement heterostructure QW design is proposed to enable lasing from tensile strained germanium ($\varepsilon$-Ge/In x Ga 1-x -Ge) in the range of 1.55 um to 4 um wavelengths as a function of QW thickness and indium (In) composition. Different recombination mechanisms were analyzed as a function of tensile strain in $\varepsilon$-Ge/In x Ga 1-x -Ge QW. Minority carrier lifetime and band alignment are key attributes of a QW laser, which were measured using microwave photoconductive decay and x-ray photoelectron spectroscopy (as a function of In composition), respectively. The transition point of Ge to a direct bandgap material is re-affirmed to be at $\varepsilon$-Ge/In x Ga 1-x = 1.6% (In ~24%) and the transition from type I to type II for $\varepsilon$-Ge/In x Ga 1-x -Ge/In x Ga 1-x As QW is found to be at In ~55%. Also, the transition to a TM mode dominant laser is identified at In ~15%. Using a tunable waveguide design to optimize confinement as a function of In composition, strain, wavelength, QW thickness, refractive index, and geometry, the $\varepsilon$-Ge/In x Ga 1-x -Ge QW laser design provided a net material gain of ~2000 cm -1 and a threshold current density of ~5 kA/cm 2 , which is an improvement over existing Ge based lasers. In conclusion, the impact of In composition and QW thickness on the band structure, polarized gain spectra, and various lasing metrics were analyzed to show $\varepsilon$-Ge/In x Ga 1-x -Ge/InGaAs QW lasers as promising for integrated photonics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Localization of fluorescent targets in deep tissue with expanded beam illumination for studies of cancer and the brain

Imaging fluorescence through millimeters or centimeters of tissue has important in vivo applications, such as guiding surgery and studying the brain. Often, the important information is the location of one of more optical reporters, rather than the specifics of the local geometry, motivating the need for a localization method that provides this information. We present an optimization approach based on a diffusion model for the fast localization of fluorescent inhomogeneities in deep tissue with expanded beam illumination that simplifies the experiment and the reconstruction. We show that the position of a fluorescent inhomogeneity can be estimated while assuming homogeneous tissue parameters and without having to model the excitation profile, reducing the computational burden and improving the utility of the method. We perform two experiments as a demonstration. First, a tumor in a mouse is localized using a near infrared folate-targeted fluorescent agent (OTL38). This result shows that localization can quickly provide tumor depth information, which could reduce damage to healthy tissue during fluorescence-guided surgery. Second, another near infrared fluorescent agent (ATTO647N) is injected into the brain of a rat, and localized through the intact skull and surface tissue. This result will enable studies of protein aggregation and neuron signaling.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

A Simulated Evaluation of Powder Flowability Through a Partially Obstructed Consumable in Blown Powder Directed Energy Deposition Systems

Abstract In the interest of continued industrialization of metal additive manufacturing in modern production environments, cost is often referenced as a primary deterrent to new adopters. Conventional economic models for additive systems, processes, and supply chains often focus on specific process applications with little generalizability, or they neglect significant costs associated with production such as machine maintenance and consumable part replacement. Compounding the latter issue are substantial knowledge gaps in consumable part wear characterization for additive and other convergent manufacturing systems. In coaxial blown powder directed energy deposition systems, gas atomized metal powder is wasted during material deposition at a rate that is partly dependent on present wear phenomena in a consumable nozzle housed in the cladding head assembly. The price and lead time required to replace the nozzle incentivizes its reuse even when visibly worn. Often this initiates a process quality decline in the form of underbuilt geometry and internal defects due to losses in powder catchment efficiency. While depositing H13 steel using a hybrid manufacturing machine tool equipped with such a deposition system, a unique partial clog with a bridge-like structure formed at the consumable nozzle exit when supporting argon gas flows failed mid-process. To further understand coaxial multi-phase powder flow in the event of support gas failure, a computational fluid dynamics simulation is tailored to relevant process parameters, H13 powder material profile, and machine operator observations collected after the incident. The resulting differences in powder flow compared to control gas flow parameters is presented and discussed. The powder flowability and performance of the clogged nozzle is then assessed by using an optical profilometer to extract the profile of the clog and recreate the clog geometry within the simulation environment. In past work this simulation has been experimentally validated for a 316L steel powder material profile and used specifically for analyzing powder stream geometry and catchment efficiency. After the initial powder flow characterization, the clog is removed, and the nozzle is reprofiled. After removing the obstructing clog, the newly unobstructed nozzle geometry, the original off the shelf nozzle geometry, and additional nozzle profiles exploring different consumable refurbishment strategies are reevaluated in the simulation. Powder catchment efficiency for all variant nozzle geometries and relevant flow variables are compared and discussed, along with potential mitigation strategies for optimizing powder flowability with worn consumables. This work expands on the known morphology of blown powder obstructions and wear defects present in consumable coaxial nozzles while discussing pragmatic simulation driven responses to unanticipated subsystem failure in hybrid manufacturing machining platforms.

DeWitte, Lisa↗
DOE OSTI2023

Magnetic design of the commutational magnet and quadrupoles for PERLE accelerator

PERLE (Powerful Energy Recovery LINAC for Experiment) is a high-power Energy Recovery LINAC (ERL) facility with 20 mA beam current and beam energy from 250 MeV to 500 MeV featuring three passes through two cryomodules. It is a hub for validation of the ERL technology development towards future energy and intensity frontier machines. Design challenges of PERLE and its beam parameters make it a test-bed to validate multi-turn high current ERL operation for the LHeC. It will be the first ERL for some pioneering experiment of the eN interaction with radioactive nuclei. In this work, design and optimization of the commutational magnet (B-com) used to spread/combine the three beams and one series of the quadrupole magnet is discussed. It gives the design parameters including: yoke geometry, pole profile, and material, and calculation of the excitation current needed to drive the magnet, the coil parameters and the number of turns. The B-com magnet is optimized for a 30° bending angle with magnetic field of 0.88 T along the magnet length and a harmonic content of 0.036%. The quadrupole magnet is designed for a gradient field of 34.15 T/m and experiences saturation above this value. Further studies to avoid saturation and achieve the maximum gradient of 44.1 T/m required by the beam dynamics is undergoing.

Abukeshek, R.↗
DOE OSTI2024

ParaStell: parametric modeling and neutronics support for stellarator fusion power plants

The three-dimensional variation inherent to stellarator geometries and fusion sources motivates three-dimensional modeling to obtain accurate results from computational modeling in support of design and analysis of first wall, blanket, and shield (FWBS) systems. Manually constructing stellarator fusion power plant geometries in computer-aided design (CAD) and defining the corresponding fusion source can be cumbersome and challenging. The open-source parametric modeling toolset ParaStell has been developed to automate construction of such geometries in low-fidelity. Low-fidelity modeling is useful during the conceptual phase of engineering design as a means of rapidly exploring the design space of a given device. The modeling capability of ParaStell includes in-vessel components and magnets, for any given stellarator configuration, using a parametric definition and plasma equilibrium data. Furthermore, the toolset automates the generation of detailed, tetrahedral neutron source definitions and DAGMC geometries for use in neutronics modeling. ParaStell assists rapid design iteration, parametric study, and design optimization of stellarator fusion cores. As a demonstration of the design iteration capability, the effect of the three-dimensional parameter space on tritium breeding and magnet shielding is investigated, using the WISTELL-D configuration as a design basis. Blanket and shield thicknesses are varied in three dimensions, using the space available between the plasma edge and magnet coils as a constraint. The corresponding effects on tritium breeding ratio and magnet heating are tallied using the open-source Monte Carlo particle transport code OpenMC. The inclusion of additional and higher-fidelity modeling capabilities is planned for ParaStell’s future, as well as its implementation in machine-driven optimization.

Moreno, Connor↗