Optimal dynamical systems of Navier-Stokes equations based on generalized helical-wave bases and the fundamental elements of turbulence

In this paper, we present the theory of constructing optimal generalized helical-wave coupling dynamical systems. Applying the helical-wave decomposition method to Navier-Stokes equations, we derive a pair of coupling dynamical systems based on optimal generalized helical-wave bases. Then with the method of multi-scale global optimization based on coarse graining analysis, a set of global optimal generalized helical-wave bases is obtained. Optimal generalized helical-wave bases retain the good properties of classical helical-wave bases. Moreover, they are optimal for the dynamical systems of Navier-Stokes equations, and suitable for complex physical and geometric boundary conditions. Then we find that the optimal generalized helical-wave vortexes fitted by a finite number of optimal generalized helical-wave bases can be used as the fundamental elements of turbulence, and have important significance for studying physical properties of complex flows and turbulent vortex structures in a deeper level.     

Local RBFs Based Collocation Methods for Unsteady Navier-Stokes Equations

The local RBFs based collocation methods (LRBFCM) is presented to solve two-dimensional incompressible Navier-Stokes equations. In avoiding the ill-conditioned problem, the weight coefficients of linear combination with respect to the function values and its derivatives can be obtained by solving low-order linear systems within local supporting domain. Then, we reformulate local matrix in the global and sparse matrix. The obtained large sparse linear systems can be directly solved instead of using more complicated iterative method. The numerical experiments have shown that the developed LRBFCM is suitable for solving the incompressible Navier-Stokes equations with high accuracy and efficiency.     

Coherence resonance in coupled chaotic oscillators

Existing works on coherence resonance, i.e., the phenomenon of noise-enhanced temporal regularity, focus on excitable dynamical systems such as those described by the FitzHugh-Nagumo equations. We extend the scope of coherence resonance to an important class of dynamical systems: coupled chaotic oscillators. In particular, we show that, when a system of coupled chaotic oscillators is under the influence of noise, the degree of temporal regularity of dynamical variables characterizing the difference among the oscillators can increase and reach a maximum value at some optimal noise level. We present numerical results illustrating the phenomenon and give a physical theory to explain it.     

A Novel Low-Dimensional Method for Analytically Solving Partial Differential Equations

This paper is concerned with a low-dimensional dynamical system model for analytically solving partial differential equations (PDEs). The model proposed is based on a posterior optimal truncated weighted residue (POT-WR) method, by which an infinite dimensional PDE is optimally truncated and analytically solved in required condition of accuracy. To end that, a POT-WR condition for PDE under consideration is used as a dynamically optimal control criterion with the solving process. A set of bases needs to be constructed without any reference database in order to establish a space to describe low-dimensional dynamical system that is required. The Lagrangian multiplier is introduced to release the constraints due to the Galerkin projection, and a penalty function is also employed to remove the orthogonal constraints. According to the extreme principle, a set of ordinary differential equations is thus obtained by taking the variational operation of the generalized optimal function. A conjugate gradient algorithm by FORTRAN code is developed to solve the ordinary differential equations. The two examples of one-dimensional heat transfer equation and nonlinear Burgers’ equation show that the analytical results on the method proposed are good agreement with the numerical simulations and analytical solutions in references, and the dominant characteristics of the dynamics are well captured in case of few bases used only.     

Energy spectra of the stable boundary layer: A theoretical model

Summary A derivation of spectral energy functions for stable atmospheric boundary layers is presented for a local similarity regime. Previous studies of this nocturnal atmospheric conditions were usually carried through purely kinematic modellizations. In this study, the atmospheric turbulence energy content is derived from a dynamical viewpoint. This is achieved through the use of the energy balance relation obtained by Rotta from Navier-Stokes equations. The result is an appropriate physical rooting of the spectral energy curves known as Kaimal’s Isopleths. The authors of this paper have agreed to not receive the proofs for correction.     

Dark Equations and Their Light Integrability

A relatively new approach to analyzing integrability, based upon differential-algebraic and symplectic techniques, is applied to some “dark equations ”of the type introduced by Boris Kupershmidt. These dark equations have unusual properties and are not particularly well-understood. In particular, dark three-component polynomial Burgers type systems are studied in detail. Their matrix Lax representations are constructed, and the related symmetry recursion operators and infinite hierarchies of integrable nonlinear dynamical systems along with their Lax representations are derived. New linear Lax spectral problems for dark integrable countable hierarchies of dynamical systems are proposed and some special cases are considered as a means of indicating that the approach presented is applicable to a far wider class of dark equations than analyzed here.     

On the Global Regularity of a Helical-Decimated Version of the 3D Navier-Stokes Equations

We study the global regularity, for all time and all initial data in H ^1/2, of a recently introduced decimated version of the incompressible 3D Navier-Stokes (dNS) equations. The model is based on a projection of the dynamical evolution of Navier-Stokes (NS) equations into the subspace where helicity (the L ²-scalar product of velocity and vorticity) is sign-definite. The presence of a second (beside energy) sign-definite inviscid conserved quadratic quantity, which is equivalent to the H ^1/2-Sobolev norm, allows us to demonstrate global existence and uniqueness, of space-periodic solutions, together with continuity with respect to the initial conditions, for this decimated 3D model. This is achieved thanks to the establishment of two new estimates, for this 3D model, which show that the H ^1/2 and the time average of the square of the H ^3/2 norms of the velocity field remain finite. Such two additional bounds are known, in the spirit of the work of H. Fujita and T. Kato (Arch. Ration. Mech. Anal. 16:269–315, 1964; Rend. Semin. Mat. Univ. Padova 32:243–260, 1962), to be sufficient for showing well-posedness for the 3D NS equations. Furthermore, they are directly linked to the helicity evolution for the dNS model, and therefore with a clear physical meaning and consequences.     

带有结构非线性的跨音速翼型颤振特性研究

以非定常N-S方程为主管方程,采用时间推进的方法,计算翼型振荡的瞬态非定常气动力,并与带有结构非线性的颤振方程耦合求解,计算了带有结构刚度非线性(间隙型,三次型刚度非线性)和结构阻尼非线性(三次型阻尼非线性)的结构响应特性和颤振特性.计算研究表明,由于同时具有结构和气动非线性,振荡极限环和气动力极为复杂.     

Bifurcations of a lattice gas flow under external forcing

We study the behavior of a Frisch-Hasslacher-Pomeau lattice gas automaton under the effect of a spatially periodic forcing. It is shown that the lattice gas dynamics reproduces the steady-state features of the bifurcation pattern predicted by a properly truncated model of the Navier-Stokes equations. In addition, we show that the dynamical evolution of the instabilities driving the bifurcation can be modeled by supplementing the truncated Navier-Stokes equation with a random force chosen on the basis of the automaton noise.     

Instabilities in the Chapman-Enskog Expansion and Hyperbolic Burnett Equations

It is well-known that the classical Chapman-Enskog procedure does not work at the level of Burnett equations (the next step after the Navier-Stokes equations). Roughly speaking, the reason is that the solutions of higher equations of hydrodynamics (Burnett's, etc.) become unstable with respect to short-wave perturbations. This problem was recently attacked by several authors who proposed different ways to deal with it. We present in this paper one of possible alternatives. First we deduce a criterion for hyperbolicity of Burnett equations for the general molecular model and show that this criterion is not fulfilled in most typical cases. Then we discuss in more detail the problem of truncation of the Chapman-Enskog expansion and show that the way of truncation is not unique. The general idea of changes of coordinates (based on analogy with the theory of dynamical systems) leads finally to nonlinear Hyperbolic Burnett Equations (HBEs) without using any information beyond the classical Burnett equations. It is proved that HBEs satisfy the linearized H-theorem. The linear version of the problem is studied in more detail, the complete Chapman-Enskog expansion is given for the linear case. A simplified proof of the Slemrod identity for Burnett coefficients is also given.     

流动数值模拟中一种并行自适应有限元算法

给出了一种流动数值模拟中的基于误差估算的并行网格自适应有限元算法.首先,以初网格上获得的当地事后误差估算值为权,应用递归谱对剖分方法划分初网格,使各子域上总体误差近似相等,以解决负载平衡问题.然后以误差值为判据对各子域内网格进行独立的自适应处理.最后应用基于粘接元的区域分裂法在非匹配的网格上求解N-S方程.区域分裂情形下N-S方程有限元解的误差估算则是广义Stokes问题误差估算方法的推广.为验证方法的可靠性,给出了不可压流经典算例的数值结果.     

大雷诺数Navier-Stokes方程的两水平亚格子模型稳定化方法

基于两重网格离散方法,提出三种求解大雷诺数定常Navier-Stokes方程的两水平亚格子模型稳定化有限元算法.其基本思想是首先在一粗网格上求解带有亚格子模型稳定项的Navier-Stokes方程,然后在细网格上分别用三种不同的校正格式求解一个亚格子模型稳定化的线性问题,以校正粗网格解.通过适当的稳定化参数和粗细网格尺寸的选取,这些算法能取得最优渐近收敛阶的有限元解.最后,用数值模拟验证三种算法的有效性.     

基于数值原子轨道基组的第一性原理计算软件ABACUS

随着超级计算机硬件和数值算法迅速发展, 使得目前利用密度泛函理论研究上千个原子体系的电子能带和结构等性质变得可行. 数值原子轨道基组由于其基组较小和局域等特性, 可以很好地与电子结构计算中的线性标度算法等的新算法结合, 用来研究较大尺寸的物理体系. 本文详细介绍了一款中国科学技术大学量子信息重点 实验室自主开发的基于数值原子轨道基组的第一性原理计算软件 Atomic-orbital Based Ab-initio Computation at UStc. 大量的测试结果表明: 该软件具有很好的准确性和较高的并行效率, 可以用于包含1000个原子左右的系统的电子结构和原子结构的研究以及分子动力学模拟计算.     

A geometrical setting for Lax equations associated to dynamical systems

We develop a geometrical framework for dealing with Lax equations associated to dynamical systems over a manifold M. We also show that this theory reproduces the global versions of Lax equations given before as well as the usual theory of reduced systems obtained from systems defined on Lie groups and with such group as a symmetry group.     

Birkhoff系统的离散最优控制及其在航天器交会对接中的应用

由于非线性,最优控制问题通常依赖于数值求解,即通过离散目标泛函和受控运动方程转化为一有限维的非线性最优化问题.最优控制问题中的受控运动方程在表示为受控Birkhoff方程的形式之后,可以利用受控Birkhoff方程的离散变分差分格式进行离散.与按照传统差分格式近似受控运动方程相比,此途径可以诱导更加真实可靠的非线性最优化问题,进而也会诱导更加精确有效的离散最优控制.应用于航天器交会对接问题,该种数值求解最优控制问题的方法在较大时间步长的情况下仍然求得了一个有效实现交会对接的离散最优控制.模拟结果验证了该方法的有效性.     

A hyperchaotic system stabilization via inverse optimal control and experimental research

In this paper, some basic dynamical properties of a four-dimensional autonomous hyperchaotic system are investigated by means of Poincar′e mapping, Lyapunov exponents and bifurcation diagram. The dynamical behaviours of this new hyperchaotic system are proved not only by performing numerical simulation and brief theoretical analysis but also by conducting an electronic circuit experiment. An efficient approaching is developed for global asymptotic stabilization of this four-dimensional hyperchaotic system. Based on the method of inverse optimal control for nonlinear systems, a linear state feedback is electronically implemented. It is remarkably simple as compared with other chaos control ways, like nonlinear state feedback.     

Poisson Involutions, Spin Calogero-Moser Systems Associated with Symmetric Lie Subalgebras and the Symmetric Space Spin Ruijsenaars-Schneider Models

We develop a general scheme to construct integrable systems starting from realizations in symmetric coboundary dynamical Lie algebroids and symmetric coboundary dynamical Poisson groupoids. The method is based on the successive use of Dirac reduction and Poisson reduction. Then we show that certain spin Calogero-Moser systems associated with symmetric Lie subalgebras can be studied in this fashion. We also consider some spin-generalized Ruisjenaars-Schneider equations which correspond to the N-soliton solutions of affine Toda field theory. In this case, we show how the equations are obtained from the Dirac reduction of some Hamiltonian system on a symmetric coboundary dynamical Poisson groupoid.     

Extracting dynamical equations from experimental data is NP hard

The behavior of any physical system is governed by its underlying dynamical equations. Much of physics is concerned with discovering these dynamical equations and understanding their consequences. In this Letter, we show that, remarkably, identifying the underlying dynamical equation from any amount of experimental data, however precise, is a provably computationally hard problem (it is NP hard), both for classical and quantum mechanical systems. As a by-product of this work, we give complexity-theoretic answers to both the quantum and classical embedding problems, two long-standing open problems in mathematics (the classical problem, in particular, dating back over 70?years).     

Discrete data assimilation in the Lorenz and 2D Navier-Stokes equations

Consider a continuous dynamical system for which partial information about its current state is observed at a sequence of discrete times. Discrete data assimilation inserts these observational measurements of the reference dynamical system into an approximate solution by means of an impulsive forcing. In this way the approximating solution is coupled to the reference solution at a discrete sequence of points in time. This paper studies discrete data assimilation for the Lorenz equations and the incompressible two-dimensional Navier-Stokes equations. In both cases we obtain bounds on the time interval h between subsequent observations which guarantee the convergence of the approximating solution obtained by discrete data assimilation to the reference solution.