A new pinning control scheme of complex networks based on data flow

In this paper, a new pinning control scheme called DF (data flow)-based pinning scheme is proposed. The new scheme can obtain the similar pinning efficiency with BC-based pinning scheme in real-world networks. Comparing with BC-based pinning scheme, DF-based pinning scheme has two main advantages. First, it just needs local information of network. Second, the new pinning scheme has a much lower time complexity than BC-based pinning scheme. In this paper, we have pinned two real-world networks (the US airline routing map network and the protein–protein network in yeast) to compare the new pinning scheme with degree-based, BC-based, LBC-based pinning schemes and we also pin a small-world network, a scale-free network to analyze DF-based pinning scheme in detail. Based on the Lyapunov stability theory, the validity of the scheme is proved. Finally, the numerical simulations are verified the effectiveness of the proposed method.     

Comparison of three solvers for viscoelastic fluid flow problems

The computational efficiency of three numerical schemes has been examined for the solution of a linearized system of equations resulting from the finite element discretization of a viscoelastic fluid flow problem. The first scheme is a modified frontal solver, which solves the linear system of equations directly. The other two, one based on a biconjugate gradient stabilized (BiCGStab) method and another based on a generalized minimal residual (GMRES) method, are iterative schemes. The stick-slip problem and the four-to-one contraction problem were analyzed and the viscoelastic fluid was assumed to obey the Oldroyd-B model. The two iterative schemes are superior to the direct scheme in terms of CPU time consumed and the BiCGStab scheme is even faster than the GMRES scheme. The range of convergence for both iterative schemes is compatible with that of the direct scheme.     

一种确定非线性振子多个共存周期解的数值方法

提出一种简洁高效数值方法，用以确定周期激励下非线性振子的多个共存周期响应。这是一种基于Ｐｏｉｎｃａｒé映射和二维自治系统奇点分类的几何方法，可在Ｐｏｉｎｃａｒé截面上得出多个周期运动的分布情况和类型。文中二个算例证实了该方法的有效性。     

随机结构静力反应概率密度演化方程的差分方法

随机结构分析的概率密度演化方法是分析随机结构静力反应的一种具有良好前景的方法。本文研究了求解随机结构静力反应概率密度演化方程的差分方法，分别探讨了单边差分格式和Lax-Wendroff格式的计算性态。二者均能满足概率相容性条件并且能够保证均值线性增长。以八层框架结构的静力随机反应为例，对两种差分格式的结果及精确解答进行了具体的比较分析。研究表明，两种差分格式均是收敛和稳定的，在不连续点处存在角点效应．单边差分格式能够保证概率非负性，而Lax-Wendroff格式具有往往更快的收敛速度。就变异系数而言，通常单边差分格式的变异系数随着区间离散数的增长而趋于稳定值，Lax-Wendroff格式则一开始就可得到恒定的值。     

A gas-kinetic theory based multidimensional high-order method for the compressible Navier–Stokes solutions

This paper presents a gas-kinetic theory based multidimensional high-order method for the compressible Naiver–Stokes solutions. In our previous study, a spatially and temporally dependent third-order flux scheme with the use of a third-order gas distribution function is employed.However, the third-order flux scheme is quite complicated and less robust than the second-order scheme. In order to reduce its complexity and improve its robustness, the secondorder flux scheme is adopted instead in this paper, while the temporal order of method is maintained by using a two stage temporal discretization. In addition, its CPU cost is relatively lower than the previous scheme. Several test cases in two and three dimensions, containing high Mach number compressible flows and low speed high Reynolds number laminar flows, are presented to demonstrate the method capacity.     

求解多维双曲守恒律方程组的四阶半离散格式

提出了求解多维双曲守恒律方程组的四阶半离散格式。该方法以中心加权基本无振荡(CWENO)重构为基础,同时考虑到在R iemann扇内波传播的局部速度,从而回避了计算过程中的网格交错,建立了数值耗散较小的介于迎风格式和中心格式之间的半离散格式。本文的四阶半离散格式是Kurganov等人的三阶半离散格式的高阶推广。大量的数值算例充分说明了本文方法的高分辨率和稳定性。     

A new high-accuracy scheme for compressible turbulent flows

Shock-capturing and broad-bandwidth scale resolutions are two main challenges of compressible turbulent flow simulation. To meet the rigorous requests, a novel fifth-order hybrid scheme based on a uniform hybrid framework is designed. With the help of a continuous weight operator, the new scheme combines an upwind compact scheme for smooth regions and a compact-reconstruction weighted essentially non-oscillatory scheme for discontinuous regions. Numerical analyses and canonical numerical tests confirm that the new scheme has high accuracy, spectral-like resolution property and shock-capturing capability. Besides, the new scheme shows high computational efficiency compared to the related shock-capturing schemes and hybrid ones.     

Particle trajectory calculations with a two‐step three‐time level semi‐Lagrangian scheme well suited for curved flows

This study proposes a new two‐step three‐time level semi‐Lagrangian scheme for calculation of particle trajectories. The scheme is intended to yield accurate determination of the particle departure position, particularly in the presence of significant flow curvature. Experiments were performed both for linear and non‐linear idealized advection problems, with different flow curvatures. Results for simulations with the proposed scheme, and with three other semi‐Lagrangian schemes, and with an Eulerian method are presented. In the linear advection problem the two‐step three‐time level scheme produced smaller root mean square errors and more accurate replication of the angular displacement of a Gaussian hill than the other schemes. In the non‐linear advection experiments the proposed scheme produced, in general, equal or better conservation of domain‐averaged quantities than the other semi‐Lagrangian schemes, especially at large Courant numbers. In idealized frontogenesis simulations the scheme performed equally or better than the other schemes in the representation of sharp gradients in a scalar field. The two‐step three‐time level scheme has some computational overhead as compared with the other three semi‐Lagrangian schemes. Nevertheless, the additional computational effort was shown to be worthwhile, due to the accuracy obtained by the scheme in the experiments with large time steps. The most remarkable feature of the scheme is its robustness, since it performs well both for small and large Courant numbers, in the presence of weak as well strong flow curvatures. Copyright © 2009 John Wiley & Sons, Ltd.     

A Non Homogeneous Riemann Solver for shallow water and two phase flows

In this work we consider a two steps finite volume scheme, recently developed to solve nonhomogeneous systems. The first step of the scheme depends on a diffusion control parameter which we modulate, using the limiters theory. Results on Shallow water equations and two phase flows are presented.     

关于对流差分格式有界性与CBC准则的讨论

通常认为CBC准则是差分格式有界性的充分条件。本文采用满足与不满足CBC准则的两种高阶差分格式对非线性问题进行了求解,重新讨论了格式有界性与CBC准则的关系,得出结论如下：在数值方法稳定的前提下,CBC准则下的有界模式是求解有界的充分条件,而非必要条件;此外,文章还分析了张涵信三阶精度格式的特点。     

Applications of a higher‐order bounded numerical scheme to turbulent flows

This paper applies the higher‐order bounded numerical scheme Weighted Average Coefficients Ensuring Boundedness (WACEB) to simulate two‐ and three‐dimensional turbulent flows. In the scheme, a weighted average formulation is used for interpolating the variables at cell faces and the weighted average coefficients are determined from a normalized variable formulation and total variation diminishing (TVD) constraints to ensure the boundedness of the solution. The scheme is applied to two turbulent flow problems: (1) two‐dimensional turbulent flow around a blunt plate; and (2) three‐dimensional turbulent flow inside a mildly curved U‐bend. In the present study, turbulence is evaluated by using a low‐Reynolds number version of the k–ω model. For the flow simulation, the QUICK scheme is applied to the momentum equations while either the WACEB scheme (Method 1) or the UPWIND scheme (Method 2) is used for the turbulence equations. The present study shows that the WACEB scheme has at least second‐order accuracy while ensuring boundedness of the solutions. The present numerical study for a pure convection problem shows that the ‘TVD’ slope ranges from 2 to 4. For the turbulent recirculating flow, two different mixed procedures (Method 1 and Method 2) produce a substantial difference for the mean velocities as well as for the turbulence kinetic energy. Method 1 predicts better results than Method 2 does, comparing the analytical solution and the experimental data. For the turbulent flow inside the mildly curved U‐bend, although the predictions of velocity distributions with two procedures are very close, a noticeable difference of turbulence kinetic energy is exhibited. It is noticed that the discrepancy exists between numerical results and the experimental data. The reason is the limit of the two‐equation turbulence model to such complex turbulent flows with extra strain‐rates. Copyright © 2001 John Wiley & Sons, Ltd.     

二维时-空守恒格式及其在激波传播计算中的应用

将改进后的一维时－空守恒格式推广到了二维情形,得到了一个新的一般形式的二维Ｅｕｌｅｒ方程时－空守恒格式,并用该格式对几个具有复杂波系的流场进行了数值模拟。结果表明,该格式保留了一维格式通用性好、结构简单的优点,其计算结果精度高,对激波等间断具有很强的分辨率。     

加权型紧致格式与加权本质无波动格式的比较

线性紧致格式和加权本质无波动格式是两种典型的高阶精度数值格式,它们各有优缺点.线性紧致格式在具有高阶精度的同时,格式的分辨率也比较高,耗散低,是计算多尺度流场结构的较好格式,但是不能计算具有强激波的流场.加权本质无波动格式是一种高阶精度捕捉激波格式,鲁棒性好,但耗散比较高,分辨率也不理想.近年来,在莱勒的线性紧致格式基础上,采用加权本质无波动格式捕捉激波思想,发展了一系列加权型紧致格式.本文较全面地比较了加权型紧致格式和加权本质无波动格式,包括构造方法、鲁棒性、分辨率、耗散特性、收敛特性以及并行计算效率.结果表明,现有的加权型紧致格式基本保持了加权本质无波动格式的性质,对于气动力等宏观量的计算,比加权本质无波动格式没有明显的优势.      

Selective lumping finite element method for shallow water flow

A finite element method for solving shallow water flow problems is presented. The standard Galerkin method is employed for spatial discretization. The numerical integration scheme for the time variation is the explicit two step scheme, which was originated by the authors and their co-workers. However, the original scheme has been improved to remove the erroneous artifical damping effect. Since the improved scheme employs a combination of lumped and unlumped coefficients, the scheme is referred to as a selective lumping scheme. Stability conditions and accuracy are investigated by considering several numerical examples. The method has been applied to the tidal flow in Osaka Bay and Yatsushiro Bay.     

非线性双曲型守恒律的两步二阶TVD差分格式

通过Mac Cormack格式和Warming-Beam的结合,构造了一种非常简单的两步二阶TVD差分格式,该差分格式更适合于使用分量形式差分计算而无须对欧拉方程组进行特征解耦。通过对流体力学方程组的大量数值试验,并与二阶ENO格式进行了比较,充分显示了该格式高精度、高分辨并且极其简单的优良特性。     

高阶多维半离散中心迎风格式及其应用

提出了求解多维对流-扩散方程的四阶半离散中心迎风格式。该格式以中心加权基本无振荡(CWENO)重构为基础,同时考虑到在Riemann扇内波传播的局部速度,从而更加准确地估计出了局部Riemann扇的宽度,最终既回避了网格的交错,又降低了格式的数值粘性,建立了介于迎风格式和中心格式之间的半离散中心迎风格式。本文还将该四阶半离散中心迎风格式与涡度-流函数方法相结合,有效地求解了二维不可压Euler方程组和Navier-Stokes方程组。     

一般坐标系下的二维Euler方程时—空守恒格式

将文「１－３」中的一维时－空守恒格式推广到了二维情形,得到了一般坐标系下的二维Ｅｕｌｅｒ方程时－空守恒格式,并用几个典型算例进行了检验计算,结果表明：本文得到的二维时－空守恒格式保留了一维格式所有的优点,格式简单,通用性强,而且对激波等间断具有很高的分辨率。     

A two‐dimensional numerical scheme of dry/wet fronts for the Saint‐Venant system of shallow water equations

We propose a new two‐dimensional numerical scheme to solve the Saint‐Venant system of shallow water equations in the presence of partially flooded cells. Our method is well balanced, positivity preserving, and handles dry states. The latter is ensured by using the draining time step technique in the time integration process, which guarantees non‐negative water depths. Unlike previous schemes, our technique does not generate high velocities at the dry/wet boundaries, which are responsible for small time step sizes and slow simulation runs. We prove that the new scheme preserves ‘lake at rest’ steady states and guarantees the positivity of the computed fluid depth in the partially flooded cells. We test the new scheme, along with another recent scheme from the literature, against the analytical solution for a parabolic basin and show the improved simulation performance of the new scheme for two real‐world scenarios. Copyright © 2014 John Wiley & Sons, Ltd.     

Development and assessment of a variable-order non-oscillatory scheme for convection term discretization

A new scheme for convection term discretization is developed, called VONOS (variable-order non-oscillatory scheme). The development of the scheme is based on the behaviour of well-known non-oscillatory schemes in the pure convection of a step profile test case. The new scheme is a combination of the QUICK and BSOU (bounded second-order upwind) schemes. These two schemes do not have the same formal order of accuracy and for that reason the formal order of accuracy of the new scheme is variable. The scheme is conservative, bounded and accurate. The performance of the new scheme was assessed in three test cases. The results showed that it is more accurate than currently used higher-order schemes, so it can be used in a general purpose algorithm in order to save computational time for the same level of accuracy. © 1998 John Wiley & Sons, Ltd.     

An efficient three-dimensional semi-implicit finite element scheme for simulation of free surface flows

An efficient semi-implicit finite element model is proposed for the simulation of three-dimensional flows in stratified seas. The body of water is divided into a number of layers and the two horizontal momentum equations for each layer of water are first integrated vertically. Nine-node Lagrangian quadratic isoparametric elements are employed for spatial discretization in the horizontal domain. The time derivatives are approximated using a second-order-accurate semi-implicit time-stepping scheme. The distinguishing feature of the proposed numerical scheme is that only nodal values on the same vertical line are coupled. Two test cases for which analytic solutions are available are employed to test the proposed scheme. The test results show that the scheme is efficient and stable. A numerical experiment is also included to compare the proposed scheme with a finite difference scheme.