Estimation of a Piecewise Constant Function Using Reparameterized Level-Set Functions
|
|
- Marshall Scott
- 6 years ago
- Views:
Transcription
1 Estimation of a iecewise Constant Function Using Reparameterized Level-Set Functions nga Berre 1,2, Martha Lien 2,1, and Trond Mannseth 2,1 1 Department of Mathematics, University of Bergen, Johs. Brunsgt. 12, N-5008 Bergen, Norway. 2 CR - Centre for ntegrated etroleum Research, University of Bergen, Realfagbygget, Allégaten 41, N-5007 Bergen, Norway. nga.berre@mi.uib.no, Martha.Lien@cipr.uib.no, and Trond.Mannseth@cipr.uib.no Abstract. n the last decade, the use of level-set functions has gained increasing popularity in solving inverse problems involving the identification of a piecewise constant function. Normally, a fine-scale representation of the level-set functions is used, yielding a high number of degrees of freedom in the estimation. n contrast, we focus on reparameterization of the level-set functions on a coarse scale. The number of coefficients in the discretized function is then reduced, providing necessary regularization for solving ill-posed problems. A coarse representation is also advantageous to reduce the computational work in solving the estimation problem. 1 Level-Set Representation of a iecewise Constant arameter Function The identification of a piecewise constant parameter function is an inverse problem arising in various applications. Eamples are image segmentation and identification of electric or fluid conductivity in reservoirs. Three features of the function can potentially be unknown: the number of regions of different constant value, the geometry of the regions, and the constant values of the parameter function. For representing piecewise constant functions in a manner suitable for identification, the level set approach [6] has become popular as it provides a fleible tool to represent region boundaries. The first approach related to inverse problems is due to Santosa [7]. Recent reviews are given by Tai and Chan [8], Burger and Osher [2], and Dorn and Lesselier [3]. The main idea is that the boundary between two regions can be represented implicitly as the zero level set of a function the level-set function. f we consider a domain Ω
2 2 nga Berre, Martha Lien, and Trond Mannseth consisting of two regions Ω 1 and Ω 2, the level-set function is defined to have the following properties: φ() > 0 for Ω 1 ; φ() < 0 for Ω 2 ; φ() = 0 for Ω 1 Ω 2. Hence, the boundary between the regions Ω 1 and Ω 2 is incorporated as the zero level-set of φ(). A parameter function p() taking different constant values c 1 in Ω 1 and c 2 in Ω 2 can now be written as p() = c 1 H(φ()) + c 2 [1 H(φ())], where H is the Heaviside function. Commonly, the level-set functions are initialized as signed distance functions. For representation of a partitioning with more than two regions, Vese and Chan [9] provide an etension of the above idea, which enables representation of up to 2 l regions with l level-set functions. 2 Coarse-Scale Level-Set Representation Usually, level-set functions are represented with one degree of freedom for each grid cell of the computational grid. We apply an approach based on a coarse-scale representation, where each level-set function is reparameterized by a few coefficients only. This approach has several advantages: the reduced number of coefficients makes sensitivity computations less demanding; the need for regularization is diminished since the possible variations in the region boundaries are more limited; and, we can achieve convergence in a low number of iterations. The latter is particularly important in solving inverse problems requiring computationally demanding forward computations and sensitivity calculations. A prime eample is the inverse problem of fluid conductivity estimation for oil reservoirs [1, 5]. A reparameterization of the level-set function is written φ() = n a i θ i (), i=1 where {θ i } denotes the set of basis functions and {a i } denotes the set of coefficients in the discretization. A coarse representation of the level-set functions enables fast identification of coarse-scale features of the parameter function with a low number of estimated coefficients. n addition, the approach provides regularization as the boundaries are confined to certain shapes based on the chosen representation. For fine-scale solutions of inverse problems, the coarse representation can enhance convergence by serving as a preconditioner for
3 arameter estimation using level-set functions 3 more fine-scale updates. Another strategy is to apply a sequential estimation, where the detail in the representation of the level-set functions is successively refined. The number of degrees of freedom in the estimation is then gradually increased, and we can achieve estimates more in correspondence with the information content of the data. When it comes to the choice of basis in the representation of the levelset functions, various alternatives are possible. However, once a set of basis functions is chosen, the reparameterized functions are restricted with respect to which shapes they can take. n the following, we discuss two different choices: a piecewise constant representation and a continuous representation. A piecewise constant representation of the level-set functions [1, 5] is well suited in combination with adaptive multiscale estimation [4], which provides a rough identification of the number and geometry of the regions of different constant value for the parameter function. n Berre et al. [1], reparameterization by characteristic basis functions in combination with a narrow-band approach is proposed. The representation of the region boundaries is gradually refined during estimation. n Lien et al. [5], a sequential approach is developed, where the number of regions of constant parameter value is sought found as part of the estimation. Through successive estimations with increasing resolution in the representation of p() identification of rather general parameter functions can be achieved. The choice of characteristic basis functions yields piecewise constant levelset functions. This leads to updates of the boundaries between the constant states of p(), where the jumps are determined by the spatial support of the basis functions for φ(), regardless of the resolution of the computational grid for the forward problem. The two plots to the left in Figure 1 illustrate this situation in 1-D with n = 2, supp θ 1 = [0, 1/2), and supp θ 2 = [1/2, 1]. A continuous representation of the level-set functions enables more gradual updates of the region boundaries, depending on the resolution of the computational grid for the forward problem. Hence, the number of possible updates of the region boundaries is greatly enhanced. The two plots to the right in Figure 1 illustrate this situation in 1-D with n = 2 and θ 1 and θ 2 being the standard linear basis on [0, 1]. As a simple eample of a continuous reparameterization of φ() with a low number of coefficients in 2-D, we consider a bilinear basis. The computational domain is partitioned into a coarse quadrilateral grid of rectangular elements. A bilinear reparameterization of a level-set function on a reference element D r = [0, 1] [0, 1] is given by φ( 1, 2 ) = a 1 ( 1 1)( 2 1) a 2 ( 1 1) 2 a 3 1 ( 2 1) + a By increasing the resolution of the quadrilateral grid, the boundaries can have more comple shapes. An illustration of the potential of applying a coarse-scale reparameterization of the level-set functions is given below; see Figures 2 and 3. We consider the inverse problem of fluid conductivity estimation for an oil reservoir based
4 4 nga Berre, Martha Lien, and Trond Mannseth φ φ a 1 a 1 a 2 a 2 p c 1 c 2 p c 1 c 2 Fig. 1. Coarse representation of φ with two coefficients, a 1 and a 2, by characteristic (top left) and linear (top right) basis functions. Optimizing with respect to a 1 cause φ to change sign in some area. The resulting updates of the parameter function are shown in the lower figures. The initial states are drawn by solid lines and the final states by dashed lines. on a level-set representation with bilinear basis functions. The forward problem describes horizontal, two-phase, immiscible, and incompressible fluid flow in porous media. The available data for the inversion consist of time series of pressure data d logged in the injection wells. There are nine injection () and four production () wells; hence, the data are very sparsely distributed. The data are obtained from a forward simulation with the reference fluid conductivity, where normally distributed errors are added to the calculated pressures. To illustrate the regularizing effect of the coarse level-set representation, we minimize a weighted least squares objective function with no additional regularizing terms: J(p) = [m(p) d] T C 1 [m(p) d]. Here m(p) denote the calculated pressures given a parameter function p(), and C denotes the (diagonal) covariance matri of the measurement errors in the data. The reference fluid conductivity describes a channel with areas of low conductivity on both sides where the coarse-scale features are contaminated with fine-scale conductivity variation; see Figure 2 (left). The initial guess for p() was a constant value for the whole reservoir. The estimation was started with a coarse level-set representation on a grid of only one element before we continued with a quadrilateral grid of four elements. With our coarse level-set approach, we are able to recover the main trends in the conductivity distribution. The data are not reconciled by the coarse estimate though the objective function is greatly reduced. Hence, in this case, the coarse-scale estimate can serve as a good starting point for more fine-scale estimations.
5 arameter estimation using level-set functions Fig. 2. Reference fluid conductivity and final estimation result Fig. 3. Level-set function at convergence together with the surface z = 0. References 1.. Berre, M. Lien, and T. Mannseth. A level-set corrector to an adaptive multiscale permeability prediction. Accepted for publication in Comput. Geosci. 2. M. Burger and S. Osher. A survey on level set methods for inverse problems and optimal design. Euro. Jnl of Applied Mathematics, 16(2):263 1, April O. Dorn and D. Lesselier. Level set methods for inverse scattering. nverse roblems, 22:67 131, A.-A. Grimstad, T. Mannseth, G. Nævdal, and H. Urkedal. Adaptive multiscale permeability estimation. Comput. Geosci., 7(1):1 25, M. Lien,. Berre, and T. Mannseth. Combined adaptive multiscale and level-set parameter estimation. Multiscale Model. Simul., 4(4): , S. Osher and J. A. Sethian. Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations. J. Comput. hys., 79(1):12 49, F. Santosa. A level-set approach for inverse problems involving obstacles. ESAM Contrôle Optim. Calc. Var., 1:17 33, 1995/ X.-C. Tai and Tony F. Chan. A survey on multiple level set methods with applications for identifying piecewise constant functions. nt. J. Numer. Anal. Model., 1(1):25 47, L. A. Vese and T. F. Chan. A multiphase level set framework for image segmentation via the Mumford and Shah model. nt. J. of Comput. Vision, 02.
Journal of Physics: Conference Series. Recent citations. To cite this article: Inga Berre et al 2008 J. Phys.: Conf. Ser.
Journal of Physics: Conference Series A multi-level strategy for regularized parameter identification using nonlinear reparameterization with sample application for reservoir characterization Recent citations
More informationReservoir description using a binary level set model
Comput Visual Sci DOI.7/s79-8-2- REGULAR ARTICLE Reservoir description using a binary level set model Lars Kristian Nielsen Hongwei Li Xue-Cheng Tai Sigurd Ivar Aanonsen Magne Espedal Received: 7 September
More informationControl Volume Finite Difference On Adaptive Meshes
Control Volume Finite Difference On Adaptive Meshes Sanjay Kumar Khattri, Gunnar E. Fladmark, Helge K. Dahle Department of Mathematics, University Bergen, Norway. sanjay@mi.uib.no Summary. In this work
More informationMultiple Motion and Occlusion Segmentation with a Multiphase Level Set Method
Multiple Motion and Occlusion Segmentation with a Multiphase Level Set Method Yonggang Shi, Janusz Konrad, W. Clem Karl Department of Electrical and Computer Engineering Boston University, Boston, MA 02215
More informationOutline. Level Set Methods. For Inverse Obstacle Problems 4. Introduction. Introduction. Martin Burger
For Inverse Obstacle Problems Martin Burger Outline Introduction Optimal Geometries Inverse Obstacle Problems & Shape Optimization Sensitivity Analysis based on Gradient Flows Numerical Methods University
More informationAutomated Segmentation Using a Fast Implementation of the Chan-Vese Models
Automated Segmentation Using a Fast Implementation of the Chan-Vese Models Huan Xu, and Xiao-Feng Wang,,3 Intelligent Computation Lab, Hefei Institute of Intelligent Machines, Chinese Academy of Science,
More informationMethod of Background Subtraction for Medical Image Segmentation
Method of Background Subtraction for Medical Image Segmentation Seongjai Kim Department of Mathematics and Statistics, Mississippi State University Mississippi State, MS 39762, USA and Hyeona Lim Department
More informationMathematical Methods in Fluid Dynamics and Simulation of Giant Oil and Gas Reservoirs. 3-5 September 2012 Swissotel The Bosphorus, Istanbul, Turkey
Mathematical Methods in Fluid Dynamics and Simulation of Giant Oil and Gas Reservoirs 3-5 September 2012 Swissotel The Bosphorus, Istanbul, Turkey Fast and robust solvers for pressure systems on the GPU
More informationParametric Level Set Methods for Inverse Problems
arxiv:1007.2467v2 [math.ap] 6 Oct 2011 Parametric Level Set Methods for Inverse Problems Alireza Aghasi, Misha Kilmer, Eric L. Miller October 7, 2011 Abstract In this paper, a parametric level set method
More informationCollege of Engineering, Trivandrum.
Analysis of CT Liver Images Using Level Sets with Bayesian Analysis-A Hybrid Approach Sajith A.G 1, Dr. Hariharan.S 2 1 Research Scholar, 2 Professor, Department of Electrical&Electronics Engineering College
More informationAn Active Contour Model without Edges
An Active Contour Model without Edges Tony Chan and Luminita Vese Department of Mathematics, University of California, Los Angeles, 520 Portola Plaza, Los Angeles, CA 90095-1555 chan,lvese@math.ucla.edu
More informationBACK AND FORTH ERROR COMPENSATION AND CORRECTION METHODS FOR REMOVING ERRORS INDUCED BY UNEVEN GRADIENTS OF THE LEVEL SET FUNCTION
BACK AND FORTH ERROR COMPENSATION AND CORRECTION METHODS FOR REMOVING ERRORS INDUCED BY UNEVEN GRADIENTS OF THE LEVEL SET FUNCTION TODD F. DUPONT AND YINGJIE LIU Abstract. We propose a method that significantly
More informationUnstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications
Unstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications Per-Olof Persson (persson@mit.edu) Department of Mathematics Massachusetts Institute of Technology http://www.mit.edu/
More informationKey words. Level set, energy minimizing, partial differential equations, segmentation.
A VARIANT OF THE LEVEL SET METHOD AND APPLICATIONS TO IMAGE SEGMENTATION JOHAN LIE, MARIUS LYSAKER, AND XUE-CHENG TAI Abstract. In this paper we propose a variant of the level set formulation for identifying
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 informationComputational Fluid Dynamics - Incompressible Flows
Computational Fluid Dynamics - Incompressible Flows March 25, 2008 Incompressible Flows Basis Functions Discrete Equations CFD - Incompressible Flows CFD is a Huge field Numerical Techniques for solving
More informationSegmentation Using Active Contour Model and Level Set Method Applied to Medical Images
Segmentation Using Active Contour Model and Level Set Method Applied to Medical Images Dr. K.Bikshalu R.Srikanth Assistant Professor, Dept. of ECE, KUCE&T, KU, Warangal, Telangana, India kalagaddaashu@gmail.com
More informationALIREZA AGHASI, MISHA KILMER, ERIC L. MILLER
PARAMETRIC LEVEL SET METHODS FOR INVERSE PROBLEMS ALIREZA AGHASI, MISHA KILMER, ERIC L. MILLER AAGHAS1@ECE.TUFTS.EDU, MISHA.KILMER@TUFTS.EDU, ELMILLER@ECE.TUTS.EDU Abstract. In this paper, a parametric
More informationThe Level Set Method. Lecture Notes, MIT J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations
The Level Set Method Lecture Notes, MIT 16.920J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations Per-Olof Persson persson@mit.edu March 7, 2005 1 Evolving Curves and Surfaces Evolving
More informationDr. Ulas Bagci
Lecture 9: Deformable Models and Segmentation CAP-Computer Vision Lecture 9-Deformable Models and Segmentation Dr. Ulas Bagci bagci@ucf.edu Lecture 9: Deformable Models and Segmentation Motivation A limitation
More informationLecture 2.2 Cubic Splines
Lecture. Cubic Splines Cubic Spline The equation for a single parametric cubic spline segment is given by 4 i t Bit t t t i (..) where t and t are the parameter values at the beginning and end of the segment.
More informationElectromagnetic migration of marine CSEM data in areas with rough bathymetry Michael S. Zhdanov and Martin Čuma*, University of Utah
Electromagnetic migration of marine CSEM data in areas with rough bathymetry Michael S. Zhdanov and Martin Čuma*, University of Utah Summary In this paper we present a new approach to the interpretation
More informationModule 1: Introduction to Finite Element Analysis. Lecture 4: Steps in Finite Element Analysis
25 Module 1: Introduction to Finite Element Analysis Lecture 4: Steps in Finite Element Analysis 1.4.1 Loading Conditions There are multiple loading conditions which may be applied to a system. The load
More informationComparison of different solvers for two-dimensional steady heat conduction equation ME 412 Project 2
Comparison of different solvers for two-dimensional steady heat conduction equation ME 412 Project 2 Jingwei Zhu March 19, 2014 Instructor: Surya Pratap Vanka 1 Project Description The purpose of this
More informationMedical Image Segmentation using Level Sets
Medical Image Segmentation using Level Sets Technical Report #CS-8-1 Tenn Francis Chen Abstract Segmentation is a vital aspect of medical imaging. It aids in the visualization of medical data and diagnostics
More informationMultigrid Pattern. I. Problem. II. Driving Forces. III. Solution
Multigrid Pattern I. Problem Problem domain is decomposed into a set of geometric grids, where each element participates in a local computation followed by data exchanges with adjacent neighbors. The grids
More informationPiecewise Constant Level Set Methods for Multiphase Motion
Piecewise Constant Level Set Methods for Multiphase Motion Hongwei Li Center for Integrated Petroleum Research, University of Bergen CIPR room 4103, Allgaten 41, N-5007 Bergen, Norway Email: hongwei.li@cipr.uib.no
More informationA Toolbox of Level Set Methods
A Toolbox of Level Set Methods Ian Mitchell Department of Computer Science University of British Columbia http://www.cs.ubc.ca/~mitchell mitchell@cs.ubc.ca research supported by the Natural Science and
More informationLevel-set and ALE Based Topology Optimization Using Nonlinear Programming
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Level-set and ALE Based Topology Optimization Using Nonlinear Programming Shintaro Yamasaki
More informationThe Immersed Interface Method
The Immersed Interface Method Numerical Solutions of PDEs Involving Interfaces and Irregular Domains Zhiiin Li Kazufumi Ito North Carolina State University Raleigh, North Carolina Society for Industrial
More informationStudies of the Continuous and Discrete Adjoint Approaches to Viscous Automatic Aerodynamic Shape Optimization
Studies of the Continuous and Discrete Adjoint Approaches to Viscous Automatic Aerodynamic Shape Optimization Siva Nadarajah Antony Jameson Stanford University 15th AIAA Computational Fluid Dynamics Conference
More informationHistory Matching, Forecasting and Updating
History Matching, Forecasting and Updating Dr. Helmy Sayyouh Petroleum Engineering Cairo University 12/26/2017 1 History Matching The most practical method for testing a reservoir model s validity and
More informationApplication of Finite Volume Method for Structural Analysis
Application of Finite Volume Method for Structural Analysis Saeed-Reza Sabbagh-Yazdi and Milad Bayatlou Associate Professor, Civil Engineering Department of KNToosi University of Technology, PostGraduate
More informationTHE level set method proposed by Osher and Sethian [1] is
VERSION NOVEMBER 30, 2004 1 A Binary Level Set Model and some Applications to Mumford-Shah Image Segmentation Johan Lie, Marius Lysaker and Xue-Cheng Tai Abstract In this work we propose a PDE based level
More informationLevel set methods for watershed image segmentation
Level set methods for watershed image segmentation Erlend Hodneland 1, Xue-Cheng Tai 2, Joachim Weickert 3, Nickolay V. Bukoreshtliev 1, Arvid Lundervold 1, and Hans-Hermann Gerdes 1 1 Department of Biomedicine,
More informationThe Finite Element Method
The Finite Element Method A Practical Course G. R. Liu and S. S. Quek Chapter 1: Computational modeling An overview 1 CONTENTS INTRODUCTION PHYSICAL PROBLEMS IN ENGINEERING COMPUTATIONAL MODELLING USING
More informationLevel set methods Formulation of Interface Propagation Boundary Value PDE Initial Value PDE Motion in an externally generated velocity field
Level Set Methods Overview Level set methods Formulation of Interface Propagation Boundary Value PDE Initial Value PDE Motion in an externally generated velocity field Convection Upwind ddifferencingi
More informationApplied Reservoir Simulation. Course. Workshop Problems
Applied Reservoir Simulation Course Workshop Problems Reservoir Simulation Application Training Course and (Eclipse) Workshop GeoQuest Training and Development, And NExT Denver and Houston March 06 Applied
More informationVariogram Inversion and Uncertainty Using Dynamic Data. Simultaneouos Inversion with Variogram Updating
Variogram Inversion and Uncertainty Using Dynamic Data Z. A. Reza (zreza@ualberta.ca) and C. V. Deutsch (cdeutsch@civil.ualberta.ca) Department of Civil & Environmental Engineering, University of Alberta
More informationContents. I The Basic Framework for Stationary Problems 1
page v Preface xiii I The Basic Framework for Stationary Problems 1 1 Some model PDEs 3 1.1 Laplace s equation; elliptic BVPs... 3 1.1.1 Physical experiments modeled by Laplace s equation... 5 1.2 Other
More informationIsophote-Based Interpolation
Brigham Young University BYU ScholarsArchive All Faculty Publications 1998-10-01 Isophote-Based Interpolation Bryan S. Morse morse@byu.edu Duane Schwartzwald Follow this and additional works at: http://scholarsarchive.byu.edu/facpub
More informationNSCT BASED LOCAL ENHANCEMENT FOR ACTIVE CONTOUR BASED IMAGE SEGMENTATION APPLICATION
DOI: 10.1917/ijivp.010.0004 NSCT BASED LOCAL ENHANCEMENT FOR ACTIVE CONTOUR BASED IMAGE SEGMENTATION APPLICATION Hiren Mewada 1 and Suprava Patnaik Department of Electronics Engineering, Sardar Vallbhbhai
More informationCS205b/CME306. Lecture 9
CS205b/CME306 Lecture 9 1 Convection Supplementary Reading: Osher and Fedkiw, Sections 3.3 and 3.5; Leveque, Sections 6.7, 8.3, 10.2, 10.4. For a reference on Newton polynomial interpolation via divided
More informationHeterogeneity preserving upscaling for heat transport in fractured geothermal reservoirs
Comput Geosci (2018) 22:451 467 https://doi.org/10.1007/s10596-017-9704-6 ORIGINAL PAPER Heterogeneity preserving upscaling for heat transport in fractured geothermal reservoirs Anna Nissen 1,2 Eirik Keilegavlen
More informationIsophote-Based Interpolation
Isophote-Based Interpolation Bryan S. Morse and Duane Schwartzwald Department of Computer Science, Brigham Young University 3361 TMCB, Provo, UT 84602 {morse,duane}@cs.byu.edu Abstract Standard methods
More informationLevel Set Methods and Fast Marching Methods
Level Set Methods and Fast Marching Methods I.Lyulina Scientific Computing Group May, 2002 Overview Existing Techniques for Tracking Interfaces Basic Ideas of Level Set Method and Fast Marching Method
More informationA Multiscale Two-Point Flux-Approximation Method
A Multiscale Two-Point Flux-Approximation Method Olav Møyner Knut Andreas Lie October 25, 2013 A large number of multiscale finite-volume methods have been developed over the past decade to compute conservative
More informationFinite element solution of multi-scale transport problems using the least squares based bubble function enrichment
Finite element solution of multi-scale transport problems using the least squares based bubble function enrichment A. Yazdani a, V. Nassehi b1 a Cranfield University, School of Applied Sciences, Cranfield,
More informationDriven Cavity Example
BMAppendixI.qxd 11/14/12 6:55 PM Page I-1 I CFD Driven Cavity Example I.1 Problem One of the classic benchmarks in CFD is the driven cavity problem. Consider steady, incompressible, viscous flow in a square
More informationAMS527: Numerical Analysis II
AMS527: Numerical Analysis II A Brief Overview of Finite Element Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 1 / 25 Overview Basic concepts Mathematical
More informationMULTIPHASE LEVEL SET EVOLUTION WITH APPLICATIONS TO AUTOMATIC GENERATIONAL TRACKING OF CELL DIVISION OF ESCHERICHIA COLI. A Thesis.
MULTIPHASE LEVEL SET EVOLUTION WITH APPLICATIONS TO AUTOMATIC GENERATIONAL TRACKING OF CELL DIVISION OF ESCHERICHIA COLI A Thesis Presented to the Faculty of San Diego State University In Partial Fulfillment
More informationA Parallel, Multiscale Approach to Reservoir Modeling. Omer Inanc Tureyen and Jef Caers Department of Petroleum Engineering Stanford University
A Parallel, Multiscale Approach to Reservoir Modeling Omer Inanc Tureyen and Jef Caers Department of Petroleum Engineering Stanford University 1 Abstract With the advance of CPU power, numerical reservoir
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More information7 Fractions. Number Sense and Numeration Measurement Geometry and Spatial Sense Patterning and Algebra Data Management and Probability
7 Fractions GRADE 7 FRACTIONS continue to develop proficiency by using fractions in mental strategies and in selecting and justifying use; develop proficiency in adding and subtracting simple fractions;
More informationA SURVEY ON MULTIPLE LEVEL SET METHODS WITH APPLICATIONS FOR IDENTIFYING PIECEWISE CONSTANT FUNCTIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume, Number, Pages 5 47 c 4 Institute for Scientific Computing and Information A SURVEY ON MULTIPLE LEVEL SET METHODS WITH APPLICATIONS FOR IDENTIFYING
More informationJune 5, Institute of Structural Analysis Graz University of Technology Lessingstr. 25/II, 8010 Graz, Austria
Higher-order meshing of implicit geometries part I: Integration and interpolation in cut elements arxiv:706.00578v [cs.na] 2 Jun 207 T.P. Fries, S. Omerović, D. Schöllhammer, J. Steidl June 5, 207 Institute
More informationWAVE PATTERNS, WAVE INDUCED FORCES AND MOMENTS FOR A GRAVITY BASED STRUCTURE PREDICTED USING CFD
Proceedings of the ASME 2011 30th International Conference on Ocean, Offshore and Arctic Engineering OMAE2011 June 19-24, 2011, Rotterdam, The Netherlands OMAE2011-49593 WAVE PATTERNS, WAVE INDUCED FORCES
More informationIntroduction to Computational Fluid Dynamics Mech 122 D. Fabris, K. Lynch, D. Rich
Introduction to Computational Fluid Dynamics Mech 122 D. Fabris, K. Lynch, D. Rich 1 Computational Fluid dynamics Computational fluid dynamics (CFD) is the analysis of systems involving fluid flow, heat
More informationComputational Study of Laminar Flowfield around a Square Cylinder using Ansys Fluent
MEGR 7090-003, Computational Fluid Dynamics :1 7 Spring 2015 Computational Study of Laminar Flowfield around a Square Cylinder using Ansys Fluent Rahul R Upadhyay Master of Science, Dept of Mechanical
More informationCHAPTER 1. Introduction
ME 475: Computer-Aided Design of Structures 1-1 CHAPTER 1 Introduction 1.1 Analysis versus Design 1.2 Basic Steps in Analysis 1.3 What is the Finite Element Method? 1.4 Geometrical Representation, Discretization
More informationPhD Student. Associate Professor, Co-Director, Center for Computational Earth and Environmental Science. Abdulrahman Manea.
Abdulrahman Manea PhD Student Hamdi Tchelepi Associate Professor, Co-Director, Center for Computational Earth and Environmental Science Energy Resources Engineering Department School of Earth Sciences
More informationSnakes, Active Contours, and Segmentation Introduction and Classical Active Contours Active Contours Without Edges
Level Sets & Snakes Snakes, Active Contours, and Segmentation Introduction and Classical Active Contours Active Contours Without Edges Scale Space and PDE methods in image analysis and processing - Arjan
More informationDigital Image Processing Fundamentals
Ioannis Pitas Digital Image Processing Fundamentals Chapter 7 Shape Description Answers to the Chapter Questions Thessaloniki 1998 Chapter 7: Shape description 7.1 Introduction 1. Why is invariance to
More informationPressure Correction Scheme for Incompressible Fluid Flow
AALTO UNIVERSITY School of Chemical Technology CHEM-E7160 Fluid Flow in Process Units Pressure Correction Scheme for Incompressible Fluid Flow Ong Chin Kai 620503 Lee De Ming Benedict 620448 Page 1 Abstract
More informationGeodesic Active Contours with Combined Shape and Appearance Priors
Geodesic Active Contours with Combined Shape and Appearance Priors Rami Ben-Ari 1 and Dror Aiger 1,2 1 Orbotech LTD, Yavneh, Israel 2 Ben Gurion University, Be er Sheva, Israel {rami-ba,dror-ai}@orbotech.com
More informationSELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND
Student Submission for the 5 th OpenFOAM User Conference 2017, Wiesbaden - Germany: SELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND TESSA UROIĆ Faculty of Mechanical Engineering and Naval Architecture, Ivana
More informationCompressible Flow in a Nozzle
SPC 407 Supersonic & Hypersonic Fluid Dynamics Ansys Fluent Tutorial 1 Compressible Flow in a Nozzle Ahmed M Nagib Elmekawy, PhD, P.E. Problem Specification Consider air flowing at high-speed through a
More informationA Fast Local Level Set Method for Inverse Gravimetry
Commun. Comput. Phys. doi: 28/cicp.7.22a Vol., No. 4, pp. 44-7 October 2 A Fast Local Level Set Method for Inverse Gravimetry Victor Isakov, Shingyu Leung 2 and Jianliang Qian 3, Department of Mathematics
More informationFast computation of arrival times in heterogeneous media
Fast computation of arrival times in heterogeneous media Inga Berre (inga.berre@mi.uib.no) Dept. of Mathematics, University of Bergen, Johs. Brunsgt. 2, N-5008 Bergen, Norway. Kenneth Hvistendal Karlsen
More informationCHAPTER 9 INPAINTING USING SPARSE REPRESENTATION AND INVERSE DCT
CHAPTER 9 INPAINTING USING SPARSE REPRESENTATION AND INVERSE DCT 9.1 Introduction In the previous chapters the inpainting was considered as an iterative algorithm. PDE based method uses iterations to converge
More informationCIPC Louis Mattar. Fekete Associates Inc. Analytical Solutions in Well Testing
CIPC 2003 Louis Mattar Fekete Associates Inc Analytical Solutions in Well Testing Well Test Equation 2 P 2 P 1 P + = x 2 y 2 α t Solutions Analytical Semi-Analytical Numerical - Finite Difference Numerical
More informationPARAMETRIC SHAPE AND TOPOLOGY OPTIMIZATION WITH RADIAL BASIS FUNCTIONS
PARAMETRIC SHAPE AND TOPOLOGY OPTIMIZATION WITH RADIAL BASIS FUNCTIONS Michael Yu Wang 1 and Shengyin Wang 1 Department of Automation and Computer-Aided Engineering The Chinese University of Hong Kong
More informationReporting Mesh Statistics
Chapter 15. Reporting Mesh Statistics The quality of a mesh is determined more effectively by looking at various statistics, such as maximum skewness, rather than just performing a visual inspection. Unlike
More informationThis document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore.
This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore. Title A binary level set model and some applications to Mumford-Shah image segmentation Author(s) Lie, Johan;
More informationIntroduction to C omputational F luid Dynamics. D. Murrin
Introduction to C omputational F luid Dynamics D. Murrin Computational fluid dynamics (CFD) is the science of predicting fluid flow, heat transfer, mass transfer, chemical reactions, and related phenomena
More informationA Framework for the Upscaling of the Electrical Conductivity in the Quasi-static Maxwell s Equations
A Framework for the Upscaling of the Electrical Conductivity in the Quasi-static Maxwell s Equations arxiv:1610.02948v1 [math.na] 7 Oct 2016 L. A. Caudillo-Mata*, E. Haber, L. J. Heagy, and C. Schwarzbach
More informationDomain Decomposition and hp-adaptive Finite Elements
Domain Decomposition and hp-adaptive Finite Elements Randolph E. Bank 1 and Hieu Nguyen 1 Department of Mathematics, University of California, San Diego, La Jolla, CA 9093-011, USA, rbank@ucsd.edu. Department
More informationToward reservoir simulation on geological grid models
1 Toward reservoir simulation on geological grid models JØRG E. AARNES and KNUT ANDREAS LIE SINTEF ICT, Dept. of Applied Mathematics, P.O. Box 124 Blindern, NO-0314 Oslo, Norway Abstract We present a reservoir
More informationMultiphase Dynamic Labeling for Variational Recognition-Driven Image Segmentation
Multiphase Dynamic Labeling for Variational Recognition-Driven Image Segmentation Daniel Cremers 1, Nir Sochen 2, and Christoph Schnörr 3 1 Department of Computer Science University of California, Los
More informationCFD Best Practice Guidelines: A process to understand CFD results and establish Simulation versus Reality
CFD Best Practice Guidelines: A process to understand CFD results and establish Simulation versus Reality Judd Kaiser ANSYS Inc. judd.kaiser@ansys.com 2005 ANSYS, Inc. 1 ANSYS, Inc. Proprietary Overview
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY. Analyzing wind flow around the square plate using ADINA Project. Ankur Bajoria
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Analyzing wind flow around the square plate using ADINA 2.094 - Project Ankur Bajoria May 1, 2008 Acknowledgement I would like to thank ADINA R & D, Inc for the full
More informationPTE 519 Lecture Note Finite Difference Approximation (Model)
PTE 519 Lecture Note 3 3.0 Finite Difference Approximation (Model) In this section of the lecture material, the focus is to define the terminology and to summarize the basic facts. The basic idea of any
More informationINTEGRATION OF LOCAL-GLOBAL UPSCALING AND GRID ADAPTIVITY FOR SIMULATION OF SUBSURFACE FLOW IN HETEROGENEOUS FORMATIONS
INTEGRATION OF LOCAL-GLOBAL UPSCALING AND GRID ADAPTIVITY FOR SIMULATION OF SUBSURFACE FLOW IN HETEROGENEOUS FORMATIONS M. GERRITSEN AND J. V. LAMBERS Abstract. We propose a methodology, called multi-level
More informationTIPS4Math Grades 4 to 6 Overview Grade 4 Grade 5 Grade 6 Collect, Organize, and Display Primary Data (4+ days)
Collect, Organize, and Display Primary Data (4+ days) Collect, Organize, Display and Interpret Categorical Data (5+ days) 4m88 Collect data by conducting a survey or an experiment to do with the 4m89 Collect
More information3D Inversion of Time-Domain Electromagnetic Data for Ground Water Aquifers
3D Inversion of Time-Domain Electromagnetic Data for Ground Water Aquifers Elliot M. Holtham 1, Mike McMillan 1 and Eldad Haber 2 (1) Computational Geosciences Inc. (2) University of British Columbia Summary
More informationcuibm A GPU Accelerated Immersed Boundary Method
cuibm A GPU Accelerated Immersed Boundary Method S. K. Layton, A. Krishnan and L. A. Barba Corresponding author: labarba@bu.edu Department of Mechanical Engineering, Boston University, Boston, MA, 225,
More informationA DISCONTINUOUS GALERKIN METHOD FOR TRANSPORT IN FRACTURED MEDIA USING UNSTRUCTURED TRIANGULAR GRIDS
A DISCONTINUOUS GALERKIN METHOD FOR TRANSPORT IN FRACTURED MEDIA USING UNSTRUCTURED TRIANGULAR GRIDS B. EIKEMO, K.-A. LIE, G. T. EIGESTAD, AND H. K. DAHLE Abstract. The possibility to couple discrete (fractures,
More informationENERGY-224 Reservoir Simulation Project Report. Ala Alzayer
ENERGY-224 Reservoir Simulation Project Report Ala Alzayer Autumn Quarter December 3, 2014 Contents 1 Objective 2 2 Governing Equations 2 3 Methodolgy 3 3.1 BlockMesh.........................................
More informationSIZE PRESERVING MESH GENERATION IN ADAPTIVITY PROCESSES
Congreso de Métodos Numéricos en Ingeniería 25-28 junio 2013, Bilbao, España c SEMNI, 2013 SIZE PRESERVING MESH GENERATION IN ADAPTIVITY PROCESSES Eloi Ruiz-Gironés 1, Xevi Roca 2 and Josep Sarrate 1 1:
More informationImplicit Active Model using Radial Basis Function Interpolated Level Sets
Implicit Active Model using Radial Basis Function Interpolated Level Sets Xianghua Xie and Majid Mirmehdi Department of Computer Science University of Bristol, Bristol BS8 1UB, England. {xie,majid}@cs.bris.ac.uk
More informationHoughton Mifflin MATHEMATICS Level 5 correlated to NCTM Standard
s 2000 Number and Operations Standard Understand numbers, ways of representing numbers, relationships among numbers, and number systems understand the place-value structure of the TE: 4 5, 8 11, 14 17,
More informationPhase field, level set method and PCLSM p.1/133
Phase field, level set method and PCLSM Xue-Cheng Tai Xue-Cheng.Tai@mi.uib.no, http://www.mi.uib.no/%7etai. Department of Mathematics, University of Bergen, Johannes Brunsgt. 12, 5007 Bergen, Norway, Phase
More informationOPTIMIZATION FOR AUTOMATIC HISTORY MATCHING
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 2, Supp, Pages 131 137 c 2005 Institute for Scientific Computing and Information OPTIMIZATION FOR AUTOMATIC HISTORY MATCHING Abstract. SHUGUANG
More informationOn-line, Incremental Learning of a Robust Active Shape Model
On-line, Incremental Learning of a Robust Active Shape Model Michael Fussenegger 1, Peter M. Roth 2, Horst Bischof 2, Axel Pinz 1 1 Institute of Electrical Measurement and Measurement Signal Processing
More informationCerebral Artery Segmentation with Level Set Methods
H. Ho, P. Bier, G. Sands, P. Hunter, Cerebral Artery Segmentation with Level Set Methods, Proceedings of Image and Vision Computing New Zealand 2007, pp. 300 304, Hamilton, New Zealand, December 2007.
More informationThe effect of irregular interfaces on the BDDC method for the Navier-Stokes equations
153 The effect of irregular interfaces on the BDDC method for the Navier-Stokes equations Martin Hanek 1, Jakub Šístek 2,3 and Pavel Burda 1 1 Introduction The Balancing Domain Decomposition based on Constraints
More informationA COUPLED FINITE VOLUME SOLVER FOR THE SOLUTION OF LAMINAR TURBULENT INCOMPRESSIBLE AND COMPRESSIBLE FLOWS
A COUPLED FINITE VOLUME SOLVER FOR THE SOLUTION OF LAMINAR TURBULENT INCOMPRESSIBLE AND COMPRESSIBLE FLOWS L. Mangani Maschinentechnik CC Fluidmechanik und Hydromaschinen Hochschule Luzern Technik& Architektur
More informationCommunications in Applied Mathematics and Computational Science
Communications in Applied Mathematics and Computational Science Volume 1 No. 1 2006 A COMPARISON OF THE EXTENDED FINITE ELEMENT METHOD WITH THE IMMERSED INTERFACE METHOD FOR ELLIPTIC EQUATIONS WITH DISCONTINUOUS
More informationNumerical schemes for Hamilton-Jacobi equations, control problems and games
Numerical schemes for Hamilton-Jacobi equations, control problems and games M. Falcone H. Zidani SADCO Spring School Applied and Numerical Optimal Control April 23-27, 2012, Paris Lecture 2/3 M. Falcone
More informationDesign and Analysis of an Euler Transformation Algorithm Applied to Full-Polarimetric ISAR Imagery
Design and Analysis of an Euler Transformation Algorithm Applied to Full-Polarimetric ISAR Imagery Christopher S. Baird Advisor: Robert Giles Submillimeter-Wave Technology Laboratory (STL) Presented in
More information