Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SCALE MODELS”

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 109 records · Page 6

High-fidelity micro-scale modeling of the thermo-visco-plastic behavior of carbon fiber polymer matrix composites

An experimentally validated micro-scale analysis of the visco-thermo-mechanical behavior of polymer matrix composites under different loads is proposed. A new constitutive law for the matrix material is developed taking into account the pressure dependence of the material as well as strain-rate and temperature dependence. Capturing the matrix behavior under multi-axial stress states is concluded to be essential to accurately predict the composite material behavior, even when considering simple load cases such as transverse compression and/or shear. Without any calibration procedure at the composite level, good agreement with the experimental data is observed for different loading conditions, including strain-rate dependency. Using this validated micro-scale model, a three-dimensional simulation of the formation of a kink band under longitudinal compression of the composite is conducted. A new evidence at micro-scale is found supporting the hypothesis that shear stresses transferred between fibers and matrix are particularly important in the formation of the kink band.

36 MATERIALS SCIENCE↗

Low length scale modeling of Cr precipitation in FeCrAl under thermal aging conditions

This report summarizes the progress on lower length scale modeling of Cr precipitation in FeCrAl alloys, which has been identified as a possible cladding material for accident tolerant fuels due to its excellent corrosion resistance via the formation of a protective, slow growing alumina surface layer. The precipitation of Cr-rich α' clusters in FeCrAl can lead to degraded mechanical properties such as the well-known “475oC embrittlement”, and should therefore be avoided to ensure cladding integrity. In FY16, a lattice kinetic Monte Carlo (LKMC) model was developed for FeCr binary alloys. Due to the assumption of constant bond energies, the previous model cannot correctly reproduce the low-temperature part of the Fe-Cr phase diagram. It is also not possible to model the anomalous ordering tendency and repulsive interactions between Cr atoms in α Fe at the dilute concentration regime using constant bond energies. To address these issues, a new local-environment-dependent bond energy formalism has been developed for Fe-Cr in FY17 that can simultaneously reproduce the negative mixing energy at low Cr concentrations and the positive mixing energy at high Cr concentrations, and the consequent large solubility of Cr in bcc Fe at low temperatures. Furthermore, the LKMC model has been extended to include Al to enable the assessment of Cr precipitation tendency in ternary FeCrAl alloys under thermal aging. Results from the present KMC simulations show good agreement with available experiments in the literature. This study is supported by the DOE NEAMS program Accident Tolerant Fuel project.

36 MATERIALS SCIENCE↗

Multi-scale modeling of the evolution of structure and properties in materials for nuclear energy applications

Nuclear energy is an important component of an overall strategy to address climate change. Idaho National Laboratory (INL) is the U.S. Department of Energy’s primary facility for research and development in nuclear science and technology for energy generation, supporting the improvement and life extension of the existing reactor fleet and the development and licensing of new reactor designs. Computational modeling is an important component of these activities, particularly in the area of materials for nuclear applications, where experimental data can be very challenging and expensive to acquire, and where data is especially scarce for new reactor designs. INL has used multi-scale modeling – linking atomistic, mesoscale, and engineering scales – to improve the ability to predict the performance of materials for nuclear energy applications. These modeling efforts make extensive of MOOSE (Multiphysics Object-Oriented Simulation Environment), a general-purpose open source finite element framework developed at INL. In this talk, I will give an overview of the approach and tools used, and several examples of application, including performance of nuclear fuels, understanding radiation-driven formation of nanoscale void and gas bubble superlattices, and powder densification through electric field assisted sintering.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Multi-scale modeling of the evolution of structure and properties in materials for nuclear energy applications

Nuclear energy is an important component of an overall strategy to address climate change. Idaho National Laboratory (INL) is the U.S. Department of Energy’s primary facility for research and development in nuclear science and technology for energy generation, supporting the improvement and life extension of the existing reactor fleet and the development and licensing of new reactor designs. Computational modeling is an important component of these activities, particularly in the area of materials for nuclear applications, where experimental data can be very challenging and expensive to acquire, and where data is especially scarce for new reactor designs. INL has used multi-scale modeling – linking atomistic, mesoscale, and engineering scales – to improve the ability to predict the performance of materials for nuclear energy applications. These modeling efforts make extensive of MOOSE (Multiphysics Object-Oriented Simulation Environment), a general-purpose open source finite element framework developed at INL. In this talk, I will give an overview of the approach and tools used, and several examples of application, including performance of nuclear fuels, understanding radiation-driven formation of nanoscale void and gas bubble superlattices, and powder densification through electric field assisted sintering.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Generalized Grain-Scale Model for the Non-Plasma and Plasma-Assisted Hydrogen Direct Reduction of Iron Ore

Direct Reduction of Iron ore using hydrogen (H-DRI) is a promising pathway towards efficient steelmaking and accurate predictive models are a necessity for scale-up and optimization of this technology. However, accurate models of this process remain limited because existing models oversimplify grain-scale phenomena, such as nonlinearity inside grain, self-sufficient porosity, surface reactions, and the role of plasma species. These phenomena are important for flash steelmaking and plasma-assisted H-DRI processes. To address this need, we present a phenomenological model for simulating H-DRI at the scale of a single micron-sized grain of the iron ore. We call this the Transient Reactive Grain Model (TRGM). TRGM incorporates key physical process: gas species transport, a chemical kinetics of material conversion, nanopore structural evolution and, adsorption-desorption surface kinetics at the reactive nanopore surface. The important contribution of this work is that the model provides a dependence on different reductant species, specifically hydrogen atoms versus molecules, so that role of hydrogen plasma reduction can be clarified compared to the use of pure hydrogen gas reduction. TRGM predictions agree well with experimental data for both molecular H2 reduction of Fe2O3 and plasma hydrogen reduction of Fe3O4. Results reveal species concentration gradients with a diffuse reaction zone, and enhanced hydrogen diffusion at the grain outer surface due to evolving porosity. These findings challenge common assumptions in existing models, including sharp reaction fronts, quasi-steady diffusion and kinetics, and the neglect of surface chemistry. As a generalized grain-scale model for H-DRI processes, TRGM has practical applications in flash steelmaking and in-flight reduction using both molecular and plasma hydrogen.

08 HYDROGEN↗

Mesa isochrones and stellar tracks (mist). I. Solar-scaled models

This is the first of a series of papers presenting the Modules for Experiments in Stellar Astrophysics (MESA) Isochrones and Stellar Tracks (MIST) project, a new comprehensive set of stellar evolutionary tracks and isochrones computed using MESA, a state-of-the-art open-source 1D stellar evolution package. In this work, we present models with solar-scaled abundance ratios covering a wide range of ages (5⩽log(Age) [year]⩽10.3), masses (0.1⩽M/M{sub ⊙}⩽300), and metallicities (−2.0⩽[Z/H]⩽0.5). The models are self-consistently and continuously evolved from the pre-main sequence (PMS) to the end of hydrogen burning, the white dwarf cooling sequence, or the end of carbon burning, depending on the initial mass. We also provide a grid of models evolved from the PMS to the end of core helium burning for −4.0⩽[Z/H]<−2.0. We showcase extensive comparisons with observational constraints as well as with some of the most widely used existing models in the literature. The evolutionary tracks and isochrones can be downloaded from the project website at http://waps.cfa.harvard.edu/MIST/.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluating WRF-Chem multi-scale model in simulating aerosol radiative properties over the tropics – A case study over India

We utilized WRF-Chem multi-scale model to simulate the regional distribution of aerosols, optical properties and its effect on radiation over India for a winter month. The model is evaluated using measurements obtained from upper-air soundings, AERONET sun photometers, various satellite instruments, and pyranometers operated by the Indian Meteorological Department. The simulated downward shortwave flux was overestimated when the effect of aerosols on radiation and clouds was neglected. Downward shortwave radiation from a simulation that included aerosol-radiation interaction processes was 5 to 25 Wm{sup -2} closer to the observations, while a simulation that included aerosol-cloud interaction processes were another 1 to 20 Wm{sup -2} closer to the observations. For the few observations available, the model usually underestimated particulate concentration. This is likely due to turbulent mixing, transport errors and the lack of secondary organic aerosol treatment in the model. The model efficiently captured the broad regional hotspots such as high aerosol optical depth over Indo-Gangetic basin as well as the northwestern and southern part of India. The regional distribution of aerosol optical depth compares well with AVHRR aerosol optical depth and the TOMS aerosol index. The magnitude and wavelength-dependence of simulated aerosol optical depth was also similar to the AERONET observations across India. Differences in surface shortwave radiation between simulations that included and neglected aerosol-radiation interactions were as high as -25 Wm{sup -2}, while differences in surface shortwave radiation between simulations that included and neglect aerosol-radiation-cloud interactions were as high as -30 Wm{sup -2}. The spatial variations of these differences were also compared with AVHRR observation. This study suggests that the model is able to qualitatively simulate the impact of aerosols on radiation over India; however, additional measurements of particulate mass and composition are needed to fully evaluate whether the aerosol precursor emissions are adequate when simulating radiative forcing in the region.

54 ENVIRONMENTAL SCIENCES↗

Pore-Scale Modeling of Electrokinetics in Geomaterials

Pore-scale finite-volume continuum models of electrokinetic processes are used to predict the Debye lengths, velocity, and potential profiles for two-dimensional arrays of circles, ellipses and squares with different orientations. The pore-scale continuum model solves the coupled Navier–Stokes, Poisson, and Nernst–Planck equations to characterize the electro-osmotic pressure and streaming potentials developed on the application of an external voltage and pressure difference, respectively. Here, this model is used to predict the macroscale permeabilities of geomaterials via the widely used Carmen–Kozeny equation and through the electrokinetic coupling coefficients. The permeability results for a two-dimensional X-ray tomography-derived sand microstructure are within the same order of magnitude as the experimentally calculated values. The effect of the particle aspect ratio and orientation on the electrokinetic coupling coefficients and subsequently the electrical and hydraulic tortuosity of the porous media has been determined. These calculations suggest a highly tortuous geomaterial can be efficient for applications like decontamination and desalination.

36 MATERIALS SCIENCE↗

SCALE Modeling of the Sodium Cooled Fast-Spectrum Advanced Burner Test Reactor

This report documents the modeling and simulation of a sodium-cooled fast reactor (SFR) as part of a U.S. Nuclear Regulatory Commission–sponsored project to assess the modeling and simulation capabilities for accident progression, source term, and consequence analysis for advanced reactor technologies with the Oak Ridge National Laboratory code SCALE and the Sandia National Laboratories (SNL) code MELCOR. Based on publicly available benchmark specifications, a fully heterogeneous 3D SCALE model of the 250 MWth Advanced Burner Test Reactor (ABTR) was developed to demonstrate SCALE’s capabilities for full-core reactivity analysis, fuel inventory prediction, and decay heat analysis of an SFR. The benchmark specifications contain modeling details for the ABTR core at the beginning of equilibrium cycle (BOEC) at operating conditions; they were derived from a 2006 preconceptual design report produced by Argonne National Laboratory. The ABTR was designed to demonstrate reactor-based transmutation of transuranics, that is, to “burn” transuranics recovered from light-water reactor (LWR) spent fuel. The ABTR’s fuel is designed to operate in 4 month cycles using uranium/transuranic (U/TRU) metallic fuel, with a TRU content of approximately 20%, at a conversion ratio of approximately 0.6. Various reactivity calculations were performed with SCALE for the ABTR and, where possible, compared with results available in the open literature. Additionally, SCALE was used to perform a full-core depletion calculation over the 4 month cycle to obtain the nuclide inventory at the end of equilibrium cycle (EOEC). These nuclide inventories, decay heat, power profiles, and reactivity feedback coefficients at EOEC represent the initial conditions for analyzing severe accident scenarios with MELCOR.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Principles of proteome allocation are revealed using proteomic data and genome-scale models

Integrating omics data to refine or make context-specific models is an active field of constraint-based modeling. Proteomics now cover over 95% of the Escherichia coli proteome by mass. Genome-scale models of Metabolism and macromolecular Expression (ME) compute proteome allocation linked to metabolism and fitness. Using proteomics data, we formulated allocation constraints for key proteome sectors in the ME model. The resulting calibrated model effectively computed the “generalist” (wild-type) E. coli proteome and phenotype across diverse growth environments. Across 15 growth conditions, prediction errors for growth rate and metabolic fluxes were 69% and 14% lower, respectively. The sector-constrained ME model thus represents a generalist ME model reflecting both growth rate maximization and “hedging” against uncertain environments and stresses, as indicated by significant enrichment of these sectors for the general stress response sigma factor σS. Finally, the sector constraints represent a general formalism for integrating omics data from any experimental condition into constraint-based ME models. The constraints can be fine-grained (individual proteins) or coarse-grained (functionally-related protein groups) as demonstrated here. Furthermore, this flexible formalism provides an accessible approach for narrowing the gap between the complexity captured by omics data and governing principles of proteome allocation described by systems-level models.

59 BASIC BIOLOGICAL SCIENCES↗

Multi-Scale Modeling of the Evolution of Structure and Properties in Materials for Nuclear Energy Applications [Slides]

Nuclear energy is an important component of an overall strategy to address climate change. Idaho National Laboratory (INL) is the U.S. Department of Energy’s primary facility for research and development in nuclear science and technology for energy generation, supporting the improvement and life extension of the existing reactor fleet and the development and licensing of new reactor designs. Computational modeling is an important component of these activities, particularly in the area of materials for nuclear applications, where experimental data can be very challenging and expensive to acquire, and where data is especially scarce for new reactor designs. INL has used multi-scale modeling – linking atomistic, mesoscale, and engineering scales – to improve the ability to predict the performance of materials for nuclear energy applications. In this talk, I will give an overview of the approach and tools used, and several examples of application, including performance of nuclear fuels, understanding radiation-driven formation of nanoscale void and gas bubble superlattices, and powder densification through electric field assisted sintering.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Multi‐Scale Model‐Informed Deep Learning for Plasma‐Nanoparticle Interaction

The Overarching Goal of this proposed research is to understand and quantitively determine the interactions between non-thermal plasma (hot electrons, reactive radicals, vibrationally excited species) and surface reactions on influencing the activity and selectivity of the desired reactions via developing multi-scale model informed deep learning algorithm. Investigating non-thermal plasma-surface interaction is feasible due to the low bulk temperature in the discharge region. To investigate the role of plasma-nanoparticle interaction on enhancing the reaction kinetics, we will focus on ammonia cracking to generate clean hydrogen over earth-abundant, non-critical metallic nanoparticles, which is of great significance for decarbonization. We hypothesize that (1) reactive radicals interacting with surface reaction species via Eley–Rideal mechanism will significantly lower the energetics of the potential rate-limiting step of nitrogen formation; (2) the surface will be charged heterogeneously under non-thermal plasma conditions and the charged site will lower the energetics of ammonia cracking through Langmuir– Hinshelwood mechanism; (3) vibrationally excited ammonia will further promote the initial N-H bond cleavage. To access the hypothesis, we will (1) reveal the surface charge effects on tunning the reaction energetics via interpretable, physics-informed deep learning accelerated density functional theory (DFT) calculations; (2) determine the reactive radicals interacting with surface reaction species on tuning the reaction energetics via DFT; (3) reveal the surface charge effects on tunning the reaction energetics via DFT and deep learning models, (4) quantify how vibrationally excited species, reactive radicals, and surface charging effects on enhancing the catalysis via developing DFT-based microkinetic modeling (MKM) and active learning. Deep and active learning of plasma-nanoparticle interactions effects on enhancing ammonia cracking to generate hydrogen represents a new paradigm for designing high performance plasma materials. The fundamental science of how plasma-nanoparticle interactions will change the plasma kinetics and will improve the energy efficiency for decarbonization and sustainability. The interpretable and physics-informed machine learning model will accelerate low temperature plasma chemistry and material discovery with physics rules and model interpretation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A simple dynamic subgrid-scale model for LES of particle-laden turbulence

In this study, a dynamic model for large-eddy simulations is proposed in order to describe the motion of small inertial particles in turbulent flows. The model is simple, involves no significant computational overhead, contains no adjustable parameters, and is flexible enough to be deployed in any type of flow solvers and grids, including unstructured setups. The approach is based on the use of elliptic differential filters to model the subgrid-scale velocity. The only model parameter, which is related to the nominal filter width, is determined dynamically by imposing consistency constraints on the estimated subgrid energetics. The performance of the model is tested in large-eddy simulations of homogeneous-isotropic turbulence laden with particles, where improved agreement with direct numerical simulation results is observed in the dispersed-phase statistics, including particle acceleration, local carrier-phase velocity, and preferential-concentration metrics.

42 ENGINEERING↗

Multi-scale modeling of inter-granular fracture in UO2

A hierarchical multi-scale approach is pursued in this work to investigate the influence of porosity, pore and grain size on the intergranular brittle fracture in UO2. In this approach, molecular dynamics simulations are performed to obtain the fracture properties for different grain boundary types. A phase-field model is then utilized to perform intergranular fracture simulations of representative microstructures with different porosities, pore and grain sizes. In these simulations the grain boundary fracture properties obtained from molecular dynamics simulations are used. The responses from the phase-field fracture simulations are then fitted with a stress-based brittle fracture model usable at the engineering scale. This approach encapsulates three different length and time scales, and allows the development of microstructurally informed engineering scale model from properties evaluated at the atomistic scale.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Reprocessing plant scale model integration

The national laboratory complex contains a number of niche modeling capabilities for fuel cycle facilities, at different scales and utilizing a range of codes or modeling platforms. The goal of this work has been to integrate modeling capabilities for reprocessing plants to develop a plant scale model. The Separations and Safeguards Performance Model (SSPM), developed in Matlab Simulink, has been used as a base to integrate more detailed codes. In the past, both the AMUSE and SEPHIS solvent extraction codes were integrated into the SSPM to provide more realistic modeling of the separations in aqueous reprocessing plants. For this work, the AMUSE integration was expanded to include time dependence and more fully integrated for both UREX and TRUEX contactor banks. The Nitron code is a new continuous dissolver model which has been examined for integration in the SSPM as well. Finally, in order to improve the model run time, a parallel processing capability was added. The integration progress along with initial results from testing are presented.

Cipiti, Benjamin B.↗

EchemAMR (electro-chemical microsctructure scale models with adaptive meshing) [SWR-23-111]

A 3D microstructure resolving electrochemical transport and interfacial chemistry solver. Electrode microstructure plays an important role in determining the performance of an electrochemical system, e.g. lithium ion battery. EchemAMR is a microstructure scale model that solves the governing equations for ion transport, electrical current continuity, interfacial chemistry and structural mechanics. Complex microstructure geometries from imaging can be directly imported into EchemAMR. A volume fraction based description of the geometry on Cartesian grid with an immersed interface formulation enables simplified meshing and large-scale simulations with millions of degrees of freedom. EchemAMR has been tested against systems with analytic solutions for numerical convergence and highly resolved lithium ion battery microstructures. EchemAMR demonstrates excellent mass conversation and efficient scaling on heterogenous High-Performance Computing (HPC) with central and graphics processing units.

Sitaraman, Hariswaran↗

Structural changes and variability of the ITCZ induced by radiation–cloud–convection–circulation interactions: inferences from the Goddard Multi-scale Modeling Framework (GMMF) experiments

In this paper, we have investigated the impact of radiation–cloud–convection–circulation interaction (RC3I) on structural changes and variability of the Inter-tropical Convergence Zone (ITCZ) using the Goddard Multi-scale Modeling Framework, where cloud processes are super-parameterized, i.e., explicitly resolved with 2-D cloud resolving models embedded in each coarse grid of the host Goddard Earth Observing System-Version 5 global climate model. Experiments have been conducted under prescribed sea surface temperature conditions for 10 years (2007–2016), with and without cloud radiation feedback in the atmosphere, respectively. Diagnostic analyses separately for January and July show that RC3I leads to an enhanced and expanded Hadley Circulation characterized by (1) a quasi-uniform warming and moistening of the tropical atmosphere and a sharpening of the ITCZ with enhanced deep convection, more intense precipitation and higher clouds, (2) extended drying of the tropical marginal convective zones, and extratropical mid- to lower troposphere, and (3) a cooling of the polar regions, with increased baroclinicity and midlatitude storm track activities. Computations based on the zonal mean thermodynamic energy balance equation show that the radiative warming and cooling are strongly balanced by local adiabatic processes associated with changes in large-scale vertical motions, as well as horizontal atmospheric heat transport. In the tropics, enhanced short-wave absorption and longwave water vapor greenhouse effects by high clouds play key roles in providing strong positive feedback to the tropospheric warming. In the extratropics, increased atmospheric heat transport associated with changes in the Hadley circulation is balanced by strong longwave cooling above, and warming below due to increased high clouds. We also find a strong positive correlation between daily and pentad heavy rain in the ITCZ core, and expansion of the drier zones coupled to a contraction of the highly convective zones in the ITCZ, indicating a strong tendency RC3I-induced convective aggregation in tropical clouds i.e., wet-regions-get-wetter and contracted, and dry-areas-get-drier and expanded.

54 ENVIRONMENTAL SCIENCES↗