Effects of Roughness Elements Distribution on Overland Flow Resistance

Roughness elements are various in a mountain area; they include gravel and ground surface vegetation that often result in surface friction drag to resist overland flows. The variation and characteristics of flow resistance strongly impact the overland flow process and watershed floods. In view of the universal existence of natural vegetation, such as Chlorophytum malayense(CM) or Ophiopogon bodinieri(OB), and the sand-gravel bed of the river channel, it is important to understand the role of different types of roughness elements in flow resistance. This study was performed to investigate and compare through flume experiments the behaviors of overland flow resistance by the reaction of multi-scale configuration of different roughness elements. The result showed that the resistance coefficient gradually reduced versus the increase of flow rate in unit width and tended to be a constant when q = 3.0 l/s.m, Fr = 1.0, and Re = 4000 for slopes of 6 to 10 degrees. The gap of the vegetated rough bed and the gravel rough bed is limited to the same as the gap of the two types of vegetation, CM and OB. It was noted that the vegetation contributed to the increase in form resistance negatively and may lead to the mean resistance on decrease. To classify the flow pattern, the laminar flows were described by DarcyWeisbach's equation. In the study the f-Re equation of vegetated bed was developed with f ?5000 Re.The friction coefficient for laminar flows can be regarded as the critical value for identifying the transformation point of the flow pattern.

CONSTRUCTION OF MODIFIED TAYLOR－GALERKIN FINITE ELEMENTS AND ITS APPLICATION IN COMPRESSIBLE FLOW COMPUTATION

In this paper we begin with Taylor-Galerkin Finite Elements,then n improve it completely and finally construct Modified Taylor-Galerkin Finite Elements.Computation is done by two methods to study two kinds of subsonic and supersonic flow.field.The results show that the new one is much better than the old one.

INFLUENCE OF SURFACE-ACTIVE ELEMENTS ON FLUID FLOW OF MIG WELDPOOL

A mathematical model describing the behavior of metal inert gas (MIG) welding isformulated in the paper. By means of numerical simulation, the influence of surface-active elements on fluid flow of MIG weldpool is studied. The calculation resultsshow that by adding surface-active elements, the fluid flow behavior is drasticallychanged and the flow fluid flows from lower to upper in vertical direction at the rearof weldpool (w>0). The physical phenomenon is explained from the viewpoint of fluidflow behavior of weldpool that the properties of weld metal is greatly improved andthe content of diffused hydrogen is reduced, thus providing a basis for developing newwelding materials.

Phase-field modeling of dendritic growth under forced flow based on adaptive finite element method

A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic growth was simulated from undercooled nickel melt under the forced flow. The simulation results show that the asymmetry behavior of the dendritic growth is caused by the forced flow. When the flow velocity is less than the critical value, the asymmetry of dendrite is little influenced by the forced flow. Once the flow velocity reaches or exceeds the critical value, the controlling factor of dendrite growth gradually changes from thermal diffusion to convection. With the increase of the flow velocity, the deflection angle towards upstream direction of the primary dendrite stem becomes larger. The effect of the dendrite growth on the flow field of the melt is apparent. With the increase of the dendrite size, the vortex is present in the downstream regions, and the vortex region is gradually enlarged. Dendrite tips appear to remelt. In addition, the adaptive finite element method can reduce CPU running time by one order of magnitude compared with uniform grid method, and the speed-up ratio is proportional to the size of computational domain.

Theoretical Analysis and Experimental Verification on Flow Field of Piezoelectric Pump with Unsymmetrical Slopes Element

Regular valveless piezoelectric pumps have rectifying elements outside their chambers to produce net flow. These rectifying elements outside the chamber will increase the overall volume of the pump and prevent its minimization. Valveless piezoelectric pump with unsymmetrical slopes elements（USE）, proposed in this paper, differs from other valveless pumps in that it is easy to be minimized by developing the chamber bottom as such a rectifying element. In this research, the working principle of the proposed pump was analyzed first. Numerical models were thereby established and numerical simulation was conducted to the chamber flow field with the method of time-dependent velocity. The effects of the USEs on the flow field in the chamber were shown clearly in simulation. And the particular feature of flow field in the chamber was discovered. It behaves a complex flow field, in which strong turbulent occurs companying a lot of vortexes in different directions and different sizes. This feature is just opposite to what regular piezoelectric pumps expect： a moderate flow field. The turbulent flow could be used to have different liquids stirred and well mixed in the chamber to produce homogeneous solution, emulsion or turbid liquid. Meanwhile, numerical simulation also presents the effect of the angles difference of the two slopes upon the flow field, and upon the flow rate of the pump, which fits to the theoretical analysis. Experiments with the proposed pump were also conducted to verify the numerical results. In these experiments, six USEs with different slope angles were used for efficiency tests, which proved the validity and reliability of the numerical analysis. The data obtained from numerical analysis agree well with that from the experiments. The errors ranged from 4.4% to 14.8% with their weighted average error being 9.7%.

COUPLED SIMULATION OF 3D ELECTRO-MAGNETO-FLOW FIELD IN HALL-HEROULT CELLS USING FINITE ELEMENT METHOD

Two full 3D steady mathematical models are developed by finite element method （FEM） to calcalate coupled physics fields. the electro-magnetic model is built and solved first and so is the fluid motion model with the acquired electromagnetic force as source body forces in Navier-Stokes equations. Effects caused by the ferromagnetic shell, busbar system around, and open boundary problem as well as inside induced current were considered in terms of the magnetic field. Furthermore, a new modeling method is found to set up solid models and then mesh them entirely with so-called structuralized grids, namely hex-mesh. Examples of 75kA prebaked cell with two kinds of busbar arrangements are presented. Results agree with those disclosed in the literature and confirm that the coupled simulation is valid. It is also concluded that the usage of these models facilitates the consistent analysis of the electric field to magnetic field and then flow motion to the greater extent, local distributions of current density and magnetic flux density are very much dependent on the cell structure, the steel shell is a shield to reduce the magnetic field and flow pattern is two dimensional in the main body of the metal pad.

Finite Element Numerical Simulation and PIV Measurement of Flow Field inside Metering-in Spool Valve

The finite element method （FEM） and particle image velocimetry （PIV） technique are utilized to get the flow field along the inlet passage, the chamber, the metering port and the outlet passage of spool valve at three different valve openings. For FEM numerical simulation, the stream function ψ-vorticity ω forms of continuity and Navier-Stokes equations are employed and FEM is applied to discrete the equations. Homemade simulation codes are executed to compute the values of stream function and vorticity at each node in the flow domain, then according to the correlation between stream function and velocity components, the velocity vectors of the whole field are calculated. For PIV experiment, pulse Nd： YAG laser is exploited to generate laser beam, cylindrical and spherical lenses are combined each other to produce 1.0 mm thickness laser sheet to illuminate the object plane, Polystyrene spherical particle with diameter of 30-50 μm is seeded in the fluid as a tracing particles, Kodak ES 1.0 CCD camera is employed to capture the images of interested, the images are processed with fast Fourier transform （FFT） cross-correlation algorithm and the processing results is displayed. Both results of numerical simulation and PIV experimental show that there are three main areas in the spool valve where vortex is formed. Numerical results also indicate that the valve opening have some effects on the flow structure of the valve. The investigation is helpful for qualitatively analyzing the energy loss, noise generating, steady state flow forces and even designing the geometry structure and flow passage.

A control volume based finite element method for simulating incompressible two-phase flow in heterogeneous porous media and its application to reservoir engineering

Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.

ANISOTROPIC MULTISTAGE FINITE ELEMENT METHOD FOR TWO DIMENSIONAL VISCOUS TRANSONIC FLOW IN TURBOMACHINERY

A new method based on the anisotropic tensor force finite element and Taylor-Galerkin finite element is presented in the present paper.Its application to two-dimensional viscous transonic flow in turbomachinery improves the conver- gence rate and stability of calculation,and the results obtained agree well with the experimental measurements.

A mixed finite element scheme for viscoelastic flows with XPP model

A mixed finite element formulation for viscoelastic flows is derived in this paper, in which the FIC （finite incremental calculus） pressure stabilization process and the DEVSS （discrete elastic viscous stress splitting） method using the Crank-Nicolson-based split are introduced within a general framework of the iterative version of the fractional step algorithm. The SU （streamline-upwind） method is particularly chosen to tackle the convective terms in constitutive equations of viscoelastic flows. Thanks to the proposed scheme the finite elements with equal low-order interpolation approximations for stress-velocity-pressure variables can be successfully used even for viscoelastic flows with high Weissenberg numbers. The XPP （extended Pom-Pom） constitutive model for describing viscoelastic behaviors is particularly integrated into the proposed scheme. The numerical results for the 4：1 sudden contraction flow problem demonstrate prominent stability, accuracy and convergence rate of the proposed scheme in both pressure and stress distributions over the flow domain within a wide range of the Weissenberg number, particularly the capability in reproducing the results, which can be used to explain the ＂die swell＂ phenomenon observed in the polymer injection molding process.

Computation of Viscous Uniform and Shear Flow over A Circular Cylinder by A Finite Element Method

The incompressible viscous uniform and shear flow past a circular cylinder is studied. The two-dimensional Navier-Stokes equations are solved by a finite element method. The governing equations are discretized by a weighted residual method in space. The stable three-step scheme is applied to the momentum equations in the time integration. The numerical model is firstly applied to the computation of the lid-driven cavity flow for its validation. The computed results agree well with the measured data and other numerical results. Then, it is used to simulate the viscous uniform and shear flow over a circular cylinder for Reynolds numbers from 100 to 1000. The transient time interval before the vortex shedding occurs is shortened considerably by introduction of artificial perturbation. The computed Strouhal number, drag and lift coefficients agree well with the experimental data. The computation shows that the finite element model can be successfully applied to the viscous flow problem.

A hybrid vertex-centered finite volume/element method for viscous incompressible flows on non-staggered unstructured meshes

This paper proposes a hybrid vertex-centered fi- nite volume/finite element method for solution of the two di- mensional （2D） incompressible Navier-Stokes equations on unstructured grids. An incremental pressure fractional step method is adopted to handle the velocity-pressure coupling. The velocity and the pressure are collocated at the node of the vertex-centered control volume which is formed by join- ing the centroid of cells sharing the common vertex. For the temporal integration of the momentum equations, an im- plicit second-order scheme is utilized to enhance the com- putational stability and eliminate the time step limit due to the diffusion term. The momentum equations are discretized by the vertex-centered finite volume method （FVM） and the pressure Poisson equation is solved by the Galerkin finite el- ement method （FEM）. The momentum interpolation is used to damp out the spurious pressure wiggles. The test case with analytical solutions demonstrates second-order accuracy of the current hybrid scheme in time and space for both veloc- ity and pressure. The classic test cases, the lid-driven cavity flow, the skew cavity flow and the backward-facing step flow, show that numerical results are in good agreement with the published benchmark solutions.

SUPG finite element method based on penalty function for lid-driven cavity flow up to Re = 27500

A streamline upwind/Petrov-Galerkin （SUPG） finite element method based on a penalty function is pro- posed for steady incompressible Navier-Stokes equations. The SUPG stabilization technique is employed for the for- mulation of momentum equations. Using the penalty function method, the continuity equation is simplified and the pres- sure of the momentum equations is eliminated. The lid-driven cavity flow problem is solved using the present model. It is shown that steady flow simulations are computable up to Re = 27500, and the present results agree well with previous solutions. Tabulated results for the properties of the primary vortex are also provided for benchmarking purposes.

A Parallel Boundary Element Formulation for Tracking Multiple Particle Trajectories in Stoke’s Flow for Microfluidic Applications

A new formulation for tracking multiple particles in slow viscous flow for microfluidic applications is presented.The method employs the manipulation of the boundary element matrices so that finally a system of equations is obtained relating the rigid body velocities of the particle to the forces applied on the particle.The formulation is specially designed for particle trajectory tracking and involves successive matrix multiplications for which SMP(Symmetric multiprocessing)parallelisation is applied.It is observed that present formulation offers an efficient numerical model to be used for particle tracking and can easily be extended for multiphysics simulations in which several physics involved.

The Finite Element Analysis for Parallel-wire Capacitance Probe in Small Diameter Two-phase Flow Pipe

This paper presents a novel capacitance probe, i.e., paraUel-wire capacitance probe （PWCP）, for two-phase flow measurement. Using finite element method （FEM）, the sensitivity field of the PWCP is investigated and the optimum sensor geometry is determiend in term of the characterisitc parameters. Then, the response of PWCP for the oil-water stratified flow is calculated, and it is found the PWCP has better linearity and sensitivity to the variation of water-layer thickness, and is almost independant of the angle between the oil-water interface and the sensor electrode. Finally, the static experiment for oil-water stratified flow is carried out and the calibration method of liquid holdup is presented.

Towards a Micromechanical Understanding of Landslides—Aiming at a Combination of Finite and Discrete Elements with Minimal Number of Degrees of Freedom

In this paper, we propose a combination of discrete elements for the soil and finite elements for the fluid flow field inside the pore space to simulate the triggering of landslides. We give the details for the implementation of third order finite elements (“P2 with bubble”) together with polygonal discrete elements, which allows the formulation with a minimal number of degrees of freedom to save computer time and memory. We verify the implementation with several standard problems from computational fluid dynamics, as well as the decay of a granular step in a fluid as test case for complex flow.

Boundary Element Analysis (Laplace Transform Solution) of Groundwater Unsteady Flow to a Multiple Well System in a Confined Aquifer

The calculations of unsteady flow to a multiple well system with the application of boundary elementmethod (BEM) are discussed. The mathematical model of unsteady well flow is a boundary value problem ofparabolic differential equation. It is changed into an elliptic one by Laplace transform to eliminate time varia-ble. The image function of water head H can be solved by BEM. We derived the boundary integral equation ofthe transformed variable H and the discretization form of it, so that there is no need to discretize the bounda-ries of well walls and it becomes easier to solve the groundwater head H by numerical inversion.

A simplified two-dimensional boundary element method with arbitrary uniform mean flow

To reduce computational costs, an improved form of the frequency domain boundary element method（BEM） is proposed for two-dimensional radiation and propagation acoustic problems in a subsonic uniform flow with arbitrary orientation. The boundary integral equation（BIE） representation solves the two-dimensional convected Helmholtz equation（CHE） and its fundamental solution, which must satisfy a new Sommerfeld radiation condition（SRC） in the physical space. In order to facilitate conventional formulations, the variables of the advanced form are expressed only in terms of the acoustic pressure as well as its normal and tangential derivatives, and their multiplication operators are based on the convected Green＇s kernel and its modified derivative. The proposed approach significantly reduces the CPU times of classical computational codes for modeling acoustic domains with arbitrary mean flow. It is validated by a comparison with the analytical solutions for the sound radiation problems of monopole,dipole and quadrupole sources in the presence of a subsonic uniform flow with arbitrary orientation.

GPU-based discrete element simulation on flow stability of flat-bottomed hopper

In this study, the flow stability of the flat-bottomed hopper was investigated via GPU-based discrete element method(DEM) simulation. With the material height inside the hopper reducing, the fluctuation of the flow rate indicates an unstable discharge. The flow regions of the unstable discharge were compared with that of the stable discharge, a key transformation zone, where the voidage showed the largest difference between unstable and stable discharge, was revealed. To identify the relevance of the key transformation zone and the hopper flow stability, the voidage variation of the key transformation zone with material height reducing was studied.A sharp increase in the voidage in the key transformation zone was considered to be the standard for judging the unstable hopper flow, and the ‘Top–Bottom effect' of the hopper was defined, which indicated the hopper flow was unstable when the hopper only had the top area and the bottom area, because the voidage of particles in the top area and the bottom area were both variables.

FINITE ELEMENT ANALYSES OF THE EXTRUSION FLOWS OF BOGER FLUIDS WITH END-DRAWING AND GRAVITY-DRAWING

The spinning flow of Boger fluids and the gravity-drawing extrusion flow of a Newtonianas well as a Boger fluid have been simulated by using the stream-line finite element method and thetechnique of matching the finite element solutions with those of one-dimensional spinning equations.The recoverable shear strain is proved not to be a basic parameter in characterising thespinning flow of Boger fluids.For Newtonian fluids this technique predicts the experimental jetshape accurately.For Boger fluids,the numerical simulation agrees with the experimental data of spin-ning flow reported by Sridhar et al.,but seems to give an insufficient swelling and over contractionof the jets when drawn by its own weight,compared with the experimental results of Trang andYeow.It implies that the Oldroyd-B model fitting the viscometric-flow data fails to describeaccurately the elasticity and extensional viscosity in the extrusion flow of Boger fluids with gravi-ty-drawing.