Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “abstraction layers”

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

Sweeping Jet Film Cooling at High Blowing Ratio on a Turbine Vane

Abstract Experimental and numerical investigations were performed to study the effects of high blowing ratios and high freestream turbulence on sweeping jet film cooling. Experiments were conducted on a nozzle guide vane suction surface in a low-speed linear cascade at a range of blowing ratios of 0.5–3.5 and freestream turbulence of 0.6% and 14.3%. Infrared thermography was used to estimate the adiabatic cooling effectiveness. Thermal field and boundary layer measurements were conducted at a cross-plane at x/D = 12 downstream of the hole exit. Results were compared with a baseline 777-shaped hole and showed that the sweeping jet hole has improved cooling effectiveness at high blowing ratios (M > 1). The thermal field data revealed that the coolant separates from the surface at high blowing ratios for the 777-shaped hole while the coolant remains attached for the sweeping jet hole. Boundary layer measurements further confirmed that due to the sweeping motion of the jet, the effective jet momentum of the sweeping jet hole remains much lower than that of a 777-shaped hole. Thus, the coolant remains closer to the wall even at high blowing ratios. Large Eddy simulations (LES) were performed for both the sweeping jet and 777-shaped hole to evaluate the interaction between the coolant jet and the freestream in the near hole regions. Results showed that 777-shaped hole has a strong jetting at high blowing ratio that originates inside the hole breakout edges thus causing the jet to blow-off from the surface. In contrast, the sweeping jet does not show this behavior due to its internal geometry and flapping motion of the jet.

Engineering↗
DOE OSTI2022

The Influence of Shallow Cloud Populations on Transitions to Deep Convection in the Amazon

Abstract In this study, a pair of convection-permitting (2-km grid spacing), month-long, wet-season Weather Research and Forecasting (WRF) Model simulations with and without the eddy-diffusivity mass-flux (EDMF) scheme are performed for a portion of the Green Ocean Amazon (GoAmazon) 2014/15 field campaign period. EDMF produces an ensemble of subgrid-scale convective plumes that evolve in response to the boundary layer meteorological conditions and can develop into shallow clouds. The objective of this study is to determine how different treatments of shallow cumulus clouds (i.e., with and without EDMF) impact the total cloud population and precipitation across the Amazonian rain forest, with emphasis on impacts on the likelihood of shallow-to-deep convection transitions. Results indicate that the large-scale synoptic conditions in the EDMF and control simulations are nearly identical; however, on the local scale their rainfall patterns diverge drastically and the biases decrease in EDMF. The EDMF scheme significantly increases the frequency of shallow clouds, but the frequencies of deep clouds are similar between the simulations. Deep convective clouds are tracked using a cloud-tracking algorithm to examine the impact of shallow cumulus on the surrounding ambient environment where deep convective clouds initiate. Results suggest that a rapid increase of low-level cloudiness acts to cool and moisten the low to midtroposphere during the day, favoring the transition to deep convection.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2003

[Results for SHEBA/FIRE]

The Atmospheric Research Laboratory's Radiation Measurement System (RAMS) was on the NCAR C-130 aircraft in May and July 1998, collecting radiometric data on the science flights conducted in the vicinity of the Surface Heat Budget of the Arctic Ocean (SHEBA) ship. These measurements were part of the FIRE Arctic Cloud Experiment (FIRE ACE). Analysis of some of the data focused on the absorption, reflection, and transmittance of Arctic clouds, especially compared to model results. In order to assess the absorption of solar radiation by the clear and cloudy atmosphere in the Arctic the measurements from the radiometers were combined in pairs of above-cloud segments and below-cloud segments. To get these pairs, the data for all sixteen of the flights (8 in May and 8 in July) were examined for occurrences of low-altitude segments in proximity to high-altitude segments. The low-altitude data are then treated as measurements of the bottom of a layer and the high-altitude data are taken as measurements of the top of the layer. With measurements of the upwelling and downwelling irradiances above and below a layer one can determine the reflectance, transmittance, and absorptance of the layer. Attachment: Doelling, D.R., P. Minnis, D.A. Spangenberg, V. Chakrapani, A. Mahesh, S.K. Pope, and F.P.J. Valero, Cloud radiative forcing at the top of the atmosphere during FIRE ACE derived from AVHRR data, J. Geophys. Res. 106, 15,279-15,296,2001. Minnis, P., D.R. Doelling, D.A. Spangenberg, A. Mahesh, S.K. Pope, and F.P.J. Valero, AVHRR-derived cloud radiative forcing over the ARM NSA and SHEBA site during FIRE ACE, abstract submitted to the ARM Science Team Meeting, San Antonio, TX, M a . 13-17,2000. Pope, S.K., and F.P.J. Valero, Measured and modeled radiometric fluxes in the Arctic during FIRE-ACE, presented as a poster at the American Geophysical Union meeting, San Francisco, CA, Dec. 13-17, 1999. Pope, S.K., and F.P.J. Valero, Measured and modeled radiometric fluxes in the Arctic during FIRE-ACEy paper presented at SHEBA/FIRE Workshop, National Center for Atmospheric Research, Boulder, Colo., Apr. 17-20,2000.

Valero, Francisco P. J.↗
DOE OSTI2023

Gate-tunable superconducting diode effect in a three-terminal Josephson device

Abstract The phenomenon of non-reciprocal critical current in a Josephson device, termed the Josephson diode effect, has garnered much recent interest. Realization of the diode effect requires inversion symmetry breaking, typically obtained by spin-orbit interactions. Here we report observation of the Josephson diode effect in a three-terminal Josephson device based upon an InAs quantum well two-dimensional electron gas proximitized by an epitaxial aluminum superconducting layer. We demonstrate that the diode efficiency in our devices can be tuned by a small out-of-plane magnetic field or by electrostatic gating. We show that the Josephson diode effect in these devices is a consequence of the artificial realization of a current-phase relation that contains higher harmonics. We also show nonlinear DC intermodulation and simultaneous two-signal rectification, enabled by the multi-terminal nature of the devices. Furthermore, we show that the diode effect is an inherent property of multi-terminal Josephson devices, establishing an immediately scalable approach by which potential applications of the Josephson diode effect can be realized, agnostic to the underlying material platform. These Josephson devices may also serve as gate-tunable building blocks in designing topologically protected qubits.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Uncovering the spin ordering in magic-angle graphene via edge state equilibration

Abstract The flat bands in magic-angle twisted bilayer graphene (MATBG) provide an especially rich arena to investigate interaction-driven ground states. While progress has been made in identifying the correlated insulators and their excitations at commensurate moiré filling factors, the spin-valley polarizations of the topological states that emerge at high magnetic field remain unknown. Here we introduce a technique based on twist-decoupled van der Waals layers that enables measurement of their electronic band structure and–by studying the backscattering between counter-propagating edge states–the determination of the relative spin polarization of their edge modes. We find that the symmetry-broken quantum Hall states that extend from the charge neutrality point in MATBG are spin unpolarized at even integer filling factors. The measurements also indicate that the correlated Chern insulator emerging from half filling of the flat valence band is spin unpolarized and suggest that its conduction band counterpart may be spin polarized.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Selective laser sintering and spark plasma sintering of (Zr,Nb,Ta,Ti,W)C compositionally complex carbide ceramics

Abstract Two advanced manufacturing processes, spark plasma sintering (SPS) and selective laser sintering (SLS), have been developed for synthesis of (Zr,Nb,Ta,Ti,W)C compositionally complex carbide (CCC) via reactive sintering of a powder mixture of constitute monocarbides. X‐ray diffraction analysis confirmed that the single‐phase CCC can be formed by both SPS and SLS. While a homogenous microstructure with uniform metal element distributions was developed during SPS, three‐layer microstructures with a thin TiC‐rich layer and two TaC‐rich layers along with a TiO 2 ‐rich surface layer containing W nanoparticles were formed during SLS. In addition, cellular structures with W, Zr, and Ti element segregation and dislocations on cell boundaries were observed in the SLS‐CCC sample, indicating the effect of nonequilibrium conditions on microstructure formation during laser melting followed by rapid cooling and solidification process. Compared to the SPS‐CCC sample, the SLS‐CCC showed enhanced hardness and reduced thermal conductivity, which may be related to their unique cellular structures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Measurements from mobile surface vehicles during the Lower Atmospheric Profiling Studies at Elevation – a Remotely-piloted Aircraft Team Experiment (LAPSE-RATE)

Abstract. Between 14 and 20 July 2018, small unmanned aircraft systems (UASs) were deployed to the San Luis Valley of Colorado (USA) alongside surface-based remote sensors, in situ sensors, and radiosonde systems as part of the Lower Atmospheric Profiling Studies at Elevation – a Remotely-piloted Aircraft Team Experiment (LAPSE-RATE). The measurements collected as part of LAPSE-RATE targeted quantities related to enhancing our understanding of boundary layer structure, cloud and aerosol properties and surface–atmosphere exchange and provide detailed information to support model evaluation and improvement work. Additionally, intensive intercomparison between the different unmanned aircraft platforms was completed. The current paper describes the observations obtained using three different types of surface-based mobile observing vehicles. These included the University of Colorado Mobile UAS Research Collaboratory (MURC), the National Oceanic and Atmospheric Administration National Severe Storms Laboratory Mobile Mesonet, and two University of Nebraska Combined Mesonet and Tracker (CoMeT) vehicles. Over the 1-week campaign, a total of 143 h of data were collected using this combination of vehicles. The data from these coordinated activities provide detailed perspectives on the spatial variability of atmospheric state parameters (air temperature, humidity, pressure, and wind) throughout the northern half of the San Luis Valley. These datasets have been checked for quality and published to the Zenodo data archive under a specific “community” setup for LAPSE-RATE (https://zenodo.org/communities/lapse-rate/, last access: 21 January 2021) and are accessible at no cost by all registered users. The primary dataset DOIs are https://doi.org/10.5281/zenodo.3814765 (CU MURC measurements; de Boer et al., 2020d), https://doi.org/10.5281/zenodo.3738175 (NSSL MM measurements; Waugh, 2020), and https://doi.org/10.5281/zenodo.3838724 (UNL CoMeT measurements; Houston and Erwin, 2020).

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

Antarctic permafrost processes and antiphase dynamics of cold-based glaciers in the McMurdo Dry Valleys inferred from 10 Be and 26 Al cosmogenic nuclides

Abstract. Soil and sediment mixing and associated permafrost processes are not widely studied or understood in the McMurdo Dry Valleys of Antarctica. In this study, we investigate the stability and depositional history of near-surface permafrost sediments to ∼ 3 m depth in the Pearse and lower Wright valleys using measured cosmogenic 10Be and 26Al depth profiles. In Pearse Valley, we estimate a minimum depositional age of ∼ 74 ka for the active layer and paleoactive-layer sediments (< 0.65 m). Combined depth profile modelling of 10Be and 26Al gives a depositional age for near-surface (< 1.65 m) permafrost in Pearse Valley of 180 +20/-40 ka, implying that the deposition of permafrost sediments predates MIS 5 advances of Taylor Glacier. Deeper permafrost sediments (> 2.09 m) in Pearse Valley are thus inferred to have a depositional age of > 180 ka. At a coastal, lower-elevation site in neighbouring lower Wright Valley, 10Be and 26Al depth profiles from a second permafrost core exhibit near-constant concentrations with depth and indicate the sediments are either vertically mixed after deposition or sufficiently young so that post-depositional nuclide production is negligible relative to inheritance. 26Al/10Be concentration ratios for both depth profiles range between 4.0 and 5.2 and are all lower than the nominal surface production rate ratio of 6.75, indicating that prior to deposition, these sediments experienced complex, yet similar, exposure–burial histories. Assuming a single-cycle exposure–burial scenario, the observed 26Al/10Be ratios are equivalent to a total minimum exposure–burial history of ∼ 1.2 Myr. In proximity to the depth profile core site, we measured cosmogenic 10Be and 26Al in three granite cobbles from thin, patchy drift (Taylor 2 Drift) in Pearse Valley to constrain the timing of retreat of Taylor Glacier. Assuming simple continuous exposure, our minimum, zero-erosion exposure ages suggest Taylor Glacier partially retreated from Pearse Valley no later than 65–74 ka. The timing of retreat after 65 ka and until the Last Glacial Maximum (LGM) when Taylor Glacier was at a minimum position remains unresolved. The surface cobble ages and permafrost processes reveal Taylor Glacier advances during MIS 5 were non-erosive or mildly erosive, preserving the underlying permafrost sediments and peppering boulders and cobbles upon an older, relict surface. Our results are consistent with U/Th ages from central Taylor Valley and suggest changes in moisture delivery over Taylor Dome during MIS 5e, 5c, and 5a appear to be associated with the extent of the Ross Ice Shelf and sea ice in the Ross Sea. These data provide further evidence of antiphase behaviour through retreat of a peripheral lobe of Taylor Glacier in Pearse Valley, a region that was glaciated during MIS 5. We suggest a causal relationship of cold-based glacier advance and retreat that is controlled by an increase in moisture availability during retreat of sea ice and perhaps the Ross Ice Shelf, as well as, conversely, a decrease during times of sea ice and Ross Ice Shelf expansion in the Ross Sea.

58 GEOSCIENCES↗
DOE OSTI2022

Large Magnetic Gap in a Designer Ferromagnet–Topological Insulator–Ferromagnet Heterostructure

Abstract Combining magnetism and nontrivial band topology gives rise to quantum anomalous Hall (QAH) insulators and exotic quantum phases such as the QAH effect where current flows without dissipation along quantized edge states. Inducing magnetic order in topological insulators via proximity to a magnetic material offers a promising pathway toward achieving the QAH effect at a high temperature for lossless transport applications. One promising architecture involves a sandwich structure comprising two single‐septuple layers (1SL) of MnBi 2 Te 4 (a 2D ferromagnetic insulator) with ultrathin few quintuple layer (QL) Bi 2 Te 3 in the middle, and it is predicted to yield a robust QAH insulator phase with a large bandgap greater than 50 meV. Here, the growth of a 1SL MnBi 2 Te 4 /4QL Bi 2 Te 3 /1SL MnBi 2 Te 4 heterostructure via molecular beam epitaxy is demonstrated and the electronic structure probed using angle‐resolved photoelectron spectroscopy. Strong hexagonally warped massive Dirac fermions and a bandgap of 75 ± 15 meV are observed. The magnetic origin of the gap is confirmed by the observation of the exchange‐Rashba effect, as well as the vanishing bandgap above the Curie temperature, in agreement with density functional theory calculations. These findings provide insights into magnetic proximity effects in topological insulators and reveal a promising platform for realizing the QAH effect at elevated temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Microstructural Passivation of Patterned Al 2 O 3 Back Contacts on CdS/CdTe Solar Cells

Abstract Patterned aluminum oxide (Al 2 O 3 ) back contact on cadmium telluride (CdTe) solar cells can effectively mitigate the loss of the majority carriers while retaining the passivation benefits for minority carriers. This study investigates the impact of the point‐contact geometry on cell performance, where the spacing of the microholes is systematically varied at a constant microhole diameter (≈16 µm). Microhole arrays are created on an evaporated Al 2 O 3 layer (≈30 nm) on CdTe using laser‐beam lithography and selective wet etching processes. Comparative analysis indicates a significant decrease in fill factor (FF) and short‐circuit current ( I SC ) when the contact fraction falls below 30%, likely due to the Al 2 O 3 barrier impeding majority carriers. Above this fraction, FF and I SC are comparable to those of baseline CdTe cells without Al 2 O 3 . Open‐circuit voltage ( V OC ) shows a gradual decrease as Al 2 O 3 coverage is reduced but exhibits a slight increase (<15 mV) when the contact fraction exceeds 30%. Cathodoluminescence (CL) characterization reveals that significantly improved radiative recombination occurs within the CdTe grain bulk with Al 2 O 3 , whereas this effect on grain boundaries is minimal. The findings imply complex local carrier dynamics in the patterned back‐contact region, closely tied to the microstructural properties of CdTe‐based solar cells.

Roy, Etee Kawna↗
DOE OSTI2020

Direct Nano‐Synthesis Methods Notably Benefit Mg‐Battery Cathode Performance

Abstract Rechargeable magnesium batteries are promising candidates for next‐generation electrochemical energy storage, but their development is severely hindered by sluggish solid‐state diffusion and significant desolvation penalties of the divalent cation. Studies suggest that nano‐sized electrode materials alleviate these issues by shortening diffusion lengths and increasing electrode/electrolyte interaction. Here, the effect of particle size and synthetic methodology on the electrochemical performance of four sulfide cathode materials in Mg batteries is investigated: layered TiS 2 , CuS, spinel Ti 2 S 4 , and CuCo 2 S 4 . In these sulfide hosts, the direct preparation of nano‐dimensional crystallites is critical to activate or improve electrochemistry. Even promising cathode materials can appear electrochemically inert when micron‐sized particles are investigated (e.g., CuCo 2 S 4 ), and mechanical milling leads to surface degradation of active material which severely limits performance. However, nano‐sized CuCo 2 S 4 prepared directly reaches a capacity nearly double that of ball‐milled material and delivers 350 mAh g −1 at 60 °C. This work provides synthetic considerations which may be crucial in the discovery and design of novel Mg cathode materials, so that promising candidates are not overlooked. By extension, in oxide materials where Mg 2+ diffusion is expected to be much more sluggish, this factor is anticipated to be even more important when screening for new hosts.

Blanc, Lauren E.↗
DOE OSTI2025

Uncertainty quantification for neural network potential foundation models

Abstract For neural network potentials (NNPs) to gain widespread use, researchers must be able to trust model outputs. However, the blackbox nature of neural networks and their inherent stochasticity are often deterrents, especially for foundation models trained over broad swaths of chemical space. Uncertainty information provided at the time of prediction can help reduce aversion to NNPs. In this work, we detail two uncertainty quantification (UQ) methods. Readout ensembling, by finetuning the readout layers of an ensemble of foundation models, provides information about model uncertainty, while quantile regression, by replacing point predictions with distributional predictions, provides information about uncertainty within the underlying training data. We demonstrate our approach with the MACE-MP-0 model, applying UQ to the foundation model and a series of finetuned models. The uncertainties produced by the readout ensemble and quantile methods are demonstrated to be distinct measures by which the quality of the NNP output can be judged.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Integral analysis of the effect of material dimension and composition on tokamak neutronics *

Abstract The neutronics performance of a tokamak has been identified as an important factor in designing a fusion power plant. The design of the tokamak should not only meet operational parameters such as sufficient tritium breeding, but also safety parameters such as low structural material activation. This paper investigates the impacts of the neutronics metrics for the ARC-class tokamak, a compact tokamak with an immersion blanket, by perturbing the first five layers of structural material—first wall, inner vacuum vessel, coolant salt channel, neutron multiplier, and outer vacuum vessel. The goal of this work is to provide insight into shaping and scaling the flux on each layer to obtain optimized operational and safety metrics through quantification of the responses from each perturbation. Results show that increased first wall thickness can increase the tritium breeding ratio (TBR) in specific configurations with high 6 Li enrichments and that vacuum vessels decrease TBR for low- 6 Li enrichment configurations. It was also found that the neutron multiplier can either increase or decrease TBR depending on the configuration. The response of metrics to the change in layer thickness and enrichment also varies depending on the vacuum vessel material. The integral impacts of 6 Li enrichment, layer thicknesses, and vacuum vessel material choice are investigated and presented in this paper.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Experimental and Computational Investigation of Integrated Internal and Film Cooling Designs Incorporating a Thermal Barrier Coating

Abstract Few studies in the open literature have studied the effect of thermal barrier coatings (TBC) when used in combination with shaped hole film cooling and enhanced internal cooling techniques. The current study presents Reynolds-averaged Navier–Stokes (RANS) conjugate heat transfer simulations that identify trends in cooling design performance as well as experimental measurements of overall effectiveness using a flat-plate matched-Biot number model with a simulated TBC layer of 0.42D thickness, where D is the film cooling hole diameter. Coolant is fed to the film cooling holes in a co-flow configuration, and the results of both the smooth and rib-turbulated channels are compared. At a constant coolant flow rate, enhanced internal cooling was found to provide a 44% increase in spatially-averaged overall effectiveness, $\bar{\bar{\phi}}$, without a TBC. The results show that the addition of a TBC can raise $\bar{\bar{\phi}}$ on a film-cooled component surface by 47%. The optimum velocity ratio was found to decrease with the addition of enhanced cooling techniques and a TBC as the film provided minimal benefit at the expense of reduced internal cooling. While the computational results closely identified trends in overall system performance without a TBC, the model over-predicted effectiveness on the metal–TBC interface. The results of this study will inform turbine component design as material science advances increase the reliability of the TBC.

Engineering↗
NASA NTRS2010

Geologic Map of the Olympia Cavi Region of Mars (MTM 85200): A Summary of Tactical Approaches

The 1:500K-scale geologic map of MTM 85200 - the Olympia Cavi region of Mars - has been submitted for peer review [1]. Physiographically, the quadrangle includes portions of Olympia Rupes, a set of sinuous scarps which elevate Planum Boreum ~800 meters above Olympia Planum. The region includes the high-standing, spiral troughs of Boreales Scopuli, the rugged and deep depressions of Olympia Cavi, and the vast dune fields of Olympia Undae. Geologically, the mapped units and landforms reflect the recent history of repeated accumulation and degradation. The widespread occurrence of both weakly and strongly stratified units implicates the drape-like accumulation of ice, dust, and sand through climatic variations. Similarly, the occurrence of layer truncations, particularly at unit boundaries, implicates punctuated periods of both localized and regional erosion and surface deflation whereby underlying units were exhumed and their material transported and re-deposited. Herein, we focus on the iterative mapping approaches that allowed not only the accommodation of the burgeoning variety and volume of data sets, but also facilitated the efficient presentation of map information. Unit characteristics and their geologic history are detailed in past abstracts [2-3].

Skinner, J. A., Jr.↗
DOE OSTI2021

Unique Features of Polarization in Ferroelectric Ionic Conductors

Abstract Ferroelectrics that are also ionic conductors offer possibilities for novel applications with high tunability, especially if the same atomic species causes both phenomena. In particular, at temperatures just below the Curie temperature, polarized states may be sustainable as the mobile species is driven in a controlled way over the energy barrier that governs ionic conduction, resulting in unique control of the polarization. This possibility is recently demonstrated in CuInP 2 S 6 , a layered ferroelectric ionic conductor in which Cu ions cause both ferroelectricity and ionic conduction. Here, it is shown that the commonly used approach to calculate the polarization of evolving atomic configurations in ferroelectrics using the modern theory of polarization, namely concerted (synchronous) migration of the displacing ions, is not well suited to describe the polarization evolution as the Cu ions cross the van der Waals gaps. An asynchronous Cu‐migration scheme is introduced, which reflects the physical process by which Cu ions migrate, resolves the difficulties, and describes the polarization evolution both for normal ferroelectric switching and for transitions across the van der Waals gaps, providing a single framework to discuss ferroelectric ionic conductors.

O'Hara, Andrew↗
DOE OSTI2024

Thermodynamic Modeling of Solid Flux Interactions with Molten Aluminum

Abstract Oxide and dross formation during aluminum melt processing results in a considerable amount of loss of metallic aluminum. The total melt loss generated during melt processing can be greatly reduced by efficient use of flux, particularly for melting aluminum scrap or secondary alloys. Effective use of cover fluxes can significantly reduce dross generation by creating a barrier between the metallic aluminum and the atmosphere. The amount of metallic aluminum trapped within the dross layer can also be reduced by up to 50% using drossing fluxes. However, reducing the industrial average melt loss has remained difficult. To evaluate flux ingredients used in the casting industry, computational thermodynamic software ThermoCalc was used to calculate the driving force for reactions between twenty-one flux ingredients and eighteen common alloying and impurity elements in foundry alloys. The thermodynamic calculations, combined with other properties, were used to provide a desirable list of cover and drossing ingredients, which are being experimentally validated.

Moodispaw, Michael P. (ORCID:0000000267635882)↗
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↗