Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “edge effects”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A Large Eddy Simulation Study of Flow Turbulence, Alumina Transport, and Bath Temperature Evolution in Conventional Aluminum-Smelting Cell Using OpenFOAM

In this study, a Large Eddy Simulation (LES) of the aluminum-smelting process is performed using OpenFOAM. To understand the coupled behavior of heat transfer, mass transfer, and flow of the smelting process, a multi-physics computational fluid dynamics (CFD) model based on the Eulerian–Eulerian multi-fluid approach is adopted. The model accounts for CO 2 bubble and magnetohydrodynamics (MHD)-driven flow, along with alumina dissolution, transport, and bath temperature evolution. The simulation predictions show small-scale turbulent vortical structures in the anode–cathode space caused by combined effect of MHD and CO 2 bubble-bath interactions and relatively large-scale asymmetric vortices in the inter-anode space caused by the CO 2 bubble-bath interactions. The vortex formation at the edges of the anodes evidently aids in transporting alumina from the central channel to the bottom of the anodes and prevents accumulation of gas bubbles in the periphery of the anode bottom. Symmetric bath cold spots are observed in the vicinity of the feeder. Cold spots are also observed in the anode–cathode distance space below the anode bottom due to the transport of undissolved solid to this region by the flow. The findings from the work are useful in developing and designing alumina-feeding strategy leading to reduced anode effects and smooth operation of the cell. The work also highlights the important flow structures in conventional aluminum-smelting cell.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Toward Drilling the Perfect Geothermal Well: An International Research Coordination Network for Geothermal Drilling Optimization Supported by Deep Machine Learning and Cloud Based Data Aggregation

The EDGE project, supported by the U.S. Department of Energy Geothermal Technologies Office under award DE-EE0008793, established a data-driven framework for improving the efficiency, cost-effectiveness, and reliability of geothermal well drilling. The project focused on developing scalable data infrastructure, advanced machine learning and probabilistic models, and integrated analytics tools to support continuous drilling optimization. A central objective was to reduce geothermal drilling costs by up to seventy percent while minimizing the risk of well failure through predictive diagnostics and adaptive planning. Over the project period, a comprehensive data repository was designed and deployed, incorporating records from over one hundred geothermal wells across varied geological settings. This repository supported both structured and unstructured data and adhered to FAIR data principles, enabling provenance tracking, quality control, and standardized metadata. The project introduced automated ingestion pipelines and a cloud-hosted platform that facilitated access to raw, processed, and derived datasets. This infrastructure served as the foundation for model development and analysis. Machine learning workflows were developed to predict key drilling metrics including rate of penetration, non-productive time, and total drilling costs. Self-organizing maps and dimensionality reduction methods were used to uncover operational patterns and outliers, while supervised learning algorithms such as random forests and deep neural networks were applied to forecast performance outcomes. The models were validated on heterogeneous datasets from both U.S. and Icelandic fields, demonstrating variable but significant predictive accuracy. The results indicated that finer temporal resolution, inclusion of lithological data, and consistency in operational annotations could substantially improve model performance. The project also implemented process mining techniques to reconstruct state-transition models from drilling event logs. These models enabled the identification of deviations from optimal workflows and provided insights into recurring failure modes. Analysis of non-productive time highlighted the impact of equipment failures, geological challenges, and human factors, offering opportunities for targeted mitigation strategies. The EDGE Dashboard was developed as a web-based expert system integrating data visualization, model outputs, and user-driven queries. It provided an accessible interface for operators to explore historical data, evaluate predicted outcomes, and compare drilling scenarios. Initial feedback from project partners suggested that the dashboard could serve as a foundation for more advanced advisory and optimization tools. Overall, the EDGE project demonstrated the feasibility and value of applying modern data science techniques to geothermal drilling. It delivered a set of interoperable tools and models that can support more efficient, lower-risk well development. The findings point toward a viable path for transitioning from advisory analytics to semi-autonomous drilling systems, contingent on continued collaboration, expanded datasets, and field validation. The project results have immediate relevance for drilling operations, data management practices, and future geothermal R&D efforts aimed at achieving reliable, cost-competitive geothermal energy at scale.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Spin Solar Cell Phenomenon on a Single-Molecule Magnet (SMM) Impacted CoFeB-Based Magnetic Tunnel Junctions

The single-molecule magnet (SMM) is demonstrated here to transform conventional magnetic tunnel junctions (MTJs), a memory device used in present-day computers, into solar cells. For the first time, we demonstrated an electronic spin-dependent solar cell effect on an SMM-transformed MTJ under illumination from unpolarized white light. We patterned cross-junction-shaped devices to form a CoFeB/MgO/CoFeB-based MTJ. The MgO barrier thickness at the intersection between the two exposed junction edges was less than the SMM extent, which enabled the SMM molecules to serve as channels to conduct spin-dependent transport. The SMM channels yielded a region of long-range magnetic ordering around these engineered molecular junctions. Our SMM possessed a hexanuclear [Mn 6 (μ 3 -O) 2 (H 2 N-sao) 6 (6-atha) 2 (EtOH) 6 ] [H 2 N-saoH = salicylamidoxime, 6-atha = 6-acetylthiohexanoate] complex and thioesters end groups to form bonds with metal films. SMM-doped MTJs were shown to exhibit a solar cell effect and yielded ≈80 mV open-circuit voltage and ≈10 mA/cm 2 saturation current density under illumination from one sun equivalent radiation dose. A room-temperature Kelvin Probe AFM (KPAFM) study provided direct evidence that the SMM transformed the electronic properties of the MTJ’s electrodes over a lateral area in excess of several thousand times larger in extent than the area spanned by the molecular junctions themselves. As a result, the decisive factor in observing this spin-photovoltaic effect is the formation of SMM spin channels between the two different ferromagnetic electrodes, which in turn is able to catalyze the long-range transformation in each electrode around the junction area.

14 SOLAR ENERGY↗
DOE OSTI2023

Neutral Atomic Beam Techniques for Diagnosis of Edge and Scrape-Off Layer Plasmas (2017-2019) [and subsequently titled] Neutral Atomic Beam Techniques for Diagnosis of Interior and Scrape-Off Layer Plasmas (2019-2022) (Final Report)

This is the Final Technical Report for the DOE Measurement Innovation grant DE-SC0017998 titled Neutral Atomic Beam Techniques for Diagnosis of Edge and Scrape-Off Layer Plasmas (2017-2019) and subsequently titled Neutral Atomic Beam Techniques for Diagnosis of Interior and Scrape-Off Layer Plasmas (2019-2022). The grant was funded by the Office of Fusion Energy Sciences (FES) for the period 1 September 2017 – 31 May 2022. Work performed through the grant has resulted in diagnostic innovations and new measurement capabilities that could improve understanding of magnetic equilibrium, edge localized modes in tokamaks, and bootstrap and Pfirsch-Schlüter currents in stellarators. We have performed simulations of a novel non-hydrogen beam-based diagnostic that would determine a local value (in the plasma) of the poloidal flux function by measuring the toroidal angle of the beam particles' velocity at a detector external to the plasma. The approach is similar to that of a Heavy Ion Beam Probe (HIBP), but would be able to operate with a smaller, more economical detection systems and, in doing so, may enable substantial extension of scenarios in which beam-based diagnostics are deployed. Key results include: Identification of the neutral and singly charged atoms that are most appropriate for use in the beam based diagnostic and the ionization cross-sections that are important for studying the plasma edge of a large device; Development of a promising geometry for operation in which the beam is injected horizontally near the torus midplane; Implementation of simulations using banana orbit trajectories of secondary beam ions to access plasma regions that are commonly in-accessible to an HIBP diagnostic, including sample volumes near the X-point of a diverted tokamak; Assessment of the ability of the diagnostic concept (originally developed under the assumption of, and reliant on, symmetry) to work in stellarators by performing simulations with both quasi-symmetric fields and non-quasi-symmetric fields; Measurement of noise levels with a prototype detector on a stellarator and exploration of techniques to minimize effects of noise on a diagnostic.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Physics-informed heterogeneous graph neural networks for DC blocker placement

The threat of geomagnetic disturbances (GMDs) to the reliable operation of the bulk energy system has spurred the development of effective strategies for mitigating their impacts. One such approach involves placing transformer neutral blocking devices, which interrupt the path of geomagnetically induced currents (GICs) to limit their impact. The high cost of these devices and the sparsity of transformers that experience high GICs during GMD events, however, calls for a sparse placement strategy that involves high computational cost. To address this challenge, we developed a physics-informed heterogeneous graph neural network (PIHGNN) for solving the graph-based dc-blocker placement problem. Our approach combines a heterogeneous graph neural network (HGNN) with a physics-informed neural network (PINN) to capture the diverse types of nodes and edges in ac/dc networks and incorporates the physical laws of the power grid. We train the PIHGNN model using a surrogate power flow model and validate it using case studies. Results demonstrate that PIHGNN can effectively and efficiently support the deployment of GIC dc-current blockers, ensuring the continued supply of electricity to meet societal demands. Furthermore, our approach has the potential to contribute to the development of more reliable and resilient power grids capable of withstanding the growing threat that GMDs pose.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Optimizing cloud motion estimation on the edge with phase correlation and optical flow

Abstract. Phase correlation (PC) is a well-known method for estimating cloud motion vectors (CMVs) from infrared and visible spectrum images. Commonly, phase shift is computed in the small blocks of the images using the fast Fourier transform. In this study, we investigate the performance and the stability of the blockwise PC method by changing the block size, the frame interval, and combinations of red, green, and blue (RGB) channels from the total sky imager (TSI) at the United States Atmospheric Radiation Measurement user facility's Southern Great Plains site. We find that shorter frame intervals, followed by larger block sizes, are responsible for stable estimates of the CMV, as suggested by the higher autocorrelations. The choice of RGB channels has a limited effect on the quality of CMVs, and the red and the grayscale images are marginally more reliable than the other combinations during rapidly evolving low-level clouds. The stability of CMVs was tested at different image resolutions with an implementation of the optimized algorithm on the Sage cyberinfrastructure test bed. We find that doubling the frame rate outperforms quadrupling the image resolution in achieving CMV stability. The correlations of CMVs with the wind data are significant in the range of 0.38–0.59 with a 95 % confidence interval, despite the uncertainties and limitations of both datasets. A comparison of the PC method with constructed data and the optical flow method suggests that the post-processing of the vector field has a significant effect on the quality of the CMV. The raindrop-contaminated images can be identified by the rotation of the TSI mirror in the motion field. The results of this study are critical to optimizing algorithms for edge-computing sensor systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Pulling thin single crystal silicon wafers from a melt: The new leading-edge solar substrate

The Floating Silicon Method (FSM) has been established as a viable, stable method for growing single crystal ribbons directly from a silicon melt. With intense helium jet cooling to drive the linear progress of a [111] facet, pulled in the <110> direction, ribbons in the 0.6 – 3.0 mm thickness range can be grown at linear growth rates from 0.3 mm/s to >6 mm/s as reported in the literature. We report on recent progress towards growing (100) oriented ribbons with a net thickness of less than 200 microns and a ribbon width up to 18 cm using a stable, continuous process in the Leading Edge prototype furnace. Here, the 3D details of the single crystal growth are explained by extending the mechanics of the Limit Cycle Theory, with novel Internal Side Effect morphology described by a proposed Facet Flow Theory. Grown-in crystalline defect distributions are described as well as values of critical impurities like oxygen, carbon, dopants, and metals that are relevant for use as wafers for solar cells.

14 SOLAR ENERGY↗
DOE OSTI2025

Fabrication and characterization of boron-terminated tetravacancies in monolayer hBN using STEM, EELS and electron ptychography

Tetravacancies in monolayer hexagonal boron nitride (hBN) with consistent edge termination (boron or nitrogen) form triangular nanopores with electrostatic potentials that can be leveraged for applications such as selective ion transport and neuromorphic computing. In order to quantitatively predict the properties of these structures, an atomic-level understanding of their local electronic and chemical environments is required. Moreover, robust methods for their precision manufacture are needed. Here we use electron irradiation in a scanning transmission electron microscope (STEM) at a high dose rate to drive the formation of boron-terminated tetravacancies in monolayer hBN. Characterization of the defects is achieved using aberration-corrected STEM, monochromated electron energy-loss spectroscopy (EELS), and electron ptychography. Z-contrast in STEM and chemical fingerprinting by core-loss EELS enable identification of the edge terminations, while electron ptychography gives insight into structural relaxation of the tetravacancies and provides evidence of enhanced electron density around the defect perimeters indicative of bonding effects.

2D hBN↗
DOE OSTI2021

Thermal quench in ITER locked mode disruptions

Simulations and theory are presented of an ITER locked mode thermal quench(TQ). In present experiments, locked mode disruptions have a long precursor phase,followed by a rapid termination and thermal quench, which can be identified with a resistive wall tearing mode (RWTM). In ITER, the RWTM will be slowed by the highly conductive vacuum vessel. Here, the rapid termination might be absent, and the plasma could remain in the precursor phase. If the edge temperature is in the collisional regime, the TQ would proceed on a long timescale, limited by the RWTM to almost100ms. This is an important self mitigating effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Investigations of radiation damaged arranged scintillating particle composites

Composites consisting of periodically arranged scintillating particles within an acrylic (poly methyl methacrylate) matrix are excellent templates for the gamma-insensitive detection of fast neutrons. In order to study the impact of radiation damage on acrylic-based scintillating particle composites, we exposed acrylic samples to different doses of fast neutrons up to 119.7 kGy and investigated the change in optical characteristics using optical absorption spectroscopy in the 200–800 nm wavelength range. The experimental results are compared to coupled MCNP6 and optical ray-tracing (FRED) simulations, and qualitative agreement is found. Neutron-induced radiation damage in the acrylic matrix manifests as a red-shift in the ultra-violet absorption edge of the acrylic and a corresponding decrease in the light propagation efficiency, which leads to detector performance degradation. We show for the first time that the effects of radiation damage can be mitigated by the use of wavelength-shifting coatings. We predict that composites with wavelength-shifting coatings can enable neutron detectors with high tolerance against neutron-induced radiation damage, a property that is particularly desired for neutron detection in high radiation environments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Direct Evidence of Multispecies Hydrodynamics in Ignition-Scale Hohlraums

A targeted experiment at the National Ignition Facility (NIF) confirms the presence of multispecies hydrodynamics in inertial confinement fusion hohlraums relevant to ignition. Here, the effects are identified by filling the gold hohlraum with a deuterium-tritium (DT) gas mixture instead of helium. As the hohlraum is heated by the NIF lasers, it implodes inward, compressing and heating the DT, which leads to fusion. The resulting DT-fusion neutrons are measured in space, time, yield, angle, and energy. A distinct, peaked, triangular shape in the radial neutron emission profile provides evidence of a “leaky piston” effect caused by the interpenetration of DT into the expanding gold. This process reduces the reversibility of DT compression and decreases neutron generation on the temporal trailing edge compared to the leading edge. These results are well described by multispecies hydrodynamics simulations, which reproduce the observed spatial and temporal features, as well as the total neutron yield and angularly resolved energy spectra. In contrast, conventional simulations that use only single-species physics fail to match the experimental data. They overpredict the neutron yield, produce a flatter-than-observed spatial profile, and show excessive emission at late times compared to the experimental measurements.

Higginson, Drew P. [Lawrence Livermore National La↗
DOE OSTI2022

Towards Compact Neural Networks via End-to-End Training: A Bayesian Tensor Approach with Automatic Rank Determination

Post-training model compression can reduce the inference costs of deep neural networks, but uncompressed training still consumes enormous hardware resources and energy. To enable low-energy training on edge devices, it is highly desirable to directly train a compact neural network from scratch with a low memory cost. Low-rank tensor decomposition is an effective approach to reduce the memory and computing costs of large neural networks. However, directly training low-rank tensorized neural networks is a very challenging task because it is hard to determine a proper tensor rank a priori, and the tensor rank controls both model complexity and accuracy. Here, this paper presents a novel end-to-end framework for low-rank tensorized training. We first develop a Bayesian model that supports various low-rank tensor formats (e.g., CANDECOMP/PARAFAC, Tucker, tensor-train, and tensor-train matrix) and reduces neural network parameters with automatic rank determination during training. Then we develop a customized Bayesian solver to train large-scale tensorized neural networks. Our training methods shows orders-of-magnitude parameter reduction and little accuracy loss (or even better accuracy) in the experiments. On a very large deep learning recommendation system with over 4.2 ×10 9 model parameters, our method can reduce the parameter number to 1.6 ×10 5 automatically in the training process (i.e., by 2.6 ×10 4 times) while achieving almost the same accuracy. Code is available at https://github.com/colehawkins/bayesian-tensor-rank-determination.

compact neural networks↗
DOE OSTI2023

Estimating Field-Level Perennial Bioenergy Grass Biomass Yields Using the Normalized Difference Red-Edge Index and Linear Regression Analysis for Central Virginia, USA

We investigated the indicative power of the normalized difference red-edge index (NDRE) for estimating field-level perennial bioenergy grass biomass yields utilizing Sentinel-2 imagery and a linear regression model as a rapid, cost-effective method for biomass yield estimations for bioenergy. We used 2019 data from three study sites containing mature perennial bioenergy grass stands in central Virginia, USA. Of the simulated daily NDRE values based on the temporally weighted averaging of two temporal neighbors, we found the strongest index–yield correlation on 11 August (R = 0.85). We estimated the perennial bioenergy grass biomass yields for (1) all sites using the data pooled from the three sites (all-site estimation) and (2) each site using the data pooled from the other two sites (cross-site estimation). The estimated field-level perennial bioenergy grass biomass yields strongly correlated with the recorded yields (average R2 = 0.76), with a root mean square error (RMSE) of 1.5 Mg/ha and a mean absolute error (MAE) of 1.2 Mg/ha for the all-site estimation. For the cross-site estimation, the site with diverse perennial grass types had the weakest correlation (R2 = 0.44) of the sites, indicating a difficulty in accounting for heterogeneous index–yield relationships in a single model. In addition to identifying a strong indicative power of the NDRE for estimating the overall perennial bioenergy grass biomass yields at a field level, the findings from this study call for an analysis across multiple perennial grasses and a comparison using multiple sites to understand (1) if the indicative power of the index shifts from the biomass of the specific perennial bioenergy grass type to the overall biomass during the growing season and (2) the level of perennial bioenergy grass heterogeneity that may hinder the remotely sensed biomass yield estimation using a single model.

09 BIOMASS FUELS↗
DOE OSTI2023

Investigation of local distortion effects on X-ray absorption of ferroelectric perovskites from first principles simulations

Understanding the role of ferroelectric polarization in modulating the electronic and structural properties of crystals is critical for advancing these materials for overcoming various technological and scientific challenges. However, due to difficulties in performing experimental methods with the required resolution, or in interpreting the results of methods therein, the nanoscale morphology and response of these surfaces to external electric fields has not been properly elaborated. Here, in this work, we investigate the effect of ferroelectric polarization and local distortions in a BaTiO 3 perovskite, using two widely used computational approaches which treat the many-body nature of X-ray excitations using different philosophies, namely the many-body, delta-self-consistent-field determinant (mb-ΔSCF) and the Bethe–Salpeter equation (BSE) approaches. We show that in agreement with our experiments, both approaches consistently predict higher excitations of the main peak in the O–K edge for the surface with upward polarization. However, the mb-ΔSCF approach mostly fails to capture the L 2,3 separations at the Ti–L edge, due to the absence of spin–orbit coupling in Kohn–Sham density functional theory (KS-DFT) at the generalized gradient approximation level. On the other hand, and most promising, we show that application of the GW/BSE approach successfully reproduces the experimental XAS, both the relative peak intensities as well as the L 2,3 separations at the Ti–L edges upon ferroelectric switching. Thus simulated XAS is shown to be a powerful method for capturing the nanoscale structure of complex materials, and we underscore the need for many-body perturbation approaches, with explicit consideration of core-hole and multiplet effects, for capturing the essential physics in these systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Tuning the interlayer coupling in La 0.7 Sr 0.3 Mn 0.95 Ru 0.05 O 3 / LaNiO 3 multilayers with perpendicular magnetic anisotropy

In ferromagnetic oxide epitaxial multilayers, magnetic properties and interlayer coupling are determined by a variety of factors. Beyond the contribution of interlayer exchange coupling, strain and interfacial effects, such as structural reconstructions or charge transfer, play significant roles, resulting in complex magnetic behavior. In this study, the interlayer coupling of ferromagnetic La 0.7 ⁢Sr 0.3⁢ Mn 0.95 ⁢Ru 0.05 ⁢O 3 (LSMRO) layers (8 nm thick) was investigated, when separated by epitaxial spacers of paramagnetic metallic LaNiO 3 (LNO), in stacks exhibiting perpendicular magnetic anisotropy. By varying the thickness of the spacer, it was found that the coupling between two LSMRO layers changes from antiferromagnetic (with a 4-unit-cell-thick LNO spacer) to ferromagnetic (with a 6-unit-cell-thick LNO spacer). For multilayers comprising five LSMRO layers and a 4-unit-cell-thick LNO spacer, the antiferromagnetic coupling was preserved. However, the effective magnetic anisotropy changed, causing the magnetization to cant more towards the in-plane direction. Further, this behavior was corroborated by X-ray magnetic circular dichroism (XMCD) investigations at the Mn and Ni L 3 edges. The XMCD results indicated that the 4-unit-cell-thick LNO spacer in the multilayer become magnetically ordered, closely following the magnetization of adjacent LSMRO layers.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Toroidal modeling of runaway electron loss due to 3D fields in ITER

Mitigation of runaway electrons (REs) by three-dimensional (3D) magnetic field perturbations is numerically investigated for the ITER 15 MA baseline D–T scenario, utilizing the MARS-F code (Liu et al Phys. Plasmas 7 3681) with a drift orbit test particle tracing module. Considered are two types of 3D fields: the n = 3 (n is the toroidal mode number) resonant magnetic perturbation (RMP) utilized for the purpose of controlling the edge localized modes in ITER, and perturbations generated by the n = 1 magneto-hydrodynamic (MHD) instabilities in a post-disruption plasma. The RMP field, applied to a pre-disruption plasma, is found to be moderately effective in mitigating the RE seeds in ITER when vacuum field model is assumed. Up to ~40% loss fraction is possible at 90 kA-turn coil current. The mitigation efficiency is however substantially reduced, down to less than 5%, when the plasma response is taken into account. This is due to strong screening of the resonant magnetic field components by the plasma response resulting in much less field line stochasticity. On the other hand, the MARS-F modeling, based on the DINA-simulated post-disruption equilibria, shows that the n = 1 resistive kink instabilities develop in these plasmas, as the edge safety factor qa evolves and drops below integer numbers. RE mitigation by these MHD instabilities is sensitive to the eigenmode structure. The best mitigation is achieved as q a drops below 3, when a global kink instability occurs that encompasses both internal and external components. This global instability is found to be capable of mitigating over 80% MeV-level passing RE orbits at a field perturbation |δB|/B 0 that is comparable to that observed in DIII-D experiments, and full mitigation if the perturbation amplitude is doubled. The 'wetted' area on the ITER limiting surface, due to MHD instability induced RE loss, generally increases with the perturbation amplitude (together with increasing loss fraction). At the highest perturbation level assumed in this study, the wetted area reaches ~60% of the total limiting surface area. Lastly, the lost RE orbits mainly strike the outer divertor region of the limiting surface, with some fraction also hitting a wide area along the inboard side of the surface.

3D fields↗
DOE OSTI2020

Drop size distribution broadening mechanisms in a bin microphysics Eulerian model

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

54 ENVIRONMENTAL SCIENCES↗