Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “drops”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Experimental Investigation of Vapor Formation in Liquid CO2 Flow Through a Converging-Diverging Nozzle

Carbon dioxide is an attractive working fluid for many cycles, including for pumped thermal energy storage (PTES). A challenge with some proposed sCO2 PTES cycles is the operation of sCO2 machinery outside the typical bounds of experience, with local phase change from the liquid state being particularly unknown. Presently, there is insufficient data in the literature regarding multiphase CO2 to adequately design a multiphase-tolerant turbine, so generation of foundational data is required. This experimental study investigates the flow characteristics of sub-sonic liquid CO2 undergoing expansion and phase change in a converging-diverging nozzle. The nozzle is instrumented to measure static pressure, unsteady pressure, temperature, and density. The static pressure transducers are located at 27 axial locations to accurately characterize the pressure profile in the nozzle. High-accuracy RTDs are located at the entrance and exit of the nozzle, and three dynamic pressure transducers are strategically located to capture any unsteady phenomena. During testing, values of mass flow and nozzle inlet pressure are swept to vary the pressure drop and fluid properties. The measured total pressure drop in the nozzle is compared to a homogenous model and the Lockhart-Martinelli correlation method, with the latter predicting loss quite closely. The resulting data set is valuable for validating multiphase numerical models in a simple geometry before implementation of these models in turbomachinery design.

25 ENERGY STORAGE↗
DOE OSTI2023

ECAR-6581 Rev 0 Hydraulic Assessment of Conceptual Design of Ring-type Spacer for MARVEL Reactor

The fuel rods in the nuclear reactor core can be deformed due to various reasons, such as thermal expansion, the effect of irradiation, hydrodynamic effects, etc. The deformation of the fuel rods can lead to contact between them, which can deteriorate the coolability of the nuclear fuel and increase the potential risk of failure. To prevent the fuel rods from coming into contact with each other, a ring-type spacer may be introduced in the MARVEL reactor. However, the ring spacers can result in an additional pressure drop, which could restrict the natural circulation operation of the MARVEL reactor. Therefore, this ECAR will examine the additional pressure drop caused by the introduction of the ring-type spacer by comparing with the bare rod condition.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2026

Next Generation Refrigerants for Domestic Heat Pump Water Heaters

To select next generation refrigerants and phase out R-134a, we conducted extensive drop-in tests of multiple new refrigerants. For laboratory investigations, we performed drop-in UEF (uniform energy factor) and FHR (first hour rating) tests on two heat pump water heaters (HPWHs), one uses a 220V compressor and 40-gallon water tank, and the other uses a 120V compressor and 50-gallon tank. The alternative refrigerants/blends include R-1234yf, R-1234ze, R-515B, R-513A, and R-516A, among which R-513A is A1, non-flammable and the others are in A2L category. All the alternative refrigerants required similar optimum system charges as the baseline R-134a. They resulted in similar or better UEFs than R-134a. R-516A led to best UEFs in both the HPWHs. The resultant discharge temperatures of the new refrigerants/blends are lower than R-134a, indicating good lubricant compatibility. The 120V HPWH delivered > 60 gallons first-hour rating using R-134a. The alternatives led to reduced FHRs due to their smaller volumetric vaporization capacities using the same compressor. Except R-515B, the other refrigerants/blends delivered FHRs > 60 gallons. R-516A produced better UEFs and FHRs in both the units than the other replacements. If considering flammability, R-513A would be the best option, which has similar UEF and FHR as R-134a.

Shen, Bo [ORNL] (ORCID:0000000336600393)↗
DOE OSTI2026

Thermal-Hydraulic Performance of a Zeotropic Refrigerant in a Brazed Plate Evaporator

In this work we report the results of a preliminary investigation focused on evaluating the thermal-hydraulic performance of R455A (CO2, R32, R1234yf – 3/21.5/75.5) in a commercial brazed plate heat exchanger. Three key parameters were measured which were refrigerant flow maldistribution, overall refrigerant side heat transfer coefficient, and total refrigerant side pressure drop. These same parameters were also measured for R134a which served as a base-line refrigerant for comparison. The experimental test matrix was set up to investigate the influence of refrigerant mass flux (3 ≤ Gref ≤ 5 kg m-2 s-1), refrigerant inlet vapor quality (0.1 ≤ xin ≤ 0.4), and refrigerant exit superheat (5 ≤ Tsh ≤ 10 K) on the above mentioned thermal-hydraulic parameters. For these experimental conditions, on average, the overall refrigerant side heat transfer coefficients for R455A were 45.2% lower than those of R134a. Refrigerant flow maldistribution for R455A, on average, was observed to be 62.6 % higher than that measured for R134a, and the total pressure drop across the heat exchanger (minor and frictional losses) for R455A was on average 26.5% higher than that measured for R134a. These initial results set the stage for a comprehensive investigation focused on further investigating the thermal-hydraulic performance of zeotropic refrigerants in brazed plate heat exchangers.

Jajja, Saad [ORNL]↗
DOE OSTI2021

Contextualizing Ancestral Pueblo Turkey (Meleagris gallopavo spp.) Management

Abstract Penning turkeys ( Meleagris gallopavo spp.) in the Ancestral Pueblo American Southwest/Mexican Northwest (SW/NW) involved the creation or use of a variety of spaces and contexts throughout AD 1–1600 and into the post-contact era. Turkey pens, or captivity, occur through simple tethering, reuse of abandoned pit houses or surface rooms, or creation of pens within villages, plazas, and elsewhere. Turkey dung, droppings, and eggshells are fundamental for determining the presence or absence of pens at archaeological sites. In this paper, I review the archaeological record for turkey pens and focus on three main questions: (1) how are turkey pens identified in the SW/NW, (2) if turkey pen construction or evidence for turkey captivity shifts through time, and (3) what the record of turkey penning informs us regarding the long-term human management of these birds and global perspectives on human–bird/human–animal management. Ancestral Pueblo peoples created an adaptive and flexible strategy for turkey penning, which successfully integrated these birds into ceremonial and socioeconomic processes for approximately 1600 years.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

A Novel Design of Guiding Stress Wave Propagation

Impact loads can induce a series of undesirable physical phenomena including vibration, acoustical shock, perforation, fracture and fragmentation, etc. The energy associated with the impact loads can lead to severe structure damage and human injuries. A design approach which effectively reduces these negative impacts through shock/stress wave diversion is highly needed. In this paper, a computational model which predicts stress wave propagation by considering different beam geometries and configurations is developed. Additionally, a novel concept of wave guide design which modifies the stress wave propagation path without disturbance is also presented. Furthermore, this design approach is not only useful for material property characterization particularly at intermediate or high strain rates, but also allows stress wave propagation in a desired direction as the shock/impact energy can be redistributed in controllable paths. The numerical results are experimentally verified through a Drop-Hopkinson bar apparatus at Sandia National Laboratories.

36 MATERIALS SCIENCE↗
DOE OSTI2020

High Speed In-situ X-ray Imaging of 3D Freeze Printing of Aerogels

3D freeze printing (3DFP) combines drop-on-demand (DOD) inkjet printing with freeze casting to fabricate lightweight and multifunctional aerogels with customized geometries. Freeze casting is an efficient and easily implemented method capable of fabricating porous, sponge-like structures for many different applications. This process enables tailoring the microstructure of the final product (i.e., pore morphology, alignment, average size distribution, etc.) by controlling the fabrication conditions and freezing kinetics. Furthermore, its combination with DOD printing provides the capability of engineering the macrostructure without relying on a mold as reported for 3D freeze-printed aerogels made from graphene, silver nanowires, and other nanocomposites. In this paper, we performed in-situ X-ray imaging to understand the inside process dynamics in 3DFP using a commercially available colloidal silica ink. We investigated the 3DFP process with the following hierarchy: first, single droplets; then, uniform lines obtained from coalescence of droplets; and finally, three consecutive lines deposited layer by layer. With the help of X-ray imaging, the importance of the balance between material deposition and freezing rates was shown in-situ by the observation inside of the freeze front following the tip of the printed line. The effects of the substrate temperature on the elimination of undesired interfacial boundaries were also shown by the observed ice crystals penetrating from lower to upper layer.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Sustainable co-production of plant lipids and cellulosic sugars from transgenic energycane at an industrially relevant scale: A proof of concept for alternative feedstocks

Development of sustainable and scalable technologies to convert lignocellulosic biomass to biofuels is critical to achieving carbon neutrality. The potential of transgenic bioenergy crops as a renewable source of sugars and lipids has been demonstrated at bench-scale. However, scaling up these processes is important for holistic analysis. Here proof-of-concept for chemical-free hydrothermal pretreatment of transgenic energycane-oilcane line L13 at an industrially relevant scale to recover vegetative lipids along with cellulosic sugars is presented. Here, pilot-scale processing of 97 kg of transgenic energycane-oilcane L13 stems and high solids pretreatment of bagasse enhanced the recovery of cellulosic glucose and xylose by 5-fold as compared to untreated bagasse and helped in the enrichment of vegetative lipids in the biomass residues which allowed its recovery at the end of the bioprocess. Palmitic and oleic acids were the predominant fatty acids (FAs) extracted from stems and leaves. The processing did not affect lipid composition. The efficiency of lipid recovery from untreated biomass was 75.9% which improved to 88.7% upon pretreatment. The vegetative tissues of transgenic energycane-oilcane L13 contained 0.42 metric tons/hectare of lipids. Processing vegetative tissues yielded 0.38 metric tons/hectare of lipids. This approaches an oil yield similar to soybean (global average 0.44 metric tons/hectare) and is almost twice as high as the oil yield from sugarcane engineered to hyperaccumulate lipids (0.20 metric tons/hectare). The study suggests that further optimization by state-of-the-art metabolic engineering and biomass processing can establish transgenic bioenergy crops for commercial drop-in fuel production.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Bio-inspired alula-based winglet design for enhanced heat transfer in high temperature fin-and-tube heat exchangers

Fin-and-tube heat exchangers (FTHEs) are widely used for high-temperature flue-gas heat recovery, but their performance is often limited by wake regions and non-uniform fin-surface temperatures. This study proposes and numerically evaluates four bio-inspired longitudinal vortex generator (VG) configurations in a high-temperature FTHE with flue-gas inlet temperature ∼1230 K: double-delta, curved double-delta, alula, and a new curved-alula geometry. The reference fin is not hydraulically plain; it already incorporates leading-edge separation columns and convex protrusions, so the alula-type winglets are assessed as downstream add-ons acting on a strongly disturbed flow. In a second step, perforations (one, two and three circular holes) are introduced into the curved-alula VGs to further tailor the flow field. Three-dimensional simulations with the Shear Stress Transpor (SST) $k - ω$ model, temperature-dependent flue-gas properties and conjugate conduction are carried out for gas-side Reynolds numbers $Re_g ≈ 8.0$ x $10^2 - 3.6$ x $10^3$ (mass flow rates 0.5 – 2.5 g/s), and the designs are compared in terms of surface heat flux, Nusselt number, friction factor and hydrothermal performance factor (HTPF). For this already-promoted fin, the additional downstream winglets provide moderate, incremental hydrothermal gains. At the highest Reynolds number, the best non-perforated design (curved-alula) increases surface heat flux from 1630.9 to 1794.7 kW/m² (∼ 10 % gain) and the Nusselt number from 227.6 to 242.6 (∼ 7 % gain), while the friction factor rises from 0.26 to about 0.30, yielding HTPF values close to unity (∼ 0.9 – 1.0). Introducing circular perforations into the curved-alula winglets acts mainly as a wake-bleeding refinement: the three-hole configuration provides a heat flux of 1824.7 kW/m² and a pressure drop of 127.9 Pa, with HTPF in the range ∼ 1.03 – 1.14 and a small (∼ 1 – 3 %) improvement over the solid curved-alula design. Flow-field analysis shows that the perforated curved-alula VGs shrink tube-wake regions, thin the thermal boundary layer and homogenize the fin-surface temperature (outlet-gas temperature ∼ 510 – 520 K and fin-surface temperature ∼ 420 – 421 K for the three-hole case). An optimal flue-gas mass flow rate of ∼ 1 g/s ($Re_g ≈ 1.5$ x $10^3$) is identified, beyond which additional heat-transfer gains are offset by rapidly increasing pressure losses. Overall, the results highlight that initial fin geometry and VG placement are as important as VG shape: alula-based winglets are expected to yield larger relative gains on simpler flat-fin layouts or when positioned closer to the fin leading edge and tube

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

The Corrected Distortion model for Lagrangian spray simulation of transcritical fuel injection

In this work, we present a detailed implementation and validation of the droplet modeling framework proposed by Dahms and Oefelein (2016) into the engine commercial CFD software CONVERGE using the User Defined Function (UDF) interface. The model accounts for the nonlinear deformation and oscillation experienced by liquid spray droplet injected into high pressure and temperature. Lagrangian spray simulations of Engine Combustion Network (ECN) Spray A are performed. Model validation against standard experimental measurements of liquid velocity, vapor mixture fraction is conducted. To perform more rigorous model validation, new experimental measurements based on Diffused Back Illumination (DBI) are introduced. The new measurements are processed for Projected Liquid Volume (PLV), which offers as close as possible one-to-one model validation for liquid penetration while offering new insights into the spray physics. Comparison with a One-D model based on adiabatic mixing theory by Siebers (1999) and Desantes et al. (2007) are also conducted. Through these model validation exercises, it is shown that the new framework improves liquid-phase penetration predictions, following a tendency for enhanced evaporation, compared to the standard approach for both Reynolds Average Navier Stokes (RANS) and Large Eddy Simulation (LES). Additionally, at the liquid length, maximum mixture fraction values predicted by the new approach are in good agreement those of an adiabatic mixing model. Qualitative analysis of the spray behaviors during the early stage of the injection process reveals that the proposed framework predicts significant increase in droplet evaporation rate with lower drop drag compared to the current standard approach.

42 ENGINEERING↗
DOE OSTI2024

Pyrolysis of high-density polyethylene: Degradation behaviors, kinetics, and product characteristics

Pyrolysis is a promising technology for converting plastic waste into valuable raw materials while offering a potential solution to the global plastic pollution crisis. In this study, the thermal pyrolysis of high-density polyethylene (HDPE) is investigated in a drop tube reactor under nearly isothermal conditions. The impact of reaction temperature and gas/volatile residence time on carbon conversion and product distribution is examined across a range of 500–900°C and 3.6–32.2s, respectively. Non-condensable gas products detected by online mass spectrometry are H 2 , CH 4 , C 2 H 4 , C 2 H 6 , C 3 H 6 , and C 3 H 8 . At elevated temperatures and prolonged residence time, H 2 yield reaches as high as 8.6 wt% of the initial HDPE mass due to intensified cracking reactions of C 2 –C 3 hydrocarbons and long-chain aliphatic compounds. Consequently, pyrolysis tars consist mainly of polycyclic aromatic hydrocarbons (PAHs) with 5–7 rings, accompanied by visible coke deposition within the reactor. HDPE decomposition to volatiles is an endothermic process and it is complete at a temperature between 492°C and 525°C, depending on the heating rate employed, from non-isothermal thermogravimetric analysis and differential scanning calorimetry (TGA-DSC) measurements. The thermal degradation of HDPE pellets follows the two-dimensional nucleation growth model for conversion levels up to 0.8 with an apparent activation energy of 259–270 kJ/mol and a pre-exponential factor of 4.83 × 10 17 –1.37 × 10 19 min -1 , determined from various isoconversional methods such as Flynn-Wall-Ozawa (FWO), Kissinger-Akahira-Sunose (KAS), and Starink, along with Criado's master plots. Further, these findings provide valuable insights into optimizing process parameters and refining reactor design for pyrolysis, which can be integrated with gasification and reforming processes to enhance hydrogen production on a larger scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Review of heat transfer enhancement techniques for single phase flows

The thermal energy exchange between a flowing fluid and its confining channel is a ubiquitous process in modern society. To enhance the fluid-to-wall or wall-to-fluid heat transfer, several techniques have been developed to maximize the contact area between the fluid and the inner wall and/or disrupt the flow to enhance circulation or induce turbulence. Deployment of channels having features capable of enhancing heat transfer enables the reduction of heat exchanger size while maintaining performance. Reduction in equipment size is critical due to the ability to minimize the required volume of costly working fluids and to mitigate potential safety concerns associated with total system fluid volume. Here, a comprehensive review of single-phase heat transfer enhancement techniques is presented. The article provides a thorough comparison by analyzing the heat transfer rate, pressure drop, and other operational aspects. Single-phase heat transfer enhancement methods are divided into active and passive techniques. Active methods such as electrohydrodynamic (EHD), magnetohydrodynamics (MHD), or mechanical motion require external power to create enhancement. Passive methods such as dimples, fins, or tape inserts do not require external input and rely only on surface modification. Although active methods are more expensive and difficult to implement compared to passive techniques, it enables active control of heat transfer augmentation. Finally, this review develops and summarizes key learning data for design optimization enabled by additive manufacturing and machine learning algorithms, helping to inform these next-generation heat exchanger design methodologies for a plethora of modern applications such as electrification of vehicles, computing, and classical industries.

42 ENGINEERING↗
DOE OSTI2024

Flow maldistribution in plate heat exchangers – Impact, analysis, and solutions

Plate heat exchangers are commonly used in various industrial applications, such as refrigeration, air conditioning, heat pumps, powerplants, and chemical industries. Plate heat exchangers are well known for their superior heat transfer performance, compactness, and low refrigerant charge. Despite offering several advantages, they suffer from flow maldistribution issues. The flow maldistribution can deteriorate both the heat transfer and pressure drop performance, ultimately resulting in a lower system efficiency where plate heat exchangers are deployed. The flow maldistribution issues become more pronounced when the heat exchanger size is relatively large and the number of plates is higher, limiting the deployment of plate heat exchangers in larger industrial systems. Consequently, analyzing and understanding the flow maldistribution behavior in plate heat exchangers and finding ways to mitigate flow maldistribution related issues become essential topics of interest. Further, this review aims to address the effect of flow maldistribution on plate heat exchanger characteristics. First, the experimental and numerical works on flow maldistribution under single-phase and two-phase conditions are detailed. Subsequently, the end-channel and end-plate effects are discussed. Then, the methods to mitigate flow maldistribution in plate heat exchangers are outlined. Finally, based on a thorough literature survey and industrial requirements, future research directions are recommended.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Machine learning-based optimization of air-cooled heat sinks

Machine learning-based models using Artificial Neural Network (ANN) and greedy search algorithm are used to optimize air-cooled parallel plate-finned heat sinks (PPFHSs) subjected to laminar flow over an extensive range of design parameters. Here, the thermal and hydraulic performances of PPFHSs are represented by heat transfer coefficient (h) and pressure drop (ΔP), respectively. Optimization objectives for PPFHS designs can vary from industry to industry depending on their design priorities. The present study proposes a novel and generalized optimization method that defines practical optimization objectives and provides an accurate optimization process to design effective PPFHSs for a wide range of industrial applications with different design requirements. Three optimization objectives are presented in this study: (i) the largest h PΔ, (ii) the largest h within a specified maximum allowed flow rate, and (iii) the lowest weight that maximizes h for operation within the maximum allowed flow rate. While the shortcoming of the first objective is demonstrated, the other two objectives are found to be suitable for designing effective heat sinks (HSs) across different applications. Results suggest a promising trend from the third objective to develop HSs with ~ 37-68% lower weight, 80-85% reduced ΔP, and negligible penalty in h compared with optimized HSs obtained from the second objective. However, since the third objective leads to HSs with thinner fins, structural analysis should be performed to ensure reliable operation of the HSs.

42 ENGINEERING↗
DOE OSTI2023

Quantifying Site Effects and Their Influence on Earthquake Source Parameter Estimations Using a Dense Array in Oklahoma

We investigate the effects of site response on source parameter estimates using earthquakes recorded by the LArge-n Seismic Survey in Oklahoma (LASSO). While it is well known that near-surface unconsolidated sediments can cause an apparent breakdown of earthquake self-similarity, the influence of laterally varying site conditions remains unclear. We analyze site conditions across the 1825-station array on a river plain within an area of 40 km by 23 km using vertical ground motions from 14 regional earthquakes. While the source radiation pattern controls P-wave ground motions below 8 Hz, the surface geology correlates with P-wave ground motions above 8 Hz and S-wave ground motions at 2–21 Hz. Stations installed in alluvial sediments have vertical ground motions that can exceed three times the array median. We use the variation of ground motion of regional earthquakes across the array as a proxy for site effects. The corner frequencies and stress drops of local earthquakes ( M L = 0.01–3) estimated using a standard single-spectra approach show negative correlations with the site-effect proxy, while the seismic moments show positive correlations. In contrast, the spectral-ratio approach effectively shows no correlation. The overall bias is small as expected for this relatively homogeneous structure; accurate estimation of site-related biases requires at least 30 stations. Correcting for site-related biases reduces the standard deviations of the source parameters by less than 13% of the total variations. Remaining variations are partially associated with source directivity and model misfits— as small earthquakes can have complex ruptures.

58 GEOSCIENCES↗
DOE OSTI2025

Integrative path modeling and QTL mapping identify maturity, stem strength, and cell wall composition driving lettuce resistance to Sclerotinia minor

Lettuce ( Lactuca sativa ) is highly vulnerable to Sclerotinia minor , the pathogen causing lettuce drop. Breeding for resistance is the most effective control strategy; however, full resistance has not been achieved, and current partial resistance sources are often linked with undesirable traits, such as early bolting. This study aimed to unravel the genetic basis of partial resistance to S. minor and its relationship with plant maturity (bolting), stem mechanical strength (SMS), and cell wall composition (CWC) using a recombinant inbred line (RIL) population derived from a cross between the susceptible iceberg cv. ‘Salinas’ and the resistant oil-seed accession PI 251246. Field evaluations indicated that resistance was linked to earlier bolting, stronger stems, and higher pentose content. Path analysis demonstrated that earlier-maturing plants exhibited increased resistance through enhanced SMS and modified CWC, particularly with higher xylose and lower arabinose levels. Further analysis indicated a significant relationship between syringyl lignin content and resistance, especially in plants with varying bolting responses. Three key quantitative trait loci (QTLs) on linkage groups (LG) 2, 6, and 7 were consistently associated with resistance, bolting, and SMS. Importantly, residual QTL analysis revealed that the resistance locus on LG7 acted independently of maturity, suggesting a distinct resistance mechanism. Callose synthase emerged as a key candidate gene within the LG7 resistance QTL, located near - but distinct from - genes associated with plant maturity and flowering. These findings provide valuable insights into decoupling resistance from early bolting, suggesting a pathway for breeding lettuce cultivars with improved disease resistance and delayed bolting.

Lactuca↗
DOE OSTI2025

High-Fidelity CFD Assessments of Flow Resistance in a 61-Pin Wire-Wrapped Assembly with Partially Blocked Channels

The examination of thermal-hydraulic behaviors in wire-wrapped rod bundles continues to be an active area of research. The sodium fast reactor, a prominent candidate in next-generation nuclear designs, utilizes a hexagonal configuration of wire-wrapped fuel pins. Here, the potential for channel blockage within this compact arrangement poses a significant safety challenge, spurring a number of recent experimental and computational investigations to evaluate its impact on flow and heat transfer. The present work aims to benchmark the high-fidelity NekRS computational fluid dynamics (CFD) solver in predicting the pressure drops associated with substantial blockages, using available experimental data as a reference. A 61-pin wire-wrapped fuel assembly with two flow blockage configurations has been simulated and investigated at a range of low to moderate Reynolds numbers (487 ≤ Re ≤ 14 600). The NekRS solver demonstrates an exponential reduction of spatial discretization error with increasing polynomial order. The high level of agreement between the numerical results and measured data confirms the accuracy and consistency of the present numerical approach. This benchmark study establishes the capability of NekRS to perform reliable hydrodynamic simulations for sodium fast reactor applications and supports its use in design, licensing, and safety analyses.

CFD Benchmarking↗
DOE OSTI2022

Limitations of Bin and Bulk Microphysics in Reproducing the Observed Spatial Structure of Light Precipitation

Abstract Coarse-gridded atmospheric models often account for subgrid-scale variability by specifying probability distribution functions (PDFs) of process rate inputs such as cloud and rainwater mixing ratios ( q c and q r , respectively). PDF parameters can be obtained from numerous sources: in situ observations, ground- or space-based remote sensing, or fine-scale modeling such as large-eddy simulation (LES). LES is appealing to constrain PDFs because it generates large sample sizes, can simulate a variety of cloud regimes/case studies, and is not subject to the ambiguities of observations. However, despite the appeal of using model output for parameterization development, it has not been demonstrated that LES satisfactorily reproduces the observed spatial structure of microphysical fields. In this study, the structure of observed and modeled microphysical fields are compared by applying bifractal analysis, an approach that quantifies variability across spatial scales, to simulations of a drizzling stratocumulus field that span a range of domain sizes, drop concentrations (a proxy for mesoscale organization), and microphysics schemes (bulk and bin). Simulated q c closely matches observed estimates of bifractal parameters that measure smoothness and intermittency. There are major discrepancies between observed and simulated q r properties, though, with bulk simulated q r consistently displaying the bifractal properties of observed clouds (smooth, minimally intermittent) rather than rain while bin simulations produce q r that is appropriately intermittent but too smooth. These results suggest fundamental limitations of bulk and bin schemes to realistically represent higher-order statistics of the observed rain structure.

54 ENVIRONMENTAL SCIENCES↗