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.

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Reynolds stress calculations of homogeneous turbulent shear flow with bounded energy states

Reynolds stress calculations of homogeneous turbulent shear flow are conducted with a second-order closure model modified to account for non-equilibrium vortex stretching in the dissipation rate transport equation, as recently proposed by Bernard and Speziale. As with the earlier reported k-epsilon model calculations incorporating this vortex stretching effect, a production-equals-dissipation equilibrium is obtained with bounded turbulent kinetic energy and dissipation. However, this equilibrium is not achieved until the dimensionless time greater than 60, an elapsed time that is at least twice as large as any of those considered in previous numerical and physical experiments on homogeneous shear flow. Direct quantitative comparisons between the model predictions and the results of experiments are quite favorable. In particular, it is shown that the inclusion of this non-equilibrium vortex stretching effect has the capability of explaining the significant range of production to dissipation ratios observed in experiments.

Speziale, Charles G.↗

NASA NTRS1993

Reynolds stress calculations of homogeneous turbulent shear flow with bounded energy states

Reynolds stress calculations of homogeneous turbulent shear flow are conducted with a second-order closure model modified to account for nonequilibrium vortex stretching in the dissipation rate transport equation as recently proposed by Bernard and Speziale (1992). As with the earlier reported K-epsilon model calculations incorporating this vortex stretching effect, a production-equals-dissipation equilibrium is obtained with bounded turbulent kinetic energy and dissipation. However, this equilibrium is now not achieved until the dimensionless time St greater than 60 - an elapsed time that is at least twice as large as any of those considered in previous numerical and physical experiments on homogeneous shear flow. Direct quantitative comparisons between the model predictions and the results of experiments are quite favorable. In particular, it is shown that the inclusion of this nonequilibrium vortex stretching effect has the capability of explaining the significant range of production to dissipation ratios observed in experiments.

Speziale, Charles G.↗

NASA NTRS1978

Numerical solutions and laser-Doppler measurements of spin-up

The spin-up flow in a cylinder of homogeneous fluid has been examined both experimentally and numerically. A series of laser-Doppler measurements was made of the zonal flow over a range of Ekman numbers and Rossby numbers at various locations in the interior of the flow. These measurements exceed previous ones in accuracy. The weak inertial modes excited by the impulsive start are detectable. The numerical simulations used the primitive equations in axisymmetric form and employed finite-difference techniques on both constant and variable grids. The number of grid points necessary to resolve the Ekman layers was determined. A thorough comparison of the simulations and the experimental measurements is made which includes the details of the amplitude and frequency of the inertial modes. Agreement to within the experimental tolerance is achieved. Analytical results for conditions identical to those in the experiments are not available but some similar linear and nonlinear theories are also compared with the experiments.

Warn-Varnas, A.↗

DOE OSTI2022

Implications of inertial subrange scaling for stably stratified mixing

In this work, we investigate the effects of the turbulent dynamic range on active scalar mixing in stably stratified turbulence by adapting the theoretical passive scalar modelling and demonstrating their usefulness through consideration of the results of direct numerical simulations of statistically stationary homogeneous stratified and sheared turbulence. By analysis of inertial and inertial–convective subrange scalings, we show that the relationship between the active scalar and turbulence time scales is predicted by the ratio of the Kolmogorov and Obukhov–Corrsin constants, provided mean flow parameters permit the two subrange scalings to be appropriate approximations. We use the resulting relationship between time scales to parameterise an appropriate turbulent mixing coefficient $\varGamma \equiv \chi /\epsilon$, defined here as the ratio of available potential energy ($E_p$) and turbulent kinetic energy ($E_k$) dissipation rates. With the analysis presented here, we show that $\varGamma$ can be estimated by $E_p,E_k$ and a universal constant provided an appropriate Reynolds number is sufficiently high. This large Reynolds number regime appears here to occur at $ {{Re_b}} \equiv \epsilon / \nu N^{2} \gtrapprox 300$ where $\nu$ is the kinematic viscosity and $N$ is the characteristic buoyancy frequency. We propose a model framework for irreversible diapycnal mixing with robust theoretical parametrisation and asymptotic behaviour in this high-$ {{Re_b}}$ limit.

42 ENGINEERING↗

DOE OSTI2024

A high-frequency, low-power resonant radio-frequency neutron spin flipper for high-resolution spectroscopy

Here, we present a resonant-mode, transverse-field, radio-frequency (rf) neutron spin flipper design that uses high-temperature superconducting films to ensure sharp transitions between uniform magnetic field regions. Resonant mode allows for low-power, high-frequency operation but requires strict homogeneity of the magnetic fields inside the device. This design was found to efficiently flip neutrons at 96.6 ± 0.6% at an effective frequency of 4 MHz in bootstrap configuration with a beam size of 2.4 × 2.5 cm 2 and a wavelength of 0.4 nm. The high frequency and efficiency enable this device to perform high-resolution neutron spectroscopy with comparable performance with currently implemented rf flipper designs. The limitation of the maximum frequency was found due to the field homogeneity of the device. We numerically analyze the maximum possible efficiency of this design using a Bloch solver simulation with magnetic fields generated from finite-element simulations. We also discuss future improvements of the efficiency and frequency to the design based on the experimental and simulation results.

47 OTHER INSTRUMENTATION↗

NASA NTRS1972

Some spatially homogeneous gas motions

A functional equation is derived for a spatially homogeneous distribution function and examples of spatially homogeneous gas motions are considered. Numerical calculation was conducted for one such motion with the use of a model equation.

Rykov, V. A.↗

NASA NTRS1984

Electromagnetic ion beam instabilities

The linear theory of electromagnetic instabilities driven by an energetic ion beam streaming parallel to a magnetic field in a homogeneous Vlasov plasma is considered. Numerical solutions of the full dispersion equation are presented. At propagation parallel to the magnetic field, there are four distinct instabilities. A sufficiently energetic beam gives rise to two unstable modes with right-hand polarization, one resonant with the beam, the other nonresonant. A beam with sufficiently large T (perpendicular to B)/T (parallel to B) gives rise to the left-hand ion cyclotron anisotropy instability at relatively small beam velocities, and a sufficiently hot beam drives unstable a left-hand beam resonant mode. The parametric dependences of the growth rates for the three high beam velocity instabilities are presented here. In addition, some properties at oblique propagation are examined. It is demonstrated that, as the beam drift velocity is increased, relative maxima in growth rates can arise at harmonics of the ion cyclotron resonance for both right and left elliptically polarized modes.

Gary, S. P.↗

NASA NTRS1985

Morphological Stability and Fluid Dynamics of Vapor Crystal Growth

A fundamental understanding of the conditions under which crystals can retain morphological stability, i.e., shape stability of the advancing interface, during growth from vapors was studied. Morphological stability (MS) is a necessary condition for the growth of homogeneous single crystals required for numerous device applications. For crystallization from melts, the MS concepts are well developed and are essentially based on heat and mass transfer conditions about the advancing interface. For crystallization from vapors, the MS requirements are more complex and not well understood. The added complexity arises from the fact that anisotropies in interfacial kinetics are typically stronger in crystallization from vapors than from melts. These pronounced anisotropies root in the distinctly lower atomic roughness of most vapor-solid interfaces.

Rosenberger, F.↗

NASA NTRS1989

The analysis and modeling of dilatational terms in compressible turbulence

It is shown that the dilatational terms that need to be modeled in compressible turbulence include not only the pressure-dilatation term but also another term - the compressible dissipation. The nature of these dilatational terms in homogeneous turbulence is explored by asymptotic analysis of the compressible Navier-Stokes equations. A non-dimensional parameter which characterizes some compressible effects in moderate Mach number, homogeneous turbulence is identified. Direct numerical simulations (DNS) of isotropic, compressible turbulence are performed, and their results are found to be in agreement with the theoretical analysis. A model for the compressible dissipation is proposed; the model is based on the asymptotic analysis and the direct numerical simulations. This model is calibrated with reference to the DNS results regarding the influence of compressibility on the decay rate of isotropic turbulence. An application of the proposed model to the compressible mixing layer has shown that the model is able to predict the dramatically reduced growth rate of the compressible mixing layer.

Sarkar, S.↗

NASA NTRS1991

The analysis and modelling of dilatational terms in compressible turbulence

Sarkar, S.↗

NASA NTRS1990

Modeling the turbulent kinetic energy equation for compressible, homogeneous turbulence

The turbulent kinetic energy transport equation, which is the basis of turbulence models, is investigated for homogeneous, compressible turbulence using direct numerical simulations performed at CTR. It is shown that the partition between dilatational and solenoidal modes is very sensitive to initial conditions for isotropic decaying turbulence but not for sheared flows. The importance of the dilatational dissipation and of the pressure-dilatation term is evidenced from simulations and a transport equation is proposed to evaluate the pressure-dilatation term evolution. This transport equation seems to work well for sheared flows but does not account for initial condition sensitivity in isotropic decay. An improved model is proposed.

Aupoix, B.↗

NASA NTRS1992

The mirror and ion cyclotron anisotropy instabilities

The linear dispersion equation for fully electromagnetic waves and instabilities at arbitrary directions of propagation relative to a background magnetic field B(0) in a homogeneous Vlasov plasma is solved numerically for bi-Maxwellian particle distributions. For isotropic plasmas the dispersion and damping of the three modes below the proton cyclotron frequency are studied as functions of Beta(i) and T(e)/T(i). The transport ratios of helicity, cross-helicity, Alfven ratio, compressibility, and parallel compressibility are defined. Under the condition that the proton temperature perpendicular to B(0) is greater than the parallel temperature, the growth rates and transport ratios of the mirror instability and the ion cyclotron anisotropy instability are examined and compared. Both the proton parallel compressibility and the proton Alfven ratio are significantly different for the two growing modes.

Gary, S. P.↗

NASA NTRS1995

A Study of Two-Equation Turbulence Models on the Elliptic Streamline Flow

Several two-equation turbulence models are compared to data from direct numerical simulations (DNS) of the homogeneous elliptic streamline flow, which combines rotation and strain. The models considered include standard two-equation models and models with corrections for rotational effects. Most of the rotational corrections modify the dissipation rate equation to account for the reduced dissipation rate in rotating turbulent flows, however, the DNS data shows that the production term in the turbulent kinetic energy equation is not modeled correctly by these models. Nonlinear relations for the Reynolds stresses are considered as a means of modifying the production term. Implications for the modeling of turbulent vortices will be discussed.

Blaisdell, Gregory A.↗

NASA NTRS1998

FU Orionis Outbursts and the Solar Nebula

Protostellar systems are variable on many timescales. The FU Orionis outburst is one of the most drastic forms of variability known to occur in low mass stellar systems. During a typical outburst lasting several decades, system luminosities may be a hundred times what is normal of the quiescent state. FU Orionis outburst events are believed to have significant impacts on the thermal structure of the protosolar nebula. Their existence has been utilized to explain features in the meteoritic record from thermally induced homogenization to chondrule formation. Recent numerical models have shown the viability of the hypothesis that the radiation observed during outburst is emitted by a luminous circumstellar disk transporting mass at a thousand times the quiescent rate. We will begin by describing what is known about the FU Orionis outburst phenomenon from recent observations and theory. We will discuss evidence that suggests that outburst radiation is emitted by a circumstellar disk rather than by the star and will briefly describe the thermal instability as a mechanism for outburst. Additional information is contained in the original extended abstract.

Bell, Robbins↗

NASA NTRS1981

Mixing of a passive scalar in isotropic and sheared homogeneous turbulence

In order to calculate the velocity and scalar fields, the three dimensional, time-dependent equations of motion and the diffusion equation were solved numerically. The following cases were treated: isotropic, homogeneous turbulence with decay of a passive scalar; and homogeneous turbulent shear flow with a passive scalar whose mean varies linearly in the spanwise direction. The solutions were obtained at relatively low Reynolds numbers so that all of the turbulent scales could be resolved without modeling. Turbulent statistics such as integral length scales, Taylor microscales, Kolmogorov length scale, one- and two-point correlations of velocity-velocity and velocity-scalar, turbulent Prandtl/Schmidt number, r.m.s. values of velocities, the scalar quantity and pressure, skewness, decay rates, and decay exponents were calculated. The results are compared with the available expermental results, and good agreement is obtained.

Shirani, E.↗

DOE OSTI2022

DNS of Hydrodynamic Instabilities of Laminar H2/O2/N2 Flames at Elevated Pressure

The atypical combustion properties of hydrogen have long provided an interest for hydrogen-enriched combustion, emphasized in recent years by increasingly stringent emission regulations. Significant work has been focused on the small hydrocarbon(C1-C3)/hydrogen blends, but much less is known about the effect of hydrogen on heavier hydrocarbon combustion, of interest for the aviation and automotive industry. In this work, we perform direct numerical simulations (DNS) of spherically expanding laminar n-dodecane/H2 flames in a constant-volume vessel at elevated pressure and a range of hydrogen seeding levels. A low Mach number model is used to represent the gas and flame dynamics allowing for a temporally varying, spatially homogeneous pressure field, while also enabling numerical evolution of the system numerically at time step sizes governed by the advective CFL limit rather than acoustic processes. We use an adaptive mesh refinement (AMR) approach to tackle the large separation of scales associated with high pressure premixed flame surfaces propagating in a laboratory-scale closed domain. We incorporate detailed chemistry and transport models for the reacting flow and focus the analysis on the effects of the pressure rise on the flame propagation characteristics and morphology as the hydrogen content is increased. The presence of fast diffusive H2 and the pressure rise results in dramatic thinning of the flame and the potential onset of hydrodyanmic and theromdiffusive instabilities. We assess the potential role of these instabilities on the ability to accurately measure laminar burning speed experimentally based on the mean flame surface propagation speed.

high pressure↗

DOE OSTI2020

Diffusional growth of cloud droplets in homogeneous isotropic turbulence: DNS, scaled-up DNS, and stochastic model

This paper presents a novel methodology to use direct numerical simulation (DNS) to study the impact of isotropic homogeneous turbulence on the condensational growth of cloud droplets. As shown by previous DNS studies, the impact of turbulence increases with the computational domain size, that is, with the Reynolds number, because larger eddies generate higher and longer-lasting supersaturation fluctuations that affect growth of individual cloud droplets. The traditional DNS can only simulate a limited range of scales because of the excessive computational cost that comes from resolving all scales involved, that is, from large scales at which the turbulent kinetic energy (TKE) is introduced down to the Kolmogorov microscale, and from following every single droplet. The novel approach is referred to as the “scaled-up DNS”. The scaling up is done in two parts, first by increasing both the computational domain and the Kolmogorov microscale and second by using super-droplets instead of real droplets. To ensure proper dissipation of TKE and scalar variance at small scales, molecular transport coefficients are appropriately scaled up with the grid length. For the scaled-up domains, say, meters and tens of meters, one needs to follow billions of real droplets. This is not computationally feasible, and so-called super-droplets are applied in scaled-up DNS simulations. Each super-droplet represents an ensemble of identical real droplets, and the number of real droplets represented by a super-droplet is referred to as the multiplicity attribute. After simple tests showing the validity of the methodology, scaled-up DNS simulations are conducted for five domains, the largest of 64 3 m 3 volume using a DNS of 256 3 grid points and various multiplicities. All simulations are carried out with vanishing mean vertical velocity and with no mean supersaturation, similarly to past DNS studies. As expected, the supersaturation fluctuations as well as the spread in droplet size distribution increase with the domain size, with the droplet radius variance increasing in time t as t 1∕2 as identified in previous DNS studies. Scaled-up simulations with different multiplicities document numerical convergence of the scaled-up solutions. Finally, we compare the scaled-up DNS results with a simple stochastic model that calculates supersaturation fluctuations based on the vertical velocity fluctuations updated using the Langevin equation. Overall, the results document similar scaling to previous small-domain DNS simulations and support the notion that the stochastic subgrid-scale model is a valuable tool for the multi-scale simulation of droplet spectral evolution applying a large-eddy simulation model.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Floquet operator dynamics and orthogonal polynomials on the unit circle

Operator spreading under stroboscopic time evolution under a unitary is studied. An operator Krylov space is constructed and mapped to orthogonal polynomials on a unit circle (OPUC), as well as to the Krylov space of the edge operator of the Floquet transverse field Ising model with inhomogeneous couplings (ITFIM). The Verblunsky coefficients in the OPUC representation are related to the Krylov angles parametrizing the ITFIM. The relations between the OPUC and spectral functions are summarized and several applications are presented. These include derivation of analytic expressions for the OPUC under persistent m-periodic dynamics, and the numerical construction of the OPUC for autocorrelations of the homogeneous Floquet-Ising model as well as the Z 3 clock model. The numerically obtained Krylov angles of the Z 3 clock model with long-lived period tripled autocorrelations show a spatial periodicity of six, and this observation is used to develop an analytically solvable model for the ITFIM that mimics this behavior.

Floquet systems↗