In modelling elastic wave propagation in a porous medium, when the ratio between the fluid viscosity and the medium permeability is comparatively large, the stiffness problem of Blot＇s poroelastic equations will be e...In modelling elastic wave propagation in a porous medium, when the ratio between the fluid viscosity and the medium permeability is comparatively large, the stiffness problem of Blot＇s poroelastic equations will be encountered. In the paper, a partition method is developed to solve the stiffness problem with a staggered high-order finite-difference. The method splits the Biot equations into two systems. One is stiff, and solved analytically, the other is nonstiff, and solved numerically by using a high-order staggered-grid finite-difference scheme. The time step is determined by the staggered finite-difference algorithm in solving the nonstiff equations, thus a coarse time step may be employed. Therefore, the computation efficiency and computational stability are improved greatly. Also a perfect by matched layer technology is used in the split method as absorbing boundary conditions. The numerical results are compared with the analytical results and those obtained from the conventional staggered-grid finite-difference method in a homogeneous model, respectively. They are in good agreement with each other. Finally, a slightly more complex model is investigated and compared with related equivalent model to illustrate the good performance of the staggered-grid finite-difference scheme in the partition method.展开更多
Using a discretized finite difference method, a numerical model was developed to study the interaction of regular waves with a perforated breakwater. Considering a non-viscous, non-rotational fluid, the governing equa...Using a discretized finite difference method, a numerical model was developed to study the interaction of regular waves with a perforated breakwater. Considering a non-viscous, non-rotational fluid, the governing equations of Laplacian velocity potential were developed, and specific conditions for every single boundary were defined. The final developed model was evaluated based on an existing experimental result. The evaluated model was used to simulate the condition for various wave periods from 0.6 to 2 s. The reflection coefficient and transmission coefficient of waves were examined with different breakwater porosities, wave steepnesses, and angular frequencies. The results show that the developed model can suitably present the effect of the structural and hydraulic parameters on the reflection and transmission coefficients. It was also found that with the increase in wave steepness, the reflection coefficient increased logarithmically, while the transmission coefficient decreased logarithmically.展开更多
A modified alternating direction implicit algorithm is proposed to solve the full-vectorial finite-difference beam propagation method formulation based on H fields. The cross-coupling terms are neglected in the first ...A modified alternating direction implicit algorithm is proposed to solve the full-vectorial finite-difference beam propagation method formulation based on H fields. The cross-coupling terms are neglected in the first sub-step, but evaluated and doubly used in the second sub-step. The order of two sub-steps is reversed for each transverse magnetic field component so that the cross-coupling terms are always expressed in implicit form, thus the calculation is very efficient and stable. Moreover, an improved six-point finite-difference scheme with high accuracy independent of specific structures of waveguide is also constructed to approximate the cross-coupling terms along the transverse directions. The imaginary-distance procedure is used to assess the validity and utility of the present method. The field patterns and the normalized propagation constants of the fundamental mode for a buried rectangular waveguide and a rib waveguide are presented. Solutions are in excellent agreement with the benchmark results from the modal transverse resonance method.展开更多
In this study,single-particle energy was examined using the finite difference method by taking 208Pb as an example.If the first derivative term in the spherical Dirac equation is discretized using a three-point formul...In this study,single-particle energy was examined using the finite difference method by taking 208Pb as an example.If the first derivative term in the spherical Dirac equation is discretized using a three-point formula,a one-to-one correspondence occurs between the physical and spurious states.Although these energies are exactly the same,the wave functions of the spurious states exhibit a much faster staggering than those of the physical states.Such spurious states can be eliminated when applying the finite difference method by introducing an extra Wilson term into the Hamiltonian.Furthermore,it was also found that the number of spurious states can be reduced if we improve the accuracy of the numerical differential formula.The Dirac equation is then solved in a momentum space in which there is no differential operator,and we found that the spurious states can be completely avoided in the momentum space,even without an extra Wilson term.展开更多
In this paper, three implicit finite difference methods are developed to solve one dimensional time fractional advection-diffusion equation. The fractional derivative is treated by applying right shifted Grünwald...In this paper, three implicit finite difference methods are developed to solve one dimensional time fractional advection-diffusion equation. The fractional derivative is treated by applying right shifted Grünwald-Letnikov formula of order α ∈(0, 1). We investigate the stability analysis by using von Neumann method with mathematical induction and prove that these three proposed methods are unconditionally stable. Numerical results are presented to demonstrate the effectiveness of the schemes mentioned in this paper.展开更多
This paper investigates the phenomenon of three-pulse photon echo in thick rare-earth ions doped crystal whose thickness is far larger than 0.002 cm which is adopted in previous works.The influence of thickness on the...This paper investigates the phenomenon of three-pulse photon echo in thick rare-earth ions doped crystal whose thickness is far larger than 0.002 cm which is adopted in previous works.The influence of thickness on the three-pulse photon echo＇s amplitude and efficiency is analyzed with the Maxwell-Bloch equations solved by finite-difference timedomain method.We demonstrate that the amplitude of three-pulse echo will increase with the increasing of thickness and the optimum thickness to generate three-pulse photon echo is 0.3 cm for Tm^（3＋）：YAG when the attenuation of the input pulse is taken into account.Meanwhile,we find the expression 0.09 exp（α＇L）,which is previously employed to describe the relationship between echo＇s efficiency and thickness,should be modified as 1.3 · 0.09 exp（2.4 ·α＇L） with the propagation of echo considered.展开更多
Coupled system of multilayer dynamics of fluids in porous media is to describe the history of oil-gas transport and accumulation in basin evolution.It is of great value in rational evaluation of prospecting and exploi...Coupled system of multilayer dynamics of fluids in porous media is to describe the history of oil-gas transport and accumulation in basin evolution.It is of great value in rational evaluation of prospecting and exploiting oil-gas resources.The mathematical model can be described as a coupled system of nonlinear partial differential equations with moving boundary values.The upwind finite difference schemes applicable to parallel arithmetic are put forward and two-dimensional and three-dimensional schemes are used to form a complete set.Some techniques,such as change of variables,calculus of variations, multiplicative commutation rule of difference operators,decomposition of high order difference operators and prior estimates,are adopted.The estimates in l~2 norm are derived to determine the error in the approximate solution.This method was already applied to the numerical simulation of migration-accumulation of oil resources.展开更多
The finite-difference time-domain method is used to simulate the optical characteristics of an in-plane switching blue phase liquid crystal display.Compared with the matrix optic methods and the refractive method,the ...The finite-difference time-domain method is used to simulate the optical characteristics of an in-plane switching blue phase liquid crystal display.Compared with the matrix optic methods and the refractive method,the finite-difference timedomain method,which is used to directly solve Maxwell＇s equations,can consider the lateral variation of the refractive index and obtain an accurate convergence effect.The simulation results show that e-rays and o-rays bend in different directions when the in-plane switching blue phase liquid crystal display is driven by the operating voltage.The finitedifference time-domain method should be used when the distribution of the liquid crystal in the liquid crystal display has a large lateral change.展开更多
An efficient conformal locally one-dimensional finite-difference time-domain（LOD-CFDTD） method is presented for solving two-dimensional（2D） electromagnetic（EM） scattering problems. The formulation for the 2D tra...An efficient conformal locally one-dimensional finite-difference time-domain（LOD-CFDTD） method is presented for solving two-dimensional（2D） electromagnetic（EM） scattering problems. The formulation for the 2D transverse-electric（TE） case is presented and its stability property and numerical dispersion relationship are theoretically investigated. It is shown that the introduction of irregular grids will not damage the numerical stability. Instead of the staircasing approximation, the conformal scheme is only employed to model the curve boundaries, whereas the standard Yee grids are used for the remaining regions. As the irregular grids account for a very small percentage of the total space grids, the conformal scheme has little effect on the numerical dispersion. Moreover, the proposed method, which requires fewer arithmetic operations than the alternating-direction-implicit（ADI） CFDTD method, leads to a further reduction of the CPU time. With the total-field/scattered-field（TF/SF） boundary and the perfectly matched layer（PML）, the radar cross section（RCS） of two2 D structures is calculated. The numerical examples verify the accuracy and efficiency of the proposed method.展开更多
An explicit,time-dependent variable grid finite difference method is introduced and analyzed for approximating the solution of a scalar conservation law in two dimension. The scheme is stable,and the numerical solutio...An explicit,time-dependent variable grid finite difference method is introduced and analyzed for approximating the solution of a scalar conservation law in two dimension. The scheme is stable,and the numerical solution is proved to converges to the relevant physical solution.展开更多
The key issue in accelerating method of characteristics(MOC)transport calculations is in obtaining a completely equivalent low-order neutron transport or diffusion equation.Herein,an equivalent low-order angular flux ...The key issue in accelerating method of characteristics(MOC)transport calculations is in obtaining a completely equivalent low-order neutron transport or diffusion equation.Herein,an equivalent low-order angular flux nonlinear finite difference equation is proposed for MOC transport calculations.This method comprises three essential features:(1)the even parity discrete ordinates method is used to build a low-order angular flux nonlinear finite difference equation,and different boundary condition treatments are proposed;(2)two new defined factors,i.e.,the even parity discontinuity factor and odd parity discontinuity factor,are strictly defined to achieve equivalence between the low-order angular flux nonlinear finite difference method and MOC transport calculation;(3)the energy group and angle are decoupled to construct a symmetric linear system that is much easier to solve.The equivalence of this low-order angular flux nonlinear finite difference equation is analyzed for two-dimensional(2D)pin,2D assembly,and 2D C5G7 benchmark problems.Numerical results demonstrate that a low-order angular flux nonlinear finite difference equation that is completely equivalent to the pin-resolved transport equation is established.展开更多
In the article, the fully discrete finite difference scheme for a type of nonlinear reaction-diffusion equation is established. Then the new function space is introduced and the stability problem for the finite differ...In the article, the fully discrete finite difference scheme for a type of nonlinear reaction-diffusion equation is established. Then the new function space is introduced and the stability problem for the finite difference scheme is discussed by means of variational approximation method in this function space. The approach used is of a simple characteristic in gaining the stability condition of the scheme.展开更多
The transmission coefficients of electromagnetic （EM） waves due to a superconductor-dielectric superlattice are numerically calculated. Shift operator finite difference time domain （SO-FDTD） method is used in the ...The transmission coefficients of electromagnetic （EM） waves due to a superconductor-dielectric superlattice are numerically calculated. Shift operator finite difference time domain （SO-FDTD） method is used in the analysis. By using the SO-FDTD method, the transmission spectrum is obtained and its characteristics are investigated for different thicknesses of superconductor layers and dielectric layers, from which a stop band starting from zero frequency can be apparently observed. The relation between this low-frequency stop band and relative temperature, and also the London penetration depth at a superconductor temperature of zero degree are discussed, separately. The low-frequency stop band properties of superconductor-dielectric superlattice thus are well disclosed.展开更多
The wave force exerted on vertical piles of offshore structures is the main criterion in designing them.In structures with more than one large pile,the influence of piles on each other is one of the most important iss...The wave force exerted on vertical piles of offshore structures is the main criterion in designing them.In structures with more than one large pile,the influence of piles on each other is one of the most important issues being concerned in past researches.An efficient method for determining the interaction of piles is introduced in present research.First the wave force is calculated by the exact method using the diffraction theory,then in the finite difference numerical method the force is calculated by adding the velocity potentials of each pile and integration of pressure on their surface.The results showed that the ratio of the wave force on each of the double piles to a single pile has a damped oscillation around unity in which the amplitude of oscillation decreases with the increase in the spacing parameter.Also different wave incident directions and diffraction parameters were used and the results showed that the numerical solution has acceptable accuracy when the diffraction parameter is larger than unity.展开更多
In this work, we first derive the one-point large deviations principle (LDP) for both the stochastic Cahn–Hilliard equation with small noise and its spatial finite difference method (FDM). Then, we focus on giving th...In this work, we first derive the one-point large deviations principle (LDP) for both the stochastic Cahn–Hilliard equation with small noise and its spatial finite difference method (FDM). Then, we focus on giving the convergence of the one-point large deviations rate function (LDRF) of the spatial FDM, which is about the asymptotical limit of a parametric variational problem. The main idea for proving the convergence of the LDRF of the spatial FDM is via the Γ-convergence of objective functions. This relies on the qualitative analysis of skeleton equations of the original equation and the numerical method. In order to overcome the difficulty that the drift coefficient is not one-sided Lipschitz continuous, we derive the equivalent characterization of the skeleton equation of the spatial FDM and the discrete interpolation inequality to obtain the uniform boundedness of the solution to the underlying skeleton equation. These play important roles in deriving the T-convergence of objective functions.展开更多
To enhance the computational efficiency of spatio-temporally discretized phase-field models,we present a high-speed solver specifically designed for the Poisson equations,a component frequently used in the numerical c...To enhance the computational efficiency of spatio-temporally discretized phase-field models,we present a high-speed solver specifically designed for the Poisson equations,a component frequently used in the numerical computation of such models.This efficient solver employs algorithms based on discrete cosine transformations(DCT)or discrete sine transformations(DST)and is not restricted by any spatio-temporal schemes.Our proposed methodology is appropriate for a variety of phase-field models and is especially efficient when combined with flow field systems.Meanwhile,this study has conducted an extensive numerical comparison and found that employing DCT and DST techniques not only yields results comparable to those obtained via the Multigrid(MG)method,a conventional approach used in the resolution of the Poisson equations,but also enhances computational efficiency by over 90%.展开更多
基金Project supported by the ＂100 Talents Project＂ of the Chinese Academy of Sciences and the Major Program of the National Natural Science Foundation of China （Grant No 10534040）.
文摘In modelling elastic wave propagation in a porous medium, when the ratio between the fluid viscosity and the medium permeability is comparatively large, the stiffness problem of Blot＇s poroelastic equations will be encountered. In the paper, a partition method is developed to solve the stiffness problem with a staggered high-order finite-difference. The method splits the Biot equations into two systems. One is stiff, and solved analytically, the other is nonstiff, and solved numerically by using a high-order staggered-grid finite-difference scheme. The time step is determined by the staggered finite-difference algorithm in solving the nonstiff equations, thus a coarse time step may be employed. Therefore, the computation efficiency and computational stability are improved greatly. Also a perfect by matched layer technology is used in the split method as absorbing boundary conditions. The numerical results are compared with the analytical results and those obtained from the conventional staggered-grid finite-difference method in a homogeneous model, respectively. They are in good agreement with each other. Finally, a slightly more complex model is investigated and compared with related equivalent model to illustrate the good performance of the staggered-grid finite-difference scheme in the partition method.
文摘Using a discretized finite difference method, a numerical model was developed to study the interaction of regular waves with a perforated breakwater. Considering a non-viscous, non-rotational fluid, the governing equations of Laplacian velocity potential were developed, and specific conditions for every single boundary were defined. The final developed model was evaluated based on an existing experimental result. The evaluated model was used to simulate the condition for various wave periods from 0.6 to 2 s. The reflection coefficient and transmission coefficient of waves were examined with different breakwater porosities, wave steepnesses, and angular frequencies. The results show that the developed model can suitably present the effect of the structural and hydraulic parameters on the reflection and transmission coefficients. It was also found that with the increase in wave steepness, the reflection coefficient increased logarithmically, while the transmission coefficient decreased logarithmically.
文摘A modified alternating direction implicit algorithm is proposed to solve the full-vectorial finite-difference beam propagation method formulation based on H fields. The cross-coupling terms are neglected in the first sub-step, but evaluated and doubly used in the second sub-step. The order of two sub-steps is reversed for each transverse magnetic field component so that the cross-coupling terms are always expressed in implicit form, thus the calculation is very efficient and stable. Moreover, an improved six-point finite-difference scheme with high accuracy independent of specific structures of waveguide is also constructed to approximate the cross-coupling terms along the transverse directions. The imaginary-distance procedure is used to assess the validity and utility of the present method. The field patterns and the normalized propagation constants of the fundamental mode for a buried rectangular waveguide and a rib waveguide are presented. Solutions are in excellent agreement with the benchmark results from the modal transverse resonance method.
基金partly supported by the National Natural Science Foundation of China(No.11875070)the Natural Science Foundation of Anhui Province(No.1908085MA16)
文摘In this study,single-particle energy was examined using the finite difference method by taking 208Pb as an example.If the first derivative term in the spherical Dirac equation is discretized using a three-point formula,a one-to-one correspondence occurs between the physical and spurious states.Although these energies are exactly the same,the wave functions of the spurious states exhibit a much faster staggering than those of the physical states.Such spurious states can be eliminated when applying the finite difference method by introducing an extra Wilson term into the Hamiltonian.Furthermore,it was also found that the number of spurious states can be reduced if we improve the accuracy of the numerical differential formula.The Dirac equation is then solved in a momentum space in which there is no differential operator,and we found that the spurious states can be completely avoided in the momentum space,even without an extra Wilson term.
文摘In this paper, three implicit finite difference methods are developed to solve one dimensional time fractional advection-diffusion equation. The fractional derivative is treated by applying right shifted Grünwald-Letnikov formula of order α ∈(0, 1). We investigate the stability analysis by using von Neumann method with mathematical induction and prove that these three proposed methods are unconditionally stable. Numerical results are presented to demonstrate the effectiveness of the schemes mentioned in this paper.
基金Project supported by Tianjin Research Program Application Foundation and Advanced Technology,China(Grant No.15JCQNJC01100)
文摘This paper investigates the phenomenon of three-pulse photon echo in thick rare-earth ions doped crystal whose thickness is far larger than 0.002 cm which is adopted in previous works.The influence of thickness on the three-pulse photon echo＇s amplitude and efficiency is analyzed with the Maxwell-Bloch equations solved by finite-difference timedomain method.We demonstrate that the amplitude of three-pulse echo will increase with the increasing of thickness and the optimum thickness to generate three-pulse photon echo is 0.3 cm for Tm^（3＋）：YAG when the attenuation of the input pulse is taken into account.Meanwhile,we find the expression 0.09 exp（α＇L）,which is previously employed to describe the relationship between echo＇s efficiency and thickness,should be modified as 1.3 · 0.09 exp（2.4 ·α＇L） with the propagation of echo considered.
基金supported by the Major State BasicResearch Program of China(19990328)the National Tackling Key Problem Programs(20050200069)+4 种基金the National Natural Science Foundation of China(1077112410372052)the Doctorate Foundation of the Ministryof Education of China(20030422047)Shandong Provance Natural Science Foundation(2R2009AQ12)the Independent Innovation Foundation of Shandong University(2010TS031)
文摘Coupled system of multilayer dynamics of fluids in porous media is to describe the history of oil-gas transport and accumulation in basin evolution.It is of great value in rational evaluation of prospecting and exploiting oil-gas resources.The mathematical model can be described as a coupled system of nonlinear partial differential equations with moving boundary values.The upwind finite difference schemes applicable to parallel arithmetic are put forward and two-dimensional and three-dimensional schemes are used to form a complete set.Some techniques,such as change of variables,calculus of variations, multiplicative commutation rule of difference operators,decomposition of high order difference operators and prior estimates,are adopted.The estimates in l~2 norm are derived to determine the error in the approximate solution.This method was already applied to the numerical simulation of migration-accumulation of oil resources.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11304074,61475042,and 11274088)the Natural Science Foundation of Hebei Province,China(Grant Nos.A2015202320 and GCC2014048)the Key Subject Construction Project of Hebei Province University,China
文摘The finite-difference time-domain method is used to simulate the optical characteristics of an in-plane switching blue phase liquid crystal display.Compared with the matrix optic methods and the refractive method,the finite-difference timedomain method,which is used to directly solve Maxwell＇s equations,can consider the lateral variation of the refractive index and obtain an accurate convergence effect.The simulation results show that e-rays and o-rays bend in different directions when the in-plane switching blue phase liquid crystal display is driven by the operating voltage.The finitedifference time-domain method should be used when the distribution of the liquid crystal in the liquid crystal display has a large lateral change.
基金supported by the National Natural Science Foundation of China(Grant Nos.61331007 and 61471105)
文摘An efficient conformal locally one-dimensional finite-difference time-domain（LOD-CFDTD） method is presented for solving two-dimensional（2D） electromagnetic（EM） scattering problems. The formulation for the 2D transverse-electric（TE） case is presented and its stability property and numerical dispersion relationship are theoretically investigated. It is shown that the introduction of irregular grids will not damage the numerical stability. Instead of the staircasing approximation, the conformal scheme is only employed to model the curve boundaries, whereas the standard Yee grids are used for the remaining regions. As the irregular grids account for a very small percentage of the total space grids, the conformal scheme has little effect on the numerical dispersion. Moreover, the proposed method, which requires fewer arithmetic operations than the alternating-direction-implicit（ADI） CFDTD method, leads to a further reduction of the CPU time. With the total-field/scattered-field（TF/SF） boundary and the perfectly matched layer（PML）, the radar cross section（RCS） of two2 D structures is calculated. The numerical examples verify the accuracy and efficiency of the proposed method.
文摘An explicit,time-dependent variable grid finite difference method is introduced and analyzed for approximating the solution of a scalar conservation law in two dimension. The scheme is stable,and the numerical solution is proved to converges to the relevant physical solution.
基金the National Key R&D Program of China(No.2018YFE0180900).
文摘The key issue in accelerating method of characteristics(MOC)transport calculations is in obtaining a completely equivalent low-order neutron transport or diffusion equation.Herein,an equivalent low-order angular flux nonlinear finite difference equation is proposed for MOC transport calculations.This method comprises three essential features:(1)the even parity discrete ordinates method is used to build a low-order angular flux nonlinear finite difference equation,and different boundary condition treatments are proposed;(2)two new defined factors,i.e.,the even parity discontinuity factor and odd parity discontinuity factor,are strictly defined to achieve equivalence between the low-order angular flux nonlinear finite difference method and MOC transport calculation;(3)the energy group and angle are decoupled to construct a symmetric linear system that is much easier to solve.The equivalence of this low-order angular flux nonlinear finite difference equation is analyzed for two-dimensional(2D)pin,2D assembly,and 2D C5G7 benchmark problems.Numerical results demonstrate that a low-order angular flux nonlinear finite difference equation that is completely equivalent to the pin-resolved transport equation is established.
文摘In the article, the fully discrete finite difference scheme for a type of nonlinear reaction-diffusion equation is established. Then the new function space is introduced and the stability problem for the finite difference scheme is discussed by means of variational approximation method in this function space. The approach used is of a simple characteristic in gaining the stability condition of the scheme.
基金Project supported partly by the Open Research Program in State Key Laboratory of Millimeter Waves of China(Grant No.K200802)partly by the National Natural Science Foundation of China(Grant No.60971122)
文摘The transmission coefficients of electromagnetic （EM） waves due to a superconductor-dielectric superlattice are numerically calculated. Shift operator finite difference time domain （SO-FDTD） method is used in the analysis. By using the SO-FDTD method, the transmission spectrum is obtained and its characteristics are investigated for different thicknesses of superconductor layers and dielectric layers, from which a stop band starting from zero frequency can be apparently observed. The relation between this low-frequency stop band and relative temperature, and also the London penetration depth at a superconductor temperature of zero degree are discussed, separately. The low-frequency stop band properties of superconductor-dielectric superlattice thus are well disclosed.
文摘The wave force exerted on vertical piles of offshore structures is the main criterion in designing them.In structures with more than one large pile,the influence of piles on each other is one of the most important issues being concerned in past researches.An efficient method for determining the interaction of piles is introduced in present research.First the wave force is calculated by the exact method using the diffraction theory,then in the finite difference numerical method the force is calculated by adding the velocity potentials of each pile and integration of pressure on their surface.The results showed that the ratio of the wave force on each of the double piles to a single pile has a damped oscillation around unity in which the amplitude of oscillation decreases with the increase in the spacing parameter.Also different wave incident directions and diffraction parameters were used and the results showed that the numerical solution has acceptable accuracy when the diffraction parameter is larger than unity.
基金supported by the National Natural Science Foundation of China(12201228,12171047)the Fundamental Research Funds for the Central Universities(3034011102)supported by National Key R&D Program of China(2020YFA0713701).
文摘In this work, we first derive the one-point large deviations principle (LDP) for both the stochastic Cahn–Hilliard equation with small noise and its spatial finite difference method (FDM). Then, we focus on giving the convergence of the one-point large deviations rate function (LDRF) of the spatial FDM, which is about the asymptotical limit of a parametric variational problem. The main idea for proving the convergence of the LDRF of the spatial FDM is via the Γ-convergence of objective functions. This relies on the qualitative analysis of skeleton equations of the original equation and the numerical method. In order to overcome the difficulty that the drift coefficient is not one-sided Lipschitz continuous, we derive the equivalent characterization of the skeleton equation of the spatial FDM and the discrete interpolation inequality to obtain the uniform boundedness of the solution to the underlying skeleton equation. These play important roles in deriving the T-convergence of objective functions.
基金Supported by Shanxi Province Natural Science Research(202203021212249)Special/Youth Foundation of Taiyuan University of Technology(2022QN101)+3 种基金National Natural Science Foundation of China(12301556)Research Project Supported by Shanxi Scholarship Council of China(2021-029)International Cooperation Base and Platform Project of Shanxi Province(202104041101019)Basic Research Plan of Shanxi Province(202203021211129)。
文摘To enhance the computational efficiency of spatio-temporally discretized phase-field models,we present a high-speed solver specifically designed for the Poisson equations,a component frequently used in the numerical computation of such models.This efficient solver employs algorithms based on discrete cosine transformations(DCT)or discrete sine transformations(DST)and is not restricted by any spatio-temporal schemes.Our proposed methodology is appropriate for a variety of phase-field models and is especially efficient when combined with flow field systems.Meanwhile,this study has conducted an extensive numerical comparison and found that employing DCT and DST techniques not only yields results comparable to those obtained via the Multigrid(MG)method,a conventional approach used in the resolution of the Poisson equations,but also enhances computational efficiency by over 90%.
