Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical homogenization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 127 records · Page 7

Exponential-fitted methods for integrating stiff systems of ordinary differential equations: Applications to homogeneous gas-phase chemical kinetics

Conventional algorithms for the numerical integration of ordinary differential equations (ODEs) are based on the use of polynomial functions as interpolants. However, the exact solutions of stiff ODEs behave like decaying exponential functions, which are poorly approximated by polynomials. An obvious choice of interpolant are the exponential functions themselves, or their low-order diagonal Pade (rational function) approximants. A number of explicit, A-stable, integration algorithms were derived from the use of a three-parameter exponential function as interpolant, and their relationship to low-order, polynomial-based and rational-function-based implicit and explicit methods were shown by examining their low-order diagonal Pade approximants. A robust implicit formula was derived by exponential fitting the trapezoidal rule. Application of these algorithms to integration of the ODEs governing homogenous, gas-phase chemical kinetics was demonstrated in a developmental code CREK1D, which compares favorably with the Gear-Hindmarsh code LSODE in spite of the use of a primitive stepsize control strategy.

Pratt, D. T.↗

Effect of the Reynolds Number on the Freestream Disturbance Environment in a Mach 6 Nozzle

To understand the impact of unit Reynolds number on the acoustic disturbance field inside a high-speed wind tunnel, we use Direct Numerical Simulations (DNS) to model the turbulent boundary layers along the walls of a quasi-two-dimensional nozzle configuration. Intended as a stepping stone to fully three-dimensional simulations of freestream noise inside the NASA 20-Inch Mach 6 Wind Tunnel, the present simulations are based on periodic boundary conditions across the spanwise width of the computational domain that corresponds to about one third of the actual tunnel width. These simulations are performed at four different unit Reynolds numbers, ranging from 3.56e6 to 14.0e6 per meter. The predominantly hydrodynamic fluctuations inside the boundary layer are shown to be nearly unaffected by the presence of freestream forcing associated with the impinging acoustic radiation from the opposite wall. Thus, the Reynolds number trends associated with boundary-layer quantities are consistent with previously published DNS of flat-plate boundary layers at similar Mach numbers and wall temperature ratios. The unsteady disturbance environment within the nozzle core region is found to be approximately spatially homogeneous and purely acoustic in nature. The numerical results are used to make comparisons with tunnel noise measurements by Chou et al. [1,2]. Unlike previous comparisons involving static-pressure fluctuations based on the DNS and pitot-pressure fluctuations measured in the wind tunnel, direct comparisons involving the fluctuations in the same physical quantity, namely, the streamwise mass flux, have been reported for the first time. The predicted decrease in the root-mean-square fluctuations in pressure and mass flux with an increasing unit Reynolds number is in agreement with the measurements in the NASA 20-Inch Mach 6 Wind Tunnel. Additional details of the acoustic radiation field are quantified and should be useful toward a digital synthesis of the tunnel disturbance environment that would enable realistic simulations of the natural transition process.

Turbulence↗

Effect of the Reynolds Number on the Freestream Disturbance Environment in a Mach 6 Nozzle

To understand the impact of unit Reynolds number on the acoustic disturbance field inside a high-speed wind tunnel, we use Direct Numerical Simulations (DNS) to model the turbulent boundary layers along the walls of a quasi-two-dimensional nozzle configuration. Intended as a stepping stone to fully three-dimensional simulations of freestream noise inside the NASA 20-Inch Mach 6 Wind Tunnel, the present simulations are based on periodic boundary conditions across the spanwise width of the computational domain that corresponds to about one third of the actual tunnel width. These simulations are performed at four different unit Reynolds numbers, ranging from 3.56e6 to 14.0e6 per meter. The predominantly hydrodynamic fluctuations inside the boundary layer are shown to be nearly unaffected by the presence of freestream forcing associated with the impinging acoustic radiation from the opposite wall. Thus, the Reynolds number trends associated with boundary-layer quantities are consistent with previously published DNS of flat-plate boundary layers at similar Mach numbers and wall temperature ratios. The unsteady disturbance environment within the nozzle core region is found to be approximately spatially homogeneous and purely acoustic in nature. The numerical results are used to make comparisons with tunnel noise measurements by Chou et al. [1,2]. Unlike previous comparisons involving static-pressure fluctuations based on the DNS and pitot-pressure fluctuations measured in the wind tunnel, direct comparisons involving the fluctuations in the same physical quantity, namely, the streamwise mass flux, have been reported for the first time. The predicted decrease in the root-mean-square fluctuations in pressure and mass flux with an increasing unit Reynolds number is in agreement with the measurements in the NASA 20-Inch Mach 6 Wind Tunnel. Additional details of the acoustic radiation field are quantified and should be useful toward a digital synthesis of the tunnel disturbance environment that would enable realistic simulations of the natural transition process.

Acoustics↗

Dipole Alignment in Rotating MHD Turbulence

We present numerical results from long-term CPU and GPU simulations of rotating, homogeneous, magnetohydrodynamic (MHD) turbulence, and discuss their connection to the spherically bounded case. We compare our numerical results with a statistical theory of geodynamo action that has evolved from the absolute equilibrium ensemble theory of ideal MHD turbulence, which is based on the ideal MHD invariants are energy, cross helicity and magnetic helicity. However, for rotating MHD turbulence, the cross helicity is no longer an exact invariant, although rms cross helicity becomes quasistationary during an ideal MHD simulation. This and the anisotropy imposed by rotation suggests an ansatz in which an effective, nonzero value of cross helicity is assigned to axisymmetric modes and zero cross helicity to non-axisymmetric modes. This hybrid statistics predicts a large-scale quasistationary magnetic field due to broken ergodicity , as well as dipole vector alignment with the rotation axis, both of which are observed numerically. We find that only a relatively small value of effective cross helicity leads to the prediction of a dipole moment vector that is closely aligned (less than 10 degrees) with the rotation axis. We also discuss the effect of initial conditions, dissipation and grid size on the numerical simulations and statistical theory.

Shebalin, John V.↗

Striving to translate shale physics across ten orders of magnitude: What have we learned?

Shales will play an important role in the successful transition of energy from fossil-based resources to renewables in the coming decades. Aside from being a significant source of low-carbon intensity fuels, like natural gas, they also serve as geologic seals of subsurface formations that may be used to isolate nuclear waste, sequester CO 2 , or store intermittent energy (e.g., solar hydrogen). Despite their importance, shales pose significant engineering and environmental challenges due to their nanoporous structure and extreme heterogeneity that spans at least ~10 orders of magnitude in spatial scale. Two challenges inhibit a system-level understanding: (1) the physics of fluid flow and phase behavior in shales are poorly understood due to the dominant molecular interactions between minerals and fluids under confinement, and (2) the apparent lack of scale separation that prevents a reliable (closed) description of the physics at any single scale of observation. In this review, we focus on the latter issue and discuss scale translation, which in its broadest sense is transforming data or simulations from one spatiotemporal scale to another. While effective scale translation is not exclusive to shales, but all geologic porous media, the need for it is especially acute in shales given their high degree of heterogeneity. Classical theories like homogenization, while indispensable, fail when scales are not separated. Other methods, like numerical upscaling, scale-translate in only one direction: small to large, but not the reverse, called downscaling. However, the confluence of advances in three areas are bringing challenging problems such as shales within reach: increased computational power and scalable algorithms; high-resolution imaging and multi-modal data acquisition; and machine learning to process massive amounts of data. While these advances equip geoscientists with a wide array of experimental and computational tools, no individual tool can probe the entire gamut of heterogeneity in shales. Their effective use, therefore, requires an ability to bridge between various data types obtained at different scales. The aim of this review is to present a coherent account of computational and experimental methods that may be used to achieve just that, i.e., to perform scale translation. We provide a broader definition of scale translation, one that transcends classical homogenization and upscaling methods, but is consistent with them and accommodates notions like downscaling and data translation. After a brief introduction to homogenization, we review hybrid methods, numerical upscaling and its recent extensions, multiscale computing, high-resolution imaging, and machine learning. We place particular emphasis on multiscale computing and propose an algorithmic framework to bridge between the pore (micro) and Darcy (macro) scales. Throughout the paper, we draw comparisons between the various methods and highlight their (often hidden) similarities, differences, benefits, and pitfalls. We finally conclude with two case studies on shales that exemplify some of the methods presented.

58 GEOSCIENCES↗

PDF turbulence modeling and DNS

The problem of time discontinuity (or jump condition) in the coalescence/dispersion (C/D) mixing model is addressed in probability density function (pdf). A C/D mixing model continuous in time is introduced. With the continuous mixing model, the process of chemical reaction can be fully coupled with mixing. In the case of homogeneous turbulence decay, the new model predicts a pdf very close to a Gaussian distribution, with finite higher moments also close to that of a Gaussian distribution. Results from the continuous mixing model are compared with both experimental data and numerical results from conventional C/D models. The effect of Coriolis forces on compressible homogeneous turbulence is studied using direct numerical simulation (DNS). The numerical method used in this study is an eight order compact difference scheme. Contrary to the conclusions reached by previous DNS studies on incompressible isotropic turbulence, the present results show that the Coriolis force increases the dissipation rate of turbulent kinetic energy, and that anisotropy develops as the Coriolis force increases. The Taylor-Proudman theory does apply since the derivatives in the direction of the rotation axis vanishes rapidly. A closer analysis reveals that the dissipation rate of the incompressible component of the turbulent kinetic energy indeed decreases with a higher rotation rate, consistent with incompressible flow simulations (Bardina), while the dissipation rate of the compressible part increases; the net gain is positive. Inertial waves are observed in the simulation results.

Hsu, A. T.↗

Direct simulation of particle dispersion in a decaying isotropic turbulence

Results of a numerical investigation of the dispersion of solid particles in decaying isotropic turbulence are presented. The 3D time-dependent velocity field of a homogeneous nonstationary turbulence is computed using the method of direct numerical simulation (DNS). The dispersion characteristics of three different solid particles (corn, copper, and glass) injected in the flow are obtained by integrating the complete equation of particle motion along the instantaneous trajectories of 22-cubed particles for each particle type, and then by performing ensemble averaging. Good agreement was achieved between the present DNS results and the measured time development of the mean-square displacement of the particles. Questions of how and why the dispersion statistics of a solid particle differ from those of its corresponding fluid point and surrounding fluid and what influences inertia and gravity have on these statistics are also discussed.

Elghobashi, S.↗

On the convergence of statistics in simulations of stationary incompressible turbulent flows

When reporting statistics from simulations of statistically stationary chaotic phenomenon, it is important to verify that the simulations are time-converged. This condition is connected with the statistical error or number of digits with which statistics can be reliably reported. In this work we consider numerical experiments of low Reynolds number incompressible homogeneous and isotropic turbulence as a model problem to investigate statistical convergence over finite simulation times. Specifically, we investigate the time integration requirements that allow meaningful reporting of the statistical error associated with finiteness of the temporal domain. We address two key questions: (1) How long should a simulation be performed in terms of large eddy time, and (2) How should the simulation time be divided among temporal windows over which a quantity of interest is estimated so that its statistical error could be reliably reported? We find that reliable reporting of statistical errors requires simulations on the order of 10 4 large eddy times, which is orders of magnitude longer than typically performed. Additionally, data post-processing should employ windows of at least ten times the large eddy time scale, with the most robust computation of statistical error of the mean requiring window sizes of an additional factor of ten. For practical simulations, we demonstrate that it is possible to estimate the statistical error within a factor of two under a less stringent condition in which a minimum of four windows with size at least ten large eddy times are used. In conclusion, our observations for homogeneous isotropic turbulence are also shown to hold in turbulent channel flow.

42 ENGINEERING↗

Coarse grained simulation of convectively driven turbulent mixing, transition, and turbulence decay

Accurate predictions with quantifiable uncertainties are essential to many practical turbulent flow applications exhibiting extreme geometrical complexity and broad ranges of length and time scales. Under-resolved computer simulations are typically unavoidable in such applications, and implicit large-eddy simulation (ILES) often becomes the effective strategy. In this work, we focus on assessing ILES initialized with well-characterized 25 6 3 homogeneous isotropic turbulence datasets generated with direct numerical simulation (DNS). ILES solutions based on the LANL xRAGE code are studied as function of resolution for 6 4 3 , 12 8 3 , 25 6 3 , and 51 2 3 grids. ILES performance of new directionally-unsplit high-order numerical hydrodynamics algorithms in xRAGE is examined in this context. Compared to the initial 25 6 3 DNS, we find longer inertial subranges and higher turbulence Reynolds number (Re) for 25 6 3 and 51 2 3 xRAGE — attributed to having linked DNS (Navier–Stokes based) solutions to nominally inviscid (higher Re) Euler based ILES solutions. For fixed 25 6 3 resolution, we find that significantly larger simulated turbulence Re can be achieved with the higher-order unsplit (vs. split) discretizations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

2024 Second Half Semi Annual Report: Modeling plasticity-mediated flow in metals with pressurized cavities

The objective is to better predict the bulk-scale mechanical behavior of porous metals that have over pressurized cavities (e.g., irradiated metals with helium bubbles) by quantifying the complex coupling among cavity aspects (e.g., size distribution, inhomogeneous overpressure values, spatial arrangement) and metal properties (e.g., rate-dependency, crystallographic lattice). This requires up-scaling local mechanical fields from the single crystal scale and will be accomplished using a homogenization approach that combines full-field numerical simulations, analytical formalisms, and physics-informed machine learning to produce symbolically-defined constitutive equations (e.g., gauge functions). These equations will satisfy the objective because they enable computationally efficient predictions that approach the accuracy of computationally expensive full-field numerical simulations, abide by theoretical requirements (e.g., conservation of energy, work conjugacy), and retain the transparency of analytical models.

36 MATERIALS SCIENCE↗

Boundary-layer effects in composite laminates: Free-edge stress singularities, part 6

A rigorous mathematical model was obtained for the boundary-layer free-edge stress singularity in angleplied and crossplied fiber composite laminates. The solution was obtained using a method consisting of complex-variable stress function potentials and eigenfunction expansions. The required order of the boundary-layer stress singularity is determined by solving the transcendental characteristic equation obtained from the homogeneous solution of the partial differential equations. Numerical results obtained show that the boundary-layer stress singularity depends only upon material elastic constants and fiber orientation of the adjacent plies. For angleplied and crossplied laminates the order of the singularity is weak in general.

Wanag, S. S.↗

Simulation and testing of digital control on a flexible beam

Large space structures are expected to have control problems due to low stiffness and damping, and control laws for these structures must deal with shape and configuration control as well as attitude and orbit maintenance. In general, these control tasks must be accomplished without adversely interacting with the lightly damped and low frequency vibration modes of the structure. Modal control schemes have been proposed to deal with these problems. A discrete time parameter adaptive control scheme which uses modal control has been proposed by Montgomery and Johnson (1978). In the present investigation the method considered by Montgomery and Johnson is applied to a homogeneous free-free beam in both numerical simulation and laboratory experimentation. Mathematical modeling of the beam is treated in a manner expedient for digital simulation and control implementation.

Williams, J. P.↗

Study of Lagrangian characteristic times using direct numerical simulation of turbulence

Direct numerical simulations and large eddy simulations of homogeneous isotropic turbulence are used to compute Lagrangian statistics of turbulence and, in particular, its time scales. The computed time scales are compared with the spectral times scales that are frequently used in Eddy Damped Quasi-Normal Markovian calculations of the spectrum. The time scale models are rather good at high wavenumber and the results point to directions for improvement of the time scales at low wavenumber.

Lee, C. H.↗

Direct numerical simulations of stably-stratified sheared turbulence: Implications for oceanic mixing

Direct numerical simulations of the time evolution of homogeneous stably stratified turbulent shear flows have been performed for several Richardson numbers Ri and Reynolds numbers R(sub lambda) in earlier works. The results show excellent agreement with length scale models developed from laboratory experiments to characterize oceanic turbulence. When the Richardson number Ri is less than the stationary value Ri(sub s), the turbulence intensity grows at all scales, and the growth rate appears to be a function of Ri. The size of the vertical density inversions also increases. On the other hand, when Ri is greater than or equal to Ri(sub s) the largest turbulent eddies become vertically constrained by buoyancy when the Ellison (turbulence) scale L(sub E) and the Ozmidov (buoyancy) scale L(sub O) are equal. At this point, the mixing efficiency is maximal and corresponds to a flux Richardson number R(sub f) = 0.20. The vertical mass flux becomes counter-gradient when epsilon = 19(nu)N(exp 2) and vertical density overturns are suppressed in less than half a Brunt-Vaisala period. The results of the simulations were also recast in terms of the Hydrodynamic Phase Diagram introduced in fossil turbulence models. The so-called point of fossilization occurs when epsilon = 4DCN(exp 2); Gibson proposed 13DCN(exp 2). This value is in agreement with indirect laboratory observations and field observations. Finally, the validity of the steady-state models to estimate vertical eddy diffusivities in the oceanic thermocline is discussed.

Itsweire, E. C.↗

A dynamic subgrid-scale model for compressible turbulence and scalar transport

The dynamic subgrid-scale (SGS) model of Germano et al. (1991) is generalized for the large eddy simulation (LES) of compressible flows and transport of a scalar. The model was applied to the LES of decaying isotropic turbulence, and the results are in excellent agreement with experimental data and direct numerical simulations. The expression for the SGS turbulent Prandtl number was evaluated using direct numerical simulation (DNS) data in isotropic turbulence, homogeneous shear flow, and turbulent channel flow. The qualitative behavior of the model for turbulent Prandtl number and its dependence on molecular Prandtl number, direction of scalar gradient, and distance from the wall are in accordance with the total turbulent Prandtl number from the DNS data.

Moin, P.↗

Two-temperature models of old supernova remnants with ion and electron thermal conduction

To investigate the potential effects thermal conduction may have on the evolution of old supernova remnants, we present the results of 1D (spherically symmetric) numerical simulations of a remnant in a homogeneous interstellar medium for four different cases: (1) without thermal conduction; (2) with both electron and ion thermal conduction assuming equal temperatures; (3) with electron thermal conduction only, following electron and ion temperatures separately; and (4) with both electron and ion thermal conduction following separate temperatures. We followed the entire evolution until the completion of the remnant bubble collapse. Our most significant result is that in remnant evolution studies concerned principally with either the shell or bubble evolution at late times, reasonable results are obtained with single-temperature models. When the electron and ion temperatures are followed separately, however, ion thermal conduction cannot safely be ignored.

Cui, Wei↗

Pseudospectral simulation of compressible turbulence using logarithmic variables

The direct numerical simulation of dissipative, highly compressible turbulent flow is performed using a pseudospectral Fourier technique. The governing equations are cast in a form where the important physical variables are the fluid velocity and the natural logarithms of the fluid density and temperature. Bulk viscosity is utilized to model polyatomic gases more accurately and to ensure numerical stability in the presence of strong shocks. Numerical examples include three-dimensional supersonic homogeneous turbulence and two-dimensional shock-turbulence interactions.

Shebalin, John V.↗

Mixing in a stratified shear flow: Energetics and sampling

Direct numerical simulations of the time evolution of homogeneous stably stratified shear flows have been performed for Richardson numbers from 0 to 1 and for Prandtl numbers between 0.1 and 2. The results indicate that mixing efficiency R(sub f) varies with turbulent Froude number in a manner consistent with laboratory experiments performed with Prandtl numbers of 0.7 and 700. However, unlike the laboratory results, for a particular Froude number, the simulations do not show a clear dependence on the magnitude of R(sub f) on Pr. The observed maximum value of R(sub f) is 0.25. When averaged over vertical length scales of an order of magnitude greater than either the overturning or Ozmidov scales of the flow, the simulations indicate that the dissipation rate epsilon is only weakly lognormally distributed with an intermittency of about 0.01 whereas estimated values in the ocean are 3 to 7.

Ivey, G. N.↗