Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large system simulation”

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 739 records · Page 41

Climate Shifts Within Major Agricultural Seasons for +1.5 and +2.0 °C Worlds: HAPPI Projections and AgMIP Modeling Scenarios

This study compares climate changes in major agricultural regions and current agricultural seasons associated with global warming of +1.5 or +2.0 °C above pre-industrial conditions. It describes the generation of climate scenarios for agricultural modeling applications conducted as part of the Agricultural Model Intercomparison and Improvement Project (AgMIP) Coordinated Global and Regional Assessments. Climate scenarios from the Half a degree Additional warming, Projections, Prognosis and Impacts project (HAPPI) are largely consistent with transient scenarios extracted from RCP4.5 simulations of the Coupled Model Intercomparison Project phase 5 (CMIP5). Focusing on food and agricultural systems and top-producing breadbaskets in particular, we distinguish maize, rice, wheat, and soy season changes from global annual mean climate changes. Many agricultural regions warm at a rate that is faster than the global mean surface temperature (including oceans) but slower than the mean land surface temperature, leading to regional warming that exceeds 0.5 °C between the +1.5 and +2.0 °C Worlds. Agricultural growing seasons warm at a pace slightly behind the annual temperature trends in most regions, while precipitation increases slightly ahead of the annual rate. Rice cultivation regions show reduced warming as they are concentrated where monsoon rainfall is projected to intensify, although projections are influenced by Asian aerosol loading in climate mitigation scenarios. Compared to CMIP5, HAPPI slightly underestimates the CO2 concentration that corresponds to the +1.5 °C World but overestimates the CO2 concentration for the +2.0 °C World, which means that HAPPI scenarios may also lead to an overestimate in the beneficial effects of CO2 on crops in the +2.0 °C World. HAPPI enables detailed analysis of the shifting distribution of extreme growing season temperatures and precipitation, highlighting widespread increases in extreme heat seasons and heightened skewness toward hot seasons in the tropics. Shifts in the probability of extreme drought seasons generally tracked median precipitation changes; however, some regions skewed toward drought conditions even where median precipitation changes were small. Together, these findings highlight unique seasonal and agricultural region changes in the +1.5 °C and +2.0 °C worlds for adaptation planning in these climate stabilization targets.

Alexander C. Ruane↗

FUN3D Manual: 13.5

This manual describes the installation and execution of FUN3D (Fully-UNstructured three-dimensional CFD (Computational Fluid Dynamics) code) version 13.5, including optional dependent packages. FUN3D is a suite of computational fluid dynamics simulation and design tools that uses mixed-element unstructured grids in a large number of formats, including structured multiblock and overset grid systems. A discretely-exact adjoint solver enables efficient gradient-based design and grid adaptation to reduce estimated discretization error. FUN3D is available with and without a reacting, real-gas capability. This generic gas option is available only for those persons that qualify for its beta release status.

Biedron, Robert T.↗

FUN3D Manual: 13.6

This manual describes the installation and execution of FUN3D version 13.6, including optional dependent packages. FUN3D is a suite of computational fluid dynamics simulation and design tools that uses mixed-element unstructured grids in a large number of formats, including structured multiblock and overset grid systems. A discretely-exact adjoint solver enables efficient gradient-based design and grid adaptation to reduce estimated discretization error. FUN3D is available with and without a reacting, real-gas capability. This generic gas option is available only for those persons that qualify for its beta release status.

Robert T Biedron↗

Implementation of A New Microwave Scattering Database and A Forward Model for Active Microwave Sensors in CRTM

Radiative transfer models are extensively used for the assimilation of satellite observations into NWP models as well as retrieving geophysical products from satellite measurements. CRTM is a community model developed by NOAA JCSDA and widely used for different purposes requiring RT calculations. CRTM requires bulk scattering lookup tables in order to perform all-sky RT calculations. However, the current CRTM lookup tables for microwave frequencies were generated based on the Mie theory by assuming spherical frozen particles. The scattering lookup tables generated using the DDA technique has shown to largely improve the RT scattering calculations in the MW region. This presentation targets (i) the implementation and validation of a DDA database that was originally developed for the ARTS RT model into CRTM, and (ii) developing the CRTM active sensor module that takes advantage of the backscattering coefficients computed using the DDA method. The DDA database only provides single scattering properties of different habits, while CRTM requires bulk scattering properties. The CRTM cloud coefficients were previously generated based on the effective radius for representing the size of the particles. However, effective radius is neither measurable nor provided by the NWP models, thus need to be estimated from other geophysical variables such as water content. Therefore, in addition to calculating the CRTM bulk scattering properties from the DDA single scattering database, the CRTM was also largely modified to use cloud water content (kg.m-3), instead of effective radius, for performing the interpolation over size/mass of the particles. CRTM already requires water content as input, thus no extra variables are required for performing scattering calculations using the new ARTS DDA database. The CRTM scattering modules search for effective radius in cloud coefficient files and will use the cloud water content if the effective radius dimension is not found in the cloud coefficient files. Figure 1 shows the CRTM simulated brightness temperatures computed using different cloud coefficients versus ATMS observed values over Hurricane Irma on September 7, 2017 at 18:00 UTC. We used all the cloud water content values included in ERA5 with default CRTM/DDA habits for water, rain, snow, ice, hail, and graupel. ERA5 does not provide separate water content values for ice, hail, and graupel, thus the ice water content values were divided between ice, hail, and graupel clouds similar to what was explained in the previous section. In channels with a frequency lower than 90 GHz, emission from water and rain clouds can compensate for cloud scattering so that cloud contaminated Tbs are larger than corresponding clear sky Tbs. The DDA simulations for channels 1-7 largely perform better than the Mie simulations. The DDA simulations show a mix of small negative and positive simulated minus observed values, while the Mie results show large negative biases. The weighting functions for some of the ATMS temperature sounding channels (channels 9-15) peak mostly above the clouds, therefore the measured Tbs become less sensitive to clouds so that the results of both Mie and DDA become very similar. The Mie lookup tables generate excessive scattering for channel 16, but not enough scattering for the water vapor channels. In the specific case of Hurricane Maria, the DDA lookup tables do not generate enough scattering for channel 16, but the DDA results are much more consistent with observations for water vapor channels than for channel 16. It should be noted that the results may vary if we use other habits to represent snow, hail, and graupel in the DDA simulations. Although these results clearly show the advantage of the DDA database over the Mie dataset, different error sources such as error in the observations, displacement of clouds in the ERA5 reanalysis, and also lack of convective clouds or in general errors in the input atmospheric and cloud profiles contribute to the differences between the simulated and observed values. Aside from the improvements in the simulations, a major advantage of the new dataset is a large number of habits that can be used to tune the data assimilation systems to perform well in different weather conditions.

Isaac Moradi↗

Development of a Hybrid RANS/LES Method for Compressible Mixing Layer Simulations

A hybrid method has been developed for simulations of compressible turbulent mixing layers. Such mixing layers dominate the flows in exhaust systems of modem day aircraft and also those of hypersonic vehicles currently under development. The hybrid method uses a Reynolds-averaged Navier-Stokes (RANS) procedure to calculate wall bounded regions entering a mixing section, and a Large Eddy Simulation (LES) procedure to calculate the mixing dominated regions. A numerical technique was developed to enable the use of the hybrid RANS/LES method on stretched, non-Cartesian grids. The hybrid RANS/LES method is applied to a benchmark compressible mixing layer experiment. Preliminary two-dimensional calculations are used to investigate the effects of axial grid density and boundary conditions. Actual LES calculations, performed in three spatial directions, indicated an initial vortex shedding followed by rapid transition to turbulence, which is in agreement with experimental observations.

Georgiadis, Nicholas J.↗

Marshall Space Flight Center Materials and Processes Laboratory

Marshall?s Materials and Processes Laboratory has been a core capability for NASA for over fifty years. MSFC has a proven heritage and recognized expertise in materials and manufacturing that are essential to enable and sustain space exploration. Marshall provides a "systems-wise" capability for applied research, flight hardware development, and sustaining engineering. Our history of leadership and achievements in materials, manufacturing, and flight experiments includes Apollo, Skylab, Mir, Spacelab, Shuttle (Space Shuttle Main Engine, External Tank, Reusable Solid Rocket Motor, and Solid Rocket Booster), Hubble, Chandra, and the International Space Station. MSFC?s National Center for Advanced Manufacturing, NCAM, facilitates major M&P advanced manufacturing partnership activities with academia, industry and other local, state and federal government agencies. The Materials and Processes Laborato ry has principal competencies in metals, composites, ceramics, additive manufacturing, materials and process modeling and simulation, space environmental effects, non-destructive evaluation, and fracture and failure analysis provide products ranging from materials research in space to fully integrated solutions for large complex systems challenges. Marshall?s materials research, development and manufacturing capabilities assure that NASA and National missions have access to cutting-edge, cost-effective engineering design and production options that are frugal in using design margins and are verified as safe and reliable. These are all critical factors in both future mission success and affordability.

Tramel, Terri L.↗

Efficient massively parallel simulation of dynamic channel assignment schemes for wireless cellular communications

Fast, efficient parallel algorithms are presented for discrete event simulations of dynamic channel assignment schemes for wireless cellular communication networks. The driving events are call arrivals and departures, in continuous time, to cells geographically distributed across the service area. A dynamic channel assignment scheme decides which call arrivals to accept, and which channels to allocate to the accepted calls, attempting to minimize call blocking while ensuring co-channel interference is tolerably low. Specifically, the scheme ensures that the same channel is used concurrently at different cells only if the pairwise distances between those cells are sufficiently large. Much of the complexity of the system comes from ensuring this separation. The network is modeled as a system of interacting continuous time automata, each corresponding to a cell. To simulate the model, conservative methods are used; i.e., methods in which no errors occur in the course of the simulation and so no rollback or relaxation is needed. Implemented on a 16K processor MasPar MP-1, an elegant and simple technique provides speedups of about 15 times over an optimized serial simulation running on a high speed workstation. A drawback of this technique, typical of conservative methods, is that processor utilization is rather low. To overcome this, new methods were developed that exploit slackness in event dependencies over short intervals of time, thereby raising the utilization to above 50 percent and the speedup over the optimized serial code to about 120 times.

Greenberg, Albert G.↗

Multifractal characterizations of nonstationarity and intermittency in geophysical fields: Observed, retrieved, or simulated

Geophysical data rarely show any smoothness at any scale, and this often makes comparison with theoretical model output difficult. However, highly fluctuating signals and fractal structures are typical of open dissipative systems with nonlinear dynamics, the focus of most geophysical research. High levels of variability are excited over a large range of scales by the combined actions of external forcing and internal instability. At very small scales we expect geophysical fields to be smooth, but these are rarely resolved with available instrumentation or simulation tools; nondifferentiable and even discontinuous models are therefore in order. We need methods of statistically analyzing geophysical data, whether measured in situ, remotely sensed or even generated by a computer model, that are adapted to these characteristics. An important preliminary task is to define statistically stationary features in generally nonstationary signals. We first discuss a simple criterion for stationarity in finite data streams that exhibit power law energy spectra and then, guided by developments in turbulence studies, we advocate the use of two ways of analyzing the scale dependence of statistical information: singular measures and qth order structure functions. In nonstationary situations, the approach based on singular measures seeks power law behavior in integrals over all possible scales of a nonnegative stationary field derived from the data, leading to a characterization of the intermittency in this (gradient-related) field. In contrast, the approach based on structure functions uses the signal itself, seeking power laws for the statistical moments of absolute increments over arbitrarily large scales, leading to a characterization of the prevailing nonstationarity in both quantitative and qualitative terms. We explain graphically, step by step, both multifractal statistics which are largely complementary to each other. The geometrical manifestations of nonstationarity and intermittency, 'roughness' and 'sparseness', respectively, are illustrated and the associated analytical (differentiability and continuity) properties are discussed. As an example, the two techniques are applied to a series of recent measurements of liquid water distributions inside marine stratocumulus decks; these are found to be multifractal over scales ranging from approximately 60 m to approximately 60 km. Finally, we define the 'mean multifractal plane' and show it to be a simple yet comprehensive tool with many applications including data intercomparison, (dynamical or stochastic) model and retrieval validations.

Davis, Anthony↗

Mesoscale Simulations of CRYSTAL-FACE 23 July 2002 Case

A key objective of the Cirrus Regional Study of Tropical Anvils and Cirrus Layers - Florida Area Cirrus Experiment (CRYSTAL-FACE) is to understand the relationships between properties of tropical convective cloud systems and the lifecycle of the extended cirrus anvils they produce. We report here on a case study of 23 July 2002 where a line of land-based convective storms was generated between Lake Okeechobee and the Florida east coast as a result of complex interactions between lake and sea breeze fronts and outflow boundaries. A central goal of this study is to develop a description of convective input to the anvil system and to quantify the ongoing dynamical forcing of anvil processes by mesoscale and large-scale dynamics. This information is then used to force high-resolution cloud simulations with a model that explicitly resolves cloud microphysical processes (bin model) for study of cirrus anvil microphysical development.

Starr, David↗

JSC-Rocknest: A Large-Scale Mojave Mars Simulant (MMS) Based Soil Simulant for In-Situ Resource Utilization Water-Extraction Studies

The Johnson Space Center-Rocknest (JSC-RN) simulant was developed in response to a need by NASA's Advanced Exploration Systems (AES) In-Situ Resource Utilization (ISRU) project for a simulant to be used in component and system testing for water extraction from Mars regolith. JSC-RN was designed to be chemically and mineralogically similar to material from the aeolian sand shadow named Rocknest in Gale Crater, particularly the 1-3 wt.% low temperature (<450 ºC) water release as measured by the Sample Analysis at Mars (SAM) instrument on the Curiosity rover. Sodium perchlorate, goethite, pyrite, ferric sulfate, regular and high capacity granular ferric oxide, and forsterite were added to a Mojave Mars Simulant (MMS) base in order to match the mineralogy, evolved gases, and elemental chemistry of Rocknest. Mineral and rock components were sent to the United States Geological Survey (USGS) in Denver for mixing. Approximately 800 kg of JSC-RN was sent back to NASA in 5 gallon buckets, which were subsampled and characterized. All samples of the USGS-produced simulants had similar evolved gas profiles as a small prototype batch of JSC-RN made in JSC laboratories, with the exception of HCl, and were similar in terms of mineralogy and total chemistry. Also, all JSC-RN subsamples were homogenous and had similar mineralogy, total chemistry, and low-temperature evolved gas profiles as the Rocknest aeolian sand shadow examined with Curiosity’s instrument suite on Mars. In particular, the low temperature water releases were similar and the amount of water evolved from JSC-RN at <450 ºC was similar to the water content of Rocknest based on SAM water peak integrations. Overall, JSC-RN is ideally suited for ISRU studies of water extraction of global martian soil due to its excellent agreement with measured properties of martian soils and its proven feasibility for large-scale production.

Simulant↗

Recent Improvements in the Dart Model for Atmosphere, Topography, Large Landscape, Chlorophyll Fluorescence, Satellite Image Inversion

Physical models simulating the radiative budget (RB) and remote sensing (RS) observation of three-dimensional (3D) landscapes are critical to better understand human and natural components of the Earth system and further develop RS technology. DART is one of the most comprehensive 3D models of Earth-atmosphere optical radiative transfer (RT), from ultraviolet (UV) to thermal infrared (TIR). It simulates the optical signal of proximal, aerial and satellite imaging spectrometers and laser scanners, the 3D RB and solar induced chlorophyll fluorescence (SIF) signal, for any urban or natural landscape and any experimental or instrument configuration. It is freely available for research and teaching activities (https://dart.omp.eu). Here, five recent advances are presented. 1) Atmosphere RT. 2) RT in non repetitive topography. 3) Monte Carlo modelling for fast RS image simulation of large landscapes. 4) SIF modelling for vegetation simulated as facets and turbid cells. 5) RS image inversion for mapping the optical properties of urban material and the urban radiative budget.

DART↗

Survey Summary of AGATE Concepts Demonstration at Annual NATA Convention March 24-26, 1997: Original Survey Data - Volume 2

An AGATE Concepts Demonstration was conducted at the Annual National Air Transportation Association (NATA) Convention in 1997. Following, a 5-minute introductory briefing, an interactive simulation of a single-pilot, single-engine aircraft was conducted. The participant was able to take off, fly a brief enroute segment, fly a Global Positioning System (GPS) approach and landing, and repeat the approach and landing segment. The participant was provided an advanced 'highway-in-the-sky' presentation on both a simulated head-up display and on a large LCD head-down display to follow throughout the flight. A single-lever power control and display concept was also provided for control of the engine throughout the flight. A second head-down, multifunction display in the instrument panel provided a moving map display for navigation purposes and monitoring of the status of the aircraft's systems.

Source record↗

Survey Summary of AGATE Concepts Demonstration at Annual NATA Convention, March 24-26, 1997: Basic Report - Volume 1

An AGATE Concepts Demonstration was conducted at the annual National Air Transportation Association (NATA) Convention in 1997. Following a 5-minute introductory briefing, an interactive simulation of a single-pilot, single-engine aircraft was conducted. The participant was able to take off, fly a brief enroute segment, fly a Global Positioning System (GPS) approach and landing, and repeat the approach and landing segment. The participant was provided an advanced "highway-in-the-sky" presentation on both a simulated head-up display and on a large LCD head-down display to follow throughout the flight. A single-lever power control and display concept was also provided for control of the engine throughout the flight. A second head-down, multifunction display in the instrument panel provided a moving map display for navigation purposes and monitoring of the status of the aircraft's systems. An estimated 100 people observed or participated in the demonstration, and 68 surveys were collected. The pilot ratings of the participants ranged from student to Air Transport Rating with an average of 6089 hours total flight time. The overwhelming response was that technologies that simplify piloting tasks are enthusiastically welcomed by pilots of all experience levels. The increase in situation awareness and use of the head-up display were universally accepted and lauded as steps in the right direction.

Source record↗

An ultrahigh-accuracy body pointing system for the Large Space Telescope

The Large Space Telescope (LST) program is aimed at placing a three-meter diffraction-limited telescope in a 270-nm orbit to perform astronomical observations that are not possible with earth-based telescopes. A complex simulation model is described which was developed to determine whether the stringent pointing stability requirements could be met. The model (programmed on a hybrid computer) included detailed dynamic representation of control moment gyros (CMGs) and reaction wheels (RWs), including their noise characteristics; dynamic sensor representation (including noise); shockmounts for the CMG actuators; detailed representation of an image motion compensation system; and a detailed flexible body vehicle model. Stability and performance studies based on the simulation model showed that the body pointing system will meet LST requirements in the presence of CMG vibrational disturbances and sensor noise. The recommended system consists of three orthogonally mounted RWs for primary short-term control, and a cluster of CMG actuators for continuous RW desaturation and vehicle maneuvering.

Rybak, S. C.↗

Opportunities for Breakthroughs in Large-Scale Computational Simulation and Design

Opportunities for breakthroughs in the large-scale computational simulation and design of aerospace vehicles are presented. Computational fluid dynamics tools to be used within multidisciplinary analysis and design methods are emphasized. The opportunities stem from speedups and robustness improvements in the underlying unit operations associated with simulation (geometry modeling, grid generation, physical modeling, analysis, etc.). Further, an improved programming environment can synergistically integrate these unit operations to leverage the gains. The speedups result from reducing the problem setup time through geometry modeling and grid generation operations, and reducing the solution time through the operation counts associated with solving the discretized equations to a sufficient accuracy. The opportunities are addressed only at a general level here, but an extensive list of references containing further details is included. The opportunities discussed are being addressed through the Fast Adaptive Aerospace Tools (FAAST) element of the Advanced Systems Concept to Test (ASCoT) and the third Generation Reusable Launch Vehicles (RLV) projects at NASA Langley Research Center. The overall goal is to enable greater inroads into the design process with large-scale simulations.

Alexandrov, Natalia↗

The Role of Atmospheric Aerosol Concentration on Deep Convective Precipitation: Cloud-resolving Model Simulations

Aerosols and especially their effect on clouds are one of the key components of the climate system and the hydrological cycle [Ramanathan et al., 20011. Yet, the aerosol effect on clouds remains largely unknown and the processes involved not well understood. A recent report published by the National Academy of Science states "The greatest uncertainty about the aerosol climate forcing - indeed, the largest of all the uncertainties about global climate forcing - is probably the indirect effect of aerosols on clouds NRC [2001]." The aerosol effect on clouds is often categorized into the traditional "first indirect (i.e., Twomey)" effect on the cloud droplet sizes for a constant liquid water path and the "semi-direct" effect on cloud coverage. The aerosol effect on precipitation processes, also known as the second type of aerosol indirect effect, is even more complex, especially for mixed-phase convective clouds. ln this paper, a cloud-resolving model (CRM) with detailed spectral-bin microphysics was used to examine the effect of aerosols on three different deep convective cloud systems that developed in different geographic locations: South Florida, Oklahoma and the Central Pacific. In all three cases, rain reaches the ground earlier for the low CCN (clean) case. Rain suppression is also evident in all three cases with high CCN (dirty) case. However, this suppression only occurs during the first hour of the simulations. During the mature stages of the simulations, the effects of increasing aerosol concentration range from rain suppression in the Oklahoma case, to almost no effect in the Florida case, to rain enhancement in the Pacific case. These results show the complexity of aerosol interactions with convection.

Tao, W.-K.↗

The Role of Atmospheric Aerosol Concentration on Deep Convective Precipitation: Cloud-Resolving Model Simulations

Aerosols and especially their effect on clouds are one of the key components of the climate system and the hydrological cycle [Ramanathan et al., 2001]. Yet, the aerosol effect on clouds remains largely unknown and the processes involved not well understood. A recent report published by the National Academy of Science states "The greatest uncertainty about the aerosol climate forcing - indeed, the largest of all the uncertainties about global climate forcing - is probably the indirect effect of aerosols on clouds NRC [2001]." The aerosol effect on Clouds is often categorized into the traditional "first indirect (i.e., Twomey)" effect on the cloud droplet sizes for a constant liquid water path and the "semi-direct" effect on cloud coverage. The aerosol effect on precipitation processes, also known as the second type of aerosol indirect effect, is even more complex, especially for mixed-phase convective clouds. In this paper, a cloud-resolving model (CRM) with detailed spectral-bin microphysics was used to examine the effect of aerosols on three different deep convective cloud systems that developed in different geographic locations: South Florida, Oklahoma and the Central Pacific, In all three cases, rain reaches the ground earlier for the low CCN (clean) case. Rain suppression is also evident in all three cases with high CCN (dirty) case. However, this suppression only occurs during the first hour of the simulations. During the mature stages of the simulations, the effects of increasing aerosol concentration range from rain suppression in the Oklahoma case, to almost no effect in the Florida case, to rain enhancement in the Pacific case. These results show the complexity of aerosol interactions with convection. The model results suggest that evaporative cooling is a key process in determining whether high CCN reduces or enhances precipitation. Stronger evaporative cooling can produce a stronger cold pool and thus stronger low-level convergence through interactions with the low-level wind shear. Consequently, precipitation processes can be more vigorous. For example,, the evaporative cooling is more than two times stronger in the lower troposphere with high CCN for the Pacific case. Sensitivity tests also suggest that ice processes are crucial for suppressing precipitation in the Oklahoma case with high CCN.

Tao, Wei-Kuo↗

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↗