Finite element solution based on fast numerical technique for large-scale electromagnetic computation

A numerical technique of the target-region locating （TRL） solver in conjunction with the wave-front method is presented for the application of the finite element method （FEM） for 3-D electromagnetic computation. First, the principle of TRL technique is described. Then, the availability of TRL solver for nonlinear application is particularly discussed demonstrating that this solver can be easily used while still remaining great efficiency. The implementation on how to apply this technique in FEM based on magnetic vector potential （MVP） is also introduced. Finally, a numerical example of 3-D magnetostatic modeling using the TRL solver and FEMLAB is given. It shows that a huge computer resource can be saved by employing the new solver.

DSMC: Fast direct simulation Monte Carlo solver for the Boltzmann equation by Multi-Chain Markov Chain and multicore programming

Direct Simulation Monte Carlo(DSMC)solves the Boltzmann equation with large Knudsen number.The Boltzmann equation generally consists of three terms:the force term,the diffusion term and the collision term.While the first two terms of the Boltzmann equation can be discretized by numerical methods such as the finite volume method,the third term can be approximated by DSMC,and DSMC simulates the physical behaviors of gas molecules.However,because of the low sampling efficiency of Monte Carlo Simulation in DSMC,this part usually occupies large portion of computational costs to solve the Boltzmann equation.In this paper,by Markov Chain Monte Carlo(MCMC)and multicore programming,we develop Direct Simulation Multi-Chain Markov Chain Monte Carlo(DSMC3):a fast solver to calculate the numerical solution for the Boltzmann equation.Computational results show that DSMC3 is significantly faster than the conventional method DSMC.

A Fast Direct Solver for a Class of 3-D Elliptic Partial Differential Equation with Variable Coefficient

We propose a direct solver for the three-dimensional Poisson equation with a variable coefficient,and an algorithm to directly solve the associated sparse linear systems that exploits the sparsity pattern of the coefficient matrix.Introducing some appropriate finite difference operators,we derive a second-order scheme for the solver,and then two suitable high-order compact schemes are also discussed.For a cube containing N nodes,the solver requires O(N^(3/2)log^(2)N)arithmetic operations and O(NlogN)memory to store the necessary information.Its efficiency is illustrated with examples,and the numerical results are analysed.

A Fast High Order Iterative Solver for the Electromagnetic Scattering by Open Cavities Filled with the Inhomogeneous Media

The scattering of the open cavity filled with the inhomogeneous media is studied.The problem is discretized with a fourth order finite difference scheme and the immersed interfacemethod,resulting in a linear system of equations with the high order accurate solutions in the whole computational domain.To solve the system of equations,we design an efficient iterative solver,which is based on the fast Fourier transformation,and provides an ideal preconditioner for Krylov subspace method.Numerical experiments demonstrate the capability of the proposed fast high order iterative solver.

Nonconforming finite element Stokes complexes in three dimensions

Two nonconforming finite element Stokes complexes starting from the conforming Lagrange element and ending with the nonconforming P_1-P_0 element for the Stokes equation in three dimensions are studied.Commutative diagrams are also shown by combining nonconforming finite element Stokes complexes and interpolation operators.The lower order H(gradcurl)-nonconforming finite element only has 14 degrees of freedom,whose basis functions are explicitly given in terms of the barycentric coordinates.The H(gradcurl)-nonconforming elements are applied to solve the quad-curl problem,and the optimal convergence is derived.By the nonconforming finite element Stokes complexes,the mixed finite element methods of the quad-curl problem are decoupled into two mixed methods of the Maxwell equation and the nonconforming P_1-P_0 element method for the Stokes equation,based on which a fast solver is discussed.Numerical results are provided to verify the theoretical convergence rates.

An Efficient Neural-Network and Finite-Difference Hybrid Method for Elliptic Interface Problems with Applications

A new and efficient neural-network and finite-difference hybrid method is developed for solving Poisson equation in a regular domain with jump discontinuities on embedded irregular interfaces.Since the solution has low regularity across the interface,when applying finite difference discretization to this problem,an additional treatment accounting for the jump discontinuities must be employed.Here,we aim to elevate such an extra effort to ease our implementation by machine learning methodology.The key idea is to decompose the solution into singular and regular parts.The neural network learning machinery incorporating the given jump conditions finds the singular solution,while the standard five-point Laplacian discretization is used to obtain the regular solution with associated boundary conditions.Regardless of the interface geometry,these two tasks only require supervised learning for function approximation and a fast direct solver for Poisson equation,making the hybrid method easy to implement and efficient.The two-and three-dimensional numerical results show that the present hybrid method preserves second-order accuracy for the solution and its derivatives,and it is comparable with the traditional immersed interface method in the literature.As an application,we solve the Stokes equations with singular forces to demonstrate the robustness of the present method.

An Efficient Cartesian Grid-Based Method for Scattering Problems with Inhomogeneous Media

Boundary integral equations provide a powerful tool for the solution of scattering problems.However,often a singular kernel arises,in which case the standard quadratures will give rise to unavoidable deteriorations in numerical precision,thus special treatment is needed to handle the singular behavior.Especially,for inhomogeneous media,it is difficult if not impossible to find out an analytical expression for Green’s function.In this paper,an efficient fourth-order accurate Cartesian grid-based method is proposed for the two-dimensional Helmholtz scattering and transmission problems with inhomogeneous media.This method provides an alternative approach to indirect integral evaluation by solving equivalent interface problems on Cartesian grid with a modified fourth-order accurate compact finite difference scheme and a fast Fourier transform preconditioned conjugate gradient(FFT-PCG)solver.A remarkable point of this method is that there is no need to know analytical expressions for Green’s function.Numerical experiments are provided to demonstrate the advantage of the current approach,including its simplicity in implementation,its high accuracy and efficiency.

An Indirect-Forcing Immersed Boundary Method for Incompressible Viscous Flows with Interfaces on Irregular Domains

An indirect-forcing immersed boundary method for solving the incompressible Navier-Stokes equations involving the interfaces and irregular domains is developed.The rigid boundaries and interfaces are represented by a number of Lagrangian control points.Stationary rigid boundaries are embedded in the Cartesian grid and singular forces at the rigid boundaries are applied to impose the prescribed velocity conditions.The singular forces at the interfaces and the rigid boundaries are then distributed to the nearby Cartesian grid points using the immersed boundary method.In the present work,the singular forces at the rigid boundaries are computed implicitly by solving a small system of equations at each time step to ensure that the prescribed velocity condition at the rigid boundary is satisfied exactly.For deformable interfaces,the forces that the interface exerts on the fluid are computed from the configuration of the elastic interface and are applied to the fluid.The Navier-Stokes equations are discretized using finite difference method on a staggered uniform Cartesian grid by a second order accurate projection method.The ability of the method to simulate viscous flows with interfaces on irregular domains is demonstrated by applying to the rotational flow problem,the relaxation of an elastic membrane and flow in a constriction with an immersed elastic membrane.

A direct solver with O（N） complexity for integral equations on one-dimensional domains

An algorithm for the direct inversion of the linear systems arising from NystrSm discretization of integral equations on one-dimensional domains is described. The method typically has O（N） complexity when applied to boundary integral equations （BIEs） in the plane with non-oscillatory kernels such as those associated with the Laplace and Stokes＇ equations. The scaling coefficient suppressed by the ＂big-O＂ notation depends logarithraically on the requested accuracy. The method can also be applied to BIEs with oscillatory kernels such as those associated with the Helmholtz and time-harmonic Maxwell equations; it is efficient at long and intermediate wave-lengths, but will eventually become prohibitively slow as the wave-length decreases. To achieve linear complexity, rank： deficiencies in the off-diagonal blocks of the coefficient matrix are exploited. The technique is conceptually related to the H- and H2-matrix arithmetic of Hackbusch and coworkers, and is closely related to previous work on Hierarchically Semi-Separable matrices.

An Efficient Numerical Solution Method for Elliptic Problems in Divergence Form

In this paper the problem−div(a(x,y)∇u)=f with Dirichlet boundary conditions on a square is solved iteratively with high accuracy for u and∇u using a new scheme called“hermitian box-scheme”.The design of the scheme is based on a“hermitian box”,combining the approximation of the gradient by the fourth order hermitian derivative,with a conservative discrete formulation on boxes of length 2h.The iterative technique is based on the repeated solution by a fast direct method of a discrete Poisson equation on a uniform rectangular mesh.The problem is suitably scaled before iteration.The numerical results obtained show the efficiency of the numerical scheme.This work is the extension to strongly elliptic problems of the hermitian box-scheme presented by Abbas and Croisille(J.Sci.Comput.,49(2011),pp.239–267).

POISSON PRECONDITIONING FOR SELF-ADJOINT ELLIPTIC PROBLEMS

In this paper, we formulate interface problem and Neumann elliptic boundary value problem into a form of linear operator equations with self-adjoint positive definite op- erators. We prove that in the discrete level the condition number of these operators is independent of the mesh size. Therefore, given a prescribed error tolerance, the classical conjugate gradient algorithm converges within a fixed number of iterations. The main computation task at each iteration is to solve a Dirichlet Poisson boundary value problem in a rectangular domain, which can be furnished with fast Poisson solver. The overall computational complexity is essentially of linear scaling.

An Augmented Approach for the Pressure Boundary Condition in a Stokes Flow

An augmented method is proposed for solving stationary incompressible Stokes equations with a Dirichlet boundary condition along parts of the boundary.In this approach,the normal derivative of the pressure along the parts of the boundary is introduced as an additional variable and it is solved by the GMRES iterative method.The dimension of the augmented variable in discretization is the number of grid points along the boundary which is O(N).Each GMRES iteration(or one matrix-vector multiplication)requires three fast Poisson solvers for the pressure and the velocity.In our numerical experiments,only a few iterations are needed.We have also combined the augmented approach for Stokes equations involving interfaces,discontinuities,and singularities.