Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA Finite Difference Method on Irregular Grids with Local Second Order Ghost Point Extension for Solving Maxwells Equations Around Curved PEC Objects
71
0
0.0
(
0
)
Added by
Haiyu Zou
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
A new finite difference method on irregular, locally perturbed rectangular grids has been developed for solving electromagnetic waves around curved perfect electric conductors (PEC). This method incorporates the back and forth error compensation and correction method (BFECC) and level set method to achieve convenience and higher order of accuracy at complicated PEC boundaries. A PDE-based local second order ghost cell extension technique is developed based on the level set framework in order to compute the boundary value to first order accuracy (cumulatively), and then BFECC is applied to further improve the accuracy while increasing the CFL number. Numerical experiments are conducted to validate the properties of the method.
rate research
Read More
At present, deep learning based methods are being employed to resolve the computational challenges of high-dimensional partial differential equations (PDEs). But the computation of the high order derivatives of neural networks is costly, and high order derivatives lack robustness for training purposes. We propose a novel approach to solving PDEs with high order derivatives by simultaneously approximating the function value and derivatives. We introduce intermediate variables to rewrite the PDEs into a system of low order differential equations as what is done in the local discontinuous Galerkin method. The intermediate variables and the solutions to the PDEs are simultaneously approximated by a multi-output deep neural network. By taking the residual of the system as a loss function, we can optimize the network parameters to approximate the solution. The whole process relies on low order derivatives. Numerous numerical examples are carried out to demonstrate that our local deep learning is efficient, robust, flexible, and is particularly well-suited for high-dimensional PDEs with high order derivatives.
In this paper, we combine the nonlinear HWENO reconstruction in cite{newhwenozq} and the fixed-point iteration with Gauss-Seidel fast sweeping strategy, to solve the static Hamilton-Jacobi equations in a novel HWENO framework recently developed in cite{mehweno1}. The proposed HWENO frameworks enjoys several advantages. First, compared with the traditional HWENO framework, the proposed methods do not need to introduce additional auxiliary equations to update the derivatives of the unknown function $phi$. They are now computed from the current value of $phi$ and the previous spatial derivatives of $phi$. This approach saves the computational storage and CPU time, which greatly improves the computational efficiency of the traditional HWENO scheme. In addition, compared with the traditional WENO method, reconstruction stencil of the HWENO methods becomes more compact, their boundary treatment is simpler, and the numerical errors are smaller on the same mesh. Second, the fixed-point fast sweeping method is used to update the numerical approximation. It is an explicit method and does not involve the inverse operation of nonlinear Hamiltonian, therefore any Hamilton-Jacobi equations with complex Hamiltonian can be solved easily. It also resolves some known issues, including that the iterative number is very sensitive to the parameter $varepsilon$ used in the nonlinear weights, as observed in previous studies. Finally, in order to further reduce the computational cost, a hybrid strategy is also presented. Extensive numerical experiments are performed on two-dimensional problems, which demonstrate the good performance of the proposed fixed-point fast sweeping HWENO methods.
We propose a controllability method for the numerical solution of time-harmonic Maxwells equations in their first-order formulation. By minimizing a quadratic cost functional, which measures the deviation from periodicity, the controllability method determines iteratively a periodic solution in the time domain. At each conjugate gradient iteration, the gradient of the cost functional is simply computed by running any time-dependent simulation code forward and backward for one period, thus leading to a non-intrusive implementation easily integrated into existing software. Moreover, the proposed algorithm automatically inherits the parallelism, scalability, and low memory footprint of the underlying time-domain solver. Since the time-periodic solution obtained by minimization is not necessarily unique, we apply a cheap post-processing filtering procedure which recovers the time-harmonic solution from any minimizer. Finally, we present a series of numerical examples which show that our algorithm greatly speeds up the convergence towards the desired time-harmonic solution when compared to simply running the time-marching code until the time-harmonic regime is eventually reached.
We develop a stable finite difference method for the elastic wave equations in bounded media, where the material properties can be discontinuous at curved interfaces. The governing equations are discretized in second order form by a fourth or sixth order accurate summation-by-parts operator. The mesh size is determined by the velocity structure of the material, resulting in nonconforming grid interfaces with hanging nodes. We use order-preserving interpolation and the ghost point technique to couple adjacent mesh blocks in an energy-conserving manner, which is supported by a fully discrete stability analysis. In numerical experiments, we demonstrate that the convergence rate is optimal, and is the same as when a globally uniform mesh is used in a single domain. In addition, with a predictor-corrector time integration method, we obtain time stepping stability with stepsize almost the same as given by the usual Courant Friedrichs Lewy condition.
The locally modified finite element method, which is introduced in [Frei, Richter: SINUM 52(2014), p. 2315-2334] is a simple fitted finite element method that is able to resolve weak discontinuities in interface problems. The method is based on a fixed structured coarse mesh, which is then refined into sub-elements to resolve an interior interface. In this work, we extend the locally modified finite element method to second order using an isoparametric approach in the interface elements. Thereby we need to take care that the resulting curved edges do not lead to degenerate sub-elements. We prove optimal a priori error estimates in the $L^2$-norm and in a modified energy norm, as well as a reduced convergence order of ${cal O}(h^{3/2})$ in the standard $H^1$-norm. Finally, we present numerical examples to substantiate the theoretical findings.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
