Atomization of ethanol jets in a combustor with oscillatory combustion-gas flow
Drop-size distribution data for ethanol jets atomized by downstream injection into combustion gas stream inside rocket combustor
SEARCH · Search NASA
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.
Drop-size distribution data for ethanol jets atomized by downstream injection into combustion gas stream inside rocket combustor
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.
The results of the December 2013 to February 2014 Icing Research Tunnel full icing cloud calibration are presented. The calibration steps included establishing a uniform cloud and conducting drop size and liquid water content calibrations. The goal of the calibration was to develop a uniform cloud, and to generate a transfer function from the inputs of air speed, spray bar atomizing air pressure and water pressure to the outputs of median volumetric drop diameter and liquid water content. This was done for both 14 CFR Parts 25 and 29, Appendix C ('typical' icing) and soon-to-be released Appendix O (supercooled large drop) conditions.
The results of the December 2013 to February 2014 Icing Research Tunnel full icing cloud calibration are being presented to the SAE AC-9C committee, as represented in the 2014 cloud calibration report. The calibration steps included establishing a uniform cloud and conducting drop size and liquid water content calibrations. The goal of the calibration was to develop a uniform cloud, and to generate a transfer function from the inputs of air speed, spray bar atomizing air pressure and water pressure to the outputs of median volumetric drop diameter and liquid water content. This was done for both 14 CFR Parts 25 and 29, Appendix C (typical icing) and soon-to-be released Appendix O (supercooled large drop) conditions.
The justification for describing a spray by treating drop velocity as a random variable on an equal statistical basis with drop size was studied experimentally. A double-exposure technique using fluorescent drop photography was used to make size and velocity measurements at selected locations in a steady ethanol spray formed by a swirl atomizer. The size-velocity data were categorized to construct bivariate spray density functions to describe the spray immediately after formation and during downstream propagation. It was found that a statistical treatment of drop velocity was supported by the data. Spray density function shapes and modal characteristics depended strongly on position and the amount of droplet-gas interaction that had occurred. Bimodal density functions were formed by environmental interaction during downstream propagation. Large differences were also found between spatial mass density and mass flux size distributions at the same location.
This study compares macrophysical and microphysical properties of single-layered, liquid-dominant MBL clouds from the Measurements of Aerosols, Radiation, and Clouds over the Southern Ocean (MARCUS) (above 60°S) and the ARM East North Atlantic (ENA) site during the Aerosol and Cloud Experiments in Eastern North Atlantic (ACE-ENA) field campaign. A total of 1,136 (16.5% of clouds) and 6,034 5-min cloud samples are selected from MARCUS and ARM ENA in this study. MARCUS clouds have higher cloud-top heights, thicker cloud layers, larger liquid water path, and colder cloud temperatures than ENA. Thinner, warmer MBL clouds at ENA can contain higher layer-mean liquid water content due to higher cloud and ocean surface temperatures along with greater precipitable water vapor (PWV). MARCUS has a higher drizzle frequency rate (71.8%) than ENA (45.1%). Retrieved cloud and drizzle microphysical properties from each field campaign show key differences. MARCUS clouds feature smaller cloud droplets, whereas ENA clouds have larger cloud droplets, especially at the upper region of the cloud. From cloud top to cloud base, drizzle drop sizes increase while number concentrations decrease. Drizzle drop radius and number concentration decrease from cloud base to drizzle base due to net evaporation, and MARCUS' lower specific humidity leads to a higher drizzle base than ENA. In conclusion, the broader surface pressure and lower tropospheric stability (LTS) distributions during MARCUS have demonstrated that there are different synoptic patterns for selected cases during MARCUS with less PWV, while ENA is dominated by high pressure systems with nearly doubled PWV.
Non-intrusive and non-invasive optical imaging techniques has generated great interest among researchers for their potential applications to biological study, device characterization, surface defect detection, and jet fuel dynamics. Non-linear optical parametric amplification gate (NLOPG) has been used to detect back-scattered images of objects hidden in diluted Intralipid solutions. To directly detect objects hidden in highly scattering media, the diffusive component of light needs to be sorted out from early arrived ballistic and snake photons. In an optical imaging system, images are collected in transmission or back-scattered geometry. The early arrival photons in the transmission approach, always carry the direct information of the hidden object embedded in the turbid medium. In the back-scattered approach, the result is not so forth coming. In the presence of a scattering host, the first arrival photons in back-scattered approach will be directly photons from the host material. In the presentation, NLOPG was applied to acquire time resolved back-scattered images under the phase matching condition. A time-gated amplified signal was obtained through this NLOPG process. The system's gain was approximately 100 times. The time-gate was achieved through phase matching condition where only coherent photons retain their phase. As a result, the diffusive photons, which were the primary contributor to the background, were removed. With a large dynamic range and high resolution, time-gated early light imaging has the potential for improving rocket/aircraft design by determining jets shape and particle sizes. Refinements to these techniques may enable drop size measurements in the highly scattering, optically dense region of multi-element rocket injectors. These types of measurements should greatly enhance the design of stable, and higher performing rocket engines.
The justification for describing a spray by treating drop velocity as a random variable on an equal statistical basis with drop size was studied experimentally. A double exposure technique using fluorescent drop photography was used to make size and velocity measurements at selected locations in a steady ethanol spray formed by a swirl atomizer. The size velocity data were categorized to construct bivariate spray density functions to describe the spray immediately after formation and during downstream propagation. Bimodal density functions were formed by environmental interaction during downstream propagation. Large differences were also found between spatial mass density and mass flux size distribution at the same location.
The head-on collision of equal sized drops is studied by full numerical simulations. The Navier-Stokes equations are solved for fluid motion both inside and outside the drops using a front tracking/finite difference technique. The drops are accelerated toward each other by a body force that is turned off before the drops collide. When the drops collide, the fluid between them is pushed outward leaving a thin later bounded by the drop surface. This layer gets progressively thinner as the drops continue to deform and in several of the calculations this double layer is artificially removed once it is thin enough, thus modeling rupture. If no rupture takes place, the drops always rebound, but if the film is ruptured the drops may coalesce permanently or coalesce temporarily and then split again.
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.
Recently, a method to retrieve rain drop size distribution (DSD) moments from X-band dual polarization radar measurements has been developed using copolar reflectivity, differential reflectivity and specific attenuation. Two reference moments are retrieved first, followed by reconstructing the DSDs and calculating other moments using a generalized gamma model to represent the underlying shape, h(x), corresponding to the pair of chosen reference moments. Here we explore a similar approach but the retrieval method in this study uses dual-frequency radar measurements. In the case of GPM-DPR, the two frequencies are 13.8and 35 GHz, and their products include attenuation corrected reflecivities at both frequencies (Z and Z) and the specific attenuation (k, k, if available). Our approach is to use these products to determine two reference moments, namely M3 andM6 representing the third and the sixth moments respectively. Then, as with the polarimetric radar retrievals, we use the most probable h(x) to reconstruct the full DSD spectra, from which other moments are calculated. The best two DPR products for estimating M3 and M6 appear to be A and Z. Figure 1 shows the retrievals versus the ‘true’ moments. 2930 three-minute DSDs were used for the (T-matrix) scattering calculations at Ku and Ka bands but only cases with A > 0.5 dB/km were chosen. For the retrieved moments, the integration was performed only up to 6 mm drop diameter. The [1:1] line is included in Fig. 1. Even with just two DPR products as input the results seem promising, although the lower order moments show somewhat more scatter, especially the zeroth moment, M0. We will quantify the retrieval errors, and additionally examine the stability of h(x). For the latter, data from (i) Greeley, Colorado, (ii) Huntsville, Alabama, and (iii)Wallops, Virginia, will be used. Finally, a GPM overpass case over Huntsville, during a widespread rain event on 11 April 2016, will be considered as an initial test case.
The Vaisala FD70 is a multi-parameter present weather and visibility sensor designed to measure precipitation type, intensity, and visibility with high accuracy in diverse environmental conditions. It uses a combination of forward-scatter measurement and optical disdrometer technologies to detect drop size, fall speeds, and optical properties, enabling the classification of various precipitation types such as rain, snow, sleet, and freezing rain along is visibility estimates. The FD70 provides quantitative estimates of liquid-equivalent precipitation rate and meteorological optical range (MOR), supporting applications in meteorological research, aviation, and road weather monitoring. These measurements are collected at the Argonne Testbed for Multiscale Observational Science (ATMOS), a 20 acre prairie site at Argonne National Lab, located in Lemont, IL. Data is available in netcdf format. Each file contains one second interval data, for approximately 24 hrs each day. File naming convention includes the project (CROCUS), location (ATMOS), instrument name, data level (raw, a1), and date (year, month, day).
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
Detonation wave propagation, drop size effects on velocity, and streak schlieren photography on drop shattering in liquid fuel-oxygen mixture
The mechanism for initiating drizzle drop remains a gap in the current understanding of warm rain formation. One prevalent hypothesis suggests that the presence of Giant Cloud Condensation Nuclei (GCCN) generates drizzle‐sized drops necessary to trigger the Collision‐Coalescence (C‐C) process. Here, in this study, this hypothesis is investigated using a novel framework that integrates in situ observations, remote sensing measurements, and idealized models. Results show that GCCN can efficiently generate drizzle drops through condensation, producing a broad Droplet Size Distribution (DSD) comparable to in situ observations. The large drizzle drop and broad DSD strongly facilitate C‐C, further accelerating drizzle initiation. To compare with observation, the model‐generated DSDs are used to generate radar Doppler spectra where radar reflectivity and Doppler skewness is estimated. The simulated radar quantities correspond well with radar observations, providing critical evidence for the GCCN‐induced drizzle initiation mechanism.
Subsonic, external flow about nonlifting bodies, lifting bodies or combinations of lifting and nonlifting bodies is calculated by a modified version of the Hess lifting code. Trajectory calculations can be performed for any atmospheric conditions and for all water drop sizes, from the smallest cloud droplet to large raindrops. Experimental water drop drag relations are used in the water drop equations of motion and effects of gravity settling are included. Inlet flow can be accommodated, and high Mach number compressibility effects are corrected for approximately. Seven codes are described: (1) a code used to debug and plot body surface description data; (2) a code that processes the body surface data to yield the potential flow field; (3) a code that computes flow velocities at arrays of points in space; (4) a code that computes water drop trajectories from an array of points in space; (5) a code that computes water drop trajectories and fluxes to arbitrary target points; (6) a code that computes water drop trajectories tangent to the body; and (7) a code that produces stereo pair plots which include both the body and trajectories. Accuracy of the calculations is discussed, and trajectory calculation results are compared with prior calculations and with experimental data.
This paper presents experimental work conducted in the Icing Research Tunnel at NASA Glenn Research Center to characterize the velocity of large drops in the test section. Some icing spray clouds with large drops were generated with Mod1 nozzles at low nozzle air pressure of 2 to 4 psig for various tunnel air speeds. Drop diameters and drop velocities were measured via high-resolution imaging with a Particle Imaging Particle Tracking Velocimetry probe developed by Artium Technologies. The probe was mounted at four different locations aligned with the centerline of the test section from near the end of the contraction to the constant height test section part of the tunnel. CFD analyses were performed. It showed that the probe head geometry affects the local air flow in the measurement area of the probe between the prongs and in front of the probe. Initial analysis of the air velocity data during the test also indicated that the probe mounting stand has blockage effect on the local tunnel air velocity measurement by a pitot-static probe affixed on the mounting plate next to the stand. Those findings were later verified in the Icing Research Tunnel using a new pitot-static probe design with a linear motion system to measure the local tunnel air velocity with and without the probe. As a result, additional drop trajectory simulations were performed with the airflow moving towards the probe head with a Langmuir-D 7-bin distribution cloud containing large drops. The simulation results helped identify a critical drop-size threshold of 300 μm above which the velocities of larger drops are negligibly affected by the adverse pressure gradient generated by the probe head due to their large drop inertia. From the dimensional analysis of the drop velocity measurement data obtained, it showed that at the tunnel test section reference location a generalized empirical correlation was developed for the non-dimensional drop velocity as a function of the non-dimensional corrected drop diameter independent of the actual spraybar pressure settings and the tunnel air speeds. The generalized curve-fit correlation showed that the drop velocity was universally asymptotic to about 86 percent of the corresponding tunnel air speed at the test section reference location for the largest drop diameter captured by the probe. Further evaluation of this correlation is recommended to assess its applicability for Supercooled Large Drop icing scaling applications in the Icing Research Tunnel.
This study uses observations from a ground-based instrument suite to investigate the rain-snow level in stratiform rainfall from January 2014 to April 2020 in the Upper Great Lakes Region. The height above the surface where ice melts to rain, the rain-snow level (RSL), influences microphysical assumptions in remote sensing precipitation retrievals and the ability of space-based radar to discriminate surface precipitation phase because of ground clutter. The instrument suite is installed at the Marquette, MI (MQT) National Weather Service station adjacent to Lake Superior. Rain events and the RSL are studied through a ground-based vertically profiling radar (Micro Rain Radar), a custom NASA-developed video disdrometer (Precipitation Imaging Package), and reanalyses products from ECMWF ERA5 and NASA MERRA-2. Further, distinct macro and micro physical characteristics are observed in precipitation events with shallow RSL ( < 1.8 km above ground level [AGL]) and intermediate RSL ( > 1.8 km AGL). Intermediate RSL correspond to rain events with relatively higher rain rates and a higher concentration of small drops in the drop size distributions (DSDs). Shallow RSL DSDs contain relatively higher concentrations of large drops with lower fall speeds suggesting that partially melted snowflakes may be reaching the surface. Reflectivity-rain-rate relationships are also impacted by microphysical differences associated with RSL regimes. Radar-detected RSL agree with reanalysis-derived melt levels- especially at wet-bulb isotherms of 0.5°C and 1°C. Seasonal differences such as shallow rain-snow levels in winter, fall, and spring have subsequent implications for satellite detectability.