On the use of a meshless solver for PDEs governing electromagnetic transients
Introduction
In the last two decades the meshless methods have known a great success in the simulation of a wide problems variety as a valid computational alternative to grid methods. They share common features such as the avoidance of the use of grids, but are different in functions approximation and computational processes.
The numerical technique known as Smoothed Particle Hydrodynamics (SPH) [4], [5], [6], [10], [11] is a meshless method and its attractiveness and popularity is due to the evaluation of unknown field functions and relative differential operators by means of an integral representation based on a suitable interpolating function. The integral representation is discretized by using a set of particles scattered in the problem domain.
The appropriate choice of the smoothing kernel function is a crucial task before performing any computation using the SPH solver. The smoothing kernel function is of remarkable importance since it not only determines the interpolating pattern, but also defines the width of the influence area of a particle determined by a parameter h called as smoothing length. The choice of this parameter is a key variable for the kernel’s worth. In fact, the smoothing kernel function should have a certain degree of consistency which can be expressed by its ability to reproduce the polynomials in both integral and discrete formulations [2], [3], [7]. For each smoothing kernel function only a set of h values, related to the interspacing particles, verifying the polynomial reproducing conditions must be taken into account. In this paper an analysis of the smoothing length h values is carried out by adopting two bell-shaped smoothing kernel functions widely used in literature [7], [9], [11]. Namely, Gaussian and cubic B-spline smoothing kernel functions are considered in one and two dimensions. Validation tests are performed by considering the SPH method suitably reformulated for solving the partial differential equations (PDEs) governing electromagnetic transients [1]. The particle expressions of the Maxwell’s curl equations are provided in one and two dimensions in free space. By working with curl equations the derivatives of the unknown field functions components must verify the consistency conditions. Various simulations are reported with different values of the smoothing length h by considering the transient propagation of a time and space variable pulse. The 2-D model is proposed for a transverse electric wave.
The paper is organized as follows. In Section 2 the background of SPH method is assumed and an analysis of the fundamental issues is reported; namely, the discrete constant and linear consistency conditions are investigated and a set of h values is determined. In Section 3 the meshless formulation of the Maxwell’s curl equations is provided in one and two dimensions. Simulation results referred to canonical case studies are reported.
Section snippets
Basic issues
In order to approximate a sufficiently regular function f (x) in a domain Ω ⊆ Rd by means of a convolution function fh ≡ f∗W the numerical technique known as SPH adopts the so-called kernel approximation [6], [10], [11]:In (1) the function W is the smoothing kernel function depending on the spatial variables and on the smoothing length parameter h:where R = ∥x-y∥/h and αd is a dimension-dependent normalization constant. The smoothing kernel function is assumed
Numerical investigations
In this section, the SPH methodology is used to numerically achieve the time-domain electric and magnetic fields. In this context a “particle” means an electromagnetic field point.
In order to better clarify the main features of the method applied to electromagnetic phenomena, let us consider the time-dependent Maxwell’s curl equations in free space:where E and H are the electric and magnetic vector fields, ε0 is the vacuum permittivity and μ0 is the vacuum permeability.
Conclusions
This paper provides a method, based on the polynomial reproducing conditions, to determine a set of values of the smoothing length h to achieve a good SPH approximation of the unknown field functions. Two bell-shaped smoothing kernel functions, i.e. the Gaussian and the cubic B-spline smoothing kernel functions, have been taken into account. At first, the criterion has been adopted to recognize proper h values verifying the constant and linear consistency conditions. Moreover, the consistency
Acknowledgement
This work has been supported by the Italian Ministero dell’Università e della Ricerca.
References (13)
- et al.
Variational and momentum preservation aspects of smooth particle hydrodynamics formulations
Computer Methods in Applied Mechanics and Engineering
(1999) - et al.
Constructing smoothing functions in smoothed particle hydrodynamics with applications
Journal of Computational and Applied Mathematics
(2003) An introduction to SPH
Computer Physics Communications
(1988)- et al.
A smoothed particle interpolation scheme for transient electromagnetic simulation
IEEE Transactions on Magnetics
(2006) - et al.
On the completeness of meshfree methods
International Journal for Numerical Methods in Engineering
(1998) - et al.
Smoothed particle hydrodynamics: theory and application to nonspherical stars
Monthly Notices of the Royal Astronomical Society
(1985)
Cited by (11)
Numerical solution of EFIE using MLPG methods
2017, Engineering Analysis with Boundary ElementsCitation Excerpt :Nevertheless, these solvers are very flexible and can properly handle complex media. Currently, two types of MFMs corresponding to such solvers are developed; i.e., time-domain PDE solvers [5–17] that are mostly based on the meshfree collocation methods (MCM) [18] and frequency-domain PDE solvers [19–24] exploiting meshless local Petrov-Galerkin (MLPG) methods [25]. These types are meshfree counterparts of, respectively, FDTD and FEM.
Numerical simulation of electric field in complex geometries for different electrode arrangements using meshless local MQ-DQ method
2011, Journal of ElectrostaticsCitation Excerpt :They compared their results with analytical methods and gained accurate results. Francomano et al. [26] developed the previous work using SPH methodology to analyze electromagnetic transients. In their work, the time-dependent Maxwell’s curl equations in free space for 1-D and 2-D formulations was considered and it was focused on the interpolating smoothing kernel function which strongly influences the computational results.
Meshfree analysis of high-frequency field-effect transistors: distributed modeling approach
2019, Journal of Computational ElectronicsAnomalous performance enhancement effects in Ruthenium-based Dye Sensitized Solar Cells
2017, 2017 6th International Conference on Clean Electrical Power: Renewable Energy Resources Impact, ICCEP 2017Renewable energy evaluation in Italy and Brazil
2017, 2017 6th International Conference on Clean Electrical Power: Renewable Energy Resources Impact, ICCEP 2017Numerical analysis of light soaking phenomenon in Ruthenium based Dye Sensitized Solar Cells
2017, 2017 6th International Conference on Clean Electrical Power: Renewable Energy Resources Impact, ICCEP 2017