Transient electromagnetic field computations for polygonal loops on layered earths

Geophysics ◽  
1987 ◽  
Vol 52 (6) ◽  
pp. 785-793 ◽  
Author(s):  
A. P. Raiche

The transient electromagnetic response (vertical and horizontal components of dB/dt) of a large polygonal transmitting loop on a layered earth is calculated using a nested interpolation scheme based on the dipole‐dipole response function. The frequency‐domain field of a vertical magnetic dipole is inverse Laplace transformed into the time domain, for selected values of the Hankel transform variable, using the Gaver‐Stehfest method. After interpolation, the result is inverse Hankel transformed (for selected values of distance) using digital filters. Interpolating over distance allows integration of the dipole response over the area of one or more transmitting loops. An interpolation over time gives the step‐current response, which in turn is convolved with the transmitter‐receiver characteristics to yield the system response. This method allows robust calculation of several transmitting loops (with different signal parameters) and several receiver positions in little more time than that required for one loop with one receiver. The computation time for the single transmitter‐receiver response can be decreased by analytical integration of the Bessel function over the transmitter loop area before performing the inverse Hankel transform. Since this procedure precludes the use of standard digital filters for the inverse Hankel transform, it is not efficient for multireceiver computations.

Geophysics ◽  
1986 ◽  
Vol 51 (7) ◽  
pp. 1450-1461 ◽  
Author(s):  
Y. Goldman ◽  
C. Hubans ◽  
S. Nicoletis ◽  
S. Spitz

We present a numerical method for solving Maxwell’s equations in the case of an arbitrary two‐dimensional resistivity distribution excited by an infinite current line. The electric field is computed directly in the time domain. The computations are carried out in the lower half‐space only because exact boundary conditions are used on the free surface. The algorithm follows the finite‐element approach, which leads (after space discretization) to an equation system with a sparse matrix. Time stepping is done with an implicit time scheme. At each time step, the solution of the equation system is provided by the fast system ICCG(0). The resulting algorithm produces good results even when large resistivity contrasts are involved. We present a test of the algorithm’s performance in the case of a homogeneous earth. With a reasonable grid, the relative error with respect to the analytical solution does not exceed 1 percent, even 2 s after the source is turned off.


2020 ◽  
Vol 10 (3) ◽  
pp. 1027 ◽  
Author(s):  
Yanju Ji ◽  
Xiangdong Meng ◽  
Weimin Huang ◽  
Yanqi Wu ◽  
Gang Li

The grounded electrical-source airborne transient electromagnetic (GREATEM) system is widely used in mineral exploration. Meanwhile, the induced polarization (IP) effect, which indicates the polarizability of the earth, is often found. In this paper, the Maxwell equations in the frequency domain are transformed into fictitious wave domain, where Maxwell equations are solved by the time domain finite difference method. Then, an integral transformation method is used to convert the calculation results back to the time domain. A three-dimensional (3D) numerical simulation in a polarizable medium is presented. The accuracy of this method is proven by comparing it with the analytical solution and the existing method, and the calculation efficiency is increased five-fold. The simulation results show that the GREATEM system has a higher response amplitude in the conductive region, while IP effects cannot be identified in the conductive area. The GREATEM system has a higher response amplitude in the low-resistance region, but IP effects cannot be identified in the low-resistance area, and the detection of IP effects is more suitable for the high-resistance area. Therefore, it is necessary to improve the detection ability of the GREATEM system in the low-resistance area.


Author(s):  
Niels Hørbye Christiansen ◽  
Per Erlend Torbergsen Voie ◽  
Jan Høgsberg ◽  
Nils Sødahl

Dynamic analyses of slender marine structures are computationally expensive. Recently it has been shown how a hybrid method which combines FEM models and artificial neural networks (ANN) can be used to reduce the computation time spend on the time domain simulations associated with fatigue analysis of mooring lines by two orders of magnitude. The present study shows how an ANN trained to perform nonlinear dynamic response simulation can be optimized using a method known as optimal brain damage (OBD) and thereby be used to rank the importance of all analysis input. Both the training and the optimization of the ANN are based on one short time domain simulation sequence generated by a FEM model of the structure. This means that it is possible to evaluate the importance of input parameters based on this single simulation only. The method is tested on a numerical model of mooring lines on a floating off-shore installation. It is shown that it is possible to estimate the cost of ignoring one or more input variables in an analysis.


Geophysics ◽  
1983 ◽  
Vol 48 (7) ◽  
pp. 934-952 ◽  
Author(s):  
P. Weidelt

An exact solution is given for the electromagnetic induction in a dipping dike of finite conductivity, represented as a thin half‐sheet in a nonconducting surrounding. The problem is formulated for arbitrary dipole or circular loop [Formula: see text] configurations. The formal solution obtained by the Wiener‐Hopf technique is cast into a rapidly convergent triple integral suitable for an effective numerical treatment. A good agreement is found between numerical results and analog measurements available for harmonic excitation. The transient response is obtained as a superposition of the half‐sheet free‐decay modes and is illustrated by some numerical examples for coincident loops, including a diagram for the approximate determination of conductance and depth of a vertical dike.


2017 ◽  
Vol 9 (6) ◽  
Author(s):  
Stephen L. Canfield ◽  
Reabetswe M. Nkhumise

This paper develops an approach to evaluate a state-space controller design for mobile manipulators using a geometric representation of the system response in tool space. The method evaluates the robot system dynamics with a control scheme and the resulting response is called the controllability ellipsoid (CE), a tool space representation of the system’s motion response given a unit input. The CE can be compared with a corresponding geometric representation of the required motion task (called the motion polyhedron) and evaluated using a quantitative measure of the degree to which the task is satisfied. The traditional control design approach views the system response in the time domain. Alternatively, the proposed CE views the system response in the domain of the input variables. In order to complete the task, the CE must fully contain the motion polyhedron. The optimal robot arrangement would minimize the total area of the CE while fully containing the motion polyhedron. This is comparable to minimizing the power requirements of robot design when applying a uniform scale to all inputs. It will be shown that changing the control parameters changes the eccentricity and orientation of the CE, implying a preferred set of control parameters to minimize the design motor power. When viewed in the time domain, the control parameters can be selected to achieve desired stability and time response. When coupled with existing control design methods, the CE approach can yield robot designs that are stable, responsive, and minimize the input power requirements.


Author(s):  
Krishnakumar Gopalakrishnan ◽  
Teng Zhang ◽  
Gregory J. Offer

Research into reduced-order models (ROM) for Lithium-ion batteries is motivated by the need for a real-time embedded model possessing the accuracy of physics-based models, while retaining computational simplicity comparable to equivalent-circuit models. The discrete-time realization algorithm (DRA) proposed by Lee et al. (2012, “One-Dimensional Physics-Based Reduced-Order Model of Lithium-Ion Dynamics,” J. Power Sources, 220, pp. 430–448) can be used to obtain a physics-based ROM in standard state-space form, the time-domain simulation of which yields the evolution of all the electrochemical variables of the standard pseudo-2D porous-electrode battery model. An unresolved issue with this approach is the high computation requirement associated with the DRA, which needs to be repeated across multiple SoC and temperatures. In this paper, we analyze the computational bottleneck in the existing DRA and propose an improved scheme. Our analysis of the existing DRA reveals that singular value decomposition (SVD) of the large Block–Hankel matrix formed by the system's Markov parameters is a key inefficient step. A streamlined DRA approach that bypasses the redundant Block–Hankel matrix formation is presented as a drop-in replacement. Comparisons with existing DRA scheme highlight the significant reduction in computation time and memory usage brought about by the new method. Improved modeling accuracy afforded by our proposed scheme when deployed in a resource-constrained computing environment is also demonstrated.


2021 ◽  
pp. 1-30
Author(s):  
Mohsen Bayani ◽  
Casper Wickman ◽  
Aswin Dhananjai Krishnaswamy ◽  
Chidambaram Sathappan ◽  
Rikard Söderberg

Abstract Avoiding quality problems in passenger cars, such as squeak and rattle (S&R), has been a remarkable cost-saving consideration. The introduction of electric engines and less engaged drivers due to autonomous driving is expected to further stress the need for quieter cabins. However, the complexity of the virtual evaluation of S&R events has obstructed the practical treatment of these quality issues in the pre-design-freeze phases of product development. In this study, new quantified frequency-domain metrics are proposed to measure the risk for the generation of S&R in subsystem assemblies. The proposed metrics measure the resonance risk and the mode shape similarity in the critical interfaces for S&R. The calculations are done based on the system response in the frequency domain. Compared to the time-domain evaluation methods, the knowledge about the system excitation levels is not essential and the calculations are more time-efficient. The proposed metrics can be used in closed-loop design optimisation processes to involve S&R attributes in the pre-design-freeze attribute trade-off activities besides other attributes. In this work, these metrics were used in a two-stage optimisation problem to optimise the connection configuration in two industrial cases. As compared to the baseline design, the risk for S&R was reduced by improving the system behaviour in terms of resonance risk and mode shape similarity. This was achieved by applying some adjustments to the location of the fasteners while maintaining the same general connection configuration concept.


2018 ◽  
Vol 2018 ◽  
pp. 1-16 ◽  
Author(s):  
Tian Hao Huang ◽  
Kuo Lung Lian

In response to the growing demand for medium- and high-power trends, multilevel voltage source converters (VSCs) have been attracting growing considerations. One of the widely used VSCs are the diode-clamped multilevel VSC (DCM-VSC). As these converters proliferate, their harmonic impact may become significant. Nevertheless, a harmonic model for the DCM-VSC is currently lacking in the literature. In this paper, the ABCD matrix, mapping the input harmonics to the output harmonics of DCM-VSC, is derived. As the matrix is formulated in the time-domain, the output harmonics are exact and do not suffer from harmonic truncation errors. As the paper will demonstrate, the derived ABCD matrix can be easily applied to a microgrid system and users can easily predict all the uncharacteristic harmonics when a microgrid is subjected to various conditions of imbalance. In addition to all the results being validated with those of PSCAD/EMTDC, the computation time of the proposed method is in contrast much shorter.


Sign in / Sign up

Export Citation Format

Share Document