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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Electromagnetic ion beam instabilities - Hot beams at interplanetary shocks

This paper considers Vlasov instabilities driven by a very hot ion beam streaming parallel to a magnetic field. The linear theory of electromagnetic instabilities driven by such a beam in a homogeneous, collisionless, nonrelativistic plasma is used. Numerical solutions of the full linear dispersion equation for bi-Maxwellian distribution functions are presented. If the thermal speed of the beam is greater than its drift speed, the two dominant modes are the right-hand and left-hand resonant ion beam instabilities. The parametric dependencies of the threshold drift speeds and growth rates for both modes are presented. In particular, the thresholds are shown to be sensitive functions of the beam temperature anisotropies. Arguments are presented that these two modes are driven by the suprathermal ion component at interplanetary shocks, and that the growth of these modes is sufficient to account for the amplitudes of MHD-like waves observed at such shocks.

Gary, S. P.↗
NASA NTRS1987

Analysis of subgrid-scale eddy viscosity with use of results from direct numerical simulations

Without resort to any modeling, subgrid-scale eddy viscosity is computed from the results of high-resolution (64 cubed and 128 cubed grid points) direct numerical simulations of three-dimensional homogeneous isotropic decaying turbulence. In these simulations the eddy viscosity peaks sharply at the cutoff wavenumber, in rough agreement with the results of Kraichnan (1976). In addition, in the low-wavenumber range the eddy viscosity may be negative, contrary to the generally accepted concept of subgrid-scale eddy viscosity. Some possible explanations of this behavior are discussed.

Domaradzki, J. Andrzej↗
NASA NTRS1989

Multifractal spectra in shear flows

Numerical simulations of three-dimensional homogeneous shear flow and fully developed channel flow, are used to calculate the associated multifractal spectra of the energy dissipation field. Only weak parameterization of the results with the nondimensional shear is found, and this only if the flow has reached its asymptotic development state. Multifractal spectra of these flows coincide with those from experiments only at the range alpha less than 1.

Keefe, L. R.↗
NASA NTRS1991

Analysis of turbulence in the orthonormal wavelet representation

The usefulness of the wavelet transform for the analysis of turbulent flow fields is explored by examining the wavelet transform properties of a decomposition of turbulent velocity fields into modes that exhibit the localization in a wavenumber and physical space. The calculations are performed on 3D fields from direct numerical simulations of isotropic flow and homogeneous shear flow, and from measurements in two laboratory wind-tunnel experimental velocity signals (boundary layer and wake behind a circular cylinder). The analysis confirmed that there is strong spatial intermittency in nonlinear quantities; their mean spectral behavior results from a delicate balance between large positive and negative excursions. The wavelet analysis is a way to quantify these observations in a standardized fashion by using 'flow-independent eddies' to decompose the velocity field.

Meneveau, Charles↗
NASA NTRS1974

Finite elasto-plastic deformation. I - Theory and numerical examples

It is demonstrated that the problem of elasto-plastic finite deformation is governed by a quasi-linear model irrespective of deformation magnitude. This feature follows from the adoption of a rate viewpoint toward finite deformation analysis in an Eulerian reference frame. Objectivity of the formulation is preserved by introduction of a frame-invariant stress rate. Equations for piece-wise linear incremental finite element analysis are developed by application of the Galerkin method to the instantaneously linear governing differential equations of the quasi-linear model. Numerical solution capability has been established for problems of plane strain and plane stress. The accuracy of the numerical analysis is demonstrated by consideration of a number of problems of homogeneous finite deformation admitting comparative analytic solution. It is shown that accurate, objective numerical solutions can be obtained for problems involving dimensional changes of an order of magnitude and rotations of a full radian.

Osias, J. R.↗
NASA NTRS1993

Homogeneous ice nucleation and supercooled liquid water in orographic wave clouds

This study investigates ice nucleation mechanisms in cold lenticular wave clouds, a cloud type characterized by quasi-steady-state air motions and microphysical properties. It is concluded that homogeneous ice nucleation is responsible for the ice production in these clouds at temperatures below about -33 C. The lack of ice nucleation observed above -33 C indicates a dearth of ice-forming nuclei, and hence heterogeneous ice nucleation, in these clouds. Aircraft measurements in the temperature range -31 to -41 C show the following complement of simultaneous and abrupt changes in cloud properties that indicate a transition from the liquid phase to ice: disappearance of liquid water; decrease in relative humidity from near water saturation to ice saturation; increase in mean particle size; change in particle concentration; and change in temperature due to the release of latent heat. A numerical model of cloud particle growth and homogeneous ice nucleation is used to aid in interpretation of our in situ measurements. The abrupt changes in observed cloud properties compare favorably, both qualitatively and quantitatively, with results from the homogeneous ice nucleation model. It is shown that the homogeneous ice nucleation rates from the measurements are consistent with the temperature-dependent rates employed by the model (within a factor of 100, corresponding to about 1 C in temperature) in the temperature range -35 deg to -38 C. Given the theoretical basis of the modeled rates, it may be reasonable to apply them throughout the -30 to -50 C temperature range considered by the theory.

Heymsfield, Andrew J.↗
NASA NTRS2006

Homogeneous Aerosol Freezing in the Tops of High-Altitude Tropical Cumulonimbus Clouds

Numerical simulations of deep, intense continental tropical convection indicate that when the cloud tops extend more than a few kilometers above the liquid water homogeneous freezing level, ice nucleation due to freezing of entrained aqueous sulfate aerosols generates large concentrations of small crystals (diameters less than approx. equal to 20 micrometers). The small crystals produced by aerosol freezing have the largest impact on cloud-top ice concentration for convective clouds with strong updrafts but relatively low aerosol concentrations. An implication of this result is that cloud-top ice concentrations in high anvil cirrus can be controlled primarily by updraft speeds in the tops of convective plumes and to a lesser extent by aerosol concentrations in the uppermost troposphere. While larger crystals precipitate out and sublimate in subsaturated air below, the population of small crystals can persist in the saturated uppermost troposphere for many hours, thereby prolonging the lifetime of remnants from anvil cirrus in the tropical tropopause layer.

Jensen, E. J.↗
NASA NTRS1984

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.↗
NASA NTRSDate not supplied

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↗
NASA NTRSDate not supplied

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↗
NASA NTRS2012

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.↗
NASA NTRS1992

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.↗
NASA NTRS1992

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.↗
NASA NTRS1981

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.↗
NASA NTRS1982

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.↗
NASA NTRS1987

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.↗
NASA NTRS1990

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.↗