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 163 records · Page 9

In General, the Total Voltammetric Current from a Mixture of Redox-Active Substances will Not be the Sum of the Currents that Each Substance would Produce Independently at the Same Concentration as in the Mixture

At the potential range where both decamethylferrocene (dMeFc) and ferrocene (Fc) are oxidized with rates controlled by linear diffusion, electrogenerated Fc(+) radicals diffusing outwards from the electrode react quantitatively (K23 C=5.8 x 10(exp 8) with dMeFc diffusing towards the electrode and produce Fc and dMeFc. That reaction replaces dMeFc with Fc, whose diffusion coefficient is higher than that of dMeFc(+), and the total mass-transfer limited current from the mixture is increased by approximately 10%. Analogous observations are made when mass-transfer is controlled by convective-diffusion as in RDE voltammetry. Similar results have been obtained with another, and for all practical purposes randomly selected pair of redox-active substances, [Co(bipy)3](2+) and N - methylphenothiazine (MePTZ); reaction of MePTZ(+) with [Co(bipy)3](2+) replaces the latter with MePTZ, which diffuses faster and the current increases by approximately 20%. The experimental voltammograms have been simulated numerically and the role of (a) the rate constant of the homogeneous reaction; (b) the relative concentrations; and, (c) the diffusion coefficients of all species involved have been studied in detail. Importantly, it was also identified that within any given redox system the dependence of the mass-transfer limited current on the bulk concentrations of the redox-active species is expected to be non-linear. These findings are discussed in terms of their electroanalytical implications.

Leventis, Nicholas↗

Hydrodynamic Modeling of Oxidizer-Rich Staged Combustion Injector Flow

The main objective of this work is to determine the unsteady hydrodynamic characteristics of coaxial swirl atomizers of interest in oxidizer-rich staged combustion (ORSC) liquid rocket engines. To this end, the pseudo-density (homogeneous flow) treatment combined with the Marker-and-Cell (MAC) numerical algorithm has been used to develop an axisymmetric with swirl, two-phase, unsteady model. The numerical model is capable of assessing the time-dependent orifice exit conditions and internal mixing for arbitrary fuel and oxidizer gas injection conditions. Parametric studies have been conducted to determine the effect of geometry, gas properties, and liquid properties on the exit massflow rate and velocity. It has been found that the frequency at which the liquid film oscillates increases as the density ratio and thickness increase, decreases as film thickness and liquid swirl velocity increase, and is unaffected by the mixing length. Additionally, it has been determined that the variation in the massflow rate increases as the liquid swirl velocity and liquid film thickness increase, and decreases as the density ratio, collar thickness, and mixing length increase.

Canino, J. V.↗

Applicability of the Effective-Medium Approximation to Heterogeneous Aerosol Particles.

The effective-medium approximation (EMA) is based on the assumption that a heterogeneous particle can have a homogeneous counterpart possessing similar scattering and absorption properties. We analyze the numerical accuracy of the EMA by comparing superposition T-matrix computations for spherical aerosol particles filled with numerous randomly distributed small inclusions and Lorenz-Mie computations based on the Maxwell-Garnett mixing rule. We verify numerically that the EMA can indeed be realized for inclusion size parameters smaller than a threshold value. The threshold size parameter depends on the refractive-index contrast between the host and inclusion materials and quite often does not exceed several tenths, especially in calculations of the scattering matrix and the absorption cross section. As the inclusion size parameter approaches the threshold value, the scattering-matrix errors of the EMA start to grow with increasing the host size parameter and or the number of inclusions. We confirm, in particular, the existence of the effective-medium regime in the important case of dust aerosols with hematite or air-bubble inclusions, but then the large refractive-index contrast necessitates inclusion size parameters of the order of a few tenths. Irrespective of the highly restricted conditions of applicability of the EMA, our results provide further evidence that the effective-medium regime must be a direct corollary of the macroscopic Maxwell equations under specific assumptions.

aerosols↗

Annual Research Briefs - 1996, Center for Turbulence Research

Partial contents include: (1) Simulation of jet impingement heat transfer with the k-e-V2 model; (2) A new methodology for turbulence modelers using Direct Numerical Simulation (DNS) database analysis; (3) A particle representation model for the deformation of homogeneous turbulence; (4) Modeling of inhomogeneous compressible turbulence using a two-scale statistical theory; (5) Lewis number effects on partially premixed flames; (6) Conditional moment closure of mixing and reaction in turbulent non-premixed combustion; (7) Dynamic models for LES of turbulent front propagation with a spectral method; (8) Jet noise using large eddy simulation; (9) Towards numerical simulations of trailing-edge aeroacoustics; (10) A model for the onset of vortex breakdown; (11) Compressing an elliptic vortex: transition to turbulence by tumble breakdown; (12) A new method for the adaptive control of vortex-wall interactions; (13) Symmetries in turbulent boundary layer flows; (14) Near-wall models in large eddy simulations of flow behind a backward-facing step; (15) Assessment of the subgrid-scale models at low and high Reynolds numbers; (16) Large-eddy simulations of flow around a NACA 4412 airfoil using unstructured grids; (17) Progress on LES of flow past a circular cylinder; (18) Experimental investigation of flow through an asymmetric plane diffuser; (19) Progress in the large-eddy simulation of an asymmetric plane diffuser; (20) A numerical study of self-similarity in a turbulent plane wake using large-eddy simulation; (21) On the generation of turbulent inflow conditions for boundary layer simulations, and (22) A new class of finite difference schemes.

COMPUTATIONAL FLUID DYNAMICS↗

Spurious Behavior of Shock-Capturing Methods: Problems Containing Stiff Source Terms and Discontinuities

The goal of this paper is to relate numerical dissipations that are inherited in high order shock-capturing schemes with the onset of wrong propagation speed of discontinuities. For pointwise evaluation of the source term, previous studies indicated that the phenomenon of wrong propagation speed of discontinuities is connected with the smearing of the discontinuity caused by the discretization of the advection term. The smearing introduces a nonequilibrium state into the calculation. Thus as soon as a nonequilibrium value is introduced in this manner, the source term turns on and immediately restores equilibrium, while at the same time shifting the discontinuity to a cell boundary. The present study is to show that the degree of wrong propagation speed of discontinuities is highly dependent on the accuracy of the numerical method. The manner in which the smearing of discontinuities is contained by the numerical method and the overall amount of numerical dissipation being employed play major roles. Moreover, employing finite time steps and grid spacings that are below the standard Courant-Friedrich-Levy (CFL) limit on shockcapturing methods for compressible Euler and Navier-Stokes equations containing stiff reacting source terms and discontinuities reveals surprising counter-intuitive results. Unlike non-reacting flows, for stiff reactions with discontinuities, employing a time step and grid spacing that are below the CFL limit (based on the homogeneous part or non-reacting part of the governing equations) does not guarantee a correct solution of the chosen governing equations. Instead, depending on the numerical method, time step and grid spacing, the numerical simulation may lead to (a) the correct solution (within the truncation error of the scheme), (b) a divergent solution, (c) a wrong propagation speed of discontinuities solution or (d) other spurious solutions that are solutions of the discretized counterparts but are not solutions of the governing equations. The present investigation for three very different stiff system cases confirms some of the findings of Lafon & Yee (1996) and LeVeque & Yee (1990) for a model scalar PDE. The findings might shed some light on the reported difficulties in numerical combustion and problems with stiff nonlinear (homogeneous) source terms and discontinuities in general.

Yee, Helen M. C.↗

Analysis of Plasma-Sprayed Thermal Barrier Coatings With Homogeneous and Heterogeneous Bond Coats Under Spatially Uniform Cyclic Thermal Loading

This report summarizes the results of a numerical investigation into the spallation mechanism in plasma-sprayed thermal barrier coatings observed under spatially-uniform cyclic thermal loading. The analysis focuses on the evolution of local stress and inelastic strain fields in the vicinity of the rough top/bond coat interface during thermal cycling, and how these fields are influenced by the presence of an oxide film and spatially uniform and graded distributions of alumina particles in the metallic bond coat aimed at reducing the top/bond coat thermal expansion mismatch. The impact of these factors on the potential growth of a local horizontal delamination at the rough interface's crest is included. The analysis is conducted using the Higher-Order Theory for Functionally Graded Materials with creep/relaxation constituent modeling capabilities. For two-phase bond coat microstructures, both the actual and homogenized properties are employed in the analysis. The results reveal the important contributions of both the normal and shear stress components to the delamination growth potential in the presence of an oxide film, and suggest mixed-mode crack propagation. The use of bond coats with uniform or graded microstructures is shown to increase the potential for delamination growth by increasing the magnitude of the crack-tip shear stress component.

Arnold, Steven M.↗

On the Susceptibility of Cold Tropical Cirrus to Ice Nuclei Abundance

Numerical simulations of cirrus formation in the tropical tropopause layer (TTL) are used to evaluate the impact heterogeneous ice nuclei (IN) abundance on cold cloud microphysical properties. The model includes homogeneous and heterogeneous ice nucleation, deposition growthsublimation, and sedimentation. Reanalysis temperature and wind fields are used to force the simulations, with addition of high-frequency wave temperature variability based on in situ measurements. The model results are constrained by recent in situ observations of TTL cirrus and relative humidity, as well as satellite measurements of cirrus bulk properties and occurrence frequencies. Temperature variability driven by high-frequency waves has a dominant influence on TTL cirrus microphysical properties and occurrence frequencies, and inclusion of these waves is required to produce agreement between the simulated and observed cloud properties as well as the abundance of TTL cirrus. With a composition-independent supersaturation threshold for homogeneous freezing of aqueous aerosols, the model produces excessive ice concentrations compared with in situ observations. Inclusion of relatively numerous heterogeneous nuclei ( 100L) in the simulations improves the agreement with observed ice concentrations. However, when IN contribute significantly to TTL cirrus ice nucleation, the occurrence frequency of large supersaturations with respect to ice is less than indicated by in situ measurements. The simulated TTL cirrus extinction statistics agree with observations (within uncertainties) from both in situ and remote-sensing measurements. We find that the sensitivity of TTL cirrus extinction and ice water content statistics to heterogeneous ice nuclei abundance is relatively weak. Likewise, the simulated occurrence frequencies of TTL cirrus are insensitive to ice nuclei abundance, both in terms of cloud frequency height distribution and regional distribution throughout the tropics.

occurences↗

The influence of line shape and band structure on temperatures in planetary atmospheres

Numerical experiments are performed to examine the effects of line shape and band structure on the radiative equilibrium temperature profile in planetary atmospheres. In order to accurately determine these effects, a method for calculating radiative terms is developed which avoids the usual approximations. It differs from the more commonly used methods in that it allows arbitrary dependence of the absorption coefficient on wave number, without requiring tedious line by line integration and without the constraints of band models. The present formulation is restricted to homogeneous atmospheres but the concept can be extended to the more general case. The numerical experiments reveal that the line shape and band structure of the absorbing gases have a large effect on temperatures in the higher layers of the atmosphere (corresponding to the stratosphere and mesosphere). The more nongrey the spectrum (that is, the higher the peaks and the deeper the troughs in the spectrum), the lower the temperature.

Arking, A.↗

Three-dimensional time dependent computation of turbulent flow

The three-dimensional, primitive equations of motion are solved numerically for the case of isotropic box turbulence and the distortion of homogeneous turbulence by irrotational plane strain at large Reynolds numbers. A Gaussian filter is applied to governing equations to define the large scale field. This gives rise to additional second order computed scale stresses (Leonard stresses). The residual stresses are simulated through an eddy viscosity. Uniform grids are used, with a fourth order differencing scheme in space and a second order Adams-Bashforth predictor for explicit time stepping. The results are compared to the experiments and statistical information extracted from the computer generated data.

Kwak, D.↗

Magnetic dynamo activity in mechanically driven compressible magnetohydrodynamic turbulence

Magnetic dynamo activity in a homogeneous, dissipative, polytropic, two-dimensional, turbulent magneto-fluid is simulated numerically. The magneto-fluid is simulated numerically. The magneto-fluid is, in a number of cases, mechanically forced so that energy input balances dissipation, thereby maintaining constant energy. In the presence of a mean magnetic field, a magneto-fluid whose initial turbulent magnetic energy is zero quickly arrives at a state of non-zero turbulent magnetic energy. If the mean magnetic field energy density is small, the turbulent magnetic field can achieve a local energy density more than four hundred times larger; if the mean magnetic field energy density is large, then equipartition between the turbulent magnetic and kinetic energy is achieved. Compared to the presence of a mean magnetic field, compressibility appears to have only a marginal effect in mediating the transfer of turbulent kinetic energy into magnetic energy.

Shebalin, John V.↗

Inelastic transient dynamic analysis of three-dimensional problems by BEM

A general direct boundary element formulation and its numerical implementation for solving transient dynamic problems of three-dimensional isotropic homogeneous or piecewise homogeneous solids involving material nonlinearities are presented. The algorithm produces accurate results for static nonlinear problems by using large time steps. When a large value of yield stress is selected, the incremental inelastic transient algorithm produces results identical to those obtained by elastodynamic analysis.

Ahmad, S.↗

Flight experiment on acousto-optic crystal growth

The physical vapor transport method was used for growing mercurous chloride crystals in different convective conditions. Optical homogeneity is found to be extremely dependent on convection levels. Results of numerical studies indicate that for a gravity level of 0.001 g or less the Stefan wind drives the flow and no recirculating cells are observed.

Singh, N. B.↗

Simulation and modeling of the elliptic streamline flow

Direct numerical simulations are performed for the elliptic streamline flow, which is a homogeneous turbulent flow that combines the effects of solid body rotation and strain. Simulations are run over a range of parameters in order to determine the effect of changing rotation and strain separately. For early times the nonlinear cascade is suppressed, but then is re-established at later times. The growth rate of turbulent kinetic energy agrees at early times with the trends from linear theory, but at later times the flow seems to approach an asymptotic state that is independent of the ratio of mean flow rotation rate to strain rate. A comparison with standard Reynolds stress turbulence models is made. It is found that for strong rotation rates, the models predict decay of the turbulence, while the simulations show exponential growth. Close examination of the simulation results shows that they are affected by excessively low Reynolds numbers. Suggestions for reducing low Reynolds number effects in future simulations is given.

Blaisdell, G. A.↗

Applying Orbital Multi-Angle Photopolarimetric Observations to Study Properties of Aerosols in the Earth's Atmosphere: Implications of Measurements in the 1.378 µm Spectral Channel to Retrieve Microphysical Characteristics and Composition of Stratospheric Aerosols

We analyze the possibilities of orbital photopolarimetric measurements to study properties of aerosols in the Earth's atmosphere. As an example, we consider the case when such measurements are performed within a narrow spectral channel centered at 1.378 µm that allows to retrieve microphysical characteristics of stratospheric aerosols separately from those of tropospheric aerosols. We consider the case of stratospheric aerosols caused by volcanic eruption, and adopt the model of the stratosphere in the form of a homogeneous plane-parallel layer composed of polydisperse spherical particles. We use numerically exact solutions of the vector radiative transfer equation to theoretically simulate measurements carried out at various numbers of scattering angles, including: (i) radiance measurements alone; (ii) polarization measurements alone; and (iii) radiance and polarization measurements together. The results of computations show that the simultaneous use of radiance and polarization measurements at a sufficiently large number of scattering angles enables one to retrieve the optical thickness, effective radius, and refractive index of aerosols with adequate accuracy. We demonstrate how the accuracy of the derived values of the optical parameters of aerosols depends on the accuracy of measurements of the intensity and polarization of the reflected light, optical thickness of aerosol layer itself, effective radius of aerosols, width of the particle size distribution, and number of viewing angles.

Aerosols↗

A New Approach for a Wider Class of Entropy Split Methods for Compressible Gas Dynamics and MHD

The high order entropy split methods of Sjögreen & Yee [1, 2] by entropy splitting of the compressible Euler (inviscid) flux derivatives for a thermally-perfect gas are based on Harten’s entropy function [3, 4, 5]. Their derivation takes advantage of the homogeneity property of Euler flux, symmetrizable Euler flux derivatives and energy-norm stability in conjunction with high order classical spatial central, DRP (dispersion relation-preserving) [6, 7, 8] or Padé (compact) spatial discretizations [9] with summation-by-parts (SBP) operators [10]. Our entropy split methods have been proven entropy conserving and stable [1, 11, 12]. Our proofs do not rely on a two-point numerical flux, but rather only a linear difference operator is required to derive these methods. To extend the entropy split method for the MHD, we used the Godunov symmetrizable non-conservative MHD form [12, 13, 14]. These high order entropy split methods not only preserve certain physical properties of the chosen governing equations but are also known to either improve numerical stability, and/or minimize aliasing errors in long time integration of turbulent flow computations without the aid of added numerical dissipation. In our previous published work, extensive error norm comparison with grid refinement was performed to show the high accuracy performance of these methods. These studies also showed how well the entropy split methods conserve the entropy, momentum and mass, and preserve the kinetic energy for long time integration of the various flows [1, 2, 12, 13, 14]. The objective of the present work is to use a new approach to obtain a wider class of entropy split methods consisting of a two-point numerical flux portion and a non-conservative portion in such a way that the homogeneity property of the compressible Euler flux is not required. For high order classical spatial central, DRP (dispersion relation-preserving) or Padé (compact) spatial discretizations, this new approach can be proven to be entropy conservative with conservative spatial dsicretizations while at the same time allowing a wider class of symmetrizable inviscid flux derivatives. We also use this generalization to derive an entropy split scheme that is entropy conserving for the equations of MHD without the homogeneity property using the Godunov symmetrizable ideal MHD formulation [15].

High Order Physical Preserving Methods↗

Direct numerical simulations of turbulent shear flows

Numerical simulations of wakes of axisymmetric bodies and of turbulent mixing layers are reported. The flows were assumed to be statistically homogeneous in the mean flow direction, in concert with experimental data and the self-similarity theorem. The nonlinear Navier-Stokes equations were solved by a pseudo-spectral numerical method using a 32 x 32 x 33 point grid and an algorithm for fast Fourier transforms and inverse transforms. Leapfrog time differencing was employed on nonlinear terms and time differencing on viscous terms. Towed wakes and wakes behind a self-propelled body were simulated, showing that the towed wakes exhibited a proper temporal behavior after an initial period of adjustment, including the development of a kurtosis near the wake edge, which is experimentally verifiable. The mixing-layer simulation displayed the laboratory demonstrated presence of large scale features such as vortex cores, while the lateral coherence was weak.

Metcalfe, R. W.↗

A numerical method for interface problems in elastodynamics

The numerical implementation of a formulation for a class of interface problems in elastodynamics is discussed. This formulation combines the use of the finite element and boundary integral methods to represent the interior and the exteriro regions, respectively. In particular, the response of a semicylindrical alluvial valley in a homogeneous halfspace to incident antiplane SH waves is considered to determine the accuracy and convergence of the numerical procedure. Numerical results are obtained from several combinations of the incidence angle, frequency of excitation, and relative stiffness between the inclusion and the surrounding halfspace. The results tend to confirm the theoretical estimates that the convergence is of the order H(2) for the piecewise linear elements used. It was also observed that the accuracy descreases as the frequency of excitation increases or as the relative stiffness of the inclusion decreases.

Mcghee, D. S.↗

Statistics of pressure and pressure gradient in homogeneous isotropic turbulence

The statistics of pressure and pressure gradient in stationary isotropic turbulence are measured within direct numerical simulations at low to moderate Reynolds numbers. It is found that the one-point pdf of the pressure is highly skewed and that the pdf of the pressure gradient is of stretched exponential form. The power spectrum of the pressure P(k) is found to be larger than the corresponding spectrum P(sub G)(k) computed from a Gaussian velocity field having the same energy spectrum as that of the DNS field. The ratio P(k)/P(sub G)(k), a measure of the pressure-field intermittence, grows with wavenumber and Reynolds number as -R(sub lambda)(exp 1/2)log(k/k(sub d)) for k less than k(sub d)/2 where k(sub d) is the Kolmogorov wavenumber. The Lagrangian correlations of pressure gradient and velocity are compared and the Lagrangian time scale of the pressure gradient is observed to be much shorter than that of the velocity.

Gotoh, T.↗