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 91 records · Page 5

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↗

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↗

Buoyancy Driven Shear Flows of Bubble Suspensions

In this work the gas volume fraction and the root-mean-squared fluid velocity are measured in buoyancy driven shear flows of bubble suspensions in a tall, inclined, rectangular channel. The experiments are performed under conditions where We << 1a nd Re >> 1, for which comparisons are made with kinetic theory and numerical simulations. Here Re = gamma(a(exp 2)/nu is the Reynolds number and We = rho(gamma(exp 2))a(exp 3)/sigma is the Weber number; gamma is the shear rate, a is the bubble radius, nu is the kinematic viscosity of the liquid, rho is the density of the liquid, and sigma is the surface tension of the gas/liquid interface. Kang et al. calculated the bubble phase pressure and velocity variance of sheared bubble suspensions under conditions where the bubbles are spherical and the liquid phase velocity field can be approximated using potential flow theory, i.e. We= 0 and Re >> 1. Such conditions can be achieved in an experiment using gas bubbles, with a radius of O(0.5mm), in water. The theory requires that there be no average relative motion of the gas and liquid phases, hence the motivation for an experimental program in microgravity. The necessity of performing preliminary, Earth based experiments, however, requires performing experiments where the gas phase rises in the liquid, which significantly complicates the comparison of experiments with theory. Rather than comparing experimental results with theory for a uniform, homogeneous shear flow, experiments can be compared directly with solutions of the averaged equations of motion for bubble suspensions. This requires accounting for the significant lift force acting on the gas phase when the bubbles rise parallel to the average velocity of the sheared suspension. Shear flows can be produced in which the bubble phase pressure gradient, arising from shear induced collisions amongst the bubbles, balances a body force (centrifugal or gravitational) on the gas phase. A steady, non-uniform gas volume fraction can be measured, from which the bubble phase pressure gradient can be obtained and compared to theory and numerical simulations. The presence of bounding walls further complicates the experiments, since the detailed interactions of the bubbles with bounding walls is not well understood, especially in the presence of gravity, where the momentum and energy exchange depends on the inclination of the wall.

Koch, D. L.↗

Energy Transfer in Rotating Turbulence

The influence or rotation on the spectral energy transfer of homogeneous turbulence is investigated in this paper. Given the fact that linear dynamics, e.g. the inertial waves regime tackled in an RDT (Rapid Distortion Theory) fashion, cannot Affect st homogeneous isotropic turbulent flow, the study of nonlinear dynamics is of prime importance in the case of rotating flows. Previous theoretical (including both weakly nonlinear and EDQNM theories), experimental and DNS (Direct Numerical Simulation) results are gathered here and compared in order to give a self-consistent picture of the nonlinear effects of rotation on tile turbulence. The inhibition of the energy cascade, which is linked to a reduction of the dissipation rate, is shown to be related to a damping due to rotation of the energy transfer. A model for this effect is quantified by a model equation for the derivative-skewness factor, which only involves a micro-Rossby number Ro(sup omega) = omega'/(2(OMEGA))-ratio of rms vorticity and background vorticity as the relevant rotation parameter, in accordance with DNS and EDQNM results fit addition, anisotropy is shown also to develop through nonlinear interactions modified by rotation, in an intermediate range of Rossby numbers (Ro(omega) = (omega)' and Ro(omega)w greater than 1), which is characterized by a marco-Rossby number Ro(sup L) less than 1 and Ro(omega) greater than 1 which is characterized by a macro-Rossby number based on an integral lengthscale L and the micro-Rossby number previously defined. This anisotropy is mainly an angular drain of spectral energy which tends to concentrate energy in tile wave-plane normal to the rotation axis, which is exactly both the slow and the two-dimensional manifold. In Addition, a polarization of the energy distribution in this slow 2D manifold enhances horizontal (normal to the rotation axis) velocity components, and underlies the anisotropic structure of the integral lengthscales. Finally is demonstrated the ability of a generalized EDQNM (Eddy Damped Quasi-Normal Markovian) model to predict the underlying spectral transfer structure and all the subsequent developments of classic anisotropy indicators in physical space, when compared to recent LES results. Even if the applications mainly concern developed strong turbulence, a particular emphasis is given to the strong formal analogy of this EDQNM2 model with recent weakly nonlinear approaches to wave-turbulence.

Cambon, Claude↗

Computational thermomechanics for crystalline rock. Part II: Chemo-damage-plasticity and healing in strongly anisotropic polycrystals

We present a thermal–mechanical–chemical-phase field model that captures the multi-physical coupling effects of precipitation creeping, crystal plasticity, anisotropic fracture, and crack healing in polycrystalline rock at various temperature and strain-rate regimes. This model is solved via a fast Fourier transfer solver with an operator-split algorithm to update displacement, temperature and phase field, and chemical concentration incrementally. In nuclear waste disposal in salt formation, brine inside the crystal salt may migrate along the grain boundary and cracks due to the gradient of interfacial energy and pressure. This migration has a significant implication on the permeability evolution, creep deformation, and crack healing within rock salt but is difficult to incorporate implicitly via effective medium theories compared with computational homogenization. As such, we introduce a thermodynamic framework and a corresponding computational implementation that explicitly captures the brine diffusion along the grain boundary and crack at the grain scale. Meanwhile, the anisotropic fracture and healing are captured via a high-order phase field that represents the regularized crack region in which a newly derived non-monotonic driving force is used to capture the fracture and healing due to the solution–precipitation. Numerical examples are presented to demonstrate the capacity of the thermodynamic framework to capture the multiphysics material behaviors of rock salt.

42 ENGINEERING↗

Electron–electrophile coupled dinitrogen reduction in a cerium– meta -tetraphenolate system: a computational study

The use of lanthanide complexes for catalytic dinitrogen reduction is a new development in homogeneous catalysis. Density functional theory calculations on our recently reported cerium phenolate catalyst [K 2 Ce 2 (sol) 4 (mTP) 2 ] (mTP = {(OC 6 H 2 - 2- t Bu-4-Me) 2 CH} 2 - 1,3-C 6 H 4 ; sol = OMe 2 here; THF in the experiment) have been undertaken to elucidate the reduction, activation and silylation steps at the bound dinitrogen molecule, in the presence of the reductant, potassium metal (K 0 ) and the electrophile Me 3 SiCl (TMSCl). Out of the total of six electron reductions required to cleave the N 2 , the first two-electron reduction step was found to be highly disfavoured unless potassium cations (K + ) are included, upon which the step is rendered strongly exergonic; N–Si bond formation at the two-electron stage is predicted to be unfavourable. The three-electron-reduced N 2 -adduct is found to be at the reductive activation limit in the absence of added electrophiles, which can form N-element bonds and lower the overall charge. Added electron density beyond three-electron reduction no longer localises on N 2 , preventing formal N 2 4− formation. A pathway in which both K 0 and Me 3 SiCl work in concert was modelled, and six sequential reduction–silylation steps were calculated, showing how the N–N bond is cleaved after the third reduction, eventually releasing two equivalents of N(SiMe 3 ) 3 , and regenerating the starting complex with the highest barrier of any step being 22 kcal mol −1 . We establish alkali metal coordination and coupled electron–electrophile transfer as key factors in the design of rare-earth-mediated dinitrogen functionalisation.

Ahmad, Shahbaz [Univ. of Manchester (United Kingdo↗

Light scattering properties of spheroidal particles

In the present paper, the light scattering characteristics of spheroidal particles are evaluated within the framework of a scattering theory developed for a homogeneous isotropic spheroid. This approach is shown to be well suited for computing the scattering quantities of spheroidal particles of fairly large sizes (up to a size parameter of 30). The effects of particle size, shape, index of refraction, and orientation on the scattering efficiency factors and the scattering intensity functions are studied and interpreted physically. It is shown that, in the case of oblique incidence, the scattering properties of a long slender prolate spheroid resemble those of an infinitely long circular cylinder.

Asano, S.↗

Thermoelastic response of metal matrix composites with large-diameter fibers subjected to thermal gradients

A new micromechanical theory is presented for the response of heterogeneous metal matrix composites subjected to thermal gradients. In contrast to existing micromechanical theories that utilize classical homogenization schemes in the course of calculating microscopic and macroscopic field quantities, in the present approach the actual microstructural details are explicitly coupled with the macrostructure of the composite. Examples are offered that illustrate limitations of the classical homogenization approach in predicting the response of thin-walled metal matrix composites with large-diameter fibers when subjected to thermal gradients. These examples include composites with a finite number of fibers in the thickness direction that may be uniformly or nonuniformly spaced, thus admitting so-called functionally gradient composites. The results illustrate that the classical approach of decoupling micromechanical and macromechanical analyses in the presence of a finite number of large-diameter fibers, finite dimensions of the composite, and temperature gradient may produce excessively conservative estimates for macroscopic field quantities, while both underestimating and overestimating the local fluctuations of the microscopic quantities in different regions of the composite. Also demonstrated is the usefulness of the present approach in generating favorable stress distributions in the presence of thermal gradients by appropriately tailoring the internal microstructure details of the composite.

Aboudi, Jacob↗

The Effects of Wave Escape on Fast Magnetosonic Wave Turbulence in Solar Flares

One of the leading models for electron acceleration in solar flares is stochastic acceleration by weakly turbulent fast magnetosonic waves ("fast waves"). In this model, large-scale flows triggered by magnetic reconnection excite large-wavelength fast waves, and fast-wave energy then cascades from large wavelengths to small wavelengths. Electron acceleration by large-wavelength fast-waves is weak, and so the model relies on the small-wavelength waves produced by the turbulent cascade. In order for the model to work, the energy cascade time for large-wavelength fast waves must be shorter than the time required for the waves to propagate out of the solar-flare acceleration region. To investigate the effects of wave escape, we solve the wave kinetic equation for fast waves in weak turbulence theory, supplemented with a homogeneous wave-loss term.We find that the amplitude of large-wavelength fast waves must exceed a minimum threshold in order for a significant fraction of the wave energy to cascade to small wavelengths before the waves leave the acceleration region.We evaluate this threshold as a function of the dominant wavelength of the fast waves that are initially excited by reconnection outflows.

Pongkitiwanichakul, Peera↗

Simultaneous Q-switching and mode-locking in the CW Nd:YAG laser

The theory of transient mode-locking for an active modulator in a laser with a homogeneously broadened line is presented. The theory is applied to simultaneously Q-switched and mode-locked Nd:YAG lasers and good agreement between theory and experiment is obtained. The main conclusion is that under usual Q-switched operating conditions the mode-locking process does not have sufficient time to build up to steady-state conditions. We also present a method to overcome this problem by allowing the laser to prelase before the Q-switch is opened. Mode-locked pulses whose width approaches the steady-state value are obtained. The transient mode-locking theory presented here also applies to actively mode-locking TEA CO2 lasers and to other types of simultaneously pulsed and mode-locked lasers.

Kuizenga, D. J.↗

Antennas in matter: Fundamentals, theory, and applications

The volume provides an introduction to antennas and probes embedded within or near material bodies such as the earth, the ocean, or a living organism. After a fundamental analysis of insulated and bare antennas, an advanced treatment of antennas in various media is presented, including a detailed study of the electromagnetic equations in homogeneous isotropic media, the complete theory of the bare dipole in a general medium, and a rigorous analysis of the insulated antenna as well as bare and insulated loop antennas. Finally, experimental models and measuring techniques related to antennas and probes in a general dissipative or dielectric medium are examined.

King, R. W. P.↗

Theory of beat-resonant coupling of electrostatic modes

A general expression is derived for the beat-resonant coupling electrostatic modes in a Vlasov plasma. The result for the coupling of two modes has a simple structure: the appropriate momentum gradient of the equilibrium particle distribution is weighted by a positive coupling coefficient and averaged over the resonance surface in momentum space. The contributions of all the resonance surfaces are then summed. This basic structure had been previously exhibited only for specific homogeneous plasma models. The present theory, which unifies and greatly simplifies these individual treatments, is based on a variational formulation of the Vlasov-Poisson equations. Using Lie transforms, the variational principle is reexpressed in oscillation-center variables, and then the nonlinear wave dynamics are obtained from the independent variations of the wave phase and the wave amplitude. The power of the method is then applied to a strongly magnetized, strongly inhomogeneous, non-neutral plasma model.

Crawford, John David↗

On the consistency of Reynolds stress turbulence closures with hydrodynamic stability theory

The consistency of second-order closure models with results from hydrodynamic stability theory is analyzed for the simplified case of homogeneous turbulence. In a recent study, Speziale, Gatski, and MacGiolla Mhuiris showed that second-order closures are capable of yielding results that are consistent with hydrodynamic stability theory for the case of homogeneous shear flow in a rotating frame. It is demonstrated in this paper that this success is due to the fact that the stability boundaries for rotating homogeneous shear flow are not dependent on the details of the spatial structure of the disturbances. For those instances where they are -- such as in the case of elliptical flows where the instability mechanism is more subtle -- the results are not so favorable. The origins and extent of this modeling problem are examined in detail along with a possible resolution based on rapid distortion theory (RDT) and its implications for turbulence modeling.

Speziale, Charles G.↗

Efficient Reformulation of the Thermoelastic Higher-order Theory for Fgms

Functionally graded materials (FGMs) are characterized by spatially variable microstructures which are introduced to satisfy given performance requirements. The microstructural gradation gives rise to continuously or discretely changing material properties which complicate FGM analysis. Various techniques have been developed during the past several decades for analyzing traditional composites and many of these have been adapted for the analysis of FGMs. Most of the available techniques use the so-called uncoupled approach in order to analyze graded structures. These techniques ignore the effect of microstructural gradation by employing specific spatial material property variations that are either assumed or obtained by local homogenization. The higher-order theory for functionally graded materials (HOTFGM) is a coupled approach developed by Aboudi et al. (1999) which takes the effect of microstructural gradation into consideration and does not ignore the local-global interaction of the spatially variable inclusion phase(s). Despite its demonstrated utility, however, the original formulation of the higher-order theory is computationally intensive. Herein, an efficient reformulation of the original higher-order theory for two-dimensional elastic problems is developed and validated. The use of the local-global conductivity and local-global stiffness matrix approach is made in order to reduce the number of equations involved. In this approach, surface-averaged quantities are the primary variables which replace volume-averaged quantities employed in the original formulation. The reformulation decreases the size of the global conductivity and stiffness matrices by approximately sixty percent. Various thermal, mechanical, and combined thermomechanical problems are analyzed in order to validate the accuracy of the reformulated theory through comparison with analytical and finite-element solutions. The presented results illustrate the efficiency of the reformulation and its advantages in analyzing functionally graded materials.

Bansal, Yogesh↗

The Langdon effect in laser plasmas: Absorption and conduction

A plasma heated by inverse bremsstrahlung absorption of laser light develops a non-Maxwellian electron distribution function, called the Langdon effect [A. B. Langdon, Phys. Rev. Lett. 44, 575 (1980)]. These non-Maxwellian distributions are sufficiently long-lived to impact the absorption processes itself as well as the transport of heat by electrons. The theory of the Langdon effect in a homogeneous plasma is reviewed to clarify some aspects of Langdon's derivation as well as to confirm that the widely used super-Gaussian approximation works fairly well to describe the shape of the distribution function and reduction of the absorption rate. The Langdon effect on thermal conduction in an inhomogeneous plasma is developed by considering perturbations in a homogeneous absorbing plasma, which develops a heat flux due to both temperature and density gradients. A practical theory of the heat flux is developed by fitting the results of Vlasov–Fokker–Planck simulations, which avoids several approximations that compromised the usefulness of past theoretical predictions, most critically, the effect of electron–electron collisions on the fluxes. The present fits parameterize the coefficients of the temperature gradient (thermal conductivity) and the density gradient for a plasma of any ionization state and for any laser intensity where the theory of the Langdon effect remains locally valid. It is expected that this generalized theory of heat flow in an absorbing plasma will improve the predictive capability of radiation-hydrodynamics simulations of laser-produced plasmas, especially those formed in inertial confinement fusion experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An evaluation of a coupled microstructural approach for the analysis of functionally graded composites via the finite-element method

A comparison is presented between the predictions of the finite-element analysis and a recently developed higher-order theory for functionally graded materials subjected to a thorough-thickness temperature gradient. In contrast to existing micromechanical theories that utilize classical (i.e., uncoupled) homogenization schemes to calculate micro-level and macro-level stress and displacement fields in materials with uniform or nonuniform fiber spacing (i.e., functionally graded materials), the new theory explicitly couples the microstructural details with the macrostructure of the composite. Previous thermo-elastic analysis has demonstrated that such coupling is necessary when: the temperature gradient is large with respect to the dimension of the reinforcement; the characteristic dimension of the reinforcement is large relative to the global dimensions of the composite and the number of reinforcing fibers or inclusions is small. In these circumstances, the standard micromechanical analyses based on the concept of the representative volume element used to determine average composite properties produce questionable results. The comparison between the predictions of the finite-element method and the higher-order theory presented herein establish the theory's accuracy in predicting thermal and stress fields within composites with a finite number of fibers in the thickness direction subjected to a thorough-thickness thermal gradient.

Pindera, Marek-Jerzy↗

Activation of Low-Valent, Multiply M–M Bonded Group VI Dimers toward Catalytic Olefin Metathesis via Surface Organometallic Chemistry

Olefin metathesis is a broadly employed reaction with applications that range from fine chemicals to materials and petrochemicals. The design and investigation of olefin metathesis catalysts have been ongoing for over half a century, with advancements made in terms of activity, stability, and selectivity. Immobilization of organometallic complexes onto solid supports such as silica or alumina is a promising strategy for catalyst heterogenization, often resulting in increased activity and stability. Consequently, a broad range of early transition metal catalysts bearing alkyl, oxide/alkoxide, and amide ligands have been grafted onto silica and their reactivities investigated. In this paper, we report a series of silica-supported tungsten and molybdenum dimers (X 3 M≡MX 3, where M = W and Mo; X = neopentyl, tert-butoxide, and dimethyl amide) and their reactivities toward catalytic olefin metathesis. Dynamic nuclear polarization (DNP)-enhanced solid-state nuclear magnetic resonance (SSNMR), diffuse reflectance infrared Fourier transform (DRIFT), UV resonance Raman, and X-ray absorption (XAS) spectroscopies suggest that upon heterogenization the dimers bind to the surface in a monopodal fashion, with the M≡M triple bond remaining intact. These structural assignments were further corroborated by density functional theory (DFT) calculations. While the homogeneous dimer counterparts are inert, the supported low-valent alkyl W and Mo dimers become active for the disproportionative self-metathesis of propylene to ethylene and butenes and 4-nonene to 4-octene and 5-decene under mild conditions. The lack of activity observed for the free and supported tert-butoxide and dimethyl amide dimers likely suggests that the neopentyl groups are necessary for the formation of a putative alkylidene active species. The difference in reactivity between the free and supported dimers could be explained either by the lowering of the activation barrier of the complex through the electronic effects of the surface or by site isolation of catalytically relevant reactive intermediates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Data-driven electron-diffraction approach reveals local short-range ordering in CrCoNi with ordering effects

Abstract The exceptional mechanical strength of medium/high-entropy alloys has been attributed to hardening in random solid solutions. Here, we evidence non-random chemical mixing in a CrCoNi alloy, resulting from short-range ordering. A data-mining approach of electron nanodiffraction enabled the study, which is assisted by neutron scattering, atom probe tomography, and diffraction simulation using first-principles theory models. Two samples, one homogenized and one heat-treated, are observed. In both samples, results reveal two types of short-range-order inside nanoclusters that minimize the Cr–Cr nearest neighbors (L1 2 ) or segregate Cr on alternating close-packed planes (L1 1 ). The L1 1 is predominant in the homogenized sample, while the L1 2 formation is promoted by heat-treatment, with the latter being accompanied by a dramatic change in dislocation-slip behavior. These findings uncover short-range order and the resulted chemical heterogeneities behind the mechanical strength in CrCoNi, providing general opportunities for atomistic-structure study in concentrated alloys for the design of strong and ductile materials.

36 MATERIALS SCIENCE↗