Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical homogenization”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Supported Metal Pair-Site Catalysts

Complexes with neighboring metal centers in isolated pairs and their analogues on surfaces are drawing increasing attention as catalysts. These include molecular homogeneous catalysts incorporating various ligands, enzymes, and solids that include pairs of metal atoms mounted on supports. Catalysts in this broad class are active for numerous reactions and offer unexplored opportunities to address challenging reactions, such as oxidation of methane and oxidation of water in artificial photosynthesis. The subject of supported metal pair-site catalysts is in its infancy, facing challenges in (a) precise synthesis, (b) structure determination at the atomic scale, and (c) stabilization in reactive atmospheres. In this Perspective, we summarize key characteristics of molecular and enzymatic catalysts that incorporate neighboring metal centers and build on this foundation to assess the emerging literature of metal pair-site catalysts on various supports. The supported catalysts include those synthesized by anchoring molecular dinuclear precursors to support surfaces and those synthesized by selective formation of dinuclear surface species from mononuclear surface species. Examples of metals in this class are rhodium and iridium, and examples of supports are MgO and Fe 2 O 3 . We summarize characterization of these materials by electron microscopy and spectroscopy, emphasizing atomic-resolution aberration-corrected scanning transmission electron microscopy and spectroscopies that provide atomic-scale structural information and allow characterization of functioning catalysts, especially X-ray absorption spectroscopy. Finally, we list some opportunities for research, including suggestions that might lead to structurally well-defined supported metal pair-sites with new catalytic properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Quantum Drude oscillators coupled with Coulomb potential as an efficient model for bonded and non-covalent interactions in atomic dimers

The quantum Drude oscillator (QDO) model has been widely used as an efficient surrogate to describe the electric response properties of matter as well as long-range interactions in molecules and materials. Most commonly, QDOs are coupled within the dipole approximation so that the Hamiltonian can be exactly diagonalized, which forms the basis for the many-body dispersion method [Phys. Rev. Lett. 108, 236402 (2012)]. The dipole coupling is efficient and allows us to study non-covalent many-body effects in systems with thousands of atoms. However, there are two limitations: (i) the need to regularize the interaction at short distances with empirical damping functions and (ii) the lack of multipolar effects in the coupling potential. In this work, we convincingly address both limitations of the dipole-coupled QDO model by presenting a numerically exact solution of the Coulomb-coupled QDO model by means of quantum Monte Carlo methods. We calculate the potential-energy surfaces of homogeneous QDO dimers, analyzing their properties as a function of the three tunable parameters: frequency, reduced mass, and charge. We study the coupled-QDO model behavior at short distances and show how to parameterize this model to enable an effective description of chemical bonds, such as the covalent bond in the H2 molecule.

Chemistry↗
DOE OSTI2023

Development and evaluation of a skeletal mechanism for EHN additized gasoline mixtures in large Eddy simulations of HCCI combustion

Advanced Low Temperature Combustion modes, such as the Sandia proposed Additive-Mixing Fuel Injection (AMFI), can unlock significant potential to boost fuel conversion efficiency and ultimately improve the energy conversion of internal combustion engines. This is a novel improved combustion process that is enabled by supplying small (<5%) variable amounts of autoignition improver to the fuel to enhance the engine operation and control. Common, diesel-fuel ignition-quality enhancing additive, 2-ethylexyl nitrate (EHN), is doped into gasoline to enable Sandia LTGC + AMFI combustion. This manuscript focuses on the development of a reduced sub-mechanism for EHN chemical kinetics at engine relevant conditions that is implemented into a skeletal mechanism for chemical kinetic studies of gasoline surrogate fuels. The mechanism validation utilized zero-dimensional numerical simulations and comparison to shock tube ignition-delay data of pure and EHN-doped n-heptane. Additional validation is presented with Homogeneous Charge Compression-Ignition (HCCI) engine data of pure and EHN-doped research-grade E10 gasoline. Then, the mechanism was deployed in a 3-D computational fluid dynamics (CFD) using Large Eddy Simulations (LES) to model the HCCI engine experiments of 0.4% vol EHN additized E10 gasoline at several equivalence ratios. Simulations showed a very good performance of the mechanism, and the model accurately reproduced (a) the ignition point, (b) combustion phasing, (c) combustion duration, and (d) the peak of the heat release rates of the engine experiments. The results show that EHN promotes Low-Temperature Heat Release, ultimately driving the gasoline to autoignite at thermodynamic conditions where the fuel would not otherwise ignite. Overall, this work demonstrates a viable reduced chemical-kinetic mechanism for EHN and shows that it can be combined with a skeletal gasoline mechanism for CFD-LES analysis of well-mixed LTGC that matches well with experimental results. The CFD-LES analysis also shows the spatial distribution of EHN-fuel interactions that control the autoignition throughout the combustion chamber.

Guleria, Gaurav↗
DOE OSTI2025

Experimental Investigation of Vapor Formation in Liquid CO2 Flow Through a Converging-Diverging Nozzle

Carbon dioxide is an attractive working fluid for many cycles, including for pumped thermal energy storage (PTES). A challenge with some proposed sCO2 PTES cycles is the operation of sCO2 machinery outside the typical bounds of experience, with local phase change from the liquid state being particularly unknown. Presently, there is insufficient data in the literature regarding multiphase CO2 to adequately design a multiphase-tolerant turbine, so generation of foundational data is required. This experimental study investigates the flow characteristics of sub-sonic liquid CO2 undergoing expansion and phase change in a converging-diverging nozzle. The nozzle is instrumented to measure static pressure, unsteady pressure, temperature, and density. The static pressure transducers are located at 27 axial locations to accurately characterize the pressure profile in the nozzle. High-accuracy RTDs are located at the entrance and exit of the nozzle, and three dynamic pressure transducers are strategically located to capture any unsteady phenomena. During testing, values of mass flow and nozzle inlet pressure are swept to vary the pressure drop and fluid properties. The measured total pressure drop in the nozzle is compared to a homogenous model and the Lockhart-Martinelli correlation method, with the latter predicting loss quite closely. The resulting data set is valuable for validating multiphase numerical models in a simple geometry before implementation of these models in turbomachinery design.

25 ENERGY STORAGE↗
DOE OSTI2020

Observation of Road Salt Aerosol Driving Inland Wintertime Atmospheric Chlorine Chemistry

Inland sources of particulate chloride for atmospheric nitryl chloride (ClNO2) formation remain unknown and unquantified, hindering air quality assessments. Globally each winter, tens of millions of tons of road salt are spread on roadways for deicing. Here, we identify road salt aerosol as the primary chloride aerosol source, accounting for 80–100% of ClNO 2 formation, at an inland urban area in the wintertime. This study provides experimental evidence of the connection between road salt and air quality through the production of this important reservoir for nitrogen oxides and chlorine radicals, which significantly impact atmospheric composition and pollutant fates. A numerical model was employed to quantify the contributions of chloride sources to ClNO 2 production. The traditional method for simulating ClNO 2 considers chloride to be homogeneously distributed across the atmospheric particle population; yet, we show that only a fraction of the particulate surface area contains chloride. Our new single-particle parametrization considers this heterogeneity, dramatically lowering overestimations of ClNO 2 levels that have been routinely reported using the prevailing methods. The identification of road salt as a ClNO 2 source links this common deicing practice to atmospheric composition and air quality in the urban wintertime environment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Impact Modeling and Testing of Pouch and Prismatic Cells

Understanding battery response under impact is critical to improve the safety of electrified vehicles. Nevertheless, predicting the impact behaviors of batteries is not straightforward since a battery cell usually contains hundreds of thin layers with dramatically different material properties and multiple physical processes occur simultaneously during cell deformation. Here we utilized both empirical tests and numerical models to capture the failure process of pouch and prismatic cells in various impact scenarios. In each test, a cell was hit once by an indenter dropped from a certain height. During which the cell penetration, loading force, voltage and temperature were monitored to characterize the cell’s response. Meanwhile, numerical models were developed to capture the coupled mechanical, electrical, electrochemical and thermal responses of batteries. In these models, the cell bulk was treated as a homogeneous part to achieve computational efficiency required by large-scale simulations, and it was represented by the geologic cap model that allows both shear and compaction deformation. Simulation results showed agreement with experimental data in essential features of cell behaviors during impact. Details of the test setup, model development and cell failure behaviors are presented in this paper. Additionally, capabilities, limitations and future improvement of the battery safety modeling are discussed.

25 ENERGY STORAGE↗
DOE OSTI2021

Transport in Stochastic Media with Random Chord Length Distributions

Thermal radiation transport computations in binary Markovian random mixtures rely almost exclusively on the Levermore-Pomraning (LP) model which is obtained by applying a heuristic closure to the ensemble averaged random medium transport equation. The validity of this model has been extensively tested by comparing numerical results over a broad parameter range (material types and mixing parameters) against benchmark solutions in planar geometry. The conditions under which the LP-model provides useful results and when it breaks down are now well established, but work to date has been largely restricted to homogeneous mixing statistics, i.e., the mean chord lengths of both materials are taken to be spatially constant. In recent work, this limitation was relaxed by allowing the mean chord lengths and, in a consistent fashion, the volume fractions in the LP-model to vary spatially and in a follow-up investigation benchmark solutions were obtained by ensemble averaging results over material realizations sampled from a nonhomogeneous Poisson process (NHPP). Numerical experiments in rod geometry with specifically linear and quadratic spatial dependence of chord lengths showed that the material averaged radiation intensities vary nonmonotonically with depth into the medium, in stark contrast to solutions obtained assuming uniform chord lengths. Moreover, depending on the local optical thickness and strength of scattering, the LP-model results showed locally more nuanced deviations from the benchmark solutions than was the case with constant chord lengths. These limited numerical investigations highlight the nontrivial qualitative and quantitative consequences of nonhomogeneous mixing statistics, in particular that closure approximations may not be uniformly valid or invalid over the problem domain.

42 ENGINEERING↗
DOE OSTI2020

Predicting mesoscale spectral thermal conductivity using advanced deterministic phonon transport techniques

We present a review, demonstration and simulation of phonon transport for the purposes of predicting materials performance at the mesoscale. We focus primarily on the development and implementation of a unified methodology to enable predictive heat transport. We report on the current state of the art as it pertains to deterministic phonon transport methodologies, discussing various topics concerning phonons. In application, we focus on the Self-Adjoint Angular Flux (SAAF) formulation of the Boltzmann transport equation for phonons, and develop the spatial, angular, and material property discretization required to accurately simulate the predictive physics of heat transport in dielectrics. We discuss thermal interfacial resistance and present our formulation of the diffuse mismatch model for simulating phonon interactions at internal boundaries. We have recently developed a deterministic, spectral phonon transport method for predicting effective thermal conductivity ($\kappa_{\textrm{eff}}$), using Bose-Einstein source terms coupled through an average material temperature. In this method, we introduce a closure term to the phonon transport system which acts as a redistribution function for the total energy of the system, and serves as a guide for the amount of non-equilibrium behavior occurring in the system. This method predicts thermal conductivity and equilibrium temperature distributions in homogeneous and heterogeneous materials using data generated by ab initio density functional theory methods. We employ polarization, density of states and full dispersion spectra to resolve thermal conductivity with numerous angular and spatial discretizations. Our implementation utilizes a Richardson iteration on a modified version of the phonon scattering source. The equations associated with this method are solved via a modification of traditional source iteration. We compare the performance of source iteration applied to an existing uncoupled, traditional SAAF method to our new method and comment on the iterative performance of each. We observe ballistic and diffusive phonon scattering as acoustic thickness of the domain changes, and are able to make comparisons between the accuracy and efficiency of both methods.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fabrication, thermal analysis, and heavy ion irradiation resistance of epoxy matrix nanocomposites loaded with silane-functionalized ceria nanoparticles

This paper describes a detailed understanding of how nanofillers function as radiation barriers within the polymer matrix, and how their effectiveness is impacted by factors such as composition, size, loading, surface chemistry, and dispersion. Here, we designed a comprehensive investigation of heavy ion irradiation resistance in epoxy matrix composites loaded with surface-modified ceria nanofillers, utilizing tandem computational and experimental methods to elucidate radiolytic damage processes and relate them to chemical and structural changes observed through thermal analysis, vibrational spectroscopy, and electron microscopy. A detailed mechanistic examination supported by FTIR spectroscopy data identified the bisphenol A moiety as a primary target for degradation reactions. Results of computational modeling by the Stopping Range of Ions in Matter (SRIM) Monte Carlo simulation were in good agreement with damage analysis from surface and cross-sectional SEM imaging. All metrics indicated that ceria nanofillers reduce the damage area in polymer nanocomposites, and that nanofiller loading and homogeneity of dispersion are key to effective damage prevention. The results of this study represent a significant pathway for engineered irradiation tolerance in a diverse array of polymer nanocomposite materials. Numerous areas of materials science can benefit from utilizing this facile and effective method to extend the reliability of polymer materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

On the numerical accuracy in finite-volume methods to accurately capture turbulence in compressible flows

The goal of the present article is to understand the impact of numerical schemes for the reconstruction of data at cell faces in finite-volume methods, and to assess their interaction with the quadrature rule used to compute the average over the cell volume. Here, third-, fifth- and seventh-order WENO-Z schemes are investigated. On a problem with a smooth solution, the theoretical order of convergence rate for each method is retrieved, and changing the order of the reconstruction at cell faces does not impact the results, whereas for a shock-driven problem all the methods collapse to first-order. Here, study of the decay of compressible homogeneous isotropic turbulence reveals that using a high-order quadrature rule to compute the average over a finite-volume cell does not improve the spectral accuracy and that all methods present a second-order convergence rate. However the choice of the numerical method to reconstruct data at cell faces is found to be critical to correctly capture turbulent spectra. In the context of simulations with finite-volume methods of practical flows encountered in engineering applications, it becomes apparent that an efficient strategy is to perform the average integration with a low-order quadrature rule on a fine mesh resolution, whereas high-order schemes should be used to reconstruct data at cell faces.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Partially polaron-transformed quantum master equation for exciton and charge transport dynamics

Polaron-transformed quantum master equation (PQME) offers a unified framework to describe the dynamics of quantum systems in both limits of weak and strong couplings to environmental degrees of freedom. Thus, the PQME serves as an efficient method to describe charge and exciton transfer/transport dynamics for a broad range of parameters in condensed or complex environments. However, in some cases, the polaron transformation (PT) being employed in the formulation invokes an over-relaxation of slow modes and results in premature suppression of important coherence terms. A formal framework to address this issue is developed in the present work by employing a partial PT that has smaller weights for low frequency bath modes. In this work, it is shown that a closed form expression of a second order time-local PQME including all the inhomogeneous terms can be derived for a general form of partial PT, although more complicated than that for the full PT. All the expressions needed for numerical calculation are derived in detail. Applications to a model of a two-level system coupled to a bath of harmonic oscillators, with test calculations focused on those due to homogeneous relaxation terms, demonstrate the feasibility and the utility of the present approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Fierro Version 2.x

FIERRO is a parallel C++ code designed to simulate fluid mechanics, heat transfer, and solid mechanics in two- and three-dimensional space. FIERRO is written to run on homogeneous (CPU) and heterogeneous (CPU+GPU) high performance computing machines. Fierro can aid a) modeling and design efforts that have historically relied on commercial implicit and explicit finite element codes, b) numerical methods research, c) manufacturing research, and d) computer science research. The code contains diverse numerical methods to solve the governing physics equations for both quasi-static and dynamic problems. Mathematical optimization solvers are coupled to the numerical methods to research topology and shape optimization that has application to additive manufacturing, and to create novel numerical approaches. Phase-field methods with micromechanical solvers are provided to simulate microstructure formation and evolution in manufacturing processes. The micromechanical solvers can also help research efforts create continuum-scale constitutive models for solids, as a function of the microstructure, in situ in a calculation or in a stand-alone manner. No physical data exists within the code.

Morgan, Nathaniel↗
DOE OSTI2023

Fierro

FIERRO is a parallel C++ code designed to simulate fluid mechanics, heat transfer, and solid mechanics in two- and three dimensional space. FIERRO is written to run on homogeneous (CPU) and heterogeneous (CPU+GPU) high performance computing machines. Fierro can aid a) modeling and design efforts that have historically relied on commercial implicit and explicit finite element codes, b) numerical methods research, c) manufacturing research, and d) computer science research. The code contains diverse numerical methods to solve the governing physics equations for both quasi-static and dynamic problems. Mathematical optimization solvers are coupled to the numerical methods to research topology and shape optimization that has application to additive manufacturing, and to create novel numerical approaches. Phase-field methods with micromechanical solvers are provided to simulate microstructure formation and evolution in manufacturing processes. The micromechanical solvers can also help research efforts create continuum-scale constitutive models for solids, as a function of the microstructure, in situ in a calculation or in a stand-alone manner. No physical data exists within the code.

Morgan, Nathaniel↗
DOE OSTI2021

Numerical and experimental evaluation of mechanical performance of the multifunctional energy storage composites

This work presents numerical simulation methods to model the mechanical behavior of the multifunctional energy storage composites (MESCs), which consist of a stack of multiple thin battery layers reinforced with through-the-hole polymer rivets and embedded inside carbon fiber composite laminates. MESC has been demonstrated through earlier experiments on its exceptional behavior as a structural element as well as a battery. However, the inherent complex infrastructure of the MESC design has created significant challenges in simulation and modeling. A novel homogenization technique was adopted to characterize the multi-layer properties of battery material using physics-based constitutive equations combined with nonlinear deformation theories to handle the interface between the battery layers. Second, mechanical damage and failure modes among battery materials, polymer reinforcements, and carbon fiber-polymer interfaces were characterized through appropriate models and experiments. The model of MESCs has been implemented in a commercial finite element code in ABAQUS. A comparison of structural response and failure modes from numerical simulations and experimental tests are presented. The results of the study showed that the predictions of elastic and damage responses of MESCs at various loading conditions agreed well with the experimental data. © 2021

Materials Science↗
DOE OSTI2020

Modeling of deformation of battery cells using thick shell element formulation

Here, we describe a new approach for modeling nonlinear deformation and stress distribution of battery cells using a new thick shell finite element formulation with a through-thickness calculation of stresses and strains that satisfy equilibrium conditions. Battery cells are transversely layered materials that contain numerous thin layers in a repeating sequence. The layers are made of materials with significantly different mechanical responses. Explicit discretization of the layers is computationally impractical except for very small domains, while homogenized material models cannot account for stress and strain variations and partition through the thickness. The new formulation allows for calculation of the transverse and interlayer stresses based on the material properties of the individual cell layers. The application to problems where the transverse stresses have strong influence is also possible.

25 ENERGY STORAGE↗
DOE OSTI2024

Quantifying local and global mass balance errors in physics-informed neural networks

Physics-informed neural networks (PINN) have recently become attractive for solving partial differential equations (PDEs) that describe physics laws. By including PDE-based loss functions, physics laws such as mass balance are enforced softly in PINN. This paper investigates how mass balance constraints are satisfied when PINN is used to solve the resulting PDEs. We investigate PINN’s ability to solve the 1D saturated groundwater flow equations (diffusion equations) for homogeneous and heterogeneous media and evaluate the local and global mass balance errors. We compare the obtained PINN’s solution and associated mass balance errors against a two-point finite volume numerical method and the corresponding analytical solution. We also evaluate the accuracy of PINN in solving the 1D saturated groundwater flow equation with and without incorporating hydraulic heads as training data. We demonstrate that PINN’s local and global mass balance errors are significant compared to the finite volume approach. Tuning the PINN’s hyperparameters, such as the number of collocation points, training data, hidden layers, nodes, epochs, and learning rate, did not improve the solution accuracy or the mass balance errors compared to the finite volume solution. Mass balance errors could considerably challenge the utility of PINN in applications where ensuring compliance with physical and mathematical properties is crucial.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Kelvin–Helmholtz instability under stabilizing parallel magnetic field in nonhomogeneous compressible MHD flows

We study the Kelvin–Helmholtz instability (KHI) for the general case of a compressible, nonhomogeneous, magnetized plasma flow. The study is limited to a vortex sheet interface with an imposed parallel magnetic field. We introduce a new formalism based on a convective Mach number M c , a convective Alfvénic Mach number M Ac , and a total convective Mach number that combines the two. We derive an analytic expression of the KHI growth rate for a homogeneous flow (i.e., zero Atwood number, A=0) that converges toward both the expression for unmagnetized compressible flow and Chandrasekhar's expression for magnetized incompressible flow. Otherwise, the dispersion relation is solved numerically and allows deriving general stability diagrams of magnetized KHI for the triplet (A, M c , β −plasma) parameters. We show these parameters uniquely define all configurations for a parallel magnetic field. We also construct diagrams with respect to the convective Alfvénic Mach number, the β − plasma parameter, or the magnetic field showing which magnetic field strength is required for stabilizing a given shear flow. The theoretical growth rates are compared with 18 simulations made with the GAMERA code, currently used for 3D magnetospheric simulations. Finally, we apply our results to the analysis of a past KHI experiment performed at the OMEGA laser facility, showing linear theory succeeds to provide accurate estimates of the growth rate at early times. We further discuss how our results can inform future experiments in the high-Mach magnetized regime at the National Ignition Facility. Possible limitations of the study due to resistive, mixing, or turbulence effects are discussed.

compressible flows↗
DOE OSTI2023

Synthetically non-Hermitian nonlinear wave-like behavior in a topological mechanical metamaterial

Topological mechanical metamaterials have enabled new ways to control stress and deformation propagation. Exemplified by Maxwell lattices, they have been studied extensively using a linearized formalism. Herein, we study a two-dimensional topological Maxwell lattice by exploring its large deformation quasi-static response using geometric numerical simulations and experiments. We observe spatial nonlinear wave-like phenomena such as harmonic generation, localized domain switching, amplification-enhanced frequency conversion, and solitary waves. We further map our linearized, homogenized system to a non-Hermitian, nonreciprocal, one-dimensional wave equation, revealing an equivalence between the deformation fields of two-dimensional topological Maxwell lattices and nonlinear dynamical phenomena in one-dimensional active systems. Our study opens a regime for topological mechanical metamaterials and expands their application potential in areas including adaptive and smart materials and mechanical logic, wherein concepts from nonlinear dynamics may be used to create intricate, tailored spatial deformation and stress fields greatly transcending conventional elasticity.

36 MATERIALS SCIENCE↗