耦合广义非线性薛定谔方程的相互作用表象龙格库塔算法及其误差分析

本文通过表象变换, 将耦合广义非线性薛定谔方程 (C-GNLSE) 变换成相互作用表象中的向量方程, 再利用向量形式的4阶龙格-库塔迭代格式, 建立了一种在频域内求解C-GNLSE的同步更新迭代算法. 通过将该向量形式的相互作用表象中的4阶龙格-库塔 (V-JH-RK4IP) 算法应用于高双折射光子晶体光纤中超连续谱产生的数值模拟, 验证了算法的有效性, 通过与现有其他典型算法的比较, 表明以V-JH-RK4IP算法求解C-GNLSE具有最高的计算精度和计算效率. 关键词： 耦合广义非线性薛定谔方程(C-GNLSE) 相互作用表象 4阶龙格-库塔算法 超连续谱产生     

一维双曲守恒律的龙格-库塔控制体积间断有限元方法

给出数值求解一维双曲守恒律方程的新方法——龙格-库塔控制体积间断有限元方法(RKCVDFEM),其中空间离散基于控制体积有限元方法,时间离散基于二阶TVB Runge-Kutta技术,有限元空间选取为分段线性函数空间.理论分析表明,格式具有总变差有界(TVB)的性质,而且空间和时间离散形式上具有二阶精度.数值算例表明,数值解收敛到熵解并且对光滑解的收敛阶是最优的,优于龙格-库塔间断Galerkin方法(RKDGM)的计算结果.     

啁啾超短激光脉冲对稠密共振介质特性的调控

采用光与物质相互作用的半经典理论,建立了稠密二能级体系中含啁啾相位项的光学Bloch方程,并用高准确度、快速和可靠的四阶龙格-库塔法数值求解了该方程.通过数值计算,研究了啁啾超短激光脉冲与稠密共振二能级体系作用下的特性,得到了局域场和啁啾参量对Bloch矢量影响的规律.研究结果表明:局域场修正系数和啁啾参量对Bloch矢量的瞬态相干过程和稳态特性都会产生明显的调制作用,可以通过调节线性啁啾参量和局域场修正系数实现稳定的粒子数布居反转的调控.     

三点强紧致加权高分辨率高阶格式及其应用

本文分别给出了采用有限体积计算时采用高阶格式遇到的网格单元体界面上三阶、四阶与五阶格式精度下的数值通量表达式,并且给出了确定权函数过程中所遇到的光滑因子表达式．文中首先对模型方程进行了格式分辨率方面的检验,然后将格式用于某涡轮级的实际流场计算．数值结果表明:所给出的高分辨率、高阶格式具有较高的激波分辨率并且具有较高的数值精度,在计算流场时,它可以用较少网格点去取代普通低阶精度格式下所采用的较密网格．     

耗散紧致格式的频谱特性研究与应用

采用Fourier分析方法,给出线性耗散紧致格式和基于m级Runge-Kutta时间积分方法的全离散格式的频谱特性,并应用五阶耗散紧致格式模拟典型高频波传播和超声速平面Couette流动的特征值问题及其稳定性边值问题,展示耗散紧致差分格式良好的频谱特性.     

等离子体电极普克尔盒电光开关单脉冲过程数值模拟

采用流体模型对等离子体电极普克尔盒(PEPC)电光开关单脉冲过程进行了数值模拟分析.模型包括带电粒子连续性方程、动量守恒方程、电子平均能量方程及空间电位泊松方程.分别采用隐式指数差分格式，超松弛迭代法(SOR)和经典四阶龙格-库塔法(R-K)对带电粒子连续性方程，泊松方程和电子平均能量方程进行数值求解.模拟分析了PEPC单脉冲过程中的带电粒子浓度、电子温度、空间电场、PEPC的放电电流、晶体两侧电压和开关效率的时间演化特性.模型得出了PEPC中气体放电等离子体的微观物理过程与PEPC宏观参量的关系，对设计 关键词： 等离子体电极普克尔盒 电光开关 数值模拟 气体放电     

基于多步高阶差分格式求解时域Maxwell方程

提出基于无穷维哈密尔顿系统及分裂算子理论的多步高阶差分格式,求解时域Maxwell方程.在时间方向上,针对Maxwell方程采用不同阶数的辛算法进行差分离散;在空间方向上,采用四阶差分格式进行差分离散.探讨多步高阶差分格式的稳定性及数值色散性,最后给出数值计算结果.结果表明,五级四阶格式为最有效的多步高阶差分格式,具有高精度、占用较少的计算机资源等优点,适用于长时间的数值模拟.     

混沌和超混沌系统中的奇怪吸引子及其分析

用四阶定步长龙格—库塔算法对几种混沌和超混沌系统进行数值求解，绘制了各种系统典型奇怪吸引子的相图，对奇怪吸引子的结构和特性进行了分析。     

啁啾超短激光脉冲对二能级体系特性的调控

研究了啁啾超短激光脉冲与二能级体系近共振作用下的特性.采用光与物质相互作用的半经典理论,建立了含啁啾相位项的修正光学Bloch方程,并用高准确度、快速和可靠的四阶龙格-库塔法数值求解了该方程.通过数值计算,得到啁啾符号和啁啾量与Bloch矢量的关系,以及共振失谐量符号和啁啾符号对Bloch矢量性质影响的规律.结果表明,啁啾的符号和大小对Bloch矢量的瞬态相干过程和稳态特性都会产生明显的调制作用,由此可以通过调节激光脉冲的啁啾特性来实现对二能级体系性质的调控.     

基于(火积)耗散理论的喷淋室内传热传质性能优化研究

基于(火积)理论,在气-水传热传质方程的基础上,建立了立式喷淋室内热学分析模型。利用四阶龙格-库塔方法,对该模型进行了数值求解,分析了不同参数对传热传质过程的影响。结果表明,(火积)耗散热阻用来优化喷淋室内传热传质效果是有效的,(火积)耗散热阻可以用来评价喷淋室参数变化对换热效果的影响,各参数对换热效果影响的显著性按下列顺序依次降低:空气进口温度、空气流量、空气相对湿度、水量、水滴直径、塔高、气流速度、水滴速度。     

Nonuniform time-step Runge–Kutta discontinuous Galerkin method for Computational Aeroacoustics

With many superior features, Runge–Kutta discontinuous Galerkin method (RKDG), which adopts Discontinuous Galerkin method (DG) for space discretization and Runge–Kutta method (RK) for time integration, has been an attractive alternative to the finite difference based high-order Computational Aeroacoustics (CAA) approaches. However, when it comes to complex physical problems, especially the ones involving irregular geometries, the time step size of an explicit RK scheme is limited by the smallest grid size in the computational domain, demanding a high computational cost for obtaining time accurate numerical solutions in CAA. For computational efficiency, high-order RK method with nonuniform time step sizes on nonuniform meshes is developed in this paper. In order to ensure correct communication of solutions on the interfaces of grids with different time step sizes, the values at intermediate-stages of the Runge–Kutta time integration on the elements neighboring such interfaces are coupled with minimal dissipation and dispersion errors. Based upon the general form of an explicit p-stage RK scheme, a linear coupling procedure is proposed, with details on the coefficient matrices and execution steps at common time-levels and intermediate time-levels. Applications of the coupling procedures to Runge–Kutta schemes frequently used in simulation of fluid flow and acoustics are given, including the third-order TVD scheme, and low-storage low dissipation and low dispersion (LDDRK) schemes. In addition, an analysis on the stability of coupling procedures on a nonuniform grid is carried out. For validation, numerical experiments on one-dimensional and two-dimensional problems are presented to illustrate the stability and accuracy of proposed nonuniform time-step RKDG scheme, as well as the computational benefits it brings. Application to a one-dimensional nonlinear problem is also investigated.     

An Eulerian–Lagrangian WENO finite volume scheme for advection problems

We develop a locally conservative Eulerian–Lagrangian finite volume scheme with the weighted essentially non-oscillatory property (EL–WENO) in one-space dimension. This method has the advantages of both WENO and Eulerian–Lagrangian schemes. It is formally high-order accurate in space (we present the fifth order version) and essentially non-oscillatory. Moreover, it is free of a CFL time step stability restriction and has small time truncation error. The scheme requires a new integral-based WENO reconstruction to handle trace-back integration. A Strang splitting algorithm is presented for higher-dimensional problems, using both the new integral-based and pointwise-based WENO reconstructions. We show formally that it maintains the fifth order accuracy. It is also locally mass conservative. Numerical results are provided to illustrate the performance of the scheme and verify its formal accuracy.     

A conservative numerical scheme for solving an autonomous Hamiltonian system

A new numerical scheme is proposed for solving Hamilton’s equations that possesses the properties of symplecticity. Just as in all symplectic schemes known to date, in this scheme the conservation laws of momentum and angular momentum are satisfied exactly. A property that distinguishes this scheme from known schemes is proved: in the new scheme, the energy conservation law is satisfied for a system of linear oscillators. The new numerical scheme is implicit and has the third order of accuracy with respect to the integration step. An algorithm is presented by which the accuracy of the scheme can be increased up to the fifth and higher orders. Exact and numerical solutions to the two-body problem, calculated by known schemes and by the scheme proposed here, are compared.     

Adaptive mesh refinement based on high order finite difference WENO scheme for multi-scale simulations

In this paper, we propose a finite difference AMR-WENO method for hyperbolic conservation laws. The proposed method combines the adaptive mesh refinement (AMR) framework  and  with the high order finite difference weighted essentially non-oscillatory (WENO) method in space and the total variation diminishing (TVD) Runge–Kutta (RK) method in time (WENO-RK)  and  by a high order coupling. Our goal is to realize mesh adaptivity in the AMR framework, while maintaining very high (higher than second) order accuracy of the WENO-RK method in the finite difference setting. The high order coupling of AMR and WENO-RK is accomplished by high order prolongation in both space (WENO interpolation) and time (Hermite interpolation) from coarse to fine grid solutions, and at ghost points. The resulting AMR-WENO method is accurate, robust and efficient, due to the mesh adaptivity and very high order spatial and temporal accuracy. We have experimented with both the third and the fifth order AMR-WENO schemes. We demonstrate the accuracy of the proposed scheme using smooth test problems, and their quality and efficiency using several 1D and 2D nonlinear hyperbolic problems with very challenging initial conditions. The AMR solutions are observed to perform as well as, and in some cases even better than, the corresponding uniform fine grid solutions. We conclude that there is significant improvement of the fifth order AMR-WENO over the third order one, not only in accuracy for smooth problems, but also in its ability in resolving complicated solution structures, due to the very low numerical diffusion of high order schemes. In our work, we found that it is difficult to design a robust AMR-WENO scheme that is both conservative and high order (higher than second order), due to the mass inconsistency of coarse and fine grid solutions at the initial stage in a finite difference scheme. Resolving these issues as well as conducting comprehensive evaluation of computational efficiency constitute our future work.     

A high-order accurate hybrid scheme using a central flux scheme and a WENO scheme for compressible flowfield analysis

A high-order accurate hybrid central-WENO scheme is proposed. The fifth order WENO scheme [G.S. Jiang, C.W. Shu, Efficient implementation of weighted ENO schemes, J. Comput. Phys. 126 (1996) 202–228] is divided into two parts, a central flux part and a numerical dissipation part, and is coupled with a central flux scheme. Two sub-schemes, the WENO scheme and the central flux scheme, are hybridized by means of a weighting function that indicates the local smoothness of the flowfields. The derived hybrid central-WENO scheme is written as a combination of the central flux scheme and the numerical dissipation of the fifth order WENO scheme, which is controlled adaptively by a weighting function. The structure of the proposed hybrid central-WENO scheme is similar to that of the YSD-type filter scheme [H.C. Yee, N.D. Sandham, M.J. Djomehri, Low-dissipative high-order shock-capturing methods using characteristic-based filters, J. Comput. Phys. 150 (1999) 199–238]. Therefore, the proposed hybrid scheme has also certain merits that the YSD-type filter scheme has. The accuracy and efficiency of the developed hybrid central-WENO scheme are investigated through numerical experiments on inviscid and viscous problems. Numerical results show that the proposed hybrid central-WENO scheme can resolve flow features extremely well.     

Polarization-free second-order nonlinear frequency conversion using the optical superlattice

We propose a universal and practical scheme to realize polarization-free second-order nonlinear frequency-conversion processes, including sum-frequency generation, difference-frequency generation, and optical parametric amplification. This scheme is based on the optical superlattice with a noncritical phase-matching condition, which is suitable for optical integration. A chirped dual-period structure is proposed to ensure the efficiency and stability of the conversions.     

Super-sensitivity measurement of tiny Doppler frequency shifts based on parametric amplification and squeezed vacuum state

The precision measurement of Doppler frequency shifts is of great significance for improving the precision of speed measurement. This paper proposes a precision measurement scheme of tiny Doppler shifts by a parametric amplification process and squeezed vacuum state. This scheme takes a parametric amplification process and squeezed vacuum state into a detection system, so that the measurement precision of tiny Doppler shifts can exceed the Cram′er–Rao bound of coherent light. Simultaneously, a simulation study is carried out on the theoretical basis, and the following results are obtained: for the signal light of Gaussian mode, when the amplification factor g = 1 and the squeezed factor r = 0.5, the measurement error of Doppler frequency shifts is 14.4% of the Cramer–Rao bound of the coherent light in our system. At the same time,when the local light mode and squeezed vacuum state mode are optimized, the measurement precision of this scheme can be further improved by ■ times, where n is the mode-order of the signal light.     

Large eddy simulation using a new set of sixth order schemes for compressible viscous terms

A set of conservative sixth order central differencing schemes is suggested for compressible flows with variable viscosity coefficient. This new set of central differencing schemes has the stencil width matching that of the seventh order weighted essentially non-oscillatory scheme (WENO). This feature is important to maintain the compactness of the seventh order WENO scheme and facilitate boundary condition treatment. As an application example, a large eddy simulation (LES) is conducted for a cylinder flow using the seventh order WENO scheme for the convective terms and the new set of sixth order central differencing scheme for the viscous terms. The results are compared with those from other research groups and those obtained using the fifth order WENO scheme and fourth order central differencing.     

一类新的差分格式优化目标函数

对差分格式进行优化处理可以提高其谱精度。与高精度(指Taylor展开精度)格式相比,优化格式放大因子的误差随波数的变化不是单调的,而是必然会出现极值点,这样就存在临界距离R_cr,在此距离内优化格式描述的数值波的积累误差小于高精度格式,而超出此距离后优化格式的误差反而大,对于非定常流及气动声学计算来说,控制差分格式的临界距离是必要的。一般的优化目标函数以每个时间推进步的误差为基础(即放大因子法),R_cr不能在优化过程中确定。对此进行分析,指出积累误差的重要性并提出以此为基础的新的优化目标函数,这样在对格式进行优化时可以直接指定临界距离,从而为控制谱精度提供方便。     

Sliding observer for synchronization of fractional order chaotic systems with mismatched parameter

In this paper, we propose an observer-based fractional order chaotic synchronization scheme. Our method concerns fractional order chaotic systems in Brunovsky canonical form. Using sliding mode theory, we achieve synchronization of fractional order response with fractional order drive system using a classical Lyapunov function, and also by fractional order differentiation and integration, i.e. differintegration formulas, state synchronization proved to be established in a finite time. To demonstrate the efficiency of the proposed scheme, fractional order version of a well-known chaotic system; Arnodo-Coullet system is considered as illustrative examples.