Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “3D model”

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 559 records · Page 31

STELLOPT Modeling of the 3D Diagnostic Response in ITER

The ITER three dimensional diagnostic response to an n=3 resonant magnetic perturbation is modeled using the STELLOPT code. The in-vessel coils apply a resonant magnetic perturbation (RMP) fi eld which generates a 4 cm edge displacement from axisymmetry as modeled by the VMEC 3D equilibrium code. Forward modeling of flux loop and magnetic probe response with the DIAGNO code indicates up to 20 % changes in measured plasma signals. Simulated LIDAR measurements of electron temperature indicate 2 cm shifts on the low field side of the plasma. This suggests that the ITER diagnostic will be able to diagnose the 3D structure of the equilibria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of three-dimensional models for WIPP performance assessment using PFLOTRAN - 15073

The Waste Isolation Pilot Plant (WIPP), located in southeastern New Mexico, has been developed by the U.S. Department of Energy (DOE) for the deep geologic disposal of transuranic (TRU) waste. WIPP performance assessment (PA) calculations estimate the probability and consequence of potential radionuclide releases from the repository to the accessible environment for a regulatory period of 10,000 years after facility closure. To date, WIPP PA calculations have employed multiple 2D numerical models requiring simplification of the mesh and processes. Introducing 3D numerical models within WIPP PA enables increasingly realistic representations of the WIPP subsurface domain and improved flexibility for incorporating relevant features. Furthermore, the increased computing capacity afforded by modern hardware architectures enables this 3D modeling. Modeling in 3D will be more flexible because less homogenization is required and better accuracy of the results can be obtained by making fewer assumptions and/or simplifications. The use of PFLOTRAN, a state-of-art massively parallel subsurface flow and reactive transport code, will enable this 3D modeling, and eliminate the need for multiple related, but decoupled 2D models. However, there are capabilities that PFLOTRAN must be demonstrated in order to fulfill requirements for WIPP PA. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Precipitation Processes developed during ARM (1997), TOGA COARE (1992), GATE (1974), SCSMEX (1998), and KWAJEX (1999), Consistent 2D, semi-3D and 3D Cloud Resolving Model Simulations

Real clouds and cloud systems are inherently three-dimensional (3D). Because of the limitations in computer resources, however, most cloud-resolving models (CRMs) today are still two-dimensional (2D). A few 3D CRMs have been used to study the response of clouds to large-scale forcing. In these 3D simulations, the model domain was small, and the integration time was 6 hours. The major objectives of this paper are: (1) to assess the performance of the super-parameterization technique (i.e. is 2D or semi-3D CRM appropriate for the super-parameterization?); (2) calculate and examine the surface energy (especially radiation) and water budgets; (3) identify the differences and similarities in the organization and entrainment rates of convection between simulated 2D and 3D cloud systems.

Tao, Wei-Kuo↗

Energy Performance Assessment of Radiant Cooling System through Modeling and Calibration at Component Level

The paper describes a case study of an information technology office building with a radiant cooling system and a conventional variable air volume (VAV) system installed side by side so that performancecan be compared. First, a 3D model of the building involving architecture, occupancy, and HVAC operation was developed in EnergyPlus, a simulation tool. Second, a different calibration methodology was applied to develop the base case for assessing the energy saving potential. This paper details the calibration of the whole building energy model to the component level, including lighting, equipment, and HVAC components such as chillers, pumps, cooling towers, fans, etc. Also a new methodology for the systematic selection of influence parameter has been developed for the calibration of a simulated model which requires large time for the execution. The error at the whole building level [measured in mean bias error (MBE)] is 0.2%, and the coefficient of variation of root mean square error (CvRMSE) is 3.2%. The total errors in HVAC at the hourly are MBE = 8.7% and CvRMSE = 23.9%, which meet the criteria of ASHRAE 14 (2002) for hourly calibration. Different suggestions have been pointed out to generalize the energy saving of radiant cooling system through the existing building system. So a base case model was developed by using the calibrated model for quantifying the energy saving potential of the radiant cooling system. It was found that a base case radiant cooling system integrated with DOAS can save 28% energy compared with the conventional VAV system.

Building Energy Simulation↗

On the Origin of Sinter‐Resistance and Catalyst Accessibility in Raspberry‐Colloid‐Templated Catalyst Design

Abstract Nanoparticle (NP) sintering is a major cause of the deactivation of supported catalysts. Raspberry‐colloid‐templated (RCT) catalysts are an emerging class of materials that show an unprecedented level of sinter‐resistance and exhibit high catalytic activity. Here a comprehensive study of the origin of NP stability and accessibility in RCT catalysts using theoretical modeling, 3D electron microscopy, and epitaxial overgrowth is reported. The approach is showcased for silica‐based RCT catalysts containing dilute Pd‐in‐Au NPs previously used in hydrogenation and oxidation catalysis. Modeling of the contact line of the silica precursor infiltrating into the assembled raspberry colloids suggests that a large part of the particles must be embedded into silica, which is confirmed by quantitative visualization of >200 individual NPs by dual‐axis electron tomography. The RCT catalysts have a unique structure in which all NPs reside at the pore wall but have >50% of their surface embedded in the matrix, giving rise to the strongly enhanced thermal and mechanical stability. Importantly, epitaxial overgrowth of Ag on the supported NPs reveals that not only the NP surface exposed to the pore but the embedded interface as well remained chemically accessible. This mechanistic understanding provides valuable guidance in the design of stable catalytic materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Validated Reactive Force Field Quantifies MXene Interfacial Properties, Mechanics, and Thermal Transport

MXenes combine rich surface chemistry, mechanical strength, and high conductivity for a multitude of emerging applications. Predictive modeling supports accelerated materials designs and has been limited by the absence of validated and transferable force fields. Here, we introduce an interpretable, reactive INTERFACE force field (IFF and IFF-R) for Ti 3 C 2 T x MXenes that is trained based on chemical knowledge and achieves quantitative agreement with experiments across lattice parameters (<0.5%), density (<0.2%), liquid contact angles, Raman spectra, and the in-plane elastic modulus (∼320 GPa). The models cover surface terminations from hydroxyl (−OH) to fluorine (−F) groups and are extensible to other chemistries. We introduce pH-resolved surface chemistry and identify dopamine adsorption mechanisms at MXene–aqueous interfaces supported by QCM-D and UV–Vis experiments. The data reveal coplanar and perpendicular binding modes and concentration-dependent multilayer assembly. We predict previously inaccessible properties, including termination-dependent cleavage energies, interlayer shear moduli and dynamic shear failure, nanoindentation and brittle fracture, anisotropic in-plane and out-of-plane thermal conductivities, including the role of defects. Agreement with available experimental data is consistently close and exceeds DFT accuracy across the benchmark properties examined. The IFF/IFF-R model is compatible with CHARMM, AMBER, OPLS, and CVFF force fields for simulations of MXenes with diverse surface terminations, electrolyte interfaces, biointerfaces, and polymer composites without additional parameters. Parameter sets, 3D models, and analysis scripts are provided for community use. The validated, reactive, and transferable IFF framework facilitates predictive design of MXene-based films, membranes, sensing interfaces, and composites.

MXene↗

A virtual simulator designed for collision prevention in proton therapy

Purpose: In proton therapy, collisions between the patient and nozzle potentially occur because of the large nozzle structure and efforts to minimize the air gap. Thus, software was developed to predict such collisions between the nozzle and patient using treatment virtual simulation. Methods: Three-dimensional (3D) modeling of a gantry inner-floor, nozzle, and robotic-couch was performed using SolidWorks based on the manufacturer’s machine data. To obtain patient body information, a 3D-scanner was utilized right before CT scanning. Using the acquired images, a 3D-image of the patient’s body contour was reconstructed. The accuracy of the image was confirmed against the CT image of a humanoid phantom. The machine components and the virtual patient were combined on the treatment-room coordinate system, resulting in a virtual simulator. The simulator simulated the motion of its components such as rotation and translation of the gantry, nozzle, and couch in real scale. A collision, if any, was examined both in static and dynamic modes. The static mode assessed collisions only at fixed positions of the machine’s components, while the dynamic mode operated any time a component was in motion. A collision was identified if any voxels of two components, e.g., the nozzle and the patient or couch, overlapped when calculating volume locations. The event and collision point were visualized, and collision volumes were reported. Results: All components were successfully assembled, and the motions were accurately controlled. The 3D-shape of the phantom agreed with CT images within a deviation of 2 mm. Collision situations were simulated within minutes, and the results were displayed and reported. Conclusions: The developed software will be useful in improving patient safety and clinical efficiency of proton therapy.

60 APPLIED LIFE SCIENCES↗

Characterizing Structural Controls of EGS-Candidate and Conventional Geothermal Reservoirs in the Great Basin: Developing Successful Exploration Strategies in Extended Terranes

We conducted a comprehensive analysis of the structural controls of geothermal systems within the Great Basin and adjacent regions. Our main objectives were to: 1) Produce a catalogue of favorable structural environments and models for geothermal systems. 2) Improve site-specific targeting of geothermal resources through detailed studies of representative sites, which included innovative techniques of slip tendency analysis of faults and 3D modeling. 3) Compare and contrast the structural controls and models in different tectonic settings. 4) Synthesize data and develop methodologies for enhancement of exploration strategies for conventional and EGS systems, reduction in the risk of drilling non-productive wells, and selecting the best EGS sites.

15 GEOTHERMAL ENERGY↗

3D Radiative Hydrodynamics Modeling of Convection of Stars to Probe Their Interiors and Photospheric Properties

The dramatic flow of data from the Kepler and K2 missions opens the opportunity to significantly improve our knowledge of stellar interiors, surface dynamics, and structure. However, interpretation of these observations is a challenging task because it depends on tiny effects that can be studied only with advanced first-principles modeling. We present results of 3D time-dependent radiative hydrodynamic simulations of stellar outer convection zones and atmospheres taking into account chemical composition, radiative transfer, turbulence effects, and a realistic equation of state for main sequence stars. We will discuss properties of convective structure and dynamics, convective overshoot, effects of magnetic fields and rotation, as well as the potential influence of turbulent surface dynamics on high-precision RV measurements.

Photospheric↗

Simulation-Aided Design and Interpretation of Results of Microgravity-Based Solidification Experiments

The preparation and interpretation of microgravity experiments conducted on the International Space Station (ISS) can be challenging because of the remote location of the equipment. Although a “ground unit” is typically available on Earth, it is not located in the laboratory of the principal investigator. Furthermore, test experiments in the ground unit cannot fully duplicate heat transfer conditions of the weightless environment, and thus be a reliable reference for experiment design. Under such conditions, numerical simulations are a valuable tool both for preparing the experiments and interpreting the results obtained under terrestrial and micro-gravity conditions. In this work, we developed a detailed 3D model of the “SUBSA” furnace located at the ISS, and demonstrate its use in our present investigation which includes melt growth and vapor growth of InI crystals on the ISS. In 2002, the SUBSA furnace was designed and built for the growth of InSb crystals in the Microgravity Science Glovebox (MSG) at the ISS. The SUBSA furnace features a transparent section without heating elements, where the axial temperature gradient is high. In our present investigation, we have used the same furnace to grow four InI crystals from the melt, and two InI crystals from the vapor phase (2018-2021). InI crystals have ~20 times lower thermal conductivity than InSb crystals. As a result, during melt growth, the InI crystals were exposed to extremely high temperature gradients. To reduce the axial temperature gradient, we added external jackets made from high thermal conductivity material to the growth ampoules. The SUBSA furnace could not be modeled as asymmetric because of its rectangular shape and right-angled transparent section, which contains a convoluted geometry. The rectangular shape of the transparent section further complicates modelling of radiative heat transfer, which dominates in the zone. At present, numerical simulations are being conducted to resolve an apparent discrepancy between the readings of the 5 thermocouples located in the SUBSA furnace, and the observed position of the solid and molten material. The 3D simulation will be presented, along with the calibration and growth experiment results conducted in microgravity and in the SUBSA ground unit. The developed model should be useful for upcoming microgravity investigations planned to be conducted in the SUBSA facility.

V Riabov↗

CFD modeling of debris melting phenomena during late phase Candu 6 severe accident

The objective of this paper was to study the phase change of the debris formed on the Candu 6 calandria bottom in a postulated accident sequence. The molten pool and crust formation were studied employing the Ansys-Fluent code. The 3D model using Large Eddy Simulation (LES) predicts the conjugate, radiative and convective heat transfer inside and from the corium pool. LES simulations require a very fine grid to capture the crust formation and the free convection flow. This aspect (fine mesh requirement) correlated with the long transient has imposed the use of a slice from the 3D calandria geometry in order not to exceed the computing resources. The preliminary results include heat transfer coefficients, temperature profiles and heat fluxes through calandria wall. From the safety point of view it is very important to maintain a heat flux through the wall below the CHF assuring the integrity of the calandria vessel. This can be achieved by proper cooling of the tank water which contains the vessel. Also, transient duration can be estimated being important in developing guidelines for severe accidents management. The debris physical structure and material properties have large uncertainties in the temperature range of interest. Thus, further sensitivity studies should be carried out in order to better understand the influence of these parameters on this complex phenomenon. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

SU-F-BRD-11: A Virtual Simulator Designed for Collision Prevention in Proton Therapy

Purpose: In proton therapy, collisions between patient and nozzle potentially occur in attaining minimal air gap due to the large nozzle structure. Thus, we developed software predicting the collisions of the nozzle and patient by simulating treatments. Methods: 3D modeling of a gantry inner-floor, nozzle and robotic-couch was done by using the SolidWorks based on the manufacturer’s machine data. To obtain patient body information, a 3D-scanner was utilized to scan a patient right before CT scanning. From the acquired images, a 3D-image of the patient’s body contour was reconstructed. The accuracy of the image was confirmed against the CT image for a humanoid phantom. The machine components and the virtual patient were combined on the treatment-room coordinate system, resulting in a virtual simulator. The simulator simulated the motion of its components such as rotation and translation of gantry, nozzle and couch, in real scale. Collision, if any, was examined both in static mode and dynamic mode. The static mode checks only at fixed positions of the machine’s components while dynamic mode examines while one component is in motion. Collision was notified if any voxel of two components, for example a nozzle and a patient or couch, overlapped when calculating volume locations. The event and collision point are visualized and colliding volumes are reported. Results: All components were successfully assembled and the motions could be accurately controlled. The 3D-shape of a phantom agreed with CT images within a deviation of 2 mm. Collision situations can be simulated within minutes and the results are displayed and reported. Conclusion: The developed software will be useful in improving patient safety and clinical efficiency for proton therapy. This work was supported by the National Research Foundation of Korea funded by Ministry of Science, ICT & Future Planning (2012M3A9B6055201, 2013M2A2A7043507), and Samsung Medical Center grant (GFO1130081)

60 APPLIED LIFE SCIENCES↗

WE-G-17A-09: Novel Magnetic Shielding Design for Inline and Perpendicular Integrated 6 MV Linac and 1.0 T MRI Systems

Purpose: The influence of fringe magnetic fields delivered by magnetic resonance imaging (MRI) on the beam generation and transportation in Linac is still a major challenge for the integration of linear accelerator and MRI (Linac-MRI). In this study, we investigated an optimal magnetic shielding design for Linac-MRI and further characterized the beam trajectory in electron gun. Methods: Both inline and perpendicular configurations were analyzed in this study. The configurations, comprising a Linac-MRI with a 100cm SAD and an open 1.0 T superconductive magnet, were simulated by the 3D finite element method (FEM). The steel shielding around the Linac was included in the 3D model, the thickness of which was varied from 1mm to 20mm, and magnetic field maps were acquired with and without additional shielding. The treatment beam trajectory in electron gun was evaluated using OPERA 3d SCALA with and without shielding cases. Results: When Linac was not shielded, the uniformity of diameter sphere volume (DSV) (30cm) was about 5 parts per million (ppm) and the fringe magnetic fields in electron gun were more than 0.3 T. With shielding, the magnetic fields in electron gun were reduced to less than 0.01 T. For the inline configuration, the radial magnetic fields in the Linac were about 0.02T. A cylinder steel shield used (5mm thick) altered the uniformity of DSV to 1000 ppm. For the perpendicular configuration, the Linac transverse magnetic fields were more than 0.3T, which altered the beam trajectory significantly. A 8mm-thick cylinder steel shield surrounding the Linac was used to compensate the output losses of Linac, which shifted the magnetic fields' uniformity of DSV to 400 ppm. Conclusion: For both configurations, the Linac shielding was used to ensure normal operation of the Linac. The effect of magnetic fields on the uniformity of DSV could be modulated by the shimming technique of the MRI magnet. NIH/NIGMS grant U54 GM104944, Lincy Endowed Assistant Professorship.

43 PARTICLE ACCELERATORS↗

Risk-Informed Systems Analysis (RISA) Dynamic Fire PRA Roadmap

Modeling and implementing fire safety for nuclear power plants is a costly activity. Because of the complexity of fire phenomenon and multiple operational procedures, it is difficult to computationally provide assurance that the mitigation methods are adequate for critical areas using current analysis methods. An economical method to provide more accurate modeling and optimize mitigation methods is needed to improve nuclear power viability. This report describes the initial investigation into modeling and simulation tools for application of fire as part of the Risk-Informed Systems Analysis (RISA), formerly Risk-Informed Safety Margin Characterization (RISMC). The report provides a framework of how 3D modeling and simulation techniques could be combined with a dynamic Probabilistic Risk Analysis (PRA) to reduce compounding conservatism present in current Fire PRA methods. Electrical Power Research Institute has analyzed current Fire PRA practices and identified the most significant contributors to risk and areas for improving analysis. This framework describes how to apply dynamic PRA and simulation methods for key contributors/scenarios, credit missing factors of manual fire suppression, and applying conditional probabilities.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Hawai‘i Play Fairway (Final Report)

Most of Hawai'i's geothermal resources are blind—their manifestations, such as hot springs and steam vents, do not appear on the ground surface because the heated water flows far below. With the exception of $K\bar{i}lauea East Rift Zone$, in most areas of Hawai'i, high lateral permeability in the first kilometer below ground surface prevents surface thermal features from developing. As a methodology for discovering these blind resources, Play Fairway Analysis (PFA) involves finding potential locations of blind hydrothermal systems and describing potential geothermal sources in rift-zone settings. Using the PFA to find Hawai'i's geothermal resources, the University of Hawai'i (UH) conducted the Hawai‘i Play Fairway Project, Hawai'i's first statewide geothermal resource assessment since 1985. Sponsored by the U.S. Department of Energy, the Hawai'i Play Fairway Project provided an updated resource assessment, a roadmap for additional exploration activities, and the identification of areas for further exploration. Benefitting from UH's core competency in earth sciences and experienced geothermal researchers, the project comprised three phases. During the first phase, the team identified, compiled, and ranked existing geologic, groundwater, and geophysical datasets relevant to subsurface heat, fluid and permeability. Using a Bayesian statistical approach, the team developed a statistical methodology to integrate these data into a resource probability map. The team evaluated the confidence in the probability value and considered development viability of areas with geothermal resources. With these analyses, the team identified 10 locations in the Hawaiian Islands for exploration activities. For the second phase, the team collected new groundwater data in 10 locations across the state and new geophysical data on $L\bar{a}na‘i, Maui$, and central Hawai'i Island and modeled topographically induced stress to better characterize subsurface permeability. Analyzing the subsurface stresses, the team evaluated the potential for fracture-induced permeability. The team inverted the MT and gravity data to produce 3D models of resistivity and density, respectively, on $L\bar{a}na‘i$, across $Haleakal\bar{a}'s$ SW rift (Maui), and surrounding Mauna Kea (Hawai‘i Island). The team developed and applied a new method for incorporating depth information about resistivity, density, and potential for fracture-induced permeability into the statistical method for computing resource probability in these three focus areas. The team incorporated the new groundwater results with the new geophysical results and the calculations of potential for fracture-induced permeability to produce updated maps of resource probability and confidence. Through combining data from the first and second phases, the team determined locations for further exploration during the third phase. For MT and gravity surveys, the team recommended $Kaua'i's$ $L\bar{i}hu'e$ $Basin$, the east rift of $Maui's$ $Haleakal\bar{a}$ volcano, and the southwest rift of Hawai'i Island's Mauna Loa volcano. The MT and gravity surveys aimed to enable improved confidence in the resource potential in these locations. For drilling deep groundwater well(s), the team recommended Southeast Mauna Kea and $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$. During the third phase, further exploration involved drilling a groundwater well in $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$ and performing more geophysical surveys. We deepened an existing water well proximal to our target area on $L\bar{a}na'i$ due to funding constraints that precluded us from spudding a new well that would exceed 1km depth. Drilling was preceded by a number of substantial elements including: writing an Environmental Assessment and the subsequent legal process, performance of deviation logging, lowering a camera down the well, coordinating site preparation with $P\bar{u}lama$ $L\bar{a}na'i$, shipping the UH-owned rig interisland, procuring supplies, and leading 3 community meetings on $L\bar{a}na'i$. Drilling occurred 24/7 the entire month of June 2019 over which time $L\bar{a}na'i$ $Well$ 10 was deepened from 427 m to 1057 m, with continuous core collected. We measured a roughly linear temperature gradient averaging 42°C/km and a maximum bottom hole temperature of 66°C. This gradient is more than twice the background for Hawai'i and within a range of gradients measured in this depth range for some exploration wells within KERZ. We consider these results encouraging for $L\bar{a}na'i's$ resource potential and recommend following with a slim hole within $L\bar{a}na'i's$ caldera (our target zone) to ~ 2 km. Further, the positive implications such results have for the island of O‘ahu are substantial - the shield stage of O'ahu's volcanoes ended 1-2 My earlier. However, O'ahu uses more electricity than the rest of the islands combined, and the utility recently called for 500-700MW of firm, dispatchable renewable electricity on O'ahu by 2033. In Phase 3, we also collected limited new encouraging groundwater data, and updated our thoughts on the probabilities of fluid and permeability at resource depths (Pr F = 1; Pr P = mostly unconstrained). Ultimately, we advocate for using our final probability of heat, and confidence in this probability, to drive the next phase of exploration. We contend further development of geothermal in Hawai‘i will enable the state to achieve its 100% renewable policy objective and Hawai'i to transition off of fossil fuels through geothermal discovery and development. The project not only produced a large amount of data and expanded the existing knowledge of Hawai'i's geothermal resources, but also produced publications, theses, presentations, core photos, datasets, media reports, television interviews, community events, and a blog. Students and new professionals benefitted from the project's hands-on research experiences and educational opportunities and earned awards and recognition.

15 GEOTHERMAL ENERGY↗

Initial 1.5D BISON Simulation Verification and Validation CASL FY17 Letter Report

This milestone report describes the current efforts to improve the usability of the recently developed 1.5 dimensional capability in BISON, the initial verification of this BISON capability, and initial validation of the 1.5 dimensional models against existing 2 dimensional axisymmetric full length light water reactor (LWR) fuel rods. The 1.5 dimensional (1.5D) capability in BISON was developed to enable the simulation of LWR fuel rods in a less computationally demanding simulation formulation while retaining key physics of interest. BISON was originally developed to solve 2D axisymmetric and 3D models, and this original development focus allows BISON to examine a wider range of fuel types and geometrically complexity than can be solved with lower dimensional frameworks. As the set of problems which BISON is used to solve expands, the software design requirements continue to evolve. Coupling physics in a higher dimensional geometry is not always required to resolve the physics of interest, such as the modeling of radial heat transfer and radial displacement. BISON, built on the dimension-agnostic MOOSE framework, is able to solve problems in 3D, 2D, and 1D with only minor modifications to the material definition code so that the appropriate geometrical representation can be applied in different problems of interest. The motivation for the development and validation of the BISON 1.5D capability is the CASL program's need for a simplified and computationally efficient version of BISON for use in the large multi-physics VERA simulations which include hundreds of individual fuel rods. The original development of the BISON 1.5D simulation capability is described in the CASL FY17 Letter Report from January 2017, and this report documents the continuing work identified in the previous report. The improvements to the usability of the 1.5D capability through input file simplifications, the results of the verification simulations, the results of the initial validation against the 2D assessment cases in BISON, and the future development needs are included within this report.

1.5D Capability↗

When Disorder Looks Like Order: A New Model to Explain Radial Magnetic Fields in Young Supernova Remnants

Radial magnetic fields are observed in all known young, shell-type supernova remnants in our Galaxy, including Cas A, Tycho, Kepler, and SN1006, and yet the nature of these radial fields has not been thoroughly explored. Using a 3D model, we consider the existence and observational implications of an intrinsically radial field. We also present a new explanation of the origin of the radial pattern observed from polarization data as resulting from a selection effect due to the distribution of cosmic-ray electrons (CREs). We show that quasi-parallel acceleration can concentrate CREs at regions where the magnetic field is radial, making a completely turbulent field appear ordered, when it is in fact disordered. We discuss observational properties that may help distinguish between an intrinsically radial magnetic field and the case where it only appears radial due to the CRE distribution. We also show that the case of an intrinsically radial field with a quasi-perpendicular CRE acceleration mechanism has intriguing similarities to the observed polarization properties of SN1006.

79 ASTRONOMY AND ASTROPHYSICS↗

Subject-specific modeling framework for particle deposition using computational fluid dynamics

Quantifying particle deposition and dose in the respiratory tract requires a physiologically realistic representation and reproducible computational workflows. However, existing modeling frameworks, such as the International Commission on Radiological Protection (ICRP) compartmental models and the Multiple Path Particle Dosimetry (MPPD) tool, lack detailed deposition profiles and subject-specific capabilities. The combination of advances in computer vision algorithms applied to the respiratory tract and Computational Fluid and Particle Dynamics (CFPD) allows high-fidelity simulations of particle behavior in anatomically accurate geometries derived from individual CT scans. The segmentation, preprocessing, and file preparation task for a CFPD simulation was often time-consuming, and no prior studies to-date have yet presented a fully automated framework. This work presents a fully automated workflow to obtain individualized particle deposition profiles in the human respiratory tract. The pipeline starts with segmenting upper and lower airway geometries using morphological and deep learning-based methods, generating three-dimensional (3D) models from CT imaging data. Next, a series of algorithms are presented to quality check and prepare the 3D geometry for a CFD or CFPD simulation. The preprocessing step includes correcting geometric artifacts, enforcing a physically consistent mesh, and automatically identifying and capping multiple outlets, which is required for CFD/CFPD simulations. These processed models are then input into open-source (OpenFOAM) or commercial (StarCCM+) CFD solvers, where flow and transient particle transport equations — including turbulence and particle–wall interactions are solved under realistic breathing conditions. Finally, the resulting particle deposition profiles can be integrated with Monte Carlo radiation transport codes and state-of-the-art computational phantoms to assess organ-specific absorbed doses in scenarios of radioactive aerosol inhalation. The presented work streamlines respiratory tract segmentation, preprocessing for CFD/CFPD simulations, and integration with dose assessment workflows, reducing manual intervention and improving access to high-fidelity, subject-specific modeling. The high precision in predicted particle deposition and dose distributions can improve personalized treatment strategies in respiratory medicine and refine dose estimates for radiation protection.

AI↗