Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Drop size distribution”

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 109 records · Page 6

Comparison of Bimodal Ice Shapes on Full-Chord and Hybrid CRM Midspan Models

An ice shape database has been created to document ice accretions on a full-chord length and hybrid versions of a 65 percent scale size of the mid-span section of the Common Research Model resulting from an exposure to a Super-cooled Large Drop (SLD) icing cloud with a bimodal drop size distribution. The ice shapes created were documented with photographs, laser scanned surface measurements over a section of the model span, and measurement of the ice mass over the same section of each ice accretion. Drop distribution effects were evaluated by using the same IRT icing conditions except with either monomodal or bimodal drop distributions. Ice shapes resulting from the bimodal distribution as well as from monomodal distributions of identical cloud droplet Median Volume Diameter (MVD) were obtained and compared. Results indicated that the ice shapes resulting from the monomodal and bimodal drop size distributions have similar 2D Maximum Combined Cross Section, or MCCS, ice shape around the model leading edge region and equivalent ice mass measurement per unit span to within the IRT SLD cloud calibration uncertainty, as expected with this approach, but the bimodal distributions have distinctly larger ice feathers that extended further back on the chord of the model. For most test conditions, the overall size and shape of the main ice accretions compared well between the full-chord and hybrid models

Icing , Simulation, Validation↗

Comparison of Bimodal Ice Shapes on Full-Chord and Hybrid CRM Midspan Models

An ice shape database has been created to document ice accretions on a full-chord length and hybrid versions of a 65 percent scale size of the mid-span section of the Common Research Model resulting from an exposure to a Super-cooled Large Drop (SLD) icing cloud with a bimodal drop size distribution. The ice shapes created were documented with photographs, laser scanned surface measurements over a section of the model span, and measurement of the ice mass over the same section of each ice accretion. Drop distribution effects were evaluated by using the same IRT icing conditions except with either monomodal or bimodal drop distributions. Ice shapes resulting from the bimodal distribution as well as from monomodal distributions of identical cloud droplet Median Volume Diameter (MVD) were obtained and compared. Results indicated that the ice shapes resulting from the monomodal and bimodal drop size distributions have similar 2D Maximum Combined Cross Section, or MCCS, ice shape around the model leading edge region and equivalent ice mass measurement per unit span to within the IRT SLD cloud calibration uncertainty, as expected with this approach, but the bimodal distributions have distinctly larger ice feathers that extended further back on the chord of the model. For most test conditions, the overall size and shape of the main ice accretions compared well between the full-chord and hybrid models.

Icing↗

Laser Calibration of an Impact Disdrometer

A practical approach to developing an operational low-cost disdrometer hinges on implementing an effective in situ adaptive calibration strategy. This calibration strategy lowers the cost of the device and provides a method to guarantee continued automatic calibration. In previous work, a collocated tipping bucket rain gauge was utilized to provide a calibration signal to the disdrometer's digital signal processing software. Rainfall rate is proportional to the 11/3 moment of the drop size distribution (a 7/2 moment can also be assumed, depending on the choice of terminal velocity relationship). In the previous case, the disdrometer calibration was characterized and weighted to the 11/3 moment of the drop size distribution (DSD). Optical extinction by rainfall is proportional to the 2nd moment of the DSD. Using visible laser light as a means to focus and generate an auxiliary calibration signal, the adaptive calibration processing is significantly improved.

calibration↗

Microwave scattering from wind- and rain-roughened seas

Microwave scattering from wind- and rain-roughened seas is studied under the assumption that, for light wind and light rain, the processes contributing to the scatterometer signal are independent. It is argued that calibrations from exclusively rain and exclusively wind conditions are all that is needed to predict scatterometer cross section for combined conditions. It is demonstrated that scatterometer response to a rain-roughened water surface is linear, but the slope of the relationship is sensitive to drop-size distribution. It is concluded that first-order wind and rain can be regarded as independent processes. Average cross section varies linearly with rainfall rate, but the slope of the relationship is a function of drop-size distribution. For combined light rain plus light wind, sigma-c is found to be well predicted by the addition of wind and rain models.

Bliven, Larry F.↗

Microphysical Retrievals Over Stratiform Rain Using Measurements from an Airborne Dual-Wavelength Radar-Radiometer

The need to understand the complementarity of the radar and radiometer is important not only to the Tropical Rain Measuring Mission (TRMM) program but to a growing number of multi-instrumented airborne experiment that combine single or dual-frequency radars with multichannel radiometers. The method of analysis used in this study begins with the derivation of dual-wavelength radar equations for the estimation of a two-parameter drop size distribution (DSD). Defining a "storm model" as the set of parameters that characterize snow density, cloud water, water vapor, and features of the melting layer, then to each storm model there will usually correspond a set of range-profiled drop size distributions that are approximate solutions of the radar equations. To test these solutions, a radiative transfer model is used to compute the brightness temperatures for the radiometric frequencies of interest. A storm model or class of storm models is considered optimum if it provides the best reproduction of the radar and radiometer measurements. Tests of the method are made for stratiform rain using simulated storm models as well as measured airborne data. Preliminary results show that the best correspondence between the measured and estimated radar profiles usually can be obtained by using a moderate snow density (0.1-0.2 g/cu cm), the Maxwell-Garnett mixing formula for partially melted hydrometeors (water matrix with snow inclusions), and low to moderate values of the integrated cloud liquid water (less than 1 kg/sq m). The storm-model parameters that yield the best reproductions of the measured radar reflectivity factors also provide brightness temperatures at 10 GHz that agree well with the measurements. On the other hand, the correspondence between the measured and modeled values usually worsens in going to the higher frequency channels at 19 and 34 GHz. In searching for possible reasons for the discrepancies, It is found that changes in the DSD parameter Mu, the radar constants, or the path-integrated attenuation can affect the high frequency channels significantly. In particular, parameters that cause only modest increases in the median mass diameter of the snow, and which have a minor effect on the radar returns or the low frequency brightness temperature, can produce a strong cooling of the 34 GHz brightness temperature.

Meneghini, Robert↗

DSD Characteristics of a Mid-Winter Tornadic Storm Using C-Band Polarimetric Radar and Two 2D-Video Disdrometers

Drop size distributions in an evolving tornadic storm are examined using C-band polarimetric radar observations and two 2D-video disdrometers. The E-F2 storm occurred in mid-winter (21 January 2010) in northern Alabama, USA, and caused widespread damage. The evolution of the storm occurred within the C-band radar coverage and moreover, several minutes prior to touch down, the storm passed over a site where several disdrometers including two 2D video disdrometers (2DVD) had been installed. One of the 2DVDs is a low profile unit and the other is a new next generation compact unit currently undergoing performance evaluation. Analyses of the radar data indicate that the main region of precipitation should be treated as a "big-drop" regime case. Even the measured differential reflectivity values (i.e. without attenuation correction) were as high as 6-7 dB within regions of high reflectivity. Standard attenuation-correction methods using differential propagation phase have been "fine tuned" to be applicable to the "big drop" regime. The corrected reflectivity and differential reflectivity data are combined with the co-polar correlation coefficient and specific differential phase to determine the mass-weighted mean diameter, Dm, and the width of the mass spectrum, (sigma)M, as well as the intercept parameter , Nw. Significant areas of high Dm (3-4 mm) were retrieved within the main precipitation areas of the tornadic storm. The "big drop" regime assumption is substantiated by the two sets of 2DVD measurements. The Dm values calculated from 1-minute drop size distributions reached nearly 4 mm, whilst the maximum drop diameters were over 6 mm. The fall velocity measurements from the 2DVD indicate almost all hydrometeors to be fully melted at ground level. Drop shapes for this event are also being investigated from the 2DVD camera data.

Thurai, M.↗

On the Use of the Log-Normal Particle Size Distribution to Characterize Global Rain

Although most parameterizations of the drop size distributions (DSD) use the gamma function, there are several advantages to the log-normal form, particularly if we want to characterize the large scale space-time variability of the DSD and rain rate. The advantages of the distribution are twofold: the logarithm of any moment can be expressed as a linear combination of the individual parameters of the distribution; the parameters of the distribution are approximately normally distributed. Since all radar and rainfall-related parameters can be written approximately as a moment of the DSD, the first property allows us to express the logarithm of any radar/rainfall variable as a linear combination of the individual DSD parameters. Another consequence is that any power law relationship between rain rate, reflectivity factor, specific attenuation or water content can be expressed in terms of the covariance matrix of the DSD parameters. The joint-normal property of the DSD parameters has applications to the description of the space-time variation of rainfall in the sense that any radar-rainfall quantity can be specified by the covariance matrix associated with the DSD parameters at two arbitrary space-time points. As such, the parameterization provides a means by which we can use the spaceborne radar-derived DSD parameters to specify in part the covariance matrices globally. However, since satellite observations have coarse temporal sampling, the specification of the temporal covariance must be derived from ancillary measurements and models. Work is presently underway to determine whether the use of instantaneous rain rate data from the TRMM Precipitation Radar can provide good estimates of the spatial correlation in rain rate from data collected in 5(sup 0)x 5(sup 0) x 1 month space-time boxes. To characterize the temporal characteristics of the DSD parameters, disdrometer data are being used from the Wallops Flight Facility site where as many as 4 disdrometers have been used to acquire data over a 2 km path. These data should help quantify the temporal form of the covariance matrix at this site.

Meneghini, Robert↗

The meteorological parameterization of specific attenuation in rain viewed at Nadir

Polynomial expressions are presented for the parameterization of the specific attenuation in rain from 9 to 38 GHz that are applicable to a wide range of naturally occurring drop size distributions and for viewing angles close to nadir. Because the temperature T affects the specific attenuation at some frequencies, expressions for the polynomial coefficients as functions of T are also provided for -10C less than or equal to T less than or equal to 30C. The advantage of this parameterization is that even without a detailed specification of the drop size distribution, useful estimates of the specific attenuation are often possible given only two out of three parameters, namely, the rainfall rate R, the rainwater content W, or D(sub m) (the mass-weighted mean drop diameter), particularly if the temperature is also specified.

Jameson, A. R.↗

Automated holographic drop-size analyzer

System analyzes drop-size distribution in liquid-droplet-spray combustion fields. Holographic camera takes "stop-motion" hologram of combustion volume; it is then viewed by vidicon camera connected to digital data-processing system that identifies particles or droplets, determining their size and count, and displays histogram of drop-size distribution in holographic field.

Feinstein, S. P.↗

Optimal Estimation Retrievals and Their Uncertainties: What Every Atmospheric Scientist Should Know

Remote sensing instruments are heavily used to provide observations for both the operational and research communities. These sensors do not provide direct observations of the desired atmospheric variables, but instead, retrieval algorithms are necessary to convert the indirect observations into the variable of interest. It is critical to be aware of the underlying assumptions made by many retrieval algorithms, including that the retrieval problem is often ill posed and that there are various sources of uncertainty that need to be treated properly. In short, the retrieval challenge is to invert a set of noisy observations to obtain estimates of atmospheric quantities. The problem is often complicated by imperfect forward models, by imperfect prior knowledge, and by the existence of nonunique solutions. Optimal estimation (OE) is a widely used physical retrieval method that combines measurements, prior information, and the corresponding uncertainties based on Bayes’s theorem to find an optimal solution for the atmospheric state. Furthermore, OE also allows the relative contributions of the different sources of error to the uncertainty in the final retrieved atmospheric state to be understood. Here, we provide a novel Python library to illustrate the use of OE for inverse problems in the atmospheric sciences. We introduce two example problems: how to retrieve drop size distribution parameters from radar observations and how to retrieve the temperature profile from ground-based microwave sensors. Using these examples, we discuss common pitfalls, how the various error sources impact the retrieval, and how the quality of the retrieval results can be quantified.

54 ENVIRONMENTAL SCIENCES↗

Application of digital image analysis techniques to antimisting fuel spray characterization

Pulsed ruby laser sheet illumination of the spray is used for the initial recording of data on very-high-resolution photographic film. The digitization of mosaic elements is effected with a vidicon and video digitizer whose output is stored in computer RAM memory for processing. Highly nonspherical elements and a broad range of drop diameters (8-2000 microns) resulting from the unusual rheological properties of the fuel-additive system are accommodated by the device configuration and algorithms. It is found that the generation of two-dimensional images by means of scattered light also eliminates errors resulting from variations in the index of refraction and from the submicron scattering sites that are often present within the modified fuel. No a priori information on the drop size distribution or on the system response to various drop sizes is required.

Fleeter, R.↗

Phase segregation due to simultaneous migration and coalescence

Ground-based modeling and experiments have been performed on the interaction and coalescence of drops leading to macroscopic phase separation. The focus has been on gravity-induced motion, with research also initiated on thermocapillary motion of drops. The drop size distribution initially shifts toward larger drops with time due to coalescence, and then a back towards smaller drops due to the larger preferentially settling out. As a consequence, the phase separation rate initially increases with time and then decreases.

Davis, Robert H.↗

Identifying and tracking bubbles and drops in simulations: A toolbox for obtaining sizes, lineages, and breakup and coalescence statistics

Knowledge of bubble and drop size distributions in two-phase flows is important for characterizing a wide range of phenomena, including combustor ignition, sonar communication, and cloud formation. The physical mechanisms driving the background flow also drive the time evolution of these distributions. Accurate and robust identification and tracking algorithms for the dispersed phase are necessary to reliably measure this evolution and thereby quantify the underlying mechanisms in interface-resolving flow simulations. The identification of individual bubbles and drops traditionally relies on an algorithm used to identify connected regions. This traditional algorithm can be sensitive to the presence of spurious structures. A cost-effective refinement is proposed to maximize volume accuracy while minimizing the identification of spurious bubbles and drops. An accurate identification scheme is crucial for distinguishing bubble and drop pairs with large size ratios. The identified bubbles and drops need to be tracked in time to obtain breakup and coalescence statistics that characterize the evolution of the size distribution, including breakup and coalescence frequencies, and the probability distributions of parent and child bubble and drop sizes. An algorithm based on mass conservation is proposed to construct bubble and drop lineages using simulation snapshots that are not necessarily from consecutive time steps. These lineages are then used to detect breakup and coalescence events, and obtain the desired statistics. Accurate identification of large-size-ratio bubble and drop pairs enables accurate detection of breakup and coalescence events over a large size range. Accurate detection of successive breakup and coalescence events requires that the snapshot interval be an order of magnitude smaller than the characteristic breakup and coalescence times to capture these successive events while minimizing the identification of repeated confounding events. Together, these algorithms serve as a toolbox for detailed analysis of two-phase simulations, and enable insights into the mechanisms behind bubble and drop formation and evolution in flows of practical importance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Evolution of Droplet Size Distributions During the Transition of an Ultraclean Stratocumulus Cloud System to Open Cell Structure: An LES Investigation Using Lagrangian Microphysics

Here, a state-of-the-art Lagrangian microphysics scheme is used in a large-eddy simulation to investigate the stratocumulus transition from closed to open cell structure. Processes controlling precipitation development, which is a key to the transition, are analyzed by leveraging unique benefits of Lagrangian microphysics, particularly the ability to track computational drops in the flow. Sufficient time is needed for coalescence growth of cloud drops to drizzle within the updraft-downdraft cycle of large eddies. This favors broad drop size distributions (DSDs) and drizzle growth in downdrafts, where drops are typically much older than in updrafts. During the closed cell stage, mean cloud drop radius is too small, and the DSDs are too narrow, so that the timescale for coalescence is much longer than the large eddy turnover time and drizzle growth is limited. The closed-to-open cell transition occurs when these timescales become comparable and the precipitation flux increases sharply.

54 ENVIRONMENTAL SCIENCES↗

Rain Fallspeeds and Rates Derived from Airborne Nadir-Pointing Doppler Radar Measurements

The use of vertical-incidence Doppler velocity in addition to radar reflectivity may yield information on drop size distribution and therefore result in better rainrate estimates. Doppler velocity can provide useful information on the raindrop size distribution. Doppler velocities from a zenith-pointing radar represent the sum of the mean reflectivity-weighted hydrometeor fallspeed and the vertical air motion. Dual-parameter rain estimation methods using the Doppler velocity, require that the latter can be removed, or is negligible. Atlas et al. (1972) derived relations between Doppler velocity, reflectivity, and rain rate assuming an exponential size distribution for rain. Ulbrich (1994) expanded on this work by deriving the relation between the Doppler velocity and the reflectivity assuming a Gamma size distribution. This distribution provides a more realistic representation of the small rain drops. To get accurate information on raindrop size distributions with the above method, the air motions must be removed from the observed Doppler velocities

Heymsfield, Gerald M.↗

An alternative approach to estimating rainfall rate by radar using propagation differential phase shift

In this work it is shown that for frequencies from 3 to 13 GHz, the ratio of the specific propagation differential phase shift phi(sub DP) to the rainfall rate can be specified essentially independently of the form of the drop size distribution by a function only of the mass-weighted mean drop size D(sub m). This significantly reduces one source of substantial bias errors common to most other techniques for measuring rain by radar. For frequencies 9 GHz and greater, the coefficient can be well estimated from the ratio of the specific differential attenuation to phi(sub DP), while at nonattenuating frequencies such as 3 GHz, the coefficient can be well estimated using the differential reflectivity. In practice it appears that this approach yields better estimates of the rainfall rate than any other current technique. The best results are most likely at 13.80 GHz, followed by those at 2.80 GHz. An optimum radar system for measuring rain should probably include components at a both frequencies so that when signals at 13.8 GHz are lost because of attenuation, good measurements are still possible at the lower frequency.

Jameson, A. R.↗

Estimating the Concentration of Large Raindrops from Polarimetric Radar and Disdrometer Observations

Estimation of rainfall integral parameters, including radar observables, and empirical relations between them are sensitive to the truncation of the drop size distribution (DSD), particularly at the large drop end. The sensitivity of rainfall integral parameters to the maximum drop diameter (D(sub max)) is exacerbated at C‐band since resonance effects are pronounced for large drops in excess of 5 mm diameter (D). Due to sampling limitations, it is often difficult to reliably estimate D(sub max) with disdrometers. The resulting uncertainties in D(sub max0 potentially increase errors in radar retrieval methods, particularly at C‐band, that rely on disdrometer observations for DSD input to radar models. In fact, D(sub max) is typically an assumed DSD parameter in the development of radar retrieval methods. Because of these very uncertainties, it is difficult to independently confirm disdrometer estimates of D(sub max) with polarimetric radar observations. A couple of approaches can be taken to reduce uncertainty in large drop measurement. Longer integration times can be used for the collection of larger disdrometer samples. However, integration periods must be consistent with a radar resolution volume (RRV) and the temporal and spatial scales of the physical processes affecting the DSD therein. Multiple co‐located disdrometers can be combined into a network to increase the sample size within a RRV. However, over a reasonable integration period, a single disdrometer sample volume is many orders of magnitudes less than a RRV so it is not practical to devise a network of disdrometers that has an equivalent volume to a typical RRV. Since knowledge of DSD heterogeneity and large drop occurrence in time and space is lacking, the specific accuracy or even general representativeness of disdrometer based D(sub max) and large drop concentration estimates within a RRV are currently unknown. To address this complex issue, we begin with a simpler question. Is the frequency of occurrence of large rain drops (D > 5 mm) in disdrometer observations, either stand alone or networked, generally representative and consistent with polarimetric radar observations? We first show from simulations that the concentration of large (D > 5 mm) rain drops (N(sub T5)) can be estimated from polarimetric observations of specific differential phase (K(sub dp)) and differential reflectivity (Z(sub dr)), N(sub T5)=F(K(sub dp),Z(sub dr)), or horizontal reflectivity (Z(sub h)) and Z(sub dr), N(sub T5)=(Z(sub h),Z(sub dr)). We assess the error associated with polarimetric retrieval of N(sub T5), including sensitivity to D(sub max) parameterization assumptions and measurement error in the radar simulations. Polarimetric measurements at S‐band and C‐band will then be used to retrieve estimates of N(sub T5) and compared to disdrometer estimates of N(sub T5). After careful consideration of retrieval error, we will check consistency between disdrometer and polarimetric radar estimates of N(sub T5) and the frequency of occurrence of large rain drops in a variety of precipitating regimes using data from NASA's Global Precipitation Measurement (GPM) Ground Validation (GV) program, including field campaigns such as MC3E (Oklahoma) and IFloodS (Iowa) and extended measurements over Huntsville, Alabama and NASA Wallops Flight Facility in coastal Virginia.

Carey, Lawrence D.↗