Peaceman and Thiem well models or how to remove a logarithmic singularity fr. your numerical solution
|
|
- Delphia Wilcox
- 6 years ago
- Views:
Transcription
1 Peaceman and Thiem well models or how to remove a logarithmic singularity from your numerical solution Department of Mathematics Oregon State University AMC, 3/2/2007 Current work supported by NSF and DOE Acknowledgement: Cristiano Garibotti
2 Outline Problem statement: what are wells 1 Problem statement: what are wells Steady state, rate specified wells Idealized and real well models 2 Not enough regularity for classical FE... Why subtracting singularities is not sufficient 3 Effective well radius idea: Peaceman correction Beyond Peaceman correction 4 Closeness of solutions: elliptic and parabolic estimates
3 What are wells Problem statement: what are wells Steady state, rate specified wells Idealized and real well models
4 Steady state, rate specified wells Idealized and real well models Well in a confined aquifer/hydrocarbon reservoir
5 Steady state, rate specified wells Idealized and real well models Pressure profiles around the well: side view
6 Steady state, rate specified wells Idealized and real well models Pressure profiles around the well: aerial view
7 Mathematical model: geometry Steady state, rate specified wells Idealized and real well models
8 Steady state, rate specified wells Idealized and real well models Problem statement for steady state, rate specified, R 2 Conservation equations: u = K p, x Ω Darcy slaw u = q, conservation of mass q(x) represents source/sink terms, and is known, K is the conductivity tensor and is known. Put it together K p = q(x), x Ω p Ω = p GIVEN = 0
9 Steady state, rate specified wells Idealized and real well models Mathematical model for wells: idealized and real Geometry Ω w B(x 0, r w ) with boundary Γ w 1 (R)Real model of injection/production wells K p = 0, x Ω \ Ω w K p ν = q, x Ω w 2πr w 2 (ID) Idealization of a true real-life situation (K p δ ) = q(x) = qδ x0 (x), x Ω where by δ x0 (x) we mean the Dirac-δ distribution δ Rx0 (x) = 0, x x 0 D δx 0 (x) = 1, and where q (total rate) is given as data.
10 Steady state, rate specified wells Idealized and real well models Solution to the mathematically idealized problem (ID) Analytical solution in 2D, steady state, with one well at 0, single-phase flow with homogeneous isotropic K = K I K ( p δ ) = qδ 0, x Ω Change variable to polar coordinates p δ = p δ (r, θ), assume radial solution p δ = p δ (r): 1 r pδ (Kr r r ) = qδ 0, r > 0 Solution has a logarithmic singularity p δ (r) = C 1 log(r) + C 2, r > 0 Fix the constants (use total input q and boundary conditions on Ω C 1 = q 2πK, C 2 = p GIVEN (easiest case C 2 = 0) to get p δ (r) = q 2πK log(r)
11 Steady state, rate specified wells Idealized and real well models Analytical solution to (ID) is not very useful... when (bottom hole) pressure in the well is needed Problem!!!! p δ (x 0 ) = when multiple wells are present: K p = i q iδ xi, x Ω
12 Steady state, rate specified wells Idealized and real well models Solution to the real single well problem (R) Known: K, r w, q. Solve K p = 0, x Ω \ Ω w K p ν = q, x Ω w 2πr w Similarly as in (ID) model we can derive p(r) p(r) = q 2πK log( r R ) where p(r) is given at some R > 0. For example, p(r) can be given at some distance away from wells (Thiem solution) [Guenther,Lee 8.1/#12,Marsilly 86] Analytical solution still not useful when multiple wells are present
13 Steady state, rate specified wells Idealized and real well models Try numerical solution for (R) real well model... Very fine grid around (multiple) wells would be necessary. Still no good way to get p(r w )
14 Not enough regularity for classical FE... Why subtracting singularities is not sufficient Numerical handling of idealized model using FE Functional spaces... notation W m,p (Ω) := {u : m u L p (Ω)} W 0,2 (Ω) L 2 (Ω) W m,2 (Ω) H m (Ω) Assume smooth q for the moment... derive weak form of (K p) = q K p w = qw, w H0 1 (Ω) Ω Ω (Classical FE Galerkin method:) Find p h V h which approximates p and solves: K p h w h = qw h, w h V h H0 1 (Ω) Ω Ω Error estimate p p h H1 (Ω) Ch p H2 Ω But q is not smooth... and p is not smooth...
15 Not enough regularity for classical FE... Why subtracting singularities is not sufficient How nonsmooth are solutions to (ID)...? Elliptic theory[giltru,dautray-lions] tells us that weak solution to (K p) = q has the following regularity q L 2 (Ω) = W 0,2 (Ω) = p H 2 (Ω) q H 1 (Ω) = W 1,2 (Ω) = p H 1 (Ω) q L 1 (Ω) = W 0,1 (Ω) W 1,p (Ω) = p W 1,p (Ω) (p < 2). But δ x0 L p (Ω)!! Note: Ω δ x 0 w = w(x 0 ) requires w C 0 (Ω) or that δ 0 (C 0 (Ω)) = M(Ω) (dual space to C 0 (Ω)). However, we have W 1,p (Ω) C 0 ( Ω) so q (C 0 ( Ω)) M(Ω) W 1,p (Ω) if p < 2 (Sobolev imbedding Thm) hence the solution p W 1,p (Ω) for any p < 2. This is not enough regularity for FE solution to converge
16 Not enough regularity for classical FE... Why subtracting singularities is not sufficient Classical FE solution may not converge. What now? We do not have p H 1 (Ω), hence, p h may fail to converge to p. Idea: subtract the singular part of the solution Consider a general problem K p = i q iδ xi, x Ω Construct p S (x) = i q i 2πK i log(r i (x)) and instead of p approximate p REGULAR := p p S, that is, solve K ph REGULAR w h = K p S νw h, w h V h Ω i Ω and so p REGULAR H 2 (Ω) (K homogeneous). If K i are different, then one can show p REGULAR H 2 ɛ (Ω). and almost optimal convergence rates for p REGULAR ph REGULAR are obtained. [Ewing 82,Wheeler/Russell 82]
17 Not enough regularity for classical FE... Why subtracting singularities is not sufficient Why subtracting singularities is not sufficient The above idea (solve for (p REGULAR ) h, add p S ) works well at a sufficiently large distance from wells. However, at wells it is STILL designed to give p(x i ) =. good for rate-specified wells (when q is given) when p(x i ) is not needed not applicable with pressure-specified wells (p(x i ) is given) Idea: consider the idea of an effective well radius... where we exploit the familiar equation from (R) real well model p(r) p(r w ) = q 2πK log( r ) r w
18 Knowing q, how to get p w? Effective well radius idea: Peaceman correction Beyond Peaceman correction Given q, we can hypothetically compute a numerical solution p h p = (p 0, p 1,... p N ) where p 0 p(x 0 ) we have a numerical method in which we solve for p h using q Ap = f but clearly p w p 0!!! but... we have an analytical formula p(r) p(r w ) = q 2πK log( r r w ) we hope this formula applies already to p 1,2,3,4 Idea: combine these two methods to get an expression for p w
19 Effective well radius idea: Peaceman correction Beyond Peaceman correction Peaceman correction [Peaceman 77] (homogeneous isotropic K, isotropic grid) Cell Centered Finite Difference Mixed Finite Element with RT0 1 1 h q h 2 ( K p 1 p 0 h + K p 2 p 0 h + K p 3 p 0 h + K p 4 p 0 ) h = 2 p 1,2,3,4 p 0 == q 2πK log( h r w ) Combining 1,2 we get p 0 = p w + q 2πK (log( h ) π r w 2 ) = p w + q 2πK (log( r 0 ) r w where r 0 = exp( π/2)h 0.208h is the effective well radius.
20 Effective well radius idea: Peaceman correction Beyond Peaceman correction But what if K is not isotropic homogeneous? K = K T (permeability, conductivity, mobility,...) is in general anisotropic and heterogeneous. And grids are not uniform and isotropic... Peaceman models allow for K = diag(k 11, K 22 ), h x h y homogenization for Darcy flow around wells [Zijl, Trykozko 01] multiscale FE method for Darcy flow near wells[z. Chen et al 03] extension to non-darcy flow, K const [Lazarov, Ewing et al 99] CG & MP... in progress
21 Effective well radius idea: Peaceman correction Beyond Peaceman correction Results: Darcy with rate-wells and BHP wells
22 Idealized model versus real model Closeness of solutions: elliptic and parabolic estimates But.. good in many situations when diam(ω) >>>>>>> diam(ω w ) How close is p δ to the true p? What if there are more wells and flow is not single phase but multiphase? use numerical solution: but how to guarantee convergence? We need to know p(r w ) (pressure in well)! in order to determine phase behavior What if K is not constant and not isotropic?
23 How close is p δ to p? Closeness of solutions: elliptic and parabolic estimates For non-stationary version of the model, as r w 0, one has (if coefficients are extended properly) [G. Chavent, J. Jaffré 86] p δ p L2 (Ω (0,T )) 0 For the original stationary problem: [Li Ta Tsien, A. Damlamian 80], [Da Qian Li, Shu Xing Chen (in Chinese) 78], Recent result [Zhiming Chen, Xinye Yue 03] max x Ω pδ p Cr w log(r w )
24 Back to recent estimate Closeness of solutions: elliptic and parabolic estimates [Zhiming Chen, Xinye Yue 03] Assume: K C 0,1 (Ω) with Lipschitz constant Λ Thus for any p > 0 the solution satisfies (p δ p) W 2,p (Ω) = (p δ p) C 1 (Ω) and we have the following estimate max x Ω pδ p CΛ(1 + log(λ)) 1/2 r w log(r w )
25 Closeness of solutions: elliptic and parabolic estimates
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part IV Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
More informationHomogenization and Multiscale Modeling
Ralph E. Showalter http://www.math.oregonstate.edu/people/view/show Department of Mathematics Oregon State University Multiscale Summer School, August, 2008 DOE 98089 Modeling, Analysis, and Simulation
More informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part II Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
More informationWeak Galerkin Finite Element Methods and Applications
Weak Galerkin Finite Element Methods and Applications Lin Mu mul1@ornl.gov Computational and Applied Mathematics Computationa Science and Mathematics Division Oak Ridge National Laboratory Georgia Institute
More informationParallel Simulation of Subsurface Fluid Flow
Parallel Simulation of Subsurface Fluid Flow Scientific Achievement A new mortar domain decomposition method was devised to compute accurate velocities of underground fluids efficiently using massively
More informationGENERALIZED PSEUDOPRESSURE WELL TREATMENT
GENERALIZED PSEUDOPRESSURE WELL TREATMENT IN RESERVOIR SIMULATION Curtis H. Whitson a,b Øivind Fevang b a Norwegian University of Science and Technology (NTNU) b PERA a/s ABSTRACT This paper presents a
More informationIntroduction to Well Hydraulics Fritz R. Fiedler
Introduction to Well Hydraulics Fritz R. Fiedler A well is a pipe placed in a drilled hole that has slots (screen) cut into it that allow water to enter the well, but keep the aquifer material out. A well
More informationMAT389 Fall 2016, Problem Set 4
MAT389 Fall 2016, Problem Set 4 Harmonic conjugates 4.1 Check that each of the functions u(x, y) below is harmonic at every (x, y) R 2, and find the unique harmonic conjugate, v(x, y), satisfying v(0,
More informationXingye Yue. Soochow University, Suzhou, China.
Relations between the multiscale methods p. 1/30 Relationships beteween multiscale methods for elliptic homogenization problems Xingye Yue Soochow University, Suzhou, China xyyue@suda.edu.cn Relations
More informationAnalysis of Two-Grid Methods for Nonlinear Parabolic Equations by Expanded Mixed Finite Element Methods
Advances in Applied athematics and echanics Adv. Appl. ath. ech., Vol. 1, No. 6, pp. 830-844 DOI: 10.408/aamm.09-m09S09 December 009 Analysis of Two-Grid ethods for Nonlinear Parabolic Equations by Expanded
More informationSeparation of Variables in Linear PDE: One-Dimensional Problems
Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,
More informationMonotonicity Conditions for Discretization of Parabolic Conservation Laws. Hilde Kristine Hvidevold
Monotonicity Conditions for Discretization of Parabolic Conservation Laws Master of Science Thesis in Applied Mathematics Hilde Kristine Hvidevold Department of Mathematics University of Bergen June 2,
More informationA generalised Ladyzhenskaya inequality and a coupled parabolic-elliptic problem
A generalised Ladyzhenskaya inequality and a coupled parabolic-elliptic problem Dave McCormick joint work with James Robinson and José Rodrigo Mathematics and Statistics Centre for Doctoral Training University
More informationIn this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,
Chapter 8 Elliptic PDEs In this chapter we study elliptical PDEs. That is, PDEs of the form 2 u = lots, where lots means lower-order terms (u x, u y,..., u, f). Here are some ways to think about the physical
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationContinuous dependence estimates for the ergodic problem with an application to homogenization
Continuous dependence estimates for the ergodic problem with an application to homogenization Claudio Marchi Bayreuth, September 12 th, 2013 C. Marchi (Università di Padova) Continuous dependence Bayreuth,
More informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
More informationNumerical Simulation of Single-Phase and Multiphase Non-Darcy Flow in Porous and Fractured Reservoirs
Transport in Porous Media 49: 209 240, 2002. 2002 Kluwer Academic Publishers. Printed in the Netherlands. 209 Numerical Simulation of Single-Phase and Multiphase Non-Darcy Flow in Porous and Fractured
More informationHomogenization and numerical Upscaling. Unsaturated flow and two-phase flow
Homogenization and numerical Upscaling Unsaturated flow and two-phase flow Insa Neuweiler Institute of Hydromechanics, University of Stuttgart Outline Block 1: Introduction and Repetition Homogenization
More informationA NUMERICAL APPROXIMATION OF NONFICKIAN FLOWS WITH MIXING LENGTH GROWTH IN POROUS MEDIA. 1. Introduction
Acta Math. Univ. Comenianae Vol. LXX, 1(21, pp. 75 84 Proceedings of Algoritmy 2 75 A NUMERICAL APPROXIMATION OF NONFICIAN FLOWS WITH MIXING LENGTH GROWTH IN POROUS MEDIA R. E. EWING, Y. LIN and J. WANG
More informationOperator Upscaling for the Wave Equation
Operator Upscaling for the Wave Equation Tetyana Vdovina Susan E. Minkoff UMBC), Oksana Korostyshevskaya Department of Computational and Applied Mathematics Rice University, Houston TX vdovina@caam.rice.edu
More informationPetroleum Engineering 613 Natural Gas Engineering. Texas A&M University. Lecture 07: Wellbore Phenomena
Petroleum Engineering 613 Natural Gas Engineering Texas A&M University Lecture 07: T.A. Blasingame, Texas A&M U. Department of Petroleum Engineering Texas A&M University College Station, TX 77843-3116
More informationPreconditioned space-time boundary element methods for the heat equation
W I S S E N T E C H N I K L E I D E N S C H A F T Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
More informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More informationLocal Time Step for a Finite Volume Scheme I.Faille F.Nataf*, F.Willien, S.Wolf**
Controlled CO 2 Diversified fuels Fuel-efficient vehicles Clean refining Extended reserves Local Time Step for a Finite Volume Scheme I.Faille F.Nataf*, F.Willien, S.Wolf** *: Laboratoire J.L.Lions **:Université
More informationResearch Article A Stabilized Low Order Finite-Volume Method for the Three-Dimensional Stationary Navier-Stokes Equations
Mathematical Problems in Engineering Volume 212, Article ID 297269, 14 pages doi:1.1155/212/297269 Research Article A Stabilized Low Order Finite-Volume Method for the Three-Dimensional Stationary Navier-Stokes
More informationE&M. 1 Capacitors. January 2009
E&M January 2009 1 Capacitors Consider a spherical capacitor which has the space between its plates filled with a dielectric of permittivity ɛ. The inner sphere has radius r 1 and the outer sphere has
More informationMATH Final Project Mean Curvature Flows
MATH 581 - Final Project Mean Curvature Flows Olivier Mercier April 30, 2012 1 Introduction The mean curvature flow is part of the bigger family of geometric flows, which are flows on a manifold associated
More informationApplications of Partial Differential Equations in Reservoir Simulation
P-32 Applications of Partial Differential Equations in Reservoir Simulation Deepak Singh Summary The solution to stochastic partial differential equations may be viewed in several manners. One can view
More informationComm. Nonlin. Sci. and Numer. Simul., 12, (2007),
Comm. Nonlin. Sci. and Numer. Simul., 12, (2007), 1390-1394. 1 A Schrödinger singular perturbation problem A.G. Ramm Mathematics Department, Kansas State University, Manhattan, KS 66506-2602, USA ramm@math.ksu.edu
More information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More informationESTIMATES FOR THE MONGE-AMPERE EQUATION
GLOBAL W 2,p ESTIMATES FOR THE MONGE-AMPERE EQUATION O. SAVIN Abstract. We use a localization property of boundary sections for solutions to the Monge-Ampere equation obtain global W 2,p estimates under
More informationRadius of single fluxon electron model identical with classical electron radius
Radius of single fluxon electron model identical with classical electron radius U. Kayser-Herold 12 February 2017 Abstract Analytical determination of the magnetic flux included in the electron s dipole
More informationA new meshless method for steady-state heat conduction problems in anisotropic and inhomogeneous media
Archive of Applied Mechanics 74 25 563--579 Springer-Verlag 25 DOI 1.17/s419-5-375-8 A new meshless method for steady-state heat conduction problems in anisotropic and inhomogeneous media H. Wang, Q.-H.
More informationAn adaptive RBF-based Semi-Lagrangian scheme for HJB equations
An adaptive RBF-based Semi-Lagrangian scheme for HJB equations Roberto Ferretti Università degli Studi Roma Tre Dipartimento di Matematica e Fisica Numerical methods for Hamilton Jacobi equations in optimal
More informationIn all of the following equations, is the coefficient of permeability in the x direction, and is the hydraulic head.
Groundwater Seepage 1 Groundwater Seepage Simplified Steady State Fluid Flow The finite element method can be used to model both steady state and transient groundwater flow, and it has been used to incorporate
More informationAMSC 663 Project Proposal: Upgrade to the GSP Gyrokinetic Code
AMSC 663 Project Proposal: Upgrade to the GSP Gyrokinetic Code George Wilkie (gwilkie@umd.edu) Supervisor: William Dorland (bdorland@umd.edu) October 11, 2011 Abstract Simulations of turbulent plasma in
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationSparse Grids. Léopold Cambier. February 17, ICME, Stanford University
Sparse Grids & "A Dynamically Adaptive Sparse Grid Method for Quasi-Optimal Interpolation of Multidimensional Analytic Functions" from MK Stoyanov, CG Webster Léopold Cambier ICME, Stanford University
More informationUpscaling non-darcy flow
Transport in Porous Media manuscript No. (will be inserted by the editor) Upscaling non-darcy flow C.R. Garibotti M.Peszyńska Received: date / Accepted: date Abstract We consider upscaling of non-darcy
More informationGraded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons
Graded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons Anna-Margarete Sändig, Miloslav Feistauer University Stuttgart, IANS Journées
More informationANISOTROPIC EQUATIONS: UNIQUENESS AND EXISTENCE RESULTS
ANISOTROPIC EQUATIONS: UNIQUENESS AND EXISTENCE RESULTS STANISLAV ANTONTSEV, MICHEL CHIPOT Abstract. We study uniqueness of weak solutions for elliptic equations of the following type ) xi (a i (x, u)
More informationLocal Mesh Refinement with the PCD Method
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 1, pp. 125 136 (2013) http://campus.mst.edu/adsa Local Mesh Refinement with the PCD Method Ahmed Tahiri Université Med Premier
More informationMATH 124B Solution Key HW 03
6.1 LAPLACE S EQUATION MATH 124B Solution Key HW 03 6.1 LAPLACE S EQUATION 4. Solve u x x + u y y + u zz = 0 in the spherical shell 0 < a < r < b with the boundary conditions u = A on r = a and u = B on
More informationKevin James. MTHSC 206 Section 16.4 Green s Theorem
MTHSC 206 Section 16.4 Green s Theorem Theorem Let C be a positively oriented, piecewise smooth, simple closed curve in R 2. Let D be the region bounded by C. If P(x, y)( and Q(x, y) have continuous partial
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationA posteriori error estimates, stopping criteria, and adaptivity for multiphase compositional Darcy flows in porous media
A posteriori error estimates, stopping criteria, and adaptivity for multiphase compositional Darcy flows in porous media D. A. Di Pietro, M. Vohraĺık, and S. Yousef Université Montpellier 2 Marseille,
More informationREGULARITY FOR INFINITY HARMONIC FUNCTIONS IN TWO DIMENSIONS
C,α REGULARITY FOR INFINITY HARMONIC FUNCTIONS IN TWO DIMENSIONS LAWRENCE C. EVANS AND OVIDIU SAVIN Abstract. We propose a new method for showing C,α regularity for solutions of the infinity Laplacian
More informationHomogenization and error estimates of free boundary velocities in periodic media
Homogenization and error estimates of free boundary velocities in periodic media Inwon C. Kim October 7, 2011 Abstract In this note I describe a recent result ([14]-[15]) on homogenization and error estimates
More informationDynamics and Control of Rotorcraft
Dynamics and Control of Rotorcraft Helicopter Aerodynamics and Dynamics Abhishek Department of Aerospace Engineering Indian Institute of Technology, Kanpur February 3, 2018 Overview Flight Dynamics Model
More informationThe Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition
The Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition Sukjung Hwang CMAC, Yonsei University Collaboration with M. Dindos and M. Mitrea The 1st Meeting of
More informationMORTAR MULTISCALE FINITE ELEMENT METHODS FOR STOKES-DARCY FLOWS
MORTAR MULTISCALE FINITE ELEMENT METHODS FOR STOKES-DARCY FLOWS VIVETTE GIRAULT, DANAIL VASSILEV, AND IVAN YOTOV Abstract. We investigate mortar multiscale numerical methods for coupled Stokes and Darcy
More informationAPPROXIMATE IDENTITIES AND YOUNG TYPE INEQUALITIES IN VARIABLE LEBESGUE ORLICZ SPACES L p( ) (log L) q( )
Annales Academiæ Scientiarum Fennicæ Mathematica Volumen 35, 200, 405 420 APPROXIMATE IDENTITIES AND YOUNG TYPE INEQUALITIES IN VARIABLE LEBESGUE ORLICZ SPACES L p( ) (log L) q( ) Fumi-Yuki Maeda, Yoshihiro
More informationFoliations of hyperbolic space by constant mean curvature surfaces sharing ideal boundary
Foliations of hyperbolic space by constant mean curvature surfaces sharing ideal boundary David Chopp and John A. Velling December 1, 2003 Abstract Let γ be a Jordan curve in S 2, considered as the ideal
More informationFinite difference method for solving Advection-Diffusion Problem in 1D
Finite difference method for solving Advection-Diffusion Problem in 1D Author : Osei K. Tweneboah MATH 5370: Final Project Outline 1 Advection-Diffusion Problem Stationary Advection-Diffusion Problem in
More informationReservoir Flow Properties Fundamentals COPYRIGHT. Introduction
Reservoir Flow Properties Fundamentals Why This Module is Important Introduction Fundamental understanding of the flow through rocks is extremely important to understand the behavior of the reservoir Permeability
More informationA Narrow-Stencil Finite Difference Method for Hamilton-Jacobi-Bellman Equations
A Narrow-Stencil Finite Difference Method for Hamilton-Jacobi-Bellman Equations Xiaobing Feng Department of Mathematics The University of Tennessee, Knoxville, U.S.A. Linz, November 23, 2016 Collaborators
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial ifferential Equations «Viktor Grigoryan 3 Green s first identity Having studied Laplace s equation in regions with simple geometry, we now start developing some tools, which will lead
More informationMEYERS TYPE ESTIMATES FOR APPROXIMATE SOLUTIONS OF NONLINEAR ELLIPTIC EQUATIONS AND THEIR APPLICATIONS. Yalchin Efendiev.
Manuscript submitted to AIMS Journals Volume X, Number 0X, XX 200X Website: http://aimsciences.org pp. X XX MEYERS TYPE ESTIMATES FOR APPROXIMATE SOLUTIONS OF NONLINEAR ELLIPTIC EUATIONS AND THEIR APPLICATIONS
More informationUpscaling Wave Computations
Upscaling Wave Computations Xin Wang 2010 TRIP Annual Meeting 2 Outline 1 Motivation of Numerical Upscaling 2 Overview of Upscaling Methods 3 Future Plan 3 Wave Equations scalar variable density acoustic
More informationFundamental Solutions and Green s functions. Simulation Methods in Acoustics
Fundamental Solutions and Green s functions Simulation Methods in Acoustics Definitions Fundamental solution The solution F (x, x 0 ) of the linear PDE L {F (x, x 0 )} = δ(x x 0 ) x R d Is called the fundamental
More informationUniformity of the Universe
Outline Universe is homogenous and isotropic Spacetime metrics Friedmann-Walker-Robertson metric Number of numbers needed to specify a physical quantity. Energy-momentum tensor Energy-momentum tensor of
More informationMMsFEM and Streamlines
MMsFEM and Streamlines Combining a mixed multiscale FEM with streamline simulation for enhanced reservoir performance prediction on large grid models. Jørg E. Aarnes, Vegard Kippe and Knut-Andreas Lie
More information1. Introduction Boundary estimates for the second derivatives of the solution to the Dirichlet problem for the Monge-Ampere equation
POINTWISE C 2,α ESTIMATES AT THE BOUNDARY FOR THE MONGE-AMPERE EQUATION O. SAVIN Abstract. We prove a localization property of boundary sections for solutions to the Monge-Ampere equation. As a consequence
More informationHölder regularity estimation by Hart Smith and Curvelet transforms
Hölder regularity estimation by Hart Smith and Curvelet transforms Jouni Sampo Lappeenranta University Of Technology Department of Mathematics and Physics Finland 18th September 2007 This research is done
More information2 A Model, Harmonic Map, Problem
ELLIPTIC SYSTEMS JOHN E. HUTCHINSON Department of Mathematics School of Mathematical Sciences, A.N.U. 1 Introduction Elliptic equations model the behaviour of scalar quantities u, such as temperature or
More informationξ,i = x nx i x 3 + δ ni + x n x = 0. x Dξ = x i ξ,i = x nx i x i x 3 Du = λ x λ 2 xh + x λ h Dξ,
1 PDE, HW 3 solutions Problem 1. No. If a sequence of harmonic polynomials on [ 1,1] n converges uniformly to a limit f then f is harmonic. Problem 2. By definition U r U for every r >. Suppose w is a
More informationOn the Stokes semigroup in some non-helmholtz domains
On the Stokes semigroup in some non-helmholtz domains Yoshikazu Giga University of Tokyo June 2014 Joint work with Ken Abe (Nagoya U.), Katharina Schade (TU Darmstadt) and Takuya Suzuki (U. Tokyo) Contents
More informationA POSTERIORI ERROR ESTIMATES OF TRIANGULAR MIXED FINITE ELEMENT METHODS FOR QUADRATIC CONVECTION DIFFUSION OPTIMAL CONTROL PROBLEMS
A POSTERIORI ERROR ESTIMATES OF TRIANGULAR MIXED FINITE ELEMENT METHODS FOR QUADRATIC CONVECTION DIFFUSION OPTIMAL CONTROL PROBLEMS Z. LU Communicated by Gabriela Marinoschi In this paper, we discuss a
More informationSolving PDEs with freefem++
Solving PDEs with freefem++ Tutorials at Basque Center BCA Olivier Pironneau 1 with Frederic Hecht, LJLL-University of Paris VI 1 March 13, 2011 Do not forget That everything about freefem++ is at www.freefem.org
More informationMixed Finite Element Methods. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016
Mixed Finite Element Methods Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016 Linear elasticity Given the load f : Ω R n, find the displacement u : Ω R n and the
More informationEvidence that the Diffie-Hellman Problem is as Hard as Computing Discrete Logs
Evidence that the Diffie-Hellman Problem is as Hard as Computing Discrete Logs Jonah Brown-Cohen 1 Introduction The Diffie-Hellman protocol was one of the first methods discovered for two people, say Alice
More information3.4.1 Distinct Real Roots
Math 334 3.4. CONSTANT COEFFICIENT EQUATIONS 34 Assume that P(x) = p (const.) and Q(x) = q (const.), so that we have y + py + qy = 0, or L[y] = 0 (where L := d2 dx 2 + p d + q). (3.7) dx Guess a solution
More informationAnalytical Coupled Axial and Radial Productivity Model for Steady-State Flow in Horizontal Wells. Thormod Johansen*, Lesley James, Jie Cao
Int. J. Petroleum Engineering, Vol. x, No. x, 1 19 1 Analytical Coupled Axial and Radial Productivity Model for Steady-State Flow in Horizontal Wells Thormod Johansen*, Lesley James, Jie Cao Engineering
More informationDiamond and Other Carbon Materials for Solar-Thermionic Energy Conversion
Diamond and Other Carbon Materials for Solar-Thermionic Energy Conversion Timothy Fisher tsfisher@purdue.edu School of Mechanical Engineering, and Birck Nanotechnology Center Purdue University October
More informationIntegro-differential equations: Regularity theory and Pohozaev identities
Integro-differential equations: Regularity theory and Pohozaev identities Xavier Ros Oton Departament Matemàtica Aplicada I, Universitat Politècnica de Catalunya PhD Thesis Advisor: Xavier Cabré Xavier
More informationSparse Tensor Galerkin Discretizations for First Order Transport Problems
Sparse Tensor Galerkin Discretizations for First Order Transport Problems Ch. Schwab R. Hiptmair, E. Fonn, K. Grella, G. Widmer ETH Zürich, Seminar for Applied Mathematics IMA WS Novel Discretization Methods
More informationHp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme
Hp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme Aravind Balan, Michael Woopen and Georg May AICES Graduate School, RWTH Aachen University, Germany 22nd AIAA Computational
More informationFachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. A note on splitting-type variational problems with subquadratic growth Dominic reit Saarbrücken
More informationDegenerate Two-Phase Flow in Porous Media including Dynamic Effects
Degenerate Two-Phase Flow in Porous Media including Dynamic Effects Xiulei Cao Supervisor: Iuliu Sorin Pop Center for Analysis, Scientific Computing and Applications (CASA) April 9, 2014 Where innovation
More informationSolving PDEs with Multigrid Methods p.1
Solving PDEs with Multigrid Methods Scott MacLachlan maclachl@colorado.edu Department of Applied Mathematics, University of Colorado at Boulder Solving PDEs with Multigrid Methods p.1 Support and Collaboration
More informationA A + B. ra + A + 1. We now want to solve the Einstein equations in the following cases:
Lecture 29: Cosmology Cosmology Reading: Weinberg, Ch A metric tensor appropriate to infalling matter In general (see, eg, Weinberg, Ch ) we may write a spherically symmetric, time-dependent metric in
More informationAN ABSTRACT OF THE THESIS OF. Cristiano Rodrigues Garibotti for the degree of Master of Science in Mathematics
AN ABSTRACT OF THE THESIS OF Cristiano Rodrigues Garibotti for the degree of Master of Science in Mathematics presented on March 8, 007. Title: Upscaling Non-Darcy Flow Using Mixed Finite Element Method.
More informationThe Poisson boundary of certain Cartan-Hadamard manifolds of unbounded curvature
The Poisson boundary of certain Cartan-Hadamard manifolds of unbounded curvature University of Luxembourg Workshop on Boundaries Graz University of Technology June 29 July 4, 2009 Reference Marc Arnaudon,
More informationOn a parabolic-hyperbolic system for contact inhibition of cell growth
On a parabolic-hyperbolic system for contact inhibition of cell growth M. Bertsch 1, D. Hilhorst 2, H. Izuhara 3, M. Mimura 3, T. Wakasa 4 1 IAC, CNR, Rome 2 University of Paris-Sud 11 3 Meiji University
More informationA New Method for Treating Wells in Reservoir Simulation
Brigham Young University BYU ScholarsArchive All Theses and Dissertations 007-06-7 A New Method for Treating Wells in Reservoir Simulation Gregory M. Gessel Brigham Young University - Provo Follow this
More informationMotion of a Cylindrical Dielectric Boundary
Motion of a Cylindrical Dielectric Boundary Bo Li Department of Mathematics and NSF Center for Theoretical Biological Physics UC San Diego Collaborators Li-Tien Cheng, Michael White, and Shenggao Zhou
More informationMethod of Frobenius. General Considerations. L. Nielsen, Ph.D. Dierential Equations, Fall Department of Mathematics, Creighton University
Method of Frobenius General Considerations L. Nielsen, Ph.D. Department of Mathematics, Creighton University Dierential Equations, Fall 2008 Outline 1 The Dierential Equation and Assumptions 2 3 Main Theorem
More informationINTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR GRIDS. 1. Introduction
Trends in Mathematics Information Center for Mathematical Sciences Volume 9 Number 2 December 2006 Pages 0 INTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR
More informationIn the expanding Universe, a comoving volume element expands along with the cosmological flow, getting physically larger over time.
Cosmological models In the expanding Universe, a comoving volume element expands along with the cosmological flow, getting physically larger over time. The expansion is described by the scale factor R(t).
More informationTD M1 EDP 2018 no 2 Elliptic equations: regularity, maximum principle
TD M EDP 08 no Elliptic equations: regularity, maximum principle Estimates in the sup-norm I Let be an open bounded subset of R d of class C. Let A = (a ij ) be a symmetric matrix of functions of class
More informationPoint Pattern Analysis
Point Pattern Analysis Nearest Neighbor Statistics Luc Anselin http://spatial.uchicago.edu principle G function F function J function Principle Terminology events and points event: observed location of
More informationCourse Description for Real Analysis, Math 156
Course Description for Real Analysis, Math 156 In this class, we will study elliptic PDE, Fourier analysis, and dispersive PDE. Here is a quick summary of the topics will study study. They re described
More informationAbsorption cross section of RN black hole
3 Absorption cross section of RN black hole 3.1 Introduction Even though the Kerr solution is the most relevant one from an astrophysical point of view, the solution of the coupled Einstein-Maxwell equation
More informationA High Order Method for Three Phase Flow in Homogeneous Porous Media
A High Order Method for Three Phase Flow in Homogeneous Porous Media Justin Dong Advisor: B eatrice Rivi`ere February 5, 2014 Abstract The modeling of three-phase fluid flow has many important applications
More informationDetermination of thin elastic inclusions from boundary measurements.
Determination of thin elastic inclusions from boundary measurements. Elena Beretta in collaboration with E. Francini, S. Vessella, E. Kim and J. Lee September 7, 2010 E. Beretta (Università di Roma La
More informationNumerical Simulation of Flows in Highly Heterogeneous Porous Media
Numerical Simulation of Flows in Highly Heterogeneous Porous Media R. Lazarov, Y. Efendiev, J. Galvis, K. Shi, J. Willems The Second International Conference on Engineering and Computational Mathematics
More informationMULTISCALE FINITE ELEMENT METHODS FOR STOCHASTIC POROUS MEDIA FLOW EQUATIONS AND APPLICATION TO UNCERTAINTY QUANTIFICATION
MULTISCALE FINITE ELEMENT METHODS FOR STOCHASTIC POROUS MEDIA FLOW EQUATIONS AND APPLICATION TO UNCERTAINTY QUANTIFICATION P. DOSTERT, Y. EFENDIEV, AND T.Y. HOU Abstract. In this paper, we study multiscale
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
More information