Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “scattering matrix simulations”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Probabilistic Modeling of Ceramic Matrix Composite Strength

Uncertainties associated with the primitive random variables such as manufacturing process (processing temperature, fiber volume ratio, void volume ratio), constituent properties (fiber, matrix and interface), and geometric parameters (ply thickness, interphase thickness) have been simulated to quantify the scatter in the first matrix cracking strength (FMCS) and the ultimate tensile strength of SCS-6/RBSN (SiC fiber (SCS-6) reinforced reaction-bonded silicon nitride composite) ceramic matrix composite laminate at room temperature. Cumulative probability distribution function for the FMCS and ultimate tensile strength at room temperature (RT) of (0)(sub 8), (0(sub 2)/90(sub 2), and (+/-45(sub 2))(sub S) laminates have been simulated and the sensitivity of primitive variables to the respective strengths have been quantified. Computationally predicted scatter of the strengths for a uniaxial laminate have been compared with those from limited experimental data. Also the experimental procedure used in the tests has been described briefly. Results show a very good agreement between the computational simulation and the experimental data. Dominating failure modes in (0)(sub 8), (0/90)(sub s) and (+/-45)(sub S) laminates have been identified. Results indicate that the first matrix cracking strength for the (0)(sub S), and (0/90)(sub S) laminates is sensitive to the thermal properties, modulus and strengths of both the fiber and matrix whereas the ultimate tensile strength is sensitive to the fiber strength and the fiber volume ratio. In the case of a (+/-45)(sub S), laminate, both the FMCS and the ultimate tensile strengths have a small scatter range and are sensitive to the fiber tensile strength as well as the fiber volume ratio.

Shan, Ashwin R.↗
NASA NTRS1991

Hypercube technology

The JPL designed MARKIII hypercube supercomputer has been in application service since June 1988 and has had successful application to a broad problem set including electromagnetic scattering, discrete event simulation, plasma transport, matrix algorithms, neural network simulation, image processing, and graphics. Currently, problems that are not homogeneous are being attempted, and, through this involvement with real world applications, the software is evolving to handle the heterogeneous class problems efficiently.

Parker, Jay W.↗
NASA NTRS1995

The use of optimal polarizations for studying the microphysics of precipitation: Nonattenuating wavelengths

The objective of this work is to explore relationships between the microphysical properties of precipitation and optimal polarizations. The dependence of three optimal polarization parameters (asymmetry ratio A, optimal tilt tau(sub op), and optimal ellipticity epsilon(sub op) on the reflectivity-weighted mean drop shape, mean canting angle, and standard deviation of a Gaussian canting angle distribution is studied. This is accomplished by using computer simulations that provide the rms scattering matrix for an ensemble of canted drops with a prescribed two-parameter canting angle distribution. Also examined are the effects of propagation on the polarization parameters for nonattenuating wavelengths. The asymmetry ratio A is simply the ratio of the maximal to minimal total backscattered energy (ratio of the largest and smallest eigenvalue of the Graves power matrix G is identically equal to S dagger S). Similar to Z(sub DR), this ratio decreases with increasing mean axial ratio, but unlike Z(sub DR), it is not affected by canting (for a single drop). The dependence of A on the reflectivity-weighted mean drop shape is examined, and a power-law relationship similar to that which exists for Z(sub DR) is established. The asymmetry ratio A can be regarded as a generalization of Z(sub DR) because it requires only a measurement of linear depolarization ratio (in addition to Z(sub DR)), is independent of the propagation phase, and is less sensitive to canting. In a similar manner, the dependence of optimal ellipticity and tilt on the microphysical parameters is studied. In particular, it appears that the rms tilt of the optimal polarization ellipse is proportional to the variance of the canting angle distribution. Several other promising relationships between optimal polarizations and the microphysical variables of an ensemble of hydrometeors are also discussed.

Kwiatkowski, John M.↗
NASA NTRS2022

Multiple Scattering Effects of Trees on L-band Microwave Using a Hybrid Method

The multiple scattering of electromagnetic waves in forests is studied using a two-step hybrid method. First, the T-matrix of a single tree is calculated based on its far-field computed with the full-wave simulations of the FEKO software. The T-matrix captures the multiple scattering caused by the tree structure. Second, the the interactions among different trees are considered using the T-matrices of the individual trees and the Foldy-Lax equations. The result of the two-step hybrid method is validated with FEKO by solving scattering from two trees directly. The multiple scattering effects are illustrated by the field solutions.

Yueh, Simon↗
NASA NTRS1977

Extending radiative transfer models by use of Bayes rule

This paper presents a procedure that extends some existing radiative transfer modeling techniques to problems in atmospheric science where curvature and layering of the medium and dynamic range and angular resolution of the signal are important. Example problems include twilight and limb scan simulations. Techniques that are extended include successive orders of scattering, matrix operator, doubling, Gauss-Seidel iteration, discrete ordinates and spherical harmonics. The procedure for extending them is based on Bayes' rule from probability theory.

Whitney, C.↗
NASA NTRS2022

Full-Wave Simulations of Scattering by Corn Fields at L-Band

In this paper, the Numerical Maxwell Model of 3D (NMM3D) full-wave simulation is performed over a corn field using a hybrid method to study the vegetation effect on the microwave. The commercial software of FEKO is used to extract T-matrix of single corn in the first step. Then the calculated T-matrix is combined with Wave Multiple Scattering Theory (W-MST) in the second step to consider the multiple scattering among different plants. The hybrid method is validated with HFSS by solving scattering from 2 corns. A corn field of 25 corn is simulated using the hybrid method and the transmission is calculated and compared with those obtained from the classical radiative transfer model.

Yueh, Simon↗
NASA NTRS1993

Probabilistic assessment of composite structures

A general computational simulation methodology for an integrated probabilistic assessment of composite structures is discussed and demonstrated using aircraft fuselage (stiffened composite cylindrical shell) structures with rectangular cutouts. The computational simulation was performed for the probabilistic assessment of the structural behavior including buckling loads, vibration frequencies, global displacements, and local stresses. The scatter in the structural response is simulated based on the inherent uncertainties in the primitive (independent random) variables at the fiber matrix constituent, ply, laminate, and structural scales that describe the composite structures. The effect of uncertainties due to fabrication process variables such as fiber volume ratio, void volume ratio, ply orientation, and ply thickness is also included. The methodology has been embedded in the computer code IPACS (Integrated Probabilistic Assessment of Composite Structures). In addition to the simulated scatter, the IPACS code also calculates the sensitivity of the composite structural behavior to all the primitive variables that influence the structural behavior. This information is useful for assessing reliability and providing guidance for improvement. The results from the probabilistic assessment for the composite structure with rectangular cutouts indicate that the uncertainty in the longitudinal ply stress is mainly caused by the uncertainty in the laminate thickness, and the large overlap of the scatter in the first four buckling loads implies that the buckling mode shape for a specific buckling load can be either of the four modes.

Shiao, Michael E.↗
NASA NTRS1997

Probabilistic Evaluation of Bolted Joints in Polymer Matrix Composites

Computational methods are described to probabilistically simulate fracture in bolted composite structures. Progressive fracture is simulated via an innovative approach independent of stress intensity factors and fracture toughness. The effect on structure damage of design variable uncertainties is quantified. The Fast Probability Integrator is used to assess the scatter in the composite structure response before and after damage. Sensitivity of the response to design variables is evaluated. The methods are demonstrated for bolted joint polymer matrix composite panels under end loads. The effects of fabrication process are included in the simulation of damage in the bolted panel. The results show that the most effective way to reduce the end displacement at fracture is to control the load and ply thickness.

Chamis, C. C.↗
NASA NTRS2011

Direct Simulation of Multiple Scattering by Discrete Random Media Illuminated by Gaussian Beams

The conventional orientation-averaging procedure developed in the framework of the superposition T-matrix approach is generalized to include the case of illumination by a Gaussian beam (GB). The resulting computer code is parallelized and used to perform extensive numerically exact calculations of electromagnetic scattering by volumes of discrete random medium consisting of monodisperse spherical particles. The size parameters of the scattering volumes are 40, 50, and 60, while their packing density is fixed at 5%. We demonstrate that all scattering patterns observed in the far-field zone of a random multisphere target and their evolution with decreasing width of the incident GB can be interpreted in terms of idealized theoretical concepts such as forward-scattering interference, coherent backscattering (CB), and diffuse multiple scattering. It is shown that the increasing violation of electromagnetic reciprocity with decreasing GB width suppresses and eventually eradicates all observable manifestations of CB. This result supplements the previous demonstration of the effects of broken reciprocity in the case of magneto-optically active particles subjected to an external magnetic field.

Mackowski, Daniel W.↗
NASA NTRS2021

Multi-Frequency NMM3D Simulations of Wave Propagation in Vegetation for Remote Sensing of Soil Moisture

To investigate the feasibility of using multi-frequencies to fur- ther improve the soil moisture retrieval, the recently devel- oped hybrid method is used to perform full-wave scattering simulations of a wheat field at three different frequencies of L-, S- and C- bands. In the hybrid method, the multiple scattering within a single wheat plant is first captured using the T-matrix of vec- tor cylindrical waves calculated from the full-wave solutions of HFSS. In the second step, the scattering among different plants is considered using the Foldy-Lax equations of multiple scattering theory (MST) and the extracted T-matrix. The transmission of microwaves through wheat field at the L-, S- and C-bands are calculated using the hybrid method. Results show that: (1) the transmission obtained from the full-wave simulations is much larger than those computed from the radiative transfer equations (RTE) model, (2) the hybrid method transmission has a weaker frequency dependence than those of the RTE, and (3) the attenuation caused by the vegetation layer would saturate with an increase in frequency.

Yueh, Simon↗
NASA NTRS2010

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↗
NASA NTRS1985

A General-applications Direct Global Matrix Algorithm for Rapid Seismo-acoustic Wavefield Computations

A new matrix method for rapid wave propagation modeling in generalized stratified media, which has recently been applied to numerical simulations in diverse areas of underwater acoustics, solid earth seismology, and nondestructive ultrasonic scattering is explained and illustrated. A portion of recent efforts jointly undertaken at NATOSACLANT and NORDA Numerical Modeling groups in developing, implementing, and testing a new fast general-applications wave propagation algorithm, SAFARI, formulated at SACLANT is summarized. The present general-applications SAFARI program uses a Direct Global Matrix Approach to multilayer Green's function calculation. A rapid and unconditionally stable solution is readily obtained via simple Gaussian ellimination on the resulting sparsely banded block system, precisely analogous to that arising in the Finite Element Method. The resulting gains in accuracy and computational speed allow consideration of much larger multilayered air/ocean/Earth/engineering material media models, for many more source-receiver configurations than previously possible. The validity and versatility of the SAFARI-DGM method is demonstrated by reviewing three practical examples of engineering interest, drawn from ocean acoustics, engineering seismology and ultrasonic scattering.

Schmidt, H.↗
NASA NTRS2001

Microwave Polarized Signatures Generated within Cloud Systems: SSM/I Observations Interpreted with Radiative Transfer Simulations

Special Sensor Microwave /Imager (SSM/I) observations in cloud systems are studied over the tropics. Over optically thick cloud systems, presence of polarized signatures at 37 and 85 GHz is evidenced and analyzed with the help of cloud top temperature and optical thickness extracted from visible and IR satellite observations. Scattering signatures at 85 GHz (TbV(85) less than or = 250 K) are associated with polarization differences greater than or = 6 K, approx. 50%, of the time over ocean and approx. 40% over land. In addition. over thick clouds the polarization difference at 37 GHz is rarely negligible. The polarization differences at 37 and 85 GHz do not stem from the surface but are generated in regions of relatively homogeneous clouds having high liquid water content. To interpret the observations, a radiative transfer model that includes the scattering by non-spherical particles is developed. based on the T-matrix approach and using the doubling and adding method. In addition to handling randomly and perfectly oriented particles, this model can also simulate the effect of partial orientation of the hydrometeors. Microwave brightness temperatures are simulated at SSM/I frequencies and are compared with the observations. Polarization differences of approx. 2 K can be simulated at 37 GHz over a rain layer, even using spherical drops. The polarization difference is larger for oriented non-spherical particles. The 85 GHz simulations are very sensitive to the ice phase of the cloud. Simulations with spherical particles or with randomly oriented non-spherical ice particles cannot replicate the observed polarization differences. However, with partially oriented non-spherical particles, the observed polarized signatures at 85 GHz are explained, and the sensitivity of the scattering characteristics to the particle size, asphericity, and orientation is analyzed. Implications on rain and ice retrievals are discussed.

Prigent, Catherine↗
NASA NTRS2018

Active and Passive Radiative Transfer Modeling of the Olympic Mountains Experiment

Sensor forward models are an important tool for interpreting remote sensing observations of geophysical phenomena. By implementing a three-dimensional framework, we can simulate and analyze observations from various sensors on disparate platforms. To demonstrate our model framework, we simulate observations from the Olympic Mountains Experiment (OLYMPEX). The use of cloud model simulations allows us to understand sensor response to cloud ice, falling snow, and other processes and features, and the application of model tools to observations allows us to quantify precipitation.MIIST 3D Forward ModelThe Multi-Instrument Inverse Solver Testbed(MIIST) uses the Atmospheric Radiative TransferSimulator (ARTS) for solving the vector radiativetransfer (RT) equation in up to three spatialdimensions within a spherical geometry• Gas absorptiono Line-by-line calculationso Fast transmittance tables• Hydrometeor scattering solverso Discrete ordinateo RT4 (Evans, 1D)o Radar Single Scattering (1D or 3D)o Monte Carlo (3D)Scattering TablesHigh-fidelity hydrometeor scatteringtables are necessary for accurateand consistent forward modeling ofmulti-frequency observations• Requires full Stokes matriceso And absorption vector• Randomly oriented particleso Discrete Dipole Approximationo Characteristic Basis Function Method(coming soon)• Horizontally-oriented plateso Invariant Imbedding T-matrix MethodCloud Resolving SimulationsCloud resolving simulations (e.g.,NU-WRF) supply output consistentwith ARTS needs• Atmospheric Informationo Temperatureo Pressure / heighto Water vapor• Hydrometeor Profileso ARTS architecture ripe for explicit binmicrophysics• Examples use Morrison 2M schemeThe Olympic Mountains Experiment (OLYMPEX)Validation for GPM of mid-latitudefrontal systems approaching nearcoastalmountains from the ocean• Large collection of ground-based andairborne sensorso Radarso Radiometerso In situ• Contemporaneous with RADEXo Two sets of radar at same frequenciesRadiometer Simulation (3 km NUWRF, 20151203, 15:00)2018.12.14 7Simulate 166 GHz polarizationdifference• Corresponds to the presence of aligned icecrystals• Look at trends for both simulations andobservations• Simulations can tolerate lower resolutiono Larger domainSimulations from Observations: OLYMPEXSimulate sensor response usinggeophysical retrievals as input• Single frequency radar retrievals• Multiple scattering enhancementapparent at W band• Spatially dependent phenomenonModeling Application: 1D Retrievals03 December 2015• DC-8 and ER-2 flightso Focus on APR-3 (DC-8)• Citationo Stacked microphysics legso Qualitative comparisonso Range of frozen habitso Presence of supercooledliquid cloudsResults• Retrievals match probeso Good qualitative match• Bands of increasedreflectivity correspond tolarge Dm and highaggregate fraction• Significant amounts ofsupercooled liquid water

Adams, Ian S.↗
NASA NTRS2012

Assessment of Polarization Effect on Efficiency of Levenberg-Marquardt Algorithm in Case of Thin Atmosphere Over Black Surface

The Levenberg-Marquardt algorithm [1, 2] provides a numerical iterative solution to the problem of minimization of a function over a space of its parameters. In our work, the Levenberg-Marquardt algorithm retrieves optical parameters of a thin (single scattering) plane parallel atmosphere irradiated by collimated infinitely wide monochromatic beam of light. Black ground surface is assumed. Computational accuracy, sensitivity to the initial guess and the presence of noise in the signal, and other properties of the algorithm are investigated in scalar (using intensity only) and vector (including polarization) modes. We consider an atmosphere that contains a mixture of coarse and fine fractions. Following [3], the fractions are simulated using Henyey-Greenstein model. Though not realistic, this assumption is very convenient for tests [4, p.354]. In our case it yields analytical evaluation of Jacobian matrix. Assuming the MISR geometry of observation [5] as an example, the average scattering cosines and the ratio of coarse and fine fractions, the atmosphere optical depth, and the single scattering albedo, are the five parameters to be determined numerically. In our implementation of the algorithm, the system of five linear equations is solved using the fast Cramer s rule [6]. A simple subroutine developed by the authors, makes the algorithm independent from external libraries. All Fortran 90/95 codes discussed in the presentation will be available immediately after the meeting from sergey.v.korkin@nasa.gov by request.

Korkin, S.↗
NASA NTRS2013

Spectrally Consistent Scattering, Absorption, and Polarization Properties of Atmospheric Ice Crystals at Wavelengths from 0.2 to 100 um

A data library is developed containing the scattering, absorption, and polarization properties of ice particles in the spectral range from 0.2 to 100 microns. The properties are computed based on a combination of the Amsterdam discrete dipole approximation (ADDA), the T-matrix method, and the improved geometric optics method (IGOM). The electromagnetic edge effect is incorporated into the extinction and absorption efficiencies computed from the IGOM. A full set of single-scattering properties is provided by considering three-dimensional random orientations for 11 ice crystal habits: droxtals, prolate spheroids, oblate spheroids, solid and hollow columns, compact aggregates composed of eight solid columns, hexagonal plates, small spatial aggregates composed of 5 plates, large spatial aggregates composed of 10 plates, and solid and hollow bullet rosettes. The maximum dimension of each habit ranges from 2 to 10,000 microns in 189 discrete sizes. For each ice crystal habit, three surface roughness conditions (i.e., smooth, moderately roughened, and severely roughened) are considered to account for the surface texture of large particles in the IGOM applicable domain. The data library contains the extinction efficiency, single-scattering albedo, asymmetry parameter, six independent nonzero elements of the phase matrix (P11, P12, P22, P33, P43, and P44), particle projected area, and particle volume to provide the basic single-scattering properties for remote sensing applications and radiative transfer simulations involving ice clouds. Furthermore, a comparison of satellite observations and theoretical simulations for the polarization characteristics of ice clouds demonstrates that ice cloud optical models assuming severely roughened ice crystals significantly outperform their counterparts assuming smooth ice crystals.

Optical properties↗
NASA NTRS2017

Detection of the Pairwise Kinematic Sunyaev-Zel'dovich Effect with BOSS DR11 and the Atacama Cosmology Telescope

We present a new measurement of the kinematic Sunyaev-Zel'dovich effect using data from the Atacama Cosmology Telescope (ACT) and the Baryon Oscillation Spectroscopic Survey (BOSS). Using 600 square degrees of overlapping sky area, we evaluate the mean pairwise baryon momentum associated with the positions of 50,000 bright galaxies in the BOSS DR11 Large Scale Structure catalog. A non-zero signal arises from the large-scale motions of halos containing the sample galaxies. The data fits an analytical signal model well, with the optical depth to microwave photon scattering as a free parameter determining the overall signal amplitude. We estimate the covariance matrix of the mean pairwise momentum as a function of galaxy separation, using microwave sky simulations, jackknife evaluation, and bootstrap estimates. The most conservative simulation-based errors give signal-to-noise estimates between 3.6 and 4.1 for varying galaxy luminosity cuts. We discuss how the other error determinations can lead to higher signal-to-noise values, and consider the impact of several possible systematic errors. Estimates of the optical depth from the average thermal Sunyaev-Zel'dovich signal at the sample galaxy positions are broadly consistent with those obtained from the mean pairwise momentum signal.

De Bernardis, F.↗
NASA NTRS2021

Polarimetric sensitivity of light-absorbing carbonaceous aerosols over ocean: A theoretical assessment

Visible-light-absorbing carbonaceous aerosols within the boundary layer affect the radiance and polarization states of the radiation at the top of the atmosphere. Remote sensing from suborbital and satellite-based platforms utilizes these radiance and polarization signals to retrieve the key properties of these aerosols. Recent retrieval algorithms have shown a progressive trend toward including multi-angular and multi-spectral polarimetric measurements to produce better retrieval accuracy in comparison to those using measurements based on a single viewing angle. Here, we perform a theoretical investigation of the top of atmosphere (TOA) radiance-related reflectance factor (bidirectional reflectance factor (BRF)) and the two types of polarimetry-related factors (polarized bidirectional reflectance factor (pBRF) and the degree of linear polarization (DoLP)) for different types of atmospheric light-absorbing carbonaceous aerosols as a function of particle size distribution. We selected three polarimetric bands corresponding to those utilized by NASA's Airborne Multiangle SpectroPolarimetric Imager (AirMSPI)—near-UV (470 nm), visible (660 nm), and near-infrared (865 nm)—for our simulations which were performed over ocean surface using the successive order of scattering (SOS) algorithm coupled to a Lorenz-Mie aerosol optics model. The analysis of particle phase matrix elements indicates a close relationship between the angular dependencies of DoLP and associated phase matrix components at the shortest polarimetric band (470 nm). Using Jacobian analysis, we find that the radiance- and polarimetry-related reflectance factors of weakly light-absorbing aerosols, such as brown carbon, are more sensitive to changes in particle size and imaginary refractive index in comparison with those of black carbon, which is strongly light-absorbing. Our results suggest that the DoLP data could be used by future retrieval algorithms for reliably estimating microphysical properties of absorbing carbonaceous aerosols with imaginary refractive index less than 0.4.

refractive index↗