High-Order Unstructured Lagrangian One-Step WENO Finite Volume Schemes for Non-Conservative Hyperbolic Systems: Applications to Compressible Multi-Phase Flows

In this article we present the first better than second order accurate unstructured Lagrangian-type one-step WENO finite volume scheme for the solution of hyperbolic partial differential equations with non-conservative products. The method achieves high order of accuracy in space together with essentially non-oscillatory behavior using a nonlinear WENO reconstruction operator on unstructured triangular meshes. High order accuracy in time is obtained via a local Lagrangian space-time Galerkin predictor method that evolves the spatial reconstruction polynomials in time within each element. The final one-step finite volume scheme is derived by integration over a moving space-time control volume, where the non-conservative products are treated by a path-conservative approach that defines the jump terms on the element boundaries. The entire method is formulated as an Arbitrary-Lagrangian-Eulerian (ALE) method, where the mesh velocity can be chosen independently of the fluid velocity. The new scheme is applied to the full seven-equation Baer-Nunziato model of compressible multi-phase flows in two space dimensions. The use of a Lagrangian approach allows an excellent resolution of the solid contact and the resolution of jumps in the volume fraction. The high order of accuracy of the scheme in space and time is confirmed via a numerical convergence study. Finally, the proposed method is also applied to a reduced version of the compressible Baer-Nunziato model for the simulation of free surface water waves in moving domains. In particular, the phenomenon of sloshing is studied in a moving water tank and comparisons with experimental data are provided

An experimental and computational study of the vortex shape in a partially baffled agitated vessel

The vortex shape in a non-standard partially baffled agitated vessel in the form of a glass-lined, under-baffled stirred vessel has been investigated using both experimental and numerical approaches for an air/water system for different rotation speeds of the agitator. A simple and flexible experimental strategy was developed for determination of the time-averaged location of the unstable free surface using a process involving superimposition of images. CFD simulations were made to predict the vortex shape by using an Eulerian–Eulerian multiphase model coupled with a homogenous turbulence model. The simplifying assumptions of a constant bubble size, a constant drag coefficient and use of the k–ε turbulence model were made. An assessment of the capability of the numerical method to predict the vortex shape was carried out through comparison between experimental data and numerical results. Considering for comparison purposes a water isosurface volume fraction equal to 0.9, to account for the existence of air/water mixture present at the interface in the experiments, instead of the classical value of 0.5, gave very good agreement with the experimental data

Impact of thixotropy on flow patterns induced in a stirred tank : numerical and experimental studies

Agitation of a thixotropic shear-thinning fluid exhibiting a yield stress is investigated both experimentally and via simulations. Steady-state experiments are conducted at three impeller rotation rates (1, 2 and 8 s−1) for a tank stirred with an axial-impeller and flow-field measurements are made using particle image velocimetry (PIV) measurements. Threedimensional numerical simulations are also performed using the commercial CFD code ANSYS CFX10.0. The viscosity of the suspension is determined experimentally and is modelled using two shear-dependant laws, one of which takes into account the flow instabilities of such fluids at low shear rates. At the highest impeller speed, the flow exhibits the familiar outward pumping action associated with axial-flow impellers. However, as the impeller speed decreases, a cavern is formed around the impeller, the flow generated in the vicinity of the agitator reorganizes and its pumping capacity vanishes. An unusual flow pattern, where the radial velocity dominates, is observed experimentally at the lowest stirring speed. It is found to result from wall slip effects. Using blades with rough surfaces prevents this peculiar behaviour and mainly resolves the discrepancies between the experimental and computational results

CFD prediction and validation of ship-bank interaction in a canal

This paper utilizes CFD (Computational Fluid Dynamics) methods to investigate the bank effects on a tanker moving straight ahead at low speed in a canal characterized by surface piercing banks. For varying water depths and ship-to-bank distances, the sinkage and trim as well as the viscous hydrodynamic forces on the hull are predicted mainly by a steady state RANS (Reynolds Averaged Navier-Stokes) solver, in which the double model approximation is adopted to simulate the flat free surface. A potential flow method is also applied to evaluate the effect of the free surface and viscosity on the solutions. In addition, focus is placed on V&V (Verification and Validation) based on a grid convergence study and comparison with EFD (Experimental Fluid Dynamics) data, as well as the exploration of the modelling error in RANS computations to enable more accurate and reliable predictions of the bank effects

Engineers handbook - Low-G propellant behavior

Engineers handbook applicable in prediction of low gravity behavior of liquids in rocket propellant tank

Autocatalytic plume pinch-off

A localized source of buoyancy flux in a non-reactive fluid medium creates a plume. The flux can be provided by either heat, a compositional difference between the fluid comprising the plume and its surroundings, or a combination of both. For autocatalytic plumes produced by the iodate-arsenous acid reaction, however, buoyancy is produced along the entire reacting interface between the plume and its surroundings. Buoyancy production at the moving interface drives fluid motion, which in turn generates flow that advects the reaction front. As a consequence of this interplay between fluid flow and chemical reaction, autocatalytic plumes exhibit a rich dynamics during their ascent through the reactant medium. One of the more interesting dynamical features is the production of an accelerating vortical plume head that in certain cases pinches-off and detaches from the upwelling conduit. After pinch-off, a new plume head forms in the conduit below, and this can lead to multiple generations of plume heads for a single plume initiation. We investigated the pinch-off process using both experimentation and simulation. Experiments were performed using various concentrations of glycerol, in which it was found that repeated pinch-off occurs exclusively in a specific concentration range. Autocatalytic plume simulations revealed that pinch-off is triggered by the appearance of accelerating flow in the plume conduit.Comment: 10 figures. Accepted for publication in Phys Rev E. See also http://www.physics.utoronto.ca/nonlinear/papers_chemwave.htm

Bubble simulations with an interface tracking technique based on a partitioned fluid-structure interaction algorithm

Numerical techniques frequently used for the simulation of one bubble can be classified as interface tracking techniques and interface capturing techniques. Most of these techniques calculate both the flow around the bubble and the shape of the interface between the gas and the liquid with one code. In this paper, a rising axisymmetric bubble is simulated with an interface tracking technique that uses separate codes to determine the position of the gas-liquid interface and to calculate the flow around the bubble. The grid converged results correspond well with experimental data. The gas-liquid interface is conceived as a zero-mass, zero-thickness structure whose position is determined by the liquid forces, a uniform gas pressure and surface tension. Iterations between the two codes are necessary to obtain the coupled solution of both problems and these iterations are stabilized with a fluid-structure interaction (FSI) algorithm. The flow around the bubble is calculated on a moving mesh in a reference frame that rises at the same speed as the bubble. The flow solver first updates the mesh throughout the liquid domain given a position of the gas-liquid interface and then calculates the flow around the bubble. It is considered as a black box with the position of the gas-liquid interface as input and the liquid forces on the interface as output. During the iterations, a reduced-order model of the flow solver is generated from the inputs and outputs of the solver. The solver that calculates the interface position uses this model to adapt the liquid forces on the gas-liquid interface during the calculation of the interface position

Nonlinear Modeling and Verification of a Heaving Point Absorber for Wave Energy Conversion

Although the heaving Point Absorber (PA) concept is well known in wave energy conversion research, few studies focus on appropriate modelling of non-linear fluid viscous and mechanical friction dynamics. Even though these concepts are known to have non-linear effects on the hydrodynamic system, most research studies consider linearity as a starting point and in so doing have a weak approach to modelling the true dynamic behaviour, particularly close to resonance. The sole use of linear modelling leads to limited ability to develop control strategies capable of true power capture optimisation and suitable device operation. Based on a 1/50 scale cylindrical heaving PA, this research focuses on a strategy for hydrodynamic model development and experimental verification. In this study, nonlinear dynamics are considered, including the lumped effect of the fluid viscous and mechanical friction forces. The excellent correspondence between the derived non-linear model and wave tank tested PA behaviours provides a strong background for wave energy tuning and control system design