Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Convection lines”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Urban Impacts on Convective Squall Lines over Chicago in the Warm Season—Part I: Observations of Multi-Scale Convective Evolution

In this study, our aim is to diagnose how two quasi-linear convective systems (QLCS) are organized so one can determine the possible role of the city of Chicago, IL, USA, in modifying convective precipitation systems. In this Part I of a two-part study, we employ large-scale analyses, radiosonde soundings, surface observations, and Doppler radar data to diagnose the precursor atmospheric circulations that organize the evolution of two mesoscale convective systems and compare those circulations to radar and precipitation. Several multi-scale processes are found that organize and modify convection over the Chicago metroplex. Two sequential quasi-linear convective systems (QLCS #1 and #2) were organized that propagated over Chicago, IL, USA, during an eight-hour period on 5–6 July 2018. The first squall line (QLCS #1) built from the southwest to the northeast while strengthening as it propagated over the city, and the second (QLCS #2) propagated southeastwards and weakened as it passed over the city in association with a polar cold front. The weak upper-level divergence associated with a diffluent flow poleward of an expansive ridge built over and strengthened a low-level trough and confluence zone, triggering QLCS #1. Convective downdrafts from QLCS #1 produced a cold pool that interacted with multiple confluent low-level jets surrounding and focused on the metroplex urban heat island, thus advecting the convection poleward over the metroplex. The heaviest precipitation occurred just south-southeast of Midway Airport, Chicago. Subsequently, a polar cold front propagated into the metroplex, which triggered QLCS #2. However, the descending air above it under the polar jet and residual cold pool from QLCS #1 rapidly dissipated the cold frontal convection. This represents a case study where very active convection built over the metroplex and was likely modified by it, as evidenced in numerical simulations to be described in Part II.

Kaplan, Michael L.↗
DOE OSTI2021

Climatology of Linear Mesoscale Convective System Morphology in the United States based on Random Forests Method

This study uses machine learning methods, specifically the random forest (RF), on a radar-based mesoscale convective system (MCS) tracking dataset to classify the five types of linear MCS morphology in the contiguous United States during the period 2004-2016. The algorithm is trained using radar- and satellite-derived spatial and morphological parameters, and reanalysis environmental information from 5-26yr manually identified nonlinear and five linear MCS modes. The algorithm is then used to automate the classification of linear MCSs over 8 years with high accuracy, providing a systematic, long-term climatology of linear MCSs. Results reveal that nearly 40% of MCSs are classified as linear MCSs, in which half of the linear events belong to the type of system having a leading convective line. The occurrence of linear MCSs shows large annual and seasonal variations. On average, 113 linear MCSs occur annually during the warm season (through March to October), with most of these events clustered from May through August in the central eastern Great Plains. MCS characteristics, including duration, propagation speed, orientation, and system cloud size, have large variability among the different linear modes. The systems having a trailing convective line and the systems having a back-building area of convection typically move more slowly and have higher precipitation rate, and thus have higher potential in producing extreme rainfall and flash flooding. Analysis of the environmental conditions associated with linear MCSs show that the storm-relative flow is of most importance in determining the organization mode of linear MCSs.

Cui, Wenjun↗
DOE OSTI2022

Dependencies of Four Mechanisms of Secondary Ice Production on Cloud-Top Temperature in a Continental Convective Storm

Various mechanisms of secondary ice production (SIP) cause multiplication of numbers of ice particle, after the onset of primary ice. A measure of SIP is the ice enhancement ratio (“IE ratio”) defined here as the ratio between number concentrations of total ice (excluding homogeneously nucleated ice) and active ice-nucleating particles (INPs). A convective line observed on 11 May 2011 over the Southern Great Plains in the Mesoscale Continental Convective Cloud Experiment (MC3E) campaign was simulated with the “Aerosol–Cloud” (AC) model. AC is validated against coincident MC3E observations by aircraft, ground-based instruments, and satellite. Four SIP mechanisms are represented in AC: the Hallett–Mossop (HM) process of rime splintering, and fragmentation during ice–ice collisions, raindrop freezing, and sublimation. The vertical profile of the IE ratio, averaged over the entire simulation, is almost uniform (10 2 to 10 3 ) because fragmentation in ice–ice collisions dominates at long time scales, driving the ice concentration toward a theoretical maximum. The IE ratio increases with both the updraft (HM process, fragmentation during raindrop freezing, and ice–ice collisions) and downdraft speed (fragmentation during ice–ice collisions and sublimation). As reported historically in aircraft sampling, IE ratios were predicted to peak near 10 3 for cloud-top temperatures close to the −12°C level, mostly due to the HM process in typically young clouds with their age less than 15 min. Here, at higher altitudes with temperatures of −20° to −30°C, the predicted IE ratios were smaller, ranging from 10 to 10 2 , and mainly resulted from fragmentation in ice–ice collisions.

Aerosols↗
DOE OSTI2020

Impacts of Stochastic Mixing in Idealized Convection-Permitting Simulations of Squall Lines

This study investigates impacts of altering subgrid-scale mixing in “convection-permitting” kilometer-scale horizontal-grid-spacing (Δ h ) simulations by applying either constant or stochastic multiplicative factors to the horizontal mixing coefficients within the Weather Research and Forecasting Model. In quasi-idealized 1-km Δ h simulations of two observationally based squall-line cases, constant enhanced mixing produces larger updraft cores that are more dilute at upper levels, weakens the cold pool, rear-inflow jet, and front-to-rear flow of the squall line, and degrades the model’s effective resolution. Reducing mixing by a constant multiplicative factor has the opposite effect on all metrics. Completely turning off parameterized horizontal mixing produces bulk updraft statistics and squall-line mesoscale structure closest to an LES “benchmark” among all 1-km simulations, although the updraft cores are too undilute. The stochastic mixing scheme, which applies a multiplicative factor to the mixing coefficients that varies stochastically in time and space, is employed at 0.5-, 1-, and 2-km Δ h . It generally reduces midlevel vertical velocities and enhances upper-level vertical velocities compared to simulations using the standard mixing scheme, with more substantial impacts at 1- and 2-km Δ h compared to 0.5-km Δ h . Further, the stochastic scheme also increases updraft dilution to better agree with the LES for one case, but has less impact on the other case. Stochastic mixing acts to weaken the cold pool but without a significant impact on squall-line propagation. It also does not affect the model’s overall effective resolution unlike applying constant multiplicative factors to the mixing coefficients.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Evaluating the Conservation of Energy Variables in Simulations of Deep Moist Convection

Abstract It is often assumed in parcel theory calculations, numerical models, and cumulus parameterizations that moist static energy (MSE) is adiabatically conserved. However, the adiabatic conservation of MSE is only approximate because of the assumption of hydrostatic balance. Two alternative variables are evaluated here: MSE − IB and MSE + KE, wherein IB is the path integral of buoyancy ( B ) and KE is kinetic energy. Both of these variables relax the hydrostatic assumption and are more precisely conserved than MSE. This article quantifies the errors that result from assuming that the aforementioned variables are conserved in large-eddy simulations (LES) of both disorganized and organized deep convection. Results show that both MSE − IB and MSE + KE better predict quantities along trajectories than MSE alone. MSE − IB is better conserved in isolated deep convection, whereas MSE − IB and MSE + KE perform comparably in squall-line simulations. These results are explained by differences between the pressure perturbation behavior of squall lines and isolated convection. Errors in updraft B diagnoses are universally minimized when MSE − IB is assumed to be adiabatically conserved, but only when moisture dependencies of heat capacity and temperature dependency of latent heating are accounted for. When less accurate latent heat and heat capacity formulae were used, MSE − IB yielded poorer B predictions than MSE due to compensating errors. Our results suggest that various applications would benefit from using either MSE − IB or MSE + KE instead of MSE with properly formulated heat capacities and latent heats.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Tracing the Rain Formation Pathways in Numerical Simulations of Deep Convection

Quantifying the microphysical process contributions to surface precipitation in numerical simulations can be challenging. This is due to the fact that many microphysical processes contribute to the formation and depletion of rain drops and there is almost always a spatial/temporal mismatch between where/when rain is formed and where/when it strikes the surface. In this work, we develop a tracing method that tracks the sources and sinks of raindrop mass and number as they are advected by the Weather Research and Forecasting model. Applying the method to an idealized squall line confirms that convective precipitation is dominated by warm rain processes (autoconversion and accretion) while stratiform precipitation is dominated by the melting of rimed and unrimed ice crystals. Sensitivity experiments in which the prescribed cloud drop number concentration is increased confirm the conventional wisdom that weakened autoconversion increases the fraction of raindrops originating from cold rain processes. The method also reveals that when applied to deep convection the Khairoutdinov and Kogan autoconversion scheme produces an excessive number of raindrops which are subsequently clipped in P3 microphysics to keep the rain size distribution within prescribed limits. This problem can mostly be mitigated by increasing the assumed radius for raindrops created by autoconversion.

58 GEOSCIENCES↗
DOE OSTI2022

A measurement of stellar surface gravity hidden in radial velocity differences of comoving stars

The gravitational redshift induced by stellar surface gravity is notoriously difficult to measure for non-degenerate stars, since its amplitude is small in comparison with the typical Doppler shift induced by stellar radial velocity. In this study, we make use of the large observational data set of the Gaia mission to achieve a significant reduction of noise caused by these random stellar motions. By measuring the differences in velocities between the components of the pairs of comoving stars and wide binaries, in this work we are able to statistically measure the combined effects of gravitational redshift and convective blueshifting of spectral lines, and nullify the effect of the peculiar motions of the stars. For the subset of stars considered in this study, we find a positive correlation between the observed differences in Gaia radial velocities and the differences in surface gravity and convective blueshift inferred from effective temperature and luminosity measurements. The results rule out a null signal at the 5σ level for our full data set. Additionally, we study the subdominant effects of binary motion, and possible systematic errors in radial velocity measurements within Gaia. Results from the technique presented in this study are expected to improve significantly with data from the next Gaia data release. Such improvements could be used to constrain the mass–luminosity relation and stellar models that predict the magnitude of convective blueshift.

(stars:) binaries: general↗
DOE OSTI2023

Observations of Boundary Layer Convergence Lines and Associated Updrafts in the U.S. Southern Great Plains

Abstract Boundary layer convergence lines (CLs) are highly effective at deep-convection initiation (DCI), suggesting that their associated updraft properties differ from those of more widespread turbulent updrafts in the planetary boundary layer (PBL). This study exploits observations at the Atmospheric Radiation Measurement Southern Great Plains (ARM SGP) observatory in Oklahoma from 2011 to 2016 to quantify CL properties and their relation to turbulent PBL eddies preceding CL arrival. Two independent methods for estimating CL properties are developed at two locations in the SGP region, both relying on the assumption of a 2D circulation in the CL-normal plane but using different combinations of instruments. The first (the radar method) relies mainly on scanning radar data and is applied to 61 CLs passing near a high-resolution scanning radar based in Nardin, Oklahoma, while the second (the surface method) relies mainly on surface wind data and is applied to 68 CLs crossing the SGP facility in nearby Lamont, Oklahoma. Mean daytime (1000–1900 LST) CL width (∼2 km) and convergence magnitude (∼0.003 s −1 ) are similar for both methods, and mean daytime CL depth is ∼0.75 km. The two methods disagree at night (0000–1000 and 1900–2400 LST), where the surface method estimates wider and weaker CLs than the radar method. This difference may stem from the radar beam overshooting the shallow, highly stable nocturnal PBL. The largest CL updrafts are slightly wider (∼20%) and stronger (∼40%) than the largest PBL updrafts in the pre-CL period, generating 50%–100% larger updraft mass fluxes over most of the PBL depth. Significance Statement Deep convection is commonly initiated by boundary layer convergence lines (CLs), which are associated with intense surface-based wind convergence and strong updrafts that may lift air to saturation. Although CLs form regularly, they are far less common than ordinary, short-lived turbulent thermals in the daytime boundary layer. To better understand why CLs are so effective at deep-convection initiation, we observationally quantify their morphologies and strengths and compare these properties to those of surrounding turbulent updrafts. Perhaps surprisingly, the CLs are found to exhibit only slightly larger scales and strengths as the turbulent updrafts. Although these marginal increases help to explain the preference for storms to initiate along CLs, they likely are not the whole story.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Hydrometeor Transport and the Stratiform Water Budget in Amazonian Mesoscale Convective Systems

A significant portion of precipitation falls within the stratiform regions of mesoscale convective systems (MCSs). These regions are formed from convective clouds that have either weakened and formed a horizontally stratified structure or have advected hydrometeors to form a trailing-stratiform region, with further precipitation growth promoted by in situ dynamical and microphysical processes. Green Ocean Amazon 2014/5 (GoAmazon2014/5) field campaign observations are used to quantify the stratiform water budget of six MCSs (three each from the wet and dry seasons) utilizing the framework from Leary and Houze, with an emphasis on the transfer of hydrometeors (C A ) from the convective to the stratiform regions. The storms were predominantly leading-line/trailing-stratiform MCSs, but two had embedded convection in the stratiform region, which can complicate water budget calculations. A pulse-like pattern with a mean interval of 21 min occurs in C A during each storm resulting from a lag between the initial convective updrafts and the subsequent hydrometeor transport and sedimentation in the stratiform region. Here, the peaks in C A were followed by peaks in hydrometeor growth associated with the mesoscale updraft (C mu ) in the stratiform region. Averaged over the six storms, the positive contributions of C A and C mu to the stratiform water budget were roughly equal, although storm-to-storm variations existed. Negative contributions from evaporation below cloud base in the mesoscale downdraft and detrainment from the stratiform region into the nonprecipitating anvil cloud and environment accounted for 53% and 21% of the sum of C A and C mu , respectively, resulting in a mean stratiform region precipitation efficiency of 26%.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

How does vertical wind shear influence entrainment in squall lines?

The influence of vertical wind shear on updraft entrainment in squall lines is not well understood. To address this knowledge gap, a suite of high-resolution idealized numerical model simulations of squall lines were run in various vertical wind shear (hereafter “shear”) environments to study the effects of shear on entrainment in deep convective updrafts. Low-level horizontal mass flux into the leading edge of the cold pool was strongest in the simulations with the strongest low-level shear. These simulations consequently displayed wider updrafts, less entrainment-driven dilution, and larger buoyancy than the simulations with comparatively weak low-level shear. An analysis of vertical accelerations along trajectories that passed through updrafts showed larger net accelerations from buoyancy in the simulations with stronger low-level shear, which demonstrates how less entrainment-driven dilution equated to stronger updrafts. Additionally, the effects of upper-level shear on entrainment and updraft vertical velocities were generally less pronounced than the effects of low-level shear. We argue that in addition to the outflow boundary-shear interactions and their effect on updraft tilt established by previous authors, decreased entrainment-driven dilution is yet another beneficial effect of strong low-level shear on squall line updraft intensity.

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

Effects of Surface Turbulence Flux Parameterizations on the MJO: The Role of Ocean Surface Waves

This study investigates the sensitivity of the Madden–Julian oscillation (MJO) to changes to the bulk flux parameterization and the role of ocean surface waves in air–sea coupling using a fully coupled ocean–atmosphere–wave model. The atmospheric and ocean model components of the Energy Exascale Earth System Model (E3SM) are coupled to a spectral wave model, WAVEWATCH III (WW3). Two experiments with wind speed–dependent bulk algorithms (NCAR and COARE3.0a) and one experiment with wave-state-dependent flux (COR3.0a-WAV) were conducted. We modify COARE3.0a to include surface roughness calculated within WW3 and also account for the buffering effect of waves on the relative difference between air-side and ocean-side momentum flux. Differences in surface fluxes, primarily caused by discrepancies in drag coefficients, result in significant differences in MJO’s properties. While COARE3.0a has better convection–circulation coupling than NCAR, it exhibits anomalous MJO convection east of the date line. The wave-state-dependent flux (COR3.0-WAV) improves the MJO representation over the default COARE3.0 algorithm. Strong easterlies over the Pacific Ocean in COARE3.0a enhance the latent heat flux (LHFLX). This is responsible for the anomalous MJO propagation after the date line. In COR3.0a-WAV, waves reduce the anomalous easterlies, leading to a decrease in LHFLX and MJO dissipation after the date line. These findings highlight the role of surface fluxes in MJO simulation fidelity. Most importantly, we show that the proper treatment of wave-induced effects in bulk flux parameterization improves the simulation of coupled climate variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Hourly dynamic line ratings for existing transmission across the contiguous United States (preliminary results) [Slides]

The thermal rating of a transmission line depends on the rate at which heat generated by current flow through the line can be dissipated into the environment by convective and radiative cooling, keeping the conductor below a maximum operating temperature. Static or seasonal line ratings assume nearly-worst-case weather conditions (high ambient air temperature, high solar irradiance, and low wind speed), resulting in a low risk of "overrating" a line (i.e., allowing more current to flow than the line is designed for) but a high fraction of time when the line is "underrated" (i.e., able to flow more current while staying below the maximum temperature). Dynamic transmission line ratings (DLR) use real-time measurements of weather conditions, line sag, or line temperature to determine the thermal rating of transmission lines at hourly or higher time resolution, allowing for increased ratings during some (cold, windy, and/or dark) periods and correcting for overrating during other (hot, calm, and/or sunny) periods. Many studies have explored the benefits of DLR, but typically for small collections of transmission lines under a limited sample of weather conditions, leaving open the question of how widely applicable the benefits of DLR might be over diverse geographies and weather conditions. Here, we combine geospatial data on tens of thousands of transmission line routes across the contiguous United States with high-resolution hourly historical weather data for 2007-2013 to estimate the impact of DLR and other weather-adjusted rating systems. We compare static and dynamic line ratings to intermediate ratings adjusted by ambient temperature, day/night irradiance, clear-sky irradiance, and/or measured irradiance, illustrating the benefits of incorporating progressively more detailed weather data into the calculated line ratings. Moving from nighttime-adjusted ratings to clear-sky-adjusted ratings provides more benefit than moving from clear-sky-adjusted ratings to measured-irradiance-adjusted ratings. Depending on the assumed static windspeed for irradiance-adjusted-ratings, DLR either corrects for overrated hours or provides substantially increased ratings in the majority of hours on the majority of lines.

14 SOLAR ENERGY↗
DOE OSTI2021

Scaling theory of three-dimensional magnetic reconnection spreading

We develop a first-principles scaling theory of the spreading of three-dimensional (3D) magnetic reconnection of finite extent in the out of plane direction. This theory addresses systems with or without an out of plane (guide) magnetic field, and with or without Hall physics. The theory reproduces known spreading speeds and directions with and without guide fields, unifying previous knowledge in a single theory. New results include: (1) Reconnection spreads in a particular direction if an x-line is induced at the interface between reconnecting and non-reconnecting regions, which is controlled by the out of plane gradient of the electric field in the outflow direction. (2) The spreading mechanism for anti-parallel collisionless reconnection is convection, as is known, but for guide field reconnection it is magnetic field bending. We confirm the theory using 3D two-fluid and resistive-magnetohydrodynamics simulations. (3) The theory explains why anti-parallel reconnection in resistive-magnetohydrodynamics does not spread. (4) The simulation domain aspect ratio, associated with the free magnetic energy, influences whether reconnection spreads or convects with a fixed x-line length. (5) We perform a simulation initiating anti-parallel collisionless reconnection with a pressure pulse instead of a magnetic perturbation, finding spreading is unchanged rather than spreading at the magnetosonic speed as previously suggested. The results provide a theoretical framework for understanding spreading beyond systems studied here, and are important for applications including two-ribbon solar flares and reconnection in Earth’s magnetosphere.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The impacts of biomass burning activities on convective systems over the Maritime Continent

Convective precipitation associated with Sumatra squall lines and diurnal rainfall over Borneo is an important weather feature of the Maritime Continent in Southeast Asia. Over the past few decades, biomass burning activities have been widespread during summertime over this region, producing massive fire aerosols. These additional aerosols, when brought into the atmosphere, besides influencing the local radiation budget through directly scattering and absorbing sunlight, can also act as cloud condensation nuclei or ice nuclei to alter convective clouds and precipitation over the Maritime Continent via so-called aerosol indirect effects. Based on 4-month simulations with or without biomass burning aerosols, conducted using theWeather Research and Forecasting model coupled with a chemistry module (WRFChem), we have investigated the aerosol–cloud interactions associated with biomass burning aerosols over the Maritime Continent. Results from selected cases of convective events have specifically shown the significant impact of fire aerosols on weak convections by their increasing of the quantities of hydrometeors and rainfall in both the Sumatra and Borneo regions. Statistical analysis over the fire season also suggests that fire aerosols have impacts on the nocturnal convections associated with the local anticyclonic circulation in western Borneo and weaken nocturnal rainfall intensity by about 9 %. Such an effect is likely to have come from the near-surface heating due to absorbing aerosols emitted from fires, which could weaken land breezes and thus the convergence of anticyclonic circulation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Heat transfer and pseudo phase transition for low-Reynolds, mixed-convection channel flow in the supercritical thermodynamic regime

Fluids at supercritical thermodynamic conditions are inherently complex due to the large variations in thermodynamic and transport properties. Multiple recent numerical and experimental investigations indicate ongoing interest for these fluids, especially supercritical CO and supercritical water, for a variety of applications. The potential heat transfer benefits in this regime are directly influenced by the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. Limited studies have been conducted for low-Reynolds/intermediate-Rayleigh numbers mixed-convection channel flow of supercritical water, which is the focus of this study. To investigate the thermally driven hydrodynamic instabilities in this regime, we use a high-order fully-implicit numerical method. As expected, in this channel configuration and flow regime, a competition arises between the forced and the natural convection. Buoyancy, driven by the thermal gradient at the wall, pushes the fluid upwards in the channel, while inertial forces drive the fluid from inlet to outlet. We observe significant differences in the heat transfer coefficient between the top and bottom walls of the channel flow, due to natural convection, and unstable plumes. We further quantify the buoyancy forces along the channel flow through evaluation of the Richardson number. Simulations show that with a higher negative Richardson number, there are more unstable structures along the channel, resulting in enhanced mixing. Furthermore, we demonstrate the existence of gas-like, liquid-like, and the possibility of pseudo-two-phase regimes in the supercritical fluid. The pseudo-two-phase regime occurs when the range of fluid temperatures is crossing the pseudo-critical line, also known as the Widom line. The density of the fluid varies significantly near the Widom line, and further influences the natural convection and heat transfer in the flow.

42 ENGINEERING↗
DOE OSTI2020

Assessing CLUBB PDF Closure Assumptions for a Continental Shallow-to-Deep Convective Transition Case Over Multiple Spatial Scales

Assumed-PDF (probability density function) higher-order turbulence closures (APHOCs) are now widely used for parameterizing boundary layer turbulence and shallow convection in Earth system models (ESMs). A better understanding of the resolution-dependent behavior of APHOCs is essential for improving the performance of next-generation ESMs with intended horizontal resolutions finer than 10 km. In this study, we evaluate the PDF family of Analytic double-Gaussian 1 implemented in Cloud Layers Unified By Binormals (CLUBB) over a range of spatial scales (Dx) from 2 to 100 km. A 120-km-wide large eddy simulation (LES) for a continental convection case during 2016 Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field campaign serves as benchmark to evaluate the PDF closure using an off-line approach. We find during the shallow convection period, the CLUBB PDF closure tends to produce positive biases of cloud properties and liquid water flux near cloud base for all scales of analysis. It produces negative biases for these variables near cloud top that are more severe for Dx larger than 25 km. Results show that replacing the CLUBB-parameterized moisture and temperature skewnesses with LES-derived ones can fix most of the biases if clipping of input moments is allowed to prevent the occurrence of unrealizable solutions. Overall, the performance of the PDF closure is better for smaller Dx = 2–5 km than for larger Dx = 50–100 km; for a given grid spacing, it is better when the convective clouds become deeper in the late afternoon. Likely causes for the resolution dependence and implications for improving the PDF closure are discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The DPG Method for the Convection-Reaction Problem, Revisited

In this work, we study both conforming and non-conforming versions of the practical DPG method for the convection-reaction problem. We determine that the most common approach for DPG stability analysis – construction of a local Fortin operator – is infeasible for the convection-reaction problem. We then develop a line of argument based on a direct proof of discrete stability; we find that employing a polynomial enrichment for the test space does not suffice for this purpose, motivating the introduction of a (two-element) subgrid mesh. The argument combines mathematical analysis with numerical experiments.

97 MATHEMATICS AND COMPUTING↗