1,822 research outputs found
An adaptive fixed-mesh ALE method for free surface flows
In this work we present a Fixed-Mesh ALE method for the numerical simulation of free surface flows capable of using an adaptive finite element mesh covering a background domain. This mesh is successively refined and unrefined at each time step in order to focus the computational effort on the spatial regions where it is required. Some of the main ingredients of the formulation are the use of an Arbitrary-Lagrangian–Eulerian formulation for computing temporal derivatives, the use of stabilization terms for stabilizing convection, stabilizing the lack of compatibility between velocity and pressure interpolation spaces, and stabilizing the ill-conditioning introduced by the cuts on the background finite element mesh, and the coupling of the algorithm with an adaptive mesh refinement procedure suitable for running on distributed memory environments. Algorithmic steps for the projection between meshes are presented together with the algebraic fractional step approach used for improving the condition number of the linear systems to be solved. The method is tested in several numerical examples. The expected convergence rates both in space and time are observed. Smooth solution fields for both velocity and pressure are obtained (as a result of the contribution of the stabilization terms). Finally, a good agreement between the numerical results and the reference experimental data is obtained.Postprint (published version
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
Exponential basis functions in solution of incompressible fluid problems with moving free surfaces
In this report, a new simple meshless method is presented for the solution of incompressible
inviscid fluid flow problems with moving boundaries. A Lagrangian formulation established on
pressure, as a potential equation, is employed. In this method, the approximate solution is
expressed by a linear combination of exponential basis functions (EBFs), with complex-valued
exponents, satisfying the governing equation. Constant coefficients of the solution series are
evaluated through point collocation on the domain boundaries via a complex discrete transformation
technique. The numerical solution is performed in a time marching approach using an implicit
algorithm. In each time step, the governing equation is solved at the beginning and the end of the
step, with the aid of an intermediate geometry. The use of EBFs helps to find boundary velocities
with high accuracy leading to a precise geometry updating. The developed Lagrangian meshless
algorithm is applied to variety of linear and nonlinear benchmark problems. Non-linear sloshing
fluids in rigid rectangular two-dimensional basins are particularly addressed
Interface Tracking and Solid-Fluid Coupling Techniques with Coastal Engineering Applications
Multi-material physics arise in an innumerable amount of engineering problems. A broadly
scoped numerical model is developed and described in this thesis to simulate the dynamic interaction
of multi-fluid and solid systems. It is particularly aimed at modelling the interaction
of two immiscible fluids with solid structures in a coastal engineering context; however it can
be extended to other similar areas of research. The Navier Stokes equations governing the
fluids are solved using a combination of finite element (FEM) and control volume finite element
(CVFE) discretisations. The sharp interface between the fluids is obtained through the
compressive transport of material properties (e.g. material concentration). This behaviour is
achieved through the CVFE method and a conveniently limited flux calculation scheme based
on the Hyper-C method by Leonard (1991). Analytical and validation test cases are provided,
consisting of steady and unsteady flows. To further enhance the method, improve accuracy, and
exploit Lagrangian benefits, a novel moving mesh method is also introduced and tested. It is
essentially an Arbitrary Lagrangian Eulerian method in which the grid velocity is defined by
semi-explicitly solving an iterative functional minimisation problem.
A multi-phase approach is used to introduce solid structure modelling. In this approach,
solution of the velocity field for the fluid phase is obtained using Model B as explained by
Gidaspow (1994, page 151). Interaction between the fluid phase and the solids is achieved
through the means of a source term included in the fluid momentum equations. The interacting
force is calculated through integration of this source term and adding a buoyancy contribution.
The resulting force is passed to an external solid-dynamics model such as the Discrete Element
Method (DEM), or the combined Finite Discrete Element Method (FEMDEM).
The versatility and novelty of this combined modelling approach stems from its ability to
capture the fluid interaction with particles of random size and shape. Each of the three main
components of this thesis: the advection scheme, the moving mesh method, and the solid interaction
are individually validated, and examples of randomly shaped and sized particles are
shown. To conclude the work, the methods are combined together in the context of coastal engineering
applications, where the complex coupled problem of waves impacting on breakwater
amour units is chosen to demonstrate the simulation possibilities. The three components developed
in this thesis significantly extend the application range of already powerful tools, such
as Fluidity, for fluids-modelling and finite discrete element solids-modelling tools by bringing
them together for the first time
Recommended from our members
A zonal hybrid approach coupling FNPT with OpenFOAM for modelling wave-structure interactions with action of current
This paper presents a hybrid numerical approach, which combines a two-phase Navier- Stokes model (NS) and the fully nonlinear potential theory (FNPT), for modelling wave-structure interaction. The former governs the computational domain near the structure, where the viscous and turbulent effects are significant, and is solved by OpenFOAM/InterDyMFoam which utilising the finite volume method (FVM) with a Volume of Fluid (VOF) for the phase identification. The latter covers the rest of the domain, where the fluid may be considered as incompressible, inviscid and irrotational, and solved by using the Quasi Arbitrary Lagrangian- Eulerian finite element method (QALE-FEM). These two models are weakly coupled using a zonal (spatially hierarchical) approach. Considering the inconsistence of the solutions at the boundaries between two different sub-domains governed by two fundamentally different models, a relaxation (transitional) zone is introduced, where the velocity, pressure and surface elevations are taken as the weighted summation of the solutions by two models. In order to tackle the challenges associated and maximise the computational efficiency, further developments of the QALE-FEM have been made. These include the derivation of an arbitrary Lagrangian- Eulerian FNPT and application of a robust gradient calculation scheme for estimating the velocity. The present hybrid model is applied to the numerical simulation of a fixed horizontal cylinder subjected to a unidirectional wave with or without following current. The convergence property, the optimisation of the relaxation zone, the accuracy and the computational efficiency are discussed. Although the idea of the weakly coupling using the zonal approach is not new, the present hybrid model is the first one to couple the QALE-FEM with OpenFOAM solver and/or to be applied to numerical simulate the wave-structure interaction with presence of current
Recommended from our members
Quasi ALE finite element method for nonlinear water waves
This paper presents a newly developed quasi arbitrary Lagrangian-Eulerian finite element method (QALE-FEM) for simulating water waves based on fully nonlinear potential theory. The main difference of this method from the conventional finite element method developed by one of authors of this paper and others (see, e.g., [11] and [22]) is that the complex mesh is generated only once at the beginning and is moved at all other time steps in order to conform to the motion of the free surface and structures. This feature allows one to use an unstructured mesh with any degree of complexity without the need of regenerating it every time step, which is generally inevitable and very costly. Due to this feature, the QALE-FEM has high potential in enhancing computational efficiency when applied to problems associated with the complex interaction between large steep waves and structures since the use of an unstructured mesh in such a case is likely to be necessary. To achieve overall high efficiency, the numerical techniques involved in the QALE-FEM are developed, including the method to move interior nodes, technique to re-distribute the nodes on the free surface, scheme to calculate velocities and so on. The model is validated by water waves generated by a wavemaker in a tank and the interaction between water waves and periodic bars on the bed of tank. Satisfactory agreement is achieved with analytical solutions, experimental data and numerical results from other methods
- …