Fourier Series Development for Solving a Steady Fluid Structure Interaction Problem
|
|
- Logan Wheeler
- 5 years ago
- Views:
Transcription
1 Int. J. Contemp. Math. Sciences, Vol. 8, 203, no. 2, HIKARI Ltd, Fourier Series Development for Solving a Steady Fluid Structure Interaction Problem Ibrahima Mbaye University of Thies, Department of Mathematics ibambaye2000@yahoo.fr Abstract This paper presents a new approximation of pressures on the interface between fluid and structure based on Fourier series development. This approximation gives easily an analytic solution to structure equation. Using the least squares method, we introduce an optimization problem to determine Fourier coefficients series and deduct the structure displacement, the fluid velocity and the fluid pressure. This present work aims to extend and improve the approximation method of pressures on the interface. Mathematics Subject Classification: 74F0, 65T40, 37M05 Keywords: Fluid-Structure, Fourier series development, Optimization Introduction Problems involving fluid structure interaction occur in a wide variety of engineering problems and therefore have attracted the interest of many investigations from different engineering disciplines. As a results, much efforts has gone into the development of general computational method for fluid structure systems [5], [6], [2], [5], [6],[3], [3], [9], [4], [0], [], [2]. We saw in a previous paper [3] that pressures on the interface was approximated by a linear combination of shape functions as well as by a development of Fourier sine series in []. In [3], fluid-structure interaction problem has been solved in the reference domain Ω F 0. Thus, this paper aims at showing a new approximation method of pressures on the interface between fluid and structure based on Fourier cosine series development. This appproximation gives easily an analytic solution to structure
2 76 I. Mbaye equation. The latter depends on Fourier coefficients series. Using the least squares method, we establish an optimization problem to compute Fourier coefficients and obtain the structure displacement, the fluid velocity and the fluid pressure. In this work, we note that the fluid structure interaction problem is solved in the moving domain Ω F u. To solve this optimization problem by the quasi-newton BFGS method we need the derivative of the cost function. Despite, it is a difficult task, the gradient of the cost function is computed in this work. In addition, the fluid is modeled by two dimensional Stoke equations for steady flow and the structure is represented by the one dimensional beam equation. 2 Fluid structure interaction problem We denote by Ω F u the two-dimensional domain occupied by the fluid, and Γ u the elastic interface between fluid and structure and Γ = Σ Σ 2 Σ 3 be the remaining external boundaries of the fluid as depicted Fig.. 2 u F u 3 2 Figure : Sets appearing to the fluid structure interaction problem Moreover, the reference domain is defined by Ω F 0 =[0,L] [0,H]. The coupled problem is to find (u, v, p) such that: ) ( (σ(p, v) n) e 2 +(u (x )) 2 for0 <x <L() EI d4 u = dx 4 u(0) = u(l) = du(0) dx = du(l) = 0 (2) dx μδv + p = f F, in Ω F u (3) v = 0, in Ω F u (4) v = g, on Σ (5) pi 2 n + μ v n = 0, on Σ 3 (6) v = 0, on Γ u (7) v 2 = 0, on Σ 2 (8) v x 2 = 0, on Σ 2 (9)
3 fluid structure interaction problem 77 Where, E is the young modulus, I is the moment of inertia, I 2 is the identity matrix, μ is the fluid viscosity, v first component of v, v 2 second component of v, the vector e 2 =(0, ), n is the unit outward normal vector, Σ 2 is the symetric axis, On Σ 2, we have the non penetration condition: v n = v 2 =0, On Σ 2, we have the continuity of Cauchy shear stress: σ n = v x 2 =0, g the velocity profil in inflow Σ. We assume that ( (σ(p, v) n) e 2 +(u (x )) 2 = p(x,h + u(x )) [3], [9]. As a consequence, equation () becomes EI d4 u = p(x dx 4,H + u(x )). 3 Fourier cosine series development We assume that on the interface between the fluid and the structure: Where ω = 2π T p(x,h+ u(x )) = α 0 + α n cos(nωx ). (0) n= is the frequency, α =(α 0,α,...,α N ) is Fourier cosine coefficients series. Now, to have numerical results, the partial sum is used: N p(x,h+ u(x )) α 0 + α n cos(nωx ). () n=
4 78 I. Mbaye Where N. Thus, this approximation gives the structure displacement u as follows: u(x )= ( x 4 N α 0 EI 24 + α n n= (nω) cos(nωx x 3 )+k k x 2 ) k 3x + k 4, (2) Thanks to equation (2) we determine constants k,k 2,k 3,k 4. Using the least squares method, we introduce the cost functional J to find Fourier coefficients. We have the following optimization problem: inf J(α) = ( L N 0 2 α n cos(nωx ) p(x,h+ u(x ))) dx, (3) n=0 subject to equations (2)-(2) and (3)-(9). To solve this optimization problem by the quasi-newton BFGS method we need the derivative of J with respect to α k, k {0,,...,N}. 4 Computation of gradient Let α k the k-th component of the vector α, we compute for each k in {0,,...,N} the partial derivative of J with respect to α k as follow: ( J L N )( = 2 a n cos(nωx ) p(x,h + u(x )) cos(kωx ) p(x ),H + u(x )) dx 0 n=0 To compute the term p(x,h+u(x )), we introduce two propositions [3], [9]. Proposition Find v (H (Ω F u ))2,v=g on Ω F u and p L2 (Ω F u )/R such that: { Ω μ v wdx F u Ω ( v)pdx = f F,w, w (H F u 0(Ω F u )) 2 Ω ( v)qdx = 0, q (4) F L2 (Ω F u u )/R have a unique solution. Thanks to the transformation T u : Ω F 0 Ω F u such that: x = x, ( x, x 2 ) Ω F 0 T u ( x, x 2 ) = (x,x 2 )= x 2 = H+u( x ) x H 2, ( x, x 2 ) Ω F 0 (5) we write the proposition in the reference domain Ω F 0 in order to partial derivatives under integral [3], [9]. After derivation, we use again the transformation T u to write equations in the moving domain Ω F u. Then, we deduct the following proposition.
5 fluid structure interaction problem 79 Proposition 2 Applications α R m v (H (Ω F u ))2 and α R m p L 2 (Ω F v u )/R are differentiable and there partial derivative (H (Ω F u ))2 and L 2 (Ω F u )/R verify the following system: p a F (α, v,w)+b F (α, w, p ) = c F (α, v, p), w (H (Ω F u ))2 Where, a F (α, v,w) = v Ω μ F u wdx b F (α, w, p ) = p Ω w F u dx c F (α, v, p) = 2 i= Ω F u b F (α, v,q) = 0 q L 2 (Ω F u )/R, a a f F i w idx (6) ( a w + μ ( a 2 a a Ω F u a x a a 2 2 ) w ) pdx x 2 2 i= Ω F u a = H + u(x ) (x,a 2 = u ) H H x 2, a = H Remark To compute the gradient, we need to solve: the structure equation ()-(2), to find u, the fluid equation (3)-(9), to find v and p, ( a v i w i μ ( a 2 a a a x x a a 2 2 ) v ) i w i dx x x 2 (7) u(x,α), a 2 = u (x,α) x 2 H the system (6), to find p(x,h+u(x )). Remark 2 In [3], author solves fluid-structure problem in reference domain Ω F 0. In the paper, we solve the fluid-structure problem in the moving domain Ω F u. 5 Numerical results We assume that the velocity on the boundary fluid domain is: { x 30( 2 2 ), (x v (x,x 2 ) = H 2,x 2 ) Σ 30, (x,x 2 ) Σ 2 v 2 (x,x 2 ) = 0, (x,x 2 ) Γ u The parameter values of the fluid and the structure are:
6 80 I. Mbaye Parameter related to fluid: The fluid viscosity is μ =0.035 g, the channel length is L =3cm, the channel width is H =0.5cm, the volume force cm s of the fluid f F =(0, 0). Paramater related to structure: The structure thickness h = 0.cm, Young s modulus is E = g, the moment of inertia is I = h3. cm s 2 2 We use the P 2 Lagrange finite element to approach velocities and P Lagrange finite element is used to approach the pressure. FreeFem++ [8] is used for the numerical tests. Table : Optimal values with BFGS method after 0 iterations, the starting point α initial =0. N ω J(α initial ) J(α optimal ) J CPU s s s s s Table 2: α optimal with BFGS method after 0 iterations, the starting point α initial =0. N α optimal ( , ) 2 ( , , ) 3 ( , , ,.27328) 4 (9.777, , ,.25828, ) 6 ( , , ,.30454, , , )
7 fluid structure interaction problem 8 25 pressure Fourier approximation 25 pressure Fourier approximation Figure 2: Pressure and approximation case Figure 3: Pressure and approximation case N=2. N=4. 25 pressure Fourier approximation Pressure N=2 Pressure N=4 Pressure N=6 Fourier N=2 Fourier N=4 Fourier N= Figure 4: Pressure and approximation case N=6. Figure 5: case N=2,4 and 6. The above figures display the pressure on the interface and the Fourier series approximation. When the number of Fourier coefficients N is increased, the approximation increases as well. It follows, therefore, that the approximation of the pressure on the interface by Fourier cosine series is becoming more accurate that greater the number of coefficients taken and tending to the pressure itself.
8 82 I. Mbaye Figure 6: velocity case N=2. Figure 7: Pressure case N=2. Figure 8: velocity case N=4. Figure 9: Pressure case N=4. Figure 0: velocity case N=6. Figure : Pressure case N=6. The above figures display the structure displacement, the fluid velocity and the pressure.
9 fluid structure interaction problem 83 6 Conclusion In this work, we introduce Fourier series development to solve fluid structure interaction problem. This approximation enable us to determine easily an analytic solution to structure equation. Using least squares method, we establish an optimization problem to compute Fourier coefficients α, the displacement u, the velocity v and the pressure p. Fourier series gives good results for N=6. It can be noted, if Fourier coefficients N is increased then J goes to zero and the coupled problem becomes equivalent to the optimization problem. This article is coming also with a view to extend and improve approximation methods of pressures on the interface. References [] S. Deparis, M.A. Fernandez, L. Formaggia. Acceleration of a fixed point algorithm for fluid-structure interaction using transpiration conditions. M2AN Math. Model. Numer. Anal., [2] L. Formaggia, J.F. Gerbeau, F. Nobile, A. Quarteroni. On the coupling of 3D and D Navier-Stokes equations for flow problems in compliant vessels. Comput. methods Appl. Mech. Engrg. 9(6-7) , 200. [3] M.A. Fernandez, M. Moubachir. An exact Block-Newton algorithm for solving fluid-structure interaction problems. C. R. Acad. Sci. Paris, Ser.I , [4] M.A. Fernandez, M. Moubachir. A Newton method using exact jacobians for solving fluid-structure coupling. Computers and Structures, vol.83, num.2-3, [5] J.F. Gerbeau, M. Vidrascu. A quasi-newton Algorithm on a reduced model for fluid-structure interaction problems in blood flows. Math. Model. Num. Anal., [6] V. Girault, P. A. Raviart. Finite element methods for Navier-Stokes equations. Theory and algorithms. Springer Series in Computational Mathematics, Springer-Verlag, Berlin, 986. [7] C. Grandmont, Y. Maday. Existence for an unsteady fluid-strucrure interaction problem. M2AN Math. Model. Numer. Anal. 34(3) , [8] F. Hecht, O. Pironneau. A finite element Software for PDE: FreeFem++. www-rocq.inria.fr/frederic.hecht, 2003.
10 84 I. Mbaye [9] I. Mbaye, C. Murea, Numerical procedure with analytic derivative for unsteady fluid-structure interaction. in Commun. Numer. Meth. Engng., 24 (2008), no., [0] I. Mbaye, Controllability approach for a fluid structure interaction problem, Applied Mathematics, 202, 3, 23-26, doi:0.4236/am [] I. Mbaye, Fourier series development resulting from an unsteady coupled problem, Int. J. Contemp. Math. Sciences, Vol. 7, 202, no. 34, [2] C. M. Murea, S. Sy, A fast method for solving fluid-structure interaction problems numerically, Int. J. Numer. Meth. Fluids, 60 (2009), no. 0, 49-72, ( DOI: 0.002/fld.93). [3] C. Murea. The BFGS algorithm for nonlinear least squares problem arising from blood flow in arteries. Comput. Math. Appl., , [4] C. Murea, C. Maday. Existence of an optimal Control for a Nonlinear Fluid-Cable Interaction Problem. Rapport de recherche CEMRACS, C.I.R.M. Lyminy, France, 996. [5] A. Osses, J.P. Puel. Approximate Controllabiity for a linear Model of Fluid Structure Interaction. ESAIM: Control, Optimisation and Calculus of Variations, vol.4, p , [6] A. Osses. A rotated multiplier applied to the controllability of waves, elasticity and tangential Stokes control. SIAM J. Control Optim.40(3) , 200. [7] Alfio Quarteroni, Massimiliano Tuveri, Alessandro Veneziani. Computational vascular fluid dynamics: Problem, models and methods. Comput Visual Sci 2: 63-97, Received: October, 202
Numerical modelling of shear-thinning non-newtonian flows in compliant vessels
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2007; 00:1 [Version: 2002/09/18 v1.01] Numerical modelling of shear-thinning non-newtonian flows in compliant vessels M.
More informationNumerical Analysis and Computation of Fluid/Fluid and Fluid/Structure Interaction Problems Jason S. Howell
J.S.Howell Research Statement 1 Numerical Analysis and Computation of Fluid/Fluid and Fluid/Structure Interaction Problems Jason S. Howell Fig. 1. Plot of the magnitude of the velocity (upper left) and
More informationA Multiphysics Strategy for Free Surface Flows
A Multiphysics Strategy for Free Surface Flows Edie Miglio, Simona Perotto, and Fausto Saleri MOX, Modeling and Scientific Computing, Department of Mathematics, Politecnico of Milano, via Bonardi 9, I-133
More informationNUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL
Proceedings of ALGORITMY 212 pp. 29 218 NUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL MARTIN HADRAVA, MILOSLAV FEISTAUER, AND PETR SVÁČEK Abstract. The present paper
More informationFREE BOUNDARY PROBLEMS IN FLUID MECHANICS
FREE BOUNDARY PROBLEMS IN FLUID MECHANICS ANA MARIA SOANE AND ROUBEN ROSTAMIAN We consider a class of free boundary problems governed by the incompressible Navier-Stokes equations. Our objective is to
More informationNumerical Simulations of the Fluid-Structure Coupling in Physiological Vessels Mini-Workshop I
Numerical Simulations of the Fluid-Structure Coupling in Physiological Vessels INRIA Rocquencourt, projet REO miguel.fernandez@inria.fr CRM Spring School, Montréal May 16, 2005 Motivations It is generally
More informationA multi-solver quasi-newton method for the partitioned simulation of fluid-structure interaction
A multi-solver quasi-newton method for the partitioned simulation of fluid-structure interaction J Degroote, S Annerel and J Vierendeels Department of Flow, Heat and Combustion Mechanics, Ghent University,
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationEXISTENCE AND REGULARITY OF SOLUTIONS FOR STOKES SYSTEMS WITH NON-SMOOTH BOUNDARY DATA IN A POLYHEDRON
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 147, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE AND REGULARITY OF SOLUTIONS FOR
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationMéthode de quasi-newton en interaction fluide-structure. J.-F. Gerbeau & M. Vidrascu
Méthode de quasi-newton en interaction fluide-structure J.-F. Gerbeau & M. Vidrascu Motivations médicales Athérosclérose : Sistole Diastole La formation des plaques est liée, entre autre, la dynamique
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationA Fully Implicit Domain Decomposition Based ALE Framework for Three-dimensional Fluid-structure Interaction with Application in Blood Flow Computation
A Fully Implicit Domain Decomposition Based ALE Framework for Three-dimensional Fluid-structure Interaction with Application in Blood Flow Computation Yuqi Wu a, Xiao-Chuan Cai b a Department of Applied
More informationL p Theory for the div-curl System
Int. Journal of Math. Analysis, Vol. 8, 2014, no. 6, 259-271 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2014.4112 L p Theory for the div-curl System Junichi Aramaki Division of Science,
More informationReduced basis method for the reliable model reduction of Navier-Stokes equations in cardiovascular modelling
Reduced basis method for the reliable model reduction of Navier-Stokes equations in cardiovascular modelling Toni Lassila, Andrea Manzoni, Gianluigi Rozza CMCS - MATHICSE - École Polytechnique Fédérale
More informationOne dimensional models for blood flow in arteries
One dimensional models for blood flow in arteries Luca Formaggia Daniele Lamponi Alfio Quarteroni Institut de Mathématiques FS/IMA) Ecole Polytecnique Fédérale de Lausanne CH-5 Lausanne Switzerland MOX
More informationNavier-Stokes Flow in Cylindrical Elastic Tubes
Navier-Stokes Flow in Cylindrical Elastic Tubes Taha Sochi University College London, Department of Physics & stronomy, Gower Street, London, WC1E 6BT Email: t.sochi@ucl.ac.uk. bstract nalytical expressions
More informationFinite volume method for nonlinear transmission problems
Finite volume method for nonlinear transmission problems Franck Boyer, Florence Hubert To cite this version: Franck Boyer, Florence Hubert. Finite volume method for nonlinear transmission problems. Proceedings
More informationGlowinski Pironneau method for the 3D ω-ψ equations
280 GUERMOND AND QUARTAPELLE Glowinski Pironneau method for the 3D ω-ψ equations Jean-Luc Guermond and Luigi Quartapelle 1 LIMSI CNRS, Orsay, France, and Dipartimento di Fisica, Politecnico di Milano,
More informationVariable Exponents Spaces and Their Applications to Fluid Dynamics
Variable Exponents Spaces and Their Applications to Fluid Dynamics Martin Rapp TU Darmstadt November 7, 213 Martin Rapp (TU Darmstadt) Variable Exponent Spaces November 7, 213 1 / 14 Overview 1 Variable
More information20. A Dual-Primal FETI Method for solving Stokes/Navier-Stokes Equations
Fourteenth International Conference on Domain Decomposition Methods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 23 DDM.org 2. A Dual-Primal FEI Method for solving Stokes/Navier-Stokes
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationLeast-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations
Least-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations by Wilhelm Heinrichs Universität Duisburg Essen, Ingenieurmathematik Universitätsstr.
More informationA Least-Squares Finite Element Approximation for the Compressible Stokes Equations
A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,
More informationEffect of Rotatory Inertia and Load Natural. Frequency on the Response of Uniform Rayleigh. Beam Resting on Pasternak Foundation
Applied Mathematical Sciences, Vol. 12, 218, no. 16, 783-795 HIKARI Ltd www.m-hikari.com https://doi.org/1.12988/ams.218.8345 Effect of Rotatory Inertia and Load Natural Frequency on the Response of Uniform
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationAnalytical Solution of Unsteady Flow of a Viscoelastic Fluid due to an Oscillating Porous Wall
Applied Mathematics 015, 5(4): 88-9 DOI: 10.593/j.am.0150504.03 Analytical Solution of Unsteady Flow of a Viscoelastic Fluid due to an Oscillating Porous Wall Abdulsalam Ya u Gital 1,*, Muhammed Abdulhameed,
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationPseudo-divergence-free element free Galerkin method for incompressible fluid flow
International Workshop on MeshFree Methods 23 Pseudo-divergence-free element free Galerkin method for incompressible fluid flow Y. Vidal () and A. Huerta (2) Abstract: Incompressible modelling in finite
More informationMathematical Models and Numerical Simulations for the Blood Flow in Large Vessels
Mathematical Models and Numerical Simulations for the Blood Flow in Large Vessels Balazs ALBERT 1 Titus PETRILA 2a Corresponding author 1 Babes-Bolyai University M. Kogalniceanu nr. 1 400084 Cluj-Napoca
More informationWeak Solutions to Nonlinear Parabolic Problems with Variable Exponent
International Journal of Mathematical Analysis Vol. 1, 216, no. 12, 553-564 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ijma.216.6223 Weak Solutions to Nonlinear Parabolic Problems with Variable
More informationThe Multiple Solutions of Laminar Flow in a. Uniformly Porous Channel with Suction/Injection
Adv. Studies Theor. Phys., Vol. 2, 28, no. 1, 473-478 The Multiple Solutions of Laminar Flow in a Uniformly Porous Channel with Suction/Injection Botong Li 1, Liancun Zheng 1, Xinxin Zhang 2, Lianxi Ma
More informationFluid-structure algorithms based on Steklov-Poincaré operators
Fluid-structure algorithms based on Steklov-Poincaré operators S. Deparis, M. Discacciati, G. Fourestey, A. Quarteroni April 25, 2005 Abstract In this paper we review some classical algorithms for fluid-structure
More informationIntroduction to Fluid Mechanics
Introduction to Fluid Mechanics Tien-Tsan Shieh April 16, 2009 What is a Fluid? The key distinction between a fluid and a solid lies in the mode of resistance to change of shape. The fluid, unlike the
More informationOn a variational inequality of Bingham and Navier-Stokes type in three dimension
PDEs for multiphase ADvanced MATerials Palazzone, Cortona (Arezzo), Italy, September 17-21, 2012 On a variational inequality of Bingham and Navier-Stokes type in three dimension Takeshi FUKAO Kyoto University
More informationSpectrum and Exact Controllability of a Hybrid System of Elasticity.
Spectrum and Exact Controllability of a Hybrid System of Elasticity. D. Mercier, January 16, 28 Abstract We consider the exact controllability of a hybrid system consisting of an elastic beam, clamped
More informationA domain decomposition framework for fluid-structure interaction problems
A domain decomposition framework for fluid-structure interaction problems Simone Deparis 1, Marco Discacciati 1, and Alfio Quarteroni 1,2 1 Institute of Analysis and Scientific Computing, CMCS-IACS, EPF
More informationOUTLET DISSIPATIVE CONDITIONS FOR AIR FLOW IN THE BRONCHIAL TREE
ESAIM: PROCEEDINGS, September 2005, Vol.14, 201-212 Eric Cancès & Jean-Frédéric Gerbeau, Editors DOI: 10.1051/proc:2005015 OUTLET DISSIPATIVE CONDITIONS FOR AIR FLOW IN THE BRONCHIAL TREE B. Maury 1, N.
More informationTWO-SCALE SIMULATION OF MAXWELL S EQUATIONS
ESAIM: PROCEEDINGS, February 27, Vol.6, 2-223 Eric Cancès & Jean-Frédéric Gerbeau, Editors DOI:.5/proc:27 TWO-SCALE SIMULATION OF MAWELL S EQUATIONS Hyam Abboud,Sébastien Jund 2,Stéphanie Salmon 2, Eric
More informationA Robust Preconditioned Iterative Method for the Navier-Stokes Equations with High Reynolds Numbers
Applied and Computational Mathematics 2017; 6(4): 202-207 http://www.sciencepublishinggroup.com/j/acm doi: 10.11648/j.acm.20170604.18 ISSN: 2328-5605 (Print); ISSN: 2328-5613 (Online) A Robust Preconditioned
More informationA. Devys 1, 2, C. Grandmont 3, B. Grec 4, B. Maury 5 and D.Yakoubi 3
ESAIM: PROCEEDINGS, August 2009, Vol. 28, p. 162-181 M. Ismail, B. Maury & J.-F. Gerbeau, Editors NUMERICAL METHOD FOR THE 2D SIMULATION OF THE RESPIRATION. A. Devys 1, 2, C. Grandmont 3, B. Grec 4, B.
More informationA Numerical Solution of Classical Van der Pol-Duffing Oscillator by He s Parameter-Expansion Method
Int. J. Contemp. Math. Sciences, Vol. 8, 2013, no. 15, 709-71 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijcms.2013.355 A Numerical Solution of Classical Van der Pol-Duffing Oscillator by
More informationPartitioned Formulation for Solving 3D Frictional Contact Problems with BEM using Localized Lagrange Multipliers
Copyright c 2007 ICCES ICCES, vol.2, no.1, pp.21-27, 2007 Partitioned Formulation for Solving 3D Frictional Contact Problems with BEM using Localized Lagrange Multipliers L. Rodríguez-Tembleque 1, J.A.
More informationCoupling techniques for partitioned fluid structure interaction simulations with black box solvers
Coupling techniques for partitioned fluid structure interaction simulations with black box solvers Joris Degroote 1, Robby Haelterman 2, Sebastiaan Annerel 1, Jan Vierendeels 1 1 Ghent University, Department
More informationNew Discretizations of Turbulent Flow Problems
New Discretizations of Turbulent Flow Problems Carolina Cardoso Manica and Songul Kaya Merdan Abstract A suitable discretization for the Zeroth Order Model in Large Eddy Simulation of turbulent flows is
More informationDesign optimization of multi-point constraints in structures
11 th World Congress on Structural and Multidisciplinary Optimization 7 th - 12 th, June 2015, Sydney Australia Design optimization of multi-point constraints in structures Daniel N. Wilke 1, Schalk Kok
More informationREE Internal Fluid Flow Sheet 2 - Solution Fundamentals of Fluid Mechanics
REE 307 - Internal Fluid Flow Sheet 2 - Solution Fundamentals of Fluid Mechanics 1. Is the following flows physically possible, that is, satisfy the continuity equation? Substitute the expressions for
More informationHIGHER-ORDER LINEARLY IMPLICIT ONE-STEP METHODS FOR THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
STUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume LIII, Number 1, March 2008 HIGHER-ORDER LINEARLY IMPLICIT ONE-STEP METHODS FOR THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IOAN TELEAGA AND JENS
More informationON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS
ON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS Abdelhafid Younsi To cite this version: Abdelhafid Younsi. ON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS. 4 pages. 212. HAL Id:
More informationA unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation
A unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation S. Bordère a and J.-P. Caltagirone b a. CNRS, Univ. Bordeaux, ICMCB,
More informationA QUASI-NEWTON ALGORITHM BASED ON A REDUCED MODEL FOR FLUID-STRUCTURE INTERACTION PROBLEMS IN BLOOD FLOWS
ESAIM: M2AN Vol. 37, N o 4, 2003, pp. 631 647 DOI: 10.1051/m2an:2003049 ESAIM: Mathematical Modelling and Numerical Analysis A QUASI-NEWTON ALGORITHM BASED ON A REDUCED MODEL OR LUID-STRUCTURE INTERACTION
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationBUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM
BUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM ERIK BURMAN AND BENJAMIN STAMM Abstract. We propose a low order discontinuous Galerkin method for incompressible flows. Stability of the
More informationNumerical control and inverse problems and Carleman estimates
Numerical control and inverse problems and Carleman estimates Enrique FERNÁNDEZ-CARA Dpto. E.D.A.N. - Univ. of Sevilla from some joint works with A. MÜNCH Lab. Mathématiques, Univ. Blaise Pascal, C-F 2,
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationVibration of Thin Beams by PIM and RPIM methods. *B. Kanber¹, and O. M. Tufik 1
APCOM & ISCM -4 th December, 23, Singapore Vibration of Thin Beams by PIM and RPIM methods *B. Kanber¹, and O. M. Tufik Mechanical Engineering Department, University of Gaziantep, Turkey. *Corresponding
More informationAnalysis of Axially Loaded Non-prismatic Beams with General End Restraints Using Differential Quadrature Method
ISBN 978-93-84422-56-1 Proceedings of International Conference on Architecture, Structure and Civil Engineering (ICASCE'15 Antalya (Turkey Sept. 7-8, 2015 pp. 1-7 Analysis of Axially Loaded Non-prismatic
More informationFLUID-STRUCTURE AND ELECTRIC INTERACTION ANALYSIS OF PIEZOELECTRIC FLAP IN A CHANNEL USING A STRONGLY COUPLED FEM SCHEME
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 2018, Glasgow, UK FLUID-STRUCTURE AND ELECTRIC INTERACTION ANALYSIS
More informationNonlinear elasticity and gels
Nonlinear elasticity and gels M. Carme Calderer School of Mathematics University of Minnesota New Mexico Analysis Seminar New Mexico State University April 4-6, 2008 1 / 23 Outline Balance laws for gels
More informationNon-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3
Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Tommaso Ruggeri Department of Mathematics and Research Center of Applied Mathematics University of Bologna January 21, 2017 ommaso
More informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationDiscontinuous Galerkin Method for interface problem of coupling different order elliptic equations
Discontinuous Galerkin Method for interface problem of coupling different order elliptic equations Igor Mozolevski, Endre Süli Federal University of Santa Catarina, Brazil Oxford University Computing Laboratory,
More informationThe Homotopy Perturbation Method for free vibration analysis of beam on elastic foundation
Shiraz University of Technology From the SelectedWorks of Habibolla Latifizadeh 2011 The Homotopy Perturbation Method for free vibration analysis of beam on elastic foundation Habibolla Latifizadeh, Shiraz
More informationControl of Interface Evolution in Multi-Phase Fluid Flows
Control of Interface Evolution in Multi-Phase Fluid Flows Markus Klein Department of Mathematics University of Tübingen Workshop on Numerical Methods for Optimal Control and Inverse Problems Garching,
More informationA Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements
W I S S E N T E C H N I K L E I D E N S C H A F T A Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements Matthias Gsell and Olaf Steinbach Institute of Computational Mathematics
More informationChristel Hohenegger A simple model for ketchup-like liquid, its numerical challenges and limitations April 7, 2011
Notes by: Andy Thaler Christel Hohenegger A simple model for ketchup-like liquid, its numerical challenges and limitations April 7, 2011 Many complex fluids are shear-thinning. Such a fluid has a shear
More informationSTEADY AND UNSTEADY 2D NUMERICAL SOLUTION OF GENERALIZED NEWTONIAN FLUIDS FLOW. Radka Keslerová, Karel Kozel
Conference Applications of Mathematics 1 in honor of the th birthday of Michal Křížek. Institute of Mathematics AS CR, Prague 1 STEADY AND UNSTEADY D NUMERICAL SOLUTION OF GENERALIZED NEWTONIAN FLUIDS
More informationComputing Time-Periodic Solutions of Nonlinear Systems of Partial Differential Equations
1 Computing ime-periodic Solutions of Nonlinear Systems of Partial Differential Equations David M. Ambrose Department of Mathematics, Drexel University, Philadelphia, PA 191 USA E-mail: ambrose@math.drexel.edu
More informationα 2 )(k x ũ(t, η) + k z W (t, η)
1 3D Poiseuille Flow Over the next two lectures we will be going over the stabilization of the 3-D Poiseuille flow. For those of you who havent had fluids, that means we have a channel of fluid that is
More informationSimulating Thin Shells with MPM
Simulating Thin Shells with MPM Biswajit Banerjee Center for the Simulation of Accidental Fires and xplosions University of Utah March 14, 2005 Outline The Problem. Shell Theory. MPM Formulation. Results.
More informationFSI with Application in Hemodynamics Analysis and Simulation
FSI with Application in Hemodynamics Analysis and Simulation Mária Lukáčová Institute of Mathematics, University of Mainz A. Hundertmark (Uni-Mainz) Š. Nečasová (Academy of Sciences, Prague) G. Rusnáková
More informationUNIT IV BOUNDARY LAYER AND FLOW THROUGH PIPES Definition of boundary layer Thickness and classification Displacement and momentum thickness Development of laminar and turbulent flows in circular pipes
More informationAnalysis of the Oseen-Viscoelastic Fluid Flow Problem
Analysis of the Oseen-Viscoelastic Fluid Flow Problem Vincent J. Ervin Hyesuk K. Lee Louis N. Ntasin Department of Mathematical Sciences Clemson University Clemson, South Carolina 9634-0975 December 30,
More informationIntegration of Sequential Quadratic Programming and Domain Decomposition Methods for Nonlinear Optimal Control Problems
Integration of Sequential Quadratic Programming and Domain Decomposition Methods for Nonlinear Optimal Control Problems Matthias Heinkenschloss 1 and Denis Ridzal 2 1 Department of Computational and Applied
More informationPSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS Jie Shen Department of Mathematics, Penn State University University Par, PA 1680, USA Abstract. We present in this
More informationA DUALITY ALGORITHM FOR THE OBSTACLE PROBLEM
Ann. Acad. Rom. Sci. Ser. Math. Appl. ISSN 2066-6594 Vol. 5, No. -2 / 203 A DUALITY ALGORITHM FOR THE OBSTACLE PROBLEM Diana Merluşcă Abstract We consider the obstacle problem in Sobolev spaces, of order
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationFluid Mechanics. Jim Emery 9/3/ Bernoulli s Law 2. 2 Torricelli s Law 3. 3 Time to Empty a Tank 3. 4 Viscosity 4. 5 The Acceleration 5
Fluid Mechanics Jim Emery 9/3/2010 Contents 1 Bernoulli s Law 2 2 Torricelli s Law 3 3 Time to Empty a Tank 3 4 Viscosity 4 5 The Acceleration 5 6 The Equation of Motion Due to Cauchy 7 7 The Deformation
More informationON EFFECTIVE IMPLICIT TIME INTEGRATION IN ANALYSIS OF FLUID-STRUCTURE PROBLEMS
SHORT COMMUNICATIONS 943 ON EFFECTIVE IMPLICIT TIME INTEGRATION IN ANALYSIS OF FLUID-STRUCTURE PROBLEMS KLAUS-JURGEN BATHE? AND VUAY SONNADS Dcpaflment of Mechanical Engineering, Massachusetts Institute
More informationSECOND-ORDER FULLY DISCRETIZED PROJECTION METHOD FOR INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
Tenth MSU Conference on Differential Equations and Computational Simulations. Electronic Journal of Differential Equations, Conference 3 (06), pp. 9 0. ISSN: 07-669. URL: http://ejde.math.txstate.edu or
More informationNon-Newtonian Fluids and Finite Elements
Non-Newtonian Fluids and Finite Elements Janice Giudice Oxford University Computing Laboratory Keble College Talk Outline Motivating Industrial Process Multiple Extrusion of Pastes Governing Equations
More informationOne-Shot Domain Decomposition Methods for Shape Optimization Problems 1 Introduction Page 565
1 One-Shot Domain Decomposition Methods for Shape 2 Optimization Problems 3 Rongliang Chen 1 and Xiao-Chuan Cai 2 4 1 College of Mathematics and Econometrics, Hunan University, Changsha, Hunan 410082,
More informationAA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow
AA210A Fundamentals of Compressible Flow Chapter 1 - Introduction to fluid flow 1 1.2 Conservation of mass Mass flux in the x-direction [ ρu ] = M L 3 L T = M L 2 T Momentum per unit volume Mass per unit
More informationarxiv: v1 [math.na] 15 Nov 2017
Noname manuscript No. (will be inserted by the editor) An HDG method with orthogonal projections in facet integrals Issei Oikawa arxiv:1711.05396v1 [math.na] 15 Nov 2017 Received: date / Accepted: date
More informationOSGOOD TYPE REGULARITY CRITERION FOR THE 3D NEWTON-BOUSSINESQ EQUATION
Electronic Journal of Differential Equations, Vol. 013 (013), No. 3, pp. 1 6. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu OSGOOD TYPE REGULARITY
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationThe Representation of Energy Equation by Laplace Transform
Int. Journal of Math. Analysis, Vol. 8, 24, no. 22, 93-97 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.2988/ijma.24.442 The Representation of Energy Equation by Laplace Transform Taehee Lee and Hwajoon
More informationMultiscale method and pseudospectral simulations for linear viscoelastic incompressible flows
Interaction and Multiscale Mechanics, Vol. 5, No. 1 (2012) 27-40 27 Multiscale method and pseudospectral simulations for linear viscoelastic incompressible flows Ling Zhang and Jie Ouyang* Department of
More informationA Remark on the Regularity of Solutions of Maxwell s Equations on Lipschitz Domains
A Remark on the Regularity of Solutions of Maxwell s Equations on Lipschitz Domains Martin Costabel Abstract Let u be a vector field on a bounded Lipschitz domain in R 3, and let u together with its divergence
More informationAn optimal control problem governed by implicit evolution quasi-variational inequalities
Annals of the University of Bucharest (mathematical series) 4 (LXII) (213), 157 166 An optimal control problem governed by implicit evolution quasi-variational inequalities Anca Capatina and Claudia Timofte
More informationOptimization and control of a separated boundary-layer flow
Optimization and control of a separated boundary-layer flow Journal: 2011 Hawaii Summer Conferences Manuscript ID: Draft lumeetingid: 2225 Date Submitted by the Author: n/a Contact Author: PASSAGGIA, Pierre-Yves
More information(Refer Slide Time: 0:35)
Fluid Dynamics And Turbo Machines. Professor Dr Shamit Bakshi. Department Of Mechanical Engineering. Indian Institute Of Technology Madras. Part A. Module-1. Lecture-4. Tutorial. (Refer Slide Time: 0:35)
More informationA Direct Proof of Caristi s Fixed Point Theorem
Applied Mathematical Sciences, Vol. 10, 2016, no. 46, 2289-2294 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2016.66190 A Direct Proof of Caristi s Fixed Point Theorem Wei-Shih Du Department
More informationSome remarks on the derivation of the Sverdrup relation.
Some remarks on the derivation of the Sverdrup relation. Didier BRESCH, Thierry COLIN Laboratoire de Mathématiques Appliquées, CNRS (UMR 662), Université Blaise Pascal (Clermont-Ferrand 2), 63177 Aubière
More informationUnsteady Flow of a Newtonian Fluid in a Contracting and Expanding Pipe
Unsteady Flow of a Newtonian Fluid in a Contracting and Expanding Pipe T S L Radhika**, M B Srinivas, T Raja Rani*, A. Karthik BITS Pilani- Hyderabad campus, Hyderabad, Telangana, India. *MTC, Muscat,
More informationA semi-implicit approach for fluid-structure interaction based on an algebraic fractional step method
A semi-implicit approach for fluid-structure interaction based on an algebraic fractional step method Annalisa Quaini, Alfio Quarteroni 3rd October 6 Chair of Modelling and Scientific Computing (CMCS),
More informationA domain decomposition algorithm for contact problems with Coulomb s friction
A domain decomposition algorithm for contact problems with Coulomb s friction J. Haslinger 1,R.Kučera 2, and T. Sassi 1 1 Introduction Contact problems of elasticity are used in many fields of science
More informationNONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS
NONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS J.-L. GUERMOND 1, Abstract. This paper analyzes a nonstandard form of the Stokes problem where the mass conservation
More information