SEARCH · Search NASA

Results for “Microscale 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.

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Assessing turbulence and mixing parameterizations in the gray-zone of multiscale simulations over mountainous terrain during the METEX21 field experiment

Multiscale numerical weather prediction models transition from mesoscale, where turbulence is fully parameterized, to microscale, where the majority of highly energetic scales of turbulence are resolved. The turbulence gray-zone is situated between these two regimes and multiscale models must downscale through these resolutions. Here, we compare three multiscale simulations which vary by the parameterization used for turbulence and mixing within the gray-zone. The three parameterizations analyzed are the Mellor-Yamada Nakanishi and Niino (MYNN) Level 2.5 planetary boundary layer scheme, the TKE-1.5 large eddy simulation (LES) closure scheme, and a recently developed three-dimensional planetary boundary layer scheme based on the Mellor-Yamada model. The simulation domain includes complex (i.e., mountainous) terrain in Nevada that was instrumented with meteorological towers, profiling and scanning lidars, a tethered balloon, and a surface flux tower. Simulations are compared to each other and to observations, with assessment of model skill at predicting wind speed, wind direction and TKE, and qualitative evaluations of transport and dispersion of smoke from controlled releases. This analysis demonstrates that microscale predictions of transport and dispersion can be significantly influenced by the choice of turbulence and mixing parameterization in the terra incognita , particularly over regions of complex terrain and with strong local forcing. This influence may not be apparent in the analysis of model skill, and motivates future field campaigns involving controlled tracer releases and corresponding modeling studies of the turbulence gray-zone.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Confronting the Challenge of Modeling Cloud and Precipitation Microphysics

In the atmosphere, microphysics refers to the microscale processes that affect cloud and precipitation particles and is a key linkage among the various components of Earth’s atmospheric water and energy cycles. The representation of microphysical processes in models continues to pose a major challenge leading to uncertainty in numerical weather forecasts and climate simulations. In this paper, the problem of treating microphysics in models is divided into two parts: i) how to represent the population of cloud and precipitation particles, given the impossibility of simulating all particles individually within a cloud, and ii) uncertainties in the microphysical process rates owing to fundamental gaps in knowledge of cloud physics. The recently-developed Lagrangian particle-based method is advocated as a way to address several conceptual and practical challenges of representing particle populations using traditional bulk and bin microphysics parameterization schemes. For addressing critical gaps in cloud physics knowledge, sustained investment for observational advances from laboratory experiments, new probe development, and next-generation instruments in space is needed. Greater emphasis on laboratory work, which has apparently declined over the past several decades relative to other areas of cloud physics research, is argued to be an essential ingredient for improving process-level understanding. More systematic use of natural cloud and precipitation observations to constrain microphysics schemes is also advocated. Because it is generally difficult to quantify individual microphysical process rates from these observations directly, this presents an inverse problem that can be viewed from the standpoint of Bayesian statistics. Following this idea, a probabilistic framework is proposed that combines elements from statistical and physical modeling. Besides providing rigorous constraint of schemes, there is an added benefit of quantifying uncertainty systematically. Finally, a broader hierarchical approach is proposed to accelerate improvements in microphysics schemes, leveraging the advances described in this paper related to process modeling (using Lagrangian particle-based schemes), laboratory experimentation, cloud and precipitation observations, and statistical methods.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Generalized thermo-mechanical framework for heterogeneous materials through asymptotic homogenization

Abstract A fundamental understanding of the interaction between microstructure and underlying physical mechanisms is essential, especially for developing more accurate multi-physics models for heterogeneous materials. Effects of microstructure on the material response at the macroscale are modeled by using the generalized thermomechanics. In this study, strain gradient theory is employed as a higher-order theory on the macroscale with thermodynamics modeled as a first-order theory on the microscale. Hence, energy depends only on the temperature such that we circumvent an extension of Fourier’s law and analyze the “simplest” thermo-mechanical model in strain gradient elasticity. Developing multiphysics models for heterogeneous materials is indeed a challenge and even this “simplest” model in generalized thermomechanics creates dozens of parameters to be determined. We develop a thermo-mechanical framework, in which microstructure is modeled as a periodic structure and through asymptotic homogenization approach, higher-order parameters at macroscopic scale are calculated. To illustrate the importance of higher-order parameters in overall thermo-mechanical response of a heterogeneous materials, finite element method (FEM) is employed with the aid of open-source codes (FEniCS). Verification example of a bulk system and several case studies of porous structures demonstrate how such numerical framework can be beneficial in the design of materials with tailored microstructures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Micro-cantilever beam experiments and modeling in porous polycrystalline UO 2

Understanding the impact of microstructure on the thermo-mechanical behavior of oxide nuclear fuels is vital to predicting their performance through multiscale models. Evaluating the mechanical properties at the sub-grain length scale is key to developing these multiscale models. In this work, 3D finite element (FE) models were constructed to simulate the micrometer-scale bending of micro-cantilever beams fabricated using porous polycrystalline uranium dioxide (UO 2 ) and tested at room temperature. Here, the results showed that the porosity and elastic anisotropy of individual grains can play a significant role in determining the effective mechanical properties of the material deduced from the tests. Specifically, the porosity had a non-negligible effect, given that the pore size was of the same order of magnitude as the dimensions of the micro-beams. Correlations between load-deflection data, pore location, and elastic properties (effective Young's modulus) were investigated using UO 2 micro-beam FE models, where pore clusters were included and placed at different locations along the length of the beam. Results indicated that the presence of pore clusters near the substrate, i.e., the clamp of the micro-cantilever beam, has the strongest effect on the load-deflection behavior, with the porosity leading to a reduction of stiffness that is the largest for any location of the pore clusters. Furthermore, it was also found that pore clusters located towards the middle of the span and close to the end of the beam have a comparatively small effect on the load-deflection behavior. Therefore, it is concluded that accurate estimates of Young's modulus can be obtained from micro-cantilever experiments after accounting for porosity on the one third of the beam length close to the clamp. This, in turn, provides an avenue to improve microscale experiments and their analysis in porous, anisotropic elastic materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

A concurrent irradiation-mechanics multiscale coupling model

The vast majority of our current knowledge regarding the basic mechanisms controlling irradiation effect on mechanical properties is based almost entirely on the results of post-irradiation experiments or theoretical models. However, the concurrent effects of irradiation, mechanical stress, and thermal damage on the failure phenomena of materials and components remain largely unexplored due to its internal multiscale-multiphysics coupling nature. Here, we present here a concurrent irradiation-mechanics multiscale coupling model. The concurrent evolutions of nanoscale irradiation defect clusters, microscale dislocation configurations, and mechanical responses are well captured based on coupling cluster dynamics, discrete dislocation dynamics, and the finite element methods using an effective time marching scheme. Model predictions of defect densities and size are in general agreement with experimental observations. Irradiation hardening is shown to take place also in samples undergoing concurrent irradiation-mechanical loading, similar to samples tested post-irradiation. However, the occurrence of plastic flow localization and dislocation channel formation is not accompanied with apparent yield drop (softening) under concurrent irradiation-mechanical loading conditions, which is different from the post-irradiation case.

42 ENGINEERING↗

DOE OSTI2023

Validity condition for the local sheath impedance boundary condition and a nonlocal generalization

ICRF sheaths can cause unwanted interactions of high-power RF waves with material surfaces in magnetic fusion devices. In previous work, a local RF sheath impedance boundary condition (BC) was derived for use in ICRF codes together with a microscale (i.e., Debye or sheath width scale) model for obtaining the sheath impedance used in that BC. This local RF sheath BC matches the normal component of current and electrostatic potential across the sheath-plasma interface. Collapsing the matching conditions at the sheath-plasma interface to a BC depends on the assumption of scale separation, which can be violated when conditions along the local radius of curvature of the surface vary sufficiently rapidly. The validity condition is explored in this contribution, with special attention to the case where the magnetic field approaches being tangent to the surface. When the local sheath BC no longer applies, a non-local sheath BC is developed under the assumption of a more relaxed scale separation assumption. It is shown that the non-local sheath BC reduces to the previous local sheath BC under appropriate conditions. Furthermore, a surface-integrated sheath admittance parameter describes the 2D physics in the new BC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Advancing material modeling in hydrocodes using a concurrent finite-element and molecular dynamics multiscale framework

We present a multiscale simulation framework that couples the finite-element method with molecular dynamics. Bypassing traditional equations of state (EOS) by using in-line atomistic simulations, the method offers the advantage of incorporating detailed microscale physics not easily represented with coarse-grained models. Coupling consistency with the continuum code is ensured through the use of lifting and restriction operators, in line with heterogeneous multiscale methods. The concurrent continuum-atomistic framework is validated through comparison with experimental results and conventional EOS models, and demonstrated in a shock-driven hydrodynamic flow simulation under extreme conditions. We further evaluate the framework's usability by comparing it to state-of-the-art EOS models of deuterium. A computational performance study reveals that the atomistic EOS evaluation is a feasible alternative to conventional approaches, and demonstrates a weak scaling of 99% efficiency. These results highlight the framework's potential for large-scale multiscale modeling across a broad range of materials and conditions.

Computer science↗

DOE OSTI2022

Stochastic Modeling in a Multimaterial Continuum Mixture Shock Physics Code

Stochastic modelling approaches are presented to capture random effects at multiple time and length scales. Random processes that occur at the microscale produce nondeterministic effects at the macroscale. Here we present three stochastic modeling approaches that describe random processes at microscopic length scales and map these processes to the macroscopic length scale. The first stochastic modeling approach is based upon a particle based numerical technique to solve a Stochastic Differential Equation (SDE) using an arbitrary diffusion process to capture random processes at the microstructural level. The second approach prescribes a Probability Density Function (PDF) for the drift and diffusion of the random variable derived using the forward and backward Kolmogorov equations. This method requires mean and drift evolution PDF transport equations. The third approach is the coupling of multiple random variables which are dependent on each other. The relationship of the PDFs and a coupling function, known as a copula, produces a Joint Probability Density Function (JPDF). These stochastic modeling approaches are implemented into a Multiple Component (MC) shock physics computational code and used to model statistical fracture and reactive flow applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Material Property Determinations for Rock Valley Legacy Samples in Support of Source Physics Experiment

The Source Physics Experiment (SPE) is a long-term NNSA research and development effort designed to improve nonproliferation verification and monitoring capabilities. The overarching goals of the SPE program are to improve understanding of prompt signals and physical signatures that develop from underground chemical explosions and associated modeling capabilities. Our work focuses on a primary factor controlling chemical explosion induced signals and signatures: the material properties of the rocks in which the chemical explosion takes place. This document reports on material property determinations of legacy core USGS Test Well F and outcrop analogs for the subsurface stratigraphy for the third phase of SPE in the Rock Valley (RV) area of the Nevada National Security Site (NNSS). The objective of this work is to establish a baseline set of lithologic descriptions and material properties expected prior to observatory borehole drilling in support of the SPE-RVDC (Rock Valley Direct Comparison) experiment. We determine for each rock type the compressional failure envelope, elastic properties as a function of stress (bulk modulus versus mean stress, shear modulus versus shear stress, Young’s modulus versus axial stress and Poisson’s ratio versus axial stress), indirect tensile strength, and porosity. Geologic characterization, both at the core-scale and microscale, provides context for using the data in modeling efforts and to inform interpretations for the material properties testing.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Life Expectancy of Evaporating Capillary Bridges Predicted by Tertiary Creep Modeling

The evaporation of capillary bridges is experimentally investigated at the microscale through a three-grain capillary cluster. This setting provides the minimum viable description of Haines jumps during evaporation, that is, capillary instabilities stemming from air entry into a saturated granular material. The displacement profile of a meniscus is obtained via digital image correlation for different grain materials, geometries, and separations. While it is well known that Haines jumps are triggered at the pore throat, we find that these instabilities are of three types depending on the separation. We also provide a temporal characterization of Haines jumps; we find that they are accurately described, as tertiary creep instabilities, by Voight’s relation, similarly to landslides and volcanic eruptions. This finding extends the description of capillary instabilities beyond their onset predicted by Laplace equilibrium. Our contribution also paves the way for a microscopically-informed description of desiccation cracks, of which Haines jumps are the precursors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Limitations of Bin and Bulk Microphysics in Reproducing the Observed Spatial Structure of Light Precipitation

Abstract Coarse-gridded atmospheric models often account for subgrid-scale variability by specifying probability distribution functions (PDFs) of process rate inputs such as cloud and rainwater mixing ratios ( q c and q r , respectively). PDF parameters can be obtained from numerous sources: in situ observations, ground- or space-based remote sensing, or fine-scale modeling such as large-eddy simulation (LES). LES is appealing to constrain PDFs because it generates large sample sizes, can simulate a variety of cloud regimes/case studies, and is not subject to the ambiguities of observations. However, despite the appeal of using model output for parameterization development, it has not been demonstrated that LES satisfactorily reproduces the observed spatial structure of microphysical fields. In this study, the structure of observed and modeled microphysical fields are compared by applying bifractal analysis, an approach that quantifies variability across spatial scales, to simulations of a drizzling stratocumulus field that span a range of domain sizes, drop concentrations (a proxy for mesoscale organization), and microphysics schemes (bulk and bin). Simulated q c closely matches observed estimates of bifractal parameters that measure smoothness and intermittency. There are major discrepancies between observed and simulated q r properties, though, with bulk simulated q r consistently displaying the bifractal properties of observed clouds (smooth, minimally intermittent) rather than rain while bin simulations produce q r that is appropriately intermittent but too smooth. These results suggest fundamental limitations of bulk and bin schemes to realistically represent higher-order statistics of the observed rain structure.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Grain2mesh: A Python and cubit mesh generator from unprocessed mesoscale images

Predicting bulk behavior from microscale features constitutes a key objective in multiscale modeling research, often involving numerical models composed of finite elements that capture the diversity of constituent phases, shapes, and orientations within the material. The Grain2mesh toolbox allows the user to input unprocessed mesoscopic images for automatic segmentation, pre-processing, quality control, and numerical mesh generation. The numerical mesh generation incorporates Cubit routines to generate robust multi-phase mesh structure for use in computational mechanics solvers. The python classes developed contain detailed documentation and examples to support standard usage and case-specific alternative options.

58 GEOSCIENCES↗

DOE OSTI2020

An anisotropic viscoplasticity model for shale based on layered microstructure homogenization

Viscoplastic deformation of shale is frequently observed in many subsurface applications. Many studies have suggested that this viscoplastic behavior is anisotropic—specifically, transversely isotropic—and closely linked to the layered composite structure at the microscale. In this work, we develop a two-scale constitutive model for shale in which anisotropic viscoplastic behavior naturally emerges from semianalytical homogenization of a bilayer microstructure. The microstructure is modeled as a composite of soft layers, representing a ductile matrix formed by clay and organics, and hard layers, corresponding to a brittle matrix composed of stiff minerals. This layered microstructure renders the macroscopic behavior anisotropic, even when the individual layers are modeled with isotropic constitutive laws. Using a common correlation between clay and organic content and magnitude of creep, we apply a viscoplastic modified Cam-Clay plasticity model to the soft layers, while treating the hard layers as a linear elastic material to minimize the number of calibration parameters. We then describe the implementation of the proposed model in a standard material update subroutine. The model is validated with laboratory creep data on samples from three gas shale formations. We also demonstrate the computational behavior of the proposed model through simulation of time-dependent borehole closure in a shale formation with different bedding plane directions.

58 GEOSCIENCES↗

DOE OSTI2024

A Multiphysics Thermoelastoviscoplastic Damage Internal State Variable Constitutive Model including Magnetism

We present a macroscale constitutive model that couples magnetism with thermal, elastic, plastic, and damage effects in an Internal State Variable (ISV) theory. Previous constitutive models did not include an interdependence between the internal magnetic (magnetostriction and magnetic flux) and mechanical fields. Although constitutive models explaining the mechanisms behind mechanical deformations caused by magnetization changes have been presented in the literature, they mainly focus on nanoscale structure–property relations. A fully coupled multiphysics macroscale ISV model presented herein admits lower length scale information from the nanoscale and microscale descriptions of the multiphysics behavior, thus capturing the effects of magnetic field forces with isotropic and anisotropic magnetization terms and moments under thermomechanical deformations. For the first time, this ISV modeling framework internally coheres to the kinematic, thermodynamic, and kinetic relationships of deformation using the evolving ISV histories. For the kinematics, a multiplicative decomposition of deformation gradient is employed including a magnetization term; hence, the Jacobian represents the conservation of mass and conservation of momentum including magnetism. The first and second laws of thermodynamics are used to constrain the appropriate constitutive relations through the Clausius–Duhem inequality. The kinetic framework employs a stress–strain relationship with a flow rule that couples the thermal, mechanical, and magnetic terms. Experimental data from the literature for three different materials (iron, nickel, and cobalt) are used to compare with the model’s results showing good correlations.

36 MATERIALS SCIENCE↗

DOE OSTI2024

A Multiphysics Thermoelastoviscoplastic Damage Internal State Variable Constitutive Model including Magnetism

36 MATERIALS SCIENCE↗

DOE OSTI2021

Characterize traction–separation relation and interfacial imperfections by data-driven machine learning models

Abstract Interfacial mechanical properties are important in composite materials and their applications, including vehicle structures, soft robotics, and aerospace. Determination of traction–separation (T–S) relations at interfaces in composites can lead to evaluations of structural reliability, mechanical robustness, and failures criteria. Accurate measurements on T–S relations remain challenging, since the interface interaction generally happens at microscale. With the emergence of machine learning (ML), data-driven model becomes an efficient method to predict the interfacial behaviors of composite materials and establish their mechanical models. Here, we combine ML, finite element analysis (FEA), and empirical experiments to develop data-driven models that characterize interfacial mechanical properties precisely. Specifically, eXtreme Gradient Boosting (XGBoost) multi-output regressions and classifier models are harnessed to investigate T–S relations and identify the imperfection locations at interface, respectively. The ML models are trained by macroscale force–displacement curves, which can be obtained from FEA and standard mechanical tests. The results show accurate predictions of T–S relations ( R 2 = 0.988) and identification of imperfection locations with 81% accuracy. Our models are experimentally validated by 3D printed double cantilever beam specimens from different materials. Furthermore, we provide a code package containing trained ML models, allowing other researchers to establish T–S relations for different material interfaces.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

The Need for Urban‐Resolving Climate Modeling Across Scales

Abstract Increasing urbanization, evolving urban landscapes, and growing populations require an accurate representation of urban areas and urban processes at microscale, regional, and global scales and their feedback across the scales. Important urban processes should be represented in climate models, and current approaches used are discussed. We discuss strategies for urban‐resolving climate research to capture the effect of climate on urban regions and the non‐linear urban processes on climate. Advancing urban‐resolving climate modeling will take us closer to capabilities that can truly capture relevant natural and human components.

Sharma, Ashish↗

DOE OSTI2021

Multi-scale Coupled Processes Modeling of Fractures as Porous, Interfacial and Granular Systems from Rock Images with the Numerical Manifold Method

Abstract The greatest challenges of rigorously modeling coupled hydro-mechanical processes in fractured rocks at different scales are associated with computational geometry. In addition, selections of continuous or discontinuous models, physical laws, and coupling priorities at different scales based on different geometric features determine the applicability of a numerical model for a certain type of problem. In this study, we present our multi-scale modeling capabilities that have been developed based on the numerical manifold method for analyzing coupled hydro-mechanical processes in fractured rocks. Based on their geometric features, the fractures are modeled as continua—finite-thickness porous zones, and discontinua—discontinuous interfaces and microscale asperities and granular systems. Different governing equations, physical laws, coupling priorities, and approaches for addressing fracture intersections and shearing are then applied to describe these. We applied these models to simulate coupled processes in fractured rocks using realistic geometry obtained from rock images at different scales. We first calculated shearing of a single fracture with different models and demonstrated the impacts of asperities on shearing. We then applied the continuous and discontinuous models to simulate a network of rough fractures, demonstrating that contact dynamics contribute significantly to the geometric, multi-physical evolution of systems where rough fractures are not mineral filled. For a discrete fracture network, our coupled processes modeling demonstrates that shearing of the discrete fractures can have a major impact on stress and pore pressure distribution. Lastly, we applied the discontinuous granular model to simulate evolution of a complex granular system with a deformation band, demonstrating that the deformation band can dominate contact dynamics, the structural and the stress evolution of the granular system.

58 GEOSCIENCES↗