Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “three phase 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 1,117 records · Page 62

Observational evidence of preferred flow regimes in the Northern Hemisphere winter stratosphere

Ten years of stratospheric geopotential height data are analyzed in an attempt to determine whether there are preferred flow regimes in the Northern Hemisphere winter stratosphere. The data are taken from Stratospheric Sounding Units on board NOAA satellites. The probability density estimate of the amplitude of the wavenumber 1 10-mb height is found to be bimodal. The density distribution is composed of a dominant large-amplitude mode and a less frequent low-amplitude mode. When the wavenumber 1 10-mb height data are projected onto the phase plane defined by the 10-mb zonal-mean winds and wavenumber 1 100-mb heights, three preferred regimes are evident. The small-amplitude mode separates into a strong zonal wind-weak wave regime and a weak zonal wind-weak wave regime. The large-amplitude mode is an intermediate zonal wind-strong wave regime. Transitions between the large-amplitude regime and the weak zonal wind-weak wave regime are found to be associated with major stratospheric warmings. The clustering of the stratospheric data into the preferred flow regimes is interpreted in light of the bifurcation properties of the Holton and Mass model. The interannual variability of the Northern Hemisphere winter stratosphere is interpreted in terms of the relative frequency of the observed preferred regimes.

Pierce, R. B.↗

A calorimetric determination of the enthalpy of formation and a description of the defect structure of the ordered beta-phase /Ni, Cu/ /1-x/ Al/x/

In order to describe thermodynamically the defect structure of an ordered B-Hume-Rothery phase, the heat of formation of (Ni,Cu)(1-x)Al(x) was measured at 1100 K as a function of concentration in the range x (sub Al) = 0.4 and 0.55 for three substitution rations x (sub Ni)/x (sub Cu) = infinity; 11; 5. The heat of formation of the NiAl beta-phase is strongly negative. For the stoichiometric composition it is -72.2 kJ/g-atom. On both the nickel-rich side and the aluminum-rich side the magnitude of the enthalpy of formation decreases linearly with concentration. Substitution of nickel for copper decreases the magnitude of the enthalpy of formation over the entire homogeneity range for the phase (Ni,Cu)(1-x)Al(x). The curve for the enthalpy of formation as well as the literature values for the chemical potential of aluminum are described with great accuracy by the disorder model of Wagner-Schottky.

Henig, E. T.↗

Instrumentation and robotic image processing using top-down model control

A top-down image processing scheme is described. A three-dimensional model of a robotic working environment, with robot manipulators, workpieces, cameras, and on-the-scene visual enhancements is employed to control and direct the image processing, so that rapid, robust algorithms act in an efficient manner to continually update the model. Only the model parameters are communicated, so that savings in bandwidth are achieved. This image compression by modeling is especially important for control of space telerobotics. The background for this scheme lies in an hypothesis of human vision put forward by the senior author and colleagues almost 20 years ago - the Scanpath Theory. Evidence was obtained that repetitive sequences of saccadic eye movements, the scanpath, acted as the checking phase of visual pattern recognition. Further evidence was obtained that the scanpaths were apparently generated by a cognitive model and not directly by the visual image. This top-down theory of human vision was generalized in some sense to the frame in artificial intelligence. Another source of the concept arose from bioengineering instrumentation for measuring the pupil and eye movements with infrared video cameras and special-purpose hardware.

Stark, Lawrence↗

A provisional RR Lyrae distance scale

Hydrodynamic pulsation models are matched with observations of globular cluster RRc stars to derive masses, luminosities and temperatures as functions of two observables: pulsation period, P1, and Fourier phase parameter, phi sub 31. We find that mean RRc masses and luminosities increase, and mean temperatures fall, with decreasing cluster metallicity. The Oosterhoff dichotomy is explained as mainly a temperature effect, while hints of a new dichotomy (in RRc mass and luminosity) are found among certain Oo II clusters. A provisional RR Lyrae distance scale emerges for the RRc stars in the form of a relation between luminosity and the two observables, P1 and phi sub 31. Finally, the scheme we propose passes three independent tests: (1) it reproduces the observed hierarchy of relative luminosity among a large sample of RRc stars in Omega Centauri; (2) it yields mean RRc masses which are fully consistent with the RRd masses derived for the clusters M68 and M15; and (3) it gives an LMC distance modulus in agreement with that obtained by other methods.

Simon, Norman R.↗

Analysis of Petal Rotation Trajectory Characteristics

In this study, the characteristics of petal rotation trajectories are explored in both the two-body and circular restricted three-body problem (CRTBP) models. Petal rotation trajectories alternate long and short resonances of different kinds to rotate the line of apsides. They are typically computed using the patched conic model, and they are used in a number of different missions and mission concepts including Cassini, JUICE, and Europa mission concepts. Petal rotation trajectories are first analyzed here using the patched conic model to quantify their characteristics and search for cases with fast rotation of the line of apsides. When they are computed in the CRTBP, they are unstable periodic orbits with corresponding stable and unstable manifolds. The characteristics of these orbits are explored from a dynamical systems perspective in the second phase of the study.

Resonance↗

Long-Term Creep and Creep Rupture Behavior of Woven Ceramic Matrix Composites

Tensile creep behavior of SiC/SiNC ceramic matrix composites at elevated temperatures and at various stress levels have been investigated for turbine engine applications. The objective of this research is to present creep behavior of SiC/SiCN composites at stress levels above and below the monotonic proportional limit strength and predict the life at creep rupture conditions. Tensile creep-rupture tests were performed on an Instron 8502 servohydraulic testing machine at constant load conditions up to a temperature limit of 1000 C. Individual creep curves indicate three stages such as primary, secondary, and tertiary. The creep rate increased linearly at an early stage and then gradually became exponential at higher strains. The stress exponent and activation energy were also obtained at 700 and 1000 C. The specimen lifetime was observed to be 55 hrs at 121 MPa and at 700 C. The life span reduced to 35 hrs at 143 MPa and at 1000 C. Scanning electron microscopy observations revealed significant changes in the crystalline phases and creep damage development. Creep failures were accompanied by extensive fiber pullout, matrix cracking, and debonding along with fiber fracture. The creep data was applied to Time-Temperature-Stress superposition model and the Manson-Haferd parametric model for long-time life prediction.

Haque, A.↗

Numerical Predictions of Wind Turbine Power and Aerodynamic Loads for the NREL Phase II and IV Combined Experiment Rotor

Accurate, reliable and robust numerical predictions of wind turbine rotor power remain a challenge to the wind energy industry. The literature reports various methods that compare predictions to experiments. The methods vary from Blade Element Momentum Theory (BEM), Vortex Lattice (VL), to variants of Reynolds-averaged Navier-Stokes (RaNS). The BEM and VL methods consistently show discrepancies in predicting rotor power at higher wind speeds mainly due to inadequacies with inboard stall and stall delay models. The RaNS methodologies show promise in predicting blade stall. However, inaccurate rotor vortex wake convection, boundary layer turbulence modeling and grid resolution has limited their accuracy. In addition, the inherently unsteady stalled flow conditions become computationally expensive for even the best endowed research labs. Although numerical power predictions have been compared to experiment. The availability of good wind turbine data sufficient for code validation experimental data that has been extracted from the IEA Annex XIV download site for the NREL Combined Experiment phase II and phase IV rotor. In addition, the comparisons will show data that has been further reduced into steady wind and zero yaw conditions suitable for comparisons to "steady wind" rotor power predictions. In summary, the paper will present and discuss the capabilities and limitations of the three numerical methods and make available a database of experimental data suitable to help other numerical methods practitioners validate their own work.

Duque, Earl P. N.↗

Flow interaction experiment. Volume 1: Aerothermal modeling, phase 2

An experimental and computational study is reported for the flow of a turbulent jet discharging into a rectangular enclosure. The experimental configurations consisting of primary jets only, annular jets only, and a combination of annular and primary jets are investigated to provide a better understanding of the flow field in an annular combustor. A laser Doppler velocimeter is used to measure mean velocity and Reynolds stress components. Major features of the flow field include recirculation, primary and annular jet interaction, and high turbulence. A significant result from this study is the effect the primary jets have on the flow field. The primary jets are seen to create statistically larger recirculation zones and higher turbulence levels. In addition, a technique called marker nephelometry is used to provide mean concentration values in the model combustor. Computations are performed using three levels of turbulence closures, namely k-epsilon model, algebraic second moment (ASM), and differential second moment (DSM) closure. Two different numerical schemes are applied. One is the lower-order power-law differencing scheme (PLDS) and the other is the higher-order flux-spline differencing scheme (FSDS). A comparison is made of the performance of these schemes. The numerical results are compared with experimental data. For the cases considered in this study, the FSDS is more accurate than the PLDS. For a prescribed accuracy, the flux-spline scheme requires a far fewer number of grid points. Thus, it has the potential for providing a numerical error-free solution, especially for three-dimensional flows, without requiring an excessively fine grid. Although qualitatively good comparison with data was obtained, the deficiencies regarding the modeled dissipation rate (epsilon) equation, pressure-strain correlation model, and the inlet epsilon profile and other critical closure issues need to be resolved before one can achieve the degree of accuracy required to analytically design combustion systems.

Nikjooy, M.↗

Flow interaction experiment. Volume 2: Aerothermal modeling, phase 2

An experimental and computational study is reported for the flow of a turbulent jet discharging into a rectangular enclosure. The experimental configurations consisting of primary jets only, annular jets only, and a combination of annular and primary jets are investigated to provide a better understanding of the flow field in an annular combustor. A laser Doppler velocimeter is used to measure mean velocity and Reynolds stress components. Major features of the flow field include recirculation, primary and annular jet interaction, and high turbulence. A significant result from this study is the effect the primary jets have on the flow field. The primary jets are seen to create statistically larger recirculation zones and higher turbulence levels. In addition, a technique called marker nephelometry is used to provide mean concentration values in the model combustor. Computations are performed using three levels of turbulence closures, namely k-epsilon model, algebraic second moment (ASM), and differential second moment (DSM) closure. Two different numerical schemes are applied. One is the lower-order power-law differencing scheme (PLDS) and the other is the higher-order flux-spline differencing scheme (FSDS). A comparison is made of the performance of these schemes. The numerical results are compared with experimental data. For the cases considered in this study, the FSDS is more accurate than the PLDS. For a prescribed accuracy, the flux-spline scheme requires a far fewer number of grid points. Thus, it has the potential for providing a numerical error-free solution, especially for three-dimensional flows, without requiring an excessively fine grid. Although qualitatively good comparison with data was obtained, the deficiencies regarding the modeled dissipation rate (epsilon) equation, pressure-strain correlation model, and the inlet epsilon profile and other critical closure issues need to be resolved before one can achieve the degree of accuracy required to analytically design combustion systems.

Nikjooy, M.↗

Effects of Thermomechanical History on the Tensile Behavior of Nitinol Ribbon

Shape memory alloys (SMAs) have enormous potential for a wide variety of applications. A large body of work exists on the characterization of the microstructure and stress-strain behavior of these alloys, Nitinol (NiTi) in particular. However, many attributes of these materials are yet to be fully understood. Previous work at NASA Langley Research Center (LaRC) has included fabrication of hybrid composite specimens with embedded Nitinol actuators and modeling of their thermomechanical behavior. An intensive characterization effort has been undertaken to facilitate fundamental understanding of this alloy and to promote implementation of Nitinol in aerospace applications. Previous work revealed attributes of the Nitinol ribbon that were not easily rationalized with existing data in the literature. In particular, tensile behavior at ambient temperature showed significant dependence on the thermomechanical history prior to testing. The present work is focused on characterizing differences in the microstructure of Nitinol ribbons exposed to four different thermomechanical histories and correlation of the microstructure with tensile properties. Differential scanning calorimetry (DSC) and x-ray diffraction (XRD) analysis were employed to rationalize the microstructures present after exposure to various thermomechanical histories. Three of the Nitinol ribbon conditions were reversible upon heating (in the DSC) through the reverse transformation temperature (A(sub f) to transform the microstructure to austenite. However, the prior thermomechanical conditioning for the Nitinol ribbon that reflected the entire fabrication procedure (4% thermal cycle condition) was found to have an irreversible effect on the microstructure, as it remained unchanged after repeated complete thermal cycles. Tensile tests were conducted to determine the effect of prior thermomechancal conditioning on both the tensile behavior of the Nitinol ribbons and the stress state of the microstructure. The stress-strain behavior of the Nitinol actuators appears to be governed by the interplay between two major variables: namely, microstructural constituents such as the R-phase and the martensite; and the stress state of these constituents (whether twinned with low residual stresses, or detwinned with high residual stresses). The most significant difference in the stress-strain behavior of the four conditions, the critical stress required to achieve an initial stress plateau, was found to depend on both the amount and stress state (twinned or detwinned) of R-phase present in the initial microstructure. Thus, the effect of prior thermomechanical processing is critical to the resulting tensile behavior of the Nitinol actuator. For numerical modeling inputs one must take into account the entire fabrication process on the Nitinol actuator.

Lach, Cynthia L.↗

Numerical investigation of the control of the secondary instability process in boundary layers

A numerical model has been developed for investigating boundary layer transition control for a flat plate boundary layer. Active control of a periodically forced boundary layer in an incompressible fluid is studied using surface heating techniques. The spatially evolving boundary layer is simulated. The Navier-Stokes and energy equations are integrated using a fully implicit finite difference/spectral method. Temperature perturbations are introduced locally along finite heater strips to directly attenuate the instability waves in the flow. A feedback control loop is employed in which a downstream sensor is used to monitor wall shear stress fluctuations. Active control of small amplitude two-dimensional and three-dimensional disturbances is numerically simulated. With proper phase control, in-phase reinforcement and out-of-phase attenuation are demonstrated. A receptivity study of the localized temperature perturbations is made. It is shown that narrow heater strips are more receptive in that they maximize the amplitude level of the disturbances in the flow. Active control of the early stages of the fundamental breakdown process is also numerically simulated. Control is achieved with either two-dimensional or three-dimensional control inputs.

Kral, L. D.↗

The ocean response to operational westerly wind bursts during the 1991 - 1992 El Nino

Numerical simulations of the remotely forced ocean response to westerly wind bursts prior to and during the 1991 - 1992 El Nino are examined; the models are forced by wind stress from the U.S. Navy's atmospheric global operational analysis/forecast system. The study focuses on (1) the relative response of the first and second internal modes to a single episode of westerly bursts; (2) the role of westerly bursts in producing the eastern Pacific sea level variations from October 1990 to February 1992; and (3) the relative importance of the remotely forced sea level responses generated by central and western Pacific wind anomalies. The simulation use the Naval Research Laboratory global multilayer formulation; the suite of experiments includes hydrodynamic simulations that use both one- and three- active-layer reduced gravity models as well as an experiment that also includes thermodynamic effects. The models are executed from January 1, 1990, to March 1, 1992, a period that includes 10 significant westerly wind bursts or burst clusters and the 1991 - 1992 El Nino event. The numerical experiments reveal an ability to accurately simulate the eastern Pacific sea level variations during this period. In response to a single burst, the three-layer hydrodynamic simulation reveals that the second internal mode Kelvin wave yields a sea level change at the eastern boundary that is approximately one third that of the first mode and a surface velocity signature that is equivalent to the first mode. During the onset of the El Nino event the inclusion of higher modes also produces a more realistic representation of the observed eastern boundary sea level signal. Furthermore, by comparing the model response to particular wind bursts with the observed sea level at Baltra, Galapagos, a value of 2.5 - 2.6 m/s is suggested as the most appropriate mean speed for the first internal mode Kelvin wave during the onset phase. A scenario for the onset of the 1991 - 1992 El Nino is given.

Kindle, John C.↗

Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows

The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.

Source record↗

X-Ray and Infrared Observations of Embedded Young Stars in NGC 2264

Images of the NGC 2264 star-forming region, which we have acquired with the XMM-Newton spacecraft, reveal strong X-ray emission from three deeply embedded (Av > 10 mag) young stellar objects in the vicinity of Allen's infrared source (AFGL 989 = IRS 1) and Castelaz & Grasdalen s infrared source (RNO-EW = IRS 2). Thermal plasma models for the brightest source in X-rays, located 11 southwest of Allen's star, yield a quasi-steady luminosity of Lx = 10 ergs s-1 and an extraordinarily high X-ray temperature of 100 MK. The high temperature is consistent with the presence of emission lines of Fe xxv and Fe xxvi at photon energies of 6.7 and 6.9 keV, respectively. An even higher temperature of nearly 140 MK was observed during the rise phase of a powerful impulsive X-ray flare of another young star in the IRS 2 region. Moderate-resolution near-infrared (1-4 um) spectra of the embedded objects, obtained at the NASA Infrared Telescope Facility, exhibit deep water ice absorption bands, as well as a variety of emission and absorption features of H I, CO, and both neutral and ionized metals.

Simon, Theordore↗

Application of Extended Messinger Models to Complex Geometries

Since, ice accretion can significantly degrade the performance and the stability of an airborne vehicle, it is imperative to be able to model it accurately. While ice accretion studies have been performed on airplane wings and helicopter blades in abundance, there are few that attempt to model the process on more complex geometries such as fuselages. This paper proposes a methodology that extends an existing in-house Extended Messinger solver to complex geometries by introducing the capability to work with unstructured grids and carry out spatial surface streamwise marching. For the work presented here commercial solvers such as STAR-CCM+ and ANSYS Fluent are used for the flow field and droplet dispersed phase computations. The ice accretion is carried out using an in-house icing solver called GT-ICE. The predictions by GT-ICE are compared to available experimental data, or to predictions by other solvers such as LEWICE and STAR-CCM+. Three different cases with varying levels of complexity are presented. The first case considered is a commercial transport airfoil, followed by a three-dimensional MS(1)-317 swept wing. Finally, ice accretion calculations performed on a Robin fuselage have been discussed. Good agreement with experimental data, where applicable, is observed. Differences between the ice accretion predictions by different solvers have been discussed.

Gupta, Avani↗

Intracranial Effects of Artificial Gravity: A 3T MRI Study

INTRODUCTION Spaceflight associated neuro-ocular syndrome (SANS) is characterized by the development of optic disc edema, posterior globe flattening, choroidal/retinal folds and hyperopic refractive errors1. SANS is hypothesized to be a result of headward fluid shifts that invariably occurs in the microgravity environment. As a countermeasure, artificial gravity (AG) through centrifugation has been proposed to reduce this headward fluid shift, however there is no current proof of benefit. The goal of this study was to determine if the application of AG can prevent or reduce known changes in brain volumetry, internal carotid artery (ICA) stroke volume and cerebral spinal fluid (CSF) flow velocity that occurs during simulated chronic headward fluid shift using head down tilt bed rest (HDTBR) methodology2 as an indicator of countermeasure efficacy. METHODS Healthy volunteers were recruited for an IRB approved HDTBR study performed at the German Aerospace Center in Cologne, Germany. Strict six-degree HDTBR was used as a spaceflight analog to induce a continuous headward fluid shift. HDTBR was carried out for 60 days for all subjects. Short-arm centrifugation was utilized to generate AG equating to ~0.3g of acceleration at the level of the eye. The subjects were divided equally into three groups: NoAG (control; n=8), daily intermittent AG (6 x 5 min iAG; n=8), and daily continuous 30 min (cAG; n=8). All studies were performed on a single dedicated 3T MRI Scanner. Pulse-gated MRI phase-contrast flow imaging was used to quantify ICA stroke volume and peak-to-peak CSF flow velocity in the mid cerebral aqueduct. 3D-SPGR was acquired for volumetric segmentation of the brain and CSF spaces. MRI acquisitions were obtained at baseline (BDC), 14 days into HDTBR (HDTBR14), 52 days into HDTBR (HDTBR52) and 3-5 days after HDTBR (recovery, R+3/5).The data were analyzed by the mixed model, which included intervention and time (BDC, HDTBR 14, HDTBR 52, R+3/5) as the fixed effects and included subject as the random effect.RESULTS24 healthy subject volunteers (16 men, 8 women, mean age = 33 years ± 9 [standard deviation] and mean BMI = 24.3 kg/m2 ± 2.0) successfully completed all phases of the study. Strict six-degree HDTBR was characterized by progressive and statistically significant (p<.01) increases in mean combined brain and CSF volumes and mean aqueductal CSF peak-to-peak flow velocity, as well as statistically significant (p<.01) progressive decrease in mean ICA stroke volume from baseline to 52 days post intervention (Figs. 1-3). Compared to baseline, only combined brain and CSF volumes did not return to baseline values in the recovery period (p=NS). Neither iAG nor cAG exerted any significant effects on the measured MRI brain parameters as compared to HDTBR alone (p=NS). CONCLUSION Our results indicate that HDTBR at 6-degrees was effective in producing alterations in ICA stroke volume, aqueductal CSF flow velocity, and combined brain and CSF volumetric change that is associated with chronic headward fluid shift. Short duration, 30-min daily exposure to either iAG or cAG appears to be insufficient in preventing or reducing the effects of chronic HDTBR and thus may not be a suitable countermeasure as currently deployed. AG protocol modifications, including increased duration and magnitude of exposure, should be considered for future research.

L A Kramer↗

Crystal growth from the vapour phase (ASTP experiment MA-085)

The positive effects of microgravity on crystal quality and fundamental properties of the vapor transport reaction were established by analyzing the results of three transport experiments on multi-component systems performed during the Apollo Soyuz Test Project (ASTP) mission. The systems employed were GeSe(0.99)Te(0.01) - GeI4 (A), GeS(0.98)Se(0.02) - GeCl4 (B), and GeS - GeCl4-Ar (C). The crystallographic analysis is based on a direct comparison of space and ground-based (prototype) crystals, employing X-ray diffraction, microprobe, microscopic, and chemical etching techniques. The results demonstrate a considerable improvement of the space-grown crystals in terms of chemical and crystalline homogeneity, surface morphology, and bulk perfection relative to ground specimens. The confirmation of greater mass transport rates than predicted for micro-gravity by present vapor transport models is of basic scientific and practical importance for the improvement of transport models and techniques on earth and for the development of a transport model for space conditions. The internal consistency and agreement between ASTP and Skylab results obtained for different compounds, transport agents, and temperature gradients strongly support the validity of these observations. The combined findings of improved crystal quality and high mass transport rates are of technological significance for space-processing applications.

Wiedemeier, H.↗

Modeling the radiant transfers of sparse vegetation canopies

The scattering dynamics of sparse vegetation canopies were studied within the framework of the three-dimensional radiative transfer model of Kimes (1984). The model was upgraded by including an algorithm for the anisotropic scattering of a soil boundary. Validation of the model was carried out using measured directional reflectance data for two canopies exhibiting typical scattering behavior with low and intermediate vegetation density. The canopies were: an orchard grass canopy; and a hard wheat canopy. A number of factors were found contributing to the final reflectance distribution of the canopies, including: (1) the strong anisotropic scattering properties of the soil; (2) the geometric effect of the vegetation probability gap function on the soil anisotropy and solar irradiance; and (3) the anisotropic scattering of vegetation which is controlled by the phase function and the layering of leaves. The application of the theoretical results to the development of earth-observing sensor systems is discussed.

Kimes, D. S.↗