边界拟合坐标法在尾水河道计算中的应用及其验证

本文讲述应用边界拟合坐标法计算柘溪与丰满水电站尾水河道流场。同时进行了水工模型试验进行验证核对两者吻合很好。说明边界拟合坐标法精度较高。

高精度分数时延滤波器设计的边界拟合Remez算法

在卫星星间测距模拟信号仿真及实际测试中,为提高测距模拟信号精度和系统可用带宽,提出基于边界拟合Remez算法的高精度分数时延滤波器的设计算法。该算法利用Farrow结构的多项式近似思想,采用多项式拟合Remez算法设计滤波器的冲激响应边界系数,通过多相分解实现分数时延滤波器组。该算法改善了当设计的滤波器阶数较高时冲激响应边界的不连续现象,进而降低了群时延误差,提高了精度。仿真结果表明,该算法设计的滤波器的分数时延精度得到了提高,同时系统可用带宽提高近一倍,实现时需使用的乘法器数目也有明显降低。

热传导方程反问题的数值解法

文章分别应用拟边界法和积分方程方法对一类热传导反问题进行了研究,并通过数值算例验证了方法的有效性.

NUMERICAL SIMULATION AND ANALYSIS OF HORSESHOE-SHAPED VORTEX IN NEAR-WALL REGION OF TURBULENT BOUNDARY LAYER

The low-Reynolds-number full developed turbulent flow in channels is simulated using large eddy simulation（LES）method with the preconditioned algorithm and the dynamic subgrid-scale model,with a given disturbance in inlet boundary,after a short development section.The inlet Reynolds number based on momentum thickness is 670.The computed results show good agreement with direct numerical simulation（DNS）,which include root mean square fluctuated velocity distribution and average velocity distribution.It is also found that the staggered phenomenon of the coherent structures is caused by sub-harmonic.The results clearly show the formation and evolution of horseshoe vortex in the turbulent boundary layer,including horseshoe vortex structure with a pair of streamwise vortexes and one-side leg of horseshoe vortex.Based on the results,the development of the horseshoe-shaped coherent structures is analyzed in turbulent boundary layer.

Numerical Simulation for Water Entry of a Wedge at Varying Speed by a High Order Boundary Element Method

A high order boundary element method was developed for the complex velocity potential problem. The method ensures not only the continuity of the potential at the nodes of each element but also the velocity. It can be applied to a variety of velocity potential problems. The present paper, however, focused on its application to the problem of water entry of a wedge with varying speed. The continuity of the velocity achieved herein is particularly important for this kind of nonlinear free surface flow problem, because when the time stepping method is used, the free surface is updated through the velocity obtained at each node and the accuracy of the velocity is therefore crucial. Calculation was made for a case when the distance S that the wedge has travelled and time t follow the relationship s=Dtα, where D and α are constants, which is found to lead to a self similar flow field when the effect due to gravity is ignored.

Unstructured Grid Immersed Boundary Method for Numerical Simulation of Fluid Structure Interaction

This paper presents an improved unstructured grid immersed boundary method.The advantages of both immersed boundary method and body fitted grids which are generated by unstructured grid technology are used to enhance the computation efficiency of fluid structure interaction in complex domain.The Navier-Stokes equation was discretized spacially with collocated finite volume method and Euler implicit method in time domain.The rigid body motion was simulated by immersed boundary method in which the fluid and rigid body interface interaction was dealt with VOS(volume of solid) method.A new VOS calculation method based on graph was presented in which both immersed boundary points and cross points were collected in arbitrary order to form a graph.The method is verified with flow past oscillating cylinder.

Multigrid fictitious boundary method for incompressible viscous flows around moving airfoils

In this paper, an efficient multigrid fictitious boundary method （MFBM） coupled with the FEM solver package FEATFLOW was used for the detailed simulation of incompressible viscous flows around one or more moving NACA0012 airfoils. The calculations were carded on a fixed multigrid finite element mesh on which fluid equations were satisfied everywhere, and the airfoils were allowed to move freely through the mesh. The MFBM was employed to treat interactions between the fluid and the airfoils The motion of the airfoils was modeled by Newton-Euler equations. Numerical results of experiments verify that this method provides an efficient way to simulate incompressible viscous flows around moving airfoils.

Numerical Simulation of the Flow around Two-dimensional Partially Cavitating Hydrofoils

In the present study, a new approach is applied to the cavity prediction for two-dimensional （2D） hydrofoils by the potential based boundary element method （BEM）. The boundary element method is treated with the source and doublet distributions on the panel surface and cavity surface by usethe of the Dirichlet type boundary conditions. An iterative solution approach is used to determine the cavity shape on partially cavitating hydrofoils. In the case of a specified cavitation number and cavity length, the iterative solution method proceeds by addition or subtraction of a displacement thickness on the cavity surface of the hydrofoil. The appropriate cavity shape is obtained by the dynamic boundary condition of the cavity surface and the kinematic boundary condition of the whole foil surface including the cavity. For a given cavitation number the cavity length of the 2D hydrofoil is determined according to the minimum error criterion among different cavity lengths, which satisfies the dynamic boundary condition on the cavity surface. The NACA 16006, NACA 16012 and NACA 16015 hydrofoil sections are investigated for two angles of attack. The results are compared with other potential based boundary element codes, the PCPAN and a commercial CFD code （FLUENT）. Consequently, it has been shown that the results obtained from the two dimensional approach are consistent with those obtained from the others.

Iterative Method on Well Bore Boundary in Numerical Modeling of Variably Saturated Flow

An iterative method was developed for incorporating the well bore boundary into the finite difference model of water flow in variably saturated porous media. Six cases were presented involving groundwater pumping or injection to demonstrate the advantages of the iterative method over the traditional method. For the iterative method, the total flux gradually approached the well discharge and the flux profile was non-uniform. And the iterative method took into account the variation of well bore water table. Compared to the traditional method, the iterative method can simulate the variably saturated flow caused by pumping or injection more realistically.

3D Numerical Modeling of Wave Forces on Tandem Fixed Cylinders Using the BEM

In this paper a 3D numerical model was developed to study the complicated interaction between waves and a set of tandem fixed cylinders.The fluid was considered to be inviscid and irrotational.Therefore,the Helmholtz equation was used as a governing equation.The boundary element method（BEM） was adopted to discretize the relevant equations.Open boundaries were used in far fields of the study domain.Linear waves were generated and propagated towards tandem fixed cylinders to estimate the forces applied on them.Special attention was paid to consideration of the effect on varying non-dimensional cylinder radius and distance between cylinders,ka and kd on forces and trapped modes.The middle cylinder wave forces and trapped modes in a set of nine tandem cylinders were validated utilizing analytical data.The comparisons confirm the accuracy of the model.The results of the inline wave force estimation on n tandem cylinders show that the critical cylinder in the row is the middle one for odd numbers of cylinders.Furthermore the results show that the critical trapped mode effect occurs for normalized cylinder radiuses close to 0.5 and 1.0.Finally the force estimation for n tandem cylinders confirms that force amplitude of the middle cylinder versus normalized separation distance fluctuates about that of a single cylinder.

An element-free Galerkin method for ground penetrating radar numerical simulation

An element-free Galerkin method(EFGM) is used to solve the two-dimensional(2D) ground penetrating radar(GPR)modelling problems, due to its simple pre-processing, the absence of elements and high accuracy. Different from element-based numerical methods, this approach makes nodes free from the elemental restraint and avoids the explicit mesh discretization. First, we derived the boundary value problem for the 2D GPR simulation problems. Second, a penalty function approach and a boundary condition truncated method were used to enforce the essential and the absorbing boundary conditions, respectively. A three-layered GPR model was used to verify our element-free approach. The numerical solutions show that our solutions have an excellent agreement with solutions of a finite element method(FEM). Then, we used the EFGM to simulate one more complex model to show its capability and limitations. Simulation results show that one obvious advantage of EFGM is the absence of element mesh, which makes the method very flexible. Due to the use of MLS fitting, a key feature of EFM, is that both the dependent variable and its gradient are continuous and have high precision.

Study of numerical simulation on dual-frequency IP method with FEM

Using the finite element method and Cole-Cole model for dual-frequency IP method to research numerical simulation, the authors introduced the fundamental principle of the dual-frequency IP method and the boundary value problem and variational equations, then replaced the complex resistivity of the model with the Cole-Cole model's parameters under ignoring the EM effect. Through solving the last linear equations, electric potentials of all the model's points were obtained. With changing model's parameters, the authors got different curves of the Fs and phases. According to the results of the simulation, the algorithm is proved to be correct and adaptable.

Calculation of Turbulent Flow with Pseudo-spectral Method

This paper deals with the numerical simulation of incompressible turbulent boundary flow of a flat plate with the pseudo-spectral matrix method. In order to appear more than 10 nodes in the turbulent base-stratum and transition of 43×43 computational grids,a coordinate transformation is put up from physical panel to computational panel. Several zero turbulent models are computed comparatively. The results are credible when comparing with the previous methods.

Numerical Investigation of an Oscillating Airfoil using Immersed Boundary Method

Based on the immersed boundary method, a numerical simulation for an oscillating airfoil is established and a preliminary analysis of the oscillating airfoil is presented with an emphasis on the physical understanding of fluid-structure interaction. In order to validate the method, two simulation cases： oscillating circular cylinder at low K-C number and two degrees of freedom oscillating cylinder are carded out first and the results are in good agreement with the previous re：searches. In the oscillating airfoil simulation, it is found that the reduced velocity U^＊. is a very sensitive factor and especially U^＊-2.8 is the critical stable boundary in the present work. The method shows the predominance of time saving in computational process for such a complicated fluid-structure interac- tion problem.

Simulation of Unsteady Flows with Fluctuating Inflow or Outflow Boundary Conditions

This paper describes the numerical simulation of unsteady flows due to incoming wakes and/or varying back pressure,The solution method is based upon the one-step finite-volume TVD Lax-Wendroff scheme.Dual time-step approach and multigrid algorithm are adopted to improve the computational efficiency of the baseline scheme.Numerical results for the transonic unsteady flow in a channel bump and the unsteady flow in a flat plate cascade and the VKI cascade are presented.