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.展开更多
The compact implicit integration factor （cIIF） method is an efficient time discretization scheme for stiff nonlinear diffusion equations in two and three spatial dimensions. In the current work, we apply the cIIF me...The compact implicit integration factor （cIIF） method is an efficient time discretization scheme for stiff nonlinear diffusion equations in two and three spatial dimensions. In the current work, we apply the cIIF method to some complex-valued nonlinear evolutionary equations such as the nonlinear SchrSdinger （NLS） equation and the complex Ginzburg-Landau （GL） equation. Detailed algorithm formulation and practical implementation of cIIF method are performed. The numerical results indicate that this method is very accurate and efficient.展开更多
Determining the venting time of a gas trunk pipeline segment provides an important basis for formulating an emergency plan in the advent of unexpected accidents.As the natural gas venting process corresponds to the tr...Determining the venting time of a gas trunk pipeline segment provides an important basis for formulating an emergency plan in the advent of unexpected accidents.As the natural gas venting process corresponds to the transient flow,it is necessary to establish a transient hydraulic-thermal simulation model in order to determine the venting time.In this paper,based on two kinds of venting scenarios in which there is only one venting point in the venting system of a gas trunk pipeline segment—namely,where the venting point is either at one of the two ends or at the junction of two gas trunk pipeline segments—transient hydraulic-thermal simulation models are established.The models consist of gas flow governing equations,the gas state equation,gas physical property equations,initial conditions,and appropriate boundary conditions.The implicit central difference method is used to discretize the gas flow partial differential equations,and the trust-region-dogleg algorithm is used to solve the equations corresponding to each time step,in order to dynamically simulate the whole venting process.The judgment condition for the end of the venting process is that the average pressure of gas trunk pipeline segment is less than 0.11 MPa(actual pressure).Comparing the simulation results of the proposed model with those of the OLGA software and real operational data,we find that the venting time error is less than 10%.On this basis,a venting valve opening control principle is proposed,which prevents the venting noise from exceeded the specified noise value(85 d B)in the venting design of domestic gas pipeline projects.The established calculation model for venting time(dynamic simulation model)for a gas trunk pipeline segment and the proposed opening control principle of venting valve provide reference for the optimal operation of gas pipeline venting systems.展开更多
In order to eliminate Courant-Friedrich-Levy(CFL) condition restraint and improvecomputational efficiency,a new finite-difference time-domain(FDTD)method based on the alternating-direction implicit(ADI) technique is i...In order to eliminate Courant-Friedrich-Levy(CFL) condition restraint and improvecomputational efficiency,a new finite-difference time-domain(FDTD)method based on the alternating-direction implicit(ADI) technique is introduced recently.In this paper,a theoretical proof of the stabilityof the three-dimensional(3-D)ADI-FDTD method is presented.It is shown that the 3-D ADI-FDTDmethod is unconditionally stable and free from the CFL condition restraint.展开更多
针对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%,且温度分布较为接近,证明了该数学模型的合理性与正确性。展开更多
文摘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.
文摘The compact implicit integration factor （cIIF） method is an efficient time discretization scheme for stiff nonlinear diffusion equations in two and three spatial dimensions. In the current work, we apply the cIIF method to some complex-valued nonlinear evolutionary equations such as the nonlinear SchrSdinger （NLS） equation and the complex Ginzburg-Landau （GL） equation. Detailed algorithm formulation and practical implementation of cIIF method are performed. The numerical results indicate that this method is very accurate and efficient.
基金supported by the National Natural Science Foundation of China(Grant No.52174064)
文摘Determining the venting time of a gas trunk pipeline segment provides an important basis for formulating an emergency plan in the advent of unexpected accidents.As the natural gas venting process corresponds to the transient flow,it is necessary to establish a transient hydraulic-thermal simulation model in order to determine the venting time.In this paper,based on two kinds of venting scenarios in which there is only one venting point in the venting system of a gas trunk pipeline segment—namely,where the venting point is either at one of the two ends or at the junction of two gas trunk pipeline segments—transient hydraulic-thermal simulation models are established.The models consist of gas flow governing equations,the gas state equation,gas physical property equations,initial conditions,and appropriate boundary conditions.The implicit central difference method is used to discretize the gas flow partial differential equations,and the trust-region-dogleg algorithm is used to solve the equations corresponding to each time step,in order to dynamically simulate the whole venting process.The judgment condition for the end of the venting process is that the average pressure of gas trunk pipeline segment is less than 0.11 MPa(actual pressure).Comparing the simulation results of the proposed model with those of the OLGA software and real operational data,we find that the venting time error is less than 10%.On this basis,a venting valve opening control principle is proposed,which prevents the venting noise from exceeded the specified noise value(85 d B)in the venting design of domestic gas pipeline projects.The established calculation model for venting time(dynamic simulation model)for a gas trunk pipeline segment and the proposed opening control principle of venting valve provide reference for the optimal operation of gas pipeline venting systems.
基金Supported by the Specialized Research Fund for the Doctoral Program of Higher Education(No.20010614003)
文摘In order to eliminate Courant-Friedrich-Levy(CFL) condition restraint and improvecomputational efficiency,a new finite-difference time-domain(FDTD)method based on the alternating-direction implicit(ADI) technique is introduced recently.In this paper,a theoretical proof of the stabilityof the three-dimensional(3-D)ADI-FDTD method is presented.It is shown that the 3-D ADI-FDTDmethod is unconditionally stable and free from the CFL condition restraint.
文摘针对COB-LED(Chip on Board-Light Emitting Diode)散热问题,文中基于二维热传导方程建立了一个可快速计算COB-LED散热器表面热分布的数学模型。为了便于模型求解,采用有限差分法求解该数学模型并选择交替方向隐格式作为其差分格式。根据模型中的边界条件和初始条件设计COB-LED常温点亮实验,并基于ANSYS有限元分析软件进行仿真分析。通过比较求解结果、仿真结果和实验结果验证该数学模型的合理性。结果表明,求解结果与实验结果中最高温度相对误差约23.57%,且两者的温度变化趋势一致。求解结果与仿真结果中最高温度相对误差约34.84%,且温度分布较为接近,证明了该数学模型的合理性与正确性。
