Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “small scale”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The Pantheon+ Analysis: Evaluating Peculiar Velocity Corrections in Cosmological Analyses with Nearby Type Ia Supernovae

Separating the components of redshift due to expansion and peculiar motion in the nearby universe (z < 0.1) is critical for using Type Ia Supernovae (SNe Ia) to measure the Hubble constant (H 0 ) and the equation-of-state parameter of dark energy (w). Here, we study the two dominant "motions" contributing to nearby peculiar velocities: large-scale, coherent-flow (CF) motions and small-scale motions due to gravitationally associated galaxies deemed to be in a galaxy group. We use a set of 584 low-z SNe from the Pantheon+ sample, and evaluate the efficacy of corrections to these motions by measuring the improvement of SN distance residuals. We study multiple methods for modeling the large and small-scale motions and show that, while group assignments and CF corrections individually contribute to small improvements in Hubble residual scatter, the greatest improvement comes from the combination of the two (relative standard deviation of the Hubble residuals, Rel. SD, improves from 0.167 to 0.157 mag). We find the optimal flow corrections derived from various local density maps significantly reduce Hubble residuals while raising H 0 by ~0.4 km s -1 Mpc -1 as compared to using CMB redshifts, disfavoring the hypothesis that unrecognized local structure could resolve the Hubble tension. We estimate that the systematic uncertainties in cosmological parameters after optimally correcting redshifts are 0.06–0.11 km s -1 Mpc -1 in H 0 and 0.02–0.03 in w which are smaller than the statistical uncertainties for these measurements: 1.5 km s -1 Mpc -1 for H 0 and 0.04 for w.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2001

Microgravity Propellant Tank Geyser Analysis and Prediction

An established correlation for geyser height prediction of an axial jet inflow into a microgravity propellant tank was analyzed and an effort to develop an improved correlation was made. The original correlation, developed using data from ethanol flow in small-scale drop tower tests, uses the jet-Weber number and the jet-Bond number to predict geyser height. A new correlation was developed from the same set of experimental data using the jet-Weber number and both the jet-Bond number and tank-Bond number to describe the geyser formation. The resulting correlation produced nearly a 40% reduction in geyser height predictive error compared to the original correlation with experimental data. Two additional tanks were computationally modeled in addition to the small-scale tank used in the drop tower testing. One of these tanks was a 50% enlarged small-scale tank and the other a full-scale 2 m radius tank. Simulations were also run for liquid oxygen and liquid hydrogen. Results indicated that the new correlation outperformed the original correlation in geyser height prediction under most circumstances. The new correlation has also shown a superior ability to recognize the difference between flow patterns II (geyser formation only) and III (pooling at opposite end of tank from the bulk fluid region).

Thornton, Randall J.↗
DOE OSTI2021

The impact of magnetic fields on momentum transport and saturation of shear-flow instability by stable modes

The Kelvin–Helmholtz (KH) instability of a shear layer with an initially uniform magnetic field in the direction of flow is studied in the framework of 2D incompressible magnetohydrodynamics with finite resistivity and viscosity using direct numerical simulations. The shear layer evolves freely, with no external forcing, and thus broadens in time as turbulent stresses transport momentum across it. As with hydrodynamic KH, the instability here features a conjugate stable mode for every unstable mode in the absence of dissipation. Stable modes are shown to transport momentum up its gradient, shrinking the layer width whenever they exceed unstable modes in amplitude. In simulations with weak magnetic fields, the linear instability is minimally affected by the field, but enhanced small-scale fluctuations relative to the hydrodynamic case are observed. These enhanced fluctuations coincide with increased energy dissipation and faster layer broadening, with these features more pronounced in simulations with stronger fields. These trends result from the magnetic field reducing the effects of stable modes relative to the transfer of energy to small scales. As field strength increases, stable modes become less excited, thus transporting less momentum against its gradient. Furthermore, the energy that would otherwise transfer back to the driving shear because of the stable modes is instead allowed to cascade to small scales, where it is lost to dissipation. Approximations of the turbulent state in terms of a reduced set of modes are explored. While the Reynolds stress is well-described using just two modes per wavenumber at large scales, the Maxwell stress is not.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2010

Cascade and Dissipation of Solar Wind Turbulence at Electron Scales: Whistlers or Kinetic Alfv\'en Waves?

Over the past few decades, large-scales solar wind (SW) turbulence has been studied extensively, both theoretically and observationally. Observed power spectra of the low frequency turbulence, which can be described in the magnetohydrodynamic (MHD) limit, are shown to obey the Kolmogorov scaling, $k"{ -5/3 }$, down the local proton gyrofrequency ($C{ci} \sim O.l$-Hz). Turbulence at frequencies above $C{ci}$ has not been thoroughly investigated and remains far less well understood. Above $C{ ci}$ the spectrum steepens to $\sim f"{ -2.5}$ and a debate exists as to whether the turbulence has become dominated by dispersive kinetic Alfven waves (KA W) or by whistler waves, before it is dissipated at small scales, In a case study Sahraoui et al., PRL (2009) have reported the first direct determination of the dissipation range of solar wind turbulence near the electron gyroscale using the high resolution Cluster magnetic and electric field data (up to $10"2$-Hz in the spacecraft reference frame). Above the Doppler-shifted proton scale $C{\rho i}$ a new inertial range with a scaling $\sim f"{ -2.3}$ has been evidenced and shown to remarkably agree with theoretical predictions of a quasi-two-dimensional cascade into KA W turbulence. Here, we use a wider sample of data sets of small scale SW turbulence under different plasma conditions, and investigate under which physical criteria the KA W (or the whistler) turbulence may be observed to carry out the cascade at small scales, These new observations/criteria are compared to the predictions on the cascade and the (kinetic) dissipation from the Vlasov theory. Implications of the results on the heating problem of the solar wind will be discussed.

Sahraoui, Fouad↗
NASA NTRS2011

Scale Dependence of Magnetic Helicity in the Solar Wind

We determine the magnetic helicity, along with the magnetic energy, at high latitudes using data from the Ulysses mission. The data set spans the time period from 1993 to 1996. The basic assumption of the analysis is that the solar wind is homogeneous. Because the solar wind speed is high, we follow the approach first pioneered by Matthaeus et al. by which, under the assumption of spatial homogeneity, one can use Fourier transforms of the magnetic field time series to construct one-dimensional spectra of the magnetic energy and magnetic helicity under the assumption that the Taylor frozen-in-flow hypothesis is valid. That is a well-satisfied assumption for the data used in this study. The magnetic helicity derives from the skew-symmetric terms of the three-dimensional magnetic correlation tensor, while the symmetric terms of the tensor are used to determine the magnetic energy spectrum. Our results show a sign change of magnetic helicity at wavenumber k approximately equal to 2AU(sup -1) (or frequency nu approximately equal to 2 microHz) at distances below 2.8AU and at k approximately equal to 30AU(sup -1) (or nu approximately equal to 25 microHz) at larger distances. At small scales the magnetic helicity is positive at northern heliographic latitudes and negative at southern latitudes. The positive magnetic helicity at small scales is argued to be the result of turbulent diffusion reversing the sign relative to what is seen at small scales at the solar surface. Furthermore, the magnetic helicity declines toward solar minimum in 1996. The magnetic helicity flux integrated separately over one hemisphere amounts to about 10(sup 45) Mx(sup 2) cycle(sup -1) at large scales and to a three times lower value at smaller scales.

Brandenburg, Axel↗
NASA NTRS2010

Downscaling Soil Moisture in the Southern Great Plains Through a Calibrated Multifractal Model for Land Surface Modeling Applications

Accounting for small-scale spatial heterogeneity of soil moisture (theta) is required to enhance the predictive skill of land surface models. In this paper, we present the results of the development, calibration, and performance evaluation of a downscaling model based on multifractal theory using aircraft!based (800 m) theta estimates collected during the southern Great Plains experiment in 1997 (SGP97).We first demonstrate the presence of scale invariance and multifractality in theta fields of nine square domains of size 25.6 x 25.6 sq km, approximately a satellite footprint. Then, we estimate the downscaling model parameters and evaluate the model performance using a set of different calibration approaches. Results reveal that small-scale theta distributions are adequately reproduced across the entire region when coarse predictors include a dynamic component (i.e., the spatial mean soil moisture ) and a stationary contribution accounting for static features (i.e., topography, soil texture, vegetation). For wet conditions, we found similar multifractal properties of soil moisture across all domains, which we ascribe to the signature of rainfall spatial variability. For drier states, the theta fields in the northern domains are more intermittent than in southern domains, likely because of differences in the distribution of vegetation coverage. Through our analyses, we propose a regional downscaling relation for coarse, satellite-based soil moisture estimates, based on ancillary information (static and dynamic landscape features), which can be used in the study area to characterize statistical properties of small-scale theta distribution required by land surface models and data assimilation systems.

Mascaro, Giuseppe↗
NASA NTRS1979

Sound production due to large-scale coherent structures

The acoustic pressure fluctuations due to large-scale finite amplitude disturbances in a free turbulent shear flow are calculated. The flow is decomposed into three component scales; the mean motion, the large-scale wave-like disturbance, and the small-scale random turbulence. The effect of the large-scale structure on the flow is isolated by applying both a spatial and phase average on the governing differential equations and by initially taking the small-scale turbulence to be in energetic equilibrium with the mean flow. The subsequent temporal evolution of the flow is computed from global energetic rate equations for the different component scales. Lighthill's theory is then applied to the region with the flowfield as the source and an observer located outside the flowfield in a region of uniform velocity. Since the time history of all flow variables is known, a minimum of simplifying assumptions for the Lighthill stress tensor is required, including no far-field approximations. A phase average is used to isolate the pressure fluctuations due to the large-scale structure, and also to isolate the dynamic process responsible. Variation of mean square pressure with distance from the source is computed to determine the acoustic far-field location and decay rate, and, in addition, spectra at various acoustic field locations are computed and analyzed. Also included are the effects of varying the growth and decay of the large-scale disturbance on the sound produced.

Gatski, T. B.↗
DOE OSTI2023

Temporal Subsampling Diminishes Small Spatial Scales in Recurrent Neural Network Emulators of Geophysical Turbulence

The immense computational cost of traditional numerical weather and climate models has sparked the development of machine learning (ML) based emulators. Because ML methods benefit from long records of training data, it is common to use data sets that are temporally subsampled relative to the time steps required for the numerical integration of differential equations. Here, we investigate how this often overlooked processing step affects the quality of an emulator's predictions. We implement two ML architectures from a class of methods called reservoir computing: (a) a form of Nonlinear Vector Autoregression (NVAR), and (b) an Echo State Network (ESN). Despite their simplicity, it is well documented that these architectures excel at predicting low dimensional chaotic dynamics. We are therefore motivated to test these architectures in an idealized setting of predicting high dimensional geophysical turbulence as represented by Surface Quasi-Geostrophic dynamics. In all cases, subsampling the training data consistently leads to an increased bias at small spatial scales that resembles numerical diffusion. Interestingly, the NVAR architecture becomes unstable when the temporal resolution is increased, indicating that the polynomial based interactions are insufficient at capturing the detailed nonlinearities of the turbulent flow. The ESN architecture is found to be more robust, suggesting a benefit to the more expensive but more general structure. Spectral errors are reduced by including a penalty on the kinetic energy density spectrum during training, although the subsampling related errors persist. Future work is warranted to understand how the temporal resolution of training data affects other ML architectures.

58 GEOSCIENCES↗
NASA NTRS2022

Atmospheric Energy Spectra in Global Kilometre-Scale Models

Eleven 40-day long integrations of five different global models with horizontal resolutions of less than 9 km are compared in terms of their global energy spectra. The method of normal-mode function decomposition is used to distinguish between balanced (Rossby wave; RW) and unbalanced (inertia-gravity wave; IGW) circulation. The simulations produce the expected canonical shape of the spectra, but their spectral slopes at mesoscales, and the zonal scale at which RW and IGW spectra intersect differ significantly. The partitioning of total wave energies into RWs an IGWs is most sensitive to the turbulence closure scheme and this partitioning is what determines the spectral crossing scale in the simulations, which differs by a factor of up to two. It implies that care must be taken when using simple spatial filtering to compare gravity wave phenomena in storm-resolving simulations, even when the model horizontal resolutions are similar. In contrast to the energy partitioning between the RWs and IGWs, changes in turbulence closure schemes do not seem to strongly affect spectral slopes, which only exhibit major differences at mesoscales. Despite their minor contribution to the global (horizontal kinetic plus potential available) energy, small scales are important for driving the global mean circulation. Our results support the conclusions of previous studies that the strength of convection is a relevant factor for explaining discrepancies in the energies at small scales. The models studied here produce the major large-scale features of tropical precipitation patterns. However, particularly at large horizontal wavenumbers, the spectra of upper tropospheric vertical velocity, which is a good indicator for the strength of deep convection, differ by factors of three or more in energy. High vertical kinetic energies at small scales are mostly found in those models that do not use any convective parameterisation.

DYAMOND↗
DOE OSTI2021

Spent Crystalline Silicotitanate Storage Study—Post AP-105DF Processing

A Tank-Side Cesium Removal (TSCR) system is under development by Washington River Protection Solutions to support initial production of immobilized low-activity waste (LAW) by feeding Hanford tank supernate from tank farms to the Hanford Waste Treatment and Immobilization Plant (WTP) LAW Facility. Tank waste supernate will be filtered to remove suspended solids and then Cs will be removed by processing it through crystalline silicotitanate (CST) ion exchange media manufactured by Honeywell UOP, LLC. The Cs-loaded CST columns will be stored indefinitely, with a goal of eventual CST removal and treatment. Thus, the spent CST needs to be recoverable after undetermined storage time. Previous testing with AP-105 simulant showed that rinsing the CST bed with 3 bed volumes (1.4 apparatus volumes [AVs]) of 0.1 M NaOH resulted in a dried bed that maintained flow characteristics indicative of ease of recovery. This study explored the intermediate conditions (between feed dried in place and the 1.4 AVs of 0.1 M NaOH rinse) to evaluate CST bed properties after: 1) stoppage with feed in place; 2) stoppage after draining feed; 3) stoppage after 0.7 AV of 0.1 M NaOH rinse through column; 4) stoppage after 1.5 AVs of 0.1 M NaOH rinse through column. Post processing, each column was heated at 50 °C for 19 days under pseudo-storage conditions to simulate the expected dried and stored CST bed conditions. Testing was conducted at the small scale (12-mL bed volume); actual, Cs-depleted, AP-105 tank waste was used as the feed. Post-dried CST bed physical properties (angle of repose and penetration depth) were measured to evaluate how CST moved and flowed. All process stop-conditions resulted in a solidified CST bed except for the final condition, 1.5 AVs of 0.1 M NaOH rinse. At this small scale, the three CST beds presented an issue for retrievability after the short storage period (19 days at 50 °C). The latter case confirmed the results from simulant testing. The testing was intended to provide a preliminary assessment of issues that may arise from desiccation of CST during storage with the indicated salt solutions in place. Since these were small-scale tests, the processing system did not scale to full scale conditions exactly; however, the tests did provide insight into the impact on the dried and stored CST bed after stopping processing at an earlier step (upset condition) than normal. These results indicate that if an upset condition occurs at TSCR, a dilute hydroxide rinse should be considered before the CST dries from internal heating.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1997

On the Computation of Sound by Large-Eddy Simulations

The effect of the small scales on the source term in Lighthill's acoustic analogy is investigated, with the objective of determining the accuracy of large-eddy simulations when applied to studies of flow-generated sound. The distribution of the turbulent quadrupole is predicted accurately, if models that take into account the trace of the SGS stresses are used. Its spatial distribution is also correct, indicating that the low-wave-number (or frequency) part of the sound spectrum can be predicted well by LES. Filtering, however, removes the small-scale fluctuations that contribute significantly to the higher derivatives in space and time of Lighthill's stress tensor T(sub ij). The rms fluctuations of the filtered derivatives are substantially lower than those of the unfiltered quantities. The small scales, however, are not strongly correlated, and are not expected to contribute significantly to the far-field sound; separate modeling of the subgrid-scale density fluctuations might, however, be required in some configurations.

Piomelli, Ugo↗
DOE OSTI2020

Spent Crystalline Silicotitanate Storage Study

Washington River Protection Solutions is working to support initial production of immobilized low-activity waste (LAW) by feeding Hanford tank supernate from tank farms to the Hanford Waste Treatment and Immobilization Plant (WTP) LAW Facility. This goal incorporates the design of a Tank-Side Cesium Removal (TSCR) system, which filters tank waste supernate to remove suspended solids and then removes Cs by processing it through crystalline silicotitanate (CST) ion exchange media manufactured by Honeywell UOP, LLC. The 137 Cs-depleted product is intended to be sent to the WTP for vitrification. The Cs-loaded CST columns will be stored indefinitely, with a goal of eventual CST removal and treatment. Thus, the spent CST needs to be recoverable. The testing described herein looks to potential upset process conditions where CST storage may be required before various rinse steps are completed. This study evaluated upset conditions at three sequential processing steps envisioned for TSCR (feed, 0.1 M NaOH rinse, water rinse), and in-column drying with compressed air (normal end step) to assess impacts on the nature of the CST bed. Testing was conducted at the small scale (12-mL bed volume); simulated AP-105 tank waste was used as the feed. Following process disruption, the CST bed was dried in place at 70 °C. Post-dried CST bed physical properties (angle of repose, penetration depth, particle morphological changes) were measured to evaluate how CST moved and flowed. The testing is intended to provide a preliminary assessment of issues that may arise from desiccation of CST with the indicated salt solutions in place. Since these were small-scale tests, the processing conditions will not match full scale conditions exactly; however, the tests do provide insight into the impact of stopping processing at an earlier step than normal. Except for the feed that was dried in-place, all other process stop-conditions showed the CST bed flowed well after drying. At this small scale, CST beds would not present an issue for retrievability. The feed that was dried in place had solidified into a rock-hard monolith with no movement possible. Samples had to be chipped from the surface.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Novel Lightweight, Low-Cost Heliostat for Concentrating Solar Power (SBIR Phase I Final Report)

The overall goal of this Department of Energy (DOE) SBIR Phase I project was to demonstrate the feasibility of the L.Garde, Inc. heliostat mirror facet for use in concentrated solar power applications. Prior to this Phase I project, L.Garde, Inc. developed a design for a lightweight and low cost mirror facet and successfully fabricated small scale (0.667 m 2 ) with a mirror surface flatness of 0.7 mrad. Therefore, the focus of the Phase I effort was to (1) demonstrate scalability by fabricating larger size mirrors while retaining the same mirror surface flatness, (2) demonstrate the design can withstand operating environment conditions through field testing, and (3) perform cost analysis to show the heliostat mirror facet design can meet DoE goals of < 12 kg/m2 and <$30/m 2 . During the Phase I project, L.Garde successfully fabricated large area (1.395 m 2 ) heliostat mirror facets, with a mass per unit area of 8 kg/m 2 and a cost of $29.80/m 2 . Surface flatness measurements (using contact probe technique) were performed on the mirror surface and results showed a flatness of 0.68 mrad, which were in agreement with previously fabricated small scale mirror facets. Slope error measurements was also performed on the mirror surface and results showed a slope error of 0.77 mrad which also was in agreement with previously fabricated small scale mirror facets. One large area mirror facet was subjected to a three (3) month field test. Surface flatness and slope error measurements were performed on the mirror surface after the field test and results were 0.663 mrad and 0.58 mrad respectively. The large scale heliostat mirror facets were also subjected to accelerated aging. However, during this exposure we observed failure due to the breaking of the mirror facets. Further investigation into the source of the failure showed that the mirror component alone showed signs of degradation during the same accelerated aging test. We believe that moisture ingress is occurring when subjected to high humidity, leading to buildup of water with the layers of protective paints and resulting in mirror breakage when water expands in freezing conditions. In order to mitigate or eliminate this issue, a new weatherproofing strategy may be implemented or a new mirror component may be used that is able to survive accelerated aging conditions.

14 SOLAR ENERGY↗
DOE OSTI2025

Multiscale Modeling Framework Using Element‐Based Galerkin Methods for Moist Atmospheric Limited‐Area Simulations

This paper presents a multiscale modeling framework (MMF) to model moist atmospheric limited-area weather. The MMF resolves large-scale convection using a coarse grid while simultaneously resolving local features through numerous fine local grids and coupling them seamlessly. Both large- and small-scale processes are modeled using the compressible Navier-Stokes equations within the Nonhydrostatic Unified Model of the Atmosphere (NUMA), and are discretized using a continuous element-based Galerkin method (spectral elements) with high-order basis functions. Consequently, the large-scale and small-scale models share the same dynamical core but have the flexibility to be adjusted individually. The proposed MMF method is tested in 2D and 3D idealized limited-area weather problems involving storm clouds produced by squall line and supercell simulations. Numerical results from the MMF showed enhanced representation of cloud processes compared to the coarse model.

Kang, Soonpil [Naval Postgraduate School, Monterey↗
DOE OSTI2023

Vibration Analysis of a Unimorph Nanobeam with a Dielectric Layer of Both Flexoelectricity and Piezoelectricity

In this study, for the first time, free and forced vibrational responses of a unimorph nanobeam consisting of a functionally graded base, along with a dielectric layer of both piezoelectricity and flexoelectricity, is investigated based on paradox-free local/nonlocal elasticity. The formulation and boundary conditions are attained by utilizing the energy method Hamilton’s principle. In order to set a comparison, the formulation of a model in the framework of differential nonlocal is first presented. An effective implementation of the generalized differential quadrature method (GDQM) is then utilized to solve higher-order partial differential equations. This method can be utilized to solve the complex equations whose analytic results are quite difficult to obtain. Lastly, the impact of various parameters is studied to characterize the vibrational behavior of the system. Additionally, the major impact of flexoelectricity compared to piezoelectricity on a small scale is exhibited. The results show that small-scale flexoelectricity, rather than piezoelectricity, is dominant in electromechanical coupling. One of the results that can be mentioned is that the beams with higher nonlocality have the higher voltage and displacement under the same excitation amplitude. The findings can be helpful for further theoretical as well as experimental studies in which dielectric material is used in smart structures.

36 MATERIALS SCIENCE↗
NASA NTRS1986

Chaotic motion in the Jovian atmosphere

Strong nonlinear interactions among unstable waves and the mean flow occur in a simplified quasigeostrophic spectral model of the upper troposphere of Jupiter. The upper boundary of the layer inhibits vertical motion while at the lower boundary perturbations of the potential temperature are not permitted. On an infinite beta plane the forced flow of alternating zones of prograde and retrograde zonal winds, decreasing with height, are linearly unstable and it is shown that the nonlinear terms stabilize the flow by bounding the growth of the eddies. Explicit viscosity terms are not needed. This does not imply that energy would not cascade to the small scale flow but suggests that the nature of the large scale flow is independent of the viscosity at small scales. Numerical time integration shows the flow to be chaotic but, in some cases, with transient propagating features and meandering zonal flow.

Pirraglia, Joseph↗
NASA NTRS1985

Combined aerodynamic and structural dynamic problem emulating routines (CASPER): Theory and implementation

The Combined Aerodynamic and Structural Dynamic Problem Emulating Routines (CASPER) is a collection of data-base modification computer routines that can be used to simulate Navier-Stokes flow through realistic, time-varying internal flow fields. The Navier-Stokes equation used involves calculations in all three dimensions and retains all viscous terms. The only term neglected in the current implementation is gravitation. The solution approach is of an interative, time-marching nature. Calculations are based on Lagrangian aerodynamic elements (aeroelements). It is assumed that the relationships between a particular aeroelement and its five nearest neighbor aeroelements are sufficient to make a valid simulation of Navier-Stokes flow on a small scale and that the collection of all small-scale simulations makes a valid simulation of a large-scale flow. In keeping with these assumptions, it must be noted that CASPER produces an imitation or simulation of Navier-Stokes flow rather than a strict numerical solution of the Navier-Stokes equation. CASPER is written to operate under the Parallel, Asynchronous Executive (PAX), which is described in a separate report.

Jones, William H.↗
NASA NTRS1993

Scaling, soil moisture and evapotranspiration in runoff models

The effects of small-scale heterogeneity in land surface characteristics on the large-scale fluxes of water and energy in the land-atmosphere system has become a central focus of many of the climatology research experiments. The acquisition of high resolution land surface data through remote sensing and intensive land-climatology field experiments (like HAPEX and FIFE) has provided data to investigate the interactions between microscale land-atmosphere interactions and macroscale models. One essential research question is how to account for the small scale heterogeneities and whether 'effective' parameters can be used in the macroscale models. To address this question of scaling, the probability distribution for evaporation is derived which illustrates the conditions for which scaling should work. A correction algorithm that may appropriate for the land parameterization of a GCM is derived using a 2nd order linearization scheme. The performance of the algorithm is evaluated.

Wood, Eric F.↗