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 199 records · Page 11

Simulation and modeling of homogeneous, compressed turbulence

Low Reynolds number homogeneous turbulence undergoing low Mach number isotropic and one-dimensional compression was simulated by numerically solving the Navier-Stokes equations. The numerical simulations were performed on a CYBER 205 computer using a 64 x 64 x 64 mesh. A spectral method was used for spatial differencing and the second-order Runge-Kutta method for time advancement. A variety of statistical information was extracted from the computed flow fields. These include three-dimensional energy and dissipation spectra, two-point velocity correlations, one-dimensional energy spectra, turbulent kinetic energy and its dissipation rate, integral length scales, Taylor microscales, and Kolmogorov length scale. Results from the simulated flow fields were used to test one-point closure, two-equation models. A new one-point-closure, three-equation turbulence model which accounts for the effect of compression is proposed. The new model accurately calculates four types of flows (isotropic decay, isotropic compression, one-dimensional compression, and axisymmetric expansion flows) for a wide range of strain rates.

Wu, C. T.↗

Rapid calculation of radiative heating rates and photodissociation rates in inhomogeneous multiple scattering atmospheres

The solution of the generalized two-stream approximation for radiative transfer in homogeneous multiple scattering atmospheres is extended to vertically inhomogeneous atmospheres in a manner which is numerically stable and computationally efficient. It is shown that solar energy deposition rates, photolysis rates, and infrared cooling rates all may be calculated with the simple modifications of a single algorithm. The accuracy of the algorithm is generally better than 10 percent, so that other uncertainties, such as in absorption coefficients, may often dominate the error in calculation of the quantities of interest to atmospheric studies.

Toon, Owen B.↗

Thermoelectric effects in decaying homogeneous magneto-gas turbulence

In the formulation of compressible MHD (i.e., magneto-gas dynamics), a 'generalized Ohm's law is required. In particular, an electron pressure term and a 'Hall effect' term may appear as non-negligible additions to the Ohm's law that is conventionally used for incompressible MHD. In 'high-beta' (i.e., relatively low magnetic energy) situations, the Hall term may be neglected (at least initially) but, as it turns out, the electron pressure term cannot be neglected. Here, three-dimensional, high-beta, homogeneous, decaying, magneto-gas turbulence is examined with regard to this additional term. Through numerical simulation, it is found that 'thermoelectric effects' are produced that significantly alter the evolution of the magnetic field and electric current strengths.

Shebalin, John V.↗

Effects of Eddy Viscosity on Time Correlations in Large Eddy Simulation

Subgrid-scale (SGS) models for large. eddy simulation (LES) have generally been evaluated by their ability to predict single-time statistics of turbulent flows such as kinetic energy and Reynolds stresses. Recent application- of large eddy simulation to the evaluation of sound sources in turbulent flows, a problem in which time, correlations determine the frequency distribution of acoustic radiation, suggest that subgrid models should also be evaluated by their ability to predict time correlations in turbulent flows. This paper compares the two-point, two-time Eulerian velocity correlation evaluated from direct numerical simulation (DNS) with that evaluated from LES, using a spectral eddy viscosity, for isotropic homogeneous turbulence. It is found that the LES fields are too coherent, in the sense that their time correlations decay more slowly than the corresponding time. correlations in the DNS fields. This observation is confirmed by theoretical estimates of time correlations using the Taylor expansion technique. Tile reason for the slower decay is that the eddy viscosity does not include the random backscatter, which decorrelates fluid motion at large scales. An effective eddy viscosity associated with time correlations is formulated, to which the eddy viscosity associated with energy transfer is a leading order approximation.

He, Guowei↗

Easing into Reality: Experimental Impacts into Slopes and Layers

Impact cratering is the dominant geo-logic process affecting the surfaces of solid bodies throughout our solar system. Because large impacts are (luckily) rare on Earth, the process is studied through experiments, observations of existing structures, numerical modeling, and theory, most of which make the simplifying assumptions that the target is homogeneous, with no substantial topography. Craters do not always form on level targets com-posed of homogeneous loose material. Rather (Fig. 1), they often form on sloped surfaces and in layered tar-gets, both of which significantly influence the excavation and ejecta deposition processes. Such craters are common on the Moon and asteroids. We are investigating crater formation in two separate suites of experiments using sloped and layered targets (Fig. 2) at the Experimental Impact Laboratory at NASA Johnson Space Center. An experiment was also performed in a flat, homogenous target to serve as a reference.

Ebel, J. M.↗

An Analytical Method to Calculate Effective Elastic Properties of General Multifunctional Honeycomb Cores in Sandwich Composites

Sandwich composite structures are ideal configurations in which to incorporate additional functionality beyond load-carrying capabilities. The inner core-walls can be layered to incorporate other functions such as power storage for a battery. In this work we investigate an assemblage of analytical tools to compute effective properties that allow complex layered core architectures to be homogenized into a single continuum layer. This provides a great increase in computational efficiency to numerically simulate the structural response of multifunctional sandwich structures under applied loads. We present a coupled analytical method including an extensive numerical verification of the accuracy of this method.

Sather, Erik↗

Estimation of Effective Elastic Properties of General Multifunctional Honeycomb Structures Using a Unit Cell Method

Sandwich composite structures are ideal configurations in which to incorporate additional functionality beyond load carrying capabilities. The inner core can be layered to facilitate other functions such as power storage for a battery. In this work we investigate an assemblage of analytical tools to compute effective properties that allow complex layered core architectures to be homogenized into a single continuum layer. This provides a great increase in computational efficiency to numerically simulate the structural response of multifunctional sandwich structures under applied loads.

Krishnamurthy, T.↗

An investigation of voids formation mechanisms and their effects on freeze and thaw processes of lithium and lithium fluoride

The mechanisms of void formation during the cooldown and freezing of lithium coolant within the primary loop of SP-100 type systems are investigated. These mechanisms are: (1) homogeneous nucleation; (2) heterogeneous nucleation; (3) normal segregation of helium gas dissolved in liquid lithium; and (4) shrinkage of lithium during freezing. To evaluate the void formation potential due to segregation, a numerical scheme that couples the freezing and mass diffusion processes in both the solid and liquid regions is developed. The results indicated that the formation of He bubbles is unlikely by either homogeneous or heterogeneous nucleation during the cooldown process. However, homogeneous nucleation of He bubbles following the segregation of dissolved He in liquid lithium ahead of the solid-liquid interface is likely to occur. Results also show that total volume of He void is insignificant when compared to that of shrinkage voids. In viewing this, the subsequent research focuses on the effects of shrinkage void forming during freezing of lithium on subsequent thaw processes are investigated using a numerical scheme that is based on a single (solid/liquid) cell approach. The cases of lithium-fluoride are also investigated to show the effect of larger volume shrinkage upon freezing on the freeze and thaw processes. Results show that a void forming at the wall appreciably reduces the solid-liquid interface velocity, during both freeze and thaw, and causes a substantial rise in the wall temperature during thaw. However, in the case of Li, the maximum wall temperature was much lower than the melting temperature of PWC-11, which is used as the structure material in the SP-100 system. Hence, it is included that a formation of hot spots is unlikely during the startup or restart of the SP-100 system.

El-Genk, Mohamed S.↗

Modeling the Dynamics of Ensemble-Averaged Linear Disturbances in Homogeneous Shear Flow

In order to expand the predictive capability of single-point turbulence closure models to account for the early-stage transition regime, a methodology for the formulation and calibration of model equations for the ensemble-averaged disturbance kinetic energy and energy dissipation rate is presented. The calibration is based on homogeneous shear flow where disturbances can be described by rapid distortion theory (RDT). The relationship between RDT and linear stability theory is exploited in order to obtain a closed set of modeled equations. The linear disturbance equations are solved directly so that the numerical simulation yields a database from which the closure coefficients in the ensemble-averaged disturbance equations can be determined. Key words: transition/turbulence modeling, homogeneous shear, rapid distortion theory.

Thacker, W. D.↗

Condensation of water vapor in rarefaction waves. I - Homogeneous nucleation

A detailed theoretical investigation has been made of the condensation of water vapor/carrier gas mixtures in the nonstationary rarefaction wave generated in a shock tube. It is assumed that condensation takes place by homogeneous nucleation. The equations of motion together with the nucleation rate and the droplet growth equations were solved numerically by the method of characteristics and Lax's method of implicit artificial viscosity. It is found that, for the case considered, the condensation wave formed by the collapse of the metastable nonequilibrium state is followed by a shock wave generated by the intersection of characteristics of the same family. The expansion is practically isentropic up to the onset of condensation. The condensation front accelerates in the x,t plane. The results of the computations for a chosen case of water vapor/nitrogen mixture are presented by plotting variations of pressure, nucleation rate, number density of critical clusters, and condensate mass-fraction along three particle paths. Some consideration is given to homogeneous condensation experiments conducted in a shock tube. Although a direct comparison of the present theoretical work and these experiments is not possible, several worthwhile interpretative features have resulted nevertheless.

Sislian, J. P.↗

A nonlocal theory of an electrostatic sinusoidal density drift instability

The stability of space plasmas in which the macroscopic variables vary in the direction x perpendicular to the magnetic field is usually studied by means of the local approximation, in which the fluctuating potential is assumed to be independent of x. To remove this approximation, a nonlocal problem is studied in which the background density, and hence the fluctuating potentials, are periodic in x. From the linear Vlasov/Poisson equations, a set of coupled, linear, homogeneous, algebraic equations relating the Fourier amplitudes of the eigensolution is derived. The equations are solved numerically and a hierarchy of exact eigenmodes characterized by different growth rates and spatial structures is found. At points where the density gradient is locally zero, the mode amplitude is generally several orders of magnitude lower than the peak amplitude. The dependence of these solutions on the parameters which define the background density perturbation is studied and a correspondence between these nonlocal results and those obtained with the local approximation is demonstrated.

Thomsen, M. F.↗

Robot arm dynamic model reduction for control

General methods are described by which the mathematical complexities of explicit and exact state equations of robot arms can be reduced to a simplified and compact state equation representation without introducing significant errors into the robot arm dynamic model. The model reduction methods are based on homogeneous coordinates and on the Langrangian algorithm for robot arm dynamics, and utilize matrix, vector and numeric analysis techniques. The derivation of differential vector representation of centripetal and Coriolis forces which has not yet been established in the literature is presented.

Bejczy, A. K.↗

Mechanisms of voids formation during cooldown and freezing of lithium in SP-100 type systems

The mechanisms of void formation during the cooldown and freezing of lithium coolant within the primary loop of SP-100 type systems are investigated. These mechanisms are: (a) homogeneous nucleation, (b) heterogeneous nucleation, (c) normal segregation of helium gas dissolved in liquid lithium, and (d) shrinkage of lithium during freezing. To evaluate the void formation potential due to segregation, a numerical scheme that couples the freezing and mass diffusion processes in both the solid and liquid regions is developed. The results indicated that the formation of He bubbles is unlikely by either homogeneous or heterogeneous nucleation during the cooldown process. However, homogeneous nucleation of He bubbles following the segregation of dissolved He in liquid Lithium ahead of the solid-liquid interface is likely to occur. Results also show that total volume of He void is insignificant when compared to that of shrinkage voids.

Yang, Jae Y.↗

Local dynamic subgrid-scale models in channel flow

The dynamic subgrid-scale (SGS) model has given good results in the large-eddy simulation (LES) of homogeneous isotropic or shear flow, and in the LES of channel flow, using averaging in two or three homogeneous directions (the DA model). In order to simulate flows in general, complex geometries (with few or no homogeneous directions), the dynamic SGS model needs to be applied at a local level in a numerically stable way. Channel flow, which is inhomogeneous and wall-bounded flow in only one direction, provides a good initial test for local SGS models. Tests of the dynamic localization model were performed previously in channel flow using a pseudospectral code and good results were obtained. Numerical instability due to persistently negative eddy viscosity was avoided by either constraining the eddy viscosity to be positive or by limiting the time that eddy viscosities could remain negative by co-evolving the SGS kinetic energy (the DLk model). The DLk model, however, was too expensive to run in the pseudospectral code due to a large near-wall term in the auxiliary SGS kinetic energy (k) equation. One objective was then to implement the DLk model in a second-order central finite difference channel code, in which the auxiliary k equation could be integrated implicitly in time at great reduction in cost, and to assess its performance in comparison with the plane-averaged dynamic model or with no model at all, and with direct numerical simulation (DNS) and/or experimental data. Other local dynamic SGS models have been proposed recently, e.g., constrained dynamic models with random backscatter, and with eddy viscosity terms that are averaged in time over material path lines rather than in space. Another objective was to incorporate and test these models in channel flow.

Cabot, William H.↗

Induced nucleation of carbon dust in red giant stars

This study quantitatively tests the proposed model of induced nucleation of carbonaceous grains in carbon-rich red giant stars. Induced nucleation is the process of grain growth initiated by the presence of reactive surfaces provided by seed particles. The numerical study was performed using a deailed chemical kinetic model of carbon deposition, grain coagulation, and homogeneous nucleation of polycyclic aromatic hydrocarbons (PAHs). The model uses a method of moments to keep track of developing grain population in the forming dust shell. We test the efficiency of grain formation for large ranges of dust shell parameters typical for carbon stars. Our model is capable of producing a range of optically thick and thin dust shells in carbon stars. Results are in accord with (IRAS) spectral classes of carbon stars. The resulting composite grains produced are consistent with those recently found in ancient meteorites. This model also provides a realistic explanation for high abundances of (PAHs) in the interstellar medium and some planetary nebulae.

Cadwell, Brian J.↗

Electromagnetic Field Penetration Studies

A numerical method is presented to determine electromagnetic shielding effectiveness of rectangular enclosure with apertures on its wall used for input and output connections, control panels, visual-access windows, ventilation panels, etc. Expressing EM fields in terms of cavity Green's function inside the enclosure and the free space Green's function outside the enclosure, integral equations with aperture tangential electric fields as unknown variables are obtained by enforcing the continuity of tangential electric and magnetic fields across the apertures. Using the Method of Moments, the integral equations are solved for unknown aperture fields. From these aperture fields, the EM field inside a rectangular enclosure due to external electromagnetic sources are determined. Numerical results on electric field shielding of a rectangular cavity with a thin rectangular slot obtained using the present method are compared with the results obtained using simple transmission line technique for code validation. The present technique is applied to determine field penetration inside a Boeing-757 by approximating its passenger cabin as a rectangular cavity filled with a homogeneous medium and its passenger windows by rectangular apertures. Preliminary results for, two windows, one on each side of fuselage were considered. Numerical results for Boeing-757 at frequencies 26 MHz, 171-175 MHz, and 428-432 MHz are presented.

Deshpande, M.D.↗

Improved methods for large-eddy simulations of turbulence

By using Fourier transforms for evaluating spatial derivatives, the accuracy of the large-eddy simulation of homogeneous isotropic turbulence is improved. Especially, the treatment of certain terms that arise in filtering the equations is considerably improved in both speed and accuracy. Use of vorticity as the principal variable is shown to be a viable and potentially useful alternative to the primitive variables. A method of deriving conservation properties of numerical schemes is given which is much simpler than previous methods and is widely applicable. The methods are applied to the computation of homogeneous isotropic turbulence, and it is found that the subgrid scale model is improved by using finite differences in place of 'exact' derivatives.

Mansour, N. N.↗

Turbulence structure at high shear rate

The structure of homogeneous turbulence in the presence of a high shear rate is studied using results obtained from three-dimensional time-dependent numerical simulations of the Navier-Stokes equations on a grid of 512 x 128 x 128 node points. It is shown that high shear rate enhances the streamwise fluctuating motion to such an extent that a highly anisotropic turbulence state with a one-dimensional velocity field and two-dimensional small-scale turbulence develops asymptotically as total shear increases. Instantaneous velocity fields show that high shear rate in homogeneous turbulent shear flow produces structures which are similar to the streaks present in the viscous sublayer of turbulent boundary layers.

Lee, Moon Joo↗