分析半透明单晶生长传热过程的边界元方法

论文发展了一个能求解带相变运动界面非定常传热和非线性热物理特性问题的双倒易边界元方法。数值模拟了半透明单晶生长中热过程的一个例子。由于方法是纯边界积分方法,计算量与计算内存都大大减少。获得了单晶生长过程瞬态温度场分布和固液相界面形状时间推进的一些结果。     

基于自动微分的线性/非线性谐波方法研究

本文运用线性/非线性谐波方法并结合自动微分开发求解一维非定常欧拉方程组的求解器.线性/非线性谐波方法求解非定常流不仅精度高,而且计算时间短,自动微分方法的应用又极大地减少了谐波方程求解器开发的工作量.在定常求解器的基础上,利用自动微分工具生成求解谐波方程的基本代码,通过手动组合和微调最终得到求解一维非定常欧拉方程的求解器.最后通过准一维拉瓦尔喷管算例和直管算例来验证求解器的准确性.     

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

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

最小二乘有限元法求解非定常应力的Navier-Stokes方程

为精确求解非定常层流问题,发展一种非定常速度-应力-压力的方法.采用牛顿法对非线性对流项进行线性化处理和预处理共轭梯度法,实现了非定常应力形式Navier-Stokes方程的求解.方腔层流流动比较发现,非定常应力形式比涡量形式与试验结果更加吻合,精度更高.该方法有效地解决亚格子应力项的问题,实现基于最小二乘有限元法的湍流求解.比较方腔湍流流动的试验与仿真结果,证明本文的方法具有可行性,为湍流大涡模拟计算打下基础.     

Navier-Stokes方程的有限元边界元耦合方法

主要研究了三维外部区域上具有Dirichlet边界条件的非定常Navier-Stokes方程的有限元边界元耦合方法,并分析了这一数值解的收敛速度.     

非线性薛定谔方程的高阶分裂改进光滑粒子动力学算法

为提高传统光滑粒子动力学(SPH)方法求解高维非线性薛定谔(nonlinear Schr?dinger/Gross-Pitaevskii equation, NLS/GP)方程的数值精度和计算效率,本文首先基于高阶时间分裂思想将非线性薛定谔方程分解成线性导数项和非线性项,其次拓展一阶对称SPH方法对复数域上线性导数部分进行显式求解,最后引入MPI并行技术,结合边界施加虚粒子方法给出一种能够准确、高效地求解高维NLS/GP方程的高阶分裂修正并行SPH方法.数值模拟中,首先对带有周期性和Dirichlet边界条件的NLS方程进行求解,并与解析解做对比,准确地得到了周期边界下孤立波的奇异性,且对提出方法的数值精度、收敛速度和计算效率进行了分析;随后,运用给出的高阶分裂粒子方法对复杂二维和三维NLS/GP问题进行了数值预测,并与其他数值结果进行比较,准确地展现了非线性孤立波传播中的奇异现象和玻色-爱因斯坦凝聚态中带外旋转项的量子涡旋变化过程.     

非定常Navier-Stokes方程基于两重网格离散的有限元并行算法

基于两重网格离散和区域分解技巧，提出三种求解非定常Navier-Stokes方程的有限元并行算法.算法的基本思想是在每一时间迭代步，在粗网格上采用Oseen迭代法求解非线性问题，在细网格上分别并行求解Oseen、Newton、Stokes线性问题以校正粗网格解.对于空间变量采用有限元离散，时间变量采用向后Euler格式离散.数值实验验证了算法的有效性.     

谱方法用于非定常流动计算的隐式求解

本文对谱方法用于周期性非定常流动的隐式求解方法进行了探讨,分析了影响计算稳定性和收敛速度的因素.提出了结合多重网格的隐式求解方法并对算法进行了验证,初步计算表明本文算法具有良好的稳定性和收敛速度.对于周期性非定常流动,结合本文提出的隐式求解的时域谱方法可以达到很高的精度且具有良好的计算效率.     

非定常空气流场中污染气体的对流与粘性扩散

基于离散涡方法求得的非定常、不稳定空气流场，数值求解了污染气体在空气中的对流湍流扩散方程．两者在每一个时间步耦合，以求非定常对流项．考虑速度脉动，改进了扩散模式．用SF6作为示踪气体，进行了现场测量实验．测得了各测点的SF6浓度随时间的变化．实验结果与数值仿真进行了对比，两者吻合．     

振动叶栅中三维振荡粘性流场的压力密度消去法

1前言三维非定常流场的求解是目前国内外的一个热点研究课题山。文献［2］完成了三维可压非定常欧拉流场的求解,这一方法在求解三维非定常欧拉流场时,应用了四阶fringe－Kutta方法对控制方程进行积分,用中心差分进行空间离散,采用了四阶人工粘性项来保证计算格式稳定,计算稳定性要求严格,时间步长不能大,计算时间长。本文从非定常三维粘性N－S方程组出发,通过合理的数学方法,消去压力及密度项,得到只包含振荡速度矢量项对空间的偏微分方程组,在已知定常速度场后,这一方程组很容易求解。2基本方程在以角速度为n作旋转的相对坐…     

On the vortex formation at the moving front of lightweight granular particles

The formation of vortices at a moving front of lightweight granular particles is investigated experimentally. The particles used in this study are made of polystyrene foam with three different diameters of nearly uniform size. Pairs of vortices are found to emerge at the moving front at regular intervals, thereby forming a wavy pattern. Once the vortices are produced, the flow velocity tends to increase. A simple analysis suggests the existence of a velocity boundary layer at the moving front, whose thickness increases with increasing particle diameter. The frontal radius of each vortex pair is about the size of this boundary layer; when the radius exceeds this size, the front tends to bifurcate into a train of vortices with the size of the boundary layer. The formation of twin vortices leads to a reduction in the air drag force exerted on the system, and thereby the system attains a higher flow velocity, i.e., a higher conversion rate of gravitational potential energy to the kinetic energy of the particle motion. The higher conversion rate of potential energy thus feeds back to the development of the vortex motion, resulting in the twin vortex formation.     

Vortex sound generation due to a flow impedance discontinuity on a flat surface

The sound generated by the unsteady motion of a vortex filament moving over a flat boundary with a sharp flow impedance discontinuity is studied theoretically. Theoretical results show that the vortex filament undergoes significant accelerating or decelerating motions and radiates sound at the instant when it moves across the plane of impedance discontinuity. The accelerations and decelerations of the vortex filament are shown to be the major mechanisms of sound generation. The sound so produced has a large low-frequency content such that the change in the flow impedance affects only the sound generation process but not the subsequent sound propagation to the far field.     

可压缩涡流场中空泡运动规律及声辐射特性研究

基于可压缩流体力学基本理论, 通过边界积分方程, 采用不同表面压力模型, 求解空泡在计及可压缩性的涡流场中的运动规律; 通过表面离散及坐标变换, 采用Kirchhoff动边界积分方程, 将空泡表面视为运动变形边界, 作为直接噪声源, 获得涡流场中空泡运动产生的时域声压分布; 分析了涡流场参数对空泡运动规律及声辐射特性的影响. 研究结果表明: 计及流场可压缩性, 空泡的脉动幅度会随时间减弱, 辐射声压幅值随之减小; 空泡在涡流场中会发生延展、 颈缩、 撕裂, 并在撕裂后子空泡中形成射流; 当流场中的压力减小时, 空泡运动过程中的最大半径与撕裂前的最大长度逐渐增加, 且当流场中压力较小时, 空泡撕裂时形成的子空泡增多; 空泡辐射声压的指向性较弱, 撕裂会使辐射声压产生突变, 形成极大峰值; 随着涡通量的增大或空泡数的减小, 空泡脉动周期及其诱导的辐射声压波动周期随之延长, 辐射声压峰值逐渐滞后并减小. 本文结果旨在为涡流场中空泡运动规律及声辐射特性的相关研究提供参考. 关键词： 可压缩 涡流场 空泡 声辐射     

Dynamical phase transition in two-dimensional fully frustrated Josephson junction arrays with resistively shunted junction dynamics

The dynamical phase transitions in two-dimensional fully frustrated Josephson junction arrays at zero temperature are investigated numerically with the resistively shunted junction model through the fluctuating twist boundary condition. The model is subjected to a driving current with nonzero orthogonal components i _x , i _y parallel to both axes of the square lattice. We find a roughly lattice size independent phase diagram with three dynamical phases: a pinned vortex lattice phase, a moving vortex lattice phase and a moving plastic phase. The phase diagram shows a direct transition from the pinned vortex to the moving vortex phase and the separation of the pinned vortex and the moving plastic phases. The time-dependent voltages v _x and v _y are periodic in the moving vortex lattice phase. But they are aperiodic in the moving plastic phase, resulting in non-monotonic characteristics and hysteresis in the current-voltage curves. It is found that the characteristic frequency is twice the time-averaged voltage in the moving vortex phase and around the time-averaged voltage in the plastic flow regime.Received: 29 May 2003, Published online: 2 October 2003PACS: 64.60.Ht Dynamic critical phenomena - 74.25.Sv Critical currents - 74.25.Fy Transport properties     

非紧凑低马赫数运动边界散射流动噪声的预测

提出了一种求解非紧凑低马赫数运动边界散射流动诱发噪声的预测方法。该方法首先基于运动坐标系下的连续性方程和动量方程推导得到该坐标系下的波动方程及其积分解,然后在该坐标系下采用边界元方法求解得到非紧凑运动边界表面的声压,最后将求解得到的壁面声压回代到静止或运动坐标系下的积分方程中实现对远场噪声的预测。推导得到的积分方程适用于流体与飞机机身、高速列车车身及旋转叶片等非紧凑结构边界作用诱发噪声的预测。     

Linear mechanism of surface gravity wave generation in horizontally sheared flow

An analysis is presented of a linear mechanism of surface gravity wave generation in a horizontally sheared flow in a fluid layer with free boundary. A free-surface flow of this type is found to be algebraically unstable. The development of instability leads to the formation of surface gravity waves whose amplitude grows with time according to a power law. Flow stability is analyzed by using a nonmodal approach in which the behavior of a spatial Fourier harmonic of a disturbance is considered in a semi-Lagrangian frame of reference moving with the flow. Shear-flow disturbances are divided into two classes (wave and vortex disturbances) depending on the value of potential vorticity. It is shown that vortex disturbances decay with time while the energy of wave disturbances increases indefinitely. Transformation of vortex disturbances into wave ones under strong shear is described.     

Étude aérodynamique dans une machine à compression rapide

This study deals with aerodynamics phenomena in a rapid compression machine. The mixture is initially at rest in the cylinder, and the piston is at bottom dead center; therefore there is no flow generated during the intake stroke. Moreover, flat and parallel piston and cylinder-head faces avoid squish flow. A corner vortex is generated when the piston surface, moving toward the top dead center, scrapes the boundary layer on the cylinder wall. This vortex is visualized by shadowgraph and then modeled to evaluate its characteristics. During the expansion, the gas trapped in the crevices between the piston, piston ring and cylinder wall flows back into the combustion chamber. Both crevice jet and corner vortex exist in real engines and their interactions with combustion are important: the corner vortex impairs flame kernel development and the crevice gases constitute a major source of unburnt gas.     

磁通密度对第Ⅱ类超导体磁通动力学的影响

计算了二维无序钉扎系统中磁通运动的平均速度、微分电阻、纵向电压噪声和静态结构因子.通过在不同磁通密度下的磁通运动形式,给出了磁通运动的动力学相图.研究表明,磁通晶格存在钉扎相、塑性流相、近晶流相,和运动玻璃相.在运动玻璃相中,随着驱动力的进一步增加,横向玻璃态和运动Bragg玻璃态相继出现.磁通密度增大有利于有序相的出现.当磁通密度增大到一定程度时,近晶流动相会消失.磁通运动随着外加驱动电流增大发生从塑性流动相到运动玻璃相的转变 关键词： Ⅱ类超导体 磁通动力学 运动玻璃     

Simulations of single and multiple swimmers with non-divergence free deforming geometries

We present a vortex particle method coupled with a penalization technique to simulate single and multiple swimmers in an incompressible, viscous flow in two and three dimensions. The proposed algorithm can handle arbitrarily deforming bodies and their corresponding non-divergence free deformation velocity fields. The method is validated on a number of benchmark problems with stationary and moving boundaries. Results include flows of tumbling objects and single and multiple self-propelled swimmers.     

用晶格玻尔兹曼方法研究颗粒在涡流中的运动

用晶格玻尔兹曼方法研究小颗粒在涡流中的运动.涡流由流经空腔的流体产生.用动量交换法和压力张量积分法计算颗粒在涡流中的运动轨迹、速度和角速度.最后用张量积分法计算两个不同半径的颗粒在涡流中的运动. 关键词： 晶格玻尔兹曼方法 涡流 颗粒