Modélisation et discrétisation par éléments finis des discontinuités d interfaces pour les problèmes couplés

Transcription

1 Modélisation et discrétisation par éléments finis des discontinuités d interfaces pour les problèmes couplés Application de XFEM pour les discontinuités mobiles Andreas Kölke Assistant Professor for Fluid-Structure Interaction Institute for Structural Analysis / Computational Sciences in Engineering Beethovenstraße 51, Braunschweig, Germany a.koelke@tu-bs.de RÉSUMÉ. Le traitement numérique des problèmes couplés multiphysiques présentant des changements topologiques d interfaces fait l objet de recherches intensives. Afin de suivre de façon Eulérienne ce type de frontières mobiles, une technique basée sur des lignes de niveau est utilisée (Osher et al., 1988). Le domaine, les lignes de niveau ainsi que les différents champs physiques inconnus sont discrétisés par des éléments finis espace-temps (Hughes et al., 1988, Hübner et al., 2004). Ceci conduit à une approche monolytique du système couplé. Les espaces d approximations des différentes variables physiques sont enrichis afin de pouvoir représenter correctement les discontinuités fortes et faibles dues à la présence des interfaces (méthode des éléments finis espace-temps enrichis). Les variables sont couplées à l interface par des multiplicateurs de Lagrange distribués (Kölke et al., 2004). ABSTRACT. The treatment of boundary-coupled multifield problems with topologically changing interfaces is subject of intensive research in computational mechanics. For these types of moving boundaries, the Eulerian interface capturing level set method (Osher et al., 1988) is applied. Space-time finite elements (Hughes et al., 1988, Hübner et al., 2004) are used for discretization of the domain, the level set function and the set of physical unknowns, promoting the monolithic analysis of the coupled system. The choice of extended space-time approximation spaces for physics enables to capture strong and weak discontinuous solutions introduced by interfacial constraints and ensured by distributed Lagrange multipliers (Kölke et al., 2004). MOTS-CLÉS : frontières mobiles, discontinuités d interfaces, éléments finis espace-temps enrichis KEYWORDS: moving boundaries, interfacial discontinuities, enriched space-time finite elements. L objet. Volume 8 n 2/2005, pages 1 à 15

2 2 L objet. Volume 8 n 2/ Introduction Numerous models in engineering are concerned with the existence of strong and weak discontinuities, occuring with bi-material contact surfaces, chemical reaction fronts, solidification processes, two-fluid flow with surface tension, or propagation of discrete cracks in elastic solids. The motion of these interfaces is governed by the dynamics of the underlying physical process. The shape of physical quantities across and nearby the interface can strongly influence the behavior of the coupled multifield system. The challenge within the numerical solution of such problems is in describing the motion of a topologically changing interface or front and in the treatment of interface induced nonlinearities. 2. General Mixed Euler-Lagrange Weighted Residuals Form The space-time finite element formulation of a general boundary-coupled system, consisting of two compressible and istothermal continua with impermeable interface for the n-th space-time slab χ = Ω χ [t n, t n1 ] is described by mass balance χ i δ (ρ χ i va i va i ) ( ρ χ i t ) χ (ρ χ i vc i ) d χ, [1] χ and conservation of momentum χ i δ (ρ χ i va i ) ( (ρ χ i va i ) t χ i ) χ (ρ χ i va i vc i ) ρχ i b i d χ [2] χ δ( χ ρ χ i va i ) :x P χ d χ, [3] on the outer domain boundaries the Dirichlet boundary conditions are defined δt χ i (ρχ vi a ρχ v i a) dp χ 2 P χ D,i P χ D,i δ(ρ χ i va i ) tχ i dp χ, [4] as well as Neumann boundary conditions P χ N,i δ(ρ χ i va i ) t χ i dp χ. [5] At the common moving coupling boundary Dirichlet-type interface conditions R χ δt χ 1 (ρχ v a 1 ρχ v a 2 ) drχ, [6]

3 Éléments des discontinuités d interfaces 3 have to be fulfilled and Neumann interface conditions δ(ρ χ 1 va 1) t χ 1 drχ δ(ρ χ 2 va 2) (t χ 1 NT t ) dr χ [7] R χ R χ The initial conditions are imposed in integral form δ (ρ χ i va i (t a)) (ρ χ i va i (t a) ρ χ ) i va i,a dω χ [8] Ω χ = 0 δρ χ i, δva i, δtχ i where P χ D is the outer Dirichlet boundary, P χ N is the outer Neumann boundary of the domain and R χ is the deformable interface in the reference configuration K χ of the arbitrary Euler-Lagrange description. Eq. 1 shows the balance of mass after Petrov-Galerkin weighting, that is performed consistent to physical units. Eqs. 2 and 3 show the weak form of the momentum balance equation after Bubnov-Galerkin weighting, wherevi a is the material velocity and vi c is the material velocity in reference form, and partial integration of the stress term with the 1st Piola-Kirchhoff stress tensor x P χ. The boundary condition for mass flux ρ χ v a on the Dirichlet boundary is introduced in integral form with Eq. 4, using boundary tractions as additional variables. The second term in this equation results from partial integration of the momentum balance equation. A specified traction boundary condition t χ i on the Neumann boundary takes Eq. 5 into account. The mass conservation aspect at the common interface R of the continua in 1 and 2 represents Eq. 6, while Eq. 7 contains the dynamic interface conditions with t as source term of interfacial tractions, e.g. surface tension in two-fluid flow. The local basis N consists of the normal and tangential vectors at the interface R. In Eq. 8 the initial conditions are introduced to the variational formulation. This formulation is the basis for further derivation and simplifictaion to specialized coupled continua, such as two-fluid flow or fluid-solid systems. 3. Evolving and Topologically Changing Interfaces For computations of multi-field physics with moving boundaries, all interfaces need to be described in space and time. The level set method (Osher et al., 1988) is applied to distinguish the two fluid domains 1 and 2 using an implicit description of the interface position in space and time. The method introduces a smooth scalar function φ(x, t), which specifies the signed smallest Euclidean distance between x and the time dependent interface R.

4 4 L objet. Volume 8 n 2/2005 n 1 n n1 y t x y t x Figure 1. Discretization of the space-time domain. The motion of the interface R in presence of the fluid velocity v is described by the solution of a convection equation for φ ( ) φ δφ t v φ d = 0 δφ. [9] Note that the solution of equation (9) will not necessarily remain the distance property of φ. The definition of the level set function φ allows the identification of each phase at any time Enriched Space-Time Finite Element Method The concept of the space-time formulation (Hughes et al., 1989, Hughes et al., 1988, Tezduyar et al., 1992) is the integration of the model equations simultaneously in space and time in a single step. The domain n is discretized by a set of nodes N and a triangulation T is formed for quadrature purposes, leading to space-time finite elements. The interpolation of physical variables within each domain, e.g. the velocity v, is chosen to be discontinuous in time. Independent degrees of freedom are used at discrete time levels t n of adjacent time slabs, leading to a time-discontinuous Galerkin formulation, which is of higher order accurate, implicit, and A-stable for integration in time (Johnson, 1993). In case of linear interpolation functions used for first-order ODE s the accuracy is of third order in time. Further, a Galerkin-/Least-squares stabilization is applied for suppression of numerical oscillations in solutions to hyperbolic differential equations (Hughes et al., 1989, Hughes et al., 1988, Tezduyar et al., 1992).

5 Éléments des discontinuités d interfaces Local partition of unity FEM The fundamental idea of XFEM (Belytschko et al., 2001) based on the local partition of unity method (Babuška et al., 1997, Melenk et al., 1996) is to extend the approximation space V std defined by standard ansatz functions N i on the support of discrete nodes i N, by problem dependent constructed enrichment functions ψ. These enriching functions incorporate an a priori knowledge about the solution behavior of the considered problem. Instead of performing an enrichment of the entire computational domain, only a sub-domain with nodes j M that needs special attention is usually enhanced by a local partition of unity V ext = V std N j (x, t)ψ. [10] j M This enrichment conception is used to construct the extended approximation for spacetime finite elements w h (x, t) = N i (x, t)w i N j (x, t)ψ j (x, t)wj [11] i N j M with ψ j the space-time node specific enrichment functions ψ j (x, t) = φ(x, t) φ(x j, t j ) to assure the nodal interpolating property of the extended approximation (Chessa et al., 2003). For the treatment of strong and also weak discontinuous solutions in space and time the enrichment function ψ(x, t) = sign φ(x, t) is constructed with help of the level set value φ used for an implicit interface description. The number of additional degrees of freedom wj depends on the number of nodes M involved in the enrichment and varies from time slab to time slab since the fluidfluid interface is moving. For two-fluid flows an extended basis (11) is used for all velocity components and the pressure field Enforcing coupling conditions In the variational form, see Eqs. 1 to 8, additional interface tractions t on R are introduced for enforcing the inter-continua coupling conditions, e.g. continuity of velocity. Since the interface is moving and changing topologically the discrete description of these tractions in an explicit manner is not easy to realize. It is promising to apply an implicit technique for that as it is done for the interface itself with help of the level set formulation. The interface traction t(x, t) is expressed by a higher-dimensional function t(x, t) and the iso-surface φ(x, t) = 0 of the level set, t(x, t) = t(x, t)δ D (φ(x, t)), [12] with δ D the Dirac delta function. The embedded Lagrangian multiplier t is defined in the vicinity e of the interface R, usually only at nodes whose support is cut by the interface. This approach results in the implicit formulation of the interfacial tractions.

6 6 L objet. Volume 8 n 2/ Bibliographie Babuška I., Melenk J. M., «The Partition of Unity Method», Int. J. Numer. Methods Eng., vol. 40, p , Belytschko T., Moës T., Usui S., Parimi C., «Arbitrary Discontinuities in Finite Elements», Int. J. Numer. Methods Eng., vol. 50(4), p , Chessa J., Belytschko T., «The Extended Finite Element Method for Two-Phase Fluids», ASME J. Appl. Mech., vol. 70(1), n 1, p , Hübner B., Walhorn E., Dinkler D., «A Monolithic Approach to Fluid-Sructure Interaction using Space-Time Finite Elements», Comput. Methods Appl. Mech. Eng., vol. 193, n 23 26, p , Hughes T. J. R., Franca L. P., Hulbert G. M., «A New Finite Element Formulation for Computational Fluid Dynamics : VIII. The Galerkin / Least-Squares Method for Advective-Diffusive Equations», Comput. Methods Appl. Mech. Eng., vol. 73, p , Hughes T. J. R., Hulbert G. M., «Space-Time Finite Element Methods for Elastodynamics : Formulation and Error Estimates», Comput. Methods Appl. Mech. Eng., vol. 66, p , Johnson G., «Discontinuous Galerkin Finite Element Methods for Second Order Hyperbolic Problems», Comput. Methods Appl. Mech. Eng., vol. 107, p , Kölke A., Dinkler D., «Extended Space-Time Finite Elements for Two-Fluid Flows in Fluid- Structure Interaction», in, M. Papadrakis,, E. Onate,, B. Schrefler (eds), Proceedings of Sixth World Congress on Computational Mechanics Beijing, WCCM, Melenk J. M., Babuška I., «The Partition of Unity Finite Element Method : Basic Theory and Applications», Comput. Methods Appl. Mech. Eng., vol. 39, p , Osher S., Sethian J. A., «Fronts Propagating with Curvature-Dependent Speed : Algorithms Based on Hamilton-Jacobi Formulations», J. Comput. Phys., vol. 79, p , Tezduyar T. E., Behr M., Liou J., «A New Strategy for Finite Element Computations Involving Moving Boundaries and Interfaces - The Deforming-Spatial-Domain/Space-Time Procedure : I. The Concept and the Preliminary Numerical Tests», Comput. Methods Appl. Mech. Eng., vol. 94, p , Article reçu le 22/09/1996. Version révisée le 03/03/2005. Rédacteur responsable : GUILLAUME LAURENT SERVICE ÉDITORIAL HERMES-LAVOISIER 14 rue de Provigny, F Cachan cedex Tél : revues@lavoisier.fr Serveur web :