A vortex level set method for the two-way coupling of an incompressible fluid with colliding rigid bodies

2016

IJESMS

De-icing systems work periodically to remove ice formed on aircraft protected surfaces. They achieve this result by reducing adhesive shear strength between ice and surfaces. In the present work, an immersed boundary method will predict ice-shedding trajectories and collisions. At the starting-point of shedding trajectory calculations, when the ice-shed particle is still close to the wing surface, there is a high probability of collision with the surface. A linear spring collision model is presented assuming the interacting bodies are connected tangentially and normally by linear springs. The verification of the model is made by comparison with literature results and study of ice trajectories around a cylinder and an aerofoil.

Section: Methodsmentioning

confidence: 99%

Mathematical model for ice wall interactions within a level set method

2016

IJESMS

“…Equation (1) is solved using a classical VIC scheme [21], detailed in [10]. The grid values of Φ , u , and are evaluated at time = Δ .…”

Section: Fluid Structure Interaction With Penalized Vic Schemementioning

confidence: 99%

Parametric Studies of Flat Plate Trajectories Using VIC and Penalization

Modelling and Simulation in Engineering

2018

Flying debris is generated in several situations: when a roof is exposed to a storm, when ice accretes on rotating wind turbines, or during inflight aircraft deicing. Four dimensionless parameters play a role in the motion of flying debris. The goal of the present paper is to investigate the relative importance of four dimensionless parameters: the Reynolds number, the Froude number, the Tachikawa number, and the mass moment of inertia parameters. Flying debris trajectories are computed with a fluid-solid interaction model formulated for an incompressible 2D laminar flow. The rigid moving solid effects are modelled in the NavierStokes equations using penalization. A VIC scheme is used to solve the flow equations. The aerodynamic forces and moments are used to compute the acceleration and the velocity of the solid. A database of 64 trajectories is built using a two-level full factorial design for the four factors. The dispersion of the plate position at a given horizontal position decreases with the Froude number. Moreover, the Tachikawa number has a significant effect on the median plate position.

“…Depending on the location in the computational domain, the velocity is either the fluid velocity u f or the solid velocity u s i . The idea is to extend the velocity field inside the solid body and to solve the flow equations with a penalization term to enforce rigid motion inside the solid [21]. Given a penalization parameter λ 1 and denoting by χ s i the characteristic function of the solid S i , the model equation is the following:…”

Section: Physicalmentioning

confidence: 99%

“…First, a vortex method [19] is proposed to simulate the interaction of an incompressible flow with rigid bodies. For our method, we followed the innovative approach proposed by Coquerelle et al [20,21]. They propose an efficient and accurate technique to simulate unsteady incompressible viscous flows using hybrid vortex methods [22].…”

Section: Introductionmentioning

confidence: 99%

Computation of Ice Shedding Trajectories Using Cartesian Grids, Penalization, and Level Sets

Modelling and Simulation in Engineering

Gallizio³

et al. 2011

We propose to model ice shedding trajectories by an innovative paradigm that is based on cartesian grids, penalization and level sets. The use of cartesian grids bypasses the meshing issue, and penalization is an efficient alternative to explicitly impose boundary conditions so that the body-fitted meshes can be avoided, making multifluid/multiphysics flows easy to set up and simulate. Level sets describe the geometry in a nonparametric way so that geometrical and topological changes due to physics and in particular shed ice pieces are straight forward to follow. The model results are verified against the case of a free falling sphere. The capabilities of the proposed model are demonstrated on ice trajectories calculations for flow around iced cylinder and airfoil.