Hybrid numerical methods in static and dynamic fracture mechanics

In this article, the concept of the hybrid numerical methods is clarified. On the basis of this concept, various hybrid numerical methods used in static and dynamic fracture mechanics are classified into five categories: (i) hybrid experimental–numerical methods, (ii) hybrid numerical–experimental methods, (iii) hybrid analytical–numerical methods, (iv) hybrid numerical–analytical methods, and (v) hybrid numerical–numerical methods. Features of each category of hybrid numerical method are presented with pertinent numerical results.     

串级结构Z扫描特性的分析

对非线性介质串级结构的Z扫描特性进行了分析，分别给出了解析解和数值解. 解析分析主要是基于高斯分解方法（GDM）和分布透镜方法的组合；数值分析采用了Crank-Nicholson有限差分方法和快速汉克尔变换方法(FHT). 同时，也给出了这两种数值分析的使用条件和方法. 将解析结果、数值结果和实验结果进行了比较，结果显示它们具有很好的一致性. 关键词： 串级结构 Z扫描 有限差分方法 快速汉克尔变换     

EM scattering analysis of a ferrite cylinder by a FD-FD technique with an effective absorbing boundary condition and MEI

Electromagnetic scattering problem of an arbitrarily shaped ferrite cylinder is analyzed based on the finite difference-frequency domain(FD-FD) method with an effective numerical absorbing boundary condition (ABC) and the measured equation of invariance (MEI) on the terminated boundary. Compared with the method of moments (MoM), both the numerical ABC presented in this paper and the MEI result in dramatic savings in computing time and memory requirement for electrically large objects due to the sparsity of the finite difference equation. The absorbing characteristic of this numerical ABC is demonstrated numerically. The accuracy, memory needs and CPU time of the FD-FD with the numerical ABC or the MEI and the MoM are compared and then result in some important conclusions. Besides, the RCS of some ferrite cylinders are presented.     

Direct Numerical Simulation of Three-Dimensional Richtmyer--Meshkov Instability

Direct numerical simulation （DNS） is used to study flow characteristics after interaction of a planar shock with a spherical media interface in each side of which the density is different. This interracial instability is known as the Richtmyer-Meshkov （R-M） instability. The compressible Navier-Stoke equations are discretized with group velocity control （GVC） modified fourth order accurate compact difference scheme. Three-dimensional numerical simulations are performed for R-M instability installed passing a shock through a spherical interface. Based on numerical results the characteristics of 3D R-M instability are analysed. The evaluation for distortion of the interface, the deformation of the incident shock wave and effects of refraction, reflection and diffraction are presented. The effects of the interracial instability on produced vorticity and mixing is discussed.     

A numerical study of three-dimensional liquid sloshing in tanks

A numerical model NEWTANK (Numerical Wave TANK) has been developed to study three-dimensional (3-D) non-linear liquid sloshing with broken free surfaces. The numerical model solves the spatially averaged Navier–Stokes equations, which are constructed on a non-inertial reference frame having arbitrary six degree-of-freedom (DOF) of motions, for two-phase flows. The large-eddy-simulation (LES) approach is adopted to model the turbulence effect by using the Smagorinsky sub-grid scale (SGS) closure model. The two-step projection method is employed in the numerical solutions, aided by the Bi-CGSTAB technique to solve the pressure Poisson equation for the filtered pressure field. The second-order accurate volume-of-fluid (VOF) method is used to track the distorted and broken free surface. Laboratory experiments are conducted for both 2-D and 3-D non-linear liquid sloshing in a rectangular tank. A linear analytical solution of 3-D liquid sloshing under the coupled surge and sway excitation is also developed in this study. The numerical model is first validated against the available analytical solution and experimental data for 2-D liquid sloshing of both inviscid and viscous fluids. The validation is further extended to 3-D liquid sloshing. The numerical results match with the analytical solution when the excitation amplitude is small. When the excitation amplitude is large where sloshing becomes highly non-linear, large discrepancies are developed between the numerical results and the analytical solutions, the former of which, however, agree well with the experimental data. Finally, as a demonstration, a violent liquid sloshing with broken free surfaces under six DOF excitations is simulated and discussed.     

On Physically Unacceptable Numerical Solutions Yielding Strong Chaotic Signals

Physically unacceptable chaotic numerical solutions of nonlinear circuits and systems are discussed in this paper. First, as an introduction, a simple example of a wrong choice of a numerical solver to deal with a second-order linear ordinary differential equation is presented. Then, the main result follows with the analysis of an ill-designed numerical approach to solve and analyze a particular nonlinear memristive circuit. The obtained trajectory of the numerical solution is unphysical (not acceptable), as it violates the presence of an invariant plane in the continuous systems. Such a poor outcome is then turned around, as we look at the unphysical numerical solution as a source of strong chaotic sequences. The 0–1 test for chaos and bifurcation diagrams are applied to prove that the unacceptable (from the continuous system point of view) numerical solutions are, in fact, useful chaotic sequences with possible applications in cryptography and the secure transmission of data.     

二维随机粗糙面上导体目标复合瞬态散射的混合算法

提出用时域积分方程法(TDIE)与时域基尔霍夫近似法(TDKA)的混合算法来求解二维导体随机粗糙面及其上方二维导体目标的复合瞬态散射,推导出了在TM波入射情形下显式及隐式格式的时间步进方程.将粗糙面与目标分别进行TDKA和TDIE计算,并考虑目标与粗糙面之间的耦合,对TDKA和TDIE进行混合迭代,既大大降低了粗糙面求解的复杂度,又保证了计算精度.数值算例中,考虑了角反射器(开放体)和圆柱(封闭体)两种目标,分别计算了目标表面电流响应和电场远场响应.计算结果表明,和单纯TDIE法相比,本文混合方法计算效率 关键词： 随机粗糙面 复合瞬态散射 时域积分方程法和时域基尔霍夫近似法 混合算法     

On the use of Fresnel boundary and interface conditions in radiative-transfer calculations for multilayered media

The ADO (analytical discrete ordinates) method is used to establish a concise and accurate solution for a multi-layer radiative-transfer problem with Fresnel boundary and interface conditions. A finite plane-parallel medium composed of a number (K) of sub-strata with different material properties is considered to be illuminated by isotropically incident radiation. While a general result is obtained, emphasis in the numerical work is given to computing accurately the currents and the intensities that exit each of the two exterior surfaces. Monochromatic forms (with anisotropic scattering) of the radiative-transfer equation are used, and numerical results are given for several specific cases. The complications introduced by the Fresnel boundary and interface conditions are well resolved, so that the numerical results obtained are thought to define a very high standard.     

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.     

A central Rankine–Hugoniot solver for hyperbolic conservation laws

A numerical method in which the Rankine–Hugoniot condition is enforced at the discrete level is developed. The simple format of central discretization in a finite volume method is used together with the jump condition to develop a simple and yet accurate numerical method free of Riemann solvers and complicated flux splittings. The steady discontinuities are captured accurately by this numerical method. The basic idea is to fix the coefficient of numerical dissipation based on the Rankine–Hugoniot (jump) condition. Several numerical examples for scalar and vector hyperbolic conservation laws representing the inviscid Burgers equation, the Euler equations of gas dynamics, shallow water equations and ideal MHD equations in one and two dimensions are presented which demonstrate the efficiency and accuracy of this numerical method in capturing the flow features.     

Numerical errors of diffraction computing using plane wave spectrum decomposition

In the paper the numerical determination of diffraction patterns using plane wave spectrum decomposition (PWS) is investigated. The simple formula for sampling selection for error-free numerical computation is proposed and its applicability is discussed. The usage of this formula presents practical difficulty for some diffraction problems due to required large memory load. A new multi-Fourier transform PWS (MPWS) method is elaborated which overcomes memory requirement of the PWS method. The performances of the PWS and MPWS methods are verified through extensive numerical simulations.     

Efficient numerical methods for computing ground states and dynamics of dipolar Bose–Einstein condensates

New efficient and accurate numerical methods are proposed to compute ground states and dynamics of dipolar Bose–Einstein condensates (BECs) described by a three-dimensional (3D) Gross–Pitaevskii equation (GPE) with a dipolar interaction potential. Due to the high singularity in the dipolar interaction potential, it brings significant difficulties in mathematical analysis and numerical simulations of dipolar BECs. In this paper, by decoupling the two-body dipolar interaction potential into short-range (or local) and long-range interactions (or repulsive and attractive interactions), the GPE for dipolar BECs is reformulated as a Gross–Pitaevskii–Poisson type system. Based on this new mathematical formulation, we prove rigorously existence and uniqueness as well as nonexistence of the ground states, and discuss the existence of global weak solution and finite time blow-up of the dynamics in different parameter regimes of dipolar BECs. In addition, a backward Euler sine pseudospectral method is presented for computing the ground states and a time-splitting sine pseudospectral method is proposed for computing the dynamics of dipolar BECs. Due to the adoption of new mathematical formulation, our new numerical methods avoid evaluating integrals with high singularity and thus they are more efficient and accurate than those numerical methods currently used in the literatures for solving the problem. Extensive numerical examples in 3D are reported to demonstrate the efficiency and accuracy of our new numerical methods for computing the ground states and dynamics of dipolar BECs.     

EXPRESSIONS FOR PREDICTING FUNDAMENTAL NATURAL FREQUENCIES OF NON-CYLINDRICAL HELICAL SPRINGS

Numerical and analytical studies are performed for the free vibration analysis of non-cylindrical (conical, barrel and hyperboloidal types) helical springs. The stiffness matrix method is used in the numerical analysis. A total of 12 degrees of freedom (six displacements and six rotations) is described for an element. The exact element stiffness matrix and the exact concentrated element inertia matrix are used in the formulation. The rotary inertia, the shear and extensional deformation effects are considered in the analysis. Comparison of the numerical results with the reported results obtained numerically and experimentally gives satisfactory values. After verification of the numerical frequencies, the non-dimensional fundamental frequencies of fixed-fixed non-cylindrical helical springs with circular section are expressed in a simple formula with a maximum absolute relative error of 5% using those numerical values for the constant helix pitch angles (5°, 10°, and 15°). These expressions restricted to the fundamental frequencies are also verified with ANSYS results.     

A lattice Bhatnagar--Gross--Krook model for a class of the generalized Burgers equations

A new lattice Bhatnagar--Gross--Krook (LBGK) model for a class of the generalized Burgers equations is proposed. It is a general LBGK model for nonlinear Burgers equations with source term in arbitrary dimensional space. The linear stability of the model is also studied. The model is numerically tested for three problems in different dimensional space, and the numerical results are compared with either analytic solutions or numerical results obtained by other methods. Satisfactory results are obtained by the numerical simulations.     

混凝土HJC模型抗侵彻参数敏感性数值模拟研究

 混凝土是一种应用广泛的复合材料,其动态性能及抗侵彻毁伤效应的研究越来越依赖于数值模拟,其中混凝土HJC模型是应用最广泛的本构模型之一。在介绍HJC模型的基础上,应用LS-DYNA数值模拟方法进行HJC本构模型抗侵彻性能的参数敏感性研究,通过大量数值模拟得出HJC模型的抗侵彻敏感参数有：无侧限单轴抗压强度f_c、压实点压力p_lock、特征化粘性强度A、特征化压力硬化因子B和压力硬化指数N。在混凝土的抗侵彻数值模拟过程中,重点是确定敏感参数的取值,以提高数值模拟的精度和质量,做到合理地选择试验方法和安排试验量。     

Design and numerical simulation of novel DBRs

In this paper, a numerical simulation of the traditional graded distributed Bragg reflector (DBR) and a design of the novel DBR with short period superlattices (SPSs DBR) used by vertical cavity surface emitting laser (VCSEL) are reported. First, the optical characteristic matrix of the graded DBRs is derived using the theories of thin film optics. Second, its reflective spectrum is numerical simulated and it is found that the simulative results are similar with the experimental data. The difference of the cavity mode position between the experimental and simulative data is discussed. Finally, based on the simulative results of graded DBR, a novel DBR with 4.5-pair GaAs/AlAs SPSs is designed, and its reflective spectrum is numerical simulated and analyzed.     

Numerical simulation of partially coherent broadband optical imaging using the finite-difference time-domain method

Rigorous numerical modeling of optical systems has attracted interest in diverse research areas ranging from biophotonics to photolithography. We report the full-vector electromagnetic numerical simulation of a broadband optical imaging system with partially coherent and unpolarized illumination. The scattering of light from the sample is calculated using the finite-difference time-domain (FDTD) numerical method. Geometrical optics principles are applied to the scattered light to obtain the intensity distribution at the image plane. Multilayered object spaces are also supported by our algorithm. For the first time, numerical FDTD calculations are directly compared to and shown to agree well with broadband experimental microscopy results.     

Numerical analysis of a CO2 SFUR laser cavity

Summary In this paper we present a numerical study of the SFUR (self-filtering unstable resonator) cavity resonant modes and numerical results are compared with the experimental data. The numerical code developed for analysing the SFUR cavity will also allow for studying a wider class of resonators, and is therefore a very useful tool in designing new laser cavities.     

Prediction-step staggered-in-time FDTD: An efficient numerical scheme to solve the linearised equations of fluid dynamics in outdoor sound propagation

The finite-difference time-domain (FDTD) method to solve the linearised equations of fluid dynamics has shown to be very powerful and useful in outdoor sound propagation. Practical applications are however limited due to the large need for computational resources. The numerical discretisation influences computational efficiency to an important degree. In this paper, some possible ways to discretise temporal derivatives are studied. Two obvious ways of time-discretisation namely staggered-in-time (SIT) and a simple collocated-in-time (CIT) scheme are compared to the prediction-step staggered-in-time (PSIT) scheme. The latter is intended to be used for the calculation of sound propagation in the typical low wind speeds encountered in the outdoor environment at low heights above the earth’s surface. It was shown that the PSIT scheme is more stable than the SIT scheme, so practical calculations are possible. Computational efficiency is increased to an important degree compared to the CIT scheme. The numerical accuracy (more precisely the amplitude error) of the PSIT scheme is an important improvement upon SIT. The CIT scheme on the other hand conserves amplitude better. The amplitude error becomes larger with increasing wind speed because of some simplifications during the numerical discretisation. In low wind speeds, the PSIT algorithm can serve as an interesting compromise between numerical accuracy and the required amount of computing power.     

An algorithm for LES of premixed compressible flows using the Conditional Moment Closure model

A novel numerical method has been developed to couple a recent high order accurate fully compressible upwind method with the Conditional Moment Closure combustion model. The governing equations, turbulence modelling and numerical methods are presented in full. The new numerical method is validated against direct numerical simulation (DNS) data for a lean premixed methane slot burner. Although the modelling approaches are based on non-premixed flames and hence not expected to be valid for a wide range of premixed flames, the predicted flame is just 10% longer than that in the DNS and excellent agreement of mean mass fractions, conditional mass fractions and temperature is demonstrated. This new numerical method provides a very useful framework for future application of CMC to premixed as well as non-premixed combustion.