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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 649 records · Page 36

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Mechanistic Insights into Adsorptive and Catalytic Reactions from Controllable Distributions of Metal Cations (Pd, Pt, Ni, Cr, Cu) as [M‐OH] +1 /1Al or M +2 /2Al in Zeolites

Anchoring divalent metal ions in the same zeolite framework with similar Si/Al ratio selectively as zeolite-bound M +2 or [M +2 -OH] +1 cationic species enables critical comparison of the species’ intrinsic reactivity for industrially and fundamentally relevant reactions. H-BEA zeolites with similar Si/Al ratios but differing framework Al siting were used to anchored multiple divalent metal cations (Ni, Pd, Pt, Cr, Cu) in the zeolite micropores. State-of-the-art infrared (IR) spectroscopy, electron paramagnetic resonance (EPR) measurements, including two-dimensional pulsed HYSCORE EPR, extended X-ray absorption fine structure (EXAFS), and density functional theory (DFT) calculations together provide unambiguous evidence for the selective formation of divalent metal cations as M +2 /2Al species (for H-BEA prepared in the conventional hydroxide media), and [M +2 OH] +1 /1Al species for H-BEA prepared in HF. Solid-state proton-decoupled triple-quantum magic-angle spinning (3Q MAS) NMR measurements confirmed contrasting Al distributions in the two H-BEA zeolites, which led to a contrasting divalent cation speciation. The reactivities of the two cationic species were explored for catalytic and adsorptive applications in both organometallic homogeneous and heterogeneous catalysis. This work demonstrates their divergent reactivity in ethylene dimerization, ethylene oxidation (Wacker process), selective catalytic reduction (SCR) of NO, NO adsorption, and methane oxidation. Both M +2 /2Al and [M +2 OH] +1 /1Al cations are both active for ethylene dimerization, but [M +2 OH] +1 /1Al species show higher reaction rates for each Pd, Ni, Pt. [M +2 OH] +1 /1Al is active for acetaldehyde formation in Wacker ethylene oxidation. A new active site for ethylene oligomerization is proposed that possesses a terminal OH group (Cr-OH) in Phillips catalysts evident by a nearly inactive isolated Cr +2 /2Al species that contrast an active Cr─OH motif.

Divalent metal cations in a zeolite↗

DOE OSTI2024

Intertwined charge and spin density waves in a topological kagome material

Using neutrons and x rays we show the topological kagome antiferromagnet Mn 3 Sn for T < 285 K forms a homogeneous spin and charge ordered state comprising a longitudinally polarized spin density wave with wave vector k β = k β c ̂ , a helical modulated version of the room temperature antichiral magnetic order with k χ = k χ c ̂ , and charge density waves with wave vectors 2 k β , 2 k χ , and k β + k χ . Though k χ and k β coincide for 200 K < T < 230 K , they exhibit distinct continuous T dependencies before locking to commensurate values of k β = 1 12 c * and k χ = 5 48 c * at low T . Density functional theory indicates this complex modulated state may be associated with the nesting of Fermi surfaces from correlated flat kagome bands, which host Weyl nodes that are annihilated as it forms. Published by the American Physical Society 2024

Chen, Y. (ORCID:0000000177169179)↗

DOE OSTI2025

Non-equilibrium simulations of hydraulic permeation: Role of mechanical boundary conditions in dense membranes

Understanding the mechanisms of water transport in reverse osmosis membranes is critical for improving membrane performance and guiding material design. While classical models describe transport as either solution-diffusion (SD)—involving concentration-driven diffusion through a homogeneous medium—or pore-flow (PF)—involving pressure-driven convection through percolated water channels—their applicability to crosslinked polyamide membranes remains debated. Here, using non-equilibrium molecular dynamics simulations, we investigate the impact of mechanical support conditions on pressure-driven water transport in polyamide membranes across varying crosslink densities and pressure differentials (1000–5000 bar). Two support conditions are considered: graphene-restrained, representing experimentally relevant supported membranes, and freeze-restrained, mimicking a self-supported structure. In graphene-restrained systems, water concentration gradients and constant pressure profiles emerge, consistent with SD theory and incompatible with PF assumptions due to the absence of percolated pores and sub-nanometer voids. In contrast, freeze-restrained systems display uniform water concentration and linearly decreasing pressure at 1000 bar, and exhibit compressibility-induced water gradients and partial percolation at 5000 bar, resembling PF-like behavior. However, the underlying assumptions of PF theory—continuous solvent pathways and pressure transmission through water-filled pores—are not met under most conditions. Our results demonstrate that accurate modeling of reverse osmosis membranes must incorporate realistic mechanical boundary conditions to distinguish between transport mechanisms. For dense polyamide membranes supported by porous substrates, graphene-restrained simulations best reflect experimental setups and support the SD model as the dominant mechanism of water permeation.

molecular dynamics↗

DOE OSTI2026

Suppressing Phase Segregation and Improving Stability in Mixed-Halide Perovskites through Spinel Oxide-Directed Epitaxy

Mixed-cation mixed-halide perovskite compositions are essential for achieving the required bandgaps for high-efficiency multijunction photovoltaics, yet their stability remains limited by interfacial defects, phase segregation, and degradation. Here, we introduce spinel oxides as a new family of lattice-matched substrates that enable crystalline, phase-pure, compositionally-uniform, bromide-rich perovskite film growth. The effect of spinel oxides is two-fold: reducing defects at the bottom interface by templating film growth and inducing beneficial compressive strain through mismatch-dependent substrate-perovskite lattice coupling. Spinel oxide substrates facilitate growth of highly crystalline films and eliminate detrimental secondary phases across thicknesses. Using grazing incidence X-ray diffraction, X-ray fluorescence, cathodoluminescence–scanning electron microscopy, cryogenic photoluminescence, and density functional theory, we reveal that Mg-halide bonds at the bottom interface induce lattice mismatch-dependent compressive strain that suppresses halide segregation and further reduces defect formation. In addition, films grown on spinel oxides maintain over 87% of the perovskite phase after 12 h under 100% relative humidity, as monitored by in situ grazing incidence wide-angle X-ray scattering (GIWAXS), compared to less than 70% for control samples. This work extends lattice matching from vapor-deposited epitaxial semiconductors to solution-processed halide perovskites to establish a broadly applicable strategy for defect suppression, phase homogenization, and long-term stability. Based on the fundamental science explored here, we set the stage for the development of lattice-matched spinel oxide charge transport layers to be integrated into perovskite solar cells and other optoelectronic devices.

36 MATERIALS SCIENCE↗

NASA NTRS1978

Nonneutral and quasi-neutral diffusion of weakly ionized multiconstituent plasma

A theory is developed for quasi-neutral and nonneutral ambipolar diffusion in a weakly ionized multiconstituent plasma. Magnetic field effects are disregarded. The effects of negative ions and various positive and negative ion diffusion coefficients are identified. It is found that the presence of negative ions leads to enhanced electron diffusion and a retarded diffusion of positive ions. In some cases electrons may be forcibly propelled from the area of perturbations in ionic concentrations. The ambipolar electric field is observed to cause negative ions to flow in a countergradient fashion where the gradient is small. Electrons tend toward a more homogeneous spatial distribution in the nonneutral case than in the quasi-neutral case.

Hill, R. J.↗

DOE OSTI2025

Onset of cavitation and vapor bubble development over hydrophilic and hydrophobic surfaces

Cavitation, the formation of vapor bubbles as the liquid pressure is reduced below the saturated vapor pressure, often requires a substantial negative relative pressure in a pure liquid. Classical nucleation theory (CNT) provides an estimate for the rate of cavitation but there is often a disconnect between the predictions at the molecular scale compared to observations at the macroscale. We report on mesoscale simulations of cavitation based on many-body dissipative particle dynamics (mDPD), a coarse-grained molecular dynamics (MD), which bridges the two scales. A liquid layer is confined between smooth planar walls at a constant temperature, while the pressure is reduced slowly by expanding the wall-bounded domain. The wetting properties of the liquid are determined by the parameters of the interaction potentials. With hydrophilic walls, homogeneous nucleation is observed in the liquid bulk. As a bubble forms and grows, it creates a strong pressure pulse and oscillations that cause other bubbles that may have formed slightly later to collapse. For a nearly neutral wall with a contact angle close to 90 ° , heterogeneous nucleation occurs at the walls at a smaller negative pressure and generates weaker pressure oscillations. With hydrophobic walls or seed particles, heterogeneous nucleation readily occurs, where fluctuations and the merger of transient surface bubbles are significant.

Science & Technology - Other Topics↗

NASA NTRS1986

Kinetic efficiencies of stellar wind bubbles

The theoretical fraction of a stellar wind's energy converted into the kinetic energy of an expanding swept-up shell is 20 percent in the classical theory of stellar wind bubbles. Observational estimates of this conversion factor based on the amount of ionized material in wind-swept shells about Wolf-Rayet stars generally yield results of 1 percent. If there is a substantial amount of neutral material in the shell, it will not be counted and the kinetic efficiency will be underestimated. Presented here is a dynamical estimate which accounts for this neutral material in deducing the kinetic efficiencies of stellar wind bubbles. Bubbles classified as wind-blown shells have kinetic efficiencies in line with theoretical expectations for energy-conserving evolution in a homogeneous medium. Ringlike nebulae have significantly lower efficiencies, probably because they have been 'poisoned' by the photoevaporation of clouds engulfed during evolution into a cloudy substrate.

Van Buren, D.↗

NASA NTRS2021

Micromechanics Modeling of Textiles for Re-Entry Parachute Applications

Recent flight test projects and NASA missions have highlighted the challenges associated with accurately and efficiently modeling the behavior of parachute deployment systems needed for parachute design. Moreover, parachute deployment has been identified as one of the higher risk components for such missions. The analysis of textile fabrics used for atmospheric entry is inherently complex due to the multiple scales present in the fabric structure, including individual fiber filaments at the microscale, yarn bundles of fibers at the mesoscale, and the overall woven fabric at the macroscale. Computational tools for simulating fabric behavior must be able to account for the different mechanisms present at each scale without sacrificing computational efficiency. This work examines the generalized multiscale method of cells micromechanics theory, which has previously been used for the analysis of reinforced composite structures, to unreinforced textile fabrics. Modifications to the existing composite multiscale framework, implemented in NASA’s Multiscale Analysis Tool (NASMAT), include the specific mechanics unique to unreinforced textile fabrics, and overcoming the assumptions of a fixed fiber angle. It looks to assess the feasibility of using the NASMAT tool for efficient prediction of the response of unreinforced fabrics to loading such that it can ultimately be applied to fluid structure interaction tools for the prediction of parachute deployment systems. In this work, fabric behavior is simulated in NASMAT through homogenization of a triply periodic repeating unit cell, where the geometry of the subcells can change as a function of loading to represent the relative rotation and uncrimping that can occur in fabric tows. Predictions from the amended NASMAT code are compared to experimental data for uniaxial and off-axis tension to verify the ability of the code to incorporate lower-scale mechanics in prediction of unreinforced fabrics under loading.

Micromechanics↗

DOE OSTI2022

A Micro-Mechanical Constitutive Model to Predict Hygrothermal Aging of Cross-Linked Polymers

Abstract A multi-physics material model is presented to describe the effects of temperature, oxygen, and humidity on the constitutive response of cross-linked polymers. The effect of hygrothermal damage on the mechanical integrity of the polymer matrix can be considered as the result of damage accumulation of two independent aging mechanism namely, i) thermo-oxidative, and ii) hydrolytic aging. In order to capture the mutual effects of thermo-oxidative and hydrolytic aging, an assumption has been made that each of the aging phenomenon can be superposed to each other. In fact, each of them works independently and as a result, they can compete with each other. Utilizing the theory of network decomposition, all phenomena and their correlation were modeled and thus, the strain energy function of the polymer matrix is written with respect to four independent mechanisms, i) the shrinking original matrix that has neither been attacked by water nor oxygen, ii) conversion of the first network to two new network due to the reduction and formation of cross-links, and iii) energy loss from network degradation due to attack of the water molecules to polymer active agents. Moreover, the proposed model is micro-mechanically based and is mainly relevant on thin samples due to our underlying assumption of homogeneous diffusion of oxygen and water throughout the matrix. The model has been validated against extensive data-sets obtained from experiments we specifically designed for concept validation.

Bahrololoumi, Amir↗

DOE OSTI2020

An LES-based airborne Doppler lidar simulator and its application to wind profiling in inhomogeneous flow conditions

Wind profiling by Doppler lidar is common practice and highly useful in a wide range of applications. Airborne Doppler lidar can provide additional insights relative to ground-based systems by allowing for spatially distributed and targeted measurements. Providing a link between theory and measurement, a first large eddy simulation (LES)-based airborne Doppler lidar simulator (ADLS) has been developed. Simulated measurements are conducted based on LES wind fields, considering the coordinate and geometric transformations applicable to real-world measurements. The ADLS provides added value as the input truth used to create the measurements is known exactly, which is nearly impossible in real-world situations. Thus, valuable insight can be gained into measurement system characteristics as well as retrieval strategies. As an example application, airborne Doppler lidar wind profiling is investigated using the ADLS. For commonly used airborne velocity azimuth display (AVAD) techniques, flow homogeneity is assumed throughout the retrieval volume, a condition which is violated in turbulent boundary layer flow. Assuming an ideal measurement system, the ADLS allows to isolate and evaluate the error in wind profiling which occurs due to the violation of the flow homogeneity assumption. Overall, the ADLS demonstrates that wind profiling is possible in turbulent wind field conditions with reasonable errors (root mean squared error of 0.36 m s –1 for wind speed when using a commonly used system setup and retrieval strategy for the conditions investigated). Nevertheless, flow inhomogeneity, e.g., due to boundary layer turbulence, can cause an important contribution to wind profiling error and is non-negligible. Results suggest that airborne Doppler lidar wind profiling at low wind speeds (<5 m s –1 ) can be biased, if conducted in regions of inhomogeneous flow conditions.

17 WIND ENERGY↗

NASA NTRS1987

A method of calculating the total flow from a given sea surface topography

Using a simple dynamical model of a wind-driven ocean circulation of the Stommel type, and an analytical basis developed to objectively analyze the sea surface height residuals from an altimeter and, in the process, to determine the total flow instead of just the near surface geostrophic component associated with the given sea surface topography. The method is based on first deriving the solution to the forced problem for a given wind stress required to develop a hypothetical true or perfect data field and to establishing the basis for the objective analysis. The stream function and the surface height field for the forced problem are developed in terms of certain characteristic functions with the same expansion coefficients for both fields. These characteristic functions are simply the solutions for a homogeneous elliptic equation for the stream function and the solutions of an inhomogeneous balance equation for the height field. For the objective analysis, using a sample of randomly selected height values from the true data field, the height field characteristic functions are used to fit the given topography in a least squares sense. The resulting expansion coefficients then permit the synthesis of the total flow field via the stream function characteristic modes and the solution is perfectly well behaved even along the equator. The method of solution is easily adaptable to realistic ocean basis by straight forward numerical methods. The analytical basis of the theory and the results for an ideal rectangular basin on a beta plane are described.

Rao, Desiraju B.↗

NASA NTRS1987

Effect of wave localization on plasma instabilities

The Anderson model of wave localization in random media is involved to study the effect of solar wind density turbulence on plasma processes associated with the solar type III radio burst. ISEE-3 satellite data indicate that a possible model for the type III process is the parametric decay of Langmuir waves excited by solar flare electron streams into daughter electromagnetic and ion acoustic waves. The threshold for this instability, however, is much higher than observed Langmuir wave levels because of rapid wave convection of the transverse electromagnetic daughter wave in the case where the solar wind is assumed homogeneous. Langmuir and transverse waves near critical density satisfy the Ioffe-Reigel criteria for wave localization in the solar wind with observed density fluctuations -1 percent. Numerical simulations of wave propagation in random media confirm the localization length predictions of Escande and Souillard for stationary density fluctations. For mobile density fluctuations localized wave packets spread at the propagation velocity of the density fluctuations rather than the group velocity of the waves. Computer simulations using a linearized hybrid code show that an electron beam will excite localized Langmuir waves in a plasma with density turbulence. An action principle approach is used to develop a theory of non-linear wave processes when waves are localized. A theory of resonant particles diffusion by localized waves is developed to explain the saturation of the beam-plasma instability. It is argued that localization of electromagnetic waves will allow the instability threshold to be exceeded for the parametric decay discussed above.

Levedahl, William Kirk↗

NASA NTRS1980

Effects of axisymmetric contractions on turbulence of various scales

Digitally acquired and processed results from an experimental investigation of grid generated turbulence of various scales through and downstream of nine matched cubic contour contractions ranging in area ratio from 2 to 36, and in length to inlet diameter ratio from 0.25 to 1.50 are reported. An additional contraction with a fifth order contour was also utilized for studying the shape effect. Thirteen homogeneous and nearly isotropic test flow conditions with a range of turbulence intensities, length scales and Reynolds numbers were generated and used to examine the sensitivity of the contractions to upstream turbulence. The extent to which the turbulence is altered by the contraction depends on the incoming turbulence scales, the total strain experienced by the fluid, as well as the contraction ratio and the strain rate. Varying the turbulence integral scale influences the transverse turbulence components more than the streamwise component. In general, the larger the turbulence scale, the lesser the reduction in the turbulence intensity of the transverse components. Best agreement with rapid distortion theory was obtained for large scale turbulence, where viscous decay over the contraction length was negligible, or when a first order correction for viscous decay was applied to the results.

Tan-Atichat, J.↗

NASA NTRS1986

Stability of semidiscrete approximations for hyperbolic initial-boundary-value problems: An eigenvalue analysis

A hyperbolic initial-boundary-value problem can be approximated by a system of ordinary differential equations (ODEs) by replacing the spatial derivatives by finite-difference approximations. The resulting system of ODEs is called a semidiscrete approximation. A complication is the fact that more boundary conditions are required for the spatially discrete approximation than are specified for the partial differential equation. Consequently, additional numerical boundary conditions are required and improper treatment of these additional conditions can lead to instability. For a linear initial-boundary-value problem (IBVP) with homogeneous analytical boundary conditions, the semidiscrete approximation results in a system of ODEs of the form du/dt = Au whose solution can be written as u(t) = exp(At)u(O). Lax-Richtmyer stability requires that the matrix norm of exp(At) be uniformly bounded for O less than or = t less than or = T independent of the spatial mesh size. Although the classical Lax-Richtmyer stability definition involves a conventional vector norm, there is no known algebraic test for the uniform boundedness of the matrix norm of exp(At) for hyperbolic IBVPs. An alternative but more complicated stability definition is used in the theory developed by Gustafsson, Kreiss, and Sundstrom (GKS). The two methods are compared.

Warming, Robert F.↗

NASA NTRS1990

Stratus: An interactive steady state mixed layer model for personal computers

A steady-state, horizontally homogeneous, cloud-topped marine boundary layer model based primarily on the work of Lilly (1968) and Schubert et al., (1979) is presented. The conservative thermodynamic variables are equivalent potential temperature, theta(sub e), and total water mixing ratio, q + l. Some of the differences between this and Lilly's (1968) model are: radiation is allowed to penetrate into the boundary layer; cloud top values of longwave radiation, equivalent potential temperature, and water vapor mixing ratio are linear functions of height derived from climatological data at California coastal stations; and the closure assumption assumes a weighted average of Lilly's (1968) maximum and minimum entrainment theories. This model was programmed in FORTRAN and will run interactively on an IBM-compatible personal computer. The program allows the user to specify the geographical location, the wind speed, the sea-surface temperature, the large scale horizontal divergence, and the initial guess for cloud top height. Output includes the steady state values of cloud top and cloud base height, mixed layer equivalent potential temperature and total water mixing ratio, and the associated convective and radiative fluxes.

Guinn, Thomas A.↗

DOE OSTI2021

Quantifying the hierarchy of structural and mechanical length scales in granular systems

Continuum modeling of granular media is made possible by the existence of a length scale at and above which grain-resolved properties can be meaningfully homogenized. Progress has been made in identifying such length scales relevant to local structural properties such as porosity. However, a systematic analysis of scales above which different mechanical properties can be homogenized has yet to emerge. Here, X-ray tomography and 3D X-ray diffraction data are examined to identify such length scales. The data was obtained in-situ in compressed granular materials with rigid and flexible confinement. The experimental data are supplemented with validated discrete element simulations which examine different system sizes and different boundary conditions. Overall, our study reveals a hierarchy in the length scales of granular solids, with lengths governing structural variables being the shortest, lengths of stress variables being intermediate, and lengths of energy dissipation being the longest. All structural and mechanical length scales obey a power law based on the theory of Geostatistics, implying that the length scales can be found by analyzing samples significantly smaller than the length scales themselves. The length scales are also found to be sensitive to boundary conditions, implying that they are extrinsic features of granular media.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Implementation of Perturbation Theory and Sensitivity Capabilities in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor Multiphysics analysis application, jointly developed by Argonne and Idaho National Laboratories under the DOE-NE NEAMS program. This fiscal year, capabilities for reactivity and sensitivity evaluation using perturbation methods were implemented and verified. The First Order Perturbation Method (FOPT) was employed to compute reactivity worth resulting from small perturbations in input parameters, while the Generalized Perturbation Theory (GPT) was used to evaluate sensitivities of a range of response types, including reaction rate ratio, k-eigenvalue, neutron generation time, and effective delayed neutron fraction. These perturbation methods enable users to quantify how response quantities change due to a perturbation in a input parameter without explicitly performing an additional transport simulation for each perturbed state. In particular, the GPT formulation accounts for indirect effects arising from flux changes by solving generalized inhomogeneous equations, for which a Neumann series-based iterative solution method was developed and implemented in Griffin. The implemented reactivity and sensitivity evaluation capabilities were verified using two test problems: an infinite homogeneous system and a two-dimensional hexagonal core. The results showed excellent agreement with reference solutions obtained by a direct method based on finite difference approximation as well as GPT-based results from the PERSENT code, confirming the accuracy of both reactivity and sensitivity evaluations. Additionally, preliminary uncertainty quantification (UQ) results were obtained by combining the sensitivity values computed using GPT and external covariance data, demonstrating that the implemented sensitivity results can be reliably used for uncertainty calculations. To further demonstrate the generality and practical strength of the implementation, the sensitivity evaluation capability was successfully applied to the Empire microreactor with a geometrically complex design that poses significant modeling challenges. The results confirm that Griffin enables sensitivity evaluations even for irregular and highly heterogeneous reactor configurations, thereby establishing a foundation for UQ applications in advanced reactor designs and analyses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1996

A Numerical Simulation of Scattering from One-Dimensional Inhomogeneous Dielectric Random Surfaces

In this paper, an efficient numerical solution for the scattering problem of inhomogeneous dielectric rough surfaces is presented. The inhomogeneous dielectric random surface represents a bare soil surface and is considered to be comprised of a large number of randomly positioned dielectric humps of different sizes, shapes, and dielectric constants above an impedance surface. Clods with nonuniform moisture content and rocks are modeled by inhomogeneous dielectric humps and the underlying smooth wet soil surface is modeled by an impedance surface. In this technique, an efficient numerical solution for the constituent dielectric humps over an impedance surface is obtained using Green's function derived by the exact image theory in conjunction with the method of moments. The scattered field from a sample of the rough surface is obtained by summing the scattered fields from all the individual humps of the surface coherently ignoring the effect of multiple scattering between the humps. The statistical behavior of the scattering coefficient sigma(sup 0) is obtained from the calculation of scattered fields of many different realizations of the surface. Numerical results are presented for several different roughnesses and dielectric constants of the random surfaces. The numerical technique is verified by comparing the numerical solution with the solution based on the small perturbation method and the physical optics model for homogeneous rough surfaces. This technique can be used to study the behavior of scattering coefficient and phase difference statistics of rough soil surfaces for which no analytical solution exists.

Sarabandi, Kamal↗