Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mixing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Comparing measurements of the mixing layer height to modeled values

The mixing layer height (MLH) is the top layer of turbulent mixing within the lower atmosphere, above the Earth’s surface. Estimates of the mixing layer height allow us to determine the volume available for the dispersion of pollutants throughout the atmosphere. Our goal is to identify the most suitable mixing layer height input for our dispersion modeling tool (AERMOD). For this project we evaluated two different methods of estimating the local mixing layer height. First, we use 3 estimates of the mixing layer height obtained from a ceilometer, each height estimate corresponds to a different gradient in aerosol backscatter which is used as a proxy for mixing layer height identification. Second, we use AERMET, our AERMOD modeling system preprocessor for meteorological data, which estimates the mixing layer height using several equations combined with measured meteorological data. We evaluate the ceilometer and AERMET estimates over a 5-year period (2015-2019) to see how well aligned the estimated mixing layer heights are. Our results suggest that our model, AERMET, is on average, aligned with ceilometer estimates during the daytime hours. However, daily maximum MLH estimates by AERMET occur earlier in the day than those estimated by the ceilometer. We believe that the ceilometer struggles with accurately measuring the MLH during nighttime hours as a result of sensor limitations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Similarity of the turbulent kinetic energy dissipation rate distribution in the upper mixed layer of the tropical Indian Ocean

Turbulence within the upper ocean mixed layer plays a key role in various physical, biological, and chemical processes. Between September and November 2011, a dataset of 570 vertical profiles of the turbulent kinetic energy (TKE) dissipation rate, as well as conventional hydrological and meteorological data, were collected in the upper layer of the tropical Indian Ocean. These data were used to statistically analyze the vertical distribution of the TKE dissipation rate in the mixed layer. The arithmetic-mean method made the statistical TKE dissipation rate profile more scattered than the median and geometric-mean methods. The statistical TKE dissipation rate were respectively scaled by the surface buoyancy flux and the TKE dissipation rate at the mixed-layer base. It was found that the TKE dissipation rate scaled by that at the mixed-layer base exhibited better similarity characteristics than that scaled by the surface buoyancy flux, whether the stability parameter D/|L MO | was greater or less than 10, indicating that the TKE dissipation rate at the mixed-layer base is a better characteristic scaling parameter for reflecting the intrinsic structure of the TKE dissipation rate in the mixed layer, where D and L MO are respectively the mixed-layer thickness and the Monin-Obukhov length scale. The parameterization of the TKE dissipation rate at the mixed-layer base on the shear-driven dissipation rate and the surface buoyancy flux was further explored. It was found that the TKE dissipation rate at the mixed-layer base could be well fitted by a linear combination of three terms: the wind-shear-driven dissipation rate, the surface buoyancy flux, and a simple nonlinear coupling term of these two.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Measurement report: Comparison of wintertime individual particles at ground level and above the mixed layer in urban Beijing

Beijing has been suffering from frequent severe air pollution events, with concentrations affected significantly by the mixed-layer height. Major efforts have been made to study the physico-chemical properties, compositions, and sources of aerosol particles at ground level. However, little is known about the morphology, elemental composition, and mixing state of aerosol particles above the mixed layer. In this work, we collected individual aerosol particles simultaneously at ground level (2 m above ground) and above the mixed layer in urban Beijing (within the Atmospheric Pollution and Human Health in a Chinese Megacity, APHH-Beijing, 2016 winter campaign). The particles were analyzed offline by transmission electron microscopy coupled with energy dispersive X-ray spectroscopy. Our results showed that the relative number contribution of mineral particles to all measured particles was much higher during non-haze periods (42.5 %) than haze periods (18.1 %); in contrast, internally mixed particles contributed more during haze periods (21.9 %) than non-haze periods (7.2 %) at ground level. In addition, more mineral particles were found at ground level than above the mixed-layer height. Around 20 % of individual particles showed core–shell structures during haze periods, whereas only a few core–shell particles were observed during non-haze periods (2 %). The results showed that the particles above the mixed layer were more aged, with a larger proportion of organic particles originating from coal combustion. Our results indicate that a large fraction of the airborne particles above the mixed layer come from surrounding areas influenced by coal combustion activities. This source contributes to the surface particle concentrations in Beijing when polluted air is mixed down to the ground level.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1993

CFD mixing analysis of axially opposed rows of jets injected into confined crossflow

A computational fluid dynamics (CFD) parametric study was performed to analyze axially opposed rows of jets mixing with crossflow in a rectangular duct. Isothermal analysis was conducted to determine the influence of lateral geometric arrangement on mixing. Two lateral arrangements were analyzed: (1) inline (jets' centerlines aligned with each other on top and bottom walls), and (2) staggered (jets' centerlines offset with each other on top and bottom walls). For a jet-to-mainstream mass flow ratio (MR) of 2.0, design parameters were systematically varied for jet-to-mainstream momentum-flux ratios (J) between 16 and 64 and orifice spacing-to-duct height ratios (S/H) between 0.125 and 1.5. Comparisons were made between geometries optimized for S/H at a specified J. Inline configurations had a unique spacing for best mixing at a specified J. In contrast, staggered configurations had two 'good mixing' spacings for each J, one corresponding to optimum inline spacing and the other corresponding to optimum non-impinging jet spacing. The inline configurations, due to their smaller orifice size at optimum S/H, produced better initial mixing characteristics. At downstream locations (e.g. x/H of 1.5), the optimum non-impinging staggered configuration produced better mixing than the optimum inline configuration for J of 64; the opposite results were observed for J of 16. Increasing J resulted in better mixing characteristics if each configuration was optimized with respect to orifice spacing. Mixing performance was shown to be similar to results from previous dilution jet mixing investigations (MR less than 0.5).

Bain, D. B.↗
NASA NTRS1994

Quantitative investigation of compressible mixing: Staged transverse injection into Mach 2 flow

Planar measurements of the injectant mole fraction distribution and the velocity field within a supersonic mixing flowfield have been made using laser-induced iodine fluorescence. The flowfield investigated in this work is staged transverse injection of air into a Mach 2 freestream. A complete three-dimensional survey of the injectant mole fraction distribution has been generated, and a single planar velocity measurement has been completed. The measurements reveal the dramatic effect of streamwise vortices on the mixing in the near field of the injectors, as well as the rapid mixing generated by staging two fuel injectors. Analysis of the downstream decay of the maximum injectant mole fraction in this and other supersonic mixing flowfields indicates that the relative rate of injectant mixing well downstream of the injectors is independent of injection geometry, freestream Mach number, and injectant molecular weight. Mixing within this region of the flowfield is dominated by small-scale turbulence within the injectant plume. The transition of the dominant mixing mechanism, from vortex-driven mixing in the near field to small-scale turbulent mixing in the far field, was found to occur in the region about 10 diameters downstream of the injectors.

Hollo, Steven D.↗
NASA NTRS2004

South Polar Ar Enhancement as a Tracer for Southern Winter Horizontal Meridional Mixing

Measurements made by the Gamma Ray Spectrometer (GRS) on Mars Odyssey during 2002 and 2003 show an obvious increase in the gamma flux of 1294 keV gamma rays resulting from the decay of (41)Ar. (41)Ar is made by the capture of thermal neutrons by atmospheric (40)Ar. The increase measured above the southern polar region has permitted calculation of the increase in mixing ratio of Ar from L(sub s) 8 to 100 between latitudes 75 S and 90 S. The peak in Ar enhancement occurs about 200 Earth days after CO2 freeze-out has begun, indicating that up to this time equatorward meridional mixing is rapid enough to move enhanced Ar from the polar regions northward. Although the CO2 frost depth continues to increase from L(sub s) 110 deg to 190 deg, the Ar enhancement steadily decreases to its baseline value reached at about L(sub s) 200 deg. Our data permit an estimate of the horizontal eddy mixing coefficient useful for constraining equatorward meridional mixing during southern winter and a characteristic mixing time for the polar southern winter atmosphere. Also, using the drop in excess Ar measured by the GRS from L(sub s) 110 deg to 200 deg, we estimate an eddy coefficient appropriate for meridional mixing of the entire Ar excess back to the baseline value. The horizontal eddy mixing coefficients are derived using Ar as a tracer much as the vertical eddy mixing coefficient for the Earth's troposphere is derived using CH4 as a minor constituent tracer. The estimation of meridional mixing for high latitudes at Mars is important for constraining parameters used in atmospheric modeling and predicting seasonal and daily behavior. The calculations are order of magnitude estimates that should improve as the data set becomes more robust and improves our models.

Sprague, A. L.↗
NASA NTRS2014

Influence of Mixed Mode I-Mode II Loading on Fatigue Delamination Growth Characteristics of a Graphite Epoxy Tape Laminate

Mixed mode I-mode II interlaminar tests were conducted on IM7/8552 tape laminates using the mixed-mode bending test. Three mixed mode ratios, G(sub II)/G(sub T) = 0.2, 0.5, and 0.8, were considered. Tests were performed at all three mixed-mode ratios under quasi-static and cyclic loading conditions, where the former static tests were used to determine initial loading levels for the latter fatigue tests. Fatigue tests at each mixed-mode ratio were performed at four loading levels, Gmax, equal to 0.5G(sub c), 0.4G(sub c), 0.3G(sub c), and 0.2G(sub c), where G(sub c) is the interlaminar fracture toughness of the corresponding mixed-mode ratio at which a test was performed. All fatigue tests were performed using constant-amplitude load control and delamination growth was automatically documented using compliance solutions obtained from the corresponding quasi-static tests. Static fracture toughness data yielded a mixed-mode delamination criterion that exhibited monotonic increase in Gc with mixed-mode ratio, G(sub II)/G(sub T). Fatigue delamination onset parameters varied monotonically with G(sub II)/G(sub T), which was expected based on the fracture toughness data. Analysis of non-normalized data yielded a monotonic change in Paris law exponent with mode ratio. This was not the case when normalized data were analyzed. Fatigue data normalized by the static R-curve were most affected in specimens tested at G(sub II)/G(sub T)=0.2 (this process has little influence on the other data). In this case, the normalized data yielded a higher delamination growth rate compared to the raw data for a given loading level. Overall, fiber bridging appeared to be the dominant mechanism, affecting delamination growth rates in specimens tested at different load levels and differing mixed-mode ratios.

Ratcliffe, James G.↗
NASA NTRS2014

Results of Propellant Mixing Variable Study Using Precise Pressure-Based Burn Rate Calculations

A designed experiment was conducted in which three mix processing variables (pre-curative addition mix temperature, pre-curative addition mixing time, and mixer speed) were varied to estimate their effects on within-mix propellant burn rate variability. The chosen discriminator for the experiment was the 2-inch diameter by 4-inch long (2x4) Center-Perforated (CP) ballistic evaluation motor. Motor nozzle throat diameters were sized to produce a common targeted chamber pressure. Initial data analysis did not show a statistically significant effect. Because propellant burn rate must be directly related to chamber pressure, a method was developed that showed statistically significant effects on chamber pressure (either maximum or average) by adjustments to the process settings. Burn rates were calculated from chamber pressures and these were then normalized to a common pressure for comparative purposes. The pressure-based method of burn rate determination showed significant reduction in error when compared to results obtained from the Brooks' modification of the propellant web-bisector burn rate determination method. Analysis of effects using burn rates calculated by the pressure-based method showed a significant correlation of within-mix burn rate dispersion to mixing duration and the quadratic of mixing duration. The findings were confirmed in a series of mixes that examined the effects of mixing time on burn rate variation, which yielded the same results.

Stefanski, Philip L.↗
NASA NTRS2016

Impact of Flight Enthalpy, Fuel Simulant, and Chemical Reactions on the Mixing Characteristics of Several Injectors at Hypervelocity Flow Conditions

The high total temperatures or total enthalpies required to duplicate the high-speed flight conditions in ground experiments often place stringent requirements on the material selection and cooling needs for the test articles and intrusive flow diagnostic equipment. Furthermore, for internal flows, these conditions often complicate the use of nonintrusive diagnostics that need optical access to the test section and interior portions of the flowpath. Because of the technical challenges and increased costs associated with experimentation at high values of total enthalpy, an attempt is often made to reduce it. This is the case for the Enhanced Injection and Mixing Project (EIMP) currently underway in the Arc-Heated Scramjet Test Facility at the NASA Langley Research Center. The EIMP aims to investigate supersonic combustion ramjet (scramjet) fuel injection and mixing physics, improve the understanding of underlying physical processes, and develop enhancement strategies and functional relationships between mixing performance and losses relevant to flight Mach numbers greater than 8. The experiments will consider a "direct-connect" approach and utilize a Mach 6 nozzle to simulate the combustor entrance flow of a scramjet engine. However, while the value of the Mach number is matched to that expected at the combustor entrance in flight, the maximum value of the total enthalpy for these experiments is limited by the thermal-structural limits of the uncooled experimental hardware. Furthermore, the fuel simulant is helium, not hydrogen. The use of "cold" flows and non-reacting mixtures of fuel simulants for mixing experiments is not new and has been extensively utilized as a screening technique for scramjet fuel injectors. In this study, Reynolds-averaged simulations are utilized (RAS) to systematically verify the implicit assumptions used by the EIMP. This is accomplished by first performing RAS of mixing for two injector configurations at planned nominal experimental conditions. The mixing parameters of interest, such as mixing efficiency and total pressure recovery, are then computed and compared to the values obtained from RAS under the true enthalpy conditions and using helium and hydrogen. Finally, the impact of combustion on mixing, often deemed small enough to neglect at hypervelocity conditions, is assessed by comparing the results obtained from the hydrogen-fueled reacting and non-reacting RAS. For reacting flows, in addition to mixing efficiency and total pressure recovery, the combustion efficiency and thrust potential are also considered. In all of the simulations, the incoming air Mach number and the fuel-to-air ratio are the same, while the total pressure, total enthalpy, and the fuel simulant vary depending on the case considered. It is found that under some conditions the "cold" flow experiments are a good approximation of the flight.

Drozda, Tomasz G.↗
NASA NTRS2016

CFD Analysis of Mixing Characteristics of Several Fuel Injectors at Hypervelocity Flow Conditions

CFD analysis is presented of the mixing characteristics and performance of three fuel injectors at hypervelocity flow conditions. The calculations were carried out using the VULCAN-CFD solver and Reynolds-Averaged Simulations (RAS). The high Mach number flow conditions match those proposed for the planned experiments conducted as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center. The EIMP aims to investigate scramjet fuel injection and mixing physics, improve the understanding of underlying physical processes, and develop enhancement strategies and functional relationships relevant to flight Mach numbers greater than eight. Because of the high Mach number flow considered, the injectors consist of a fuel placement device, a strut; and a fluidic vortical mixer, a ramp. These devices accomplish the necessary task of distributing and mixing fuel into the supersonic cross-flow albeit via different strategies. Both of these devices were previously studied at lower flight Mach numbers where they exhibited promising performance in terms of mixing efficiency and total pressure recovery. For comparison, a flush-wall injector is also included. This type of injector generally represents the simplest method of introducing fuel into a scramjet combustor, however, at high flight Mach number conditions, the dynamic pressure needed to induce sufficient fuel penetration may be difficult to achieve along with other requirements such as achieving desired levels of fuel-to-air mixing at the required equivalence ratio. The three injectors represent the baseline configurations planned for the experiments. The current work discusses the mixing flow field behavior and differences among the three fuel injectors, mixing performance as described by the mixing efficiency and the total pressure recovery, and performance considerations based on the thrust potential.

Drozda, Tomasz G.↗
DOE OSTI2022

Fluid mixing apparatus using liquid fuel and high- and low-pressure fluid streams

A fluid mixing apparatus includes mixing conduits that extend through a fluid plenum. The fluid plenum, which surrounds a first wall defining a main passage fluidly coupled to a low-pressure fluid source, is surrounded itself by a second wall defining a high-pressure plenum fluidly coupled to a high-pressure fluid source. An insulated tube disposed at the inlet of the first wall delivers a third fluid. The mixing conduits fluidly couple the high-pressure plenum to the main passage, where the high-pressure fluid is mixed with low-pressure fluid and the third fluid. Optionally, the fluid plenum may house a fourth fluid that is injected through injection holes in the mixing conduits. The fluid mixing apparatus may be used to mix one or more fuels with high- and low-pressure air in a gas turbine combustor. Alternately, the fluid mixing apparatus may mix a fluid with high- and low-pressure water streams.

Berry, Joseph Jonathan↗
DOE OSTI2020

On the competition between mixing rate and uniformity in a coaxial hydrodynamic focusing mixer

Fast microfluidic mixers for use with line-of-sight integrating detection schemes pose unique challenges. Such detectors typically cannot discriminate signal from slow moving (e.g. near internal walls) and fast-moving portions of the fluid stream. This convolves reaction rate dynamics with fluid flow residence time dynamics. Further, the small cross sections of typical three-dimensional hydrodynamic focusing devices lead to lower detection signals. The current study focuses on achieving both small time scales of mixing and homogenous residence times. This is achieved by injecting sample through a center capillary and hydrodynamically focusing using a sheath flow within a tapered second capillary. The current design also features a third, larger coaxial capillary. The mixed stream flows into the large cross-section of this third capillary to decelerate and expand the stream by up to 14-fold to improve line-of-sight signal strength of reaction products. Hydrodynamic focusing, mixing, and expansion are studied using analytical and numerical models and also studied experimentally using a fluorescein-iodide quenching reaction. The experimentally validated models are used to explore trade-offs between mixing rate and uniformity. For the first time, this work presents detailed analysis of the Lagrangian time history of species transport during mixing inside coaxial capillaries to measure mixing nonuniformity. The mixing region enables order 100 μs mixing times and residence time widths of the same order (140 μs).

42 ENGINEERING↗
DOE OSTI2025

Evaluation of probability density function descriptions for three-component Rayleigh–Taylor mixing

Results from simulations of a three-component Rayleigh–Taylor (RT) mixing problem are presented. These simulations are conducted in heavy–light–heavy and heavy–intermediate–light configurations, and each of these configurations are further considered in high- and low-Reynolds-number regimes. This results in RT-unstable flow with one or both interfaces initially unstable, permitting the influence of problem configuration on the statistical description of three-component RT-driven mixing to be considered. Mass fraction covariances are observed to undergo a sign change through the mixing layer in all four configurations considered. This appears to be unique to the multi-component case and represents another way in which multi-component RT mixing differs from the two-component case. Qualitative and quantitative comparisons of joint and marginal probability density function (PDF) descriptions of species concentration are made. Three-, five-, and six-parameter model PDFs are compared against simulation data to assess how accurately they describe the mixing, and it is found that three-component mixing requires at least a five-parameter model PDF to accurately describe the mixing. Notably, the marginal distributions of three-component mixing do not appear to conform to a beta distribution, representing a departure from the classical two-component RT case. In conclusion, statistical neutrality also appears to influence the optimal choice of model PDF, which is found to be a function of problem configuration.

Large-eddy simulation↗
DOE OSTI2023

Effects of Lower Troposphere Vertical Mixing on Simulated Clouds and Precipitation Over the Amazon During the Wet Season

Planetary boundary layer (PBL) schemes parameterize unresolved turbulent mixing within the PBL and free troposphere (FT). Previous studies reported that precipitation simulation over the Amazon in South America is quite sensitive to PBL schemes and the exact relationship between the turbulent mixing and precipitation processes is, however, not disentangled. In this study, regional climate simulations over the Amazon in January–February 2019 are examined at process level to understand the precipitation sensitivity to PBL scheme. The focus is on two PBL schemes, the Yonsei University (YSU) scheme, and the asymmetric convective model v2 (ACM2) scheme, which show the largest difference in the simulated precipitation. During daytime, while the FT clouds simulated by YSU dissipate, clouds simulated by ACM2 maintain because of enhanced moisture supply due to the enhanced vertical moisture relay transport process: (a) vertical mixing within PBL transports surface moisture to the PBL top, and (b) FT mixing feeds the moisture into the FT cloud deck. Due to the thick cloud deck over Amazon simulated by ACM2, surface radiative heating is reduced and consequently the convective available potential energy is reduced. As a result, precipitation is weaker from ACM2. Two key parameters dictating the vertical mixing are identified, p , an exponent determining boundary layer mixing and λ , a scale dictating FT mixing. Sensitivity simulations with altered p , λ , and other treatments within YSU and ACM2 confirm the precipitation sensitivity. The FT mixing in the presence of clouds appears most critical to explain the sensitivity between YSU and ACM2.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Shock-driven three-fluid mixing with various chevron interface configurations

When a shock wave crosses a density interface, the Richtmyer–Meshkov instability causes perturbations to grow. Richtmyer–Meshkov instabilities arise from the deposition of vorticity from the misaligned density and pressure gradients at the shock front. In many engineering applications, microscopic surface roughness will grow into multi-mode perturbations, inducing mixing between the fluid on either side of an initial interface. Applications often have multiple interfaces, some of which are close enough to interact in the later stages of instability growth. In this study, we numerically investigate the mixing of a three-layer system with periodic zigzag (or chevron) interfaces, calculating the dependence of the width and mass of mixed material on properties such as the shock timing, chevron amplitude, multi-mode perturbation spectrum, density ratio, and shock mach number. The multi-mode case is also compared with a single-mode perturbation. The Flash hydrodynamic code is used to solve the Euler equations in three dimensions with adaptive grid refinement. Key results include a significant increase in mixed mass when changing from a single-mode to a multi-mode perturbation on one of the interfaces. The mixed width is mainly sensitive to the density ratio and chevron amplitude, whereas the mixed mass also depends on the multi-mode spectrum. In conclusion, steeper initial perturbation spectra have lower mixed mass at early times but a greater mixed mass after the reflected shock transits back across the layer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Cloud microphysical response to entrainment and mixing is locally inhomogeneous and globally homogeneous: Evidence from the lab

Entrainment of dry air into clouds strongly influences cloud optical and precipitation properties and the response of clouds to aerosol perturbations. The response of cloud droplet size distributions to entrainment–mixing is examined in the Pi convection-cloud chamber that creates a turbulent, steady-state cloud. The experiments are conducted by injecting dry air with temperature (T e ) and flow rate (Q e ) through a flange in the top boundary, into the otherwise well-mixed cloud, to mimic the entrainment–mixing process. Due to the large-scale circulation, the downwind region is directly affected by entrained dry air, whereas the upwind region is representative of the background conditions. Droplet concentration (Cn) and liquid water content (L) decrease in the downwind region, but the difference in the mean diameter of droplets (D m ) is small. The shape of cloud droplet size distributions relative to the injection point is unchanged, to within statistical uncertainty, resulting in a signature of inhomogeneous mixing, as expected for droplet evaporation times small compared to mixing time scales. As T e and Q e of entrained air increase, however, Cn, L, and D m of the whole cloud system decrease, resulting in a signature of homogeneous mixing. The apparent contradiction is understood as the cloud microphysical responses to entrainment and mixing differing on local and global scales: locally inhomogeneous and globally homogeneous. This implies that global versus local sampling of clouds can lead to seemingly contradictory results for mixing, which informs the long-standing debate about the microphysical response to entrainment and the parameterization of this process for coarse-resolution models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Mixing Cell Model: A One-Dimensional Numerical Model for Assessment of Water Flow and Contaminant Transport in the Unsaturated Zone

This report describes the Mixing Cell Model code, a one dimensional model for water flow and solute transport in the unsaturated zone under steady state or transient flow conditions. The model is based on the principles and assumptions underlying mixing-cell model formulations. The unsaturated zone is discretized into a series of independent mixing cells. Each cell may have unique hydrologic, lithologic, and sorptive properties. Ordinary differential equations describe the material (water and solute) balance within each cell. Water-flow equations are derived from the continuity equation, assuming that unit gradient conditions exist at all times in each cell. Pressure gradients are considered implicitly through model discretization. Unsaturated hydraulic conductivity and moisture contents are determined by the material specific moisture-characteristic curves. Solute-transport processes include explicit treatment of advective processes, first order chain decay, and linear sorption reactions. Dispersion is addressed through implicit and explicit dispersion. Implicit dispersion is an inherent feature of all mixing-cell models and originates from the formulation of the problem in terms of mass balance around fully mixed volume elements. Expressions are provided that relate implicit dispersion to the physical dispersion of the system. Two FORTRAN codes were developed to solve the water flow and solute-transport equations: (1) the Mixing Cell Model for Flow (MCMF) solves transient water-flow problems and (2) the Mixing Cell Model for Transport (MCMT) solves the solute-transport problem. The transient water-flow problem is typically solved first by estimating the water flux through each cell in the model domain as a function of time using the MCMF code. These data are stored in either ASCII or binary files that are later read by the solute transport code (MCMT). Code output includes solute pore water concentrations, water and solute inventories in each cell and at each specified output time, and water and solute fluxes through each cell and specified output time. Computer run times for coupled transient water flow and solute transport were typically several seconds on a 2 GHz Intel Pentium IV desktop computer. The model was benchmarked against analytical solutions and finite element approximations to the partial differential equations (PDE) describing unsaturated flow and transport. Differences between the maximum solute flux estimated by the mixing cell model and the PDE models were typically less than two percent. This revision includes an option for a fixed concentration lower boundary condition for diffusive fluxes for versions 020321 and later.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Quantifying the Turbulent Entrainment‐Mixing Processes Based on Z ‐LWC Relationships of Cloud Droplets

Abstract Turbulent entrainment‐mixing processes profoundly influence the relationship between radar reflectivity factor and liquid water content ( Z ‐LWC) of cloud droplets. However, quantification of the entrainment‐mixing mechanisms based on the Z ‐LWC relationship is still lacking. To address this gap, 12,218 entrainment‐mixing cases are simulated using the Explicit Mixing Parcel Model. We examine the variations of the parameters in the power‐law relationship Z = a LWC b , and the relationship between parameter b and homogeneous mixing degree ( ψ ), a measure quantifying entrainment‐mixing processes. The results indicate that parameter b distributes within the range of 1–2, with a positive correlation between parameter b and ψ . The b ‐ ψ relationship is fitted, which connects the Z ‐LWC relationship for various entrainment‐mixing types. The results suggest the potential for employing a remote sensing approach to investigate the entrainment‐mixing mechanisms of non‐precipitating small cumulus/stratocumulus clouds, thereby overcoming the limitations of traditional observational studies that rely solely on aircraft observations.

54 ENVIRONMENTAL SCIENCES↗