Coupling Linear Sloshing with Six Degrees of Freedom Rigid Body Dynamics

Fluid motion in tanks is usually described in space industry with the so-called Lomen hypothesis which assumes the vorticity is null in the moving frame. We establish in this contribution that this hypothesis is valid only for uniform rotational motions. We give a more general formulation of this coupling problem, with a compact formulation. We consider the mechanical modeling of a rigid body with a motion of small amplitude, containing an incompressible fluid in the linearized regime. We first establish that the fluid motion remains irrotational in a Galilean referential if it is true at the initial time. When continuity of normal velocity and pressure are prescribed on the free surface, we establish that the global coupled problem conserves an energy functional composed by three terms. We introduce the Stokes - Zhukovsky vector fields, solving Neumann problems for the Laplace operator in the fluid in order to represent the rotational rigid motion with irrotational vector fields. Then we have a good framework to consider the coupled problem between the fluid and the rigid motion. The coupling between the free surface and the ad hoc component of the velocity potential introduces a "Neumann to Dirichlet" operator that allows to write the coupled system in a very compact form. The final expression of a Lagrangian for the coupled system is derived and the Euler-Lagrange equations of the coupled motion are presented.Comment: 23 page

Boundary Element and Finite Element Coupling for Aeroacoustics Simulations

We consider the scattering of acoustic perturbations in a presence of a flow. We suppose that the space can be split into a zone where the flow is uniform and a zone where the flow is potential. In the first zone, we apply a Prandtl-Glauert transformation to recover the Helmholtz equation. The well-known setting of boundary element method for the Helmholtz equation is available. In the second zone, the flow quantities are space dependent, we have to consider a local resolution, namely the finite element method. Herein, we carry out the coupling of these two methods and present various applications and validation test cases. The source term is given through the decomposition of an incident acoustic field on a section of the computational domain's boundary.Comment: 25 page

Lorentz Transform and Staggered Finite Differences for Advective Acoustics

38 pagesWe study acoustic wave propagation in a uniform stationary flow. We develop a method founded on the Lorentz transform and a hypothesis of irrotationality of the acoustic perturbation. After a transformation of the space-time and of the unknown fields, we derive a system of partial differential equations that eliminates the external flow and deals with the classical case of non advective acoustics. A sequel of the analysis is a new set of perfectly matched layers equations in the spirit of the work of Berenger and Collino. The numerical implementation of the previous ideas is presented with the finite differences method HaWAY on cartesian staggered grids. Relevant numerical tests are proposed

Coupling discontinuous Galerkin methods and retarded potentials for transient wave propagation on unbounded domains

International audienceThis work deals with the numerical simulation of wave propagation on unbounded domains with localized heterogeneities. To do so, we propose to combine a discretization based on a discontinuous Galerkin method in space and explicit finite differences in time on the regions containing heterogeneities with the retarded potential method to account the unbounded nature of the computational domain. The coupling formula enforces a discrete energy identity ensuring the stability under the usual CFL condition in the interior. Moreover, the scheme allows to use a smaller time step in the interior domain yielding to quasi-optimal discretization parameters for both methods. The aliasing phenomena introduced by the local time stepping are reduced by a post-processing by averaging in time obtaining a stable and second order consistent (in time) coupling algorithm. We compute the numerical rate of convergence of the method for an academic problem. The numerical results show the feasibility of the whole discretization procedure. © 2011 Elsevier Inc

Boundary Element and Finite Element Coupling for Aeroacoustics Simulations

25 pagesWe consider the scattering of acoustic perturbations in a presence of a flow. We suppose that the space can be split into a zone where the flow is uniform and a zone where the flow is potential. In the first zone, we apply a Prandtl-Glauert transformation to recover the Helmholtz equation. The well-known setting of boundary element method for the Helmholtz equation is available. In the second zone, the flow quantities are space dependent, we have to consider a local resolution, namely the finite element method. Herein, we carry out the coupling of these two methods and present various applications and validation test cases. The source term is given through the decomposition of an incident acoustic field on a section of the computational domain's boundary

Wave concept iterative procedure applied to cylinders

International audienceThe principles and advantages of an integral iterative method based on a wave concept known as the wave concept iterative procedure are outlined. The spectral operators used in the formulation are expressed in cylindrical coordinates so that they can be used to solve problems that involve cylinders. In a first step, the validity and robustness of the procedure are verified by simulations of classical structures and in a second step, simulations are successfully compared to measurements for coupled slot antennas on a perfectly conducting cylinder