Auxiliary Grids Research Articles

Efficient implementation of equivalent medium parameterization in finite-difference seismic wave simulation methods

Summary Grid point discretization of the model has a significant impact on the accuracy of finite-difference seismic waveform simulations. Discretizing the discontinuous velocity model using local point medium parameters can lead to artifact diffraction caused by the stair-step representation and inaccuracies in calculated waveforms due to interface errors, particularly evident when employing coarse grids. To accurately represent model interfaces and reduce interface errors in finite-difference calculations, various equivalent medium parametrization methods have been developed in recent years. Most of these methods require volume-integrated averaging calculations of the medium parameter values within grid cells. The simplest way to achieve this volume averaging is to apply numerical integration averaging to all grid cells. However, this approach demands considerable computational time. To address this computational challenge, we propose employing a set of auxiliary grids to identify which grid cells intersected by the welded interface and perform volume averaging only on these specific cells, thereby reducing unnecessary computational overhead. Additionally, we present a three-dimensional tilted transversely isotropic equivalent medium parameterization method, which effectively suppresses interface errors and artefact diffraction under the application of coarse grids. We also provide an approach for computing the normal direction of the interface, which is essential for the tilted transversely isotropic equivalent medium parameterization. Numerical tests validate the accuracy of the tilted transversely isotropic equivalent medium parameterization method and demonstrate the practicality of the implementation proposed in this paper for complex models.

Semantic Complete Scene Forecasting from a 4D Dynamic Point Cloud Sequence

We study a new problem of semantic complete scene forecasting (SCSF) in this work. Given a 4D dynamic point cloud sequence, our goal is to forecast the complete scene corresponding to the future next frame along with its semantic labels. To tackle this challenging problem, we properly model the synergetic relationship between future forecasting and semantic scene completion through a novel network named SCSFNet. SCSFNet leverages a hybrid geometric representation for high-resolution complete scene forecasting. To leverage multi-frame observation as well as the understanding of scene dynamics to ease the completion task, SCSFNet introduces an attention-based skip connection scheme. To ease the need to model occlusion variations and to better focus on the occluded part, SCSFNet utilizes auxiliary visibility grids to guide the forecasting task. To evaluate the effectiveness of SCSFNet, we conduct experiments on various benchmarks including two large-scale indoor benchmarks we contributed and the outdoor SemanticKITTI benchmark. Extensive experiments show SCSFNet outperforms baseline methods on multiple metrics by a large margin, and also prove the synergy between future forecasting and semantic scene completion.The project page with code is available at scsfnet.github.io.

Analysis of the calculation methods for shell circulating current and grounding current in 500 kV HGIS and the optimization of grounding scheme

The strong electromagnetic coupling between the wire and the shell of Hybrid Gas Insulated Switchgear (HGIS) will lead to a large induced current on the shell, which will further cause local heating of the outside shell. Moreover, the utilization of grounding wire to release this current may also lead to the problem of its heating and ground potential rise. This study is intended to further promote the calculation accuracy of the shell circulating current and grounding current of the 500 kV HGIS and also propose the proper grounding schemes. The simulation results show that the ground grid must be considered in the calculation model and the phase difference of the conductor current could be ignored when the three-phase spacing exceeded 7.5 m. Furthermore, the influences of shorting bars and ground grid on the simulation results were compared and analyzed, and the optimization of grounding scheme and corresponding suggestions were put forward. Generally, it was not recommended to erect shorting bars and auxiliary ground grids. But if the temperature rise exceeds the limit, the shorting bars can be added at the two ends of a HGIS branch and the auxiliary ground grid should be laid using materials with higher flow capacity.

A grid-assisted 3D printing method for magnetically driven micro soft robot

The magnetically driven micro soft robot is one of the research hotspots in the field of micro robots. A 3D printing method for manufacturing micro soft robots with the help of auxiliary grids is proposed. The shape flexibility of the micro soft robots manufactured by the micro 3D printing method can be improved, and internal structures composed of different materials inside the micro robots can be printed using this method. The technical details of the grid-assisted 3D printing process are introduced in this article. Micro soft robots with specific patterns were 3D printed, and the deformation and movement capabilities of the robots were verified in magnetically driven motion experiments. The experimental results prove that the grid-assisted 3D printing technology can not only manufacture magnetically driven micro soft robots, but also adjust the internal structure and motion performance of the robots.

Research and application of large diameter steel cylinder island construction technology in Hong Kong-Zhuhai-Macao Bridge project

ABSTRACT The large diameter deep inserted steel cylinder is used for the first time in the artificial island project of Hong Kong-Zhuhai-Macao Bridge. Through the ingenious conception of the scheme of large diameter deep inserted steel cylinder, this paper analyzes the working mechanism, design key points and construction key points of large diameter deep inserted steel cylinder structure, and introduces the research and application of large diameter deep inserted steel cylinder technology in foundation pit support engineering of East–West artificial island of Hong Kong Zhuhai Macao Bridge. Two arc-shaped steel plate auxiliary grids are used to connect and seal the adjacent cylinders. The arc steel plate can be manufactured in the factory and transported to the site for construction, which greatly reduces the construction difficulty and material cost. The formed cofferdam Island wall has good water stop performance and good structural stability.

Global Bathymetry and Topography at 15 Arc Sec: SRTM15+

An updated global bathymetry and topography grid is presented using a spatial sampling interval of 15 arc sec. The bathymetry is produced using a combination of shipboard soundings and depths predicted using satellite altimetry. New data consists of >33.6 million multibeam and singlebeam measurements collated by several institutions, namely, the National Geospatial‐Intelligence Agency, Japan Agency for Marine‐Earth Science and Technology, Geoscience Australia, Center for Coastal and Ocean Mapping, and Scripps Institution of Oceanography. New altimetry data consists of 48, 14, and 12 months of retracked range measurements from Cryosat‐2, SARAL/AltiKa, and Jason‐2, respectively. With respect to SRTM15_PLUS (Olson et al.,), the inclusion of these new data results in a ∼1.4‐km improvement in the minimum wavelength recovered for sea surface free‐air gravity anomalies, a small increase in the accuracy of altimetrically derived predicted depths, and a 1.24% increase, from 9.60% to 10.84%, in the total area of ocean floor that is constrained by shipboard soundings at 15‐arc sec resolution. Bathymetric grid cells constrained by satellite altimetry have estimated uncertainties of ±150 m in the deep oceans and ±180 m between coastlines and the continental rise. Onshore, topography data are sourced from previously published digital elevation models, predominately SRTM‐CGIAR V4.1 between 60°N and 60°S. ArcticDEM is used above 60°N, while Reference Elevation Model of Antarctica is used below 62°S. Auxiliary grids illustrating shipboard data coverage, marine free‐air gravity anomalies, and vertical gradient gradients are also provided in common data formats.

Adiabatic fission barriers in superheavy nuclei

Using the microscopic-macroscopic model based on the deformed Woods-Saxon single-particle potential and the Yukawa-plus-exponential macroscopic energy we calculated static fission barriers $B_{f}$ for 1305 heavy and superheavy nuclei $98\leq Z \leq 126$, including even - even, odd - even, even - odd and odd - odd systems. For odd and odd-odd nuclei, adiabatic potential energy surfaces were calculated by a minimization over configurations with one blocked neutron or/and proton on a level from the 10-th below to the 10-th above the Fermi level. The parameters of the model that have been fixed previously by a fit to masses of even-even heavy nuclei were kept unchanged. A search for saddle points has been performed by the "Imaginary Water Flow" method on a basic five-dimensional deformation grid, including triaxiality. Two auxiliary grids were used for checking the effects of the mass asymmetry and hexadecapole non-axiallity. The ground states were found by energy minimization over configurations and deformations. We find that the non-axiallity significantly changes first and second fission barrier in many nuclei. The effect of the mass - asymmetry, known to lower the second, very deformed barriers in actinides, in the heaviest nuclei appears at the less deformed saddles in more than 100 nuclei. It happens for those saddles in which the triaxiallity does not play any role, what suggests a decoupling between effects of the mass-asymmetry and triaxiality. We studied also the influence of the pairing interaction strength on the staggering of $B_f$ for odd- and even-particle numbers. Finally, we provide a comparison of our results with other theoretical fission barrier evaluations and with available experimental estimates.

IE-DDM with a novel multiple-grid p-FFT for analyzing multiscale structures in half space

We present an integral equation domain decomposition method accelerated by a novel multiple-grid precorrected fast Fourier transform (MG-p-FFT) for the efficient analysis of multiscale structures in a half space. Based on the philosophy of DDM, the original computational domain is partitioned into several non-overlapping sub-domains. By employing non-conformal discretizations to each domain boundaries, combined field integral equation with half-space dyadic Green’s function is proposed for each individual sub-domain. Subsequently, the MG-p-FFT with auxiliary Cartesian grids with different size, order, location, and spacing, is adopted in each sub-domain independently to account for the self-interactions. Here, the proposed MG-p-FFT scheme outperforms the existing single-grid p-FFT scheme for multiscale problems by reducing the computational time and memory consumption. The proposed method can also be viewed as an effective preconditioning scheme for multiscale problems in a half space. The validity and advantages of the proposed method are illustrated by several representative numerical examples.

Erosion protection Phytoreinforcement of SCARP steep slopes of the holy virgin’s DITCH

Erosion protection landscaping embedment of steep subsoil slopes is a time-sensitive issue of road construction and planning of recreational area that are often fit on a challenging picturesque terrain unsuitable for site development. The article provides the results of a 4-year experiment on landscaping and plant fixing of up to 4.5 m soil slopes with 1:1 and 2:1 grades; the experiment was carried out by the MGSU on the territory of a convent in the south of the Nizhniy Novgorod region. The site has slopes oriented towards all cardinals. At some places the slopes are bedimmed by trees. All these factors create a wide range of geo-ecological conditions for lawns. All the slopes are fixed with geo-fibrefill grids; slopes with 2:1 grade are strengthened by auxiliary grids made of reinforced metal bars, anchors and braces on the bottom of the Holy Moat. The paper recommends composition of grass plants as well as techniques to build up lawns suitable for various micro-climate conditions. It also advises the structure of multi-tier plant entity. The suggested methods are tested during a 3-year maintenance of slopes built for constant use.

Power particles

This paper introduces a new particle-based approach to incompressible fluid simulation. We depart from previous Lagrangian methods by considering fluid particles no longer purely as material points, but also as volumetric parcels that partition the fluid domain. The fluid motion is described as a time series of well-shaped power diagrams (hence the name power particles ), offering evenly spaced particles and accurate pressure computations. As a result, we circumvent the typical excess damping arising from kernel-based evaluations of internal forces or density without having recourse to auxiliary Eulerian grids. The versatility of our solver is demonstrated by the simulation of multiphase flows and free surfaces.

Ferroresonance Phenomenon in Power Plant Auxiliary Grids

The paper discusses the problem of the ferroresonance phenomenon occurring in a medium voltage grid with an ungrounded neutral point of the supply transformer. Results of a study of ferroresonance in the auxiliary grid at Żarnowiec Pumped Storage Power Plant are presented. Analysis of the study results has led to the development and implementation of an effective method of preventing the adverse ferroresonance phenomenon. DOI: 10.12736/issn.2300-3022.2015105

Surface relief and domain structure of ferromagnetic shape memory alloys

A study is made of the surface corrugation during thermal cycling of ferromagnetic shape memory alloys (FSMA). This specific feature is a property of FSMA alloys and is a consequence of martensitic phase transformations not necessarily connected with surface defects of the material. The surface relief structure was studied together with martensite and magnetic domain structure changes during thermal cycling of the samples with the aid of differential polarized light microscopy. The analysis was facilitated making use of auxiliary reference grids applied to the surface of the samples.

An auxilliary grid method for the calculation of electrostatic terms in density functional theory on a real-space grid.

In this work we show the implementation of a linear scaling algorithm for the calculation of the Poisson integral. We use domain decomposition and non-uniform auxiliary grids (NGs) to calculate the electrostatic interaction. We demonstrate the approach within the PARSEC density functional theory code and perform calculations of long 1D carbon chains and other long molecules. Finally, we discuss possible applications to additional problems and geometries.

Robust and efficient overset grid assembly for partitioned unstructured meshes

This paper presents a method to perform efficient and automated Overset Grid Assembly (OGA) on a system of overlapping unstructured meshes in a parallel computing environment where all meshes are partitioned into multiple mesh-blocks and processed on multiple cores. The main task of the overset grid assembler is to identify, in parallel, among all points in the overlapping mesh system, at which points the flow solution should be computed (field points), interpolated (receptor points), or ignored (hole points). Point containment search or donor search, an algorithm to efficiently determine the cell that contains a given point, is the core procedure necessary for accomplishing this task. Donor search is particularly challenging for partitioned unstructured meshes because of the complex irregular boundaries that are often created during partitioning.Another challenge arises because of the large variation in the type of mesh-block overlap and the resulting large load imbalance on multiple processors. Desirable traits for the grid assembly method are efficiency (requiring only a small fraction of the solver time), robustness (correct identification of all point types), and full automation (no user input required other than the mesh system). Additionally, the method should be scalable, which is an important challenge due to the inherent load imbalance. This paper describes a fully-automated grid assembly method, which can use two different donor search algorithms. One is based on the use of auxiliary grids and Exact Inverse Maps (EIM), and the other is based on the use of Alternating Digital Trees (ADT). The EIM method is demonstrated to be more efficient than the ADT method, while retaining robustness. An adaptive load re-balance algorithm is also designed and implemented, which considerably improves the scalability of the method.

WAVE PROPAGATION MODELING IN HIGHLY HETEROGENEOUS MEDIA BY A POLY-GRID CHEBYSHEV SPECTRAL ELEMENT METHOD

A wide range of applications requires the modeling of wave propagation phenomena in media with variable physical properties in the domain of interest, while highly accurate algorithms are needed to avoid unphysical effects. Spectral element methods (SEM), based on either a Chebyshev or a Legendre polynomial basis, have excellent properties of accuracy and flexibility in describing complex models, outperforming other techniques. In the standard SEM approach the computational domain is discretized by using very coarse meshes and constant-property elements, but in some cases the accuracy and the computational efficiency may be seriously reduced. For instance, a finely heterogeneous medium requires grid resolution down to the finest scales, leading to an extremely large problem dimension. In such problems the wavelength scale of interest is much larger but cannot be exploited in order to reduce the problem size. A poly-grid Chebyshev spectral element method (PG-CSEM) can overcome this limitation. In order to accurately deal with continuous variation in the properties, or even with small scale fluctuations, temporary auxiliary grids are introduced which avoid the need of using any finer global grid, and at the macroscopic level the wave field propagation is solved maintaining the SEM accuracy and computational efficiency.

A numerical approach for modeling near-fault ground motion and its application in the 1994 Northridge earthquake

An approach for simulating near-fault ground motion was presented by combining the finite fault model with a numerical algorithm, named investigated lump method presented in this paper for wave propagation. The investigated lumps are constructed from the auxiliary quadrilateral grids. The dynamic equilibrium equations of a typical investigated lump have been derived and obtained by integrating the stresses along the contour of the investigated lump. The stresses are calculated using the constitutive relations and the interpolation techniques. The investigated lump method is then implemented using the equilibrium equations of investigated lumps and the calculations of stresses alternately in time domain. The stability criterion of the algorithm has been given. Comparisons with the discrete wave-number method solutions for predicting the ground motions at the Pacoima Dam during the San Fernando earthquake show the validity of the method presented in this paper for simulating near-field ground motions. A finite fault source model has been implemented in the algorithm here. The source parameters given by Wald et al. (1996) [18] are applied to synthesize the ground motions at three stations during the 1994 Northridge earthquake. The simulating results qualitatively match to the corresponding ground motion records. The studies demonstrated that the approach presented in this paper is an effective tool for the numerical simulation of near-fault ground motion.

A Multiple-Grid Adaptive Integral Method for Multi-Region Problems

A multiple-grid extension of the adaptive integral method (AIM) is presented for fast analysis of scattering from piecewise homogeneous structures. The proposed scheme accelerates the iterative method-of-moments solution of the pertinent surface integral equations by employing multiple auxiliary Cartesian grids: If the structure of interest is composed of <i xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">K</i> homogeneous regions, it introduces <i xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">K</i> different auxiliary grids. It uses the <i xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">k</i> <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">th</sup> auxiliary grid first to determine near-zones for the basis functions and then to execute AIM projection, propagation, interpolation, and near-zone pre-correction stages in the <i xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">k</i> <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">th</sup> region. Thus, the AIM stages are executed a total of <i xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">K</i> times using different grids and different groups of basis functions. The proposed multiple-grid AIM scheme requires a total of <i xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">O</i> ( <i xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">N</i> <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">nz,near</sup> +? <sub xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">k</sub> <i xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">N</i> <sub xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">k</sub> <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">C</sup> log <i xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">N</i> <sub xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">k</sub> <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">C</sup> ) operations per iteration, where <i xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">N</i> <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">nz,near</sup> denotes the total number of near-zone interactions in all regions and <i xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">N</i> <sub xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">k</sub> <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">C</sup> denotes the number of nodes of the <i xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">k</i> <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">th</sup> Cartesian grid. Numerical results validate the method's accuracy and reduced complexity for large-scale canonical structures with large numbers of regions (up to ~10 <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">6</sup> degrees of freedom and ~10 <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">3</sup> regions). Moreover, an investigation of HF-band wave propagation in a loblolly pine forest model demonstrates the method's generality and practical applicability. Multiple-grid AIM accelerated simulations with various tree models show that higher fidelity models for the trunk material and branch geometry are needed for accurate calculation of horizontally-polarized field propagation while lower fidelity models can be satisfactory for analyzing vertically-polarized field propagation.

Triangular grid method for stress-wave propagation in 2-D orthotropic materials

A triangular grid method is presented to calculate propagation problems of elastic stress waves in 2-D orthotropic materials. This method is based on the dynamic equilibrium equations of the computational cells formed among the auxiliary triangular grids. The solution is obtained by calculating alternately the nodal displacements and the central point stresses of the spatial grids. The numerical results are compared with the corresponding solutions of the finite element method. Comparisons show that the triangular grid method yields a higher calculational speed than the finite element method. The stress concentrations are investigated from wave-field analyses when the stress wave propagates within an orthotropic plate with a hole. Finally, the presented numerical method is used to study the features of wave propagation and diffraction in a square orthotropic plate with a hole when an impact load is applied to the top of the plate.

Unstructured grid method for stress wave propagation in elastic media

A new numerical method, named unstructured grid method, is presented to calculate two-dimensional stress wave propagation problems in elastic media. This method is based on the dynamic equilibrium equations of the investigated lumps that are formed among the auxiliary unstructured grids. The approach to the solution is obtained by calculating the node displacements and the central point stresses of the spatial grids alternately. The stability criterion of the method is given and tested numerically. Analysis of dispersion for numerical phase velocity is studied. The rule is provided to determine the spatial step size of discretization mesh so as to ensure the numerical accuracy of the unstructured grid method. Numerical results are compared with the available exact solutions and with the solutions of central difference method for demonstrating the good accuracy of the new method. The unstructured grid method possesses higher calculating speed than the central difference method. Finally, the dynamic behavior of a square plate with a hole under an impact loading is researched.

Protection of 6-kV TPP Auxiliary Grids with Two Neutral Ground Modes from Single-Phase Ground Faults

A 6-kV power plant auxiliary grid with low ground fault currents is considered, which can operate in modes with insulated or compensated neutral depending on the connection of the working or backup auxiliary supply element. A high-sensitivity zero-phase sequence current protection from single-phase ground faults is suggested, which operates on commercial-frequency capacitive ground fault current and is suitable for the two kinds of auxiliary grid.