To speed up three-dimensional (3D) DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method (AGMG). We first discretize the differential equation of the secondar...To speed up three-dimensional (3D) DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method (AGMG). We first discretize the differential equation of the secondary potential field with mixed boundary conditions by using a seven-point finite-difference method to obtain a large sparse system of linear equations. Then, we introduce the theory behind the pairwise aggregation algorithms for AGMG and use the conjugate-gradient method with the V-cycle AGMG preconditioner (AGMG-CG) to solve the linear equations. We use typical geoelectrical models to test the proposed AGMG-CG method and compare the results with analytical solutions and the 3DDCXH algorithm for 3D DC modeling (3DDCXH). In addition, we apply the AGMG-CG method to different grid sizes and geoelectrical models and compare it to different iterative methods, such as ILU-BICGSTAB, ILU-GCR, and SSOR-CG. The AGMG-CG method yields nearly linearly decreasing errors, whereas the number of iterations increases slowly with increasing grid size. The AGMG-CG method is precise and converges fast, and thus can improve the computational efficiency in forward modeling of three-dimensional DC resistivity.展开更多
To deal with the problem of low computational precision at the nodes near the source and satisfy the requirements for computational efficiency in inversion imaging and finite-element numerical simulations of the direc...To deal with the problem of low computational precision at the nodes near the source and satisfy the requirements for computational efficiency in inversion imaging and finite-element numerical simulations of the direct current method, we propose a new mesh refinement and recoarsement method for a two-dimensional point source. We introduce the mesh refinement and mesh recoarsement into the traditional structured mesh subdivision. By refining the horizontal grids, the singularity owing to the point source is minimized and the topography is simulated. By recoarsening the horizontal grids, the number of grid cells is reduced significantly and computational efficiency is improved. Model tests show that the proposed method solves the singularity problem and reduces the number of grid cells by 80% compared to the uniform grid refinement.展开更多
The structure, separation principle and feasibility research for a new type of vehicle air filter called the high speed rotary positive air filter were described. The analysis of the experimental data showed that the ...The structure, separation principle and feasibility research for a new type of vehicle air filter called the high speed rotary positive air filter were described. The analysis of the experimental data showed that the principle and structure of it were feasible and it possessed high separation efficiency and great self cleaning ability. Compared with the conventional air filter it also has lower air intake loss. So it is worth further practical research.展开更多
As a fundamental study on recovery of valuable metals from nonferrous metallurgical slags,electrical conductivity values of MO(MO=FeO,NiO)-containing CaO-MgO-SiO2-Al2O3 slag with a low basicity were measured at diff...As a fundamental study on recovery of valuable metals from nonferrous metallurgical slags,electrical conductivity values of MO(MO=FeO,NiO)-containing CaO-MgO-SiO2-Al2O3 slag with a low basicity were measured at different temperatures using AC impedance spectroscopy.The result shows that the electrical conductivity increased from 1.4 S/m to 14.4 S/m with the increase of the temperature from 1 573 to 1 773 K and the content of MO which is less than 12% under the constant mass ratio of (CaO+MgO) to (SiO2+Al2O3) of 0.47.Moreover,the increase magnitude of the electrical conductivity was also promoted with the increase of the content of MO.The electrical conductivity of FeO-containing slags was close to that of NiO-containing slags when the content was less than 8%;however,it was obviously larger than that of NiO-containing slags when the content was 12%.The activation energy of the electrical conductivity decreased with the increase of MO content.展开更多
Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid...Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid finite element method was proposed. In this algorithm, a-posteriori error estimator was employed to generate adaptively refined mesh on a given initial mesh. On these iterative meshes, V-cycle based multigrid method was adopted to fast solve each linear equation with each initial iterative term interpolated from last mesh. With this error estimator, the unknowns were nearly optimally distributed on the final mesh which guaranteed the accuracy. The numerical results show that the multigrid solver is faster and more stable compared with ICCG solver. Meanwhile, the numerical results obtained from the final model discretization approximate the analytical solutions with maximal relative errors less than 1%, which remarkably validates this algorithm.展开更多
Based on an analytical solution for the current point source in an anisotropic half-space,we study the apparent resistivity and apparent chargeability of a transversely isotropic medium with vertical and horizontal ax...Based on an analytical solution for the current point source in an anisotropic half-space,we study the apparent resistivity and apparent chargeability of a transversely isotropic medium with vertical and horizontal axes symmetry,respectively.We then provide a simple derivation of the anisotropy paradoxes in direct current resistivity and time-domain induced polarization methods.Analogous to the mean resistivity,we propose a formulation for deriving the mean polarizability.We also present a three-dimensional finite element algorithm for modeling the direct current resistivity and time-domain induced polarization using an unstructured tetrahedral grid.Finally,we provide the apparent resistivity and apparent chargeability curves of a tilted,transversely isotropic medium with diff erent angles,respectively.The subsequent results illustrate the anisotropy paradoxes of direct current resistivity and time-domain induced polarization.展开更多
The strategies that minimize the overall solution time of multiple linear systems in 3D finite element method (FEM) modeling of direct current (DC) resistivity were discussed. A global stiff matrix is assembled and st...The strategies that minimize the overall solution time of multiple linear systems in 3D finite element method (FEM) modeling of direct current (DC) resistivity were discussed. A global stiff matrix is assembled and stored in two parts separately. One part is associated with the volume integral and the other is associated with the subsurface boundary integral. The equivalent multiple linear systems with closer right-hand sides than the original systems were constructed. A recycling Krylov subspace technique was employed to solve the multiple linear systems. The solution of the seed system was used as an initial guess for the subsequent systems. The results of two numerical experiments show that the improved algorithm reduces the iterations and CPU time by almost 50%, compared with the classical preconditioned conjugate gradient method.展开更多
A fixed artificial source(greater than 200 kW) was used and the source location was selected at a high resistivity region to ensure high emission efficiency. Some publications used the "earth-ionosphere" mod...A fixed artificial source(greater than 200 kW) was used and the source location was selected at a high resistivity region to ensure high emission efficiency. Some publications used the "earth-ionosphere" mode in modeling the electromagnetic(EM) fields with the offset up to a thousand kilometer, and such EM fields still have a signal/noise ratio of 10-20 dB. This means that a new EM method with fixed source is feasible, but in their calculation, the displacement in air was neglected. In this work, some three-layer modeling results were presented to illustrate the basic EM fields' characteristics in the near, far and waveguide areas under "earth-ionosphere" mode, and a standard is given to distinguish the boundary of near, far and waveguide areas. Due to the influence of the ionosphere and displacement current in the air, the "earth-ionosphere" mode EM fields have an extra waveguide zone, where the fields' behavior is very different from that of the far field zone.展开更多
In order to explore the influence of soil resistivity on stray current in power supply system of urban rail transit,we establish an equivalent circuit model of the rail-to-ground structure based on resistance network ...In order to explore the influence of soil resistivity on stray current in power supply system of urban rail transit,we establish an equivalent circuit model of the rail-to-ground structure based on resistance network method first.After measuring the soil resistivity of a real subway system,a simulation model is established in Matlab to obtain the stray currents at different soil resistivities.Then the influence of soil resistivity on stray current is analyzed.Afterwards,to verify the rationality and reliability of the simulation model,we design a test circuit to measure the stray current and rail-to-ground voltage in a real subway system,and a comparison of the experimentally measured results and simulation results is presented.The results show that the stray current is the maximum when the soil resistivity is 211.57Ω·m;when the soil resistivity is 768.47Ω·m,the stray current is the minimum,that is,the smaller the soil resistivity,the greater the stray current.Therefore,the resistivity should be increased as much as possible when ramming the track foundation in urban rail transit system.展开更多
基金supported by the Natural Science Foundation of China(Nos.41404057,41674077 and 411640034)the Nuclear Energy Development Project of China,and the‘555’Project of Gan Po Excellent People
文摘To speed up three-dimensional (3D) DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method (AGMG). We first discretize the differential equation of the secondary potential field with mixed boundary conditions by using a seven-point finite-difference method to obtain a large sparse system of linear equations. Then, we introduce the theory behind the pairwise aggregation algorithms for AGMG and use the conjugate-gradient method with the V-cycle AGMG preconditioner (AGMG-CG) to solve the linear equations. We use typical geoelectrical models to test the proposed AGMG-CG method and compare the results with analytical solutions and the 3DDCXH algorithm for 3D DC modeling (3DDCXH). In addition, we apply the AGMG-CG method to different grid sizes and geoelectrical models and compare it to different iterative methods, such as ILU-BICGSTAB, ILU-GCR, and SSOR-CG. The AGMG-CG method yields nearly linearly decreasing errors, whereas the number of iterations increases slowly with increasing grid size. The AGMG-CG method is precise and converges fast, and thus can improve the computational efficiency in forward modeling of three-dimensional DC resistivity.
基金financially supported by the National Natural Science Foundation of China(No.41574127 and 41174104)the National Key Technology R&D Program for the 13th five-year plan(No.2016ZX05018006-006)
文摘To deal with the problem of low computational precision at the nodes near the source and satisfy the requirements for computational efficiency in inversion imaging and finite-element numerical simulations of the direct current method, we propose a new mesh refinement and recoarsement method for a two-dimensional point source. We introduce the mesh refinement and mesh recoarsement into the traditional structured mesh subdivision. By refining the horizontal grids, the singularity owing to the point source is minimized and the topography is simulated. By recoarsening the horizontal grids, the number of grid cells is reduced significantly and computational efficiency is improved. Model tests show that the proposed method solves the singularity problem and reduces the number of grid cells by 80% compared to the uniform grid refinement.
文摘The structure, separation principle and feasibility research for a new type of vehicle air filter called the high speed rotary positive air filter were described. The analysis of the experimental data showed that the principle and structure of it were feasible and it possessed high separation efficiency and great self cleaning ability. Compared with the conventional air filter it also has lower air intake loss. So it is worth further practical research.
基金Project(50574011) supported by the National Natural Science Foundation of China
文摘As a fundamental study on recovery of valuable metals from nonferrous metallurgical slags,electrical conductivity values of MO(MO=FeO,NiO)-containing CaO-MgO-SiO2-Al2O3 slag with a low basicity were measured at different temperatures using AC impedance spectroscopy.The result shows that the electrical conductivity increased from 1.4 S/m to 14.4 S/m with the increase of the temperature from 1 573 to 1 773 K and the content of MO which is less than 12% under the constant mass ratio of (CaO+MgO) to (SiO2+Al2O3) of 0.47.Moreover,the increase magnitude of the electrical conductivity was also promoted with the increase of the content of MO.The electrical conductivity of FeO-containing slags was close to that of NiO-containing slags when the content was less than 8%;however,it was obviously larger than that of NiO-containing slags when the content was 12%.The activation energy of the electrical conductivity decreased with the increase of MO content.
基金Projects(2006AA06Z105, 2007AA06Z134) supported by the National High-Tech Research and Development Program of ChinaProjects(2007, 2008) supported by China Scholarship Council (CSC)
文摘Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid finite element method was proposed. In this algorithm, a-posteriori error estimator was employed to generate adaptively refined mesh on a given initial mesh. On these iterative meshes, V-cycle based multigrid method was adopted to fast solve each linear equation with each initial iterative term interpolated from last mesh. With this error estimator, the unknowns were nearly optimally distributed on the final mesh which guaranteed the accuracy. The numerical results show that the multigrid solver is faster and more stable compared with ICCG solver. Meanwhile, the numerical results obtained from the final model discretization approximate the analytical solutions with maximal relative errors less than 1%, which remarkably validates this algorithm.
基金the special funding of Guiyang science and technology bureau and Guiyang University[GYUKY-[2021]]the National Key Research and Development Program of China-Geophysical Comprehensive Exploration and Information Extraction of Deep Mineral Resources(2016YFC0600505)the National K&D Program(2018YFC1504901,2018YFC1504904).
文摘Based on an analytical solution for the current point source in an anisotropic half-space,we study the apparent resistivity and apparent chargeability of a transversely isotropic medium with vertical and horizontal axes symmetry,respectively.We then provide a simple derivation of the anisotropy paradoxes in direct current resistivity and time-domain induced polarization methods.Analogous to the mean resistivity,we propose a formulation for deriving the mean polarizability.We also present a three-dimensional finite element algorithm for modeling the direct current resistivity and time-domain induced polarization using an unstructured tetrahedral grid.Finally,we provide the apparent resistivity and apparent chargeability curves of a tilted,transversely isotropic medium with diff erent angles,respectively.The subsequent results illustrate the anisotropy paradoxes of direct current resistivity and time-domain induced polarization.
基金Projects(40974077,41164004)supported by the National Natural Science Foundation of ChinaProject(2007AA06Z134)supported by the National High Technology Research and Development Program of China+2 种基金Projects(2011GXNSFA018003,0832263)supported by the Natural Science Foundation of Guangxi Province,ChinaProject supported by Program for Excellent Talents in Guangxi Higher Education Institution,ChinaProject supported by the Foundation of Guilin University of Technology,China
文摘The strategies that minimize the overall solution time of multiple linear systems in 3D finite element method (FEM) modeling of direct current (DC) resistivity were discussed. A global stiff matrix is assembled and stored in two parts separately. One part is associated with the volume integral and the other is associated with the subsurface boundary integral. The equivalent multiple linear systems with closer right-hand sides than the original systems were constructed. A recycling Krylov subspace technique was employed to solve the multiple linear systems. The solution of the seed system was used as an initial guess for the subsequent systems. The results of two numerical experiments show that the improved algorithm reduces the iterations and CPU time by almost 50%, compared with the classical preconditioned conjugate gradient method.
基金Projects(41204054,41541036,41604111)supported by the National Natural Science Foundation of China
文摘A fixed artificial source(greater than 200 kW) was used and the source location was selected at a high resistivity region to ensure high emission efficiency. Some publications used the "earth-ionosphere" mode in modeling the electromagnetic(EM) fields with the offset up to a thousand kilometer, and such EM fields still have a signal/noise ratio of 10-20 dB. This means that a new EM method with fixed source is feasible, but in their calculation, the displacement in air was neglected. In this work, some three-layer modeling results were presented to illustrate the basic EM fields' characteristics in the near, far and waveguide areas under "earth-ionosphere" mode, and a standard is given to distinguish the boundary of near, far and waveguide areas. Due to the influence of the ionosphere and displacement current in the air, the "earth-ionosphere" mode EM fields have an extra waveguide zone, where the fields' behavior is very different from that of the far field zone.
基金National Natural Science Foundation of China(Nos.51476073,51266004)Natural Science Foundation of Gansu Province(No.1308RJZA199)。
文摘In order to explore the influence of soil resistivity on stray current in power supply system of urban rail transit,we establish an equivalent circuit model of the rail-to-ground structure based on resistance network method first.After measuring the soil resistivity of a real subway system,a simulation model is established in Matlab to obtain the stray currents at different soil resistivities.Then the influence of soil resistivity on stray current is analyzed.Afterwards,to verify the rationality and reliability of the simulation model,we design a test circuit to measure the stray current and rail-to-ground voltage in a real subway system,and a comparison of the experimentally measured results and simulation results is presented.The results show that the stray current is the maximum when the soil resistivity is 211.57Ω·m;when the soil resistivity is 768.47Ω·m,the stray current is the minimum,that is,the smaller the soil resistivity,the greater the stray current.Therefore,the resistivity should be increased as much as possible when ramming the track foundation in urban rail transit system.