Search NASA⌕ Search

SEARCH · Search NASA

Results for “Discrete dipole”

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 37 records · Page 2

Fractal dust grains around R Coronae Borealis stars

Discrete dipole approximation calculations of the optical properties of random fractal aggregates of graphite spheroids show a UV absorption feature that is too wide and centered at too long a wavelength to fit the observed interstellar 2200-A feature, but which is a good match to the 2400-A feature seen in the hydrogen-deficient R CrB stars reported by Hecht et al. (1984). Graphite fractal grains also match the UV bump and large long-wavelenvth extinction seen in laboratory studies of carbon smoke published by Bussoletti et al. (1987), which are usually attributed to amorphous carbon.

Wright, Edward L.↗

The effects of compositional inhomogeneities and fractal dimension on the optical properties of astrophysical dust

Mathis and Whiffen (1989) have recently suggested that interstellar dust particles are fluffy aggregates of submicron-size particles composed of various astronomical minerals. These dust particles should exhibit optical properties that are quite different from standard dust, characterized by spherical particles of various homogeneous mineral composition. In this paper, the discrete dipole approximation (DDA) method is used to examine the effects of chemical inhomogeneities and spatial structure on the optical properties of interstellar Mathis-Whiffen-type dust particles. The spatial structure of the dust is represented by its fractal dimension, and the chemical inhomogeneities are simulated by randomly assigning the composition of the occupied sites in the structure to be either carbon or silicate. It is found that compositional inhomogeneities are the dominant parameter affecting the shape of the 9.7 and 18 micron silicate bands. Some bands-shape variations can be attributed to the fractal dimension of the dust. The results derived here can be used to explain or constrain variations in these parameters among various astronomical objects.

Bazell, David↗

Radiative properties of visible and subvisible Cirrus: Scattering on hexagonal ice crystals

One of the main objectives of the First International Satellite Cloud Climatology Project (ISCCP) Regional Experiment (FIRE) is to provide a better understanding of the physics of upper level clouds. The focus is on just one specific aspect of cirrus physics, namely on characterizing the radiative properties of single, nonspherical ice particles. The basis for further more extensive studies of the radiative transfer through upper level clouds is provided. Radiation provides a potential mechanism for strong feedback between the divergence of in-cloud radiative flux and the cloud microphysics and ultimately on the dynamics of the cloud. Some aspects of ice cloud microphysics that are relevant to the radiation calculations are described. Next, the Discrete Dipole Approximation (DDA) is introduced and some new results of scattering by irregular crystals are presented. The Anomalous Diffraction Theory (ADT) was adopted to investigate the scattering properties of even larger crystals. In this way the scattering properties of nonspherical particles were determined over a range of particle sizes.

Flatau, Piotr J.↗

Millimeter wave radiative transfer studies for precipitation measurements

Scattering calculations using the discrete dipole approximation and vector radiative transfer calculations were performed to model multiparameter radar return and passive microwave emission for a simple model of a winter storm. The issue of dendrite riming was addressed by computing scattering properties of thin ice disks with varying bulk density. It was shown that C-band multiparameter radar contains information about particle density and the number concentration of the ice particles. The radiative transfer modeling indicated that polarized multifrequency passive microwave emission may be used to infer some properties of ice hydrometers. Detailed radar modeling and vector radiative transfer modeling is in progress to enhance the understanding of simultaneous radar and radiometer measurements, as in the case of the proposed TRMM field program. A one-dimensional cloud model will be used to simulate the storm structure in detail and study the microphysics, such as size and density. Multifrequency polarized radiometer measurements from the SSMI satellite instrument will be analyzed in relation to dual-frequency and dual-polarization radar measurements.

Vivekanandan, J.↗

Modeling composite and fluffy grains: The effects of porosity

Recent studies of interplanetary and interstellar dust provide evidence that cosmic dust grains are fluffy, composite objects, highlighting the need for models of the electromagnetic scattering by these grains. Effective medium theory (EMT) with Mie-type series solutions has been used to explore the effects of porosity which would be important in composite dust particles. While this indirect approach is both flexible and computationally efficient, it is not necessarily a good approximation. The need for EMT and its rather restrictive assumptions may be circumvented through a direct computation of the scattering properties via finite element methods, such as the discrete dipole approximation (DDA). Recently, the utility of the DDA method has been advanced significantly through improvements in theory, in numerical algorithms, and in computer hardware. Extensive calculations with the DDA method are used here to examine more directly the effects of porosity. A particular emphasis is placed upon developing a valid methodology. For both solid and porous targets we establish both numerical and physical convergence properties over the range of size parameter that is required for our study. DDA cross sections for grains with a range of porosity are compared to those computed by the EMT/series expansion technique to examine the applicability of several mixing rules, including two extensions of the Bruggeman rule. We show that for particles with Rayleigh vacuum inclusions, the extension proposed by Rouleau & Martin is quite successful. We also investigate the effects of larger, non-Rayleigh vacuum inclusions for various levels of porosity and find that they can be significant.

Wolff, Michael J.↗

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

Snow Crystal Orientation Effects on the Scattering of Passive Microwave Radiation

For this study, consideration is given to the role crystal orientation plays in scattering and absorbing microwave radiation. A discrete dipole scattering model is used to measure the passive microwave radiation, at two polarizations (horizontal and vertical), scattered by snow crystals oriented in random and non random positions, having various sizes (ranging between 1 micrometers to 10,000 micrometers in radius), and shapes (including spheroids, cylinders, hexagons). The model results demonstrate that for the crystal sizes typically found in a snowpack, crystal orientation is insignificant compared to crystal size in terms of scattering microwave energy in the 8,100 gm (37 GHz) region of the spectrum. Therefore, the assumption used in radiative transfer approaches, where snow crystals are modeled as randomly oriented spheres, is adequate to account for the transfer of microwave energy emanating from the ground and passing through a snowpack.

Foster, J. L.↗

Studies of Radiation and Microphysics in Cirrus and Marine Stratocumulus Clouds

Two tasks were completed during this period. In the first, we examined the polarization of millimeter-wavelength radar beams scattered by ice crystals. Because of their non-spherical shape and size, ice crystals depolarize the incident polarized radar beam. In principle, this depolarization can be used to identify ice from liquid water, as well as provide some information on size. However, the amount of de-polarization is small, producing only a weak signal at the receiver. Our task was to determine the magnitude of such a signal and decide if our radar would be capable of measuring it under typical cirrus conditions. The theoretical study was carried out by Henrietta Lemke, a visiting graduate student from Germany. She had prior experience using a discrete dipole code to compute scattering depolarization. Dr. Kultegin Aydin of the Penn State Electrical Engineering Department, who is also expert in this area, consulted with us on this project at no cost to the project. Our conclusion was that the depolarization signal is too weak to be usefully measured by our system. Therefore we proceeded no further in this study. The second task involved the study of the effect of stratus microphysics on surface cloud forcing. Manajit Sengupta, a graduate student, and the project PI jointly carried out this task. The study used data culled from over a year of continuous radar and radiometer observations at the Atmospheric Radiation Measurement (ARM) site in Oklahoma. The study compared solar radiation calculations made using constant microphysics with calculations made using a retrieved mean particle size. The results showed that on average the constant microphysics produced the correct forcing when compared with the observed forcing. We conclude, therefore, that there is little impetus on radiation grounds alone to include explicit microphysics in climate models. The question of pollutant particle emission impacts on microphysics remains to be resolved. A manuscript is in preparation and will be submitted this year.

Source record↗

A New Method for Estimating the Single Scattering Phase Functions of Regolith Grains

Hapke's photometric model has been widely used in solar system remote sensing applications for nearly two decades. Recently, Hapke extended his model to describe the coherent-backscatter opposition effect and multiple-scattering by particles with anisotropic single particle phase functions (SPPF's). A practical difficulty for retrieving Hapke's model parameters from typical planet, satellite, and asteroid photometry data sets is that the model employs a large number of adjustable parameters (at least eight) that can be reliably constrained only for a small number of planetary data sets in which both disk-resolved and whole-disk observations are available from opposition to very large phase angles. The present work aims to reduce the number of adjustable parameters and preserve (or even enhance) the model's accuracy and usefulness by expressing Hapke's parameters in terms of more fundamental physical properties on which they mutually depend. The most difficult part of this task, described here, is to develop a simple method for computing the effective SPPF for structurally complex regolith grains from optical constants, grain-size distribution, and average regolith porosity. The development of light-scattering models for irregularly shaped particles is a large, complex subject and many sophisticated methods, such as Discrete Dipole Approximation (DDA) and Monte-Carlo simulations, have been explored elsewhere. Many of these methods are computationally intensive and probably impractical for routine substitution in Hapke's model. Here, progress is reported in developing a practical, semi-empirical method for estimating the directional scattering behavior (i.e. SPPF) of irregular regolith grains. The method employs Optical Transfer Function (OTF) techniques to model how the structural complexity of regolith particles broaden and attenuate the angular distribution of scattered light relative to that expected from ideal spherical particles of equivalent size and composition.

Helfenstein, P.↗

Scattering Properties and Brightness Temperatures Associated with Solid Precipitation

In the past few years, early solid precipitation detection and retrieval algorithms have been developed and shown to be applicable for snowing clouds and blizzards. NOAA has an operational snow versus rain classifier based on AMSU-B observations. Solid precipitation retrieval algorithms reported in the literature over the past two years include those that rely on neural nets, statistics, or physical relationships. All of the algorithms require the use of millimeter-wave radiometer observations. The millimeter-wave frequencies are especially sensitive to the scattering and emission properties of frozen particles due to the ice particle refractive index. Passive radiometric channels respond to both the integrated particle mass throughout the volume and field of view, and to the amount, location, and size distribution of the frozen (and liquid) particles with the sensitivity varying for different frequencies and hydrometeor types. This investigation probes the sensitivity of scattering and absorption coefficients, and hence computed brightness temperatures, resulting from variations in solid precipitation cloud profiles. The first study compares the single scattering, absorption, and asymmetry parameters associated with snow particles in clouds. Several methodologies are used to convert the physical characteristics (e.g., shape, size distributions, ice-air-water ratios) of ice particles to electromagnetic properties (e.g., absorption, scattering, and asymmetry factors). These methodologies include: conversion to solid ice particles, homogeneous dielectric mixing, or discrete dipole approximation. Changes in the conversion methodology can produce computed brightness temperature differences greater than 50 Kelvin.

Skofronick-Jackson, Gail M.↗

Comet Mineralogy as Inferred from Infrared Spectra of Comets

For most comets, infrared (IR) spectroscopy (remote sensing) is the method through which we diagnose the mineralogy and size distribution of dust in their comae. The shape and contrast of the IR spectral features depend on the particle size: optically active minerals (absorbing of visible and near-IR solar photons) and submicron solid grains or highly porous (> 90% vacuum) grains primarily contribute to the shapes of the observed resonances. Comet mineralogies typically are determined by fitting thermal emission models of ensembles of discrete mineral grains to observed IR spectral energy distributions. The absorptivities (Q-abs) and scattering efficiencies (Q-scat) of the discrete mineral grains are computed using Mie scattering, Maxwell-Garnet mixing, Discrete Dipole Approximation, and Multi-Layered Sphere codes. These techniques when applied to crystalline minerals, specifically olivine (Mg_x, Fe_1-x)2 Si04, x>0.9, require the use of ellipsoidal shaped particles with elongated axial ratios or hollow spheres to produce the shapes of the resonances observed both from comet comae and laboratory samples. The wavelength positions of the distinct resonances from submicron-radii crystalline silicates, as well as their thermal equilibrium temperatures, constrain the crystalline olivine to have a relatively high Mg-content (x>0.9, or Fo>90). Only resonances computed for submicron Mg-rich crystalline olivine and crystalline orthopyroxene match the observed IR spectral features. However, this has led to the interpretation that micron-radii and larger crystals are absent from comet comae. Furthermore, the mass fraction of silicate crystals is dependent upon whether just the submicron portion of the size distribution is being compared or the submicron crystals compare to the aggregates of porous amorphous silicates that are computationally tractable as porous spheres. We will discuss the Deep Impact results as examples of these challenges to interpreting mid-IR spectra of short-period comets.

Wooden, Diane H.↗

A Physical Model to Estimate Snowfall over Land using AMSU-B Observations

In this study, we present an improved physical model to retrieve snowfall rate over land using brightness temperature observations from the National Oceanic and Atmospheric Administration's (NOAA) Advanced Microwave Sounder Unit-B (AMSU-B) at 89 GHz, 150 GHz, 183.3 +/- 1 GHz, 183.3 +/- 3 GHz, and 183.3 +/- 7 GHz. The retrieval model is applied to the New England blizzard of March 5, 2001 which deposited about 75 cm of snow over much of Vermont, New Hampshire, and northern New York. In this improved physical model, prior retrieval assumptions about snowflake shape, particle size distributions, environmental conditions, and optimization methodology have been updated. Here, single scattering parameters for snow particles are calculated with the Discrete-Dipole Approximation (DDA) method instead of assuming spherical shapes. Five different snow particle models (hexagonal columns, hexagonal plates, and three different kinds of aggregates) are considered. Snow particle size distributions are assumed to vary with air temperature and to follow aircraft measurements described by previous studies. Brightness temperatures at AMSU-B frequencies for the New England blizzard are calculated using these DDA calculated single scattering parameters and particle size distributions. The vertical profiles of pressure, temperature, relative humidity and hydrometeors are provided by MM5 model simulations. These profiles are treated as the a priori data base in the Bayesian retrieval algorithm. In algorithm applications to the blizzard data, calculated brightness temperatures associated with selected database profiles agree with AMSU-B observations to within about +/- 5 K at all five frequencies. Retrieved snowfall rates compare favorably with the near-concurrent National Weather Service (NWS) radar reflectivity measurements. The relationships between the NWS radar measured reflectivities Z(sub e) and retrieved snowfall rate R for a given snow particle model are derived by a histogram matching technique. All of these Z(sub e)-R relationships fall in the range of previously established Z(sub e)-R relationships for snowfall. This suggests that the current physical model developed in this study can reliably estimate the snowfall rate over land using the AMSU-B measured brightness temperatures.

Kim, Min-Jeong↗

Exploring Alternate Parameterizations for Snowfall with Validation from Satellite and Terrestrial Radars

Increases in computational resources have allowed operational forecast centers to pursue experimental, high resolution simulations that resolve the microphysical characteristics of clouds and precipitation. These experiments are motivated by a desire to improve the representation of weather and climate, but will also benefit current and future satellite campaigns, which often use forecast model output to guide the retrieval process. Aircraft, surface and radar data from the Canadian CloudSat/CALIPSO Validation Project are used to check the validity of size distribution and density characteristics for snowfall simulated by the NASA Goddard six-class, single-moment bulk water microphysics scheme, currently available within the Weather Research and Forecast (WRF) Model. Widespread snowfall developed across the region on January 22, 2007, forced by the passing of a midlatitude cyclone, and was observed by the dual-polarimetric, C-band radar King City, Ontario, as well as the NASA 94 GHz CloudSat Cloud Profiling Radar. Combined, these data sets provide key metrics for validating model output: estimates of size distribution parameters fit to the inverse-exponential equations prescribed within the model, bulk density and crystal habit characteristics sampled by the aircraft, and representation of size characteristics as inferred by the radar reflectivity at C- and W-band. Specified constants for distribution intercept and density differ significantly from observations throughout much of the cloud depth. Alternate parameterizations are explored, using column-integrated values of vapor excess to avoid problems encountered with temperature-based parameterizations in an environment where inversions and isothermal layers are present. Simulation of CloudSat reflectivity is performed by adopting the discrete-dipole parameterizations and databases provided in literature, and demonstrate an improved capability in simulating radar reflectivity at W-band versus Mie scattering assumptions.

Molthan, Andrew L.↗

Exploring Alternative Parameterizations for Snowfall with Validation from Satellite and Terrestrial Radars

Increases in computational resources have allowed operational forecast centers to pursue experimental, high resolution simulations that resolve the microphysical characteristics of clouds and precipitation. These experiments are motivated by a desire to improve the representation of weather and climate, but will also benefit current and future satellite campaigns, which often use forecast model output to guide the retrieval process. The combination of reliable cloud microphysics and radar reflectivity may constrain radiative transfer models used in satellite simulators during future missions, including EarthCARE and the NASA Global Precipitation Measurement. Aircraft, surface and radar data from the Canadian CloudSat/CALIPSO Validation Project are used to check the validity of size distribution and density characteristics for snowfall simulated by the NASA Goddard six-class, single moment bulk water microphysics scheme, currently available within the Weather Research and Forecast (WRF) Model. Widespread snowfall developed across the region on January 22, 2007, forced by the passing of a mid latitude cyclone, and was observed by the dual-polarimetric, C-band radar King City, Ontario, as well as the NASA 94 GHz CloudSat Cloud Profiling Radar. Combined, these data sets provide key metrics for validating model output: estimates of size distribution parameters fit to the inverse-exponential equations prescribed within the model, bulk density and crystal habit characteristics sampled by the aircraft, and representation of size characteristics as inferred by the radar reflectivity at C- and W-band. Specified constants for distribution intercept and density differ significantly from observations throughout much of the cloud depth. Alternate parameterizations are explored, using column-integrated values of vapor excess to avoid problems encountered with temperature-based parameterizations in an environment where inversions and isothermal layers are present. Simulation of CloudSat reflectivity is performed by adopting the discrete-dipole parameterizations and databases provided in literature, and demonstrate an improved capability in simulating radar reflectivity at W-band versus Mie scattering assumptions.

Molthan, Andrew L.↗

Modeling Optical Properties of Mineral Aerosol Particles by Using Nonsymmetric Hexahedra

We explore the use of nonsymmetric geometries to simulate the single-scattering properties of airborne dust particles with complicated morphologies. Specifically, the shapes of irregular dust particles are assumed to be nonsymmetric hexahedra defined by using the Monte Carlo method. A combination of the discrete dipole approximation method and an improved geometric optics method is employed to compute the single-scattering properties of dust particles for size parameters ranging from 0.5 to 3000. The primary optical effect of eliminating the geometric symmetry of regular hexahedra is to smooth the scattering features in the phase function and to decrease the backscatter. The optical properties of the nonsymmetric hexahedra are used to mimic the laboratory measurements. It is demonstrated that a relatively close agreement can be achieved by using only one shape of nonsymmetric hexahedra. The agreement between the theoretical results and their measurement counterparts can be further improved by using a mixture of nonsymmetric hexahedra. It is also shown that the hexahedron model is much more appropriate than the "equivalent sphere" model for simulating the optical properties of dust particles, particularly, in the case of the elements of the phase matrix that associated with the polarization state of scattered light.

Bi, Lei↗

Collisional Histories of Comets and Trojan Asteroids: Diopside, Magnesite, and Fayalite Impact Studies

Comets and asteroids have weathered dynamic histories, as evidenced by their rough surfaces. The Nice model describes a violent reshuffling of small bodies during the Late Heavy Bombardment, with collisions acting to grind these planetesimals away. This creates an additional source of impact material that can re-work the surfaces of the larger bodies over the lifetime of the solar system. Here, we investigate the possibility that signatures due to impacts (e.g. from micrometeoroids or meteoroids) could be detected in their spectra, and how that can be explained by the physical manifestation of shock in the crystalline structure of minerals. All impact experiments were conducted in the Johnson Space Center Experimental Impact Laboratory using the vertical gun. Impact speeds ranged from approx.2.0 km/s to approx.2.8 km/s. All experiments were conducted at room temperature. Minerals found in comets and asteroids were chosen as targets, including diopside (MgCaSi2O6, monoclinic pyroxene), magnesite (MgCO3, carbonate), and fayalite (FeSiO4, olivine). Impacted samples were analyzed using a Fourier Transform Infrared Spectrometer (FTIR) and a Transmission Electron Microscope (TEM). Absorbance features in the 8-13 m spectral region demonstrate relative amplitude changes as well as wavelength shifts. Corresponding TEM images exhibit planar shock dislocations in the crystalline structure, attributed to deformation at high strain and low temperatures. Elongating or shortening the axes of the crystalline structure of forsterite (Mg2SiO4, olivine) using a discrete dipole approximation model (Lindsay et al., submitted) yields changes in spectral features similar to those observed in our impacted laboratory minerals.

Lederer, S. M.↗

A Sensitivity Study on the Effects of Particle Chemistry, Asphericity and Size on the Mass Extinction Efficiency of Mineral Dust in the Earth's Atmosphere: From the Near to Thermal IR

To determine a plausible range of mass extinction efficiencies (MEE) of terrestrial atmospheric dust from the near to thermal IR, sensitivity analyses are performed over an extended range of dust microphysical and chemistry perturbations. The IR values are subsequently compared to those in the near-IR, to evaluate spectral relationships in their optical properties. Synthesized size distributions consistent with measurements, model particle size, while composition is defined by the refractive indices of minerals routinely observed in dust, including the widely used OPAC/Hess parameterization. Single-scattering properties of representative dust particle shapes are calculated using the T-matrix, Discrete Dipole Approximation and Lorenz-Mie light-scattering codes. For the parameterizations examined, MEE ranges from nearly zero to 1.2 square meters per gram, with the higher values associated with non-spheres composed of quartz and gypsum. At near-IR wavelengths, MEE for non-spheres generally exceeds those for spheres, while in the thermal IR, shape-induced changes in MEE strongly depend on volume median diameter (VMD) and wavelength, particularly for MEE evaluated at the mineral resonant frequencies. MEE spectral distributions appear to follow particle geometry and are evidence for shape dependency in the optical properties. It is also shown that non-spheres best reproduce the positions of prominent absorption peaks found in silicates. Generally, angular particles exhibit wider and more symmetric MEE spectral distribution patterns from 8-10 micrometers than those with smooth surfaces, likely due to their edge-effects. Lastly, MEE ratios allow for inferring dust optical properties across the visible-IR spectrum. We conclude the MEE of dust aerosol are significant for the parameter space investigated, and are a key component for remote sensing applications and the study of direct aerosol radiative effects.

Hansell, R. A., Jr.↗

Efficient Implementation of the Invariant Imbedding T-Matrix Method and the Separation of Variables Method Applied to Large Nonspherical Inhomogeneous Particles

Three terms, ''Waterman's T-matrix method'', ''extended boundary condition method (EBCM)'', and ''null field method'', have been interchangeable in the literature to indicate a method based on surface integral equations to calculate the T-matrix. Unlike the previous method, the invariant imbedding method (IIM) calculates the T-matrix by the use of a volume integral equation. In addition, the standard separation of variables method (SOV) can be applied to compute the T-matrix of a sphere centered at the origin of the coordinate system and having a maximal radius such that the sphere remains inscribed within a nonspherical particle. This study explores the feasibility of a numerical combination of the IIM and the SOV, hereafter referred to as the IIMþSOV method, for computing the single-scattering properties of nonspherical dielectric particles, which are, in general, inhomogeneous. The IIMþSOV method is shown to be capable of solving light-scattering problems for large nonspherical particles where the standard EBCM fails to converge. The IIMþSOV method is flexible and applicable to inhomogeneous particles and aggregated nonspherical particles (overlapped circumscribed spheres) representing a challenge to the standard superposition T-matrix method. The IIMþSOV computational program, developed in this study, is validated against EBCM simulated spheroid and cylinder cases with excellent numerical agreement (up to four decimal places). In addition, solutions for cylinders with large aspect ratios, inhomogeneous particles, and two-particle systems are compared with results from discrete dipole approximation (DDA) computations, and comparisons with the improved geometric-optics method (IGOM) are found to be quite encouraging.

Invariant imbedding method↗