Parallel finite element algorithm with domain decomposition for three-dimensional magnetotelluric modelling

We present a new finite element (FE) method for magnetotelluric modelling of three-dimensional conductivity structures. Maxwell's equations are treated as a system of first-order partial differential equations for the secondary fields. Absorbing boundary conditions are introduced, minimizing undesired boundary effects and allowing the use of small computational domains. The numerical algorithm presented here is an iterative, domain decomposition procedure employing a nonconforming FE space. It does not use global matrices, therefore allowing the modellization of large and complicated structures. The algorithm is naturally parallellizable, and we show results obtained in the IBM SP2 parallel supercomputer at Purdue University. The accuracy of the numerical method is verified by checking the computed solutions with the results of COMMEMI, the international project on the comparison of modelling methods for electromagnetic induction.     

Numerical simulations of very large impacts on the Earth

Vertical impacts on the Earth of asteroids 500-3000 km in diameter at 15 km/s have been numerically modelled using the hydrodynamic SOVA code. This code has been modified for the spherical system of coordinates well suited for simulations of very large impacts when the entire Earth is involved in motion. The simulations include cratering process, upward motion of deep mantle layers, fall of ejecta on the Earth, escape of matter to space, and formation of rock vapour atmospheres. The calculations were made for the period preceding disappearance of rock vapour atmospheres caused by radiation several years after the largest impacts. For very large vertical impacts at 15 km/s, escaping masses proved to be negligibly small. Quantities of kinetic, internal, potential, and radiated away energies are obtained as functions of time and space. After the impacts, a global layer of condensed ejecta covers the whole of the Earth's surface and the ejecta energy is sufficient to vaporise an ocean 3 km deep. The mass of rock vapour atmosphere is 10-23% of the impactor mass. This atmosphere has a greater mass than the water atmosphere if impactor is 2000 km in diameter or larger.     

吸积盘谱讨论

本文用数值方法讨论了吸积盘谱与黑体谱的形状和谱指数。结果表明:对于吸积盘用F_v∝v~(1/3)作为其特征谱是不合适的;同时Frank等(1985)所得到的解析式F_v∝v~2也是不正确的。     

Round-off error in long-term orbital integrations using multistep methods

Techniques for reducing roundoff error are compared by testing them on high-order Störmer and summetric multistep methods. The best technique for most applications is to write the equation in summed, function-evaluation form and to store the coefficients as rational numbers. A larger error reduction can be achieved by writing the equation in backward-difference form and performing some of the additions in extended precision, but this entails a larger cpu cost.     

Advances in Numerical Modeling of Astrophysical and Space Plasmas

Advances in the simulation of astrophysical and cosmic plasmas are the direct result of advances in computational capabilities, today consisting of new techniques such as multilevel concurrent simulation, multi-teraflop computational platforms and experimental facilities for producing and diagnosing plasmas under extreme conditions for the benchmarking of simulations. Examples of these are the treatment of mesoscalic plasma and the scaling to astrophysical and cosmic dimensions and the Accelerated Strategic Computing Initiative whose goal is to construct petaflop (10¹⁵ floating operations per second) computers, and pulsed power and laser inertial confinement plasmas where megajoules of energy are delivered to highly-diagnosed plasmas. This paper concentrates on the achievements to date in simulating and experimentally producing plasmas scaled to both astrophysical and cosmic plasma dimensions. A previous paper (Part I, Peratt, 1997) outlines the algorithms and computational growth.     

Numerical determination of pressure-dependent effective thermal conductivity in Berea sandstone

In this paper, the methods of digital rock physics are applied to determine pressure-dependent effective thermal conductivity in rock samples. Simulations are performed with an in-house three-dimensional finite volume code. In the first step, four numerical models are derived from a given tomographic scan of Berea sandstone. Consequently, simulations of the thermal conductivity at ambient conditions are performed and validated with experimental data. In a second step, a new workflow for the determination of the pressure-dependent thermal conductivity in rock samples is elaborated, tested and calibrated. Results originating from the derived workflow show very good agreement with experimental data.     

Experimental and numerical investigation of groundwater head losses on and nearby short intersections between disc-shaped fractures

Discrete fracture models are used for investigating precise processes of groundwater flow in fractured rocks,while a disc-shaped parallel-plates model for a single fracture is more reasonable and efficient for computational treatments.The flow velocity has a large spatial differentiation which is more likely to produce non-linear flow and additional head losses on and nearby intersections in such shaped fractures,therefore it is necessary to understand and quantify them.In this study,both laboratory experiments and numerical simulations were performed to investigate the total head loss on and nearby the intersections as well as the local head loss exactly on the intersections,which were not usually paid sufficient attention or even ignored.The investigation results show that these two losses account for 29.17%-84.97%and 0-73.57%of the entire total head loss in a fracture,respectively.As a result,they should be necessarily considered for groundwater modeling in fractured rocks.Furthermore,both head losses become larger when aperture and flow rate increase and intersection length decreases.Particularly,the ratios of these two head losses to the entire total head loss in a fracture could be well statistically explained by power regression equations with variables of aperture,intersection length,and flow rates,both of which achieved high coefficients of determination.It could be feasible through this type of study to provide a way on how to adjust the groundwater head from those obtained by numerical simulations based on the traditional linear flow model.Finally,it is practicable and effective to implement the investigation approach combining laboratory experiments with numerical simulations for quantifying the head losses on and nearby the intersections between disc-shaped fractures.     

板钉连接CFST柱-RC梁节点梁端塑性铰区受力性能数值模拟

建立竖板-栓钉连接钢管混凝土（CFST）柱-钢筋混凝土（RC）梁节点试件（SSJD）拟静力加载试验有限元模型,并在节点损伤情况、梁端荷载-位移曲线等数值模拟结果与试验结果吻合较好的基础上,进一步开展了RC梁混凝土强度、配筋率ρ_s和连接竖板长度L_b及界面连接情况等对CFST柱-RC梁节点梁端塑性铰区域力学性能的影响。研究结果表明,RC梁混凝土强度对试件SSJD塑性铰区域受力性能的影响较小;适筋范围内RC梁配筋率增加可适当提高试件SSJD承载力和延性;随着连接竖板长度的增加,梁端塑性铰区域外移,梁破坏荷载增大;本研究给出的RC梁与CFST柱之间的界面抗剪承载力模拟值与计算值吻合较好,可用于界面抗剪设计。     

Numerical simulation of deep-hole resistivity anomaly caused by drilling construction in Xinfengjiang geoelectric station

The Xinfengjiang reservoir in Guangdong Province, is one of the large reservoirs that have triggered earthquakes of magnitude greater than 6. Numerous earthquakes have occurred since the impoundment of the reservoir, making it one of the most active seismic zones in south China. In 2015, a set of deep-hole resistivity anomalies was observed in the Heping geoelectric station in Dongyuan county, located near the Xinfengjiang reservoir. After a field investigation, we found that a planned well drilling construction of new measuring channels was being carried out during that corresponding period of time. After careful comparison and analysis on the basis of the collected raw data, we had a reason to believe that drilling construction, rather than the inducement of the Xinfengjiang reservoir, was the main culprit for those unusual georesistivity values. So as to verify the above conjecture, we constructed a series of 3D finite element models based on the geological and hydrological information around Heping station and analyzed the drilling disturbances, respectively. Some significant conclusions were finally drawn according to the precise numerical simulation. This study gives a good example by combining numerical simulation with engineering practice as a way to understand the root cause of georesistivity anomalies in reality.     

Seismic performance of two-story subway station structures with different isolating systems

When the sliding bearing is fixed only at the top of the middle column of the underground structure, the cracks at the side end of the middle plate should be aggravated while the seismic damage of the mid-column should be alleviated. To enhance the seismic performance of the mid-plate, a new isolation design method has been mentioned while the elastic sliding bearings are set at the top of the mid-columns and between the side end of the mid-plate and the side wall at the same time. By establishing a nonlinear finite element analysis model for the static-dynamic coupling interaction system, the seismic response characteristics of the cast-in-place station structure without a sliding bearing have been analyzed and compared with those of the station structure with the sliding bearing fixed only at the top of the middle columns, and those of the station structure with sliding bearing be fixed between the mid-plate and the sidewall at the same time. The results show that the new isolation station structures suffer fewer earthquake damages at the mid-plate and mid-columns at the same time, which can improve the overall seismic performance of the subway station structure.