Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Homogenization theory”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Automated Symbolic Upscaling: 2. Model Generation for Extended Applicability Regimes

Abstract In this second part of the two paper series, we detail an algorithmic procedure for systematically implementing the generalized closure form strategy presented in Part 1. This strategy extends the applicability of homogenized models with respect to classical homogenization theory, as demonstrated in Part 1 where upscaled models are rigorously derived in moderately reactive physical regimes. After encoding the algorithm into Symbolica, an automated upscaling framework, we upscale two reactive mass transport problems and numerically validate the resulting nonlinear homogenized models by showing the absolute error estimates predicted by homogenization theory are satisfied. In both problems, nontrivial closure forms and closure problems are automatically formulated using the encoded strategy with no human interaction, nor prior knowledge regarding the closure required for the systems. We hope these demonstrations spark further interest in automated analytical frameworks for multiscale modeling, as such capabilities are invaluable for generating rigorous multiscale models of complex phenomena in porous media.

Pietrzyk, Kyle↗
DOE OSTI2023

Automated Symbolic Upscaling: 1. Model Generation for Extended Applicability Regimes

Abstract In porous media theory, upscaling techniques are fundamental to deriving rigorous Darcy‐scale models for flow and reactive transport in subsurface systems. Due to limitations in classical techniques, a number of ad hoc approaches have been proposed to address physical regimes in which reactive time scales are similar to, or faster than, diffusive time scales. In Part 1 of this two part series, we present a strategy for expanding the applicability of classical homogenization theory by generalizing the assumed closure form. We detail the implementation of this strategy on two reactive mass transport problems with moderately reactive physics. The strategy produces nontrivial homogenized models with emergent terms and effective parameters that couple reactive, diffusive, and advective transport. The differences in equation forms between the macroscopic and pore‐scale descriptions advise caution to further studies where the forms of macroscopic equations are assumed, as opposed to rigorously derived. Numerical validation is provided for each problem to show that the error estimates of homogenization theory are satisfied, and to justify the implemented strategy. In Part 2, the presented strategy is automated using symbolic computing to expedite its implementation.

Pietrzyk, Kyle↗
DOE OSTI2024

Theories of homogeneous and electrochemical electron transfer in complex media and interfaces (Final Technical Report)

This project makes the next step in establishing practical theories of charge transfer in complex media. The development of formal models is supported by extensive atomistic simulations, quantum calculations of force-field parameters, and direct measurements of charge-transfer spectra. All theory development is supported by experiment, extensive numerical simulations, and through external collaborations.

14 SOLAR ENERGY↗
DOE OSTI2024

Enforcing global constraints for the dispersion closure problem: τ 2 -SIMPLE algorithm

Permeability and effective dispersion tensors are critical parameters to characterize flow and transport in porous media at the continuum scale. Homogenization theory defines a framework in which such effective properties are first computed from solving a closure problem in a repeating unit cell of the periodic microstructure and then used in a macroscopic formulation for efficient computation. The closure problem is formulated as a local boundary value problem subjected to global constraints, which guarantee the uniqueness of the solution and can be difficult to satisfy for complex geometries and at high flow conditions. These constraints also ensure that pore-scale pressure, velocity, and concentration fields can be accurately reconstructed from the closure variable. Building on a previous work, here we present a framework that allows to satisfy global constraints associated to both the permeability and the dispersion closure problems by introducing two artificial time scales. The algorithm, called τ 2 -SIMPLE, computes both permeability and effective dispersion given an arbitrarily complex geometry and flow condition. Furthermore, this algorithm is demonstrated to be accurate for both 2D and 3D geometries across varying flow conditions, and thus it can be used to quickly characterize effective properties from porous media images in many applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

The Thermo-Elastic Properties and Damping of U-6wt%Nb

While thermal expansion data exists for quenched (as well as aged) U-6wt%Nb, there is wide variation in the reported room temperature elastic moduli. To better understand the room temperature elastic behavior and to address the complete absence of data on the temperature dependence of the elastic moduli, room temperature and in-situ elevated temperature resonant ultrasonic spectroscopy (RUS) was performed along with x-ray diffraction (XRD) and dilatometry using a thermomechanical analyzer (TMA). An in-situ small- and wide-angle x-ray scattering (SAXS/WAXS) experiment was performed on a companion sample to help interpret the results. The room temperature polycrystalline dynamic moduli were measured to be E = 95.4 ± 4.5 GPa and G = 35.2 ± 0.2 GPa, respectively. As the temperature is raised, the stiffness slowly decreases, consistent with the empirical rule proposed by Varshni, up to the point of precipitation of the equilibrium α-U phase, which is stiffer than the martensitic phases. Homogenization theory of composites can be used to rationalize the observed response at high temperatures and after high-temperature exposures. Further, aging of the material at low temperatures (≤200°C) does not affect the linear elastic stiffness, but does impact the damping behavior which can be measured through RUS. This change in damping provides another perspective on the microstructure changes induced by low-temperature aging which also result in significant strengthening.

36 MATERIALS SCIENCE↗
DOE OSTI2021

An FFT-based approach for Bloch wave analysis: application to polycrystals

A method based on the Fast Fourier Transform is proposed to obtain the dispersion relation of acoustic waves in heterogeneous periodic media with arbitrary microstructures. The microstructure is explicitly considered using a voxelized Representative Volume Element (RVE). The dispersion diagram is obtained solving an eigenvalue problem for Bloch waves in Fourier space. To this aim, two linear operators representing stiffness and mass are defined through the use of differential operators in Fourier space. The smallest eigenvalues are obtained using the implicitly restarted Lanczos and the subspace iteration methods, and the required inverse of the stiffness operator is done using the conjugate gradient with a preconditioner. The method is used to study the propagation of acoustic waves in elastic polycrystals, showing the strong effect of crystal anistropy and polycrystaline texture on the propagation. It is shown that the method combines the simplicity of classical Fourier series analysis with the versatility of Finite Elements to account for complex geometries proving an efficient and general approach which allows the use of large RVEs in 3D.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A deep learning upscaling framework: Reactive transport and mineral precipitation in fracture-matrix systems

Pore-scale modeling has limited applicability at large scales due to its high computational cost. One common approach to upscale pore-scale models is the use of effective medium theories, which homogenize small-scale features in a porous structure and characterize the medium by macroscale properties (e.g., permeability) and equations (e.g., Darcy’s law). However, there are classes of physical processes for which effective medium approximations may become inaccurate, e.g., mineral precipitation and clogging during reactive transport. We have developed a deep learning upscaling framework, in which pore-scale modeling is directly employed in macroscale systems, without relying on effective medium approximations. The upscaling framework is first developed for general multiscale systems and then applied to modeling reactive transport with mineral precipitation in the altered layer in fracture-matrix structures. Solute transport from the fractures to the matrix is modeled as a wall boundary condition for the fractures, which, in turn, is predicted by recurrent neural networks using the concentration histories at the fracture-matrix boundary. Specifically, we consider a meter-scale fracture network embedded in sandstones, where the smallest feature is at the micron scale. Here the proposed framework allows us to span five orders of magnitude in length scales by capturing mineral precipitation in the altered layer of the rock matrix at the pore scale across the entire meter-scale fracture network.

42 ENGINEERING↗
DOE OSTI2021

Striving to translate shale physics across ten orders of magnitude: What have we learned?

Shales will play an important role in the successful transition of energy from fossil-based resources to renewables in the coming decades. Aside from being a significant source of low-carbon intensity fuels, like natural gas, they also serve as geologic seals of subsurface formations that may be used to isolate nuclear waste, sequester CO 2 , or store intermittent energy (e.g., solar hydrogen). Despite their importance, shales pose significant engineering and environmental challenges due to their nanoporous structure and extreme heterogeneity that spans at least ~10 orders of magnitude in spatial scale. Two challenges inhibit a system-level understanding: (1) the physics of fluid flow and phase behavior in shales are poorly understood due to the dominant molecular interactions between minerals and fluids under confinement, and (2) the apparent lack of scale separation that prevents a reliable (closed) description of the physics at any single scale of observation. In this review, we focus on the latter issue and discuss scale translation, which in its broadest sense is transforming data or simulations from one spatiotemporal scale to another. While effective scale translation is not exclusive to shales, but all geologic porous media, the need for it is especially acute in shales given their high degree of heterogeneity. Classical theories like homogenization, while indispensable, fail when scales are not separated. Other methods, like numerical upscaling, scale-translate in only one direction: small to large, but not the reverse, called downscaling. However, the confluence of advances in three areas are bringing challenging problems such as shales within reach: increased computational power and scalable algorithms; high-resolution imaging and multi-modal data acquisition; and machine learning to process massive amounts of data. While these advances equip geoscientists with a wide array of experimental and computational tools, no individual tool can probe the entire gamut of heterogeneity in shales. Their effective use, therefore, requires an ability to bridge between various data types obtained at different scales. The aim of this review is to present a coherent account of computational and experimental methods that may be used to achieve just that, i.e., to perform scale translation. We provide a broader definition of scale translation, one that transcends classical homogenization and upscaling methods, but is consistent with them and accommodates notions like downscaling and data translation. After a brief introduction to homogenization, we review hybrid methods, numerical upscaling and its recent extensions, multiscale computing, high-resolution imaging, and machine learning. We place particular emphasis on multiscale computing and propose an algorithmic framework to bridge between the pore (micro) and Darcy (macro) scales. Throughout the paper, we draw comparisons between the various methods and highlight their (often hidden) similarities, differences, benefits, and pitfalls. We finally conclude with two case studies on shales that exemplify some of the methods presented.

58 GEOSCIENCES↗
DOE OSTI2023

Learning macroscopic internal variables and history dependence from microscopic models

This paper concerns the study of history dependent phenomena in heterogeneous materials in a two-scale setting where the material is specified at a fine microscopic scale of heterogeneities that is much smaller than the coarse macroscopic scale of application. Here, we specifically study a polycrystalline medium where each grain is governed by crystal plasticity while the solid is subjected to macroscopic dynamic loads. The theory of homogenization allows us to solve the macroscale problem directly with a constitutive relation that is defined implicitly by the solution of the microscale problem. However, the homogenization leads to a highly complex history dependence at the macroscale, one that can be quite different from that at the microscale. In this paper, we examine the use of machine-learning, and especially deep neural networks, to harness data generated by repeatedly solving the finer scale model to: (i) gain insights into the history dependence and the macroscopic internal variables that govern the overall response; and (ii) to create a computationally efficient surrogate of its solution operator, that can directly be used at the coarser scale with no further modeling. We do so by introducing a recurrent neural operator (RNO), and show that: (i) the architecture and the learned internal variables can provide insight into the physics of the macroscopic problem; and (ii) that the RNO can provide multiscale, specifically FE 2 , accuracy at a cost comparable to a conventional empirical constitutive relation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Calculation of Hydraulic Conductivity for a Porous Medium with Staggered Array of Solid Grains - 20087

For a porous medium with periodic array of solid grains. the flow of an incompressible viscous fluid in the pore space is considered with an aim of calculating the permeability of the medium. The solid grains are of ellipsoidal shape. For this purpose, a unit-cell boundary-value problem (Stokes problem) is solved numerically. With the assumption of periodicity of the medium structure and the fluid properties (in fact constant), the flow field and pressure distributions are determined numerically with inertial effects ignored. By applying the theory of homogenization to the viscous flow through pore space, a microscale boundary-value problem is defined. After solving the problem, the permeability is determined by taking the micro-cell average of the velocity components, Various porous structures are chosen to allow a wide range of porosity values. The permeability increases mildly for porosity values in the lower range. But it increases sharply for porosity values in the upper range of the porosity. The permeability of the present study agrees with the Kozeny-Carman relation excellently for smaller porosity. Certain discrepancy between the calculated permeability in this study and the Kozeny-Carman relation appears and increases with porosity for larger porosity. It is believed that the discrepancy is due to the replacement of the pores by a bundle of capillaries, which is unrealistic, in the Kozeny-Carman relation. (authors)

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

Analysis of multi-electron, multi-step homogeneous catalysis by rotating disc electrode voltammetry: theory, application, and obstacles

Rotating disc electrode (RDE) voltammetry has been widely adopted for the study of heterogenized molecular electrocatalysts for multi-step fuel-forming reactions but this tool has never been comprehensively applied to their homogeneous analogues. Here, the utility and limitations of RDE techniques for mechanistic and kinetic analysis of homogeneous molecular catalysts that mediate multi-electron, multi-substrate redox transformations are explored. Using the ECEC' reaction mechanism as a case study, two theoretical models are derived based on the Nernst diffusion layer model and the Hale transformation. Current–potential curves generated by these computational strategies are compared under a variety of limiting conditions to identify conditions under which the more minimalist Nernst Diffusion Layer approach can be applied. Based on this theoretical treatment, strategies for extracting kinetic information from the plateau current and the foot of the catalytic wave are derived. RDEV is applied to a cobaloxime hydrogen evolution reaction (HER) catalyst under non-aqueous conditions in order to experimentally validate this theoretical framework and explore the feasibility of RDE as a tool for studying homogeneous catalysts. Importantly, analysis of the foot-of-the-wave via this theoretical framework provides rate constants for elementary reaction steps that agree with those extracted from stationary voltammetric methods, supporting the application of RDE to study homogeneous fuel-forming catalysts. Lastly, obstacles encountered during the kinetic analysis of cobaloxime, along with the voltammetric signatures used to diagnose this reactivity, are discussed with the goal of guiding groups working to improve RDE set-ups and help researchers avoid misinterpretation of RDE data.

08 HYDROGEN↗
DOE OSTI2022

Nucleation of Grain Boundary Phases

Here, we derive a theory that describes homogeneous nucleation of grain boundary (GB) phases. Our analysis takes account of the energy resulting from the GB phase junction, the line defect separating two different GB structures, which is necessarily a dislocation as well as an elastic line force due to the jump in GB stresses. The theory provides analytic forms for the elastic interactions and the core energy of the GB phase junction that, along with the change in GB energy, determines the nucleation barrier. We apply the resulting nucleation model to simulations of GB phase transformations in tungsten. Our theory explains why under certain conditions GBs cannot spontaneously change their structure even to a lower energy state.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Effect of Carbon on Void Nucleation in Iron

The study reports the significance of carbon presence in affecting void nucleation in Fe. Without carbon, void nucleation rates decrease gradually at high temperatures but remain significantly high and almost saturated at low temperatures. With carbon present, even at 1 atomic parts per million, void nucleation rates show a low-temperature cutoff. With higher carbon levels, the nucleation temperature window becomes narrower, the maximum nucleation rate becomes lower, and the temperature of maximum void nucleation shifts to a higher temperature. Fundamentally, this is caused by the change in effective vacancy diffusivity due to the formation of carbon-vacancy complexes. The high sensitivity of void nucleation to carbon comes from the high sensitivity of void nucleation to the vacancy arrival rate in a void. The void nucleation is calculated by first obtaining the effective vacancy diffusivity considering the carbon effect, then calculating the defect concentration and defect flux change considering both carbon effects and pre-existing dislocations, and finally calculating the void nucleation rate based on the recently corrected homogeneous void nucleation theory. The study is important not only in the fundamental understanding of impurity effects in ion/neutron irradiation but also in alloy engineering for judiciously introducing impurities to increase swelling resistance, as well as in the development of simulation and modeling methodologies applicable to other metals.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A cavitation and dynamic void growth model for a general class of strain-softening amorphous materials

In this work, we derive a homogenized framework for cavitation and dynamic void growth in a general class of strain-softening materials with particular emphasis on amorphous materials. Analytic solutions are derived by considering an assemblage of spherical shells subject to dynamic hydrostatic tensile pressure. The isochoric response of the matrix material is captured by a free volume theory based amorphous viscoplasticity model. The framework accounts for the coupled rate-dependent, strain-softening behavior typically exhibited in amorphous materials. Validation of the homogenized theory is carried out against numerous direct numerical simulations with excellent agreement. Furthermore, we find that a reduction in cavitation strength and significant acceleration of dynamic void growth rates is induced as a consequence of strain-softening in the amorphous matrix material. Lastly, the model is utilized to interpret some non-intuitive experimental observations of amorphous materials subject to shock compression and subsequent spall failure.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Orbital-selective metal skin induced by alkali-metal-dosing Mott-insulating Ca2RuO4

Abstract Doped Mott insulators are the starting point for interesting physics such as high temperature superconductivity and quantum spin liquids. For multi-band Mott insulators, orbital selective ground states have been envisioned. However, orbital selective metals and Mott insulators have been difficult to realize experimentally. Here we demonstrate by photoemission spectroscopy how Ca 2 RuO 4 , upon alkali-metal surface doping, develops a single-band metal skin. Our dynamical mean field theory calculations reveal that homogeneous electron doping of Ca 2 RuO 4 results in a multi-band metal. All together, our results provide evidence for an orbital-selective Mott insulator breakdown, which is unachievable via simple electron doping. Supported by a cluster model and cluster perturbation theory calculations, we demonstrate a type of skin metal-insulator transition induced by surface dopants that orbital-selectively hybridize with the bulk Mott state and in turn produce coherent in-gap states.

Physics↗
DOE OSTI2023

On the Cauchy problem for Boltzmann equation modeling a polyatomic gas

In the present article, we consider the Boltzmann equation that models a polyatomic gas by introducing one additional continuous variable, referred to as microscopic internal energy. We establish existence and uniqueness theory in the space homogeneous setting for the full non-linear case, under an extended Grad-type assumption on transition probability rates, which comprises hard potentials for both the relative speed and internal energy with the rate in the interval 0,2, multiplied by an integrable angular part and integrable partition functions. The Cauchy problem is resolved by means of an abstract ordinary differential equation (ODE) theory in Banach spaces for the initial data with finite and strictly positive gas mass and energy, finite momentum, and additionally finite K* polynomial moment, with K* depending on the rate of the transition probability and the structure of a polyatomic molecule or its internal degrees of freedom. Moreover, we prove that polynomially and exponentially weighted Banach space norms associated with the solution are both generated and propagated uniformly in time.

Physics↗
DOE OSTI2021

Machine learning dielectric screening for the simulation of excited state properties of molecules and materials

Accurate and efficient calculations of absorption spectra of molecules and materials are essential for the understanding and rational design of broad classes of systems. Solving the Bethe–Salpeter equation (BSE) for electron–hole pairs usually yields accurate predictions of absorption spectra, but it is computationally expensive, especially if thermal averages of spectra computed for multiple configurations are required. We present a method based on machine learning to evaluate a key quantity entering the definition of absorption spectra: the dielectric screening. We show that our approach yields a model for the screening that is transferable between multiple configurations sampled during first principles molecular dynamics simulations; hence it leads to a substantial improvement in the efficiency of calculations of finite temperature spectra. We obtained computational gains of one to two orders of magnitude for systems with 50 to 500 atoms, including liquids, solids, nanostructures, and solid/liquid interfaces. Importantly, the models of dielectric screening derived here may be used not only in the solution of the BSE but also in developing functionals for time-dependent density functional theory (TDDFT) calculations of homogeneous and heterogeneous systems. Overall, our work provides a strategy to combine machine learning with electronic structure calculations to accelerate first principles simulations of excited-state properties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The Cauchy Problem for Boltzmann Bi-linear Systems: The Mixing of Monatomic and Polyatomic Gases

Abstract From a unified vision of vector valued solutions in weighted Banach spaces, this paper establishes the existence and uniqueness for space homogeneous Boltzmann bi-linear systems with conservative collisional forms arising in complex gas dynamical structures. This broader vision is directly applied to dilute multi-component gas mixtures composed of both monatomic and polyatomic gases. Such models can be viewed as extensions of scalar Boltzmann binary elastic flows, as much as monatomic gas mixtures with disparate masses and single polyatomic gases, providing a unified approach for vector valued solutions in weighted Banach spaces. Novel aspects of this work include developing the extension of a general ODE theory in vector valued weighted Banach spaces, precise lower bounds for the collision frequency in terms of the weighted Banach norm, energy identities, angular or compact manifold averaging lemmas which provide coerciveness resulting into global in time stability, a new combinatorics estimate for p -binomial forms producing sharper estimates for the k -moments of bi-linear collisional forms. These techniques enable the Cauchy problem improvement that resolves the model with initial data corresponding to strictly positive and bounded initial vector valued mass and total energy, in addition to only a $$2^+$$ 2 + moment determined by the hard potential rates discrepancy, a result comparable in generality to the classical Cauchy theory of the scalar homogeneous Boltzmann equation.

Physics↗