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. F...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.展开更多
在导弹类金属-介质复合目标电磁散射特性求解过程中,采用常规迭代求解方法存在难以收敛以及内迭代边界积分区域重复求解的问题。针对该问题,在传统有限元边界积分区域分解法(finite element boundary integral domain decomposition met...在导弹类金属-介质复合目标电磁散射特性求解过程中,采用常规迭代求解方法存在难以收敛以及内迭代边界积分区域重复求解的问题。针对该问题,在传统有限元边界积分区域分解法(finite element boundary integral domain decomposition method,FE-BI-DDM)的基础上,采用了更为灵活的多区多求解器的方法(multi domain multi solver method,MDMSM)。该方法对导弹类金属-介质复合目标中难以收敛的金属区域,使用快速直接求逆的方法求解,由于可以使用独立的网格模型进行电磁建模,避免了内迭代部分的模型重复建立过程,从而大幅减少了整体模型求解时间。实验结果表明:所提方法可以在相同计算精度的条件下,以不过多增加内存空间为前提,大幅缩短了导弹类目标的金属-介质复合模型的电磁求解时间。该方法为开展导弹类目标特性分析提供了一条可行的技术途径。展开更多
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 fi...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.展开更多
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 coeff...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.展开更多
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...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.展开更多
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...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.展开更多
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 dia...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.展开更多
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 lo...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.展开更多
基金Open Funds of State Key Laboratory of MillimeterWaves,China (No. K200401), Outstanding Teaching and ResearchAwards for Young Teachers of Nanjing Normal University (No.1320BL51)
文摘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.
文摘在导弹类金属-介质复合目标电磁散射特性求解过程中,采用常规迭代求解方法存在难以收敛以及内迭代边界积分区域重复求解的问题。针对该问题,在传统有限元边界积分区域分解法(finite element boundary integral domain decomposition method,FE-BI-DDM)的基础上,采用了更为灵活的多区多求解器的方法(multi domain multi solver method,MDMSM)。该方法对导弹类金属-介质复合目标中难以收敛的金属区域,使用快速直接求逆的方法求解,由于可以使用独立的网格模型进行电磁建模,避免了内迭代部分的模型重复建立过程,从而大幅减少了整体模型求解时间。实验结果表明:所提方法可以在相同计算精度的条件下,以不过多增加内存空间为前提,大幅缩短了导弹类目标的金属-介质复合模型的电磁求解时间。该方法为开展导弹类目标特性分析提供了一条可行的技术途径。
文摘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.
基金supported by NFS No.11001257,was stimulated by Per-Gunnar Martinsson’s paper”A Fast Direct Solver for a Class of Elliptic Partial Differential Equations”.Professor Jingfang Huang suggested solving the Poisson equation with variable coefficient as a test case.We are very grateful to both of them for their selfless help.
文摘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.
基金The author is grateful for Professor Tao Tang and Dr.Zhonghua Qiao for many helpful and fruitful discussions,and would like to thank Professor Weiwei Sun for constructive suggestions。
文摘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.
文摘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.
基金supported by National Natural Science Foundation of China (Grant Nos.12171300 and 11771338)the Natural Science Foundation of Shanghai (Grant No.21ZR1480500)the Fundamental Research Funds for the Central Universities (Grant No.2019110066)。
文摘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.
基金the supports by National Science and Technology Council,Taiwan,under the research grants 111-2115-M-008-009-MY3,111-2628-M-A49-008-MY4,111-2115-M-390-002,and 110-2115-M-A49-011-MY3,respectivelythe supports by National Center for Theoretical Sciences,Taiwan.
文摘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.
