Development of SPH variable resolution using dynamic particle coalescing and splitting
B.D. Rogers, P.K. Stansby, P. Mignosa, J. Feldman
In this paper a novel variable resolution method using particle splitting and coalescing for the SPH numerical schemes... more
In this paper a novel variable resolution method using particle splitting and coalescing for the SPH numerical schemes of the Navier-Stokes equations is presented. The key idea of the scheme is to dynamically modify the particle sizes by means of splitting and coalescing individual particles. The SPH scheme adopted is variationally derived and this guaratees that both mass and momentum are conserved including particles with different smoothing lengths. To prevent highly anisotropic distributions of the particles, a generalized shifting procedure which can address also domains discretized with variable mass particles is included. The algorithm has been tested against Poiseuille flow showing that the error introduced by the splitting and coalescing is negligible.
The capability of the numerical scheme for increasing efficiency is also shown: the SPHERIC test case of the moving square in a box has shown that the particle refinement procedure is able to increase the efficiency while maintaining the same level of accuracy, as a uniform distribution with the most refined
resolution.
A correction for balancing discontinuous bed slopes in two-dimensional smoothed particle hydrodynamics shallow water modeling
R. Vacondio, B.D. Rogers, P.K. Stansby, P. Mignosa
In this paper, a smoothed particle hydrodynamics (SPH) numerical model for the shallow water equations (SWEs) with bed... more In this paper, a smoothed particle hydrodynamics (SPH) numerical model for the shallow water equations (SWEs) with bed slope source term balancing is presented. The solution of the SWEs using SPH is attractive being a conservative, mesh-free, automatically adaptive method without special treatment for wet-dry interfaces. Recently, the capability of the SPH–SWEs numerical scheme with shock capturing and general boundary conditions has been used for predicting practical flooding problems. The balance between the bed slope source term and fluxes in shallow water models is desirable for reliable simulations of flooding over bathymetries where discontinuities are present and has received some attention in the framework of Finite Volume Eulerian models. The imbalance because of the source term resulting from the calculation of the the water depth is eradicated by means of a corrected mass, which is able to remove the error introduced by a bottom discontinuity. Two different discretizations of the momentum equation are presented herein: the first one is based on the variational formulation of the SWEs in order to obtain a fully conservative formulation, whereas the second one is obtained using a non-conservative form of the free-surface elevation gradient. In both formulations, a variable smoothing length is considered. Results are presented demonstrating the corrections preserve still water in the vicinity of either 1D or 2D bed discontinuities and provide close agreement with 1D analytical solutions for rapidly varying flows over step changes in the bed. The method is finally applied to 2D dam break flow over a square obstacle where the balanced formulation improves the agreement with experimental measurements of the free surface
A high order Discontinuous Galerkin - Fourier incompressible 3D Navier-Stokes solver with rotating sliding meshes
Journal of Computational Physics
Received 28 July 2011. Revised 23 April 2012. Accepted 27 April 2012. Available online 9 May 2012.
We present the development of a sliding mesh capability for an unsteady high order (order>3) h/p Discontinuous... more
We present the development of a sliding mesh capability for an unsteady high order (order>3) h/p Discontinuous Galerkin solver for the three-dimensional incompressible Navier-Stokes equations. A high order sliding mesh method is developed and implemented for flow simulation with relative rotational motion of an inner mesh with respect to an outer static mesh, through the use of curved boundary elements and mixed triangular-quadrilateral meshes.
A second order stiffly stable method is used to discretise in time the Arbitrary Lagrangian-Eulerian form of the incompressible Navier-Stokes equations. Spatial discretisation is provided by the Symmetric Interior Penalty Galerkin formulation with modal basis functions in the it x-y plane, allowing hanging nodes and sliding meshes without the requirement to use mortar type techniques. Spatial discretisation in the z-direction is provided by a purely spectral method that uses Fourier series and allows computation of spanwise periodic three-dimensional flows. The developed solver is shown to provide high order solutions, second order in time convergence rates and spectral convergence when solving the incompressible Navier-Stokes equations on meshes where fixed and rotating elements coexist.
In addition, an exact implementation of the no-slip boundary condition is included for curved edges; circular arcs and NACA 4-digit airfoils, where analytic expressions for the geometry are used to compute the required metrics.
The solver capabilities are tested for a number of two dimensional problems governed by the incompressible Navier-Stokes equations on static and rotating meshes: the Taylor vortex problem, a static and rotating symmetric NACA0015 airfoil and flows through three bladed cross-flow turbines. In addition, three dimensional flow solutions are demonstrated for a three bladed cross-flow turbine and a circular cylinder shadowed by a pitching NACA0012 airfoil.
A Theoretical Analysis of the Observed Variability of the Geomagnetic Dipole Field
P. Hoyng, D. Schmitt, M. Ossendrijver, Phys. of the Earth and Planet. Interiors, 130, 143-157 (2002)
A mass conservative 3-D numerical model for predicting solute fluxes in estuarine waters
by Roger Falconer - Cardiff University
Paper 52: Wu, Y. and Falconer, R. A. 2000. A mass conservative 3-D numerical model for predicting solute fluxes in estuarine waters. Advances in Water Resources. 23(5), 531-543.
Abstract
A refined three-dimensional layer-integrated model to predict accurately salt and cohesive sediment... more
Abstract
A refined three-dimensional layer-integrated model to predict accurately salt and cohesive sediment transport in estuarine waters is described herein. A splitting algorithm has been used to split the three-dimensional transport equation into a horizontal two-dimensional equation and a vertical one-dimensional equation due to the different length scales. An additional source term associated with the layer average of the free-surface flow is introduced in the conservative form of the layer-integrated pollutant transport equation. The one-dimensional QUICKEST scheme has been extended to two dimensions and included in the layer-integrated advective–diffusion equation. A modified one-dimensional ULTIMATE algorithm has also been added to avoid unphysical numerical oscillations. Numerical tests for discontinuities have been carried out to study the performance of the ULTIMATE QUICKEST scheme used in the present model. The model has also been used to simulate solute transport in an idealized harbor. It has been found that the additional source term was crucial for the mass conservation of pollutant. Finally the refined model has been applied to simulate salt and cohesive sediment transport in the Humber Estuary, UK. Good agreement has been obtained with the field measured data.
An acoustic analogy formulation for moving sources in uniformly moving media
Najafi-Yazdi, A. and Bres, G.A. and Mongeau, L., (2011) Proceedings of the Royal Society A: Mathematical, Physical and Engineering Science (467)
Acoustic analogy methods are used as post-processing tools to predict aerodynamically generated sound from numerical... more Acoustic analogy methods are used as post-processing tools to predict aerodynamically generated sound from numerical solutions of unsteady flow. The Ffowcs Williams–Hawkings (FW–H) equation and related formulations, such as Farassat’s Formulations 1 and 1A, are among the commonly used analogies because of their relative low computation cost and their robustness. These formulations assume the propagation of sound waves in a medium at rest. The present paper describes a surface integral formulation based on the convective wave equation, which takes into account the presence of a mean flow. The formulation was derived to be easy to implement as a numerical post-processing tool for computational fluid dynamics codes. The new formulation constitutes one possible extension of Farassat’s Formulation 1 and 1A based on the convective form of the FW–H equation.
Robust WAF-HLL Scheme for Compressible Two-Pressure Two-Velocity Multiphase Flow Model
by K.W. Chau
Geum-Su Yeom, Keun-Shik Chang and Seung Wook Baek, Engineering Applications of Computational Fluid Mechanics, March 2012, Volume 6, Number 1, pp. 144-162.
Mathematical Modelling of Flow and Heat/Mass Transfer during Reactive Spraying Deposition Technology (RSDT) Process for High Temperature Fuel Cells
by K.W. Chau
Guilin Hu, R. Neagu, Qianpu Wang, Zhiguo Zhang, Guoneng Li and Youqu Zheng, Engineering Applications of Computational Fluid Mechanics, March 2012, Volume 6, Number 1, pp. 134-143.
Computational Fluid Dynamics Study of Large-Scale Mixing Systems with Side-Entering Impellers
by K.W. Chau
Binxin Wu, Engineering Applications of Computational Fluid Mechanics, March 2012, Volume 6, Number 1, pp. 123-133.
Numerical Simulation of Airflow and Particle Collection by Vegetative Barriers
by K.W. Chau
Li Guo and Ronaldo G. Maghirang, Engineering Applications of Computational Fluid Mechanics, March 2012, Volume 6, Number 1, pp. 110-122.
13 views
Simulating Flow over Circular Spillways by Using Different Turbulence Models
by K.W. Chau
H. Rahimzadeh, R. Maghsoodi, H. Sarkardeh and S. Tavakkol, Engineering Applications of Computational Fluid Mechanics, March 2012, Volume 6, Number 1, pp. 100-109.
11 views
Causes and Conditions for Reamer Blade Balling during Hole Enlargement while Drilling
by K.W. Chau
Xiao Hua Zhu, Jun Jing and Tong Hua, Engineering Applications of Computational Fluid Mechanics, March 2012, Volume 6, Number 1, pp. 87-99.
Discontinuous Galerkin Method for 1D Shallow Water Flows in Natural Rivers
by K.W. Chau
W. Lai and A. A. Khan, Engineering Applications of Computational Fluid Mechanics, March 2012, Volume 6, Number 1, pp. 74-86.
Numerical Analysis of Recirculation Bubble Sizes of Turbulent Co-Flowing Jet
by K.W. Chau
H. Mahmoud, W. Kriaa, H. Mhiri, G. Le Palec and P. Bournot, Engineering Applications of Computational Fluid Mechanics, March 2012, Volume 6, Number 1, pp. 58-73.
Semi-Explicit Modelling of Watersheds with Urban Drainage Systems
by K.W. Chau
Benjamin J. Dewals, Pierre Archambeau, Bruno Khuat Duy, Sébastien Erpicum and Michel Pirotton, Engineering Applications of Computational Fluid Mechanics, March 2012, Volume 6, Number 1, pp. 46-57.
Numerical Study of Flow around Iced Wind Turbine Airfoil
by K.W. Chau
Fernando Villalpando, Marcelo Reggio and Adrian Ilinca, Engineering Applications of Computational Fluid Mechanics, March 2012, Volume 6, Number 1, pp. 39-45.
22 views
Calculation of Solitary Wave Shoaling on Plane Beaches by Extended Boussinesq Equations
by K.W. Chau
Parviz Ghadimi, Mohammad Hadi Jabbari and Arsham Reisinezhad, Engineering Applications of Computational Fluid Mechanics, March 2012, Volume 6, Number 1, pp.25-38.
Numerical Modelling of Cavitation: Validation and Parametric Studies
by K.W. Chau
X. Margot, S. Hoyas, A. Gil and S. Patouna, Engineering Applications of Computational Fluid Mechanics, March 2012, Volume 6, Number 1, pp.15-24.
