Results for “effective medium approximations” — Page 12

NASA NTRS1996

The Ice and Silicate Spectral Features for Dust Aggregates

The optical properties of inhomogeneous aggregates of dust particles are calculated. The Discrete Dipole Approximation (DDA) is applied to the calculation of light scattering by the dust aggregates. The mixtures of ices and silicates are considered. The IR profiles near ice and silicate spectral features (3 micron and 10 micron) are constructed. The influence of grain topology, chemical composition and porosity have been investigated. The comparison of exact results for inhomogeneous aggregate and the effective medium theory (the rules of Maxwell-Garnett and Bruggeman) is made.

Ilin, A. E.↗

NASA NTRSDate not supplied

Enhancing Radiative Transfer Models for Optimized Assimilation of Microwave and Radar Observations

Radiative transfer models are extensively used for the assimilation of satellite observations into NWP models as well as retrieving geophysical products from satellite measurements. The Community Radiative Transfer Model (CRTM) is a community model developed by NOAA JCSDA and widely used for different purposes requiring RT calculations. CRTM requires bulk optical properties of hydrometeors in the form of lookup tables in order to perform all-sky RT calculations. However, the current cloud scattering lookup tables in CRTM assume spherical shapes for all frozen hydrometeors, whereas actual clouds contain frozen particles with diverse shapes. The first part of this talk presents the implementation and validation of a comprehensive Discrete Dipole Approximation (DDA) cloud scattering database into CRTM, specifically targeting microwave frequencies. The DDA technique proves effective in simulating the optical properties of non-spherical hydrometeors in the microwave region. The original DDA database assumes total random orientation in calculating single scattering properties. To generate the required mass scattering parameters for CRTM, the single scattering properties and water content dependent particle size distributions are used. The evaluation of results involved a collocated dataset comprising short-term forecasts from the Integrated Forecast System of the European Center for Medium-Range Weather Forecasts and satellite microwave data. The findings demonstrate that the DDA lookup tables significantly reduce discrepancies between simulated and observed values when compared to the Mie tables. Passive instruments lack the ability to provide vertically resolved measurements of clouds and precipitation, which can be obtained by active radar instruments. However, incorporating these active measurements into data assimilation systems presents challenges due to the absence of fast forward radiative transfer models and difficulties in error modeling. The second part of the talk focuses on the development, evaluation, and sensitivity analysis of a forward radar model integrated into CRTM. The forward radar model utilizes scattering properties obtained from hydrometeor lookup tables generated using the discrete dipole approximation. By utilizing CRTM instrument-specific coefficients, the model can calculate both reflectivity and attenuated reflectivity for any given radar instrument and zenith angles. Evaluation using CloudSat measurements demonstrates a strong agreement between simulations and observations when the input profiles of hydrometeors align with the measured reflectivity profiles.

Isaac Moradi↗

NASA NTRSDate not supplied

Optimizing Assimilation of Microwave and Radar Observations in the NWP Models

Radiative transfer models are extensively used for the assimilation of satellite observations into NWP models as well as retrieving geophysical products from satellite measurements. The Community Radiative Transfer Model (CRTM) is a community model developed by NOAA JCSDA and widely used for different purposes requiring RT calculations. CRTM requires bulk optical properties of hydrometeors in the form of lookup tables in order to perform all-sky RT calculations. However, the current cloud scattering lookup tables in CRTM assume spherical shapes for all frozen hydrometeors, whereas actual clouds contain frozen particles with diverse shapes. The first part of this talk presents the implementation and validation of a comprehensive Discrete Dipole Approximation (DDA) cloud scattering database into CRTM, specifically targeting microwave frequencies. The DDA technique proves effective in simulating the optical properties of non-spherical hydrometeors in the microwave region. The original DDA database assumes total random orientation in calculating single scattering properties. To generate the required mass scattering parameters for CRTM, the single scattering properties and water content dependent particle size distributions are used. The evaluation of results involved a collocated dataset comprising short-term forecasts from the Integrated Forecast System of the European Center for Medium-Range Weather Forecasts and satellite microwave data. The findings demonstrate that the DDA lookup tables significantly reduce discrepancies between simulated and observed values when compared to the Mie tables. Passive instruments lack the ability to provide vertically resolved measurements of clouds and precipitation, which can be obtained by active radar instruments. However, incorporating these active measurements into data assimilation systems presents challenges due to the absence of fast forward radiative transfer models and difficulties in error modeling. The second part of the talk focuses on the development, evaluation, and sensitivity analysis of a forward radar model integrated into CRTM. The forward radar model utilizes scattering properties obtained from hydrometeor lookup tables generated using the discrete dipole approximation. By utilizing CRTM instrument-specific coefficients, the model can calculate both reflectivity and attenuated reflectivity for any given radar instrument and zenith angles. Evaluation using CloudSat measurements demonstrates a strong agreement between simulations and observations when the input profiles of hydrometeors align with the measured reflectivity profiles.

Isaac Moradi↗

NASA NTRSDate not supplied

Assimilation of Spaceborne Microwave and Radar Observations

Radiative transfer models are extensively used for the assimilation of satellite observations into NWP models as well as retrieving geophysical products from satellite measurements. The Community Radiative Transfer Model (CRTM) is a community model developed by NOAA JCSDA and widely used for different purposes requiring RT calculations. CRTM requires bulk optical properties of hydrometeors in the form of lookup tables in order to perform all-sky RT calculations. However, the current cloud scattering lookup tables in CRTM assume spherical shapes for all frozen hydrometeors, whereas actual clouds contain frozen particles with diverse shapes. The first part of this talk presents the implementation and validation of a comprehensive Discrete Dipole Approximation (DDA) cloud scattering database into CRTM, specifically targeting microwave frequencies. The DDA technique proves effective in simulating the optical properties of non-spherical hydrometeors in the microwave region. The original DDA database assumes total random orientation in calculating single scattering properties. To generate the required mass scattering parameters for CRTM, the single scattering properties and water content dependent particle size distributions are used. The evaluation of results involved a collocated dataset comprising short-term forecasts from the Integrated Forecast System of the European Center for Medium-Range Weather Forecasts and satellite microwave data. The findings demonstrate that the DDA lookup tables significantly reduce discrepancies between simulated and observed values when compared to the Mie tables. Passive instruments lack the ability to provide vertically resolved measurements of clouds and precipitation, which can be obtained by active radar instruments. However, incorporating these active measurements into data assimilation systems presents challenges due to the absence of fast forward radiative transfer models and difficulties in error modeling. The second part of the talk focuses on the development, evaluation, and sensitivity analysis of a forward radar model integrated into CRTM. The forward radar model utilizes scattering properties obtained from hydrometeor lookup tables generated using the discrete dipole approximation. By utilizing CRTM instrument-specific coefficients, the model can calculate both reflectivity and attenuated reflectivity for any given radar instrument and zenith angles. Evaluation using CloudSat measurements demonstrates a strong agreement between simulations and observations when the input profiles of hydrometeors align with the measured reflectivity profiles.

Isaac Moradi↗

NASA NTRS1988

Effective permittivity of dielectric mixtures

General mixing formulas are derived for discrete scatterers immersed in a host medium. The inclusion particles are assumed to be ellipsoidal. The electric field inside the scatterers is determined by quasi-static analysis, assuming the diameter of the inclusion particles to be much smaller than one wavelength. The results are applicable to general multiphase mixtures, and the scattering ellipsoids of the different phases can have different sizes and arbitrary ellipticity distribution and axis orientation, i.e., the mixture may be isotropic or anisotropic. The resulting mixing formula is nonlinear and is suitable for iterative solutions. The formula contains a quantity called the apparent permittivity, and with different choices of this quantity, the result leads to the generalized Lorentz-Lorenz formula, the generalized Polder-van Santen formula, and the generalized coherent potential-quasicrystalline approximation formula. The results are applied to calculating the complex effective permittivity of dry and wet snow, and sea ice.

Sihvola, Ari H.↗

NASA NTRS1976

The effect of adhesive layer on crack propagation in laminates

The effect of the adhesive layer on crack propagation in composite materials is investigated. The composite medium consists of parallel load carrying laminates and buffer strips arranged periodically and bonded with thin adhesive layers. The strips, assumed to be isotropic and linearly elastic, contain symmetric cracks of arbitrary lengths located normal to the interfaces. Two problems are considered: (1) thin adhesive layers are approximated by uncoupled tension and shear springs distributed along the interfaces of the strips for which only the case of internal cracks can be treated rigorously; (2) broken laminates and the true singular behavior in the presence of the adhesive layer are studied. The adhesive is then treated as an isotropic, linearly elastic continuum. General expressions for field quantities are obtained in terms of infinite Fourier integrals. These expressions give a system of singular integral equations in terms of the crack surface displacement derivatives. By using appropriate quadrature formulas, the integral equations reduce to a system of linear algebraic equations which are solved numerically.

Gecit, M. R.↗

NASA NTRS2013

Direct Simulation of Extinction in a Slab of Spherical Particles

The exact multiple sphere superposition method is used to calculate the coherent and incoherent contributions to the ensemble-averaged electric field amplitude and Poynting vector in systems of randomly positioned nonabsorbing spherical particles. The target systems consist of cylindrical volumes, with radius several times larger than length, containing spheres with positional configurations generated by a Monte Carlo sampling method. Spatially dependent values for coherent electric field amplitude, coherent energy flux, and diffuse energy flux, are calculated by averaging of exact local field and flux values over multiple configurations and over spatially independent directions for fixed target geometry, sphere properties, and sphere volume fraction. Our results reveal exponential attenuation of the coherent field and the coherent energy flux inside the particulate layer and thereby further corroborate the general methodology of the microphysical radiative transfer theory. An effective medium model based on plane wave transmission and reflection by a plane layer is used to model the dependence of the coherent electric field on particle packing density. The effective attenuation coefficient of the random medium, computed from the direct simulations, is found to agree closely with effective medium theories and with measurements. In addition, the simulation results reveal the presence of a counter-propagating component to the coherent field, which arises due to the internal reflection of the main coherent field component by the target boundary. The characteristics of the diffuse flux are compared to, and found to be consistent with, a model based on the diffusion approximation of the radiative transfer theory.

spheres↗

DOE OSTI2020

Medium-induced gluon radiation with full resummation of multiple scatterings for realistic parton-medium interactions

The new precision era of jet quenching observables at both RHIC and the LHC calls for an improved and more precise description of in-medium gluon emissions. The development of new theoretical tools and analytical calculations to tackle this challenge has been hampered by the inability to include the effects of multiple scatterings with the medium using a realistic model for the parton-medium interactions. In this paper, we show how the analytical expressions for the full in-medium spectrum, including the resummation of all multiple scatterings, can be written in a form where the numerical evaluation can be easily performed without the need of the usually employed harmonic or single hard approximations. We present the transverse momentum and energy-dependent medium- induced gluon emission distributions for known realistic interaction models to illustrate how our framework can be applied beyond the limited kinematic regions of previous calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1985

On the accuracy of the pseudocompressibility method in solving the incompressible Navier-Stokes equations

The method of pseudocompressibility is being tested for its accuracy in solving the incompressible Navier-Stokes equations. An implicit, finite-difference computer code is used to solve the equations in a three-dimensional, curvilinear coordinate system. The code employs artificial compressibility for solving the pressure field, coupled with an implicit, approximate-factorization scheme. This coupling is known as the pseudocompressibility method. The pseudocompressibility method introduces pressure waves of finite speed into the fluid medium that would otherwise have an infinite sound speed. The waves die out as the solution converges, and the steady state solution approaches a divergence-free condition. However, these waves limit the time-accuracy of the computations. The effects of these waves are analyzed and criteria are set for choosing the pseudocompressibility parameters that govern the pressure wave speed. Test cases are presented that verify these criteria. The code is tested by computing laminar flow over a two-dimensional, backward-facing step and over a two-dimensional, impulsively started circular cylinder.

Rogers, S. E.↗

NASA NTRS1993

A Fourier-Riccati approach to radiative transfer. I - Foundations

The three-dimensional equation of radiative transfer is formally solved using a Fourier-Riccati approach while calculations are performed on cloudy media embedded in a two-dimensional space. An extension to Stephens' work this study addresses the coupling between space and angle asserted by the equation of transfer. In particular, the accuracy of the computed radiation field as it is influenced by the angular resolution of the phase function and spatial discretization of the cloudy medium is discussed. The necessity of using a large number of quadrature points to calculate fluxes even when the phase function is isotropic for media exhibiting vertical and horizontal inhomogeneities is demonstrated. Effects of incorrect spatial sampling on both radiance and flux fields are also quantified by example. Radiance and flux comparisons obtained by the Fourier-Riccati model and the independent pixel approximation for inhomogeneous cloudy media illustrate the inadequacy of the latter even for tenuous clouds.

Gabriel, P. M.↗

NASA NTRS1995

Solar wind/local interstellar medium interaction including charge exchange with neural hydrogen

We present results from a hydrodynamic model of the interaction of the solar wind with the local interstellar medium (LISM), self-consistently taking into account the effects of charge exchange between the plasma component and the interstellar neutrals. The simulation is fully time dependent, and is carried out in two or three dimensions, depending on whether the helio-latitudinal dependence of the solar wind speed and number density (both giving rise to three dimensional effects) are included. As a first approximation it is assumed that the neutral component of the flow can be described by a single, isotropic fluid. Clearly, this is not the actual situation, since charge exchange with the supersonic solar wind plasma in the region of the nose results in a 'second' neutral fluid propagating in the opposite direction as that of the LISM neutrals.

Pauls, H. Louis↗

NASA NTRS1995

Locations of termination shock and heliopause based on Voyager plasma and magnetic field data

The locations of the termination shock and the heliopause are studied taking into account the effects of pickup protons. The study uses available plasma and magnetic field data from Voyagers over a 14-year period (1978-1991) and Voyager observation of the 1992-93 radio emission event. Outside 30 AU, pickup protons have a significant influence on dynamical structures of the outer heliosphere. The solar wind is treated as a mixture of electrons, solar wind protons, and interstellar pickup protons. If the magnitude of the interstellar magnetic field B(sub int) is given, one can quantitatively study the motion and location of the termination shock. The location is anti-correlated with the sun spot number and the shock has an average speed of approx. 24 km/s. Because B(sub int) is poorly known, additional information is needed in studying the termination shock. Cummings, et al. have used observations of anomalous cosmic rays to estimate the location of the shock. The observations of the 1991 GMIR and GMIR shock and the 1992-93 radio emission event provide another handle for the study of the termination shock and the heliopause. After its penetration through the termination shock, the GMIR shock continued to propagate in the subsonic region of the solar wind and eventually interacted with the heliopause. This interaction produces a transmitted shock propagating outward in the interstellar medium and a reflected shock propagating inward toward the sun in the subsonic solar wind. The plasma frequencies behind the reflected and the transmitted shock can be, respectively, responsible for the 2- and 3-kHz radio emissions. Taking into account the effects of pickup protons we found that the average locations of the termination shock and the heliopause in 1991-92 are at approximately 66 AU and 150 AU, respectively.

Whang, Y. C.↗

NASA NTRS2006

Observations of High-Redshift X-Ray Selected Clusters with the Sunyaev-Zel'dovich Array

We report measurements of the Sunyaev-Zel'dovich (SZ) effect in three highredshift (0.89 less than or equal to z less than or equal to 1.03), X-ray selected galaxy clusters. The observations were obtained at 30 GHz during the commissioning period of a new, eight-element interferometer - the Sunyaev-Zel dovich Array (SZA) - built for dedicated SZ effect observations. The SZA observations are sensitive to angular scales larger than those subtended by the virial radii of the clusters. Assuming isothermality and hydrostatic equilibrium for the intracluster medium, and gas-mass fractions consistent with those for clusters at moderate redshift, we calculate electron temperatures, gas masses, and total cluster masses from the SZ data. The SZ-derived masses, integrated approximately to the virial radii, are 1.9(sup +0.5)(sub -0.4) x 10(sup 14) solar mass for ClJ1415.1+3612, 3.4 (sup +0.6)(sup -0.5) x 10(sup 14) solar mass for ClJ1429.0+4241 and 7.2(sup +1.3)(sub -0.9) x 10(sup 14) solar mass for ClJ1226.9+3332. The SZ-derived quantities are in good agreement with the cluster properties derived from X-ray measurements.

Muchovej, Stephen↗

NASA NTRS2006

Observations of High-Redshift X-Ray Selected Clusters with the Sunyaev-Zel'dovich Array

We report measurements of the Sunyaev-Zel'dovich (SZ) effect in three high redshift (0.89 less than or equal to z less than or equal to 1.03), X-ray selected galaxy clusters. The observations were obtained at 30 GHz during the commissioning period of a new, eight-element interferometer - the Sunyaev-Zel'dovich Array (SZA) - built for dedicated SZ effect observations. The SZA observations are sensitive to angular scales larger than those subtended by the virial radii of the clusters. Assuming isothermality and hydrostatic equilibrium for the intracluster medium, and gas-mass fractions consistent with those for clusters at moderate redshift, we calculate electron temperatures, gas masses, and total cluster masses from the SZ data. The SZ-derived masses, integrated approximately to the virial radii, are 1.9 (sup +0.5)(sub -0.4) x 10(exp 14) solar mass for Cl J1415.1+3612, 3.4 (sup +0.6)(sub -0.5) x 10(exp 14) solar mass for Cl J1429.0+4241 and 7.2 (sup +1.3)(sub -0.9) x 10(exp 14) solar mass for Cl J1226.9+3332. The SZ-derived quantities are in good agreement with the cluster properties derived from X-ray measurements.

Muchovej, Stephen↗

DOE OSTI2026

A comprehensive first-principles study of the effects of the exchange-correlation functional and magnetism on defect and diffusion properties of the CoCrNi medium-entropy alloy

The present work is a novel, systematic study of the effect of density functional theory input parameters on the vacancy formation energy (VFE), migration barrier for diffusion, and electronic structure for each element in the CoCrNi medium-entropy alloy (MEA). In particular, the novelties include: (1) calculating the aforementioned properties of Co, Cr, or Ni, in the CoCrNi MEA using magnetic and non-magnetic states, and two versions of the generalized gradient approximation: Perdew, Burke, and Ernzerhof (PBE) and the PBE version for solids (PBEsol), and (2) a detailed comparison of 0 K activation energy to experimental creep activation energies. First-principles calculations at 0 K are performed using the Vienna ab-initio simulation package. Special quasirandom structures (SQS) and Widom-type substitution are employed. For each element, Co, Cr, or Ni, non-magnetic calculations result in a higher VFE and larger range of calculated values for the configurations studied. The averaged migration barrier is the highest for Co in the CoCrNi for three of four sets of calculation parameters in the configurations studied. Finally, the results indicate that the average 0 K activation energy for diffusion makes up 70–80% of the experimental creep activation energy, depending on the exchange-correlation functional employed.

36 MATERIALS SCIENCE↗

NASA NTRS2022

Numerical Results for Polarized Light Scattering in a Spherical Atmosphere

We report numerical results for polarized light reflection from the top of a Rayleigh scattering spherical atmosphere with height-dependent single scattering albedo over a dark surface. Michael Mishchenko considered this scenario back in the 1990’s, for a plane-parallel atmosphere of unit optical thickness (OT), for which radiance errors arising from neglecting polarization reaches their highest values. To further extend Mishchenko’s results, we consider a value of OT = 0.25, for which the effect of atmospheric curvature is pronounced. New results are generated using three state-of-the art radiative transfer (RT) codes. These are: the MYSTIC and MCSSA models, which simulate light scattering in a true- spherical atmosphere using Monte Carlo methods; and the discrete ordinate code VLIDORT, operating with a new multiple-scatter spherical correction designed to deliver reasonable approximations to spherical-medium scattering. In this work, we report results for both single and multiple scattering; this will help to support the validation of existing and future polarized spherical RT codes, especially those using approximative methods to deal with sphericity.

S Korkin↗

NASA NTRS1990

Heating and melting of small icy satellites by the decay of 26Al

We study the effect of radiogenic heating due to 26Al on the thermal evolution of small icy satellites. Our object is to find the extent of internal melting as a function of the satellite radius and of the initial 26Al abundance. The implicit assumption, based on observations of young stars, is that planet and satellite accretion occurred on a time scale of approximately 10(6) yr (comparable with the lifetime of 26Al). The icy satellites are modeled as spheres of initially amorphous ice, with chondritic abundances of 40K, 232Th, 235U, 238U, corresponding to an ice/dust mass ratio of 1. Evolutionary calculations are carried out, spanning 4.5 x 10(9) yr, for different combinations of the two free parameters. Heat transfer by subsolidus convection is neglected for these small satellites. Our main conclusion is that the initial 26Al abundance capable of melting icy bodies of satellite size to a significant extent is more than 10 times lower than that prevailing in the interstellar medium (or that inferred from the Ca-Al rich inclusions of the Allende meteorite, approximately 7 x 10(-7) by mass). We find, for example, that an initial 26Al mass fraction of approximately 4 x 10(-8) is sufficient for melting almost completely icy spheres with radii of 800 km, typical of the larger icy planetary satellites. We also find that for any given 26Al abundance, there is a narrow range of radii below which only marginal melting occurs and above which most of the ice melts (and refreezes later). Since extensive melting may have important consequences, such as differentiation, gas release, and volcanic activity, the effect of 26Al should be included in future studies of satellite interiors.

Non-NASA Center↗

DOE OSTI2025

The role of specimen size and grain boundary characteristics in the yield strength of tungsten in microtensile tests

To effectively use the measured properties from small-scale tensile tests for bulk material performance predictions, it is essential to understand the threshold of specimen size-effect strengthening and the interaction between dislocations and microstructures within miniaturized specimens. This study uses pure tungsten to investigate the size effect in terms of specimen size, grain size, and grain boundary characteristics relative to the yield strength of tungsten at room temperature. We evaluate the transition from miniaturized specimen properties to bulk properties and the deformation behavior through small-scale tensile tests of three specimen sizes (large: 80 × 100 × 233 µm³; medium: 7 × 7 × 18 µm³; and small: 2 × 2 × 5 µm³). The testing results reveal that the small and medium specimens exhibit high yield strength with ductile behavior, while the large specimens exhibit brittle failure, consistent with the room temperature strength of tungsten, indicating bulk behavior. We further explore the specimen size-effect sensitivity to yield stress and the scaling relationship between yield strength and the number of grains involved in the deformation. A power-law relationship with the exponent value of approximately -0.5 was found in the yield strength–grain number scaling, implying the Hall-Petch like behavior. A minimum of 7–17 effective grain boundaries across the tensile gauge dimension is required to accurately measure bulk properties.

36 - MATERIALS SCIENCE↗