Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “production simulation codes”

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,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Importance of 3-D S{sub N} depletion in non-proliferation using BSOLVE

We present a single Pressurized Water Reactor (PWR) 3-D fuel rod design for depletion analysis using BSOLVE, our newly developed Runge-Kutta-Fehlberg based depletion code. BSOLVE is coupled with the deterministic 3-D S{sub N} particle transport code, PENTRAN, applied here with a 4-neutron energy group comparison to Continuous Energy (C/E) SERPENT2 Monte Carlo results. Differences are expected, as PENTRAN+BSOLVE retains full (multi-group) energy information for reactions, nuclide specific fission contributions, and energy dependent fission yields, using the latest available ENDF-BVIII data, important to retain accurate burned fuel inventories; SERPENT2 collapses burnup reactions to a single energy value. For depletion times up to ∼ 700 days and typical PWR power densities, relative differences between multigroup 3-D S{sub N} with full energy data and Monte Carlo one group burnup for trans-uranium nuclide concentrations and fission products are up to ∼20%. System eigenvalues are consistent, but with differences early and late in the cycle attributed to multigroup vs. C/E Monte Carlo cross sections. This work highlights the importance of low variance transport driven burnup for non-proliferation concerns, since plutonium quality varies significantly along axial lengths, and is more challenging to converge using Monte Carlo; details of depletion steps with spatial/zone dependent plutonium quality are provided. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

Processing of Ice Cloud In Situ Data Collected by Bulk Water, Scattering, and Imaging Probes: Fundamentals, Uncertainties, and Efforts toward Consistency

Abstract In situ observations of cloud properties made by airborne probes play a critical role in ice cloud research through their role in process studies, parameterization development, and evaluation of simulations and remote sensing retrievals. To determine how cloud properties vary with environmental conditions, in situ data collected during different field projects processed by different groups must be used. However, because of the diverse algorithms and codes that are used to process measurements, it can be challenging to compare the results. Therefore it is vital to understand both the limitations of specific probes and uncertainties introduced by processing algorithms. Since there is currently no universally accepted framework regarding how in situ measurements should be processed, there is a need for a general reference that describes the most commonly applied algorithms along with their strengths and weaknesses. Methods used to process data from bulk water probes, single-particle light-scattering spectrometers and cloud-imaging probes are reviewed herein, with emphasis on measurements of the ice phase. Particular attention is paid to how uncertainties, caveats, and assumptions in processing algorithms affect derived products since there is currently no consensus on the optimal way of analyzing data. Recommendations for improving the analysis and interpretation of in situ data include the following: establishment of a common reference library of individual processing algorithms, better documentation of assumptions used in these algorithms, development and maintenance of sustainable community software for processing in situ observations, and more studies that compare different algorithms with the same benchmark datasets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Particle-in-cell modeling for MJ scale dense plasma focus with varied anode shape

Megajoule scale dense plasma focus (DPF) Z-pinches with deuterium gas fill are compact devices capable of producing 10{sup 12} neutrons per shot but past predictive models of large-scale DPF have not included kinetic effects such as ion beam formation or anomalous resistivity. We report on progress of developing a predictive DPF model by extending our 2D axisymmetric collisional kinetic particle-in-cell (PIC) simulations from the 4 kJ, 200 kA LLNL DPF to 1 MJ, 2 MA Gemini DPF using the PIC code LSP. These new simulations incorporate electrodes, an external pulsed-power driver circuit, and model the plasma from insulator lift-off through the pinch phase. To accommodate the vast range of relevant spatial and temporal scales involved in the Gemini DPF within the available computational resources, the simulations were performed using a new hybrid fluid-to-kinetic model. This new approach allows single simulations to begin in an electron/ion fluid mode from insulator lift-off through the 5-6 μs run-down of the 50+ cm anode, then transition to a fully kinetic PIC description during the run-in phase, when the current sheath is 2-3 mm from the central axis of the anode. Simulations are advanced through the final pinch phase using an adaptive variable time-step to capture the fs and sub-mm scales of the kinetic instabilities involved in the ion beam formation and neutron production. Validation assessments are being performed using a variety of different anode shapes, comparing against experimental measurements of neutron yield, neutron anisotropy and ion beam production.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Hydrodynamic Simulations of Oxygen–Neon Classical Novae as Galactic 7 Li Producers and Potential Accretion-induced Collapse Progenitors*

Abstract We report on studies of classical nova (CN) explosions where we follow the evolution of thermonuclear runaways (TNRs) on oxygen–neon (ONe) white dwarfs (WDs). Using NOVA, a 1D hydrodynamic computer code, we accrete solar matter until the TNR is ongoing and then switch to a mixed composition. This approach is guided by the results of multidimensional studies of TNRs in WDs, which find that sufficient mixing with WD core material occurs after the TNR is well underway, and levels of enrichment of the CNONeMg elements are reached that agree with observations of CN ejecta abundances. Because the amount of accreted material is inversely proportional to the oxygen abundance, by first accreting solar matter, the amount of accreted material is larger than in those simulations with an initially enriched composition. We vary the mass of the WD (from 0.6 M ⊙ to 1.35 M ⊙ ) and the composition of the mixed materials. Our results show large enrichments of 7 Be in the ejected gases, implying that ONe CNe and CO CNe may be responsible for a significant fraction (∼100 M ⊙ ) of the Galactic 7 Li (∼1000 M ⊙ ). The production of 22 Na and 26 Al in CN explosions and the γ -ray emission predicted by our simulations are discussed. The WDs in all our simulations eject less material than they accrete and we predict that the WD is growing in mass as a consequence of the CN outburst. ONe CNe, therefore, may be an important channel for accretion-induced collapse events.

Starrfield, Sumner (ORCID:0000000213596312)↗
DOE OSTI2020

Failure investigation of nuclear grade POCO graphite target in high energy neutrino physics through numerical simulation

In many high energy physics neutrino beamlines, targets are made up of isotropic POCO graphite grade for production of neutrinos for high energy particle physics research and are bombarded with highly energetic pulsed proton beam. The pulsed proton beam with a small beam size creates thermal stress waves as well as radiation damage such as displacement damage, void formation, swelling, and gas formation amongst others. As a result, after a few years of operation one of the longest operated targets appears to have undergone bulk swelling and fractured along the beam centerline. A complex interaction of dynamic loading due to beam, material degradation due to irradiation are thought to have caused such a failure. Here, we present a numerical simulation to capture the combined effects of swelling and the dynamic effect of thermal stress waves due to beam heating, to explain the failure of target. An empirical formula has been developed to take into account swelling as a function of temperature and proton fluence and has been implemented in a commercial finite element code. Extensive X-ray diffraction (XRD) were performed on graphite samples from the failed target to understand the lattice parameter changes due to irradiation and build a valid empirical model. Simulation results show that swelling plays a major role in elevating the stress state in the material exceeding the failure strength while the pulsed beam heating also introduces fatigue loading that would explain brittle fracture emanating from inside of the material and propagating outward. Such empirical formulations will help in identifying the threshold values of critical parameters to extend the service life of future multi-megawatt neutrino targets and targets in other accelerator environments.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Achieving Higher Order Accuracy in Space in Hydrodynamic Simulations of Self-Gravitating Gas

Modern astrophysical simulation codes employ a variety of numerical algorithms capable of achieving higher-order accuracy in both space and time. Albeit they succeed in achieving an effective higher spatial resolution and in suppressing the numerical damping of waves, to our knowledge, all current astrophysical simulations invoking self-gravity are limited to second-order accuracy in space. If we can devise an algorithm to evaluate self-gravity with a higher-order spatial accuracy, we can better the evaluation of the gravitational acceleration and gravitational energy release which dictate the evolution of many astrophysical systems. Herein, we present a numerical algorithm for self-gravitating hydrodynamics capable of achieving fourth order accuracy for a given density distribution on a Cartesian uniform grid. First, we derive the cell-averaged gravitational potential at fourth-order accuracy from the cell-averaged density by solving the Poisson equation. Next, we obtain the cell average of the product of the density and gravitational acceleration, which differs from the cell-averaged density multiplied by the cell-averaged gravitational acceleration. We then show the verification of the algorithm by applying it to critical test problems: (1) maintaining equilibria of self-gravitating slabs, even upon advection, (2) evolving a polytropic sphere with a massive power-law envelope, and (3) conservation of specific entropy during the propagation of a sound wave.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

CFD flowfield simulation of Delta Launch Vehicles in a power-on configuration

This paper summarizes recent work at McDonnell Douglas Aerospace (MDA) to develop and validate computational fluid dynamic (CFD) simulations of under expanded rocket plume external flowfields for multibody expendable launch vehicles (ELVs). Multi engine reacting gas flowfield predictions of ELV base pressures are needed to define vehicle base drag and base heating rates for sizing external nozzle and base region insulation thicknesses. Previous ELV design programs used expensive multibody power-on wind tunnel tests that employed chamber/nozzle injected high pressure cold or hot-air. Base heating and pressure measurements were belatedly made during the first flights of past ELV's to correct estimates from semi-empirical engineering models or scale model tests. Presently, CFD methods for use in ELV design are being jointly developed at the Space Transportation Division (MDA-STD) and New Aircraft Missiles Division (MDA-NAMD). An explicit three dimensional, zonal, finite-volume, full Navier-Stokes (FNS) solver with finite rate hydrocarbon/air and aluminum combustion kinetics was developed to accurately compute ELV power-on flowfields. Mississippi State University's GENIE++ general purpose interactive grid generation code was chosen to create zonal, finite volume viscous grids. Axisymmetric, time dependent, turbulent CFD simulations of a Delta DSV-2A vehicle with a MB-3 liquid main engine burning RJ-1/LOX were first completed. Hydrocarbon chemical kinetics and a k-epsilon turbulence model were employed and predictions were validated with flight measurements of base pressure and temperature. Zonal internal/external grids were created for a Delta DSV-2C vehicle with a MB-3 and three Castor-1 solid motors burning and a Delta-2 with an RS-27 main engine (LOX/RP-1) and 9 GEM's attached/6 burning. Cold air, time dependent FNS calculations were performed for DSV-2C during 1992. Single phase simulations that employ finite rate hydrocarbon and aluminum (solid fuel) combustion chemistry are currently in progress. Reliable and efficient Eulerian algorithms are needed to model two phase (solid-gas) momentum and energy transfer mechanisms for solid motor fuel combustion products.

Pavish, D. L.↗
DOE OSTI2024

Achieving Higher Order Accuracy in Space in Hydrodynamic Simulations of Self-Gravitating Gas

Modern astrophysical simulation codes employ a variety of numerical algorithms capable of achieving higher-order accuracy in both space and time. Albeit they succeed in achieving an effective higher spatial resolution and in suppressing the numerical damping of waves, to our knowledge, all current astrophysical simulations invoking self-gravity are limited to second-order accuracy in space. If we can devise an algorithm to evaluate self-gravity with a higher-order spatial accuracy, we can better the evaluation of the gravitational acceleration and gravitational energy release which dictate the evolution of many astrophysical systems. Herein, we present a numerical algorithm for self-gravitating hydrodynamics capable of achieving fourth-order accuracy for a given density distribution on a Cartesian uniform grid. First, we derive the cell-averaged gravitational potential at fourth-order accuracy from the cell-averaged density by solving the Poisson equation. Next, we obtain the cell average of the product of the density and gravitational acceleration, which differs from the cell-averaged density multiplied by the cell-averaged gravitational acceleration. We then show the verification of the algorithm by applying it to critical test problems: (1) maintaining equilibria of self-gravitating slabs, even upon advection, (2) evolving a polytropic sphere with a massive power-law envelope, and (3) conservation of specific entropy during the propagation of a sound wave.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Advances in ArborX to support exascale applications

ArborX is a performance portable geometric search library developed as part of the Exascale Computing Project (ECP). In this paper, we explore a collaboration between ArborX and a cosmological simulation code HACC. Large cosmological simulations on exascale platforms encounter a bottleneck due to the in-situ analysis requirements of halo finding, a problem of identifying dense clusters of dark matter (halos). This problem is solved by using a density-based DBSCAN clustering algorithm. With each MPI rank handling hundreds of millions of particles, it is imperative for the DBSCAN implementation to be efficient. In addition, the requirement to support exascale supercomputers from different vendors necessitates performance portability of the algorithm. We describe how this challenge problem guided ArborX development, and enhanced the performance and the scope of the library. We explore the improvements in the basic algorithms for the underlying search index to improve the performance, and describe several implementations of DBSCAN in ArborX. Further, we report the history of the changes in ArborX and their effect on the time to solve a representative benchmark problem, as well as demonstrate the real world impact on production end-to-end cosmology simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Criticality Evaluation of Low Density Product Oxides for Air Transport for the ES-3100 Package

This document presents the criticality safety evaluation of the Y-12 National Security Complex (Y-12) Model ES-3100 package for the addition of lower bulk density limits for HEU product oxides transported by air. The current authorized oxide density range for air transport is 2.0 to 6.54 g/cm 3 . The analysis in this report looks at lower densities, specifically in the range of 1.0 to 2.0 g/cm 3 . For reasons presented in this report, it was decided to present the entire range of oxide densities (1.0 to 6.54 g/cm 3 ) as a series of sub-ranges. The lowest density in each range was evaluated. Thus calculations were performed for the following oxide densities: 5.0 g/cm 3 , 4.0 g/cm 3 , 3.0 g/cm 3 , 2.0 g/cm 3 , 1.5 g/cm 3 and 1.0 g/cm 3 . Calculations presented herein were performed using a third generation computing platform; the Microway Parallel Linux Cluster ENGHPC workstation using the Monte Carlo Criticality Safety Analysis Sequence (CSAS5) criticality code sequence and 238-group ENDF/B-VII cross-section library of Standardized Computer Analysis for Licensing Evaluation (SCALE) 6.1.3 (ORNL/TM-2005/39, SCALE: A Comprehensive Modeling and Simulation Suite for Nuclear Safety Analysis and Design).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Heating effects on jack pine pyrogenic organic matter properties from a pyrocosm study in 2022

This dataset contains data associated with the preprint “Fire removes preexisting pyrogenic organic matter from the ecosystem through the mechanisms of both direct combustion and increasing mineralizability” (Luo et al., 2025b), which is the complementary study to the published paper “Reburning pyrogenic organic matter: a laboratory method for dosing dynamic heat fluxes from above” (Luo et al., 2025a). We designed a full-factorial experiment with different burial depths of jack pine (Pinus banksiana Lamb) pyrogenic organic matter (PyOM) (Surface, 1 cm, and 5 cm) and different heat-flux profiles (High, Low, and Control) to examine how subsequent fires affect the properties of preexisting PyOM. We measured total carbon (C), pH, dissolved organic carbon (DOC), dissolved inorganic carbon (DIC), and mineralized C (as CO₂-C, from a 12-week incubation).We found that high heat flux and/or surface placement resulted in substantial direct C losses through combustion. Intermediate heat exposure produced both combustion losses and increases in DOC and mineralizability, which may have complex long-term implications: an increased dissolved fraction of PyOM may promote downward transport into mineral soils and potentially contribute to deeper, longer-term C storage, but it may also make PyOM more susceptible to microbial decomposition. Under the lowest heat flux and deepest burial, most PyOM was retained, and changes in DOC and C mineralization were minimal. Finally, PyOM pH, an important chemical property, decreased under low-temperature heating but increased under higher temperatures.We uploaded pH data for all samples (“pH_of_all_samples.csv”); pH and temperature-related data (peak temperature and degree hours) for samples in High and Low heat-flux treatments (“pH_vs_peakT_and_degree_hours_only_for_heated_samples.csv”); total C data (“CN_pct_C_stock_C_loss_in_samples.csv”); DOC and DIC data (“doc_dic.csv”); and mineralized C (CO₂-C) data (“CO2-C_all_original.csv”). Additional details can be found in the Methods & Sampling section.All datasets uploaded to ESS-DIVE are clearly labeled, cleaned, and include both raw and derived data, ready for reuse in other analyses. All analysis code and raw datasets are also available on GitHub: https://github.com/MengmengLuo/Fire-removes-preexisting-pyrogenic-organic-matter-from-the-ecosystem.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Impact of increased smoothing by spectral dispersion bandwidth on stimulated Brillouin scattering in laser driven Hohlraums

Experiments were conducted at the National Ignition Facility (NIF) to investigate the impact of increased smoothing by spectral dispersion (SSD) bandwidth on the production of stimulated Brillouin scattering (SBS) within an indirect-drive inertial confinement fusion (ICF) Hohlraum. This was done in a subscale gold Hohlraum driven by 192 laser beams depositing 1.1 MJ of energy. The laser bandwidth was increased from 45 to 118 GHz (before frequency tripling) on the 30° and 50° cones, where backscatter could be measured using the Full Aperture Backscatter Station (FABS). It was expected that this 2.6-fold increase in bandwidth would suppress SBS generated within the Hohlraum plasma and lower the backscattered SBS energy on the 50° cones by over a factor of four. Experimental results, however, show that this SSD change only reduced the 50° cone SBS during the main capsule drive by -18±31% and -4.5±7.8% over the entire pulse. This is small compared to expected shot-to-shot SBS reproducibility (∼30%), such that the result can be considered within normal performance fluctuations. New 3D parallel paraxial code (pF3D) simulations, accounting for beam refractive intensification reproduce this result, suggesting that closer to 300 GHz of bandwidth would have been required to mitigate SBS to the expected level. Delivering such a high bandwidth is not feasible at NIF when operating at high peak power and would potentially prevent NIF’s ability to use cross-beam energy transfer (CBET) for implosion symmetry tuning.

Physics↗
NASA NTRS1995

NAS Experiences of Porting CM Fortran Codes to HPF on IBM SP2 and SGI Power Challenge

Current Connection Machine (CM) Fortran codes developed for the CM-2 and the CM-5 represent an important class of parallel applications. Several users have employed CM Fortran codes in production mode on the CM-2 and the CM-5 for the last five to six years, constituting a heavy investment in terms of cost and time. With Thinking Machines Corporation's decision to withdraw from the hardware business and with the decommissioning of many CM-2 and CM-5 machines, the best way to protect the substantial investment in CM Fortran codes is to port the codes to High Performance Fortran (HPF) on highly parallel systems. HPF is very similar to CM Fortran and thus represents a natural transition. Conversion issues involved in porting CM Fortran codes on the CM-5 to HPF are presented. In particular, the differences between data distribution directives and the CM Fortran Utility Routines Library, as well as the equivalent functionality in the HPF Library are discussed. Several CM Fortran codes (Cannon algorithm for matrix-matrix multiplication, Linear solver Ax=b, 1-D convolution for 2-D datasets, Laplace's Equation solver, and Direct Simulation Monte Carlo (DSMC) codes have been ported to Subset HPF on the IBM SP2 and the SGI Power Challenge. Speedup ratios versus number of processors for the Linear solver and DSMC code are presented.

Saini, Subhash↗
DOE OSTI2016

DEVELOPMENT AND TESTING OF FAULT-DIAGNOSIS ALGORITHMS FOR REACTOR PLANT SYSTEMS

Argonne National Laboratory is further developing fault diagnosis algorithms for use by the operator of a nuclear plant to aid in improved monitoring of overall plant condition and performance. The objective is better management of plant upsets through more timely, informed decisions on control actions with the ultimate goal of improved plant safety, production, and cost management. Integration of these algorithms with visual aids for operators is taking place through a collaboration under the concept of an operator advisory system. This is a software entity whose purpose is to manage and distill the enormous amount of information an operator must process to understand the plant state, particularly in off-normal situations, and how the state trajectory will unfold in time. The fault diagnosis algorithms were exhaustively tested using computer simulations of twenty different faults introduced into the chemical and volume control system (CVCS) of a pressurized water reactor (PWR). The algorithms are unique in that each new application to a facility requires providing only the piping and instrumentation diagram (PID) and no other plant-specific information; a subject-matter expert is not needed to install and maintain each instance of an application. The testing approach followed accepted procedures for verifying and validating software. It was shown that the code satisfies its functional requirement which is to accept sensor information, identify process variable trends based on this sensor information, and then to return an accurate diagnosis based on chains of rules related to these trends. The validation and verification exercise made use of GPASS, a one-dimensional systems code, for simulating CVCS operation. Plant components were failed and the code generated the resulting plant response. Parametric studies with respect to the severity of the fault, the richness of the plant sensor set, and the accuracy of sensors were performed as part of the validation exercise. The background and overview of the software will be presented to give an overview of the approach. Following, the verification and validation effort using the GPASS code for simulation of plant transients including a sensitivity study on important parameters will be presented

Fault Diagnosis Algorithm↗
DOE OSTI2011

Status of NIO1 construction

The NIO1 (Negative Ion Optimization phase 1) project consists of a multiaperture negative ion source mounted on a 60 kV accelerating column; up to 9 beamlets (15 mA H{sup -} each one) arranged in 3x3 matrix with 14 mm spacing can be extracted. The moderate size and the modular concept make some relative rotation of the source magnetic filter and the electrodes possible, so that the effect of crossed and aligned field can be easily compared. Other goals of source experimental program are emittance (and beam profile) measurement at several distances (for simulation code validation), testing of diagnostic components and of radiofrequency coupling. A full set of construction drawing was completed; also the Fast Emittance Scanner (FES) and its vacuum chambers were built (four mounting positions are reserved to FES). Some low power rf matching boxes were developed for a test plasma, approximately half the source size. A cesium oven compatible with NIO1 is being also developed; by using some industrial standard 100 W heaters (with a proper driver) a careful control of temperature is planned.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

PPPL Report on Reduced Modeling of Fusion Alpha Transport in ARC Burning Plasmas

We are reporting on the modeling of fusion alpha particle transport in the planned ARC fusion device being designed by the CFS (Commonwealth Fusion Systems: https://cfs.energy). The ARC tokamak is designed to operate in a burning-plasma regime characterized by a substantial population of fusion-born alpha particles. Alfvén eigenmode (AE) stability is assessed both analytically and numerically, incorporating alpha-particle drive, ion Landau, and radiative damping from thermal species and collisional damping from trapped electrons. Regions of unstable and near-threshold AE activity are mapped across ARC’s operational parameter space. Linear stability analysis with NOVA indicates multiple, often marginally unstable AEs, extending to toroidal mode numbers up to n= 30. The present report focuses on the ARC flat-top operating point prior to the sawtooth event. Alpha-particle transport on timescales exceeding the neoclassical slowing-down time is assessed using the NUBEAM module [1][2] of the TRANSP code [3], employing transport coefficients derived from the RBQ quasilinear modeling (cf. Appendix B). These global simulations identify favorable and unfavorable operating regimes with respect to alpha confinement, pressure redistribution, and overall alpha-heating efficiency. We also evaluate additional transport mechanisms—including neoclassical tearing mode (TM)–induced stochasticity, sawtooth-driven redistribution, and toroidal-field ripple using the kick model (cf. Appendix C) which makes use of the guiding-center code ORBIT, see Section 5. The kick model is integrated into TRANSP to enable self-consistent predictions of alpha-driven current formation and sustainment within the ARC scenario. Sensitivity scans are performed over the mode frequency, rational-surface alignment, island width, mode amplitude, and proximity of the limiter to the plasma. Our study provides an initial, physics-based guidance for machine design, operational planning, and equilibrium control, ensuring adequate alpha confinement and robust self-heating performance in ARC. Our simulations mostly targeted worst case scenarios, e.g. for TMs and sawteeth. Overall, we expect benign effects for the ARC scenario investigated in this work on fusion alpha confinement and losses in the presence of AEs, tearing modes and sawteeth. This report addresses three thrusts identified at the outset. The first thrust focuses on analytic estimates of the parametric dependencies of EP relaxation based on local AE stability simulations (Section 3). The second thrust involves global evaluations of AE stability using the NOVA, RBQ, and NUBEAM codes (Section 4). Finally, we investigate alpha-particle transport driven by low-frequency instabilities associated with sawteeth and tearing modes (Section 5).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Multiple Degradation Mechanisms in Reinforced Concrete Structures, Modeling and Risk Analysis

The overarching goal of this project is to complement ongoing Department of Energy (DOE) Light Water Reactor Sustainability (LWRS)-funded Grizzly concrete modeling development efforts by providing improved multi-physics models for incorporation into the Grizzly and BlackBear codes. Objectives of this work identified at the outset of this project include: 1. Coupling between mechanical damage and transport processes. The mechanical damage will be characterized by nonlinear mechanical constitutive models appropriate for concrete used in nuclear power plants (NPPs). 2. Improved representation of the coupling among the transport process, building on the current moisture and thermal transport models in Grizzly. The parameters for coupling among heat and mass transport in the multi-physics model will be experimentally determined. 3. Coupling among different scale levels of concrete constituents, fine and coarse aggregates, cement paste, hydration products, and pore structure. 4. Probabilistic analysis of the random nature of the heterogeneous concrete structure and the environmental factors (ambient temperature and humidity) and their multiple effects on concrete deterioration. 5. Benchmark verification, validation and uncertainty quantification of the thermo-hygro- chemo-mechanical (THCM) concrete formulation. The end goal of this work was to provide a comprehensive and robust simulation capability that can be applied at the engineering scale for analysis of realistic deterioration scenarios in NPP concrete structures.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Comparative modelling of the coupled thermal–hydraulic-mechanical (THM) processes in a heated bentonite pellet column with hydration

For the deep geological disposal of high-level radioactive waste in argillaceous rocks, the heat production of the waste is an important driver for thermal–hydraulic-mechanical (THM)-coupled processes. These THM processes influence the properties and conditions of the near field that in many repository designs contains bentonite as a clay buffer. One task in the DECOVALEX-2015 (DEvelopment of COupled models and their VALidation against Experiments) project was the modelling of a heated bentonite column (Villar et al. in Long-term THM tests reports: THM cells for the HE-E test: update of results until February 2014. Deliverable-no: D2.2-7.3. CIEMAT Technical Report IEMAT/DMA/2G210/03/2014, 2014) in preparation for the in situ heater experiment HE-E at the underground rock laboratory Mont Terri. DECOVALEX is an international cooperative project that focuses on the development and validation of mathematical models for simulating such coupled processes associated with disposal in deep geological repositories. Eight modelling teams developed their own THM-coupled models for the bentonite column experiment, using six different simulation codes. Each of the teams individually calibrated the THM parameters for the bentonite material. The eight resulting parameter sets agree well and allow a satisfactory reproduction of the TH measurements by all models. The modelling results for the evolution of temperature and relative humidity over time at three sensors in the bentonite column are in good agreement between the teams and with the measured data. Also, changes of the temperature due to modifications of the insulation and the adjustment of the heating power during the course of the experiment are well reproduced. The models were thus able to reproduce the main physical processes of the experiment, both for vapour-dominated diffusion during the heating phase and combined liquid and vapour transport during a subsequent heating and hydration phase. Lastly, based on the parameter sets, the teams predict a penetration of the water infiltration front in the 48-cm column filled with bentonite pellets to a depth between 25 and 35 cm over the 15,000 h (i.e. over 20 months) of the hydration phase of the experiment.

58 GEOSCIENCES↗