Wave Propagation in Poroelastic Materials
|
|
- Jody Ross
- 5 years ago
- Views:
Transcription
1 Wave Propagation in Poroelastic Materials M. Yvonne Ou Department of Mathematical Sciences University of Delaware Newark, DE ructure Interaction Probs, May Acknowledgement: Research funded by NSF-DMS Math. Biology Grant and NIH CBER PILOT Grant , and collaborators E. Cherkaev (UUtah), R. Lemoine & R. LeVeque (UW), S. Nintcheu(ORNL), M. Mak & M. Xu (Chinese U. of HK)
2 Outline Motivation Biot s equation for poroelasticity Balance form and finite volume method Numerical results of a plane strain problem using CLAWPACK Application to cancellous bones in literature Numerical simulation for effective elastic moduli by boundary element method (BEM) Conclusion and Future Work
3 Images of cancellous bone and cortical bone Figure:
4 Main interest: Two classes of inverse problems for composite materials To infer statistical information from measurement of effective properties: Dehomogenization To infer effective properties from waves requires building an effective wave solver for the governing equation and a better inversion scheme [Buchanan-Gilbert-MYO-211].
5 Effective permittivity with isotropic constituents [Golden-Papanicolaou-83] The complex permittivity of the medium is modeled by a (spatially) stationary random field (x ; ), x 2 R d and 2, where is the set of all realizations of the random medium, (x ; ) = 1 1 (x ; ) (x ; ): (1) E (x ; ): electric field, D(x ; ): displacement field e k : unit vector in direction k, < >: ensemble average over (PDE) r D = ; r E = ; < E >= e k : (2) Effective complex permittivity < E > def = < D >
6 Integral representation formula (IRF) Theorem (Golden-Papanicolaou-83) Let be one of the diagonal component of the effective permittivity tensor and s def = 1=(1 h); h def = 1 = 2. Then F (s) def = 1 (h)= 2 is analytic off [; 1] in the s plane with the integral representation F (s) = 1 Z 1 = 2 dz s z ; (3) where the positive measure on [; 1] is the spectral measure of the self-adjoint operator 1, where = r( ) 1 r (Helmholtz projector). Separation of contrast (h) and information of microstructure (). This provides a characterization of microstructure (Cherkaev and Zhang). If the is known, then computation of the effective property amounts to only plug in the contrast s into IRF.
7 Numerical results for measure Uniqueness of from data on an arc on the!-plane [E. Cherkaev 21] p=.11 p=.25 p= Why Moments?.25 p=.43 p=.25 p= Figure: Reconstructions of the spectral function of Maxwell-Garnett composite using quadratic regularization (left figure) and nonnegativity constraint (right figure). [Cherkaev-MYO-8]
8 First few moments Table: n calculated for Maxwell-Garnett composite solving regularized problem for the spectral function. Volume fractions: 11%; 25%; 43%. [Cherkaev-MYO-8] Exact N.neg Tik Exact N.neg Tik Exact N.neg Tik Note: Tik: Tikhonov; n-neg: Non-negativity
9 Relations between the moments of and the microstructure [Golden-Papanicolaou-83] m(h) := = 2 ; h := 1 = 2, analytic outside ( 1; ]. Expanding m at h = 1 and F (s) at s = 1 Z ( 1) n 1 1 m (n) (1) = z n 1 (dz ) =: (n 1) : n! Differentiating m at h = 1 m (1) = volume fraction of region occupied by 1 : (n 1) is related to the integral of n-point correlation functions, n = 1; 2; 3; : : :. If the permittivity 1 or 2 depends on parameters such as frequency or temperature, we may reconstruct the moments from information of 1, 2 and.
10 Reconstruction of Moments [Cherkaev-MYO-8, MYO-1] F (s) = Z 1 dz s z = s + 1 s s 3 + :::; n = Z 1 z n dz Reconstruction of coefficients of analytic functions in Hardy space H 2 in the unit circle jj < 1, = 1=s. 1 Polynomial interpolation of the kernel function s k z using n nodes z i, i ; k = 1; : : : ; n on [; 1] this results in an inversion scheme for the weights of a finite sum of Dirac measures supported at z i whose moments are close to those of. [L=M ] Padé approximation of the effective function F (s) this results in an inversion scheme for the coefficients of Padé approximants whose first L + M + 1 Taylor coefficients agree with those of F (s) at s = 1. Analytic continuation of F (s) into jsj < 1 with simple poles on [; 1].
11 Biot s equations for waves in poroelastic materials M. A. Biot, J. Acous. Socie. Amer. (1956a, 1956b). Three pieces of physics 1 Strain-stress relation 2 Kinetic energy 3 Energy dissipation u: displacement vector for solid U : displacement vector for fluid * Stress acting on the solid part of each face: ij, i ; j = 1; 2; 3 * Stress tensor acting on the fluid part of each face: s ij, s = p, : porosity. * Strain tensor in the solid: e ij = 1 2 j * Strain tensor in the fluid: =r U i ), i ; j = 1; 2; 3 Coupled stress-strain relation ij = 2N e ij + ij (Ae + Q), s = Q + Re. Here e def = r u.
12 Kinetic energy T = j _uj ( _u _U ) j U _ j 2 Dissipation function D = bj _u _U j 2 = b[( u_ 1 U_ 1 ) 2 + ( u_ 2 U_ 2 ) 2 + ( u_ 3 U_ 3 ) 2 ] *^b is frequency-dependent for high frequency! >! c Balance P of force in x 1 -direction Solid j = j U _ U _ _ u 1 Dynamic equations N r 2 u + r[(a + N )e + 2 ( 11 u + 12 U ) (u U ) r[qe + 2 ( 12 u + 22 U (u U ) Two longitudinal waves by taking r on both sides to get wave equations r 2 (Pe + 2 ( 11 e + 12 ) (e ) r 2 (Qe + 2 ( 12 e (e ) *P=A+2N
13 Biot s 1962 formulation in JAP Stress-strain relations = C u ɛ = [ 11 ; 22 ; 33 ; 23 ; 13 ; 12 ; p] t ɛ = [ 11 ; 22 ; 33 ; 2 23 ; 2 13 ; 2 12 ; ] t = div [(U u)] : variation of Fluid content Equations of j ( ij ) t v i + t q i + f i ; i = 1; 2; 3 v solid velocity, q = (@ t U v): fluid velocity rel. to solid, = (1 ) s + f Dynamic Darcy s i p = t v i + i t q i (t) i (t): viscodynamic operator
14 J. Carcione, JASA 1996! c := min i m, m i i := f T i i For! <! c : i (t) = m i (t) + H (t) i p = t v i + m t q i + i q i For!! c, approximate by the generalized Zener kernel 2 XL i 3 i (t) = 1 e t il 5 H i p = t v i + i L i l=1 il il il ; il : relaxation times; il il 1 i q i i L i XL i 1 il A qi + il t e 1 il = e 1 il + il L i il 1 XL i i e il l=1 il q i il
15 A transversely isotropic case: plane strain problem Stress-strain t xx = c u x v x + c u z v z + 1 M (@ x q x z q z ) t s 11 t zz = c u x v x + c u z v z + 3 M (@ x q x z q z ) t s 33 t xz = c u 55 (@ z v x x v z ) t s 55 t p = 1 x v x 3 z v z M (@ x q x z q z ) t s f (7) Equations of motion Darcy s law for low frequency t v x + t q x x xx z xz t v z + t q z x xz z zz x p = t v x + m t q x z p = t v z + m t q z + q x (1) 1 q z (11) 3
16 Balance Law where A t g + A@ x g + B@ z g = Dg t s (12) t g = xx zz xz v x v z p q x q z t s t s t s t s t s f c11 u 1 M c13 u 3 M c55 u m 1 41 f 41 m M M f f 43 1 C A (13)
17 B = c13 u 1 M c33 u 3 M c55 u m 1 41 m 3 43 f 43 3 M M f 41 f f 4 D = 11 f C A 1 C A (14) (15)
18 Finite Volume Method Normal and transverse Riemann problem solver Godunov operator splitting with 4t = O(1 5 ) MC limiter Preliminary results 1. Materials with different porosities. Left: :2, Right: :5 2. AMRCLAW zz and p excited at point source by Ricker pulse with peak frequency of 3135Hz, sandstone and shale
19 e228 A. Hosokawa / Ultrasonics 44 (26) e227 e231 Outline Dehomogenization Ultrasound Image processing Computation Parameters of bovine cancellous bone models in the viscoelastic FDTD Ultrasound in cancellous bone and method relevant effective 2.1. Viscoelastic FDTD method Trabecular thickness a parameters s.15 mm Pore space between trabeculae af 1.35 mm e s - r - - t e F ðjþ ¼ 1 jt ðjþ ; 2. FDTD methods The wave propagation in a 2D viscoelastic medium is Trabecular length ls 4.5 mm governed by the following equations: sonics Hosakawa 44 (26) e227 e231 26, Ultrasonics e229 q o_u i ot ¼ os ii oi þ os These values were uniformly distributed in the range of ±25%. xy oj ; ð1aþ Table 3 os ii Þ Biot s parameters of bovine cancellous ot þ c iis bone ii ¼ðkþ2lÞ o_u i models oi þ k o_u j oj with ; parallel ð1bþ and perpendicular trabecular orientations os xy Þ ot þ c Parallel s xy xy ¼ l o_u x oy þ o_u y ; ð1cþ ox Perpendicular Young s modulus of solid bone E s where i, 22. j = x, y GPa (i 5 j), u i are the 22. particle GPa displacement Þ Poisson s ratio of solid bone m s components.32 of the medium in the.23 i-direction (an overdot Bulk modulus of bone marrow K f denotes 2. a timegpa derivative), s 2. GPa Density of solid bone q s 196 kg/m 3 ii are the normal stress components in the i-direction, s 196 kg/m 3 Þ Density of bone marrow q f 93 kg/m 3 xy is the shear stress component on the x y plane, c 93 kg/m 3 ii and c xy are, respectively, the resistance Poisson s ratio of trabecular frame m b coefficients for the normal and shear stress components related to the attenuation coefficients of the longitudinal and Variable r Þ Viscosity of bone marrow g 1.5 P s 1.5 P s shear waves, q is the medium density and k, l are Lamé Parameter n coefficients. For Eqs. (1a) (1c), central differencing with Pore size a 1.35 mm 1.35 mm staggered grids is Permeability k used to constitute a leap-frog system [8] for _u Porosity b i and s.83 ii, s xy..83 Resistance coefficients c ii (i = x,y) A schematic 1. drawing 1 6 of the cancellous 1. bone 1 model 6 in the Resistance coefficient c viscoelastic xy 1. FDTD 1 method 7 [9] is shown 1. in Fig The skeletal frame of cancellous bone was made of numerous trabecular rods. The porosity b was adjusted by changing the trabecular thickness a s and pore space between trabeculae a f, and the trabecular orientation was adjusted by changing M. Y. Ou (UDel) the trabecular length Poroelasticity l s and angle h. In order to form an Table 1 Fig. 2. Geometric drawings of bovine cancellous bone models in the viscoelastic FDTD method, with (a) parallel and (b) perpendicular trabecular orientations. Table 2
20 to the trabecular alignment, and Fig. 4(b) shows the waveforms in the perpendicular direction. The experimental waveforms in the parallel and perpendicular directions Hosakawa s result Outline Fig. Dehomogenization 4 shows the viscoelastic FDTD (upper) Ultrasound and Biot s it is shown in Fig. Image 4(b) that, processing using Biot s FDTD method, Computation FDTD (middle) results of the pulsed waveforms propagating the simulated waveform in the perpendicular direction is through bovine cancellous bone, together with the not in good agreement with the experimental waveform. experimentally observed results (lower). Fig. 4(a) shows The slower speed of the fast wave and the higher attenuation the waveforms for the propagation in the direction parallel of the slow wave, which are the propagation characteristics in the perpendicular direction, could be analyzed by increasing value of the parameter n and by decreasing value of the permeability k because n and k are inversely propor- Amplitude (arb. units) Amplitude (arb. units) Time (2 μs/div.) (a) Time (2 μs/div.) (b) Fig. 4. Viscoelastic FDTD (upper), Biot s FDTD (middle) and experimental (lower) results for observations of pulsed waveforms propagating through bovine cancellous bone in the directions (a) parallel and (b) perpendicular to the trabecular alignment.
21 The fluid water characteristics are: bulk modulus K f 2.28 GPa, density f 1 kg m 3, viscosity 1 3 kgms 1. The incident experimental signal generated by the transducer is shown in Fig. 3 and its spectrum in Fig. 4. Fellah et. al. 24, 26,JASA Outline Dehomogenization Ultrasound Image processing Computation Ultrasound in cancellous bone -continue- TABLE I. Biot s model parameters of cancellous bone. Biot s model parameters of cancellous bone M1 M2 M3 Thickness L cm Density s (kg/m 3 ) Porosity Tortuosity Viscous characteristic length m Bulk modulus of the elastic solid Ks (GPa) Bulk modulus of bone skeletal frame Kb (GPa) Shear modulus of the frame N (GPa) FIG. 17. Comparison between experimental transmitted signal solid lin and simulated FIG. 15. transmitted Experimental signal setup dashed for ultrasonic line for bone measurements. sample M2. 7 J. Acoust. Soc. Am., Vol. 116, No. 1, July 24 fluid phases of Fellah theet porous al.: Ultrasonic material propagation and thus in cancellous on the oth bo parameters whichncts were kept constant Wkspduring on this Fluid-St study. In Fellah 26,, 1,, K b and N are sought after by solving the inverse problem.
22 VI. CONCLUSION Outline a fast wave, and adehomogenization part is transmitted as a slow wave. When Ultrasound Image processing Computation any of these components, traveling at different speeds, hit the In this paper, analytical calculus of the reflection and second surface, a similar effect takes place: part is transmitted into the fluid, part is reflected as a fast or slow wave. elastic frame is established. This calculus is based on Biot s transmission coefficient of a slab of cancellous bone with an Fellah et al s results Sample characteristics are measured using standard theory modified by the Johnson et al. model to describe the methods 11,16,17,19,25 27 Appendix B and given in Table I. viscous interaction between fluid and structure. Numerical Solid The line: fluid water Experimental characteristics transmitted are: bulk signal; modulus Dashed K line: Simulated transmitted signal f simulations of transmitted waves in the time domain were 2.28 GPa, density f 1 kg m 3, viscosity 1 3 kgms 1. The incident experimental signal generated by the transducer is shown in Fig. 3 and its spectrum in Fig. 4. TABLE I. Biot s model parameters of cancellous bone. Biot s model parameters of cancellous bone M1 M2 M3 Thickness L cm Density s (kg/m 3 ) Porosity Tortuosity Viscous characteristic length m nts. Bulk modulus of the elastic solid Ks (GPa) Bulk modulus of bone skeletal frame Kb the other(gpa) tudy. Shear modulus FIG. 16. of the Comparison frame N (GPa) between experimental 2.6 transmitted 1.7 signal.35 solid line FIG. 17. Comparison between experimental transmitted For longitudinal signal solid waves, line the relation A5 and simulated transmitted signal dashed line for bone sample M1. and simulated transmitted signal dashed line for bone sample M2. s x, P 2N 2 s f 7 J. Acoust. Figures Soc. 16, Am., 17, Vol. and 116, 18 No. show 1, Julythe 24comparison between Fellah et al.: Ultrasonic propagation in cancellous bone x 2 x 2 cal simuperiments signals dashed line for bone samples M1, M2, and M3, experimental transmitted signals solid line and simulated Using relations 12 14, metrics A respectively. The experimental transmitted wave forms are s x, Z c uency of traveling through the cancellous bone in the same direction a 558PR as the trabecular alignment x direction. For some situations, Z re amplihigh fre- Z fast and slow waves are superimposed, depending on the coupling between the two porous medium phases, so that this le in this is due to the modified Biot s parameters for cancellous bone, Z removed which as mentioned in Sec. IV, play an important role in the setup is arrival time of the two waves. were ma- M. Y. The Ou experimental (UDel) data and FIG. theoretical 18. Comparison prediction between Poroelasticity experimental are transmitted signal solid line where the coefficients Z 1 ( ) and Z 2 (
23 Translation between microstructure and effective parameters Numerical simulation from images of cancellous bone to some effective properties. Ongoing numerical projects focusing on cancellous bones. Step 1 Image segmentation and 3D reconstruction from -CT scans (a) raw -CT scan (b) digitalized data (c) 3-D structure Figure: -CT scans are provided by Prof. L. Cardoso from CUNY Med. Center
24 Numerical projects for effective parameters of cancellous bone, SCBLAP with Dr. S. Nincheu at ORNL Step 2 Meshing and solving (BVP) by BEM Advantage Meshing only the interface and outer surface, not the entire domain. (a) meshed trabeculae (b) electric potential Figure: Our preliminary results
25 Effective dielectric tensor ɛ (BVP) ( Div [(ɛ1 1 (x ) + ɛ 2 2 (x ))r] = ; = P d j =1 E j x j (16) with the constant vector E = e j, the j -th unit vector in Cartesian coordinates ((e j ) i = ij ) followed by evaluating the integral Z ɛ def 1 ik = [ 1 1 (x ) (x )]@ i dx (17) jj 1 Z = jj [ 1 1 n 1 ds n 2 ds ] (18) 2
26 Illustration of ideas To solve Laplace equation: 4u = in a simply connected bounded domain with boundary and constant boundary elements. Key: Green s identity ( R R u(x); x 2 G(x; y)t(y)d y H(x; y) n(y)u(y)d y = ; x 2 R 3 n R G(x; y) = 1 x y 4kx yk, H(x; y) = 4kx yk 3, 1 Triangulate into N non-overlapping elements j, j = 1; ; N. 2 Define g j (x) := G(x; y)d y and h j (x) := H(x; y) n(y)d y. j j 3 Discretized system is Gftg = Hfug, with u = (u 1 ; ; u N ) T and u = (t 1 ; ; t N ) T G ij = lim x!x j g j (x ), H ij = lim x!x j h j (x ), x j 2 j, j = 1; ; N. 4 Substitute in boundary data and rearrange to obtain R Afz g = fbg with fz g containing those u i and t j which are unknown, and fbg with the boundary data. 5 Solve the linear system to obtain u j and t j, j = 1; ; N. 6 For any x 2, we have u(x) = j=1 j=1 NX t j g j (x) NX u j h j (x)
27 Outline Dehomogenization 3. Real Trabecula Marrow Ultrasound Image processing Computation and A small piece of trabecula is used in this test case. The following are the results Figure (viii) Trabecula and Marrow Number of Iteration jj 11 jj 12 jj 13 Energy Norm Figure (ix) Cross-sectional View at z = 15 = Trabecula, Red = Marrow) (Blue
28 Future Work To generalize dehomogenization schemes for viscoelastic composites by using the 2-parameter IRF [MYO-211] To develop a better image processing scheme in order to handle bone sample of realistic size. To derive integral representation formula for permeability tensor. To compare the simulated results with lab measurement. To build highly parallel numerical codes for effective parameters computation for bones Numerical analysis of the FVM and BEM codes.
29 Thank you.
Poroelasticity simulation in CLAWPACK for geophysical and biological applications
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
More information13th International Symposium on Nondestructive Characterization of Materials (NDCM-XIII), May 2013, Le Mans, France
3th International Symposium on Nondestructive Characterization of Materials (NDCM-XIII), 2-24 May 23, Le Mans, France www.ndt.net/?id=5532 Biot waves in porous ceramic plates : influence of boundary conditions
More informationDETERMINATION OF THE PARAMETERS OF CANCELLOUS BONE USING LOW FREQUENCY ACOUSTIC MEASUREMENTS
Journal of Computational Acoustics, Vol. 12, No. 2 (2004) 99 126 c IMACS DETERMINATION OF THE PARAMETERS OF CANCELLOUS BONE USING LOW FREQUENCY ACOUSTIC MEASUREMENTS JAMES L. BUCHANAN Mathematics Department,
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 informationMacroscopic theory Rock as 'elastic continuum'
Elasticity and Seismic Waves Macroscopic theory Rock as 'elastic continuum' Elastic body is deformed in response to stress Two types of deformation: Change in volume and shape Equations of motion Wave
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 informationElements of Rock Mechanics
Elements of Rock Mechanics Stress and strain Creep Constitutive equation Hooke's law Empirical relations Effects of porosity and fluids Anelasticity and viscoelasticity Reading: Shearer, 3 Stress Consider
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationNDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16.
CAVITY INSPECTION NDT&E Methods: UT VJ Technologies NDT&E Methods: UT 6. NDT&E: Introduction to Methods 6.1. Ultrasonic Testing: Basics of Elasto-Dynamics 6.2. Principles of Measurement 6.3. The Pulse-Echo
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
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 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 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 informationProceedings, International Snow Science Workshop, Banff, 2014
SNOW SLAB FAILURE DUE TO BIOT-TYPE ACOUSTIC WAVE PROPAGATION Rolf Sidler 2 Simon Fraser University, 8888 University Drive, Burnaby V5A1S6, Canada ABSTRACT: Even though seismic methods are among the most
More informationAn alternative Biot s displacement formulation for porous materials
An alternative Biot s displacement formulation for porous materials Olivier Dazel, a Bruno Brouard, Claude Depollier, and Stéphane Griffiths Laboratoire d Acoustique de l Université du Maine - UMR CNRS
More informationAcoustic Detection of Buried Objects in 3-D Fluid Saturated Porous Media: Numerical Modeling
IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING, VOL. 39, NO. 6, JUNE 2001 1165 Acoustic Detection of Buried Objects in 3-D Fluid Saturated Porous Media: Numerical Modeling Yan Qing Zeng, Student Member,
More informationHigh-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 informationPhysics of Continuous media
Physics of Continuous media Sourendu Gupta TIFR, Mumbai, India Classical Mechanics 2012 October 26, 2012 Deformations of continuous media If a body is deformed, we say that the point which originally had
More informationFig. 1. Circular fiber and interphase between the fiber and the matrix.
Finite element unit cell model based on ABAQUS for fiber reinforced composites Tian Tang Composites Manufacturing & Simulation Center, Purdue University West Lafayette, IN 47906 1. Problem Statement In
More informationPhysical and Biological Properties of Agricultural Products Acoustic, Electrical and Optical Properties and Biochemical Property
Physical and Biological Properties of Agricultural Products Acoustic, Electrical and Optical Properties and Biochemical Property 1. Acoustic and Vibrational Properties 1.1 Acoustics and Vibration Engineering
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 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 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 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 informationAcoustic wave reflection from the transition layer of surficial marine sediment
Acoust. Sci. & Tech. 25, 3 (2004) PAPER Acoustic wave reflection from the transition layer of surficial marine sediment Masao Kimura and Takuya Tsurumi School of Marine Science and Technology, Tokai University
More informationModeling Scattering from Rough Poroelastic Surfaces Using COMSOL Multiphysics
Modeling Scattering from Rough Poroelastic Surfaces Using COMSOL Multiphysics Anthony L. Bonomo *1 Marcia J. Isakson 1 and Nicholas P. Chotiros 1 1 Applied Research Laboratories The University of Texas
More informationInterface conditions for Biot s equations of poroelasticity Boris Gurevich The Geophysical Institute of Israel, P.O. Box 2286, Holon 58122, Israel
Interface conditions for Biot s equations of poroelasticity Boris Gurevich The Geophysical Institute of Israel, P.O. Box 2286, Holon 58122, Israel Michael Schoenberg Schlumberger-Doll Research, Old Quarry
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 informationOn the study of elastic wave scattering and Rayleigh wave velocity measurement of concrete with steel bar
NDT&E International 33 (2000) 401 407 www.elsevier.com/locate/ndteint On the study of elastic wave scattering and Rayleigh wave velocity measurement of concrete with steel bar T.-T. Wu*, J.-H. Sun, J.-H.
More information6th NDT in Progress Lamb waves in an anisotropic plate of a single crystal silicon wafer
6th NDT in Progress 2011 International Workshop of NDT Experts, Prague, 10-12 Oct 2011 Lamb waves in an anisotropic plate of a single crystal silicon wafer Young-Kyu PARK 1, Young H. KIM 1 1 Applied Acoustics
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 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 informationPorous Materials for Sound Absorption and Transmission Control
Purdue e-pubs Publications of the Ray W. School of Mechanical Engineering 8-2005 Porous Materials for Sound Absorption and Transmission Control J Stuart Bolton, bolton@purdue.edu Follow this and additional
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 informationSource Free Surface x
Finite-dierence time-domain model for elastic waves in the ground Christoph T. Schroeder and Waymond R. Scott, Jr. School of Electrical and Computer Engineering Georgia Institute of Technology Atlanta,
More informationMechanics of Biomaterials
Mechanics of Biomaterials Lecture 7 Presented by Andrian Sue AMME498/998 Semester, 206 The University of Sydney Slide Mechanics Models The University of Sydney Slide 2 Last Week Using motion to find forces
More informationAnalytical Solution for Acoustic/Poroelastic Interaction
Analytical Solution for Acoustic/Poroelastic Interaction Julien Diaz,2 Abdelâaziz Ezziani 2, INRIA Bordeaux Sud-Ouest, EPI Magique 3D 2 LMA, UMR 542, Université de Pau 7 April 2009 Acoustic/Poroelastic
More informationMecanum. Acoustic Materials: Characterization. We build silence. Mecanum Inc.
ecanum We build silence Acoustic aterials: Characterization ecanum Inc. info@mecanum.com www.mecanum.com otivation Sound quality in vehicles starts at the design stage odels are used to simulate the acoustics
More informationContinuum Mechanics. Continuum Mechanics and Constitutive Equations
Continuum Mechanics Continuum Mechanics and Constitutive Equations Continuum mechanics pertains to the description of mechanical behavior of materials under the assumption that the material is a uniform
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationHandbook of Radiation and Scattering of Waves:
Handbook of Radiation and Scattering of Waves: Acoustic Waves in Fluids Elastic Waves in Solids Electromagnetic Waves Adrianus T. de Hoop Professor of Electromagnetic Theory and Applied Mathematics Delft
More information1 f. result from periodic disturbance same period (frequency) as source Longitudinal or Transverse Waves Characterized by
result from periodic disturbance same period (frequency) as source Longitudinal or Transverse Waves Characterized by amplitude (how far do the bits move from their equilibrium positions? Amplitude of MEDIUM)
More informationAbstract. 1 Introduction
Efficiency of absorbing boundary conditions forfluid-saturatedporous media T. Akiyoshi, K. Fuchida, H.L. Fang Department of Civil and Environmental Engineering, Kumamoto University, 2-39-1 Kurokami, Kumamoto
More informationCharacterization of discretely graded materials using acoustic wave propagation
Computational Materials Science 37 (006) 0 8 www.elsevier.com/locate/commatsci Characterization of discretely graded materials using acoustic wave propagation Ritesh Samadhiya a, Abhijit Mukherjee a, *,
More information3D Elasticity Theory
3D lasticity Theory Many structural analysis problems are analysed using the theory of elasticity in which Hooke s law is used to enforce proportionality between stress and strain at any deformation level.
More informationLaboratory experiments and numerical simulation on Bitumen Saturated Carbonates: A Rock Physics Study for 4D Seismology
Laboratory experiments and numerical simulation on Bitumen Saturated Carbonates: A Rock Physics Study for 4D Seismology Arif Rabbani Jason Nycz* Zizhen Wong Doug Schmitt Ken Gray Department of Physics
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 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 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 informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Standard Solids and Fracture Fluids: Mechanical, Chemical Effects Effective Stress Dilatancy Hardening and Stability Mead, 1925
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 informationFractional acoustic wave equations from mechanical and. Sverre Holm
Fractional acoustic wave equations from mechanical and thermal constitutive equations Sverre Holm Outline 2 Perspectives Non fractional wave equations Low and high frequency regions Fractional viscoelastic
More informationSound Propagation in Porous Media
Final Project Report for ENGN34 Sound Propagation in Porous Media ---Numerical simulation based on MATLAB Name: Siyuan Song Department: Engineering Date: Dec.15 17 1 Name: Siyuan Song Department: Engineering
More informationEfficient outdoor sound propagation modelling in time-domain
Efficient outdoor sound propagation modelling in time-domain Timothy Van Renterghem Acoustics Group Department of Information Technology Ghent University/iMinds CEAS, X-Noise : Atm. and ground effects
More informationBorehole Geophysics. Acoustic logging measurements
Acoustic logging measurements - Review of basic physics background - Concept of P- and S-wave measurements and logging tools - Tube waves - Seismic imaging - Synthetic seismograms - Field application examples
More informationElements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004
Elements of Continuum Elasticity David M. Parks Mechanics and Materials II 2.002 February 25, 2004 Solid Mechanics in 3 Dimensions: stress/equilibrium, strain/displacement, and intro to linear elastic
More informationVelocity Measurements of Pore Fluids at Pressure and Temperature: Application to bitumen
Velocity Measurements of Pore Fluids at Pressure and Temperature: Application to bitumen Arif Rabbani 1*, Douglas R Schmitt 1, Jason Nycz 2, and Ken Gray 3 1 Institute for Geophysical Research, Department
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 informationOne-Dimensional Numerical Solution of the Maxwell-Minkowski Equations
Tamkang Journal of Science and Engineering, Vol. 12, No. 2, pp. 161168 (2009) 161 One-Dimensional Numerical Solution of the Maxwell-Minkowski Equations Mingtsu Ho 1 and Yao-Han Chen 2 1 Department of Electronic
More informationReceiver. Johana Brokešová Charles University in Prague
Propagation of seismic waves - theoretical background Receiver Johana Brokešová Charles University in Prague Seismic waves = waves in elastic continuum a model of the medium through which the waves propagate
More informationSummary PHY101 ( 2 ) T / Hanadi Al Harbi
الكمية Physical Quantity القانون Low التعريف Definition الوحدة SI Unit Linear Momentum P = mθ be equal to the mass of an object times its velocity. Kg. m/s vector quantity Stress F \ A the external force
More informationarxiv: v1 [physics.geo-ph] 4 Feb 2015
1 arxiv:152.1284v1 [physics.geo-ph] 4 Feb 215 Acoustic waves in a Biot-type porous snow model: The fast slow wave in light snow Rolf Sidler Department of Earth Sciences, Simon Fraser University, 8888 University
More informationFinite Element Modeling of Ultrasonic Transducers for Polymer Characterization
Excerpt from the Proceedings of the COMSOL Conference 2009 Milan Finite Element Modeling of Ultrasonic Transducers for Polymer Characterization Serena De Paolis *, Francesca Lionetto and Alfonso Maffezzoli
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 informationSENSITIVITY ANALYSIS OF THE PETROPHYSICAL PROPERTIES VARIATIONS ON THE SEISMIC RESPONSE OF A CO2 STORAGE SITE. Juan E. Santos
SENSITIVITY ANALYSIS OF THE PETROPHYSICAL PROPERTIES VARIATIONS ON THE SEISMIC RESPONSE OF A CO2 STORAGE SITE Juan E. Santos Instituto del Gas y del Petróleo, Facultad de Ingeniería UBA and Department
More informationSaturation Effects of Soils on Ground Motion at Free Surface Due to Incident SV Waves
Saturation Effects of Soils on Ground Motion at Free Surface Due to Incident SV Waves Jun Yang, M.ASCE 1 Abstract: A study is presented of saturation effects of subsoil on seismic motions at the free surface
More informationElasticity: Term Paper. Danielle Harper. University of Central Florida
Elasticity: Term Paper Danielle Harper University of Central Florida I. Abstract This research was conducted in order to experimentally test certain components of the theory of elasticity. The theory was
More informationOPAC102. The Acoustic Wave Equation
OPAC102 The Acoustic Wave Equation Acoustic waves in fluid Acoustic waves constitute one kind of pressure fluctuation that can exist in a compressible fluid. The restoring forces responsible for propagating
More informationNumerical simulation of the Biot slow wave in water-saturated Nivelsteiner Sandstone
GEOPHYSICS, VOL. 66, NO. 3 (MAY-JUNE 2001; P. 890 896, 4 FIGS., 2 TABLES. Short Note Numerical simulation of the Biot slow wave in water-saturated Nivelsteiner Sstone Børge Arntsen José M. Carcione INTRODUCTION
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 informationIDENTIFICATION OF VISCOELASTIC MODULI OF COMPOSrrn MATERIALS FROM TIlE PLATE TRANSMISSION COEFFICIENTS
IDENTIFICATION OF VISCOELASTIC MODULI OF COMPOSrrn MATERIALS FROM TIlE PLATE TRANSMISSION COEFFICIENTS Bernard Hosten and Michel Castaings Laboratoire de Mecanique Physique, University of Bordeaux I, URA
More informationLECTURE 5 - Wave Equation Hrvoje Tkalčić " 2 # & 2 #
LECTURE 5 - Wave Equation Hrvoje Tkalčić " 2 # "t = ( $ + 2µ ) & 2 # 2 % " 2 (& ' u r ) = µ "t 2 % & 2 (& ' u r ) *** N.B. The material presented in these lectures is from the principal textbooks, other
More information20. Rheology & Linear Elasticity
I Main Topics A Rheology: Macroscopic deformation behavior B Linear elasticity for homogeneous isotropic materials 10/29/18 GG303 1 Viscous (fluid) Behavior http://manoa.hawaii.edu/graduate/content/slide-lava
More informationModelling Anisotropic, Hyperelastic Materials in ABAQUS
Modelling Anisotropic, Hyperelastic Materials in ABAQUS Salvatore Federico and Walter Herzog Human Performance Laboratory, Faculty of Kinesiology, The University of Calgary 2500 University Drive NW, Calgary,
More informationFlow and Transport. c(s, t)s ds,
Flow and Transport 1. The Transport Equation We shall describe the transport of a dissolved chemical by water that is traveling with uniform velocity ν through a long thin tube G with uniform cross section
More informationChapter 16 Waves in One Dimension
Chapter 16 Waves in One Dimension Slide 16-1 Reading Quiz 16.05 f = c Slide 16-2 Reading Quiz 16.06 Slide 16-3 Reading Quiz 16.07 Heavier portion looks like a fixed end, pulse is inverted on reflection.
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationWenyong Pan and Lianjie Huang. Los Alamos National Laboratory, Geophysics Group, MS D452, Los Alamos, NM 87545, USA
PROCEEDINGS, 44th Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 019 SGP-TR-14 Adaptive Viscoelastic-Waveform Inversion Using the Local Wavelet
More informationON THE MODELING AND SIMULATION OF AN ACOUSTIC CLOAK
5 th International Conference Advanced Composite Materials Engineering COMAT 014 16-17 October 014, Braşov, Romania ON THE MODELING AND SIMULATION OF AN ACOUSTIC CLOAK Veturia Chiroiu 1, Rodica Ioan 1,,
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 informationA unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation
A unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation S. Bordère a and J.-P. Caltagirone b a. CNRS, Univ. Bordeaux, ICMCB,
More informationValidity of the limp model for porous materials: A criterion based on the Biot theory
Validity of the limp model for porous materials: A criterion based on the Biot theory Olivier Doutres, a Nicolas Dauchez, Jean-Michel Génevaux, and Olivier Dazel Laboratoire d Acoustique UMR CNRS 6613,
More informationToday s menu. Last lecture. Ultrasonic measurement systems. What is Ultrasound (cont d...)? What is ultrasound?
Last lecture Measurement of volume flow rate Differential pressure flowmeters Mechanical flowmeters Vortex flowmeters Measurement of mass flow Measurement of tricky flows" Today s menu Ultrasonic measurement
More informationABHELSINKI UNIVERSITY OF TECHNOLOGY
ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for the Morley plate element Jarkko Niiranen Department of Structural
More informationA new second order absorbing boundary layer formulation for anisotropic-elastic wavefeld simulation
A new second order absorbing boundary layer formulation for anisotropic-elastic wavefeld simulation Junxiao Li, Kris Innanen and Bing Wang University of Calgary China University of Petroleum-Beijing Summary
More informationSemeniuk, B., Göransson, P. (2018) Modelling the Dynamic Viscous and Thermal Dissipation Mechanisms in a Fibrous Porous Material In:
http://www.diva-portal.org This is the published version of a paper presented at COMSOL Conference 2018 Lausanne. Citation for the original published paper: Semeniuk, B., Göransson, P. (2018) Modelling
More informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More informationNon-Newtonian Fluids and Finite Elements
Non-Newtonian Fluids and Finite Elements Janice Giudice Oxford University Computing Laboratory Keble College Talk Outline Motivating Industrial Process Multiple Extrusion of Pastes Governing Equations
More informationThe acoustic characterization of porous media and its standards
The acoustic characterization of porous media and its standards Luc JAOUEN 1, François-Xavier BECOT, Fabien CHEVILLOTTE Matelys, France ABSTRACT While there is a growing number of methods for the acoustic
More informationDevelopment of an analytical solution of modified Biot s equations for the optimization of lightweight acoustic protection
Development of an analytical solution of modified Biot s equations for the optimization of lightweight acoustic protection Jamil Kanfoud a and Mohamed Ali Hamdi Laboratoire Roberval, Université de Technologie
More informationRESEARCH PROPOSAL. Effects of scales and extracting methods on quantifying quality factor Q. Yi Shen
RESEARCH PROPOSAL Effects of scales and extracting methods on quantifying quality factor Q Yi Shen 2:30 P.M., Wednesday, November 28th, 2012 Shen 2 Ph.D. Proposal ABSTRACT The attenuation values obtained
More informationIntroduction to Seismology Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 12.510 Introduction to Seismology Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Stress and Strain
More informationPhysics in Faculty of
Why we study Physics in Faculty of Engineering? Dimensional analysis Scalars and vector analysis Rotational of a rigid body about a fixed axis Rotational kinematics 1. Dimensional analysis The ward dimension
More informationInversion of seismic AVA data for porosity and saturation
Inversion of seismic AVA data for porosity and saturation Brikt Apeland Thesis for the degree Master of Science Department of Earth Science University of Bergen 27th June 213 2 Abstract Rock physics serves
More informationSeismic behaviour of CO 2 saturated Fontainebleau sandstone under in situ conditions
Seismic behaviour of CO 2 urated Fontainebleau sandstone under in situ conditions Md Mizanul Huq Chowdhury*, University of Alberta, Edmonton, AB, Canada, mhchowdh@ualberta.ca and Douglas R. Schmitt, University
More informationLINK BETWEEN ATTENUATION AND VELOCITY DISPERSION
LINK BETWEEN ATTENUATION AND VELOCITY DISPERSION Jack Dvorkin Stanford University and Rock Solid Images April 25, 2005 SUMMARY In a viscoelastic sample, the causality principle links the attenuation of
More informationComputational Astrophysics
16 th Chris Engelbrecht Summer School, January 2005 3: 1 Computational Astrophysics Lecture 3: Magnetic fields Paul Ricker University of Illinois at Urbana-Champaign National Center for Supercomputing
More informationChapter 1 Fluid Characteristics
Chapter 1 Fluid Characteristics 1.1 Introduction 1.1.1 Phases Solid increasing increasing spacing and intermolecular liquid latitude of cohesive Fluid gas (vapor) molecular force plasma motion 1.1.2 Fluidity
More informationLecture 20: Modes in a crystal and continuum
Physics 16a, Caltech 11 December, 218 The material in this lecture (which we had no time for) will NOT be on the Final exam. Lecture 2: Modes in a crystal and continuum The vibrational modes of the periodic
More informationESTIMATION OF LAYER THICKNESS BY THE COST FUNCTION OPTIMIZATION: PHANTOM STUDY
ESTIMATION OF LAYER THICKNESS BY THE COST FUNCTION OPTIMIZATION: PHANTOM STUDY JURIJ TASINKIEWICZ, JERZY PODHAJECKI, JANUSZ WOJCIK, KATARZYNA FALIŃSKA, JERZY LITNIEWSKI Ultrasound Department Institute
More information