A boundary element method for small-amplitude viscous fluid sloshing in couple with structural vibration

A boundary element method is presented for the coupled motion analysis of structural vibration with small-amplitude fluid sloshing in two-dimensional space. The linearized Navier-Stokes equations are considered in frequency domain and transformed into boundary integral equations. An appropriate fundamental solution for the Helmholtz equation with pure imaginary constant is found. The condition of zero-stress is imposed on the free surface, and non-slip condition of fluid particles is imposed on the walls of the container. For rigid motion models, the expressions for added mass and added damping to the structural motion equations are obtained. Some typical numerical examples are presented.     

Finite element study of vortex‐induced cross‐flow and in‐line oscillations of a circular cylinder at low Reynolds numbers

Vortex‐induced vibrations of a circular cylinder placed in a uniform flow at Reynolds number 325 are investigated using a stabilized space–time finite element formulation. The Navier–Stokes equations for incompressible fluid flow are solved for a two‐dimensional case along with the equations of motion of the cylinder that is mounted on lightly damped spring supports. The cylinder is allowed to vibrate, both in the in‐line and in the cross‐flow directions. Results of the computations are presented for various values of the structural frequency of the oscillator, including those that are sub and superharmonics of the vortex‐shedding frequency for a stationary cylinder. In most of the cases, the trajectory of the cylinder corresponds to a Lissajou figure of 8. Lock‐in is observed for a range of values of the structural frequency. Over a certain range of structural frequency (F_s), the vortex‐shedding frequency of the oscillating cylinder does not match F_s exactly; there is a slight detuning. This phenomenon is referred to as soft‐lock‐in. Computations show that this detuning disappears when the mass of the cylinder is significantly larger than the mass of the surrounding fluid it displaces. A self‐limiting nature of the oscillator with respect to cross‐flow vibration amplitude is observed. It is believed that the detuning of the vortex‐shedding frequency from the structural frequency is a mechanism of the oscillator to self‐limit its vibration amplitude. The dependence of the unsteady solution on the spatial resolution of the finite element mesh is also investigated. Copyright © 1999 John Wiley & Sons, Ltd.     

A hybrid immersed boundary and material point method for simulating 3D fluid–structure interaction problems

A numerical method is developed for solving the 3D, unsteady, incompressible Navier–Stokes equations in curvilinear coordinates containing immersed boundaries (IBs) of arbitrary geometrical complexity moving and deforming under forces acting on the body. Since simulations of flow in complex geometries with deformable surfaces require special treatment, the present approach combines a hybrid immersed boundary method (HIBM) for handling complex moving boundaries and a material point method (MPM) for resolving structural stresses and movement. This combined HIBM & MPM approach is presented as an effective approach for solving fluid–structure interaction (FSI) problems. In the HIBM, a curvilinear grid is defined and the variable values at grid points adjacent to a boundary are forced or interpolated to satisfy the boundary conditions. The MPM is used for solving the equations of solid structure and communicates with the fluid through appropriate interface‐boundary conditions. The governing flow equations are discretized on a non‐staggered grid layout using second‐order accurate finite‐difference formulas. The discrete equations are integrated in time via a second‐order accurate dual time stepping, artificial compressibility scheme. Unstructured, triangular meshes are employed to discretize the complex surface of the IBs. The nodes of the surface mesh constitute a set of Lagrangian control points used for tracking the motion of the flexible body. The equations of the solid body are integrated in time via the MPM. At every instant in time, the influence of the body on the flow is accounted for by applying boundary conditions at stationary curvilinear grid nodes located in the exterior but in the immediate vicinity of the body by reconstructing the solution along the local normal to the body surface. The influence of the fluid on the body is defined through pressure and shear stresses acting on the surface of the body. The HIBM & MPM approach is validated for FSI problems by solving for a falling rigid and flexible sphere in a fluid‐filled channel. The behavior of a capsule in a shear flow was also examined. Agreement with the published results is excellent. Copyright © 2007 John Wiley & Sons, Ltd.     

A fully coupled Navier–Stokes solver for calculation of turbulent incompressible free surface flow past a ship hull

This paper deals with the calculation of free surface flow of viscous incompressible fluid around the hull of a boat moving with rectilinear motion. An original method used to avoid a large part of the theoretical problems connected with free surface boundary conditions in three‐dimensional Navier–Stokes–Reynolds equations is proposed here. The linearised system of convective equations for velocities, pressure and free surface elevation unknowns is discretised by finite differences and two methods to solve the fully coupled resulting matrix are presented here. The non‐linear convergence of fully coupled algorithm is compared with the velocity–pressure weakly coupled algorithm SIMPLER. Turbulence is taken into account through Reynolds decomposition and k–ε or k–ω model to close the equations. These two models are implemented without wall function and numerical calculations are performed up to the viscous sub‐layer. Numerical results and comparisons with experiments are presented on the Series 60 CB=0.60 ship model for a Reynolds number Rn=4.5×10⁶ and a Froude number Fn=0.316. Copyright © 1999 John Wiley & Sons, Ltd.     

RETRACTED ARTICLE: Flow of fractional Maxwell fluid between coaxial cylinders

This paper deals with the study of unsteady flow of a Maxwell fluid with fractional derivative model, between two infinite coaxial circular cylinders, using Laplace and finite Hankel transforms. The motion of the fluid is produced by the inner cylinder that, at time t = 0⁺, is subject to a time-dependent longitudinal shear stress. Velocity field and the adequate shear stress are presented under series form in terms of the generalized G and R functions. The solutions that have been obtained satisfy all imposed initial and boundary conditions. The corresponding solutions for ordinary Maxwell and Newtonian fluids are obtained as limiting cases of general solutions. Finally, the influence of the pertinent parameters on the fluid motion as well as a comparison between the three models is underlined by graphical illustrations.     

Exact analytical solutions for axial flow of a fractional second grade fluid between two coaxial cylinders

The velocity field and the adequate shear stress corresponding to the longitudinal flow of a fractional second grade fluid, between two infinite coaxial circular cylinders, are determined by applying the Laplace and finite Hankel transforms. Initially the fluid is at rest, and at time t = 0⁺, the inner cylinder suddenly begins to translate along the common axis with constant acceleration. The solutions that have been obtained are presented in terms of generalized G functions. Moreover, these solutions satisfy both the governing differential equations and all imposed initial and boundary conditions. The corresponding solutions for ordinary second grade and Newtonian fluids are obtained as limiting cases of the general solutions. Finally, some characteristics of the motion, as well as the influences of the material and fractional parameters on the fluid motion and a comparison between models, are underlined by graphical illustrations.     

A spectral‐element discontinuous Galerkin thermal lattice Boltzmann method for conjugate heat transfer applications

We present a spectral‐element discontinuous Galerkin thermal lattice Boltzmann method for fluid–solid conjugate heat transfer applications. Using the discrete Boltzmann equation, we propose a numerical scheme for conjugate heat transfer applications on unstructured, non‐uniform grids. We employ a double‐distribution thermal lattice Boltzmann model to resolve flows with variable Prandtl (Pr) number. Based upon its finite element heritage, the spectral‐element discontinuous Galerkin discretization provides an effective means to model and investigate thermal transport in applications with complex geometries. Our solutions are represented by the tensor product basis of the one‐dimensional Legendre–Lagrange interpolation polynomials. A high‐order discretization is employed on body‐conforming hexahedral elements with Gauss–Lobatto–Legendre quadrature nodes. Thermal and hydrodynamic bounce‐back boundary conditions are imposed via the numerical flux formulation that arises because of the discontinuous Galerkin approach. As a result, our scheme does not require tedious extrapolation at the boundaries, which may cause loss of mass conservation. We compare solutions of the proposed scheme with an analytical solution for a solid–solid conjugate heat transfer problem in a 2D annulus and illustrate the capture of temperature continuities across interfaces for conductivity ratio γ > 1. We also investigate the effect of Reynolds (Re) and Grashof (Gr) number on the conjugate heat transfer between a heat‐generating solid and a surrounding fluid. Steady‐state results are presented for Re = 5?40 and Gr = 10⁵?10⁶. In each case, we discuss the effect of Re and Gr on the heat flux (i.e. Nusselt number Nu) at the fluid–solid interface. Our results are validated against previous studies that employ finite‐difference and continuous spectral‐element methods to solve the Navier–Stokes equations. Copyright © 2016 John Wiley & Sons, Ltd.     

Stability analysis of pipes conveying fluid with fractional viscoelastic model

Divergence and flutter instabilities of pipes conveying fluid with fractional viscoelastic model has been investigated in the present work. Attention is concentrated on the boundaries of the stability. Based on the Euler–Bernoulli beam theory for structural dynamics, viscoelastic fractional model for damping and, plug flow model for fluid flow, equation of motion has been derived. The effects of gravity, and distributed follower forces are also considered. By transferring the equation of motion to the Laplace domain and using the Galerkin method, the characteristic equations are obtained. By solving the eigenvalue problem, frequencies and dampings of the system have been obtained versus flow velocity. Some numerical test cases have been studied with viscoelastic fractional model and the effect of the fractional derivative order and the retardation time is investigated for various boundary conditions.

     

Hamilton principle for an inviscid compressible fluid model in Eulerian coordinates

It is shown that the well-known variational principles for the ideal compressible fluid model in Eulerian coordinates have the following deficiencies:

1. They are not related to the corresponding variational principles in Lagrangian coordinates;
2. The variation procedure in these variational problems does not lead to the equations of motion themselves in the Euler form; rather it leads to relations which correspond to definite classes of solutions of the Euler equations. Here allowance for the equations of the constraints imposed by the adiabaticity and continuity conditions limits the region of application of these variational principles to only potential flows;
3. More general results, involving flows other than potential, are achieved by artificial selection of certain additional constraint conditions imposed on the quantities being varied, and in this case additional clarification is required to ascertain whether any inviscid compressible fluid flow is the extremum of the corresponding variational problem.

A new formulation of the Hamilton principle for the inviscid compressible fluid in Eulerian coordinates is suggested which is free from these deficiencies.     

A numerical method to solve the m‐terms of a submerged body with forward speed

To model mathematically the problem of a rigid body moving below the free surface, a control surface surrounding the body is introduced. The linear free surface condition of the steady waves created by the moving body is satisfied. To describe the fluid flow outside this surface a potential integral equation is constructed using the Kelvin wave Green function whereas inside the surface, a source integral equation is developed adopting a simple Green function. Source strengths are determined by matching the two integral equations through continuity conditions applied to velocity potential and its normal derivatives along the control surface. After solving for the induced fluid velocity on the body surface and the control surface, an integral equation is derived involving a mixed distribution of sources and dipoles using a simple Green function and one component of the fluid velocity. The normal derivatives of the fluid velocity on the body surface, namely the m‐terms, are then solved by this matching integral equation method (MIEM). Numerical results are presented for two elliptical sections moving at a prescribed Froude number and submerged depth and a sensitivity analysis undertaken to assess the influence of these parameters. Furthermore, comparisons are performed to analyse the impact of different assumptions adopted in the derivation of the m‐terms. It is found that the present method is easy to use in a panel method with satisfactory numerical precision. Copyright © 2002 John Wiley & Sons, Ltd.     

A control-volume based finite element method for three-dimensional incompressible turbulent fluid flow,heat transfer,and related phenomena

A control-volume based finite element method of equal-order type for three-dimensional incompressible turbulent fluid flow, heat transfer, and related phenomena is presented. The discretization equations are based mainly on the physics of the phenomena under consideration, more than on mathematical arguments. Special emphasis is devoted to the discretization of the convective terms and the continuity equation, and to the treatment of the boundary conditions imposed by the use of a high Reynolds k-?, type turbulence model. The pressure-velocity coupling in the fluid flow calculation is made from a derivative of the original SIMPLER method, without pressure correction. The discretized equations are solved in a sequential, rather than a coupled, form with significant advantage in the required computer time and storage. The method is an extension of a former version proposed by us for two-dimensional, laminar problems, and is here successfully applied to the following situations: three-dimensional deflected turbulent jet, and flows in 90° and 45° junctions of ducts with rectangular cross sections. The calculated results are in very good agreement with the experimental and numerical (obtained with the well established finite difference method) data available in the literature.     

On the unsteady rotational flow of fractional Oldroyd-B fluid in cylindrical domains

This paper concerned with the unsteady rotational flow of fractional Oldroyd-B fluid, between two infinite coaxial circular cylinders. To solve the problem we used the finite Hankel and Laplace transforms. The motion is produced by the inner cylinder that, at time t=0⁺, is subject to a time-dependent rotational shear. The solutions that have been obtained, presented under series form in terms of the generalized G functions, satisfy all imposed initial and boundary conditions. The corresponding solutions for ordinary Oldroyd-B, fractional and ordinary Maxwell, fractional and ordinary second grade, and Newtonian fluids, performing the same motion, are obtained as limiting cases of general solutions.     

Exact solutions for the flow of second grade fluid in annulus between torsionally oscillating cylinders

The velocity field and the associated shear stress corresponding to the torsional oscillatory flow of a second grade fluid, between two infinite coaxial circular cylinders, are determined by means of the Laplace and Hankel transforms. At time t = 0, the fluid and both the cylinders are at rest and at t = 0 + , cylinders suddenly begin to oscillate around their common axis in a simple harmonic way having angular frequencies ω 1 and ω 2 . The obtained solutions satisfy the governing differential equation and all imposed initial and boundary conditions. The solutions for the motion between the cylinders, when one of them is at rest, can be obtained from our general solutions. Furthermore, the corresponding solutions for Newtonian fluid are also obtained as limiting cases of our general solutions.     

Instabilities during batch sedimentation in geometries containing obstacles: A numerical and experimental study

Batch sedimentation of non‐colloidal particle suspensions is studied with nuclear magnetic resonance flow visualization and continuum‐level numerical modelling of particle migration. The experimental method gives particle volume fraction as a function of time and position, which then provides validation data for the numerical model. A finite element method is used to discretize the equations of motion, including an evolution equation for the particle volume fraction and a generalized Newtonian viscosity dependent on local particle concentration. The diffusive‐flux equation is based on the Phillips model (Phys. Fluids A 1992; 4 :30–40) and includes sedimentation terms described by Zhang and Acrivos (Int. J. Multiphase Flow 1994; 20 :579–591). The model and experiments are utilized in three distinct geometries with particles that are heavier and lighter than the suspending fluid, depending on the experiment: (1) sedimentation in a cylinder with a contraction; (2) particle flotation in a horizontal cylinder with a horizontal rod; and (3) flotation around a rectangular inclusion. Secondary flows appear in both the experiments and the simulations when a region of higher density fluid is above a lower density fluid. The secondary flows result in particle inhomogeneities, Rayleigh–Taylor‐like instabilities, and remixing, though the effect in the simulations is more pronounced than in the experiments. Published in 2007 by John Wiley & Sons, Ltd.     

Modelling a recirculating density-driven turbulent flow

A mathematical model of turbulent density-driven flows is presented and is solved numerically. A form of the k–? turbulence model is used to characterize the turbulent transport, and both this non-linear model and a sediment transport equation are coupled with the mean-flow fluid motion equations. A partitioned, Newton–Raphson-based solution scheme is used to effect a solution. The model is applied to the study of flow through a circular secondary sedimentation basin.     

A local‐analytic‐based discretization procedure for the numerical solution of incompressible flows

An Erratum has been published for this article in International Journal for Numerical Methods in Fluids 2005, 49(8): 933. We present a local‐analytic‐based discretization procedure for the numerical solution of viscous fluid flows governed by the incompressible Navier–Stokes equations. The general procedure consists of building local interpolants obtained from local analytic solutions of the linear multi‐dimensional advection–diffusion equation, prototypical of the linearized momentum equations. In view of the local analytic behaviour, the resulting computational stencil and coefficient values are functions of the local flow conditions. The velocity–pressure coupling is achieved by a discrete projection method. Numerical examples in the form of well‐established verification and validation benchmarks are presented to demonstrate the capabilities of the formulation. The discretization procedure is implemented alongside the ability to treat embedded and non‐matching grids with relative motion. Of interest are flows at high Reynolds number, ??(10⁵)–??(10⁷), for which the formulation is found to be robust. Applications include flow past a circular cylinder undergoing vortex‐induced vibrations (VIV) at high Reynolds number. Copyright © 2005 John Wiley & Sons, Ltd.     

Simulating 3D deformable particle suspensions using lattice Boltzmann method with discrete external boundary force

A method for direct numerical analysis of three‐dimensional deformable particles suspended in fluid is presented. The flow is computed on a fixed regular ‘lattice’ using the lattice Boltzmann method (LBM), where each solid particle is mapped onto a Lagrangian frame moving continuously through the domain. Instead of the bounce‐back method, an external boundary force (EBF) is used to impose the no‐slip boundary condition at the fluid–solid interface for stationary or moving boundaries. The EBF is added directly to the lattice Boltzmann equation. The motion and orientation of the particles are obtained from Newtonian dynamics equations. The advantage of this approach is outlined in comparison with the standard and higher‐order interpolated bounce‐back methods as well as the LBM immersed‐boundary and the volume‐of‐fluid methods. Although the EBF method is general, in this application, it is used in conjunction with the lattice–spring model for deformable particles. The methodology is validated by comparing with experimental and theoretical results. Copyright © 2009 John Wiley & Sons, Ltd.     

Analysis of non-axisymmetric wave propagation in a homogeneous piezoelectric solid circular cylinder of transversely isotropic material

A study concerning the propagation of free non-axisymmetric waves in a homogeneous piezoelectric cylinder of transversely isotropic material with axial polarization is carried out on the basis of the linear theory of elasticity and linear electro-mechanical coupling. The solution of the three dimensional equations of motion and quasi-electrostatic equation is given in terms of seven mechanical and three electric potentials. The characteristic equations are obtained by the application of the mechanical and two types of electric boundary conditions at the surface of the piezoelectric cylinder. A novel method of displaying dispersion curves is described in the paper and the resulting dispersion curves are presented for propagating and evanescent waves for PZT-4 and PZT-7A piezoelectric ceramics for circumferential wave numbers m = 1, 2, and 3. It is observed that the dispersion curves are sensitive to the type of the imposed boundary conditions as well as to the measure of the electro-mechanical coupling of the material.