Time-Domain Numerical Solutions of Maxwell Interface Problems with Discontinuous Electromagnetic Waves

This paper is devoted to time domain numerical solutions of two-dimensional (2D) material interface problems governed by the transverse magnetic (TM) and transverse electric (TE) Maxwell's equations with discontinuous electromagnetic solutions. Due to the discontinuity in wave solutions across the interface, the usual numerical methods will converge slowly or even fail to converge. This calls for the development of advanced interface treatments for popular Maxwell solvers. We will investigate such interface treatments by considering two typical Maxwell solvers – one based on collocation formulation and the other based on Galerkin formulation. To restore the accuracy reduction of the collocation finite-difference time-domain (FDTD) algorithm near an interface, the physical jump conditions relating discontinuous wave solutions on both sides of the interface must be rigorously enforced. For this purpose, a novel matched interface and boundary (MIB) scheme is proposed in this work, in which new jump conditions are derived so that the discontinuous and staggered features of electric and magnetic field components can be accommodated. The resulting MIB time-domain (MIBTD) scheme satisfies the jump conditions locally and suppresses the staircase approximation errors completely over the Yee lattices. In the discontinuous Galerkin time-domain (DGTD) algorithm – a popular Galerkin Maxwell solver, a proper numerical flux can be designed to accurately capture the jumps in the electromagnetic waves across the interface and automatically preserves the discontinuity in the explicit time integration. The DGTD solution to Maxwell interface problems is explored in this work, by considering a nodal based high order discontinuous Galerkin method. In benchmark TM and TE tests with analytical solutions, both MIBTD and DGTD schemes achieve the second order of accuracy in solving circular interfaces. In comparison, the numerical convergence of the MIBTD method is slightly more uniform, while the DGTD method is more flexible and robust.     

初始纵横向荷载效应对薄板自振频率的影响分析The analysis for the effect of both initial in-plane and transverse loads on the natural frequency of plates

首先基于能量变分原理,给出了同时考虑初始纵向、横向荷载效应情况下板的应变能表达、动力平衡微分方程及单元刚度矩阵。再以动力平衡微分方程为基础,运用伽辽金法解得四种典型板同时考虑初始纵向、横向荷载效应的基频近似解,并运用瑞利法解得简支矩形板考虑初始纵向、横向荷载效应的前三阶频率近似解。然后相互验证了考虑初始纵向、横向荷载效应的板单元刚度矩阵和频率近似解的正确性,并进一步分析了初始纵向、横向荷载及相关因素对板自振频率的影响。结果表明,考虑初始纵向、横向荷载效应后,板的自振频率主要受初始纵向荷载、初始横向荷载、板的厚度及边界条件等因素的影响;初始纵向、横向荷载效应对板的基频影响明显于高阶频率;初始纵向、横向荷载对板的自振频率影响分别呈线性和抛物线规律。     

横观各向同性油气藏水力压裂裂纹扩展规律研究

针对横观各向同性与各向同性油气藏水力压裂裂纹扩展的差异性,基于扩展有限元法建立水力压裂力学模型,通过ABAQUS子程序开发了各向同性和横观各向同性岩体的起裂判据。在各向同性岩体数值模拟结果与解析解以及现场压裂典型曲线对比吻合的基础上,得到了包含层理构造的横观各向同性岩体水力压裂过程中裂纹扩展规律。层理类岩体水力压裂的裂缝扩展方向由地应力状态、层理方向以及岩体与层理界面抗拉强度共同决定;水力压裂过程中,注水压力在裂纹尖端产生应力集中,层理面法向分量先达到界面抗拉强度时,裂纹沿层理方向开裂,反之裂纹沿垂直最小地应力的方向扩展;裂纹扩展速度随层理抗拉强度的增加而降低;由于地层的滤失,随压裂液的注入,裂纹长宽尺度增长速率降低。     

Effect of rotary inertia on stability of axially accelerating viscoelastic Rayleigh beams

The dynamic stability of axially moving viscoelastic Rayleigh beams is presented. The governing equation and simple support boundary condition are derived with the extended Hamilton’s principle. The viscoelastic material of the beams is described as the Kelvin constitutive relationship involving the total time derivative. The axial tension is considered to vary longitudinally. The natural frequencies and solvability condition are obtained in the multi-scale process. It is of interest to investigate the summation parametric resonance and principal parametric resonance by using the Routh-Hurwitz criterion to obtain the stability condition. Numerical examples show the effects of viscosity coefficients, mean speed, beam stiffness, and rotary inertia factor on the summation parametric resonance and principle parametric resonance. The differential quadrature method (DQM) is used to validate the value of the stability boundary in the principle parametric resonance for the first two modes.     

Valley Hall effect and Nernst effect in strain engineered graphene

We theoretically predict the existence of tunneling valley Hall effect and Nernst effect in the normal/strain/normal graphene junctions, where a strained graphene is sandwiched by two normal graphene electrodes. By applying an electric bias a pure transverse valley Hall current with longitudinal charge current is generated. If the system is driven by a temperature bias, a valley Nernst effect is observed, where a pure transverse valley current without charge current propagates. Furthermore, the transverse valley current can be modulated by the Fermi energy and crystallographic orientation. When the magnetic field is further considered, we obtain a fully valley-polarized current. It is expected these features may be helpful in the design of the controllable valleytronic devices.     

不同横纵牵引角度下拖轮艏拖桩结构强度分析

以3675kW拖轮为例, 利用有限元软件MSC.Patran/Nastran对其艏拖桩结构进行结构强度分析. 通过计算得出该艏拖桩结构强度能够满足法国船级社规范要求, 并得到了不同牵引角度下的应力分布规律, 为今后拖轮结构设计及营运提供了参考.     

Fraction‐Dimensional Accessible Solitons in a Parity‐Time Symmetric Potential

By using the modified Snyder‐Mitchell (MSM) model, which can describe the propagation of a paraxial beam in fractional dimensions (FDs), we find the exact "accessible soliton” solutions in the strongly nonlocal nonlinear media with a self‐consistent parity‐time (PT) symmetric complex potential. The exact solutions are constructed with the help of two special functions: the complex Gegenbauer and the generalized Laguerre polynomials in polar coordinates, parametrized by two nonnegative integer indices ‐ the radial and azimuthal mode numbers (n,m), and the beam modulation depth. By the choice of different soliton parameters, the intensity and angular profiles display symmetric and asymmetric structures. We believe that it is important to explore the MSM model in FDs and PT‐symmetric potentials, for a better understanding of nonlinear FD physical phenomena. Different physical systems in which the model might be of relevance are briefly discussed.     

Analytical solution of transverse oscillation in cyclotron using LP method

We have carried out an approximate analytical solution to precisely consider the influence of magnetic field on the transverse oscillation of particles in a cyclotron.The differential equations of transverse oscillation are solved from the Lindstedt-Poincare method.After careful deduction,accurate first-order analytic solutions are obtained.The analytical solutions are applied to the magnetic field from an isochronous cyclotron with four spiral sectors.The accuracy of these analytical solutions is verified and confirmed from comparison with a numerical method.Finally,we discussed the transverse oscillation at v_0=N/2,using the same analytical solution.     

分流技术强化超声速混合及燃烧的数值模拟

本文提出一个利用分流技术强化超声速横侧射流混合及燃烧的新方法,即通过一个细小通道连通燃料喷口上下游,上游高压流体部分分流到下游,从底部抬高燃料穿透深度.应用修正可压缩性的 k-ε湍流模型,对不同喷射压力下冷态及燃烧流场进行数值模拟.有无采用分流技术的结果对比表明,分流技术能有效增大混合效率和燃烧效率,分流孔位置是一重要影响因素.观察到两类强化混合及燃烧的机制,合理布置分流孔位置可使燃烧效率相对提高 25% 以上.     

与质点连接的杆横振动的频率特性

推导了一端固定、一端与质点连接的弹性杆横振动的频率特性方程,并详细讨论了几种质量比情形下的特征频率。