Flow simulation on moving boundary‐fitted grids and application to fluid–structure interaction problems

We present a method for the parallel numerical simulation of transient three‐dimensional fluid–structure interaction problems. Here, we consider the interaction of incompressible flow in the fluid domain and linear elastic deformation in the solid domain. The coupled problem is tackled by an approach based on the classical alternating Schwarz method with non‐overlapping subdomains, the subproblems are solved alternatingly and the coupling conditions are realized via the exchange of boundary conditions. The elasticity problem is solved by a standard linear finite element method. A main issue is that the flow solver has to be able to handle time‐dependent domains. To this end, we present a technique to solve the incompressible Navier–Stokes equation in three‐dimensional domains with moving boundaries. This numerical method is a generalization of a finite volume discretization using curvilinear coordinates to time‐dependent coordinate transformations. It corresponds to a discretization of the arbitrary Lagrangian–Eulerian formulation of the Navier–Stokes equations. Here the grid velocity is treated in such a way that the so‐called Geometric Conservation Law is implicitly satisfied. Altogether, our approach results in a scheme which is an extension of the well‐known MAC‐method to a staggered mesh in moving boundary‐fitted coordinates which uses grid‐dependent velocity components as the primary variables. To validate our method, we present some numerical results which show that second‐order convergence in space is obtained on moving grids. Finally, we give the results of a fully coupled fluid–structure interaction problem. It turns out that already a simple explicit coupling with one iteration of the Schwarz method, i.e. one solution of the fluid problem and one solution of the elasticity problem per time step, yields a convergent, simple, yet efficient overall method for fluid–structure interaction problems. Copyright © 2005 John Wiley & Sons, Ltd.     

Residual‐based turbulence models for moving boundary flows: hierarchical application of variational multiscale method and three‐level scale separation

This paper presents residual‐based turbulence models for problems with moving boundaries and interfaces. The method is developed via a hierarchical application of variational multiscale ideas and the models are cast in an arbitrary Lagrangian–Eulerian (ALE) frame to accommodate the deformation of domain boundaries. An overlapping additive decomposition of velocity and pressure fields into coarse and fine scale components leads to coarse and fine scale mixed‐field problems. The problem governing fine scales is subjected to a further decomposition of the fine scale velocity into overlapping components termed as fine scales level I and level II. In turn, in the bottom‐up integration of scales, the model for level II fine scales serves to stabilize the problem governing level I fine scales, and model for level I fields yields the turbulence models. From the computational perspective, the coarse scales are represented in terms of the standard Lagrange shape functions, whereas level I and level II scales are represented via quadratic and fourth order polynomial bubbles, respectively. Because of the bubble functions approach employed in the consistently derived fine scale models, the resulting method is free of any embedded or tunable parameters. The proposed turbulence models share a common feature with the LES models in that the largest scales in the flow are numerically resolved, whereas the subgrid scales are modeled. The method is applied to flow around a plunging airfoil at Re = 40,000, and results are compared with experimental and numerical data published in the literature. Also presented are the results for the plunging airfoil at Re = 60,000 to show the robustness and range of applicability of the method. Copyright © 2013 John Wiley & Sons, Ltd.     

CFD development for macro particle simulations

Numerous industrial operations involve fluid-particle systems, in which both phases display very complex behaviour. Some examples include fluidisation technology in catalytic reactors, pneumatic transportation of grain or powder materials, carbon nanotube alignment in the nano-devices and circuit integration and so on. In this paper, a macro particle method is developed to model the fluid-particle flows. The macro particle is formed by a collection of micro-sized particles so that the number of macro particles to be tracked is much less than the number of smaller particles. Unlike the calculations of instantaneous point variables of fluid phase with moving discrete boundaries of the smaller particles with direct numerical simulation, the boundary of each macro particle is just dealt with the blocked-off approach. On the other hand, the flow fields based on the present method is solved by original Navier–Stokes, rather than the modified ones based on the locally averaged theorem. The flow fields are solved on the length scale of computational cells, while the resolutions of solid particles are the size of macro particle, which is determined as needed in specific applications. The macro particle method is validated by several selected cases, which demonstrate that the macro particle method could accurately resolve fluid-particle systems in an efficient, robust and flexible fashion.     

Improved MPS method with variable‐size particles

Using variable‐size particles in the moving particle semi‐implicit method (MPS) could lead to inaccurate predictions and/or numerical instability. In this paper, a variable‐size particle moving particle semi‐implicit method (VSP‐MPS) scheme is proposed for the MPS method to achieve more reliable simulations with variable‐size particles. To improve stability and accuracy, a new gradient model is developed based on a previously developed MPS scheme that requires no surface detection MPS. The dynamic particle coalescing and splitting algorithm is revised to achieve dynamic multi‐resolution. A cubic spline function with additional function is employed as the kernel function. The effectiveness of the VSP‐MPS method is demonstrated by three verification examples, that is, a hydrostatic pressure problem, a complicated free surface flow problem with large deformation, and a dynamic impact problem. The new VSP‐MPS scheme with variable‐size particles is found to have balanced efficiency and accuracy that is suitable for simulating large systems with complex flow patterns. Copyright © 2015 John Wiley & Sons, Ltd.     

Coupling methods between finite element–based Boussinesq‐type wave and particle‐based free‐surface flow models

We develop one‐way coupling methods between a Boussinesq‐type wave model based on the discontinuous Galerkin finite element method and a free‐surface flow model based on a mesh‐free particle method to strike a balance between accuracy and computational cost. In our proposed model, computation of the wave model in the global domain is conducted first, and the nonconstant velocity profiles in the vertical direction are reproduced by using its results. Computation of the free‐surface flow is performed in a local domain included within the global domain with interface boundaries that move along the reproduced velocity field in a Lagrangian fashion. To represent the moving interfaces, we used a polygon wall boundary model for mesh‐free particle methods. Verification and validation tests of our proposed model are performed, and results obtained by the model are compared with theoretical values and experimental results to show its accuracy and applicability.     

Coupled Navier–Stokes—Molecular dynamics simulations using a multi‐physics flow simulation framework

Simulation of nano‐scale channel flows using a coupled Navier–Stokes/Molecular Dynamics (MD) method is presented. The flow cases serve as examples of the application of a multi‐physics computational framework put forward in this work. The framework employs a set of (partially) overlapping sub‐domains in which different levels of physical modelling are used to describe the flow. This way, numerical simulations based on the Navier–Stokes equations can be extended to flows in which the continuum and/or Newtonian flow assumptions break down in regions of the domain, by locally increasing the level of detail in the model. Then, the use of multiple levels of physical modelling can reduce the overall computational cost for a given level of fidelity. The present work describes the structure of a parallel computational framework for such simulations, including details of a Navier–Stokes/MD coupling, the convergence behaviour of coupled simulations as well as the parallel implementation. For the cases considered here, micro‐scale MD problems are constructed to provide viscous stresses for the Navier–Stokes equations. The first problem is the planar Poiseuille flow, for which the viscous fluxes on each cell face in the finite‐volume discretization are evaluated using MD. The second example deals with fully developed three‐dimensional channel flow, with molecular level modelling of the shear stresses in a group of cells in the domain corners. An important aspect in using shear stresses evaluated with MD in Navier–Stokes simulations is the scatter in the data due to the sampling of a finite ensemble over a limited interval. In the coupled simulations, this prevents the convergence of the system in terms of the reduction of the norm of the residual vector of the finite‐volume discretization of the macro‐domain. Solutions to this problem are discussed in the present work, along with an analysis of the effect of number of realizations and sample duration. The averaging of the apparent viscosity for each cell face, i.e. the ratio of the shear stress predicted from MD and the imposed velocity gradient, over a number of macro‐scale time steps is shown to be a simple but effective method to reach a good level of convergence of the coupled system. Finally, the parallel efficiency of the developed method is demonstrated. Copyright © 2009 John Wiley & Sons, Ltd.     

Direct numerical simulation of decaying two‐dimensional turbulence in a no‐slip square box using smoothed particle hydrodynamics

This paper explores the application of SPH to a DNS of decaying turbulence in a two‐dimensional no‐slip wall‐bounded domain. In this bounded domain, the inverse energy cascade, and a net torque exerted by the boundary, results in a spontaneous spin‐up of the fluid, leading to a typical end state of a large monopole vortex that fills the domain. The SPH simulations were compared against published results using a high‐accuracy pseudo‐spectral code. Ensemble averages of the kinetic energy, enstrophy and average vortex wavenumber compared well against the pseudo‐spectral results, as did the evolution of the total angular momentum of the fluid. However, although the pseudo‐spectral results emphasised the importance of the no‐slip boundaries as generators of long‐lived coherent vortices in the flow, no such generation was seen in the SPH results. Vorticity filaments produced at the boundary were always dissipated by the flow shortly after separating from the boundary layer. The kinetic energy spectrum of the SPH results was calculated using an SPH Fourier transform that operates directly on the disordered particles. The ensemble kinetic energy spectrum showed the expected k^?3 scaling over most of the inertial range. However, the spectrum flattened at smaller length scales (initially less than 7.5 particle spacings and growing in size over time), indicating an excess of small‐scale kinetic energy.Copyright © 2011 John Wiley & Sons, Ltd.     

Projection‐based variational multiscale method for incompressible Navier–Stokes equations in time‐dependent domains

A variational multiscale method for computations of incompressible Navier–Stokes equations in time‐dependent domains is presented. The proposed scheme is a three‐scale variational multiscale method with a projection‐based scale separation that uses an additional tensor valued space for the large scales. The resolved large and small scales are computed in a coupled way with the effects of unresolved scales confined to the resolved small scales. In particular, the Smagorinsky eddy viscosity model is used to model the effects of unresolved scales. The deforming domain is handled by the arbitrary Lagrangian–Eulerian approach and by using an elastic mesh update technique with a mesh‐dependent stiffness. Further, the choice of orthogonal finite element basis function for the resolved large scale leads to a computationally efficient scheme. Simulations of flow around a static beam attached to a square base, around an oscillating beam and around a plunging aerofoil are presented. Copyright © 2016 John Wiley & Sons, Ltd.     

Non‐self‐overlapping structured grid generation on an n‐sided surface

Most existing meshing algorithms for a 2D or shell figure requires the figure to have exactly four sides. Generating structured grids in the n‐sided parametric region of a trimmed surface thus usually requires to first partition the region into four‐sided sub‐regions. We address the automatic structured grid generation problem in an n‐sided region by fitting a planar Gregory patch so that the partition requirement is naturally avoided. However, self‐overlapping may occur in some portions of the algebraically generated grid; this severely limits its usage in most of engineering and scientific applications where a grid system with no self‐intersecting is strictly required. To solve the problem, we use a functional optimization approach to move grid nodes in the u?v domain of the trimmed surface to eliminate the self‐overlapping. The derivatives of a Gregory patch, which are extremely difficult to compute analytically, are not required in our method. Thus, our optimization algorithm compares favourably at least in terms of speed with some other mesh optimization algorithms, such as the elliptic PDE method. In addition, to overcome the difficulty of guessing a good initial position of every grid node for the conjugate gradient method, a progressive optimization algorithm is incorporated in our optimization. Experiment results are given to illustrate the usefulness and effectiveness of the presented method. Copyright © 2004 John Wiley & Sons, Ltd.     

Three‐dimensional immersed boundary conditions for moving solids in the lattice‐Boltzmann method

This paper establishes the range of validity for a previously published three‐dimensional moving solid boundary condition for the lattice‐Boltzmann method. This method was reasonably formulated from a mass and momentum balance perspective, but was only verified for a small range of (primarily two‐dimensional) problems. One of the advantages of this boundary condition is that it offers resolution at the sub‐grid scale, allowing for accurate and stable calculation of the force and torque for solids which are moving through a lattice, even for small solid sizes relative to the computational grid size. We verify the boundary condition for creeping flows by comparison to analytical solutions that include both the force and the torque on fixed and moving spheres, and then follow this with comparisons to experimental and empirical results for both fixed as well moving spheres in inertial flows. Finally, we compare simulation results to numerical results of other investigators for the settling of an offset sphere and the drafting–kissing–tumbling of two sedimenting spheres. We found that an accurate calculation of the collision‐operator weighting used to obtain sub‐grid‐scale resolution was necessary in order to prevent spikes in the velocities, forces, and moments when solid objects cross‐computational cells. The wide range of comparisons collected and presented in this paper can be used to establish the validity of other numerical models, in addition to the one examined here. Published in 2007 by John Wiley & Sons, Ltd.     

The effect of solutal undercooling on double‐diffusive convection and macrosegregation during binary alloy solidification: a numerical investigation

In this paper, we present a macroscopic numerical model that is capable of capturing the interaction between the double‐diffusive convective field and a localized fluid flow on account of solutal undercooling during non‐equilibrium solidification of binary alloys. The model is essentially based on a fixed‐grid enthalpy based control volume approach. In the present model, microscopic features pertaining to non‐equilibrium effects on account of solutal undercooling are incorporated through the formulation of a modified partition‐coefficient. The effective partition‐coefficient is numerically modelled by means of a number of macroscopically observable parameters related to the solidifying domain. This feature has made the present treatment different from micro‐macro modelling of alloy solidification, which involves certain parameters that may not be macroscopically resolvable. Numerical simulations are performed for the case of two‐dimensional transient solidification of Pb–Sn alloys (both hypoeutectic and hypereutectic) in a rectangular cavity, employing the present model. The simulation results are also compared with the corresponding experimental results quoted in the literature, and the agreement is excellent. From the results, it can be concluded that non‐equilibrium effects on account of solutal undercooling result in a more enhanced macrosegregation. Copyright © 2002 John Wiley & Sons, Ltd.     

A new multi‐domain/multi‐resolution method for large‐eddy simulation

A new multi‐domain/multi‐resolution method is presented in the framework of the large‐eddy simulation (LES). The proposed treatment at the interfaces is conceived to deal with the problem of discontinuities on the characteristic length scales met in the case of two domains having different resolutions. It gives rise to an original approach taking into account not only the discontinuous aspect on the flow fields values but also, consequently, the non‐conservative aspect of transfer of fluxes through the interfaces. This new treatment at the interfaces has been assessed successfully in the case of a subsonic compressible channel flow. Copyright © 2001 John Wiley & Sons, Ltd.     

OVERLAPPING GRIDS FOR THE SIMULATION OF PARTICLE INTERACTION AT THE MICRO SCALE

This paper describes the use of overlapping grids for the calculation of flow around single and multipleparticle configurations at the micro scale. The basic equations for calculation are those for conservation of mass and momentum which are solved using a common Finite-Volume formulation. The hydrodynamic particle-particle and particle-wall interaction can be calculated by using an overlapping or Chimera grid scheme. With the grid structuring procedure it is possible to use simple and structured grids around the particles and the overall main grid geometry. The particle grids are lapped over the main grid such that they can move independently after each time step without remeshing the whole geometry. The paper gives results for the validation of the code developed for general test cases, for a rotating ellipsoid in simple shear flow, the flow around particles attached to a wall, the motion of a particle in the vicinity of a wall and some results for the flow through a packed bed configuration.     

On the immersed boundary‐lattice Boltzmann simulations of incompressible flows with freely moving objects

For simulating freely moving problems, conventional immersed boundary‐lattice Boltzmann methods encounter two major difficulties of an extremely large flow domain and the incompressible limit. To remove these two difficulties, this work proposes an immersed boundary‐lattice Boltzmann flux solver (IB‐LBFS) in the arbitrary Lagragian–Eulerian (ALE) coordinates and establishes a dynamic similarity theory. In the ALE‐based IB‐LBFS, the flow filed is obtained by using the LBFS on a moving Cartesian mesh, and the no‐slip boundary condition is implemented by using the boundary condition‐enforced immersed boundary method. The velocity of the Cartesian mesh is set the same as the translational velocity of the freely moving object so that there is no relative motion between the plate center and the mesh. This enables the ALE‐based IB‐LBFS to study flows with a freely moving object in a large open flow domain. By normalizing the governing equations for the flow domain and the motion of rigid body, six non‐dimensional parameters are derived and maintained to be the same in both physical systems and the lattice Boltzmann framework. This similarity algorithm enables the lattice Boltzmann equation‐based solver to study a general freely moving problem within the incompressible limit. The proposed solver and dynamic similarity theory have been successfully validated by simulating the flow around an in‐line oscillating cylinder, single particle sedimentation, and flows with a freely falling plate. The obtained results agree well with both numerical and experimental data. Copyright © 2016 John Wiley & Sons, Ltd.     

TOWARDS A NEW MODEL-FREE SIMULATION OF HIGH-REYNOLDS-FLOWS: LOCAL AVERAGE DIRECT NUMERICAL SIMULATION

Studies on the numerical simulation of high-Reynolds-number flows encounter difficulties due to the wide range of characteristic length and time scales existing in the flow field. These are often much smaller than the computational grid size. A new approach based on a ‘model-free ’ local average direct numerical simulation is presented which incorporates a strategy to filter the non-resolvable scales by means of an integration over the domain and to recover the contribution of the subgrid scales by using an integral formulation developed for them. The resulting weak formulation allows us to define a numerical flux that, thanks to the filtering operation, is highly accurate. Several computation test-cases concerning theoretical accuracy and the Navier–Stokes equations at high Reynolds number are carried out without using any turbulence model. The obtained accuracy for all computations confirms that this approach can be considered a valid contribution in the field of direct numerical simulation.     

A structured but non‐uniform Cartesian grid‐based model for the shallow water equations

In large‐scale shallow flow simulations, local high‐resolution predictions are often required in order to reduce the computational cost without losing the accuracy of the solution. This is normally achieved by solving the governing equations on grids refined only to those areas of interest. Grids with varying resolution can be generated by different approaches, e.g. nesting methods, patching algorithms and adaptive unstructured or quadtree gridding techniques. This work presents a new structured but non‐uniform Cartesian grid system as an alternative to the existing approaches to provide local high‐resolution mesh. On generating a structured but non‐uniform Cartesian grid, the whole computational domain is first discretized using a coarse background grid. Local refinement is then achieved by directly allocating a specific subdivision level to each background grid cell. The neighbour information is specified by simple mathematical relationships and no explicit storage is needed. Hence, the structured property of the uniform grid is maintained. After employing some simple interpolation formulae, the governing shallow water equations are solved using a second‐order finite volume Godunov‐type scheme in a similar way as that on a uniform grid. Copyright © 2010 John Wiley & Sons, Ltd.     

Constrained moving least‐squares immersed boundary method for fluid‐structure interaction analysis

A numerical method is presented for the analysis of interactions of inviscid and compressible flows with arbitrarily shaped stationary or moving rigid solids. The fluid equations are solved on a fixed rectangular Cartesian grid by using a higher‐order finite difference method based on the fifth‐order WENO scheme. A constrained moving least‐squares sharp interface method is proposed to enforce the Neumann‐type boundary conditions on the fluid‐solid interface by using a penalty term, while the Dirichlet boundary conditions are directly enforced. The solution of the fluid flow and the solid motion equations is advanced in time by staggerly using, respectively, the third‐order Runge‐Kutta and the implicit Newmark integration schemes. The stability and the robustness of the proposed method have been demonstrated by analyzing 5 challenging problems. For these problems, the numerical results have been found to agree well with their analytical and numerical solutions available in the literature. Effects of the support domain size and values assigned to the penalty parameter on the stability and the accuracy of the present method are also discussed.     

Multi‐scale time integration for transient conjugate heat transfer

The quasi‐steady assumption is commonly adopted in existing transient fluid–solid‐coupled convection–conduction (conjugate) heat transfer simulations, which may cause non‐negligible errors in certain cases of practical interest. In the present work, we adopt a new multi‐scale framework for the fluid domain formulated in a triple‐timing form. The slow‐varying temporal gradient corresponding to the time scales in the solid domain has been effectively included in the fluid equations as a source term, whilst short‐scale unsteadiness of the fluid domain is captured by a local time integration at a given ‘frozen’ large scale time instant. For concept proof, validation and demonstration purposes, the proposed methodology has been implemented in a loosely coupled procedure in conjunction with a hybrid interfacing treatment for coupling efficiency and accuracy. The present results indicate that a much enhanced applicability can be achieved with relatively small modifications of existing transient conjugate heat transfer methods at little extra cost. Copyright © 2016 John Wiley & Sons, Ltd.