SEARCH · Engineering Papers

Results for “Droplet 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 469 records · Page 26

All sources NASA NTRS DOE OSTI

DOE OSTI2023

CPC data from TAMU TRACER campaign in the Houston, TX region from July to September 2022

During TRACER, the Texas A&M Rapid Onsite Atmospheric Measurements Van (ROAM-V) was deployed to capture airmasses behind (maritime) and ahead (continental) of the passage of the sea-breeze front through Houston. On select sampling days, ROAM-V sampled in the morning/mid-day on the coast and then transited to a second inland site for the afternoon/evening. The suite of instruments deployed on ROAM-V included a Condensation Particle Counter (CPC; GRIMM Model 5.403 CPC), Scanning Mobility Particle Sizer (SMPS; TSI 3750 detector, TSI 3082 classifier, TSI 3088 neutralizer, TSI 3081A Differential Mobility Analyzer), Cloud Condensation Nuclei counter (Droplet Measurement Technologies CCN Counter), micro pulse lidar (Droplet Measurement Technologies Micro Pulse LiDAR (miniMPL)), and a Davis Rotating Uniform size-cut Monitor (DRUM; DRUMAir 4-DRUM). Before sampling at each location, the latitude and longitude were recorded using the GPS on the phone application “My Altitude”. Onboard the ROAM-V, aerosol samples are drawn through a shared isokinetic inlet at a flow rate ranging from 3.5 to 7.0 LPM. A portion of this flow is directed through a cyclone impactor (Brechtel, Inc. Model SCC 0.732) and 0.3 LPM is directed to the CPC. To correct the data for particle losses, we used a two-step method. First, the measured SMPS size distributions were used to calculate total particle loss through the inlet for every SMPS scan at a single deployment location. Second, the CPC data was corrected for particle losses using the average of the total losses per scan. This calculation was done separately for each deployment location due to changes in the measured size distributions between locations. Particle loss from diffusion (based on Kesten, 1991 and Gormley, 1949), inertial impaction in 90-degree bends (based on Aerosol Measurement, 2011 and Crane, 1977), and cyclone impactor efficiency (based on Dirgo, 1985) were included in the loss calculation. When the SMPS was not sampling at a location (in the case of an instrument malfunction or operator error), the reported CPC data was corrected with an average of the total losses for the entire campaign at the specified deployment location (e.g., if we needed to correct Galveston data, then the average of all calculated losses at Galveston was taken). These flatline corrections were used for all data on 22/07/13, 22/07/20, 22/07/22, and the data from Galveston on 22/08/09. An additional flatline correction of 14% was applied to all CPC data based on laboratory calibration prior to and after the campaign. In laboratory calibrations, we identified that for the same sample of air the field CPC (GRIMM Model 5.403 CPC) undercounted the total concentration of particles by 14% compared to the ground truth laboratory CPC (TSI Model 3750). This data was collected for ARM Field Campaign AFC07055 and supported by DOE ASR grant DE-SC0021047. For any further questions, please feel free to contact the instrument PI, Sarah D. Brooks, sbrooks@tamu.edu. Kesten et. al. Calibration of a TSI Model 3025 Ultrafine Condensation Particle Counter. Aerosol Science and Technology, 15:2, 107-111, 1991. Gormley et. al. Diffusion from a Stream Flowing through a Cylindrical Tube. Proceedings of the Royal Irish Academy, Vol 52, 163-169, 1948. Aerosol Measurement: Principles, Techniques, and Applications, Third Edition. John Wiley & Sons, Inc, 2011. Crane et. al. Inertial Deposition of Particles in a Bent Pipe. Journal of Aerosol Science, Vol 8, 161-170, 1977. Dirgo et. al. Cyclone Collection Efficiency: Comparison of Experimental Results with Theoretical Predictions. Aerosol Science and Technology, 4:4, 401-415, 1985.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1993

LRAT: Lightning Radiative Transfer

In this report, we extend to cloud physics the work done for single and multiple scattering of electromagnetic waves. We consider the scattering of light, visible or infrared, by a spherical cloud represented by a statistically homogeneous ensemble of configurations of N identical spherical water droplets whose centers are uniformly distributed in its volume V. The ensemble is specified by the average number rho of scatterers in unit volume and by rho f(R) with f(R) as the distribution function for separations R of pairs. The incident light, vector-phi(sub 0) a plane electromagnetic wave with harmonic time dependence, is from outside the cloud. The propagation parameter kappa(sub 0) and the index of refraction eta(sub 0) determine physically the medium outside the distribution of scatterers. We solve the interior problem separately to obtain the bulk parameters for the scatterer equivalent to the ensemble of spherical droplets. With the interior solution or the equivalent medium approach, the multiple scattering problem is reduced to that of an equivalent single scatterer excited from outside illumination. A dispersion relation which determines the bulk propagation parameter K and the bulk index of refraction eta of the cloud is given in terms of the vector equivalent scattering amplitude vector-G and the dyadic scattering amplitude tilde-g of the single object in isolation. Based on this transfer model we will have the ability to consider clouds composed of inhomogeneous distribution of water and/or ice particles and we will be able to take into account particle size distributions within the cloud. We will also be able to study the effects of cloud composition (i.e., particle shape, size, composition, orientation, location) on the polarization of the single or the multiple scattered waves. Finally, this study will provide a new starting point for studying the problem of lightning radiative transfer.

Phanord, Dieudonne D.↗

NASA NTRS2019

Observation of Absorbing Aerosols Above Clouds over the South-East Atlantic Ocean from the Geostationary Satellite SEVIRI – Part 1: Method Description and Sensitivity

High temporal resolution observations from satellites have a great potential for studying the impact of biomass burning aerosols and clouds over the South East Atlantic Ocean (SEAO). This paper presents a method developed to retrieve simultaneously aerosol and cloud properties in aerosol above cloud conditions from the geostationary instrument Meteosat Second Generation/Spinning Enhanced Visible and Infrared Imager (MSG/SEVIRI). The above-cloud Aerosol Optical Thickness (AOT), the Cloud Optical Thickness (COT) and the Cloud droplet Effective Radius (CER) are derived from the spectral contrast and the magnitude of the signal measured in three channels in the visible to shortwave infrared region. The impact of the absorption from atmospheric gases on the satellite signal is corrected by applying transmittances calculated using the water vapour profiles from a Met Office forecast model. The sensitivity analysis shows that a 10% error on the humidity profile leads to an 18.5% bias on the above-cloud AOT, which highlights the importance of an accurate atmospheric correction scheme. In situ measurements from the CLARIFY-2017 airborne field campaign are used to constrain the aerosol size distribution and refractive index that is assumed for the aforementioned retrieval algorithm. The sensitivities in the retrieved AOT, COT and CER to the aerosol model assumptions are assessed. Between 09:00-15:00 UTC, an uncertainty of 40% is estimated on the above-cloud AOT, which is dominated by the sensitivity of the retrieval to the single scattering albedo. The absorption AOT is less sensitive to the aerosol assumptions with an uncertainty generally lower than 17% between 09:00-15:00 UTC. Outside of that time range, as the scattering angle decreases, the sensitivity of the AOT and the absorption AOT to the aerosol model increases. The retrieved cloud properties are only weakly sensitive to the aerosol model assumptions throughout the day, with biases lower than 6% on the COT and 3% on the CER. The stability of the retrieval over time is analysed. For observations outside of the backscattering glory region, the time-series of the aerosol and cloud properties are physically consistent, which confirms the ability of the retrieval to monitor the temporal evolution of aerosol above cloud events over the SEAO.

Peers, Fanny↗

NASA NTRS2019

Observation of absorbing aerosols above clouds over the south-east Atlantic Ocean from the geostationary satellite SEVIRI – Part 1: Method description and sensitivity

High-temporal-resolution observations from satellites have a great potential for studying the impact of biomass burning aerosols and clouds over the south-east Atlantic Ocean (SEAO). This paper presents a method developed to simultaneously retrieve aerosol and cloud properties in aerosol above-cloud conditions from the geostationary instrument Meteosat Second Generation/Spinning Enhanced Visible and Infrared Imager (MSG/SEVIRI). The above-cloud aerosol optical thickness (AOT), the cloud optical thickness (COT) and the cloud droplet effective radius (CER) are derived from the spectral contrast and the magnitude of the signal measured in three channels in the visible to shortwave infrared region. The impact of the absorption from atmospheric gases on the satellite signal is corrected by applying transmittances calculated using the water vapour profiles from a Met Office forecast model. The sensitivity analysis shows that a 10 % error on the humidity profile leads to an 18.5 % bias on the above-cloud AOT, which highlights the importance of an accurate atmospheric correction scheme. In situ measurements from the CLARIFY-2017 airborne field campaign are used to constrain the aerosol size distribution and refractive index that is assumed for the aforementioned retrieval algorithm. The sensitivities in the retrieved AOT, COT and CER to the aerosol model assumptions are assessed. Between 09:00 and 15:00 UTC, an uncertainty of 40 % is estimated on the above-cloud AOT, which is dominated by the sensitivity of the retrieval to the single-scattering albedo. The absorption AOT is less sensitive to the aerosol assumptions with an uncertainty generally lower than 17 % between 09:00 and 15:00 UTC. Outside of that time range, as the scattering angle decreases, the sensitivity of the AOT and the absorption AOT to the aerosol model increases. The retrieved cloud properties are only weakly sensitive to the aerosol model assumptions throughout the day, with biases lower than 6 % on the COT and 3 % on the CER. The stability of the retrieval over time is analysed. For observations outside of the backscattering glory region, the time series of the aerosol and cloud properties are physically consistent, which confirms the ability of the retrieval to monitor the temporal evolution of aerosol above-cloud events over the SEAO.

cloud properties↗

DOE OSTI2023

SMPS data from TAMU TRACER campaign in the Houston TX region from July to September 2022

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Aerosol Microphysics and Chemical Measurements at Mt. Soledad and Scripps Pier during the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE) from February 2023 to February 2024 UCSD Library Collection

This dataset includes guest instrument measurements and other PI products for aerosol microphysics and chemical measurements collected at Mt. Soledad and Scripps Pier during the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE) from February 2023 to February 2024. The measurements include the following instruments at Mt. Soledad: High-Resolution Time-of-Flight Aerosol Mass Spectrometer (HR-ToF-AMS, Aerodyne), Scanning Electrical Mobility Spectrometer (SEMS, Brechtel Manufacturing Inc.), Aerodynamic Particle Sizer (APS, Droplet Measurements Technologies), Single Particle Soot Photometer (SP2, Drople Measurements Technologies), Meteorological Station (WXT520, Vaisala), Ozone (Teco), and trace gas proxies (Teledyne). In addition, the analyses of particle filters collected at Mt. Soledad for three dry-diameter size cuts (<1 micron, <0.5 micron, <0.18 micron) and at Scripps Pier for one dry-diameter size cut (<1 micron) by Fourier Transform Infrared (FTIR) and X-ray Fluorescence (XRF) are reported. A differential mobility analyzer operated as a scanning mobility particle sizer (SMPS, TSI Inc.), a printed particle optical spectrometer (POPS, Grimm), and a continuous flow diffusion cloud condensation nuclei (CCN, DMT) counter provide the mobility aerosol size distribution (30-360 nm), optical size distribution (150 - 6000 nm), size-resolved CCN distribution (30-360 nm) at 0.2, 0.4, 0.6, 0.8, and 1.0% supersaturation. Measurements are reported for both sampling from an isokinetic aerosol inlet and from a Counterflow Virtual Impactor (CVI, Brechtel Manufacturing Inc.). The data are available at the following link: https://library.ucsd.edu/dc/collection/bb0898306q

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1999

Comparison of Techniques for Non-Intrusive Fuel Drop Size Measurements in a Subscale Gas Turbine Combustor

In aviation gas turbine combustors, many factors, such as the degree and extent of fuel/air mixing and fuel vaporization achieved prior to combustion, influence the formation of undesirable pollutants. To assist in analyzing the extent of fuel/air mixing, flow visualization techniques have been used to interrogate the fuel distributions during subcomponent tests of lean-burning fuel injectors. Discrimination between liquid and vapor phases of the fuel was determined by comparing planar laser-induced fluorescence (PLIF) images, elastically-scattered light images, and phase/Doppler interferometer measurements. Estimates of Sauter mean diameters are made by ratioing PLIF and Mie scattered intensities for various sprays, and factors affecting the accuracy of these estimates are discussed. Mie calculations of absorption coefficients indicate that the fluorescence intensities of individual droplets are proportional to their surface areas, instead of their volumes, due to the high absorbance of the liquid fuel for the selected excitation wavelengths.

Zaller, Michelle↗

NASA NTRS1993

Optical properties of marine stratocumulus clouds modified by ships

Results are presented of an application of the diffusion domain method to multispectral solar radiation measurements obtained deep within a marine stratocumulus cloud layer modified by pollution from ships. In situ airborne measurements of the relative angular distribution of scattered radiation are compared to known asymptotic expressions for the intensity field deep within an optically thick cloud layer. Analytical expressions relating the ratio of the nadir-to-zenith intensities to surface reflectance, similarity parameter, and scaled optical depth beneath the aircraft flight level are used to analyze measurements obtained with the cloud absorption radiometer mounted on the University of Washington's C-131A research aircraft. It is shown that the total optical thickness of the cloud layer increased in the ship tracks, in contrast to the similarity parameter, which decreased. The decrease in absorption was a direct consequence of the reduction in cloud droplet size that occurred within the ship tracks.

King, Michael D.↗

NASA NTRS1996

Independent Pixel and Two Dimensional Estimates of LANDSAT-Derived Cloud Field Albedo

A theoretical study has been conducted on the effects of cloud horizontal inhomogeneity on cloud albedo bias. A two-dimensional (2D) version of the Spherical Harmonic Discrete Ordinate Method (SHDOM) is used to estimate the albedo bias of the plane parallel (PP-IPA) and independent pixel (IPA-2D) approximations for a wide range of 2D cloud fields obtained from LANDSAT. They include single layer trade cumulus, open and closed cell broken stratocumulus, and solid stratocumulus boundary layer cloud fields over ocean. Findings are presented on a variety of averaging scales and are summarized as a function of cloud fraction, mean cloud optical depth, cloud aspect ratio, standard deviation of optical depth, and the gamma function parameter Y (a measure of the width of the optical depth distribution). Biases are found to be small for small cloud fraction or mean optical depth, where the cloud fields under study behave linearly. They are large (up to 0.20 for PP-IPA bias, -0.12 for IPA-2D bias) for large v. On a scene average basis PP-IPA bias can reach 0.30, while IPA-2D bias reaches its largest magnitude at -0.07. Biases due to horizontal transport (IPA-2D) are much smaller than PP-IPA biases but account for 20% RMS of the bias overall. Limitations of this work include the particular cloud field set used, assumptions of conservative scattering, constant cloud droplet size, no gas absorption or surface reflectance, and restriction to 2D radiative transport. The LANDSAT data used may also be affected by radiative smoothing.

Chambers, L. H.↗

DOE OSTI2022

Sensitivity of puff dynamics and airborne droplet nuclei distribution to variations in violent expiration events

Here, we use large eddy simulations to investigate the puff and droplet dynamics from violent expiratory events such as coughs and sneezes in the first few seconds following an ejection. For each of the eleven simulations considered, over 60,000 droplets are ejected and individually tracked using the point-particle Euler–Lagrange approach. We test the sensitivity of the puff and droplet dynamics to various parameters including the ejection volume, momentum, and orientation. We also explore the effect of the mouth shape on the aforementioned dynamics by considering elliptical and circular inlet cross-sections. The results from the simulations compare favorably with a recent theoretical framework put forth by Balachandar et al. in terms of the puff size and propagation velocity. More importantly however, the theory is able to accurately predict the number and size spectra of the potentially virus-laden droplet nuclei that remain airborne within the puff. We observe that the ejection angle and mouth shape do not significantly affect the puff and droplet dynamics. Additionally, we quantify the carrying capacity of the detached puff portions in terms of the number and size spectra of droplets/droplet nuclei suspended within.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1999

Dynamics of Single Chains of Suspended Ferrofluid Particles

We present an experimental study of the dynamics of isolated chains made of super-paramagnetic particles under the influence of a magnetic field. The motivation of this work is to understand if the chain fluctuations exist and, if it does, how does the fluctuation affect chain aggregation. We find that single chains strongly fluctuate and that the characteristic frequency of their fluctuations is inversely proportional to the magnetic field strength. The higher the field the lower the characteristic frequency of the chain fluctuations. In the high magnetic field limit, chains behave like rigid rods without any internal motions. In this work, we used ferrofluid particles suspended in water. These particles do not have any intrinsic magnetization. Once a magnetic field is applied, a dipole moment is induced in each particle, proportional to the magnetic field. A dipolar magnetic interaction then occurs between particles. If dipole-dipole magnetic energy is higher than the thermal energy, the result is a structure change inside the dipolar fluid. The ratio of these two energies is expressed by a coupling constant lambda as: lambda = (pi(a(exp 3))(chi(exp 2))(mu(sub 0))(H(sub 0))(exp 2))/18kT Where a is the particle radius, mu(sub 0) is the vacuum magnetic permeability, H(sub 0) the applied magnetic field, k the Boltzmann constant and T the absolute temperature. If lambda > 1, magnetic particles form chains along the field direction. The lateral coalescence of several chains may form bigger aggregates especially if the particle volume fraction is high. While many studies and applications deal with the rheological properties and the structural changes of these dipolar fluids, this work focuses on the understanding of the chain dynamics. In order to probe the chain dynamics, we used dynamic light scattering (DLS) in self-beating mode as our experimental technique. The experimental geometry is such that the scattering plane is perpendicular to the magnetic field. Therefore, only motions in this plane are probed. A very dilute sample of a ferrofluid emulsion with a particle volume fraction of 10(exp -5) is used in this experiment. We chose such a low volume fraction to avoid multiple light scattering as well as lateral chain-chain aggregation. DLS measures the dynamic structure factor S(q,t) of the sample (q is the scattering wave vector, t is the time). In the absence of the magnetic field, identical particles of ferrofluid droplets are randomly distributed and S(q,t) reduces to exp(-q(exp 2)2D(sub 0)t). D(sub 0)=(kT/(6(pi)(eta)(a)) is the diffusion coefficient of Brownian particles (where Xi = (6(pi)(eta)(a)) is the Stokes frictional coefficient of a spherical particle in a fluid of viscosity eta). If interactions or polydispersity can not be ignored, an effective diffusion coefficient is introduced. Formally, D(sub eff) is defined as: D(sub eff) = - q(exp -2) partial derivative of (ln(S(q,t)) with respect to time, as t goes to 0. D(sub eff) reduces to D(sub 0) if no interactions and only a few particles size are present. Therefore, we can use DLS to measure particle size. The particle radius was found to be a=0.23 mu m with 7% of polydispersity. In this case, if we vary the scattering angle theta (and so q) we do not have any change in the measured diffusion coefficient: it is q-independent. When a magnetic field is applied, particles aggregate into chains if lambda > 1. We first studied the kinetics of the chain formation when lambda = 406. At a fixed scattering angle, we measured diffusion coefficient D(sub eff) as a function of time. Experimentally, we find that D(sub eff) decreases monotonously with time. Physically, this means that chains are becoming longer and longer. Since we are only sensitive to motions in the scattering plane and since chains have their main axis perpendicular to this plane, the measured diffusion coefficient is the trans-verse diffusion coefficient. We can relate D(sub eff) to the mean number of particles per chain N(t) at a given time and to the diffusion coefficient of an isolated particle D(sub 0) as D(sub eff)=f(N(t))D(sub 0). Since f(N) is known from other recent work, N can be expressed as a function of the time. We found a square root dependency: N(t) proportional to the square root of t. As expected for very low volume fraction, this behavior is characteristic of a diffusion-limited aggregation as suggested by several authors and by our previous work. In this study, we focus on the dependence of the effective diffusion coefficient on the scattering angle and the magnetic field strength. After the magnetic field is applied (lambda = 406) for a long time, typically 6 hours, kinetics of chain formation becomes very slow. Chain size does not vary much over the next hour period. Thus, we can perform different interesting experiments. First, at a fixed magnetic field, we measure the effective diffusion coefficient as a function of the scattering angle (from 5 to 130 deg). Our results show that the measured diffusion coefficient increases linearly with the scattering angle: D(sub eff) proportional to q. If we do the same experiment for different lambda values, D(sub eff) depends on lambda as D(sub eff) proportional to lambda(exp -1/2). We also find for different lambda values that the same asymptotic D(sub eff) value is obtained when q approaches zero. The angle dependency of D(sub eff) suggests that an additional motion exists besides chain drifting. Chain size is constant during experiment, which was verified by measuring the same diffusion coefficient at the beginning and at the end of the angle switching. If chains are rigid, D(sub eff) is independent of q. Therefore, we found that D(sub eff) not only measures the motion of the entire chain but also its internal fluctuations. These internal motions are the fluctuations of the particles in the chain. To understand the q dependency of D(sub eff), let us look at the probing length used. In our study, the characteristic length scale probed is l=2pi/q which is in the range of 0.9<l/a<20. When l is much larger than the particle radius (l/a going to infinity) we are mainly sensitive to the center of mass diffusion of the chain. D(sub eff) is then the diffusion of the entire chain and depends only on N but not on q and lambda. The value D(sub eff) allows us to obtain the number of particles per chain N. In the opposite limit (i.e. l/a < 1), we are sensitive to motions on the size of individual particles. The main contributions to the measured diffusion coefficient come from internal motions of the chains (i.e. particles' fluctuations inside the chain). We investigated also the effect of the chain size on D(sub eff). For the same value of the magnetic interaction lambda and different chain sizes, we found that the slope of D(sub eff) versus q decreases when N increases as a power law D(sub eff) proportional to N(exp -0.7). This is an important behavior because this means that the longer the chain the slower its internal fluctuations besides the whole chain motion.

Cutillas, S.↗

DOE OSTI2024

Glaciation of mixed-phase clouds: insights from bulk model and bin-microphysics large-eddy simulation informed by laboratory experiment

Abstract. Mixed-phase clouds affect precipitation and radiation differently from liquid and ice clouds, posing greater challenges to their representation in numerical simulations. Recent laboratory experiments using the Pi Cloud Chamber explored cloud glaciation conditions based on increased injection of ice-nucleating particles. In this study, we use two approaches to reproduce the results of the laboratory experiments: a bulk scalar mixing model and large-eddy simulation (LES) with bin microphysics. The first approach assumes a well-mixed domain to provide an efficient assessment of the mean cloud properties for a wide range of conditions. The second approach resolves the energy-carrying turbulence, the particle size distribution, and their spatial distribution to provide more details. These modeling approaches enable a separate and detailed examination of liquid and ice properties, which is challenging in the laboratory. Both approaches demonstrate that, with an increased ice number concentration, the flow and microphysical properties exhibit the same changes in trends. Additionally, both approaches show that the ice integral radius reaches the theoretical glaciation threshold when the cloud is subsaturated with respect to liquid water. The main difference between the results of the two approaches is that the bulk model allows for the complete glaciation of the cloud. However, LES reveals that, in a dynamic system, the cloud is not completely glaciated as liquid water droplets are continuously produced near the warm lower boundary and subsequently mixed into the chamber interior. These results highlight the importance of the ice mass fraction in distinguishing the mixed-phase clouds and ice clouds.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1992

Numerical modeling for dilute and dense sprays

We have successfully implemented a numerical model for spray-combustion calculations. In this model, the governing gas-phase equations in Eulerian coordinate are solved by a time-marching multiple pressure correction procedure based on the operator-splitting technique. The droplet-phase equations in Lagrangian coordinate are solved by a stochastic discrete particle technique. In order to simplify the calculation procedure for the circulating droplets, the effective conductivity model is utilized. The k-epsilon models are utilized to characterize the time and length scales of the gas phase in conjunction with turbulent modulation by droplets and droplet dispersion by turbulence. This method entails random sampling of instantaneous gas flow properties and the stochastic process requires a large number of computational parcels to produce the satisfactory dispersion distributions even for rather dilute sprays. Two major improvements in spray combustion modelings were made. Firstly, we have developed a probability density function approach in multidimensional space to represent a specific computational particle. Secondly, we incorporate the Taylor Analogy Breakup (TAB) model for handling the dense spray effects. This breakup model is based on the reasonable assumption that atomization and drop breakup are indistinguishable processes within a dense spray near the nozzle exit. Accordingly, atomization is prescribed by injecting drops which have a characteristic size equal to the nozzle exit diameter. Example problems include the nearly homogeneous and inhomogeneous turbulent particle dispersion, and the non-evaporating, evaporating, and burning dense sprays. Comparison with experimental data will be discussed in detail.

Chen, C. P.↗

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

Direct numerical simulation of turbulent dispersion of evaporative aerosol clouds produced by an intense expiratory event

Airborne particles are a major route for transmission of COVID-19 and many other infectious diseases. When a person talks, sings, coughs, or sneezes, nasal and throat secretions are spewed into the air. After a short initial fragmentation stage, the expelled material is mostly composed of spherical particles of different sizes. While the dynamics of the largest droplets are dominated by gravitational effects, the smaller aerosol particles, mostly transported by means of hydrodynamic drag, form clouds that can remain afloat for long times. In subsaturated air environments, the dependence of pathogen-laden particle dispersion on their size is complicated due to evaporation of the aqueous fraction. Particle dynamics can significantly change when ambient conditions favor rapid evaporation rates that result in a transition from buoyancyto-drag dominated dispersion regimes. To investigate the effect of particle size and evaporation on pathogen-laden cloud evolution, a direct numerical simulation of a mild cough was coupled with an evaporative Lagrangian particle advection model. The results suggest that while the dispersion of cough particles in the tails of the size distribution are unlikely to be disrupted by evaporative effects, preferential aerosol diameters (30-40 μm) may exhibit significant increases in the residence time and horizontal range under typical ambient conditions. Using estimations of the viral concentration in the spewed fluid and the number of ejected particles in a typical respiratory event, we obtained a map of viral load per volume of air at the end of the cough and the number of virus copies per inhalation in the emitter vicinity.

60 APPLIED LIFE SCIENCES↗

NASA NTRS1979

Development of a Model and Computer Code to Describe Solar Grade Silicon Production Processes

Mathematical models and computer codes based on these models, which allow prediction of the product distribution in chemical reactors for converting gaseous silicon compounds to condensed-phase silicon were developed. The following tasks were accomplished: (1) formulation of a model for silicon vapor separation/collection from the developing turbulent flow stream within reactors of the Westinghouse (2) modification of an available general parabolic code to achieve solutions to the governing partial differential equations (boundary layer type) which describe migration of the vapor to the reactor walls, (3) a parametric study using the boundary layer code to optimize the performance characteristics of the Westinghouse reactor, (4) calculations relating to the collection efficiency of the new AeroChem reactor, and (5) final testing of the modified LAPP code for use as a method of predicting Si(1) droplet sizes in these reactors.

Srivastava, R.↗

NASA NTRS2017

Numerical Analysis of Mixed-Phase Icing Cloud Simulations in the NASA Propulsion Systems Laboratory

This presentation describes the development of a numerical model that couples the thermal interaction between ice particles, water droplets, and the flowing gas of an icing wind tunnel for simulation of NASA Glenn Research Centers Propulsion Systems Laboratory (PSL). The ultimate goal of the model is to better understand the complex interactions between the test parameters and have greater confidence in the conditions at the test section of the PSL tunnel. The model attempts to explain the observed changes in test conditions by coupling the conservation of mass and energy equations for both the cloud particles and flowing gas mass. Model predictions were compared to measurements taken during May 2015 testing at PSL, where test conditions varied gas temperature, pressure, velocity and humidity levels, as well as the cloud total water content, particle initial temperature, and particle size distribution.

NASA PSL↗

NASA NTRS2017

Numerical Analysis of Mixed-Phase Icing Cloud Simulations in the NASA Propulsion Systems Laboratory

This paper describes the development of a numerical model that couples the thermal interaction between ice particles, water droplets, and the flowing gas of an icing wind tunnel for simulation of NASA Glenn Research Centers Propulsion Systems Laboratory (PSL). The ultimate goal of the model is to better understand the complex interactions between the test parameters and have greater confidence in the conditions at the test section of the PSL tunnel. The model attempts to explain the observed changes in test conditions by coupling the conservation of mass and energy equations for both the cloud particles and flowing gas mass. Model predictions were compared to measurements taken during May 2015 testing at PSL, where test conditions varied gas temperature, pressure, velocity and humidity levels, as well as the cloud total water content, particle initial temperature, and particle size distribution.

NASA PSL↗