Revised Three-Dimensional Navier-Stokes Characteristic Boundary Conditions for Intense Reactive Turbulence

The three-dimensional Navier-Stokes characteristic boundary conditions(3D-NSCBC), although physically reasonable and popular in many applications, may encounter the instability problem in simulating complex flows, especially for large Reynolds number reactive turbulence where locally the strong reversed flow appears at the outflow boundary surfaces. In the present work, a revised 3D-NSCBC strategy is proposed based on the kinematic relation in different moving coordinate systems. Following this strategy, a systematic formulation is presented for the outflow surface with local reversed flow and can be easily extended to the coupled edge and corner boundaries. Direct numerical simulation(DNS) tests of flow with different turbulence intensities are carried out. Compared with the conventional 3D-NSCBC, the newly proposed method exhibits satisfactory performance to confine numerical instability in the strong reversed flow region. The results confirm the robustness and effectiveness of this newly proposed algorithm.     

某集装箱船推进性能的数值预报

分析了船体兴波作用对推进性能的影响,提出将兴波计算从自航船模流场计算中分离的简化计算模型,并研究了船模预报实船自航点和推进因子的数值方法,同时对KCS缩尺船模推进性能进行了预报,预报结果与试验吻合较好.通过对6条不同缩尺比的KCS自航船模的计算,分析了尺度效应对推进性能预报的影响.     

Numerical simulation of viscous flow over a grooved surface by the lattice Boltzmann method

The motivation of this work is to investigate a grooved surface’s drag reduction. The viscous flow through a two-dimensional microchannel with the grooved surface is analyzed by the lattice Boltzmann method (LBM). The effects of the grooved surface on the streamline patterns, the velocity distributions near wall region and the fluid shear stress distributions on the walls at different Reynolds numbers are studied. In addition, the influences of the groove’s geometrical parameters on the grooved surface’s drag reduction are discussed. The numerical results confirm the grooved surface’s drag reduction and present the drag reduction law of the grooved surface.     

横向风作用下高速铁路车桥系统绕流特性分析

基于不可压缩二维雷诺平均N-S方程，采用RING（renormalization group）k-ε两方程湍流模型，对德国ICE（intercity-express）列车和标准简支梁组成的车桥系统在横向风作用下的二维绕流特性进行了分析．计算网格为由三角形非结构网格和四边形结构网格组成的混合网格，采用有限体积法对微分方程进行离散，应用SIMPLE算法解决压强-速度耦合关系．根据高速铁路桥梁带挡风墙的特点，引入空腔流动研究成果对车桥系统的绕流特性进行了理论分析．结果表明，高速铁路桥梁挡风墙的挡风机理在于充分运用了空腔流的旋涡特性．     

A new algorithm for non-newtonian flow and its application in mould filling process

The simulation of injection molding process requires a stable algorithm to model the molten polymer with non-isothermal non-Newtonian property. In this paper, a staggered and iterative scheme is particularly designed to solve the velocity-pressure-temperature variables. In consideration of the polymer characteristic of high viscosity and low thermal conductivity, the non-Newtonian momentum-mass conservation equations are solved by the Crank-Nicolson method based split (CNBS) scheme, and the energy conservation equation with convective character is discretized by the characteristic Galerkin (CG) method. In addition, an arbitrary Lagrangian Eulerian (ALE) free surface tracking and mesh generation method is introduced to catch the front of the fluid flow. The efficiency of the proposed scheme is demonstrated by numerical experiments including a lid-driven cavity flow problem and an injection molding problem.     

三维波浪作用下大直径圆柱绕流的数值模拟

为探讨三维波浪与结构物的相互作用,以两相流概念、大涡模拟的不可压缩粘性流体运动方程和自由水面追踪分段线性近似的流体体积(VOF)法为基础,建立了三维波浪与结构物相互作用的数学模型;对三维波浪作用下大直径圆柱绕流进行了数值模拟,用两步边界定位法和虚拟边界力法确定波浪与结构物接触面.结果表明:大直径圆柱绕流系数的数值计算结果与理论解吻合,所建立的数学模型能够很好地模拟三维波浪与结构物的相互作用.     

模腔注射填充的全矢量显式算法

为满足注射成形模拟工业应用的实际需要，基于已有的模拟注填充显式算法，提出并实现了一种新的全矢量化显式算法，对于各类注射填充的模拟问题，此算法避免了压力场的全局求解，矩阵操作仅在单元一级进行，通过反馈修正，各时间增量步满足不可压缩条件，由于避免了所有的全局耦合求解，计算代价与节点的自由度近似成正比，可实现填充模拟过程的高效运算，容易实现并行计算，通过新算法的填充模拟结果与MINI单元的三维模拟结果比较，证实了该算法的有效性。     

智能交通最短路径Dijkstra模糊动态方法分析

提高道路网络上的动态路由选择效率是当前路网建设中的难题,现有解决路由问题的方法都是基于静态方向的研究。提出一种基于模糊的Dijkstra最短路径动态算法,以被量化的限定用户参量为语言变量,用模糊集理论提供模糊数输出在网络上的最短路径预测。通过处理模糊参数,提出用模糊数表示不同路径与边长间距离的两种比较方法。采用梯级平均综合表示法改进Dijkstra算法,并以数值的道路网络实例来说明该方法的有效性。     

An Iterative Algorithm for Computed Tomography Image Reconstruction from Limited-Angle Projections

In application of tomography imaging, limited-angle problem is a quite practical and important issue.In this paper, an iterative reprojection-reconstruction(IRR) algorithm using a modified Papoulis-Gerchberg(PG)iterative scheme is developed for reconstruction from limited-angle projections which contain noise. The proposed algorithm has two iterative update processes, one is the extrapolation of unknown data, and the other is the modification of the known noisy observation data. And the algorithm introduces scaling factors to control the two processes, respectively. The convergence of the algorithm is guaranteed, and the method of choosing the scaling factors is given with energy constraints. The simulation result demonstrates our conclusions and indicates that the algorithm proposed in this paper can obviously improve the reconstruction quality.     

基于Capton最小功率估计的方位角估计算法

在麦克风阵列信号处理中,通常假定信号源发出的信号经过传播到达不同麦克风阵元时,信号满足平面波传播模型.在此模型基础上,为解决基于麦克风阵列语音增强系统中方位角的估计问题,提出了一种基于Capton最小功率估计的方位角估计算法.仿真结果表明与传统的波束形成算法比较,该算法具有较高的方位角分辨率.     

基于Capton最小功率估计的方位角估计算法

在麦克风阵列信号处理中,通常假定信号源发出的信号经过传播到达不同麦克风阵元时,信号满足平面波传播模型.在此模型基础上,为解决基于麦克风阵列语音增强系统中方位角的估计问题,提出了一种基于Capton最小功率估计的方位角估计算法.仿真结果表明与传统的波束形成算法比较,该算法具有较高的方位角分辨率.     

注射成形填充模拟的修正算法

为提高注射成形过程数值模拟结果的准确性,采用类似迎风法的概念,通过对平流方程作为填充状态控制方程而导致模拟失真原因的分析,提出了相应的数值修正算法.该方法在高效显式数值模拟算法和有限元软件基础上,通过对模型内流体速度场的系统性操作,削弱空气流动速度场对喂料填充的不真实影响,使得填充状态的输送更基于流动前沿面后部的流场作用,并对⊥型和L型模腔进行了数值仿真.结果表明:该方法能有效抑制填充流的失真现象,实现了在流动过程中流向发生大角度变化模腔的正确填充过程.      

Numerical study for turbulent flow of drag and flow noise with wavy wall

In this paper, a numerical simulation of flow-induced noise by the low Mach number turbulent flow with a sinusoidal wavy wall was presented based on the unsteady incompressible Navier-Stokes equations and Lighthill’s acoustic analogy. Large eddy simulation (LES) was used to investigate the space-time flow field and the Smagorinsky sub-grid scale (SGS) model was introduced for turbulence model. Using Lighthill’s acoustics analogy, the flow field simulated by LES was taken as near-field sound sources and radiated sound from turbulent flow was computed by the Curle’s integral formulation under the low Mach number approximation. Both spanwise wavy wall and streamwise wavy wall with various wall wave amplitudes were discussed to investigate their effects on reducing the drag and flow noise. The relationship between flow noise and drag on the wavy wall is also studied.     

Numerical study on flow structure of a shallow laminar round jet

The evolution mechanism and characteristics of the submerged laminar round jet in a viscous homogenous shallow water layer are investigated through computational modeling. The laminar mode is used to solve the Navier-Stokes equations. In order to visualize the formation and evolution of the flow pattern, the volume of fluid (VOF) method is adopted to simulate the free surface of the water layer below the air and to trace the jet fluid. The results show that the jet forms a class of quasi-two-dimensional (Q2D) vortex structures in the ambient fluid with unequal influence from the bottom wall and free surface. The time dependence of three parameters, defined for the flow pattern as jet length, spiral radius and pattern length, is investigated quantitatively in their non-dimensional forms. Three different Reynolds numbers and two injection durations are further considered to discuss their influence on the flow pattern.     

考虑桩土作用的铁路曲线梁桥车桥耦合振动分析

以某铁路曲线多跨简支梁桥为例,讨论了考虑桩土作用的铁路曲线梁桥车桥耦合振动. 在对曲线梁桥车桥耦合振动的分析中,建立了具有35个自由度的铁路车辆曲线通过模型和动力方程,建立了曲线梁的动力模型及其动力方程;建立了一种基于激励非线性振动的数值方法,并在Windows9X/2000/XP工作环境下利用PowerStation和VisualC++完成了计算程序的编制,取得了较好的计算结果. 分析中将曲线通过的车辆和曲线梁桥分为两个由非线性轮轨接触力所联系的振动子系统,通过迭代法进行求解这两个子系统;轨道不平顺采用在给定轨道条件下的人工模拟不平顺,在分析过程中计入了不同车速对曲线通过的车辆及曲线连续梁桥振动的影响,得到一些有益的结论.     

时空有限元方法的高效计算

为提高时空有限元方法求解二维瞬态不可压缩的Navier-Stokes方程的计算效率并降低对计算机内存的需求,用行格式存储法存储大型稀疏矩阵,用Newton-Raphson迭代法求解非线性代数方程组,用无填充不完全分解预处理方法以及重启型GMRES方法求解子迭代步的线性方程组.为验证该方法的可行性,对Reynolds数为100的圆柱绕流问题进行数值模拟.采用行格式存储法的存储空间仅为等带宽存储法的3.68%.     

高雷诺数下方柱绕流特性的数值模拟

为研究高雷诺数为22 000下方柱周围流场形态及气动力特性,基于开源计算流体动力学（computational fluid dynamic,CFD）软件OpenFoam平台,采用基于动态亚格子模型的大涡模拟（large eddy simulation,LES）方法,对均匀来流作用下的方柱绕流进行了三维数值模拟.首先,通过对基于时间积分的平均积分分量的比较,验证了本数值计算的准确性;其次,深入分析了方柱周围及其尾流区的流场结构,给出了流场结构的平均和湍流特征,并在此基础上,研究了其气动力特性;最后,分析了两种长径比下表面压力的展向空间相关性.研究结果表明:雷诺数为22 000下方柱尾流区回转长度为1.37倍方柱宽度,Strouhal数为0.121,脉动升力系数为1.40;展向长度取8倍方柱宽度可更准确地获得周围湍流特性.      

IMPROVED GENETIC ALGORITHM TO OPTIMAL PORTFOLIO WITH RISK CONTROL

1IntroductionAlog-optimalinvestimentportfoliowithoutriskconstraintshasbeenintroducedandsystem-aticallystudiedinChapter15ofCoverandl'homas[11.Themodelcanbediscribedasfol-.lows:Astockmarketisrepresentedasavectorofstocksx=(X,,X,,..',X.)',X,>o,i=1,2,.-',m,wheremisthenumberofstocksandthepricerelativeXirepresentstheratiooftheclos-ingpriceofthecurrentbusinessdaytotheprevi-ousday.LetF(x),x=(x,,x2,..',x-)'eR"bethejointdistributionofvectorX.Aportfoliob=(b,,bs,..',b,)',b,>o,i=l,2,..',m,Zb,=1istheall…     

模腔注射模拟显式算法中对流项作用的评估

为研究Navier-Stokes方程中对流项对粘性流注射模拟结果的影响,在基于矢量化显式算法模拟软件中加入对流作用项.应用修正后的软件模拟两个典型的注射填充过程的算例,并与相应的实验比较.对模拟和实验结果的比较表明,对流项对粘性流注射填充过程和的其中的自由表面形状的模拟结果影响可以忽略.     

Joint Tomlinson-Harashima source and linear relay precoder design in amplify-and-forward multiple-input multiple-output two-way relay systems

Existing minimum-mean-squared-error (MMSE) transceiver designs in amplified-and-forward (AF) multiple-input multiple-output (MIMO) two-way relay systems all assume a linear precoder at the sources. Nonlinear source precoders in such a system have not been considered yet. In this paper, we study the joint design of source Tomlinson-Harashima precoders (THPs), relay linear precoder and MMSE receivers in two-way relay systems. This joint design problem is a highly nonconvex optimization problem. By dividing the original problem into three sub-problems, we propose an iterative algorithm to optimize precoders and receivers. The convergence of the algorithm is ensured since the updated solution is optimal to each sub-problem. Numerical simulation results show that the proposed iterative algorithm outperforms other algorithms in the high signal-to-noise ratio (SNR) region.