An Immersed Interface Method for the Incompressible Navier-Stokes Equations in Irregular Domains
|
|
- Blaise Cooper
- 5 years ago
- Views:
Transcription
1 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 University o Singapore Department o Aeronautics and Astronautics, Massachusetts Institute o Technology Abstract We present an immersed interace method or the incompressible Navier Stokes equations capable o handling rigid immersed boundaries. The immersed boundary is represented by a set o Lagrangian control points. In order to guarantee that the no-slip condition on the boundary is satisied, singular orces are applied on the luid at the immersed boundary. The orces are related to the jumps in pressure and the jumps in the derivatives o both pressure and velocity, and are interpolated using cubic splines. The strength o singular orces is determined by solving a small system o equations at each time step. The Navier-Stokes equations are discretized on a staggered Cartesian grid by a second order accurate projection method or pressure and velocity. Keywords: Immersed interace method, Navier-Stokes equations, Cartesian grid method, inite dierence, ast Poisson solvers, irregular domains. I. INTRODUCTION This paper considers the immersed interace method (IIM) or the incompressible Navier-Stokes equations in general domains involving rigid boundaries. In a 2-dimensional bounded domain Ω that contains a rigid interace Γ, we consider the incompressible Navier-Stokes equations, written as u t + (u )u + p = µ u + F (1) with boundary and initial conditions u = (2) u Ω = u b (3) u(x, ) = u (4) where u is the luid velocity, p the pressure, and µ the viscosity o the luid. Here, we simply assume that the density, ρ 1, and the viscosity, µ, are constant. The singular orce F has the orm F (x, t) = (s, t)δ(x X(s, t))ds (5) Γ where X(s, t) is the arc-length parameterization o Γ, s is the arc-length, x = (x, y) is spatial position, and (s, t) is the orce density. The Navier-Stokes equations are discretized using inite dierences on a staggered Cartesian grid. The main eatures o our method are: It is a Cartesian grid method; the method does not require complex mesh generation. It is second order accurate or velocities. The Poisson-like equations resulting at each time step are solved using a cyclic reduction algorithm which has a complexity O(N logn), where N is the number o degrees o reedom. Methods utilizing a Cartesian grid or solving interace problems or problems with complex geometry have become popular in recent years. One o the most successul Cartesian grid methods is Peskin s immersed boundary (IB) method ( [1], [11], [15]). In order to deal with rigid boundaries, Lai and Peskin [11] propose to evaluate the orce density using an expression o the orm, (s, t) = κ(x e (s) X(s, t)), (6) where κ is a constant, κ 1, and X e is the arc-length parameterization o the required boundary position. The orcing term in Eq (6) is a particular case o the eedback orcing ormulation proposed by Goldstein et al. [12] with β =. In [12], the orce is expressed as (s, t) = α t U(s, t )dt + βu(s, t) (7) where U is the velocity at the control points, and α and β are chosen to be negative and large enough so that U will stay close to zero. Lima E Silva et al. [15] proposed an alternative model to compute the orce density based upon the evaluation o the various terms in the momentum equation (1) at the control points. The orce density is calculated by computing all the Navier-Stokes terms at the control points. Once the orce density is obtained at the boundary, the immersed boundary method uses a discrete delta unction to spread the orce density to the nearby Cartesian grid points. Since the IB method uses the discrete delta unction approach, it smears out sharp interace to a thickness o order o the meshwidth and it is only irst-order accurate or problems with non-smooth but continuous solutions. In contrast, the immersed interace method (IIM) can avoid this smearing and maintains a second-order accuracy by incorporating the known jumps into the inite dierence
2 scheme near the interace. The IIM was originally proposed by LeVeque and Li [3] or solving elliptic equations, and later extended to Stokes lows by LeVeque et al. [4]. The method was developed urther or the Navier-Stokes equations in Li et al. [6], Lee [9] and Le et al. [16] or problems with lexible boundaries. The method was also used by Calhoun [8] and Li et al. [7] or solving the two-dimensional streamunctionvorticity equations in irregular domains. In [7], [8] the no-slip boundary conditions are imposed directly by determining the correct jump conditions or streamunction and vorticity. Another Cartesian grid approach has been presented by Ye et al. [13] and Udaykumar et al. [14] using inite volume techniques. They reshaped the immersed boundary cells and use a polynomial interpolating unction to approximate the luxes and gradients on the aces o the boundary cells while preserving second-order accuracy. In this paper, we extend our earlier works presented in Le et al. [16] or problems with deormable boundaries, to solve problems with rigid immersed boundaries. Our approach uses the immersed interace method to solve the incompressible Navier-Stokes equations ormulated in primitive variables. In [16], the singular orce is computed based on the coniguration o the interace, i.e., the interace is assumed to be governed by either surace tension, or by an elastic membrane. In present work, the singular orce at the immersed boundary is determined to impose the no-slip boundary condition at the rigid boundary. At each time step the singular orce is computed implicitly by solving a small, dense linear system o equations. Having computed the singular orce, we then compute the jump in pressure and jumps in the derivatives o both pressure and velocity. The jumps in the solution and its derivatives are incorporated in the inite dierence discretization to obtain sharp interace resolution. Fast solvers rom Fishpack sotware library [18] are used to solve the resulting discrete systems o equations. The remainder o the paper is organized as ollows. In section II, we present the relations that must be satisied along the immersed boundary between the singular orce and the jumps in the velocity and pressure and their derivatives. In section III, we describe the generalized inite dierence approximations to the solution derivatives, which incorporate solution jumps. In section IV, we present details o the numerical algorithm. In section V, some numerical examples are presented and inally, some conclusions and suggestions or uture work are given in section VI. II. JUMP CONDITIONS ACROSS THE INTERFACE Let n and τ be the unit outward normal and tangential vectors to the interace, respectively. The normal, 1 = (s, t) n, and tangential, 2 = (s, t) τ, components o the orce density, can be related to the jump conditions or pressure and velocity as ollows (see [4], [6], [9] or details), [u] =, [µu ξ ] = 2 τ, [u η ] = (8) [p] = 1, [p ξ ] = 2 s, [p η] = 1 s. (9) [µu ηη ] = κ 2 τ, [µu ξη ] = 2 η τ κ 2n, [µu ξξ ] = [µu ηη ] + [p ξ ]n + [p η ]τ + [u ξ ]u n (1) The jump, [ ], denotes the dierence between the value o its argument outside and inside the interace, and (ξ, η) are the coordinates associated with the directions o n and τ, respectively. Here, κ is the signed valued o the curvature o the interace (i.e. we assume that n τ = k constant, so that n can point either towards, or outwards rom, the center o curvature). From expressions (8) (1) the values o the jumps o the irst and second derivatives o velocity and pressure with respect to the (x, y) coordinates are easily obtained by a simple coordinate transormation. For instance, we write, [u x ] = [u ξ ]n 1 + [u η ]τ 1 [u yy ] = [u ξξ ]n [u ξη ]n 2 τ 2 + [u ξη ]τ2 2, where n = (n 1, n 2 ) and τ = (τ 1, τ 2 ) are the Cartesian components o the normal and tangential vectors to the interace at the point considered. III. GENERALIZED FINITE DIFFERENCE FORMULAS From Taylor series expansions, it is possible to show that i the interace cuts a grid line between two grid points at x = α, x i α < x i+1, then, the ollowing approximations hold or a piecewise twice dierentiable unction v(x): v x (x i ) = v i+1 v i 1 v x (x i+1 ) = v i+2 v i (h + ) m [v (m) ] + O(h 2 ) m= 2 (h ) m [v (m) ] + O(h 2 ) m= v xx (x i ) = v i+1 2v i + v i 1 h 2 1 h 2 v xx (x i+1 ) = v i+2 2v i+1 + v i h h 2 2 (h + ) m [v (m) ] + O(h) m= 2 (h ) m [v (m) ]+O(h) m= where v (m), denotes the m-th derivative o v, v i = v(x i ), h + = x i+1 α, h = x i α, and h, is the mesh width in the x direction. The jumps in v and its derivatives are deined as [v (m) ] α = lim x α + v(m) (x) lim x α v(m) (11) in short, [ ] = [ ] α, and v () Bube [5] or more details. A. Projection method IV. NUMERICAL ALGORITHM = v. See Wiegmann and We employ a pressure-increment projection algorithm or the discretization o the Navier-Stokes equations. This projection algorithm is analogous to that presented in Brown et al. [1]. It leads to a second order accuracy or both velocity and pressure provided all the spatial derivatives are approximated to second order accuracy. The spatial discretization is carried out on a
3 standard MAC staggered grid analogous to that in Kim and Moin [2]. The ENO second-order upwind scheme is used or the advective terms (Shu and Osher [17]). With the MAC mesh, the pressure ield is deined at the cell center where the continuity equation is enorced. The velocity ields u and v are deined at the vertical edges and horizontal edges, respectively. Given the velocity u n, and the pressure p n 1/2, we compute the velocity u n+1 and pressure p n+1/2 in three steps: Step 1: Compute an intermediate velocity ield u by solving u u n = (u u) n+ 1 2 p n µ 2 u n+ 1 2 (12) u Ω = u n+1 b where, the advective term is extrapolated using the ormula, (u. u) n+ 1 2 = 3 2 (u u)n 1 2 (u u)n 1 (13) the diusion term is approximated implicitly as, 2 u n+1/2 = 1 2 ( 2 hu + 2 hu n ) + C 1. (14) and the pressure gradient term is given by p n 1 2 = G MAC p n C2 (15) The MAC gradient operators are deined as (G MAC x p) i+ 1 2,j = p i+1,j p i,j x (G MAC y p) i,j+ 1 2 = p i,j+1 p i,j y Step 2: Compute a pressure update φ n+1 by solving the Poisson equation 2 φ n+1 = u, n φ n+1 Ω = (16) This is accomplished by solving the discrete system 2 hφ n+1 = DMAC u + C 3 (17) where the MAC divergence operator is deined as ollows (D MAC u) i,j = u i+ 1 2,j u i 1 2,j + v i,j+ 1 v 2 i,j 1 2 x y Step 3: Update pressure and velocity ield according to, u n+1 = u tg MAC φ n+1 + C 4 (18) p n+1/2 = p n 1/2 + φ n+1 µ ( D MAC u ) + C 5 (19) 2 The operators h and 2 h are the standard three point central dierence operators and C i, i = 1,..., 5, are the correction terms which are only non-zero at the points near the interace and are calculated using generalized inite dierence ormulas o the type introduced in the previous section. This method requires solving two Helmholtz equations or u in Eq (12) and one Poisson equation or φ n+1 in Eq (17). Since the correction terms only aect the right-hand sides o the discrete systems, we can take advantage o the ast solvers rom Fishpack [18] to solve these equations. B. Singular orce evaluation Having solved or u n+1 at the grid points, we now compute the velocity at the interace. In our method, we use a set o control points to represent the interace. The velocity at the control points, U k, is interpolated rom the velocity at the grid points. Thus, we can write U k = U(X k ) = B(u n+1 ) (2) where B is the bilinear interpolation operator which includes the appropriate correction terms which are required to guarantee second order accuracy when the derivatives o the velocity are discontinuous. In summary, the equations that need to be solved in order to calculate u n+1 and U k, can be written symbolically as, Eq (12) Hu = C + B 1 Eq (17) Lφ n+1 = Du + B 2 Eq (18) u n+1 = u Gφ n+1 + B 3 Eq (2) U k = Mu n+1 + B 4 Eliminating u, φ n+1 and u n+1 rom the above equations, we can compute the velocity U k at the control points as ollows, U k = M ( H 1 C GL 1 DH 1 C ) + (M ( H 1 B 1 GH 1 B 1 GL 1 ) ) B 2 + B 3 + B4 For convenience, we can write the above equation as U k = U k + A (21) where U k is simply the velocity at the control points obtained by solving Eqs (12) (2) with =, given u n and p n 1/2. A is a 2N b 2N b matrix, where N b is the number o control points. The vector A is the velocity at the control points obtained by solving the ollowing equations u = µ 2 2 u, u Ω = (22) 2 φ n+1 = u, n φ n+1 Ω = (23) u n+1 = u t φ n+1 (24) A = B(u n+1 ) (25) with being the singular orce at the immersed boundary. Eq (21) can be used to determine the singular orce i we know the prescribed velocity U p at the immersed boundary. Thus, the singular orce at the control points can be computed by solving A = U p U k (26) The matrix A is computed once and stored. We solve Eqs (22) (25) 2N b times, i.e., one or each column. Each time, the orce strength is set to zero except or the entry in the column we want to calculate which is set to one. Once the matrix A has been calculated, only the right hand side U p U k, needs to be computed at each timestep. The resulting small system o equations (26) is then solved at each timestep or the singular orce. Finally, we solve Eqs (12) (19) to obtain u n+1 and p n+1/2.
4 Z U ield at t = 1 N N b E(u) order E(u) 2 order N N b E(p) order E(p) 2 order TABLE I: The grid reinement analysis or the rotational low problem with µ =.2, = x/4, at t = Y Y.5 1 (a) 1 Velocity ield at t = 1.5 X X (b) Fig. 1: Velocity ield at time t = 1 with a grid, µ =.2, = x/4. The immersed boundary rotates with angular velocity ω = 2. Fig. 1(a) Plot o the x component o velocity ield. Fig. 1(b) Plot o velocity vector ield. V. NUMERICAL RESULTS In this section we present the numerical results or three problems which involve immersed boundaries. A. Rotational low In this problem, the interace is a circle with radius r =.3 embedded in square domain [ 1, 1] [ 1, 1]. We prescribe the interace to rotate with angular velocity ω = 2. We set µ =.2 and consider the solution when t = 1. The velocity ield is shown in Figure 1. We carried out a grid reinement analysis, using a reerence grid o , to determine the order o convergence o the algorithm. The results in Table I show that the velocity is second order accurate and pressure is nearly second order accurate. B. Flow past a circular cylinder In this example, we simulate an unsteady low past a circular cylinder immersed in a rectangular domain Ω = [, 3] [, 1.5]. 1 The cylinder has a diameter d =.1 and its center is located at (1.6,.75). The luid density is ρ = 1. and the reestream velocity is set to unity, U = 1. The viscosity is determined by the Reynolds number, Re d. Simulations have been perormed at Re d = 2, 4, 8, 1, 2 and 3 on computational mesh. We use 4 control points to represent the circular cylinder. At the inlow boundary we speciy the velocity corresponding to the reestream velocity, and a homogeneous Neumann boundary condition is applied at the top, bottom and exit boundaries. The pressure is set to zero at the exit boundary. The pressure ield plots are shown in Figure 2. The results appear to be in good agreement with the other numerical simulations and experimental results. It is important to note that the matrix A, or a closed immersed boundary, is singular. This happens because the pressure inside the close boundary is not uniquely determined. We choose the pressure inside the cylinder such that there is no jump in pressure at one o the control points, i.e., the normal orce at that point is set to zero. Thereore, we can eliminate one column and row o the matrix A corresponding to that control point, thus making the problem solvable. C. Flow past a lat plate In this example, we simulate an unsteady low past a lat plate immersed in a rectangular domain Ω = [, 3] [.75,.75]. The lat plate whose length is L =.1 is oriented in the crosslow direction and located at x = 1.4. Simulations have been perormed at Re L = 2, 5, 1, 1 and 5 on a computational mesh. We use 8 control points to represent the lat plate. The same boundary conditions as or the low past a cylinder problem are applied in this problem. The pressure ield plots are shown in Figure 3. VI. CONCLUSION We have presented a ormally second order accurate immersed interace method or the solution o the incompressible Navier- Stokes equations in irregular domain. The implementation has been tested with three examples involving rotational low and unsteady lows over a circular cylinder and over a lat plate. Numerical experiments have shown that our method can handle problems with rigid boundaries. We plan to incorporate the current approach with our earlier works [16] or problems with deormable interaces to deal with deormable boundaries in irregular domain problems and contact problems.
5 (a) (b) (c) (d) (e) () Fig. 2: Pressure ield at t = 1, = x/4. (a) Re = 2, (b) Re = 4, (c) Re = 8, (d) Re = 1, (e) Re = 2, () Re = (a) (b) (c) (d) (e) Fig. 3: Pressure ield at t = 1, = x/4. (a) Re = 2, (b) Re = 5, (c) Re = 1, (d) Re = 1, (e) Re = 5.
6 REFERENCES [1] Brown DL, Cortez R, and Minion ML. Accurate projection methods or the incompressible Navier-Stokes equations, J. Comput. Phys. 168, 464 (21). [2] Kim J and Moin P. Application o a ractional step method to incompressible Navier-Stokes equations, J. Comput. Phys. 59, 38 (1985). [3] LeVeque RJ and Li Z. The immersed interace method or elliptic equations with discontinuous coeicients and singular sources, SIAM J. Numer. Anal. 31, 119 (1994). [4] LeVeque RJ and Li Z. Immersed interace method or Stokes low with elastic boundaries or surace tension, SIAM J. Sci. Comput. 18, 79 (1997). [5] Wiegmann A and Bube KP. The explicit-jump immersed interace method: Finite dierence methods or PDEs with piecewise smooth solutions, SIAM J. Numer. Anal. 37, No. 3, 827 (1997). [6] Li and Lai MC. The immersed interace method or the Navier-Stokes equations with singular orces, J. Comput. Phys. 171, 822 (21). [7] Li Z and Wang C. A ast inite dierence method or solving Navier- Stokes equations on irregular domains, Comm. Math. Sci. Vol 1, No 1, pp [8] Calhoun D. A Cartesian grid method or solving the two-dimensional streamunction-vorticity equations in irregular regions, J. Comput. Phys. 176, (22). [9] Lee L. An immersed interace method or the incompressible Navier- Stokes equations, Ph.D. thesis (University o Washington, 22). [1] Peskin CS. The immersed boundary method, Acta Numerica (22), pp [11] Lai MC and Peskin CS. An immersed boundary method with ormal second order accuracy and reduced numerical viscosity, J. Comput. Phys. 16, (2). [12] Goldstein D, Handler R and Sirovich L. Modeling a no-slip low with an external orce ield, J. Comput. Phys. 15, 354 (1993). [13] Ye T, Mittal R, Udaykumar HS and Shyy W. An accurate Cartesian grid method or viscous incompressible lows with complex immersed boundary, J. Comput. Phys. 156, (1999). [14] Udaykumar HS, Mittal R, Rampunggoon P and Khanna A. A sharp interace Cartesian grid method or simulating lows with complex moving boundaries, J. Comput. Phys. 174, (21). [15] Lima E Silva ALF, Silveira-Neto A and Damasceno JJR. Numerical simulation o two-dimensional lows over a circular cylinder using the immersed boundary method, J. Comput. Phys. 189, (23). [16] Le DV, Khoo BC, and Peraire J. An immersed interace method or the incompressible Navier-Stokes equations, presented at the SMA Symposium, Singapore 24. [17] Shu CW and Osher S. Eicient implementation o essentially nonoscillatory shock capturing scheme, II, J. Comput. Phys. 83, 32 (1989). [18] Adams J, Swarztrauber P and Sweet R. FISHPACK: Eicient FOR- TRAN Subprograms or the Solution o Separable Eliptic Partial Dierential Equations, National Center or Athmospheric Research, 1999,
Introduction to immersed boundary method
Introduction to immersed boundary method Ming-Chih Lai mclai@math.nctu.edu.tw Department of Applied Mathematics Center of Mathematical Modeling and Scientific Computing National Chiao Tung University 1001,
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 informationANALYSIS AND NUMERICAL METHODS FOR SOME CRACK PROBLEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING, SERIES B Volume 2, Number 2-3, Pages 155 166 c 2011 Institute for Scientific Computing and Information ANALYSIS AND NUMERICAL METHODS FOR SOME
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 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 informationExplicit Jump Immersed Interface Method: Documentation for 2D Poisson Code
Eplicit Jump Immersed Interface Method: Documentation for 2D Poisson Code V. Rutka A. Wiegmann November 25, 2005 Abstract The Eplicit Jump Immersed Interface method is a powerful tool to solve elliptic
More information2 Z. Li and C. Wang where u = (u; v) is the velocity, p is the pressure, ν is the viscosity. We assume that the boundary of the is piecewise
A fast finite difference method for solving Navier-Stokes Equations on irregular domains Zhilin Li Λ Cheng Wang y April 17, 2002 keywords: Navier-Stokes equations, irregular domains, vorticity stream-function
More informationc 2003 International Press
COMMUNICATIONS IN MATH SCI Vol. 1, No. 1, pp. 181 197 c 2003 International Press A FAST FINITE DIFFERENCE METHOD FOR SOLVING NAVIER-STOKES EQUATIONS ON IRREGULAR DOMAINS ZHILIN LI AND CHENG WANG Abstract.
More informationA Cartesian adaptive level set method for two-phase flows
Center or Turbulence Research Annual Research Bries 2003 227 A Cartesian adaptive level set method or two-phase lows By F. Ham AND Y.-N. Young 1. Motivation and objectives Simulations o lows involving
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 informationA finite difference Poisson solver for irregular geometries
ANZIAM J. 45 (E) ppc713 C728, 2004 C713 A finite difference Poisson solver for irregular geometries Z. Jomaa C. Macaskill (Received 8 August 2003, revised 21 January 2004) Abstract The motivation for this
More informationA Simple Compact Fourth-Order Poisson Solver on Polar Geometry
Journal of Computational Physics 182, 337 345 (2002) doi:10.1006/jcph.2002.7172 A Simple Compact Fourth-Order Poisson Solver on Polar Geometry Ming-Chih Lai Department of Applied Mathematics, National
More informationPhysics 5153 Classical Mechanics. Solution by Quadrature-1
October 14, 003 11:47:49 1 Introduction Physics 5153 Classical Mechanics Solution by Quadrature In the previous lectures, we have reduced the number o eective degrees o reedom that are needed to solve
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 informationAN AUGMENTED IIM-LEVEL SET METHOD FOR STOKES EQUATIONS WITH DISCONTINUOUS VISCOSITY
Sixth Mississippi State Conference on Differential Equations and Computational Simulations, Electronic Journal of Differential Equations, Conference 15 (2007), pp. 193 210. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu
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 informationAn Immersed Boundary Method for Restricted Diffusion with Permeable Interfaces
An Immersed Boundary Method for Restricted Diffusion with Permeable Interfaces Huaxiong Huang Kazuyasu Sugiyama Shu Takagi April 6, 009 Keywords: Restricted diffusion; Permeable interface; Immersed boundary
More informationIterative Methods for Stokes/Darcy Coupling
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
More information(1:1) 1. The gauge formulation of the Navier-Stokes equation We start with the incompressible Navier-Stokes equation 8 >< >: u t +(u r)u + rp = 1 Re 4
Gauge Finite Element Method for Incompressible Flows Weinan E 1 Courant Institute of Mathematical Sciences New York, NY 10012 Jian-Guo Liu 2 Temple University Philadelphia, PA 19122 Abstract: We present
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 information1 Relative degree and local normal forms
THE ZERO DYNAMICS OF A NONLINEAR SYSTEM 1 Relative degree and local normal orms The purpose o this Section is to show how single-input single-output nonlinear systems can be locally given, by means o a
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
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 informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
More informationSolving Partial Differential Equations Numerically. Miklós Bergou with: Gary Miller, David Cardoze, Todd Phillips, Mark Olah
Solving Partial Dierential Equations Numerically Miklós Bergou with: Gary Miller, David Cardoze, Todd Phillips, Mark Olah Overview What are partial dierential equations? How do we solve them? (Example)
More informationROAD MAP... D-1: Aerodynamics of 3-D Wings D-2: Boundary Layer and Viscous Effects D-3: XFLR (Aerodynamics Analysis Tool)
AE301 Aerodynamics I UNIT D: Applied Aerodynamics ROAD MAP... D-1: Aerodynamics o 3-D Wings D-2: Boundary Layer and Viscous Eects D-3: XFLR (Aerodynamics Analysis Tool) AE301 Aerodynamics I : List o Subjects
More informationABSTRACT. coefficients. The main motivation is to get not only a second order accurate solution but
ABSTRACT CHEN, GUANYU. Accurate Gradient Computation for Elliptic Interface Problems with Discontinuous and Variable Coefficients. (Under the direction of Dr. Zhilin Li.) A new numerical method is proposed
More informationA Cartesian Grid Method for Solving the Two-Dimensional Streamfunction-Vorticity Equations in Irregular Regions
Journal of Computational Physics 176, 231 275 (2002) doi:10.1006/jcph.2001.6970, available online at http://www.idealibrary.com on A Cartesian Grid Method for Solving the Two-Dimensional Streamfunction-Vorticity
More informationExperience with DNS of particulate flow using a variant of the immersed boundary method
Experience with DNS of particulate flow using a variant of the immersed boundary method Markus Uhlmann Numerical Simulation and Modeling Unit CIEMAT Madrid, Spain ECCOMAS CFD 2006 Motivation wide range
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 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 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 informationNew Finite Difference Methods Based on IIM for Inextensible Interfaces in Incompressible Flows
East Asian Journal of Applied Mathematics Vol. x, No. x, pp. -7 doi:.48/eajam xxx x New Finite Difference Methods Based on IIM for Inextensible Interfaces in Incompressible Flows Zhilin Li, and Ming-Chih
More informationA Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
More informationAn Immersed Boundary Method with Formal Second-Order Accuracy and Reduced Numerical Viscosity 1
Journal of Computational Physics 160, 705 719 (2000) doi:10.1006/jcph.2000.6483, available online at http://www.idealibrary.com on An Immersed Boundary Method with Formal Second-Order Accuracy and Reduced
More informationRunge-Kutta-Chebyshev Projection Method
Runge-Kutta-Chebyshev Projection Method Zheming Zheng Linda Petzold Department of Mechanical Engineering, University of California Santa Barbara, Santa Barbara, CA 93106, USA June 8, 2006 This work was
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 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 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 informationNUMERICAL COMPUTATION OF A VISCOUS FLOW AROUND A CIRCULAR CYLINDER ON A CARTESIAN GRID
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2000 Barcelona, 11-14 September 2000 c ECCOMAS NUMERICAL COMPUTATION OF A VISCOUS FLOW AROUND A CIRCULAR CYLINDER
More informationGauge finite element method for incompressible flows
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2000; 34: 701 710 Gauge finite element method for incompressible flows Weinan E a, *,1 and Jian-Guo Liu b,2 a Courant Institute
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 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 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 informationLecture : Feedback Linearization
ecture : Feedbac inearization Niola Misovic, dipl ing and Pro Zoran Vuic June 29 Summary: This document ollows the lectures on eedbac linearization tought at the University o Zagreb, Faculty o Electrical
More informationComputational Methods for Domains with! Complex Boundaries-I!
http://www.nd.edu/~gtrggva/cfd-course/ Computational Methods or Domains with Comple Boundaries-I Grétar Trggvason Spring For most engineering problems it is necessar to deal with comple geometries, consisting
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
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 informationA Momentum Exchange-based Immersed Boundary-Lattice. Boltzmann Method for Fluid Structure Interaction
APCOM & ISCM -4 th December, 03, Singapore A Momentum Exchange-based Immersed Boundary-Lattice Boltzmann Method for Fluid Structure Interaction Jianfei Yang,,3, Zhengdao Wang,,3, and *Yuehong Qian,,3,4
More informationMathematical Notation Math Calculus & Analytic Geometry III
Mathematical Notation Math 221 - alculus & Analytic Geometry III Use Word or WordPerect to recreate the ollowing documents. Each article is worth 10 points and should be emailed to the instructor at james@richland.edu.
More informationEfficient Fluid-Solid Coupling Algorithm for Single-Phase Flows
Efficient Fluid-Solid Coupling Algorithm for Single-Phase Flows Yen Ting Ng Chohong Min Frédéric Gibou April 11, 009 Abstract We present a simple and efficient fluid/solid coupling algorithm in two-spatial
More informationSecond Order Slip Flow of Cu-Water Nanofluid Over a Stretching Sheet With Heat Transfer
Second Order Slip Flow o Cu-Water Nanoluid Over a Stretching Sheet With Heat Transer RAJESH SHARMA AND ANUAR ISHAK School o Mathematical Sciences, Faculty o Science and Technology Universiti Kebangsaan
More informationThe Embedded Boundary Integral Method (EBI) for the Incompressible Navier-Stokes equations
The Embedded Boundary Integral Method (EBI) for the Incompressible Navier-Stokes equations George Biros Lexing Ying and Denis Zorin Courant Institute of Mathematical Sciences, New York University, NY 10012,
More informationMathematical Notation Math Calculus & Analytic Geometry III
Name : Mathematical Notation Math 221 - alculus & Analytic Geometry III Use Word or WordPerect to recreate the ollowing documents. Each article is worth 10 points and can e printed and given to the instructor
More informationOn the Solution of the Elliptic Interface Problems by Difference Potentials Method
On the Solution of the Elliptic Interface Problems by Difference Potentials Method Yekaterina Epshteyn and Michael Medvinsky Abstract Designing numerical methods with high-order accuracy for problems in
More informationOPTIMALLY STAGGERED FINNED CIRCULAR AND ELLIPTIC TUBES IN FORCED CONVECTION
OPTIMALLY STAGGERED FINNED CIRCULAR AND ELLIPTIC TUBES IN FORCED CONVECTION R. S. Matos a, T. A. Laursen b, J. V. C. Vargas a, and A. Bejan c, a Universidade Federal do Paraná Departamento de Engenharia
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 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 informationSpatial discretization scheme for incompressible viscous flows
Spatial discretization scheme for incompressible viscous flows N. Kumar Supervisors: J.H.M. ten Thije Boonkkamp and B. Koren CASA-day 2015 1/29 Challenges in CFD Accuracy a primary concern with all CFD
More informationFeedback Linearization Lectures delivered at IIT-Kanpur, TEQIP program, September 2016.
Feedback Linearization Lectures delivered at IIT-Kanpur, TEQIP program, September 216 Ravi N Banavar banavar@iitbacin September 24, 216 These notes are based on my readings o the two books Nonlinear Control
More informationFlow Field and Oscillation Frequency of a Rotating Liquid Droplet
Flow Field and Oscillation Frequency of a Rotating Liquid Droplet TADASHI WATANABE Center for Computational Science and e-systems Japan Atomic Energy Agency (JAEA) Tokai-mura, Naka-gun, Ibaraki-ken, 319-1195
More informationarxiv: v2 [math.na] 29 Sep 2011
A Correction Function Method for Poisson Problems with Interface Jump Conditions. Alexandre Noll Marques a, Jean-Christophe Nave b, Rodolfo Ruben Rosales c arxiv:1010.0652v2 [math.na] 29 Sep 2011 a Department
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
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 informationOn the Development of Implicit Solvers for Time-Dependent Systems
School o something FACULTY School OF OTHER o Computing On the Development o Implicit Solvers or Time-Dependent Systems Peter Jimack School o Computing, University o Leeds In collaboration with: P.H. Gaskell,
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 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 informationManoj Kumar 1, Pratibha Joshi 2. (Received 31 August 2011, accepted 26 December 2011)
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.14(2012) No.1,pp.11-22 A Mathematical Model and Numerical Solution of a One Dimensional Steady State Heat Conduction
More informationSyllabus Objective: 2.9 The student will sketch the graph of a polynomial, radical, or rational function.
Precalculus Notes: Unit Polynomial Functions Syllabus Objective:.9 The student will sketch the graph o a polynomial, radical, or rational unction. Polynomial Function: a unction that can be written in
More informationA REPORT ON PERFORMANCE OF ANNULAR FINS HAVING VARYING THICKNESS
VOL., NO. 8, APRIL 6 ISSN 89-668 ARPN Journal o Engineering and Applied Sciences 6-6 Asian Research Publishing Networ (ARPN). All rights reserved. A REPORT ON PERFORMANCE OF ANNULAR FINS HAVING VARYING
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 informationThe Relationship between Discrete Calculus Methods and other Mimetic Approaches
The Relationship between Discrete Calculus Methods and other Mimetic Approaches Blair Perot Theoretical and Computational Fluid Dynamics Laboratory Woudschoten Conerence Oct. 2011 October 7, 2011 Background
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 informationFundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics
Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/
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 informationNumerical Methods of Applied Mathematics -- II Spring 2009
MIT OpenCourseWare http://ocw.mit.edu 18.336 Numerical Methods of Applied Mathematics -- II Spring 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Discretization of Boundary Conditions Discretization of Boundary Conditions On
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 informationLevel Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows
Level Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows Abstract Maged Ismail Claremont Graduate University Level Set and Phase Field methods are well-known interface-capturing
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 informationTransport Properties: Momentum Transport, Viscosity
Transport Properties: Momentum Transport, Viscosity 13th February 2011 1 Introduction Much as mass(material) is transported within luids (gases and liquids), linear momentum is also associated with transport,
More informationIntroduction to Simulation - Lecture 2. Equation Formulation Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy
Introduction to Simulation - Lecture Equation Formulation Methods Jacob White Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Outline Formulating Equations rom Schematics Struts and Joints
More informationBoundary-Fitted Coordinates!
Computational Fluid Dnamics http:wwwndedu~gtrggvacfdcourse Computational Fluid Dnamics Computational Methods or Domains with Comple BoundariesI Grétar Trggvason Spring For most engineering problems it
More informationDUE: WEDS MARCH 26TH 2018
HOMEWORK # 2: FINITE DIFFERENCES MAPPING AND TWO DIMENSIONAL PROBLEMS DUE: WEDS MARCH 26TH 2018 NOTE: In this homework, you will choose ONE of the following three questions to perform and hand-in. Each
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 informationApplied Numerical Mathematics. High-order numerical schemes based on difference potentials for 2D elliptic problems with material interfaces
Applied Numerical Mathematics 111 (2017) 64 91 Contents lists available at ScienceDirect Applied Numerical Mathematics www.elsevier.com/locate/apnum High-order numerical schemes based on difference potentials
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 informationFeedback Linearization
Feedback Linearization Peter Al Hokayem and Eduardo Gallestey May 14, 2015 1 Introduction Consider a class o single-input-single-output (SISO) nonlinear systems o the orm ẋ = (x) + g(x)u (1) y = h(x) (2)
More informationPartially fluidized shear granular flows: Continuum theory and molecular dynamics simulations
Partially luidized shear granular lows: Continuum theory and molecular dynamics simulations Dmitri Volson, 1 Lev S. Tsimring, 1 and Igor S. Aranson 2 1 Institute or Nonlinear Science, University o Caliornia,
More informationAn adaptive mesh redistribution method for nonlinear Hamilton-Jacobi equations in two- and three-dimensions
An adaptive mesh redistribution method for nonlinear Hamilton-Jacobi equations in two- and three-dimensions H.-Z. Tang, Tao Tang and Pingwen Zhang School of Mathematical Sciences, Peking University, Beijing
More informationChapter 3 Water Flow in Pipes
The Islamic University o Gaza Faculty o Engineering Civil Engineering Department Hydraulics - ECI 33 Chapter 3 Water Flow in Pipes 3. Description o A Pipe Flow Water pipes in our homes and the distribution
More informationScattered Data Approximation of Noisy Data via Iterated Moving Least Squares
Scattered Data Approximation o Noisy Data via Iterated Moving Least Squares Gregory E. Fasshauer and Jack G. Zhang Abstract. In this paper we ocus on two methods or multivariate approximation problems
More informationAN OVERVIEW OF THE IMMERSED INTERFACE METHOD AND ITS APPLICATIONS
. AN OVERVIEW OF THE IMMERSED INTERFACE METHOD AND ITS APPLICATIONS ZHILIN LI Abstract. Interface problems have many applications. Mathematically, interface problems usually lead to differential equations
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 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 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 information1 Introduction. J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods
J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods 1 Introduction Achieving high order time-accuracy in the approximation of the incompressible Navier Stokes equations by means of fractional-step
More informationDefinition: Let f(x) be a function of one variable with continuous derivatives of all orders at a the point x 0, then the series.
2.4 Local properties o unctions o several variables In this section we will learn how to address three kinds o problems which are o great importance in the ield o applied mathematics: how to obtain the
More informationAn immersed boundary method for interfacial flows with insoluble surfactant
An immersed boundary method for interfacial flows with insoluble surfactant Ming-Chih Lai Yu-Hau Tseng Huaxiong Huang Abstract In this paper, an immersed boundary method is proposed for the simulation
More informationTRANSIENT SIMULATION OF LIQUID ROCKET ENGINES: A STEP TOWARDS A MORE EDUCATED PROPELLANT CHOICE BETWEEN KEROSENE AND METHANE.
TRANSIENT SIMULATION OF LIQUID ROCKET ENGINES: A STEP TOWARDS A MORE EDUCATED PROPELLANT CHOICE BETWEEN KEROSENE AND METHANE Chiara Manletti Space Launcher Systems Analysis (SART), DLR, Cologne, Germany,
More information