IMPROVING THE NUMERICAL ACCURACY OF HYDROTHERMAL RESERVOIR SIMULATIONS USING THE CIP SCHEME WITH THIRD-ORDER ACCURACY
|
|
- Garry Ramsey
- 6 years ago
- Views:
Transcription
1 PROCEEDINGS, Thirty-Seventh Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 30 - February 1, 2012 SGP-TR-194 IMPROVING THE NUMERICAL ACCURACY OF HYDROTHERMAL RESERVOIR SIMULATIONS USING THE CIP SCHEME WITH THIRD-ORDER ACCURACY Mitsuo Matsumoto Energy & Mineral Resources Department, Idemitsu Kosan Co., Ltd , Toranomon, Minato-ku, Tokyo , Japan mitsuo.matsumoto@si.idemitsu.co.jp ABSTRACT We have developed a numerical algorithm using the CIP scheme, or its enhanced version, the RCIP scheme, for improving the numerical accuracy when simulating fluid flows with a high gradient of temperature or chemical concentration in hydrothermal reservoirs. The CIP and RCIP schemes achieve third-order accuracy, and they can effectively control the numerical diffusion that often causes serious numerical errors when using the conventional first-order upstream-difference scheme. The developed algorithm consists of the following two procedures: 1. Time integration of the equations governing the conservation of mass and enthalpy; 2. Determining the thermodynamic functions and their derivatives by solving the systems of nonlinear and linear equations. This paper presents practical numerical techniques by describing the numerical algorithms using not only the CIP or RCIP scheme but also the conventional first-order upstream-difference scheme, for reference purposes. Using these algorithms, several numerical simulations of one-dimensional cold sweep are demonstrated. From the numerical solutions, we conclude that: 1. the algorithm using the CIP or RCIP scheme can effectively control the numerical diffusion; 2. the algorithm using the RCIP scheme yields more realistic solutions by controlling the overshoot that is caused when using the CIP scheme; 3. spatial steps that are four to eight times longer can be used without compromising numerical accuracy, by replacing the conventional firstorder upstream-difference scheme with the RCIP scheme. INTRODUCTION Subsurface fluid flows with a high gradient of temperature or chemical concentration often appear in a hydrothermal reservoir. One example is a cooling front propagated by the injection of cold water into the reservoir cold sweep. Accurate simulations of this effect are important for predicting the productivity of reservoirs during the operation of geothermal power plants. Another example is the flow of a tracer injected into the reservoir. Because the tracer did not exist originally, the injected tracer forms a sharp concentration peak in the reservoir. Simulating the tracer flow accurately and matching it with the observed changes in tracer concentration at the observation wells provide us with useful constraints when constructing a reliable numerical model of the reservoir. Although these high-gradient flows play an important role, attempts to simulate them often entail overcoming the drawback of insufficient numerical accuracy. The transportation of enthalpy or a chemical component due to advection is governed by the advection equation. This equation is known as one of the hyperbolic-type partial differential equations that are difficult to solve accurately using conventional numerical techniques such as the firstorder upstream-difference scheme. The major origin of the numerical error is the numerical diffusion that severely blunts the sharp gradient and yields failed numerical solutions. Improving the numerical accuracy is essential for simulating high-gradient flows. We have improved the numerical accuracy when simulating the high-gradient flow in hydrothermal reservoirs by using the constrained interpolation profile CIP scheme Takewaki et al., 1985 and its enhanced version, the rational function CIP RCIP scheme Xiao et al., Compared to the conventional first-order upstream-difference scheme,
2 In an idealized hydrothermal system composed of pure water and rocks, the mass of fluid per unit volume M is expressed as follows: 1 Figure 1: Solutions of the advection equation. The red, yellow, and blue lines denote the analytic solution AS, numerical solution obtained using the CIP scheme CIP, and that obtained using the firstorder upstream-difference scheme USD, respectively. The time and spatial step lengths are 0.01 and 0.025, respectively. these schemes can control the numerical diffusion effectively by achieving third-order accuracy Figure. 1. This paper presents practical techniques for applying the CIP and RCIP schemes to hydrothermal reservoir simulations. First, we describe the equations governing the subsurface fluid flow and physical properties of water. Second, we describe numerical algorithms using not only the CIP or RCIP scheme but also the conventional first-order upstreamdifference scheme for reference purposes. Finally, we demonstrate several numerical simulations of a onedimensional problem involving cold sweep, and we show how the CIP and RCIP schemes effectively improve the numerical accuracy. GOVERNING EQUATIONS Conservation of Mass and Enthalpy We adopt the mathematical model proposed by Faust and Mercer 1979 for describing the subsurface fluid flow. This model is also adopted in the reservoir simulator HYDROTHERM Kipp et al., 2008, which is a well-known multi-purpose simulator. In this model, we regard the rocks as porous media. We assume that fluid exists in the pores of rocks. The fluid flow is governed by Darcy s law. In this paper, we consider only one-dimensional problems in the Cartesian coordinate system, and for simplicity, we assume that the physical properties of rocks are constant. The numerical techniques described in this paper, however, can be immediately applied with minor modifications to multi-dimensional problems entailing variable physical properties of rocks. where ϕ and ρ f represent the porosity of rocks and density of fluid, respectively. The total enthalpy per unit volume H, including the enthalpy of fluid and rocks, is defined as follows: 2 where h, ρ r, c r, and T represent the specific enthalpy of fluid, density of rocks, specific heat capacity of rocks, and temperature, respectively. The first and second terms of the right-hand side of Eq. 2 represent the enthalpy of fluid and that of rocks, respectively. The principles of conservation of mass and enthalpy govern the changes in the mass of fluid M and total enthalpy H, respectively, as follows: 3 4 where t, x, u, λ, q M, and h src represent the time, coordinate, volumetric flux of fluid, thermal conductivity, mass flow rate of a source, and specific enthalpy of the fluid flowing from a source into the reservoir, respectively. The volumetric flux of fluid obeys Darcy s law. Physical Properties of Water For determining the physical properties of water, we have referred to the empirical equation of state IAPWS-95 IAPWS, 1996 and the empirical equation for calculating the viscosity IAPWS, 1997, developed by the International Association for the Properties of Water and Steam IAPWS. Although these empirical equations can simulate the physical properties of water with high accuracy, the computational load is too large to calculate them in numerical simulations. To speed up the calculation of the physical properties of water, we have developed the electronic steam table GANSEKI Matsumoto, GANSEKI provides a module that can be imported to an arbitrary numerical code written in Fortran 90/95, and it approximates the pressure, specific enthalpy, and viscosity of water on the basis of the abovementioned empirical equations. A comparison of the calculations using GANSEKI and the original empirical equations is shown in Figure 2. As shown
3 Figure 2: Changes in pressure, specific enthalpy, and viscosity of water depending on density and temperature calculated using the electronic steam table GANSEKI and the empirical equations developed by IAPWS. in Figure 2, the independent variables of GANSEKI are density and temperature. When developing GANSEKI, we divided the density-temperature space into a number of rectangular cells, whose lengths are 0.1 to 46.5 kg/m 3 and 0.8 to 20.0 deg. C. At the four corners of each cell, the values of the functions to be approximated and their derivatives are memorized in the database. Then, changes in the functions in each cell are approximated using bicubic interpolation. The relative error of this approximation is at most %. The total data size of GANSEKI is 7.8 MB. The calculation load of GANSEKI is less than one-tenth of that of the original empirical equations. ALGORITHM 1: FIRST-ORDER UPSTREAM- DIFFERENCE SCHEME In this section, we describe the numerical algorithm, using the conventional first-order upstreamdifference scheme for reference purposes. As shown in a later section, the numerical solutions using this algorithm coincide with those obtained using the well-known reservoir simulator HYDROTHERM. This algorithm consists of the following two procedures: 1. Integrating the discretized forms of Eqs. 3 and 4 explicitly with time, for yielding the mass of fluid M i n+1 and total enthalpy H i n+1 at the next time step, where i and n are the indexes of time and spatial steps, respectively; n+1 2. Substituting M and H for M i and H n+1 i, respectively, in Eqs. 1 and 2 and solving the system of nonlinear equations formed by these two equations numerically for determining the thermodynamic functions of water at the next time step. Time Integration In the first procedure, we adopt the time-splitting method and integrate the advection and nonadvection terms of Eqs. 3 and 4 separately as follows: For the advection terms: 5
4 6 18 For the non-advection terms: 19 where where Δt, Δx, k, and μ represent the time step length, spatial step length, permeability of rocks, and viscosity of fluid, respectively. The viscosity of fluid is a function of fluid density and temperature, and it can be computed using the electronic steam table GANSEKI. The non-advection terms Eqs. 7 to 9 are integrated with time using the finite difference method as follows: When integrating the advection terms Eqs. 5 and 6 with time, we adopt the first-order upstreamdifference scheme as follows: with The total enthalpy H i n+1 at the next time step can be computed by summing up as follows: 24 Thus, the time integration of Eqs. 3 and 4 can be reduced to that of the advection Eqs. 5 and 6 and non-advection Eqs. 7 to 9 equations. The source of the serious numerical error due to numerical diffusion is the time integration of Eqs. 5 and 6. As shown later, we intend to overcome this problem by improving the numerical accuracy when integrating these two equations by using the CIP or RCIP scheme. We have adopted explicit time integration to simplify the numerical algorithms, although the available time step lengths are shorter by several orders of magnitude, as compared to implicit time integration. We have preferred the simplicity of the numerical algorithms to shorter calculation time for the first stage of our attempt to apply the CIP and RCIP schemes. We intend to adopt implicit time integration to reduce the calculation time after verifying the adoptability of these schemes.
5 Determining the Thermodynamic Functions In the second procedure, we determine the thermodynamic functions, such as pressure P i n+1 and temperature T i n+1, by using the mass of fluid M i n+1 and total enthalpy H i n+1 determined in the previous procedure. We assume that the fluid density ρ f,i n+1 and temperature T i n+1 are the independent unknowns and that the other variables are functions depending on these unknowns. Since the porosity of rocks ϕ is assumed to be constant, we can immediately obtain the fluid density ρ f,i n+1 from Eq. 1. Thus, we can reduce Eqs. 1 and 2 to a single nonlinear equation whose unknown is solely temperature T i n+1, as follows: where where we regard the enthalpy of fluid h as the function of the already obtained fluid density ρ f,i n+1 and temperature T i n+1. We can compute this function using the electronic steam table GANSEKI. After solving Eq. 25 using the Newton-Raphson method, pressure P i n+1 is determined using GANSEKI again. Finally, the enthalpy of fluid H f,i n+1 and that of rocks H r,i n+1 defined in Eqs. 11 and 12 are computed using the determined thermodynamic functions. ALGORITHM 2: CIP AND RCIP SCHEMES Time Integration Next, we describe the numerical algorithm using the CIP or RCIP scheme by modifying the algorithm presented in the previous section. In the first procedure, we consider not only Eqs. 3 and 4 but also the derivatives of Eqs. 3 and 4 as follows: where By using the CIP scheme, the pair of advection equations 5 and 30 can be integrated with time as follows: where By using the time-splitting method, the time integration of Eqs. 26 and 27 can be reduced to that of the following advection and non-advection equations, 45
6 46 47 equations 32 to 34 are integrated using the finite difference method as follows: The pair of advection equations 6 and 31 can be integrated in the same manner. We also consider the RCIP scheme that can control over- and undershoot effectively. To pdopt the RCIP scheme, we replace Eqs. 43, 44, 46, and 47 as follows: 60 with Note that the finite differences of Eqs. 21 to 23 in space are used for calculating α, β f, and β r in Eqs. 58 to 60. The derivative of the total enthalpy n+1 Ψ i at the next time step can be obtained by summing up as follows: Determining the Thermodynamic Functions and Their Derivatives In the second procedure, we consider not only Eqs. 1 and 2 but also the derivatives of Eqs. 1 and 2 as follows: where if Otherwise,. The non-advection equations 7 to 9 are integrated as shown in Eqs. 21 to 23. The non-advection 63 The derivatives of the enthalpy of fluid are determined using the electronic steam table GANSEKI. Since all the thermodynamic functions are determined by solving Eq. 25, we can solve the system of linear equations formed by Eqs. 62 and n+1 63 for δ i and τ n+1 i. Finally, the derivatives Ψ n+1 f,i and Ψ n+1 r,i defined in Eqs. 35 and 36 are computed using the determined thermodynamic functions and their derivatives.
7 NUMERICAL SIMULATIONS Conditions We present several numerical simulations using the abovementioned numerical algorithms. In these simulations, we considered a one-dimensional problem of cold sweep in a hydrothermal reservoir, as shown in Figure 3. We simulated the change in pressure and temperature in a one-dimensional space whose length was 2.0 km. We assumed that the physical properties of rocks are constant, as shown in Table 1. The initial pressure and temperature were uniformly 20 MPa and 250 deg. C, respectively. Cold water with a specific enthalpy of 550 kj/kg the specific enthalpy of saturated water at 131 deg. C was injected at a source located 800 m away from the left boundary. The flow rate of the source was kg/s/m 2. The pressure and temperature at the two boundaries remained constant at the initial values. Numerical Solutions and Discussions Figure 4 shows the numerical solutions of the abovementioned problem obtained using the algorithms with the first-order upstream-difference, CIP, and RCIP schemes. For reference purposes, each of these solutions overlays the numerical solution obtained using the well-known reservoir simulator HYDROTHERM under the same conditions. The common time and spatial step lengths were years and 100 m, respectively. Since the injection of cold water at the source, the cooling fronts at which the temperature begins to decrease from the initial value propagate to both sides of the one-dimensional space. Most of the space has been cooled down over 30 years. The numerical solution obtained using the algorithm with the first-order upstream-difference scheme is compared with that obtained using HYDROTHERM Figure 4a. Both solutions precisely coincide with each other. This ensures the validity of the numerical solutions using the first-order upstream-difference scheme. Figures 4b and 4c show the numerical solutions obtained using the algorithm with the CIP and RCIP schemes, respectively, with third-order accuracy. Both solutions simulate the temperature profiles at the cooling fronts more sharply than the solution with the first-order upstream-difference scheme Figure 4a by controlling the numerical diffusion effectively. The RCIP scheme yields a more realistic solution than the CIP scheme by controlling the overshoot at the cooling fronts. We have confirmed the validity of the numerical solution obtained using the RCIP scheme by studying the convergence of the numerical solutions. Figures 5 and 6 show the simulated temperature profiles using the RCIP scheme with a spatial step length of 100 m and the first-order upstream-difference scheme with Figure 3: One-dimensional problem for the simulations. Table 1: Physical properties of rocks for the simulations. Property Value Porosity 0.10 Density kg/m 3 Permeability m 2 Thermal conductivity Specific heat capacity Table 2: 1.0 W/m/K J/kg/K Pairs of time and spatial step lengths for the simulations. Time step length [years] Spatial step length [m] spatial step lengths of 12.5 to 100 m. Refer to Table 2 for the time step lengths used for various spatial step lengths. We find that the spatial and temporal profiles of temperature simulated using the first-order upstream-difference scheme converge with that simulated by using the RCIP scheme as the time and spatial step lengths decrease. From this convergence, we can conclude that the numerical solution obtained using the RCIP scheme is more accurate than that obtained using the first-order upstream-difference scheme with the common spatial step length. Figures 5 and 6 show that the RCIP scheme can yield a simulation with accuracy that is comparable with that of the first-order upstream-difference scheme, even if we adopt spatial steps that are four to eight times longer. This advantage becomes the critical factor, especially when performing three-dimensional simulations, because the number of grid points can be reduced by a factor of several tens to hundreds by adopting the RCIP scheme instead of the first-order upstream-difference scheme. The lack of numerical accuracy possibly results in seriously incorrect predictions of the change in the reservoir temperature. As shown in Figure 5, the location of the cooling front differs by a few hundred
8 meters in the profiles simulated using the RCIP and first-order upstream-difference schemes with the same spatial step length of 100 m. The temperature difference between the two profiles is up to around 20 deg. C. The time when the temperature begins to decrease i.e., when the cooling front passes also differs by around ten years, as shown in Figure 6. Such differences can be a critical factor in predicting the change in the productivity of hydrothermal reservoirs and in studying the feasibility of operating a geothermal power plant based on the numerical reservoir model. CONCLUSIONS From the numerical solutions presented in this paper, Figure 4: Changes in pressure and temperature simulated using the first-order upstream-difference a, CIP b, and RCIP c schemes. The numerical solution under the same condition obtained using the reservoir simulation HYDROTHERM is also shown for reference purposes. The time and spatial step lengths are years and 100 m, respectively.
9 Figure 5: Temperature profiles 15 years after the initial state obtained using the RCIP scheme RCIP with a spatial step length of 100 m and the first-order upstream-difference scheme USD with spatial step lengths of 12.5 to 100 m. we can conclude that: 1. the developed algorithm using the CIP and RCIP schemes can effectively control numerical diffusion; 2. the algorithm using the RCIP scheme yields more realistic solutions by controlling the overshoot that appears using the CIP scheme; 3. we can use spatial steps that are four to eight times longer without compromising numerical accuracy, by replacing the conventional firstorder upstream-difference scheme with the RCIP scheme. REFERENCES Faust, C. R., and Mercer, J. W Geothermal Reservoir Simulation 1 Mathematical Models for Liquid- and Vapor-Dominated Hydrothermal Systems, Water Resources Res., 15, IAPWS 1996 Release on the IAPWS Formulation 1995 for the Thermodynamic Properties of Ordinary Water Substance of General and Scientific Use, 18p. IAPWS 1997 Revised Release on the IAPS Formulation 1985 for the Viscosity of Ordinary Substance, 15p. Kipp, K. L., Jr., Hsieh, P. A., and Charlton, S. R Guide to the Revised Ground-Water Flow and Heat Transport Simulator: Figure 6: Changes in temperature 1,500 m away from the left boundary obtained using the RCIP scheme RCIP with a spatial step length of 100 m and the first-order upstream-difference scheme USD with spatial step lengths of 12.5 to 100 m. HYDROTHERM Version 3, USGS Techniques and Methods, 6-A25, 160p. Matsumoto, M Numerical Simulation of Two-phase Hydrothermal System Using a Simple Explicit Scheme and Parallel Computing Development and Application of the Electronic Steam Table GANSEKI, J. Geotherm. Res. Soc. Japan, 33, in Japanese with an English abstract. Takewaki, H., Nishiguchi, A., and Yabe, T The Cubic-Interpolated Pseudo-Particle CIP Method for Solving Hyperbolic-Type Equations, Journal of Computational Physics, 61, Xiao, F., Yabe, T., and Ito, T Constructing Oscillation Preventing Scheme for Advection Equation by Rational Function, Comput. Phys. Commn., 93, 1-12.
Study on the Numerical Accuracy for the CFD
Study on the Numerical Accuracy for the CFD T.Yamanashi 1, H.Uchida, and M.Morita 1 Department of Mathematics, Master s Research Course of Faculty of Science, Tokyo University of Science,1-3 Kagurazaka,
More informationNumerical Analysis of Shock Tube Problem by using TVD and ACM Schemes
Numerical Analysis of Shock Tube Problem by using TVD and Schemes Dr. Mukkarum Husain, Dr. M. Nauman Qureshi, Syed Zaid Hasany IST Karachi, Email: mrmukkarum@yahoo.com Abstract Computational Fluid Dynamics
More informationDevelopment of Numerical Methods for Estimating Fluid Flow Path in Fractured Geothermal Reservoir
PROCEEDINGS, 43rd Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 12-14, 2018 SGP-TR-213 Development of Numerical Methods for Estimating Fluid Flow Path
More informationPartial Differential Equations
Simulation in Computer Graphics Partial Differential Equations Matthias Teschner Computer Science Department University of Freiburg Motivation various dynamic effects and physical processes are described
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 informationMODELING OF FLOW AND TRANSPORT IN ENHANCED GEOTHERMAL SYSTEMS
PROCEEDINGS, Thirty-Sixth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 31 - February 2, 2011 SGP-TR-191 MODELING OF FLOW AND TRANSPORT IN ENHANCED GEOTHERMAL
More informationAnimation of Fluids. Animating Fluid is Hard
Animation of Fluids Animating Fluid is Hard Too complex to animate by hand Surface is changing very quickly Lots of small details In short, a nightmare! Need automatic simulations AdHoc Methods Some simple
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 informationAdarsh Krishnamurthy (cs184-bb) Bela Stepanova (cs184-bs)
OBJECTIVE FLUID SIMULATIONS Adarsh Krishnamurthy (cs184-bb) Bela Stepanova (cs184-bs) The basic objective of the project is the implementation of the paper Stable Fluids (Jos Stam, SIGGRAPH 99). The final
More informationAn explicit and conservative remapping strategy for semi-lagrangian advection
An explicit and conservative remapping strategy for semi-lagrangian advection Sebastian Reich Universität Potsdam, Potsdam, Germany January 17, 2007 Abstract A conservative semi-lagrangian advection scheme
More informationMicrowell Mixing with Surface Tension
Microwell Mixing with Surface Tension Nick Cox Supervised by Professor Bruce Finlayson University of Washington Department of Chemical Engineering June 6, 2007 Abstract For many applications in the pharmaceutical
More informationGALAXY ADVANCED ENGINEERING, INC. P.O. BOX 614 BURLINGAME, CALIFORNIA Tel: (650) Fax: (650)
GALAXY ADVANCED ENGINEERING, INC. P.O. BOX 614 BURLINGAME, CALIFORNIA 94011 Tel: (650) 740-3244 Fax: (650) 347-4234 E-mail: bahmanz@aol.com PUFF-TFT/PC A Material Response Computer Code for PC Computer
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 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 informationPorePy: A Python Simulation Tool for Fractured and Deformable Porous Media. Eirik Keilegavlen, Alessio Fumagalli, Runar Berge, Ivar Stefansson
PorePy: A Python Simulation Tool for Fractured and Deformable Porous Media Eirik Keilegavlen, Alessio Fumagalli, Runar Berge, Ivar Stefansson Summary Meshing and discretization of dynamics in fractured
More informationlecture 8 Groundwater Modelling -1
The Islamic University of Gaza Faculty of Engineering Civil Engineering Department Water Resources Msc. Groundwater Hydrology- ENGC 6301 lecture 8 Groundwater Modelling -1 Instructor: Dr. Yunes Mogheir
More informationC. A. D. Fraga Filho 1,2, D. F. Pezzin 1 & J. T. A. Chacaltana 1. Abstract
Advanced Computational Methods and Experiments in Heat Transfer XIII 15 A numerical study of heat diffusion using the Lagrangian particle SPH method and the Eulerian Finite-Volume method: analysis of convergence,
More informationA brief description of the particle finite element method (PFEM2). Extensions to free surface
A brief description of the particle finite element method (PFEM2). Extensions to free surface flows. Juan M. Gimenez, L.M. González, CIMEC Universidad Nacional del Litoral (UNL) Santa Fe, Argentina Universidad
More informationChapter 1 - Basic Equations
2.20 Marine Hydrodynamics, Fall 2017 Lecture 2 Copyright c 2017 MIT - Department of Mechanical Engineering, All rights reserved. 2.20 Marine Hydrodynamics Lecture 2 Chapter 1 - Basic Equations 1.1 Description
More informationNUMERICAL INVESTIGATION OF THE FLOW BEHAVIOR INTO THE INLET GUIDE VANE SYSTEM (IGV)
University of West Bohemia» Department of Power System Engineering NUMERICAL INVESTIGATION OF THE FLOW BEHAVIOR INTO THE INLET GUIDE VANE SYSTEM (IGV) Publication was supported by project: Budování excelentního
More informationTOUGHREACT Example: Aqueous Transport with Adsorption and Decay
403 Poyntz Avenue, Suite B Manhattan, KS 66502 USA +1.785.770.8511 www.thunderheadeng.com TOUGHREACT Example: Aqueous Transport with Adsorption and Decay PetraSim 2016.1 Table of Contents Overview...1
More informationOn the thickness of discontinuities computed by THINC and RK schemes
The 9th Computational Fluid Dynamics Symposium B7- On the thickness of discontinuities computed by THINC and RK schemes Taku Nonomura, ISAS, JAXA, Sagamihara, Kanagawa, Japan, E-mail:nonomura@flab.isas.jaxa.jp
More informationLiterature Report. Daniël Pols. 23 May 2018
Literature Report Daniël Pols 23 May 2018 Applications Two-phase flow model The evolution of the momentum field in a two phase flow problem is given by the Navier-Stokes equations: u t + u u = 1 ρ p +
More informationQUASI-3D SOLVER OF MEANDERING RIVER FLOWS BY CIP-SOROBAN SCHEME IN CYLINDRICAL COORDINATES WITH SUPPORT OF BOUNDARY FITTED COORDINATE METHOD
QUASI-3D SOLVER OF MEANDERING RIVER FLOWS BY CIP-SOROBAN SCHEME IN CYLINDRICAL COORDINATES WITH SUPPORT OF BOUNDARY FITTED COORDINATE METHOD Keisuke Yoshida, Tadaharu Ishikawa Dr. Eng., Tokyo Institute
More informationA mass-conservative version of the semi- Lagrangian semi-implicit HIRLAM using Lagrangian vertical coordinates
A mass-conservative version of the semi- Lagrangian semi-implicit HIRLAM using Lagrangian vertical coordinates Peter Hjort Lauritzen Atmospheric Modeling & Predictability Section National Center for Atmospheric
More informationPorous Reactor with Injection Needle
Porous Reactor with Injection Needle Introduction This model treats the flow field and species distribution in an experimental reactor for studies of heterogeneous catalysis. The model exemplifies the
More informationTOUGH2 Example: Initial State and Production Analyses in a Layered Model
403 Poyntz Avenue, Suite B Manhattan, KS 66502 USA +1.785.770.8511 www.thunderheadeng.com TOUGH2 Example: Initial State and Production Analyses in a Layered Model PetraSim 5 Table of Contents Acknowledgements...
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 informationFaculty of Mechanical and Manufacturing Engineering, University Tun Hussein Onn Malaysia (UTHM), Parit Raja, Batu Pahat, Johor, Malaysia
Applied Mechanics and Materials Vol. 393 (2013) pp 305-310 (2013) Trans Tech Publications, Switzerland doi:10.4028/www.scientific.net/amm.393.305 The Implementation of Cell-Centred Finite Volume Method
More informationCHAPTER 3 MODELING OF DEAERATOR AND SIMULATION OF FAULTS
27 CHAPTER 3 MODELING OF DEAERATOR AND SIMULATION OF FAULTS 3.1 INTRODUCTION Modeling plays an important role in the prediction and assessment of plant performance. There are two ways of getting the model
More informationImaging of flow in porous media - from optimal transport to prediction
Imaging of flow in porous media - from optimal transport to prediction Eldad Haber Dept of EOS and Mathematics, UBC October 15, 2013 With Rowan Lars Jenn Cocket Ruthotto Fohring Outline Prediction is very
More informationCorrected/Updated References
K. Kashiyama, H. Ito, M. Behr and T. Tezduyar, "Massively Parallel Finite Element Strategies for Large-Scale Computation of Shallow Water Flows and Contaminant Transport", Extended Abstracts of the Second
More informationGeneric framework for taking geological models as input for reservoir simulation
Generic framework for taking geological models as input for reservoir simulation Collaborators: SINTEF: Texas A&M: NTNU: Stanford Stein Krogstad, Knut-Andreas Lie, Vera L. Hauge Yalchin Efendiev and Akhil
More informationSEAWAT 4 Tutorial - Example Problem
SEAWAT 4 Tutorial - Example Problem Introduction About SEAWAT-4 SEAWAT-4 is the latest release of SEAWAT computer program for simulation of threedimensional, variable-density, transient ground-water flow
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 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 informationINTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 3, 2012
INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 3, 2012 Copyright 2010 All rights reserved Integrated Publishing services Research article ISSN 0976 4399 Efficiency and performances
More informationPhase-field simulation of two-phase micro-flows in a Hele-Shaw cell
Computational Methods in Multiphase Flow III 7 Phase-field simulation of two-phase micro-flows in a Hele-Shaw cell Y. Sun & C. Beckermann Department of Mechanical and Industrial Engineering, University
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 informationCFD-1. Introduction: What is CFD? T. J. Craft. Msc CFD-1. CFD: Computational Fluid Dynamics
School of Mechanical Aerospace and Civil Engineering CFD-1 T. J. Craft George Begg Building, C41 Msc CFD-1 Reading: J. Ferziger, M. Peric, Computational Methods for Fluid Dynamics H.K. Versteeg, W. Malalasekara,
More informationModule 1 Lecture Notes 2. Optimization Problem and Model Formulation
Optimization Methods: Introduction and Basic concepts 1 Module 1 Lecture Notes 2 Optimization Problem and Model Formulation Introduction In the previous lecture we studied the evolution of optimization
More informationMid-Year Report. Discontinuous Galerkin Euler Equation Solver. Friday, December 14, Andrey Andreyev. Advisor: Dr.
Mid-Year Report Discontinuous Galerkin Euler Equation Solver Friday, December 14, 2012 Andrey Andreyev Advisor: Dr. James Baeder Abstract: The focus of this effort is to produce a two dimensional inviscid,
More informationSolve Non-Linear Parabolic Partial Differential Equation by Spline Collocation Method
Solve Non-Linear Parabolic Partial Differential Equation by Spline Collocation Method P.B. Choksi 1 and A.K. Pathak 2 1 Research Scholar, Rai University,Ahemdabad. Email:pinalchoksey@gmail.com 2 H.O.D.
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 informationExperimental Validation of the Computation Method for Strongly Nonlinear Wave-Body Interactions
Experimental Validation of the Computation Method for Strongly Nonlinear Wave-Body Interactions by Changhong HU and Masashi KASHIWAGI Research Institute for Applied Mechanics, Kyushu University Kasuga
More informationRealistic Animation of Fluids
Realistic Animation of Fluids p. 1/2 Realistic Animation of Fluids Nick Foster and Dimitri Metaxas Realistic Animation of Fluids p. 2/2 Overview Problem Statement Previous Work Navier-Stokes Equations
More informationAn investigation of spatial and temporal weighting schemes for use in unstructured mesh control volume finite element methods
ANZIAM J. 44 (E) ppc759 C779, 23 C759 An investigation of spatial and temporal weighting schemes for use in unstructured mesh control volume finite element methods S. L. Truscott I. W. Turner (Received
More informationOptimizing Bio-Inspired Flow Channel Design on Bipolar Plates of PEM Fuel Cells
Excerpt from the Proceedings of the COMSOL Conference 2010 Boston Optimizing Bio-Inspired Flow Channel Design on Bipolar Plates of PEM Fuel Cells James A. Peitzmeier *1, Steven Kapturowski 2 and Xia Wang
More informationApplication of Simulation Method for Analysis of Heat Exchange in Electronic Cooling Applications
Application of Simulation Method for Analysis of Heat Exchange in Electronic Cooling Applications D.Ramya 1, Vivek Pandey 2 Department of Mechanical Engineering, Raghu Engineering College 1 Associate Professor,
More informationDevelopment of a Maxwell Equation Solver for Application to Two Fluid Plasma Models. C. Aberle, A. Hakim, and U. Shumlak
Development of a Maxwell Equation Solver for Application to Two Fluid Plasma Models C. Aberle, A. Hakim, and U. Shumlak Aerospace and Astronautics University of Washington, Seattle American Physical Society
More informationPARAMETRIC STUDY WITH GEOFRAC: A THREE-DIMENSIONAL STOCHASTIC FRACTURE FLOW MODEL. Alessandra Vecchiarelli, Rita Sousa, Herbert H.
PROCEEDINGS, Thirty-Eighth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 3, 23 SGP-TR98 PARAMETRIC STUDY WITH GEOFRAC: A THREE-DIMENSIONAL STOCHASTIC
More informationThe Spherical Harmonics Discrete Ordinate Method for Atmospheric Radiative Transfer
The Spherical Harmonics Discrete Ordinate Method for Atmospheric Radiative Transfer K. Franklin Evans Program in Atmospheric and Oceanic Sciences University of Colorado, Boulder Computational Methods in
More informationA Novel Approach to High Speed Collision
A Novel Approach to High Speed Collision Avril Slone University of Greenwich Motivation High Speed Impact Currently a very active research area. Generic projectile- target collision 11 th September 2001.
More informationTutorial 1. Introduction to Using FLUENT: Fluid Flow and Heat Transfer in a Mixing Elbow
Tutorial 1. Introduction to Using FLUENT: Fluid Flow and Heat Transfer in a Mixing Elbow Introduction This tutorial illustrates the setup and solution of the two-dimensional turbulent fluid flow and heat
More informationUsing a Single Rotating Reference Frame
Tutorial 9. Using a Single Rotating Reference Frame Introduction This tutorial considers the flow within a 2D, axisymmetric, co-rotating disk cavity system. Understanding the behavior of such flows is
More informationInvestigating The Stability of The Balance-force Continuum Surface Force Model of Surface Tension In Interfacial Flow
Investigating The Stability of The Balance-force Continuum Surface Force Model of Surface Tension In Interfacial Flow Vinh The Nguyen University of Massachusetts Dartmouth Computational Science Training
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 informationThe 3D DSC in Fluid Simulation
The 3D DSC in Fluid Simulation Marek K. Misztal Informatics and Mathematical Modelling, Technical University of Denmark mkm@imm.dtu.dk DSC 2011 Workshop Kgs. Lyngby, 26th August 2011 Governing Equations
More informationCS-184: Computer Graphics Lecture #21: Fluid Simulation II
CS-184: Computer Graphics Lecture #21: Fluid Simulation II Rahul Narain University of California, Berkeley Nov. 18 19, 2013 Grid-based fluid simulation Recap: Eulerian viewpoint Grid is fixed, fluid moves
More informationDownscaling saturations for modeling 4D seismic data
Downscaling saturations for modeling 4D seismic data Scarlet A. Castro and Jef Caers Stanford Center for Reservoir Forecasting May 2005 Abstract 4D seismic data is used to monitor the movement of fluids
More informationImplicit surface tension model for stimulation of interfacial flows
surface tension model for stimulation of interfacial flows Vinh The University of Massachusetts Dartmouth March 3rd, 2011 About my project Project Advisor Dr. Mehdi Raessi Department of Mechanical Engineering
More informationContents Foreword...v Acknowledgments...vii 1. Introduction to Simulation and History Matching...1
Foreword Mathematical simulation of reservoir behavior may be used to help understand reservoir processes and predict reservoir behavior in addition simulation can be used as a tool for reservoir description
More informationFluid Simulation. [Thürey 10] [Pfaff 10] [Chentanez 11]
Fluid Simulation [Thürey 10] [Pfaff 10] [Chentanez 11] 1 Computational Fluid Dynamics 3 Graphics Why don t we just take existing models from CFD for Computer Graphics applications? 4 Graphics Why don t
More informationContinuum-Microscopic Models
Scientific Computing and Numerical Analysis Seminar October 1, 2010 Outline Heterogeneous Multiscale Method Adaptive Mesh ad Algorithm Refinement Equation-Free Method Incorporates two scales (length, time
More informationAn explicit and conservative remapping strategy for semi-lagrangian advection
ATMOSPHERIC SCIENCE LETTERS Atmos. Sci. Let. 8: 58 63 (2007) Published online 22 May 2007 in Wiley InterScience (www.interscience.wiley.com).151 An explicit and conservative remapping strategy for semi-lagrangian
More information99 International Journal of Engineering, Science and Mathematics
Journal Homepage: Applications of cubic splines in the numerical solution of polynomials Najmuddin Ahmad 1 and Khan Farah Deeba 2 Department of Mathematics Integral University Lucknow Abstract: In this
More informationThe Development of a Navier-Stokes Flow Solver with Preconditioning Method on Unstructured Grids
Proceedings of the International MultiConference of Engineers and Computer Scientists 213 Vol II, IMECS 213, March 13-15, 213, Hong Kong The Development of a Navier-Stokes Flow Solver with Preconditioning
More informationUnderstanding and Using MINC
Understanding and Using MINC Background In TOUGH2, the MINC (Multiple Interacting Continua) approach is used to model flow in fractured media. It is a generalization of the classical double-porosity concept
More informationDevelopment of the Compliant Mooring Line Model for FLOW-3D
Flow Science Report 08-15 Development of the Compliant Mooring Line Model for FLOW-3D Gengsheng Wei Flow Science, Inc. October 2015 1. Introduction Mooring systems are common in offshore structures, ship
More informationNumerical Modeling Study for Fish Screen at River Intake Channel ; PH (505) ; FAX (505) ;
Numerical Modeling Study for Fish Screen at River Intake Channel Jungseok Ho 1, Leslie Hanna 2, Brent Mefford 3, and Julie Coonrod 4 1 Department of Civil Engineering, University of New Mexico, Albuquerque,
More informationEFFICIENT PRODUCTION OPTIMIZATION USING FLOW NETWORK MODELS. A Thesis PONGSATHORN LERLERTPAKDEE
EFFICIENT PRODUCTION OPTIMIZATION USING FLOW NETWORK MODELS A Thesis by PONGSATHORN LERLERTPAKDEE Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the requirements
More informationNEW DATA AND FILE REQUIREMENTS FOR TET -1
PROCEEDINGS, Twenty-Eighth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 27-29, 2003 SGP-TR-173 NEW DATA AND FILE REQUIREMENTS FOR TET -1 G. Michael Shook
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 informationCFD modelling of thickened tailings Final project report
26.11.2018 RESEM Remote sensing supporting surveillance and operation of mines CFD modelling of thickened tailings Final project report Lic.Sc.(Tech.) Reeta Tolonen and Docent Esa Muurinen University of
More informationLab 9: FLUENT: Transient Natural Convection Between Concentric Cylinders
Lab 9: FLUENT: Transient Natural Convection Between Concentric Cylinders Objective: The objective of this laboratory is to introduce how to use FLUENT to solve both transient and natural convection problems.
More informationThermal Reservoir Simulator for Microsoft Windows TM
Thermal Reservoir Simulator for Microsoft Windows TM Designed for use with the Microsoft Windows XP TM, and later, operating systems: Simulator Features Simulator Applications Radial and Cartesian Coordinates
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 informationCoping with the Ice Accumulation Problems on Power Transmission Lines
Coping with the Ice Accumulation Problems on Power Transmission Lines P.N. Shivakumar 1, J.F.Peters 2, R.Thulasiram 3, and S.H.Lui 1 1 Department of Mathematics 2 Department of Electrical & Computer Engineering
More informationDevelopment of a Computational Framework for Block-Based AMR Simulations
Procedia Computer Science Volume 29, 2014, Pages 2351 2359 ICCS 2014. 14th International Conference on Computational Science Development of a Computational Framework for Block-Based AMR Simulations Hideyuki
More informationSPH: Why and what for?
SPH: Why and what for? 4 th SPHERIC training day David Le Touzé, Fluid Mechanics Laboratory, Ecole Centrale de Nantes / CNRS SPH What for and why? How it works? Why not for everything? Duality of SPH SPH
More informationSolving Partial Differential Equations on Overlapping Grids
**FULL TITLE** ASP Conference Series, Vol. **VOLUME**, **YEAR OF PUBLICATION** **NAMES OF EDITORS** Solving Partial Differential Equations on Overlapping Grids William D. Henshaw Centre for Applied Scientific
More informationVerification and Validation in CFD and Heat Transfer: ANSYS Practice and the New ASME Standard
Verification and Validation in CFD and Heat Transfer: ANSYS Practice and the New ASME Standard Dimitri P. Tselepidakis & Lewis Collins ASME 2012 Verification and Validation Symposium May 3 rd, 2012 1 Outline
More informationThree dimensional meshless point generation technique for complex geometry
Three dimensional meshless point generation technique for complex geometry *Jae-Sang Rhee 1), Jinyoung Huh 2), Kyu Hong Kim 3), Suk Young Jung 4) 1),2) Department of Mechanical & Aerospace Engineering,
More informationJournal 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 informationNUMERICAL ANALYSIS OF WIND EFFECT ON HIGH-DENSITY BUILDING AERAS
NUMERICAL ANALYSIS OF WIND EFFECT ON HIGH-DENSITY BUILDING AERAS Bin ZHAO, Ying LI, Xianting LI and Qisen YAN Department of Thermal Engineering, Tsinghua University Beijing, 100084, P.R. China ABSTRACT
More informationTMVOC Buckley-Leverett Flow
403 Poyntz Avenue, Suite B Manhattan, KS 66502 USA +1.785.770.8511 www.thunderheadeng.com TMVOC Buckley-Leverett Flow PetraSim 2016.1 Table of Contents 1. Buckley-Leverett Flow...1 Description... 1 Create
More informationModeling Evaporating Liquid Spray
Tutorial 17. Modeling Evaporating Liquid Spray Introduction In this tutorial, the air-blast atomizer model in ANSYS FLUENT is used to predict the behavior of an evaporating methanol spray. Initially, the
More informationLagrangian and Eulerian Representations of Fluid Flow: Kinematics and the Equations of Motion
Lagrangian and Eulerian Representations of Fluid Flow: Kinematics and the Equations of Motion James F. Price Woods Hole Oceanographic Institution Woods Hole, MA, 02543 July 31, 2006 Summary: This essay
More informationThe WENO Method in the Context of Earlier Methods To approximate, in a physically correct way, [3] the solution to a conservation law of the form u t
An implicit WENO scheme for steady-state computation of scalar hyperbolic equations Sigal Gottlieb Mathematics Department University of Massachusetts at Dartmouth 85 Old Westport Road North Dartmouth,
More informationTutorial 2. Modeling Periodic Flow and Heat Transfer
Tutorial 2. Modeling Periodic Flow and Heat Transfer Introduction: Many industrial applications, such as steam generation in a boiler or air cooling in the coil of an air conditioner, can be modeled as
More informationCHAPTER 5 STUDY OF THERMAL COMFORT IN A ROOM WITH INSECT PROOF SCREEN
146 CHAPTER 5 STUDY OF THERMAL COMFORT IN A ROOM WITH INSECT PROOF SCREEN 5.1 INTRODUCTION In recent days, most of the buildings are equipped with insect proof screens to keep the insect not to enter inside
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 informationNumerical and theoretical analysis of shock waves interaction and reflection
Fluid Structure Interaction and Moving Boundary Problems IV 299 Numerical and theoretical analysis of shock waves interaction and reflection K. Alhussan Space Research Institute, King Abdulaziz City for
More informationAbstract. 1 Introduction
On the Use of Pore-Scale Computational Models for Two-Phase Porous-media Flows M.A. Celia, P.C Reeves, H.K. Dahle Environmental Engineering and Water Resources Program, Department of Civil Engineering
More informationNUMERICAL DIFFUSION AND DISSIPATION IN HYDROSTATIC MODELS OF INTERNAL WAVES
NUMERICAL DIFFUSION AND DISSIPATION IN HYDROSTATIC MODELS OF INTERNAL WAVES Ben R. Hodges 1 (A.M., ASCE) and Sarah Kelly Delavan 2 ABSTRACT Analysis of numerical diffusion and dissipation rates in a hydrostatic
More informationHigh-fidelity CFD modeling of particle-to-fluid heat transfer in packed bed reactors
High-fidelity CFD modeling of particle-to-fluid heat transfer in packed bed reactors Proceedings of European Congress of Chemical Engineering (ECCE-6) Copenhagen, 6-20 September 2007 High-fidelity CFD
More informationModeling Evaporating Liquid Spray
Tutorial 16. Modeling Evaporating Liquid Spray Introduction In this tutorial, FLUENT s air-blast atomizer model is used to predict the behavior of an evaporating methanol spray. Initially, the air flow
More informationThe International Association for the Properties of Water and Steam
IAPWS G13-15 The International Association for the Properties of Water and Steam Stockholm, Sweden July 015 Guideline on the Fast Calculation of Steam and Water Properties with the Spline-Based Table Look-Up
More informationALE Seamless Immersed Boundary Method with Overset Grid System for Multiple Moving Objects
Tenth International Conference on Computational Fluid Dynamics (ICCFD10), Barcelona,Spain, July 9-13, 2018 ICCFD10-047 ALE Seamless Immersed Boundary Method with Overset Grid System for Multiple Moving
More informationTutorial: Heat and Mass Transfer with the Mixture Model
Tutorial: Heat and Mass Transfer with the Mixture Model Purpose The purpose of this tutorial is to demonstrate the use of mixture model in FLUENT 6.0 to solve a mixture multiphase problem involving heat
More information