Progress with the 5D full-F continuum gyrokinetic code COGENT

COGENT is an Eulerian gyrokinetic code being developed for edge plasma modelling. The code is distinguished by the use of a high-order finite-volume (conservative) discretization combined with mapped multi-block grid technology. Our recent work is focused on the development of a 5D full-F COGENT version. A numerical algorithm utilizing locally a field-aligned multi-block coordinate system is implemented to facilitate simulations of highly anisotropic microturbulence in the presence of a strong magnetic shear. In this approach, the toroidal direction is divided into blocks such that, within each block, the cells are field-aligned and a non-matching (non-conformal) grid interface is allowed at the block boundaries. Here we report on details of the numerical implementation and present preliminary results of verification studies performed for the case of the ion temperature gradient (ITG) instability in a sheared toroidal annulus geometry.     

结构型网格分块生成技术

构造了一种适用于模拟具有复杂几何特征的流动问题的结构型网格生成方法,该方法应用CAD技术中非均匀有理B样条(NURBS)曲面表达形式,根据给定物体的典型站位截面,首先对物体进行分片表面网格的构造与重构,同时根据物体几何特征及流场特性对绕流流场进行区域分解并构造其块边界,在此基础上再采用数值方法进行相应的分块网格生成,并给出了应用所建立的方法生成的几种天地往返运输系统流场网格.     

High-order compact finite-difference methods on general overset grids

This work investigates the coupling of a very high-order finite-difference algorithm for the solution of conservation laws on general curvilinear meshes with overset-grid techniques originally developed to address complex geometric configurations. The solver portion of the algorithm is based on Padé-type compact finite-differences of up to sixth-order, with up to 10th-order filters employed to remove spurious waves generated by grid non-uniformities, boundary conditions and flow non-linearities. The overset-grid approach is utilized as both a domain-decomposition paradigm for implementation of the algorithm on massively parallel machines and as a means for handling geometric complexity in the computational domain. Two key features have been implemented in the current work; the ability of the high-order algorithm to accommodate holes cut in grids by the overset-grid approach, and the use of high-order interpolation at non-coincident grid overlaps. Several high-order/high-accuracy interpolation methods were considered, and a high-order, explicit, non-optimized Lagrangian method was found to be the most accurate and robust for this application. Several two-dimensional benchmark problems were examined to validate the interpolation methods and the overall algorithm. These included grid-to-grid interpolation of analytic test functions, the inviscid convection of a vortex, laminar flow over single- and double-cylinder configurations, and the scattering of acoustic waves from one- and three-cylinder configurations. The employment of the overset-grid techniques, coupled with high-order interpolation at overset boundaries, was found to be an effective way of employing the high-order algorithm for more complex geometries than was previously possible.     

A high-order cut-cell method for numerical simulation of hypersonic boundary-layer instability with surface roughness

Laminar-turbulent transition of hypersonic boundary layers can be affected significantly by the existence of surface roughness. Currently many important mechanisms of roughness-induced transition are not well understood. In recent years, direct numerical simulation (DNS) has been extensively applied for investigating instability and transition mechanisms of hypersonic boundary layers. Most of the past DNS studies, however, have been based on body-fitted grids for smooth surfaces without roughness. Due to complex geometry of embedded roughness, the use of body-fitted grids can be very difficult for flow with arbitrary surface roughness. In this paper, we present a new high-order cut-cell method to overcome the natural complexities in grid generation around arbitrary surface of roughness. The new method combines a non-uniform-grid finite-difference method for discrete grid points near the solid boundary and a shock-fitting method for the treatment of the bow shock. The non-uniform-grid finite-difference formulas are expressed in a general explicit form so that they can be applied to different multi-dimensional problems without any modification. The computational accuracy of new algorithms of up to O(h⁴) are tested on several one- and two-dimensional elliptic equations in irregular domains. In addition, the new method is applied to the simulation of the receptivity process of a Mach 5.92 flow over a flat plate under the combined effect of an isolated surface roughness element and surface blow and suction. A good agreement is found between the unsteady flow results and those obtained by a Linear Stability Theory (LST).     

A discrete action principle for electrodynamics and the construction of explicit symplectic integrators for linear,non-dispersive media

In this work, we derive a discrete action principle for electrodynamics that can be used to construct explicit symplectic integrators for Maxwell’s equations. Different integrators are constructed depending on the choice of discrete Lagrangian used to approximate the action. By combining discrete Lagrangians in an explicit symplectic partitioned Runge–Kutta method, an integrator capable of achieving any order of accuracy is obtained. Using the von Neumann stability analysis, we show that the integrators greatly increase the numerical stability and reduce the numerical dispersion compared to other methods. For practical purposes, we demonstrate how to implement the integrators using many features of the finite-difference time-domain method. However, our approach is also applicable to other spatial discretizations, such as those used in finite element methods. Using this implementation, numerical examples are presented that demonstrate the ability of the integrators to efficiently reduce and maintain a minimal amount of numerical dispersion, particularly when the time-step is less than the stability limit. The integrators are therefore advantageous for modeling large, inhomogeneous computational domains.     

最有效的相容剂是两嵌段还是多嵌段聚合物？

采用Monte Carlo模拟方法研究了多嵌段聚合物在A／B／嵌段聚合物三组份体系作为相容剂使用的有效性．占总体积19％的A组份在体系中为分散相．模拟结果显示了两嵌段和多嵌段聚合物在界面上的聚集行为,以及如何影响这个不相容体系的相形为．两嵌段聚合物趋于直立在相界面上,而多嵌段聚合物更容易横跨在相界面上并占据较大的界面积．从而导致多嵌段聚合物更有效的阻止体系相分离的发生．     

A Spectral Time-Domain Method for Computational Electrodynamics

Ever since its introduction by Kane Yee over forty years ago, the finite-difference time-domain (FDTD) method has been a widely-used technique for solving the time-dependent Maxwell's equations that has also inspired many other methods. This paper presents an alternative approach to these equations in the case of spatially-varying electric permittivity and/or magnetic permeability, based on Krylov subspace spectral (KSS) methods. These methods have previously been applied to the variable-coefficient heat equation and wave equation, and have demonstrated high-order accuracy, as well as stability characteristic of implicit time-stepping schemes, even though KSS methods are explicit. KSS methods for scalar equations compute each Fourier coefficient of the solution using techniques developed by Golub and Meurant for approximating elements of functions of matrices by Gaussian quadrature in the spectral, rather than physical, domain. We show how they can be generalized to coupled systems of equations, such as Maxwell's equations, by choosing appropriate basis functions that, while induced by this coupling, still allow efficient and robust computation of the Fourier coefficients of each spatial component of the electric and magnetic fields. We also discuss the application of block KSS methods to problems involving non-self-adjoint spatial differential operators, which requires a generalization of the block Lanczos algorithm of Golub and Underwood to unsymmetric matrices.     

Output error estimation for summation-by-parts finite-difference schemes

The paper develops a posteriori error estimates of integral output functionals for summation-by-parts finite-difference methods. The error estimates are based on the adjoint-weighted residual method and take advantage of a variational interpretation of summation-by-parts discretizations. The estimates are computed on a fixed grid and do not require an embedded grid or explicit interpolation operators. For smooth boundary-value problems containing first and second derivatives the error estimates converge to the exact error as the mesh is refined. The theory is verified using linear boundary-value problems and the Euler equations.     

Multi-block Lattice Boltzmann method: Extension to 3D and validation in turbulence

In this paper, a multi-block method is applied to 3D problems. In this application, interpolation schemes are carefully addressed to avoid the inconsistency when information is transferred from coarse to fine blocks. Two test cases are employed to assess information transfer scheme and accuracy improvement with respect to grid refinement.     

Split step non-paraxial finite difference method for 3D scalar wave propagation

We describe a simple and efficient wide-angle, split-step finite-difference based explicit transfer matrix method to solve the 3D scalar wave equation. The formulation is completely analytic and does not involve any numerical matrix inversion or diagonalization. It also does not use the ADI scheme. The PML boundary condition can be easily implemented with only a marginal increase in computational effort.     

Resolution estimates for simple one-dimensional flames

The accurate resolution of flame structure is critically important to the direct numerical simulation of reacting flows. While grid-resolution estimates are readily available for cold flows, premixed flames appear to have received relatively little attention. In this paper, a premixed flame characterised by single-step chemistry at moderate-to-high Zel'dovitch numbers (β) is analysed, and its structure is used to provide estimates for sufficient grid resolution. It is found that the critical region of the flame exhibits a weak inverse Zel'dovitch number dependence, and that heuristic methods of resolution estimation based on flame thickness grow relatively less meaningful as β → ∞. Resolution estimates for second- and fourth-order finite-difference schemes are provided.     

Efficient computation of compressible and incompressible flows

The combination of explicit Runge–Kutta time integration with the solution of an implicit system of equations, which in earlier work demonstrated increased efficiency in computing compressible flow on highly stretched meshes, is extended toward conditions where the free stream Mach number approaches zero. Expressing the inviscid flux Jacobians in terms of Mach number, an artificial speed of sound as in low Mach number preconditioning is introduced into the Jacobians, leading to a consistent formulation of the implicit and explicit parts of the discrete equations. Besides extension to low Mach number flows, the augmented Runge–Kutta/Implicit method allowed the admissible Courant–Friedrichs–Lewy number to be increased from O(1 0 0) to O(1 0 0 0). The implicit step introduced into the Runge–Kutta framework acts as a preconditioner which now addresses both, the stiffness in the discrete equations associated with highly stretched meshes, and the stiffness in the analytical equations associated with the disparity in the eigenvalues of the inviscid flux Jacobians. Integrated into a multigrid algorithm, the method is applied to efficiently compute different cases of inviscid flow around airfoils at various Mach numbers, and viscous turbulent airfoil flow with varying Mach and Reynolds number. Compared to well tuned conventional methods, computation times are reduced by half an order of magnitude.     

磁化等离子体光子晶体缺陷态的研究

采用磁化等离子体的分段线性电流密度卷积时域有限差分(PLCDRC-FDTD)算法研究具有单一缺陷层的一维磁化等离子体光子晶体的缺陷模特性. 从频域角度分析得到微分高斯脉冲的透射率,并讨论该光子晶体的缺陷层厚度、位置、周期常数和等离子体参数对其缺陷模的影响. 结果表明,改变位置和周期常数不会影响缺陷模的频率,改变缺陷层的厚度可以增加缺陷模数,改变等离子体参数能同时影响缺陷模的频率和峰值. 关键词： 磁化等离子体光子晶体 光子晶体 缺陷模 时域有限差分法     

A multi-block ADI finite-volume method for incompressible Navier–Stokes equations in complex geometries

An efficient second-order accurate finite-volume method is developed for a solution of the incompressible Navier–Stokes equations on complex multi-block structured curvilinear grids. Unlike in the finite-volume or finite-difference-based alternating-direction-implicit (ADI) methods, where factorization of the coordinate transformed governing equations is performed along generalized coordinate directions, in the proposed method, the discretized Cartesian form Navier–Stokes equations are factored along curvilinear grid lines. The new ADI finite-volume method is also extended for simulations on multi-block structured curvilinear grids with which complex geometries can be efficiently resolved. The numerical method is first developed for an unsteady convection–diffusion equation, then is extended for the incompressible Navier–Stokes equations. The order of accuracy and stability characteristics of the present method are analyzed in simulations of an unsteady convection–diffusion problem, decaying vortices, flow in a lid-driven cavity, flow over a circular cylinder, and turbulent flow through a planar channel. Numerical solutions predicted by the proposed ADI finite-volume method are found to be in good agreement with experimental and other numerical data, while the solutions are obtained at much lower computational cost than those required by other iterative methods without factorization. For a simulation on a grid with O(10⁵) cells, the computational time required by the present ADI-based method for a solution of momentum equations is found to be less than 20% of that required by a method employing a biconjugate-gradient-stabilized scheme.     

Compact-2D FDTD for waveguides including materials with negative dielectric permittivity, magnetic permeability and refractive index

An efficient compact-2D finite-difference time-domain method is presented for the numerical analysis of guided modes in waveguides that may include negative dielectric permittivity, negative magnetic permeability and negative refractive index materials. Both complex variable and real variable methods are given. The method is demonstrated for the analysis of channel-plasmon-polariton guided modes in triangular groves on a metal surface. The presented method can be used for a range of waveguide problems that were previously unsolvable analytically, due to complex geometries, or numerically, due to computational requirements of conventional three-dimensional finite-difference time-domain methods. A three-dimensional finite-difference time-domain algorithm that also allows analysis in the presence of bound or free electric and equivalent magnetic charges is presented and an example negative refraction demonstrates the method.     

Auxiliary differential equation: efficient implementation in the finite-difference time-domain method

An efficient algorithm for modeling dispersive media in finite-difference time-domain methods is presented. It is based on the auxiliary differential equation method for treatment of Lorentz media with an arbitrary number of relaxations. The algorithm shows excellent accuracy of second order in time and space and is efficient in both memory requirements and computational effort.     

Grid-based photogrammetry system for large scale sheet metal strain measurement

This paper presents a new approach to fast strain measurement with high accuracy for large scale sheet metal based on the surface circular grid and digital close range photogrammetry. A multi-block measuring method of discretization is implemented to archive large scale measurement. The sheet metal is separated into several blocks for respective calculating and joined together by common reference points. A surface circular searching method is presented for fast and robust 3D grid generation. A flexible bundle adjustment method is proposed for large amount 3D grid nodes reconstruction, which employs the conception of sampling points and is proved to be efficient. Furthermore, a multi-stage grid registration method is introduced to improve the accuracy of strain field by correcting the true deformation gradient tensor. A novel system is developed and performances well in actual large scale sheet metal strain measurement. Two accuracy tests confirm that the system strain measurement error is less than 0.2%.     

使用类别形状函数的多目标气动外形优化设计

在飞行器的气动外形优化设计中, 参数化方法和优化算法具有十分重要的作用, 对优化的计算时间设计空间的数学特性有着深刻的影响.类别形状函数(class and shape transformation, CST)方法是一种简洁高效的参数化方法, 但对于复杂曲面很难使用统一的CST方法进行拟合.文章首先介绍了CST方法的三维实现, 分析了其数学性质, 提出了分块CST参数化方法, 保留CST方法的特性, 实现了分块曲面之间的光滑连接.针对气动外形优化设计的复杂情况, 需要根据具体的飞行任务提出设计目标, 并处理不同目标的矛盾问题.其次采用Pareto策略自动寻找最优方案集, 并基于分块CST参数化方法遗传算法和气动力快速计算方法, 对类乘波翼身组合飞行器进行了优化设计, 并改变原有问题的设定条件优化得到了全新外形.研究结果表明分块CST方法参数少, 精度高, Pareto策略处理多目标准确有效, 是气动外形优化设计中非常有用的工具.      

基于有限差分束传播法的Mach-Zehnder型波导电光调制器的模拟设计

针对掺钛铌酸锂Mach-Zehnder型波导电光调制器的结构及折射率分布,运用有限差分束传播方法(FD BPM)进行模拟计算,取得了二维情况下的最佳设计尺寸.在此基础上,把所得的设计结果与原有的用快速傅里叶束传播方法(FFT BPM)所得的结果进行比较和分析.结果显示,运用有限差分束传播方法不仅取得了预期的结果,而且,这种方法比原有的快速傅里叶束传播方法更快捷和方便,是器件设计的理想和现实的数值模拟工具.     

Generation and circuit implementation of multi-block multidirectional grid multi-scroll chaotic attractors

Due to the dynamic characteristics of the Chua's system, multi-scroll chaotic attractors are still confined in a single block and fail to break the limit. This paper proposes an approach for generating novel multi-block multidirectional grid multi-scroll chaotic attractors that can break the limit via novel nonlinear modulating functions. According to this method, the recursion rules used to generate multi-block multidirectional grid multi-scroll attractors are mathematically obtained. The new system is autonomous; the effectiveness of this method has been verified by theoretical analysis, numerical simulation, and circuit implementation.