Search NASA⌕ Search

SEARCH · Search NASA

Results for “non-linearly consistent implementation”

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 37 records · Page 2

Electron and ion kinetic effects on non-linearly driven electron plasma and ion acoustic waves

Fully non-linear kinetic simulations of electron plasma and ion acoustic waves (IAWs) have been carried out with a new multi-species, parallelized Vlasov code. The numerical implementation of the Vlasov model and the methods used to compute the wave frequency are described in detail. For the first time, the nonlinear frequency of IAWs, combining the contributions from electron and ion kinetic effects and from harmonic generation, has been calculated and compared to Vlasov results. Excellent agreement of theory with simulation results is shown at all amplitudes, harmonic generation being an essential component at large amplitudes. For IAWs, the positive frequency shift from trapped electrons is confirmed and is dominant for the effective electron-to-ion temperature ratio, Z T{sub e}/T{sub i} Greater-Than-Or-Equivalent-To 10 with Z as the charge state. Furthermore, numerical results demonstrate unambiguously the dependence [R. L. Dewar, Phys. Fluids 15, 712 (1972)] of the kinetic shifts on details of the distribution of the trapped particles, which depends in turn on the conditions under which the waves were generated. The trapped particle fractions and energy distributions are derived and, upon inclusion of harmonic effects, shown to agree with the simulation results, completing a consistent picture. Fluid models of the wave evolution are considered but prove unable to capture essential details of the kinetic simulations. Detrapping by collisions and sideloss is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Improving initialization and evolution accuracy of cosmological neutrino simulations

Neutrino mass constraints are a primary focus of current and future large-scale structure (LSS) surveys. Non-linear LSS models rely heavily on cosmological simulations — the impact of massive neutrinos should therefore be included in these simulations in a realistic, computationally tractable, and controlled manner. Here, a recent proposal to reduce the related computational cost employs a symmetric neutrino momentum sampling strategy in the initial conditions. We implement a modified version of this strategy into the Hardware/Hybrid Accelerated Cosmology Code (HACC) and perform convergence tests on its internal parameters. We illustrate that this method can impart Script $\mathcal{O}$(1%) numerical artifacts on the total matter field on small scales, similar to previous findings, and present a method to remove these artifacts using Fourier-space filtering of the neutrino density field. Moreover, we show that the converged neutrino power spectrum does not follow linear theory predictions on relatively large scales at early times at the 15% level, prompting a more careful study of systematics in particle-based neutrino simulations. We also present an improved method for backscaling linear transfer functions for initial conditions in massive neutrino cosmologies that is based on achieving the same relative neutrino growth as computed with Boltzmann solvers. Our self-consistent backscaling method yields sub-percent accuracy in the total matter growth function. Comparisons for the non-linear power spectrum with the Mira-Titan emulator at a neutrino mass of m ν = 0.15 eV are in very good agreement with the expected level of errors in the emulator and in the direct N-body simulation.

79 ASTRONOMY AND ASTROPHYSICS↗

An overview of recent physics results from NSTX

Currently, the National Spherical Torus Experiment (NSTX) is being upgraded to operate at twice the toroidal field and plasma current (up to 1 T and 2 MA), with a second, more tangentially aimed neutral beam (NB) for current and rotation control, allowing for pulse lengths up to 5 s. Recent NSTX physics analyses have addressed topics that will allow NSTX-Upgrade to achieve the research goals critical to a Fusion Nuclear Science Facility. These include producing stable, 100% non-inductive operation in high-performance plasmas, assessing plasma-material interface (PMI) solutions to handle the high heat loads expected in the next-step devices and exploring the unique spherical torus (ST) parameter regimes to advance predictive capability. Non-inductive operation and current profile control in NSTX-U will be facilitated by co-axial helicity injection (CHI) as well as radio frequency (RF) and NB heating. CHI studies using NIMROD indicate that the reconnection process is consistent with the 2D Sweet-Parker theory. Full-wave AORSA simulations show that RF power losses in the scrape-off layer (SOL) increase significantly for both NSTX and NSTX-U when the launched waves propagate in the SOL. Moreover, Toroidal Alfven eigenmode avalanches and higher frequency Alfven eigenmodes can affect NB-driven current through energy loss and redistribution of fast ions. The inclusion of rotation and kinetic resonances, which depend on collisionality, is necessary for predicting experimental stability thresholds of fast growing ideal wall and resistive wall modes. Neutral beams and neoclassical toroidal viscosity generated from applied 3D fields can be used as actuators to produce rotation profiles optimized for global stability. DEGAS-2 has been used to study the dependence of gas penetration on SOL temperatures and densities for the MGI system being implemented on the Upgrade for disruption mitigation. PMI studies have focused on the effect of ELMs and 3D fields on plasma detachment and heat flux handling. Simulations indicate that snowflake and impurity seeded radiative divertors are candidates for heat flux mitigation in NSTX-U. Studies of lithium evaporation on graphite surfaces indicate that lithium increases oxygen surface concentrations on graphite, and deuterium-oxygen affinity, which increases deuterium pumping and reduces recycling. In situ and test-stand experiments of lithiated graphite and molybdenum indicate temperature-enhanced sputtering, although that test-stand studies also show the potential for heat flux reduction through lithium vapour shielding. Non-linear gyro kinetic simulations have indicated that ion transport can be enhanced by a shear-flow instability, and that non-local effects are necessary to explain the observed rapid changes in plasma turbulence. Predictive simulations have shown agreement between a microtearing-based reduced transport model and the measured electron temperatures in a microtearing unstable regime. Finally, two Alfven eigenmode-driven fast ion transport models have been developed and successfully benchmarked against NSTX data. Upgrade construction is moving on schedule with initial physics research operation of NSTX-U planned for mid-2015.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Accelerating phase field simulations through a hybrid adaptive Fourier neural operator with U-net backbone

Prolonged contact between a corrosive liquid and metal alloys can cause progressive dealloying. For one such process as liquid-metal dealloying (LMD), phase field models have been developed to understand the mechanisms leading to complex morphologies. However, the LMD governing equations in these models often involve coupled non-linear partial differential equations (PDE), which are challenging to solve numerically. In particular, numerical stiffness in the PDEs requires an extremely refined time step size (on the order of 10 -12 s or smaller). This computational bottleneck is especially problematic when running LMD simulation until a late time horizon is required. This motivates the development of surrogate models capable of leaping forward in time, by skipping several consecutive time steps at-once. In this paper, we propose a U-shaped adaptive Fourier neural operator (U-AFNO), a machine learning (ML) based model inspired by recent advances in neural operator learning. U-AFNO employs U-Nets for extracting and reconstructing local features within the physical fields, and passes the latent space through a vision transformer (ViT) implemented in the Fourier space (AFNO). We use U-AFNOs to learn the dynamics of mapping the field at a current time step into a later time step. We also identify global quantities of interest (QoI) describing the corrosion process (e.g., the deformation of the liquid-metal interface, lost metal, etc.) and show that our proposed U-AFNO model is able to accurately predict the field dynamics, in spite of the chaotic nature of LMD. Most notably, our model reproduces the key microstructure statistics and QoIs with a level of accuracy on par with the high-fidelity numerical solver, while achieving a significant 11, 200 × speed-up on a high-resolution grid when comparing the computational expense per time step. Finally, we also investigate the opportunity of using hybrid simulations, in which we alternate forward leaps in time using the U-AFNO with high-fidelity time stepping. We demonstrate that while advantageous for some surrogate model design choices, our proposed U-AFNO model in fully auto-regressive settings consistently outperforms hybrid schemes.

36 MATERIALS SCIENCE↗

Predicting concrete compressive strength using hybrid ensembling of surrogate machine learning models

This study aims to implement a hybrid ensemble surrogate machine learning technique in predicting the compressive strength (CS) of concrete, an important parameter used for durability design and service life prediction of concrete structures in civil engineering projects. For this purpose, an experimental database consisting of 1030 records has been compiled from the machine learning repository of the University of California, Irvine. The database was used to train and validate four conventional machine learning (CML) models, namely Artificial Neural Network (ANN), Linear and Non-Linear Multivariate Adaptive Regression Splines (MARS-L and MARS-C), Gaussian Process Regression (GPR), and Minimax Probability Machine Regression (MPMR). Subsequently, the predicted outputs of CML models were combined and trained using ANN to construct the Hybrid Ensemble Model (HENSM). It is observed that the proposed HENSM produces higher predictive accuracy compared to the CML models used in the present study. The predictive performance of all models for CS prediction was compared using the testing dataset and it is found that the HENSM model attained the highest predictive accuracy in both phases. Based on the experimental results, the newly constructed HENSM model is very potential to be a new alternative in handling the overfitting issues of CML models and hence, can be used to predict the concrete CS, including the design of less polluting and more sustainable concrete constructions.

36 MATERIALS SCIENCE↗

Estimating the Electric Fields Driving Lightning Dart Leader Development With BIMAP‐3D Observations

Abstract In this paper, a numerical dart leader model has been implemented to understand the leader's development and the corresponding electric field changes observed by the 3D broadband mapping and polarization (BIMAP‐3D) system. The model assumes the extending leader channel is equipotential and has a linear charge distribution induced by an ambient electric field. The charge distribution induced by the ambient field can be used to model the electric field change at the ground. We then find the ambient electric field which best fits the field change measurements at the two BIMAP stations. The estimated ambient electric field decreases in the direction of dart leader propagation. Our observations and modeling results are consistent with our earlier hypothesis that dart leader speed is proportional to the electric field at the leader tip. The model also supports our earlier analysis that leader speed variations near branch junctions were due to previous charge deposits near the junctions. The modeled tip electric field is generally lower than the breakdown field unless the pre‐dart‐leader channel has a significant temperature of 3,000 K. This is consistent with the fact that dart leaders typically do not form new branches into the virgin air. Furthermore, the tip field is generally close to the negative streamer stability field at ambient temperatures, explaining the nature of the narrow and well‐defined channel structure. In addition to the charge distribution and the ambient and tip electric field, the development of the channel potential and current distribution are also presented.

3D lightning mapping↗

Using the Gerchberg-Saxton algorithm to reconstruct nonmodulated pyramid wavefront sensor measurements

Adaptive optics (AO) is a technique to improve the resolution of ground-based telescopes by correcting, in real-time, optical aberrations due to atmospheric turbulence and the telescope itself. With the rise of Giant Segmented Mirror Telescopes (GSMT), AO is needed more than ever to reach the full potential of these future observatories. One of the main performance drivers of an AO system is the wavefront sensing operation, consisting of measuring the shape of the above mentioned optical aberrations. Aims. The non-modulated pyramid wavefront sensor (nPWFS) is a wavefront sensor with high sensitivity, allowing the limits of AO systems to be pushed. The high sensitivity comes at the expense of its dynamic range, which makes it a highly non-linear sensor. We propose here a novel way to invert nPWFS signals by using the principle of reciprocity of light propagation and the Gerchberg-Saxton (GS) algorithm. We test the performance of this reconstructor in two steps: the technique is first implemented in simulations, where some of its basic properties are studied. Then, the GS reconstructor is tested on the Santa Cruz Extreme Adaptive optics Laboratory (SEAL) testbed located at the University of California Santa Cruz. This new way to invert the nPWFS measurements allows us to drastically increase the dynamic range of the reconstruction for the nPWFS, pushing the dynamics close to a modulated PWFS. The reconstructor is an iterative algorithm requiring heavy computational burden, which could be an issue for real-time purposes in its current implementation. However, this new reconstructor could still be helpful in the case of many wavefront control operations. In conclusion, this reconstruction technique has also been successfully tested on the Santa Cruz Extreme AO Laboratory (SEAL) bench where it is now used as the standard way to invert nPWFS signal.

47 OTHER INSTRUMENTATION↗

A review of pressure anisotropy caused by electron trapping in collisionless plasma, and its implications for magnetic reconnection

From spacecraft data, it is evident that electron pressure anisotropy develops in collisionless plasmas. This is in contrast to the results of theoretical investigations, which suggest this anisotropy should be limited. Common for such theoretical studies is that the effects of electron trapping are not included; simply speaking, electron trapping is a non-linear effect and is, therefore, eliminated when utilizing the standard methods for linearizing the underlying kinetic equations. Here, we review our recent work on the anisotropy that develops when retaining the effects of electron trapping. A general analytic model is derived for the electron guiding center distribution f(v{sub ∥},v{sub ⊥}) of an expanding flux tube. The model is consistent with anisotropic distributions observed by spacecraft, and is applied as a fluid closure yielding anisotropic equations of state for the parallel and perpendicular components (relative to the local magnetic field direction) of the electron pressure. In the context of reconnection, the new closure accounts for the strong pressure anisotropy that develops in the reconnection regions. It is shown that for generic reconnection in a collisionless plasma nearly all thermal electrons are trapped, and dominate the properties of the electron fluid. A new numerical code is developed implementing the anisotropic closure within the standard two-fluid framework. The code accurately reproduces the detailed structure of the reconnection region observed in fully kinetic simulations. These results emphasize the important role of pressure anisotropy for the reconnection process. In particular, for reconnection geometries characterized by small values of the normalized upstream electron pressure, β{sub e∞}, the pressure anisotropy becomes large with p{sub ∥}≫p{sub ⊥} and strong parallel electric fields develop in conjunction with this anisotropy. The parallel electric fields can be sustained over large spatial scales and, therefore, become important for electron acceleration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A generalized Poisson and Poisson-Boltzmann solver for electrostatic environments

The computational study of chemical reactions in complex, wet environments is critical for applications in many fields. It is often essential to study chemical reactions in the presence of applied electrochemical potentials, taking into account the non-trivial electrostatic screening coming from the solvent and the electrolytes. As a consequence, the electrostatic potential has to be found by solving the generalized Poisson and the Poisson-Boltzmann equations for neutral and ionic solutions, respectively. In the present work, solvers for both problems have been developed. A preconditioned conjugate gradient method has been implemented for the solution of the generalized Poisson equation and the linear regime of the Poisson-Boltzmann, allowing to solve iteratively the minimization problem with some ten iterations of the ordinary Poisson equation solver. In addition, a self-consistent procedure enables us to solve the non-linear Poisson-Boltzmann problem. Both solvers exhibit very high accuracy and parallel efficiency and allow for the treatment of periodic, free, and slab boundary conditions. The solver has been integrated into the BigDFT and Quantum-ESPRESSO electronic-structure packages and will be released as an independent program, suitable for integration in other codes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

GPAW: An open Python package for electronic structure calculations

We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually independent and can be connected by transformations via the real-space grid. This multi-basis feature renders GPAW highly versatile and unique among similar codes. By virtue of its modular structure, the GPAW code constitutes an ideal platform for the implementation of new features and methodologies. Moreover, it is well integrated with the Atomic Simulation Environment (ASE), providing a flexible and dynamic user interface. In addition to ground-state DFT calculations, GPAW supports many-body GW band structures, optical excitations from the Bethe–Salpeter Equation, variational calculations of excited states in molecules and solids via direct optimization, and real-time propagation of the Kohn–Sham equations within time-dependent DFT. A range of more advanced methods to describe magnetic excitations and non-collinear magnetism in solids are also now available. In addition, GPAW can calculate non-linear optical tensors of solids, charged crystal point defects, and much more. Recently, support for graphics processing unit (GPU) acceleration has been achieved with minor modifications to the GPAW code thanks to the CuPy library. We end the review with an outlook, describing some future plans for GPAW.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗

A new method to measure galaxy bias by combining the density and weak lensing fields

We present a new method to measure redshift-dependent galaxy bias by combining information from the galaxy density field and the weak lensing field. This method is based on the work of Amara et al., who use the galaxy density field to construct a bias-weighted convergence field κg. The main difference between Amara et al.'s work and our new implementation is that here we present another way to measure galaxy bias, using tomography instead of bias parametrizations. The correlation between κg and the true lensing field κ allows us to measure galaxy bias using different zero-lag correlations, such as <κgκ>/<κκ> or <κgκg>/<κgκ>. Our method measures the linear bias factor on linear scales, under the assumption of no stochasticity between galaxies and matter. We use the Marenostrum Institut de Ciències de l'Espai (MICE) simulation to measure the linear galaxy bias for a flux-limited sample (i < 22.5) in tomographic redshift bins using this method. This article is the first that studies the accuracy and systematic uncertainties associated with the implementation of the method and the regime in which it is consistent with the linear galaxy bias defined by projected two-point correlation functions (2PCF). We find that our method is consistent with a linear bias at the per cent level for scales larger than 30 arcmin, while non-linearities appear at smaller scales. This measurement is a good complement to other measurements of bias, since it does not depend strongly on σ8 as do the 2PCF measurements. We will apply this method to the Dark Energy Survey Science Verification data in a follow-up article.

79 ASTRONOMY AND ASTROPHYSICS↗

Energy occupation of waves and structures in 3D compressive MHD turbulence

Abstract Structures and propagating waves are often observed in solar wind turbulence. Their origins and features remain to be uncovered. In this work, we use 3D driven, compressible MHD turbulence simulations to investigate the global signatures of the driven fluctuations in whole spatial and temporal domain. With four-dimensional spatial-temporal (x, y, z, t) Fourier transformations implemented, we have identified two distinct main populations: waves, which satisfy the $\omega -\boldsymbol {k}$ dispersion relations and are propagating; and structures, which satisfy the polarization relations but non-propagating (ω = 0). Whereas the overall turbulent energy spectrum is still consistent with k−5/3, the contributions from waves and structures show very different behaviour in $\boldsymbol {k}$ space, with structures dominating at small k but waves becomes comparable to structures at large k. Overall, the fluctuations in the directions perpendicular to the large-scale mean field $\boldsymbol {B_0}$ are a manifestation of structures, while along the parallel direction, the fluctuations are dominated by waves. Also, a significant portion of the incompressible structures are the Alfvénic nature, and with imbalanced increased, the waves predominantly propagate in one direction and nearly perpendicular to $\boldsymbol {B_0}$. Differentiating the relative contributions from waves and structures could have important implications for understanding the non-linear cascade processes in the inertial range as well as particle-fluctuation interactions at small scales.

Yang, L. P.↗

Mechanical Behavior of UO 2 at Sub-grain Length Scales: Quantification of Elastic, Plastic and Creep Properties via Microscale Testing

Techniques were developed to measure properties at sub-grain scales using depleted Uranium Oxide (d-UO2) samples heat-treated to obtain different grain sizes and oxygen stoichiometries, through three main tasks: 1) sample processing and characterization, 2) microscale and conventional testing and 3) modeling. Grain size and crystallography were characterized using Scanning Electron Microscopy (SEM), in conjunction with Electron Backscattering Diffraction (EBSD) and Electron Channeling Contrast Imaging (ECCI). Grains were then carefully selected based on their crystallographic orientations to perform ex-situ micromechanical tests with samples machined via Focused Ion Beam (FIB), with emphasis on micro-cantilever bending. These experiments were performed under controlled atmospheres, to insure stoichiometry control, at temperatures up to 700 °C and allowed measurements involving elastic (effective Young’s modulus), plastic (critical resolved shear stresses) and creep (creep strain rates) behavior. Conventional compression experiments were performed simultaneously to compare with the ex-situ measurements and study potential size effects. Modeling was implemented using anisotropic elasticity and inelastic constitutive relations for plasticity and creep based on kinematics and kinetics of dislocation glide that account for the effects of crystal orientation, and stress. The models will be calibrated and validated using the experimental data. This project provided insight on correlations among stoichiometry, crystallography and mechanical behavior in advanced oxide fuels, provided valuable experimental data to validate and calibrate mesoscale fuel performance codes and also a framework to measure sub-grain scale mechanical properties that should be suitable for use with irradiated samples due to small volumes required. The goals and metrics of the ongoing study of thermo-mechanical behavior in depleted uranium dioxide (d-UO 2 ) outlined in this project have been concluded successfully, resulting in: 1) the successful fabrication, processing, and characterization of large-grained samples with various orientations (up to and including single crystals) having stoichiometric and hyper-stoichiometric O/U ratios; 2) formulation, calibration, and validation of a crystal plasticity constitutive model to describe the creep deformation of UO 2 at the sub-grain length scale (single crystal level) at intermediate temperatures; 3) the successful calibration of a crystal plasticity constitutive model to describe the elasto-plastic deformation of microcantilever beams, also at moderate temperatures. Samples were prepared from natural uranium oxide powder of production-quality provided by Areva. The powder was pressed in a die to a pressure of 100 MPa to produce green pellets with no sintering aids, lubricants, or any other additives. The green pellets were then heated up to 1700 °C under ultra-high purity argon atmosphere (~1 ppm O2). The atmosphere was then changed to 79% Argon, 21% O 2 and the temperature was held at 1700 °C for 2 hours to sinter the pellets under oxidative conditions [1] that are known to increase grain growth kinetics in UO 2 [2]. Samples were then cooled down under Ar-4%H2 atmosphere to reduce the samples back to stoichiometric UO 2 . For macro-scale procedures, testing of UO 2 samples with large grains was performed at 1200 °C using a modified load frame capable of applying dead-weight loads to ensure constant stress conditions, while displacement of the sample produced by the applied load was measured with high precision micrometers to obtain strains. Stress steps were used during testing and the strains were monitored to measured creep strain rates under steady state for each level of stress used, so that stress exponents could be obtained. The results of the mechanical testing, along with sample geometry and crystal orientation of the grains in the samples, as well as post-test sample characterization were used to formulate a viscoplastic model to account for steady state (stage II) creep behavior, along with basic assumptions from crystal plasticity and kinematic constraints due to testing fixtures. In the micro-scale, testing of microcantilever beams at temperatures ranging from 25 to 570 °C was performed in-situ with a scanning electron microscope with a special attachment to apply load and measure displacement while the samples were at temperature. The load-displacement curves obtained showed linear behavior before fracture for all temperatures attempted except 570 °C, where clear deviations from non-linearity were observed before fracture. These deviations were consistently observed for all samples tested for a given orientation. A viscoplastic model was used to account for the presence of inelastic strain, along with basic assumptions from crystal plasticity and beam theory. These models were kept as simple as possible, and results from tests performed in a set of samples with a given crystal orientation were used to calibrate the material constants for the model, while results from a different sample set were then used for validation, thus satisfying the conditions of all main tasks within the parameters of this project. Details of these efforts are outlined in this report.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

SAM Theory Manual

The System Analysis Module (SAM) is an advanced and modern system analysis tool under development at Argonne National Laboratory for advanced non-LWR reactor safety analysis. It aims to provide fast-running, modest-fidelity, whole-plant transient analyses capabilities, which are essential for fast turnaround design scoping and engineering analyses of advanced reactor concepts. While SAM is being developed as a system-level modeling and simulation tool, advanced modeling techniques being implemented include a reduced-order three-dimensional module, pseudo 3-D conjugate heat transfer modeling in reactor core, flexible and multi-scale modeling of heat transfer between fluid and structures, in addition to the advances in software environments and design, and numerical methods. SAM aims to be a generic system-level safety analysis tool for advanced non-LWRs, including Liquid-Metal-cooled fast Reactors (LMR), Molten Salt Reactors (MSR), Fluoride-salt-cooled High- temperature Reactors (FHR), and High-Temperature Gas-cooled Reactors (HTGR). SAM takes ad- vantage of advances in physical modeling, numerical methods, and software engineering to enhance its user experience and usability. It utilizes an object-oriented computational framework (MOOSE), and its underlying meshing and finite-element library and linear and non-linear solvers, to leverage the modern advanced software environments and numerical methods. This document provides the theoretical and technical basis of the code to help users understand the underlying physical models (such as governing equations, closure models, and component models), system modeling approaches, numerical discretization and solution methods, and the overall capabilities in SAM. As new code capabilities and features are added, the SAM Theory Manual will be updated periodically to keep it consistent with the state of the development.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

SAM Theory Manual

The System Analysis Module (SAM) is an advanced and modern system analysis tool under development at Argonne National Laboratory for advanced non-LWR reactor safety analysis. It aims to provide fast-running, modest-fidelity, whole-plant transient analyses capabilities, which are essential for fast turnaround design scoping and engineering analyses of advanced reactor concepts. While SAM is being developed as a system-level modeling and simulation tool, advanced modeling techniques being implemented include a reduced-order three-dimensional module, pseudo 3-D conjugate heat transfer modeling in reactor core, flexible and multi-scale modeling of heat transfer between fluid and structures, in addition to the advances in software environments and design, and numerical methods. SAM aims to be a generic system-level safety analysis tool for advanced non-LWRs, including Liquid-Metal-cooled fast Reactors (LMR), Molten Salt Reactors (MSR), Fluoride-salt-cooled High-temperature Reactors (FHR), and High-Temperature Gas-cooled Reactors (HTGR). SAM takes advantage of advances in physical modeling, numerical methods, and software engineering to enhance its user experience and usability. It utilizes an object-oriented computational framework (MOOSE), and its underlying meshing and finite-element library and linear and non-linear solvers, to leverage the modern advanced software environments and numerical methods. This document provides the theoretical and technical basis of the code to help users understand the underlying physical models (such as governing equations, closure models, and component models), system modeling approaches, numerical discretization and solution methods, and the overall capabilities in SAM. As new code capabilities and features are added, the SAM Theory Manual will be updated periodically to keep it consistent with the state of the development.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Smoothed Particle Hydrodynamics Model for Ice Sheet and Ice Shelf Dynamics

Mathematical modeling of ice sheets is complicated by the non-linearity of the governing equations and boundary conditions. Standard grid-based methods require complex front tracking techniques and have limited capability to handle large material deformations and abrupt changes in bottom topography. As a consequence, numerical methods are usually restricted to shallow ice sheet and ice shelf approximations. We propose a new smoothed particle hydrodynamics (SPH) model for coupled ice sheet and ice shelf dynamics. SPH is a fully Lagrangian particle method. It is highly scalable and its Lagrangian nature and meshless discretization are well suited to the simulation of free surface flows, large material deformation, and material fragmentation. In this paper SPH is used to study ice sheet/ice shelf behavior, and the dynamics of the grounding line. The steady state position of the grounding line obtained from the SPH simulations is in good agreement with laboratory observations for a wide range of simulated bedrock slopes, and density ratios similar to those of ice and sea water. The numerical accuracy of the SPH algorithm is further verified by simulating the plane shear flow of two immiscible fluids and the propagation of a highly viscous blob of fluid along a horizontal surface. In the experiment, the ice was represented with a viscous newtonian fluid. For consistency, in the described SPH model the ice is also modeled as a viscous newtonian fluid. Typically, ice sheets are modeled as a non-Newtonian fluid, accounting for the changes in the mechanical properties of ice. Implementation of a non-Newtonian rheology in the SPH model is the subject of our ongoing research.

54 ENVIRONMENTAL SCIENCES↗