International Journal for Numerical Methods in Fluids

Hauke, Guillermo; Doweidar, Mohamed H.; Fuentes, Sergio

doi: 10.1002/fld.2666pmid: N/A

Recently, we developed an explicit a posteriori error estimator especially suited for fluid dynamics problems solved with a stabilized method. The technology is based upon the theory that inspired stabilized methods, namely, the variational multiscale theory. The salient features of the formulation are that it can be readily implemented in existing codes, it is a very economical procedure, and it yields very accurate local error estimates uniformly from the diffusive to the advective regime.

Plumerault, L.‐R.; Astruc, D.; Villedieu, P.; Maron, P.

doi: 10.1002/fld.2667pmid: N/A

This work presents a numerical model designed for the simulation of water‐wave impacts on a structure when aeration of the liquid phase is considered. The model is based on a multifluid Navier–Stokes approach in which all fluids are assumed compressible. The numerical method is based on a finite volume algorithm in space and a second order Runge–Kutta method in time. A validation of this model is performed. It shows a good accuracy for acoustic and shock wave propagation in a bubbly liquid and for wave breaking. Copyright © 2011 John Wiley & Sons, Ltd.

Storti, Mario A.; Garelli, Luciano; Paz, Rodrigo R.

doi: 10.1002/fld.2669pmid: N/A

In this article, a new methodology for developing discrete geometric conservation law (DGCL) compliant formulations is presented. It is carried out in the context of the finite element method for general advective–diffusive systems on moving domains using an ALE scheme. There is an extensive literature about the impact of DGCL compliance on the stability and precision of time integration methods. In those articles, it has been proved that satisfying the DGCL is a necessary and sufficient condition for any ALE scheme to maintain on moving grids the nonlinear stability properties of its fixed‐grid counterpart. However, only a few works proposed a methodology for obtaining a compliant scheme. In this work, a DGCL compliant scheme based on an averaged ALE Jacobians formulation is obtained. This new formulation is applied to the θ family of time integration methods. In addition, an extension to the three‐point backward difference formula is given. With the aim to validate the averaged ALE Jacobians formulation, a set of numerical tests are performed. These tests include 2D and 3D diffusion problems with different mesh movements and the 2D compressible Navier–Stokes equations. Copyright © 2011 John Wiley & Sons, Ltd.

Zhou, C.H.; Shu, C.

doi: 10.1002/fld.2670pmid: N/A

The local domain‐free discretization method is extended in this work to simulate fluid–structure interaction problems, the class of which is exemplified by the self‐propelled anguilliform swimming of deforming bodies in a fluid medium. Given the deformation of the fish body in its own reference frame, the translational and rotational motions of the body governed by Newton's Law are solved together with the surrounding flow field governed by Navier–Stokes equations. When the body is deforming and moving, no mesh regeneration is required in the computation. The loose coupling strategy is employed to simulate the fluid–structure interaction involved in the self‐propelled swimming. The local domain‐free discretization method and an efficient algorithm for classifying the Eulerian mesh points are described in brief. To validate the fluid–structure interaction solver, we simulate the ‘lock‐in’ phenomena associated with the vortex‐induced vibrations of an elastically mounted cylinder. Finally, we demonstrate applications of the method to two‐dimensional and three‐dimensional anguilliform‐swimming fish. The kinematics and dynamics associated with the center of mass are shown and the rotational movement is also presented via the angular position of the body axis. The wake structure is visualized in terms of vorticity contours. All the obtained numerical results show good agreement with available data in the literature. Copyright © 2011 John Wiley & Sons, Ltd.

Khan, M.; Shaheen, N.; Shahzad, A.

doi: 10.1002/fld.2673pmid: N/A

In this paper, the steady flow and heat transfer of a magnetohydrodynamic fluid is studied. The fluid is assumed to be electrically conducting in the presence of a uniform magnetic field and occupies the porous space in annular pipe. The governing nonlinear equations are modeled by introducing the modified Darcy's law obeying the Sisko model. The system is solved using the homotopy analysis method (HAM), which yields analytical solutions in the form of a rapidly convergent infinite series. Also, HAM is used to obtain analytical solutions of the problem for noninteger values of the power index. The resulting problem for velocity field is then numerically solved using an iterative method to show the accuracy of the analytic solutions. The obtained solutions for the velocity and temperature fields are graphically sketched and the salient features of these solutions are discussed for various values of the power index parameter. We also present a comparison between Sisko and Newtonian fluids. Copyright © 2011 John Wiley & Sons, Ltd.

Choi, Doo Yong; Yuan, Hengliang

doi: 10.1002/fld.2676pmid: N/A

A horizontally curvilinear non‐hydrostatic free surface model that embeds the second‐order projection method, the so‐called θ scheme, in fractional time stepping is developed to simulate nonlinear wave motion in curved boundaries. The model solves the unsteady, Navier–Stokes equations in a three‐dimensional curvilinear domain by incorporating the kinematic free surface boundary condition with a top‐layer boundary condition, which has been developed to improve the numerical accuracy and efficiency of the non‐hydrostatic model in the standard staggered grid layout. The second‐order Adams–Bashforth scheme with the third‐order spatial upwind method is implemented in discretizing advection terms. Numerical accuracy in terms of nonlinear phase speed and amplitude is verified against the nonlinear Stokes wave theory over varying wave steepness in a two‐dimensional numerical wave tank. The model is then applied to investigate the nonlinear wave characteristics in the presence of dispersion caused by reflection and diffraction in a semicircular channel. The model results agree quantitatively with superimposed analytical solutions. Finally, the model is applied to simulate nonlinear wave run‐ups caused by wave‐body interaction around a bottom‐mounted cylinder. The numerical results exhibit good agreement with experimental data and the second‐order diffraction theory. Overall, it is shown that the developed model, with only three vertical layers, is capable of accurately simulating nonlinear waves interacting within curved boundaries. Copyright © 2011 John Wiley & Sons, Ltd.

Luo, X.L.; Gu, Z.L.; Lei, K.B.; Wang, S.; Kase, K.

doi: 10.1002/fld.2678pmid: N/A

A three‐dimensional Cartesion cut cell method is presented for the simulations of incompressible viscous flows with irregular domains. A new model (referred to as ‘6+N’ model) is proposed to describe arbitrarily shaped cut cells and treat all the cells as polyhedrons with 6+N faces. The finite volume discretization of the Navier–Stokes equation is then implemented by using the ‘6+N’ model to separate the surface flux integrals into two parts, that is, the fluxes through the basic face of the hexahedron and those through the cutting surfaces. The previously proposed Kitta Cube algorithm and volume computer‐aided design platform (J. Comput. Aided. Des. 2005; 37(4): 1509–1520. Doi:10.1016/j.cad.2005.03.006) are adopted to generate cut cells and provide shape data and physical attributes for the numerical analysis. A modified SIMPLE‐based smoothing pressure correction scheme is applied to suppress checkerboard pressure oscillations caused by the collocated arrangement of velocities and pressure. The calculation accuracy of the numerical method expressed by L1 and L ∞  norm errors is first demonstrated by the simulation of a pipe flow. Then its feasibility, efficiency, and potential in engineering applications are verified by applying it to solve natural convections between concentric spheres and between eccentric spheres. The heat transfer patterns in eccentric spheres are also obtained by using the numerical method. Copyright © 2011 John Wiley & Sons, Ltd.