Fully implicit higher-order schemes applied to polymer flooding. K.-A. Lie(SINTEF), T.S. Mykkeltvedt(IRIS), X. Raynaud (SINTEF)
|
|
- Marshall Weaver
- 5 years ago
- Views:
Transcription
1 Fully implicit higher-order schemes applied to polymer flooding K.-A. Lie(SINTEF), T.S. Mykkeltvedt(IRIS), X. Raynaud (SINTEF)
2 Motivation Challenging and complex to simulate water-based EOR methods, unresolved simulation can lead to misleading predictions Polymer flooding enhances the water s ability to push the oil through the reservoir, and can reduce channeling through high flow zones Crucial to capture polymer fronts sharply to resolve the nonlinear displacement mechanisms correctly Polymer fronts will in the worst case be linear waves and generally have significantly less self-sharpening effects than water fronts challenge when using standard low-order methods
3 General flow equations for two-phase flow Gathering the equations, we have (φρ α S α ) t + (ρ α u α ) = qα, α = {w, n} u α = Kk rα µ α ( pα ρ α g z ) p c = p n p w, S w + S n = 1
4 Model for polymer flooding Simple model: introduce extra immiscible component and mixture law Polymer transported in water: u w = k rw(s) µ w,eff (c) R k (c) K( p w ρ ) wg z u p = krw(s) µ p,eff (c)r k (c) K( p w ρ ) wg z Conservation of polymer component: [ t φ(1 Sipv)cb ws + ρ r C a (1 φ) ] + (cb ) w u p = qp
5 Model for polymer flooding Simple model: introduce extra immiscible component and mixture law Polymer transported in water: u w = viscosity enhancement k rw(s) µ w,eff (c) R k (c) K( p w ρ ) wg z Todd Longstaff mixing: 1 µ = 1 c/cm w,eff µm(c) ω µ 1 ω w + c/cm µm(c) ω µ 1 ω p u p = krw(s) µ p,eff (c)r k (c) K( p w ρ ) wg z Conservation of polymer component: t [ φ(1 Sipv)cb ws + ρ r C a (1 φ) ] + (cb w u p ) = qp
6 Model for polymer flooding Simple model: introduce extra immiscible component and mixture law Polymer transported in water: u w = viscosity enhancement k rw(s) µ w,eff (c) R k (c) K( p w ρ ) wg z permeability reduction Todd Longstaff mixing: 1 µ = 1 c/cm w,eff µm(c) ω µ 1 ω w + c/cm µm(c) ω µ 1 ω p u p = krw(s) µ p,eff (c)r k (c) K( p w ρ ) wg z Conservation of polymer component: t [ φ(1 Sipv)cb ws + ρ r C a (1 φ) ] + (cb w u p ) = qp
7 Model for polymer flooding Simple model: introduce extra immiscible component and mixture law Polymer transported in water: u w = viscosity enhancement k rw(s) µ w,eff (c) R k (c) K( p w ρ ) wg z permeability reduction Todd Longstaff mixing: 1 µ = 1 c/cm w,eff µm(c) ω µ 1 ω w + c/cm µm(c) ω µ 1 ω p u p = krw(s) µ p,eff (c)r k (c) K( p w ρ ) wg z Conservation of polymer component: t [ φ(1 Sipv)cb ws + ρ r C a (1 φ) ] + (cb w u p ) = qp inaccessible pore space
8 Model for polymer flooding Simple model: introduce extra immiscible component and mixture law Polymer transported in water: u w = viscosity enhancement k rw(s) µ w,eff (c) R k (c) K( p w ρ ) wg z permeability reduction Todd Longstaff mixing: 1 µ = 1 c/cm w,eff µm(c) ω µ 1 ω w + c/cm µm(c) ω µ 1 ω p u p = krw(s) µ p,eff (c)r k (c) K( p w ρ ) wg z Conservation of polymer component: t [ φ(1 Sipv)cb ws + ρ r C a (1 φ) ] + (cb w u p ) = qp inaccessible pore space adsorption
9 Numerical framework Dicrete cell-averages ( q = p, s, c ) q i (t) = 1 ZZ q(x, y, t)dx i i Discretized eq. for water: [ w (p i ) (p i )s i ] n+1 =[ w (p i ) (p i )s i ] n t X Z [ w (p) w (s, c)] m ij i (vm w n) ij ds i j=1 ij i j
10 Numerical framework Dicrete cell-averages ( q = p, s, c ) q i (t) = 1 ZZ q(x, y, t)dx i i Discretized eq. for water: [ w (p i ) (p i )s i ] n+1 =[ w (p i ) (p i )s i ] n t X Z [ w (p) w (s, c)] m ij i (vm w n) ij ds i j=1 ij i j How to compute the mass flux evaluated at the interface? Z ij (p) ij (s, c)(v n) ij ds 1 2 [ (p i)+ (p j )] v ij Z (s, c)ds ij ij
11 Which quadrature rule to use for the integral? How to reconstruct the necessary one-sided point values from the cell-averages? How to approximate the mobility given point values that generally are different on each side of? ij
12 Which quadrature rule to use for the integral? ØMidpoint, Simpsons rule, 4 th order Gauss How to reconstruct the necessary one-sided point values from the cell-averages? How to approximate the mobility given point values that generally are different on each side of ij? ( (sij,c ij ), if v ij 0, ij(s, c) = (s + ij,c+ ij ), otherwise.
13 First-order scheme: constant reconstruction q(x) =q i Second-order slope-limiter scheme: linear reconstruction q(x) =q i + i x (x x i )+ y i (y y i) with slopes estimated from cell-averages x i x = q i q i (1,0),q i+(1,0) q i
14 First-order scheme: constant reconstruction q(x) =q i Second-order slope-limiter scheme: linear reconstruction q(x) =q i + i x (x x i )+ y i (y y i) with slopes estimated from cell-averages x i x = q i q i (1,0),q i+(1,0) q i 3 vanleer superbee minmod 5 (1,b) = b
15 WENO schemes WENO builds on ENO, which picks the least oscillatory stencil, uses a convex combination of local stencils Here: simplified version, with four linear reconstructions Polynomial: q NE (x) =q i + E i (x x i )+ N i (y y i ) Smoothness indicator: NE i = 1 h i x E 2 i + y N 2 l 4 i + Weights: w NE = NE / Linear combination: q i (x) = X v=ne,nw,se,sw NE + NW + SE + SW w v q v (x) E i = q i+(1,0) q i x N i = q i+(0,1) q i y
16 Solving the discretized systems System on residual form: F(y) =0 y =[p, s, c] Use Newton-Raphson method: y = y 0 + y 0 F(y 0 )+J y, J = df dy Structure of the Jacobian matrix J = R R R R R R R R R
17 Structure of the Jacobian matrix for the implicit WENO w Second order: WENO J = R R R R R R R R R 3 7
18 Structure of the Jacobian matrix for the implicit WENO w Second order: WENO When using lagged evaluation of the WENO weights, the structure of the Jacobian
19 A motivating example Consider φq t + q x = 0
20 A motivating example Consider φq t + q x = 0 τ=xφ transformed q t + q τ = 0 φ M
21 A motivating example Consider φq t + q x = 0 τ=xφ transformed q t + q τ = 0 φ M Modified equation (implicit/explicit first order) q t + q τ = 1 2 ( τ ± t)q ττ Smearing of discontinuity O ( t( τ ± t))
22 A motivating example Consider φq t + q x = 0 τ=xφ transformed q t + q τ = 0 φ M Modified equation (implicit/explicit first order) q t + q τ = 1 2 ( τ ± t)q ττ Smearing of discontinuity O ( t( τ ± t)) Hence: 9φ( ) xφ ± t + φ } 10 {{} 10M high-porosity region ( xφ M ± t ) } {{ } low-porosity region
23 A motivating example Consider φq t + q x = 0 τ=xφ transformed q t + q τ = 0 φ M Modified equation (implicit/explicit first order) q t + q τ = 1 2 ( τ ± t)q ττ Smearing of discontinuity O ( t( τ ± t)) Hence: 9φ( ) xφ ± t + φ } 10 {{} 10M high-porosity region ( xφ M ± t ) } {{ } low-porosity region = 9 xφ ( ν ) φ 2 x ( ) 1± + 1±ν 10 M 10M 2
24 1D: pure waterflooding with regions of different porosity = 0.01, if x 2 [0.1L, 0.2L], 0.2, otherwise. t = 8 < : 0.01, - 0.1, -.- 1, (exceeds the CFL condition in the low-porosity region) Explicit first-order scheme Implicit first-order scheme t = T t =8T t =8T t =3T t =3T t = T 0.1L 0.2L 0.1L 0.2L Explicit: oscillations (disappear as time evolves) Implicit: stable solutions for both time steps
25 1D: polymer flooding with regions of different porosity 0.8 = 0.01, if x 2 [0.1L, 0.2L], 0.2, otherwise. Saturation Explicit 1 st - order Implicit 1 st - order Implicit 2 nd - order / / Concentration t =0.01 t =0.1/ 1 1 t =5t t =5t t = T t = T 0 0.1L 0.2L 0.1L 0.2L the high-order implicit scheme resolves the displacement fronts as good as the explicit scheme the low-order implicitscheme fails to sharply resolve the structure of the oil bank that arises 0
26 1D: polymer flooding with regions of different porosity Comparison of the computational costs:
27 1D: polymer flooding with regions of different porosity Comparison of the computational costs: oscillations in the low porosity region the second-order implicit scheme is at least as effective as the explicit schemes
28 Implicit (50 control steps) Approximated solution with different reconstructions, on a uniform Cartesian grid with 50x50 cells Explicit (1600 control steps) constant: 93.2 sec minmod: 98.1 sec vanleer: sec superbee: sec weno: 111 sec constant: 9.1 sec minmod: 23.5 sec vanleer: 31.5 sec superbee: 71.2 sec weno: 42.2 sec constant: 9.1 sec minmod: 13.5 sec vanleer: 13.5 sec superbee: 13.4 sec weno: 16.6 sec
29 Implicit (50 control steps) Approximated solution with different reconstructions, on a uniform Cartesian grid with 50x50 cells Explicit (1600 control steps) constant: 93.2 sec minmod: 98.1 sec vanleer: sec superbee: sec weno: 111 sec constant: 9.1 sec minmod: 23.5 sec vanleer: 31.5 sec superbee: 71.2 sec weno: 42.2 sec constant: 9.1 sec minmod: 13.5 sec vanleer: 13.5 sec superbee: 13.4 sec weno: 16.6 sec
30 Approximated solution with different reconstructions, on a uniform Cartesian grid with 50x50 cells Explicit (1600 control steps) constant: 93.2 sec minmod: 98.1 sec vanleer: sec superbee: sec weno: 111 sec constant: 9.1 sec minmod: 23.5 sec vanleer: 31.5 sec superbee: 71.2 sec weno: 42.2 sec Implicit (50 control steps) steeper slopes constant: 9.1 sec minmod: 13.5 sec vanleer: 13.5 sec superbee: 13.4 sec weno: 16.6 sec The overall system more nonlinear and coupled
31 Implicit (50 control steps) Approximated solution with different reconstructions, on a uniform Cartesian grid with 50x50 cells Explicit (1600 control steps) constant: 93.2 sec minmod: 98.1 sec vanleer: sec superbee: sec weno: 111 sec constant: 9.1 sec minmod: 23.5 sec vanleer: 31.5 sec superbee: 71.2 sec weno: 42.2 sec constant: 9.1 sec minmod: 13.5 sec vanleer: 13.5 sec superbee: 13.4 sec weno: 16.6 sec
32 Implicit (50 control steps) Approximated solution with different reconstructions, on a uniform Cartesian grid with 50x50 cells Explicit (1600 control steps) constant: 93.2 sec minmod: 98.1 sec vanleer: sec superbee: sec weno: 111 sec constant: 9.1 sec minmod: 23.5 sec vanleer: 31.5 sec superbee: 71.2 sec weno: 42.2 sec constant: 9.1 sec minmod: 13.5 sec vanleer: 13.5 sec superbee: 13.4 sec weno: 16.6 sec Using a second-order reconstruction and improved spatial quadrature gives more accurate solution profiles both for the explicit and the implicit schemes (as expected)
33 Implicit (50 control steps) Implicit* (50 control steps) Approximated solution with different reconstructions, on a uniform Cartesian grid with 50x50 cells constant: 9.1 sec minmod: 23.5 sec vanleer: 31.5 sec superbee: 71.2 sec weno: 42.2 sec constant: 9.1 sec minmod: 13.5 sec vanleer: 13.5 sec superbee: 13.4 sec weno: 16.6 sec
34 Implicit (50 control steps) Implicit* (50 control steps) Approximated solution with different reconstructions, on a uniform Cartesian grid with 50x50 cells constant: 9.1 sec minmod: 23.5 sec vanleer: 31.5 sec superbee: 71.2 sec weno: 42.2 sec constant: 9.1 sec minmod: 13.5 sec vanleer: 13.5 sec superbee: 13.4 sec weno: 16.6 sec Implicit *: lagged evaluation of slope limiters and WENO weights reduces the number of nonlinear iterations no adverse effect on the stability and only reduces accuracy slightly
35 Number of iterations required by the explicit schemes constant minmod vanleer superbee weno
36 Number of iterations required by the explicit schemes constant minmod 0 steeper slopes vanleer The overall system more nonlinear and coupled 1 superbee weno
37 Number of iterations required by the implicit schemes 4 3 constant 2 1 minmod van Leer superbee WENO Control step Ok iterations Failed iterations
38 Number of iterations required by the implicit schemes constant minmod vanleer superbee weno Time [days]
39 Number of iterations required by the implicit* schemes constant minmod vanleer superbee weno Time [days]
40 Test of large-time-step capability constant constant t 8 7 minmod Implicit vanleer superbee minmod vanleer superbee d r e weno weno o constant constant CFL number t CFL number 8 7 minmod Implicit* vanleer superbee minmod vanleer superbee d r e weno weno o CFL number CFL number Simulating a single time-step, starting from a well-established displacement profile, 1 day
41 Test of large-time-step capability constant constant t 8 7 minmod Implicit vanleer superbee minmod vanleer superbee d r e weno weno o constant constant CFL number t CFL number 8 7 minmod Implicit* vanleer superbee minmod vanleer superbee d r e weno weno o CFL number CFL number Simulating a single time-step, starting from a well-established displacement profile, 1 day (left) and 5 days (right)
42 Spatial convergence 20 x x x 80 Implicit first-order constant: 1.6 sec constant: 5.9 sec constant: 26.7 sec Implicit second-order minmod: 7.9 sec minmod: 15.8 sec minmod: 74.5 sec
43 Spatial convergence 20 x x x 80 Implicit first-order constant: 1.6 sec constant: 5.9 sec constant: 26.7 sec Implicit second-order minmod: 7.9 sec minmod: 15.8 sec minmod: 74.5 sec
44 Spatial convergence s Explicit CFL=1 CFL=2 CFL=4 CFL=8 CFL=16 CFL=32 CFL=64 CFL= x
45 Temporal convergence s Explicit CFL=1 CFL=2 CFL=4 CFL=8 CFL=16 CFL=32 CFL=64 CFL= x
46 Temporal convergence s Explicit CFL=1 CFL=2 CFL=4 CFL=8 CFL=16 CFL=32 CFL=64 CFL= x
47 Temporal convergence
48 Temporal convergence Observe: can safely increase the CFL number for the implicit scheme to one order of magnitude beyond the stability limit before the numerical dissipation causes significant smearing
49 Test of grid-orientation errors Producer Injector Quarter five-spot Rotated five-spot
50 Test of grid-orientation errors Producer Injector Quarter five-spot Rotated five-spot Original setup: overestimates the frontal movement into stagnant regions, and underestimates the movement in the high-flow zones along the diagonal by smearing the tip of the finger Rotated setup: overestimates the movement of the front in the high-flow zone, and underestimates the movement towards stagnant zones
51 Test of grid-orientation errors 1 st order constant Colors: rotated setup setup Solid lines: original
52 Test of grid-orientation errors 1 st order constant 2 nd order minmod 2 nd order van Leer 2 nd order WENO Colors: rotated setup setup Solid lines: original
53 Test of grid-orientation errors 1 st order constant 2 nd order minmod 2 nd order van Leer 2 nd order WENO Colors: rotated setup setup Solid lines: original
54 Test of grid-orientation errors st order: rotated minmod: rotated van Leer: rotated st order: original minmod: original 0.2 van Leer: original Time [days] Water saturation in cell containing the producer
55 Channelized reservoir (K, )= ( (200 md, 0.25), for x 2 [100, 200] and x 2 [300, 400], (100 md, 0.20), otherwise. x
56 Channelized reservoir (K, )= ( (200 md, 0.25), for x 2 [100, 200] and x 2 [300, 400], (100 md, 0.20), otherwise. x
57 Channelized reservoir Explicit: constant Explicit: vanleer Implicit: constant Implicit: vanleer Explicit: constant Explicit: vanleer Implicit: constant Implicit: vanleer ) x = 150 2) x = 250
58 Channelized reservoir ) y = 700 4)y = 300
59 Including buoyancy Does including buoyancy make the nonlinear convergence more challenging? t =0days t = 110 days t = 1000 days 0
60 Including buoyancy Cummulative numberof iterations required by the nonlinear solver: implicit (solid) and implicit* (dashed) constant minmod vanleer superbee weno Time [days]
61 Summary Studied explicit and fully implicit schemes with second-order formal spatial accuracy applied to polymer flooding The use of high-resolution spatial stencils improves the accuracy and reduces grid-orientation effects Demonstrated that implicit time discretization is more suitable than explicit time integration through numerous numerical experiments To make high-resolution methods amenable for implicit discretizations, preference should be given to spatial stencils and nonlinear limiter functions that are as smooth as possible to avoid exacerbating the nonlinearity of the implicit flow equations
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationSindre Tonning Hilden. Upscaling of Water-Flooding Scenarios and Modeling of Polymer Flow. Doctoral theses at NTNU, 2016:12. Sindre Tonning Hilden
Doctoral theses at NTNU, 216:12 Sindre Tonning Hilden Sindre Tonning Hilden Upscaling of Water-Flooding Scenarios and Modeling of Polymer Flow ISBN 978-82-326-1368-7 (printed version) ISBN 978-82-326-1369-4
More informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationA Finite Volume Code for 1D Gas Dynamics
A Finite Volume Code for 1D Gas Dynamics Michael Lavell Department of Applied Mathematics and Statistics 1 Introduction A finite volume code is constructed to solve conservative systems, such as Euler
More informationField Scale Modeling of Local Capillary Trapping during CO 2 Injection into the Saline Aquifer. Bo Ren, Larry Lake, Steven Bryant
Field Scale Modeling of Local Capillary Trapping during CO 2 Injection into the Saline Aquifer Bo Ren, Larry Lake, Steven Bryant 2 nd Biennial CO 2 for EOR as CCUS Conference Houston, TX October 4-6, 2015
More informationOn efficient implicit upwind schemes
On efficient implicit upwind schemes Jostein R. Natvig Knut Andreas Lie June 30, 2008 Abstract Although many advanced methods have been devised for the hyperbolic transport problems that arise in reservoir
More informationDerivation of the fractional flow equation for a one-dimensional oil-water system. Consider displacement of oil by water in a system of dip angle α
TPG45 Reservoir Recovery Techniques 27 /9 BUCKLEY-LEVERETT ANALYSIS Derivation of the fractional flow equation for a one-dimensional oil-water system Consider displacement of oil by water in a system of
More informationTask 6: Reservoir simulation tools.
Task 6: Reservoir simulation tools Randi.Valestrand@iris.no Theme 2: Tasks and Task Leaders Task6 Robert Köfkorn Senior Researcher IRIS Robert.Kloefkorn@iris.no DEVELOPMENT & IMPROVED METHOLOGY FOR FULL
More informationORDERING-BASED NONLINEAR SOLVER WITH ADAPTIVE COUPLING FOR MULTIPHASE FLOW AND TRANSPORT
ORDERING-BASED NONLINEAR SOLVER WITH ADAPTIVE COUPLING FOR MULTIPHASE FLOW AND TRANSPORT A DISSERTATION SUBMITTED TO THE DEPARTMENT OF ENERGY RESOURCES ENGINEERING AND THE COMMITTEE ON GRADUATE STUDIES
More informationThe RAMSES code and related techniques I. Hydro solvers
The RAMSES code and related techniques I. Hydro solvers Outline - The Euler equations - Systems of conservation laws - The Riemann problem - The Godunov Method - Riemann solvers - 2D Godunov schemes -
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 informationHistory matching of experimental and CMG STARS results
https://doi.org/1.17/s13-1-55- ORIGINAL PAPER - PRODUCTION ENGINEERING History matching of experimental and CMG STARS results Ahmed Tunnish 1 Ezeddin Shirif 1 Amr Henni Received: 1 February 17 / Accepted:
More informationIMPERIAL COLLEGE LONDON. Upscaling polymer flooding in heterogeneous reservoirs
IMPERIAL COLLEGE LONDON Department of Earth Science and Engineering Centre for Petroleum Studies Upscaling polymer flooding in heterogeneous reservoirs By Jun Lee A report submitted in partial fulfilment
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
More informationFlow an open source research tool for reservoir simulation
Flow an open source research tool for reservoir simulation Robert Klöfkorn IRIS & The National IOR Centre of Norway April 27, 2016 Task 6 Reservoir Simulation Tools Pål (PostDoc UiS) Trine (PostDoc IRIS)
More informationOn the Comparison of the Finite Volume and Discontinuous Galerkin Methods
Diploma Thesis Institute for Numerical Simulation, TUHH On the Comparison of the Finite Volume and Discontinuous Galerkin Methods Corrected version Katja Baumbach August 17, 2006 Supervisor: Prof. Dr.
More informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationIntroduction 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 informationChapter 6. Finite Element Method. Literature: (tiny selection from an enormous number of publications)
Chapter 6 Finite Element Method Literature: (tiny selection from an enormous number of publications) K.J. Bathe, Finite Element procedures, 2nd edition, Pearson 2014 (1043 pages, comprehensive). Available
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationSolution of Two-Dimensional Riemann Problems for Gas Dynamics without Riemann Problem Solvers
Solution of Two-Dimensional Riemann Problems for Gas Dynamics without Riemann Problem Solvers Alexander Kurganov, 1, * Eitan Tadmor 2 1 Department of Mathematics, University of Michigan, Ann Arbor, Michigan
More informationGODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS
GODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS ARTURO S. LEON Dept. of Civil and Envir. Engng., Univ. of Illinois at Urbana-Champaign, 2519 Hydrosystems Lab., MC-250. 205 North Mathews Av., Urbana,
More informationEfficient and Accurate simulation of nonlinearly coupled multiphase flow in porous media. 25 August 2015 Meeneesh Kardale
Efficient and Accurate simulation of nonlinearly coupled multiphase flow in porous media 25 August 2015 Meeneesh Kardale Title : Efficient and accurate simulation of nonlinearly coupled multiphase flow
More informationNon-linear Methods for Scalar Equations
Non-linear Methods for Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 3, 04 / 56 Abstract
More informationA PSEUDO FUNCTION APPROACH IN RESERVOIR SIMULATION
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 2, Supp, Pages 58 67 c 2005 Institute for Scientific Computing and Information A PSEUDO FUNCTION APPROACH IN RESERVOIR SIMULATION ZHANGXIN
More informationModeling of two-phase flow in fractured porous media on unstructured non-uniform coarse grids
Modeling of two-phase flow in fractured porous media on unstructured non-uniform coarse grids Jørg Espen Aarnes and Vera Louise Hauge SINTEF ICT, Deptartment of Applied Mathematics Applied Mathematics
More informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationTwo Phase Transport in Porous Media
Two Phase Transport in Porous Media Lloyd Bridge Iain Moyles Brian Wetton Mathematics Department University of British Columbia www.math.ubc.ca/ wetton CRM CAMP Seminar, October 19, 2011 Overview Two phase
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationChapter Seven. For ideal gases, the ideal gas law provides a precise relationship between density and pressure:
Chapter Seven Horizontal, steady-state flow of an ideal gas This case is presented for compressible gases, and their properties, especially density, vary appreciably with pressure. The conditions of the
More informationIntroduction to Reservoir Simulation as Practiced in Industry
Introduction to Reservoir Simulation as Practiced in Industry Knut Andreas Lie SINTEF Digital, Mathematics and Cybernetics / Department of Mathematical Sciences, NTNU, Norway Computational Issues in Oil
More informationScalable Non-Linear Compact Schemes
Scalable Non-Linear Compact Schemes Debojyoti Ghosh Emil M. Constantinescu Jed Brown Mathematics Computer Science Argonne National Laboratory International Conference on Spectral and High Order Methods
More informationFDM for wave equations
FDM for wave equations Consider the second order wave equation Some properties Existence & Uniqueness Wave speed finite!!! Dependence region Analytical solution in 1D Finite difference discretization Finite
More informationA Higher-Order Finite Volume Method for Non-Conservative Equations on 2D Unstructured Quadrilateral Grids
Diese Arbeit wurde vorgelegt am Lehrstuhl für Mathematik (MathCCES) A Higher-Order Finite Volume Method for Non-Conservative Equations on 2D Unstructured Quadrilateral Grids CES Seminar Computational Engineering
More informationMultigrid solvers for equations arising in implicit MHD simulations
Multigrid solvers for equations arising in implicit MHD simulations smoothing Finest Grid Mark F. Adams Department of Applied Physics & Applied Mathematics Columbia University Ravi Samtaney PPPL Achi Brandt
More informationComparison of cell-centered and node-centered formulations of a high-resolution well-balanced finite volume scheme: application to shallow water flows
Comparison of cell-centered and node-centered formulations of a high-resolution well-balanced finite volume scheme: application to shallow water flows Dr Argiris I. Delis Dr Ioannis K. Nikolos (TUC) Maria
More informationHyemin Park, Jinju Han, Wonmo Sung*
Experimental Investigation of Polymer Adsorption-Induced Permeability Reduction in Low Permeability Reservoirs 2014.10.28 Hyemin Park, Jinju Han, Wonmo Sung* Hanyang Univ., Seoul, Rep. of Korea 1 Research
More informationConservation Laws & Applications
Rocky Mountain Mathematics Consortium Summer School Conservation Laws & Applications Lecture V: Discontinuous Galerkin Methods James A. Rossmanith Department of Mathematics University of Wisconsin Madison
More informationIMPERIAL COLLEGE LONDON. Department of Earth Science and Engineering. Centre for Petroleum Studies
IMPERIAL COLLEGE LONDON Department of Earth Science and Engineering Centre for Petroleum Studies Taking Chemical EOR Data From Lab Results to Simulation By Laurent Libert A report submitted in partial
More informationIMPERIAL COLLEGE LONDON. Department of Earth Science and Engineering. Centre for Petroleum Studies
IMPERIAL COLLEGE LONDON Department of Earth Science and Engineering Centre for Petroleum Studies Upscaling of Relative Permeability to Minimise Numerical Dispersion By Anthony Tobechukwu Afoaku A report
More informationHigh Order Semi-Lagrangian WENO scheme for Vlasov Equations
High Order WENO scheme for Equations Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Andrew Christlieb Supported by AFOSR. Computational Mathematics Seminar, UC Boulder
More informationStreamline calculations. Lecture note 2
Streamline calculations. Lecture note 2 February 26, 2007 1 Recapitulation from previous lecture Definition of a streamline x(τ) = s(τ), dx(τ) dτ = v(x,t), x(0) = x 0 (1) Divergence free, irrotational
More informationA Numerical Simulation Framework for the Design, Management and Optimization of CO 2 Sequestration in Subsurface Formations
A Numerical Simulation Framework for the Design, Management and Optimization of CO 2 Sequestration in Subsurface Formations Investigators Hamdi Tchelepi, Associate Professor of Energy Resources Engineering;
More informationThe effect of heterogeneity on unsteady-state displacements
The effect of heterogeneity on unsteady-state displacements Abstract Darryl Fenwick, Nicole Doerler, and Roland Lenormand, Institut Français du Pétrole In this paper, we discuss the effect of heterogeneity
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 informationMulti-rate mass transfer modeling of two-phase flow in highly heterogeneous fractured and porous media
Multi-rate mass transfer modeling of two-phase flow in highly heterogeneous fractured and porous media Jan Tecklenburg a,, Insa Neuweiler a, Jesus Carrera b, Marco Dentz b a Institute of Fluid Mechanics
More informationA Fourth-Order Central Runge-Kutta Scheme for Hyperbolic Conservation Laws
A Fourth-Order Central Runge-Kutta Scheme for Hyperbolic Conservation Laws Mehdi Dehghan, Rooholah Jazlanian Department of Applied Mathematics, Faculty of Mathematics and Computer Science, Amirkabir University
More informationA FRONT-TRACKING METHOD FOR HYPERBOLIC THREE-PHASE MODELS
A FRONT-TRACKING METHOD FOR HYPERBOLIC THREE-PHASE MODELS Ruben Juanes 1 and Knut-Andreas Lie 2 1 Stanford University, Dept. Petroleum Engineering, USA 2 SINTEF IKT, Dept., Norway ECMOR IX, August 30 September
More informationPartial Saturation Fluid Substitution with Partial Saturation
Fluid Substitution with 261 5 4.5 ρ fluid S w ρ w + S o ρ o + S g ρ g Vp (km/s) 4 3.5 K fluid S w K w + S o K o + S g K g Patchy Saturation Drainage 3 2.5 2 Fine-scale mixing 1 = S w + S o + S g K fluid
More informationOn limiting for higher order discontinuous Galerkin method for 2D Euler equations
On limiting for higher order discontinuous Galerkin method for 2D Euler equations Juan Pablo Gallego-Valencia, Christian Klingenberg, Praveen Chandrashekar October 6, 205 Abstract We present an implementation
More informationCONSTRUCTION OF HIGH-ORDER ADAPTIVE IMPLICIT METHODS FOR RESERVOIR SIMULATION
CONSTRUCTION OF HIGH-ORDER ADAPTIVE IMPLICIT METHODS FOR RESERVOIR SIMULATION A REPORT SUBMITTED TO THE DEPARTMENT OF ENERGY RESOURCES ENGINEERING OF STANFORD UNIVERSITY IN PARTIAL FULFILLMENT OF THE REQUIREMENTS
More informationNUMERICAL TREATMENT OF TWO-PHASE FLOW IN CAPILLARY HETEROGENEOUS POROUS MEDIA BY FINITE-VOLUME APPROXIMATIONS
NUMERICAL TREATMENT OF TWO-PHASE FLOW IN CAPILLARY HETEROGENEOUS POROUS MEDIA BY FINITE-VOLUME APPROXIMATIONS HELMER ANDRÉ FRIISA AND STEINAR EVJE B,C Abstract. This paper examines two-phase flow in porous
More informationICES REPORT Adaptive Numerical Homogenization for Non-Linear Multiphase Flow and Transport
ICES REPORT 17-13 June 217 Adaptive Numerical Homogenization for Non-Linear Multiphase Flow and Transport by Gurpreet Singh, Yerlan Amanbek, and Mary F. Wheeler The Institute for Computational Engineering
More informationarxiv: v3 [math.na] 10 Mar 2015
A Posteriori Subcell Limiting of the Discontinuous Galerkin Finite Element Method for Hyperbolic Conservation Laws Michael Dumbser a, Olindo Zanotti a, Raphaël Loubère b, Steven Diot c a Department of
More informationTodd Arbogast. Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and Sciences (ICES)
CONSERVATIVE CHARACTERISTIC METHODS FOR LINEAR TRANSPORT PROBLEMS Todd Arbogast Department of Mathematics and, (ICES) The University of Texas at Austin Chieh-Sen (Jason) Huang Department of Applied Mathematics
More informationProceedings of the ASME nd International Conference on Ocean, Offshore and Arctic Engineering OMAE2013 June 9-14, 2013, Nantes, France
Proceedings of the ASME 2 2nd International Conference on Ocean, Offshore and Arctic Engineering OMAE2 June 9-4, 2, Nantes, France OMAE2-5 FLUID DYNAMICAL AND MODELING ISSUES OF CHEMICAL FLOODING FOR ENHANCED
More informationGodunov methods in GANDALF
Godunov methods in GANDALF Stefan Heigl David Hubber Judith Ngoumou USM, LMU, München 28th October 2015 Why not just stick with SPH? SPH is perfectly adequate in many scenarios but can fail, or at least
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 informationA Third Order Fast Sweeping Method with Linear Computational Complexity for Eikonal Equations
DOI.7/s95-4-9856-7 A Third Order Fast Sweeping Method with Linear Computational Complexity for Eikonal Equations Liang Wu Yong-Tao Zhang Received: 7 September 23 / Revised: 23 February 24 / Accepted: 5
More informationShock Capturing for Discontinuous Galerkin Methods using Finite Volume Sub-cells
Abstract We present a shock capturing procedure for high order Discontinuous Galerkin methods, by which shock regions are refined in sub-cells and treated by finite volume techniques Hence, our approach
More informationA BENCHMARK CALCULATION OF 3D HORIZONTAL WELL SIMULATIONS
INTERNATINAL JURNAL F NUMERICAL ANALYSIS AND MDELING Volume 1, Number 2, Pages 189 201 c 2004 Institute for Scientific Computing and Information A BENCHMARK CALCULATIN F 3D HRIZNTAL WELL SIMULATINS ZHANGIN
More informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
More informationA Fifth Order Flux Implicit WENO Method
A Fifth Order Flux Implicit WENO Method Sigal Gottlieb and Julia S. Mullen and Steven J. Ruuth April 3, 25 Keywords: implicit, weighted essentially non-oscillatory, time-discretizations. Abstract The weighted
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationA high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows
A high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows F. Bassi, A. Crivellini, D. A. Di Pietro, S. Rebay Dipartimento di Ingegneria Industriale, Università di Bergamo CERMICS-ENPC
More informationCENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION
CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. The central scheme of
More informationDownloaded 08/04/16 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 38, No. 4, pp. B593 B618 c 2016 Society for Industrial and Applied Mathematics ACTIVE-SET REDUCED-SPACE METHODS WITH NONLINEAR ELIMINATION FOR TWO-PHASE FLOW PROBLEMS IN POROUS
More informationLES Simulations of Quiet Sun Magnetism
LES Simulations of Quiet Sun Magnetism Matthias Rempel HAO/NCAR Quiet sun magnetism Origin and spatial distribution of quiet sun field Small scale dynamo? Remnant field from large scale dynamo? Vögler,
More informationUpscaling for Two-phase Flows in Porous Media
Upscaling for Two-phase Flows in Porous Media Thesis by Andrew Westhead In Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy California Institute of Technology Pasadena, California
More informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationWaterflooding Performance of Communicating Stratified Reservoirs With Log-Normal Permeability Distribution
Waterflooding Performance of Communicating Stratified Reservoirs With Log-Normal Permeability Distribution Noaman El-Khatib, SPE, King Saud U. Summary An analytical solution is developed for waterflooding
More informationSimulation of Carbonated Water Injection (CWI), Challenges and Solution
Simulation of Carbonated Water Injection (CWI), Challenges and Solution Foroozesh J., Jamiolahmady M. (Jami), Sohrabi M. Presented by: Jalal Foroozesh 1 Outline Introduction Objectives Our modelling approach
More informationarxiv: v4 [physics.flu-dyn] 23 Mar 2019
High-Order Localized Dissipation Weighted Compact Nonlinear Scheme for Shock- and Interface-Capturing in Compressible Flows Man Long Wong a, Sanjiva K. Lele a,b arxiv:7.895v4 [physics.flu-dyn] 3 Mar 9
More informationHydraulic Modelling for Drilling Automation
Hydraulic Modelling for Drilling Automation CASA Day Harshit Bansal April 19, 2017 Where innovation starts Team: Supervisors at TU/e : W.H.A. Schilders, N. van de Wouw, B. Koren, L. Iapichino Collaborators:
More informationThe WENO Method for Non-Equidistant Meshes
The WENO Method for Non-Equidistant Meshes Philip Rupp September 11, 01, Berlin Contents 1 Introduction 1.1 Settings and Conditions...................... The WENO Schemes 4.1 The Interpolation Problem.....................
More informationWeighted ENO Schemes
Xiaolei Chen Advisor: Prof. Xiaolin Li Department of Applied Mathematics and Statistics Stony Brook University, The State University of New York February 7, 014 1 3 Mapped WENO-Z Scheme 1D Scalar Hyperbolic
More informationComplexity of Two-Phase Flow in Porous Media
1 Complexity of Two-Phase Flow in Porous Media Rennes September 16, 2009 Eyvind Aker Morten Grøva Henning Arendt Knudsen Thomas Ramstad Bo-Sture Skagerstam Glenn Tørå Alex Hansen 2 Declining oil production,
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 informationPolymer flooding improved sweep efficiency for utilizing IOR potential Force seminar April April 2016
Polymer flooding improved sweep efficiency for utilizing IOR potential Force seminar April 2016 Classic polymer screening Viscosifying effect Solution preparation Bulk rheology Flow properties in porous
More informationDiscretisation of the pseudo-incompressible Equations with implicit LES for Gravity Wave Simulation DACH2010. Felix Rieper, Ulrich Achatz
Discretisation of the pseudo-incompressible Equations with implicit LES for Gravity Wave Simulation DACH2010, Ulrich Achatz Institut für Atmosphäre und Umwelt Aim of the project Analysis: Propagation and
More informationICES REPORT An Approximate Jacobian Nonlinear Solver for Multiphase Flow and Transport
ICES REPORT 7-22 September 207 An Approximate Jacobian Nonlinear Solver for Multiphase Flow and Transport by Gurpreet Singh, Gergina Pencheva, Mary F. Wheeler The Institute for Computational Engineering
More informationHybrid LES RANS Method Based on an Explicit Algebraic Reynolds Stress Model
Hybrid RANS Method Based on an Explicit Algebraic Reynolds Stress Model Benoit Jaffrézic, Michael Breuer and Antonio Delgado Institute of Fluid Mechanics, LSTM University of Nürnberg bjaffrez/breuer@lstm.uni-erlangen.de
More informationModeling Multiphase Flow in Porous Media with Complementary Constraints
Background Modeling Multiphase Flow in Porous Media with Complementary Constraints Applied Math, Statistics, and Scientific Computation, University of Maryland - College Park October 07, 2014 Advisors
More informationSECOND ORDER ALL SPEED METHOD FOR THE ISENTROPIC EULER EQUATIONS. Min Tang. (Communicated by the associate editor name)
Manuscript submitted to AIMS Journals Volume X, Number X, XX 2X Website: http://aimsciences.org pp. X XX SECOND ORDER ALL SPEED METHOD FOR THE ISENTROPIC EULER EQUATIONS Min Tang Department of Mathematics
More informationA central discontinuous Galerkin method for Hamilton-Jacobi equations
A central discontinuous Galerkin method for Hamilton-Jacobi equations Fengyan Li and Sergey Yakovlev Abstract In this paper, a central discontinuous Galerkin method is proposed to solve for the viscosity
More informationHybrid semi-lagrangian finite element-finite difference methods for the Vlasov equation
Numerical Analysis and Scientific Computing Preprint Seria Hybrid semi-lagrangian finite element-finite difference methods for the Vlasov equation W. Guo J. Qiu Preprint #21 Department of Mathematics University
More informationICES REPORT A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws
ICES REPORT 7- August 7 A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws by Todd Arbogast, Chieh-Sen Huang, and Xikai Zhao The Institute for Computational Engineering and Sciences The
More informationc 2003 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol., No., pp. 57 7 c 00 Society for Industrial and Applied Mathematics ON THE ARTIFICIAL COMPRESSION METHOD FOR SECOND-ORDER NONOSCILLATORY CENTRAL DIFFERENCE SCHEMES FOR SYSTEMS
More informationAn asymptotic preserving unified gas kinetic scheme for the grey radiative transfer equations
An asymptotic preserving unified gas kinetic scheme for the grey radiative transfer equations Institute of Applied Physics and Computational Mathematics, Beijing NUS, Singapore, March 2-6, 2015 (joint
More informationA parametrized maximum principle preserving flux limiter for finite difference RK-WENO schemes with applications in incompressible flows.
A parametrized maximum principle preserving flux limiter for finite difference RK-WENO schemes with applications in incompressible flows Tao Xiong Jing-ei Qiu Zhengfu Xu 3 Abstract In Xu [] a class of
More informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
More informationDetonation Simulations with a Fifth-Order TENO Scheme
Commun. Comput. Phys. doi: 10.4208/cicp.OA-2018-0008 Vol. x, No. x, pp. 1-37 xxx 201x Detonation Simulations with a Fifth-Order TENO Scheme Haibo Dong 1, Lin Fu 2, Fan Zhang 3,, Yu Liu 4 and Jun Liu 1
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 informationImproved Seventh-Order WENO Scheme
48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Eposition 4-7 January 2, Orlando, Florida AIAA 2-45 Improved Seventh-Order WENO Scheme Yiqing Shen Gecheng Zha Dept.
More informationSpace-time domain decomposition for a nonlinear parabolic equation with discontinuous capillary pressure
Space-time domain decomposition for a nonlinear parabolic equation with discontinuous capillary pressure Elyes Ahmed, Caroline Japhet, Michel Kern INRIA Paris Maison de la Simulation ENPC University Paris
More informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
More information30 crete maximum principle, which all imply the bound-preserving property. But most
3 4 7 8 9 3 4 7 A HIGH ORDER ACCURATE BOUND-PRESERVING COMPACT FINITE DIFFERENCE SCHEME FOR SCALAR CONVECTION DIFFUSION EQUATIONS HAO LI, SHUSEN XIE, AND XIANGXIONG ZHANG Abstract We show that the classical
More informationIMPROVED STEAMFLOOD ANALYTICAL MODEL. A Thesis SUANDY CHANDRA
IMPROVED STEAMFLOOD ANALYTICAL MODEL A Thesis by SUANDY CHANDRA Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the requirements for the degree of MASTER OF
More information