Iterative Methods for Stokes/Darcy Coupling
|
|
- Abner Wilkinson
- 5 years ago
- Views:
Transcription
1 Iterative Methods or Stokes/Darcy Coupling Marco Discacciati Ecole Polytechnique Fédérale de Lausanne Institut d Analyse et Calcul Scientiique Chair o Modelling and Scientiic Computing marco.discacciati@epl.ch Summary. We present iterative subdomain methods based on a domain decomposition approach to solve the coupled Stokes/Darcy problem using inite elements. The dependence o the convergence rate on the grid parameter h and on the physical data is discussed; some diiculties encountered when applying the algorithms are indicated together with possible improvement strategies. 1 Introduction and problem setting The simulation o incompressible lows in heterogenous media is an interesting topic with many applications: considering the particular case o ree luids which can iltrate through porous media, we recall or example the hydrological environmental applications and mass transer in biomechanics. The Stokes/Darcy coupled system provides a linear model to describe such phenomena. We consider a bounded domain Ω R d (d = 2, 3) ormed by two non-overlapping subdomains Ω and Ω p separated by a surace Γ = Ω Ω p. Ω is the region occupied by the luid whose motion is described by the Stokes equations which can be written in adimensional orm as: Re 1 u + p = u = 0 in Ω, (1) where u and p are the adimensional velocity and pressure, respectively, while Re is the Reynolds number deined as Re = L U /ν, ν > 0 being the luid kinematic viscosity and L, U a characteristic length and velocity, respectively. The iltration through the porous region Ω p is modeled using Darcy s equations, whose adimensional orm reads: u p = εre p p p u p = 0 in Ω p, (2)
2 564 Marco Discacciati where u p and p p are the adimensional luid velocity and pressure, respectively, Re p is the Reynolds number Re p = δ p U p /ν, U p being a characteristic velocity through the porous medium and δ p a characteristic pore size. Finally, ε = δ p /L p is the adimensional ratio between the micro and the macro scales in Ω p. Across the interace Γ the continuity o normal stresses and luxes is required; precisely, we impose: u n = u p n n T(u, p ) n = p p on Γ (3) τ j T(u, p ) n = αu τ j where T(u, p ) is the stress tensor, n is the unit normal outward vector to Ω and τ j n = 0 (j = 1,...,d 1); α is a dimensionless coeicient depending essentially on ν and the hydraulic conductivity o the porous medium. For an extensive discussion about these coupling conditions we reer to Discacciati et al. [2002], Jäger and Mikelić [1996], Layton et al. [2003], Payne and Straughan [1998]. The mathematical analysis o the coupled problem has been addressed in previous works concerning both the continuous case and the inite element approximation (see Discacciati and Quarteroni [2003], Layton et al. [2003]). In order to solve the coupled problem, an iterative substructuring method was proposed and analyzed in Discacciati and Quarteroni [2004]. Here, we test it on a model problem in order to investigate its eectiveness and robustness, with particular emphasis on the role that physical and grid parameters play on the convergence properties. We consider the computational domain Ω R 2, with Ω = (0, 1) (1, 2), Ω p = (0, 1) 2 and interace Γ = (0, 1) {1}. We rewrite Darcy s equation as (εre p p p ) = 0 in Ω p, and consider the ollowing analytic solution: u = (y 2 2y+1, x 2 x) T, p = 2(x+y 1)/Re +1/(3εRe p ) and p p = (x(1 x)(y 1) + y 3 /3 y 2 + y)/(εre p ) + 2x/Re. Concerning the inite element discretization, P 2 P 1 Taylor-Hood elements have been used or Stokes equations, while P 2 Lagrangian elements have been adopted or Darcy s problem. All the computational meshes are conorming on Γ. 2 Dirichlet-Neumann (DN) methods Considering the interace conditions (3) 1 and (3) 2, we can choose as scalar interace variable λ either λ = u n on Γ or λ = p p on Γ. These two choices deine two dierent DN-type methods, which can be outlined as ollows, respectively: Algorithm DN 1 0. choose λ = u n on Γ and an initial guess λ (0) on Γ ;
3 For k = 0, 1,... until convergence, Do Iterative Methods or Stokes/Darcy Coupling solve Darcy s equation with b.c. εre p p p (k+1) n = λ (k) on Γ ; 2. solve Stokes problem with b.c. n T(u (k+1), p (k+1) ) n = p p (k+1) and (3) 3 on Γ; 3. λ (k+1) = θu (k+1) n + (1 θ)λ (k) on Γ, θ (0, 1) ; End For Algorithm DN 2 0. choose λ = p p on Γ and an initial guess λ (0) on Γ ; For k = 0, 1,... until convergence, Do 1. solve Stokes problem with b.c. n T(u (k+1), p (k+1) ) n = λ (k) and (3) 3 on Γ; 2. solve Darcy s equation with b.c. εre p p p (k+1) n = u (k+1) n on Γ ; 3. λ (k+1) = θp (k+1) p + (1 θ)λ (k) on Γ, θ (0, 1) ; End For The two DN methods are equivalent to preconditioned Richardson methods to solve the symmetric Steklov-Poincaré equations associated to the coupled problem, and they allow to characterize optimal preconditioners or Krylov type methods (e.g. the Conjugate Gradient) or the corresponding interace problems (see Discacciati and Quarteroni [2004]). 2.1 Numerical results We consider Re = 1, εre p = 1 and tol = 10 5 ; in Table 1 we report the number o iterations or both the Richardson and the Preconditioned Conjugate Gradient (PCG) method. These convergence results are satisactory as they show the optimality o the preconditioners with respect to the grid parameter h. Number o DN 1 PCG DN 2 PCG mesh elements (θ = 0.7) λ = u n (θ = 0.7) λ = p p Table 1. Number o iterations on dierent grids with Re = 1 and εre p = 1 However, i the luid viscosity and the hydraulic conductivity decrease, small relaxation parameters θ must be adopted to guarantee convergence, in
4 566 Marco Discacciati accordance with the theoretical estimate o the upper bound θ max given in Discacciati and Quarteroni [2004]. Unortunately, in some cases θ should be so small that in practice it prevents the numerical scheme rom converging. To quote an example, i Re = 10 3 and εre p = 10 2, then θ should be unreasonably small (smaller than 10 4!) in DN 1 to prevent divergence. This diiculty should not be ascribed to the non-optimal choice o the relaxation parameter θ. In act, i we apply the PCG method which embeds the choice o the optimal acceleration parameter (see, e.g., Quarteroni et al. [2000] p. 150), the iterative algorithm converges, but the optimal properties o the preconditioners are lost, since the number o iterations depends on the mesh parameter h, as reported in Table 2. Mesh elements PCG iterations (λ = u n) Table 2. Number o iterations on dierent grids with Re = 10 3 and εre p = 10 2 On the basis o the numerical results we have obtained we can conclude that DN methods are eective only when the ratio Re /(εre p ) is suiciently small, while dealing with large values causes some diiculties. However, the latter are the very values o interest in real-lie applications and, thereore, a robust numerical method is required. 3 Dirichlet-Neumann or a time-dependent problem We introduce a ormal argument to better understand the results obtained in Sect. 2.1 and to set up a more eective numerical scheme. This approach will be treated rom a precise mathematical viewpoint in a orthcoming work Discacciati [2004]. The underlying idea is that our diiculties in solving the Stokes/Darcy problem may come rom the dierent structure o equations (1) 1 and (2) 1, which become even more dissimilar when Re 1 and εre p 1. In act, in that case, under the physically reasonable hypothesis that u and p p are suiciently small, (1) 1 reduces almost to C I + p =, while (2) 1 becomes u p +C p I = 0, where C and C p denote two positive constants 1. We rewrite (2) 1 as (εre p ) 1 u p + p p = 0 in Ω p, (4) and ormally comparing (4) to (1) 1, we are led to modiy the latter by adding a mass term like (εre p ) 1 u p as ollows:
5 Iterative Methods or Stokes/Darcy Coupling 567 β(εre p ) 1 u Re 1 u + p =, β R +, (5) possibly with a consequent modiication o the right hand side (see Remark 1) that we have denoted by. In this way we obtain a generalized Stokes momentum equation, and note that now (5) has the same behaviour o (4) in the cases o our interest, that is when Re 1 and εre p 1. We expect that the mass term β(εre p ) 1 u would help improving the positivity o the discrete Steklov-Poincaré operator which acts as a preconditioner in the DN 1 method. With this aim, we have carried out some numerical tests using the PCG algorithm with λ = u n as interace variable to solve the modiied problem (2), (5). The convergence results reported in Table 3 show that the numerical scheme has really improved. Re εre p β Iterations on the mesh with 688 el el el Table 3. Number o iterations to solve problem (2), (5) or dierent values o Re, εre p and β Remark 1. Equation (5) can be regarded as a discretization in time o the time-dependent Stokes momentum equation t u Re 1 u + p = in Ω. Precisely, i we consider β(εre p ) 1 u,n+1 Re 1 u,n+1 + p,n+1 = n+1 n 0 with n+1 = (x, t n+1 )+β(εre p ) 1 u,n, we have a backward Euler discretization in time with β(εre p ) 1 playing the role o the inverse o a time step. From the physical viewpoint, since the luid velocities in Ω are much higher than the ones through the porous medium (see Ene and Sanchez- Palencia [1975]), a time-dependent model better represents the phenomena occurring during the iltration process. 3.1 The tdn algorithm Let [0, T] be a characteristic time interval; using or the sake o simplicity the irst-order backward Euler scheme, denoting by t > 0 the time step and N = T/ t, the iterative method that we propose to solve the time-dependent coupled problem reads (the subscript n reers to the nth time level):
6 568 Marco Discacciati Algorithm tdn For n = 0,...,N 1 Do 0. choose an initial guess λ (0) n+1 or the normal velocity on Γ at the (n + 1)th time level; For k = 0, 1,... until convergence, Do 1. solve Darcy s equation with b.c. εre p p (k+1) p,n+1 n = λ(k) n+1 on Γ ; 2. solve Stokes problem ( t) 1 u (k+1),n+1 Re 1 u(k+1),n+1 + p(k+1),n+1 = ( t) 1 u,n + n+1 u (k+1),n+1 = 0 in Ω with b.c. n T(u (k+1),n+1, p(k+1),n+1 ) n = p(k+1) 3. λ (k+1) n+1 = θu (k+1),n+1 End For End For n + (1 θ)λ(k) n+1 p,n+1 and (3) 3 on Γ ; on Γ, θ (0, 1) ; 3.2 Numerical tests We consider the horizontal section o a channel 12 m long and 8 m wide which is partially occupied by a porous medium with discontinuous conductivity, as represented in Fig. 1 (let). A parabolic inlow proile is imposed on the let hand side boundary with maximal velocity equal to 0.1m/s. On the right an outlow condition is imposed. The time interval is t [0, 0.5] and the time step t = 10 3 s; or space discretization three dierent computational meshes have been adopted. In a irst case we have considered Re = and a discontinuous coeicient εre p = 10 3 in Ω p (1), εre p = 10 7 in Ω p (2). In Fig. 1 (right) we have represented the computed solution at time t = 0.05 s, while in Fig. 2 a zoom o the velocity ield through the porous medium is shown; it can be seen that the velocity is almost null in the less permeable areas o the porous medium. Finally, Table 4 (let) reports the number o iterations obtained or three computational grids at dierent time levels, showing that the number o iterations is low and independent o h. The same test has been perormed considering dierent values o the parameters: Re = , εre p = 10 1 in Ω p (1) and εre p = 10 5 in the less permeable part o the porous medium Ω p (2). The convergence results show that the number o iterations is essentially independent o these parameters, as it can be seen comparing the previous convergence results with those reported in Table 4 (right).
7 Iterative Methods or Stokes/Darcy Coupling 569 Ω Ω (2) p Outlow Inlow Γ 8m m Ω p (1) Outlow Fig. 1. Computational domain (let) and computed velocity ield at t = 0.05 s (right) Fig. 2. Zoom o the velocity ield through the porous medium Time Iterations on the mesh with level 232 el. 928 el el Time Iterations on the mesh with level 232 el. 928 el el Table 4. Number o iterations on dierent grids with Re = , εre p = 10 3 and 10 7 (let); with Re = , εre p = 10 1 and 10 5 (right) 4 Conclusions and perspectives Numerical results show that considering a time-dependent problem allows to set up a ar more eicient DN algorithm or problems with parameters in a range o physical interest. However, as we have shown, the value o t generally depends on εre p and Re, and in some cases we could be orced to consider very small time steps t 1. This could be quite annoying since one might be interested in considering long time scales, or example in modeling the iltration o pollutants in groundwater.
8 570 Marco Discacciati This limitation on t drives us to reconsider the steady coupled model. In act, should we ind an algorithm whose behaviour were as much as possible independent o the physical parameters, then not only we would be able to solve the steady problem itsel, but we could also use it in the ramework o the time-dependent model where t would be chosen under the sole requirements o stability and accuracy. A possible approach we are currently considering is a Robin-Robin type method ollowing the ideas presented in Lube et al. [2001] or Oseen equations; its analysis and numerical results will be presented in a uture work. Acknowledgement. I wish to thank Pro. A. Quarteroni, Pro. A. Valli and Dr. N. Neuss or their suggestions and the helpul discussions concerning the subject o this work. Reerences M. Discacciati. PhD thesis, Ecole Polytechnique Fédérale de Lausanne, Switzerland, In preparation. M. Discacciati, E. Miglio, and A. Quarteroni. Mathematical and numerical models or coupling surace and groundwater lows. Appl. Numer. Math., 43:57 74, M. Discacciati and A. Quarteroni. Analysis o a domain decomposition method or the coupling o Stokes and Darcy equations. In F. Brezzi, A. Bua, S. Corsaro, and A. Murli, editors, Numerical Mathematics and Advanced Applications, ENUMATH 2001, pages Springer. Milan, M. Discacciati and A. Quarteroni. Convergence analysis o a subdomain iterative method or the inite element approximation o the coupling o Stokes and Darcy equations. Comput. Visual. Sci., To be published, H. I. Ene and E. Sanchez-Palencia. Equations et phénomenès de surace pour l écoulement dans un modèle de milieu poreux. J. Mécanique, 14(1):73 108, W. Jäger and A. Mikelić. On the boundary conditions at the contact interace between a porous medium and a ree luid. Ann. Scuola Norm. Sup. Pisa Cl. Sci., 23: , W. L. Layton, F. Schieweck, and I. Yotov. Coupling luid low with porous media low. SIAM J. Num. Anal., 40: , G. Lube, L. Müller, and F. C. Otto. A nonoverlapping domain decomposition method or stabilized inite element approximations o the Oseen equations. J. Comput. Appl. Math., 132: , L. E. Payne and B. Straughan. Analysis o the boundary condition at the interace between a viscous luid and a porous medium and related modelling questions. J. Math. Pures Appl., 77: , A. Quarteroni, R. Sacco, and F. Saleri. Numerical Mathematics. Springer- Verlag, New York Berlin Heidelberg, 2000.
A Two-grid Method for Coupled Free Flow with Porous Media Flow
A Two-grid Method for Coupled Free Flow with Porous Media Flow Prince Chidyagwai a and Béatrice Rivière a, a Department of Computational and Applied Mathematics, Rice University, 600 Main Street, Houston,
More informationCoupling Stokes and Darcy equations: modeling and numerical methods
Coupling Stokes and Darcy equations: modeling and numerical methods Marco Discacciati NM2PorousMedia 24 Dubrovnik, September 29, 24 Acknowledgment: European Union Seventh Framework Programme (FP7/27-23),
More informationNumerical analysis of the Navier-Stokes/Darcy coupling
Numerical analysis o the Navier-Stokes/Darcy coupling Lori Badea 1, Marco Discacciati 2, Alio Quarteroni 2,3 1. Institute o Mathematics o the Romanian Academy, P.O. Box 1-764, 014700 Bucharest, Romania.
More informationarxiv: v1 [math.na] 29 Oct 2018
A Multirate Approach or Fluid-Structure Interaction Computation with Decoupled Methods Lian Zhang 1, Mingchao Cai 2 Mo Mu 3, arxiv:1810.11933v1 [math.na] 29 Oct 2018 Abstract We investigate a multirate
More informationConstruction of a New Domain Decomposition Method for the Stokes Equations
Construction of a New Domain Decomposition Method for the Stokes Equations Frédéric Nataf 1 and Gerd Rapin 2 1 CMAP, CNRS; UMR7641, Ecole Polytechnique, 91128 Palaiseau Cedex, France 2 Math. Dep., NAM,
More informationControlling the Heat Flux Distribution by Changing the Thickness of Heated Wall
J. Basic. Appl. Sci. Res., 2(7)7270-7275, 2012 2012, TextRoad Publication ISSN 2090-4304 Journal o Basic and Applied Scientiic Research www.textroad.com Controlling the Heat Flux Distribution by Changing
More informationCONVECTIVE HEAT TRANSFER CHARACTERISTICS OF NANOFLUIDS. Convective heat transfer analysis of nanofluid flowing inside a
Chapter 4 CONVECTIVE HEAT TRANSFER CHARACTERISTICS OF NANOFLUIDS Convective heat transer analysis o nanoluid lowing inside a straight tube o circular cross-section under laminar and turbulent conditions
More informationA solid-fluid mixture theory of porous media
A solid-luid mixture theory o porous media I-Shih Liu Instituto de Matemática, Universidade Federal do Rio de Janeiro Abstract The theories o mixtures in the ramework o continuum mechanics have been developed
More informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN
Weierstraß-Institut ür Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 0946 8633 On iterative subdomain methods or the Stokes Darcy problem Alonso Caiazzo
More informationA Simple Explanation of the Sobolev Gradient Method
A Simple Explanation o the Sobolev Gradient Method R. J. Renka July 3, 2006 Abstract We have observed that the term Sobolev gradient is used more oten than it is understood. Also, the term is oten used
More informationGlobal Weak Solution of Planetary Geostrophic Equations with Inviscid Geostrophic Balance
Global Weak Solution o Planetary Geostrophic Equations with Inviscid Geostrophic Balance Jian-Guo Liu 1, Roger Samelson 2, Cheng Wang 3 Communicated by R. Temam) Abstract. A reormulation o the planetary
More information39.1 Gradually Varied Unsteady Flow
39.1 Gradually Varied Unsteady Flow Gradually varied unsteady low occurs when the low variables such as the low depth and velocity do not change rapidly in time and space. Such lows are very common in
More informationLecture 8 Optimization
4/9/015 Lecture 8 Optimization EE 4386/5301 Computational Methods in EE Spring 015 Optimization 1 Outline Introduction 1D Optimization Parabolic interpolation Golden section search Newton s method Multidimensional
More informationTelescoping Decomposition Method for Solving First Order Nonlinear Differential Equations
Telescoping Decomposition Method or Solving First Order Nonlinear Dierential Equations 1 Mohammed Al-Reai 2 Maysem Abu-Dalu 3 Ahmed Al-Rawashdeh Abstract The Telescoping Decomposition Method TDM is a new
More informationMAGNETOHYDRODYNAMIC GO-WATER NANOFLUID FLOW AND HEAT TRANSFER BETWEEN TWO PARALLEL MOVING DISKS
THERMAL SCIENCE: Year 8, Vol., No. B, pp. 383-39 383 MAGNETOHYDRODYNAMIC GO-WATER NANOFLUID FLOW AND HEAT TRANSFER BETWEEN TWO PARALLEL MOVING DISKS Introduction by Mohammadreza AZIMI and Rouzbeh RIAZI
More informationNew constructions of domain decomposition methods for systems of PDEs
New constructions of domain decomposition methods for systems of PDEs Nouvelles constructions de méthodes de décomposition de domaine pour des systèmes d équations aux dérivées partielles V. Dolean?? F.
More informationOn Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities
On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities Heiko Berninger, Ralf Kornhuber, and Oliver Sander FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)
More informationAnalysis of Non-Thermal Equilibrium in Porous Media
Analysis o Non-Thermal Equilibrium in Porous Media A. Nouri-Borujerdi, M. Nazari 1 School o Mechanical Engineering, Shari University o Technology P.O Box 11365-9567, Tehran, Iran E-mail: anouri@shari.edu
More informationBOUNDARY LAYER ANALYSIS ALONG A STRETCHING WEDGE SURFACE WITH MAGNETIC FIELD IN A NANOFLUID
Proceedings o the International Conerence on Mechanical Engineering and Reneable Energy 7 (ICMERE7) 8 December, 7, Chittagong, Bangladesh ICMERE7-PI- BOUNDARY LAYER ANALYSIS ALONG A STRETCHING WEDGE SURFACE
More informationСИБИРСКИЕ ЭЛЕКТРОННЫЕ МАТЕМАТИЧЕСКИЕ ИЗВЕСТИЯ
... S e MR ISSN 1813-3304 СИБИРСКИЕ ЭЛЕКТРОННЫЕ МАТЕМАТИЧЕСКИЕ ИЗВЕСТИЯ Siberian Electronic Mathematical Reports http://semr.math.nsc.ru Том 13, стр. 75 88 (2016) УДК 517.95;534.21 DOI 10.17377/semi.2016.13.006
More informationRobin-Robin preconditioned Krylov methods for fluid-structure interaction problems
MOX Report No. 19/2008 Robin-Robin preconditioned Krylov methods or luid-structure interaction problems Santiago Badia, Fabio Nobile, Christian Vergara MOX, Dipartimento di Matematica F. Brioschi Politecnico
More informationA Critical Investigation of High-Order Flux Limiters In Multiphase Flow Problems
A Critical Investigation o High-Order Flux Limiters In Multiphase Flow Problems Chris Guenther Fluent In., 3647 Collins Ferry Rd., Morgantown, WV 26505, USA cpg@luent.com ABSTRACT. In recent years inite
More informationAn Immersed Interface Method for the Incompressible Navier-Stokes Equations in Irregular Domains
An Immersed Interace Method or the Incompressible Navier-Stokes Equations in Irregular Domains Duc-Vinh Le, Boo Cheong Khoo, Jaime Peraire Singapore-MIT Alliance Department o Mechanical Engineering, National
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 Semi-Analytical Solution for a Porous Channel Flow of a Non-Newtonian Fluid
Journal o Applied Fluid Mechanics, Vol. 9, No. 6, pp. 77-76, 6. Available online at www.jamonline.net, ISSN 735-357, EISSN 735-3645. A Semi-Analytical Solution or a Porous Channel Flow o a Non-Newtonian
More informationNON-SIMILAR SOLUTIONS FOR NATURAL CONVECTION FROM A MOVING VERTICAL PLATE WITH A CONVECTIVE THERMAL BOUNDARY CONDITION
NON-SIMILAR SOLUTIONS FOR NATURAL CONVECTION FROM A MOVING VERTICAL PLATE WITH A CONVECTIVE THERMAL BOUNDARY CONDITION by Asterios Pantokratoras School o Engineering, Democritus University o Thrace, 67100
More informationConstantine, Algeria. Received Accepted Keywords: Copper nanoparticles; heat transfer; circular cylinder; steady regime.
Metallurgical and Materials Engineering Association o Metallurgical Engineers o Serbia AMES Scientiic paper UDC: 669.245 NUMERICAL INVESTIGATION OF FLUID FLOW AND HEAT TRANSFER AROUND A CIRCULAR CYLINDER
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationYURI LEVIN AND ADI BEN-ISRAEL
Pp. 1447-1457 in Progress in Analysis, Vol. Heinrich G W Begehr. Robert P Gilbert and Man Wah Wong, Editors, World Scientiic, Singapore, 003, ISBN 981-38-967-9 AN INVERSE-FREE DIRECTIONAL NEWTON METHOD
More informationPressure-correction algorithm to solve Poisson system with constant coefficients for fast two-phase simulations
Center or Turbulence Research Annual Research Bries 2009 367 Pressure-correction algorithm to solve Poisson system with constant coeicients or ast two-phase simulations By D. Kim AND P. Moin 1. Motivation
More informationOE4625 Dredge Pumps and Slurry Transport. Vaclav Matousek October 13, 2004
OE465 Vaclav Matousek October 13, 004 1 Dredge Vermelding Pumps onderdeel and Slurry organisatie Transport OE465 Vaclav Matousek October 13, 004 Dredge Vermelding Pumps onderdeel and Slurry organisatie
More informationMechanical Engineering Research Journal BUOYANT FLOW OF NANOFLUID FOR HEAT-MASS TRANSFER THROUGH A THIN LAYER
Dept. o Mech. Eng. CUET Published Online March 2015 (http://www.cuet.ac.bd/merj/index.html) Mechanical Engineering Research Journal Vol. 9, pp. 712, 2013 M E R J ISSN: 1990-5491 BUOYANT FLOW OF NANOFLUID
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationHeat and Fluid Flow of Gases in Porous Media with Micropores: Slip Flow Regime
16 Heat and Fluid Flow o Gases in Porous Media with Micropores: Slip Flow Regime Moghtada Mobedi, Murat Barisik, and A. Nakayama CONTENTS 16.1 Introduction...407 16. VAM Motion and Heat Transer Equations...409
More informationDescription of a One-Dimensional Numerical Model of an Active Magnetic Regenerator Refrigerator
This is a 1D model o an active magnetic regenerative rerigerator (AMRR) that was developed in MATLAB. The model uses cycle inputs such as the luid mass low and magnetic ield proiles, luid and regenerator
More informationSIMULATION OF A FLOW THROUGH A RIGID POROUS MEDIUM SUBJECTED TO A CONSTRAINT
SIMULATION O A LOW THROUGH A RIGID POROUS MEDIUM SUBJECTED TO A CONSTRAINT R. M. Saldanha da Gama a, J. J. Pedrosa ilho a, and M. L. Martins-Costa b a Universidade do Estado do Rio de Janeiro Departamento
More informationCOMPARISON OF THERMAL CHARACTERISTICS BETWEEN THE PLATE-FIN AND PIN-FIN HEAT SINKS IN NATURAL CONVECTION
HEFAT014 10 th International Conerence on Heat Transer, Fluid Mechanics and Thermodynamics 14 6 July 014 Orlando, Florida COMPARISON OF THERMA CHARACTERISTICS BETWEEN THE PATE-FIN AND PIN-FIN HEAT SINKS
More informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationHYDROMAGNETIC DIVERGENT CHANNEL FLOW OF A VISCO- ELASTIC ELECTRICALLY CONDUCTING FLUID
Rita Choudhury et al. / International Journal o Engineering Science and Technology (IJEST) HYDROAGNETIC DIVERGENT CHANNEL FLOW OF A VISCO- ELASTIC ELECTRICALLY CONDUCTING FLUID RITA CHOUDHURY Department
More informationStokes-Darcy coupling for periodically curved interfaces
Stokes-Darcy coupling for periodically curved interfaces Sören Dobberschütz sdobber@nano.ku.dk February 21, 2014 Abstract We investigate the boundary condition between a free fluid and a porous medium,
More informationStochastic multiscale modeling of subsurface and surface flows. Part III: Multiscale mortar finite elements for coupled Stokes-Darcy flows
Stochastic multiscale modeling of subsurface and surface flows. Part III: Multiscale mortar finite elements for coupled Stokes-Darcy flows Ivan otov Department of Mathematics, University of Pittsburgh
More information3D Numerical Modelling of Convective Heat Transfer through Two-sided Vertical Channel Symmetrically Filled with Metal Foams
P Periodica Polytechnica Mechanical Engineering P 60(4), pp. 193-202, 2016 DOI: 10.3311/PPme.8511 Creative Commons Attribution b 3D Numerical Modelling o Convective Heat Transer through Two-sided Vertical
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 informationAlgorithms for fluid-structure interaction problems
9 Algorithms or luid-structure interaction problems Miguel A. Fernández, Jean-Frédéric Gerbeau Fluid-structure coupled problems and their analysis have been addressed in the previous chapters particularly
More informationGlobal Weak Solution of Planetary Geostrophic Equations with Inviscid Geostrophic Balance
Global Weak Solution o Planetary Geostrophic Equations with Inviscid Geostrophic Balance Jian-Guo Liu 1, Roger Samelson 2, Cheng Wang 3 Communicated by R. Temam Abstract. A reormulation o the planetary
More informationIOSR Journal of Mathematics (IOSR-JM) e-issn: , p-issn: X.Volume12,Issue 1 Ver. III (Jan.-Feb.2016)PP
IOSR Journal o Mathematics (IOSR-JM) e-issn:78-578, p-issn: 39-765X.Volume,Issue Ver. III (Jan.-Feb.6)PP 88- www.iosrjournals.org Eect o Chemical Reaction on MHD Boundary Layer Flow o Williamson Nanoluid
More informationOn Limiting Behavior of Contaminant Transport Models in Coupled Surface and Groundwater Flows
Axioms 25, xx, -x; doi:.339/ Article OPEN ACCESS axioms ISSN 275-68 www.mdpi.com/journal/axioms On Limiting Behavior of Contaminant Transport Models in Coupled Surface and Groundwater Flows Vincent J.
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 informationComments on Magnetohydrodynamic Unsteady Flow of A Non- Newtonian Fluid Through A Porous Medium
Comments on Magnetohydrodynamic Unsteady Flow o A Non- Newtonian Fluid Through A Porous Medium Mostaa A.A.Mahmoud Department o Mathematics, Faculty o Science, Benha University (358), Egypt Abstract The
More informationON THE QUASISTATIC APPROXIMATION IN THE STOKES-DARCY MODEL OF GROUNDWATER-SURFACE WATER FLOWS
ON THE QUASSTATC APPROXMATON N THE STOKES-DARCY MODEL OF GROUNDWATER-SURFACE WATER FLOWS Marina Moraiti Department of Mathematics University of Pittsburgh Pittsburgh, PA 156, U.S.A. Tel: +1 (41) 889-9945
More informationTwo-scale Dirichlet-Neumann preconditioners for boundary refinements
Two-scale Dirichlet-Neumann preconditioners for boundary refinements Patrice Hauret 1 and Patrick Le Tallec 2 1 Graduate Aeronautical Laboratories, MS 25-45, California Institute of Technology Pasadena,
More informationCoupling Non-Linear Stokes and Darcy Flow using Mortar Finite Elements
Coupling Non-Linear Stokes and Darcy Flow using Mortar Finite Elements V.J. Ervin E.W. Jenkins S. Sun Department of Mathematical Sciences Clemson University, Clemson, SC, 29634-0975, USA. Abstract We study
More informationc 2007 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 45, No. 5, pp. 1801 1813 c 2007 Society for Industrial and Applied Mathematics A TWO-GRID METHOD OF A MIXED STOKES DARCY MODEL FOR COUPLING FLUID FLOW WITH POROUS MEDIA FLOW MO
More information2015 American Journal of Engineering Research (AJER)
American Journal o Engineering Research (AJER) 2015 American Journal o Engineering Research (AJER) e-issn: 2320-0847 p-issn : 2320-0936 Volume-4, Issue-7, pp-33-40.ajer.org Research Paper Open Access The
More informationA posteriori error estimates applied to flow in a channel with corners
Mathematics and Computers in Simulation 61 (2003) 375 383 A posteriori error estimates applied to flow in a channel with corners Pavel Burda a,, Jaroslav Novotný b, Bedřich Sousedík a a Department of Mathematics,
More informationSimulation of electromagnetic devices using coupled models
Manuscript received Aug. 2, 2007; revised Nov. 30, 2007 Simulation o electromagnetic devices using coupled models ION CÂRSEA Department o Computer Engineering and Communication University o Craiova Str.
More informationICES REPORT April Saumik Dana and Mary F. Wheeler
ICES REPORT 18-05 April 2018 Convergence analysis o ixed stress split iterative scheme or small strain anisotropic poroelastoplasticity: A primer by Saumik Dana and Mary F. Wheeler The Institute or Computational
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationUniform inf-sup condition for the Brinkman problem in highly heterogeneous media
Uniform inf-sup condition for the Brinkman problem in highly heterogeneous media Raytcho Lazarov & Aziz Takhirov Texas A&M May 3-4, 2016 R. Lazarov & A.T. (Texas A&M) Brinkman May 3-4, 2016 1 / 30 Outline
More informationSpace-Time Nonconforming Optimized Schwarz Waveform Relaxation for Heterogeneous Problems and General Geometries
Space-Time Nonconforming Optimized Schwarz Waveform Relaxation for Heterogeneous Problems and General Geometries Laurence Halpern, Caroline Japhet, and Jérémie Szeftel 3 LAGA, Université Paris XIII, Villetaneuse
More informationA Review of Preconditioning Techniques for Steady Incompressible Flow
Zeist 2009 p. 1/43 A Review of Preconditioning Techniques for Steady Incompressible Flow David Silvester School of Mathematics University of Manchester Zeist 2009 p. 2/43 PDEs Review : 1984 2005 Update
More informationA Nonoverlapping Subdomain Algorithm with Lagrange Multipliers and its Object Oriented Implementation for Interface Problems
Contemporary Mathematics Volume 8, 998 B 0-88-0988--03030- A Nonoverlapping Subdomain Algorithm with Lagrange Multipliers and its Object Oriented Implementation for Interface Problems Daoqi Yang. Introduction
More informationNumerical Calculation of Effective Density and Compressibility Tensors in Periodic Porous Media: A Multi-Scale Asymptotic Method
Ecerpt rom the Proceedings o the COMSOL Conerence 8 Boston Numerical Calculation o Eective Density and Compressibility Tensors in Periodic Porous Media: A Multi-Scale Asymptotic Method Chang-Yong Lee,
More information1 General description of the OpenFOAM suite
Santiago Márquez Damián-Final Work-Computational Fluid Dynamics Description and utilization o interfoam multiphase solver 1 General description o the OpenFOAM suite OpenFOAM Open Field Operation and Manipulation
More informationSTABILIZED HP-DGFEM FOR INCOMPRESSIBLE FLOW. Math. Models Methods Appl. Sci., Vol. 13, pp , 2003
Mathematical Models and Methods in Applied Sciences c World Scientiic Publishing Company STABILIZED HP-DGFEM FOR INCOMPRESSIBLE FLOW D. SCHÖTZAU1 C. SCHWAB 2 A. TOSELLI 2 1 Department o Mathematics, University
More informationStrain and Stress Measurements with a Two-Dimensional Detector
Copyright ISSN (C) 97-, JCPDS-International Advances in X-ray Centre Analysis, or Volume Diraction 4 Data 999 5 Strain and Stress Measurements with a Two-Dimensional Detector Baoping Bob He and Kingsley
More informationApproximation of the Stokes-Darcy system by optimization
Approximation of the Stokes-Darcy system by optimization V.J. Ervin E.W. Jenkins Hyesuk Lee June 2, 202 Abstract A solution algorithm for the linear/nonlinear Stokes-Darcy coupled problem is proposed and
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 informationOn some numerical convergence studies of mixed finite element methods for flow in porous media
On some numerical convergence studies of mixed finite element methods for flow in porous media Gergina Pencheva Abstract We consider an expanded mixed finite element method for solving second-order elliptic
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 informationNatural convection in a vertical strip immersed in a porous medium
European Journal o Mechanics B/Fluids 22 (2003) 545 553 Natural convection in a vertical strip immersed in a porous medium L. Martínez-Suástegui a,c.treviño b,,f.méndez a a Facultad de Ingeniería, UNAM,
More informationStokes and Darcy equations, Robin-Robin domain decomposition, geometrical
A PARALLEL ROBIN-ROBIN DOMAIN DECOMPOSITION METHOD FOR THE STOKES-DARCY SYSTEM WENBIN CHEN, MAX GUNZBURGER, FEI HUA, AND XIAOMING WANG Abstract. We propose a new parallel Robin-Robin domain decomposition
More informationFILM STACKING IMPREGNATION MODEL FOR THERMOPLASTIC COMPOSITES APPLIED TO A NOVEL NET-SHAPE PREFORMING PROCESS
The 8 th International Conerence on Flow Processes in Composite Materials (FPCM8) FILM STACKING IMPREGNATION MODEL FOR THERMOPLASTIC COMPOSITES APPLIED TO A NOVEL NET-SHAPE PREFORMING PROCESS S.T.Jespersen
More informationTermination criteria for inexact fixed point methods
Termination criteria for inexact fixed point methods Philipp Birken 1 October 1, 2013 1 Institute of Mathematics, University of Kassel, Heinrich-Plett-Str. 40, D-34132 Kassel, Germany Department of Mathematics/Computer
More informationBuoyancy Driven Heat Transfer of Water-Based CuO Nanofluids in a Tilted Enclosure with a Heat Conducting Solid Cylinder on its Center
July 4-6 2012 London U.K. Buoyancy Driven Heat Transer o Water-Based CuO Nanoluids in a Tilted Enclosure with a Heat Conducting Solid Cylinder on its Center Ahmet Cihan Kamil Kahveci and Çiğdem Susantez
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 informationKeywords Perforated pinned heat sinks, Conjugate heat transfer, Electronic component cooling.
Eect o Dierent Perorations Shapes on the Thermal-hydraulic Perormance o Perorated Pinned Heat Sinks Amer Al-Damook 1,, J.L. Summers 1, N. Kapur 1, H. Thompson 1 mnajs@leeds.ac.uk, j.l.summers@leeds.ac.uk,
More informationAn Improved Expression for a Classical Type of Explicit Approximation of the Colebrook White Equation with Only One Internal Iteration
International Journal o Hydraulic Engineering 06, 5(): 9-3 DOI: 0.593/j.ijhe.06050.03 An Improved Expression or a Classical Type o Explicit Approximation o the Colebrook White Equation with Only One Internal
More informationRESOLUTION MSC.362(92) (Adopted on 14 June 2013) REVISED RECOMMENDATION ON A STANDARD METHOD FOR EVALUATING CROSS-FLOODING ARRANGEMENTS
(Adopted on 4 June 203) (Adopted on 4 June 203) ANNEX 8 (Adopted on 4 June 203) MSC 92/26/Add. Annex 8, page THE MARITIME SAFETY COMMITTEE, RECALLING Article 28(b) o the Convention on the International
More informationA New Algorithm Based SLM to Find Critical Eigenvalue in the Rayleigh-Bénard-Brinkman Convection Problem with Boundary Conditions
Advances in Computational Sciences and Technology ISSN 0973-6107 Volume 10, Number 6 (017) pp. 1557-1570 Research India Publications http://www.ripublication.com A New Algorithm Based SLM to Find Critical
More informationNumerical calculation of the electron mobility in ZnS and ZnSe semiconductors using the iterative method
International Journal o the Physical Sciences Vol. 5(11), pp. 1752-1756, 18 September, 21 Available online at http://www.academicjournals.org/ijps ISSN 1992-195 21 Academic Journals Full Length Research
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 informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
More informationAn Iterative Substructuring Method for Coupled Fluid-Solid Acoustic Problems 1
An Iterative Substructuring Method or Coupled Fluid-Solid Acoustic Problems 1 Jan Mandel Department o Mathematics and Center or Computational Mathematics, University o Colorado at Denver, Denver, CO 8217-3364,
More informationFluctuationlessness Theorem and its Application to Boundary Value Problems of ODEs
Fluctuationlessness Theorem and its Application to Boundary Value Problems o ODEs NEJLA ALTAY İstanbul Technical University Inormatics Institute Maslak, 34469, İstanbul TÜRKİYE TURKEY) nejla@be.itu.edu.tr
More informationDevelopment and Testing of a Robust Free-Surface Finite Volume Method. Marwan Darwish, Ph.D.
Development and Testing o a Robust Free-Surace Finite Volume Method Marwan Darwish, Ph.D. Faculty o Engineering and Architecture, American University o Beirut July 2003 Summary A robust numerical technique
More informationSpace-Time Domain Decomposition Methods for Transport Problems in Porous Media
1 / 49 Space-Time Domain Decomposition Methods for Transport Problems in Porous Media Thi-Thao-Phuong Hoang, Elyes Ahmed, Jérôme Jaffré, Caroline Japhet, Michel Kern, Jean Roberts INRIA Paris-Rocquencourt
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 informationNon-newtonian Rabinowitsch Fluid Effects on the Lubrication Performances of Sine Film Thrust Bearings
International Journal o Mechanical Engineering and Applications 7; 5(): 6-67 http://www.sciencepublishinggroup.com/j/ijmea doi:.648/j.ijmea.75.4 ISSN: -X (Print); ISSN: -48 (Online) Non-newtonian Rabinowitsch
More informationFLUID MECHANICS. Lecture 7 Exact solutions
FLID MECHANICS Lecture 7 Eact solutions 1 Scope o Lecture To present solutions or a ew representative laminar boundary layers where the boundary conditions enable eact analytical solutions to be obtained.
More informationWind-Driven Circulation: Stommel s gyre & Sverdrup s balance
Wind-Driven Circulation: Stommel s gyre & Sverdrup s balance We begin by returning to our system o equations or low o a layer o uniorm density on a rotating earth. du dv h + [ u( H + h)] + [ v( H t y d
More informationNon-Isothermal Displacements with Step-Profile Time Dependent Injections in Homogeneous Porous Media
Proceedings o the World Congress on Mechanical, Chemical, and Material Engineering (MCM 015) Barcelona, Spain July 0-1, 015 Paper No. 93 Non-Isothermal Displacements with Step-Proile Time Dependent Injections
More informationNUMERICAL ANALYSIS OF FORCED CONVECTION HEAT TRANSFER FROM TWO TANDEM CIRCULAR CYLINDERS EMBEDDED IN A POROUS MEDIUM
THERMAL SCIENCE, Year 2017, Vol. 21, No. 5, pp. 2117-2128 2117 Introduction NUMERICAL ANALYSIS OF FORCED CONVECTION HEAT TRANSFER FROM TWO TANDEM CIRCULAR CYLINDERS EMBEDDED IN A POROUS MEDIUM by Habib-Ollah
More informationNUMERICAL STUDY ON THE EFFECT OF INCLINATION ANGLE ON HEAT TRANSFER PERFORMANCE IN BACK-WARD FACING STEP UTILIZING NANOFLUID
NUMERICAL STUDY ON THE EFFECT OF INCLINATION ANGLE ON HEAT TRANSFER PERFORMANCE IN BACK-WARD FACING STEP UTILIZING NANOFLUID Saleh Etaig*, Etaig.Mahmoud@Northumbria.ac.uk Reaz Hasan, Reaz.Hasan@Northumria.ac.uk
More informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
More informationVector and scalar penalty-projection methods
Numerical Flow Models for Controlled Fusion - April 2007 Vector and scalar penalty-projection methods for incompressible and variable density flows Philippe Angot Université de Provence, LATP - Marseille
More informationDirichlet-Neumann and Neumann-Neumann Methods
Dirichlet-Neumann and Neumann-Neumann Methods Felix Kwok Hong Kong Baptist University Introductory Domain Decomposition Short Course DD25, Memorial University of Newfoundland July 22, 2018 Outline Methods
More informationThe achievable limits of operational modal analysis. * Siu-Kui Au 1)
The achievable limits o operational modal analysis * Siu-Kui Au 1) 1) Center or Engineering Dynamics and Institute or Risk and Uncertainty, University o Liverpool, Liverpool L69 3GH, United Kingdom 1)
More informationAdv. Theor. Appl. Mech., Vol. 7, 2014, no. 1, 1-20 HIKARI Ltd,
Adv. Theor. Appl. Mech., Vol. 7, 2014, no. 1, 1-20 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/atam.2014.31124 Magneto-Hydrodynamic Eect with Temperature Dependent Viscosity on Natural Convection
More informationRotating Flow of Magnetite-Water Nanofluid over a Stretching Surface Inspired By Non-Linear Thermal Radiation and Mass Transfer
International Journal o Mathematics Research. ISSN 0976-5840 Volume 9, Number (017), pp. 89-97 International Research Publication House http://www.irphouse.com Rotating Flow o Magnetite-Water Nanoluid
More information