Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sawtooth Oscillations”

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.

16 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Core electron temperature turbulence and transport during sawtooth oscillations in the DIII-D tokamak

Abstract Sawteeth are one of the concerning instabilities in ITER and future burning plasma experiments. Sawtooth dynamics and its interaction with broadband plasma turbulence has been a challenge for predictive simulations of core transport in future fusion devices. This study provides new observations of core turbulence behavior during sawtooth oscillations in DIII-D hydrogen L-mode neutral beam injection heated plasmas in an inner wall limited configuration. A strong correlation of electron temperature and density turbulence levels with the sawtooth oscillation phase has been observed at locations inside the T e inversion radius and/or safety factor q = 1 magnetic surface. The T e turbulence amplitude in the core during the sawtooth ramp exhibits a critical T e gradient behavior inside but not near the T e inversion radius/ q = 1 magnetic surface. The most unstable mode calculated from the trapped gyro-landau fluid turbulence simulations reveal a change from low- k ion-type to low-k electron-type modes from pre- to post- sawtooth crash time periods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Power-transfer and fixed-point analysis of sawtooth simulations of a current-carrying stellarator

Power-transfer and fixed-point analysis of previous NIMROD simulations (Roberds et al. , Phys. Plasmas , vol. 23, issue 9, 2016, 092513) improved the understanding of the effect of 3D (non-axisymmetric equilibrium) magnetic fields on sawtooth oscillations in the Compact Toroidal Hybrid (CTH) experiment. Computing the locations of order-1 fixed points, their Greene's residues, and local values for the rotational transform results in a description of CTH sawteeth consistent with Kadomtsev's model. A power-transfer analysis quantifies the distribution of energy among toroidal Fourier modes and their nonlinear interactions. The Lorentz power transfer drives sawtooth growth, and it is unambiguously interpreted as the flow of energy from toroidal mode $n'$ to mode $n$ , catalysed by $\boldsymbol {B}_{n-n'}$ . It has been reported previously that the CTH sawtooth frequency increases with the 3D field strength. This is attributed to an increased growth rate of the internal kink that drives sawtooth oscillations. Here, 3D fields remove energy from the kink, eliminating the possibility that these fields are an additional energy source that drives growth. Instead, 3D fields catalyse energy transfer from large-to-small scales, where magnetic reconnection is stronger. It is proposed that this energy transfer increases the reconnection rate at small scales, which is consistent with the increased growth rate observed at higher 3D field strengths.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

3D Equilibrium Reconstruction and MHD Studies in Stellarators (Final Report)

This is the final technical report for the grant Three Dimensional Equilibrium Reconstruction and MHD Studies in Stellarators. The major goals of this project were to: further develop V3FIT code capabilities, continue to support experiments using the V3FIT code, continue theoretical support for the CTH experiment at Auburn, particularly in the area of NIMROD simulations of CTH discharges. Dr. Lopez Ortiz has made significant progress on his analysis of NIMROD simulations of CTH MHD phenomena. In particular, his study has focused on the fixed point dynamics and energy flow amongst toroidal Fourier modes in sawtooth oscillations simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Velocity-space distribution function of fast ions in a sawtoothing plasma

This study explores the influence of sawtooth oscillations on the velocity space distribution of fast ions in tokamak plasma discharges. The relevant Fokker–Planck equation for fast ions is solved analytically. Two distinct effects arising from the temperature drop associated with a sawtooth crash and their impact on the distribution function of fast ions are considered. The first effect involves the modulation of the fusion alpha particle source on the timescale of the sawtooth period, linked to the drop in fusion yield resulting from the sawtooth temperature relaxations. The second effect is tied to the increase of the slowing-down time during the sawtooth ramp, causing particles born later in the sawtooth cycle to experience reduced slowing down compared to those born right after the crash, creating an accumulation-like mechanism at higher energies. In regimes where the sawtooth period is shorter than the fast ion slowing-down time, the combined influence of these effects gives rise to fast ion distribution functions that transiently exhibit positive slopes in velocity space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Statistical and temporal characteristics of sawtooth events

Magnetospheric sawtooth events are characterized by periodic particle injections and magnetic dipolarizations spread quasi-simultaneously across a wide range of magnetic local times. We present a comprehensive statistical study of magnetospheric sawtooth events (STEs) during solar cycle 24 (2008–2016), extending previous catalogs and enabling solar cycle comparisons. Our results confirm that STEs predominantly occur during the rising and declining phases of the solar cycle, and are strongly associated with geomagnetic storms. Superposed epoch analysis reveals near-simultaneous particle injections across all magnetic local time sectors, but magnetic field dipolarization confined to the midnight region. These results support a scenario in which nightside tail reconnection and enhanced convection are the primary drivers of sawtooth oscillations. The localization of magnetic dipolarizations during STEs challenges global instability interpretations and suggests that STEs represent a stormtime substorm mode triggered under specific solar wind and magnetotail conditions. Superposed epoch analyses also show enhanced oxygen content in the magnetosphere during sawtooth events, but do not show a significant difference from geomagnetic storms that do not exhibit periodic behavior.

58 GEOSCIENCES↗
DOE OSTI2022

Computational study of runaway electrons in MST tokamak discharges with applied resonant magnetic perturbation

A numerical study of magnetohydrodynamics (MHD) and tracer-particle evolution investigates the effects of resonant magnetic perturbations (RMPs) on the confinement of runaway electrons (REs) in tokamak discharges conducted in the Madison Symmetric Torus. In computational results of applying RMPs having a broad toroidal spectrum but a single poloidal harmonic, m = 1 RMP does not suppress REs, whereas m = 3 RMP achieves significant deconfinement but not the complete suppression obtained in the experiment. MHD simulations with the NIMROD code produce sawtooth oscillations, and the associated magnetic reconnection can affect the trajectory of REs starting in the core region. Simulations with m = 3 RMP produce chaotic magnetic topology over the outer region, but the m = 1 RMP produces negligible changes in field topology, relative to applying no RMP. Using snapshots of the MHD simulation fields, full-orbit relativistic electron test particle computations with KORC show ≈50% loss from the m = 3 RMP compared to the 10%–15% loss from the m = 1 RMP. Here, test particle computations of the m = 3 RMP in the time-evolving MHD simulation fields show correlation between MHD activity and late-time particle losses, but total electron confinement is similar to computations using magnetic-field snapshots.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Impact of suprathermal ions on neutron yield in the pre-DT phase of ITER operation

An assessment of neutron production during the pre-fusion power-operation (PFPO) phase has been carried out for a representative set of plasma scenarios predicted by the ITER Research Plan. A range of heating systems, namely neutral beam injection (NBI) (hydrogen), electron cyclotron resonance heating (ECRH), and ion cyclotron resonance heating (ICRH) are planned to be used for PFPO studies in helium, hydrogen, and mixed hydrogen–helium plasmas. Fast ions (protons and 3He) originating from NBI and ICRH systems can increase neutron production in PFPO plasmas by directly interacting with intrinsic Be impurities or through secondary processes, as also evidenced at JET. The generation of fast ions in ITER PFPO scenarios has been modelled using the ASTRA-NBI and TORIC-SSFPQL codes. A significant impact of the synergy between hydrogen NBI and hydrogen-minority ICRH on neutron production in helium plasmas is reported. In addition, the stability of the toroidicity-induced Alfvén eigenmodes (TAE) is analyzed for PFPO plasmas with a high pressure of suprathermal ions and a weak reversed shear. The possible impact of sawtooth oscillations and TAEs on neutron production is discussed, based on a linear stability analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

On FIRE mode in KSTAR

We report on the status of Fast Ion Regulated Enhancement (FIRE) mode experiments in the Korea Superconducting Tokamak Advanced Research. This regime is being developed for high-performance, steady-state operation which features a stationary ion internal transport barrier, enabling a central ion temperature approaching 10 keV to be sustained for up to 50 s, without the need for delicate profile control and with no significant impurity accumulation. As its key novelty lies in the significant contribution of fast ions that stabilize core turbulence, the regime has been named FIRE mode. To achieve this regime, neutral beam injection is applied at moderate power levels near the L–H power threshold, while maintaining low plasma density to avoid the L–H transition. The scenario is typically established in diverted magnetic configurations. The core features of FIRE mode were investigated through power balance analysis and fluctuation measurements, revealing a clear transport bifurcation in the ion channel. At the plasma edge, FIRE mode occasionally exhibits I-mode characteristics, particularly in unfavorable magnetic null configurations with q 95 ∼ 4, including the presence of weakly coherent modes. In terms of MHD activity, sawtooth oscillations are observed but appear to be stabilized during the high-performance phase. Fast ion-driven Alfvénic eigenmodes (AEs), indicated by strong frequency chirping near 200 kHz, are also observed. Additionally, lower-frequency MHD activities, distinct from the fast ion-driven AEs, are present and have some influence on plasma performance. The enhancement of core confinement is primarily attributed to fast ion effects, which were evaluated with respect to dilution, alpha stabilization, and resonant interactions with turbulence using gyrokinetic analyses. Among these, the dilution effect was found to be the most dominant. The characteristics of FIRE mode were compared with those of other hot ion plasma scenarios, such as supershot and hot ion mode. While they share many similarities, FIRE mode is distinguished by the accessibility to conditions with T i ≈ T e , presence of I-mode edge features, and its long-duration sustainment. Predictive simulations of FIRE mode were performed using integrated transport modeling with TRIASSIC, employing the TGLF anomalous transport model. These simulations confirmed the critical role of fast ions in achieving this regime. The future prospect of FIRE mode for application in fusion reactors is discussed, with an emphasis on possibly extending the regime to higher density operation.

FIRE mode↗
DOE OSTI2021

Validation of the EDGES Low-band Antenna Beam Model

The response of the antenna is a source of uncertainty in measurements with the Experiment to Detect the Global Epoch of Reionization Signature (EDGES). We aim to validate the electromagnetic beam model of the low-band (50–100 MHz) dipole antenna with comparisons between models and against data. We find that simulations of a simplified model of the antenna over an infinite perfectly conducting ground plane are, with one exception, robust to changes in the numerical electromagnetic solver code or algorithm. For simulations of the antenna with the actual finite ground plane and realistic soil properties, we find that two out of three numerical solvers agree well. Applying our analysis pipeline to a simulated drift-scan observation from an early EDGES low-band instrument that had a 10 m × 10 m ground plane, we find residual levels after fitting and removing a five-term foreground model from the simulated data binned in local sidereal time (LST) average about 250 mK with ±40 mK variation between numerical solvers. A similar analysis of the primary 30 m × 30 m sawtooth ground plane reduced the LST-averaged residuals to about 90 mK with ±10 mK between the two viable solvers. More broadly we show that larger ground planes generally perform better than smaller ground planes. Simulated data have a power that is within 4% of real observations, a limitation of net accuracy of the sky and beam models. We observe that residual spectral structures after foreground model fits match qualitatively between simulated data and observations, suggesting that the frequency dependence of the beam is reasonably represented by the models. We find that a soil conductivity of 0.02 S m{sup -1} and relative permittivity of 3.5 yield good agreement between simulated spectra and observations. This is consistent with the soil properties reported by Sutinjo et al. for the Murchison Radio-astronomy Observatory, where EDGES is located.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

A New Plasma Radar Concept for Simultaneous Magnetic and Density Measurements

An innovative, compact 288GHz interferometer has been fabricated, tested, installed and successfully demonstrated on the LAPD-U magnetized plasma at UCLA. The system takes advantage of frequency modulated (FM) radar techniques to deliver a compact heterodyne system. In addition, the reflected power from the source is taken advantage of to eliminate the need for additional quasi-optical components. Electron density in LAPD-U plasma has recently been increased substantially thereby requiring a higher frequency/shorter wavelength interferometer to avoid deleterious refractive effects. This system satisfies those needs. The system uses a 96GHz varactor tuned Gunn oscillator which passes to a passive tripler. This tripler has ~3% conversion efficiency. The 288GHz radiation is then coupled to free space using a so-called dual-mode or Pickett horn. The output 288GHz beam is then coupled to an aspheric lens manufactured from low-loss, high-density polyethylene. This lens is employed to collimate the emerging beam. Small axial adjustment of the lens position can also be used to create a slowly focusing beam so as to optimize the measured signal. In addition, up-down or side-to-side adjustment of the lens can be utilized to steer the beam vertically or horizontally – again to optimize alignment. The propagating beam passes through a beam splitter and then through a water-free, bubble-free fused quartz window into the LAPD-U vacuum vessel. The beam-splitter is a thin sheet of G10 which reflects a small fraction of the incident power (~5 %) towards a zero-bias detector optimized for the frequency range from 220 to 300GHz. Note that waveguides at this frequency have dimensions of ~0.9mm x 0.45mm and so have very large conductive losses. This drives the use of quasi-optical propagation. The detector requires no DC bias and is very responsive (> 1V/mW into 1MΩ). Radiation is coupled to the detector via a similar lens-horn arrangement used for the launch. This reflected beam acts as the local oscillator or reference millimeter-wave beam for the detector. The remainder of the launched source beam then enters the LAPD-U vacuum vessel and passes through the plasma at the mid-plane until reaching the opposing port which is closed off with an aluminum flange. This flange is used as a mirror to retroreflect the incident 288GHz beam back along its path. The retroreflected beam exits the input port but does NOT couple directly into the zero-bias detector. Instead, the majority of the return power continues towards the 288GHz source. As mentioned above the transmitted beam enters the source a second time. This would appear undesirable. However, at these frequencies multipliers are highly non-linear elements which results in a significant portion of the return beam (~20%) re-emerging from the multiplier and horn and then coupling via the G10 beam-splitter to the zero-bias detector. This approach eliminated the need for a second quasi-optical beam-splitter. The system is extremely compact measuring approximately 28 inches x 20 inches. The above did not explain how heterodyne operation was achieved. As mentioned above the Gunn oscillator is able to be varactor tuned. This allows a low voltage to be applied to control the operating frequency of the Gunn oscillator. During heterodyne operation a sawtooth shaped voltage is applied to the varactor at 750kHz using an 80MHz Arbitrary Waveform Generator (AWG). This voltage changes the Gunn frequency linearly during the up-sweep which is then reset abruptly at the sawtooth crash to be immediately followed by another linear sweep. Passage through the 288GHz multiplier triples the frequency change experienced by the electromagnetic wave. These frequency changes are small – tens of megahertz. This FM radar approach results in the launched electromagnetic wave frequencies at the detector for the reference and plasma wave to be different. The approximately 10 ns delay propagation delay for the plasma beam results in the local oscillator and plasma beams NOT having an identical frequency – there is in fact a fixed difference frequency. The frequency tuning level of the Gunn oscillator is then adjusted so that there is ONE cycle of this difference frequency during each linear ramp. During the sawtooth crash or downward re-sweep this one cycle replays in reverse but on a very fast timescale. The process then repeats. Low-pass filtering eliminates the fast re-sweep to leave a pure sine wave heterodyne signal. When the plasma is present it introduces a phase delay in the sine wave (caused by the extremely small Doppler shift resulting from the optical path length change). Of course, to measure this phase change we need a reference. This is simply obtained from an arbitrary waveform generator which provides a synchronized output pulse train which again is low pass filtered to obtain a 750kHz sinusoidal voltage reference for the interferometer. The interferometer was installed on LAPD-U where it has worked reliably and has established that electron densities exceeding 1x10 13 cm -3 are routinely achieved. In addition, the system sensitivity was able to easily observe density fluctuation at frequencies up to 50kHz. FM Radar techniques have enabled a full demonstration of a compact, sensitive, high frequency (288GHz/1mm) heterodyne interferometer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

On discriminating tokamak sawtooth crash models via localized density and temperature measurements

The core electron temperature drops rapidly during the sawtooth crash in tokamak plasmas, which causes heat loss and may lead to fast particle losses or even a disruption. Several models have been proposed for the periodic crash, including the Kadomtsev model with magnetic reconnection and the quasi-interchange model with the growth of higher-mode-number pressure-driven instabilities. 3D MHD simulations were performed for these two models with a goal to develop intuition and to predict qualitatively how different types of sawtooth will appear in various diagnostics. The structures of electron density ne and electron temperature Te show a dominant (1, 1) mode for the Kadomtsev case and a dominant (4, 4) mode for the quasi-interchange case. The oscillations of ne and Te have a positive correlation near the inversion layer for both cases, while their frequencies and amplitudes are different depending on the dominant modes. Particularly, for the Kadomtsev case, we find a relation between the amount of flux reconnected during a sawtooth event and ne or Te oscillations. Therefore, we connect recently developed measurement capabilities for ne and Te to the internal sawtooth behavior. We propose that this method of analysis can help in identifying the type of sawtooth in future experiments augmented by simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Two-dimensional plasma density evolution local to the inversion layer during sawtooth crash events using Beam Emission Spectroscopy

Here, we present methods for analyzing Beam Emission Spectroscopy (BES) data to obtain the plasma density evolution associated with rapid sawtooth crash events at the DIII-D tokamak. BES allows coverage over a 2D spatial plane, inherently local measurements, with fast time responses, and, therefore, provides a valuable new channel for data during sawtooth events. A method is developed to remove sawtooth-induced edge-light pulses contained in the BES data. The edge light pulses appear to be from the D α emission produced by edge recycling during sawtooth events, and are large enough that traditional spectroscopic filtering and data analysis techniques are insufficient to deduce physically meaningful quantities. A cross-calibration of 64 BES channels is performed by using a novel method to ensure accurate measurements. For the large-amplitude density oscillations observed, we discuss and use the non-linear relationship between the BES signal δI/ I 0 and the plasma density variation δn e / n e0 . The 2D BES images cover an 8 × 20 cm 2 region around the sawtooth inversion layer and show large-amplitude density oscillations, with additional significant spatial variations across the inversion layer that grows and peaks near the time of the temperature crash. The edge light removal technique and method of converting large-amplitude δI/ I 0 to δn e / n e0 presented here may help analyze other impulsive MHD phenomena in tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Overview of physics results from the ADITYA-U tokamak and future experiments

The ADITYA upgrade (ADITYA-U), a medium-sized (R 0 = 75 cm, a = 25 cm) conventional tokamak facility in India, has been consistently producing experiments findings by usingcircular and shaped-plasmas. Recognizing the plasma parameters aligning closely with the design parameters of circular limited plasmas, ADITYA-U shifted its focus toward exploring the operational regime for experimentation on saw-tooth and MHD phenomena. Moreover, ADITYA-U has made consistent advancements toward conducting preliminary plasma shaping experiments through the activation of top and bottom divertor coils utilizing hydrogen as well as deuterium fuels. Confinement is improved by a factor of ~1.5 in D 2 plasmas when compared to H 2 plasmas of ADITYA-U. Further, ADITYA-U operations emphasize preventing disruptions and runaway electrons (REs) to ensure safe operations for future fusion devices. Significant suppression of REs has been achieved in ADITYA-U with the application of pulsed localized vertical magnetic field (LVF) perturbation, thereby establishing the technique's independence from the tokamak device. The successful RE mitigation requires a critical threshold of LVF pulse magnitude, which is approximately 1% of the toroidal magnetic field, and a minimum duration of ~5 ms. Apart from this, several novel findings have been achieved in the ADITYA-U experiments, including the modification of sawtooth duration through gas-puff, the emergence of MHD-induced geodesic acoustic mode-like oscillations, the propagation of fast heat pulses induced by MHD activity, the control of RE dynamics through Gas-puffs, the propagation of pinch-driven cold-pulses, the transport and core accumulations of argon impurities, the mass dependency of plasma toroidal rotation and the detection of 'RICE' scaling, as well as the characterization of edge plasma using wall conditioning methods, such as glow discharge cleaning using a combination of Ar-H 2 mixture, localized wall cleaning by electron cyclotron resonant plasma, and the development of machine learning-based disruption predictions, will be discussed in this paper.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Synthesis of Quantum-Confined Borophene Nanoribbons

Here, borophene nanoribbons (BNRs) are one-dimensional strips of atomically thin boron expected to exhibit quantum-confined electronic properties that are not present in extended two-dimensional borophene. While the parent material borophene has been experimentally shown to possess anisotropic metallicity and diverse polymorphic structures, the atomically precise synthesis of nanometer-wide BNRs has not yet been achieved. Here, we demonstrate the synthesis of multiple BNR polymorphs with well-defined edge configurations within the nanometer-scale terraces of vicinal Ag(977). Through atomic-scale imaging, spectroscopy, and first-principles calculations, the synthesized BNR polymorphs are characterized and found to possess distinct edge structures and electronic properties. For single-phase BNRs, v 1/6 -BNRs and v 1/5 -BNRs adopt reconstructed armchair edges and sawtooth edges, respectively. In addition, the electronic properties of single-phase v 1/6 -BNRs and v 1/5 -BNRs are dominated by Friedel oscillations and striped moiré patterns, respectively. On the other hand, mixed-phase BNRs possess quantum-confined states with increasing nodes in the electronic density of states at elevated biases. Overall, the high degree of polymorphism and diverse edge topologies in borophene nanoribbons provide a rich quantum platform for studying one-dimensional electronic states.

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

Dependence of ion-cyclotron range of frequencies instabilities on species mix and fast-ion distribution: II. Time evolution

Dedicated DIII-D experiments studied instabilities below the ion cyclotron frequency f ci and at cyclotron harmonics by injecting thirteen different neutral beam populations into plasmas with different thermal compositions of H and D at five different values of toroidal field. The resulting database of 2529 conditions is described in Part 1: Stability. This paper focuses on the temporal evolution. For some beam geometries, the ion cyclotron emission (ICE) is much stronger at the beginning of injection, suggesting the importance of ‘bump-on-tail’ energy drive. Sawtooth crashes cause the emission to cease for some beam types and emission bands, but the opposite effect occurs for others. When unstable, for virtually all conditions and frequency bands, the emission bursts in a cycle reminiscent of predator-prey oscillations. Evidence for nonlinear interaction between ICE harmonics is observed. Strongly unstable modes typically appear within 1 ms from the start of beam injection and often decay within 0.1 ms when the beam turns off, but both the initial rise and final decay of the mode power exhibit a wide variety of behavior.

electromagnetic ion cyclotron instabilities↗
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↗