Full-vectorial finite-difference beam propagation method based on the modified alternating direction implicit method

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.

Analysis of regular and chaotic dynamics of the Euler-Bernoulli beams using finite difference and finite element methods

Chaotic vibrations of flexible non-linear Euler-Bernoulli beams subjected to harmonic load and with various boundary conditions（symmetric and non-symmetric）are studied in this work.Reliability of the obtained results is verified by the finite difference method（FDM）and the finite element method（FEM）with the Bubnov-Galerkin approximation for various boundary conditions and various dynamic regimes（regular and non-regular）.The influence of boundary conditions on the Euler-Bernoulli beams dynamics is studied mainly,dynamic behavior vs.control parameters { ωp,q0 } is reported,and scenarios of the system transition into chaos are illustrated.

An efficient locally one-dimensional finite-difference time-domain method based on the conformal scheme

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.

The dimension split element-free Galerkin method for three-dimensional potential problems

This paper presents the dimension split element-free Galerkin （DSEFG） method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-dimensional potential problem can be transformed into a series of two-dimensional problems. For these two-dimensional problems, the improved moving least-squares （IMLS） approximation is applied to construct the shape function, which uses an orthogonal function system with a weight function as the basis functions. The Galerkin weak form is applied to obtain a discretized system equation, and the penalty method is employed to impose the essential boundary condition. The finite difference method is selected in the splitting direction. For the purposes of demonstration, some selected numerical examples are solved using the DSEFG method. The convergence study and error analysis of the DSEFG method are presented. The numerical examples show that the DSEFG method has greater computational precision and computational efficiency than the IEFG method.

Numerical simulation of standing wave with 3D predictor-corrector finite difference method for potential flow equations

A three-dimensional （3D） predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is mapped onto a fixed cubic tank through the proper coordinate transform schemes. The cubic tank is distributed by the staggered meshgrid, and the staggered meshgrid is used to denote the variables of the flow field. The predictor-corrector finite difference method is given to develop the difference equa- tions of the dynamic boundary equation and kinematic boundary equation. Experimental results show that, using the finite difference method of the predictor-corrector scheme, the numerical solutions agree well with the published results. The wave profiles of the standing wave with different amplitudes and wave lengths are studied. The numerical solutions are also analyzed and presented graphically.

An improved three-dimensional full-vectorial finite-difference imaginary-distance beam propagation method

A modified alternating direction implicit approach is proposed to discretize the three-dimensional full-vectorial beam propagation method （3D-FV-BPM） formulation along the longitudinal direction. The cross-coupling terms （CCTs） are neglected at the first substep, and then double used at the second substep. The order of two substeps is reversed for each transverse electric field component so that the CCTs are always expressed in an implicit form, thus the calculation is efficient and stable. Based on the multinomial interpolation, a universal finite difference scheme with a high accuracy is developed to approximate the 3D-FV-BPM formulation along the transverse directions, in which the discontinuities of the normal components of the electric field across the abrupt dielectric interfaces are taken into account and can be applied to both uniform and non-uniform grids. The corresponding imaginary-distance procedure is first applied to a buried rectangular and a GaAs-based deeply-etched rib waveguide. The field patterns and the normalized propagation constants of the fundamental and the first order modes are presented and the hybrid nature of the full-vectorial guided-modes is demonstrated, which shows the validity and utility of the present approach. Then the modal characteristics of the deeply- and shallow-etched rib waveguides based on the InGaAsp/InGaAsP strained multiple quantum wells in InP substrate are investigated in detail. The results are necessary for modeling and the design of the planar lightwave circuits or photonic integrated circuits based on these waveguides.

DYNAMIC ANALYSIS OF ARREST OF BUCKLE PROPAGATION ON A BEAM ON A NONLINEAR ELASTIC FOUNDATION BY FEM

Based on the dynamic governing equation of propagating buckle on a beam on a nonlinear elastic foundation, this paper deals with an important problem of buckle arrest by combining the FEM with a time integration technique. A new conclusion completely different from that by the quasi-static analysis about the buckle arrestor design is drawn. This shows that the inertia of the beam cannot be ignored in the analysis under consideration, especially when the buckle propagation is suddenly stopped by the arrestors.

Three-Dimensional Numerical Simulation of Stably Stratified Flows over a Two-Dimensional Hill

Stably stratified flows over a two-dimensional hill are investigated in a channel of finite depth using a three-dimensional direct numerical simulation (DNS). The present study follows onto our previous two-dimensional DNS studies of stably stratified flows over a hill in a channel of finite depth and provides a more realistic simulation of atmospheric flows than our previous studies. A hill with a constant cross-section in the spanwise (y) direction is placed in a 3-D computational domain. As in the previous 2-D simulations, to avoid the effect of the ground boundary layer that develops upstream of the hill, no-slip conditions are imposed only on the hill surface and the surface downstream of the hill;slip conditions are imposed on the surface upstream of the hill. The simulated 3-D flows are discussed by comparing them to the simulated 2-D flows with a focus on the effect of the stable stratification on the non-periodic separation and reattachment of the flow behind the hill. In neutral (K = 0, where K is a non-dimensional stability parameter) and weakly stable (K = 0.8) conditions, 3-D flows over a hill differ clearly from 2-D flows over a hill mainly because of the three-dimensionality of the flow, that is the development of a spanwise flow component in the 3-D flows. In highly stable conditions (K = 1, 1.3), long-wavelength lee waves develop downstream of the hill in both 2-D and 3-D flows, and the behaviors of the 2-D and 3-D flows are similar in the vicinity of the hill. In other words, the spanwise component of the 3-D flows is strongly suppressed in highly stable conditions, and the flow in the vicinity of the hill becomes approximately two-dimensional in the x and z directions.

Alternating Group Explicit Iterative Methods for One-Dimensional Advection-Diffusion Equation

The finite difference method such as alternating group iterative methods is useful in numerical method for evolutionary equations and this is the standard approach taken in this paper. Alternating group explicit (AGE) iterative methods for one-dimensional convection diffusion equations problems are given. The stability and convergence are analyzed by the linear method. Numerical results of the model problem are taken. Known test problems have been studied to demonstrate the accuracy of the method. Numerical results show that the behavior of the method with emphasis on treatment of boundary conditions is valuable.

基于DNN的煤矿富水区探测反演方法研究

提出了一种基于DNN的煤矿富水区探测反演算法,该算法可以快速准确地实现煤矿富水区二维分布模型的重建。首先,利用时域有限差分方法(FDTD)获得不同分布模型的数值解;随后,依据数据样本搭建网络框架,网络的输入主要为电场分量,输出为相应的模型电导率参数。通过对神经网络进行训练,得到网络的最优系数,随后对富水区分布进行反演预测;结果表明:DNN算法在单个小目标异常体反演中,可以有效克服BP神经网络模型失效的问题,且对于多目标异常体的反演效果更加准确。另外,相同数据集下,DNN的训练耗时与预测耗时也少于BP神经网络。实验结果表明,该算法可以有效提高煤矿富水区探测效率。

Three-dimensional beam propagation method based on the variable transformed Galerkin's method

A novel three-dimensional beam propagation method (BPM) based on the variable transformed Galerkin's method is introduced for simulating optical field propagation in three-dimensional dielectric structures. The infinite Cartesian x-y plane is mapped into a unit square by a tangent-type function transformation. Consequently, the infinite region problem is converted into the finite region problem. Thus, the boundary truncation is eliminated and the calculation accuracy is promoted. The three-dimensional BPM basic equation is reduced to a set of first-order ordinary differential equations through sinusoidal basis function, which fits arbitrary cladding optical waveguide, then direct solution of the resulting equations by means of the Runge-Kutta method. In addition, the calculation is efficient due to the small matrix derived from the present technique. Both z-invariant and z-variant examples are considered to test both the accuracy and utility of this approach.

二维复杂弹性空腔的边光滑有限元建模及分析

在弹性空腔结构上敷设被动约束层阻尼(passive constraint layer damping, PCLD)可达到减振降噪的效果,针对这类复杂结构建立了二维复合弹性空腔的边光滑有限元耦合动力学模型。其中,PCLD结构采用两节点四自由度的PCLD梁单元,声场采用边光滑有限元模型。以二维全敷设复合矩形空腔模型为数值算例,以精细网格下的有限元法结果作为参考解,对比研究了在相同背景网格下,边光滑有限元法和有限元法的频响结果,发现前者更接近参考解,说明同样的计算成本下,边光滑有限元法具有更高的精确性,特别是在中频计算中。最后,分析了PCLD结构对某汽车驾驶舱的降噪效果,以及黏弹层和约束层厚度参数的影响规律,发现黏弹层厚度增大,可一定程度上降低空腔噪声,而约束层厚度增大,并不能在整个频段得到很好的降噪效果。

三维表面裂纹扩展轨迹与数值仿真研究

针对三维表面裂纹扩展形态和轨迹难以预测的特点,基于ANSYS有限元分析软件,利用三维裂纹扩展仿真方法,开发三维裂纹扩展程序,研究典型的三维表面单裂纹与三维非等大共面表面双裂纹扩展轨迹,实现了任意三维多裂纹扩展轨迹的数值模拟。主要研究内容与结论如下:针对三维表面单裂纹模型,当初始裂纹形状c/a>1时,最深处的应力强度因子值大于自由表面处应力强度因子值,随着裂纹的不断扩展,前缘会渐渐趋于稳定的圆形。而对于三维非等大共面表面双裂纹,较大的裂纹扩展速率大于较小的裂纹。开始时2条裂纹均沿光滑的样条曲线扩展,后来受到另一条裂纹的影响,在彼此接近处,由于应力放大作用,此部位的应力强度因子变大,扩展速率也会高于裂纹前缘其他部位。

基于电流密度模型的自然电场二维正反演研究

为推动自然电位数据解释,本文基于二维电流密度模型开展了自然电位二维正反演研究。首先,基于有限差分法实现了二维自然电位正演计算,并通过与解析解对比,其最大误差不超过2%,验证了正演算法的有效性。在此基础上,建立了关于电流密度的线性正则化反演算法,该正则化算法通过合理选用正则化因子对数据噪声具有一定适应性。通过含5%噪声的合成数据反演测试结果显示,预测数据与合成数据拟合较好。实测数据的反演测试结果显示,反演预测数据与实测数据拟合较好,验证了反演算法的有效性。该反演算法能有效恢复地下电流密度分布,对地下场源的进一步定位提供了有效信息,将有利于进一步提升自然电位数据的解释效果。

Dynamic Analysis and Numerical Simulation of a Kind of Arrestor Arresting Buckle Propagation by FEM

For exact estimation of efficiency of a buckle arrestor, it is necessary to take the effect of structural inertia into account in the analysis of buckle propagation on elastic structures after meeting arrestors. Under this consideration, this paper deals with the dynamics of buckle arrest and its numerical simulation on the basis of the beam system model used by Chater and Hutchinson (1983). The FEM combined with an improving are-length control method is adopted to solve the dynamic equations describing the arresting of buckle propagation. A new group of parameters for arrestor design which differs greatly from that by the quasi-static analysis is obtained. The present results support the conclusion that the inertia of the beam cannot be neglected in such analysis.

Conjugate Gradient Method to Solve Fluid Structure Interaction Problem

In this paper, we propose a method to solve coupled problem. Our computational method is mainly based on conjugate gradient algorithm. We use finite difference method for the structure and finite element method for the fluid. Conjugate gradient method gives suitable numerical results according to some papers.

A Fast Element-Free Galerkin Method for 3D Elasticity Problems

In this paper,a fast element-free Galerkin(FEFG)method for three-dimensional(3D)elasticity problems is established.The FEFG method is a combination of the improved element-free Galerkin(IEFG)method and the dimension splitting method(DSM).By using the DSM,a 3D problem is converted to a series of 2D ones,and the IEFG method with a weighted orthogonal function as the basis function and the cubic spline function as the weight function is applied to simulate these 2D problems.The essential boundary conditions are treated by the penalty method.The splitting direction uses the finite difference method(FDM),which can combine these 2D problems into a discrete system.Finally,the system equation of the 3D elasticity problem is obtained.Some specific numerical problems are provided to illustrate the effectiveness and advantages of the FEFG method for 3D elasticity by comparing the results of the FEFG method with those of the IEFG method.The convergence and relative error norm of the FEFG method for elasticity are also studied.

ON THE FINITE DIFFERENCE METHODS FOR TWO DIMENSIONAL TURBULENCE BOUNDARY LAYER

In this paper,the author first establishes the general finite difference formula for the governing equations of the turbulent average velocities in a steady two dimensional incompressible fluid boundary layer-inner layer.Next, three key parameters of the difference scheme are determined respectively by several simple flow models with known analytical solutions.Finally a special five points difference system is given and its application value is showed by a numerical example for the vertical velocity distribution in an Ekman's layer.

FDTD Simulation of Three Photon Absorption and Realization of NAND Gate with GaAs Wire Waveguide

GaAs has high three photon absorption (3PA) co-efficient at mid-infrared wavelength like2.2mm and waveguides can be formed with this material like silicon nano-wires. It is shown that three-photon-absorption in GaAs wire waveguide can be utilized to form NAND gate. Three-photon-absorption is incorporated in one-dimensional Finite Difference Time Domain (FDTD) equations. The evolution of a probe pulse under the influence of a pump pulse through crossabsorption in a waveguide is investigated using FDTD simulation, where the dominant process is nonlinear three-photon-absorption. Output probe power dependence on input pump power shows that GaAs waveguide NAND gate has higher extinction ratio in comparison to NAND gate using two-photon-absorption in silicon waveguide.

径向井立体压裂裂缝扩展数值模拟

基于有限元-无网格法,考虑压裂液流动与岩石基质变形的耦合作用,建立了径向井压裂裂缝扩展数值模型,模拟了径向井压裂裂缝形态,定量表征了径向井对裂缝的诱导作用,揭示了方位夹角、水平主应力差与岩石基质渗透率等对裂缝形态演化的影响机制。研究表明,同层双分支径向井压裂时径向井间存在挤压作用,当径向井沿最小水平主应力方向对称分布且方位夹角大于15°时,挤压作用减小了径向井的裂缝引导长度,当径向井沿最大水平主应力方向对称分布时,挤压作用增大了径向井的裂缝引导长度。裂缝形态由径向井的纠偏作用、径向井间的挤压作用与最大水平主应力的偏转作用共同控制。当径向井沿最小水平主应力方向对称分布时,径向井的裂缝引导长度随着方位夹角的增大逐渐减小,当径向井沿最大水平主应力方向对称分布时,径向井的裂缝引导长度随着方位夹角的增大先减小后增大。径向井的裂缝引导长度随水平主应力差的增大而减小。岩石基质渗透率增大,径向井周围基质的孔隙压力增大,对裂缝的诱导作用增强。