Search NASA⌕ Search

SEARCH · Search NASA

Results for “wave modeling”

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 541 records · Page 30

Chondrites and the Protoplanetary Disk, Part 3

Contents include the following: Ca-, Al-Rich Inclusions and Ameoboid Olivine Aggregates: What We Know and Don t Know About Their Origin. Aluminium-26 and Oxygen Isotopic Distributions of Ca-Al-rich Inclusions from Acfer 214 CH Chondrite. The Trapping Efficiency of Helium in Fullerene and Its Implicatiion to the Planetary Science. Constraints on the Origin of Chondritic Components from Oxygen Isotopic Compositions. Role of Planetary Impacts in Thermal Processing of Chondrite Materials. Formation of the Melilite Mantle of the Type B1 CAIs: Flash Heating or Transport? The Iodine-Xenon System in Outer and Inner Portions of Chondrules from the Unnamed Antarctic LL3 Chondrite. Nucleosynthesis of Short-lived Radioactivities in Massive Stars. The Two-Fluid Analysis of the Kelvin-Helmholtz Instability in the Dust Layer of a Protoplanetary Disk: A Possible Path to the Planetesimal Formation Through the Gravitational Instability. Shock-Wave Heating Model for Chonodrule Formation: Heating Rate and Cooling Rate Constraints. Glycine Amide Hydrolysis with Water and OH Radical: A Comparative DFT Study. Micron-sized Sample Preparation for AFM and SEM. AFM, FE-SEM and Optical Imaging of a Shocked L/LL Chondrite: Implications for Martensite Formation and Wave Propagation. Infrared Spectroscopy of Chondrites and Their Components: A Link Between Meteoritics and Astronomy? Mid-Infrared Spectroscopy of CAI and Their Mineral Components. The Origin of Iron Isotope Fractionation in Chondrules, CAIs and Matrix from Allende (CV3) and Chainpur (LL3) Chondrites. Protoplanetary Disk Evolution: Early Results from Spitzer. Kinetics of Evaporation-Condensation in a Melt-Solid System and Its Role on the Chemical Composition and Evolution of Chondrules. Oxygen Isotope Exchange Recorded Within Anorthite Single Crystal in Vigarano CAI: Evidence for Remelting by High Temperature Process in the Solar Nebula. Chondrule Forming Shock Waves in Solar Nebula by X-Ray Flares. Organic Globules with Anormalous Nitrogen Isotopic Compositions in the Tagish Lake Meteorite: Products of Primitive Organic Reactions. Yet Another Chondrule Formation Scenario. CAIs are Not Supernova Condensates. Microcrystals and Amorphous Material in Comets and Primitive Meteorites: Keys to Understanding Processes in the Early Solar System. A Nearby Supernova Injected Short-lived Radionuclides into Our Protoplanetary Disk. REE+Y Systematics in CC and UOC Chondrules. Meteoritic Constraints on Temperatures, Pressures, Cooling Rates, Chemical Compositions, and Modes of Condensation in the Solar Nebula. The I-Xe Record of Long Equilibration in Chondrules from the Unnamed Antarctic Meteorite L3/LL3. Early Stellar Evolution.

Source record↗

Mean Flow Augmented Acoustics in Rocket Systems

Oscillatory motion in solid rocket motors and liquid engines has long been a subject of concern. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process and gas dynamics. The customary approach to modeling acoustic waves inside a rocket chamber is to apply the classical inhomogeneous wave equation to the combustion gas. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while the acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The converging section of a rocket nozzle, where gradients in pressure, density, and velocity become large, is a notable region where this approach is not applicable. The expulsion of unsteady energy through the nozzle of a rocket is identified as the predominate source of acoustic damping for most rocket systems. An accurate model of the acoustic behavior within this region where acoustic modes are influenced by the presence of a steady mean flow is required for reliable stability predictions. Recently, an approach to address nozzle damping with mean flow effects was implemented by French [1]. This new approach extends the work originated by Sigman and Zinn [2] by solving the acoustic velocity potential equation (AVPE) formulated by perturbing the Euler equations [3]. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, (del squared)(psi) − (lambda/c)(exp 2)(psi) − M(dot)[M(dot)(del)(del(psi))] − 2(lambda(M/c) + (M(dot)del(M))(dot)del(psi)−2(lambda)(psi)[M(dot)del(1/c)]=0 (1) with M as the Mach vector, c as the speed of sound, and lambda as the complex eigenvalue. French apply the finite volume method to solve the steady flow field within the combustion chamber and nozzle with inviscid walls. The complex eigenvalues and eigenvector are determined with the use of the ARPACK eigensolver. The present study employs the COMSOL Multphysics framework to solve the coupled eigenvalue problem using the finite element approach. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluated the gas flow inside of a solid rocket motor using St. Robert's law modeling solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are used by the coefficient form of the mathematics module to determine the eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. Pertinent results from these analyses are the complex valued eigenvalue and eigenvectors. Comparisons are made to the French results to evaluate the modeling approach. A comparison of the French results with that of the present analysis is displayed in figures 1 and 2, respectively. The graphic shows the first tangential eigenvector's real (a) and imaginary (b) values.

Fischbach, Sean R.↗

Wave Setup at the Minamitorishima Tide Gauge

A tide gauge situated on the western shore of Minamitorishima, a small atoll in the western Pacific, has measured sea levels significantly impacted by wave setup. Evidence for this rests with (a) strong correlations between sea level anomalies and high swell and (b) sea-level differences with satellite altimetry that display near-linear dependence on offshore swell heights (regression coefficient 30 cm per meter of swell). Setup is primarily induced by swell from the northwest. We develop models of wave setup which lead to corrected sea levels better reflecting the surrounding ocean, and thus more readily useful to studies of regional sea level and practical applications such as altimetry calibration. The wave setup is also evidently affecting measurements of the tide, with suppressed tide amplitudes in winter when swell is generally largest. In February 2020 the tide gauge was relocated to the southern shore of the island, and the wave setup is now markedly reduced.

Wave setup↗

Modeling a Large Submillimeter-Wave Observatory

The 25 meter aperture Cornell Caltech Atacama Telescope (CCAT) will provide an enormous increase in sensitivity in the submillimeter bands compared to existing observatories, provided it can establish and maintain excellent image quality. To accomplish this at a very low cost, it is necessary to conduct accurate engineering trades, including the most effective segment and wavefront sensing and control approach, to determine the best method for continuously maintaining wavefront quality in the operational environment. We describe an integrated structural/optical/controls model that provides accurate performance prediction. We also detail the analysis methods used to quantify critical design trades.

Segment Sensing↗

Seismically induced modification of lunar surface features

The formation of large impact craters and basins produced not only large volumes of ejecta but also catastrophic seismic waves. Theoretical models, extrapolation of terrestrial explosion data, and extrapolation of lunar impact data suggest surface displacements of 1-10 m at four crater radii from a Copernicus-size impact (4 x 10 to the 30th power ergs). Independent estimates of impact-generated stresses indicate that regions within 4 crater radii of an impact will receive shock waves with pressures exceeding 1 kbar. Greater areas also may receive shock effects as elastic waves following deep ray paths transform into shock waves near the surface. Additional seismic energy will be generated by impacts from secondary ejecta. The resulting large and long-lasting surface movement may contribute to the subdued appearance of the continuous ejecta blanket and secondary craters around large impact craters and to the degradation of small craters outside these zones.

Schultz, P. H.↗

Models of the Cartwheel ring galaxy: Spokes and starbursts

Recent observations of this famous ring galaxy, including optical and near-infrared CCD surface photometry, and VLA radio continuum and 21 cm line mapping (Higdon 1992b, in prep.), have inspired a renewed modeling effort. Toomre's (1978, in The Large-scale Structure of the Universe, eds. Longair and Einasto) series of restricted three-body simulations demonstrated how the multiple rings could be produced in a nearly head-on galaxy collision. New models with a halo-dominated potential based on the 21 cm rotation curve are able to reproduce such details as the spacing between rings, ring widths, offset of the nucleus, and several kinematical features, thus providing strong support for the collisional theory. The new observations have shown there are little or no old stars in Cartwheel; it may consist almost entirely of gas and stars produced as a result of compression in the ring wave. To model this process Smooth Particle Hydrodynamics (SPH) simulations of the Cartwheel disk have been performed. Fixed gravitational potentials were used to represent the Cartwheel and a roughly 30 percent mass collision partner. The interaction dynamics was treated as in the usual restricted three-body approximation, and the effects of local self-gravity between disk particles were calculated. We are particularly interested in testing the theory that enhanced star formation in waves is the result of gravitational instability in the compressed region (see e.g. Kennicutt 1989, ApJ 344, 685). The gas surface density in a number of simulations was initialized to a value slightly below the threshold for local gravitational instability throughout most of the disk. The first ring wave produces relatively modest compressions (a factor of order a few), triggering instability in a narrow range of wavelengths. Self-gravity in the disk is calculated over a comparable range of scales. Simulations were run with isothermal, adiabatic, and adiabatic with radiative cooling characterized by a relatively short timescale. The isothermal approximation is good except in the vicinity of the strong second (inner) ring, and several snapshots from one case are shown in the figure below. Flocculent spiral segments are present before the collision, and these are compressed into dense knots in the ring wave. These knots are likely to be sites of vigorous star formation. In the strong rarefaction behind the outer ring most of the knots are radially stretched and sheared, giving rise to spoke-like features. A few dense knots are evidently very tightly bound, because they retain their coherence and are stretched relatively little through the rarefaction. This is in accord with evidence for continuing star formation in some spokes (Marcum, Appleton and Higdon 1992). The number and spacing of spokes is a direct function of the scale of the gravitational instability in the disk. Thus, the gravitational instability theory, together with the hypothesis that massive stars are only formed in dense knots of gas, can account for most of the distinct morphology of the Cartwheel.

Struck-Marcell, Curtis↗

Improved passive microwave sounding of the atmosphere

The effort this year focused primarily on 118-GHz transmittance experiments. The data analyzed here was collected with the Microwave Temperature Sounder (MTS) radiometer package during the CAMEX deployment of 1993 with the aim of validating current models of atmospheric microwave absorption in the O2 bands near 54 and 118 GHz. Particular attention has been paid to data collected during four flights when the MTS scanned zenith while profiles of downwelling radiances were collected through ascents and descents. These radiances, in conjunction with radiosonde temperature data, permit the retrieval of band-averaged absorption profiles for each channel. The Millimeter-wave Propagation Model (MPM92) provides theoretical expressions for the absorption of microwaves by oxygen and water vapor and accounts for the interference of pressure-broadened spectral lines'. This model is a good fit to laboratory measurements at temperatures ranging from 279-327 K, but it has been suggested that extrapolation to the conditions of the atmospheric tropopause may result in underestimation of absorption by as much as 15 percent. Preliminary results of the analysis of MTS data appear to be in general agreement with the predictions of the MPM model to within the accuracy of the measurements, which through the coldest parts of the atmosphere ranges from less than plus or minus 5 percent in the most opaque channels to greater than plus or minus 10 percent in the most transparent channels. At those altitudes where each channel is most sensitive to changes in absorption, there is some indication that the modeled absorption may be biased low relative to the observations. Accurate instrument calibration provided challenges, particularly when observed radiances were as much as 260 K below the temperatures of the cold calibration load.

Staelin, D. H.↗

The propagation of spark-produced N waves through turbulence

A model experiment was designed and built to simulate the propagation of sonic booms through atmospheric turbulence. The setup of the model experiment is described briefly. Measurements of the N waves after they propagated across the turbulent velocity field reveal the same waveform distortion and change in rise time as for sonic booms. The data from the model experiment is used to test sonic boom models. Some models yield predictions for the waveform distortion, while others give estimates of the rise time of the sonic booms. A new theoretical model for the propagation of plane N waves through a turbulent medium is described.

Lipkens, Bart↗

The Role of Gravity Waves in Generating Equatorial Oscillations in Modulating Atmospheric Tides

We discuss a Numerical Spectral Mode (NSM) that extends from the ground up into the thermosphere and incorporates Hines' Doppler spread parameterization (DSP) for small scale gravity waves (GW). This model is applied to describe the seasonal variations in the mean zonal circulation, the semi-annual and quasi-biennial oscillations (SAO and QBO), as well as the tides and planetary waves in the middle atmosphere. Initial results showed that this model can reproduce the salient features observed, including the QBO extending into the upper mesosphere inferred from UARS measurements. The model has now been extended to simulate also: (a) the zonal circulation of the lower stratosphere and upper troposphere, and (b) the upwelling at equatorial latitudes associated with the Brewer Dobsen circulation that affects the dynamics significantly as pointed out by Dunkerton. Upward vertical winds increase the period of the QBO observed from the ground. To compensate for that, one needs to increase in the model the eddy diffusivity and the GW momentum flux, bringing the latter closer to values recommended in the DSP. This development is conducive to extending the QBO and SAO to higher latitudes through global scale momentum redistribution. Multi-year interannual oscillations are generated through wave filtering by the solar driven annual oscillation in the zonal circulation. In a 3D version of the model, wave momentum is absorbed and dissipated by tides and planetary waves. A somewhat larger GW source (well within the DSP range) is then required to generate realistic QBO and SAO amplitudes. Since GW momentum is deposited in the altitude regime of increasing winds, the amplitude of the diurnal tide is amplified and its vertical wavelength is reduced at altitudes between 70 and 120 km. Wave filtering by the mean zonal circulation causes the GW flux to peak during equinox, and this produces a large semi-annual variation in the tide that has been observed on UARS. Without the diurnal tide, the semidiurnal tide would also be modulated in this way. But the diurnal tide filters out the GW preferentially during equinox, so that the semidiurnal tide tends to peak during solstice. Under the influence of GW, the tides are modulated significantly by planetary waves that are generated preferentially during solstice in part due to baroclinic instability.

Mayr, H. G.↗

On truncated expansion self-tuning regulation of a large AMCD

The one-dimensional wave equation models out of plane displacement of a flexible ring and with appropriate boundary conditions the solution is one of uncoupled modes. Identification and control of the behavior of these modes by means of a self-tuning regulator requires precise measurements of the ring shape. Aliasing of modes occurs from a finite number of sensors along the ring causing interference in the modal identification. Nonuniform sensors spacing allows for reduced aliasing effects but introduces spreading of the modes and provides inconsistent measurements. The homogeneous wave equation solution for modal time behavior aids in the separation of the aliasing effects, but difficulty arises when the ring is spinning past the sensors and control forces are being applied.

Hamm, A. L.↗

The penetration of mountain waves into the middle atmosphere

A linear nonhydrostatic model of gravity waves forced by a bell-shaped ridge is used to investigate the penetration of mountain waves into the stratosphere and mesosphere during winter and fall. Gravity waves with horizontal scales less than 30 km are found to be trapped near the tropopause and the stratopause in regions of strong winds. The effect of trapping these modes produces a disturbance whose structure broadens with height. In the mesosphere the disturbance appears 20-40 km downstream from the forcing depending on the strength of the intervening winds. Wavebreaking associated with the mountain wave is predicted in the lower stratosphere as a result of wave superposition; no individual harmonic reaches breaking amplitude. In the mesosphere, wave breakdown is more prevalent, and the disturbance spectrum is relatively more monochromatic as a result of the filtering of the shorter scale modes by the lower atmosphere.

Schoeberl, M. R.↗

On the Causes of and Long Term Changes in Eurasian Heat Waves

The MERRA reanalysis, other observations, and the GEOS-S model have been used to diagnose the causes of Eurasian heat waves including the recent extreme events that occurred in Europe during 2003 and in Russia during 2010. The results show that such extreme events are an amplification of natural patterns of atmospheric variability (in this case a particular large-scale atmospheric planetary wave) that develop over the Eurasian continent as a result of internal atmospheric forcing. The amplification occurs when the wave occasionally becomes locked in place for several weeks to months resulting in extreme heat and drying with the location depending on the phase of the upper atmospheric wave. Model experiments suggest that forcing from both the ocean (SST) and land playa role phase-locking the waves. An ensemble of very long GEOS-S model simulations (spanning the 20th century) forced with observed SST and greenhouse gases show that the model is capable of generating very similar heat waves, and that they have become more extreme in the last thirty years as a result of the overall warming of the Asian continent.

Schubert, Siegfried↗

The QBO in Two GISS Global Climate Models: 1. Generation of the QBO

The adjustment of parameterized gravity waves associated with model convection and finer vertical resolution has made possible the generation of the quasi-biennial oscillation (QBO) in two Goddard Institute for Space Studies (GISS) models, GISS Middle Atmosphere Global Climate Model III and a climate/middle atmosphere version of Model E2. Both extend from the surface to 0.002 hPa, with 2deg × 2.5deg resolution and 102 layers. Many realistic features of the QBO are simulated, including magnitude and variability of its period and amplitude. The period itself is affected by the magnitude of parameterized convective gravity wave momentum fluxes and interactive ozone (which also affects the QBO amplitude and variability), among other forcings. Although varying sea surface temperatures affect the parameterized momentum fluxes, neither aspect is responsible for the modeled variation in QBO period. Both the parameterized and resolved waves act to produce the respective easterly and westerly wind descent, although their effect is offset in altitude at each level. The modeled and observed QBO influences on tracers in the stratosphere, such as ozone, methane, and water vapor are also discussed. Due to the link between the gravity wave parameterization and the models' convection, and the dependence on the ozone field, the models may also be used to investigate how the QBO may vary with climate change.

middle atmosphere dynamics↗

Self-Consistent Model of Magnetospheric Ring Current and Electromagnetic Ion Cyclotron Waves: The 2-7 May 1998 Storm

A complete description of a self-consistent model of magnetospheric ring current interacting with electromagnetic ion cyclotron waves is presented. The model is based on the system of two kinetic equations; one equation describes the ring current ion dynamics, and another equation describes the wave evolution. The effects on ring current ions interacting with electromagnetic ion cyclotron waves and back on waves are considered self-consistently by solving both equations on a global magnetospheric scale under nonsteady state conditions. The developed model is employed to simulate the entire 2-7 May 1998 storm period. First, the trapped number fluxes of the ring current protons are calculated and presented along with comparison with the data measured by the three- dimensional hot plasma instrument Polar/HYDRA. Incorporating in the model the wave-particle interaction leads to much better agreement between the experimental data and the model results. Second, examining of the wave (MLT, L shell) distributions produced by the model during the storm progress reveals an essential intensification of the wave emission about 2 days after the main phase of the storm. This result is well consistent with the earlier ground-based observations. Finally, the theoretical shapes and the occurrence rates of the wave power spectral densities are studied. It is found that about 2 days after the storm s main phase on 4 May, mainly non-Gaussian shapes of power spectral densities are produced.

Khazanov, G. V.↗

Earth's compressional foreshock boundary revisited Observations by the ISEE 1 magnetometer

A 'solar foreshock coordinate' (SFC) system is introduced in which the positions of foreshock components can be collated with a minimum of assumptions about the physical processes involved. Location behind the interplanetary magnetic field (IMF) tangent surface to the bow shock and rotational symmetry around the solar wind flow (X) axis are the only presuppositions. The SFC system has been applied to over 300 observations of the boundary of the ULF compressional waves recorded by ISEE 1's magnetometer in 1978 and 1979. The boundary locations form a coherent pattern in the SFC frame. A selection of those cases for which the cone angle of the IMF was between 40 and 50 deg, corresponding to the average stream angle, yields a least square line whose mapping back to the solar ecliptic coordinates frame has slope of about 85 deg, very close to that of the tangent ULF boundary deduced earlier from more primitive methods with entirely different data sets. The line, being parallel to neither the IMF, the typical reflected beam, nor the shock, cannot be compatible with any model of wave production by beam-solar wind interaction that depends on uniform beam distribution or fixed growth rate. Rather, its tangency suggests the influence of a separate, escaping ion population.

Greenstadt, E. W.↗

The study of coronal plasma structures and fluctuations with Faraday rotation measurements

We report dual-frequency, polarimetric measurements of Faraday rotation of extragalactic radio sources viewed through the solar corona. The observations were made at the Very Large Array in 1990 during solar maximum. Of the nine observed, an excess rotation measure of -12.6 rad/sq m was detected for one source (0010+005), which was observed at an elongation of about 9 solar radii. This measurement is in fair agreement with an a priori model rotation measure of -8.6 rad/sq m estimated from coronal potential field models and the electron density model of Paetzold et al. (1992). Our measurement provides a value for the coronal magnetic field strength at 9 solar radii given a knowledge of the magnetic field sector structure, of 12.5 +/- 2.3 mG. Rotation measurements of 0010+005 were made approximately once per hour over an 11 hr period. During this interval, a slow change of about 1 rad/sq m/hr in rotation measure was detected. Although we are not absolutely certain that this drift is not unremoved ionospheric Faraday rotation, extensive analysis of data from the other sources suggests that this is not the case (Sakurai & Spangler 1994). The very long timescale for this variation argues against the agency of magnetohydrodynamics (MHD) waves, and we suggest occultation of 0010+005 by relatively static plasma structures in the corona. We filtered our rotation measure time series to search for variations on an hourly timescale, such as those reported by Hollweg et al. (1992), which could be attributed to coronal MHD waves. We were unable to detect such fluctuations and can report only an upper limit to the rms variation of 1.6 rad/sq m. This upper limit is of the same order, but slightly larger than the values typically reported by Hollweg et al. (1982). This upper limit to the rotation measure fluctuations limits the dimensionless wave amplitude (delta B)/B in the corona to be less than 0.7. Using the number, we estimate the MHD wave flux at the coronal base to be less than 1.6 x 10(exp 5) ergs/sq cm/s. This is less than the amount of wave energy flux required by wave-driven models of the solar wind. Finally, we discuss a number of ways in which such observations could be improved in the future.

Sakurai, Takayuki↗

Laser Beam Propagation Through Inhomogeneous Media with Shock-Like Profiles: Modeling and Computing

Wave propagation in inhomogeneous media has been studied for such diverse applications as propagation of radiowaves in atmosphere, light propagation through thin films and in inhomogeneous waveguides, flow visualization, and others. In recent years an increased interest has been developed in wave propagation through shocks in supersonic flows. Results of experiments conducted in the past few years has shown such interesting phenomena as a laser beam splitting and spreading. The paper describes a model constructed to propagate a laser beam through shock-like inhomogeneous media. Numerical techniques are presented to compute the beam through such media. The results of computation are presented, discussed, and compared with experimental data.

Adamovsky, Grigory↗

Polarization and wavelength diversities of Gulf Stream fronts imaged by AIRSAR

During the 1990 Gulf Stream Experiment, NASA/JPL AIRSAR imaged the north edge of the Gulf Stream near the coast of Virginia. Simultaneous in-situ measurements of currents, temperatures, salinities, etc. were made for several crossings of the north edge by the R/V Cape Henlopen. Measurements identified two fronts with shearing and converging flows. The polarimetric SAR images from the fronts showed two bright linear features. One of them corresponds to the temperature front, which separated the warm Gulf Stream water to the south from a cool, freshwater filament to the north. The other line, located about 8 km north of the temperature front, is believed to correspond to the velocity front between the filament and the slope water. At these fronts, wave-current interactions produced narrow bands of steep and breaking waves manifesting higher radar returns in polarimetric SAR images. In general, our AIRSAR imagery shows that the signal-to-clutter ratio of radar cross sections for the temperature front is higher than that of the velocity front. In this paper, we study the polarization and wavelength diversities of radar response of these two fronts using the P-, L-, and C-Band Polarimetric SAR data. The north-south flight path of the AIRSAR crossed the temperature front several times and provided valuable data for analysis. Three individual passes are investigated. We found that for the temperature front, the cross-pol (HV) responses are much higher than co-pol responses (VV and HH), and that P-Band HV has the highest signal to clutter ratio. For the velocity front, the ratio is the strongest in P-Band VV, and it is indistinguishable for all polarizations in C-Band. The radar cross sections for all three polarization (HH, HV, and VV) and for all three bands are modelled using an ocean wave model and a composite Bragg scattering model. In our initial investigations, the theoretical model agrees qualitatively with the AIRSAR observations.

Lee, J. S.↗