Search NASA⌕ Search

SEARCH · Search NASA

Results for “artifacts”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Effects of X-rays on Perovskite Solar Cells

Synchrotron micro- and nanoprobe beamlines have demonstrated great potential to advance photovoltaic devices. Most importantly, their small X-ray spot size has enabled the direct correlation of electrical performance with elemental composition at subgrain resolution for a variety of polycrystalline solar cells. Whereas the bulk of most inorganic semiconductors is stable under the high X-ray flux of focused X-ray beams, semiconductors with organic components are prone to a variety of degradation mechanisms. This is particularly critical to evaluate for the emerging organometal halide perovskite solar cells. In this study, we investigate the effects of hard X-rays on the nanoscale performance and elemental distribution of these solar cells. We show that their composition does not change during common operando and in situ measurements at synchrotron nanoprobes. However, we found a significant X-ray-induced electronic degradation of solar cells with methylammonium lead iodide absorbers. Time- and dose-dependent measurements unveiled two characteristic degradation time constants on the order of 12 and 200 s that are independent of the X-ray flux. On the basis of heat and dose simulations, we attribute the fast decay to the dose-driven creation of recombination centers, while the slow decay is compatible with the observation of compositional changes. Finally, we detail how degradation-induced measurement artifacts can be outrun and showcase the high correlation of the X-ray-beam-induced current with the iodine and lead distribution.

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

Machine Learning Automated Analysis of Enormous Synchrotron X-ray Diffraction Datasets

X-ray diffraction (XRD) data analysis can be a time-consuming and laborious task. Deep neural network (DNN) based models trained with synthetic XRD patterns have been proven to be a highly efficient, accurate, and automated method for analyzing common XRD data collected from solid samples in ambient environments. However, it remains unclear whether synthetic XRD-based models can be effective in solving micro(μ)-XRD mapping data for in situ experiments involving liquid phases, which always have lower quality and significant artifacts. In this study, we collected μ-XRD mapping data from a LaCl 3 -calcite hydrothermal fluid system and trained two categories of models to analyze the experimental XRD patterns. Here, the models trained solely with synthetic XRD patterns showed low accuracy (as low as 64%) when solving experimental μ-XRD mapping data. However, the accuracy of the DNN models significantly improved (90% or above) when we trained them with a data set containing both synthetic and a small number of labeled experimental μ-XRD patterns. This study highlights the importance of labeled experimental patterns in training DNN models to solve μ-XRD mapping data from in situ experiments involving liquid phases.

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

X-ray Photoelectron Spectroscopy Analysis of Nafion-Containing Samples: Pitfalls, Protocols, and Perceptions of Physicochemical Properties

X-ray photoelectron spectroscopy (XPS) is one of the most common techniques used to analyze the surface composition of catalysts and support materials used in polymer electrolyte membrane (PEM) fuel cells and electrolyzers, providing important insights for further improvement of their properties. Characterization of catalyst layers (CLs) is more challenging, which can be at least partially attributed to the instability of ionomer materials such as Nafion during measurements. This work explores the stability of Nafion during XPS measurements, illuminating and addressing Nafion degradation concerns. The extent of Nafion damage as a function of XPS instrumentation, measurement conditions, and sample properties was evaluated across multiple instruments. Results revealed that significant Nafion damage to the ion-conducting sulfonic acid species (>50% loss in sulfur signal) may occur in a relatively short time frame (tens of minutes) depending on the exact nature of the sample and XPS instrument. This motivated the development and validation of a multipoint XPS data acquisition protocol that minimizes Nafion damage, resulting in reliable data acquisition by avoiding significant artifacts from Nafion instability. The developed protocol was then used to analyze both thin film ionomer samples and Pt/C-based CLs. Comparison of PEM fuel cell CLs to Nafion thin films revealed several changes in Nafion spectral features attributed to charge transfer due to interaction with conductive catalyst and support species. This study provides a method to reliably characterize ionomer-containing samples, facilitating fundamental studies of the catalyst-ionomer interface and more applied investigations of structure-processing-performance correlations in PEM fuel cell and electrolyzer CLs.

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

Universal Size-Dependent Stokes Shifts in Lead Halide Perovskite Nanocrystals

Size-dependent photoluminescence Stokes shifts (Δ E s ) universally exist in CsPbX 3 (X = Cl – , Br – , or I – ) perovskite nanocrystals (NCs). Δ E s values, which range from ~15 to 100 meV for NCs with average edge lengths ( l ) from approximately 13 to 3 nm, are halide-dependent such that Δ E s (CsPbI 3 ) > Δ E s (CsPbBr 3 ) ≳ Δ E s (CsPbCl 3 ). Observed size-dependent Stokes shifts are not artifacts of ensemble size distributions as demonstrated through measurements of single CsPbBr 3 NC Stokes shifts ((Δ E s ) = 42 ± 5 meV), which are in near quantitative agreement with associated ensemble ( l = 6.8 ± 0.8 nm) Δ E s values (Δ E s ≈ 50 meV). Transient differential absorption measurements additionally illustrate no significant spectral dynamics on the picosecond time scale that would contribute to Δ E s . This excludes polaron formation as being responsible for Δ E s . Overall, the results point to an origin for Δ E s , intrinsic to the size-dependent electronic properties of individual perovskite NCs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Dynamic Effects and Hydrogen Bonding in Mixed-Halide Perovskite Solar Cell Absorbers

The organic component (methylammonium) of CH 3 NH 3 PbI 3- x Cl x -based perovskites shows electronic hybridization with the inorganic framework via H-bonding between N and I sites. Femtosecond dynamics induced by core excitation are shown to strongly influence the measured X-ray emission spectra and the resonant inelastic soft X-ray scattering of the organic components. The N K core excitation leads to a greatly increased N-H bond length that modifies and strengthens the interaction with the inorganic framework compared to that in the ground state. The study indicates that excited-state dynamics must be accounted for in spectroscopic studies of this perovskite solar cell material, and the organic-inorganic hybridization interaction suggests new avenues for probing the electronic structure of this class of materials. Furthermore, it is incidentally shown that beam damage to the methylamine component can be avoided by moving the sample under the soft X-ray beam to minimize exposure and that this procedure is necessary to prevent the creation of experimental artifacts.

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

4D Imaging of ZnO-Coated Nanoporous Al 2 O 3 Aerogels by Chemically Sensitive Ptychographic Tomography: Implications for Designer Catalysts

The 3D chemical structure (4D spectromicroscopy) of nanoporous Al 2 O 3 aerogels coated with ZnO by atomic layer deposition (ALD) was studied by multienergy scanning transmission X-ray microscopy. These materials are representative of a class of designer catalysts in which the nanoporous support is prepared separately from the active catalyst material, which is subsequently introduced by ALD, thereby allowing independent optimization of the morphology, chemistry, and spatial distribution of the support and catalyst. The samples studied were prepared by Ga ion and Xe plasma focused ion beam (FIB) milling as well as drop casting from water suspension. Zn L and Al K edge spectra of six samples with three different ZnO loadings were measured to investigate how loading and different sample preparation methods affect the 3D distribution of the ZnO and Al 2 O 3 . Scanning transmission X-ray microscopy (STXM) and ptychographic imaging at two energies each at the Zn L 3 and Al K edge were measured. The ptychography data were analyzed by using the SHARP reconstruction code to generate quantitative 2D chemical maps of the Al 2 O 3 and the ZnO. The STXM and ptychography maps were then measured at a sequence of tilt angles, covering up to 160° of rotation. The 3D structure of the ZnO and Al 2 O 3 was derived from the tilt series data by tomographic reconstruction using a compressed sensing algorithm. A twodimensional spatial resolution (half-period) of 6 nm, measured by Fourier ring correlation, and a 3D spatial resolution (half-period) of 9 nm, measured by Fourier shell correlation, were achieved when using the COSMIC beamline at the Advanced Light Source (ALS). The results show that for all of the ZnO loadings studied there is nonuniform coverage of the ZnO on the Al 2 O 3 aerogel framework. In addition, we found that both FIB methods create sample artifacts, although the distortion was less with Xe plasma than Ga ion FIB.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Simply Fabricatable Reference Electrode for Studying Li Metal Interfaces Operando

Reference electrodes are essential for studying Li metal interfaces, yet their integration into coin cells remains challenging due to geometric constraints and assembly complexity. Here, we present a simple three-electrode coin cell design in which a thin Cu film is evaporated directly onto a trilayer PP/PE/PP separator, forming a conductive and flexible in situ current collector. The microporous Cu preserves separator porosity and pressure uniformity, avoiding artifacts introduced by wire- or mesh-type references. It is easily activated by Li plating and provides stable potentials for over 10 days in ether-based electrolytes without affecting Coulombic efficiency (>99%) or Li deposition morphology. Operando measurements using this design reveal electrode-specific impedance evolution and voltage responses otherwise obscured in two-electrode setups, uncovering a correlation between the SEI resistance and electrolyte ionic conductivity. This microporous Cu reference electrode offers a robust and accessible platform for characterizing reactive interfaces under realistic cycling conditions.

Sacci, Robert [ORNL] (ORCID:0000000200735221)↗
DOE OSTI2025

Nanoscale Polarization-Dependent Young’s Modulus of Ferroelectric BaTiO 3 Single Crystals

Flexoelectric fields produced by strong strain gradients at the nanoscale couple to ferroelectric polarization, promoting changes in the mechanical properties of ferroelectric materials as a function of the direction of the ferroelectric polarization. In this work, we calculate the asymmetry in the Young’s modulus found in oppositely polarized out-of-plane domains of BaTiO 3 single crystals by means of contact resonance AFM, and we evaluate its impact on the electromechanical response as measured by piezoresponse force microscopy, both using band excitation modes. Furthermore, we analyze the electromechanical response of the different BaTiO 3 domains using k-means to build up their mean elastic and electromechanical features and quantify the flexoelectrically induced modulation of the Young’s modulus. Finally, we discuss the use of contact Kelvin probe force microscope measurements to decouple electrostatic artifacts from purely piezoelectric ones when flexoelectric coupling appears.

BaTiO 3 single crystals↗
DOE OSTI2026

Adsorption-Induced Surface Magnetism

We report the emergence of adsorption-induced magnetism from heterohelicene molecules on a nonmagnetic Cu(100) surface. Spin-polarized low-energy electron microscopy measurements reveal spin-dependent electron reflectivity for enantiopure 7,12,17-trioxa[11]helicene (TO[11]H) monolayers, indicating the formation of a spin-polarized state localized in the topmost copper layer. Control experiments on clean Cu(100) and TO[11]H on highly oriented pyrolytic graphite show no such effect, excluding artifacts and chirality-induced spin selectivity as origins. Spin-polarized density functional theory calculations with hybrid functionals attribute the magnetism to strong chemisorption, which induces hybridization between the molecular HOMO and copper s- and d-states, driving asymmetric spin-polarized charge redistribution at the interface. An extended Newns-Anderson-Grimley model incorporating on-site Coulomb repulsion in Cu d-orbitals reproduces the emergence of interfacial spin polarization above a threshold interaction strength, highlighting the key roles of hybridization parameters and Coulomb correlation. These findings reveal a mechanism for inducing magnetism at molecule-metal interfaces without inherently magnetic components, offering avenues for engineering spin-polarized states in organic-inorganic hybrid systems.

Hubbard model↗
DOE OSTI2026

Ca Metal Batteries and Sn Anode Alloying: Resolving Misconceptions in Ca–Sn Alloy Formation

Ca metal anode suffers from surface passivation and struggles to effectively plate and strip in conventional Ca electrolytes, making the development of alloy anode for calcium metal battery essential. This work systematically and carefully evaluates the electrochemical formation and feasibility of the Sn anode and its possibility of alloying with Ca using a combination of cyclic voltammetry, scanning electron microscopy with energy-dispersive X-ray spectroscopy (SEM/EDS), X-ray diffraction (XRD), and synchrotron transmission X-ray microscopy (TXM). Electrochemical measurements and SEM/EDS analyses indicate apparent Ca alloy formation. However, detailed XRD and TXM characterization reveal substantial lithium-ion involvement when lithium metal is employed as the reference electrode, even in electrolytes with high Ca salt concentrations. These findings highlight a significant risk of misinterpretation: Although limited Ca–Sn alloying may occur, the faster kinetics of Li + can dominate the alloying process, leading to a Li-driven reactions at the Sn anode. In contrast, cells constructed with Ca metal as both reference and counter electrodes exhibit distinct redox features only during the initial cycles, followed by rapid cell failure. The high oxidation currents induce severe Ca passivation, occurring nearly concurrently with Sn dealloying. Overall, this study provides critical insight into the practical limitations and common experimental artifacts associated with Ca metal battery systems and underscores the necessity for rigorous cell design and cautious interpretation when evaluating alloy chemistry in multivalent Ca-based batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dual-Color Optical Recording of Bioelectric Potentials by Polymer Electrochromism

Optical recording based on voltage-sensitive fluorescent reporters allows for spatial flexibility of measuring from desired cells, but photobleaching and phototoxicity of the fluorescent labels often limit their sensitivity and recording duration. Voltage-dependent optical absorption, rather than fluorescence, of electrochromic materials, would overcome these limitations to achieve long-term optical recording of bioelectrical signals. Electrochromic materials such as PEDOT:PSS possess the property that an applied voltage can either increase or decrease the light absorption depending on the wavelength. In this work, we harness this anticorrelated light absorption at two different wavelengths to significantly improve the signal detection. With dual-color detection, electrical activity from cells produces signals of opposite polarity, while artifacts, mechanical motions, and technical noises are uncorrelated or positively correlated. Using this technique, we are able to optically record cardiac action potentials with a high signal-to-noise ratio, 10 kHz sampling rate, >15 min recording duration, and no time-dependent degradation of the signal. Furthermore, we can reliably perform multiple recording sessions from the same culture for over 25 days.

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

Solvent Dynamics Are Critical to Understanding Carbon Dioxide Dissolution and Hydration in Water

Simulations of carbon dioxide (CO 2 ) in water may aid in understanding the impact of its accumulation in aquatic environments and help advance technologies for carbon capture and utilization (via, e.g., mineralization). Quantum mechanical (QM) simulations based on static molecular models with polarizable continuum solvation poorly reproduce the energetics of CO 2 hydration to form carbonic acid in water, independent of the level of QM theory employed. Only with density-functional-theory-based molecular dynamics and rare-event sampling, followed by energy corrections based on embedded correlated wavefunction theory (in conjunction with density functional embedding theory), can a close agreement between theory and experiment be achieved. Such multilevel simulations can serve as benchmarks for simpler, less costly models, giving insight into potential errors of the latter. The strong influence of sampling/averaging over dynamical solvent configurations on the energetics stems from the difference in polarity of both the transition state and product (both polar) versus the reactant (nonpolar). When a solute undergoes a change in polarity during reaction, affecting its interaction with the solvent, careful assessment of the energetic contribution of the solvent response to this change is critical. Here, we show that static models (without structural sampling) that incorporate three explicit water molecules can yield far superior results than models with more explicit water molecules because fewer water molecules yield less configurational artifacts. Static models intelligently incorporating both explicit (molecules directly participating in the reaction) and implicit solvation, along with a proper QM theory, e.g., CCSD(T) for closed-shell systems, can close the accuracy gap between static and dynamic models.

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

Reply to Comment by Peterie Et Al. on “Accelerated Fill‐Up of the Arbuckle Group Aquifer and Links to U.S. Midcontinent Seismicity”

Abstract Peterie et al. question one observation in our paper: associating pressure increases to injection volumes at distances of up to 25 km from an injection well. In this reply, we show that the comment misunderstands our analysis and the evidence that led to this conclusion. We also show that gauge‐depth‐corrected pressures, used by the authors to produce statewide pressure maps, are discrepant with the static fluid level data, provided in our original compilation and analysis. The discrepancies are a result of the pressure correction method employed, which naïvely substitutes formation pressure for bottomhole pressure to calculate wellbore fluid density. Their linearly interpolated pressure maps, based on sparse data, contain interpolation and extrapolation artifacts that contradict injection trends in the state, the Theis solution, and the superposition principle. We reiterate that pressure and static fluid level increases in Class I wells existed prior to 2013, most notably in central Kansas, where recent earthquakes are cited in the comment as evidence of a pressure plume emanating from the Kansas‐Oklahoma border, 90 km away. We show that the space‐time pattern of seismicity in this area is inconsistent with a northward propagating pressure plume and, instead, seismicity appears to be centered on and near a cluster of high‐rate injection wells, two of which are among the highest rate wells in the state. These observations, along with recent M4 + earthquakes during continued decreases in wastewater injection in southern Kansas and northern Oklahoma, question the usefulness of the comment for understanding and managing societally significant earthquakes.

Ansari, Esmail↗
DOE OSTI2020

Using Ambient Noise Seismic Interferometry and Local and Teleseismic Earthquakes to Determine Crustal Thickness and Moho Structure of the Northwestern Gulf of Mexico Margin

Abstract The northwestern part of the Gulf of Mexico has undergone two episodes of continental rifting and collision and produced structural artifacts that are now buried under many kilometers of sediments, complicating investigations of the region. The deep sedimentary package precludes outcrops and points to a need for the application of seismic techniques, but low rates of seismicity in the region and sparse seismic monitoring limit the utility of traditional seismic methods. We therefore use diverse data to perform two‐dimensional seismic tomography across the dry land portion of the margin. Data are gleaned from teleseismic and regional earthquakes, postcritical SsPmp arrivals, and direct P wave energy identified with seismic interferometry that were recorded by a broadband, three‐component array and partially overlapping short‐period, vertical‐component array. The Pn and postcritical SsPmp phase help constrain the Moho discontinuity, which a previous receiver function study suggested was absent beneath the seaward portion of this transect. A high‐velocity body is observed in the crust at the same location as the Houston Magnetic Anomaly, possibly marking rocks from the Alleghenian continental assembly. The crust thins from NW to SE, indicating that extension occurred mostly to the south of the Ouachita orogeny. Our model indicates that the margin's sediment package reaches a maximum thickness of ~15 km at the coast and becomes unresolvably thin near the Llano Uplift.

58 GEOSCIENCES↗
DOE OSTI2022

Idealized Simulations of the Tropical Climate and Variability in the Single Column Atmosphere Model (SCAM): Radiative-Convective Equilibrium

To explore the interactions among column processes in the Community Atmosphere Model (CAM), the single-column version of CAM (SCAM) is integrated for 1000 days in radiative-convective equilibrium (RCE) with tropical values of boundary conditions, spanning a parameter or configuration space of model physics versions (v5 vs. v6), vertical resolution (standard and 60 levels), sea surface temperature (SST), and some interpretation-driven experiments. The simulated time-mean climate is reasonable, near observations and RCE of a cyclic cloud-resolving model. Updraft detrainment in the deep convection scheme produces distinctive grid-scale structures in humidity and cloud, which also interact with radiative transfer processes. These grid artifacts average out in multi-column RCE results reported elsewhere, illustrating the nuts-and-bolts interpretability that SCAM adds to the hierarchy of model configurations. Multi-day oscillations of precipitation arise from descent of warm convection-capping layers starting near the tropopause, eventually reset by a burst of convective deepening. Experiments reveal how these oscillations depend critically on an internal parameter that controls the number of neutral buoyancy levels allowed for determining cloud top and computing dilute convective available potential energy in the deep convection scheme, and merely modified a little by disabling cloud-base radiation (heating of cloud base). This strong dependence of transient behavior in 1D on this parameter will be tested in the second part of this work, in which SCAM is coupled to a parameterized dynamics of two-dimensional, linearized gravity wave, and in the 3D simulations in future study.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Load‐Balancing Intense Physics Calculations to Embed Regionalized High‐Resolution Cloud Resolving Models in the E3SM and CESM Climate Models

Abstract We design a new strategy to load‐balance high‐intensity sub‐grid atmospheric physics calculations restricted to a small fraction of a global climate simulation's domain. We show why the current parallel load balancing infrastructure of Community Earth System Model (CESM) and Energy Exascale Earth Model (E3SM) cannot efficiently handle this scenario at large core counts. As an example, we study an unusual configuration of the E3SM Multiscale Modeling Framework (MMF) that embeds a binary mixture of two separate cloud‐resolving model grid structures that is attractive for low cloud feedback studies. Less than a third of the planet uses high‐resolution (MMF‐HR; sub‐km horizontal grid spacing) relative to standard low‐resolution (MMF‐LR) cloud superparameterization elsewhere. To enable MMF runs with Multi‐Domain cloud resolving models (CRMs), our load balancing theory predicts the most efficient computational scale as a function of the high‐intensity work's relative overhead and its fractional coverage. The scheme successfully maximizes model throughput and minimizes model cost relative to precursor infrastructure, effectively by devoting the vast majority of the processor pool to operate on the few high‐intensity (and rate‐limiting) high‐resolution (HR) grid columns. Two examples prove the concept, showing that minor artifacts can be introduced near the HR/low‐resolution CRM grid transition boundary on idealized aquaplanets, but are minimal in operationally relevant real‐geography settings. As intended, within the high (low) resolution area, our Multi‐Domain CRM simulations exhibit cloud fraction and shortwave reflection convergent to standard baseline tests that use globally homogenous MMF‐LR and MMF‐HR. We suggest this approach can open up a range of creative multi‐resolution climate experiments without requiring unduly large allocations of computational resources.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Analysis of Narrow Bipolar Events Using Mode Decomposition Methods

Multi-resolution analysis methods can reveal the underlying physical dynamics of nonstationary signals, such as those from lightning. In this paper we demonstrate the application of two multi-resolution analysis methods: Ensemble Empirical Mode Decomposition (EEMD) and Variational Mode Decomposition (VMD) in a comparative way in the analysis of electric field change waveforms from lightning. EEMD and VMD decompose signals into a set of Intrinsic Mode Functions (IMFs). The IMFs can be combined using distance and divergence metrics to obtain noise reduction or to obtain new waveforms that isolate the physical processes of interest while removing irrelevant components of the original signal. We apply the EEMD and VMD methods to the observations of three close Narrow Bipolar Events (NBEs) that were reported by Rison et al. The ΔE observations reveal the occurrence of complex oscillatory processes after the main NBE sferic. We show that both EEMD and VMD are able to isolate the oscillations from the main NBE, with VMD being more effective of the two methods since it requires the least user supervision. The oscillations are found to begin at the end of the NBEs' downward fast positive breakdown and appear to be produced by a half-wavelength standing wave within a weakly-conducting resonant ionization cavity left behind in the wake of the streamer-based NBE event. In conclusion, additional analysis shows that one of the NBEs was likely initiated by an energetic cosmic ray shower, and also corrects a misinterpretation in the literature that fast breakdown is an artifact of NBE-like events in interferometer observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Removing Numerical Pathologies in a Turbulence Parameterization Through Convergence Testing

Abstract Discretized numerical models of the atmosphere are usually intended to faithfully represent an underlying set of continuous equations, but this necessary condition is violated sometimes by subtle pathologies that have crept into the discretized equations. Such pathologies can introduce undesirable artifacts, such as sawtooth noise, into the model solutions. The presence of these pathologies can be detected by numerical convergence testing. This study employs convergence testing to verify the discretization of the Cloud Layers Unified By Binormals (CLUBB) model of clouds and turbulence. That convergence testing identifies two aspects of CLUBB's equation set that contribute to undesirable noise in the solutions. First, numerical limiters (i.e., clipping) used by CLUBB introduce discontinuities or slope discontinuities in model fields. Second, nonlinear artificial diffusion employed for improving numerical stability can introduce unintended small‐scale features into the solution of the model equations. Smoothing the limiters and using linear artificial diffusion reduces the noise and restores the expected first‐order convergence in CLUBB's solutions. These model reformulations enhance our confidence in the trustworthiness of solutions from CLUBB by eliminating the unphysical oscillations in high‐resolution simulations. The improvements in the results at coarser, near‐operational grid spacing and timestep are also seen in cumulus cloud and dry turbulence tests. In addition, convergence testing is proven to be a valuable tool for detecting pathologies, including unintended discontinuities and grid dependence, in the model equation set.

58 GEOSCIENCES↗