2,229 research outputs found
Simulating water-entry/exit problems using Eulerian-Lagrangian and fully-Eulerian fictitious domain methods within the open-source IBAMR library
In this paper we employ two implementations of the fictitious domain (FD)
method to simulate water-entry and water-exit problems and demonstrate their
ability to simulate practical marine engineering problems. In FD methods, the
fluid momentum equation is extended within the solid domain using an additional
body force that constrains the structure velocity to be that of a rigid body.
Using this formulation, a single set of equations is solved over the entire
computational domain. The constraint force is calculated in two distinct ways:
one using an Eulerian-Lagrangian framework of the immersed boundary (IB) method
and another using a fully-Eulerian approach of the Brinkman penalization (BP)
method. Both FSI strategies use the same multiphase flow algorithm that solves
the discrete incompressible Navier-Stokes system in conservative form. A
consistent transport scheme is employed to advect mass and momentum in the
domain, which ensures numerical stability of high density ratio multiphase
flows involved in practical marine engineering applications. Example cases of a
free falling wedge (straight and inclined) and cylinder are simulated, and the
numerical results are compared against benchmark cases in literature.Comment: The current paper builds on arXiv:1901.07892 and re-explains some
parts of it for the reader's convenienc
A unified hyperbolic formulation for viscous fluids and elastoplastic solids
We discuss a unified flow theory which in a single system of hyperbolic
partial differential equations (PDEs) can describe the two main branches of
continuum mechanics, fluid dynamics, and solid dynamics. The fundamental
difference from the classical continuum models, such as the Navier-Stokes for
example, is that the finite length scale of the continuum particles is not
ignored but kept in the model in order to semi-explicitly describe the essence
of any flows, that is the process of continuum particles rearrangements. To
allow the continuum particle rearrangements, we admit the deformability of
particle which is described by the distortion field. The ability of media to
flow is characterized by the strain dissipation time which is a characteristic
time necessary for a continuum particle to rearrange with one of its
neighboring particles. It is shown that the continuum particle length scale is
intimately connected with the dissipation time. The governing equations are
represented by a system of first order hyperbolic PDEs with source terms
modeling the dissipation due to particle rearrangements. Numerical examples
justifying the reliability of the proposed approach are demonstrated.Comment: 6 figure
A conservative coupling algorithm between a compressible flow and a rigid body using an Embedded Boundary method
This paper deals with a new solid-fluid coupling algorithm between a rigid
body and an unsteady compressible fluid flow, using an Embedded Boundary
method. The coupling with a rigid body is a first step towards the coupling
with a Discrete Element method. The flow is computed using a Finite Volume
approach on a Cartesian grid. The expression of numerical fluxes does not
affect the general coupling algorithm and we use a one-step high-order scheme
proposed by Daru and Tenaud [Daru V,Tenaud C., J. Comput. Phys. 2004]. The
Embedded Boundary method is used to integrate the presence of a solid boundary
in the fluid. The coupling algorithm is totally explicit and ensures exact mass
conservation and a balance of momentum and energy between the fluid and the
solid. It is shown that the scheme preserves uniform movement of both fluid and
solid and introduces no numerical boundary roughness. The effciency of the
method is demonstrated on challenging one- and two-dimensional benchmarks
Theoretical and numerical comparison of hyperelastic and hypoelastic formulations for Eulerian non-linear elastoplasticity
The aim of this paper is to compare a hyperelastic with a hypoelastic model
describing the Eulerian dynamics of solids in the context of non-linear
elastoplastic deformations. Specifically, we consider the well-known
hypoelastic Wilkins model, which is compared against a hyperelastic model based
on the work of Godunov and Romenski. First, we discuss some general conceptual
differences between the two approaches. Second, a detailed study of both models
is proposed, where differences are made evident at the aid of deriving a
hypoelastic-type model corresponding to the hyperelastic model and a particular
equation of state used in this paper. Third, using the same high order ADER
Finite Volume and Discontinuous Galerkin methods on fixed and moving
unstructured meshes for both models, a wide range of numerical benchmark test
problems has been solved. The numerical solutions obtained for the two
different models are directly compared with each other. For small elastic
deformations, the two models produce very similar solutions that are close to
each other. However, if large elastic or elastoplastic deformations occur, the
solutions present larger differences.Comment: 14 figure
Numerical methods for the modelling of chip formation
The modeling of metal cutting has proved to be particularly complex due to the diversity of physical phenomena involved, including thermo-mechanical coupling, contact/friction and material failure. During the last few decades, there has been significant progress in the development of numerical methods for modeling machining operations. Furthermore, the most relevant techniques have been implemented in the the relevant commercial codes creating tools for the engineers working in
the design of processes and cutting devices. This paper presents a review on the numerical modeling methods and techniques used for the simulation of machining processes. The main purpose is to identify the strengths and weaknesses of each method and strategy developed up-to-now. Moreover the review covers the classical Finite Element Method covering mesh-less methods, particle-based methods
and different possibilities of Eulerian and Lagrangian approaches.Postprint (author's final draft
Recommended from our members
Finite volume modelling of low speed structural impact problems
Two investigations are described in this thesis on the common theme of applying finite volume methods to simulate structural impact problems. The first investigation is the application of the Eulerian Finite Volume Method (EFVM) to simulate the low-speed impact of ductile materials. Simulation results are validated against experiment showing that it is possible to accurately predict crater deformation profiles over the low speed speed impact regime for different projectile and substrate materials. We demonstrate how the rate dependent Johnson-Cook plasticity model is crucial to ensure correspondence to experiment.
The second investigation is concerned with the application of EFVM to simulate impact damage to thin polymeric coatings applied to the surface of metals. The aim of this work is to demonstrate how new simulation methods can help understand coating damage due stone impact. We simulate the debonding phenomenon of single layer coatings under impact by setting boundary conditions at the plate and paint interface. We show how EFVM can capture two limits of interface behaviour, sliding and separation 'slip' at one extreme and zero sliding 'welded' at the other. Results compare well to previously published experimental and simulation work, and our own finite element simulations in Abaqus. We also demonstrate how EFVM brings greater robustness and stability compared to FEM when modelling adhesive failure and higher energy impact penetration
A unified steady and unsteady formulation for hydrodynamic potential flow simulations with fully nonlinear free surface boundary conditions
This work discusses the correct modeling of the fully nonlinear free surface
boundary conditions to be prescribed in water waves flow simulations based on
potential flow theory. The main goal of such a discussion is that of
identifying a mathematical formulation and a numerical treatment that can be
used both to carry out transient simulations, and to compute steady solutions
-- for any flow admitting them. In the literature on numerical towing tank in
fact, steady and unsteady fully nonlinear potential flow solvers are
characterized by different mathematical formulations. The kinematic and dynamic
fully nonlinear free surface boundary conditions are discussed, and in
particular it is proven that the kinematic free surface boundary condition,
written in semi-Lagrangian form, can be manipulated to derive an alternative
non penetration boundary condition by all means identical to the one used on
the surface of floating bodies or on the basin bottom. The simplified
mathematical problem obtained is discretized over space and time via Boundary
Element Method (BEM) and Implicit Backward Difference Formula (BDF) scheme,
respectively. The results confirm that the solver implemented is able to solve
steady potential flow problems just by eliminating null time derivatives in the
unsteady formulation. Numerical results obtained confirm that the solver
implemented is able to accurately reproduce results of classical steady flow
solvers available in the literature.Comment: The final version of the present paper has been accepted for
publication on Applied Mathematical Modellin
Numerical Methods for the Modelling of Chip Formation
The modeling of metal cutting has proved to be particularly complex due to the diversity of physical phenomena involved, including thermo-mechanical coupling, contact/friction and material failure. During the last few decades, there has been significant progress in the development of numerical methods for modeling machining operations. Furthermore, the most relevant techniques have been implemented in the relevant commercial codes creating tools for the engineers working in the design of processes and cutting devices. This paper presents a review on the numerical modeling methods and techniques used for the simulation of machining processes. The main purpose is to identify the strengths and weaknesses of each method and strategy developed up-to-now. Moreover the review covers the classical Finite Element Method covering mesh-less methods, particle-based methods and different possibilities of Eulerian and Lagrangian approaches
- âŠ