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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Evaluating Entrainment–Mixing Characteristics through Direct Comparisons of Drop Size Distributions Using In Situ Observations from ACE-ENA

Abstract Constraining the impacts of entrainment and associated mixing (i.e., entrainment–mixing) on cloud properties continues to be difficult, partly due to observational uncertainties as well as a lacking consensus of which methodologies for diagnosing entrainment–mixing are most appropriate. This study introduces a novel method to evaluate the presence and degree of inhomogeneous and homogeneous mixing using ∼100 h of in situ observations from a research aircraft over the northeastern Atlantic. Specifically, drop size distributions are compared between regions containing negligible and significant entrainment for select flight legs, making a direct characterization of the degree of homogeneous and inhomogeneous mixing possible. A measure of drop concentration variance is used as a proxy variable to diagnose entrainment–mixing. Results correspond well with entrainment–mixing metrics, showing lower Damköhler numbers where drop size distributions shift toward smaller drop sizes (i.e., inhomogeneous mixing) and greater transition length scales where drop size distributions do not (i.e., homogeneous mixing). Inhomogeneous mixing occurs in most samples from Aerosol and Cloud Experiment in the Eastern North Atlantic (ACE-ENA) (regardless of homogeneous mixing frequencies increasing with increasing spatial resolution from ∼100 to ∼10 m) and is associated with decreased drop size relative dispersion and both greater aerosol and drop concentrations compared with homogeneous mixing. Precipitating clouds have a greater frequency of homogeneous mixing compared with nonprecipitating clouds. The proposed methodology is similarly applied to in situ observations of southeast Pacific stratocumulus, shallow convective clouds over central Oklahoma and low-level clouds over the Southern Ocean. All four locations are primarily dominated by inhomogeneous mixing with minimal variability among each region.

Clouds↗

DOE OSTI2020

Drop size distribution broadening mechanisms in a bin microphysics Eulerian model

In this study, processes that broaden drop size distributions (DSDs) in Eulerian models with two-moment bin microphysics are analyzed. Numerous tests are performed to isolate the effects of different physical mechanisms that broaden DSDs in two- and three-dimensional Weather Research and Forecasting model simulations of an idealized ice-free cumulus cloud. Additionally, sensitivity of these effects to modifying horizontal and vertical model grid spacings is also examined. As expected, collision-coalescence is a key process broadening the modeled DSDs. In-cloud droplet activation also contributes substantially to DSD broadening, whereas evaporation has only a minor effect and sedimentation has little effect. Cloud dilution (mixing of cloud-free and cloudy air) also broadens the DSDs considerably, whether or not it is accompanied by evaporation. This mechanism involves the reduction of droplet concentration from dilution along the cloud’s lateral edges, leading to locally high supersaturation and enhanced drop growth when this air is subsequently lifted in the updraft. DSD broadening ensues when the DSDs are mixed with those from the cloud core. Decreasing the horizontal and vertical model grid spacings from 100 m to 30 m has limited impact on the DSDs. However, when these physical broadening mechanisms (in-cloud activation, collision-coalescence, dilution, etc.) are turned off, there is a reduction of DSD width by up to ~20-50% when the vertical grid spacing is decreased from 100 m to 30 m, consistent with effects of artificial broadening from vertical numerical diffusion. Nonetheless, this artificial numerical broadening appears to be relatively unimportant overall for DSD broadening when physically based broadening mechanisms in the model are included for this cumulus case.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Rain Drop Size Distributions Estimated from NOAA Snow-Level Radar Data

Using NOAA’s S-band High-Power Snow-Level Radar (HPSLR), a technique for estimating the rain drop size distribution (DSD) above the radar is presented. This technique assumes the DSD can be described by a four parameter, generalized gamma distribution (GGD). Using the radar’s measured average Doppler velocity spectrum and a value (assumed, measured, or estimated) of the vertical air motion w, an estimate of the GGD is obtained. Four different methods can be used to obtain w. One method that estimates a mean mass-weighted raindrop diameter D m from the measured reflectivity Z produces realistic DSDs compared to prior literature examples. These estimated DSDs provide evidence that the radar can retrieve the smaller drop sizes constituting the “drizzle” mode part of the DSD. Here, this estimation technique was applied to 19 h of observations from Hankins, North Carolina. Results support the concept that DSDs can be modeled using GGDs with a limited range of parameters. Further work is needed to validate the described technique for estimating DSDs in more varied precipitation types and to verify the vertical air motion estimates.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Phase Doppler Interferometry for Efficient Cloud Drop Size Distribution, Number Density, and LWC Measurements

Threats to aviation safety as a result of super-cooled large drops (SLD) has been addressed by the FAA rules change (14 CFR Part 25) with the additional icing certification requirement. SLD clouds often consist of bi-modal drop size spectra leading to significant problems in simulating and characterizing these conditions in situ and in icing wind tunnels. Legacy instrumentation for measuring drop size distributions and liquid water content are challenged under these conditions. The large size range measurement problem is addressed with the development of the Phase Doppler Interferometer, Flight Probe Dual-Range (PDI FPDR). The method is described in this report along with the measurement capabilities including the dynamic measurement range and overall working size range. The PDI instrument bases drop size measurements on the light wavelength as the measurement length scale. The light wavelength is a much more robust scale, especially as compared to the light scattering intensity. Additionally, methods for accurately characterizing the sample volume in situ based on measured drop velocity and transit time are reviewed, given the importance of this parameter for merging results and measuring LWC. Droplet coincidence in the sample volume can be problematic so this condition is treated with an innovative signal parsing approach. Measurement examples acquired in the NASA IRT are provided. Measurements of LWC showed good agreement with the Artium Particle Imaging (PI) instrument but diverged from the tunnel calibration results for larger MVD values.

42 ENGINEERING↗

DOE OSTI2021

Drop Size Distribution Variability in Central Argentina during RELAMPAGO-CACTI

The Remote sensing of Electrification, Lightning, And Meso-scale/micro-scale Processes with Adaptive Ground Observations (RELAMPAGO) and the Cloud, Aerosol, and Complex Terrain Interactions Experiment Proposal (CACTI) field campaigns provided an unprecedented thirteen-disdrometer dataset in Central Argentina during the Intensive (IOP, 15 November to 15 December 2018) and Extended (EOP, 15 October 2018 to 30 April 2019) Observational Periods. The drop size distribution (DSD) parameters and their variability were analyzed across the region of interest, which was divided into three subregions characterized by the differing proximity to the Sierras de Córdoba (SDC), in order to assess the impact of complex terrain on the DSD parameters. A rigorous quality control of the data was first performed. The frequency distributions of DSD-derived parameters were analyzed, including the normalized intercept parameter (logN w ), the mean volume diameter (D 0 ), the mean mass diameter (D m ), the shape parameter (μ), the liquid water content (LWC), and the rain rate (R). The region closest to the SDC presented higher values of logN w , lower D 0 , and higher μ, while the opposite occurred in the farthest region, i.e., the concentration of small drops decreased while the concentration of bigger drops increased with the distance to the east of the SDC. Furthermore, the region closest to the SDC showed a bimodal distribution of D 0 : the lower values of D 0 were associated with higher values of logN w and were found more frequently during the afternoon, while the higher D 0 were associated with lower logN w and occurred more frequently during the night. The data were analyzed in comparison to the statistical analysis of Dolan et al. 2018 and sorted according to the classification proposed in the cited study. The logN w -D 0 and LWC-D 0 two-dimensional distributions allowed further discussion around the applicability of other mid-latitude and global precipitation classification schemes (startiform/convection) in the region of interest. Finally, three precipitation case studies were analyzed with supporting polarimetric radar data in order to relate the DSD characteristics to the precipitation type and the microphysical processes involved in each case.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Observed and Bin Model Simulated Evolution of Drop Size Distributions in High-Based Cumulus Congestus Over the United Arab Emirates

This study examines microphysical processes in developing high-based cumulus congestus over the United Arab Emirates using aircraft observations and a large-eddy-simulation model with bin microphysics. A notable feature of this case is the lack of mm-sized drops despite having a liquid cloud layer >3.5 km deep, contrasting with copious large drops observed in maritime tropical cumulus congestus having a similar vertical extent. Modeled drop size distributions are similar to observations at various temperatures between 9.5 and –12°C, including the lack of mm-sized drops. Cloud dilution leads to low-to-moderate liquid water contents (~0.5–1.5 g m –3 ) in most of the cloud core, several times smaller than adiabatic values. Dilution is enhanced in the inflowing branch of the toroidal circulations associated with individual cloud thermals, which are favored regions for secondary droplet activation (activation above cloud base). Secondary activation in general contributes substantially to the droplet population. Turning it off leads to a sharp decrease in droplet concentration and increase in mean size aloft, but does little to increase rain drop production. Warm rain generation (or lack thereof) in this case is therefore determined more by the sub-cloud aerosol and the cloud base droplet size distribution (DSD) than DSD evolution aloft from secondary droplet activation. Decreasing the aerosol concentration by a factor of 10 greatly increases production of large drops via collision-coalescence. Thus, despite its high base (low temperatures) and substantial dilution, the simulated cloud is thermodynamically and dynamically capable of rapidly producing copious mm-sized drops from collision-coalescence under pristine aerosol conditions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Warm-phase microphysical evolution in large-eddy simulations of tropical cumulus congestus: evaluating drop size distribution evolution using polarimetry retrievals, in situ measurements, and a thermal-based framework

Owing to uncertainties in convective microphysics processes, improving parameterizations in Earth system models (ESMs) can benefit from observationally constrained cases suitable for scaling between cloud-resolving models and ESMs. We propose a benchmark large-eddy simulation (LES) cumulus congestus case study from the NASA Cloud, Aerosol, and Monsoon Processes Philippines Experiment (CAMP 2 Ex) for evaluating and improving ESMs in single-column model (SCM) mode. We seek observational constraints using novel polarimetric retrievals and in situ cloud microphysics measurements. Simulations using bulk and bin microphysics initialized with observed aerosol profiles are compared to cloud-top retrievals of cloud droplet effective radius (R eff ), effective variance (ν eff ), and number concentration (N d ) from the airborne Research Scanning Polarimeter (RSP). Both schemes reproduce characteristics of cloud-top N d and R eff that increase and decrease with altitude, respectively. Cloud-top N d is low-biased relative to RSP retrievals in both schemes, potentially due to limitations in both simulations and retrieval assumptions. Cloud-top R eff is low-biased in the bulk scheme but reasonably reproduced by the bin scheme. Profiles of N d and R eff are sensitive to the collision–coalescence process and the vertical variation in aerosol size distribution. Comparison of simulated and in situ droplet size distributions (DSDs) shows that, to first order, integrated moments are always sensitive to sizes < ~ 30 µm and can also be sensitive to larger sizes if the DSDs are sufficiently broad, with implications for the assumed maximum observed size retrieved by the RSP. The bin scheme captures the observed extended tail of the DSD, while the bulk scheme is unable to due to parametric constraints. Differences in expected relationships between in situ measurements of cloud cores and cloud-top retrievals by RSP demonstrate difficulty in constraining well the case presented herein. Finally, a thermal-tracking framework demonstrates that the dilution of N d throughout a thermal's lifetime is heavily determined by collision–coalescence and the height-varying aerosol distribution and that, in the absence of these, the impact of entrainment on diluting N d is largely offset by secondary activation. Implications for evaluating warm-phase convective microphysics schemes in ESMs and translating results for use on global, space-based polarimetry platforms are discussed.

Stanford, McKenna Wallace [Columbia Univ., New Yor↗

DOE OSTI2026

Drop clustering and drop size correlations from holographic imagery suggest cloud droplet spectral broadening via entrainment-mixing

The question of how droplets rapidly grow large enough to initiate collision-coalescence has persisted for decades. Many theories explaining the production of sufficiently large drops (i.e., those in the “bottleneck” size range; ∼ 25–50 µm diameters) involve drop clustering on millimeter scales. A novel method is introduced to evaluate drop clustering trends particle-by-particle (i.e., the number/proximity of neighboring drops for given droplets; defined as drop clustering fields) which are diagnosed relative to drops within their shared drop environments – in contrast to previous studies which diagnose drop clustering of defined sample volumes, or in terms of absolute length scales. Specifically, this study evaluates the statistical likelihood that drops of a given size are associated with either a significant number of neighboring drops, or are significantly isolated from neighboring drops. Observations are acquired from the HOLODEC during the Cloud System Evolution in the Trades campaign, which sampled subtropical marine clouds. The HOLODEC measures drop size distributions and the 3D spatial coordinates of droplets. Results show drops within the bottleneck size range (diameters of ∼ 25–50 µm) are most likely to be significantly isolated from neighboring drops. This “isolated large drop trend” is primarily observed at subsaturated conditions, suggesting entrainment is the contributing factor. Holograms associated with this trend are more likely to have broader drop size distributions, larger maximum drop sizes and overly regions where precipitation reaches the lowest altitudes from the sampled cloud, suggesting entrainment-mixing drop size distribution broadening is a relevant precipitation-initiation mechanism.

D'Alessandro, John J. [Univ. of Washington, Seattl↗

DOE OSTI2020

Laser Disdrometer Quantities (LDQUANTS) and Video Disdrometer Quantities (VDISQUANTS) Value-Added Products Report

Disdrometers are useful instruments that measure the distribution of raindrop sizes (drop size distribution: DSD) and the associated rainfall rates/accumulation as those raindrops fall to the ground. These DSDs are a quantity of interest to members of the modeler and observational communities. However, use of disdrometer data for model evaluation, radar monitoring, or other activities requires careful quality control and processing for key DSD properties of interest (e.g., the number concentration of drops) to ensure appropriate physical (scattering, fall speed) assumptions. The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) Laser Disdrometer Quantities and Video Disdrometer Quantities Value-Added Product (LDQUANTS/VDISQUANTS VAP) uses standard methods from Tokay et al. (e.g., 2013, 2014) to filter drops with unrealistic fall speeds. Further, it estimates several microphysical/geophysical quantities of parameterized DSDs (gamma or exponential assumption type fitting methods) as in previous disdrometer studies and ARM long-term efforts (e.g., Testud et al. 2001, Giangrande et al. 2014, Thompson et al. 2015). Disdrometers are also beneficial for cross-checks with other instrumentation, including rain gauges and radars. To support research interests and related radar monitoring activities, this product calculates radar equivalent quantities, including dual polarization radar quantities (e.g., Z, Differential Reflectivity ZDR, etc.), using T-Matrix scattering and additional wavelength, temperature, and drop shape assumptions (e.g., Thurai et al. 2007). All of these efforts allow disdrometer data sets to become more useful and easily handled by modeling and observational studies, or routine instrumentation checks.

47 OTHER INSTRUMENTATION↗

DOE OSTI2021

Development of a Turbulent Liquid Spray Atomization Model for Diesel Engine Simulations (Final Technical Report)

This project addresses the systematic lack of predictive capabilities by spray models within engine CFD codes. We develop a new modeling approach to predict the breakup of diesel sprays based on recent literature showing that liquid turbulence plays a fundamental role in spray atomization. A new body of quantitative validation data is also developed as a critical element of the project, leveraging the joint capabilities of Georgia Tech’s high-pressure continuous-flow spray chamber and Argonne National Lab’s near-nozzle x-ray diagnostics at the Advanced Photon Source. This project contributes spatially-resolved measurements of drop size distribution within well-characterized diesel injectors, Spray A and D, from the Engine Combustion Network (ECN) to the engine combustion community for the first time. Utilizing this new body of measurements, we validate and demonstrate a new spray model for diesel sprays, termed the KH-Faeth model, that predicts global and local spray characteristic more accurately than the widely adopted and employed KH model. Predicted drop size distributions are seen to predict measured drops sizes both quantitatively and predictively, with accurate response in droplet size distributions over a wide range of ambient density, injection pressure, and injector nozzle size (Spray A and D) without model tuning. The KH-Faeth model can reduce error in the predicted centerline droplet size profile by up to 80% for ECN Spray D simulations when compared to use of the widely employed KH model.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2022

Limitations of Separate Cloud and Rain Categories in Parameterizing Collision-Coalescence for Bulk Microphysics Schemes

Warm rain collision-coalescence has been persistently difficult to parameterize in bulk microphysics schemes. We use a flexible bulk microphysics scheme with bin scheme process parameterizations, called AMP, to investigate reasons for the difficulty. AMP is configured in a variety of ways to mimic bulk schemes and is compared to simulations with the bin scheme upon which AMP is built. We find that an important limitation in traditional bulk schemes is the use of separate cloud and rain categories. When the drop size distribution is instead represented by a continuous distribution, the simulation of cloud-to-rain conversion is substantially improved. We also find large sensitivity to the threshold size to distinguish cloud and rain in traditional schemes; substantial improvement is found by decreasing the threshold from 40 to 25 μm. Neither the use of an assumed functional form for the size distribution nor the choice of predicted distribution moments has a large impact on the ability of AMP to simulate rain production. When predicting four total moments of the liquid drop size distribution, either with a traditional two-category, two-moment scheme with a reduced size threshold, or a four-moment single-category scheme, errors in the evolution of mass and the cloud size distribution are similar, but the single-category scheme has a substantially better representation of the rain size distribution. Optimal moment combinations for the single-category approach are investigated and appear to be linked more to the information content they provide for constraining the size distributions than to their correlation with collision-coalescence rates.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

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↗

DOE OSTI2021

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↗

DOE OSTI2022

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↗

DOE OSTI2021

Raindrop Size Spectrum in Deep Convective Regions of the Americas

This study compared drop size distribution (DSD) measurements on the surfaces, the corresponding properties, and the precipitation modes among three deep convective regions within the Americas. The measurement compilation corresponded to two sites in the midlatitudes: the U.S. Southern Great Plains and Córdoba Province in subtropical South America, as well as to one site in the tropics: Manacapuru in central Amazonia; these are all areas where intense rain-producing systems contribute to the majority of rainfall in the Americas’ largest river basins. This compilation included two types of disdrometers (Parsivel and 2D-Video Disdrometer) that were used at the midlatitude sites and one type of disdrometer (Parsivel) that was deployed at the tropical site. The distributions of physical parameters (such as rain rate R, mass-weighted mean diameter D m , and normalized droplet concentration N w ) for the raindrop spectra without rainfall mode classification seemed similar, except for the much broader N w distributions in Córdoba. The raindrop spectra were then classified into a light precipitation mode and a precipitation mode by using a cutoff at 0.5 mm h -1 based on previous studies that characterized the full drop size spectra. These segregated rain modes are potentially unique relative to previously studied terrain-influenced sites. In the light precipitation and precipitation modes, the dominant higher frequency observed in a broad distribution of N w in both types of disdrometers and the identification of shallow light precipitation in vertically pointing cloud radar data represent unique characteristics of the Córdoba site relative to the others. As a result, the co-variability between the physical parameters of the DSD indicates that the precipitation observed in Córdoba may confound existing methods of determining the rain type by using the drop size distribution.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Study on Emulsification/Demulsification Behavior and Mechanism of Produced Liquid from Polymer Flooding on Alaska North Slope

Heavy oil reservoirs on Alaska North Slope (ANS) are unconsolidated and contain abundant clay minerals, where the first-ever field pilot is currently implemented to validate the use of polymer floods for heavy oil enhanced oil recovery (EOR). The polymer molecules and/or fine clay particles carried with the produced liquid could potentially affect the oil/water separation, which is one of the major concerns for field operators. This dissertation aims to investigate the emulsification behavior of produced liquid, understand the emulsifying mechanism, and seek an adaptive and cost-effective method to treat the produced liquid from polymer flooding. Emulsions were prepared by mechanically mixing the actual heavy oil and the produced water from the pilot site, of which the stability was investigated by bottle test method or multiple light scattering method. Drop size distribution and interfacial properties were measured via microscope and pendant drop technique to probe the stability mechanism further. Results showed that oil-continuous or water-continuous emulsion could be generated depending on the water cut, clay types, clay concentration, and polymer concentration. In the crude oil/water system, the increasing water cut triggered the phase inversion of oil-continuous emulsion to water-continuous emulsion, resulting in faster separation and lower emulsion stability. Whereas, clay particles, no matter added to the oil or water phase, resulted in an unfavorable phase inversion from the loose water-continuous emulsion to the tight oil-continuous emulsion as clay concentration increased. For all four types of clay except Ca-montmorillonite, clay particles added to water led to an earlier phase inversion and higher emulsion stability than that added to the oil. The dual function of polymer on emulsion stability was observed. On the one hand, both sheared and unsheared polymer tended to convert the oil-continuous emulsion formed in either crude oil/water system or complex crude oil/water/clay particle system to the water-continuous emulsion, acting as a weak demulsifier to accelerate the oil/water separation. On the other hand, the addition of polymer to the water-continuous emulsion could result in enhanced emulsion stability, which is primarily attributed to the increased viscosity of the continuous phase and the decreased drop size of the dispersed oil phase. Particularly, the sheared polymer had a weaker ability to stabilize the o/w emulsion than the unsheared polymer due to the lower viscosity of the sheared polymer solution resulting from the breakdown of the polymer macromolecules. As for the chemical demulsification tests, the performance of demulsifiers showed a complex dependency upon the water cut, the shearing intensity, demulsifier type and dosage, and the polymer concentration. A compound emulsion breaker, E12+E18, exhibited the most satisfactory demulsification performance despite the varied test conditions. For severe water-continuous emulsions that might require a multi-fold dosage of demulsifier, a less expensive electrolyte, KCl, was proposed to be used in combination with demulsifier E12+E18 to improve the demulsification performance. In the proposed demulsifier formula, the effectiveness of the commercial demulsifier relied on its destructive effect on the interfacial film, while the efficacy of KCl was mainly dependent on its viscosity reduction effect on the continuous phase. This dissertation illustrates that intermediate layer elimination and water clarification are the major challenges for produced liquid treatment from polymer flooding. It also provides practical and theoretical guidance in advance for the demulsification strategy of the produced liquid from the ongoing first-ever polymer flooding pilot on ANS.

Chang, Hongli↗

DOE OSTI2023

Multi–year analysis of rain–snow levels at Marquette, Michigan

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.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

A Statistical Framework for Evaluating Rain Microphysics in Model Simulations and Disdrometer Observations

Abstract Statistical analyses of a large disdrometer data set and a diverse set of model simulations for convection using the Regional Atmospheric Modeling System were conducted, with the mutual goal of providing insights into precipitation formation and microphysical processes. We demonstrate that a two‐moment bulk microphysical model successfully captures the dominant observed modes of variability in rainfall related to rainfall intensity and raindrop size distributions. The model reproduced the general distribution of observed precipitation groups (PGs) derived from Principal Component Analysis. The multi‐variable analysis also uncovered some shortcomings in the model as well as limitations of the disdrometer data. The model solutions were constrained in their predicted drop size distributions (DSDs) due to the fixed DSD parameters assumed in a two‐moment microphysics scheme. A case study from the Mid‐latitude Continental Clouds and Convection Experiment field project demonstrated how model results can be used to contextualize the disdrometer observations which are limited in sample size, spatial coherence, and detection of small drops and low drop concentrations. The case study also showed that the spatial patterns of the statistically derived PGs revealed by the model are consistent with the hypothesized microphysical processes that determine surface rain DSDs. This work demonstrates how leveraging the strengths of observations and models together can improve our understanding and representation of rain microphysical processes.

54 ENVIRONMENTAL SCIENCES↗