Poroelasticity simulation in CLAWPACK for geophysical and biological applications
|
|
- Rosemary Fox
- 6 years ago
- Views:
Transcription
1 Poroelasticity simulation in CLAWPACK for geophysical and biological applications October 27, 2012 Joint work with Randall J. LeVeque (University of Washington) and Yvonne Ou (University of Delaware) Supported in part by grants from the NIH and NSF
2 Outline 1 Poroelasticity 2 Numerical solution 3 Results Convergence studies Qualitative checks Mapped grid results 4 Summary and future work
3 Outline 1 Poroelasticity 2 Numerical solution 3 Results Convergence studies Qualitative checks Mapped grid results 4 Summary and future work
4 Poroelasticity theory Poroelasticity: study of mechanics of fluid-filled porous solids Originally developed by Maurice Biot in 1930s-1960s for soil and rock Major early interest from oil industry Recent interest for monitoring underground fluid injection (e.g. carbon sequestration) Recently applied to bone as well Understanding wave propagation in bone is original motivator for this work Pumice stone. Transverse section of cortical bone. Images courtesy Wikimedia Commons.
5 Equations of poroelasticity Glossing over a lot of details, can model poroelasticity as a first-order linear system of PDEs: t Q + A x Q + B z Q = DQ, Q = [ p σ xx σ zz σ xz v x v z q x q z ] T p is fluid pressure; σ is solid stress tensor; v is solid velocity; q is fluid flow rate WRT solid Left side is classic hyperbolic system Right side introduces dissipation, through viscous drag as fluid flows through pores Many important poroelastic materials are anisotropic A and B depend on axes orientation Show system matrices
6 Poroelastic waves Poroelasticity equations look like elastodynamics + acoustics. Three families of propagating waves: 1 Fast P-waves, where fluid and solid move (roughly) parallel to propagation direction and in phase 2 S-waves, where fluid and solid move transverse to propagation direction 3 Slow P-waves, where fluid and solid move (roughly) parallel to propation direction but 180 degrees out of phase Slow P-waves involve large motions of fluid relative to solid heavily damped by viscosity System also supports non-wave-like diffusive slow mode where fluid seeps through pores due to pressure gradient
7 Wave structure of viscous orthotropic poroelasticity Source term DQ causes dissipation and dispersion Anisotropy also causes wave speeds to differ depending on propagation direction Plots below: phase velocity vs. frequency and direction for orthotropic layered sandstone Phase velocity and decay rate for waves in the X direction Dashed lines are speeds without dissipation Frequency (Hz) Phase velocity (m/s) Decay rate (1/m) Phase velocity (m/s) versus direction, all wave families (no dissipation) Fast P-wave S-wave315 Slow P-wave 0
8 Stiffness of relaxation term t Q + A x Q + B z Q = DQ Source term DQ has its own intrinsic time scales May be stiff depending on problem solved Can expect difficulties with solution Source term is of relaxation type, so expect solution to be close to that of reduced system, t u + A r x u + B r z u = 0 obtained by restricting to null space N (D) A r = ΠAG, where G maps reduced variables to full variables and Π maps full to reduced Reduced system satisfies subcharacteristic condition wave speeds not faster than full system Can at least expect qualitatively correct solution
9 Outline 1 Poroelasticity 2 Numerical solution 3 Results Convergence studies Qualitative checks Mapped grid results 4 Summary and future work
10 Finite volume method We use high-resolution wave propagation finite volume methods (FVM) to solve the poroelasticity equations Overview of approach for a hyperbolic problem: Discretize domain into grid of cells Assume we know average value of solution on each cell At each time step: 1 Reconstruct approximate solution from cell averages (usually discontinuous between cells) 2 Solve Riemann problem (Cauchy problem with piecewise constant initial data) at each cell interface 3 Evolve Riemann solution for time step t 4 Average new solution onto grid
11 Advantages of high-resolution FVM Riemann solution explicitly reproduces full gamut of wave behavior Can greatly reduce numerical dispersion through wave limiters Limiters also improve accuracy (in typical computations) at low cost Recent advanced limiters can even increase order of accuracy Easy to incorporate non-cartesian grids Handles low-quality meshes gracefully
12 Fluid-poroelastic interface Riemann problem approach allows fairly natural handling of material discontinuities, special interface conditions Example of interest: interface between fluid and poroelastic medium Have set of interface conditions: Fluid continuity equation: Newton s third law: q fluid n = (v poro + q poro ) n p fluid n = τ poro n Hydraulic permeability condition: p poro p fluid = 1 K q poro n Can incorporate these as constraints on Riemann solution
13 Fluid-poroelastic interface Some algorithmic complexity introduced at fluid-poroelastic interface High-resolution FVM originates with conservation laws, e.g. q t + F(q) x + G(q) z = 0 Some ideas most easily written in terms of fluxes F and G Can t define single flux between materials with different physics! Higher-order correction fluxes harder to formulate, some accuracy lost Not a serious problem in practice nonsmooth solution at interface anyway
14 Source term Source term handled via operator splitting Q t = DQ solved exactly with matrix exponential Best accuracy available for this part of system No stability restriction Used Strang splitting scheme: 1 Take 1 2 t step of source term 2 Take full t step of wave propagation terms 3 Take another 1 2 t step of source term Formally second-order accurate as t 0, but will see what we actually get for realistic step size
15 CLAWPACK Solved poroelasticity equations using CLAWPACK CLAWPACK: Conservation LAWs PACKage Can really handle any hyperbolic system, not just conservation laws High-resolution finite volume package for wave propagation Low memory overhead, parallel (multicore, PETSc) Operator splitting built in for source terms Adaptive mesh refinement available too (Berger-Colella-Oliger approach, AMRCLAW) Handles code that is common across all high-resolution FVM solvers User only needs to write routines for Riemann solve, source terms, grid mapping (if used)
16 Outline 1 Poroelasticity 2 Numerical solution 3 Results Convergence studies Qualitative checks Mapped grid results 4 Summary and future work
17 Overview of results Two general classes of test problems run so far: 1 Convergence studies comparing against known analytic solutions Advantage: Can quantitatively measure accuracy, convergence rate Disadvantage: Limited library of solutions to compare against 2 Qualitative sanity-check problems and eyeball norm comparisons to published solutions Advantage: Useful when analytic solution not available, good for ruling out some types of bug Disadvantage: Not very precise Most results are for Cartesian grids, but also have mapped grid results with more interesting geometry
18 Convergence studies Conducted convergence studies against analytic solutions Solutions used: plane waves of the form Q(x, z, t) = V exp(i(k x x + k z z ωt)) with some real ω specified. Important to test with waves propagating in variety of directions θ wave Also need to test variety of material principal directions θ mat For each (θ wave, θ mat ) pair, need to sweep over grid size to check convergence z (meters) v z (m/s) at time t = z θ mat θ wave 1 x Wave direction x (meters) Snapshot of plane wave solution showing θ wave, principal axes, and θ mat
19 Convergence results: inviscid First convergence study: ignore viscous dissipation, validate hyperbolic solver by itself Results are generally good Slow P wave was underresolved on coarse grids worse apparent performance Error measured using energy max-norm, relative to amplitude in energy norm of true solution Convergence rate Error on finest grid Wave family Best Worst Best Worst Fast P S Slow P
20 Convergence results: viscous high-frequency Ran more detailed convergence studies at lower and upper frequency ranges of Biot system High frequency chosen: 10 khz Ran slow P waves on different grid from others because of extremely rapid damping Results are good, though not as good as inviscid for the fast P and S waves Convergence rate Error on finest grid Wave family Best Worst Best Worst Fast P S Slow P
21 Convergence results: viscous low-frequency Also examined results at low frequency: 10 Hz Again ran slow P waves on different grid because of rapid damping Results are not so good Convergence rate Error on finest grid Wave family Best Worst Best Worst Fast P S Slow P
22 Convergence results: discussion What s going on here? Trouble only for timesteps comparable to relaxation time or longer Asymptotic error estimates only good as t 0 not inconsistent with bad results for Strang at large t. Literature suggests hyperbolic systems with stiff relaxation terms are hard to model accurately Dissipation causes change in Riemann solution structure at longer times Solution structure approaches that of reduced system Reduced system waves are blurred into erf-shapes, so reduced wave structure is not quite right either Possible ways to improve accuracy: Modeling full Riemann solution structure explicitly Using method of lines approach with appropriate integrator
23 Poroelastic code validation: two-material test case Test case: wave reflections from a material interface, from de la Puente et al. (2008) Excitation is a point source with a Ricker wavelet profile in time Forcing acts on σ z and fluid pressure terms with equal magnitude and opposite sign Source is located in shale overlying sandstone bed Material properties isotropic for this case, viscosity ignored AMR used to capture fine details
24 Poroelastic code validation: two-material test case v z (m/s) at time t = Shale a) b) z (meters) Sandstone x (meters) Figure 9. a Snapshot of the y-direction of the solid particle velocity at t 0.25 s. ers by empty circles. b Zoomed region showing the fine mesh required to capture z component of matrix velocity field. Left: CLAWPACK. Right: de la Puente. Rectangular outlines indicate boundaries of AMR grids. Note: Different value-to-shade maps on each plot. Downloaded 02 Jan 2011 to Redistribution subject to SEG lice
25 Poroelastic code validation: two-material test case DG methods for poroelastic media T93 15 a) u (m/s) v z (m/s) at gauge 1 8 Receiver 1 (y = 750m) time Time (s) c) ADER DG O5 FD Residuals d) ADER DG O5 FD 6 Residuals Left: CLAWPACK, right: 4 de la Puente Time (s) Time-history of matrix z velocity at topmost gauge. b) v (m/s) ADER DG O5 FD Residuals Receiver 1 (y = 750m) ADER DG O5 FD Residuals 1 2 u (m/s) 0 v (m/s)
26 Cylindrical scatterer More challenging test case: poroelastic cylinder in fluid bath Time-harmonic behavior of cylinder scattering acoustic waves Problem is not in stiff regime Tests mapped grid, fluid-poro interface Analytic solution provided by Prof. Yvonne Ou, University of Delaware Results: Observed 2nd-order convergence in 1- and 2-norms, 1st-order in max-norm Accuracy degraded by fluid-poro interface, highly distorted grid map
27 Cylindrical scatterer results 0.06 grids at time t = e z (meters) x (meters) Results for cylindrical scatterer in biological fluid. Incoming wavelength is 1/4 cylinder diameter. Left: fluid pressure ( grid); right: grid illustrating mapping Logically rectangular circle map from Calhoun, Helzel, and LeVeque (2007)
28 Simple bone model Have some quick results with simplified model of human femur Model uses mapped grid for concentric circular cross-sections Could do more complex geometry as long as mapping subroutine is available Have multiple materials, fluid-poroelastic interface: Inner core of cancellous bone Outer annulus of cortical bone Surrounding domain is water Model is not in stiff regime relaxation time is much longer than time step
29 Simple bone model p (Pa) at time t = e-05 grids at time t = z (meters) z (meters) x (meters) x (meters) Results for simplified model of human femur in water, struck by acoustic wave with Gaussian profile. Left: fluid pressure ( grid); right: grid illustrating mapping
30 Outline 1 Poroelasticity 2 Numerical solution 3 Results Convergence studies Qualitative checks Mapped grid results 4 Summary and future work
31 Summary Poroelasticity is a rich and complex system, with a wide variety of behaviors Have developed poroelasticity solver, validated against known solutions Solver capabilities: Multiple materials Fluid-poroelastic interfaces Mapped grids for moderately complex geometries Parallel execution (thanks to CLAWPACK framework) Convergence and accuracy are good away from stiff regime Convergence degraded in stiff regime, but magnitude of error is generally not bad Good results obtained for mapped grid models
32 Future work Deal with convergence problems in stiff regime More advanced Riemann solver? Exponential integrator for method of lines scheme Implement property averaging across material boundaries for geometry too complex for mapped grids? Look at micro-scale modeling to obtain poroelastic properties Examine cylindrical scatterer in more detail resonance, surface waves Extend to 3D
33 Poroelasticity system matrices [ ] [ ] 0 Ăsv 0 0 Ă = n xa + n zb =, D = Ă vs 0 0 D v n xα 1 M n zα 3 M n xm n zm Ă sv = n xc u 11 n zc u 13 n xα 1 M n zα 1 M n xc u 13 n zc u 33 n xα 3 M n zα 3 M n zc u 55 n xc u ρ n f m x n 1 m x 0 n 1 1 z 1 ρ 1 n f m z 0 n 3 m z n 3 Ă vs = 3 x 3 3 ρ ρ n x n f ρ x 0 n f 1 z 1 1 ρ ρ n z 0 n f ρ z n f 3 x 3 3 ρ 0 0 f η 0 1 κ 1 ρ f η D v = 3 κ ρη 0 1 κ ρη 3 κ 3 Note: these are in principal axes of orthotropic material Back to main
High-resolution finite volume methods for hyperbolic PDEs on manifolds
High-resolution finite volume methods for hyperbolic PDEs on manifolds Randall J. LeVeque Department of Applied Mathematics University of Washington Supported in part by NSF, DOE Overview High-resolution
More informationWeight-adjusted DG methods for elastic wave propagation in arbitrary heterogeneous media
Weight-adjusted DG methods for elastic wave propagation in arbitrary heterogeneous media Jesse Chan Department of Computational and Applied Math, Rice University ICOSAHOM 2018 July 12, 2018 Chan (CAAM)
More informationModeling seismic wave propagation during fluid injection in a fractured network: Effects of pore fluid pressure on time-lapse seismic signatures
Modeling seismic wave propagation during fluid injection in a fractured network: Effects of pore fluid pressure on time-lapse seismic signatures ENRU LIU, SERAFEIM VLASTOS, and XIANG-YANG LI, Edinburgh
More informationClawpack Tutorial Part I
Clawpack Tutorial Part I Randall J. LeVeque Applied Mathematics University of Washington Conservation Laws Package www.clawpack.org (pdf s will be posted and green links can be clicked) Some collaborators
More informationNonlinear Wave Theory for Transport Phenomena
JOSO 2016 March 9-11 2015 Nonlinear Wave Theory for Transport Phenomena ILYA PESHKOV CHLOE, University of Pau, France EVGENIY ROMENSKI Sobolev Institute of Mathematics, Novosibirsk, Russia MICHAEL DUMBSER
More informationHigh-Resolution Finite Volume Methods and Adaptive Mesh Refinement
High-Resolution Finite Volume Methods and Adaptive Mesh Refinement Randall J. LeVeque Department of Applied Mathematics University of Washington CLAWPACK and TsunamiClaw Software http://www.amath.washington.edu/~claw
More informationWave Propagation in Poroelastic Materials
Wave Propagation in Poroelastic Materials M. Yvonne Ou Department of Mathematical Sciences University of Delaware Newark, DE 19711 mou@math.udel.edu ructure Interaction Probs, May 28 211 Acknowledgement:
More informationCalculating Storm Surge and Other Coastal Hazards Using Geoclaw
Calculating Storm Surge and Other Coastal Hazards Using Geoclaw Kyle T. Mandli Department of Applied Mathematics University of Washington Seattle, WA, USA Modeling and Computations of Shallow-Water Coastal
More informationWe briefly discuss two examples for solving wave propagation type problems with finite differences, the acoustic and the seismic problem.
Excerpt from GEOL557 Numerical Modeling of Earth Systems by Becker and Kaus 2016 1 Wave propagation Figure 1: Finite difference discretization of the 2D acoustic problem. We briefly discuss two examples
More informationDifferential Acoustic Resonance Spectroscopy Analysis of Fluids in Porous Media
http://ijopaar.com; 2016 Vol. 2(1); pp. 22-30 Differential Acoustic Resonance Spectroscopy Analysis of Fluids in Porous Media Dr.Mohammad Miyan Associate Professor, Department of Mathematics, Shia P.G.College,
More information6298 Stress induced azimuthally anisotropic reservoir - AVO modeling
6298 Stress induced azimuthally anisotropic reservoir - AVO modeling M. Brajanovski* (Curtin University of Technology), B. Gurevich (Curtin University of Technology), D. Nadri (CSIRO) & M. Urosevic (Curtin
More informationEfficient numerical solution of the Biot poroelasticity system in multilayered domains
Efficient numerical solution of the Biot poroelasticity system Anna Naumovich Oleg Iliev Francisco Gaspar Fraunhofer Institute for Industrial Mathematics Kaiserslautern, Germany, Spain Workshop on Model
More informationDiffusion / Parabolic Equations. PHY 688: Numerical Methods for (Astro)Physics
Diffusion / Parabolic Equations Summary of PDEs (so far...) Hyperbolic Think: advection Real, finite speed(s) at which information propagates carries changes in the solution Second-order explicit methods
More informationHp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme
Hp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme Aravind Balan, Michael Woopen and Georg May AICES Graduate School, RWTH Aachen University, Germany 22nd AIAA Computational
More informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationApplying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models
0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle,
More informationAdvanced numerical methods for nonlinear advectiondiffusion-reaction. Peter Frolkovič, University of Heidelberg
Advanced numerical methods for nonlinear advectiondiffusion-reaction equations Peter Frolkovič, University of Heidelberg Content Motivation and background R 3 T Numerical modelling advection advection
More informationHydrogeophysics - Seismics
Hydrogeophysics - Seismics Matthias Zillmer EOST-ULP p. 1 Table of contents SH polarized shear waves: Seismic source Case study: porosity of an aquifer Seismic velocities for porous media: The Frenkel-Biot-Gassmann
More informationHydraulic properties of porous media
PART 5 Hydraulic properties of porous media Porosity Definition: Void space: n V void /V total total porosity e V void /V solid Primary porosity - between grains Secondary porosity - fracture or solution
More informationConservation Laws and Finite Volume Methods
Conservation Laws and Finite Volume Methods AMath 574 Winter Quarter, 2011 Randall J. LeVeque Applied Mathematics University of Washington January 3, 2011 R.J. LeVeque, University of Washington AMath 574,
More informationIs My CFD Mesh Adequate? A Quantitative Answer
Is My CFD Mesh Adequate? A Quantitative Answer Krzysztof J. Fidkowski Gas Dynamics Research Colloqium Aerospace Engineering Department University of Michigan January 26, 2011 K.J. Fidkowski (UM) GDRC 2011
More informationDownloaded 08/30/13 to Redistribution subject to SEG license or copyright; see Terms of Use at
Modeling the effect of pores and cracks interactions on the effective elastic properties of fractured porous rocks Luanxiao Zhao*, De-hua Han, Qiuliang Yao and Fuyong Yan, University of Houston; Mosab
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationA finite difference Poisson solver for irregular geometries
ANZIAM J. 45 (E) ppc713 C728, 2004 C713 A finite difference Poisson solver for irregular geometries Z. Jomaa C. Macaskill (Received 8 August 2003, revised 21 January 2004) Abstract The motivation for this
More informationAttenuation compensation in least-squares reverse time migration using the visco-acoustic wave equation
Attenuation compensation in least-squares reverse time migration using the visco-acoustic wave equation Gaurav Dutta, Kai Lu, Xin Wang and Gerard T. Schuster, King Abdullah University of Science and Technology
More informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
More informationAcoustics in Porous Media.
Acoustics in Porous Media. JUAN E. SANTOS work in collaboration with J. M. Carcione, S. Picotti, P. M. Gauzellino and R. Martinez Corredor. Department of Mathematics, Purdue University, W. Lafayette, Indiana,
More informationWave Propagation in Anisotropic, Anelastic, Porous and Electromagnetic Media
SECTION I. SEISMIC EXPLORATION Volume 38 WAVE FIELDS IN REAL MEDIA: Wave Propagation in Anisotropic, Anelastic, Porous and Electromagnetic Media (SECOND EDITION, REVISED AND EXTENDED) by Jose M. CARCIONE
More informationIn all of the following equations, is the coefficient of permeability in the x direction, and is the hydraulic head.
Groundwater Seepage 1 Groundwater Seepage Simplified Steady State Fluid Flow The finite element method can be used to model both steady state and transient groundwater flow, and it has been used to incorporate
More informationComputational Analysis of an Imploding Gas:
1/ 31 Direct Numerical Simulation of Navier-Stokes Equations Stephen Voelkel University of Notre Dame October 19, 2011 2/ 31 Acknowledges Christopher M. Romick, Ph.D. Student, U. Notre Dame Dr. Joseph
More informationFormation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas )
Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas ) Yasutomo ISHII and Andrei SMOLYAKOV 1) Japan Atomic Energy Agency, Ibaraki 311-0102, Japan 1) University
More informationTowards Modelling Elastic and Viscoelastic Seismic Wave Propagation in Boreholes
Towards Modelling Elastic and Viscoelastic Seismic Wave Propagation in Boreholes NA WANG, DONG SHI, BERND MILKEREIT Department of Physics, University of Toronto, Toronto, Canada M5S 1A7 Summary We are
More informationTheory - SBRC fundamentals
Chapter 2 Theory - SBRC fundamentals The approach proposed to estimate hydraulic properties of rocks using microseismicity is called Seismicity Based Reservoir Characterization (SBRC). It uses a spatio-temporal
More informationABSTRACT 1. INTRODUCTION
ABSTRACT In this paper the three-dimensional transient wave propagation is investigated due to a point force applied at the interface of a fluid and a poroelastic solid. Using the total response, it is
More informationAn efficient wave extrapolation method for tilted orthorhombic media using effective ellipsoidal models
An efficient wave extrapolation method for tilted orthorhombic media using effective ellipsoidal models Item Type Conference Paper Authors Waheed Umair bin; Alkhalifah Tariq Ali Eprint version Pre-print
More informationinter.noise 2000 The 29th International Congress and Exhibition on Noise Control Engineering August 2000, Nice, FRANCE
Copyright SFA - InterNoise 2000 1 inter.noise 2000 The 29th International Congress and Exhibition on Noise Control Engineering 27-30 August 2000, Nice, FRANCE I-INCE Classification: 2.0 TLM MODEL FOR SOUND
More informationLawrence Berkeley National Laboratory
Lawrence Berkeley National Laboratory Peer Reviewed Title: Fracture permeability and seismic wave scattering--poroelastic linear-slip interface model for heterogeneous fractures Author: Nakagawa, S. Publication
More informationCrosswell tomography imaging of the permeability structure within a sandstone oil field.
Crosswell tomography imaging of the permeability structure within a sandstone oil field. Tokuo Yamamoto (1), and Junichi Sakakibara (2) (1) University of Miami and Yamamoto Engineering Corporation, (2)
More informationUpscaling Wave Computations
Upscaling Wave Computations Xin Wang 2010 TRIP Annual Meeting 2 Outline 1 Motivation of Numerical Upscaling 2 Overview of Upscaling Methods 3 Future Plan 3 Wave Equations scalar variable density acoustic
More informationA Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials
Dublin, October 2010 A Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials FracMan Technology Group Dr Mark Cottrell Presentation Outline Some Physical
More informationEffects of Fracture Parameters in an Anisotropy Model on P-Wave Azimuthal Amplitude Responses
PROC. ITB Eng. Science Vol. 38 B, No. 2, 2006, 159-170 159 Effects of Fracture Parameters in an Anisotropy Model on P-Wave Azimuthal Amplitude Responses Fatkhan Program Studi Teknik Geofisika FIKTM-ITB
More informationTENSOR RELATIVE PERMEABILITIES: ORIGINS, MODELING AND NUMERICAL DISCRETIZATION
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 9, Number 3, Pages 701 724 c 2012 Institute for Scientific Computing and Information TENSOR RELATIVE PERMEABILITIES: ORIGINS, MODELING AND
More informationChapter 1 Direct Modeling for Computational Fluid Dynamics
Chapter 1 Direct Modeling for Computational Fluid Dynamics Computational fluid dynamics (CFD) is a scientific discipline, which aims to capture fluid motion in a discretized space. The description of the
More informationPEAT SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity
PEAT8002 - SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity Nick Rawlinson Research School of Earth Sciences Australian National University Anisotropy Introduction Most of the theoretical
More informationExercise: concepts from chapter 10
Reading:, Ch 10 1) The flow of magma with a viscosity as great as 10 10 Pa s, let alone that of rock with a viscosity of 10 20 Pa s, is difficult to comprehend because our common eperience is with s like
More informationThermal Analysis Contents - 1
Thermal Analysis Contents - 1 TABLE OF CONTENTS 1 THERMAL ANALYSIS 1.1 Introduction... 1-1 1.2 Mathematical Model Description... 1-3 1.2.1 Conventions and Definitions... 1-3 1.2.2 Conduction... 1-4 1.2.2.1
More informationSoil Damping Ratio: Theoretical Aspect and Measurement
Ratio: Theoretical Aspect and Measurement Sri Atmaja P. Rosyidi, Ph.D. Assistant Professor, Universitas Muhammadiyah Yogyakarta A Two Day Workshop on SASW for Practicing Engineer 17-18 February 2011, Faculty
More information3D simulations of an injection test done into an unsaturated porous and fractured limestone
3D simulations of an injection test done into an unsaturated porous and fractured limestone A. Thoraval *, Y. Guglielmi, F. Cappa INERIS, Ecole des Mines de Nancy, FRANCE *Corresponding author: Ecole des
More informationJuan E. Santos a,b,c, Gabriela B. Savioli a and Robiel Martínez Corredor c a
Juan E. Santos a,b,c, Gabriela B. Savioli a and Robiel Martínez Corredor c a Universidad de Buenos Aires, Fac. Ing., IGPUBA, ARGENTINA b Department of Mathematics, Purdue University, USA c Universidad
More informationTheoretical and Experimental Studies of Seismoelectric Conversions in Boreholes
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 3, No. 1, pp. 109-120 Commun. Comput. Phys. January 2008 Theoretical and Experimental Studies of Seismoelectric Conversions in Boreholes Zhenya Zhu 1,, Shihong
More informationDefmod, an earthquake simulator that adaptively switches between quasi-static and dynamic states
Defmod, an earthquake simulator that adaptively switches between quasi-static and dynamic states C. Meng (cmeng@mit.edu) B. Hager ERL, MIT May 16, 2016 Reservoir production/injection induced seismicity
More informationAn acoustic wave equation for orthorhombic anisotropy
Stanford Exploration Project, Report 98, August 1, 1998, pages 6?? An acoustic wave equation for orthorhombic anisotropy Tariq Alkhalifah 1 keywords: Anisotropy, finite difference, modeling ABSTRACT Using
More informationScienceDirect. Model-based assessment of seismic monitoring of CO 2 in a CCS project in Alberta, Canada, including a poroelastic approach
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 63 (2014 ) 4305 4312 GHGT-12 Model-based assessment of seismic monitoring of CO 2 in a CCS project in Alberta, Canada, including
More information7.2.1 Seismic waves. Waves in a mass- spring system
7..1 Seismic waves Waves in a mass- spring system Acoustic waves in a liquid or gas Seismic waves in a solid Surface waves Wavefronts, rays and geometrical attenuation Amplitude and energy Waves in a mass-
More informationRecovery-Based A Posteriori Error Estimation
Recovery-Based A Posteriori Error Estimation Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 2, 2011 Outline Introduction Diffusion Problems Higher Order Elements
More informationThe failure of the sounding assumption in electroseismic investigations
The failure of the sounding assumption in electroseismic investigations F.D. Fourie, J.F. Botha Institute for Groundwater Studies, University of the Free State, PO Box 339, Bloemfontein, 9300, South Africa
More informationOn propagation of Love waves in an infinite transversely isotropic poroelastic layer
Journal of Physics: Conference Series PAPER OPEN ACCESS On propagation of Love waves in an infinite transversely isotropic poroelastic layer To cite this article: C Nageswara Nath et al 2015 J. Phys.:
More informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
More informationWave Propagation in Fractured Poroelastic Media
Wave Propagation in Fractured Poroelastic Media WCCM, MS170: Advanced Computational Techniques in Geophysical Sciences, Barcelona, Spain, July 2014 Juan E. Santos Instituto del Gas y del Petróleo (IGPUBA),
More informationEffect of Anisotropic Permeability on Band Growth
Effect of Anisotropic Permeability on Band Growth Jesse Taylor-West October 7, 24 Introduction The mechanics of partially molten regions of the mantle are not well understood- in part due to the inaccessibility
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationDG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement
HONOM 2011 in Trento DG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement Institute of Aerodynamics and Flow Technology DLR Braunschweig 11. April 2011 1 / 35 Research
More informationMAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT
MAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT ABSTRACT A. G. Tarditi and J. V. Shebalin Advanced Space Propulsion Laboratory NASA Johnson Space Center Houston, TX
More informationNotes: Outline. Diffusive flux. Notes: Notes: Advection-diffusion
Outline This lecture Diffusion and advection-diffusion Riemann problem for advection Diagonalization of hyperbolic system, reduction to advection equations Characteristics and Riemann problem for acoustics
More informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,
More informationEffects of tool eccentricity on wave dispersion properties in borehole acoustic logging while drilling
Effects of tool eccentricity on wave dispersion properties in borehole acoustic logging while drilling Yibing Zheng and M. Nafi Toksöz Earth Resources Laboratory Dept. of Earth, Atmospheric, and Planetary
More informationFinite Volume Method for Scalar Advection in 2D
Chapter 9 Finite Volume Method for Scalar Advection in D 9. Introduction The purpose of this exercise is to code a program to integrate the scalar advection equation in two-dimensional flows. 9. The D
More informationPhysics 202 Laboratory 5. Linear Algebra 1. Laboratory 5. Physics 202 Laboratory
Physics 202 Laboratory 5 Linear Algebra Laboratory 5 Physics 202 Laboratory We close our whirlwind tour of numerical methods by advertising some elements of (numerical) linear algebra. There are three
More informationSplitting methods in the design of coupled flow and mechanics simulators
Splitting methods in the design of coupled flow and mechanics simulators Sílvia Barbeiro CMUC, Department of Mathematics, University of Coimbra PhD Program in Mathematics Coimbra, November 17, 2010 Sílvia
More information19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, 2-7 SEPTEMBER 2007
19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, 2-7 SEPTEMBER 2007 FREQUENCY DEPENDENCY AND ANISOTROPY OF THE ELASTIC CONSTANTS OF (NON-)POROUS MATERIALS AND THEIR INFLUENCE ON THE USAGE IN BUILDING
More informationTwo-dimensional flow in a porous medium with general anisotropy
Two-dimensional flow in a porous medium with general anisotropy P.A. Tyvand & A.R.F. Storhaug Norwegian University of Life Sciences 143 Ås Norway peder.tyvand@umb.no 1 Darcy s law for flow in an isotropic
More informationOn Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities
On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities Heiko Berninger, Ralf Kornhuber, and Oliver Sander FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Partial differential equations arise in a number of physical problems, such as fluid flow, heat transfer, solid mechanics and biological processes. These
More informationThe Lattice Boltzmann method for hyperbolic systems. Benjamin Graille. October 19, 2016
The Lattice Boltzmann method for hyperbolic systems Benjamin Graille October 19, 2016 Framework The Lattice Boltzmann method 1 Description of the lattice Boltzmann method Link with the kinetic theory Classical
More informationADAPTIVE MESH REFINEMENT USING WAVE-PROPAGATION ALGORITHMS FOR HYPERBOLIC SYSTEMS
SIAM J. NUMER. ANAL. c 998 Society for Industrial and Applied Mathematics Vol. 35, No. 6, pp. 2298 236, December 998 9 ADAPTIVE MESH REFINEMENT USING WAVE-PROPAGATION ALGORITHMS FOR HYPERBOLIC SYSTEMS
More informationSound Propagation through Media. Nachiketa Tiwari Indian Institute of Technology Kanpur
Sound Propagation through Media Nachiketa Tiwari Indian Institute of Technology Kanpur LECTURE-13 WAVE PROPAGATION IN SOLIDS Longitudinal Vibrations In Thin Plates Unlike 3-D solids, thin plates have surfaces
More informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationObservation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story
Observation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story Petr Kolinsky 1, Leo Eisner 1, Vladimir Grechka 2, Dana Jurick 3, Peter Duncan 1 Summary Shear
More informationTypical anisotropies introduced by geometry (not everything is spherically symmetric) temperature gradients magnetic fields electrical fields
Lecture 6: Polarimetry 1 Outline 1 Polarized Light in the Universe 2 Fundamentals of Polarized Light 3 Descriptions of Polarized Light Polarized Light in the Universe Polarization indicates anisotropy
More informationAcoustic landmine detection: a 3D poroelastic model
Acoustic landmine detection: a 3D poroelastic model Y. Q. Zeng a and Q. H. Liu a a Department of Electrical and Computer Engineering, Duke University, Durham, North Carolina 2778, U.S.A. ABSTRACT Acoustic
More informationA 3 D finite element model for sound transmission through a double plate system with isotropic elastic porous materials
Acoustics and Vibrations Group Université de Sherbrooke, QC CANADA Département génie mécanique Université de Sherbrooke Sherbrooke, QC CANADA Tel.: (819) 821-7157 Fax: (819) 821-7163 A 3 D finite element
More informationPartial Saturation Fluid Substitution with Partial Saturation
Fluid Substitution with 261 5 4.5 ρ fluid S w ρ w + S o ρ o + S g ρ g Vp (km/s) 4 3.5 K fluid S w K w + S o K o + S g K g Patchy Saturation Drainage 3 2.5 2 Fine-scale mixing 1 = S w + S o + S g K fluid
More informationEstimation of fractured weaknesses and attenuation factors from azimuthal seismic data
Estimation of fractured weaknesses and attenuation factors from azimuthal seismic data Huaizhen Chen and Kris Innanen CREWES project Department of Geoscience University of Calgary Summary Seismic wave
More informationTurbulence Modeling I!
Outline! Turbulence Modeling I! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions!
More informationCalculation of Sound Fields in Flowing Media Using CAPA and Diffpack
Calculation of Sound Fields in Flowing Media Using CAPA and Diffpack H. Landes 1, M. Kaltenbacher 2, W. Rathmann 3, F. Vogel 3 1 WisSoft, 2 Univ. Erlangen 3 inutech GmbH Outline Introduction Sound in Flowing
More informationCompensating visco-acoustic effects in anisotropic resverse-time migration Sang Suh, Kwangjin Yoon, James Cai, and Bin Wang, TGS
Compensating visco-acoustic effects in anisotropic resverse-time migration Sang Suh, Kwangjin Yoon, James Cai, and Bin Wang, TGS SUMMARY Anelastic properties of the earth cause frequency dependent energy
More informationMechanics of Materials and Structures
Journal of Mechanics of Materials and Structures ON TORSIONAL VIBRATIONS OF INFINITE HOLLOW POROELASTIC CYLINDERS M. Tajuddin and S. Ahmed Shah Volume 2, Nº 1 January 27 mathematical sciences publishers
More informationAnisotropy of Shale Properties: A Multi-Scale and Multi-Physics Characterization
Observation Scale Wavelength 10 0 10 4 10 6 10 8 10 12 10 16 10 18 10 20 Frequency (Hz) Anisotropy of Shale Properties: A Multi-Scale and Multi-Physics Characterization Elastic, Mechanical, Petrophysical
More informationD scattering of obliquely incident Rayleigh waves by a saturated alluvial valley in a layered half-space
1842. 3-D scattering of obliquely incident Rayleigh waves by a saturated alluvial valley in a layered half-space Zhenning Ba 1, Jianwen Liang 2 Department of Civil Engineering, Tianjin University, Tianjin
More informationEdinburgh Anisotropy Project, British Geological Survey, Murchison House, West Mains
Frequency-dependent AVO attribute: theory and example Xiaoyang Wu, 1* Mark Chapman 1,2 and Xiang-Yang Li 1 1 Edinburgh Anisotropy Project, British Geological Survey, Murchison House, West Mains Road, Edinburgh
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationHigh-resolution finite volume methods for extracorporeal shock wave therapy
High-resolution finite volume methods for extracorporeal shock wave therapy Kirsten Fagnan 1, Randall J. LeVeque 1, Thomas J. Matula 2 and Brian MacConaghy 2 1 Department of Applied Mathematics, University
More informationConservation Laws and Finite Volume Methods
Conservation Laws and Finite Volume Methods AMath 574 Winter Quarter, 2017 Randall J. LeVeque Applied Mathematics University of Washington January 4, 2017 http://faculty.washington.edu/rjl/classes/am574w2017
More information16 Rainfall on a Slope
Rainfall on a Slope 16-1 16 Rainfall on a Slope 16.1 Problem Statement In this example, the stability of a generic slope is analyzed for two successive rainfall events of increasing intensity and decreasing
More informationFinite Element Solver for Flux-Source Equations
Finite Element Solver for Flux-Source Equations Weston B. Lowrie A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Aeronautics Astronautics University
More informationPlasma Processes. m v = ee. (2)
Plasma Processes In the preceding few lectures, we ve focused on specific microphysical processes. In doing so, we have ignored the effect of other matter. In fact, we ve implicitly or explicitly assumed
More informationTable 17 1 Some general field equation terms. Heat Power. Current Source. 0 0 Boundary Current Porous Media Flow. Flow Source
17 Related Analogies 17.1 Basic Concepts The differential equation used in a finite element study in one discipline often appears in a different discipline, but with a different physical meaning for the
More informationA Multi-Continuum Multi-Component Model for Simultaneous Enhanced Gas Recovery and CO 2 Storage in Stimulated Fractured Shale Gas Reservoirs Jiamin
A Multi-Continuum Multi-Component Model for Simultaneous Enhanced Gas Recovery and CO 2 Storage in Stimulated Fractured Shale Gas Reservoirs Jiamin Jiang M.S. Candidate Joined Fall 2013 1 Main Points Advanced
More informationSound Pressure Generated by a Bubble
Sound Pressure Generated by a Bubble Adrian Secord Dept. of Computer Science University of British Columbia ajsecord@cs.ubc.ca October 22, 2001 This report summarises the analytical expression for the
More information