The method of splitting a plane-wave finite-difference time-domain （SP-FDTD） algorithm is presented for the initiation of plane-wave source in the total-field / scattered-field （TF/SF） formulation of high-order sy...The method of splitting a plane-wave finite-difference time-domain （SP-FDTD） algorithm is presented for the initiation of plane-wave source in the total-field / scattered-field （TF/SF） formulation of high-order symplectic finite- difference time-domain （SFDTD） scheme for the first time. By splitting the fields on one-dimensional grid and using the nature of numerical plane-wave in finite-difference time-domain （FDTD）, the identical dispersion relation can be obtained and proved between the one-dimensional and three-dimensional grids. An efficient plane-wave source is simulated on one-dimensional grid and a perfect match can be achieved for a plane-wave propagating at any angle forming an integer grid cell ratio. Numerical simulations show that the method is valid for SFDTD and the residual field in SF region is shrinked down to -300 dB.展开更多
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.展开更多
In this paper,a implicit difference scheme is proposed for solving the equation of one_dimension parabolic type by undetermined paameters.The stability condition is r=αΔt/Δx 2 1/2 and the truncation error is o(...In this paper,a implicit difference scheme is proposed for solving the equation of one_dimension parabolic type by undetermined paameters.The stability condition is r=αΔt/Δx 2 1/2 and the truncation error is o(Δt 4+Δx 4) It can be easily solved by double sweeping method.展开更多
Time fractional diffusion equation is usually used to describe the problems involving non-Markovian random walks. This kind of equation is obtained from the standard diffusion equation by replacing the first-order tim...Time fractional diffusion equation is usually used to describe the problems involving non-Markovian random walks. This kind of equation is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α∈(0, 1). In this paper, an implicit finite difference scheme for solving the time fractional diffusion equation with source term is presented and analyzed, where the fractional derivative is described in the Caputo sense. Stability and convergence of this scheme are rigorously established by a Fourier analysis. And using numerical experiments illustrates the accuracy and effectiveness of the scheme mentioned in this paper.展开更多
Finite-difference(FD)method is the most extensively employed numerical modeling technique.Nevertheless,when using the FD method to simulate the seismic wave propagation,the large spatial or temporal sampling interval ...Finite-difference(FD)method is the most extensively employed numerical modeling technique.Nevertheless,when using the FD method to simulate the seismic wave propagation,the large spatial or temporal sampling interval can lead to dispersion errors and numerical instability.In the FD scheme,the key factor in determining both dispersion errors and stability is the selection of the FD weights.Thus,How to obtain appropriate FD weights to guarantee a stable numerical modeling process with minimum dispersion error is critical.The FD weights computation strategies can be classified into three types based on different computational ideologies,window function strategy,optimization strategy,and Taylor expansion strategy.In this paper,we provide a comprehensive overview of these three strategies by presenting their fundamental theories.We conduct a set of comparative analyses of their strengths and weaknesses through various analysis tests and numerical modelings.According to these comparisons,we provide two potential research directions of this field:Firstly,the development of a computational strategy for FD weights that enhances stability;Secondly,obtaining FD weights that exhibit a wide bandwidth while minimizing dispersion errors.展开更多
针对COB-LED(Chip on Board-Light Emitting Diode)散热问题,文中基于二维热传导方程建立了一个可快速计算COB-LED散热器表面热分布的数学模型。为了便于模型求解,采用有限差分法求解该数学模型并选择交替方向隐格式作为其差分格式。根...针对COB-LED(Chip on Board-Light Emitting Diode)散热问题,文中基于二维热传导方程建立了一个可快速计算COB-LED散热器表面热分布的数学模型。为了便于模型求解,采用有限差分法求解该数学模型并选择交替方向隐格式作为其差分格式。根据模型中的边界条件和初始条件设计COB-LED常温点亮实验,并基于ANSYS有限元分析软件进行仿真分析。通过比较求解结果、仿真结果和实验结果验证该数学模型的合理性。结果表明,求解结果与实验结果中最高温度相对误差约23.57%,且两者的温度变化趋势一致。求解结果与仿真结果中最高温度相对误差约34.84%,且温度分布较为接近,证明了该数学模型的合理性与正确性。展开更多
目前的气体动理学统一算法(unified gas kinetic scheme,简称UGKS)在求解高速流动问题时的计算效率,难以满足求解复杂工程问题的需求.为了提高该算法的计算效率,本文对模型方程的对流项和碰撞项进行了隐式处理,并针对UGKS界面通量与演...目前的气体动理学统一算法(unified gas kinetic scheme,简称UGKS)在求解高速流动问题时的计算效率,难以满足求解复杂工程问题的需求.为了提高该算法的计算效率,本文对模型方程的对流项和碰撞项进行了隐式处理,并针对UGKS界面通量与演化时间相关的特点,引入了演化时间平均界面通量,通过对控制方程矩阵进行近似LU分解(lower-upper decomposition),实现了隐式UGKS.不同来流马赫数的圆柱绕流算例测试表明,只要演化时间选取得当,隐式方法可以得到与显式方法完全相同的结果,且计算效率可以提高1～2个量级.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.60931002 and 61101064)the Universities Natural Science Foundation of Anhui Province,China(Grant Nos.KJ2011A002 and 1108085J01)
文摘The method of splitting a plane-wave finite-difference time-domain （SP-FDTD） algorithm is presented for the initiation of plane-wave source in the total-field / scattered-field （TF/SF） formulation of high-order symplectic finite- difference time-domain （SFDTD） scheme for the first time. By splitting the fields on one-dimensional grid and using the nature of numerical plane-wave in finite-difference time-domain （FDTD）, the identical dispersion relation can be obtained and proved between the one-dimensional and three-dimensional grids. An efficient plane-wave source is simulated on one-dimensional grid and a perfect match can be achieved for a plane-wave propagating at any angle forming an integer grid cell ratio. Numerical simulations show that the method is valid for SFDTD and the residual field in SF region is shrinked down to -300 dB.
基金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.
文摘In this paper,a implicit difference scheme is proposed for solving the equation of one_dimension parabolic type by undetermined paameters.The stability condition is r=αΔt/Δx 2 1/2 and the truncation error is o(Δt 4+Δx 4) It can be easily solved by double sweeping method.
基金Supported by the Discipline Construction and Teaching Research Fund of LUTcte(20140089)
文摘Time fractional diffusion equation is usually used to describe the problems involving non-Markovian random walks. This kind of equation is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α∈(0, 1). In this paper, an implicit finite difference scheme for solving the time fractional diffusion equation with source term is presented and analyzed, where the fractional derivative is described in the Caputo sense. Stability and convergence of this scheme are rigorously established by a Fourier analysis. And using numerical experiments illustrates the accuracy and effectiveness of the scheme mentioned in this paper.
基金supported by the Marine S&T Fund of Shandong Province for Pilot National Laboratory for Marine Science and Technology (Qingdao) (No.2021QNLM020001)the Major Scientific and Technological Projects of Shandong Energy Group (No.SNKJ2022A06-R23)+2 种基金the Funds of Creative Research Groups of China (No.41821002)National Natural Science Foundation of China Outstanding Youth Science Fund Project (Overseas) (No.ZX20230152)the Major Scientific and Technological Projects of CNPC (No.ZD2019-183-003)。
文摘Finite-difference(FD)method is the most extensively employed numerical modeling technique.Nevertheless,when using the FD method to simulate the seismic wave propagation,the large spatial or temporal sampling interval can lead to dispersion errors and numerical instability.In the FD scheme,the key factor in determining both dispersion errors and stability is the selection of the FD weights.Thus,How to obtain appropriate FD weights to guarantee a stable numerical modeling process with minimum dispersion error is critical.The FD weights computation strategies can be classified into three types based on different computational ideologies,window function strategy,optimization strategy,and Taylor expansion strategy.In this paper,we provide a comprehensive overview of these three strategies by presenting their fundamental theories.We conduct a set of comparative analyses of their strengths and weaknesses through various analysis tests and numerical modelings.According to these comparisons,we provide two potential research directions of this field:Firstly,the development of a computational strategy for FD weights that enhances stability;Secondly,obtaining FD weights that exhibit a wide bandwidth while minimizing dispersion errors.
文摘针对COB-LED(Chip on Board-Light Emitting Diode)散热问题,文中基于二维热传导方程建立了一个可快速计算COB-LED散热器表面热分布的数学模型。为了便于模型求解,采用有限差分法求解该数学模型并选择交替方向隐格式作为其差分格式。根据模型中的边界条件和初始条件设计COB-LED常温点亮实验,并基于ANSYS有限元分析软件进行仿真分析。通过比较求解结果、仿真结果和实验结果验证该数学模型的合理性。结果表明,求解结果与实验结果中最高温度相对误差约23.57%,且两者的温度变化趋势一致。求解结果与仿真结果中最高温度相对误差约34.84%,且温度分布较为接近,证明了该数学模型的合理性与正确性。
文摘目前的气体动理学统一算法(unified gas kinetic scheme,简称UGKS)在求解高速流动问题时的计算效率,难以满足求解复杂工程问题的需求.为了提高该算法的计算效率,本文对模型方程的对流项和碰撞项进行了隐式处理,并针对UGKS界面通量与演化时间相关的特点,引入了演化时间平均界面通量,通过对控制方程矩阵进行近似LU分解(lower-upper decomposition),实现了隐式UGKS.不同来流马赫数的圆柱绕流算例测试表明,只要演化时间选取得当,隐式方法可以得到与显式方法完全相同的结果,且计算效率可以提高1～2个量级.
