Search results

In this paper, a modified meshless method, as one of the numerical techniques that has recently emerged in the area of computational electromagnetics, is extended to solving time-domain wave equation. The paper aims to discuss these issues.

In space domain, the fields at the collocation points are expanded into a series of new Shepard's functions which have been suggested recently and are treated with a meshless method procedure. For time discretization of the second-order time-derivative, two finite-difference schemes, i.e. backward difference and Newmark-β techniques, are proposed.

Both schemes are implicit and always stable and have unconditional stability with different orders of accuracy and numerical dispersion. The unconditional stability of the proposed methods is analytically proven and numerically verified. Moreover, two numerical examples for electromagnetic field computation are also presented to investigate characteristics of the proposed methods.

The paper presents two unconditionally stable schemes for meshless methods in time-domain electromagnetic problems.

The FDTD approach is rapidly becoming one of the most widely used computational methods in electromagnetics. It is a marching‐in‐time procedure which simulates the continuous actual waves by sampled data numerical analogues propagating in a data space stored in a computer. Thus it leads to a complete understanding of near fields and transient effects. In this paper, the study of the propagating modes and the calculation of the energy forced in the interior of a terminated waveguide is performed with an iterative FDTD technique. Results are obtained and compared with the theoretical ones. Their agreement is almost perfect.

The purpose of this paper is to report the implementation of an alternative time integration procedure for the dynamic non-linear analysis of structures.

The time integration algorithm discussed in this work corresponds to a spectral decomposition technique implemented in the time domain. As in the case of the modal decomposition in space, the numerical efficiency of the resulting integration scheme depends on the possibility of uncoupling the equations of motion. This is achieved by solving an eigenvalue problem in the time domain that only depends on the approximation basis being implemented. Complete sets of orthogonal Legendre polynomials are used to define the time approximation basis required by the model.

A classical example with known analytical solution is presented to validate the model, in linear and non-linear analysis. The efficiency of the numerical technique is assessed. Comparisons are made with the classical Newmark method applied to the solution of both linear and non-linear dynamics. The mixed time integration technique presents some interesting features making very attractive its application to the analysis of non-linear dynamic systems. It corresponds in essence to a modal decomposition technique implemented in the time domain. As in the case of the modal decomposition in space, the numerical efficiency of the resulting integration scheme depends on the possibility of uncoupling the equations of motion.

One of the main advantages of this technique is the possibility of considering relatively large time step increments which enhances the computational efficiency of the numerical procedure. Due to its characteristics, this method is well suited to parallel processing, one of the features that have to be conveniently explored in the near future.

The finite element method (FEM) is used to calculate the two-dimensional anti-plane dynamic response of structure embedded in D’Alembert viscoelastic multilayered soil on the rigid bedrock. This paper aims to research a time-domain absorbing boundary condition (ABC), which should be imposed on the truncation boundary of the finite domain to represent the dynamic interaction between the truncated infinite domain and the finite domain.

A high-order ABC for scalar wave propagation in the D’Alembert viscoelastic multilayered media is proposed. A new operator separation method and the mode reduction are adopted to construct the time-domain ABC.

The derivation of the ABC is accurate for the single layer but less accurate for the multilayer. To achieve high accuracy, therefore, the distance from the truncation boundary to the region of interest can be zero for the single layer but need to be about 0.5 times of the total layer height of the infinite domain for the multilayer. Both single-layered and multilayered numerical examples verify that the accuracy of the ABC is almost the same for both cases of only using the modal number excited by dynamic load and using the full modal number of infinite domain. Using the ABC with reduced modes can not only reduce the computation cost but also be more friendly to the stability. Numerical examples demonstrate the superior properties of the proposed ABC with stability, high accuracy and remarkable coupling with the FEM.

A high-order time-domain ABC for scalar wave propagation in the D’Alembert viscoelastic multilayered media is proposed. The proposed ABC is suitable for both linear elastic and D’Alembert viscoelastic media, and it can be coupled seamlessly with the FEM. A new operator separation method combining mode reduction is presented with better stability than the existing methods.

This paper aims to present a novel hybrid time integration approach for efficient numerical simulations of multiscale problems involving interactions of electromagnetic fields with fine structures.

The entire computational domain is discretized with a coarse grid and a locally refined subgrid containing the tiny objects. On the coarse grid, the time integration of Maxwell’s equations is realized by the conventional finite-difference technique, while on the subgrid, the unconditionally stable Krylov-subspace-exponential method is adopted to breakthrough the Courant–Friedrichs–Lewy stability condition.

It is shown that in contrast with the conventional finite-difference time-domain method, the proposed approach significantly reduces the memory costs and computation time while providing comparative results.

An efficient hybrid time integration approach for numerical simulations of multiscale electromagnetic problems is presented.

The aim of this paper is the experimental validation of an original time‐domain thin‐shell formulation. The numerical results of a three‐dimensional thin‐shell model are compared with the measurements performed on a heating device at different working frequencies.

A time‐domain extension of the classical frequency‐domain thin‐shell approach is used for the finite‐element analysis of a shielded pulse‐current induction heater. The time‐domain interface conditions at the shell surface are expressed in terms of the average flux density vector in the shell, as well as in terms of a limited number of higher‐order components.

A very good agreement between measurements and simulations is observed. A clear advantage of the proposed thin‐shell approach is that the mesh of the computation domain does not depend on the working frequency anymore. It provides a good compromise between computational cost and accuracy. Indeed, adding a sufficient number of induction components, a very high accuracy can be achieved.

The method is based on the coupling of a time‐domain 1D thin‐shell model with a magnetic vector potential formulation via the surface integral term. A limited number of additional unknowns for the magnetic flux density are incorporated on the shell boundary.

This paper aims to omit the difficulties of directly finding the periodic steady-state solutions for electromagnetic devices described by circuit models.

Determine the discrete integral operator of periodic functions and develop an iterative algorithm determining steady-state solutions by a multiplication of matrices only.

An alternative method to creating finite-difference relations directly determining steady-state solutions in the time domain.

Reduction of software and hardware requirements for determining steady-states of electromagnetic.

A unified approach for directly finding steady-state solutions for ordinary nonlinear differential equations presented in the normal form.

Eliminate the necessity of solving high-order finite-difference equations for steady-state analysis of electromagnetic devices described by circuit models.

The purpose of this paper is to present sensitivity analysis of electromagnetic fields in the time‐domain.

The method utilizing adjoint models is commonly used to evaluate sensitivity. In connection with widely applied finite element method, the time‐stepping scheme for discretization of time functions is used.

The proposed semi‐discrete method allows us to obtain time‐domain solution without time‐stepping. For space discretization, the authors use finite elements, as usual. The semi‐discrete method delivers analytical and continuous solution for any given time of analysis, which has a form of exponential functions. In order to obtain an analytical formula, there is necessary the integration of sensitivity equation. The paper finds possible solutions of this problem, either the application of Zassenhaus formula or improvement of commutation properties of two matrices.

Drawback of this method is matrices which are losing their symmetry and are no more banded. All calculations in this work were carried out with fully assigned matrices. Comparison of the efficiency of the semi‐discrete method with classical method shows that, despite the high demand for memory, this method can compete in relation to finite elements with the time‐stepping.

The resultant gradient information may be used for solving inverse problems, such as optimization of magnetic circuits and identification of material conductivity distributions.

The paper offers compact formula for sensitivity evaluation.

This paper gives a bibliographical review of the finite element and boundary element parallel processing techniques from the theoretical and application points of view. Topics include: theory – domain decomposition/partitioning, load balancing, parallel solvers/algorithms, parallel mesh generation, adaptive methods, and visualization/graphics; applications – structural mechanics problems, dynamic problems, material/geometrical non‐linear problems, contact problems, fracture mechanics, field problems, coupled problems, sensitivity and optimization, and other problems; hardware and software environments – hardware environments, programming techniques, and software development and presentations. The bibliography at the end of this paper contains 850 references to papers, conference proceedings and theses/dissertations dealing with presented subjects that were published between 1996 and 2002.

Recently, the original Wavenumber approach was introduced for dynamic analysis of dam-reservoir systems in frequency domain in the context of pure finite element programming. But its main disadvantages are that it cannot be implemented in time domain. The purpose of this paper is to propose an approximation to the original approach which enables one to carry out this effective method in time domain as well as in frequency domain. Based on the present investigation, it is proven that the Approximate Wavenumber approach has inherent characteristics, which allows it to be envisaged as an effective technique for calculating the response of concrete gravity dam-reservoir systems in time domain.

The method is described initially. Subsequently, the response of an idealized triangular dam-reservoir system is obtained by the proposed approach as well as by applying two other well-known absorbing conditions which are widely utilized in practice. The results are also controlled against the corresponding exact responses. It should be emphasized that all results presented herein are obtained by the FE-FE method under different absorbing conditions applied on the truncation boundary. These include two well-known absorbing conditions referred to as Sommerfeld and Sharan as well as the proposed approach of the present study (i.e. Approximate Wavenumber condition).

It is concluded that the maximum error for the Approximate Wavenumber approach is in the range of 10 percent at the major peaks of the response. This occurs mainly for the very low reservoir lengths under full reflective reservoir base condition and vertical excitation. This is a remarkable result for any kind of robust truncation boundary simulation that one may expect. The fundamental frequency of the system is captured correctly for the Approximate Wavenumber approach, even in cases of low reservoir length.

Based on this investigation, it is proven that the Approximate Wavenumber approach has inherent characteristics, which allows it to be envisaged as an effective technique for calculating the response of concrete gravity dam-reservoir systems in time domain. It is concluded that the maximum error for the Approximate Wavenumber approach is in the range of 10 percent at the major peaks of the response. This occurs mainly for the very low reservoir lengths under full reflective reservoir base condition and vertical excitation. This is a remarkable result for any kind of robust truncation boundary simulation that one may expect.