The Experts below are selected from a list of 5298 Experts worldwide ranked by ideXlab platform
Etienne Emmrich - One of the best experts on this subject based on the ideXlab platform.
-
convergence of the variable two step bdf time discretisation of nonlinear evolution problems governed by a monotone potential operator
Bit Numerical Mathematics, 2009Co-Authors: Etienne EmmrichAbstract:The initial-value problem for a first-order evolution equation is discretised in time by means of the two-step backward Differentiation Formula (BDF) on a variable time grid. The evolution equation is governed by a monotone and coercive potential operator. On a suitable sequence of time grids, the piecewise constant interpolation and a piecewise linear prolongation of the time discrete solution are shown to converge towards the weak solution if the ratios of adjacent step sizes are close to 1 and do not vary too much.
-
two step bdf time discretisation of nonlinear evolution problems governed by monotone operators with strongly continuous perturbations
Computational methods in applied mathematics, 2009Co-Authors: Etienne EmmrichAbstract:The time discretisation of the initial-value problem for a first-order evolution equation by the two-step backward Differentiation Formula (BDF) on a uniform grid is analysed. The evolution equation is governed by a time-dependent monotone operator that might be perturbed by a time-dependent strongly continuous operator. Well-posedness of the numerical scheme, a priori estimates, convergence of a piecewise polynomial prolongation, stability as well as smooth-data error estimates are provided relying essentially on an algebraic relation that implies the G-stability of the two-step BDF with constant time steps. 2000 Mathematics Subject Classification: 65M12, 65M15, 47J35, 35K55, 47H05.
-
error of the two step bdf for the incompressible navier stokes problem
Mathematical Modelling and Numerical Analysis, 2004Co-Authors: Etienne EmmrichAbstract:The incompressible Navier-Stokes problem is discretized in time by the two-step backward Differentiation Formula. Error estimates are proved under feasible assumptions on the regularity of the exact solution avoiding hardly fulfillable compatibility conditions. Whereas the time-weighted velocity error is of optimal second order, the time-weighted error in the pressure is of first order. Suboptimal estimates are shown for a linearisation. The results cover both the two- and three-dimensional case.
Cheng Wang - One of the best experts on this subject based on the ideXlab platform.
-
a second order accurate scalar auxiliary variable sav numerical method for the square phase field crystal equation
arXiv: Numerical Analysis, 2020Co-Authors: Min Wang, Qiumei Huang, Cheng WangAbstract:In this paper we propose and analyze a second order accurate (in time) numerical scheme for the square phase field crystal (SPFC) equation, a gradient flow modeling crystal dynamics at the atomic scale in space but on diffusive scales in time. Its primary difference with the standard phase field crystal model is an introduction of the 4-Laplacian term in the free energy potential, which in turn leads to a much higher degree of nonlinearity. To make the numerical scheme linear while preserving the nonlinear energy stability, we make use of the scalar auxiliary variable (SAV) approach, in which a second order backward Differentiation Formula (BDF) is applied in the temporal stencil. Meanwhile, a direct application of the SAV method faces certain difficulties, due to the involvement of the 4-Laplacian term, combined with a derivation of the lower bound of the nonlinear energy functional. In the proposed numerical method, an appropriate decomposition for the physical energy functional is Formulated, so that the nonlinear energy part has a well-established global lower bound, and the rest terms lead to constant-coefficient diffusion terms with positive eigenvalues. In turn, the numerical scheme could be very efficiently implemented by constant-coefficient Poisson-like type solvers (via FFT), and energy stability is established by introducing an auxiliary variable, and an optimal rate convergence analysis is provided for the proposed SAV method. A few numerical experiments are also presented, which confirm the efficiency and accuracy of the proposed scheme.
-
second order semi implicit projection methods for micromagnetics simulations
Journal of Computational Physics, 2020Co-Authors: Changjian Xie, Cheng Wang, Carlos J Garciacervera, Zhennan Zhou, Jingrun ChenAbstract:Abstract Micromagnetics simulations require accurate approximation of the magnetization dynamics described by the Landau-Lifshitz-Gilbert equation, which is nonlinear, nonlocal, and has a non-convex constraint, posing interesting challenges in developing numerical methods. In this paper, we propose two second-order semi-implicit projection methods based on the second-order backward Differentiation Formula and the second-order interpolation Formula using the information at previous two temporal steps. Unconditional unique solvability of both methods is proved, with their second-order accuracy verified through numerical examples in both 1D and 3D. The efficiency of both methods is compared to that of another two popular methods. In addition, we test the robustness of both methods for the first benchmark problem with a ferromagnetic thin film material from National Institute of Standards and Technology.
-
a second order energy stable backward Differentiation Formula method for the epitaxial thin film equation with slope selection
Numerical Methods for Partial Differential Equations, 2018Co-Authors: Wenqiang Feng, Cheng Wang, Steven M Wise, Zhengru ZhangAbstract:In this paper, we study a novel second-order energy stable Backward Differentiation Formula (BDF) finite difference scheme for the epitaxial thin film equation with slope selection (SS). One major challenge for the higher oder in time temporal discretization is how to ensure an unconditional energy stability and an efficient numerical implementation. We propose a general framework for designing the higher order in time numerical scheme with unconditional energy stability by using the BDF method with constant coefficient stabilized terms. Based on the unconditional energy stability property, we derive an $L^\infty_h (0,T; H_{h}^2)$ stability for the numerical solution and provide an optimal the convergence analysis. To deal with the 4-Laplacian solver in an $L^{2}$ gradient flow at each time step, we apply an efficient preconditioned steepest descent algorithm and preconditioned nonlinear conjugate gradient algorithm to solve the corresponding nonlinear system. Various numerical simulations are present to demonstrate the stability and efficiency of the proposed schemes and slovers.
Jingrun Chen - One of the best experts on this subject based on the ideXlab platform.
-
second order semi implicit projection methods for micromagnetics simulations
Journal of Computational Physics, 2020Co-Authors: Changjian Xie, Cheng Wang, Carlos J Garciacervera, Zhennan Zhou, Jingrun ChenAbstract:Abstract Micromagnetics simulations require accurate approximation of the magnetization dynamics described by the Landau-Lifshitz-Gilbert equation, which is nonlinear, nonlocal, and has a non-convex constraint, posing interesting challenges in developing numerical methods. In this paper, we propose two second-order semi-implicit projection methods based on the second-order backward Differentiation Formula and the second-order interpolation Formula using the information at previous two temporal steps. Unconditional unique solvability of both methods is proved, with their second-order accuracy verified through numerical examples in both 1D and 3D. The efficiency of both methods is compared to that of another two popular methods. In addition, we test the robustness of both methods for the first benchmark problem with a ferromagnetic thin film material from National Institute of Standards and Technology.
Carlos J Garciacervera - One of the best experts on this subject based on the ideXlab platform.
-
second order semi implicit projection methods for micromagnetics simulations
Journal of Computational Physics, 2020Co-Authors: Changjian Xie, Cheng Wang, Carlos J Garciacervera, Zhennan Zhou, Jingrun ChenAbstract:Abstract Micromagnetics simulations require accurate approximation of the magnetization dynamics described by the Landau-Lifshitz-Gilbert equation, which is nonlinear, nonlocal, and has a non-convex constraint, posing interesting challenges in developing numerical methods. In this paper, we propose two second-order semi-implicit projection methods based on the second-order backward Differentiation Formula and the second-order interpolation Formula using the information at previous two temporal steps. Unconditional unique solvability of both methods is proved, with their second-order accuracy verified through numerical examples in both 1D and 3D. The efficiency of both methods is compared to that of another two popular methods. In addition, we test the robustness of both methods for the first benchmark problem with a ferromagnetic thin film material from National Institute of Standards and Technology.
Yves Bourgault - One of the best experts on this subject based on the ideXlab platform.
-
efficient second order semi implicit finite element method for fourth order nonlinear diffusion equations
Computer Physics Communications, 2021Co-Authors: Sana Keita, Abdelaziz Beljadid, Yves BourgaultAbstract:Abstract We focus here on a class of fourth-order parabolic equations that can be written as a system of second-order equations by introducing an auxiliary variable. We design a novel second-order fully discrete mixed finite element method to approximate these equations. In our approach, we propose new techniques using the second-order backward Differentiation Formula for the time derivative and a special technique for the approximation of nonlinear terms. The use of the proposed technique for nonlinear terms makes the developed numerical scheme efficient in terms of computational cost since the proposed method only deals with a linear system at each time step and no iterative resolution is needed. A numerical convergence study is performed using the method of manufactured and analytical solutions of the system where we investigate different boundary conditions. With respect to the spatial discretization, convergence rates are found to at least match a priori error estimates available for linear problems. The convergence analysis is completed with an investigation of the temporal discretization where we numerically demonstrate the second-order time-accuracy of the proposed scheme using the method of reference solution. We present a series of numerical tests to demonstrate the efficiency and robustness of the proposed scheme.