浅水回流的混合有限分析解

混合有限分析法是一种在局部矩形单元上进行离散的数值格式,为了适应非规则边界,建立了Sigma坐标系的浅水回流数学模型。采用1)风引起的回流,2)密度驱动的回流,3)假潮,来检验数学模型和数值方法。计算结果和相应的分析解的比较表明模型和方法是可行的有效的。该技术可用于近海水域的水流和水质的数值模拟。     

用时变雷诺方程模型模拟孤立波与半圆型防波堤的相互作用

以时变雷诺方程为控制方程,用k-ε模型封闭该方程,采用体积函数(VOF)方法来跟踪波动自由表面,建立了二维垂向波浪数学模型,并用已有的实验资料进行了验证．随后用该模型模拟了半圆型防波堤与孤立波在淹没、平顶水位、完全露顶且不越浪3种典型工况下的相互作用过程．得到了半圆堤附近的流场、压强场和波面的变形过程．结果表明,在淹没状态下,半圆堤背浪面的底部会产生涡旋;平顶水位时,由于越浪的冲击作用,在半圆堤的背浪面将逐渐形成一对较大的涡旋,而半圆堤背浪面的底部,速度始终相对较小;而在露顶不越浪时,半圆堤的迎浪面会出现波浪的二次爬升的现象．为进一步研究结构物附近的污染物的输移扩散和泥沙运动提供基础．     

Depth-averaged modeling of free surface flows in open channels with emerged and submerged vegetation

The resistance induced by vegetation on the flow in a watercourse should be considered in projects of watercourse management and river restoration. Depth-averaged numerical model is an efficient tool to study this problem. In this study, a depth-averaged model using the finite volume method on a staggered curvilinear grid and the SIMPLEC algorithm for numerical solution is developed for simulating the hydrodynamics of free surface flows in watercourses with vegetation. For the model formulation the vegetation resistance is treated as a momentum sink and represented by a Manning type equation, and turbulence is parameterized by the k–ε equations. An analytical equation is derived to represent the resistance induced by submerged vegetation by an equivalent Manning roughness coefficient. Numerical simulation is carried out for the flow in an open channel with a 180° bend, and the flow in a curved open channel partly covered by emerged vegetation, as well as the flow in a straight trapezoidal channel with submerged vegetation. The agreement between the computed results and the measured data is generally good, showing that the resistance due to emerged or submerged vegetation can be represented accurately by the Manning roughness equation. The computed results demonstrate that the depth-averaged modeling is a reasonable and efficient tool to study flows in watercourses with vegetations.     

Modelling gas injection of a Peirce–Smith-converter

A commercial CFD-code PHOENICS was used to solve isothermal flow field of gas and liquid in a Peirce–Smith-converter. An Euler–Euler based algorithm was chosen for modelling fluid dynamics and evaluating controlling forces of a submerged gas injection. Predictions were made with a k–ε turbulence model in the body fitted coordinate system. The model has been verified with a 1/4 scale water model, and a parametric study with the mathematical model of submerged gas injection was made for the PS-process and the ladle injection processes. Limits of the modelling technique used were recognised, but calculated results indicate that the present model predicts the general flow field with reasonable accuracy. Predicted bubble distribution, pattern of the flow field and magnitude of flow velocities were used to evaluate scaling factors of physical models and general flow conditions of an industrial PS-converter.     

A numerical procedure for viscous free surface flows

We present a numerical procedure for two-dimensional unsteady viscous free surface flow problems with surface tension. The procedure is based on a finite difference approach to a primitive variable formulation; a coordinate transformation is used to transform the irregularly shaped flow domain onto a fixed rectangular domain. The procedure is tested on a standing wave problem and a cavity flow problem with a free surface. Satisfactory numerical solutions are obtained for both problems for Reynolds numbers up to 200.     

Finite element modelling of surface waves in a channel

A variational formulation of the vertically-integrated differential equations for free surface wave motion is presented. A finite element model is derived for solving this nonlinear system of hydrodynamic equations. The time integration scheme employed is discussed and the results obtained demonstrate its good stability and accuracy.Several applications of the model are considered: the first problem is an open channel of uniform depth and the second an open channel of linearly varying depth. The ‘inflow’ boundary condition is prescribed in terms of the velocity which represents a wavemaker and/or a flow source, while the ‘outflow’ boundary condition is specified in terms of the water elevation. The outflow condition is adjusted for two cases, a reflecting boundary (finite channel) and a non-reflecting boundary (open-ended channel). The latter boundary condition is examined in some detail and the results obtained show that the numerical model can produce the non-reflecting boundary that is similar to the analytical radiation condition for waves. Computational results for a third problem, involving wave reflection from a submerged cylinder, are also presented and compared with both experimental data and analytical predictions.The simplicity and the performance of the computational model suggest that free surface waves can be simulated without excessively complicated numerical schemes. The ability of the model to simulate outflow boundary conditions properly is of special importance since these conditions present serious problems for many numerical algorithms.     

Modelling free surfaces in oscillating pipe flows

An oscillating pipe flow with a free surface is investigated numerically and experimentally. The pipe diameter is 12mm. Due to this small diameter capillary forces play an important role. Therefore special attention has to be paid to the flow field near the free surface. The numerical model is based on the fundamental flow equations. The free surface is resolved according to the volume-of-fluid method. The model equations are solved on a moving grid. In the experiment, pictures of the flow field are taken near the free surface. The effects occuring near the interface will be presented here. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Improving multi-block sigma-coordinate for 3D simulation of sediment transport and steep slope bed evolution

This paper aims at developing a multi-block sigma-coordinate to simulate morphological evolutions. The developed multi-block sigma-coordinate can represent a steep slope topography smoothly with different numbers of horizontal layers, without producing any truncation error and artificial flux (PGFE). The multi-block sigma–coordinate can easily increase the depth-direction resolution of the sediment transport module in the sub-regions, without the aggregation of computational points in the shallow areas which may be caused by using high resolution over the entire domain. The model is beneficial for long-term simulation of morphological evaluations in lakes where the bed slope near delta region (the sedimentary area which forms where a river enters a lake/ocean) is mostly steep, and delta keeps advancing down the lake.The multi-block sigma approach at the block interface, where the number of horizontal layers varies, allocates the flux to the neighboring cells according to two essential factors, namely “the common border length” and “satisfying continuity” which lead to defining the virtual cells. A series of numerical tests have been performed. Comparison between the numerical results, the analytical solutions, and the experimental data demonstrated an appreciable accuracy, a satisfactory performance and the efficiency of this scheme.     

2D numerical simulation of submerged hydraulic jumps

Hydraulic jumps are usually used to dissipate energy in hydraulic engineering. In this paper, the turbulent submerged hydraulic jumps are simulated by solving the unsteady Reynolds averaged Navier–Stokes equations along with the continuity equation and the standard k–? equations for turbulence modeling. The Lagrangian moving grid method is employed for the simulation of the free surface. In the developed model, kinematic free-surface boundary condition is solved simultaneously with the momentum and continuity equations, so that the water elevation can be obtained along with velocity and pressure fields as part of the solution. Computational results are presented for Froude numbers ranging from 3.2 to 8.2 and submergence factors ranging from 0.24 to 0.85. Comparisons with experimental measurements show that numerical model can simulate the velocity field, variation of free surface, maximum velocity, Reynolds shear and normal stresses at various stations with reasonable accuracy.     

Rayleigh wave scattering by a plane strip in a deep ocean

A theoretical formulation to study the problem of scattering of Rayleigh waves due to the presence of a rigid plane strip in a deep ocean is presented. A rigid plane strip (0 ≤z ≤ H, 0 ≤x ≤ l) is fixed in the surface of the ocean occupyingz ≥ 0. Fourier transformation and Wiener-Hopf technique are used to arrive at the solution. The scattered Rayleigh waves behave as cylindrical waves emerging out of the corner of the strip and its image in the free surface of the ocean. The scattered waves are obtained in terms of Bessel functions whose behaviour near and far from the strip is well-known. The numerical calculations for the scattered waves show that their amplitude increases rapidly for a small increase in the value of the wave number. Scattering of Rayleigh waves due to a thin plane vertical barrier and a thin barrier in the free surface of the ocean has been considered as the special cases.     

A numerical investigation of the nonlinear wave stability of vertical thermal boundary layer flow in a porous medium

A two-dimensional, nonlinear, time-dependent, elliptic, numerical method coupled with an appropriate coordinate transformation is used to investigate the stability of free convection induced by an isothermally heated semi-infinite surface embedded in a fluid-saturated porous medium. It is found that the basic boundary layer flow is stable even to large amplitude disturbance for nondimensional distances of up to 1024 from the leading edge of the heated surface.     

Three-dimensional non-Darcy flow analysis using a hybrid numerical model

A hybrid numerical model is developed for the simulation of three-dimensional, unsteady non-Darcy flow through an unconfined aquifer. The major problem in analysing flow through unconfined aquifers is that they involve two boundaries, namely a surface of seepage and a free surface, the location of which is not known beforehand. The model that is presented here determines these boundaries via a two stage modelling technique. In the first stage a one-dimensional finite difference model is used to estimate the surface of seepage height whereas in the second stage a vertically integrated finite element model determines the free surface solution within the flow domain. A comparison between numerical and experimental results is included which indicates the sensitivity of the numerical solution to the selected aquifer parameters, particularly to those associated with the determination of the height of the surface of seepage.     

Numerical simulation of the filling behavior of a Non-Newtonian Fluid

A numerical model for free surface flows of non-newtonian liquids which are injected into a cavity is presented. These flows are regarded as a basic model of injection molding. Model experiments of the injection process are performed with a water-based gel. The flow equations are integrated according to the finite-volume-method. The volume of fluid method (VoF) is employed in order to describe the free surface flow of two incompressible phases, the phase interface is resolved by the method of geometric reconstruction. The Herschel-Bulkley model is used in order to describe shear-thinning behavior of the molding material and the effects of a yielding point. Different patterns of the filling flow depending on the injection parameters are evident in the experiment and the simulation. They are characterized and arranged with respect to the similarity parameters of the flow. Again, the results of the simulation are found to agree well with the experimental observations. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Axi symmetric 2D simulation and numerical heat transfer characteristics for the calibration furnace in a rectangular enclosure

This paper presents axi symmetric 2D numerical investigation of the spherical thermocouple calibration furnace in a rectangular enclosure. The focus is on the flow structure inside the Saturn (a hollow spherical cavity), external flow behavior due to annulus block heating and the surface temperature uniformity. Mesh sensitivity analysis is adopted to extract the mesh with minimum number of nodes but with fast convergent finite element solution. The basic strategy here is that temperature perturbation error at a single point instead of a single element contributed to the total perturbation error qualitatively remains the same. Agreement between numerical simulation results and the experiment results is good with a maximum temperature deviation 10 °C for the cavity temperature 400 °C. Finally, standard numerical temperature uncertainty due to variation in thermal conductivity is computed through the sensitivity coefficient using uncertainty analysis.     

Generalized Lagrange Jacobi-Gauss-Lobatto vs Jacobi-Gauss-Lobatto collocation approximations for solving (2 + 1)-dimensional Sine-Gordon equations

The Sine-Gordon (SG) equations are very important in that they can accurately model many essential physical phenomena. In this paper, the Jacobi-Gauss-Lobatto collocation (JGL-C) and Generalized Lagrange Jacobi-Gauss-Lobatto collocation (GLJGL-C) methods are adopted and compared to simulate the (2 + 1)-dimensional nonlinear SG equations. In order to discretize the time variable t, the Crank-Nicolson method is employed. For the space variables, two numerical methods based on the aforementioned collocation methods are applied. Furthermore, error estimation for both methods is provided. The present numerical method is truly effective, free of integration and derivative, and easy to implement. The given examples and the results assert that the GLJGL-C method outperforms the JGL-C method in terms of computation speed. Also, the presented methods are very valid, effective, and reliable.     

Modeling of Free Convection Boundary Layer Flow on a Solid Sphere with Newtonian Heating

In this paper, the mathematical model of free convection boundary layer flow on a solid sphere with Newtonian heating, in which the heat transfer from the surface is proportional to the local surface temperature, is considered. The transformed boundary layer equations are solved numerically using an efficient numerical scheme known as the Keller-box method. Numerical solutions are obtained for the local wall temperature, the local skin friction coefficient, as well as the velocity and temperature profiles. The features of the flow and heat transfer characteristics for different values of the Prandtl number Pr and conjugate parameter γ are analyzed and discussed.     

A variable boundary method for modelling two dimensional free surface flows with moving boundaries

A coordinate transformation method is proposed for modelling unsteady, depth-averaged shallow water equations for a open channel flow with moving lateral boundaries. The transformation technique which maps the changing domain onto a fixed domain and solves the governing equations in the mapped domain, facilitates the numerical treatment of an irregular boundary configuration. The transformed equations are solved on a staggered grid with a conditionally stable, explicit finite difference scheme. Several numerical experiments are carried out corresponding to different situations, viz., flow with constant discharge, flow with constant discharge and a closed boundary at the downstream, flow in a converging channel with constant discharge and finally flow with varying discharge. The experiments are used to verify the model ability to predict free surface elevation, circulatory pattern and displacement of the boundaries. The simulated results such as displaced area, depth, displacement–time and flow-field are used to evaluate the effects of excess discharge at the upstream on the movement of lateral boundaries.     

Modeling of van der Waals force with smoothed particle hydrodynamics: Application to the rupture of thin liquid films

The rupture of thin liquid films driven by the van der Waals force is of significance in many engineering processes, and most previous studies have relied on the lubrication approximation. In this paper, we develop a smoothed particle hydrodynamics (SPH) representation for the van der Waals force and simulate the rupture of thin liquid films without resort to lubrication theory. The van der Waals force in SPH is only imposed on one layer, i.e., the outermost layer of fluid particles, where a weighting function is deployed to evaluate the contributions of particles on or near the interface. However, to obtain an accurate hydrostatic pressure in reaction to the van der Waals force, a smaller smoothing length is used for the calculation of the weighting function than that used for SPH discretizations of the bulk fluid. The same surface particles are also used to model the surface tension. To deal with the rupture of a thin liquid film with a very small aspect ratio ε (ε = thickness/length), a coordinate transformation is introduced to shrink the length of the liquid film to achieve accurate numerical resolution with a manageable number of particles. As verifications of our physical model and numerical algorithm, we simulate the hydrostatic pressure in a stationary film and the relaxation of an initially square droplet and compare the SPH results with the analytical solutions. The method is then applied to simulate the rupture of thin liquid films with moderate and small aspect ratios (ε = 0.5 and 0.005). The convergence of the method is verified by refining particle spacing to four different levels. The effect of the capillary number on the rupture process is analyzed.     

Numerical solution of a Rayleigh-Taylor instability problem

The numerical solution of a Rayleigh-Taylor instability problem where an inviscid liquid of finite depth is accelerated into a gas of semi-infinite extent is obtained by transforming the irregular flow domain into a rectangular domain by a coordinate transformation. The free surface equation is solved by a Crank-Nicolson procedure. The boundary condition at the free surface for the velocity potential ø which contains the time derivative of ø is also treated by an implicit scheme. The numerical results agree well with those obtained by higher order perturbation analysis.     

A combined BIE-FE method for the Stokes equations

A numerical algorithm for the biharmonic equation in domainswith piecewise smooth boundaries is presented. It is intendedfor problems describing the Stokes flow in the situations whereone has corners or cusps formed by parts of the domain boundaryand, due to the nature of the boundary conditions on these partsof the boundary, these regions have a global effect on the shapeof the whole domain and hence have to be resolved with sufficientaccuracy. The algorithm combines the boundary integral equationmethod for the main part of the flow domain and the finite-elementmethod which is used to resolve the corner/cusp regions. Twoparts of the solution are matched along a numerical ‘internalinterface’ or, as a variant, two interfaces, and theyare determined simultaneously by inverting a combined matrixin the course of iterations. The algorithm is illustrated byconsidering the flow configuration of ‘curtain coating’,a flow where a sheet of liquid impinges onto a moving solidsubstrate, which is particularly sensitive to what happens inthe corner region formed, physically, by the free surface andthe solid boundary. The ‘moving contact line problem’is addressed in the framework of an earlier developed interfaceformation model which treats the dynamic contact angle as partof the solution, as opposed to it being a prescribed functionof the contact line speed, as in the so-called ‘slip models’.