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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Optimizing group structures using hierarchical division

Creating group structures with few groups that give low errors is a difficult problem in reactor analysis. In recent years, automated optimization techniques have been applied to this task. We continue this trend by applying the hierarchical division algorithm to generate optimized group structures that minimize a cost function. At each stage, the algorithm adds a single group boundary to an existing group structure, dividing one group into two to increase the resolution of the group structure. The location of the added boundary is the one that gives the lowest error over all possible new boundary locations. Our implementation requires a beginning group structure, a set of candidate new boundary locations, and a set of reference reaction rates. As a proof of concept, we used WIMS-69 as the initial group structure, XMAS-172 as the ending group structures, and a 344-group reference group structure. Testing on two simple, homogenized reactor problems, we found that hierarchical division was able to reduce the error by a factor of around 5 with an increase of only 15% in the number of groups. Because hierarchical division can get stuck in local minima, it often reaches a plateau in its error reduction capability as many groups are added. Nevertheless, we find hierarchical division has strong potential to make good group structures into great group structures at a modest increase in computational cost. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

VARI3D & PERSENT: Perturbation and Sensitivity Analysis

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code. This manuscript serves as a single manual for two separate codes: VARI3D and PERSENT. The VARI3D code (VARIational 3D) is based upon the classic finite difference diffusion theory solver available in DIF3D. The PERSENT code (PERturbation and SENitivity for Transport) is based upon the variational nodal method employed in DIF3D termed VARIANT. The VARIANT solver was added to DIF3D in 1995 and has seen continued development and use for the last 18 years. Because VARI3D primarily uses deprecated coding practices, rather than incorporating the perturbation and sensitivity treatments for transport within VARI3D, a new coding development was built using modern Fortran coding. The primary purpose of this manual is to describe the theory behind PERSENT (and by convenience, that of VARI3D) and discuss the input and output of PERSENT along with giving potential users an idea of how to use it. While this manuscript does describe the input and output of VARI3D, the PERSENT code is intended to be the replacement capability of VARI3D as PERSENT can generate nearly identical (if not superior) diffusion theory results. In this manuscript, the relevant aspects of generalized perturbation theory and exact perturbation theory that apply to both VARI3D and PERSENT are covered. The input and output of VARI3D is displayed by excerpting several of the example problems. Similarly, the input and output of PERSENT is displayed along with tips on how best to use the code. Note that the input and output of the inhomogeneous solver wrapped around DIF3D (DIF3D_IFS) is also discussed as it is needed to carry out some of the sensitivities in PERSENT such as reaction rate ratios. This manuscript describes several perturbation and sensitivity problems, and the results computed using PERSENT. From these sections, potential users should find that PERSENT provides not only the typical tables of numbers desired in perturbation and sensitivity analysis work, but also can visually plot the result for a more thorough understanding of the space and energy distribution (Section 5). Overall, PERSENT is observed to produce accurate reactivity worths and sensitivities for the displayed set of test problems and clearly demonstrates the need to have a transport-based sensitivity capability as evident from the thousands of percent errors observed in the 21-group hexagonal fast reactor problem (covered in Section 7). The uncertainty calculation capability is described in Section 3 and demonstrated in Section 7.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Energy Barriers for H 2 Uptake and Release by Defects in Boron Nitride

Density functional theory simulations (DFT) have been used in combination with the climbing image nudged elastic band method (CI-NEB) to determine the energy barriers required to dissociate H 2 at four types of h-BN defects and evaluate the potential for reacted H species to hop to adjacent available sites. The defects were selected based on previous work in which ab initio thermodynamics was used to identify specific defects with optimal thermodynamics for reversible H 2 storage, that is, the defects would dissociate molecular H 2 while not binding H species too strongly so that hydrogen release is favorable. The four defects investigated are the N and B monovacancies, which only offer either B or N defective sites respectively, as well as larger cluster vacancy defects named 3V(1B2N) and hexagonal 6V(3B3N), which present a mixture of B and N defect sites available for H binding. While homogeneous H 2 dissociation by B and N monovacancies have energy barriers larger than 0.4 eV, we found that mixed BN terminated defects have energy barrier lower than 0.4 eV for H 2 dissociation. In particular, H 2 dissociation by the hexagonal 6V(3B3N) defect, exclusively made of Frustrated Lewis Pairs (FLPs), is almost barrierless. The energy barrier for subsequent H hopping to available adjacent B or N sites of the defect were also calculated. It was found that the hopping of H species generally involves high energy barriers ranging from 0.26 eV to 2.57 eV. In conclusion, this work completes previous thermodynamic theoretical investigations and provides a more comprehensive picture of H 2 reaction with defective h-BN for potential H 2 storage and recovery.

Sassi, Michel J. [Pacific Northwest National Labor↗
DOE OSTI2020

A radiative transfer model for solar induced fluorescence using spectral invariants theory

Solar Induced chlorophyll Fluorescence (SIF) shows promise as an approach for estimating gross primary production (GPP) remotely. However, sun-target-sensor geometry and within-canopy absorption of SIF can alter the relationship between measured SIF and GPP, because sensors can only retrieve some unknown fraction of the total emitted SIF. Radiative transfer models that allow for variation in canopy structure and sensor angles are therefore needed to properly interpret SIF measurements. Spectral invariants allow decoupling of the wavelength-independent canopy structure and the wavelength-dependent leaf and soil spectrum in the radiative transfer process. Here we develop a simple analytical Fluorescence Radiative Transfer model based on Escape and Recollision probability (FluorRTER) to investigate the impact of canopy structure and sun-target-sensor geometry on SIF emissions. SIF simulations using the FluorRTER model agreed well the one-dimensional Soil-Canopy Observation of Photochemistry and Energy balance (SCOPE) model and the three-dimensional Fluorescence model with Weighted Photon Spread (FluorWPS) model. The fractional vegetation cover (FVC) and clumping effect have a large influence the SIF emission of 3D discontinuous canopies. For a moderate solar zenith angle (30°) and a clumped canopy (FVC?=?0.6), the difference between the directional observed SIF of a 3D discontinuous canopy and a 1D homogeneous canopy was as large as 43.2% and 38.4% for Photosystem I?+?II fluorescence at 685?nm and at 740?nm, respectively. By bridging the gap between observed SIF and total emitted SIF over 3D heterogeneous vegetation canopies, the FluorRTER model can assist with the angular normalization of SIF measurements and enable the more robust interpretation of how variations in SIF from directional and hemispherical in-situ, airborne and satellite observations relate to leaf and whole-canopy physiological processes.

Zeng, Yelu↗
DOE OSTI2022

Confinement and Collective Escape of Active Particles

Active matter broadly covers the dynamics of self-propelled particles. While the onset of collective behavior in homogenous active systems is relatively well understood, the effect of inhomogeneities such as obstacles and traps lacks overall clarity. Here, we study how interacting, self-propelled particles become trapped and released from a trap. We have found that captured particles aggregate into an orbiting condensate with a crystalline structure. As more particles are added, the trapped condensates escape as a whole. Our results shed light on the effects of confinement and quenched disorder in active matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Formation of pyrophosphates across grain boundaries induces the formation of mismatched but oriented interfaces in silver phosphate polypods

The interfaces and their misfit defects determine the materials properties for a wide range of applications, such as electronic devices, photocatalysis, and mechanical engineering, etc. However, current understanding of atomic interfacial structures is limited. Here we discover a special interfacial structure, mismatched but oriented interface via two distinct facets. Transmission and scanning electron microscopy results suggest that, in Ag 3 PO 4 polypods structures, interfaces of {1 0 0} and {1 1 0}, ({1 0 0}/(1 1 0}), {1 0 0}/{1 1 1}, {1 1 0}/{1 1 1}, and {1 0 0}/{1 0 0}, etc., have a certain orientation relationship, corresponding to the energy minima and coincident site lattice of interfacial atoms as demonstrated by molecular dynamics simulations. Density functional theory demonstrates that the formation of pyrophosphate and/or phosphates rotation to bond across the interface compensate the lattice mismatch at the interfaces, as well as deformations of Ag-O bonds. Furthermore, our work opens up a new avenue for a much wider range of interfacial structures, allow for a higher diversity of structures, and shine light on tailoring crystal structures, morphologies, and the resulting properties.

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

Dark Energy Survey Year 3 results: Measurement of the baryon acoustic oscillations with three-dimensional clustering

The three-dimensional correlation function offers an effective way to summarize the correlation of the large-scale structure even for imaging galaxy surveys. We have applied the projected three-dimensional correlation function, ξ p to measure the baryonic acoustic oscillations (BAO) scale on the first-three years Dark Energy Survey data. The sample consists of about 7 million galaxies in the redshift range 0.6< z p < 1.1 over a footprint of 4108 deg 2 . Our theory modeling includes the impact of realistic true redshift distributions beyond Gaussian photo-z approximation. ξ p is obtained by projecting the three-dimensional correlation to the transverse direction. To increase the signal-to-noise of the measurements, we have considered a Gaussian stacking window function in place of the commonly used top-hat. ξ p is sensitive to D M (z eff )/r s , the ratio between the comoving angular diameter distance and the sound horizon. Using the full sample, D M (z eff )/r s is constrained to be19.00 ± 0.67 (top-hat) and 19.15 ± 0.58 (Gaussian) at z eff = 0.835. The constraint is weaker than the angular correlation w constraint 18.84 ± 0.50), and we trace this to the fact that the BAO signals are heterogeneous across redshift. While ξ p responds to the heterogeneous signals by enlarging the error bar, w can still give a tight bound on D M s in this case. When a homogeneous BAO-signal subsample in the range 0.7 < z p <1.0 (z eff = 0.845) is considered, ξp yields 19.80 ± 0.67 (top-hat) and 19.84 ± 0.53 (Gaussian). The latter is mildly stronger than the w constraint (19.86 ± 0.55). We find that the ξ p results are more sensitive to photo-z because ξ p keeps the three-dimensional clustering information causing it to be more prone to photo-z noise. The Gaussian window gives more robust results than the top-hat as the former is designed to suppress the low signal modes. ξ p and the angular statistics such as w have their own pros and cons, and they serve an important crosscheck with each other.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Assessing factors that determine adatom migration and clustering on a thin film oxide; Pt 1 and Rh 1 on the “29” Cu x O/Cu(1 1 1) surface

The dynamical nature of single-site catalysts under non-equilibrium conditions poses significant challenges in the full characterization of their active sites. The “29” oxide is a thin Cu x O film grown on Cu(111) which provides a well-defined surface on which to study the structure and chemistry of atomically dispersed precious metals. A combination of experimental and first-principles approaches is used to investigate the factors that influence the mobility of metal adatoms and their clustering tendencies. Our work elucidates how a homogenous coordination environment of neutral single-atoms can be obtained when its tendency to break oxide bonds within its vicinity is low. Herein we demonstrate that this in turn affects the chemical atomically dispersed species on thin-film oxides, as defects on the oxide would in fact allow such species to be more accessible to CO. We thus highlight the importance of understanding the relationship between defects and atomically dispersed active sites upon designing single-site catalysts.

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

Influence of the Turing instability on the motion of domain boundaries

Turing's theory of pattern formation has provided crucial insights into the behavior of various biological, geographical, and chemical systems over the last few decades. Existing studies have focused on moving-boundary Turing systems for which the motion of the boundary is prescribed by an external agent. In this paper, we present an extension of this theory to a class of systems in which the front motion is governed by the physical processes that occur within the domain. Biological systems exhibiting apically dominant growth and corrosion of metals and alloys highlight some of the noteworthy examples of such systems. In this study, we characterize the nature of interaction between the moving front and the Turing-instability for both an activator-inhibitor and an activator-substrate model. Behavioral regimes of periodic, as well as nonperiodic (nonconstant), growth rates are obtained. Furthermore, the trends in the first show striking similarities with the cyclic-boundary-kinetics observed in experimental systems. In general, a stationary, periodic structure is also left behind the moving front. If the periodicity of the boundary kinetics agrees with the allowed range of the stable-periodic solutions, the pattern formed tends to persist. Lastly, it evolves to a nearby energy-minimum either by peak-splitting, peak-decay, or by settling down to a spatially homogeneous state.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Catalytic, Spectroscopic, and Theoretical Studies of Fe 4 S 4 -Based Coordination Polymers as Heterogenous Coupled Proton–Electron Transfer Mediators for Electrocatalysis

Iron–sulfur clusters play essential roles in biological systems, and thus synthetic [Fe 4 S 4 ] clusters have been an area of active research. Recent studies have demonstrated that soluble [Fe 4 S 4 ] clusters can serve as net H atom transfer mediators, improving the activity and selectivity of a homogeneous Mn CO 2 reduction catalyst. Here, in this work, we demonstrate that incorporating these [Fe 4 S 4 ] clusters into a coordination polymer enables heterogeneous H atom transfer from an electrode surface to a Mn complex dissolved in solution. A previously reported solution-processable Fe 4 S 4 -based coordination polymer was successfully deposited on the surfaces of different electrodes. The coated electrodes serve as H atom transfer mediators to a soluble Mn CO 2 reduction catalyst displaying good product selectivity for formic acid. Furthermore, these electrodes are recyclable with a minimal decrease in activity after multiple catalytic cycles. The heterogenization of the mediator also enables the characterization of solution-phase and electrode surface species separately. Surface enhanced infrared absorption spectroscopy (SEIRAS) reveals spectroscopic signatures for an in situ generated active Mn–H species, providing a more complete mechanistic picture for this system. The active species, reaction mechanism, and the protonation sites on the [Fe 4 S 4 ] clusters were further confirmed by density functional theory calculations. The observed H atom transfer reactivity of these coordination polymer-coated electrodes motivates additional applications of this composite material in reductive H atom transfer electrocatalysis.

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

Modeling dust mineralogical composition: sensitivity to soil mineralogy atlases and their expected climate impacts

Soil dust aerosols are a key component of the climate system, as they interact with short- and long-wave radiation, alter cloud formation processes, affect atmospheric chemistry and play a role in biogeochemical cycles by providing nutrient inputs such as iron and phosphorus. The influence of dust on these processes depends on its physicochemical properties, which, far from being homogeneous, are shaped by its regionally varying mineral composition. The relative amount of minerals in dust depends on the source region and shows a large geographical variability. However, many state-of-the-art Earth system models (ESMs), upon which climate analyses and projections rely, still consider dust mineralogy to be invariant. The explicit representation of minerals in ESMs is more hindered by our limited knowledge of the global soil composition along with the resulting size-resolved airborne mineralogy than by computational constraints. In this work we introduce an explicit mineralogy representation within the state-of-the-art Multiscale Online Nonhydrostatic AtmospheRe CHemistry (MONARCH) model. We review and compare two existing soil mineralogy datasets, which remain a source of uncertainty for dust mineralogy modeling and provide an evaluation of multiannual simulations against available mineralogy observations. Soil mineralogy datasets are based on measurements performed after wet sieving, which breaks the aggregates found in the parent soil. Our model predicts the emitted particle size distribution (PSD) in terms of its constituent minerals based on brittle fragmentation theory (BFT), which reconstructs the emitted mineral aggregates destroyed by wet sieving. Our simulations broadly reproduce the most abundant mineral fractions independently of the soil composition data used. Feldspars and calcite are highly sensitive to the soil mineralogy map, mainly due to the different assumptions made in each soil dataset to extrapolate a handful of soil measurements to arid and semi-arid regions worldwide. For the least abundant or more difficult-to-determine minerals, such as iron oxides, uncertainties in soil mineralogy yield differences in annual mean aerosol mass fractions of up to ~ 100 %. Although BFT restores coarse aggregates including phyllosilicates that usually break during soil analysis, we still identify an overestimation of coarse quartz mass fractions (above 2 µm in diameter). In a dedicated experiment, we estimate the fraction of dust with undetermined composition as given by a soil map, which makes up ~ 10 % of the emitted dust mass at the global scale and can be regionally larger. Changes in the underlying soil mineralogy impact our estimates of climate-relevant variables, particularly affecting the regional variability of the single-scattering albedo at solar wavelengths or the total iron deposited over oceans. All in all, this assessment represents a baseline for future model experiments including new mineralogical maps constrained by high-quality spaceborne hyperspectral measurements, such as those arising from the NASA Earth Surface Mineral Dust Source Investigation (EMIT) mission.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A Rhenium Bis -tetramethylphenanthroline Catalyst for CO 2 Reduction to Formate

Catalytic CO 2 reduction reactions featuring high selectivity toward formate are relatively rare. In some homogeneous molecular CO 2 -reducing electrocatalysis, using triethylamine (TEA) and isopropanol (IPA) as additives improves catalytic performance in producing formate. In this paper, we investigate whether the rhenium(I) bis-diimine dicarbonyl complexes, cis-[Re(N^N) 2 (CO) 2 ] + , where N^N is 2,2’-bipyridine ([1] + ) or 3,4,7,8-tetramethyl-1,10-phenanthroline ([2] + ), are capable of electrocatalytically reducing CO 2 to formate in acetonitrile containing TEA and IPA. Catalyst [1] + was ineffective at CO 2 reduction, yielding formate quantities comparable to those produced in experiments without the catalyst. Catalyst [2] + , however, is a promising electrocatalyst for the CO 2 reduction reaction in the presence of TEA and IPA, with formate being produced in millimolar concentrations (10.5 mM), as detected by 1 H NMR spectroscopy after 6 h electrolysis (formate Faradaic efficiency = 11%, with the major balance going to H 2 ). Upon more detailed examination, [2] + exhibited a turnover frequency (TOF) of 12 s –1 for formate, comparable to other leading molecular catalysts that competently execute this reduction. Combinations of spectroscopy, electrochemistry, and theory were used to better understand the mechanism of CO 2 reduction by [2] + . Fourier transform infrared spectroelectrochemical (FTIR-SEC) data provided no evidence for CO ligand dissociation or substitution upon one- and two-electron reduction of [2] + , suggesting that a mechanism distinct from one that is metal-hydride-based is operative in catalysis. Computational studies guide mechanistic investigations toward the proposed formation of a hydrophenanthroline-based intermediate responsible for hydride transfer to CO 2 and electrocatalytic formate production from [2] + .

Beverages↗
DOE OSTI2021

Mechanistic and Electronic Insights into a Working NiAu Single-Atom Alloy Ethanol Dehydrogenation Catalyst

Elucidation of reaction mechanisms and the geometric and electronic structure of the active sites themselves is a challenging, yet essential task in the design of new heterogeneous catalysts. Such investigations are best implemented via a multi-pronged approach that comprises ambient pressure catalysis, surface science, and theory. Herein, we employ this strategy to understand the workings of NiAu single-atom alloy (SAA) catalysts for the selective non-oxidative dehydrogenation of ethanol to acetaldehyde and hydrogen. The atomic dispersion of Ni is paramount for selective ethanol to acetaldehyde conversion, and we show that even the presence of small Ni ensembles in the Au surface results in the formation of undesirable byproducts via C-C scission. Spectroscopic, kinetic, and theoretical investigations of the reaction mechanism reveal that both C-H and O-H bond cleavage steps are kinetically relevant and single Ni atoms are confirmed as the active sites. X-ray absorption spectroscopy studies allow us to follow the charge of the Ni atoms in the Au host before, under, and after a reaction cycle. Specifically, in the pristine state the Ni atoms carry a partial positive charge which increases upon coordination to the electronegative oxygen in ethanol and decreases upon desorption. Furthermore, this type of oxidation state cycling during reaction is similar to the behavior of single-site homogenous catalysts. Given the unique electronic structure of many single-site catalysts, such a combined approach in which the atomic-scale catalyst structure and charge state of the single atom dopant can be monitored as a function of its reactive environment is a key step towards developing structure function relationships that inform the design of new catalysts.

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

Kibble mechanism for electroweak magnetic monopoles and magnetic fields

The vacuum manifold of the standard electroweak model is a three-sphere when one considers homogeneous Higgs field configurations. For inhomogeneous configurations we argue that the vacuum manifold is the Hopf fibered three sphere and that this viewpoint leads to general criteria to detect electroweak monopoles and Z-strings. We extend the Kibble mechanism to study the formation of electroweak monopoles and strings during electroweak symmetry breaking. The distribution of magnetic monopoles produces magnetic fields that have a spectrum B λ ∝ λ –2 , where λ is a smearing length scale. Even as the magnetic monopoles annihilate due to the confining Z-strings, the magnetic field evolves with the turbulent plasma and may be relevant for cosmological observations.

nonperturbative effects↗
DOE OSTI2020

Relativistic Effects in Magnetic Circular Dichroism: Restricted Magnetic Balance and Temperature Dependence

Magnetic circular dichroism of transition metal complexes and open-shell systems are challenging to simulate and analyze, mainly due to the interplay of spin–orbit couplings and finite-magnetic-field induced Zeeman effects with the complex selection rules dictated by the circularly polarized light. In this study, we introduce an ab initio relativistic two-component formalism based on the restricted magnetic-balanced Hamiltonian for simulating MCD spectra. Both homogeneous finite magnetic field and relativistic effects are included variationally in the ground state reference. Finite-field London orbitals are used to enforce the constrained gauge-origin independence in the calculation using localized atomic orbitals. Through benchmark studies of AuCl 4 – , Pt(CN) 4 2– , and Mo(CN) 8 3– , we discuss how relativistic effects are manifested in MCD for both closed-shell and open-shell molecular complexes and how the interplay between spin–orbit coupling and magnetic field modulates the MCD selection rules. Finally, an investigation on temperature-dependent MCD is carried out and compared to experiment.

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

A Micromorphic Length-Scale Coupling Framework for the Determination of Higher-Order Constitutive Models and the Multi-Scale Simulation of Heterogeneous Materials [Thesis]

Heterogeneous materials and materials with complex microstructures pose a unique challenge in the development of accurate models of their response to external stimuli. These difficulties principally arise due to the difficulty in characterizing and modeling the constituents and their interactions. It is usually possible, though non-trivial, to construct an explicit representation of the microstructure (a direct numeric simulation or DNS) but the method by which the complex modes of deformation and other processes can be homogenized to a reduced-order approximation is frequently unclear. Many homogenization approaches are ad-hoc and lack a strong justification beyond ease of computation. Furthermore, the homogenization approach can, in some cases, not utilize the full breadth of information available in the computation of the macro-scale stresses and deformation measures. It is also noteworthy that homogenization, by its very nature, will tend to obfuscate details of processes occurring at the lower length-scale. It is therefore of interest to include as much information as possible in the construction of the reduced-order model so as to be predictive in a variety of loading environments. We here present a length-scale bridging technique based upon the micromorphic continuum mechanics of Eringen which incorporates volume and surface area averages as a part of its construction. This approach enforces the balance equations at the micro-scale and then studies the effect of the spatially varying nature on the macro-scale. This leads, naturally, to further balance equations which are solved at the macro-scale. This work details a homogenization framework which arises naturally from the micromorphic construction of Eringen attempting to introduce no definitions beyond which are justifiable from micro-structural considerations. One of the results of this effort is the so-called “micromorphic filter” which has been developed and applied to several DNS to demonstrate its effectively. In order to determine the macroscopic degrees of freedom we utilize the special case of an overlap coupling technique where the macro-scale is fully constrained to the micro-scale. This enables us to study the resulting material properties as expected but also allows us to further study the boundary conditions on the additional degrees of freedom at the macro-scale.

42 ENGINEERING↗
DOE OSTI2022

The role of breakup and coalescence in fine-scale bubble-induced turbulence. I. Dynamics

We study the effect of bubble breakup and coalescence on fine-scale dynamics of bubbly turbulent flows using direct numerical simulations. We perform two different simulations of dilute bubbly flows of void fraction 0.5%: one with bubbles breaking up and coalescing and the other without these physical processes. The volume of the fluid method is used for simulating bubbles undergoing breakup and coalescence while the bubbles are treated as rigid spheres in the immersed boundary method simulation. The energy spectrum in both types of simulation, consistent with previous studies, exhibits a -3 slope. We follow a single infinitesimal fluid element as it evolves to understand velocity gradient dynamics using conditional mean trajectories. We note finite-time divergence when the fluid element evolves under the action of inertial and pressure dynamics. The inertial, pressure, and viscous velocity gradient dynamics, when considered individually, produce the same results for bubble-induced turbulence (BIT) as with the classical homogeneous isotropic turbulence (HIT). Yet when the overall velocity gradient dynamics is considered, BIT results in non-cyclic trajectories moving toward stable node and unstable saddle while classical HIT shows cyclic behavior in their trajectories that move toward the origin. Interestingly, both the volume of fluid and immersed boundary simulations produce similar results. Therefore, there are two main takeaways from this research. First, new velocity gradient models are needed for BIT as their velocity gradients behavior is entirely different from the HIT. Second, we can neglect the bubble topology, breakup, and coalescence while studying or modeling the fine-scale dynamics of BIT.

42 ENGINEERING↗
DOE OSTI2022

Approximating the linear response of physical chaos

Abstract Parametric derivatives of statistics are highly desired quantities in prediction, design optimization and uncertainty quantification. In the presence of chaos, the rigorous computation of these quantities is certainly possible, but mathematically complicated and computationally expensive. Based on Ruelle’s formalism, this paper shows that the sophisticated linear response algorithm can be dramatically simplified in higher-dimensional systems featuring statistical homogeneity in the physical space. We argue that the contribution of the SRB (Sinai–Ruelle–Bowen) measure gradient, which is an integral yet the most cumbersome part of the full algorithm, is negligible if the objective function is appropriately aligned with unstable manifolds. This abstract condition could potentially be satisfied by a vast family of real-world chaotic systems, regardless of the physical meaning and mathematical form of the objective function and perturbed parameter. We demonstrate several numerical examples that support these conclusions and that present the use and performance of a simplified linear response algorithm. In the numerical experiments, we consider physical models described by differential equations, including Lorenz 96 and Kuramoto–Sivashinsky.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗