Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Homogenization theory”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Aerosol Optical Properties Measured Onboard the Ronald H. Brown During ACE Asia as a Function of Aerosol Chemical Composition and Source Region

During the ACE Asia intensive field campaign conducted in the spring of 2001 aerosol properties were measured onboard the R/V Ronald H. Brown to study the effects of the Asian aerosol on atmospheric chemistry and climate in downwind regions. Aerosol properties measured in the marine boundary layer included chemical composition; number size distribution; and light scattering, hemispheric backscattering, and absorption coefficients. In addition, optical depth and vertical profiles of aerosol 180 deg backscatter were measured. Aerosol within the ACE Asia study region was found to be a complex mixture resulting from marine, pollution, volcanic, and dust sources. Presented here as a function of air mass source region are the mass fractions of the dominant aerosol chemical components, the fraction of the scattering measured at the surface due to each component, mass scattering efficiencies of the individual components, aerosol scattering and absorption coefficients, single scattering albedo, Angstrom exponents, optical depth, and vertical profiles of aerosol extinction. All results except aerosol optical depth and the vertical profiles of aerosol extinction are reported at a relative humidity of 55 +/- 5%. An over-determined data set was collected so that measured and calculated aerosol properties could be compared, internal consistency in the data set could be assessed, and sources of uncertainty could be identified. By taking into account non-sphericity of the dust aerosol, calculated and measured aerosol mass and scattering coefficients agreed within overall experimental uncertainties. Differences between measured and calculated aerosol absorption coefficients were not within reasonable uncertainty limits, however, and may indicate the inability of Mie theory and the assumption of internally mixed homogeneous spheres to predict absorption by the ACE Asia aerosol. Mass scattering efficiencies of non-sea salt sulfate aerosol, sea salt, submicron particulate organic matter, and dust found for the ACE Asia aerosol are comparable to values estimated for ACE 1, Aerosols99, and INDOEX. Unique to the ACE Asia aerosol was the large mass fractions of dust, the dominance of dust in controlling the aerosol optical properties, and the interaction of dust with soot aerosol.

Quinn, P. K.↗
NASA NTRS2004

Theory and Simulation of Real and Ideal Magnetohydrodynamic Turbulence

Incompressible, homogeneous magnetohydrodynamic (MHD) turbulence consists of fluctuating vorticity and magnetic fields, which are represented in terms of their Fourier coefficients. Here, a set of five Fourier spectral transform method numerical simulations of two-dimensional (2-D) MHD turbulence on a 512(sup 2) grid is described. Each simulation is a numerically realized dynamical system consisting of Fourier modes associated with wave vectors k, with integer components, such that k = |k| less than or equal to k(sub max). The simulation set consists of one ideal (non-dissipative) case and four real (dissipative) cases. All five runs had equivalent initial conditions. The dimensions of the dynamical systems associated with these cases are the numbers of independent real and imaginary parts of the Fourier modes. The ideal simulation has a dimension of 366104, while each real simulation has a dimension of 411712. The real runs vary in magnetic Prandtl number P(sub M), with P(sub M) is a member of {0.1, 0.25, 1, 4}. In the results presented here, all runs have been taken to a simulation time of t = 25. Although ideal and real Fourier spectra are quite different at high k, they are similar at low values of k. Their low k behavior indicates the existence of broken symmetry and coherent structure in real MHD turbulence, similar to what exists in ideal MHD turbulence. The value of PM strongly affects the ratio of kinetic to magnetic energy and energy dissipation (which is mostly ohmic). The relevance of these results to 3-D Navier-Stokes and MHD turbulence is discussed.

Shebalin, John V.↗
NASA NTRS2012

Magnetic Helicity and Planetary Dynamos

A model planetary dynamo based on the Boussinesq approximation along with homogeneous boundary conditions is considered. A statistical theory describing a large-scale MHD dynamo is found, in which magnetic helicity is the critical parameter

Shebalin, John V.↗
NASA NTRS2012

Characterization and Analysis of Triaxially Braided Polymer Composites under Static and Impact Loads

In order to design impact resistant aerospace components made of triaxially-braided polymer matrix composite materials, a need exists to have reliable impact simulation methods and a detailed understanding of the material behavior. Traditional test methods and specimen designs have yielded unrealistic material property data due to material mechanisms such as edge damage. To overcome these deficiencies, various alternative testing geometries such as notched flat coupons have been examined to alleviate difficulties observed with standard test methods. The results from the coupon level tests have been used to characterize and validate a macro level finite element-based model which can be used to simulate the mechanical and impact response of the braided composites. In the analytical model, the triaxial braid unit cell is approximated by using four parallel laminated composites, each with a different fiber layup, which roughly simulates the braid architecture. In the analysis, each of these laminated composites is modeled as a shell element. Currently, each shell element is considered to be a smeared homogeneous material. Simplified micromechanics techniques and lamination theory are used to determine the equivalent stiffness properties of each shell element, and results from the coupon level tests on the braided composite are used to back out the strength properties of each shell element. Recent improvements to the model include the incorporation of strain rate effects into the model. Simulations of ballistic impact tests have been carried out to investigate and verify the analysis approach.

Goldberg, Robert K.↗
NASA NTRS1989

Composite fracture using the normal stress ratio theory

Experiments were performed on center-notched specimens of unidirectional AS4/3501-6 graphite-epoxy subjected to varied far-field loading conditions. Specimens tested include on-axis and off-axis tensile coupons and Iosipescu shear specimens, allowing the study of normal, shear, and combined far-field loading conditions. From the tests, the direction of initial crack extension and the far-field stresses at crack initiation and fracture are obtained. The normal stress ratio theory is subsequently applied within a macroscopic-level homogeneous anisotropic stress solution to predict the direction of crack extension, critical stresses, and other crack growth behavior for the experiments. Comparison between experimental results and theoretical predictions indicates that application of the normal stress ratio theory within a macroscopic-level stress analysis has significant potential as an approach to modeling composite crack growth.

Beuth, J. L., Jr.↗
DOE OSTI2020

Polycrystal thermo-elasticity revisited: theory and applications

The self-consistent (SC) theory is the most commonly used mean-field homogenization method to estimate the mechanical response behavior of polycrystals based on the knowledge of the properties and orientation distribution of constituent single-crystal grains. The original elastic SC method can be extended to thermo-elasticity by adding a stress-free strain to an elastic constitutive relation that expresses stress as a linear function of strain. With the addition of this independent term, the problem remains linear. Although the thermo-elastic self-consistent (TESC) model has important theoretical implications for the development of self-consistent homogenization of non-linear polycrystals, in this paper, we focus on TESC applications to actual thermo-elastic problems involving non-cubic (i.e. thermally anisotropic) materials. To achieve this aim, here we provide a thorough description of the TESC theory, which is followed by illustrative examples involving cooling of polycrystalline non-cubic metals. The TESC model allows studying the effect of crystallographic texture and single-crystal elastic and thermal anisotropy on the effective thermo-elastic response of the aggregate and on the internal stresses that develop at the local level.

36 MATERIALS SCIENCE↗
NASA NTRS2008

Understanding Ice Supersaturation, Particle Growth, and Number Concentration in Cirrus Clouds

Many factors control the ice supersaturation and microphysical properties in cirrus clouds. We explore the effects of dynamic forcing, ice nucleation mechanisms, and ice crystal growth rate on the evolution and distribution of water vapor and cloud properties in nighttime cirrus clouds using a one-dimensional cloud model with bin microphysics and remote sensing measurements obtained at the Department of Energy's Atmospheric Radiation Measurement (ARM) Climate Research Facility located near Lamont, OK. We forced the model using both large-scale vertical ascent and, for the first time, mean mesoscale velocity derived from radar Doppler velocity measurements. Both heterogeneous and homogeneous nucleation processes are explored, where a classical theory heterogeneous scheme is compared with empirical representations. We evaluated model simulations by examining both bulk cloud properties and distributions of measured radar reflectivity, lidar extinction, and water vapor profiles, as well as retrieved cloud microphysical properties. Our results suggest that mesoscale variability is the primary mechanism needed to reproduce observed quantities. Model sensitivity to the ice growth rate is also investigated. The most realistic simulations as compared with observations are forced using mesoscale waves, include fast ice crystal growth, and initiate ice by either homogeneous or heterogeneous nucleation. Simulated ice crystal number concentrations (tens to hundreds particles per liter) are typically two orders of magnitude smaller than previously published results based on aircraft measurements in cirrus clouds, although higher concentrations are possible in isolated pockets within the nucleation zone.

Comstock, Jennifer M.↗
NASA NTRS1994

Mode properties of low-frequency waves: Kinetic theory versus Hall-MHD

In fluid theory, the ordering of low-frequency modes in a homogeneous plasma is based on the phase velocity, since modes do not intersect each other in dispersion diagrams as a function of wavenumber or other parameters. In linear kinetic theory, modes cross each other. Thus a consistent and useful classification should be based on the physical properties of the modes instead. This paper attempts such a classification by documeting the dispersion and general mode properties of the low-frequency waves (omega much less than (OMEGA(sub ci) OMEGA(sub ce) (exp 1/2)), where OMEGA(sub ci), OMEGA(sub ce) are the cyclotron frequencies of the ions and electrons, respectively) in kinetic theory, and by comparing them to the results of two-fluid theory. Kinetic theory gives a seperate Alfven/ion-cyclotron (A/IC) wave with phase speed Omega/k approximately = v(sub A) cos theta for omega much less than OMEGA(sub ci), where v(sub A) is the Alfven velocity and theta the angle of propagation between wave vector k and background magnetic field B(sub o). For a given wavenumber, the magnetosonic mode is a double-valued solution with a singular point in theta, beta parameter space, where beta is the ratio of thermal pressure to magnetic pressure. It is shown that a branch cut starting at the singular point theta approximately 30 deg, beta approximately 3 and leading to larger beta gives a practical and consitent seperation of this double-valued magnetosonic solution. Selection of this branch cut results in a moderately damped fast/magnetos onic and a heavily damped slow/sound wave. A comprehensive review of the polarization, compressibility and other mode properties is given and shown to be consistent with the selected branch cut. At small wavenumbers, the kinetic mode properties typically start to deviate significantly from their fluid counterparts at beta approximately 0.5. At larger beta, there is no longer a consistent correspondence between the fluid and kinetic modes. Kinetic theory also dictates the use of different mode properties to distinguish between them in observational data. For example, the phase between the density and magnetic field perturbation may become useless at high beta, whereas the direction of the magnetic field perturbations with respect to k and B(sub o) remains a useful characteristic. Two quantities based on this characteristic are suggested and are shown to be useful also to distinguish between the mirror mode and A/IC waves in a plasma with temperature anisotropy.

Krauss-Varban, D.↗
DOE OSTI2023

Static Self-Energy and Effective Mass of the Homogeneous Electron Gas from Quantum Monte Carlo Calculations

We discuss the methodology of quantum Monte Carlo calculations of the effective mass based on the static self-energy Σ(k,0). We then use variational Monte Carlo calculations of Σ(k,0) of the homogeneous electron gas at various densities to obtain results very close to perturbative G 0 W 0 calculations for values of the density parameter 1 ≤ r s ≤ 10. Here, the obtained values for the effective mass are close to diagrammatic Monte Carlo results and disagree with previous quantum Monte Carlo calculations based on a heuristic mapping of excitation energies to those of an ideal gas.

36 MATERIALS SCIENCE↗
NASA NTRS1987

Evaluation of a theory for pressure-strain rate

A theoretical expression for the slow part (the nonlinear fluctuation part) of the pressure-strain rate is compared with simulations of anisotropic homogeneous flows. The objective is to determine the quantitative accuracy of the theory and to test its prediction that the generalized Rotta coefficient, a non-dimensionalized ratio of slow term to the Reynolds stress anisotropy, varies with direction and can be negative. Comparisons are made between theoretical and simulation values of the slow term itself and of the generalized Rotta coefficients. The implications of the comparison for two-point closure theories and for Reynolds stress modeling are pointed out.

Weinstock, J.↗
DOE OSTI2026

Is there an exact magnetic moment for charged particle motion in a time-dependent, homogeneous magnetic field?

The non-perturbative guiding-centre model provides an exact alternative to full-orbit simulations of charged particle dynamics in situations where traditional guiding-centre theory may fail. We demonstrate that the charged particle motion in a homogeneous, time-varying magnetic field is a solvable example of the non-perturbative guiding-centre model. This entails showing that the exact magnetic moment of Qin and Davidson can be constructed to be asymptotic to the adiabatic invariant series of Kruskal. In contrast to the perturbative invariant, the exact invariant contains information about parametric resonances. These resonances destroy the conservation of the usual magnetic moment over very long times. This refutes some previous claims about the all-time invariance of the magnetic moment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Revisiting Multi-Material Composite Structures with Homogenized Composite Properties

Composite structures inherently develop residual stresses during their curing process. Driven predominately by mismatched thermal strains between differing materials or ply orientations, but also affected by curing process phenomena like polymer shrinkage, these residual stresses can lead to failure within composite structures. There are several methods varying in complexity that can be used to model the development of residual stresses, all of which are capable of capturing sufficient detail to understand the residual stress state at the ply level. However, explicitly modeling all plies of a layup in a composite structure can be prohibitively expensive based on the number of plies, structure size, and required element size. The computational cost can be reduced through the homogenization of the composite layup without losing much fidelity of the overall response of the structure. The homogenization process reduces the many plies of a laminate to a single lamina that reduces complexity and increases the mesh size where a single element can span multiple plies. This report focuses on verification and validation efforts for a homogenization process using a suite of finite element simulations rather than an analytic solution derived from classical laminate theory. Initial verification using representative element volumes indicated there was minimal error in the homogenization process; however, this compounded to a small, but acceptable error in strip and split ring experimental composite structures. The error does under predict the residual stress state in the strip and split ring and should be accounted for when simulating composite structures with homogenized properties.

36 MATERIALS SCIENCE↗
NASA NTRS1982

Atmospheric Correction for Remote Sensor Data

Mathematical model based on radiative transfer theory as applied to plane parallel, horizontally spatially homogeneous, isotropic atmosphere. Model includes effects of multiple scattering by gases and aerosols, and effects of absorption by ozone in near ultraviolet and visible parts of spectrum.

Turner, R. E.↗
NASA NTRS1990

Low-frequency scattering from cylindrical structures at oblique incidence

Classical Rayleigh scattering theory is extended to the case of a homogeneous dielectric cylinder of arbitrary cross section whose transverse dimensions are much smaller than the wavelength. By assuming that the surface fields can be approximated by those of the infinite cylinder, the far zone scattered field is expressed in terms of polarizability tensors, whose properties are discussed. Numerical results are presented for circular, semicircular, triangular, and square cylinders. The results are applicable to the remote sensing of twigs, stalks, and vegetation needles at centimeter and millimeter wavelengths.

Sarabandi, Kamal↗
NASA NTRS1994

Size distribution and scattering phase function of aerosol particles retrieved from sky brightness measurements

Ground-based measurements of the solar transmission and sky radiance in a horizontal plane through the Sun are taken in several geographical regions and aerosol types: dust in a desert transition zone in Israel, sulfate particles in Eastern and Western Europe, tropical aerosol in Brazil, and mixed continental/maritime aerosol in California. Stratospheric aerosol was introduced after the eruption of Mount Pinatubo in June 1991. Therefore measurements taken before the eruption are used to analyze the properties of tropospheric aerosol; measurements from 1992 are also used to detect the particle size and concentration of stratospheric aerosol. The measurements are used to retrieve the size distribution and the scattering phase function at large scattering angles of the undisturbed aerosol particles. The retrieved properties represent an average on the entire atmospheric column. A comparison between the retrieved phase function for a scattering angle of 120 deg, with phase function predicted from the retrieved size distribution, is used to test the assumption of particle homogeneity and sphericity in radiative transfer models (Mie theory). The effect was found to be small (20% +/- 15%). For the stratospheric aerosol (sulfates), as expected, the phase function was very well predicted using the Mie theory. A model with a power law distribution, based on the spectral dependence of the optical thickness, alpha, cannot estimate accurately the phase function (up to 50% error for lambda = 0.87 microns). Before the Pinatubo eruption the ratio between the volumes of sulfate and coarse particles was very well correlated with alpha. The Pinatubo stratospheric aerosol destroyed this correlation. The aerosol optical properties are compared with analysis of the size, shape, and composition of the individual particles by electron microscopy of in situ samples. The measured volume size distribution before the injection of stratospheric aerosol consistently show two modes, sulfate particles with r(sub m) less than 0.2 microns and coarse paritcles with r(sub m) greater than 0.7 microns. The 'window' in the tropospheric aerosol in this radius range was used to observe a stable stratospheric aerosol in 1992, with r(sub m) approximately 0.5 microns. A combination of such optical thickness and sky measurements can be used to assess the direct forcing and the climatic impact of aerosol. Systematic inversion for the key aerosol types (sulfates, smoke, dust, and maritime aerosol) of the size distribution and phase function can give the relationship between the aerosol physical and optical properties that can be used to compute the radiative forcing. This forcing can be validated in dedicated field experiments.

Kaufman, Y. J.↗
NASA NTRS2000

How do Polar Stratospheric Clouds Form?

SOLVE measurements have been compared with results from a microphysical model to understand the composition and formation of the polar stratospheric clouds (PSCs) observed during SOLVE. Evidence that the majority of the particles remain liquid throughout the winter will be presented. However, a small fraction of the particles do freeze, and the presence of these frozen particles can not be explained by current theories, in which the only freezing mechanism is homogeneous freezing to ice below the ice frost point. Alternative formation mechanisms, in particular homogeneous freezing above the ice frost point and heterogeneous freezing, have been explored using the microphysical model. Both nitric acid trihydrate (NAT) and nitric acid dihydrate (NAD) have been considered as possible compositions for the solid-phase nitric acid aerosols. Comparisons between the model results and the SOLVE measurements will be used to constrain the possible formation mechanisms. Other effects of these frozen particles will also be discussed, in particular denitrification.

Drdla, Katja↗
NASA NTRS2016

HST-STIS Spectra of Saturn's Rings and Implications for Their Reddening Agent

We obtained HST-STIS spectra of Saturn's main rings in May 2011, using the G230L (and G430L) gratings, with final averaged radial resolution of 160 (and 330) km/pixel. The dataset filled a previous 200-330nm "spectral gap" between Cassini and ground-based spectra. The data provide radial profiles as a function of wavelength, but our most basic product at this point is a set of very low-noise spectra, radially averaged over broad regions of the rings (A, B, C, and Cassini Division). The raw spectra required special processing to remove artifacts due to extended-source grating scatter. We have modeled the spectra using a new particle surface model, which corrects for on-surface shadowing due to the likely very rough ring particle surfaces, and avoids overestimation of intra-mixed "neutral absorber". We correct for non-classical layer effects and finite ring optical depth, and relate our observed reflectivities to the spherical albedos of individual smooth particles. We model these smooth particle albedos using standard Hapke theory for regolith grain mixtures that are either homogeneous and "intramixed" (nonicy absorbers dispersed in water ice regolith grains) or heterogeneous "intimate" mixtures. As candidates for the nonicy contaminants we have considered amorphous carbon, aromatic-rich and aliphatic-rich organic tholins, silicates, hematite and iron metal. For the A and B rings, we find that iron metal (including a new theoretical estimate of the refractive indices of nanometer-sized grains of iron) is not spectrally steep enough in the 200-300nm range, and that aliphatic-rich tholins are either too steep at short wavelengths or too flat at long wavelengths. However, less than 1% by mass of aromatic-rich tholins provides a very good fit across the entire spectral range with no gratuitous "neutral absorber" needed, and a minimum of additional free parameters. The best fits require forward-scattering regolith grains. For the C Ring and Cassini Division, additional absorbers are needed (updated results will be given).

Saturn's Rings↗