不可压缩Navier-Stokes方程数值模拟

建立不可压缩Navier-Stokes方程的Crank-Nicolson有限差分方法,数值模拟了在初始正弦波下的二维水槽内流体受到倾斜激励时流场和涡的演变机制.数值结果表明,初始波为正弦波时,流场内出现一个单涡,单涡下沉变成了两个小涡,两个小涡消失后流场内部出现三条规则的流场带,最后这三条流场带演变成一个尖涡,尖涡在周围流体的作用下演变成一个单涡,最后单涡在自由面消失,当耗散系数和Reynolds数增大时,流场和涡演化的周期变小.     

用DES分离涡方法数值模拟串列双圆柱绕流问题

主要开发了SST-DES和SST-DDES两种分离涡方法,并集成到基于开源代码平台Open-FOAM开发的CFD求解器naoe-FOAM-SJTU中.选用高Reynolds(雷诺)数下串列双圆柱绕流问题作为标准算例来验证所开发的分离涡方法.该标准算例此前在美国国家航空航天局兰利研究中心的两个不同风洞做过物理试验.该研究将数值模拟得到的时均流场信息和一些其他物理量同物理试验结果比较,同时讨论分析了三维瞬态流场结构.结果表明该文开发的SST-DES和SST-DDES分离涡方法能够解决高Reynolds数下有大量流动分离的复杂流动问题.     

用DES分离涡方法数值模拟串列双圆柱绕流问题（英文）

主要开发了SST-DES和SST-DDES两种分离涡方法,并集成到基于开源代码平台Open-FOAM开发的CFD求解器naoe-FOAM-SJTU中.选用高Reynolds(雷诺)数下串列双圆柱绕流问题作为标准算例来验证所开发的分离涡方法.该标准算例此前在美国国家航空航天局兰利研究中心的两个不同风洞做过物理试验.该研究将数值模拟得到的时均流场信息和一些其他物理量同物理试验结果比较,同时讨论分析了三维瞬态流场结构.结果表明该文开发的SST-DES和SST-DDES分离涡方法能够解决高Reynolds数下有大量流动分离的复杂流动问题.     

黏性分层流中椭球体俯仰振荡研究

以黏性密度分层流下椭球体自由俯仰振荡衰减过程为研究内容,建立了密度连续分层流数值计算模型.通过对经典小球黏性绕流场的数值模拟和增阻系数的计算验证了数值模型的正确性.以初始45°攻角下的椭球体俯仰振荡过程为研究对象,采用基于Aitken亚松弛适应算法的双向流固耦合方法,数值模拟了不同内Froude(弗汝德)数Fri下椭球体俯仰衰减振荡的动态绕流场.数值研究结果表明,俯仰振荡将上下搅动周围流体,在椭球体上下两侧对称形成四个密度涡环,密度的垂向分层限制了涡环的垂向传播,也加速了涡环的消失,这种限制助长了水平运动的发展,远场尾涡流场将以水平波动的形式传播.在较高的内Froude数Fri和Reynolds(雷诺)数Re下,双向耦合抑制了数值震荡.研究还发现,随着来流速度的增加,阻力系数不增反降,这说明,对于自由俯仰振荡的椭球体,负阻尼现象仍然会出现.     

超音速后掠椭圆柱横流定常涡的不稳定特性及转捩预测

实验发现,一定条件下超声速后掠圆柱的前半圆上可能发生定常涡引起的转捩.为了模拟高空超音速飞行器后掠翼前缘,以无限展向长度后掠椭圆柱为模型,基于eN方法,根据积分的N值,分析了迎风轴长度、Reynolds(雷诺)数、后掠角和Mach(马赫)数等参数的变化对超音速后掠椭圆柱横流定常涡不稳定性的影响.研究结果表明,后掠椭圆柱的迎风轴长度增加,会使横流定常涡的不稳定性增强;Reynolds数增大,使横流定常涡模态的不稳定性增强;不稳定性的强度与Reynolds数的大小为近似线性关系;飞行高度增加,Mach数变大,使横流定常涡的不稳定性变弱;后掠角在一定范围内的变化对横流定常涡的不稳定性影响不大.这些结果有助于提高对高空超音速飞行器翼前缘转捩机理的认识,为横流转捩预测提供理论指导     

电磁力控制下圆柱绕流的涡度拟能

Okubo-Weiss函数与流元的容变、畸变以及涡量相关,可以用来评估流场的涡结构．经该文数学证明,对于边界无滑移的,低Reynolds数的二维不可压缩流动,Okubo-Weiss函数的全流场积分为0．还以电磁控制的圆柱绕流为例,通过数值计算,对该结论进行了验证．根据计算结果,依据Okubo-Weiss函数值,对流场进行了划分,讨论了总涡度拟能、总变形率和Okubo-Weiss函数在流场中的分布规律,以及电磁力对分布的影响．     

槽道湍流中拟涡能输运对反向控制的瞬时响应

利用直接数值模拟研究了槽道湍流中脉动拟涡能输运对反向控制的瞬时响应． 发现流向和展向拟涡能的衰减首先由拉伸产生项的抑制引起,而法向拟涡能的减小是因为控制阻碍了平均剪切的倾斜．在控制的初始阶段,流向拟涡能的演化远远落后于其它两个分量的变化．法向涡量快速单调减小,并对其它两个分量的减弱起到了重要作用．     

转捩边界层中次生流向涡演化的数值研究

采用高精度直接数值模拟方法和高效的特征无反射边界条件,进行可压缩流转捩边界层中出现的次生流向涡演化的数值研究.精细的数值模拟结果清楚地揭示了转捩边界层的复杂流场中次生流向涡的形成和演化过程,探讨它对转捩至关重要的环状涡生成的影响,发现次生流向涡和主流向涡的共同作用形成环状涡的一种新机理.     

T型管中冷热流体混合过程热波动大涡模拟

在FLUENT平台上,运用大涡模拟湍流模型及Smagorinsky-Lilly亚格子模型,对T型混合结构中的温度波动进行了模拟.获得了分别用于描述时均温度和温度波动强度的无量纲时均温度和无量纲均方根温度.分析了用于热疲劳分析和生命周期评价关键参数,即温度波动的功率谱密度函数.对不同条件下的混合过程进行了模拟,分析不同Reynolds数Re和Richardson数Ri对混合过程和热波动的影响.     

基于有限体积法的非结构网格大涡模拟离散方法研究

非结构网格下的大涡模拟是解决复杂几何体高Reynolds(雷诺)数流动的有效途径.首先,基于有限体积法,研究了对流项和扩散项非结构网格下的离散方法.研究结果表明:基于TVD(total variation diminishing)限制器的限制中心差分格式保证了对流项的二阶精度并抑制了非物理振荡,同时,线性迎风格式虽然稳定,但数值耗散过大,且不能保证有界,中心差分格式引起了周期性非物理振荡;扩散项的超松弛非正交修正减小了网格非正交带来的离散误差,但修正系数须根据网格非正交的程度进行合理选取.为验证所述离散方法对大涡模拟的适用性,数值计算了Re=1.14×10~6下的非定常三维小球绕流,计算方法包括:计算网格用基于Delaunay三角剖分和Netgen前沿推进算法的四面体非结构网格;湍流模型用改进的延迟分离涡大涡模型;在离散格式的选取上,对流项用限制中心差分,扩散项加入非正交修正,插值格式用最小二乘法,时间项用二阶后向差分.计算结果表明,所用离散方法稳定收敛并且与实验数据基本吻合.     

Qualitative Analysis of Bifurcating Solutions in the Lengyel-Epstein Model

One of the most fundamental problems in theoretical biology is to explain the mechanisms by which patterns and forms are created in the＇living world. In his seminal paper ＂The Chemical Basis of Morphogenesis＂, Turing showed that a system of coupled reaction-diffusion equations can be used to describe patterns and forms in biological systems. However, the first experimental evidence to the Turing patterns was observed by De Kepper and her associates（1990） on the CIMA reaction in an open unstirred reactor, almost 40 years after Turing＇s prediction. Lengyel and Epstein characterized this famous experiment using a system of reaction-diffusion equations. The Lengyel-Epstein model is in the form as follows     

(0,1 ;0)-INTERPOLATION ON INFINITE INTERVAL (-∞, +∞)

In this paper, we study the explicit representation and convergence of (0, 1; 0)-interpolation on infinite interval, which means to determine a polynomial of degree ≤ 3n - 2 when the function values are prescribed at two set of points namely the zeros of Hn(x) and H′n(x) and the first derivatives at the zeros of H′n(x).     

RECENT PROGRESS ON SPHERICAL HARMONIC APPROXIMATION MADE BY BNU RESEARCH GROUP ——In memory of Professor Sun Yongsheng

As early as in 1990, Professor Sun Yongsheng, suggested his students at Beijing Normal University to consider research problems on the unit sphere. Under his guidance and encouragement his students started the research on spherical harmonic analysis and approximation. In this paper, we incompletely introduce the main achievements in this area obtained by our group and relative researchers during recent 5 years （2001-2005）. The main topics are： convergence of Cesaro summability, a.e. and strong summability of Fourier-Laplace series; smoothness and K-functionals; Kolmogorov and linear widths.     

On a class of self-similar processes with stationary increments in higher order Wiener chaoses

We study a class of self-similar processes with stationary increments belonging to higher order Wiener chaoses which are similar to Hermite processes. We obtain an almost sure wavelet-like expansion of these processes. This allows us to compute the pointwise and local Hölder regularity of sample paths and to analyse their behaviour at infinity. We also provide some results on the Hausdorff dimension of the range and graphs of multidimensional anisotropic self-similar processes with stationary increments defined by multiple Wiener–Itô integrals.     

A REPRESENTATION FORMULA RELATED TO SCHR(O)DINGER OPERATORS

Schr(o)dinger operator is a central subject in the mathematical study of quantum mechanics.Consider the Schrodinger operator H = -△ V on R, where △ = d2/dx2 and the potential function V is real valued. In Fourier analysis, it is well-known that a square integrable function admits an expansion with exponentials as eigenfunctions of -△. A natural conjecture is that an L2 function admits a similar expansion in terms of "eigenfunctions" of H, a perturbation of the Laplacian (see [7], Ch. Ⅺ and the notes), under certain condition on V.     

On a class of Riemann surfaces

It is considered the class of Riemann surfaces with dimT1 = 0, where T1 is a subclass of exact harmonic forms which is one of the factors in the orthogonal decomposition of the spaceΩH of harmonic forms of the surface, namely The surfaces in the class OHD and the class of planar surfaces satisfy dimT1 = 0. A.Pfluger posed the question whether there might exist other surfaces outside those two classes. Here it is shown that in the case of finite genus g, we should look for a surface S with dimT1 = 0 among the surfaces of the form Sg\K , where Sg is a closed surface of genus g and K a compact set of positive harmonic measure with perfect components and very irregular boundary.     

Journal of the Operations Research Society of China Special issue on Operations Research in Healthcare Management

正Guest Editors:Hong Chen,Shanghai Jiao Tong University,Shanghai,China Guohua Wan,Shanghai Jiao Tong University,Shanghai,China David Yao,Columbia University,New York,USA Scope:Healthcare delivery worldwide has been fraught with high cost,low efficiency and poor quality of patient care service.For the field of operations research(OR),healthcare offers some of the biggest challenges as well as best opportunities in     

ANALYSIS IN THEORY AND APPLICATIONS Vol. 30, 2014 CONTENTS

正~~     

An alternative mathematical algorithm for the photo- and videokeratoscope

Due to the resolution of current laser technology, the accuracy of corneal topography as measured by the videokeratoscope is no longer adequate to provide precise enough data for refractive surgery or for the fitting of customized contact lenses. We present an algorithm for recovering corneal topography that makes use of modern differential geometric techniques and numerical descent in Sobolev spaces. We believe this algorithm may be used with the photo- and videokeratoscope to increase the accuracy of the recovered corneal topography.