Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA Spacetime Finite Elements Method to Solve the Dirac Equation
75
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
In this work, a fully implicit numerical approach based on space-time finite element method is presented to solve the Dirac equation in 1 (space) + 1 (time), 2 + 1, and 3 + 1 dimensions. We utilize PETSc/Tao library to implement our linear system and for using Krylov subspace based solvers such as GMRES. We demonstrate our method by analyzing several different cases including plane wave solution, Zitterbewegung, and Klein paradox. Parallel performance of this implementation is also presented.
rate research
Read More
In this paper, we provide a procedure to solve the eigen solutions of Dirac equation with complicated potential approximately. At first, we solve the eigen solutions of a linear Dirac equation with complete eigen system, which approximately equals to the original equation. Take the eigen functions as base of Hilbert space, and expand the spinor on the bases, we convert the original problem into solution of extremum of an algebraic function on the unit sphere of the coefficients. Then the problem can be easily solved. This is a standard finite element method with strict theory for convergence and effectiveness.
In this paper, we develop a simplified hybrid weighted essentially non-oscillatory (WENO) method combined with the modified ghost fluid method (MGFM) [28] to simulate the compressible two-medium flow problems. The MGFM can turn the two-medium flow problems into two single-medium cases by defining the ghost fluids status in terms of the predicted the interface status, which makes the material interface invisible. For the single medium flow case, we adapt between the linear upwind scheme and the WENO scheme automatically by identifying the regions of the extreme points for the reconstruction polynomial as same as the hybrid WENO scheme [50]. Instead of calculating their exact locations, we only need to know the regions of the extreme points based on the zero point existence theorem, which is simpler for implementation and saves computation time. Meanwhile, it still keeps the robustness and has high efficiency. Extensive numerical results for both one and two dimensional two-medium flow problems are performed to demonstrate the good performances of the proposed method.
This research was mainly conducted to explore the possibility of formulating an efficient algorithm to find roots of nonlinear equations without using the derivative of the function. The Weerakoon-Fernando method had been taken as the base in this project to find a new method without the derivative since Weerakoon-Fernando method gives 3rd order convergence. After several unsuccessful attempts we were able to formulate the Finite Difference Weerakoon-Fernando Method (FDWFM) presented here. We noticed that the FDWFM approaches the root faster than any other existing method in the absence of the derivatives as an example, the popular nonlinear equation solver such as secant method (order of convergence is 1.618) in the absence of the derivative. And the FDWFM had three function evaluations and secant method had two function evaluations. By implementing FDWFM on nonlinear equations with complex roots and also on systems of nonlinear equations, we received very encouraging results. When applying the FDWFM to systems of nonlinear equations, we resolved the involvement of the Jacobian problem by following the procedure in the Broydens method. The computational order of convergence of the FDWFM was close to 2.5 for all these cases. This will undoubtedly provide scientists the efficient numerical algorithm, that doesnt need the derivative of the function to solve nonlinear equations, that they were searching for over centuries.
In this paper, authors focus effort on improving the conventional discrete velocity method (DVM) into a multiscale scheme in finite volume framework for gas flow in all flow regimes. Unlike the typical multiscale kinetic methods unified gas-kinetic scheme (UGKS) and discrete unified gas-kinetic scheme (DUGKS), which concentrate on the evolution of the distribution function at the cell interface, in the present scheme the flux for macroscopic variables is split into the equilibrium part and the nonequilibrium part, and the nonequilibrium flux is calculated by integrating the discrete distribution function at the cell center, which overcomes the excess numerical dissipation of the conventional DVM in the continuum flow regime. Afterwards, the macroscopic variables are finally updated by simply integrating the discrete distribution function at the cell center, or by a blend of the increments based on the macroscopic and the microscopic systems, and the multiscale property is achieved. Several test cases, involving unsteady, steady, high speed, low speed gas flows in all flow regimes, have been performed, demonstrating the good performance of the multiscale DVM from free molecule to continuum Navier-Stokes solutions and the multiscale property of the scheme is proved.
We consider the null controllability problem for the wave equation, and analyse a stabilized finite element method formulated on a global, unstructured spacetime mesh. We prove error estimates for the approximate control given by the computational method. The proofs are based on the regularity properties of the control given by the Hilbert Uniqueness Method, together with the stability properties of the numerical scheme. Numerical experiments illustrate the results.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
