Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PARTICLE MODELS”

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 1,765 records · Page 98

Triaxial-band structures, chirality, and magnetic rotation in La 133

The structure of 133 La has been investigated using the 116 Cd( 22 Ne,4pn) reaction and the Gammasphere array. Three new bands of quadrupole transitions and one band of dipole transitions are identified and the previously reported level scheme is revised and extended to higher spins. The observed structures are discussed using the cranked Nilsson-Strutinsky formalism, covariant density functional theory, and the particle-rotor model. Triaxial configurations are assigned to all observed bands. For the high-spin bands it is found that rotations around different axes can occur, depending on the configuration. The orientation of the angular momenta of the core and of the active particles is investigated, suggesting chiral rotation for two nearly degenerate dipole bands and magnetic rotation for one dipole band. As a result, it is shown that the h 11/2 neutron holes present in the configuration of the nearly degenerate dipole bands have significant angular momentum components not only along the long axis but also along the short axis, contributing to the balance of the angular momentum components along the short and long axes and thus giving rise to a chiral geometry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Multiscale and Multiphysics Modeling of Additive Manufacturing of Advanced Materials

The objective of this proposed project is to research and develop a prediction tool for advanced additive manufacturing (AAM) processes for advanced materials and develop experimental methods to provide fundamental properties and establish validation data. Aircraft structures and engines demand materials that are stronger, useable at much higher temperatures, provide less acoustic transmission, and enable more aeroelastic tailoring than those currently used. Significant improvements in properties can only be achieved by processing the materials under nonequilibrium conditions, such as AAM processes. AAM processes encompass a class of processes that use a focused heat source to create a melt pool on a substrate. Examples include Electron Beam Freeform Fabrication and Direct Metal Deposition. These types of additive processes enable fabrication of parts directly from CAD drawings. To achieve the desired material properties and geometries of the final structure, assessing the impact of process parameters and predicting optimized conditions with numerical modeling as an effective prediction tool is necessary. The targets for the processing are multiple and at different spatial scales, and the physical phenomena associated occur in multiphysics and multiscale. In this project, the research work has been developed to model AAM processes in a multiscale and multiphysics approach. A macroscale model was developed to investigate the residual stresses and distortion in AAM processes. A sequentially coupled, thermomechanical, finite element model was developed and validated experimentally. The results showed the temperature distribution, residual stress, and deformation within the formed deposits and substrates. A mesoscale model was developed to include heat transfer, phase change with mushy zone, incompressible free surface flow, solute redistribution, and surface tension. Because of excessive computing time needed, a parallel computing approach was also tested. In addition, after investigating various methods, a Smoothed Particle Hydrodynamics Model (SPH Model) was developed to model wire feeding process. Its computational efficiency and simple architecture makes it more robust and flexible than other models. More research on material properties may be needed to realistically model the AAM processes. A microscale model was developed to investigate heterogeneous nucleation, dendritic grain growth, epitaxial growth of columnar grains, columnar-to-equiaxed transition, grain transport in melt, and other properties. The orientations of the columnar grains were almost perpendicular to the laser motion's direction. Compared to the similar studies in the literature, the multiple grain morphology modeling result is in the same order of magnitude as optical morphologies in the experiment. Experimental work was conducted to validate different models. An infrared camera was incorporated as a process monitoring and validating tool to identify the solidus and mushy zones during deposition. The images were successfully processed to identify these regions. This research project has investigated multiscale and multiphysics of the complex AAM processes thus leading to advanced understanding of these processes. The project has also developed several modeling tools and experimental validation tools that will be very critical in the future of AAM process qualification and certification.

Liou, Frank↗

Developing New Tools for Modeling Rocket Plume-Surface Interactions

With NASA’s goal to land the next human on the lunar surface in the next few years, it has become vitally important that we have a better understanding of how future landing spacecraft will interact with the unique properties of regolith¬¬––the layer of loose, unconsolidated dust and rock on the lunar surface¬¬––which can cause hazards like visual obstructions, particulate clouds, and cratering of the landing zone. Researchers from the Fluid Dynamics Branch at NASA’s Marshall Space Flight Center are performing plume-surface interaction (PSI) simulations between lander engine plumes and unprepared regolith surfaces, and have developed new tools to provide predictive PSI environments for various NASA projects and missions, including the Human Lander System (HLS), Commercial Lunar Payload Services (CLPS), and future Mars landers. These tools allow the researchers to determine how to best meet the simulation and time requirements for each project by varying model fidelity. The highest fidelity tool is the Gas Granular Flow Solver (Loci/GGFS) that models gas-particle multi-phase interactions to predict regolith cratering and ejection of particles into the immediate surroundings of the lander. At its highest fidelity, it can model microscopic regolith particle interactions with a particle size/shape distribution that statistically replicates actual regolith, however, to be most effective with today’s computing resources, it is currently run using only one to three equivalent particle sizes/shapes. The team also incorporated engineering models into their software suite to create production-ready hybrid tools with reduced fidelity. At the lowest fidelity, the computational fluid dynamics (CFD) code Loci/CHEM+DIGGEM can predict crater depth over time by relating local CFD-predicted surface shear stresses to a model of erosion mass flux.

plume surface interaction↗

Mixed-Material Scintillators (MMSS Quarterly Report FY20Q2)

This project aims to invent, model, and prototype mixed-material scintillator systems (MMSSs), a new class of radiation detectors that use heterogenous internal structures of different scintillating materials to detect additional properties of radiation. These internal structures can be produced using additive manufacture (3D printing) of scintillator, a currently emerging application of additive manufacture technology. MMSSs combine the low cost and complexity of conventional scintillation detectors with capabilities currently only available in more expensive and complex detectors. By identifying promising MMSS designs, this project will enable a new class of detectors to meet DNN’s mission needs for SNM detection. This quarter, we delivered one the MMSS Particle-ID Modeling report, a major deliverable of this project. This report describes our results on simulations of the particle-ID (PID) class of detectors. We found that these detectors outperform competing approaches for neutron-gamma discrimination, neutron source pointing, and neutron spectroscopy. Further highlights of the report are described below. During this quarter, We alerted HQ to a change in schedule for the remaining simulations planned for this project: we plan to put off further simulation and reporting of the position-resolving (PR) detector class until July. This change in schedule has allowed us to go into additional detail in the PID work that is both very promising and motivates our upcoming proposal. The other active task in this project is prototyping 3D-printed scintillators with the aim of realizing the designs we’ve simulated. On that front, we had three focuses: measuring the light output, increasing the size of printed parts, and combining blue and green scintillators. We succeeded in measuring the light output, showing an output of 30% of commercial standard. This result leaves room for improvement but is within striking distance of where we need to be. We were on track to demonstrate a key requirement for increasing the size of parts and also to show combined green/blue prints before the laboratory moved to a Minimum Safe (MinSafe) Operations posture as a result of local Shelter In Place orders due to Covid-19.

42 ENGINEERING↗

Scattering Properties of Large Irregular Cosmic Dust Particles at Visible Wavelengths

The effect of internal inhomogeneities and surface roughness on the scattering behavior of large cosmic dust particles is studied by comparing model simulations with laboratory measurements. The present work shows the results of an attempt to model a dust sample measured in the laboratory with simulations performed by a ray-optics model code. We consider this dust sample as a good analogue for interplanetary and interstellar dust as it shares its refractive index with known materials in these media. Several sensitivity tests have been performed for both structural cases (internal inclusions and surface roughness). Three different samples have been selected to mimic inclusion/coating inhomogeneities: two measured scattering matrices of hematite and white clay, and a simulated matrix for water ice. These three matrices are selected to cover a wide range of imaginary refractive indices. The selection of these materials also seeks to study astrophysical environments of interest such as Mars, where hematite and clays have been detected, and comets. Based on the results of the sensitivity tests shown in this work, we perform calculations for a size distribution of a silicate-type host particle model with inclusions and surface roughness to reproduce the experimental measurements of a dust sample. The model fits the measurements quite well, proving that surface roughness and internal structure play a role in the scattering pattern of irregular cosmic dust particles.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigating the Effects of Anisotropic Mass Transport on Dendrite Growth in High Energy Density Lithium Batteries

Dendrite formation on the electrode surface of high energy density lithium (Li) batteries causes safety problems and limits their applications. Suppressing dendrite growth could significantly improve Li battery performance. Dendrite growth and morphology is a function of the mixing in the electrolyte near the anode interface. Most research into dendrites in batteries focuses on dendrite formation in isotropic electrolytes (i.e., electrolytes with isotropic diffusion coefficient). In this work, an anisotropic diffusion reaction model is developed to study the anisotropic mixing effect on dendrite growth in Li batteries. The model uses a Lagrangian particle-based method to model dendrite growth in an anisotropic electrolyte solution. The model is verified by comparing the numerical simulation results with analytical solutions, and its accuracy is shown to be better than previous particle-based anisotropic diffusion models. Several parametric studies of dendrite growth in an anisotropic electrolyte are performed and the results demonstrate the effects of anisotropic transport on dendrite growth and morphology, and show the possible advantages of anisotropic electrolytes for dendrite suppression.

SPH, high energy density lithium batteries, precip↗

GRAVITATIONAL MODEL OF HIGH-ENERGY PARTICLES IN A COLLIMATED JET

Observations suggest that relativistic particles play a fundamental role in the dynamics of jets emerging from active galactic nuclei as well as in their interaction with the intracluster medium. However, no general consensus exists concerning the acceleration mechanism of those high-energy particles. A gravitational acceleration mechanism is proposed here in which particles leaving precise regions within the ergosphere of a rotating supermassive black hole (BH) produce a highly collimated flow. These particles follow unbound geodesics which are asymptotically parallel to the spin axis of the BH and are characterized by the energy E, the Carter constant Q, and zero angular momentum of the component L{sub z} . If environmental effects are neglected, the present model predicts the presence of electrons with energies around 9.4 GeV at distances of about 140 kpc from the ergosphere. The present mechanism can also accelerate protons up to the highest energies observed in cosmic rays by the present experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

Complementarity of Galactic radio and collider data in constraining WIMP dark matter models

In this work we confront dark matter models to constraints that may be derived from radio synchrotron radiation from the Galaxy, taking into account the astrophysical uncertainties and we compare these to bounds set by accelerator and complementary indirect dark matter searches. Specifically we apply our analysis to three popular particle physics models. First, a generic effective operator approach, in which case we set bounds on the corresponding mass scale, and then, two specific UV completions, the Z' and Higgs portals. We show that for many candidates, the radio synchrotron limits are competitive with the other searches, and could even give the strongest constraints (as of today) with some reasonable assumptions regarding the astrophysical uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗

Derivation of an Analytical Solution to a Reaction-Diffusion Model for Autocatalytic Degradation and Erosion in Polymer Microspheres

A mathematical reaction-diffusion model is defined to describe the gradual decomposition of polymer microspheres composed of poly(D,L-lactic-co-glycolic acid) (PLGA) that are used for pharmaceutical drug delivery over extended periods of time. The partial differential equation (PDE) model treats simultaneous first-order generation due to chemical reaction and diffusion of reaction products in spherical geometry to capture the microsphere-size dependent effects of autocatalysis on PLGA erosion that occurs when the microspheres are exposed to aqueous media such as biological fluids. The model is solved analytically for the concentration of the autocatalytic carboxylic acid end groups of the polymer chains that comprise the microspheres as a function of radial position and time. The analytical solution for the reaction and transport of the autocatalytic chemical species is useful for predicting the conditions under which drug release from PLGA microspheres transitions from diffusion-controlled to erosion-controlled release, for understanding the dynamic coupling between the PLGA degradation and erosion mechanisms, and for designing drug release particles. The model is the first to provide an analytical prediction for the dynamics and spatial heterogeneities of PLGA degradation and erosion within a spherical particle. The analytical solution is applicable to other spherical systems with simultaneous diffusive transport and first-order generation by reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Use of Satellite-Measured Aerosol Optical Depth to Constrain Biomass Burning Emissions Source Strength in a Global Model GOCART

Small particles in the atmosphere, called "atmospheric aerosol" have a direct effect on Earth climate through scattering and absorbing sunlight, and also an indirect effect by changing the properties of clouds, as they interact with solar radiation as well. Aerosol typically stays in the atmosphere for several days, and can be transported long distances, affecting air quality, visibility, and human health not only near the source, but also far downwind. Smoke from vegetation fires is one of the main sources of atmospheric aerosol; other sources include anthropogenic pollution, dust, and sea salt. Chemistry transport models (CTMs) are among the major tools for studying the atmospheric and climate effects of aerosol. Due to the considerable variation of aerosol concentrations and properties on many temporal and spatial scales, and the complexity of the processes involved, the uncertainties in aerosol effects on climate are large, as is featured in the latest report of Intergovernmental Panel on Climate Change (IPCC) in 2007. Reducing this uncertainty in the models is very important both for predicting future climate scenarios and for regional air quality forecasting and mitigation. During vegetation fires, also called biomass burning (BB) events, complex mixture of gases and particles is emitted. The amount of BB emissions is usually estimated taking into account the intensity and size of the fire and the properties of burning vegetation. These estimates are input into CTMs to simulate BB aerosol. Unfortunately, due to large variability of fire and vegetation properties, the quantity of BB emissions is very difficult to estimate and BB emission inventories provide numbers that can differ by up to the order of magnitude in some regions. Larger uncertainties in data input make uncertainties in model output larger as well. A powerful way to narrow the range of possible model estimates is to compare model output to observations. We use satellite observations of aerosol properties, specifically aerosol optical depth, which is directly proportional to the amount of aerosol in the atmosphere, and compare it to the model output. Assuming the model represents aerosol transport and particle properties correctly, the amount of BB emissions determines the simulated aerosol optical depth. In this study, we explore the regional performance of 13 commonly used emission estimates. These are each input to global Goddard Chemistry Aerosol Radiation and Transport (GOCART) model. We then evaluate how well each emission estimate reproduces the smoke aerosol optical depth measured by the MODIS instrument. We compared GOCART-simulate aerosol optical depth with that measured from the satellite for 124 fire cases around the world during 2006 and 2007. We summarize the regional performance of each emission inventory and discuss reasons for their differences by considering the assumptions made during their development. We also show that because stronger wind disperses smoke plumes more readily, in cases with stronger wind, a larger increase in emission amount is needed to increase aerosol optical depth. In quiet, low-wind-speed environments, BB emissions produce a more significant increase in aerosol optical depth, other things being equal. Using the region-specific, quantitative relationships derived in our paper, together with the wind speed obtained from another source for a given fire case, we can constrain the amount of emission required in the model to reproduce the observations. The results of this paper are useful to the developers of BB emission inventories, as they show the strengths and weaknesses of individual emission inventories in different regions of the globe, and also for modelers who use these inventories and wish to improve their model results.

Petrenko, Mariya↗

Search for $ZW/ZZ \to \ell^+ \ell^-$ + Jets Production in $p\bar{p}$ Collisions at CDF

The Standard Model of particle physics describes weak interactions mediated by massive gauge bosons that interact with each other in well-defined ways. Observations of the production and decay of WW, WZ, and ZZ boson pairs are an opportunity to check that these self-interactions agree with the Standard Model predictions. Furthermore, final states that include quarks are very similar to the most prominent final state of Higgs bosons produced in association with a W or Z boson. Diboson production where WW is a significant component has been observed at the Tevatron collider in semi-hadronic decay modes. We present a search for ZW and ZZ production in a final state containing two charged leptons and two jets using 8.9 fb -1 of data recorded with the CDF detector at the Tevatron. We select events by identifying those that contain two charged leptons, two hadronic jets, and low transverse missing energy (E T ). We increase our acceptance by using a wide suite of high-p T lepton triggers and by relaxing many lepton identification requirements. We develop a new method for calculating corrections to jet energies based on whether the originating parton was a quark or gluon to improve the agreement between data and the Monte Carlo simulations used to model our diboson signal and dominant backgrounds. We also make use of neural-network-based discriminants that are trained to pick out jets originating from b quarks and light-flavor quarks, thereby increasing our sensitivity to Z → b$\bar{b}$ and W=Z → q$\bar{p'}$0 decays, respectively. The number of signal events is extracted through a simultaneous fit to the dijet mass spectrum in three channels: a heavy-flavor tagged channel, a light-flavor tagged channel, and an untagged channel. We measure σ ZW/ZZ = 2.5 +2.0 -1.0 pb, which is consistent with the SM cross section of 5.1 pb. We establish an upper limit on the cross section of σ ZW/ZZ < 6.1 pb at 95% CL.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Glass shell manufacturing in space

Highly-uniform, hollow glass spheres, which are used for inertial-confinement fusion targets, are formed from metal-organic gel powder feedstock in a drop-tower furnace. The modelling of this gel-to-sphere transformation has consisted of three phases: gel thermochemistry, furnance-to-gel heat transfer, and gravity-driven degradation of the concentricity of the molten shell. The heat transfer from the furnace to the free-falling gel particle was modelled with forced convection. The gel mass, dimensions, and specific heat as well as furnace temperature profile and furnace gas conductivity, were controlled variables. This model has been experimentally verified. In the third phase, a mathematical model was developed to describe the gravity-driven degradation of concentricity in molten glass shells.

Nolen, R. L.↗

A direct comparison of particle-resolved and point-particle methods in decaying turbulence

We use particle-resolved direct numerical simulation (PR-DNS) as a model-free physics-based numerical approach to validate particle acceleration modelling in gas-solid suspensions. To isolate the effect of the particle acceleration model, we focus on point-particle direct numerical simulation (PP-DNS) of a collision-free dilute suspension with solid-phase volume fraction $\unicode[STIX]{x1D719}=0.001$ in a decaying isotropic turbulent particle-laden flow. The particle diameter $d_{p}$ in the suspension is chosen to be the same as the initial Kolmogorov length scale $\unicode[STIX]{x1D702}_{0}$ ( $d_{p}/\unicode[STIX]{x1D702}_{0}=1$ ) in order to overlap with the regime where PP-DNS is valid. We assess the point-particle acceleration model for two different particle Stokes numbers, $St_{\unicode[STIX]{x1D702}}=1$ and 100. For the high Stokes number case, the Stokes drag model for particle acceleration under-predicts the true particle acceleration. In addition, second moment quantities which play key roles in the physical evolution of the gas–solid suspension are not correctly captured. Considering finite Reynolds number corrections to the acceleration model improves the prediction of the particle acceleration probability density function and second moment statistics of the point-particle model compared with the particle-resolved simulation. We also find that accounting for the undisturbed fluid velocity in the acceleration model can be of greater importance than using the most appropriate acceleration model for a given physical problem.

Mechanics↗

Detection of a particle shower at the Glashow resonance with IceCube

The Glashow resonance describes the resonant formation of a W - boson during the interaction of a high-energy electron antineutrino with an electron, peaking at an antineutrino energy of 6.3 petaelectronvolts (PeV) in the rest frame of the electron. Whereas this energy scale is out of reach for currently operating and future planned particle accelerators, natural astrophysical phenomena are expected to produce antineutrinos with energies beyond the PeV scale. Here we report the detection by the IceCube neutrino observatory of a cascade of high-energy particles (a particle shower) consistent with being created at the Glashow resonance. A shower with an energy of 6.05 ± 0.72 PeV (determined from Cherenkov radiation in the Antarctic Ice Sheet) was measured. Features consistent with the production of secondary muons in the particle shower indicate the hadronic decay of a resonant W - boson, confirm that the source is astrophysical and provide improved directional localization. The evidence of the Glashow resonance suggests the presence of electron antineutrinos in the astrophysical flux, while also providing further validation of the standard model of particle physics. Its unique signature indicates a method of distinguishing neutrinos from antineutrinos, thus providing a way to identify astronomical accelerators that produce neutrinos via hadronuclear or photohadronic interactions, with or without strong magnetic fields. As such, knowledge of both the flavour (that is, electron, muon or tau neutrinos) and charge (neutrino or antineutrino) will facilitate the advancement of neutrino astronomy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Constraints on secret neutrino interactions after Planck

Neutrino interactions beyond the standard model of particle physics may affect the cosmological evolution and can be constrained through observations. We consider the possibility that neutrinos possess secret scalar or pseudoscalar interactions mediated by the Nambu-Goldstone boson of a still unknown spontaneously broken global U(1) symmetry, as in, e.g., Majoron models. In such scenarios, neutrinos still decouple at T≅ 1 MeV, but become tightly coupled again (''recouple'') at later stages of the cosmological evolution. We use available observations of the cosmic microwave background (CMB) anisotropies, including Planck 2013 and the joint BICEP2/Planck 2015 data, to derive constraints on the quantity γ{sub νν}{sup 4}, parameterizing the neutrino collision rate due to scalar or pseudoscalar interactions. We consider both a minimal extension of the standard ΛCDM model, and more complicated scenarios with extra relativistic degrees of freedom or non-vanishing tensor amplitude. For a wide range of dataset and model combinations, we find a typical constraint γ{sub νν}{sup 4} ∼< 0.9× 10{sup −27} (95% C.L.), implying an upper limit on the redshift z{sub νrec} of neutrino recoupling 0∼< 850, leaving open the possibility that the latter occured well before hydrogen recombination. In the framework of Majoron models, the upper limit on γ{sub νν} roughly translates on a constraint g ∼< 8.2× 10{sup −7} on the Majoron-neutrino coupling constant g. In general, the data show a weak (∼ 1σ) but intriguing preference for non-zero values of γ{sub νν}{sup 4}, with best fits in the range γ{sub νν}{sup 4} = (0.15–0.35)× 10{sup −27}, depending on the particular dataset. This is more evident when either high-resolution CMB observations from the ACT and SPT experiments are included, or the possibility of non-vanishing tensor modes is considered. In particular, for the minimal model ΛCDM+γ{sub νν} and including the Planck 2013, ACT and SPT data, we report γ{sub νν}{sup 4}=(0.44{sup +0.17}{sub −0.36})×10{sup −27} (0300 ∼< z{sub νrec} ∼< 550) at 68% confidence level.

79 ASTRONOMY AND ASTROPHYSICS↗

THE POWER-LAW SPECTRA OF ENERGETIC PARTICLES DURING MULTI-ISLAND MAGNETIC RECONNECTION

Power-law distributions are a near-universal feature of energetic particle spectra in the heliosphere. Anomalous cosmic rays (ACRs), super-Alfvenic ions in the solar wind, and the hardest energetic electron spectra in flares all have energy fluxes with power laws that depend on energy E approximately as E {sup -1.5}. We present a new model of particle acceleration in systems with a bath of merging magnetic islands that self-consistently describes the development of velocity-space anisotropy parallel and perpendicular to the local magnetic field and includes the self-consistent feedback of pressure anisotropy on the merging dynamics. By including pitch-angle scattering we obtain an equation for the omnidirectional particle distribution f (v, t) that is solved in closed form to reveal v {sup -5} (corresponding to an energy flux varying as E {sup -1.5}) as a near-universal solution as long as the characteristic acceleration time is short compared with the characteristic loss time. In such a state, the total energy in the energetic particles reaches parity with the remaining magnetic free energy. More generally, the resulting transport equation can serve as the basis for calculating the distribution of energetic particles resulting from reconnection in large-scale inhomogeneous systems.

79 ASTRONOMY AND ASTROPHYSICS↗

Coriolis coupling effects in proton-pickup spectroscopic factors from 12 B

Spectroscopic factors to low-lying negative-parity states in Be 11 extracted from the B 12 ( d , He 3 ) Be 11 proton-removal reaction are interpreted within the rotational model. Earlier predictions of the p -wave proton-removal strengths in the strong-coupling limit of the Nilsson model underestimated the spectroscopic factors to the 3 / 2 1 - and 5 / 2 1 - states and suggested that deviations in the 1 + ground state of the odd-odd B 12 due to Coriolis coupling should be further explored. Furthermore, in this work we use the particle rotor model to take into account these effects and obtain a good description of the level scheme in B 11 , with a moderate K mixing of the proton Nilsson levels [110]1/2 and [101]3/2. This mixing, present in the 1 + bandhead of B 12 , is key to explaining the proton-pickup data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

An evaluation of size-resolved cloud microphysics scheme numerics for use with radar observations Part II: Condensation and evaporation

Accurate numerical modeling of clouds and precipitation is essential for weather forecasting and climate change research. While size-resolved (bin) cloud microphysics models predict particle size distributions without imposing shapes, results are subject to artificial size distribution broadening owing to numerical diffusion associated with various processes. Whereas Part I of this study addressed collision–coalescence, here we investigate numerical diffusion that occurs in solving condensation and evaporation. In a parcel model framework, all of the numerical schemes examined converge to one solution of condensation and evaporation as the mass grid is refined, and the advection-based schemes are recommended over the reassigning schemes. Including Eulerian vertical advection in a column limits the convergence to some extent, but that limitation occurs at a sufficiently fine mass grid, and the number of iterations in solving vertical advection should be minimized to reduce numerical diffusion. Insubstantial numerical diffusion in solving condensation can be amplified if collision–coalescence is also active, which in turn can be substantially diminished if turbulence effects on collision are incorporated. Here, large-eddy simulations of a drizzling stratocumulus field reveal that changes in moments of Doppler spectra obtained using different mass grids are consistent with those obtained from the simplified framework, and that spectral moments obtained using a mass grid designed to effectively reduce numerical diffusion are generally closer to observations. Notable differences between the simulations and observations still exist, and our results suggest a need to consider whether factors other than numerical diffusion in the fundamental process schemes employed can cause such differences.

54 ENVIRONMENTAL SCIENCES↗