SEARCH · Search NASA

Results for “Code Robustness”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 667 records · Page 37

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Enhancement of risk informed validation framework for external hazard scenario

In recent years, the U.S. Nuclear Regulatory Commission (USNRC) and the International Atomic Energy Agency (IAEA) have developed methodologies to assess the vulnerabilities of nuclear plants against site specific extreme hazards. In many cases, advanced simulation tools are being considered to simulate multi-physics, multi-scale phenomena and to evaluate vulnerability of nuclear facilities. The credibility of advanced simulation tools is assessed based on a formal verification, validation, and uncertainty quantification procedure. One of the key limitations in validation is the lack of relevant experimental data at system-level. This limitation leads to a decrease in the confidence of system-level risk predictions. Therefore, a robust validation framework is needed to formalize the confidence in predictive capability of advanced simulation results. Additionally, this study enhances the existing risk informed validation methodology, originally proposed by Kwag et al. [1] and Bodda et al. [2], by developing additional attributes and a new set of validation indicies for a complete and wider applicability of the framework. In this manuscript, the methodology to identify the critical path that leads to the system-level failure is illustrated. The overall validation is checked for completeness and consistency by comparing the critical path for both the system-level simulation and experimental models. The applicability of the code for an intended application is represented in terms of various maturity levels and helps in the process of decision making.

42 ENGINEERING↗

DOE OSTI2020

TEAL

TEAL is a financial performance calculator plugin for the RAVEN code, framework, resolving around the computation of Net Present Value and associated financial metrics. TEAL can make use of inflation rates, taxation, escalation factors, capital expenditure economy of scale scaling factors. The unique feature of TEAL is the capability to be linked with RAVEN external models and build corresponding cash flows using the variables computed by those external models. In addition to be able to use the capability to generate cash flows derived from complex physical models generated by RAVEN, another distinctive feature of TEAL is the capability to provide financial risk/probabilistic metrics that can empower RAVEN to perform optimization/analysis driven by financial risk augmentations. Optimization, robust optimization, parametric studies, large parallel simulations, sensitivity analysis, data mining, etc. are just some of the capabilities that can be leveraged.

Alfonsi, Andrea↗

DOE OSTI2018

(Machine) learning to do more with less

Determining the best method for training a machine learning algorithm is critical to maximizing its ability to classify data. In this paper, we compare the standard “fully supervised” approach (which relies on knowledge of event-by-event truth-level labels) with a recent proposal that instead utilizes class ratios as the only discriminating information provided during training. This so-called “weakly supervised” technique has access to less information than the fully supervised method and yet is still able to yield impressive discriminating power. In addition, weak supervision seems particularly well suited to particle physics since quantum mechanics is incompatible with the notion of mapping an individual event onto any single Feynman diagram. We examine the technique in detail — both analytically and numerically — with a focus on the robustness to issues of mischaracterizing the training samples. Weakly supervised networks turn out to be remarkably insensitive to a class of systematic mismodeling. Furthermore, we demonstrate that the event level outputs for weakly versus fully supervised networks are probing different kinematics, even though the numerical quality metrics are essentially identical. This implies that it should be possible to improve the overall classification ability by combining the output from the two types of networks. For concreteness, we apply this technology to a signature of beyond the Standard Model physics to demonstrate that all these impressive features continue to hold in a scenario of relevance to the LHC. Example code is provided on GitHub.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Investigation into Scalable and Detection-Enhanced Satellite Conjunction Assessment

Imaging opportunities (viewable conjunctions) of Resident Space Objects (RSOs) by satellites are not continuously discovered. We propose to continuously produce and report viewable conjunctions among objects in orbit. Viewable conjunctions are events in space and time when a satellite may favorably view a Resident Space Object (RSO). Favorability is defined by a set of constraints, e.g., solar illumination, distance between observer and target, orbital location for viewable event. Computing viewable conjunctions requires calculation of orbital propagation while considering constraints based on the state vectors of position, velocity, with covariance for both satellite and RSO. We propose two parallel lanes of effort: acceleration and research. The objective of acceleration is to avoid missed opportunities and reduce latency for satellite maneuver requests through continuous prediction and reporting of viewable conjunctions. The effort will begin by deploying currently available software on dedicated systems and continue with optimizing the code for high performance computing hardware. The research lane aims to expand RSO inspection and modeling capabilities. Among our current research ideas are spectral characterization of RSO materials and planning multiple observations to recover RSO 3D form. Computing resources at Oak Ridge National Laboratory (ORNL) are available for the acceleration work. Laika, Maxar conjunction prediction dashboard software, and Bluesim, Maxar orbital propagation software, are expected to be the first software in the acceleration lane. Laike and Bluesim are to be provided by the sponsor, and output will be made accessible through its dashboard. Deliverables will follow a gated schedule to the sponsor. ORNL will provide progressively more robust viewable conjunction assessments from both modelled and actual ephemerides.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Smart Preprocessing & Robust Integration Emulator

To achieve the desired particle size of biomass feedstocks during preprocessing for trouble-free handling and conversion to produce biofuels and bioproducts, the raw materials must undergo a crucial milling process. The particle size of biomass plays a critical role in subsequent biofuel manufacturing, where a larger area-to-volume ratio facilitates efficient synthesis while balancing the impact of moisture on biomass storage. To optimize biofuel production efficiency and overcome these challenges, it is imperative to accurately predict the particle size distribution (PSD) of the biomass in the design of efficient preprocessing systems. The population balance model (PBM), upon empirical calibration and validation, can provide rapid prediction of post-milling PSD of granular biomass. However, PSD has limitations related to mass conservation and the absence of moisture considerations. To overcome these drawbacks, a deep learning model called the enhanced deep neural operator (DNO+) is implemented in the code. This model not only retains the capabilities of the PBM in handling complex mapping functions but also incorporates additional factors influencing the system. By considering various experimental conditions such as sieve size and moisture content, the trained DNO+ model can effectively predict the PSD after milling for any given feed PSD. To further reduce the reliance on experimental data, the PBM is integrated into the DNO+ model, resulting in a physics-informed DNO+ (PIDNO+). The PIDNO+ model addresses the non-conservation of quality exhibited by the PBM while inheriting the advantages of the DNO+ model in considering multiple influencing factors. Moreover, the PIDNO+ model significantly reduces the amount of data required for model training. Both deep learning models, i.e., DNO+ and PIDNO+, are excellent in predictive performance, offering swift and accurate machine learning-based predictions. The use of this code that contains these models will assist in guiding the proper milling equipment selection and operational conditions to achieve the desired biomass particle sizes, ensuring the efficiency of subsequent biofuel and bioproduct production processes.

Xia, Yidong [Idaho National Laboratory (INL), Idah↗

DOE OSTI2015

Gyrokinetic studies of core turbulence features in ASDEX Upgrade H-mode plasmas

Gyrokinetic validation studies are crucial for developing confidence in the model incorporated in numerical simulations and thus improving their predictive capabilities. As one step in this direction, we simulate an ASDEX Upgrade discharge with the GENE code, and analyze various fluctuating quantities and compare them to experimental measurements. The approach taken is the following. First, linear simulations are performed in order to determine the turbulence regime. Second, the heat fluxes in nonlinear simulations are matched to experimental fluxes by varying the logarithmic ion temperature gradient within the expected experimental error bars. Finally, the dependence of various quantities with respect to the ion temperature gradient is analyzed in detail. It is found that density and temperature fluctuations can vary significantly with small changes in this parameter, thus making comparisons with experiments very sensitive to uncertainties in the experimental profiles. However, cross-phases are more robust, indicating that they are better observables for comparisons between gyrokinetic simulations and experimental measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

SU-D-BRD-02: Auto Weekly - An Automated Online Weekly Chart Check System for Medical Physics

Purpose: Radiation therapy has evolved to become not only more precise and potent, but also more complicated to monitor and deliver. More rigorous and comprehensive quality assurance is needed to safeguard ever advancing radiation therapy. ICRU standards dictate that an ever growing set of treatment parameters are manually checked weekly by medical physicists. This “weekly chart check” procedure is laborious and subject to human errors or other factors. A computer-assisted chart checking process will enable more complete and accurate human review of critical parameters, reduce the risk of medical errors, and improve the efficiency. Methods: We developed a web-based software system that enables a thorough weekly quality assurance checks. In the backend, the software retrieves all machine parameters from a Treatment Management System (TMS) and compares them against the corresponding ones from the treatment planning system. They are also checked for validity against preset rules. The results are displayed as a web page in the front-end for physicists to review. Then a summary report is generated and uploaded automatically to the TMS as a record for weekly chart checking. Results: The software system has been deployed on a web server in our department’s intranet, and has been tested thoroughly by our clinical physicists. A plan parameter would be highlighted when it is off the preset limit. The developed system has changed the way of checking charts with significantly improved accuracy, efficiency, and completeness. It has been shown to be robust, fast, and easy to use. Conclusion: A computer-assisted system has been developed for efficient, accurate, and comprehensive weekly chart checking. The system has been extensively validated and is being implemented for routine clinical use.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Anderson acceleration stability in NDA-accelerated k-eigenvalue problems

Anderson acceleration (AA) has been used to improve the stability and convergence rate of multiphysics iterative methods for reactor analysis. Most applications studied assume a tightly converged solution for the different physics problems, and AA is usually applied to state variables like temperature, density, and heat generation rate. In this paper, we study the theoretical performance of AA in NDA-accelerated k-eigenvalue problems. The problems and algorithms studied are simplified from the coupled iteration scheme adopted by MPACT and many other high-fidelity whole-core reactor codes. Compared to previous analyses of AA for these iteration schemes, we study the case with a partially converged neutronics solution and possibly partially converged nonlinear diffusion acceleration (NDA)/coarse mesh finite difference (CMFD) solutions. We observe that the performance of the iteration scheme with AA is very sensitive to the initial guess and is affected by the partially converged CMFD solutions. When the NDA solution is fully converged, using AA cannot achieve the optimal convergence rate in large-sized problems. Conversely, if the NDA solution is partially converged, the iteration scheme with AA can diverge or converge extremely slowly. It is found that the loss of robustness for AA is due to the fact that it is applied to the iterative subspace of state variables rather than the fundamental unknowns of the governing equations. To improve the robustness, the scalar flux should also be considered in the implementation of AA. After considering the residuals of flux, we observe that the stability is regardless of the partial convergence of NDA solutions. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Verification and validation of the open-source plasma fluid code: Zapdos

Zapdos is an open-source finite element plasma fluid solver based on the MOOSE multiphysics framework. This paper outlines Zapdos verification, benchmarking, and validation efforts for 1D and 2D RF capacitively coupled plasma discharge models for mid-range pressures (0.1 - 1 Torr). The verification process involved using the method of manufactured solutions to assess Zapdos spatial and temporal error convergence. L2 errors ranged from 10 -2 to 10 -4 , while the convergence’s slope were in agreement with the predicted slopes for the tested variable and time integration orders. The benchmarking process involved comparisons to previously results from the validated finite element code, LSODI. These works included 1D and 2D simulations for a range of plasma parameters (densities, temperatures, voltage, etc.). For the 1D cases, Zapdos and LSODI results were in very good agreement. In the 2D cases, variable behaviors matched, with slight discrepancies in peak values. The validation process involved comparisons to experimental works including electron density measurements by microwave interferometry and metastable density measurements by planar laser-induced fluorescence imaging. Results shown reasonable agreement at higher pressure, with results starting to diverge at low pressures. Furthermore, probable causes for this diverges are the limitation of the fluid assumption for plasmas at low pressure, or the need for more robust boundary conditions. Overall, Zapdos shown reasonable results for the verification, benchmarking, and validation efforts, and Zapdos can be downloaded at https://github.com/shannon-lab/zapdos.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

Numerical simulations of global Alfvén eigenmodes excitation and stabilization in NSTX-U

Three-dimensional nonlinear simulations of Alfvén eigenmodes in the subcyclotron frequency range demonstrate a robust physical stabilizing mechanism via modest off-axis beam injection, in agreement with experimental observations from the National Spherical Torus Experiment (NSTX-U). Experimental results from NSTX-U have demonstrated that neutral beam injection from the new beam sources with large tangency radii deposits beam ions with large pitch, which can very effectively stabilize all unstable Global Alfvén Eigenmodes (GAEs). Beam-driven GAEs have been linked to enhanced electron transport in NSTX, and the ability to control these modes will have significant implications for NSTX-U, ITER, and other fusion devices where super-Alfvénic fast ions might be present. Nonlinear simulations using the HYM code have been performed to study the excitation and stabilization of GAEs in the NSTX-U right before and shortly after the additional off-axis beam injection. The simulations reproduce the experimental finding, namely, it is shown that off-axis neutral beam injection reliably and strongly suppresses all unstable GAEs. Before additional beam injection, the simulations show unstable counter-rotating GAEs with toroidal mode numbers and frequencies that match the experimentally determined modes. Additional off-axis beam injection has been modeled by adding beam ions with large pitch and varying density. The complete stabilization occurs at less than 7% of the total beam ion inventory. New analytical theory of GAE (de)stabilization has also been derived, indicating a different interpretation for the GAE stabilization mechanism compared to previous publications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Energy Paths for the Yurok People

This Yurok Tribe energy strategic planning study examined electricity usage, costs and infrastructure, as well as electricity needs on the Reservation. It also assessed the Tribe’s energy resources. Equipped with this information, the project team conducted an energy options analysis that examined different project opportunities and identified key projects that are both feasible to implement and will help the Tribe reach its stated Energy Vision. A key finding of the energy needs assessment is that there are still numerous households that are without grid electric power, and it is likely that most of them will remain disconnected from the grid because of their remoteness. To provide these community members with reliable, affordable, modern energy services their homes will need to be electrified via a robust off-grid energy services program that installs and maintains off-grid energy systems. This is the key need that is explored and addressed in this plan. The approach taken is to establish a new business entity, a Tribal Energy Authority or Utility, that will operate a Tribal Energy Program. The Energy Strategic Plan identifies a number of potential funding sources that could be used to support both the On-grid and Off-grid Energy Programs. This plan also promotes continued efforts to reduce electricity consumption through the implementation of energy efficiency measures in existing facilities and the enforcement of a strong energy-efficient building code for all new construction on the Reservation. Finally, the plan promotes added resilience at critical Tribal facilities on the Reservation by recommending that battery energy storage systems and microgrid controls be coupled with rooftop solar PV systems to create cutting edge microgrid facilities that can provide critical power services to Tribal facilities when the main electric grid is shutdown.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2016

M DWARF FLARE CONTINUUM VARIATIONS ON ONE-SECOND TIMESCALES: CALIBRATING AND MODELING OF ULTRACAM FLARE COLOR INDICES

We present a large data set of high-cadence dMe flare light curves obtained with custom continuum filters on the triple-beam, high-speed camera system ULTRACAM. The measurements provide constraints for models of the near-ultraviolet (NUV) and optical continuum spectral evolution on timescales of ≈1 s. We provide a robust interpretation of the flare emission in the ULTRACAM filters using simultaneously obtained low-resolution spectra during two moderate-sized flares in the dM4.5e star YZ CMi. By avoiding the spectral complexity within the broadband Johnson filters, the ULTRACAM filters are shown to characterize bona fide continuum emission in the NUV, blue, and red wavelength regimes. The NUV/blue flux ratio in flares is equivalent to a Balmer jump ratio, and the blue/red flux ratio provides an estimate for the color temperature of the optical continuum emission. We present a new “color–color” relationship for these continuum flux ratios at the peaks of the flares. Using the RADYN and RH codes, we interpret the ULTRACAM filter emission using the dominant emission processes from a radiative-hydrodynamic flare model with a high nonthermal electron beam flux, which explains a hot, T ≈ 10{sup 4} K, color temperature at blue-to-red optical wavelengths and a small Balmer jump ratio as observed in moderate-sized and large flares alike. We also discuss the high time resolution, high signal-to-noise continuum color variations observed in YZ CMi during a giant flare, which increased the NUV flux from this star by over a factor of 100.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

TH-CD-202-12: Online Inter-Beam Replanning Based On Real-Time Dose Reconstruction

Purpose: This work provides a proof-of-concept study for online replanning during treatment delivery for step-and-shoot prostate SBRT, based on real-time dose reconstruction. Online replanning is expected to improve the trade-off between target coverage and organ-at-risk dose in the presence of intra-fractional motion. Methods: We have implemented an online replanning workflow on top of our previously reported real-time dose reconstruction software which connects to an Elekta research linac. The treatment planning system DynaPlan was extended to (1) re-optimize and sequence treatment plans (in clockwise beam order) before each beam, based on actual delivered dose, in a timeframe limited by the gantry rotation between subsequent beams, and (2) send the respective segments to the delivery control software DynaTrack which starts/continues treatment immediately.To investigate the impact of a reduced safety margin, we have created and delivered (on a linac emulator) a conventional CTV+5/3mm (I) and a reduced CTV+1mm margin (II) treatment plan for a prostate patient. We have assessed CTV coverage with and without inter-beam replanning, all exposed to a gradual target shift of 0–5mm in posterior and inferior direction from start until the end of delivery. Results: For the reconstructed conventional plan (I), D98 for CTV was 100% of D98 of the planned dose. For the reconstructed margin-reduced plan (II), D98 for CTV was 95% of the planned D98 without replanning, but could be recovered to 99% by replanning for each beam. Plan (II) with replanning resulted in a decrease for bladder V90% by 88% and an increase to rectum V90% by 9% compared to the conventional plan (I). Dose calculation/accumulation was performed in <15ms per MLC aperture, replanning in <15s per beam. Conclusion: We have shown that online inter-beam replanning is technically feasible and potentially allows for a margin reduction. Future investigation considering motion-robust replanning optimization parameters is in progress. We acknowledge support of the MLC research from Elekta AB. This work is supported by Cancer Research UK under Programme C33589/A19908. Research at ICR is also supported by Cancer Research UK under Programme C33589/A19727 and NHS funding to the NIHR Biomedical Research Centre at RMH and ICR.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2019

Investigation of light ion fusion reactions with plasma discharges

The scaling of reaction yields in light ion fusion to low reaction energies is important for our understanding of stellar fuel chains and the development of future energy technologies. Experiments become progressively more challenging at lower reaction energies due to the exponential drop of fusion cross sections below the Coulomb barrier. We report on experiments where deuterium-deuterium (D-D) fusion reactions are studied in a pulsed plasma in the glow discharge regime using a benchtop apparatus. We model plasma conditions using particle-in-cell codes. Advantages of this approach are relatively high peak ion currents and current densities (0.1 to several A/cm 2 ) that can be applied to metal wire cathodes for several days. We detect neutrons from D-D reactions with scintillator-based detectors. For palladium targets, we find neutron yields as a function of cathode voltage that are over 100 times higher than yields expected for bare nuclei fusion at ion energies below 2 keV (center of mass frame). A possible explanation is a correction to the ion energy due to an electron screening potential of 1000 ± 250 eV, which increases the probability for tunneling through the repulsive Coulomb barrier. Finally, our compact, robust setup enables parametric studies of this effect at relatively low reaction energies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Shell-model studies of the astrophysical r p -process reactions S 34 ( p , γ ) Cl 35 and Cl 34 g , m ( p , γ ) Ar 35

Background: Dust grains condensed in the outflows of presolar classical novae should have been present in the protosolar nebula. Candidates for such presolar nova grains have been found in primitive meteorites and can in principle be identified by their isotopic ratios, but the ratios predicted by state-of-the-art one-dimensional hydrodynamic models are uncertain due to nuclear-physics uncertainties. Purpose: To theoretically calculate the thermonuclear rates and uncertainties of the S 34 ( p , γ ) Cl 35 and Cl 34 g , m ( p , γ ) Ar 35 reactions and investigate their impacts on the predicted 34 S/ 32 S isotopic ratio for presolar nova grains. Method: A shell-model approach in a (0 +1 ) ℏ ω model space was used to calculate the properties of resonances in the 34 S(p ,γ ) 35 Cl and 34 g,m Cl (p ,γ ) 35 Ar reactions and their thermonuclear rates. Uncertainties were estimated using a Monte Carlo method. The implications of these rates and their uncertainties on sulfur isotopic nova yields were investigated using a postprocessing nucleosynthesis code. The rates for transitions from the ground state of 34 Cl as well as from the isomeric first excited state of 34 Cl were explicitly calculated. Results: At energies in the resonance region near the proton-emission threshold, many negative-parity states appear. Energies, spectroscopic factors, and proton-decay widths are reported. The resulting thermonuclear rates are compared with previous determinations. Conclusions: The shell-model calculations alone are sufficient to constrain the variation of the 34 S/ 32 S ratios to within about 30%. Uncertainties associated with other reactions must also be considered, but in general we find that the 34 S/ 32 S ratios are not a robust diagnostic to clearly identify presolar grains made from nova ejecta.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Fast Solution of Fully Implicit Runge--Kutta and Discontinuous Galerkin in Time for Numerical PDEs, Part II: Nonlinearities and DAEs

Fully implicit Runge--Kutta (IRK) methods have many desirable accuracy and stability properties as time integration schemes, but high-order IRK methods are not commonly used in practice with large-scale numerical PDEs because of the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic framework for solving the nonlinear equations that arise from IRK methods (and discontinuous Galerkin discretizations in time) applied to nonlinear numerical PDEs, including PDEs with algebraic constraints. Several new linearizations of the nonlinear IRK equations are developed, offering faster and more robust convergence than the often-considered simplified Newton, as well as an effective preconditioner for the true Jacobian if exact Newton iterations are desired. Inverting these linearizations requires solving a set of block 2 x 2 systems. Under quite general assumptions, it is proven that the preconditioned 2 x 2 operator's condition number is bounded by a small constant close to one, independent of the spatial discretization, spatial mesh, and time step, and with only weak dependence on the number of stages or integration accuracy. Moreover, the new method is built using the same preconditioners needed for backward Euler-type time stepping schemes, so can be readily added to existing codes. The new methods are applied to several challenging fluid flow problems, including the compressible Euler and Navier--Stokes equations, and the vorticity-streamfunction formulation of the incompressible Euler and Navier--Stokes equations. Up to 10th-order accuracy is demonstrated using Gauss IRK, while in all cases fourth-order Gauss IRK requires roughly half the number of preconditioner applications as required by standard Singly diagonally implicit Runge--Kutta methods.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

A Novel Workflow for Coupled Simulation of Hydraulic Stimulation with Simultaneous Injection of Proppant

Enhanced Geothermal Systems (EGS) are located in impermeable host rock. Therefore, to achieve promising circulation, hydraulic stimulation is a key factor. Stimulation involves injecting a combination of fluids and proppant particles into existing fractures under high pressure. The effectiveness of stimulation depends on a variety of different factors, including the pressurization scheme and the optimal choice of carrying fluid and proppant particles to be used. For example, high-density particles can settle out of the injected slurry quickly, whereas lighter particles as well as a dense carrying fluid allow for a long transport into the reservoir. A numerical tool is therefore essential to take all these factors into account and to make robust predictions about the performance of the hydraulic stimulation. We present a novel computational application for coupled numerical simulation of hydraulic stimulation with concurrent injection of proppant in a fractured geothermal reservoir. This application is based on the open-source MOOSE framework that has been enriched with a set of equations to account for the circulation of proppant-laden slurries. It accounts for both the fluid and proppant components of the mixture, as well as particle-driven processes that include gravity settling, particle-particle interactions, and strong density and viscosity contrasts. A particular challenge is the mutual coupling with fracture mechanical processes in a propagating fracture using the FE code Defmod. Aperture and fluid pressure on one side and slurry density and viscosity changes on the other side interact and affect each other. Their basic understanding is therefore essential for the success of the stimulation. We compare our application to common numerical benchmark problems for both slurry flow and hydraulic stimulation.

15 GEOTHERMAL ENERGY↗

DOE OSTI2022

Layered CAD/CSG geometry for spatially complex radiation transport scenarios

Many spatially complex fission, fusion, and national security Monte Carlo (MC) radiation transport scenarios involve combining computer-aided design (CAD) models with constructive solid geometry (CSG) models. A layered geometry method has been implemented in the Shift MC code to address this need. With layered geometry, multiple CAD and/or CSG models can be clipped, translated, rotated, and placed in overlapping layers to form transport-ready geometries. Here, the utility of this method is demonstrated with two problems: (1) a fixed-source simulation with a layered geometry consisting of a LiDAR-generated CAD model of the Combined Arms Collective Training Facility urban environment overlaid with CSG models of a mock hotel and a detector apparatus, and (2) a k-eigenvalue calculation using a layered geometry model of the Transformational Challenge Reactor consisting of CAD fuel elements placed in a CSG core. Tallied particle flux distributions match expectations, but tracking robustness must be improved prior to general-purpose use.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗