Grid Reconstruction Research Articles

Thermal effect of the Southeast Asian low-latitude highlands on interannual variability in the date of the Bay of Bengal summer monsoon onset

The timing of the Bay of Bengal summer monsoon (BOBSM) onset has implications for the evolution of the Asian summer monsoon and associated precipitation. This study employs Japanese 55-year Reanalysis and NOAA's Gridded Precipitation Reconstruction over Land data to explore the spatiotemporal variation in the thermal influence of the Southeast Asian low-latitude highlands (SEALLH) and its effect on the BOBSM onset date. There is a significant correlation between pre-monsoon season (February–April) shortwave radiative heating (SWH) in the SEALLH and BOBSM onset in the interannual time scale. During the enhanced SWH, an anomalous vertical circulation, which converges and ascends in the lower troposphere over the SEALLH and converges and descends in the upper troposphere north of the Bay of Bengal (BOB), develops over the SEALLH–north of BOB during the pre-monsoon season with anticlockwise rotation. The warmer air mass resulting from the adiabatic heating is accumulated around the north of BOB by the insulation effect related to the anomalous vertical circulation. By facilitating the reversal of the land-sea thermal contrast, the warmer upper troposphere over the north of the BOB causes a 16-day earlier onset of the BOBSM than the case of weakened SWH. The findings shed light on subsequent research into the thermal effects of the SEALLH on Asian summer monsoon variability.

Compressed sensing reconstruction for high-SNR, rapid dissolved 129Xe gas exchange MRI.

Three-dimensional hyperpolarized 129Xe gas exchange imaging suffers from low SNR and long breath-holds, which could be improved using compressed sensing (CS). The purpose of this work was to assess whether gas exchange ratio maps are quantitatively preserved in CS-accelerated dissolved-phase 129Xe imaging and to investigate the feasibility of CS-dissolved 129Xe imaging with reduced-cost natural abundance (NA) xenon. 129Xe gas exchange imaging was performed at 1.5 T with a multi-echo spectroscopic imaging sequence. A CS reconstruction with an acceleration factor of 2 was compared retrospectively with conventional gridding reconstruction in a cohort of 16 healthy volunteers, 5 chronic obstructive pulmonary disease patients, and 23 patients who were hospitalized following COVID-19 infection. Metrics of comparison included normalized mean absolute error, mean gas exchange ratio, and red blood cell (RBC) image SNR. Dissolved 129Xe CS imaging with NA xenon was assessed in 4 healthy volunteers. CS reconstruction enabled acquisition time to be halved, and it reduced background noise. Median RBC SNR increased from 6 (2-18) to 11 (2-100) with CS, and there was strong agreement between CS and gridding mean ratio map values (R2 = 0.99). Image fidelity was maintained for gridding RBC SNR > 5, but below this, normalized mean absolute error increased nonlinearly with decreasing SNR. CS increased the mean SNR of NA 129Xe images 3-fold. CS reconstruction of dissolved 129Xe imaging improved image quality with decreased scan time, while preserving key gas exchange metrics. This will benefit patients with breathlessness and/or low gas transfer and shows promise for NA-dissolved 129Xe imaging.

3D textured mesh scene data fusion method considering multiple scales and resolutions

ABSTRACT The fusion of 3D textured mesh scene data can effectively and seamlessly integrate multiple spatially related 3D geographical scene models into larger and more comprehensive datasets. Existing methods usually have different degrees of texture distortion and unreduced texture problems in the splicing area, and these methods do not consider the current situation of Level of Detail(LoD) data organisation in the 3D textured mesh model. We propose a 3D textured mesh scene data fusion method considering multiple scales and resolutions. First, seamless extraction of splicing areas was achieved based on spatial relationship constraints and orthogonal inverse projection principles. Second, the sparse scene structure concept was employed to generate sparse summary scene images. Finally, a dual simplification operator for both mesh and texture was proposed to construct multi-resolution 3D textured mesh models, simultaneously reconstructing triangular facets and textures within the splicing regions. We used multiple 3D textured mesh data constructed from aerial obliques to conduct data fusion experimental verification. The results indicate that our method exhibits minimal texture distortion during grid reconstruction in splicing areas and model simplification, effectively reducing the amount of simplified texture data. Additionally, our method adeptly addresses multilevel detailed 3D textured mesh models, demonstrating good feasibility and advancement.

A 561-yr (1461-2022 CE) summer temperature reconstruction for Mid-Atlantic-Northeast USA shows connections to volcanic forcing and atmospheric circulation

Contextualizing current increases in Northern Hemisphere temperatures is precluded by the short instrumental record of the past ca. 120 years and the dearth of temperature-sensitive proxy records, particularly at lower latitudes south of <50 °N. We develop a network of 29 blue intensity chronologies derived from tree rings of Tsuga canadensis (L.) Carrière and Picea rubens Sarg. trees distributed across the Mid-Atlantic and Northeast USA (MANE)—a region underrepresented by multi-centennial temperature records. We use this network to reconstruct mean March-September air temperatures back to 1461 CE based on a model that explains 62% of the instrumental temperature variance from 1901−1976 CE. Since 1998 CE, MANE summer temperatures are consistently the warmest within the context of the past 561 years exceeding the 1951−1980 mean of +1.3 °C. Cool summers across MANE were frequently volcanically forced, with significant (p<0.05) temperature departures associated with 80% of the largest tropical (n=13) and extratropical (n=15) eruptions since 1461 CE. Yet, we find that more of the identified cool events in the record were likely unforced by volcanism and either related to stochastic variability or atmospheric circulation via significant associations (p<0.05) to regional, coastal sea-surface temperatures, 500-hpa geopotential height, and 300-hpa meridional and zonal wind vectors. Expanding the MANE network to the west and south and combining it with existing temperature-sensitive proxies across North America is an important next step toward producing a gridded temperature reconstruction field for North America.

Three-dimensional radial echo-planar spectroscopic imaging for hyperpolarized 13C MRSI in vivo.

To demonstrate the feasibility of 3D echo-planar spectroscopic imaging (EPSI) technique with rapid volumetric radial k-space sampling for hyperpolarized (HP) 13C magnetic resonance spectroscopic imaging (MRSI) in vivo. A radial EPSI (rEPSI) was implemented on a 3 T clinical PET/MR system. To enable volumetric coverage, the sinusoidal shaped readout gradients per k-t-spoke were rotated along the three spatial dimensions in a golden-angle like manner. A distance-weighted, density-compensated gridding reconstruction was used, also in cases with undersampling of spokes in k-space. Measurements without and with HP 13C-labeled substances were performed in phantoms and rats using a double-resonant 13C/1H volume resonator with 72 mm inner diameter. Phantom measurements demonstrated the feasibility of the implemented rEPSI sequence, as well as the robustness to undersampling in k-space up to a factor of 5 without advanced reconstruction techniques. Applied to measurements with HP [1-13C]pyruvate in a tumor-bearing rat, we obtained well-resolved MRSI datasets with a large matrix size of 123 voxels covering the whole imaging FOV of (180 mm)3 within 6.3 s, enabling to observe metabolism in dynamic acquisitions. After further optimization, the proposed rEPSI method may be useful in applications of HP 13C-tracers where unknown or varying metabolite resonances are expected, and the acquisition of dynamic, volumetric MRSI datasets with an adequate temporal resolution is a challenge.

Smart grid data compression and reconstruction by wavelet packet transform

A smart grid is a power network from generation to consumers and provides beneficial, steady, and safe electricity. It utilizes smart meters for billing, phasor measurement units to check the system's health., etc. As a result, it contains enormous volumes of real-time data that may be shared and stored by users, control centers, and services in a smart grid. It weakens the smart grid's communication networks. The size of the data in a smart grid will grow extremely in the future. As a result, it must reduce distortion in data compression and denoise while minimizing the demand on storage and communication networks. The goal of data compression and denoising should be to maximally conserve the useful data while accurately reflecting the state of the system and providing sufficient data regeneration at the receiving end. This paper has used lower-order different wavelets to represent a design to compress and reconstruct data at level three using wavelet Packet Transform. It works on the phasor measurement unit's current magnitude and voltage sag signals.•The proposed design has a better compression ratio.•Low reconstruction error.•This design is easy to access, systematic, profitable, and not time-consuming.

A decay-modeled compressed sensing reconstruction approach for non-Cartesian hyperpolarized 129Xe MRI.

Hyperpolarized 129Xe MRI benefits from non-Cartesian acquisitions that sample k-space efficiently and rapidly. However, their reconstructions are complex and burdened by decay processes unique to hyperpolarized gas. Currently used gridded reconstructions are prone to artifacts caused by magnetization decay and are ill-suited for undersampling. We present a compressed sensing (CS) reconstruction approach that incorporates magnetization decay in the forward model, thereby producing images with increased sharpness and contrast, even in undersampled data. Radio-frequency, T1, and decay processes were incorporated into the forward model and solved using iterative methods including CS. The decay-modeled reconstruction was validated in simulations and then tested in 2D/3D-spiral ventilation and 3D-radial gas-exchange MRI. Quantitative metrics including apparent-SNR and sharpness were compared between gridded, CS, and twofold undersampled CS reconstructions. Observations were validated in gas-exchange data collected from 15 healthy and 25 post-hematopoietic-stem-cell-transplant participants. CS reconstructions in simulations yielded images with threefold increases in accuracy. CS increased sharpness and contrast for ventilation in vivo imaging and showed greater accuracy for undersampled acquisitions. CS improved gas-exchange imaging, particularly in the dissolved-phase where apparent-SNR improved, and structure was made discernable. Finally, CS showed repeatability in important global gas-exchange metrics including median dissolved-gas signal ratio and median angle between real/imaginary components. A non-Cartesian CS reconstruction approach that incorporates hyperpolarized 129Xe decay processes is presented. This approach enables improved image sharpness, contrast, and overall image quality in addition to up-to threefold undersampling. This contribution benefits all hyperpolarized gas MRI through improved accuracy and decreased scan durations.

Reinterpretation of Marine Magnetic Anomalies in the Bay of Bengal: Implications for Plate Tectonics During the Mid‐Cretaceous

AbstractThe breakup of India from Antarctica in the Early Cretaceous and subsequent seafloor spreading in the NW‐SE direction initiated the formation of the Bay of Bengal and its conjugate Enderby Basin. Two hypotheses are proposed to explain the evolutionary tectonics between India and Antarctica. One involves a single breakup, while the second invokes another breakup episode. Further, a major reorganization of the plate boundaries in the mid‐Cretaceous caused the spreading center to align E‐W by Late Cretaceous. The tectonics involved in this spreading direction change is unknown and difficult to understand as it occurred in the mid‐Cretaceous, coincident with the Cretaceous Normal Superchron. In the present study, we have generated a revised magnetic isochron map of the Bay of Bengal and used the modulus of its analytic signal to understand the evolutionary tectonics between the two plates, India and Antarctica. Mesozoic magnetic anomalies M12n to M0, and mid‐Cretaceous incursion Q2 (108 Ma) are identified in the western basin of the Bay of Bengal, while Q2, Q1 (92 Ma) and Late Cretaceous isochron A34 are inferred in its eastern basin. Gravity grid reconstructions demonstrate asymmetrical crustal accretion in the conjugate corridors. Fan‐shaped spreading and a southward ridge jump between 108 and 92 Ma in the eastern basin, aided the ∼40° clockwise rotation of the spreading center. The presence of volcanic basement beneath the additional crust suggests on‐axis mode of formation for the Ninetyeast Ridge by 100 Ma. The present study favors the two‐phase breakup model for the evolution of these plates.

High-resolution reconstruction of April–September precipitation and major extreme droughts in China over the past ∼530 years

Extreme drought events have increased, causing serious losses and damage to the social economy under current warming conditions. However, short-term meteorological data limit our understanding and projection of these extremes. With the accumulation of proxy data, especially tree-ring data, large-scale precipitation field reconstruction has provided opportunities to explore underlying mechanisms further. Using point-by-point regression, we reconstructed the April–September precipitation field in China for the past ∼530 years on the basis of 590 proxy records, including 470 tree-ring width chronologies and 120 drought/flood indices. Our regression models explained average 50% of the variance in precipitation. In the statistical test on calibration and verification, our models passed the significance level that assured reconstruction quality. The reconstruction data performed well, showing consistency and better quality than previously reported reconstructions. The first three leading modes of variability in the reconstruction revealed the main distribution modes of precipitation over China. Wet/drought and extremely wet/drought years accounted for 12.81%/10.92% (68 years/58 years) and 1.69%/3.20% (9 years/17 years) of the past ∼530 years in China, respectively. Major extreme drought events can be identified explicitly in our reconstruction. The detailed features of the Chongzhen Great Drought (1637–1643), the Wanli Great Drought (1585–1590), and the Ding-Wu Great Famine (1874–1879), indicated the existence of potentially different underlying mechanisms that need further exploration. Although further improvements can be made for remote uninhabited areas and large deserts, our gridded reconstruction of April–September precipitation in China over the past ∼530 years can provide a solid database for studies on the attribution of climate change and the mechanism of extreme drought events.

Spatial utilization of historical topographic map and its application in land reconstruction of ancient Chinese urban land use

The historical topographic map preserves rich geographic information and can provide direct assistance for the reconstruction of various geographic elements. Based on the historical data of cities throughout the Qing Dynasty, the land use scale data of cities across the country was obtained using GIS and urban perimeter conversion models. This study combines city information and city circumference records from the historical maps and archives of the late Qing Dynasty to quantitatively reconstruct the use patterns of ancient China’s urban land at a spatial resolution of 1° × 1°. Uncertainty analysis of the reconstruction results was conducted using modern remote sensing image data as the validation data set. The results showed the following. (1) During the late Qing Dynasty, the total area of urban land in the various provinces and regions was 1456.015 km2. The maximum value was 208.691 km2 in Beijing-Tianjin-Hebei, the minimum value was 1.713 km2 in Qinghai, and the average value was 56.001 km2. (2) The results of grid reconstruction show that among the 398 grids with urban land distribution, the maximum value is 64.099 km2/grid, the minimum value is 0.013 km2/grid, and the average value is 3.658 km2/grid. (3) Of all the grids with urban land, the urban land grid to the west of the Hu Line accounts for 12.5% and the east to 87.5%. (4) During the late Qing Dynasty, urban land use in China was primarily concentrated in agriculturally developed areas such as the North China Plain, the Central Plains, Jiangnan, and the Sichuan-Chongqing region. (6) The results of a kernel density estimation showed that there were obviously three core areas of urban land agglomeration in China during the late Qing Dynasty: the North China Plain-Central Plains, the Jiangsu-Shanghai-Zhejiang-Anhui area, and the Sichuan-Chongqing urban core area. This study provides basic data for urban land use during historical periods and provides a basis for the quantitative reconstruction of relevant urban land data for historical archives.

A past and present perspective on the European summer vapor pressure deficit

Abstract. The response of evapotranspiration to anthropogenic warming is of critical importance for the water and carbon cycle. Contradictory conclusions about evapotranspiration changes are caused primarily by their brevity in time and sparsity in space, as well as the strong influence of internal variability. Here, we present the first gridded reconstruction of the summer (June, July, and August) vapor pressure deficit (VPD) for the past 4 centuries at the European level. This gridded reconstruction is based on 26 European tree ring oxygen isotope records and is obtained using a random forest approach. According to validation scores obtained with the Nash–Sutcliffe model efficiency, our reconstruction is robust over large parts of Europe since 1600, in particular for the westernmost and northernmost regions, where most tree ring records are located. Based on our reconstruction, we show that from the mid-1700s a trend towards higher summer VPD occurred in central Europe and the Mediterranean region that is related to a simultaneous increase in temperature and decrease in precipitation. This increasing summer VPD trend continues throughout the observational period and in recent times. Moreover, our summer VPD reconstruction helps to visualize the local and regional impacts of the current climate change, as well as to minimize statistical uncertainties of historical VPD variability. This paper provides also new insights into the relationship between summer VPD and large-scale atmospheric circulation, and we show that summer VPD has two preferred modes of variability, namely a NW–SE dipole-like mode and a N–S dipole-like mode. Furthermore, the interdisciplinary use of the data should be emphasized, as summer VPD is a crucial parameter for many climatological feedback processes in the Earth's surface system. The reconstructed summer VPD gridded data over the last 400 years are available at the following link: https://doi.org/10.5281/zenodo.5958836 (Balting et al., 2022).

Reconfiguration of low-voltage distributed power sources within electric power's distribution network based on improved particle swarm-fish swarm fusibility algorithm

With the development of distributed power sources in the distribution network, the algorithm of distribution network reconfiguration is gaining attention from experts and scholars. Its goal is to reduce the power loss during power transmission, so as to reduce the power grid loss during power transmission. And weaken the electric heating effect in the process of electric energy transmission, thus maintaining the safety of the surrounding residents. Due to the wire impedance effect, a lot of electric energy of the circuit is lost to electric heating, which is easy to cause local overheating and lead to fire. This will not only cause power loss, but also endanger the safety of surrounding residents. To address the issue, experiments on distribution grid reconstruction are performed using the enhanced particle swarm-fish swarm algorithm with the Elecgrid self-constructed dataset. Initially, low-voltage distributed power sources in parallel are connected to the circuit, thereby decreasing internal resistance and electrical heat. Then, by controlling the circuit in the system, the double separation relay adjusts the inductance and capacitance of the conductor, thus reducing the reactance length. Additionally, particle swarm particles are mutated to enable them to jump out of the local optimum, and elite fish approach is used to expand the search area. Finally, the proposed fusion algorithm is applied to the self-built data set of Elecgrid and compared with the other three algorithms. The fusion algorithm serves as the standard test system for this comparison. The active power loss of the hybrid algorithm is 63 kW at an operating voltage of 0.74 V. The loss work of the other three algorithms is 74 kW, 97 kW and 109 kW respectively. The mixed algorithm has the lowest loss among the four algorithms. The experiments are repeated for six times, and the linear fitting degrees of the four algorithms are 0.9804, 0.9527, 0.9612 and 0.9503, respectively. The experimental results show that the application of this algorithm can effectively reduce the active loss in the process of distribution network reconfiguration, thus reducing energy consumption; At the same time, it can reduce the electric heating in the process of electric energy transmission, and then prevent the occurrence of fire. There are three main contributions of this study. Firstly, the resistance in the transmission path is reduced by using this algorithm, so that the power transmission efficiency can be analyzed more accurately. Secondly, the new algorithm enriches the power safety maintenance method; Finally, the fire caused by local overheating of the line is reduced by fusion algorithm.

Dynamics of settling disks in a cross-flow

Physical properties of rocks depend strongly on its grain packing which is determined by the natural process involving fluid-structure interaction during deposition. Here, we address the dynamics of settling solids with different sizes in a cross-flow numerically. This is an important aspect in rock formation process. Computational fluid-structure interaction usually involves the use of body-conformed grid. For the case of solid moving in a viscous fluid, the body-conformed grid has to be reconstructed in every time iteration which can be time consuming. The so-called Immersed Boundary method can be used to eliminate this grid reconstruction process. In this study, we use the couple of Immersed Boundary and Lattice Boltzmann Method (usually termed as IBLBM). In lattice Boltzmann, one does not require solving Poisson equation for pressure and a computation of intermediate field as in the usual fractional step procedure since velocity and pressure field can be obtained from the locally solved distribution function. These advantages make IBLBM can be easily employed. In this study, we perform simulations on the dynamics of a cluster of disks with two different values of diameter in a cross-flow. We compute the velocity distribution of the disks to study the sedimentation dynamics and we compare the results with a sedimentation process in a quiescent fluid. For the solid-fluid density ratio of 5 and incoming flow Reynolds number of 12, the velocity distribution is close to the normal distribution as also observed in the case of sedimentation in quiescent fluid.

Research on relay protection technology of high temperature superconducting cable system

High-temperature superconducting cables have many advantages such as low line loss, large transmission capacity, small space occupation, and environmental friendliness, which provide a highly efficient, compact, reliable, and green way of power transmission for power grids, and have great potential for application in large-capacity power transmission and urban power grid expansion and reconstruction. However, high-temperature superconducting cables have a complex structure, which is quite different from traditional cable lines in terms of operating conditions and fault mechanisms, and the requirements for the secondary system are also more stringent. Relying on the project of China Power Grid, “Research and Demonstration Application of Superconducting Transmission Technology in Urban Distribution Networks”, the thesis carries out an in-depth research on relay protection technology adapted to the access of high-temperature superconducting cables.

Model-based reconstruction for looping-star MRI.

The aim of this study was to develop a reconstruction method that more fully models the signals and reconstructs gradient echo (GRE) images without sacrificing the signal to noise ratio and spatial resolution, compared to conventional gridding and model-based image reconstruction method. By modeling the trajectories for every spoke and simplifying the scenario to only echo-in and echo-out mixture, the approach explicitly models the overlapping echoes. After modeling the overlapping echoes with two system matrices, we use the conjugate gradient algorithm (CG-SENSE) with the nonuniform FFT (NUFFT) to optimize the image reconstruction cost function. The proposed method is demonstrated in phantoms and in-vivo volunteer experiments for three-dimensional, high-resolution T2*-weighted imaging and functional MRI tasks. Compared to the gridding method, the high resolution protocol exhibits improved spatial resolution and reduced signal loss as a result of less intra-voxel dephasing. The fMRI task shows that the proposed model-based method produced images with reduced artifacts and blurring as well as more stable and prominent time courses. The proposed model-based reconstruction results shows improved spatial resolution and reduced artifacts. The fMRI task shows improved time series and activation map due to the reduced overlapping echoes and under-sampling artifacts.

Increasing prevalence of hot drought across western North America since the 16th century.

Across western North America (WNA), 20th-21st century anthropogenic warming has increased the prevalence and severity of concurrent drought and heat events, also termed hot droughts. However, the lack of independent spatial reconstructions of both soil moisture and temperature limits the potential to identify these events in the past and to place them in a long-term context. We develop the Western North American Temperature Atlas (WNATA), a data-independent 0.5° gridded reconstruction of summer maximum temperatures back to the 16th century. Our evaluation of the WNATA with existing hydroclimate reconstructions reveals an increasing association between maximum temperature and drought severity in recent decades, relative to the past five centuries. The synthesis of these paleo-reconstructions indicates that the amplification of the modern WNA megadrought by increased temperatures and the frequency and spatial extent of compound hot and dry conditions in the 21st century are likely unprecedented since at least the 16th century.

Long-term gridded land evapotranspiration reconstruction using Deep Forest with high generalizability

Previous datasets have limitations in generalizing evapotranspiration (ET) across various land cover types due to the scarcity and spatial heterogeneity of observations, along with the incomplete understanding of underlying physical mechanisms as a deeper contributing factor. To fill in these gaps, here we developed a global Highly Generalized Land (HG-Land) ET dataset at 0.5° spatial resolution with monthly values covering the satellite era (1982–2018). Our approach leverages the power of a Deep Forest machine-learning algorithm, which ensures good generalizability and mitigates overfitting by minimizing hyper-parameterization. Model explanations are further provided to enhance model transparency and gain new insights into the ET process. Validation conducted at both the site and basin scales attests to the dataset’s satisfactory accuracy, with a pronounced emphasis on the Northern Hemisphere. Furthermore, we find that the primary driver of ET predictions varies across different climatic regions. Overall, the HG-Land ET, underpinned by the interpretability of the machine-learning model, emerges as a validated and generalized resource catering to scientific research and various applications.

A cost-effective and high-efficient EV shared fast charging scheme with hierarchical coordinated operation strategy for addressing difficult-to-charge issue in old residential communities

Charging infrastructures are being deployed widely in urban and rural areas to solve the so-called “last one-kilometer” charging problem. However, for old residential communities, the existing difficult-to-charge issue is not easily solved solely by the widespread deployment of charging infrastructure owing to limited spaces, complex reconstruction, and huge costs. In this context, this paper introduces a cost-effective and high-efficient shared fast charging scheme for old residential communities. The comprised registration-based queuing charging mechanism is proposed to schedule electric vehicle (EV) charging. An integrated and multi-stage power electronic device called smart interlinking unit (SIU) is introduced for grid reconstruction. The SIU-based charging station architecture features flexible control capability, enabling DC power supply for terminal EV chargers. For this architecture, a novel hierarchical coordinated operation strategy is designed, comprising the upper-level real-time vehicle-to-grid (V2G) dispatching and lower-level coordination strategies. Case studies demonstrate that the proposed scheme could accommodate 40 % EV penetration, while exhibiting a significant decline in system investment and occupied area compared to the slow charging scheme. Additionally, a 1.3 % improvement in average transfer efficiency is verified under this architecture. The proposed strategy is proven to realize a prominent power peak-shaving and valley-filling, effective transformer load balance, and DC voltage restoration.

Gridded reconstruction of cropland cover changes in Northeast China from AD 1000 to 1200

Gridded reconstruction of cropland cover changes in Northeast China from AD 1000 to 1200

Investigating Inherent Numerical Stabilization for the Moist, Compressible, Non‐Hydrostatic Euler Equations on Collocated Grids

AbstractThis study investigates inherent numerical dissipation due to upwind fluxes and reconstruction strategies for collocated Finite‐Volume integration of the Euler equations. Idealized supercell simulations are used without any explicit dissipation. Flux terms are split into: mass flux, pressure, and advected quantities. They are computed with the following upwind strategies: central, advectively upwind, and acoustically upwind. This is performed for third and ninth‐order‐accurate reconstructions with and without Weighted Essentially Non‐Oscillatory limiting. Acoustic‐only upwinding for pressure and mass flux terms and advective‐only upwinding for advected quantities is the most flexible simplification found. It reduces data movement and computations. Assuming a constant speed of sound in acoustic upwinding gives similar results to using the true speed of sound. Dissipation from upwind adapts automatically to grid spacing, time step, reconstruction accuracy, and flow smoothness. While stability is maintained even at 21st‐order spatial accuracy, there is a limit to the spatial order of accuracy for which upwinding alone can create a realizable solution in the conditions of this study. Convex combinations of upwind and central solutions for flux terms also reduced dissipation, but as the central proportion grows, solutions become physically unrealizable. The range of length scales of the kinetic energy spectra can be extended along k−5/3 to smaller spatial scales by reducing dissipation either with higher‐order reconstructions or using convex combinations of upwind and central fluxes. However, not all extensions of the length scale range along k−5/3 exhibit physically realizable solutions, even though the spectra appear to be physical.