Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum modeling”

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 505 records · Page 28

Quantitative study of the effect of grain boundary parameters on the slip system level Hall-Petch slope for basal slip system in Mg-4Al

Several theoretical studies have reported that the geometry and structure of grain boundaries in polycrystalline materials could impose a significant effect on the Hall-Petch slope. However, experimental observations are primarily limited by the ability of the techniques to accurately quantify the grain boundary strength and validate these theoretical models. Using high-resolution electron backscatter diffraction (HR-EBSD), the local stress tensor ahead of a slip band blocked by a grain boundary was quantified and coupled with a continuum dislocation pile-up model to assess the barrier strength of specific grain boundaries to specific slip systems, referred to as micro-Hall-Petch coefficient. For basal slip system in a deformed Mg-4Al alloy, the micro-Hall-Petch coefficient ($k^{basal}_{μ}$) varied significantly, from 0.054 to 0.184 MPa - m 1/2 for nine different grain boundaries. These results were correlated with geometric descriptors of the respective grain boundaries, with three-dimensional GB profile additionally measured via focused ion beam milling. It was found that the angle between the two slip plane traces on the grain boundary plane was the most sensitive parameter affecting $k^{basal}_{μ}$, followed by the angle between the slip directions. A functional form for calculation of $k^{basal}_{μ}$ depending on these two angles is proposed to augment crystal plasticity constitutive models with slip resistance dependent on some measure of the grain size. The method allows a new pathway to calibrate grain size strengthening parameters in crystal plasticity models, allowing further computational investigations of the interrelationship between texture, grain morphology, and the Hall Petch effect.

36 MATERIALS SCIENCE↗

exoALMA. XI. ALMA Observations and Hydrodynamic Models of LkCa 15: Implications for Planetary Mass Companions in the Dust Continuum Cavity

In the past decade, the Atacama Large Millimeter/submillimeter Array (ALMA) has revealed a plethora of substructures in the disks surrounding young stars. These substructures have several proposed formation mechanisms, with one leading theory being the interaction between the disk and newly formed planets. In this Letter, we present high angular resolution ALMA observations of LkCa 15’s disk that reveal a striking difference in dust and CO emission morphology. The dust continuum emission shows a ringlike structure characterized by a dust-depleted inner region of ~40 au in radius. Conversely, the CO emission is radially smoother and shows no sign of gas depletion within the dust cavity. We compare the observations with models for the disk–planet interaction, including radiative transfer calculation in the dust and CO emission. This source is particularly interesting, as the presence of massive planets within the dust cavity has been suggested based on previous near-IR observations. We find that the level of CO emission observed within the dust cavity is inconsistent with the presence of planets more massive than Jupiter orbiting between 10 and 40 au. Instead, we argue that the LkCa 15 innermost dust cavity might be created either by a chain of low-mass planets or by other processes that do not require the presence of planets.

79 ASTRONOMY AND ASTROPHYSICS↗

Development of multi-scale computational frameworks to solve fusion materials science challenges

Over the past two decades, the US-DOE has funded multiple projects that rely on high-performance computing and exascale computing platforms to accelerate scientific discoveries and address grand scientific challenges, such as harnessing fusion energy. In this article, we review in detail one of these efforts aimed at enhancing our capability to model plasma-facing materials subject to plasma and high-energy ion/neutron irradiation. The plasma surface interactions project has built a multi-scale modeling framework where many of the plasma- and high-energy ion/neutron irradiation-induced effects occurring in tungsten are explored. Here, this knowledge is used to develop atomistically-informed, high-fidelity continuum and meso-scale models that can be validated against experiments. We review the developments within this project, with attention to experimental validation efforts, and specifically highlight activities associated with: helium bubble bursting and equation of state, and hydrogen-helium interactions in tungsten; atomistically-informed model development for beryllium-tungsten material mixing; coupling of scrape-of-layer plasma, sheath and material models; and coupling of stochastic cluster-dynamics and crystal plasticity models to address radiation effects in tungsten under stress. Finally, we present how the project is preparing for future computational architectures, for instance through efforts to adapt atomistic methods to exascale computing.

36 MATERIALS SCIENCE↗

Final Technical Report - Center for Simulation of Fusion Relevant RF Actuators

We have developed a suite of 3D electromagnetic field solvers, both FEM and FDTD based, that account for the RF antenna and vacuum vessel geometries with unprecedented accuracy. Workflows were developed that make it possible to translate CAD models for the antenna and vacuum vessel to physics meshes for RF wave simulation. Nonlinear RF sheath formation has been incorporated self-consistently as a boundary condition in these solvers. We have also carried out extensive studies of the impact of RF sheaths on the ion energy angle distribution at plasma-material interfaces, using high fidelity particle-in-cell codes. Comprehensive simulation models were developed to assess the impact of blob-like edge turbulence on RF wave propagation and the impact of the RF ponderomotive force on the plasma scrape-off layer (SOL). A fluid transport solver for the far-SOL was also developed which accounts for the high parallel to perpendicular heat anisotropy on an unstructured mesh, thus making it possible to precisely represent an antenna structure in the presence of edge transport. Finally we have developed a hierarchy of core wave propagation and absorption models that self-consistently combine continuum Fokker Planck and Monte Carlo treatments of fast ion evolution with ICRF full-wave field solvers and continuum Fokker Planck treatments of fast electron evolution with both full-wave and ray tracing models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Precision test of gauge/gravity duality in D0-brane matrix model at low temperature

We test the gauge/gravity duality between the matrix model and type IIA string theory at low temperatures with unprecedented accuracy. To this end, we perform lattice Monte Carlo simulations of the Berenstein-Maldacena-Nastase (BMN) matrix model, which is the one-parameter deformation of the Banks-Fischler-Shenker-Susskind (BFSS) matrix model, taking both the large N and continuum limits. We leverage the fact that sufficiently small flux parameters in the BMN matrix model have a negligible impact on the energy of the system while stabilizing the flat directions so that simulations at smaller N than in the BFSS matrix model are possible. Hence, we can perform a precision measurement of the large N continuum energy at the lowest temperatures to date. The energy is in perfect agreement with supergravity predictions including estimations of α'-corrections from previous simulations. At the lowest temperature where we can simulate efficiently (T = 0.25λ 1/3 , where λ is the ’t Hooft coupling), the difference in energy to the pure supergravity prediction is less than 10%. Furthermore, we can extract the coefficient of the 1/N 4 corrections at a fixed temperature with good accuracy, which was previously unknown.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Modeling support for the development of material surveillance specimens and procedures

This report describes modeling and simulation activities performed at Argonne National Laboratory supporting the development of passively actuated mechanical test articles for material surveillance in Molten Salt Reactors (MSRs). These test articles are a critical technology in formulating material surveillance programs for future MSRs to monitor the degradation in the structural properties of the materials in critical plant components. The main activity described in this report is the development of a method for inferring the amount of mechanical damage a test article has experienced during some duration of exposure to plant thermal and environmental conditions, using only mechanical test data collected from the test articles before and after exposure. The basic approach is to develop a model of the test article, including a description of mechanical degradation through continuum damage mechanics, and then use this model to cast the problem of inferring mechanical degradation in the test article materials into a shooting problem for a set of ordinary differential equations. We can then solve the shooting problem to determine the amount of damage accumulated in the sample. The report also details a few miscellaneous simulation studies completed at Argonne to support the development of the test articles themselves at Idaho National Laboratory.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Active Galactic Nuclei Continuum Reverberation Mapping Based on Zwicky Transient Facility Light Curves

We perform a systematic survey of active galactic nuclei (AGNs) continuum lags using ~3 days cadence gri-band light curves from the Zwicky Transient Facility. We select a sample of 94 type 1 AGNs at z < 0.8 with significant and consistent inter-band lags based on the interpolated cross-correlation function method and the Bayesian method JAVELIN. Within the framework of the "lamp-post" reprocessing model, our findings are: (1) The continuum emission (CE) sizes inferred from the data are larger than the disk sizes predicted by the standard thin-disk model. (2) For a subset of the sample, the CE size exceeds the theoretical limit of the self-gravity radius (12 lt-days) for geometrically thin disks. (3) The CE size scales with continuum luminosity as R CE ∝ L 0.48±0.04 with a scatter of 0.2 dex, analogous to the well-known radius–luminosity relation of broad Hβ. These findings suggest a significant contribution of diffuse continuum emission from the broad-line region (BLR) to AGN continuum lags. We find that the R CE –L relation can be explained by a photoionization model that assumes ~23% of the total flux comes from the diffuse BLR emission. In addition, the ratio of the CE size and model-predicted disk size anticorrelates with the continuum luminosity, which is indicative of a potential nondisk BLR lag contribution evolving with the luminosity. Finally, a robust positive correlation between the CE size and black hole mass is detected.

79 ASTRONOMY AND ASTROPHYSICS↗

Low-energy limit of some exotic lattice theories and UV/IR mixing

We continue our exploration of exotic, gapless lattice and continuum field theories with subsystem global symmetries. In an earlier paper, we presented free lattice models enjoying all the global symmetries (except continuous translations), dualities, and anomalies of the continuum theories. Here, we study in detail the relation between the lattice models and the corresponding continuum theories. We do that by analyzing the spectrum of the theories and several correlation functions. These lead us to uncover interesting subtleties in the way the continuum limit can be taken. In particular, in some cases, the infinite volume limit and the continuum limit do not commute. Furthermore, this signals a surprising UV/IR mixing, i.e., long distance sensitivity to short distance details.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modeling the Nonlinear Rheology of Polymer Additive Manufacturing

This report summarizes molecular and continuum simulation studies focused on developing physics - based predictive models for the evolution of polymer molecular order during the nonlinear processing flows of additive manufacturing. Our molecular simulations of polymer elongation flows identified novel mechanisms of fluid dissipation for various polymer architectures that might be harnessed to enhance material processability. In order to predict the complex thermal and flow history of polymer realistic additive manufacturing processes, we have developed and deployed a high - performance mesh - free hydrodynamics module in Sandia's LAMMPS software. This module called RHEO – short for Reproducing Hydrodynamics and Elastic Objects – hybridizes an updated - Lagrange reproducing - kernel method for complex fluids with a bonded particle method (BPM) to capture solidification and solid objects in multiphase flows. In combination, our two methods allow rapid, multiscale characterization of the hydrodynamics and molecular evolution of polymers in realistic processing geometries.

36 MATERIALS SCIENCE↗

Modeling Chemo-Mechanics with Electrolyte Infiltration to Quantify Degradation of Cathode Particles

One of the main goals in modeling lithium-ion batteries is to improve/predict longevity and resilience of new chemistries. To that end, this talk investigates the formation of stress-induced fracture within polycrystalline cathode particles and the impact on capacity loss. Physically based cathode aging dynamics is simulated in a 3D, continuum-level chemo-mechanical model. The model captures anisotropic Li diffusion within a single polycrystalline particle comprised of hundreds to thousands of randomly oriented grains. A recent addition to this model includes electrolyte infiltration, which occurs when the electrolyte seeps into surface cracks within the particle. The model predicts that particle fracture is primarily due to non-ideal grain interactions with slight dependence on high-rate charge demands. Essentially, when neighboring grains are misaligned, they expand a different rates relative to one another leading to high stresses and ultimately the formation of intraparticle cracks. The model predicts that small particles with large grains develop significantly less damage than larger particles with small grains. Finally, the model predicts most of the chemo-mechanical damage accumulates in the first charge after formation. This chemo-mechanical "damage saturation" effect indicates that initial particle fracture occurs within the first few cycles, while long-term cathode degradation is not solely chemo-mechanically induced. The principle contribution of this research is the use of an anisotropic chemo-mechanical model to test how particle geometry affect capacity fade, which predicts that particle size has a stronger effect on capacity fade than grain size and ultimately that small particles with large grains have the least capacity fade.

cathode↗

The charge reduction rate for multiply charged polymer ions via ion–ion recombination at atmospheric pressure

The charge reduction of multiply charged macromolecular ions via recombination with small ions in the gas phase is commonly employed to modulate the charge on macromolecules prior to mass spectrometric and mobility analyses. We employ a recently developed continuum-Molecular Dynamics (MD) calculation approach to determine the recombination rate coefficient of multiply charged (1 to 7 excess positive charged) polyethylene glycol ions (mass of 4600 Da) with smaller singly charged anions, modeled as NO 2 – ions. The continuum-MD approach accounts explicitly for the influence of the background gas on the recombination process, accounts explicitly for ion translational, vibrational, and rotational motion, and enables recombination rate coefficient calculation in nitrogen near atmospheric pressure, wherein neither low pressure nor high pressure recombination theories are strictly applicable. Continuum-MD simulations yield recombination rate coefficients near 3.9 × 10 –14 m 3 s –1 for singly charged ions, increasing to 3.0 × 10 –11 m 3 s –1 for the +7 ion. Pre-existing collision rate coefficient expressions for rigid ions are found to be within a factor of 2–5 of calculations for all charge states, but their use requires knowledge of an appropriate collision distance, which is not well-defined for flexible polymer ions. Continuum-MD-inferred rate coefficients are incorporated into a model of charge reduction, and the charge state distribution versus anion concentration determined with it is compared to charge reduction measurements made via atmospheric pressure differential mobility analysis. Furthermore, good agreement is observed between simulations and experiments; although results are extremely sensitive to the recombination rate coefficients, experimental results are bound by models utilizing rates within a factor of 2 (0.5–2.0×) of the continuum-MD rates.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

A superconfiguration calculation of opacity with consistent bound and continuum electron treatments using green’s functions

One of the challenges in calculating the opacity of dense plasmas is the difficulty in consistently modeling electrons bound to nuclei and those that exist within the continuum of free states in electronic structure models. We address this issue by adapting the green's function approach, originally developed for use in average atom calculations, to the determination of superconfiguration (SC) electronic structure. The spectra created using these SCs indicate that a consistent treatment of continuum electronic structure is important for phenomena involving electrons near ionization thresholds, such as the pressure ionization of bound states and the opacity due to transitions near bound-free edges. Though important for dense plasmas, the detailed incorporation of continuum electrons into structure calculations does not have significant impact on the recent discrepancies between the predicted and measured opacity of hot, dense iron (Bailey et al 2015 Nature517 56). Here we find that the inclusion of plasma effects through an ion-sphere model along with our treatment of continuum electronic states gives a description of pressure ionization in hot, dense aluminum that is in better agreement with experiment than methods that rely on perturbative descriptions of the plasma environment (Hoarty et al 2013 Phys. Rev. Lett.110 265003).

74 ATOMIC AND MOLECULAR PHYSICS↗

Tardigrade-examples V0.1.0

Tardigrade-examples is a repository of computational workflows that exercise the Tardigrade software package. The Tardigrade software package is an implementation of Eringen’s micromorphic continuum theory with capabilities to support multiscale material modeling. These capabilities include homogenization through the Micromorphic Filter, calibration of micromorphic material models, and macroscale simulation in Tardigrade-MOOSE. This repository investigates continuum upscaling of various direct numerical simulations (DNS) conducted in Abaqus finite element (FE), Ratel FE, and GEOS material point method (MPM) software. Verification of the upscaling workflow is first investigated by considering DNS of trivial stress states for homogeneous materials, results of which indicate that classical continuum behavior is recovered as expected. DNS of heterogeneous materials are then considered.

Allard, Thomas↗

Lattice regularization of reduced Kähler-Dirac fermions and connections to chiral fermions

We show how a path integral for reduced Kähler-Dirac fermions suffers from a phase ambiguity associated with the fermion measure that is an analog of the measure problem seen for chiral fermions. However, unlike the case of chiral symmetry, a doubler free lattice action exists which is invariant under the corresponding onsite symmetry. This allows for a clear diagnosis and solution to the problem using mirror fermions resulting in a unique gauge invariant measure. By introducing an appropriate set of Yukawa interactions which are consistent with ’t Hooft anomaly cancellation we conjecture the mirrors can be decoupled from low energy physics. Moreover, the minimal such Kähler-Dirac mirror model yields a light sector which corresponds, in the flat space continuum limit, to the Pati-Salam GUT model.

Physics↗

Climate-Relevant Gas Absorption Properties from AWARE and Other ARM Spectral Measurements (Final Report)

The objective of this project has been to use measurements from the Atmospheric Emitted Radiance Interferometer (AERI) deployed at ARM sites to improve our knowledge of uncertain infrared spectroscopic parameters of importance to climate, remote sensing, and data assimilation. A primary focus has been the water vapor continuum in the infrared atmospheric window. Radiation codes used to predict climate and weather base their representation of the water vapor continuum in this window on the MT_CKD model, but in this spectral region the MT_CKD water vapor continuum absorption coefficients were derived almost two decades ago by an analysis of AERI measurements for conditions with a limited range of precipitable water vapor (PWV) and temperature values, leading to uncertainty in the derived coefficients. This project has included a new, comprehensive analysis of all aspects of the water vapor continuum in the atmospheric window -- the self continuum, the foreign continuum, and the self continuum temperature dependence, all resolved spectrally. Having a wide range of PWV amounts and temperatures from deployments of the ARM Mobile Facility for GoAmazon as well as the more than a decade of AERI measurements at SGP has provided an appropriate foundation for such a comprehensive analysis.

54 ENVIRONMENTAL SCIENCES↗

A model for discrete fracture-clay rock interaction incorporating electrostatic effects on transport

Abstract A model based on the code CrunchClay is presented for a fracture-clay matrix system that takes electrostatic effects on transport into account. The electrostatic effects on transport include those associated with the development of a diffusion potential as captured by the Nernst-Planck equation, and the formation of a diffuse layer bordering negatively charged clay particles within which partial anion exclusion occurs. The model is based on a dual continuum formulation that accounts for diffuse layer and bulk water pore space, providing a more flexible framework than is found in the classical mean electrostatic potential models. The diffuse layer model is obtained by volume averaging ion concentrations in the Poisson-Boltzmann equation, but also includes the treatment of longitudinal transport within this continuum. The calculation of transport within the bulk and diffuse layer porosity is based on a new formulation for the Nernst-Planck equation that considers averaging of diffusion coefficients and accumulation factors at grid cell interfaces. Equations for function residuals and the associated Jacobian matrix are presented such that the system of nonlinear differential-algebraic equations can be solved with Newton’s method. As an example, we consider a 2D system with a single discrete fracture within which flow and advective transport occurs that is coupled to diffusion in the clay-rich matrix. The simulation results demonstrate the lack of retardation for anions (e.g., 36 Cl − ) of the contaminant plume within the fracture flow system because they are largely excluded from the charged clay rock, while the migration of cations (e.g., 90 Sr ++ ) is more strongly attenuated. The diffusive loss of divalent cations in particular from the fracture is accentuated by their accumulation in the diffuse layer within the clay-rich matrix.

58 GEOSCIENCES↗

Scale Invariance of Hot Spot Formation in TATB High Explosives

Shock-induced detonation of insensitive high explosives based on 1,3,5-triamino-2,4,6-trinitrobenzene starts with formation of hot spots at microstructural defects but has eluded atomistic modeling treatment at micron length scales. To this end, we performed multimicron scale all-atom molecular dynamics (MD) simulations of hot spots that form during the collapse of cylindrical pores with diameters between 10 and 300 nm. Our MD simulations show that hot spots formed at pores larger than 20 nm exhibit temperature fields with scale-invariant features for sizes up to at least 300 nm. Through a continuum-based grain-scale modeling framework, we span and extend beyond the size scales currently accessible to MD and find that hot spot scale invariance is a general feature that arises when the mechanical strength is insensitive to strain rate. Finally, our results demonstrate the applicability of all-atom MD to simulate the complicated dynamical evolution of micron-sized systems and bolster confidence in insights from MD simulations of materials that exhibit strength with negligible rate dependence over the relevant intervals.

36 MATERIALS SCIENCE↗

Utah FORGE: 2023 Large Upscaled Discrete Fracture Network Models

This dataset includes the data and a report on the large upscaled discrete fracture network modeling done for the Utah FORGE project in 2023. The FORGE modeling team is making five discrete fracture network (DFN) realizations of a large reservoir model available to researchers. These models have been upscaled to a continuum mesh or grid at resolutions of 10 meters and 20 meters providing reservoir properties for fracture porosity, permeability, and compressibility. The models are available in both the reference global coordinate frame and a local coordinate frame aligned with principal stress directions.

15 GEOTHERMAL ENERGY↗