Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “source code”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Development of mathematical model and simulation tool for the high capacity production of carbon fiber

Weight reduction in the transportation is the most effective way for reducing fuel consumption. Carbon fiber composites have been increasingly used by automotive manufacturers in vehicle nonstructural and structural parts. Larger penetration of carbon fiber-based materials requires very large quantities of fiber at a low price, and of consistently high quality. The advancements in carbon fiber technology that are required to achieve this production scale-up and cost reduction, are very challenging. The most expensive production step is the conversion of textile-grade polyacrylonitrile (PAN) precursor to carbon fiber. The conversion takes place in three main stages: stabilization, low and high temperature carbonization. Most of the research to date has focused on the stabilization of the PAN precursors and oxidized PAN fibers. Comparingly little effort has been dedicated to the investigation of processing conditions during low temperature (LT) carbonization. During this stage of fiber conversion, the largest fiber mass loss occurs, and the fundamental carbon fiber structure and material properties are established. The understanding of the conversion processes in the LT furnace, and the effects of spatial and temporal variations of the LT process parameters on the carbon fiber properties are necessary to identify possible process improvements and cost savings. The objective of this research was to develop computational models for gas flow and energy transport in the gas and fiber phases of the LT furnace, and to couple them with the fiber conversion process models. The model is implemented in the three-dimensional framework that enables analysis of process condition variations in all directions. The formulation is based on modified transport equation to simulate steady state conditions and thus significantly reduce computational effort associated with transient methods. The model is implemented in open source computational fluid dynamics (CFD) code OpenFOAM. The model has been demonstrated on a prototype carbon fiber LT furnace and can be extended to other geometries and process conditions. Extensions to other physical effects, such as accounting evolution of multiple gas species, and fiber tensioning are straightforward. The model uses High Performance Computing (HPC) capabilities of OpenFOAM and can be implemented on multi-processor workstations or mid-size computer clusters common in the industrial setting.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Solar Performance Insight (Final Report)

The PV Operations and Maintenance (O&M) service industry lacks an affordable, well-documented, intuitive PV modeling and analytics tool to calculate modeled performance from actual data from multiple data acquisition systems (DAS). We envision a performance modeling and analytics platform built on open-source, extensible, community-maintained code. The key innovation is the community-driven development of pvlib python delivered through a lightweight web service to provide configurable, consistent and reproducible PV modeling for O&M providers.

14 SOLAR ENERGY↗
DOE OSTI2026

Modeling The Nucleosynthetic Imprint of Stellar Merger Phenomena - Final Technical Report

This DOE CAREER project developed a comprehensive, multi-physics framework for modeling stellar mergers and their observable consequences. The work successfully integrated analytical models, three-dimensional hydrodynamic simulations, stellar evolution calculations, nucleosynthesis, and radiation transport into a unified pipeline. Key scientific advances include demonstrating stellar mergers as a unifying explanation for systems such as Betelgeuse and R Coronae Borealis stars, quantifying merger-driven chemical signatures, and establishing the connection between merger physics and circumstellar environments. A major outcome of the project is the development of SuperLite, an open-source Monte Carlo radiation transport code that enables the generation of synthetic spectra for astrophysical transients. This work bridges dynamical, thermal, and radiative timescales and provides a powerful, DOE-relevant computational framework for interpreting observations of supernovae and related transient phenomena, while also contributing to workforce development through the training of postdoctoral researchers and graduate students.

Chatzopoulos, Emmanouil [Louisiana State Universit↗
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↗
DOE OSTI2022

Monitoring Offshore CO 2 Sequestration Using Marine CSEM Methods; Constraints Inferred from Field- and Laboratory-Based Gas Hydrate Studies

Offshore geological sequestration of CO 2 offers a viable approach for reducing greenhouse gas emissions into the atmosphere. Strategies include injection of CO 2 into the deep-ocean or ocean-floor sediments, whereby depending on pressure–temperature conditions, CO 2 can be trapped physically, gravitationally, or converted to CO 2 hydrate. Energy-driven research continues to also advance CO 2 -for-CH 4 replacement strategies in the gas hydrate stability zone (GHSZ), producing methane for natural gas needs while sequestering CO 2 . In all cases, safe storage of CO 2 requires reliable monitoring of the targeted CO 2 injection sites and the integrity of the repository over time, including possible leakage. Electromagnetic technologies used for oil and gas exploration, sensitive to electrical conductivity, have long been considered an optimal monitoring method, as CO 2 , similar to hydrocarbons, typically exhibits lower conductivity than the surrounding medium. We apply 3D controlled-source electromagnetic (CSEM) forward modeling code to simulate an evolving CO 2 reservoir in deep-ocean sediments, demonstrating sufficient sensitivity and resolution of CSEM data to detect reservoir changes even before sophisticated inversion of data. Laboratory measurements place further constraints on evaluating certain systems within the GHSZ; notably, CO 2 hydrate is measurably weaker than methane hydrate, and >1 order of magnitude more conductive, properties that may affect site selection, stability, and modeling considerations.

58 GEOSCIENCES↗
DOE OSTI2016

Time stretching of the GeV emission of GRBs: Fermi-LAT data versus geometrical model

It is known that the high-energy (>100 MeV) emission of gamma-ray bursts is delayed with respect to the low-energy emission. However, the dependence of light curves on energy has not been studied for the high-energy bands. In this paper, we consider the bursts observed by Fermi LAT from 2008 August 4 to 2011 August 1, for which at least 10 photons were observed with energy greater than 1 GeV. These include four bursts: GRB 080916C, GRB 090510, GRB 090902B, and GRB 090926A. We use the Kolmogorov–Smirnov test to compare the light curves in the two bands 100 MeV 1 GeV. For GRB 080916C and GRB 090510 the light curves in the two bands are statistically compatible. However, for GRB 090926A, the higher-energy light curve is stretched compared to the lower-energy one with a statistical significance of 3.3σ and, for GRB 090902B, on the contrary, the lower-energy curve is stretched with 2.3σ significance. We argue that the observed diversity of stretching factors may be explained using a simple geometrical model. The model assumes that the jet opening angle depends on the emission energy in a way that the most energetic photons are radiated near the axis of the jet. All of the bursts are considered equivalent in their rest frames, and the observed light curves differ only due to different redshifts and view directions. The model conforms to the total burst energy constraint and matches the Fermi-LAT observations of the fraction of GRBs visible in the 100 MeV<E<1 GeV band, which may be observed at higher energies. The model predicts the distribution of observable stretching factors, which may be tested in future data. Finally, we propose a method to estimate the observer's off-axis angle based on the stretching factor and the fraction of the high-energy photons. The code for modeling is open source and is publicly available on GitHub (https://github.com/maxitg/GammaRays).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

LITHIUM DEPLETION IS A STRONG TEST OF CORE-ENVELOPE RECOUPLING

Rotational mixing is a prime candidate for explaining the gradual depletion of lithium from the photospheres of cool stars during the main sequence. However, previous mixing calculations have relied primarily on treatments of angular momentum transport in stellar interiors incompatible with solar and stellar data in the sense that they overestimate the internal differential rotation. Instead, recent studies suggest that stars are strongly differentially rotating at young ages but approach a solid body rotation during their lifetimes. We modify our rotating stellar evolution code to include an additional source of angular momentum transport, a necessary ingredient for explaining the open cluster rotation pattern, and examine the consequences for mixing. We confirm that core-envelope recoupling with a ∼20 Myr timescale is required to explain the evolution of the mean rotation pattern along the main sequence, and demonstrate that it also provides a more accurate description of the Li depletion pattern seen in open clusters. Recoupling produces a characteristic pattern of efficient mixing at early ages and little mixing at late ages, thus predicting a flattening of Li depletion at a few Gyr, in agreement with the observed late-time evolution. Using Li abundances we argue that the timescale for core-envelope recoupling during the main sequence decreases sharply with increasing mass. We discuss the implications of this finding for stellar physics, including the viability of gravity waves and magnetic fields as agents of angular momentum transport. We also raise the possibility of intrinsic differences in initial conditions in star clusters using M67 as an example.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

SYGMA: Stellar Yields for Galactic Modeling Applications

The Stellar Yields for Galactic Modeling Applications (SYGMA) code is an open-source module that models the chemical ejecta and feedback of simple stellar populations (SSPs). It is intended for use in hydrodynamical simulations and semi-analytic models of galactic chemical evolution. The module includes the enrichment from asymptotic giant branch (AGB) stars, massive stars, Type Ia supernovae (SNe Ia), and compact binary mergers. An extensive and extendable stellar yields library includes the NuGrid yields with all elements and many isotopes up to Bi. Stellar feedback from mechanic and frequency-dependent radiative luminosities are computed based on NuGrid stellar models and their synthetic spectra. The module further allows for customizable initial mass functions and SN Ia delay-time distributions to calculate time-dependent ejecta based on stellar yield input. A variety of r-process sites can be included. A comparison of SSP ejecta based on NuGrid yields with those from Portinari et al. and Marigo reveals up to factors of 3.5 and 4.8 less C and N enrichment from AGB stars at low metallicity, a result we attribute to NuGrid's modeling of hot-bottom burning. Different core-collapse supernova explosions and fallback prescriptions may lead to substantial variations for the accumulated ejecta of C, O and Si in the first 107 years at Z = 0.001. Finally, an online interface of the open-source SYGMA module enables interactive simulations, analysis, and data extraction of the evolution of all species formed by the evolution of simple stellar populations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Lambda-PFLOTRAN 1.0: a workflow for incorporating organic matter chemistry informed by ultra high resolution mass spectrometry into biogeochemical modeling

Abstract. Organic matter (OM) composition plays a central role in microbial respiration of dissolved organic matter and subsequent biogeochemical reactions. Here, a direct connection of organic matter chemistry and thermodynamics to reactive transport simulators has been achieved through the newly developed Lambda-PFLOTRAN workflow tool that succinctly incorporates carbon chemistry data generated from Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) into reaction networks to simulate organic matter degradation and the resulting biogeochemistry. Lambda-PFLOTRAN is a Python-based workflow, executed through a Jupyter notebook interface, that digests raw FTICR-MS data, develops a representative reaction network based on substrate-explicit thermodynamic modeling (also termed lambda modeling due to its key thermodynamic parameter λ used therein), and completes a biogeochemical simulation with the open source, reactive flow and transport code PFLOTRAN. The workflow consists of the following five steps: configuration, thermodynamic (lambda) analysis, sensitivity analysis, parameter estimation, and simulation output and visualization. Two test cases are provided to demonstrate the functionality of the Lambda-PFLOTRAN workflow. The first test case uses laboratory incubation data of temporal oxygen depletion to fit lambda parameters (i.e., maximum utilization rate and microbial carrying capacity). A slightly more complex second test case fits multiple lambda formulation and soil organic matter release parameters to temporal greenhouse gas generation measured during a soil incubation. Overall, the Lambda-PFLOTRAN workflow facilitates upscaling by using molecular-scale characterization to inform biogeochemical processes occurring at larger scales.

58 GEOSCIENCES↗
DOE OSTI2026

ACDC (Automated Campbell Diagram Code) [SWR-26-042]

This application provides a web-based graphical user interface to generating Campbell Diagrams and visualizing mode shapes for OpenFAST turbine models. Determining the aeroelastic stability and dynamic characteristics of wind turbines is a critical step in turbine design and analysis. Historically, extracting natural frequencies and mode shapes from OpenFAST—the industry-standard whole-turbine simulation code—has been a fragmented and tedious process. It required manual model configuration, command-line linearization execution, and complex post-processing via proprietary scripts to handle rotating-frame dynamics. To address these workflow bottlenecks, we present the Automated Campbell Diagram Code (ACDC), an open-source graphical software tool developed by the National Laboratory of the Rockies (NLR) under the DOE-funded Distributed Wind Aeroelastic Modeling (dWAM) project. ACDC streamlines the end-to-end linearization and stability analysis workflow into a single, intuitive cross-platform application. The software guides users through OpenFAST model configuration, definition of operating points, and the automated execution of steady-state trim and linearization simulations. Under the hood, ACDC automates the complex mathematical post-processing steps required for rotating systems, including Multi-Blade Coordinate (MBC) transformations, eigenanalysis, and advanced modal tracking utilizing the Modal Assurance Criterion (MAC) and spectral clustering. Finally, ACDC processes these results to automatically generate Campbell diagrams and features a robust 3D visualization engine to animate full-system mode shapes. By eliminating the reliance on external post-processing environments and manual data manipulation, ACDC significantly accelerates dynamic analysis and lowers the barrier to entry for wind energy researchers and engineers.

Summerville, Brent [National Laboratory of the Roc↗
DOE OSTI2016

Generation of Low Fidelity Experimental Covariance Matrices for ICSBEP Cases

In validation, statistical methods routinely incorporate information from integral experiments to improve estimations of code biases and uncertainties. One source of high quality integral experiments is the ICSBEP Handbook. While uncertainties for criticality benchmarks are well characterized in the Handbook, information pertaining to the shared uncertainty/covariance between the cases is sparse, yet required for the aforementioned methods; correlations can result from using the same fuel for example. To remedy the lack of available covariance data in the Handbook, recent efforts have been undertaken, notably, within the framework of the NEA Expert Group on Uncertainty Analysis in Criticality Safety Assessment (UACSA) under the Working Party on Nuclear Criticality Safety (WPNCS) . The UACSA activities have focused primarily on comparing methods used to generate high quality covariance matrices for a small number of evaluations. Despite these efforts, within the Database for ICSBEP DICE only about 2% of the nearly 5000 experiments in the ICSBEP Handbook have correlation coefficients, so the amount of work required to generate the missing data appears daunting. And so, to fill in the gaps where no correlations are available, it is proposed to generate low fidelity correlations by isolating the largest components of uncertainty and determining the fraction shared between cases within an evaluation. Uncertainties can be shared when they are not independently measured, or the source of uncertainty is unchanged between cases. The degree to which an uncertainty is shared is linked to the experimental procedure; an approximation is made that in the absence of additional information, similar experimental procedures will lead to a similar fraction of shared uncertainty. Between experimental procedures there is some commonality. One method, proposed here, to help elucidate similar procedures is to encode each procedure into a decision tree grouped by uncertainty type. End nodes of the tree require judgment to assign the fraction of shared uncertainty between cases associated with the node. This can also be thought of as a set of rules, such that if experimental procedures and conditions occurred, then assign the shared uncertainty as some number between 0.0 and 1.0. One advantage of such a procedure is the formation of a transparent logical link between the judgment applied and the amount of shared uncertainty. This makes the results both more reproducible and easy to update. Generating correlation coefficients is reduced to developing decision trees that categorize how the experiments were done, and apply judgment of the degree of shared uncertainty for the end node of the tree. Trees need to be constructed after reading section 1 and 2 of the evaluation and some evaluations certainly do not contain enough information to create informative decision trees nor to assign an amount of shared uncertainty. These evaluations will not be considered here. (authors)

42 ENGINEERING↗
DOE OSTI2024

Accelerating Thermochemical Equilibrium Calculations for Nuclear Reactor Applications

Thermochemical properties play a key role in modeling and simulation of several key phenomena in nuclear reactors. There has been an increasing interest in incorporating CALPHAD-based formulations in multiphysics simulations including for Molten Salt Reactors where knowledge of phase evolution of the salt and the chemical potentials of various elements are of utmost importance in source term analyses and redox control. However, the size of such simulations is often limited by the high computational cost of full thermodynamic equilibrium calculations. This work discusses the current efforts aimed at accelerating thermochemical equilibrium calculations for multiphysics simulations performed using the open-source finite element / finite volume code Multiphysics Object Oriented Simulation Environment (MOOSE) [1]. While several methods have been proposed for accelerating phase equilibrium calculations [2], most focus on relatively small systems and often rely on a- priori knowledge of the state-space of the system. Nuclear materials, however, are often multi-component systems owing to the evolution of composition under irradiation and an approach based on a-priori mapping of phase diagram is often not enough. This work is aimed at demonstrating an on-the-fly surrogate modeling framework that uses active learning to reduce the number of full equilibrium calculations that must be performed. By combining with efficient coupling approaches, the surrogate framework helps in reducing the computational cost of thermodynamic equilibrium informed multiphysics simulations of nuclear materials. The performance is benchmarked against full coupling with the thermochemistry library Thermochimica [3]. This work uses a machine learning based approach for constructing surrogate models to predict the stable phases in a multicomponent system. The surrogates were constructed using neural networks and Gaussian process classification. In this work, we compare the relative performance of the two methods. We also demonstrate the use of caching previous calculations by interpolating the values from nearest neighbors. References [1] Lindsay, A.D., et al. "2.0 – MOOSE: Enabling massively parallel multiphysics simulation", SoftwareX, 20 (2022): 101202. [2] Roos, W.A. and Zietsman J.H. "Accelerating complex chemical equilibrium calculations – A Review", Calphad, 77 (2022): 102380. [3] Piro, M.H.A., et al. "The thermochemistry library Thermochimica", Computational Materials Science, 67 (2013): 266-272.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Tiered Energy Use in Buildings (TEB)

This tool creates a set of building energy simulations and organizes them according to a tier-based model representing circuits that will be on or off depending on a microgrid controller (not included in this software). This Python based code uses the open-source RCBuildingSimulator that is under an amended MIT license. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525. SAND2021-2528 O

Villa, Daniel↗
DOE OSTI2024

Fusion Neutron Generator

The proposed code, named FROG (Fusion neutron Generator) is built upon the open-source particle transport Monte Carlo toolkit Geant4. Geant4 provides C++ classes that can be leveraged to build application-specific codes dealing with the transport of particles through matter. Geant4-based codes are applied in high-energy particle physics experiments, medical applications, shielding, and space applications for example. The FROG code allows the user to define the geometry of a neutron converter device shaped as a hollow cylinder, where a neutron breeding material such as lithium deuteride (LiD) is cladded by two concentric cylinders. Such neutron converter is then placed inside a regular nuclear fission reactor, where thermal neutrons will react with the neutron breeder material (typically, Lithium 6), and through a series of reactions, will generate high-energy neutrons – neutrons whose kinetic energy are around 14 MeV. The hollowed central portion can hold a specimen that will be bombarded by high-energy neutrons created inside the neutron breeding material. Figuratively speaking, this type of device transforms neutrons from thermal (~0.625 eV) to fusion (~14 MeV) energies and is sometimes termed “fusion-to-thermal neutron converters” in the literature. The code consists of C++ source file compiled and linked to generate an executable. The user can select the dimensions of the converter (radius, length, and thickness of the breeder material), the breeder material type, the cladding material, and the specimen material that will be activated or irradiated. As input, the neutron flux for a specific location inside a reactor, for instance, positions in ATR, is required. As output, the code predicts the number of high-energy neutrons produced, the total neutron flux and fluence as well as its detailed spectrum. The physics involved in such device is very complex, as it requires modeling neutron transport, light-ion (tritons) transport, as well as fusion reactions. The Geant4 toolkit provides the required physical models.

Martin, NicholasP. [Idaho National Laboratory (INL↗
DOE OSTI2023

Applications of “Synchrotron Radiation Workshop” Code

The “Synchrotron Radiation Workshop” (SRW) computer code was started in the group of Insertion Devices of the European Synchrotron Radiation Facility in 1997. The idea behind it was to combine in one software library electrodynamics methods for calculating electric fields (and derived characteristics) of synchrotron radiation (SR) generated by relativistic electrons moving in magnetic fields of nearly arbitrary configuration with wave-optics methods for simulating propagation of radiation through optical elements of beamlines. Further, the use of this general approach was initially driven by applications in electron beam diagnostics using SR. However, with time it started to be used for more complicated beamlines of synchrotron light sources, and the use of the code for these applications continues to increase, with the new light sources offering brighter and more coherent x ray beams for wide ranges of experiments in different areas of science.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

What can we learn from in-soil imaging of a live plant: X-ray Computed Tomography and 3D numerical simulation of root-soil system

Plant roots play a critical role in plant-soil-microbe interactions that occur in the rhizosphere, as well as processes with important implications to farming, forest management and climate change. X-ray computed tomography (XCT) has been proven to be an effective tool for non-invasive root imaging and analysis. A combination of XCT, open-source software, and our own code was used to noninvasively image a prairie dropseed (Sporobolus heterolepis) specimen, segment the root data to obtain a 3D image of the root structure at 31µm resolution, and extract quantitative information (root volume and surface area) from the 3D data, respectively. Based on the mesh generated from the root structure, computational fluid dynamics (CFD) simulations were applied to numerically investigate the root-soil-groundwater system. The plant root conductivity, soil hydraulic conductivity and transpiration rate were shown to control the groundwater distribution. The flow variability and soil water distributions under different scenarios were investigated. Parameterizations were evaluated to show their impacts on the average conductivity. The pore-scale modeling approach provides realistic simulations of rhizosphere flow processes and provides useful information that can be linked to upscaled models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

SU-E-T-605: A New Design for a Rotating Gamma Knife. Monte Carlo Simulation

Purpose: to determine the characteristics of the 60Co beam emerging from a new design of a rotating Gamma Knife system and to calculate 3D dose distributions at the isocenter for different source configurations and collimator openings. Methods: We employed the BEAM-Monte Carlo code to realistically model the geometry design, including 30 60Co source capsules, two circular primary collimators (diameter of 6.6mm and 6.1mm) and four different changeable collimators. The shielding of the head design was also simulated. The sources (2.8mm diameter) are distributed in six groups in the spherical geometry. Each source is individually collimated to obtained four different circular fields at the isocenter (3mm, 3.5mm, 6mm and 8mm).The phase-space particles reaching the scoring plane below the primary collimation assembly were recorded and the BEAMDP code was used to determine the fluence and energy spectra of the particles emerging from each source-collimator configuration. The dose distributions at the isocenter plane (397.6mm from the source) were calculated in a spherical component module for the circular field sizes studied. Results: The energy spectra below the head assembly and primary collimators have been obtained, which exhibited the typical 60Co peaks and a small low-energy tail due to scattered photons (from about 200keV to 1MeV). The scattered component of the spectra represents about 8 % of the total number of photons reaching the scoring plane. The radial photon fluence does not vary significantly inside the collimator openings. The spectra of particles from different source groups are compared. Conclusion: The 60Co beam emerging from each source configuration was characterized, which can be used to establish a generic source model for all the sources for fast MC dose calculation. Further investigations are needed to determine the dose variations as a result of partial switching on/off different groups of sources for advanced Gamma Knife SRS/SBRT planning.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

Manufactured solutions for an electromagnetic slot model

The accurate modeling of electromagnetic penetration is an important topic in computational electromagnetics. Electromagnetic penetration occurs through intentional or inadvertent openings in an otherwise closed electromagnetic scatterer, which prevent the contents from being fully shielded from external fields. To efficiently model electromagnetic penetration, aperture or slot models can be used with surface integral equations to solve Maxwell's equations. A necessary step towards establishing the credibility of these models is to assess the correctness of the implementation of the underlying numerical methods through code verification. Surface integral equations and slot models yield multiple interacting sources of numerical error and other challenges, which render traditional code-verification approaches ineffective. In this paper, we provide approaches to separately measure the numerical errors arising from these different error sources for the method-of-moments implementation of the electric-field integral equation with a slot model. Finally, we demonstrate the effectiveness of these approaches for a variety of cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗