اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدAdvantages of a semi-implicit scheme over a fully implicit scheme for Landau-Lifshitz-Gilbert equation
202
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
Magnetization dynamics in magnetic materials is modeled by the Landau-Lifshitz-Gilbert (LLG) equation. In the LLG equation, the length of magnetization is conserved and the system energy is dissipative. Implicit and semi-implicit schemes have been used in micromagnetics simulations due to their unconditional numerical stability. In more details, implicit schemes preserve the properties of the LLG equation, but solve a nonlinear system of equations per time step. In contrast, semi-implicit schemes only solve a linear system of equations, while additional operations are needed to preserve the length of magnetization. It still remains unclear which one shall be used if both implicit and semi-implicit schemes are available. In this work, using the implicit Crank-Nicolson (ICN) scheme as a benchmark, we propose to make this implicit scheme semi-implicit. It can be proved that both schemes are second-order accurate in space and time. For the unique solvability of nonlinear systems of equations in the ICN scheme, we require that the temporal step size scales quadratically with the spatial mesh size. It is numerically verified that the convergence of the nonlinear solver becomes slower for larger temporal step size and multiple magnetization profiles are obtained for different initial guesses. The linear systems of equations in the semi-implicit CN (SICN) scheme are unconditionally uniquely solvable, and the condition that the temporal step size scales linearly with the spatial mesh size is needed in the convergence of the SICN scheme. In terms of numerical efficiency, the SICN scheme achieves the same accuracy as the ICN scheme with less computational time. Based on these results, we conclude that a semi-implicit scheme is superior to its implicit analog both theoretically and numerically, and we recommend the semi-implicit scheme in micromagnetics simulations if both methods are available.
قيم البحث
اقرأ أيضاً
Recent theoretical and experimental advances show that the inertia of magnetization emerges at sub-picoseconds and contributes to the ultrafast magnetization dynamics which cannot be captured intrinsically by the LLG equation. Therefore, as a general
ization, the inertial Landau-Lifshitz-Gilbert (iLLG) equation is proposed to model the ultrafast magnetization dynamics. Mathematically, the LLG equation is a nonlinear system of parabolic type with (possible) degeneracy. However, the iLLG equation is a nonlinear system of mixed hyperbolic-parabolic type with degeneracy, and exhibits more complicated structures. It behaves like a hyperbolic system at the sub-picosecond scale while behaves like a parabolic system at larger timescales. Such hybrid behaviors impose additional difficulties on designing numerical methods for the iLLG equation. In this work, we propose a second-order semi-implicit scheme to solve the iLLG equation. The second temporal derivative of magnetization is approximated by the standard centered difference scheme and the first derivative is approximated by the midpoint scheme involving three time steps. The nonlinear terms are treated semi-implicitly using one-sided interpolation with the second-order accuracy. At each step, the unconditionally unique solvability of the unsymmetric linear system of equations in the proposed method is proved with a detailed discussion on the condition number. Numerically, the second-order accuracy in both time and space is verified. Using the proposed method, the inertial effect of ferromagnetics is observed in micromagnetics simulations at small timescales, in consistency with the hyperbolic property of the model at sub-picoseconds. For long time simulations, the results of the iLLG model are in nice agreements with those of the LLG model, in consistency with the parabolic feature of the iLLG model at larger timescales.
We propose new semi-implicit numerical methods for the integration of the stochastic Landau-Lifshitz equation with built-in angular momentum conservation. The performance of the proposed integrators is tested on the 1D Heisenberg chain. For this syst
em, our schemes show better stability properties and allow us to use considerably larger time steps than standard explicit methods. At the same time, these semi-implicit schemes are also of comparable accuracy to and computationally much cheaper than the standard midpoint implicit method. The results are of key importance for atomistic spin dynamics simulations and the study of spin dynamics beyond the macro spin approximation.
The Reynolds-Averaged Navier-Stokes equations and the Large-Eddy Simulation equations can be coupled using a transition function to switch from a set of equations applied in some areas of a domain to the other set in the other part of the domain. Fol
lowing this idea, different time integration schemes can be coupled. In this context, we developed a hybrid time integration scheme that spatially couples the explicit scheme of Heun and the implicit scheme of Crank and Nicolson using a dedicated transition function. This scheme is linearly stable and second-order accurate. In this paper, an extension of this hybrid scheme is introduced to deal with a temporal adaptive procedure. The idea is to treat the time integration procedure with unstructured grids as it is performed with Cartesian grids with local mesh refinement. Depending on its characteristic size, each mesh cell is assigned a rank. And for two cells from two consecutive ranks, the ratio of the associated time steps for time marching the solutions is $2$. As a consequence, the cells with the lowest rank iterate more than the other ones to reach the same physical time. In a finite-volume context, a key ingredient is to keep the conservation property for the interfaces that separate two cells of different ranks. After introducing the different schemes, the paper recalls briefly the coupling procedure, and details the extension to the temporal adaptive procedure. The new time integration scheme is validated with the propagation of 1D wave packet, the Sods tube, and the transport of a bi-dimensional vortex in an uniform flow.
An implicit multiscale method with multiple macroscopic prediction for steady state solutions of gas flow in all flow regimes is presented. The method is based on the finite volume discrete velocity method (DVM) framework. At the cell interface a mul
tiscale flux with a construction similar to discrete unified gas-kinetic scheme (DUGKS) is adopted. The idea of the macroscopic variable prediction is further developed and a multiple prediction structure is formed. A prediction scheme is constructed to give a predicted macroscopic variable based on the macroscopic residual, and the convergence is accelerated greatly in the continuum flow regime. Test cases show the present method is one order of magnitude faster than the previous implicit multiscale scheme in the continuum flow regime.
An efficient implicit kinetic scheme is developed to solve the stationary phonon Boltzmann transport equation (BTE) based on the non-gray model including the phonon dispersion and polarization. Due to the wide range of the dispersed phonon mean free
paths, the phonon transport under the non-gray model is essentially multiscale, and has to be solved differently and appropriately for varied phonon frequencies and branches. The proposed implicit kinetic scheme is composed of a microscopic iteration and a macroscopic iteration. The microscopic iteration is capable of automatically adapting with varied phonon mean free path of each phonon frequency and branch through solving the phonon BTE. The energy transfer of all phonons is gathered together by the microscopic iteration to evaluate the heat flux. The temperature field is predicted through a macroscopic heat transfer equation according to the heat flux, and the equilibrium state in the phonon BTE is also updated. The combination of the phonon BTE solver and the macroscopic equation makes the present method very efficient in a wide length scale. Three numerical tests, including the cross-plane, in-plane and nano-porous heat transfer in silicon, validate that the present scheme can handle with the phonon dispersion and polarization correctly and predict the multiscale heat transfer phenomena efficiently in a wide range. The present method could be tens of times faster than the typical implicit DOM and keeps the same amount of the memory requirements as the Fourier solver for multiscale heat transfer problem.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
