A New Method for the Rapid Calculation of Finely-Gridded Reservoir Simulation Pressures

Size: px
Start display at page:

Download "A New Method for the Rapid Calculation of Finely-Gridded Reservoir Simulation Pressures"

Transcription

1 Brigham Young University BYU ScholarsArchive All Theses and Dissertations A New Method for the Rapid Calculation of Finely-Gridded Reservoir Simulation Pressures Benjamin Arik Hardy Brigham Young University - Provo Follow this and additional works at: Part of the Chemical Engineering Commons BYU ScholarsArchive Citation Hardy, Benjamin Arik, "A New Method for the Rapid Calculation of Finely-Gridded Reservoir Simulation Pressures" (005). All Theses and Dissertations This Thesis is brought to you for free and open access by BYU ScholarsArchive. It has been accepted for inclusion in All Theses and Dissertations by an authorized administrator of BYU ScholarsArchive. For more information, please contact scholarsarchive@byu.edu, ellen_amatangelo@byu.edu.

2 A NEW METHOD FOR THE RAPID CALCULATION OF FINELY-GRIDDED RESERVOIR SIMULATION PRESSURES by Benjamin A. Hardy A thesis submitted to the faculty of Brigham Young University in partial fulfillment of the requirements for the degree of Master of Science Department of Chemical Engineering Brigham Young University December 005

3 Copyright 005 Benjamin A. Hardy All Rights Reserved

4 BRIGHAM YOUNG UNIVERSITY GRADUATE COMMITTEE APPROVAL of a thesis submitted by Benjamin A. Hardy This thesis has been read by each member of the following graduate committee and by a majority vote has been found satisfactory. Date Hugh B. Hales, Chair Date Larry L. Baxter Date Merrill W. Beckstead

5 BRIGHAM YOUNG UNIVERSITY As chair of the candidate s graduate committee, I have read the thesis of Benjamin A. Hardy in its final form and have found that (1) its format, citations, and bibliographical style are consistent and acceptable and fulfill university and department style requirements; () its illustrative materials including figures, tables, and charts are in place; and (3) the final manuscript is satisfactory to the graduate committee and is ready for submission to the university library. Date Hugh B. Hales Chair, Graduate Committee Accepted for the Department William G. Pitt Graduate Coordinator Accepted for the College Alan R. Parkinson Dean, Ira A. Fulton College of Engineering and Technology

6 ABSTRACT A NEW METHOD FOR THE RAPID CALCULATION OF FINELY-GRIDDED RESERVOIR SIMULATION PRESSURES Benjamin A. Hardy Department of Chemical Engineering Master of Science A new method for the determination of finely-gridded reservoir simulation pressures has been developed. It is estimated to be as much as hundreds to thousands of times faster than other methods for very large reservoir simulation grids. The method extends the work of Weber et al. 7 Weber demonstrated accuracies for the pressure solution normally requiring millions of cells using traditional finite-difference equations with only hundreds of cells. This was accomplished through the use of finite-difference equations that incorporate the physics of the flow. Although these coarse-grid solutions achieve accuracies normally requiring orders of magnitude more resolution, their coarse resolution does not resolve local pressure variations resulting from fine-grid permeability variations. Many oil reservoir simulation models require fine grids to adequately represent the reservoir properties. Weber s coarse grids are of little value. This study

7 takes advantage of the accurate coarse-grid solutions of Weber, by nesting them in the requisite fine grids to achieve much faster solutions of the large systems. Application of the nested-grid method involved calculating an accurate solution on a coarse grid, nesting the coarse-grid solution as fixed points into a finer grid and solving. Best results were obtained when an optimal number of coarse-grid pressure points were nested into the fine grid and when an optimal number of nested-grid systems were used. Speed Increase Factor for Weber's Results Grid Size One Optimum Nested-Grid Two Optimum Nested-Grids Three Optimum Nested Grids

8 ACKNOWLEDGEMENTS I would like to express appreciation to my advisor Dr. Hugh B. Hales for his untiring support and kindness. I would also like to thank the members of my graduate committee, Dr. Larry L. Baxter and Dr. Merrill W. Beckstead for their input and willingness to be on my committee. I am very grateful for the financial support that IRSRI has provided during my time as a graduate student. I consider myself fortunate to have had the privilege of completing my Master of Science degree at Brigham Young University among such outstanding students and faculty.

9 TABLE OF CONTENTS List of Tables...xi List of Figures xiii Nomenclature xvi Chapter 1 Introduction Technologies That Depend on Reservoir Simulation.. 1. Reservoir Descriptions Reservoir Simulation Approaches Review of Multiscale Methods Summary.. 8 Chapter Background Weber s Equations and Nested-Grid Description 9. Comparison to Multigrid Method Reservoir Description Solution of the Pressure Equation Summary Chapter 3 Testing the Performance of Various Solvers as a Function of Grid Size Size Direct Methods Iterative Methods MATLAB 7.0 Iterative Methods Summary... Chapter 4 Performance of Various Solvers Direct Methods: Results Stationary and Nonstationary Iterative Methods: Results Summary 30 Chapter 5 Determination of Best Solver viii

10 5.1 Summary 3 Chapter 6 Implementation of Solver in Nested-Grid Method Grid Geometries and General Setup Calculation of Coarse-Grid Pressure Calculation of Fine-Grid Pressure: LaPlace SOR Results Nested-Grid Method: GMRES Results Data Analysis Summary 41 Chapter 7 Interpolation Considering Interpolation Summary 46 Chapter 8 Weber s Coefficients Summary 56 Chapter 9 Study with Weber s Coefficients Applied Grid Geometries and General Setup Nested-Grid Results Regression and Dimensionless Time Analysis Summary 63 Chapter 10 Regression of All Results Regression of LaPlace and Weber SOR GMRES Regression Summary 70 Chapter 11 Dimensionless Time and Optimum Coarse-Grid Size Dimensionless Time LaPlace SOR: Optimization of the Nested-Grid Method GMRES: Optimization of the Nested-Grid Method Weber SOR: Optimization of the Nested-Grid Method Comparison of the Performance of the Different Methods Optimal Number of Fixed Points Summary Chapter 1 Error Anlaysis Error Analysis Results...8 ix

11 1. Summary 86 Chapter 13 Multiple Nested-Grid Analysis Results Practical Multiple Nested-Grid Study Summary 93 Chapter 14 Conclusion and Future Work Conclusions Future Work.100 Appendix Appendix A Appendix B Appendix C Appendix D Appendix E References..101 Tables Referred to in Body of Thesis..103 Original LaPlace and GMRES Regression and Dimensionless Time Analysis Original Weber Regression and Dimensionless Time Analysis..13 Optimal Over-Relaxation Factor Finder..19 Summary of Computer Programs x

12 LIST OF TABLES Table Page I. General Objectives of Study II. III. IV. Grid Sizes and Dimensions..18 Fine-Grid Sizes Used in Nested-Grid Study Course-Grid Sizes 34 V. Grid Dimensions and Size VI. VII. VIII. IX. Nested-Grid Points Take From a 5x5x10 Grid 84 Nested-Grid Points Taken From a 15x15x30 Grid..84 Nested-Grid Points Taken From a 5x5x10 Grid..84 Nested-Grid Points Taken From a 15x15x30 Grid..84 X. Size of Coarse and Fine Grids in Nested-Grid Method...88 XI. A.XII A.XIII A.XIV A.XV. A.XVI. A.XVII. Dimensionless Time Required to Solve Nested-Grid Method Nested-Grid Configurations LaPlace SOR: Fine and Nested-Grid Results GMRES: Fine and Nested-Grid Results Interpolation Study Interpolation Study: Comparison with No Interpolation Weber SOR: Fine and Nested-Grid Results with Weber s Coefficients..106 A.XVIII. LaPlace SOR: Dimensionless Time Results of Nested-Grid Method..107 xi

13 A.XIX A.XX. A.XXI A.XXII B.XXIII B.XXIV B.XXV B.XXVI GMRES: Dimensionless Time Results of Nested-Grid Method..107 Weber SOR: Dimensionless Time Results of Nested-Grid Method. 108 Absolute Difference from a 10-9 Solution Percent Difference from 10-9 Solution Laplace SOR: Fine and Nested-Grid Results Fine and Nested-Grid Results for GMRES Dimensionless Time Results of Nested-Grid Method Dimensionless Time Results of Nested-Grid Method...10 C.XXVII Weber SOR: Fine and Nested-Grid Results with Weber s Coefficients..13 C.XXVIII Performance of Nested-Grid Method Using Weber s Coefficients xii

14 LIST OF FIGURES Figure Page 1. Weber et al. Results Depiction of Smoothing Property Hypothetical Reservoir Being Simulated Determination of the Optimal Over-Relaxation Factor Direct Methods: Time required for Convergence as a Function of Grid Size Direct Methods: Memory (RAM) Requirements as a Function of Grid Size Iterative Methods: Time Required for Convergence as a Function of Grid Size Iterative Methods: Memory (RAM) Requirements as a Function of Grid Size Top and Side Plane Views of Fixed Points in a 5x5x10 Grid Nested-Grid Configurations: Percentage of Fixed Points LaPlace SOR: Time Required for Convergence of Various Nested-Grids GMRES: Time Required for Convergence as a Function of Grid Size Comparison of Time Required to Obtain Full-Fine Grid Pressure Solution With and Without Interpolation Application of Weber s Coefficients Coarse Grid Dimensions in Feet x5x10 Coarse-Grid Pressures Nested Into a 15x15x30 Grid xiii

15 17. Weber SOR: Time Required for Convergence of Various Nested-Grid Configurations LaPlace and Weber SOR Combined Data Set Regression LaPlace and Weber SORCombined ORF Regression GMRES Regression Time/Iteration as a Function of Total Grid Size Laplace SOR: Improvement Obtained by the Nested-Grid Method GMRES: Improvement Obtained by the Nested-Grid Method Weber SOR: Improvement Obtained by the Nested-Grid Method Ratio of Improvement Dimensionless Time Summary Optimal Number of Fixed Points Multiple Nested-Grid Analysis: 10 3, 10 5, and 10 6 Grid Sizes Practical Nested-Grid Configuration A Look at Optimal, Practical and No Nested Grid Configurations B.31. LaPlace SOR Time/Iteration as a Function of Total Grid B.3. Improvement Obtained by the Nested-Grid Method B.33. Optimal Number of Coarse-Grid Points Nested Into Fine Grid B.34. GMRES: Time/Iteration as a Function of Total Grid Size B.35. Optimized Nested-Grid Method Using SOR Compared with GMRES B.36. Comparison of Improvement Obtained by the Nested-Grid Method for GMRES and SOR....1 C.37. Weber Time/Iteration as a Function of Total Grid Size C.38. Improvement Obtained by Nested-Grid Method xiv

16 C.39. Comparison of Dimensionless Time for all Solution Algorithms D.40. ORF D.41. ORF D.4. ORF D.43. ORF D.44. ORF D.45. ORF D.46. ORF D.47 ORF D.48. ORF D.49. ORF D.50. ORF D.51. ORF D.5. ORF xv

17 NOMENCLATURE a c i,j,k n p q r x,y,z Proportionality constant Integration constant Grid point location in x,y,z respectively Number of grid refinements; number of wells Pressure Flux of fluid through porous media Distance to the wells Cartesian coordinates A, B, C, D Variables determined by regression Fx, Fy, Fz Fraction of x, y, and z-distances that the fine-grid points are from point (1) of the course-grid block K Permeability K Pseudo-permeability N Parameter used to describe effectiveness of well, determined by regression P Pressure N ITER N FG N CG N LAPLACE N WEBER N CG` N FG` Q n r n t d t total t coarse-grid t nested-grid Number of Iterations Number of fine-grid points Number of coarse-grid points Parameter used to describe effectiveness of well for LaPlace SOR Parameter used to describe effectiveness of well for Weber SOR N CG N N FG N Volumetric injection rate of well n Distance to well n Dimensionless time Total dimensionless time Dimensionless time required to calculate coarse-grid solution Dimensionless time required to calculate a given nested-grid solution λ Ratio of pseudo-permeability to actual permeability µ Fluid viscosity ω Over-relaxation factor ω opt Optimal Over-relaxation factor Ω Solid angle Positive cell face - Negative cell face xvi

18 CHAPTER 1 INTRODUCTION This thesis proposes a new method of linear algebra that provides a rapid solution of very large systems of equations. This work was motivated by the need for such a technology in reservoir simulation. However, the method is likely to be valuable in other disciplines where large systems of linear equations are encountered, such as in computational fluid mechanics. Oil is the lifeblood of America s economy. It presently supplies more than 40% of the nation s total energy demands. 1 Recently the Department of Energy (DOE) has expressed two key concerns over oil in America: (1) maintaining an immediate readiness to respond to oil supply disruptions and () keeping America s oil fields producing in the future. The DOE has determined that one way to prevent an oil supply disruption is to make certain that domestic production of oil is maintained. Remaining U.S. oil fields are becoming progressively more costly to produce because much of the easily produced oil has already been recovered. Better technology is needed to locate and produce the residual oil. 1 Many of the emerging technologies that will keep U.S. oil fields producing long into the future do and will employ reservoir simulators. 1

19 1.1 Technologies That Depend on Reservoir Simulation Reservoir simulation became feasible some thirty years ago with the dawn of the computer age. Since that time, computer technology has increased greatly. Computer hardware and software have become continuously more powerful at an astonishing rate. However, throughout the history of computers reservoir simulators have always taxed the very fastest machines, and despite these enormous advances in computer technology, reservoir simulators continue to strain the largest and fastest computer systems. The need for faster reservoir simulation remains as critical today as it has ever been. Emerging technologies in the oil industry that depend on reservoir simulators would be greatly enhanced by more accurate and faster reservoir simulations. Brief descriptions of some of the emerging technologies key to optimal oil recovery are described as follows. Geostatistics: Geostatistics involves gathering large quantities of data from exposed rock out-crops, where measurements of physical properties can be made with relative ease. The statistical variations in physical properties observed in the out-crops are applied to data obtained from well cores, logs, and from surface seismic data for subterranean reservoirs. Many possible subsurface models of the same field, each with some probability of being correct, can then be generated. Simulation of all the geostatistical reservoir models is practical only when simulations can be made rapidly. Results of these reservoir simulations are valuable as they can provide a measure of a field s potential profitability as well as the associated risks. Automated History Matching: History matching plays a critical role in monitoring displacement processes, constructing good reservoir models, and predicting future reservoir performance. Production data are the most common type of reservoir data. Matching these data allows reservoir engineers to better characterize reservoir properties such as permeability and porosity. This process has been automated, and is relatively simple when the number of data values to be adjusted is small and when simulations can be accomplished rapidly. For more extensive applications, history matching remains a difficult and time-consuming task, as many simulations of the same field must be made to determine a data set that matches the production history of the field. Optimization: Optimization of reservoir production involves repeated simulations of a particular field to determine optimal reservoir design variables

20 and design functions such as well location, geometry, completion intervals, well rates, well remedial treatments, etc. Optimization software is available, yet fast, accurate and robust reservoir models are essential to enhance the optimization process. Smart Wells: Smart wells are equipped with down-hole measurement equipment, packers, and control valves. They use real-time simulations to control the production rate from different well segments through the use of packers, which isolate the various production intervals in a well, and valves, that control the amount of flow from each interval. Smart wells have the potential to significantly increase production (up to 65%). 3 Optimal real-time control theory, used to determine proper packer and valve settings, requires rapid reservoir simulation. 1. Reservoir Descriptions Modern reservoir imaging techniques and advances in geological modeling are allowing very detailed reservoir descriptions. Often ten million cells or more may be used to describe reservoir rock properties. With current computers, most oil companies cannot afford to run routine fluid-flow simulations with more than about 100,000 grid blocks. 4 Typical reservoir models generally contain 10, ,000 cells. 5 The alternative to running multi-million cell simulation models that incorporate the geologic models, is Upscaling. Upscaling is the practice by which a fine geologic model containing a detailed description of reservoir rock properties is replaced by a coarser scale description of equivalent reservoir rock properties which is more suitable for reservoir simulation. The coarsening of the geologic model can lead to reservoir simulation runs that can be completed in hours, as opposed to days. The challenge in applying this process is to retain an accurate representation of the physical characteristics of the geologic model critical to fluid flow. 6 A variety of different methods has been proposed to upscale single and multiphase flow properties. Many studies 7,8,9,10,11,1 and reviews 4,5,13,14,15,16 have been conducted, yet all upscaling methods 3

21 suffer from the problem that either, they make assumptions about the large-scale boundary conditions, which may significantly affect the results, or they require a finegrid solution to derive coarse grid properties. 17 Other unsettled issues include grouping of upscaled relative permeabilities, robustness and process independence. There are still unresolved issues dealing with extending multi-phase scales up to three-phase flow, compositional flow, and flow in naturally occurring systems. 5 Consequently, despite a huge literature on multiphase upscaling, the best approach is still very much an open issue. In fact, the current industry practice is to limit upscaling to single-phase properties only. 17 In a recent review, M. A. Christie stated that, The most promising methods of the last few years may be those that regard upscaling as an integral part of the solution of the flow equations, rather than as an external process which have the correct boundary conditions to provide the correct answer Reservoir Simulation Approaches In the May 17 th, 004, Oil & Gas Journal, 6 Scott Evans summarized four basic approaches that can be used for reservoir simulation in today s high performance computing environments. Four basic approaches that can be successfully used in reservoir simulation are available for subsurface modeling: 1. Traditional. A geologic model upscaled to a second model and used for reservoir simulation. The result is two models, only one of which is used for simulation.. Geologic. A geologic model run without coarsening in the reservoir simulator. The result is one model but with potentially very slow simulation runs 4

22 3. Hybrid. One model generated, but with varying scale, in which the model has more detail where needed and less where it is not. The result is one model with potentially more acceptable simulation runs. 4. Multiscale. Two or more models, one fine and one coarse, linked and used simultaneously in the reservoir simulator. This is an area of industry research and would have each of the models used in the simulation as appropriate. The multiscale approach is a current area of research and development. This methodology will hopefully provide a way to more accurately represent the static and dynamic properties of a reservoir with minimal computational power. Recent studies have shown how to use both coarse and fine grid information during reservoir simulations. 1.4 Review of Multiscale Methods In 1991, Ramé and Killough 18 presented the first implementation of a multiscale simulation technique. The method used the implicit-pressure, explicit saturation (IMPES) procedure to decouple the pressure equation from the conservation equations numerically. A fourth-order finite element method was used to solve the pressure equation on each coarse grid. Fine-scale information was interpolated from the coarse grid using a splines-under-tension technique, and the conservation equations for fluid transport were solved on the fine grid. Time stepping was performed on the fine grid and after several timesteps the current mobilities on the fine grid were passed to the coarse grid to update the pressure field. D examples for miscible flow were presented. In 1995 Guérillot and Verdière 19 proposed a dual mesh method where the pressure field was first computed on a coarse grid, but where the saturation was moved on the fine grid. Using approximate boundary conditions, the velocity field was estimated within 5

23 each coarse block by solving for the pressure. This technique was applied to a D singlephase model and in 1996 they extended their work to multiphase flow. 0 They showed the results in D for two-phase flow models with simplified injection-production boundary conditions. Comparison of the time of calculation spent for the full fine-grid calculation and the dual mesh method gave a speed-up factor from 5 to 7. In 1996 Hou and Wu 1 presented the derivation of a mathematically rigorous multiscale finite-element method for solving the class of elliptical problems that arise from composite material (steady-heat conduction through a composite material with tubular fiber reinforcement in a matrix) and flows in porous media. By constructing multiscale finite-element base functions that were adaptive to the local property of the differential operator, this method was able to capture efficiently the large-scale behavior of the solution without resolving the small-scale features. Results for flow in porous media were provided in D without gravity and capillary effects. In 1999 Guedes and Schoiozer implemented the same methodology as Guérillot and Verdière, using an upscaling method developed by Hermitte and Guerillot. 3 They provided results in D and included gravity effects. Well boundary conditions were considered as sources or sinks applied in one grid block. Also in 1999, Gautier et al. 4 presented a similar approach using streamline-based simulation (an IMPES method). The pressure solve method (psm) was used to upscale the transmissivities for each coarse-grid block from petrophysical properties defined on a fine grid. Gravity effects and wells where included. The well pressures were determined using Peaceman s well model, and special attention was given to keep equivalent fluxes on both coarse and fine scales. The method was tested on a series of waterflood problems 6

24 and it was demonstrated that the method could give accurate estimates of oil production for large 3D models up to 8.5 times faster than direct simulation using streamlines on the fine grid. The results were very efficient in terms of CPU time and memory management. Arbogast and Bryant 5 introduced a slightly different dual grid approach in 001 by using Green function methods to upscale transmissivities and a mixed finite-element method to solve the pressure field. Gravity and capillary effects were considered in their methodology. Similar to Guatier et al. 4, they observed a reduction ratio for the time of calculation of about to 10 compared to fine-grid simulation. Audigane and Blunt (004) 17 present an extension of the dual-mesh method of Guérillot and Verdière to 3D cases and included gravity and wells. Using an IMPES method, the pressure field was solved on the coarse mesh with a conventional finitedifference scheme. Transmissivities were upscaled either with the psm method or with a simple geometric average (ga). The pressure field was reconstructed on the fine mesh using flux boundary conditions from the coarse-grid simulation. Multiscale methods continue to develop and make advances in reducing the cpu time of reservoir simulations. The solution of the pressure equation is the most timeconsuming step in any reservoir simulation 5, and it has become more evident that improved mathematics, which will allow faster and more accurate simulations, is fundamental to improved reservoir simulation capabilities. In reservoir simulation, the primary concern is movement of gas, oil, and water in the reservoir. 6 These fluids flow as a result of pressure variations in the reservoir. Hence, the accurate prediction of reservoir pressures is essential to a good reservoir simulator. This work proposes a new linear algebra technique for the solving of finely-gridded 7

25 reservoir pressures. The new method is multiscale, in that it involves the calculation of the reservoir pressures on a coarse grid and then uses the coarse-grid solution in a nested grid to calculate fine-grid pressures. Improved mathematics are key to the method and permit faster and more accurate solutions of the pressure equation. 1.5 Summary Better technology is needed to produce oil and gas reserves in a more effective manner. Many developing technologies that improve oil and gas production would be greatly enhanced by faster and more accurate reservoir simulators. In approaching reservoir simulation, the multiscale method is an area of industry research which has been shown by various researchers to reduce the computational time required to complete reservoir simulations. Calculation of reservoir pressures is the most time consuming step in any reservoir simulator; this thesis proposes a new linear algebraic method that incorporates the multiscale approach to calculate reservoir simulation pressures in an accurate and fast manner. 8

26 CHAPTER BACKGROUND.1 Weber s Equations and Nested-Grid Description The new pressure-solution method developed in this thesis is based on the work of Weber et al. 7 They proposed that finite-difference equations, used to represent the pressure equation, be based on mathematical expressions that incorporate the physics of the process instead of on traditional polynomial expressions. In modeling reservoir pressures, equations incorporating the physically realistic ln(r) dependence on pressure for reservoirs with straight line wells, and a 1/r dependence for reservoirs with more complex well geometries were used (r is the distance to the wells). Weber investigated formulations in which ln(r) s and 1/r s were summed over all the wells in the reservoir, and in which only the single, closest well value was used. The results of Weber s study that incorporated a 1/r dependence into the finite-difference equations are shown in Figure 1. The figure compares the accuracy of the pressures calculated by various methods for an 11x11x grid of a rectangular reservoir of similar geometry to that used in this study. The different methods were compared with an analytical solution generated by Weber. The new finite-difference equations showed a four-order-of-magnitude improvement in accuracy compared with traditional polynomial based finite-difference equations and approximately three-orders-of-magnitude improvement over Peaceman s 9

27 Correction. This dramatic improvement in accuracy was the catalyst for the development of a new method to calculate finely-gridded pressures Figure 4: Reservoir Pressure Error Summary Error (psi) Average Maximum 0.01 Traditional Finite Difference Method Traditional Method w/ Peaceman Correction Inverse-r Finite Difference Method Inverse-r Finite Volume Method Single Well Method Figure 1. Weber et al. Results 7 The new method involves two steps: (1) The creation of a course-grid solution using Weber s finite-difference equations that incorporate the physics of the flow to obtain an accurate pressure solution on a small, coarse grid, and () nesting the coarsegrid solutions into a desired fine grid and solving the system of equations to obtain detailed pressures that honor the nested-course-grid pressures. The improved mathematics of Weber et al. is vital to this method, as it permits the generation of accurate solutions on coarse grids. 10

28 . Comparison to Multigrid Method This method is not a multigrid method, yet it does utilize some of the basic ideas of the multigrid method, namely coarse-grid relaxation and nested iteration. Multigrid methods are built on the fact that many standard iterative methods possess a smoothing property. The smoothing property describes the fact that as iterative methods progress to convergence the reduction in error decreases and becomes smooth. Figure is taken from the book A Multigrid Tutorial by William L. Briggs 8 and shows this smoothing effect for the weighted Jacobi method in a plot of the absolute error versus iteration number. The error decreases quickly within the first five iterations and then decreases slowly. Absolute Error Iteration Number Figure. Depiction of Smoothing Property 8 The initial rapid decrease in error is associated with the quick elimination of highfrequency modes, and the slow decrease is due to the presence of low-frequency modes. The modes are Fourier modes where a small wave number, k, is associated with long 11

29 smooth waves, and large values of k correspond to highly oscillatory waves. As stated earlier, many standard iterative methods possess this smoothing property which makes them very effective at eliminating high-frequency or oscillatory components of the error, while leaving the low-frequency or smooth components relatively unchanged. Through the use of coarse grids, the multigrid method puts the smoothing property to good use as the smooth error is relatively more oscillatory on coarser grids. Hence relaxation becomes more effective on coarser grids. Coarse grids can be used to compute an improved initial guess for fine-grid relaxation, and a reliable way to improve the relaxation scheme on the fine grid is to use a good initial guess. This well-known technique of using a coarse grid to generate improved initial guesses on a fine grid is called nested iteration. In a nut shell, multigrid methods relax on a fine grid until the smooth error is reached and the high frequency oscillations have been eliminated. The grid is then restricted to a coarse grid where the error is relatively more oscillatory. The grid is again relaxed until smooth error is reached. The coarse grid points become a new and improved initial guess for the fine grid where the coarse grid points are mapped to the fine grid through interpolation. The fine grid is relaxed again and the process continued to convergence. 8 The research of this thesis investigates the potential value of implementing, in a nested-grid/multiscale manner, the new finite-difference equations developed by Weber et al. However, unlike the multigrid method described previously, the new method can potentially be completed in only two steps, one using a course grid and a second using a fine grid. Only two grids are necessary because of the very high accuracy of the coarsegrid solution using Weber s solution. 1

30 This thesis describes development of this new method. The steps involved are outlined in Table I. Table I. General Objectives of Study Objectives of Study Test performance of various solvers as a function of grid size (Chapters 3,4) Analyze results to determine the best solver(s) for coarse and fine grids. Selected solver(s) will be used in nested-grid solution method (Chapter 5) Implement best solver(s) in the nested-grid solution method (Chapter 6) Consider potential value of applying interpolated values to initial solution guesses (Chapter 7) Apply Weber s coefficients in nested-grid solution method (Chapters 8,9) Regress all results to determine potential correlations that fit data (Chapter 10) Conduct a dimensionless time analysis using correlations determined from regressions and optimize the nested-grid configuration (Chapter 11) Analyze the error of the nested-grid solution method (Chapter 1) Investigate the use of multiple nested-grids to improve method further (Chapter 13).3 Reservoir Description In this work, the reservoir geometry considered is shown in Figure 3. Injection at 1,500 psi occurs in one well; production at -1,500 psi occurs in the other well. The two wells are centered, one in each of the two cubic elements comprising the threedimensional rectangular reservoir. The boundary conditions of the reservoir are no flow, i.e. the pressure gradients in the direction normal to the boundaries are zero. 13

31 Injection Well Production Well Figure 3. Hypothetical Reservoir Being Simulated 7 The dimensions of the reservoir are in the following ratio: 1x1x. The actual reservoir dimensions were factored in by scaling the well radius with the size of the grid block. For grid blocks of 100 foot dimensions, the well radius of three inches becomes Although gravity effects are assumed negligible, the reservoir was considered to lie with its largest dimension in the horizontal plane..4 Solution of the Pressure Equation The pressure equation is written in terms of average pressure for the conservation of mass flowing through porous material. The incompressible, steady-state, threedimensional pressure in a reservoir for which the mobility is everywhere uniform, is given by the Laplace equation: P P P x y z = 0 (-1) 14

32 15 Replacing each of the second derivatives by second-order, centered-difference approximations at grid point (i,j,k) yields 0 1,,,, 1,, 1,,,, 1,,, 1,,,, 1, = z P P P y P P P x P P P k j i k j i k j i k j i k j i k j i k j i k j i k j i (-) In the special case where x = y = z the grid aspect ratio is unity and Equation (-) becomes 6 1,, 1,, 1,, 1,,, 1,, 1,,, = k j i k j i k j i k j i k j i k j i k j i P P P P P P P (-3) This is the traditional, finite-difference formulation used to solve the pressure equation in three dimensions..5 Summary Weber et. al 7 developed finite-difference equations that incorporate physics of flow and allow for pressure solutions, four-orders-of-magnitude more accurate, to be generated on coarse grids. The proposed new calculation method developed in this thesis takes advantage of Weber s accurate coarse-grid solutions in a nested-grid calculation method to determine finely-gridded reservoir simulation pressures. Although this method incorporates some of the basic ideas of multigrid methods, namely coarse-grid relaxation and nested iteration, it is not a multigrid method as the coarse-grid pressure solutions are nested and fixed throughout the entire relaxation process of the fine grid. The reservoir

33 being simulated is of an idealized nature. The simplify assumptions include: (1) homogeneous permeability, () neglected gravity effects, (3) incompressibility (4) steady-state, and (5) dimensions of reservoir and placement of wells generate a symmetric pressure solution. 16

34 CHAPTER 3 TESTING THE PERFORMANCE OF VARIOUS SOLVERS AS A FUNCTION OF GRID SIZE Reservoir simulations can be made at varying degrees of grid refinement ranging, as mentioned previously, from hundreds to millions of grid blocks. For this study, it was important to know how various solvers performed as a function of grid size so that the best solver(s) could be used for all grid sizes in the nested-grid method. The study was initialized by considering the performance of various numerical methods for the solution of systems of linear algebraic equations as a function of grid size on a standard desktop computer. The computer used four Intel (R) Pentium (R) processors running at 1.80GHz, and 655 MB of RAM. In the entire study, swapping RAM data to the hard drive paging was not apparent. Hence the same performance increase would be expected on any computer (workstation, supercomputer) with sufficient memory to avoid paging. Both direct and iterative methods were considered and compared. The performance metrics included convergence time, iterations required to converge, and run-time memory requirements. To study the performance of the various linear algebraic solvers as a function of grid size, a range of test grid sizes was selected. The grids were constructed so that the wells were always in the center of their respective cells. Table II summarizes the various grid sizes used. 17

35 Table II. Grid Sizes and Dimensions Grid Size Grid Dimensions 50 5X5X X7X X9X X11X X17X X1X X33X X37X X51X X65X Direct Methods Direct solution methods perform very well for small grids, but can require excessive computational effort and computer memory as the number of grid points increases. 9 Direct methods might be preferable for the course-grid solution, while iterative methods might be required for the fine grid. The direct elimination methods analyzed were Gauss Elimination 9, the Doolittle LU factorization method 9, and a band solver, DGBSV, from the LAPACK library. 30 The programs were developed using Compaq Visual Fortran 6.5 (Fortran 90). Default compiler options were used throughout so that others could duplicate results. 3. Iterative Methods Iterative methods can be divided into two general categories, stationary and nonstationary. 31 Three standard stationary iterative methods were considered, namely Jacobi, Gauss-Seidel, and Successive-Over-Relaxation (SOR). Five advanced nonstationary iterative methods from MATLAB and a hybrid of direct and iterative 18

36 methods, Line Successive-Over-Relaxation (LSOR), were also considered. Iterative methods should be better suited for the large systems of equations required in this study. The computer programs for these solvers, other than the five found in MATLAB, were developed in-house with the incorporation of the Thomas algorithm (DGTSV) from the LAPACK library 30 for LSOR. The in-house programs were developed in Compaq Visual Fortran 6.5 (Fortran 90). Some iterative methods require diagonal dominance to guarantee convergence and in general are less robust than direct solvers. The system of equations arising from the seven-point, second order approximation of the Laplace equation used in the simulations of this study is always diagonally dominant. Hence, such iterative solvers work well. For any iterative method, an initial approximation must be made for the solution to start the process. Several choices are available: (1) Simply let the solution (in this case the value of the pressure P i,j,k at the various grid points) equal zero at all non-specified points; () Approximate P i,j,k by some average of the well pressures and or boundary conditions; or (3) Construct a solution on a course grid, and map the course-grid values onto the fine grid through interpolation. 9 For the determination of the best solver to be used in the nested-grid method, the initial P i,j,k solution vector was set to 0.0 everywhere except at the wells, which was set to the average of the well pressures, 1500 and However, later in the study (See Chapters 6 and 9), coarse-grid points would be nested into the fine grid, and three-dimensional linear interpolation would also be considered to establish good starting values for iterative methods on the fine grid. 19

37 Iterative methods only generate an approximate solution after a finite number of iterative steps; the iterative process terminates when it meets a specified convergence criterion. In general, the number of iterations required to satisfy the convergence criterion is influenced by diagonal dominance, method of iteration, initial solution vector, and the convergence criterion itself. 9 The convergence criterion used by the in-house iterative methods in this study is described by the following equation. P 1,1,1 P Imax,Jmax,Kmax 10-6 (3-1) P 1,1,1 and P Imax,Jmax,Kmax are the values of the pressures at opposite corner points of the grid furthest from each other in three dimensions. The symmetry of the problem indicates that P 1,1,1 equals P Imax,Jmax,Kmax in a properly converged solution. P 1,1,1 is the first to be calculated and hence is the least converged. P Imax,Jma,Kmax is the last to be calculated and would be expected to be the most converged. Since the pressures are initially at zero and relax to their solution values asymptotically, with points near the wells changing most rapidly, this convergence criteria should represent the maximum error in the solution after many iterations. For the iterative method of SOR and the hybrid method of LSOR, another key factor in the determination of the number of iterations required for convergence was the value of the over-relaxation factor, ω. When ω equals one, SOR yields the Gauss-Seidel method. When ω is greater than one, but less than two, the system is over-relaxed; when the ω factor is equal two or greater than two, the system becomes unstable. Figure 4 shows a plot of the iterations required for convergence as a function of the over- 0

38 relaxation factor for a 17x17x34 grid solved by SOR. In this case it is apparent that by using the optimal over-relaxation factor one reduces the number of iterations required for convergence by approximately two orders of magnitude in comparison with the Gauss- Seidel iteration method. Optimal Over-Relaxation Factor: 17x17x34 1.E05 1.E04 Iterations 1.E03 1.E Over-Relaxation Factor Figure 4. Determination of the Optimal Over-Relaxation Factor The relaxation factor does not change the final solution since it multiplies the residual, which is zero when the final solution is reached. The major difficulty with the overrelaxation method is the determination of the best value for ω. In general, the optimal value of the over-relaxation factor, ω opt, depends on the size of the system of equations and the nature of the equations. As a general rule, larger values of ω opt are associated with larger systems of equations. 9 The optimum value of ω was determined by 1

39 experimentation for the various grid sizes considered in this study. The process of finding ω opt was straight forward but often time consuming, especially for large grids. 3.3 MATLAB 7.0 Iterative Methods MATLAB is a high-performance language for technical computing; the name stands for matrix laboratory. MATLAB incorporates LAPACK and BLAS libraries in its software for matrix computation. Nine functions are available in MATLAB that implement advanced iterative methods for sparse systems of simultaneous linear systems. Of these nine, five where considered: Biconjugate Gradient (BICG), Biconjugate Gradient stabilized (BICGSTAB), LSQR implementation of Conjugate Gradients on the Normal Equations (LSQR), Generalized Minimum Residual (GMRES), and Quasiminimal Residual (QMR). GMRES is a Krylove Subspace Method and was recommended for use in this study by simulation developers from ConocoPhillips. All of the MATAB algorithms were implemented without preconditioners, and the magnitude of the convergence tolerance was Summary A standard desktop computer was used for the study of the performance of various solvers as a function of gird size. Direct solution methods considered were Gauss Elimination 9, the Doolittle LU factorization method 9, and a band solver, DGBSV, from the LAPACK library. 30 Stationary iterative methods considered for the study were Jacobi, Gauss-Seidel, and Successive-Over-Relaxation (SOR). Five nonstationary iterative methods from MATLAB 7.0: Biconjugate Gradient (BICG), Biconjugate Gradient

40 stabilized (BICGSTAB), LSQR implementation of Conjugate Gradients on the Normal Equations (LSQR), Generalized Minimum Residual (GMRES), and Quasiminimal Residual (QMR) were considered. A hybrid of direct and iterative methods, Line Successive-Over-Relaxation (LSOR), was also considered. SOR and LSOR require the determination of an optimal over-relaxation factor. 3

41 4

42 CHAPTER 4 PERFORMANCE OF VARIOUS SOLVERS 4.1 Direct Methods: Results Figure 5 shows a plot of convergence time as a function of grid size for the direct methods of Gauss Elimination (GE), Doolittle LU factorization (LU) and the band solver, DGBSV, from the LAPACK library. Due to the computational limits of the computer, the solvers performance was only considered on smaller grid sizes. Direct Methods: Time vs Grid Size 1.E04 1.E03 1.E0 y = 8E-11x R = CPU Time (s) 1.E01 1.E00 1.E-01 y = 8E-09x.4374 R = GE LU DGBSV 1.E-0 1.E-03 1.E0 1.E03 1.E04 1.E05 Grid Size Figure 5. Direct Methods: Time Required for Convergence as a Function of Grid Size 5

43 A power law fit of the relationship between time and grid size is included for the DGBSV solver. GE and LU data points are nearly coincident and the same power law fit is used to describe both. Indeed, these direct methods proved to work reasonably well on small grids, yet bigger grid sizes required large amounts of memory and numerous calculation steps. The direct methods eventually became impractical for the larger grids. For example, an 11x11x grid results in =,6 equations in,6 unknowns and a full coefficient matrix with (,6) = 7,086,44 array elements. For the Gauss Elimination and Doolittle LU factorization programs, the next largest grid size considered, 17x17x34, has an input array of 96,550,76 elements and would have taken an estimated four and a half days to compute using 378,544 K of RAM. The banded solver from the LAPACK library, DGBSV, showed some improvement in both time and memory performance but still struggled on larger grids. The memory requirements for the direct methods are shown in Figure 6 with power law fits of the data for the DGBSV and LU/GE solver. These general trends were expected, yet the study was conducted for comparative purposes and to aid in understanding the solution process. Direct solutions could be preferable for obtaining the course-grid solutions for the new method. Again, since GE and LU memory requirements are so similar, the plotted points are essentially coincident and one trend line describes both. The left-most points are influenced by the finite amount of RAM that would be required even for a zero grid size, hence these points were excluded in the regression. 6

44 Direct Methods: RAM Requirements vs Grid Size 1.E06 1.E05 y = x R = RAM(K) 1.E04 y = x R = GE LU DGBSV 1.E03 1.E0 1.E0 1.E03 1.E04 1.E05 Grid Size Figure 6. Direct Methods: Memory (RAM) Requirements as a Function of Grid Size 4. Stationary and Nonstationary Iterative Methods: Results The iterative methods performance was tested at various grid sizes. Of all the iterative methods described in the previous chapter, LSOR was found to have the best convergence-time performance with SOR following closely. Of the five nonstationary MATLAB iterative methods, GMRES and LSQR were shown to have the best convergence-time performance. Figure 7 shows the convergence time performance of the various solvers as a function of grid size. Once again a power law fit is included for the best performer, LSOR. 7

45 Itertative Methods: Time vs Grid Size 1.E05 CPU Time (s) 1.E04 1.E03 1.E0 1.E01 1.E00 1.E-01 1.E-0 y = 1E-06x R = LSOR SOR GAUSS SEIDEL JACOBI LSQR GMRES BICG BICGSTAB QMR 1.E-03 1.E0 1.E03 1.E04 1.E05 1.E06 Grid Size Figure 7. Iterative Methods: Time Required for Convergence as a Function of Grid Size Extrapolation of the GMRES data from this plot suggests that GMRES could eventually perform better than any of the iterative methods at larger grid sizes. This can be seen by using the four data points for GMRES, which give a slope that indicates better performance for GMRES at larger grids; however GMRES memory requirements prohibited actual determination of results at larger grids. The potential of GMRES performing better than the other iterative methods at larger grids is somewhat diminished if the slope from the two largest grid sizes is used. Given the trends of the other iterative methods, decisions for GMRES were based on the data for the two largest grid sizes. QMR and BICG were stable for only the smallest grid size and had comparable convergence times to the other MATLAB iterative methods at the respective gird size. The use of a preconditioner at larger grid sizes could possibly help QMR and BICG to 8

46 stabilize and generate reasonable answers, but this was not done as this study was conducted without preconditioners. BICGSTAB was shown to be the slowest of the five MATLAB iterative methods. As expected, Jacobi and Gauss Seidel iterative methods performed poorly with respect to the other methods when comparing convergence time. SOR and LSOR performed better than all other iterative methods; LSOR had slightly better convergence times than SOR, but demanded more RAM than the other inhouse iterative methods for all grid sizes. It was difficult to determine how much memory was used for small grid sizes in MATLAB while the routines were in progress, yet for grid sizes larger than 17x17x34 the desktop computer did not have enough memory to complete the solution for any of the MATLAB methods. The results of the run-time memory performance of the various iterative methods are summarized in Figure 8. The values shown for SOR and LSOR are at the optimum over-relaxation factor, ω opt. Runtime memory requirements for MATLAB iterative methods at a grid size of 17x17x34 were determined for the three methods that worked at this grid size. These three points, which show similar memory requirements for MATLAB methods, are shown in Figure 8 and indicate memory requirement two orders of magnitude larger than the stationary iterative methods. This completes the first objective of this thesis, which was to test the performance of various solvers as a function of grid size. From these results, the best solver to be used on coarse and fine grids can be determined. 9

Contents. I The Basic Framework for Stationary Problems 1

Contents. I The Basic Framework for Stationary Problems 1 page v Preface xiii I The Basic Framework for Stationary Problems 1 1 Some model PDEs 3 1.1 Laplace s equation; elliptic BVPs... 3 1.1.1 Physical experiments modeled by Laplace s equation... 5 1.2 Other

More information

An Investigation into Iterative Methods for Solving Elliptic PDE s Andrew M Brown Computer Science/Maths Session (2000/2001)

An Investigation into Iterative Methods for Solving Elliptic PDE s Andrew M Brown Computer Science/Maths Session (2000/2001) An Investigation into Iterative Methods for Solving Elliptic PDE s Andrew M Brown Computer Science/Maths Session (000/001) Summary The objectives of this project were as follows: 1) Investigate iterative

More information

SELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND

SELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND Student Submission for the 5 th OpenFOAM User Conference 2017, Wiesbaden - Germany: SELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND TESSA UROIĆ Faculty of Mechanical Engineering and Naval Architecture, Ivana

More information

Homogenization and numerical Upscaling. Unsaturated flow and two-phase flow

Homogenization 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 information

PTE 519 Lecture Note Finite Difference Approximation (Model)

PTE 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 information

Hexahedral Mesh Refinement Using an Error Sizing Function

Hexahedral Mesh Refinement Using an Error Sizing Function Brigham Young University BYU ScholarsArchive All Theses and Dissertations 2011-06-01 Hexahedral Mesh Refinement Using an Error Sizing Function Gaurab Paudel Brigham Young University - Provo Follow this

More information

Application of Finite Volume Method for Structural Analysis

Application of Finite Volume Method for Structural Analysis Application of Finite Volume Method for Structural Analysis Saeed-Reza Sabbagh-Yazdi and Milad Bayatlou Associate Professor, Civil Engineering Department of KNToosi University of Technology, PostGraduate

More information

Multigrid Pattern. I. Problem. II. Driving Forces. III. Solution

Multigrid Pattern. I. Problem. II. Driving Forces. III. Solution Multigrid Pattern I. Problem Problem domain is decomposed into a set of geometric grids, where each element participates in a local computation followed by data exchanges with adjacent neighbors. The grids

More information

PhD Student. Associate Professor, Co-Director, Center for Computational Earth and Environmental Science. Abdulrahman Manea.

PhD Student. Associate Professor, Co-Director, Center for Computational Earth and Environmental Science. Abdulrahman Manea. Abdulrahman Manea PhD Student Hamdi Tchelepi Associate Professor, Co-Director, Center for Computational Earth and Environmental Science Energy Resources Engineering Department School of Earth Sciences

More information

Driven Cavity Example

Driven 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 information

CFD MODELING FOR PNEUMATIC CONVEYING

CFD MODELING FOR PNEUMATIC CONVEYING CFD MODELING FOR PNEUMATIC CONVEYING Arvind Kumar 1, D.R. Kaushal 2, Navneet Kumar 3 1 Associate Professor YMCAUST, Faridabad 2 Associate Professor, IIT, Delhi 3 Research Scholar IIT, Delhi e-mail: arvindeem@yahoo.co.in

More information

SPE Introduction Modern reservoir characterization tools can routinely generate very refined reservoir models containing millions of

SPE Introduction Modern reservoir characterization tools can routinely generate very refined reservoir models containing millions of SPE 79681 Dual Mesh Method in Upscaling Pascal Audigane, SPE, and Martin J Blunt, SPE, Centre for Petroleum Studies, Department of Earth Science and Engineering, Imperial College, London, UK. Copyright

More information

Introduction to Multigrid and its Parallelization

Introduction to Multigrid and its Parallelization Introduction to Multigrid and its Parallelization! Thomas D. Economon Lecture 14a May 28, 2014 Announcements 2 HW 1 & 2 have been returned. Any questions? Final projects are due June 11, 5 pm. If you are

More information

Generic framework for taking geological models as input for reservoir simulation

Generic 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 information

D036 Accelerating Reservoir Simulation with GPUs

D036 Accelerating Reservoir Simulation with GPUs D036 Accelerating Reservoir Simulation with GPUs K.P. Esler* (Stone Ridge Technology), S. Atan (Marathon Oil Corp.), B. Ramirez (Marathon Oil Corp.) & V. Natoli (Stone Ridge Technology) SUMMARY Over the

More information

Abstract. 1 Introduction

Abstract. 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 information

AMS526: Numerical Analysis I (Numerical Linear Algebra)

AMS526: Numerical Analysis I (Numerical Linear Algebra) AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 20: Sparse Linear Systems; Direct Methods vs. Iterative Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 26

More information

INTEGRATION OF LOCAL-GLOBAL UPSCALING AND GRID ADAPTIVITY FOR SIMULATION OF SUBSURFACE FLOW IN HETEROGENEOUS FORMATIONS

INTEGRATION OF LOCAL-GLOBAL UPSCALING AND GRID ADAPTIVITY FOR SIMULATION OF SUBSURFACE FLOW IN HETEROGENEOUS FORMATIONS INTEGRATION OF LOCAL-GLOBAL UPSCALING AND GRID ADAPTIVITY FOR SIMULATION OF SUBSURFACE FLOW IN HETEROGENEOUS FORMATIONS M. GERRITSEN AND J. V. LAMBERS Abstract. We propose a methodology, called multi-level

More information

1.2 Numerical Solutions of Flow Problems

1.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 information

Mid-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, 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 information

Fluent User Services Center

Fluent User Services Center Solver Settings 5-1 Using the Solver Setting Solver Parameters Convergence Definition Monitoring Stability Accelerating Convergence Accuracy Grid Independence Adaption Appendix: Background Finite Volume

More information

Variogram Inversion and Uncertainty Using Dynamic Data. Simultaneouos Inversion with Variogram Updating

Variogram Inversion and Uncertainty Using Dynamic Data. Simultaneouos Inversion with Variogram Updating Variogram Inversion and Uncertainty Using Dynamic Data Z. A. Reza (zreza@ualberta.ca) and C. V. Deutsch (cdeutsch@civil.ualberta.ca) Department of Civil & Environmental Engineering, University of Alberta

More information

S0432 NEW IDEAS FOR MASSIVELY PARALLEL PRECONDITIONERS

S0432 NEW IDEAS FOR MASSIVELY PARALLEL PRECONDITIONERS S0432 NEW IDEAS FOR MASSIVELY PARALLEL PRECONDITIONERS John R Appleyard Jeremy D Appleyard Polyhedron Software with acknowledgements to Mark A Wakefield Garf Bowen Schlumberger Outline of Talk Reservoir

More information

EFFICIENT PRODUCTION OPTIMIZATION USING FLOW NETWORK MODELS. A Thesis PONGSATHORN LERLERTPAKDEE

EFFICIENT 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 information

Toward reservoir simulation on geological grid models

Toward 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 information

ITERATIVE COLLISION RESOLUTION IN WIRELESS NETWORKS

ITERATIVE COLLISION RESOLUTION IN WIRELESS NETWORKS ITERATIVE COLLISION RESOLUTION IN WIRELESS NETWORKS An Undergraduate Research Scholars Thesis by KATHERINE CHRISTINE STUCKMAN Submitted to Honors and Undergraduate Research Texas A&M University in partial

More information

IMPROVING THE NUMERICAL ACCURACY OF HYDROTHERMAL RESERVOIR SIMULATIONS USING THE CIP SCHEME WITH THIRD-ORDER ACCURACY

IMPROVING THE NUMERICAL ACCURACY OF HYDROTHERMAL RESERVOIR SIMULATIONS USING THE CIP SCHEME WITH THIRD-ORDER ACCURACY 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

More information

Introduction to C omputational F luid Dynamics. D. Murrin

Introduction to C omputational F luid Dynamics. D. Murrin Introduction to C omputational F luid Dynamics D. Murrin Computational fluid dynamics (CFD) is the science of predicting fluid flow, heat transfer, mass transfer, chemical reactions, and related phenomena

More information

Example 24 Spring-back

Example 24 Spring-back Example 24 Spring-back Summary The spring-back simulation of sheet metal bent into a hat-shape is studied. The problem is one of the famous tests from the Numisheet 93. As spring-back is generally a quasi-static

More information

EXTENDING THE PETREL MODEL BUILDER FOR EDUCATIONAL AND RESEARCH PURPOSES. A Thesis OBIAJULU CHUKWUDI NWOSA

EXTENDING THE PETREL MODEL BUILDER FOR EDUCATIONAL AND RESEARCH PURPOSES. A Thesis OBIAJULU CHUKWUDI NWOSA EXTENDING THE PETREL MODEL BUILDER FOR EDUCATIONAL AND RESEARCH PURPOSES A Thesis by OBIAJULU CHUKWUDI NWOSA Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of

More information

What is Multigrid? They have been extended to solve a wide variety of other problems, linear and nonlinear.

What is Multigrid? They have been extended to solve a wide variety of other problems, linear and nonlinear. AMSC 600/CMSC 760 Fall 2007 Solution of Sparse Linear Systems Multigrid, Part 1 Dianne P. O Leary c 2006, 2007 What is Multigrid? Originally, multigrid algorithms were proposed as an iterative method to

More information

Application of Parametric NURBS Geometry to Mode Shape Identification and the Modal Assurance Criterion

Application of Parametric NURBS Geometry to Mode Shape Identification and the Modal Assurance Criterion Brigham Young University BYU ScholarsArchive All Theses and Dissertations 2012-04-12 Application of Parametric NURBS Geometry to Mode Shape Identification and the Modal Assurance Criterion Evan D. Selin

More information

Applied Reservoir Simulation. Course. Workshop Problems

Applied 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 information

THE DEVELOPMENT OF THE POTENTIAL AND ACADMIC PROGRAMMES OF WROCLAW UNIVERISTY OF TECH- NOLOGY ITERATIVE LINEAR SOLVERS

THE DEVELOPMENT OF THE POTENTIAL AND ACADMIC PROGRAMMES OF WROCLAW UNIVERISTY OF TECH- NOLOGY ITERATIVE LINEAR SOLVERS ITERATIVE LIEAR SOLVERS. Objectives The goals of the laboratory workshop are as follows: to learn basic properties of iterative methods for solving linear least squares problems, to study the properties

More information

FOUR WHAT S NEW IN THIS VERSION? 4.1 FLOW-3D Usability CHAPTER

FOUR WHAT S NEW IN THIS VERSION? 4.1 FLOW-3D Usability CHAPTER CHAPTER FOUR WHAT S NEW IN THIS VERSION? FLOW-3D v11.2.0 continues to streamline engineers simulation workflows by enabling them to more quickly set up simulations, avoid common errors, identify and enter

More information

Iterative methods for use with the Fast Multipole Method

Iterative methods for use with the Fast Multipole Method Iterative methods for use with the Fast Multipole Method Ramani Duraiswami Perceptual Interfaces and Reality Lab. Computer Science & UMIACS University of Maryland, College Park, MD Joint work with Nail

More information

Introduction to ANSYS CFX

Introduction to ANSYS CFX Workshop 03 Fluid flow around the NACA0012 Airfoil 16.0 Release Introduction to ANSYS CFX 2015 ANSYS, Inc. March 13, 2015 1 Release 16.0 Workshop Description: The flow simulated is an external aerodynamics

More information

Development 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 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 information

ENERGY-224 Reservoir Simulation Project Report. Ala Alzayer

ENERGY-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 information

CFD Project Workflow Guide

CFD Project Workflow Guide CFD Project Workflow Guide Contents Select a problem with known results for proof-of-concept testing... 1 Set up and run a coarse test case... 2 Select and calibrate numerical methods... 3 Minimize & quantify

More information

Simulating Sinkage & Trim for Planing Boat Hulls. A Fluent Dynamic Mesh 6DOF Tutorial

Simulating Sinkage & Trim for Planing Boat Hulls. A Fluent Dynamic Mesh 6DOF Tutorial Simulating Sinkage & Trim for Planing Boat Hulls A Fluent Dynamic Mesh 6DOF Tutorial 1 Introduction Workshop Description This workshop describes how to perform a transient 2DOF simulation of a planing

More information

Finite Element Method. Chapter 7. Practical considerations in FEM modeling

Finite Element Method. Chapter 7. Practical considerations in FEM modeling Finite Element Method Chapter 7 Practical considerations in FEM modeling Finite Element Modeling General Consideration The following are some of the difficult tasks (or decisions) that face the engineer

More information

STABILIZED FINITE ELEMENT METHODS FOR INCOMPRESSIBLE FLOWS WITH EMPHASIS ON MOVING BOUNDARIES AND INTERFACES

STABILIZED FINITE ELEMENT METHODS FOR INCOMPRESSIBLE FLOWS WITH EMPHASIS ON MOVING BOUNDARIES AND INTERFACES STABILIZED FINITE ELEMENT METHODS FOR INCOMPRESSIBLE FLOWS WITH EMPHASIS ON MOVING BOUNDARIES AND INTERFACES A THESIS SUBMITTED TO THE FACULTY OF THE GRADUATE SCHOOL OF THE UNIVERSITY OF MINNESOTA BY Marek

More information

Modeling Unsteady Compressible Flow

Modeling Unsteady Compressible Flow Tutorial 4. Modeling Unsteady Compressible Flow Introduction In this tutorial, FLUENT s density-based implicit solver is used to predict the timedependent flow through a two-dimensional nozzle. As an initial

More information

Modeling Strategies for Dynamic Finite Element Cask Analyses

Modeling Strategies for Dynamic Finite Element Cask Analyses Session A Package Analysis: Structural Analysis - Modeling Modeling Strategies for Dynamic Finite Element Cask Analyses Uwe Zencker, Günter Wieser, Linan Qiao, Christian Protz BAM Federal Institute for

More information

Integration of Geostatistical Modeling with History Matching: Global and Regional Perturbation

Integration of Geostatistical Modeling with History Matching: Global and Regional Perturbation Integration of Geostatistical Modeling with History Matching: Global and Regional Perturbation Oliveira, Gonçalo Soares Soares, Amílcar Oliveira (CERENA/IST) Schiozer, Denis José (UNISIM/UNICAMP) Introduction

More information

Calculate a solution using the pressure-based coupled solver.

Calculate a solution using the pressure-based coupled solver. Tutorial 19. Modeling Cavitation Introduction This tutorial examines the pressure-driven cavitating flow of water through a sharpedged orifice. This is a typical configuration in fuel injectors, and brings

More information

COMPUTATIONAL FLUID DYNAMICS USED IN THE DESIGN OF WATERBLAST TOOLING

COMPUTATIONAL FLUID DYNAMICS USED IN THE DESIGN OF WATERBLAST TOOLING 2015 WJTA-IMCA Conference and Expo November 2-4 New Orleans, Louisiana Paper COMPUTATIONAL FLUID DYNAMICS USED IN THE DESIGN OF WATERBLAST TOOLING J. Schneider StoneAge, Inc. Durango, Colorado, U.S.A.

More information

MATLAB. Advanced Mathematics and Mechanics Applications Using. Third Edition. David Halpern University of Alabama CHAPMAN & HALL/CRC

MATLAB. Advanced Mathematics and Mechanics Applications Using. Third Edition. David Halpern University of Alabama CHAPMAN & HALL/CRC Advanced Mathematics and Mechanics Applications Using MATLAB Third Edition Howard B. Wilson University of Alabama Louis H. Turcotte Rose-Hulman Institute of Technology David Halpern University of Alabama

More information

1 Exercise: 1-D heat conduction with finite elements

1 Exercise: 1-D heat conduction with finite elements 1 Exercise: 1-D heat conduction with finite elements Reading This finite element example is based on Hughes (2000, sec. 1.1-1.15. 1.1 Implementation of the 1-D heat equation example In the previous two

More information

Sweeping Jet Flushing Method for EDM

Sweeping Jet Flushing Method for EDM International Conference on Multidisciplinary Research & Practice P a g e 192 Sweeping Jet Flushing Method for EDM Kishor K Kulkarni # Production Engineering, ShivajiUniversity/Ashokrao Mane Group of Institution

More information

Optimizing Bio-Inspired Flow Channel Design on Bipolar Plates of PEM Fuel Cells

Optimizing 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 information

The Immersed Interface Method

The 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 information

SMART WELL MODELLING. Design, Scenarios and Optimisation

SMART WELL MODELLING. Design, Scenarios and Optimisation Page 1 Introduction Smart or complex wells are in increasing use by operators as reservoir environments become more challenging. The wells include a number of smart devices installed to achieve a variety

More information

Reducing Communication Costs Associated with Parallel Algebraic Multigrid

Reducing Communication Costs Associated with Parallel Algebraic Multigrid Reducing Communication Costs Associated with Parallel Algebraic Multigrid Amanda Bienz, Luke Olson (Advisor) University of Illinois at Urbana-Champaign Urbana, IL 11 I. PROBLEM AND MOTIVATION Algebraic

More information

Isotropic Porous Media Tutorial

Isotropic Porous Media Tutorial STAR-CCM+ User Guide 3927 Isotropic Porous Media Tutorial This tutorial models flow through the catalyst geometry described in the introductory section. In the porous region, the theoretical pressure drop

More information

A NURBS-BASED APPROACH FOR SHAPE AND TOPOLOGY OPTIMIZATION OF FLOW DOMAINS

A NURBS-BASED APPROACH FOR SHAPE AND TOPOLOGY OPTIMIZATION OF FLOW DOMAINS 6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 11 15 June 2018, Glasgow, UK A NURBS-BASED APPROACH FOR SHAPE AND TOPOLOGY OPTIMIZATION

More information

Estimation of a Piecewise Constant Function Using Reparameterized Level-Set Functions

Estimation of a Piecewise Constant Function Using Reparameterized Level-Set Functions Estimation of a iecewise Constant Function Using Reparameterized Level-Set Functions nga Berre 1,2, Martha Lien 2,1, and Trond Mannseth 2,1 1 Department of Mathematics, University of Bergen, Johs. Brunsgt.

More information

Module 1 Lecture Notes 2. Optimization Problem and Model Formulation

Module 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 information

COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF ORIFICE PLATE METERING SITUATIONS UNDER ABNORMAL CONFIGURATIONS

COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF ORIFICE PLATE METERING SITUATIONS UNDER ABNORMAL CONFIGURATIONS COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF ORIFICE PLATE METERING SITUATIONS UNDER ABNORMAL CONFIGURATIONS Dr W. Malalasekera Version 3.0 August 2013 1 COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF ORIFICE PLATE

More information

OPTIMIZATION FOR AUTOMATIC HISTORY MATCHING

OPTIMIZATION FOR AUTOMATIC HISTORY MATCHING INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 2, Supp, Pages 131 137 c 2005 Institute for Scientific Computing and Information OPTIMIZATION FOR AUTOMATIC HISTORY MATCHING Abstract. SHUGUANG

More information

Moving Interface Problems: Methods & Applications Tutorial Lecture II

Moving Interface Problems: Methods & Applications Tutorial Lecture II Moving Interface Problems: Methods & Applications Tutorial Lecture II Grétar Tryggvason Worcester Polytechnic Institute Moving Interface Problems and Applications in Fluid Dynamics Singapore National University,

More information

Introduction to Computational Fluid Dynamics Mech 122 D. Fabris, K. Lynch, D. Rich

Introduction to Computational Fluid Dynamics Mech 122 D. Fabris, K. Lynch, D. Rich Introduction to Computational Fluid Dynamics Mech 122 D. Fabris, K. Lynch, D. Rich 1 Computational Fluid dynamics Computational fluid dynamics (CFD) is the analysis of systems involving fluid flow, heat

More information

Chapter 7 Practical Considerations in Modeling. Chapter 7 Practical Considerations in Modeling

Chapter 7 Practical Considerations in Modeling. Chapter 7 Practical Considerations in Modeling CIVL 7/8117 1/43 Chapter 7 Learning Objectives To present concepts that should be considered when modeling for a situation by the finite element method, such as aspect ratio, symmetry, natural subdivisions,

More information

FEKO Mesh Optimization Study of the EDGES Antenna Panels with Side Lips using a Wire Port and an Infinite Ground Plane

FEKO Mesh Optimization Study of the EDGES Antenna Panels with Side Lips using a Wire Port and an Infinite Ground Plane FEKO Mesh Optimization Study of the EDGES Antenna Panels with Side Lips using a Wire Port and an Infinite Ground Plane Tom Mozdzen 12/08/2013 Summary This study evaluated adaptive mesh refinement in the

More information

Parallel High-Order Geometric Multigrid Methods on Adaptive Meshes for Highly Heterogeneous Nonlinear Stokes Flow Simulations of Earth s Mantle

Parallel High-Order Geometric Multigrid Methods on Adaptive Meshes for Highly Heterogeneous Nonlinear Stokes Flow Simulations of Earth s Mantle ICES Student Forum The University of Texas at Austin, USA November 4, 204 Parallel High-Order Geometric Multigrid Methods on Adaptive Meshes for Highly Heterogeneous Nonlinear Stokes Flow Simulations of

More information

Advanced Surface Based MoM Techniques for Packaging and Interconnect Analysis

Advanced Surface Based MoM Techniques for Packaging and Interconnect Analysis Electrical Interconnect and Packaging Advanced Surface Based MoM Techniques for Packaging and Interconnect Analysis Jason Morsey Barry Rubin, Lijun Jiang, Lon Eisenberg, Alina Deutsch Introduction Fast

More information

P257 Transform-domain Sparsity Regularization in Reconstruction of Channelized Facies

P257 Transform-domain Sparsity Regularization in Reconstruction of Channelized Facies P257 Transform-domain Sparsity Regularization in Reconstruction of Channelized Facies. azemi* (University of Alberta) & H.R. Siahkoohi (University of Tehran) SUMMARY Petrophysical reservoir properties,

More information

Improving Steering Module Efficiency for Incremental Loading Finite Element Numeric Models

Improving Steering Module Efficiency for Incremental Loading Finite Element Numeric Models Brigham Young University BYU ScholarsArchive All Theses and Dissertations 2006-03-22 Improving Steering Module Efficiency for Incremental Loading Finite Element Numeric Models Ryan L. Kitchen Brigham Young

More information

Revision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering. Introduction

Revision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering. Introduction Revision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering Introduction A SolidWorks simulation tutorial is just intended to illustrate where to

More information

Finite Element Convergence for Time-Dependent PDEs with a Point Source in COMSOL 4.2

Finite Element Convergence for Time-Dependent PDEs with a Point Source in COMSOL 4.2 Finite Element Convergence for Time-Dependent PDEs with a Point Source in COMSOL 4.2 David W. Trott and Matthias K. Gobbert Department of Mathematics and Statistics, University of Maryland, Baltimore County,

More information

Joe Wingbermuehle, (A paper written under the guidance of Prof. Raj Jain)

Joe Wingbermuehle, (A paper written under the guidance of Prof. Raj Jain) 1 of 11 5/4/2011 4:49 PM Joe Wingbermuehle, wingbej@wustl.edu (A paper written under the guidance of Prof. Raj Jain) Download The Auto-Pipe system allows one to evaluate various resource mappings and topologies

More information

2 T. x + 2 T. , T( x, y = 0) = T 1

2 T. x + 2 T. , T( x, y = 0) = T 1 LAB 2: Conduction with Finite Difference Method Objective: The objective of this laboratory is to introduce the basic steps needed to numerically solve a steady state two-dimensional conduction problem

More information

PARALLELIZATION OF POTENTIAL FLOW SOLVER USING PC CLUSTERS

PARALLELIZATION OF POTENTIAL FLOW SOLVER USING PC CLUSTERS Proceedings of FEDSM 2000: ASME Fluids Engineering Division Summer Meeting June 11-15,2000, Boston, MA FEDSM2000-11223 PARALLELIZATION OF POTENTIAL FLOW SOLVER USING PC CLUSTERS Prof. Blair.J.Perot Manjunatha.N.

More information

TPG4160 Reservoir simulation, Building a reservoir model

TPG4160 Reservoir simulation, Building a reservoir model TPG4160 Reservoir simulation, Building a reservoir model Per Arne Slotte Week 8 2018 Overview Plan for the lectures The main goal for these lectures is to present the main concepts of reservoir models

More information

CHAPTER 4. Numerical Models. descriptions of the boundary conditions, element types, validation, and the force

CHAPTER 4. Numerical Models. descriptions of the boundary conditions, element types, validation, and the force CHAPTER 4 Numerical Models This chapter presents the development of numerical models for sandwich beams/plates subjected to four-point bending and the hydromat test system. Detailed descriptions of the

More information

A Geostatistical and Flow Simulation Study on a Real Training Image

A Geostatistical and Flow Simulation Study on a Real Training Image A Geostatistical and Flow Simulation Study on a Real Training Image Weishan Ren (wren@ualberta.ca) Department of Civil & Environmental Engineering, University of Alberta Abstract A 12 cm by 18 cm slab

More information

McNair Scholars Research Journal

McNair Scholars Research Journal McNair Scholars Research Journal Volume 2 Article 1 2015 Benchmarking of Computational Models against Experimental Data for Velocity Profile Effects on CFD Analysis of Adiabatic Film-Cooling Effectiveness

More information

INVERSE PROBLEMS IN GROUNDWATER MODELING

INVERSE PROBLEMS IN GROUNDWATER MODELING INVERSE PROBLEMS IN GROUNDWATER MODELING Theory and Applications of Transport in Porous Media Series Editor: Jacob Bear, Technion - Israel Institute of Technology, Haifa, Israel Volume 6 The titles published

More information

Using a Single Rotating Reference Frame

Using 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 information

NIA CFD Seminar, October 4, 2011 Hyperbolic Seminar, NASA Langley, October 17, 2011

NIA CFD Seminar, October 4, 2011 Hyperbolic Seminar, NASA Langley, October 17, 2011 NIA CFD Seminar, October 4, 2011 Hyperbolic Seminar, NASA Langley, October 17, 2011 First-Order Hyperbolic System Method If you have a CFD book for hyperbolic problems, you have a CFD book for all problems.

More information

The Development of a 3D Self-Adaptive Goal-Oriented hp-finite Element Software for Simulations of Resistivity Logging Instruments

The Development of a 3D Self-Adaptive Goal-Oriented hp-finite Element Software for Simulations of Resistivity Logging Instruments A Progress Report on: The Development of a 3D Self-Adaptive Goal-Oriented hp-finite Element Software for Simulations of Resistivity Logging Instruments David Pardo (dzubiaur@yahoo.es), L. Demkowicz, C.

More information

The 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

The 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 information

Recent Developments in the Design and Optimization of Constant Force Electrical Contacts

Recent Developments in the Design and Optimization of Constant Force Electrical Contacts Recent Developments in the Design and Optimization of Constant orce Electrical Contacts John C. Meaders * Stephen P. Harston Christopher A. Mattson Brigham Young University Provo, UT, 84602, USA Abstract

More information

The determination of the correct

The determination of the correct SPECIAL High-performance SECTION: H i gh-performance computing computing MARK NOBLE, Mines ParisTech PHILIPPE THIERRY, Intel CEDRIC TAILLANDIER, CGGVeritas (formerly Mines ParisTech) HENRI CALANDRA, Total

More information

1 2 (3 + x 3) x 2 = 1 3 (3 + x 1 2x 3 ) 1. 3 ( 1 x 2) (3 + x(0) 3 ) = 1 2 (3 + 0) = 3. 2 (3 + x(0) 1 2x (0) ( ) = 1 ( 1 x(0) 2 ) = 1 3 ) = 1 3

1 2 (3 + x 3) x 2 = 1 3 (3 + x 1 2x 3 ) 1. 3 ( 1 x 2) (3 + x(0) 3 ) = 1 2 (3 + 0) = 3. 2 (3 + x(0) 1 2x (0) ( ) = 1 ( 1 x(0) 2 ) = 1 3 ) = 1 3 6 Iterative Solvers Lab Objective: Many real-world problems of the form Ax = b have tens of thousands of parameters Solving such systems with Gaussian elimination or matrix factorizations could require

More information

Recent Advances in MSC/PATRAN Pre-Processing Software Allows Modeling of Complex Automotive Lamp Designs

Recent Advances in MSC/PATRAN Pre-Processing Software Allows Modeling of Complex Automotive Lamp Designs Recent Advances in MSC/PATRAN Pre-Processing Software Allows Modeling of Complex Automotive Lamp Designs William I. Moore, Eric S. Donovan and Christopher R. Powers Delphi Interior and Lighting Systems

More information

CFD Post-Processing of Rampressor Rotor Compressor

CFD Post-Processing of Rampressor Rotor Compressor Gas Turbine Industrial Fellowship Program 2006 CFD Post-Processing of Rampressor Rotor Compressor Curtis Memory, Brigham Young niversity Ramgen Power Systems Mentor: Rob Steele I. Introduction Recent movements

More information

Fluid flow modelling with seismic cluster analysis

Fluid flow modelling with seismic cluster analysis Fluid flow modelling with seismic cluster analysis Fluid flow modelling with seismic cluster analysis Laurence R. Bentley, Xuri Huang 1 and Claude Laflamme 2 ABSTRACT Cluster analysis is used to construct

More information

FLUENT Secondary flow in a teacup Author: John M. Cimbala, Penn State University Latest revision: 26 January 2016

FLUENT Secondary flow in a teacup Author: John M. Cimbala, Penn State University Latest revision: 26 January 2016 FLUENT Secondary flow in a teacup Author: John M. Cimbala, Penn State University Latest revision: 26 January 2016 Note: These instructions are based on an older version of FLUENT, and some of the instructions

More information

ETNA Kent State University

ETNA Kent State University Electronic Transactions on Numerical Analysis. Volume, 2, pp. 92. Copyright 2,. ISSN 68-963. ETNA BEHAVIOR OF PLANE RELAXATION METHODS AS MULTIGRID SMOOTHERS IGNACIO M. LLORENTE AND N. DUANE MELSON Abstract.

More information

D-Optimal Designs. Chapter 888. Introduction. D-Optimal Design Overview

D-Optimal Designs. Chapter 888. Introduction. D-Optimal Design Overview Chapter 888 Introduction This procedure generates D-optimal designs for multi-factor experiments with both quantitative and qualitative factors. The factors can have a mixed number of levels. For example,

More information

NUMERICAL INVESTIGATION OF THE FLOW BEHAVIOR INTO THE INLET GUIDE VANE SYSTEM (IGV)

NUMERICAL 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 information

Simulation of In-Cylinder Flow Phenomena with ANSYS Piston Grid An Improved Meshing and Simulation Approach

Simulation of In-Cylinder Flow Phenomena with ANSYS Piston Grid An Improved Meshing and Simulation Approach Simulation of In-Cylinder Flow Phenomena with ANSYS Piston Grid An Improved Meshing and Simulation Approach Dipl.-Ing. (FH) Günther Lang, CFDnetwork Engineering Dipl.-Ing. Burkhard Lewerich, CFDnetwork

More information

Automated All-Quadrilateral Mesh Adaptation through Refinement and Coarsening

Automated All-Quadrilateral Mesh Adaptation through Refinement and Coarsening Brigham Young University BYU ScholarsArchive All Theses and Dissertations 2009-06-05 Automated All-Quadrilateral Mesh Adaptation through Refinement and Coarsening Bret D. Anderson Brigham Young University

More information

Challenge Problem 5 - The Solution Dynamic Characteristics of a Truss Structure

Challenge Problem 5 - The Solution Dynamic Characteristics of a Truss Structure Challenge Problem 5 - The Solution Dynamic Characteristics of a Truss Structure In the final year of his engineering degree course a student was introduced to finite element analysis and conducted an assessment

More information

f xx + f yy = F (x, y)

f xx + f yy = F (x, y) Application of the 2D finite element method to Laplace (Poisson) equation; f xx + f yy = F (x, y) M. R. Hadizadeh Computer Club, Department of Physics and Astronomy, Ohio University 4 Nov. 2013 Domain

More information

C. A. D. Fraga Filho 1,2, D. F. Pezzin 1 & J. T. A. Chacaltana 1. Abstract

C. 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 information

Thermal Reservoir Simulator for Microsoft Windows TM

Thermal 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 information

Electromagnetic migration of marine CSEM data in areas with rough bathymetry Michael S. Zhdanov and Martin Čuma*, University of Utah

Electromagnetic migration of marine CSEM data in areas with rough bathymetry Michael S. Zhdanov and Martin Čuma*, University of Utah Electromagnetic migration of marine CSEM data in areas with rough bathymetry Michael S. Zhdanov and Martin Čuma*, University of Utah Summary In this paper we present a new approach to the interpretation

More information