Numerical modeling of 3D DC resistivity method in the mixed space-wavenumber domain

Forward modeling is the basis of inversion imaging and quantitative interpretation for DC resistivity exploration.Currently,a numerical model of the DC resistivity method must be finely divided to obtain a highly accurate solution under complex conditions,resulting in a long calculation time and large storage.Therefore,we propose a 3D numerical simulation method in a mixed space-wavenumber domain to overcome this challenge.The partial differential equation about abnormal potential is transformed into many independent ordinary differential equations with different wavenumbers using a 2D Fourier transform along the x axis and y axis direction.In this way,a large-scale 3D numerical simulation problem is decomposed into several 1D numerical simulation problems,which significantly reduces the computational and storage requirements.In addition,these ordinary 1D differential equations with different wavenumbers are independent of each other and high parallelelism of the algorithm.They are solved using a finite-element algorithm combined with a chasing method,and the obtained solution is modified using a contraction operator.In this method,the vertical direction is reserved as the spatial domain,then grid size can be determined flexibly based on the underground current density distribution,which considers the solution accuracy and calculation efficiency.In addition,for the first time,we use the contraction operator in the integral equation method to iterate the algorithm.The algorithm takes advantage of the high efficiency of the standard Fourier transform and chasing method,as well as the fast convergence of the contraction operator.We verified the accuracy of the algorithm and the convergence of the contraction operator.Compared with a volume integral method and goal-oriented adaptive finite-element method,the proposed algorithm has lower memory requirements and high computational efficiency,making it suitable for calculating a model with large-scale nodes.Moreover,different examples are used to verify the high adaptability and parallelism of the proposed algorithm.The findings show that the 3D numerical simulation method of DC resistivity method in a mixed space-wavenumber domain is highly efficient,precise,and parallel.

Regional Multi-Agent Cooperative Reinforcement Learning for City-Level Traffic Grid Signal

This article studies the effective traffic signal control problem of multiple intersections in a city-level traffic system.A novel regional multi-agent cooperative reinforcement learning algorithm called RegionSTLight is proposed to improve the traffic efficiency.Firstly a regional multi-agent Q-learning framework is proposed,which can equivalently decompose the global Q value of the traffic system into the local values of several regions Based on the framework and the idea of human-machine cooperation,a dynamic zoning method is designed to divide the traffic network into several strong-coupled regions according to realtime traffic flow densities.In order to achieve better cooperation inside each region,a lightweight spatio-temporal fusion feature extraction network is designed.The experiments in synthetic real-world and city-level scenarios show that the proposed RegionS TLight converges more quickly,is more stable,and obtains better asymptotic performance compared to state-of-theart models.

Diagnostic Study of Global Energy Cycle of the GRAPES Global Model in the Mixed Space-Time Domain

Some important diagnostic characteristics for a model’s physical background are reflected in the model’s energy transport, conversion, and cycle. Diagnosing the atmospheric energy cycle is a suitable way towards understanding and improving numerical models. In this study, formulations of the “Mixed Space-Time Domain”energy cycle are calculated and the roles of stationary and transient waves within the atmospheric energy cycle of the Global-Regional Assimilation and Prediction System （GRAPES） model are diagnosed and compared with the NCEP analysis data for July 2011. Contributions of the zonal-mean components of the energy cycle are investigated to explain the performance of numerical models. The results show that the GRAPES model has the capability to reproduce the main features of the global energy cycle as compared with the NCEP analysis. Zonal available potential energy （AZ） is converted into stationary eddy available potential energy （ASE） and transient eddy available potential energy （ATE）, and ASE and ATE have similar values. The nonlinear conversion between the two eddy energy terms is directed from the stationary to the transient. AZ becomes larger with increased forecast lead time, reflecting an enhancement of the meridional temperature gradient, which strengthens the zonal baroclinic processes and makes the conversion from AZ to eddy potential energy larger, especially for CAT （conversion from AZ to ATE）. The zonal kinetic energy （KZ） has a similar value to the sum of the stationary and transient eddy kinetic energy. Barotropic conversions are directed from eddy to zonal kinetic energy. The zonal conversion from AZ to KZ in GRAPES is around 1.5 times larger than in the NCEP analysis. The contributions of zonal energy cycle components show that transient eddy kinetic energy （KTE） is associated with the Southern Hemisphere subtropical jet and the conversion from KZ to KTE reduces in the upper tropopause near 30？S. The nonlinear barotropic conversion between stationary and transient kinetic energy terms （CKTE） is reduced predominantly by the weaker KTE.

AN ITERATIVE PROCEDURE FOR DOMAIN DECOMPOSITION METHOD OF SECOND ORDER ELLIPTIC PROBLEM WITH MIXED BOUNDARY CONDITIONS

This paper is devoted to study of an iterative procedure for domain decomposition method of second order elliptic problem with mixed boundary conditions (i.e., Dirichlet condition on a part of boundary and Neumann condition on the another part of boundary). For the pure Dirichlet problem, Marini and Quarteroni [3], [4] considered a similar approach, which is extended to more complex problem in this paper.

NONOVERLAPPING DOMAIN DECOMPOSITION METHOD WITH MIXED ELEMENT FOR ELLIPTIC PROBLEMS

In this paper we consider the nonoverlapping domain decomposition method based on mixed element approximation for elliptic problems in two dimentional space. We give a kind of discrete domain decomposition iterative algorithm using mixed finite element, the subdomain problems of which can be implemented parallelly. We also give the existence, uniqueness and convergence of the approximate solution.