High-order parabolic beam approximation for aero-optics

The parabolic beam equations are solved using high-order compact differences for the Laplacians and Runge–Kutta integration along the beam path. The solution method is verified by comparison to analytical solutions for apertured beams and both constant and complex index of refraction. An adaptive 4th-order Runge–Kutta using an embedded 2nd-order method is presented that has demonstrated itself to be very robust. For apertured beams, the results show that the method fails to capture near aperture effects due to a violation of the paraxial approximation in that region. Initial results indicate that the problem appears to be correctable by successive approximations. A preliminary assessment of the effect of turbulent scales is undertaken using high-order Lagrangian interpolation. The results show that while high fidelity methods are necessary to accurately capture the large scale flow structure, the method may not require the same level of fidelity in sampling the density for the index of refraction. The solution is used to calculate a phase difference that is directly compared with that commonly calculated via the optical path difference. Propagation through a supersonic boundary layer shows that for longer wavelengths, the traditional method to calculate the optical path is less accurate than for shorter wavelengths. While unlikely to supplant more traditional methods for most aero-optics applications, the current method can be used to give a quantitative assessment of the other methods as well as being amenable to the addition of more physics.     

High-order FDTD methods for transverse electromagnetic systems in dispersive inhomogeneous media

This Letter introduces a novel finite-difference time-domain (FDTD) formulation for solving transverse electromagnetic systems in dispersive media. Based on the auxiliary differential equation approach, the Debye dispersion model is coupled with Maxwell's equations to derive a supplementary ordinary differential equation for describing the regularity changes in electromagnetic fields at the dispersive interface. The resulting time-dependent jump conditions are rigorously enforced in the FDTD discretization by means of the matched interface and boundary scheme. High-order convergences are numerically achieved for the first time in the literature in the FDTD simulations of dispersive inhomogeneous media.     

High-order convergent expansions for quantum many particle systems

We review recent advances in the calculations of high-order convergent expansions for quantum many-particle systems. Calculations for ground state properties, including correlation functions and static susceptibilities, for spin models as well as for models of many fermions, such as Hubbard and Kondo models, are discussed. A historical perspective to the subject is provided. Recently important technical advances have been made in perturbative calculations of the excitation spectra of quantum many-particle systems, which enable the calculation of these spectra to high orders. The method, along with its applications, are explained. Fairly comprehensive, though simplified, algorithms for generating lists of relevant clusters, their lattice embeddings and subclusters are presented. The perturbative recursion relations and their computer implementation are also discussed in detail. A compilation is made of various series expansion studies that have been carried out for condensed matter problems. The scope and limitations of these methods are explained, and several open problems are noted.     

High-order explicit-implicit numerical methods for nonlinear anomalous diffusion equations

In this paper, the high-order finite difference/element methods for the nonlinear anomalous diffusion equations of subdiffusion and superdiffusion are developed, where the high-order finite difference methods are used to approximate the time-fractional derivatives and the finite element methods are used in the spatial domain. The stability and error estimates are proved for both cases of superdiffusion and subdiffusion. Numerical examples are provided to confirm the theoretical analysis.     

ADI spectral collocation methods for parabolic problems

We discuss the Crank–Nicolson and Laplace modified alternating direction implicit Legendre and Chebyshev spectral collocation methods for a linear, variable coefficient, parabolic initial-boundary value problem on a rectangular domain with the solution subject to non-zero Dirichlet boundary conditions. The discretization of the problems by the above methods yields matrices which possess banded structures. This along with the use of fast Fourier transforms makes the cost of one step of each of the Chebyshev spectral collocation methods proportional, except for a logarithmic term, to the number of the unknowns. We present the convergence analysis for the Legendre spectral collocation methods in the special case of the heat equation. Using numerical tests, we demonstrate the second order accuracy in time of the Chebyshev spectral collocation methods for general linear variable coefficient parabolic problems.     

High-order implicit hybridizable discontinuous Galerkin methods for acoustics and elastodynamics

We present a class of hybridizable discontinuous Galerkin (HDG) methods for the numerical simulation of wave phenomena in acoustics and elastodynamics. The methods are fully implicit and high-order accurate in both space and time, yet computationally attractive owing to their following distinctive features. First, they reduce the globally coupled unknowns to the approximate trace of the velocity, which is defined on the element faces and single-valued, thereby leading to a significant saving in the computational cost. In addition, all the approximate variables (including the approximate velocity and gradient) converge with the optimal order of k + 1 in the L²-norm, when polynomials of degree k ? 0 are used to represent the numerical solution and when the time-stepping method is accurate with order k + 1. When the time-stepping method is of order k + 2, superconvergence properties allows us, by means of local postprocessing, to obtain better, yet inexpensive approximations of the displacement and velocity at any time levels for which an enhanced accuracy is required. In particular, the new approximations converge with order k + 2 in the L²-norm when k ? 1 for both acoustics and elastodynamics. Extensive numerical results are provided to illustrate these distinctive features.     

High-order finite-volume methods for the shallow-water equations on the sphere

This paper presents a third-order and fourth-order finite-volume method for solving the shallow-water equations on a non-orthogonal equiangular cubed-sphere grid. Such a grid is built upon an inflated cube placed inside a sphere and provides an almost uniform grid point distribution. The numerical schemes are based on a high-order variant of the Monotone Upstream-centered Schemes for Conservation Laws (MUSCL) pioneered by van Leer. In each cell the reconstructed left and right states are either obtained via a dimension-split piecewise-parabolic method or a piecewise-cubic reconstruction. The reconstructed states then serve as input to an approximate Riemann solver that determines the numerical fluxes at two Gaussian quadrature points along the cell boundary. The use of multiple quadrature points renders the resulting flux high-order. Three types of approximate Riemann solvers are compared, including the widely used solver of Rusanov, the solver of Roe and the new AUSM⁺-up solver of Liou that has been designed for low-Mach number flows. Spatial discretizations are paired with either a third-order or fourth-order total-variation-diminishing Runge–Kutta timestepping scheme to match the order of the spatial discretization. The numerical schemes are evaluated with several standard shallow-water test cases that emphasize accuracy and conservation properties. These tests show that the AUSM⁺-up flux provides the best overall accuracy, followed closely by the Roe solver. The Rusanov flux, with its simplicity, provides significantly larger errors by comparison. A brief discussion on extending the method to arbitrary order-of-accuracy is included.     

High-order,finite-volume methods in mapped coordinates

We present an approach for constructing finite-volume methods for flux-divergence forms to any order of accuracy defined as the image of a smooth mapping from a rectangular discretization of an abstract coordinate space. Our approach is based on two ideas. The first is that of using higher-order quadrature rules to compute the flux averages over faces that generalize a method developed for Cartesian grids to the case of mapped grids. The second is a method for computing the averages of the metric terms on faces such that freestream preservation is automatically satisfied. We derive detailed formulas for the cases of fourth-order accurate discretizations of linear elliptic and hyperbolic partial differential equations. For the latter case, we combine the method so derived with Runge–Kutta time discretization and demonstrate how to incorporate a high-order accurate limiter with the goal of obtaining a method that is robust in the presence of discontinuities and underresolved gradients. For both elliptic and hyperbolic problems, we demonstrate that the resulting methods are fourth-order accurate for smooth solutions.     

High-order time-splitting Hermite and Fourier spectral methods

In this paper, we are concerned with the numerical solution of the time-dependent Gross–Pitaevskii Equation (GPE) involving a quasi-harmonic potential. Primarily, we consider discretisations that are based on spectral methods in space and higher-order exponential operator splitting methods in time. The resulting methods are favourable in view of accuracy and efficiency; moreover, geometric properties of the equation such as particle number and energy conservation are well captured.Regarding the spatial discretisation of the GPE, we consider two approaches. In the unbounded domain, we employ a spectral decomposition of the solution into Hermite basis functions; on the other hand, restricting the equation to a sufficiently large bounded domain, Fourier techniques are applicable. For the time integration of the GPE, we study various exponential operator splitting methods of convergence orders two, four, and six.Our main objective is to provide accuracy and efficiency comparisons of exponential operator splitting Fourier and Hermite pseudospectral methods for the time evolution of the GPE. Furthermore, we illustrate the effectiveness of higher-order time-splitting methods compared to standard integrators in a long-term integration.     

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.     

High-order commutator-free exponential time-propagation of driven quantum systems

We discuss the numerical solution of the Schrödinger equation with a time-dependent Hamilton operator using commutator-free time-propagators. These propagators are constructed as products of exponentials of simple weighted sums of the Hamilton operator. Owing to their exponential form they strictly preserve the unitarity of time-propagation. The absence of commutators or other computationally involved operations allows for straightforward implementation and application also to large scale and sparse matrix problems. We explain the derivation of commutator-free exponential time-propagators in the context of the Magnus expansion, and provide optimized propagators up to order eight. An extensive theoretical error analysis is presented together with practical efficiency tests for different problems. Issues of practical implementation, in particular the use of the Krylov technique for the calculation of exponentials, are discussed. We demonstrate for two advanced examples, the hydrogen atom in an electric field and pumped systems of multiple interacting two-level systems or spins that this approach enables fast and accurate computations.     

High-order compact exponential finite difference methods for convection–diffusion type problems

A class of high-order compact (HOC) exponential finite difference (FD) methods is proposed for solving one- and two-dimensional steady-state convection–diffusion problems. The newly proposed HOC exponential FD schemes have nonoscillation property and yield high accuracy approximation solution as well as are suitable for convection-dominated problems. The O(h⁴) compact exponential FD schemes developed for the one-dimensional (1D) problems produce diagonally dominant tri-diagonal system of equations which can be solved by applying the tridiagonal Thomas algorithm. For the two-dimensional (2D) problems, O(h⁴ + k⁴) compact exponential FD schemes are formulated on the nine-point 2D stencil and the line iterative approach with alternating direction implicit (ADI) procedure enables us to deal with diagonally dominant tridiagonal matrix equations which can be solved by application of the one-dimensional tridiagonal Thomas algorithm with a considerable saving in computing time. To validate the present HOC exponential FD methods, three linear and nonlinear problems, mostly with boundary or internal layers where sharp gradients may appear due to high Peclet or Reynolds numbers, are numerically solved. Comparisons are made between analytical solutions and numerical results for the currently proposed HOC exponential FD methods and some previously published HOC methods. The present HOC exponential FD methods produce excellent results for all test problems. It is shown that, besides including the excellent performances in computational accuracy, efficiency and stability, the present method has the advantage of better scale resolution. The method developed in this article is easy to implement and has been applied to obtain the numerical solutions of the lid driven cavity flow problem governed by the 2D incompressible Navier–Stokes equations using the stream function-vorticity formulation.     

自适应光学系统的实时模式复原算法

分析了自适应光学系统中实时模式复原算法的基本原理,建立了一种新型的传感器本征模复原算法。与常用的直接斜率法相比,这种模式复原算法可以有效减小探测噪声对复原计算过程的影响,提高系统的闭环稳定性和校正效果。在61单元自适应光学系统上实现了这种模式复原算法,并在实际大气湍流中对传感器本征模复原算法和直接斜率法进行了实验对比研究。     

自适应光学系统的实时模式复原算法

 分析了自适应光学系统中实时模式复原算法的基本原理，建立了一种新型的传感器本征模复原算法。与常用的直接斜率法相比，这种模式复原算法可以有效减小探测噪声对复原计算过程的影响，提高系统的闭环稳定性和校正效果。在61单元自适应光学系统上实现了这种模式复原算法，并在实际大气湍流中对传感器本征模复原算法和直接斜率法进行了实验对比研究。     

用于开环液晶自适应光学系统的模式预测技术研究

时间延迟误差是液晶自适应光学系统的一个最主要的误差源. 本文提出了一种利用智能模式预测-----迭代最小二乘(RLS)模式预测算法来克服其对成像分辨率的影响. 首先, 介绍了具有RLS模式预测能力的开环液晶自适应光学系统的结构和工作原理. 其次, 详细讨论了RLS模式预测算法的实现过程. 再次, 设计和搭建了一套带有液晶湍流模拟器的开环液晶自适应光学系统, 对RLS模式预测算法的预测效果进行了分析, 并和直接开环校正做了比较. 分析结果表明: 当系统处于中等强度湍流条件(大气相干长度r₀=6 cm, Greenwood频率f_G=35 Hz)和只有时间延迟误差情况下, 经过RLS预测后, 残差波面的RMS值由直接校正的0.26波长(1波长=785 nm)降低到了0.15波长, 校正效果提高了42%. 最后, 对预测前后自适应光学系统的成像效果进行了对比试验. 实验结果显示, 经过预测以后, 系统的成像分辨率由直接开环校正的25.4 cycles/mm提高到了32.0 cycles/mm, 成像分辨率提高了26%, 达到了0.9倍的衍射极限分辨率. 因此, RLS模式预测技术可以有效的提高开环液晶自适应系统的成像分辨率.     

具连续分布时滞的抛物型系统的指数镇定

针对一类具分布时滞的抛物型控制系统，提出一种新的镇定设计方法.利用推广的向量Hanalay 微分不等式、Dini导数、结合Green公式及不等式分析技术，在线性反馈控制律的作用下，导出了具分布时滞的抛物型控制系统的镇定的新判据.该镇定性条件不依赖于时滞.最后给了一个算例说明所得结果的可行性.此外，该方法一个明显的优点是所得的镇定性条件容易验证，因而便于应用. 关键词： 抛物型系统 Dini导数 分布时滞 全局指数镇定     

High-order shock-fitting methods for direct numerical simulation of hypersonic flow with chemical and thermal nonequilibrium

In recent years, much progress has been made in the direct numerical simulation of laminar-turbulent transition of hypersonic boundary layer flow. However, most of the efforts at the direct numerical simulation of transition previously have been focused on the idealized perfect gas flow or “cold” hypersonic flows. For practical problems in hypersonic flows, high-temperature effects of thermal and chemical nonequilibrium are important and cannot be modeled by a perfect gas model. Therefore, it is necessary to include the real gas models in the numerical simulation of hypersonic boundary layer transition in order to accurately predict flow field parameters. Currently most numerical methods for hypersonic flow with thermo-chemical nonequilibrium are based on shock-capturing approach at relatively low order of accuracy. Shock capturing schemes reduce to first-order accuracy near the shock and have been shown to produce spurious oscillations behind curved strong shocks. There is a need to develop new methods capable of simulating nonequilibrium hypersonic flow fields with uniformly high-order accuracy and avoid spurious oscillations near the shock. This paper presents a fifth-order shock-fitting method for numerical simulation of thermal and chemical nonequilibrium in hypersonic flows. The method is developed based on the state-of-the-art real gas models for thermo-chemical nonequilibrium and transport phenomena. Shock-fitting approach is used because it has the advantage of capturing the entire flow field with high-order accuracy and without any oscillations near the shock. The new method has been tested and validated for a number of test cases over a wide span of free stream conditions. The developed method is applied for the study of receptivity of free stream acoustic waves over a blunt cone for hypervelocity flow. Some preliminary results of the computations of the high order shock fitting method for the above mentioned study have also been presented.