Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random fields”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Conspicuity of target lights: The influence of flash rate and brightness

The stimulus characteristics of lights that might aid a pilot to see and avoid, by alerting him to a potential threat were studied. The relative conspicuity of foveally equated, point-source, steady and flashing lights of several brightnesses, seen against a star background was examined. From the subject's viewpoint, these target lights could appear anywhere within a large (40 deg horizontal by 35 deg vertical) field of view. The lights appeared at random time intervals while the subject was periodically distracted by a simulated cockpit task. The results indicate that correct target detection increases and reaction time decreases with increased target intensity. Steady lights are missed more frequently and acquired more slowly than flashing lights, but no significant differences are found among the wide range of flash rates employed. The intensity of the light has a greater effect on both detection and reaction time to steady lights than to flashing lights. These results are compared with results of other researchers who used targets which appeared at fixed locations. The longest reaction times were recorded to lights which appeared either at the extremes or at the very center of the visual field.

Connors, M. M.↗
NASA NTRS2016

Electromagnetic Scattering by Fully Ordered and Quasi-Random Rigid Particulate Samples

In this paper we have analyzed circumstances under which a rigid particulate sample can behave optically as a true discrete random medium consisting of particles randomly moving relative to each other during measurement. To this end, we applied the numerically exact superposition T-matrix method to model far-field scattering characteristics of fully ordered and quasi-randomly arranged rigid multiparticle groups in fixed and random orientations. We have shown that, in and of itself, averaging optical observables over movements of a rigid sample as a whole is insufficient unless it is combined with a quasi-random arrangement of the constituent particles in the sample. Otherwise, certain scattering effects typical of discrete random media (including some manifestations of coherent backscattering) may not be accurately replicated.

Radiative transfer↗
NASA NTRS1973

Effect of prolonged bedrest and plus Gz acceleration on peripheral visual response time

Peripheral visual response time changes during +G sub z acceleration following fourteen days of bedrest are considered as well as what effect prolonged bedrest has upon this response. Eighteen test lights, placed 10 deg are apart along the horizontal meridian of the subject's field of view, were presented in a random sequence. The subject was instructed to press a button as soon as a light appeared. Response time testing occurred periodically during bedrest and continuously during centrifugation testing. The results indicate that: (1) mean response time is significantly longer to stimuli imaged in the far periphery than to stimuli imaged closer to the line of sight; (2) mean response time at each stimulus position tends to be longer at plateau g than during the preacceleration baseline period; (3) mean response time tends to lengthen as the g level is increased; (4) peripheral visual response time during +G sub x acceleration at 2, 3.2, and 3.8 g was not a reliable advanced indicator that blackout was going to occur; and (5) the subject's field of view collapsed rapidly just before blackout. Bedrest data showed that the distribution of response times to stimuli imaged across the subject's horizontal retinal meridian remained remarkably constant from day to day during both the bedrest and recovery periods.

Haines, R. F.↗
NASA NTRS1981

Correlation function studies for snow and ice

The random medium model is used to characterize snow and ice fields in the interpretation of active and passive microwave remote sensing data. A correlation function is used to describe the random permittivity fluctuations with the associated mean and variance and correlation lengths; and several samples are investigated to determine typical correlation functions for snow and ice. It is shown that correlation functions are extracted directly from appropriate ground truth data, and an exponential correlation function is observed for snow and ice with lengths corresponding to the actual size of ice particles or air bubbles. Thus, given that a medium has spatially stationary statistics and a small medium, the random medium model can interpret remote sensing data where theoretical parameters correspond to actual physical parameters of the terrain.

Vallese, F.↗
NASA NTRS1990

Comparison of GSC photometry in three southern fields

Entries in the Guide Star Catalog (GSC) within three southern declination fields surveyed by Ratnatunga in 1983 are studied. All three fields have high galactic latitude. About 98 percent of the 13,540 GSC entries within the total survey area of 60 deg sq are cross identified. Using over 9000 stars in the magnitude range J between 12.5 and 15.5, the random rms precision in photometric magnitude within these fields is estimated to be 0.1-0.2 mag, much lower than the total rms accuracy of 0.4 mag estimated in GSC. Some of the scanned images of the remaining unmatched entries were visually inspected to identify the cause of discrepancy.

Ratnatunga, Kavan U.↗
NASA NTRS1992

A unified perturbation expansion for surface scattering

Starting with the extinction theorem, a perturbation expansion which, to first and second orders, converges over a wider domain than the small perturbation expansion and the momentum transfer expansion is presented. It is shown that, in the appropriate limits, both of these theories, as well as the two-scale expansion, are recovered. There is no adjustable parameter, such as a spectral split, in the theory. This theory is applied to random rough surfaces and derive analytic expressions for the coherent field and the bistatic cross section. Finally, a numerical test of the theory against method of moments results for Gaussian random rough surfaces with a power law spectrum is given. These results show that the expansion is ramarkably accurate over a large range of surface heights and slopes for both horizontal and vertical polarization.

Rodriguez, Ernesto↗
NASA NTRS1999

Phase Diagrams of Electric-Fduced Aggregation in Conducting Colloids

Under the application of a sufficiently strong electric field, a suspension may undergo reversible phase transitions from a homogeneous random arrangement of particles into a variety of ordered aggregation patterns. The surprising fact about electric-field driven phase transitions is that the aggregation patterns, that are observed in very diverse systems of colloids, display a number of common structural features and modes of evolution thereby implying that a universal mechanism may exist to account for these phenomena. It is now generally believed that this mechanism emanates from the presence of the long-range anisotropic interactions between colloidal particles due to their polarization in an applied field. But, in spite of numerous applications of the electric-field-driven phenomena in biotechnology, separation, materials engineering, chemical analysis, etc. our understanding of these phenomena is far from complete. Thus, it is the purpose of the proposed research to develop a theory and then test experimentally, under normal- and low-gravity conditions, the accuracy of the theoretical predictions regarding the effect of the synergism of the interparticle electric and hydrodynamic interactions on the phase diagram of a suspension. The main results from our theoretical studies performed to-date enable one to trace how the variations of the electrical properties of the constituent materials influence the topology of the suspension phase diagram and then, by using an appropriate phase diagram, to evaluate how the electric-field-induced transformations will depend on the frequency and the strength of the applied field.

Khusid, B.↗
DOE OSTI2025

Karhunen–Loève deep learning method for surrogate modeling and approximate Bayesian parameter estimation

We evaluate the performance of the Karhunen-Loève Deep Neural Network (KL-DNN) framework for surrogate modeling and approximate Bayesian parameter estimation in partial differential equation models. In the surrogate model, the Karhunen-Loève (KL) expansions are used for the dimensionality reduction of the number of unknown parameters and variables, and a deep neural network is employed to relate the reduced space of parameters to that of the state variables. The KL-DNN surrogate model is used to formulate a maximum-a-posteriori-like least-squares problem, which is randomized to draw samples of the posterior distribution of the parameters. We test the proposed framework for a hypothetical unconfined aquifer via comparison with the forward MODFLOW and inverse PEST++ iterative ensemble smoother (IES) solutions as well as the state-of-the-art Fourier neural operator (FNO) and deep operator networks (DeepONets) operator learning surrogate models. Our results show that the KL-DNN surrogate model outperforms FNO and DeepONet for forward predictions. For solving inverse problems, the randomized algorithm provides the same or more accurate Bayesian predictions of the parameters than IES as evidenced by the higher log-predictive probability of both the estimated parameter field and the forecast hydraulic head. The posterior mean obtained from the randomized algorithm is closer to the reference parameter field than that obtained with FNO as the maximum a posteriori estimate.

Approximate Bayesian inference↗
NASA NTRS1999

Dynamics of Single Chains of Suspended Ferrofluid Particles

We present an experimental study of the dynamics of isolated chains made of super-paramagnetic particles under the influence of a magnetic field. The motivation of this work is to understand if the chain fluctuations exist and, if it does, how does the fluctuation affect chain aggregation. We find that single chains strongly fluctuate and that the characteristic frequency of their fluctuations is inversely proportional to the magnetic field strength. The higher the field the lower the characteristic frequency of the chain fluctuations. In the high magnetic field limit, chains behave like rigid rods without any internal motions. In this work, we used ferrofluid particles suspended in water. These particles do not have any intrinsic magnetization. Once a magnetic field is applied, a dipole moment is induced in each particle, proportional to the magnetic field. A dipolar magnetic interaction then occurs between particles. If dipole-dipole magnetic energy is higher than the thermal energy, the result is a structure change inside the dipolar fluid. The ratio of these two energies is expressed by a coupling constant lambda as: lambda = (pi(a(exp 3))(chi(exp 2))(mu(sub 0))(H(sub 0))(exp 2))/18kT Where a is the particle radius, mu(sub 0) is the vacuum magnetic permeability, H(sub 0) the applied magnetic field, k the Boltzmann constant and T the absolute temperature. If lambda > 1, magnetic particles form chains along the field direction. The lateral coalescence of several chains may form bigger aggregates especially if the particle volume fraction is high. While many studies and applications deal with the rheological properties and the structural changes of these dipolar fluids, this work focuses on the understanding of the chain dynamics. In order to probe the chain dynamics, we used dynamic light scattering (DLS) in self-beating mode as our experimental technique. The experimental geometry is such that the scattering plane is perpendicular to the magnetic field. Therefore, only motions in this plane are probed. A very dilute sample of a ferrofluid emulsion with a particle volume fraction of 10(exp -5) is used in this experiment. We chose such a low volume fraction to avoid multiple light scattering as well as lateral chain-chain aggregation. DLS measures the dynamic structure factor S(q,t) of the sample (q is the scattering wave vector, t is the time). In the absence of the magnetic field, identical particles of ferrofluid droplets are randomly distributed and S(q,t) reduces to exp(-q(exp 2)2D(sub 0)t). D(sub 0)=(kT/(6(pi)(eta)(a)) is the diffusion coefficient of Brownian particles (where Xi = (6(pi)(eta)(a)) is the Stokes frictional coefficient of a spherical particle in a fluid of viscosity eta). If interactions or polydispersity can not be ignored, an effective diffusion coefficient is introduced. Formally, D(sub eff) is defined as: D(sub eff) = - q(exp -2) partial derivative of (ln(S(q,t)) with respect to time, as t goes to 0. D(sub eff) reduces to D(sub 0) if no interactions and only a few particles size are present. Therefore, we can use DLS to measure particle size. The particle radius was found to be a=0.23 mu m with 7% of polydispersity. In this case, if we vary the scattering angle theta (and so q) we do not have any change in the measured diffusion coefficient: it is q-independent. When a magnetic field is applied, particles aggregate into chains if lambda > 1. We first studied the kinetics of the chain formation when lambda = 406. At a fixed scattering angle, we measured diffusion coefficient D(sub eff) as a function of time. Experimentally, we find that D(sub eff) decreases monotonously with time. Physically, this means that chains are becoming longer and longer. Since we are only sensitive to motions in the scattering plane and since chains have their main axis perpendicular to this plane, the measured diffusion coefficient is the trans-verse diffusion coefficient. We can relate D(sub eff) to the mean number of particles per chain N(t) at a given time and to the diffusion coefficient of an isolated particle D(sub 0) as D(sub eff)=f(N(t))D(sub 0). Since f(N) is known from other recent work, N can be expressed as a function of the time. We found a square root dependency: N(t) proportional to the square root of t. As expected for very low volume fraction, this behavior is characteristic of a diffusion-limited aggregation as suggested by several authors and by our previous work. In this study, we focus on the dependence of the effective diffusion coefficient on the scattering angle and the magnetic field strength. After the magnetic field is applied (lambda = 406) for a long time, typically 6 hours, kinetics of chain formation becomes very slow. Chain size does not vary much over the next hour period. Thus, we can perform different interesting experiments. First, at a fixed magnetic field, we measure the effective diffusion coefficient as a function of the scattering angle (from 5 to 130 deg). Our results show that the measured diffusion coefficient increases linearly with the scattering angle: D(sub eff) proportional to q. If we do the same experiment for different lambda values, D(sub eff) depends on lambda as D(sub eff) proportional to lambda(exp -1/2). We also find for different lambda values that the same asymptotic D(sub eff) value is obtained when q approaches zero. The angle dependency of D(sub eff) suggests that an additional motion exists besides chain drifting. Chain size is constant during experiment, which was verified by measuring the same diffusion coefficient at the beginning and at the end of the angle switching. If chains are rigid, D(sub eff) is independent of q. Therefore, we found that D(sub eff) not only measures the motion of the entire chain but also its internal fluctuations. These internal motions are the fluctuations of the particles in the chain. To understand the q dependency of D(sub eff), let us look at the probing length used. In our study, the characteristic length scale probed is l=2pi/q which is in the range of 0.9<l/a<20. When l is much larger than the particle radius (l/a going to infinity) we are mainly sensitive to the center of mass diffusion of the chain. D(sub eff) is then the diffusion of the entire chain and depends only on N but not on q and lambda. The value D(sub eff) allows us to obtain the number of particles per chain N. In the opposite limit (i.e. l/a < 1), we are sensitive to motions on the size of individual particles. The main contributions to the measured diffusion coefficient come from internal motions of the chains (i.e. particles' fluctuations inside the chain). We investigated also the effect of the chain size on D(sub eff). For the same value of the magnetic interaction lambda and different chain sizes, we found that the slope of D(sub eff) versus q decreases when N increases as a power law D(sub eff) proportional to N(exp -0.7). This is an important behavior because this means that the longer the chain the slower its internal fluctuations besides the whole chain motion.

Cutillas, S.↗
DOE OSTI2024

A new ductile, tougher resin for impregnation of superconducting magnets

A major remaining challenge for Nb 3 Sn high field magnets is their training due to random temperature variations in the coils. The main objective of our research is to reduce or eliminate it by finding novel impregnation materials in replacement of the epoxies currently used. An organic olefin-based thermosetting dicyclopentadiene resin, C 10 H 12 , commercially available in Japan as TELENE® by RIMTEC, was used to impregnate a short Nb 3 Sn undulator coil developed by ANL and FNAL. This magnet reached short sample limit after only two quenches, compared with ∼100 when CTD-101K® was used. Ductility, i.e. the ability to accept large strains, and toughness were identified as key properties to achieve these results. In addition, we have been investigating whether mixing TELENE with high heat capacity ceramic powders such as Gd 2 O 3 , Gd 2 O 2 S, and HoCu 2 , increases the specific heat (C p ) of impregnated Nb3Sn superconducting magnets. The viscosity, heat capacity, thermal conductivity, and other physical properties of TELENE with high-C p powder fillers were measured in this study as a function of temperature and magnetic field. The TELENE-87 wt%Gd 2 O 2 S had a peak in C p between 4.3 K and 5.3 K at fields between 0 and 8 T. We have also investigated the effect on the mechanical properties of pure and mixed TELENE under 10 MGy of gamma ray irradiation at the Takasaki Advanced Radiation Research Institute in Takasaki, Japan. TELENE-87 wt%Gd 2 O 2 S exhibited exceptional radiation resistance. Impregnating an undulator coil with TELENE mixed with Gd 2 O 2 S powder will verify whether the coils' thermal stability further improves, or whether its low diffusivity will require engineering the material with high-thermal conductivity components. Short magnet training will lead to better magnet reliability, lower magnet margins, lower risk and substantial saving in accelerators' commissioning costs. Part of this study is supported by the U.S.-Japan Science and Technology Cooperation Program in high energy physics operated by MEXT in Japan and DOE in the U.S.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

On the Relationship Between Shallow Cumulus Cloud Field Properties and Surface Solar Irradiance

Abstract Shallow cumulus clouds exhibit highly three‐dimensional (3‐D) spatial structure leading to complex variability in the surface solar irradiance (SSI) beneath. This variability is captured by the typically bimodal shape of the SSI probability density function (PDF). Using large eddy simulation to generate well‐resolved cloud fields and Monte Carlo 3‐D radiative transfer to reproduce realistic associated SSI PDFs, we seek direct relationships between the cloud field properties and the SSI PDF shape. Applying both random forest and artificial neural network algorithms, we find variations in the two modes of the SSI PDF are well predicted by just a handful of cloud field properties. The two algorithms utilize cloud properties similarly, with indistinguishable performance despite their different architectures. These results offer a marked improvement in realism relative to one‐dimensional radiative transfer while bypassing computationally expensive 3‐D radiative transfer, with immediate application to renewable energy assessments, and potential for several other geophysical applications.

58 GEOSCIENCES↗
NASA NTRS2013

Comparison of Precision of Biomass Estimates in Regional Field Sample Surveys and Airborne LiDAR-Assisted Surveys in Hedmark County, Norway

Airborne scanning LiDAR (Light Detection and Ranging) has emerged as a promising tool to provide auxiliary data for sample surveys aiming at estimation of above-ground tree biomass (AGB), with potential applications in REDD forest monitoring. For larger geographical regions such as counties, states or nations, it is not feasible to collect airborne LiDAR data continuously ("wall-to-wall") over the entire area of interest. Two-stage cluster survey designs have therefore been demonstrated by which LiDAR data are collected along selected individual flight-lines treated as clusters and with ground plots sampled along these LiDAR swaths. Recently, analytical AGB estimators and associated variance estimators that quantify the sampling variability have been proposed. Empirical studies employing these estimators have shown a seemingly equal or even larger uncertainty of the AGB estimates obtained with extensive use of LiDAR data to support the estimation as compared to pure field-based estimates employing estimators appropriate under simple random sampling (SRS). However, comparison of uncertainty estimates under SRS and sophisticated two-stage designs is complicated by large differences in the designs and assumptions. In this study, probability-based principles to estimation and inference were followed. We assumed designs of a field sample and a LiDAR-assisted survey of Hedmark County (HC) (27,390 km2), Norway, considered to be more comparable than those assumed in previous studies. The field sample consisted of 659 systematically distributed National Forest Inventory (NFI) plots and the airborne scanning LiDAR data were collected along 53 parallel flight-lines flown over the NFI plots. We compared AGB estimates based on the field survey only assuming SRS against corresponding estimates assuming two-phase (double) sampling with LiDAR and employing model-assisted estimators. We also compared AGB estimates based on the field survey only assuming two-stage sampling (the NFI plots being grouped in clusters) against corresponding estimates assuming two-stage sampling with the LiDAR and employing model-assisted estimators. For each of the two comparisons, the standard errors of the AGB estimates were consistently lower for the LiDAR-assisted designs. The overall reduction of the standard errors in the LiDAR-assisted estimation was around 40-60% compared to the pure field survey. We conclude that the previously proposed two-stage model-assisted estimators are inappropriate for surveys with unequal lengths of the LiDAR flight-lines and new estimators are needed. Some options for design of LiDAR-assisted sample surveys under REDD are also discussed, which capitalize on the flexibility offered when the field survey is designed as an integrated part of the overall survey design as opposed to previous LiDAR-assisted sample surveys in the boreal and temperate zones which have been restricted by the current design of an existing NFI.

Comparison↗
DOE OSTI2022

Pressure Drop Correlation Improvement for the Near-Wall Region of Pebble-Bed Reactors

Packed beds play an important role in several engineering fields, with their applications in nuclear energy being driven by the development of next-generation reactors utilizing pebble fuel. The random nature of a packed pebble bed creates a flow field that is complex and difficult to predict. Porous media models are an attractive option for modeling pebble-bed reactors (PBRs), as they provide intermediate fidelity results and are computationally efficient. Porous media models, however, rely on the use of correlations to estimate the effect of complicated flow features on the pressure drop and heat transfer in the system. Existing correlations were developed to predict the average behavior of the bed, but they are inaccurate in the near-wall region where the presence of the wall affects the pebble packing. This work aims to investigate the accuracy of a porous media model using the Kerntechnischer Ausschuss (KTA) correlation, the most common pressure drop correlation for PBRs compared to the high-fidelity large eddy simulation (LES). A bed of 1568 pebbles is investigated at Reynolds numbers from 625 to 10 000. The bed is divided into five concentric subdomains to compare the average velocity, friction losses, and form losses between the porous media and LES codes. The comparison between the LES simulation and the KTA correlation revealed that the KTA correlation largely underpredicts the form losses in the near-wall region, leading to an overprediction of the velocity near the wall by nearly 30%. An investigation of the form losses across the range of Reynolds numbers in the LES results provided additional insight into how the KTA correlation may be improved to better predict these spatial effects in a pebble bed. These data suggest that the form coefficient near the wall must be increased by 48% while decreasing the form coefficient of the inner bulk region of the bed by 15%. The implementation of these improvements to the KTA correlation in a porous media model produced a radial velocity profile that saw significantly improved agreement with the LES results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1987

Weak solar fields and their connection to the solar cycle

Videomagnetograph observations of the weak solar magnetic fields are discussed. Intranetwork (IN) elements are found to emerge in bipolar form, but to flow in a random pattern rather than toward the network elements. Merging of IN fields is found to be more important than the ephemeral regions as a source of new network elements. Arguments for and against the hypothesis that all solar fields are intrinsically strong are considered. It is noted that reconnection and merging of magnetic fields take place continually under the conditions presently examined.

Zirin, Harold↗
NASA NTRS1986

Radar scene generation for tactical decision aids

The Mueller matrix and polarization covariance matrix for polarimetric radar systems was studied. The clutter was modeled by a layer of random permittivity, described by a three-dimensional correlation function, with variance, and horizontal and vertical correlation lengths. This model was applied, using the wave theory with Born approximations carried to the second order, to find the backscattering elements of the polarimetric matrices. Theoretical predictions are matched with experimental data for vegetation fields. The strong fluctuation theory was used to derive the backscattering cross sections. A two-layer anisotropic random medium model was developed for the active and passive microwave remote sensing of ice fields. A three-layer random medium model was adopted to study the volume scattering effects for the active and passive microwave remote sensing of snow-covered ice fields.

Kong, J. A.↗
NASA NTRS1969

Solar protons E greater than 100 Mev incident over Antarctica during January- February 1967.

Commencing at 0825 +a -1 UT on January 28, 1967, a large and prolonged increase in the intensity of penetrating charged particles was observed by balloon-borne instruments floating over Byrd Station, Antarctica. (80°S, 120°W). A peak intensity of approximately 50 protons per cm 2 -sec-steradian with E> 100 MeV occurred at about 1230 UT on the 28th. The event was under observation almost continuously over a period of about 100 hours until the intensity decayed below cosmic-ray background on February 1. The initial decay was rapid but, some 40 hours after onset, went over into a slow exponential decay characterized by a 20 hour time-constant. The decay phase of an additional, though considerably less intense, event was observed on February 3 and 4. Presumably both events had their origins in major disturbances on the far side of the sun since nether event has been definitely linked to any feature which existed on the visible disk within an appropriate time interval. Results pertaining to the time-intensity profile and to the energy spectrum for protons E> 100 MeV are presented for the January 29 event. Comparison of the balloon results with neutron-monitor and satellite measurements and with models of interplanetary diffusion has led to some conclusions regarding the role of small-angle scattering by irregularities and by the random walk of magnetic lines of force relative to the mean interplanetary field within the orbit of earth.

Energy spectrum↗
NASA NTRS1992

The kinematics of intermediate-redshift gaseous galaxy halos

High-resolution spectroscopy of the absorption lines associated with four intermediate-redshift galaxies are presented. These observations are combined with existing spectroscopy of the redshifted 21-cm absorption lines toward 0235 + 164 and 1229 - 021. Evidence is presented to the effect that the lines of sight through the absorbers toward 0235 + 164 and 1229 - 021 penetrate material for which the velocity fields are dominated by systematic rather than by random motions, and the line profiles are found to be consistent with what is expected for rotating ensembles of clouds. It is suggested that for these two absorbers the absorption arises in individual clouds within single galaxy halos rather than, for example, in collections of coalescing subgalactic fragments, for which the expected velocity fields are more turbulent, or as a result of galaxy-galaxy clustering. It is concluded that lines of sight that impact absorbing regions at small impact parameters intercept many clouds, sample large velocity intervals, and are likely to show redshifted 21-cm absorption, whereas lines of sight that impact absorbing regions at large impact parameters intercept few clouds, sample small velocity intervals, and are unlikely to show redshifted 21-cm absorption.

Lanzetta, Kenneth M.↗
NASA NTRS2016

First-Principles Modeling Of Electromagnetic Scattering By Discrete and Discretely Heterogeneous Random Media

A discrete random medium is an object in the form of a finite volume of a vacuum or a homogeneous material medium filled with quasi-randomly and quasi-uniformly distributed discrete macroscopic impurities called small particles. Such objects are ubiquitous in natural and artificial environments. They are often characterized by analyzing theoretically the results of laboratory, in situ, or remote-sensing measurements of the scattering of light and other electromagnetic radiation. Electromagnetic scattering and absorption by particles can also affect the energy budget of a discrete random medium and hence various ambient physical and chemical processes. In either case electromagnetic scattering must be modeled in terms of appropriate optical observables, i.e., quadratic or bilinear forms in the field that quantify the reading of a relevant optical instrument or the electromagnetic energy budget. It is generally believed that time-harmonic Maxwell's equations can accurately describe elastic electromagnetic scattering by macroscopic particulate media that change in time much more slowly than the incident electromagnetic field. However, direct solutions of these equations for discrete random media had been impracticable until quite recently. This has led to a widespread use of various phenomenological approaches in situations when their very applicability can be questioned. Recently, however, a new branch of physical optics has emerged wherein electromagnetic scattering by discrete and discretely heterogeneous random media is modeled directly by using analytical or numerically exact computer solutions of the Maxwell equations. Therefore, the main objective of this Report is to formulate the general theoretical framework of electromagnetic scattering by discrete random media rooted in the Maxwell- Lorentz electromagnetics and discuss its immediate analytical and numerical consequences. Starting from the microscopic Maxwell-Lorentz equations, we trace the development of the first principles formalism enabling accurate calculations of monochromatic and quasi-monochromatic scattering by static and randomly varying multiparticle groups. We illustrate how this general framework can be coupled with state-of-the-art computer solvers of the Maxwell equations and applied to direct modeling of electromagnetic scattering by representative random multi-particle groups with arbitrary packing densities. This first-principles modeling yields general physical insights unavailable with phenomenological approaches. We discuss how the first-order-scattering approximation, the radiative transfer theory, and the theory of weak localization of electromagnetic waves can be derived as immediate corollaries of the Maxwell equations for very specific and well-defined kinds of particulate medium. These recent developments confirm the mesoscopic origin of the radiative transfer, weak localization, and effective-medium regimes and help evaluate the numerical accuracy of widely used approximate modeling methodologies.

Michael I Mishchenko↗