A Petrov–Galerkin finite element model for one‐dimensional fully non‐linear and weakly dispersive wave propagation

A new finite element method is presented to solve one‐dimensional depth‐integrated equations for fully non‐linear and weakly dispersive waves. For spatial integration, the Petrov–Galerkin weighted residual method is used. The weak forms of the governing equations are arranged in such a way that the shape functions can be piecewise linear, while the weighting functions are piecewise cubic with C²‐continuity. For the time integration an implicit predictor–corrector iterative scheme is employed. Within the framework of linear theory, the accuracy of the scheme is discussed by considering the truncation error at a node. The leading truncation error is fourth‐order in terms of element size. Numerical stability of the scheme is also investigated. If the Courant number is less than 0.5, the scheme is unconditionally stable. By increasing the number of iterations and/or decreasing the element size, the stability characteristics are improved significantly. Both Dirichlet boundary condition (for incident waves) and Neumann boundary condition (for a reflecting wall) are implemented. Several examples are presented to demonstrate the range of applicabilities and the accuracy of the model. Copyright © 2001 John Wiley & Sons, Ltd.     

Development of new finite volume schemes on unstructured triangular grid for simulating the gas–liquid two‐phase flow

This paper presents a coupled finite volume inner doubly iterative efficient algorithm for linked equations (IDEAL) with level set method to simulate the incompressible gas–liquid two‐phase flows with moving interfaces on unstructured triangular grid. The finite volume IDEAL method on a collocated grid is employed to solve the incompressible two‐phase Navier–Stokes equations, and the level set method is used to capture the moving interfaces. For the sake of mass conservation, an effective second‐order accurate finite volume scheme is developed to solve the level set equation on triangular grid, which can be implemented much easier than the classical high‐order level set solvers. In this scheme, the value of level set function on the boundary of control volume is approximated using a linear combination of a high‐order Larangian interpolation and a second‐order upwind interpolation. By the rotating slotted disk and stretching and shrinking of a circular fluid element benchmark cases, the mass conservation and accuracy of the new scheme is verified. Then the coupled method is applied to two‐phase flows, including a 2D bubble rising problem and a 2D dam breaking problem. The computational results agree well with those reported in literatures and experimental data. Copyright © 2015 John Wiley & Sons, Ltd.     

Fluid–solid interaction problems with thermal convection using the immersed element‐free Galerkin method

In this work, the immersed element‐free Galerkin method (IEFGM) is proposed for the solution of fluid–structure interaction (FSI) problems. In this technique, the FSI is represented as a volumetric force in the momentum equations. In IEFGM, a Lagrangian solid domain moves on top of an Eulerian fluid domain that spans over the entire computational region. The fluid domain is modeled using the finite element method and the solid domain is modeled using the element‐free Galerkin method. The continuity between the solid and fluid domains is satisfied by means of a local approximation, in the vicinity of the solid domain, of the velocity field and the FSI force. Such an approximation is achieved using the moving least‐squares technique. The method was applied to simulate the motion of a deformable disk moving in a viscous fluid due to the action of the gravitational force and the thermal convection of the fluid. An analysis of the main factors affecting the shape and trajectory of the solid body is presented. The method shows a distinct advantage for simulating FSI problems with highly deformable solids. Copyright © 2009 John Wiley & Sons, Ltd.     

饱和多孔介质大变形分析的一种有限元-有限体积混合方法

建立了饱和多孔介质大变形分析的一种有限元-有限体积混合计算方法.将饱和多孔介质视为由固体骨架和孔隙水组成的两相体,其基本方程包括动力平衡方程和渗流连续方程.基于u-p假定和更新的Lagrange方法,饱和多孔介质的动力平衡方程在空间域内采用有限元方法进行离散,而渗流连续方程在空阃域内则采用有限体积法进行离散.通过两个数值算例,一维有限弹性固结和动力荷载作用下堤坝动力响应的计算,验证了该方法的有效性.     

A depth‐integrated non‐hydrostatic finite element model for wave propagation

This paper presents a two‐dimensional finite element model for simulating dynamic propagation of weakly dispersive waves. Shallow water equations including extra non‐hydrostatic pressure terms and a depth‐integrated vertical momentum equation are solved with linear distributions assumed in the vertical direction for the non‐hydrostatic pressure and the vertical velocity. The model is developed based on the platform of a finite element model, CCHE2D. A physically bounded upwind scheme for the advection term discretization is developed, and the quasi second‐order differential operators of this scheme result in no oscillation and little numerical diffusion. The depth‐integrated non‐hydrostatic wave model is solved semi‐implicitly: the provisional flow velocity is first implicitly solved using the shallow water equations; the non‐hydrostatic pressure, which is implicitly obtained by ensuring a divergence‐free velocity field, is used to correct the provisional velocity, and finally the depth‐integrated continuity equation is explicitly solved to satisfy global mass conservation. The developed wave model is verified by an analytical solution and validated by laboratory experiments, and the computed results show that the wave model can properly handle linear and nonlinear dispersive waves, wave shoaling, diffraction, refraction and focusing. Copyright © 2013 John Wiley & Sons, Ltd.     

A high‐order Petrov–Galerkin finite element method for the classical Boussinesq wave model

A high‐order Petrov–Galerkin finite element scheme is presented to solve the one‐dimensional depth‐integrated classical Boussinesq equations for weakly non‐linear and weakly dispersive waves. Finite elements are used both in the space and the time domains. The shape functions are bilinear in space–time, whereas the weighting functions are linear in space and quadratic in time, with C⁰‐continuity. Dispersion correction and a highly selective dissipation mechanism are introduced through additional streamline upwind terms in the weighting functions. An implicit, conditionally stable, one‐step predictor–corrector time integration scheme results. The accuracy and stability of the non‐linear discrete equations are investigated by means of a local Taylor series expansion. A linear spectral analysis is used for the full characterization of the predictor–corrector inner iterations. Based on the order of the analytical terms of the Boussinesq model and on the order of the numerical discretization, it is concluded that the scheme is fourth‐order accurate in terms of phase velocity. The dissipation term is third order only affecting the shortest wavelengths. A numerical convergence analysis showed a second‐order convergence rate in terms of both element size and time step. Four numerical experiments are addressed and their results are compared with analytical solutions or experimental data available in the literature: the propagation of a solitary wave, the oscillation of a flat bottom closed basin, the oscillation of a non‐flat bottom closed basin, and the propagation of a periodic wave over a submerged bar. Copyright © 2008 John Wiley & Sons, Ltd.     

Similarities between the quasi‐bubble and the generalized wave continuity equation solutions to the shallow water equations

Two common strategies for solving the shallow water equations in the finite element community are the generalized wave continuity equation (GWCE) reformulation and the quasi‐bubble velocity approximation. The GWCE approach has been widely analysed in the literature. In this work, the quasi‐bubble equations are analysed and comparisons are made between the quasi‐bubble approximation of the primitive form of the shallow water equations and a linear finite element approximation of the GWCE reformulation of the shallow water equations. The discrete condensed quasi‐bubble continuity equation is shown to be identical to a discrete wave equation for a specific GWCE weighting parameter value. The discrete momentum equations are slightly different due to the bubble function. In addition, the dispersion relationships are shown to be almost identical and numerical experiments confirm that the two schemes compute almost identical results. Analysis of the quasi‐bubble formulation suggests a relationship that may guide selection of the optimal GWCE weighting parameter. Copyright © 2004 John Wiley & Sons, Ltd.     

Radial integral boundary element method for simulating phase change problem with mushy zone

A radial integral boundary element method(BEM) is used to simulate the phase change problem with a mushy zone in this paper. Three phases, including the solid phase, the liquid phase, and the mushy zone, are considered in the phase change problem. First, according to the continuity conditions of temperature and its gradient on the liquid-mushy interface, the mushy zone and the liquid phase in the simulation can be considered as a whole part, namely, the non-solid phase, and the change of latent heat is approximated by heat source which is dependent on temperature. Then, the precise integration BEM is used to obtain the differential equations in the solid phase zone and the non-solid phase zone, respectively. Moreover, an iterative predictor-corrector precise integration method(PIM) is needed to solve the differential equations and obtain the temperature field and the heat flux on the boundary. According to an energy balance equation and the velocity of the interface between the solid phase and the mushy zone, the front-tracking method is used to track the move of the interface. The interface between the liquid phase and the mushy zone is obtained by interpolation of the temperature field.Finally, four numerical examples are provided to assess the performance of the proposed numerical method.     

Estimation of particle dynamics in 2-D fluidized beds using particle tracking velocimetry

The experimental characterization of particle dynamics in fluidized beds is of great importance in fostering an understanding of solid phase motion and its effect on particle properties in granulation processes. Commonly used techniques such as particle image velocimetry rely on the cross-correlation of illumination intensity and averaging procedures. It is not possible to obtain single particle velocities with such techniques. Moreover, the estimated velocities may not accurately represent the local particle velocities in regions with high velocity gradients. Consequently, there is a need for devices and methods that are capable of acquiring individual particle velocities. This paper describes how particle tracking velocimetry can be adapted to dense particulate flows. The approach presented in this paper couples high-speed imaging with an innovative segmentation algorithm for particle detection, and employs the Voronoi method to solve the assignment problem usually encountered in densely seeded flows. Lagrangian particle tracks are obtained as primary information, and these serve as the basis for calculating sophisticated quantities such as the solid-phase flow field, granular temperature, and solid volume fraction. We show that the consistency of individual trajectories is sufficient to recognize collision events.     

Numerical investigation of the flow front behaviour in the co‐injection moulding process

This paper presents a three‐dimensional (3D) solution algorithm for solving the sequential co‐injection moulding process. The flow of skin and core materials inside a rectangular cavity is investigated both numerically and experimentally. A 3D finite element flow analysis code is used to solve the governing equations of the non‐isothermal sequential co‐injection moulding. The predicted flow front behaviour is compared to the experimental observations for various skin/core volume ratio, injection speed, injection temperature, and core injection delay. Simulation results are in good agreement with experimental data and indicate correctly the trends in solution change when processing parameters are changing. Solutions are also shown for the filling of a spiral‐flow mould. The numerical approach is shown to predict the core expansion phase during which the flow front of core and skin materials advance together without breakthrough. Breakthrough phenomena is also predicted and the numerical solution is in good agreement with the experiment. Copyright © 2005 Crown in the right of Canada. Published by John Wiley & Sons, Ltd.     

Estimation of particle dynamics in 2-D fluidized beds using particle tracking velocimetry

The experimental characterization of particle dynamics in fluidized beds is of great importance in fostering an understanding of solid phase motion and its effect on particle properties in granulation processes.Commonly used techniques such as particle image velocimetry rely on the cross-correlation of illumination intensity and averaging procedures.It is not possible to obtain single particle velocities with such techniques.Moreover,the estimated velocities may not accurately represent the local particle velocities in regions with high velocity gradients.Consequently,there is a need for devices and methods that are capable of acquiring individual particle velocities.This paper describes how particle tracking velocimetry can be adapted to dense particulate flows.The approach presented in this paper couples high-speed imaging with an innovative segmentation algorithm for particle detection,and employs the Voronoi method to solve the assignment problem usually encountered in densely seeded flows.Lagrangian particle tracks are obtained as primary information,and these serve as the basis for calculating sophisticated quantities such as the solid-phase flow field,granular temperature,and solid volume fraction.We show that the consistency of individual trajectories is sufficient to recognize collision events.     

The improved space–time conservation element and solution element scheme for two‐dimensional dam‐break flow simulation

A new numerical scheme, namely space–time conservation element and solution element (CE/SE) method, has been used for the solution of the two‐dimensional (2D) dam‐break problem. Distinguishing from the well‐established traditional numerical methods (such as characteristics, finite difference, finite element, and finite‐volume methods), the CE/SE scheme has many non‐traditional features in both concept and methodology: space and time are treated in a unified way, which is the most important characteristic for the CE/SE method; the CEs and SEs are introduced, both local and global flux conservations in space and time rather than space only are enforced; an explicit scheme with a stagger grid is adopted. Furthermore, this scheme is robust and easy to implement. In this paper, an improved CE/SE scheme is extended to solve the 2D shallow water equations with the source terms, which usually plays a critical role in dam‐break flows. To demonstrate the accuracy, robustness and efficiency of the improved CE/SE method, both 1D and 2D dam‐break problems are simulated numerically, and the results are consistent with either the analytical solutions or experimental results. Copyright © 2011 John Wiley & Sons, Ltd.     

A coupled surface‐subsurface model for hydrostatic flows under saturated and variably saturated conditions

In this paper, the governing differential equations for hydrostatic surface‐subsurface flows are derived from the Richards and from the Navier‐Stokes equations. A vertically integrated continuity equation is formulated to account for both surface and subsurface flows under saturated and variable saturated conditions. Numerically, the horizontal domain is covered by an unstructured orthogonal grid that may include subgrid specifications. Along the vertical direction, a simple z‐layer discretization is adopted. Semi‐implicit finite difference equations for velocities, and a finite volume approximation for the vertically integrated continuity equation, are derived in such a fashion that, after simple manipulation, the resulting discrete pressure equation can be assembled into a single, two‐dimensional, mildly nonlinear system. This system is solved by a nested Newton‐type method, which yields simultaneously the (hydrostatic) pressure and a nonnegative fluid volume throughout the computational grid. The resulting algorithm is relatively simple, extremely efficient, and very accurate. Stability, convergence, and exact mass conservation are assured throughout also in presence of wetting and drying, in variable saturated conditions, and during flow transition through the soil interface. A few examples illustrate the model applicability and demonstrate the effectiveness of the proposed algorithm.     

Semi-empirical prediction of bubble diameter in gas fluidized beds

Theoretical expressions for bubble diameter in both small and large particle fluidized beds are derived by the application of two phase theory and gas flow continuity. Comparison with experimental data suggests that the numerical and analytical solution of these expressions, combined with empirical bubble frequency relations, can provide an accurate prediction of bubble size and its parametric trends.Several commonly employed empirical correlations of bubble diameter are shown to be derivable from a common theory, with differences among the correlations ascribed to variations in flow regime and bubble frequency.     

Numerical simulations of high pressure carbon dioxide fluid fluidized beds

Hydrodynamics of carbon dioxide fluid-particle mixtures are predicted using a low density ratio-based kinetic theory of granular flow in high pressure carbon dioxide fluid fluidized beds. A coexistence of particle waves and particle aggregates exists along bed height. The threshold to identify the occurrence of particle aggregates is suggested based on standard deviation of solid volume fractions in aggregative fluidization. The existence time fractions and frequencies of particle aggregates are predicted along axial direction. The effect of carbon dioxide fluid temperature and pressure on volume fraction and velocity distributions are analyzed at different inlet carbon dioxide velocities and particle densities in high pressure carbon dioxide fluidized beds. Simulated results indicate that the carbon dioxide-particles fluidization transits from particulate to aggregative states with the increase of inlet carbon dioxide velocities. The computed fluid volume fractions and heterogeneity indexes are close to the measurements in a high pressure carbon dioxide fluidized bed.     

Lagrangian–Eulerian simulation of slugging fluidized bed

This work studies gas–solid slugging fluidized beds with Type-D particles, using two-dimensional simulations based on discrete element model (DEM). DEM performance is quantitatively validated by two commonly accepted correlations for determining slugging behavior. The voidage profiles simulated with bed height corresponding to Baeyens and Geldart (1974) correlation for onset of slugging demonstrate a transitional flow pattern from free bubbling to slugging. The present calculated values for the maximum slugging bed height are in good agreement with the correlation from Matsen et al. (1969). Simulations show that fluidized beds with Type-D particles can operate in the round-nosed slugging regime and also shows that wall slugs and square-nosed slugs tend to be formed with increase in superficial gas velocity and in bed height, respectively.     

A parallel monolithic algorithm for the numerical simulation of large‐scale fluid structure interaction problems

A novel parallel monolithic algorithm has been developed for the numerical simulation of large‐scale fluid structure interaction problems. The governing incompressible Navier–Stokes equations for the fluid domain are discretized using the arbitrary Lagrangian–Eulerian formulation‐based side‐centered unstructured finite volume method. The deformation of the solid domain is governed by the constitutive laws for the nonlinear Saint Venant–Kirchhoff material, and the classical Galerkin finite element method is used to discretize the governing equations in a Lagrangian frame. A special attention is given to construct an algorithm with exact total fluid volume conservation while obeying both the global and the local discrete geometric conservation law. The resulting large‐scale algebraic nonlinear equations are multiplied with an upper triangular right preconditioner that results in a scaled discrete Laplacian instead of a zero block in the original system. Then, a one‐level restricted additive Schwarz preconditioner with a block‐incomplete factorization within each partitioned sub‐domains is utilized for the modified system. The accuracy and performance of the proposed algorithm are verified for the several benchmark problems including a pressure pulse in a flexible circular tube, a flag interacting with an incompressible viscous flow, and so on. John Wiley & Sons, Ltd.     

A C2‐continuous high‐resolution upwind convection scheme

A bounded upwinding scheme for numerical solution of hyperbolic conservation laws and Navier–Stokes equations is presented. The scheme is based on convection boundedness criterion and total variation diminishing stability criteria and developed by employing continuously differentiable functions. The accuracy of the scheme is verified by assessing the error and observed convergence rate on 1‐D benchmark test cases. A comparative study between the new scheme and conventional total variation diminishing/convection boundedness criterion‐based upwind schemes to solve standard nonlinear hyperbolic conservation laws is also accomplished. The scheme is then examined in the simulation of Newtonian and non‐Newtonian fluid flows of increasing complexity; a satisfactory agreement has been observed in terms of the overall behavior. Finally, the scheme is used to study the hydrodynamics of a gas‐solid flow in a bubbling fluidized bed. Copyright © 2013 John Wiley & Sons, Ltd.