Search NASA⌕ Search

SEARCH · Search NASA

Results for “effective medium approximations”

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 OSTI2021

Conduction Heat Transfer through Solid in Porous Materials: A Comparative Study by Finite-Element Simulations and Effective Medium Approximations

The development of insulation materials with low effective thermal conductivity is essential for energy savings in various applications, including buildings, food services, pipe insulation, and refrigeration. Such materials can be developed by using micro- or nanoporous structures, as well as low-thermal conductivity gases and/or reduced pressure. A variety of effective medium approximation (EMA) models have been developed to study thermal transport through solid in porous structures. However, in many cases, the impacts of porosity on solid conduction are not well predicted by EMA models because of the assumptions made for simplification. Furthermore, the results can vary by a factor of up to 1.5, depending on the morphology of the pores. Hence, proper guidance is needed to choose the appropriate EMA model for a given morphology. This work presents a finite element method study using COMSOL Multiphysics software for various geometries, including hollow cubes, overlapping hollow spheres (normally stacked and tightly packed), and Voronoi structures, and compares the results to those obtained from various EMA models. Simulation results for the hollow cubes and Voronoi structures matched well with the Maxwell-Eucken and Russell model. The results for hollow spheres closely matched the Bauer model and the Glicksman model with fitting parameters. This work provides guidance on analyzing and designing insulation materials for energy savings in the future.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Effective medium approximation for the refractive index of stratified metal oxide composites synthesized by atomic layer deposition

Atomic layer deposition (ALD) is a unique method for synthesizing conformal layers with precise composition. It is especially useful for the synthesis of mixed metal oxides for functional materials. One application of interest is the use of ALD for tailoring the refractive index of coatings. In homogeneously distributed composites of metal oxides, the refractive properties can be approximated as the average of the indices of the components. This is known as an effective medium approximation (EMA) and can be used to design the macroscopic properties of composites. ALD produces layered, anisotropic films, so the validity of an EMA in describing these films is not clear. Here, we use optical simulations and experimental characterization of stratified composites of TiO 2 and Al 2 O 3 to study the application of an EMA to the transmission and reflection behavior of ALD-prepared mixed metal oxide thin films. We found that when the characteristic layer thickness is smaller than roughly 10 nm, the optical spectra of theoretical and experimental ALD films match the equivalent theoretical spectra for a material with a refractive index calculated from a simple, compositionally weighted EMA. Hence, this EMA could describe the transmission and reflectance spectra in ALD-derived TiO2–Al 2 O 3 nanolaminates when the films were sufficiently stratified even without thorough characterization of the real layers in the films. As a result, we demonstrated that ALD can be used to prepare effectively homogenous mixed metal oxide films with a predictable, tailorable refractive index of any value between those of the TiO 2 and Al 2 O 3 component materials.

42 ENGINEERING↗
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 OSTI2023

Computational design of isotropic and anisotropic ultralow thermal conductivity polymer foams

Current state-of-the-art commercial polymer thermal insulation foam exhibits a thermal conductivity of 24 mW·m -1 ·K -1 (equivalently thermal resistivity of R-6/in.), similar to that of static air. To further optimize building energy efficiency, achieving even lower thermal conductivity is needed, which is, however, highly challenging. This paper presents computational evidence that demonstrates the feasibility of achieving an ultra-low thermal conductivity of less than 14.4 mW·m -1 ·K -1 (equivalently R-10/in.) using isotropic and anisotropic foam cell designs. For the isotropic design, we have identified analytical effective medium approximation (EMA) models within the accuracy of ±5% as finite element analysis (FEA) in predicting the effective thermal conductivity of foams with various porosities and filler gases. For the anisotropic design, we have developed and validated new EMA models against FEA in predicting the effective thermal conductivity of general anisotropic cuboids and Voronoi foams. For both isotropic and anisotropic designs, the design spaces for 18, 16, and 14.4 mW·m -1 ·K -1 (equivalently R-8, R-9 and R-10/in.) using various filler gases are obtained. It is found that polymer foams can be improved to achieve ultralow thermal conductivity by reducing CO 2 concentration, reducing radiation, increasing porosity, and using anisotropic pore geometry. In conclusion, these findings contribute to the development of highly efficient thermal insulation materials, enhancing building energy efficiency and promoting sustainable construction practices.

42 ENGINEERING↗
DOE OSTI2021

Upscaling Reactive Transport and Clogging in Shale Microcracks by Deep Learning

Fracture networks in shales exhibit multiscale features. A rock system may contain a few main fractures and thousands of microcracks, whose length and aperture are orders of magnitude smaller than the former. It is computationally prohibitive to resolve all the fractures explicitly for such multiscale fracture networks. One traditional approach is to model the small-scale features (e.g., microcracks in shales) as an effective medium. Although this fracture-matrix conceptualization significantly reduces the problem complexity, there are classes of physical processes that cannot be accurately upscaled by effective medium approximations, e.g., microcrack clogging during mineral reactions. Here, we employ deep learning in place of effective medium theory to upscale physical processes in small-scale features. Specifically, we consider reactive transport in a fracture-microcrack network where microcracks can be clogged by precipitation. A deep learning multiscale algorithm is developed, in which the microcracks are upscaled as a wall boundary condition of the main fractures. The wall boundary condition is constructed by recurrent neural networks, which take concentration histories as input and predict the solute transport from main fractures to microcracks. The deep learning multiscale algorithm is firstly employed in specific scenarios, then a general model is developed which can work under various conditions. The new approach is validated against fully resolved simulations and an analytical solution, providing a reliable and efficient solution for problems that cannot be upscaled by effective medium models.

58 GEOSCIENCES↗
DOE OSTI2023

Implications of Electron Transport Layer and Back Metal Contact Variations in Tin–Lead Perovskite Solar Cells Assessed by Spectroscopic Ellipsometry and External Quantum Efficiency

The structural and optical properties of hybrid organic–inorganic metal halide perovskite solar cells are measured by spectroscopic ellipsometry to reveal an optically distinct interfacial layer among the back contact metal, charge transport, and absorber layers. Understanding how this interfacial layer impacts performance is essential for developing higher performing solar cells. Here, this interfacial layer is modeled by Bruggeman effective medium approximations (EMAs) to contain perovskite, C 60 , BCP, and metal. External quantum efficiency (EQE) simulations that consider scattering, electronic losses, and the formation of nonparallel interfaces are created with input derived from ellipsometry structural-optical models and compared with experimental EQE to estimate optical losses. This nonplanar interface causes optical losses in short circuit current density (JSC) of up to 1.2 mA cm –2 . A study of glass/C 60 /SnO 2 /Ag or Cu and glass/C 60 /BCP/Ag film stacks shows that C 60 and BCP mix, but replacing BCP with SnO 2 can prevent mixing between the ETLs to prevent contact between C 60 and back contact metal and enable the formation of a planar interface between ETLs and back contact metals.

14 SOLAR ENERGY↗
DOE OSTI2023

COWALKER:EFFECTIVE TRANSPORT PROPERTIES OF COMPOSITE MATERIALS

SF-23-026 This software computes effective transport properties of composite materials involving fibers and nanoparticles using a random-walk algorithm that efficiently scales to an arbitrary number of processes and cores. Effective transport properties (thermal, electrical) are key to bridge the microstructure of complex materials with its macroscopic behavior. Traditional approaches either use effective medium approximations (closed mathematical expressions that are approximation for certain conditions) or continuum simulation models such as finite element or finite volume, which require the generation of a mesh for each configuration explored. cowalker leverages the equivalence between laplacian or heat equation-based models and random walks to compute the asymptotic transport properties from an ensemble of first sojourn times of a random walker moving through the composite material. This allows us to directly define a composite material as a collection of particles and use algorithms developed for molecular dynamics to quickly compute the intersection of the walker with the different interfaces in the material. cowalker is developed in C++, and it relies on the GNU Scientific Library for random generation. cowalker is currently delivered as source code, so the GSL library is not included in cowalker's distribution. A more userfriendly version, cowalker.jl is currently in development and will be released as part of cowalker.

YANGUAS-GIL, ANGEL↗
DOE OSTI2021

Propagation of radio frequency waves through turbulent plasmas

The practical and economic viability of tokamak fusion reactors depends, in a significant way, on the efficiency of radio frequency (RF) waves to deliver energy and momentum to the plasma in the core of the reactor. The RF electromagnetic waves, excited by antenna structures placed near the wall of a tokamak, have to propagate through the turbulent edge plasma along their path to the core of the fusion device. In present day experiments, the radial width of the edge region and scrape-off layer is of the order of a few centimeters. In ITER, and in future fusion reactors, this width will be of the order of tens of centimeters. Any effects on RF waves due to plasma turbulence have to be properly understood in order to optimize the delivery of RF energy and momentum into the core. This paper is on a multi-pronged, theoretical and computational, approach that is being pursued to quantify the effect of edge plasma turbulence on the propagation of RF waves. The theoretical and analytical models are based on solutions of the Faraday-Ampere equation in a magnetized plasma and on the Kirchhoff tangent plane approximation. An effective medium approach has been developed so as to approximate the permittivity of a turbulent plasma by analytical expressions. The computations are being carried out with a newly developed code ScaRF that is based on the finite difference finite domain technique for solving Maxwell's equations. The plasma permittivity can be assigned as desired. The code is being used to validate the analytical and theoretical models and to evaluate their limitations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Deciphering the influence of high-flux helium plasma exposure on tungsten nanostructure growth using real-time, in situ spectroscopic ellipsometry

In situ analysis of surfaces during high-flux plasma exposure represents a long-standing challenge in the study of plasma–material interactions. While post-mortem microscopy can provide a detailed picture of structural and compositional changes, in situ techniques can capture the dynamic evolution of the surface. In this study, we demonstrate how spectroscopic ellipsometry can be applied to the real-time characterization of W nanostructure (also known as “fuzz”) growth during exposure to low temperature, high-flux He plasmas. Strikingly, over a wide range of sample temperatures and helium fluences, the measured ellipsometric parameters (Ψ,Δ) collapse onto a single curve that can be directly correlated with surface morphologies characterized by ex situ helium ion microscopy. The initial variation in the (Ψ,Δ) parameters appears to be governed by small changes in surface roughness (<50 nm) produced by helium bubble nucleation and growth, followed by the emergence of 50 nm diameter W tendrils. This basic behavior appears to be reproducible over a wide parameter space, indicating that the spectroscopic ellipsometry may be of general practical use as a diagnostic to study surface morphologies produced by high-flux He implantation in refractory metals. An advantage of the methods outlined here is that they are applicable at low incident ion energies, even below the sputtering threshold. As an example of this application, we apply in situ ellipsometry to examine how W fuzz growth is affected both by varying ion energy and the temperature of the surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Employing ternary fission of Pu 242 as a probe of very neutron-rich matter

Detailed assessments of the ability of recent theoretical approaches to modeling existing experimental data for ternary fission confirm earlier indications that the dominant mode of cluster formation in ternary fission is clusterization in very neutron rich, very low density, essentially chemically equilibrated, nucleonic matter. An extended study and comparison of these approaches applied to ternary fission yields in the thermal neutron induced reaction 241 Pu(n th ,f) has been undertaken to refine the characterization of the source matter. Here, the resonance gas approximation has been improved taking in-medium effects on the binding energies into account. A temperature of 1.29 MeV, density of 6.7 × 10 -5 nucleons/fm 3 and proton fraction Y p = 0.035 are found to provide a good representation of yields of the ternary emitted light particles and clusters. In particular, results for Z = 1 and 2 isotopes are presented. Isotopes with larger Z are discussed, and the roles of medium and continuum effects, even at very low density are illustrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

First-principles residual resistivity using a locally self-consistent multiple scattering method

The locally self-consistent multiple scattering (LSMS) method can perform efficient first-principles calculations of systems with a large number of atoms. Here, in this paper, we combine the Kubo-Greenwood equation with LSMS, enabling us to calculate the first-principles residual resistivity of large systems. This has been implemented in the open-source code LSMS . We apply this method to selected pure elements and binary random alloys. The results compare well with experiment, and with values obtained from a first-principles effective medium technique (the Korringa-Kohn-Rostoker coherent potential approximation). We discuss future applications of this method to complex systems where other methods are not applicable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Large-scale simulation-based parametric analysis of an optimal precooling strategy for demand flexibility in a commercial office building

Achieving success with grid-interactive efficient buildings (GEBs) is closely tied to the utilization of flexible loads. A valuable strategy involves the implementation of precooling techniques before high-demand events, such as peak hours, by adjusting zone air temperature setpoints. This leads to a reduction in thermal loads and peak electricity demand during these times, as the building’s thermal mass stores and subsequently releases thermal energy. However, the effectiveness of the pre-cooling optimization is highly contingent on specific conditions such as building thermal properties, weather conditions, utility rate structure, HVAC equipment sizing, etc. Therefore, investigating the impacts of these condition-specific factors is crucial, especially when considering precooling strategies that utilize thermal mass in commercial buildings. In this paper, we first devised a novel heuristic control approach that incorporates parameterized optimal precooling thermostat schedules to enhance demand flexibility in a commercial office building. Subsequently, we conducted a thorough performance evaluation of this control strategy. Here, the optimal thermostat schedule was parameterized using three optimization variables: the precooling start time, the precooling end time, and the precooling temperature setpoint. Utilizing the DOE medium-sized office building as the virtual testbed, we showed that the parameterized schedule effectively approximates model predictive control and requires drastically reduced computational overhead. In addition, we investigated the impact of different influencing factors on the optimal precooling strategy. These factors include building thermal mass, outdoor air conditions, and energy price profiles. Using high-performance computing, we simulated a total of 225 scenarios, consisting of three levels of thermal mass, five typical outdoor air temperature profiles, and fifteen time-of-use price plans. The results demonstrate that optimal thermostat scheduling could save substantial energy cost in medium-sized office buildings with heavy thermal mass but with some energy penalty. Although the potential for cost savings is lower in buildings with low and medium thermal mass, the energy penalty remains consistent in all three thermal mass scenarios. The study also highlights the need to account for zone diversity and recognize that a one-size-fits-all-zone setpoint schedule may not be suitable for all zones and can lead to unnecessary energy wastage. Furthermore, the results highlight that while outdoor air conditions play a role in cost and energy performance, the cooling load exerts a more immediate and substantial influence on cost savings in precooling strategies. Although cost savings are comparable under certain conditions with the same cooling load, observed deviations in energy penalty indicate potential disparities in the efficiency of the HVAC system during the load-shifting process. In addition, the duration of peak pricing and the ratio between peak and off-peak times exhibit clear correlations with cost savings and energy consumption, aligning with intuitive expectations. These findings offer valuable insights for optimizing precooling strategies in office buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

QED medium effects in (anti)neutrino-nucleus and electron-nucleus scattering: Elastic scattering on nucleons

Interpretation of current and future neutrino oscillation and electron scattering experiments requires knowledge of lepton-nucleon and lepton-nucleus interactions at the percent level. We study the exchange of photons between charged particles and the nuclear medium for (anti)neutrino-, electron-, and muon-induced reactions inside a large nucleus. While quantum electrodynamics (QED)-medium contributions are formally suppressed by two powers of the electromagnetic coupling constant α when compared to the leading-order cross sections, low-energy modes and the nuclear size enhance the effect by orders of magnitude. They require a proper infrared regularization, which we implement as a screening of the electromagnetic interactions at atomic length scales or above. We provide approximate analytic expressions for the distortion of (anti)neutrino-nucleus and charged lepton-nucleus cross sections and evaluate the QED-medium effects for realistic values of the screening scale on the example of elastic scattering with nucleons inside the nucleus. We find new permille- to percent-level effects, which were not considered in either (anti)neutrino-nucleus or electron-nucleus scattering.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Projections of two-particle correlations onto transverse rapidity in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV at STAR

Two-particle correlation measurements projected onto two-dimensional, transverse rapidity coordinates (y T1 ,y T2 ), allow access to dynamical properties of the QCD medium produced in relativistic heavy-ion collisions that angular correlation measurements are not sensitive to. Here, we report non-identified charged-particle correlations for Au + Au minimum-bias collisions at $\sqrt{s_{NN}}$ = 200 GeV taken by the STAR experiment at the Relativistic Heavy-Ion Collider (RHIC). Correlations are presented as 2D functions of transverse rapidity for like-sign, unlike-sign and all charged-particle pairs, as well as for particle pairs whose relative azimuthal angles lie on the near-side, the away-side, or at all relative azimuth. The correlations are constructed using charged particles with transverse momentum p T ≥ 0.15 GeV/c, pseudorapidity from –1 to 1, and azimuthal angles from –π to π. The significant correlation structures that are observed evolve smoothly with collision centrality. The major correlation features include a saddle shape plus a broad peak with maximum near y T ≈ 3, corresponding to p T ≈ 1.5 GeV/c. The broad peak is observed in both like- and unlike-sign charge combinations and in near- and away-side relative azimuthal angles. The all-charge, all-azimuth correlation measurements are compared with the theoretical predictions of HIJING and EPOS. The results indicate that the correlations for peripheral to mid-central collisions can be approximately described as a superposition of nucleon + nucleon collisions with minimal effects from the QCD medium. Strong medium effects are indicated in mid- to most-central collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

The Problem of overlapping formation times: In-medium virtual corrections for QCD

High energy particles traversing through medium primarily lose energy by showering through hard bremsstrahlung and pair production. These splitting processes are coherent over large distances in the very high energy limit, leading to suppression from the Landau-Pomeranchuk-Migdal (LPM) effect. Avoiding soft-emission approximations, we study the cases where the coherence lengths of two consecutive splittings overlap (which is important for calculating corrections to LPM effect in QCD). Here, we will show how to compute in-medium virtual corrections to the leading order LPM emission rates for QCD. These loop corrections will be necessary for calculating properties of in-medium high energy parton showers. To simplify the calculations, we will focus on the all-gluon case and work in the large-Nc limit of QCD, where Nc is the number of colors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The LPM effect in sequential bremsstrahlung: incorporation of “instantaneous” interactions for QCD

The splitting processes of bremsstrahlung and pair production in a medium are coherent over large distances in the very high energy limit, which leads to a suppression known as the Landau-Pomeranchuk-Migdal (LPM) effect. We continue study of the case when the coherence lengths (formation lengths) of two consecutive splitting processes overlap, avoiding soft-emission approximations. Previous work made a “nearly-complete” calculation of the effect of overlapping formation times on gluonic splittings such as g → gg → ggg (with simplifying assumptions such as an infinite QCD medium and the large-N c limit). In this paper, we extend those previous rate calculations from nearly-complete to complete by including processes involving the exchange of longitudinally-polarized gluons. In the context of Lightcone Pertubation Theory, used earlier for the “nearly-complete” calculation, such exchanges are instantaneous in lightcone time and have their own diagrammatic representation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Covariant classical charged particle motion in a plane electromagnetic wave revisited

We revisit the covariant formulation of a charged particle’s motion under the influence of a plane electromagnetic wave in classical picture. In the first case, we reproduce Meyer’s result by solving relativistic equation of motion ignoring radiation reaction in vacuum assuming zero mass approximation or null vector case ($k^µk_µ = k^2 = 0$). Then, we extend this approach to include the medium effects as perturbations in the propagation of electromagnetic wave (i. e. $k^2 ≠ 0$) and explore the effects of these perturbations in mass shift of the particle.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Incorporating the effects of objects in an approximate model of light transport in scattering media

A computationally efficient radiative transport model is presented that predicts a camera measurement and accounts for the light reflected and blocked by an object in a scattering medium. The model is in good agreement with experimental data acquired at the Sandia National Laboratory Fog Chamber Facility (SNLFC). The model is applicable in computational imaging to detect, localize, and image objects hidden in scattering media. Here, a statistical approach was implemented to study object detection limits in fog.

42 ENGINEERING↗