3D forward modeling of magnetotelluric fields in general anisotropic media and its numerical implementation in Julia

Geophysics ◽  
2018 ◽  
Vol 83 (4) ◽  
pp. F29-F40 ◽  
Author(s):  
Bo Han ◽  
Yuguo Li ◽  
Gang Li

We have developed a finite volume (FV) algorithm for magnetotelluric (MT) forward modeling in 3D conductivity structures with general anisotropy. The electric and magnetic fields are discretized on a conventional staggered grid, which cannot directly address the full-tensor conductivity. To overcome this difficulty, an interpolation scheme is used to average different components of the electric field to the same position. We formulate the algorithm in pure matrix form and implement it in a new language, Julia, making the programming process highly efficient and leading to a code with excellent readability, maintainability, and extendability. The validity of the FV Julia code is demonstrated using a layered 1D anisotropic model. For this model, the FV code provides accurate results, and the computational cost is reasonable. Being preconditioned with the electromagnetic potential ([Formula: see text]) system, the iterative solvers including quasi-minimal residual and biconjugate gradient stabilized exhibit a good convergence rate for a wide range of periods. The direct solvers MUMPS and PARDISO are highly efficient for small model sizes. For a relatively large model size with 2.18 millions unknowns, the linear system of one period can be solved by MUMPS within 360 s with multiple threads involved in the computation, and the memory usage is only 11.6 GB in the “out-of-core” mode. We further calculated MT responses of a 3D model with dipping and horizontal anisotropy, respectively. The results suggest that the electrical anisotropy can have significant influence on the MT response.

SPE Journal ◽  
2021 ◽  
pp. 1-19
Author(s):  
Yingnan Wang ◽  
Nadia Shardt ◽  
Janet A. W. Elliott ◽  
Zhehui Jin

Summary Gas-alkane interfacial tension (IFT) is an important parameter in the enhanced oil recovery (EOR) process. Thus, it is imperative to obtain an accurate gas-alkane mixture IFT for both chemical and petroleum engineering applications. Various empirical correlations have been developed in the past several decades. Although these models are often easy to implement, their accuracy is inconsistent over a wide range of temperatures, pressures, and compositions. Although statistical mechanics-based models and molecular simulations can accurately predict gas-alkane IFT, they usually come with an extensive computational cost. The Shardt-Elliott (SE) model is a highly accurate IFT model that for subcritical fluids is analytic in terms of temperature T and composition x. In applications, it is desirable to obtain IFT in terms of temperature T and pressure P, which requires time-consuming flash calculations, and for mixtures that contain a gas component greater than its pure species critical point, additional critical composition calculations are required. In this work, the SE model is combined with a machine learning (ML) approach to obtain highly efficient and highly accurate gas-alkane binary mixture IFT equations directly in terms of temperature, pressure, and alkane molar weights. The SE model is used to build an IFT database (more than 36,000 points) for ML training to obtain IFT equations. The ML-based IFT equations are evaluated in comparison with the available experimental data (888 points) and with the SE model, as well as with the less accurate parachor model. Overall, the ML-based IFT equations show excellent agreement with experimental data for gas-alkane binary mixtures over a wide range of T and P, and they outperform the widely used parachor model. The developed highly efficient and highly accurate IFT functions can serve as a basis for modeling gas-alkane binary mixtures for a broad range of T, P, and x.


CALCOLO ◽  
2021 ◽  
Vol 58 (4) ◽  
Author(s):  
Marco Donatelli ◽  
Rolf Krause ◽  
Mariarosa Mazza ◽  
Ken Trotti

AbstractWe focus on a time-dependent one-dimensional space-fractional diffusion equation with constant diffusion coefficients. An all-at-once rephrasing of the discretized problem, obtained by considering the time as an additional dimension, yields a large block linear system and paves the way for parallelization. In particular, in case of uniform space–time meshes, the coefficient matrix shows a two-level Toeplitz structure, and such structure can be leveraged to build ad-hoc iterative solvers that aim at ensuring an overall computational cost independent of time. In this direction, we study the behavior of certain multigrid strategies with both semi- and full-coarsening that properly take into account the sources of anisotropy of the problem caused by the grid choice and the diffusion coefficients. The performances of the aforementioned multigrid methods reveal sensitive to the choice of the time discretization scheme. Many tests show that Crank–Nicolson prevents the multigrid to yield good convergence results, while second-order backward-difference scheme is shown to be unconditionally stable and that it allows good convergence under certain conditions on the grid and the diffusion coefficients. The effectiveness of our proposal is numerically confirmed in the case of variable coefficients too and a two-dimensional example is given.


2021 ◽  
Author(s):  
Qiu-Hong Huang ◽  
Qian-Yi Zhou ◽  
Chen Yang ◽  
Li Chen ◽  
Jin-Pei Cheng ◽  
...  

A highly efficient desymmetrizing asymmetric bromination of bisphenol phosphine oxides was developed, providing a wide range of chiral bisphenol phosphine oxides and bisphenol phosphinates with high yields and enantioselectivities.


2021 ◽  
Vol 12 (1) ◽  
Author(s):  
Alexander D. Taylor ◽  
Qing Sun ◽  
Katelyn P. Goetz ◽  
Qingzhi An ◽  
Tim Schramm ◽  
...  

AbstractDeposition of perovskite films by antisolvent engineering is a highly common method employed in perovskite photovoltaics research. Herein, we report on a general method that allows for the fabrication of highly efficient perovskite solar cells by any antisolvent via manipulation of the antisolvent application rate. Through detailed structural, compositional, and microstructural characterization of perovskite layers fabricated by 14 different antisolvents, we identify two key factors that influence the quality of the perovskite layer: the solubility of the organic precursors in the antisolvent and its miscibility with the host solvent(s) of the perovskite precursor solution, which combine to produce rate-dependent behavior during the antisolvent application step. Leveraging this, we produce devices with power conversion efficiencies (PCEs) that exceed 21% using a wide range of antisolvents. Moreover, we demonstrate that employing the optimal antisolvent application procedure allows for highly efficient solar cells to be fabricated from a broad range of precursor stoichiometries.


Geophysics ◽  
2021 ◽  
pp. 1-78
Author(s):  
Zhiyuan Li ◽  
Youshan Liu ◽  
Guanghe Liang ◽  
Guoqiang Xue ◽  
Runjie Wang

The separation of P- and S-wavefields is considered to be an effective approach for eliminating wave-mode cross-talk in elastic reverse-time migration. At present, the Helmholtz decomposition method is widely used for isotropic media. However, it tends to change the amplitudes and phases of the separated wavefields compared with the original wavefields. Other methods used to obtain pure P- and S-wavefields include the application of the elastic wave equations of the decoupled wavefields. To achieve a high computational accuracy, staggered-grid finite-difference (FD) schemes are usually used to numerically solve the equations by introducing an additional stress variable. However, the computational cost of this method is high because a conventional hybrid wavefield (P- and S-wavefields are mixed together) simulation must be created before the P- and S-wavefields can be calculated. We developed the first-order particle velocity equations to reduce the computational cost. The equations can describe four types of particle velocity wavefields: the vector P-wavefield, the scalar P-wavefield, the vector S-wavefield, and the vector S-wavefield rotated in the direction of the curl factor. Without introducing the stress variable, only the four types of particle velocity variables are used to construct the staggered-grid FD schemes, so the computational cost is reduced. We also present an algorithm to calculate the P and S propagation vectors using the four particle velocities, which is simpler than the Poynting vector. Finally, we applied the velocity equations and propagation vectors to elastic reverse-time migration and angle-domain common-image gather computations. These numerical examples illustrate the efficiency of the proposed methods.


2018 ◽  
Author(s):  
Fabien Maussion ◽  
Anton Butenko ◽  
Julia Eis ◽  
Kévin Fourteau ◽  
Alexander H. Jarosch ◽  
...  

Abstract. Despite of their importance for sea-level rise, seasonal water availability, and as source of geohazards, mountain glaciers are one of the few remaining sub-systems of the global climate system for which no globally applicable, open source, community-driven model exists. Here we present the Open Global Glacier Model (OGGM, http://www.oggm.org), developed to provide a modular and open source numerical model framework for simulating past and future change of any glacier in the world. The modelling chain comprises data downloading tools (glacier outlines, topography, climate, validation data), a preprocessing module, a mass-balance model, a distributed ice thickness estimation model, and an ice flow model. The monthly mass-balance is obtained from gridded climate data and a temperature index melt model. To our knowledge, OGGM is the first global model explicitly simulating glacier dynamics: the model relies on the shallow ice approximation to compute the depth-integrated flux of ice along multiple connected flowlines. In this paper, we describe and illustrate each processing step by applying the model to a selection of glaciers before running global simulations under idealized climate forcings. Even without an in-depth calibration, the model shows a very realistic behaviour. We are able to reproduce earlier estimates of global glacier volume by varying the ice dynamical parameters within a range of plausible values. At the same time, the increased complexity of OGGM compared to other prevalent global glacier models comes at a reasonable computational cost: several dozens of glaciers can be simulated on a personal computer, while global simulations realized in a supercomputing environment take up to a few hours per century. Thanks to the modular framework, modules of various complexity can be added to the codebase, allowing to run new kinds of model intercomparisons in a controlled environment. Future developments will add new physical processes to the model as well as tools to calibrate the model in a more comprehensive way. OGGM spans a wide range of applications, from ice-climate interaction studies at millenial time scales to estimates of the contribution of glaciers to past and future sea-level change. It has the potential to become a self-sustained, community driven model for global and regional glacier evolution.


Geophysics ◽  
2014 ◽  
Vol 79 (6) ◽  
pp. T313-T321 ◽  
Author(s):  
Hanming Chen ◽  
Hui Zhou ◽  
Yanqi Li

A classical split perfectly matched layer (PML) method has recently been applied to the scalar arbitrarily wide-angle wave equation (AWWE) in terms of displacement. However, the classical split PML obviously increases computational cost and cannot efficiently absorb waves propagating into the absorbing layer at grazing incidence. Our goal was to improve the computational efficiency of AWWE and to enhance the suppression of edge reflections by applying a convolutional PML (CPML). We reformulated the original AWWE as a first-order formulation and incorporated the CPML with a general complex frequency shifted stretching operator into the renewed formulation. A staggered-grid finite-difference (FD) method was adopted to discretize the first-order equation system. For wavefield depth continuation, the first-order AWWE with the CPML saved memory compared with the original second-order AWWE with the conventional split PML. With the help of numerical examples, we verified the correctness of the staggered-grid FD method and concluded that the CPML can efficiently absorb evanescent and propagating waves.


2016 ◽  
Vol 2016 ◽  
pp. 1-15 ◽  
Author(s):  
Hongjie Guo ◽  
Guojun Dai ◽  
Jin Fan ◽  
Yifan Wu ◽  
Fangyao Shen ◽  
...  

This paper develops a mobile sensing system, the first system used in adaptive resolution urban air quality monitoring. In this system, we employ several taxis as sensor carries to collect originalPM2.5data and collect a variety of datasets, including meteorological data, traffic status data, and geographical data in the city. This paper also presents a novel method AG-PCEM (Adaptive Grid-Probabilistic Concentration Estimation Method) to infer thePM2.5concentration for undetected grids using dynamic adaptive grids. We gradually collect the measurements throughout a year using a prototype system in Xiasha District of Hangzhou City, China. Experimental data has verified that the proposed system can achieve good performance in terms of computational cost and accuracy. The computational cost of AG-PCEM is reduced by about 40.2% compared with a static grid method PCEM under the condition of reaching the close accuracy, and the accuracy of AG-PCEM is far superior as widely used artificial neural network (ANN) and Gaussian process (GP), enhanced by 38.8% and 14.6%, respectively. The system can be expanded to wide-range air quality monitor by adjusting the initial grid resolution, and our findings can tell citizens actual air quality and help official management find pollution sources.


2019 ◽  
Vol 99 (2) ◽  
pp. 1105-1130 ◽  
Author(s):  
Kun Yang ◽  
Vladimir Paramygin ◽  
Y. Peter Sheng

Abstract The joint probability method (JPM) is the traditional way to determine the base flood elevation due to storm surge, and it usually requires simulation of storm surge response from tens of thousands of synthetic storms. The simulated storm surge is combined with probabilistic storm rates to create flood maps with various return periods. However, the map production requires enormous computational cost if state-of-the-art hydrodynamic models with high-resolution numerical grids are used; hence, optimal sampling (JPM-OS) with a small number of (~ 100–200) optimal (representative) storms is preferred. This paper presents a significantly improved JPM-OS, where a small number of optimal storms are objectively selected, and simulated storm surge responses of tens of thousands of storms are accurately interpolated from those for the optimal storms using a highly efficient kriging surrogate model. This study focuses on Southwest Florida and considers ~ 150 optimal storms that are selected based on simulations using either the low fidelity (with low resolution and simple physics) SLOSH model or the high fidelity (with high resolution and comprehensive physics) CH3D model. Surge responses to the optimal storms are simulated using both SLOSH and CH3D, and the flood elevations are calculated using JPM-OS with highly efficient kriging interpolations. For verification, the probabilistic inundation maps are compared to those obtained by the traditional JPM and variations of JPM-OS that employ different interpolation schemes, and computed probabilistic water levels are compared to those calculated by historical storm methods. The inundation maps obtained with the JPM-OS differ less than 10% from those obtained with JPM for 20,625 storms, with only 4% of the computational time.


Sign in / Sign up

Export Citation Format

Share Document