Pressure derivatives of bulk modulus for materials at extreme compression

Similar documents
Canadian Journal of Physics. Higher Order derivatives of bulk modulus of materials at infinite pressure

NEW RELATION BETWEEN VOLUME THERMAL EXPANSION AND THERMAL PRESSURE FOR SOLIDS

Pressure Volume Temperature (P-V-T) Relationships and Thermo elastic Properties of Geophysical Minerals

EOS-FIT V6.0 R.J. ANGEL

Equations of State. Tiziana Boffa Ballaran

Volume Thermal Expansion for Ionic Solids at High Temperatures

! " k x 2k$1 # $ k x 2k. " # p $ 1! px! p " p 1 # !"#$%&'"()'*"+$",&-('./&-/. !"#$%&'()"*#%+!'",' -./#")'.,&'+.0#.1)2,'!%)2%! !"#$%&'"%(")*$+&#,*$,#

Evaluation of pressure and bulk modulus for alkali halides under high pressure and temperature using different EOS

APPENDIX 1 ERROR ESTIMATION

Fluctuationlessness Theorem and its Application to Boundary Value Problems of ODEs

Micro-canonical ensemble model of particles obeying Bose-Einstein and Fermi-Dirac statistics

Curve Sketching. The process of curve sketching can be performed in the following steps:

Extreme Values of Functions

Differentiation. The main problem of differential calculus deals with finding the slope of the tangent line at a point on a curve.

Temperature dependent study of volume and thermal expansivity of solids based on equation of state

Two Universal Equations of State for Solids

Theorem Let J and f be as in the previous theorem. Then for any w 0 Int(J), f(z) (z w 0 ) n+1

Roberto s Notes on Differential Calculus Chapter 8: Graphical analysis Section 1. Extreme points

Math Review and Lessons in Calculus

2. ETA EVALUATIONS USING WEBER FUNCTIONS. Introduction

Exponential and. Logarithmic Functions. Exponential Functions. Logarithmic Functions

Numerical Methods - Lecture 2. Numerical Methods. Lecture 2. Analysis of errors in numerical methods

Non-newtonian Rabinowitsch Fluid Effects on the Lubrication Performances of Sine Film Thrust Bearings

Numerical Solution of Ordinary Differential Equations in Fluctuationlessness Theorem Perspective

Grüneisen parameters and isothermal equations of state

Analysis of volume expansion data for periclase, lime, corundum and spinel at high temperatures

Comments on Magnetohydrodynamic Unsteady Flow of A Non- Newtonian Fluid Through A Porous Medium

CONVECTIVE HEAT TRANSFER CHARACTERISTICS OF NANOFLUIDS. Convective heat transfer analysis of nanofluid flowing inside a

Received: 30 July 2017; Accepted: 29 September 2017; Published: 8 October 2017

Consider the function f(x, y). Recall that we can approximate f(x, y) with a linear function in x and y:

Products and Convolutions of Gaussian Probability Density Functions

Accuracy of free-energy perturbation calculations in molecular simulation. I. Modeling

. This is the Basic Chain Rule. x dt y dt z dt Chain Rule in this context.

Updated: January 16, 2016 Calculus II 7.4. Math 230. Calculus II. Brian Veitch Fall 2015 Northern Illinois University

Limits of algebraic functions *

Problem Set. Problems on Unordered Summation. Math 5323, Fall Februray 15, 2001 ANSWERS

Function Operations. I. Ever since basic math, you have been combining numbers by using addition, subtraction, multiplication, and division.

Partial Fraction Decomposition Honors Precalculus Mr. Velazquez Rm. 254

Math 216A. A gluing construction of Proj(S)

Analysis of the regularity, pointwise completeness and pointwise generacy of descriptor linear electrical circuits

Natural convection in a vertical strip immersed in a porous medium

Lecture 2 The First Law of Thermodynamics (Ch.1)

Physics 5153 Classical Mechanics. Solution by Quadrature-1

Strong Lyapunov Functions for Systems Satisfying the Conditions of La Salle

Sound Attenuation at High Temperatures in Pt

RATIONAL FUNCTIONS. Finding Asymptotes..347 The Domain Finding Intercepts Graphing Rational Functions

8.4 Inverse Functions

Objectives. By the time the student is finished with this section of the workbook, he/she should be able

CHAPTER 5 Reactor Dynamics. Table of Contents

Exponential and Logarithmic. Functions CHAPTER The Algebra of Functions; Composite

y,z the subscript y, z indicating that the variables y and z are kept constant. The second partial differential with respect to x is written x 2 y,z

Radiation Effects on MHD Free Convective Heat and Mass Transfer Flow Past a Vertical Porous Flat Plate with Suction

Integration by partial fractions

Seismic velocity decrement ratios for regions of partial melt near the core-mantle boundary

8. THEOREM If the partial derivatives f x. and f y exist near (a, b) and are continuous at (a, b), then f is differentiable at (a, b).

PARTIAL FRACTION DECOMPOSITION. Mr. Velazquez Honors Precalculus

( x) f = where P and Q are polynomials.

Basic properties of limits

Controlling the Heat Flux Distribution by Changing the Thickness of Heated Wall

Math 2412 Activity 1(Due by EOC Sep. 17)

3.5 Graphs of Rational Functions

Core Mathematics 3 Algebra

Grade 9 Full Year 9th Grade Review

A REPORT ON PERFORMANCE OF ANNULAR FINS HAVING VARYING THICKNESS

Review of Integration Techniques

The Ascent Trajectory Optimization of Two-Stage-To-Orbit Aerospace Plane Based on Pseudospectral Method

Convective effects in air layers bound by cellular honeycomb arrays *

Basic mathematics of economic models. 3. Maximization

Modern Ab-initio Calculations Based on Tomas-Fermi-Dirac Theory with High-Pressure Environment

y2 = 0. Show that u = e2xsin(2y) satisfies Laplace's equation.

The concept of limit

Calculators are NOT permitted.

Bottom Ekman Pumping with Stress-Dependent Eddy Viscosity

Design criteria for Fiber Reinforced Rubber Bearings

Math 106 Fall 2014 Exam 2.1 October 31, ln(x) x 3 dx = 1. 2 x 2 ln(x) + = 1 2 x 2 ln(x) + 1. = 1 2 x 2 ln(x) 1 4 x 2 + C

Pressure Dependent Compressibility of Single Carbon Nanotubes and Graphite

CHAPTER 4 Reactor Statics. Table of Contents

Acoustic forcing of flexural waves and acoustic fields for a thin plate in a fluid

ROBUST STABILITY AND PERFORMANCE ANALYSIS OF UNSTABLE PROCESS WITH DEAD TIME USING Mu SYNTHESIS

Properties of solids under the effect of high temperature and pressure- NaCI as an example

Applications of local fractional calculus to engineering in fractal time-space:

Lecture 25: Heat and The 1st Law of Thermodynamics Prof. WAN, Xin

Mathematics 136 Calculus 2 Everything You Need Or Want To Know About Partial Fractions (and maybe more!) October 19 and 21, 2016

4.5 Integration of Rational Functions by Partial Fractions

OE4625 Dredge Pumps and Slurry Transport. Vaclav Matousek October 13, 2004

Logarithmic and Exponential Equations and Change-of-Base

Relating axial motion of optical elements to focal shift

0,0 B 5,0 C 0, 4 3,5. y x. Recitation Worksheet 1A. 1. Plot these points in the xy plane: A

Root Arrangements of Hyperbolic Polynomial-like Functions

Optimal Control. with. Aerospace Applications. James M. Longuski. Jose J. Guzman. John E. Prussing

On a Closed Formula for the Derivatives of e f(x) and Related Financial Applications

Lab on Taylor Polynomials. This Lab is accompanied by an Answer Sheet that you are to complete and turn in to your instructor.

ELLIPTIC PROBLEMS INVOLVING THE 1 LAPLACIAN AND A SINGULAR LOWER ORDER TERM

Topic02_PDE. 8/29/2006 topic02_pde 1. Computational Fluid Dynamics (AE/ME 339) MAE Dept., UMR

UNDERSTANDING of unsteady combustion of solid

Analysis of Non-Thermal Equilibrium in Porous Media

THE GAMMA FUNCTION THU NGỌC DƯƠNG

In the index (pages ), reduce all page numbers by 2.

Scattering of Solitons of Modified KdV Equation with Self-consistent Sources

ELLIPTIC PROBLEMS INVOLVING THE 1 LAPLACIAN AND A SINGULAR LOWER ORDER TERM

Transcription:

Indian Journal o ure & Applied hysics Vol. 5, October, pp. 734-738 ressure derivatives o bulk modulus or materials at extreme compression Singh* & A Dwivedi Department o hysics, Institute o Basic Sciences, Dr B R Ambedkar University, hari Campus, Agra, India *E-mail: pramod@yahoo.com Received December ; revised 4 March ; accepted 7 June The method based on the calculus o indeterminates or demonstrating that all the physically acceptable equations o state satisy the identities or the pressure derivatives o bulk modulus o materials at extreme compression, has been developed. The speciic examples o the Birch-Murnaghan inite strain equation, the oirier-tarantola logarithmic equation, the Rydberg-Vinet potential energy equation, the eane -primed equation the Stacey reciprocal -primed equation, have been considered. Expressions or the bulk modulus its pressure derivatives have been derived reduced to the limit o ininite pressure. The expressions thus obtained are useul or urther analysis o higher derivative thermoelastic properties. eywords: ressure derivatives, Bulk modulus, Equations o state, Extreme compression behaviour Introduction For investigating high-pressure properties o materials, we need equations o state representing the relationships between pressure volume V at a given temperature T,. Some important equations o state such as the Birch-Murnaghan inite strain equation 3, the oirier-tarantola logarithmic equation 4, the Rydberg-Vinet potential energy equation 5,6, the eane -primed equation 7 the Stacey reciprocal -primed equation 8 have been widely used or various materials. Theory All o these equations reveal that pressure bulk modulus both increase rapidly with the decreasing volume. both become ininite in the limit o extreme compression (V), but their ratio remains inite such that : = () where is the value o = d d, the pressure derivative o bulk modulus at ininite pressure. It should be mentioned that represents the rate o increase o bulk modulus with the increase in pressure. Value o has been ound to decrease continuously with the increase in pressure attains a constant positive value equal to in the limit o extreme compression. The basic condition or the physical acceptability o an equation o state (EOS) is that in the limit o extreme compression (V), the pressure must approach ininity, must remain inite positive. In the limit V tends to zero, changes with volume as V so that the bulk modulus = V d dv changes as V. This is consistent with Eq. (). All equations o state with greater than zero satisy Eq. () as demonstrated by Stacey. In the present study, we demonstrate that various equations o state satisying Eq. (), which has the status o an algebraic identity, also satisy the identities or the second order third order pressure derivatives o bulk modulus at ininitely large pressures. It should be emphasized that the boundary conditions at extreme compression or ininite pressure are equally important as they are at zero pressure. These boundary conditions are to be satisied by all physically acceptable equations o state. The Brich-Murnaghan EOS, the oirier-tarantola logarithmic EOS, the generalized Rydberg-Vinet EOS, all can be represented by the ollowing common ormula: ( x ) = + () where ( is a unction o x=v/v, V is the value o volume V at =. Values o ( are dierent

SINGH & DWIVEDI: RESSURE DERIVATIVES OF BUL MODULUS 735 or dierent EOS. Thus or the Birch-Murnaghan ourth order EOS, =/3, 3 4 ( 4 3) Bx + ( 3) 3 Ax Cx = 4 3 Ax Bx + Cx 3 D where the constants A,B,C,D are given below: (3) 5 A = + ( 4)( 5) + (4) B = 3 + ( 4)( 3 3) + (5) 5 C = 3 + ( 4)( 3 ) + (6) 35 D = + ( 4)( 3) + (7) where, are the values o, = d d, all at =. In case o the third-order Birch-Murnaghan EOS, we have = 3 D= (Eq.7), so that: = ( 4 3) Ax ( 3) Bx 3 4 3 4 3 3 Ax Bx + C (8) For the oirier-tarantola logarithmic ourth order EOS we have =, where A ln x + 3 A (ln = () ln x + A (ln + A (ln 3 A = ( ) () A = ( + 3 + 3) () 6 In case o the third order oirier-tarantola logarithmic EOS, we have = A =, so that: In case o the generalized Rydberg-Vinet EOS, is a material-dependent parameter related to the zeropressure parameters as ollows : 5 = + 4 + + 36 (3) For this EOS, we have : x x = +η 3 3 3 3( x ) 3 (4) η = ( 3 ) 3 + (5) It is interesting to note that in the limit o extreme compression, V or x, ( becomes zero, Eq. () reduces to Eq. (). This is true or all the equations o state considered above, demonstrating the universality o Eq. (). Expressions or pressure derivatives o bulk modulus up to third order are obtained here by dierentiating Eq. () successively with respect to. Thus, we ind: = x ( (6) x x + = + (7) + 4 + 6 3 = x + 3 x + x (8) where = d dx, = d dx 3 3 =. In the derivation o Eqs (6-8), d dx we have used the ollowing relationship: d d = V = x dv dx Eqs (6-8) reveal that in the limit x: () A ln = () ln x + A (ln = ()

736 INDIAN J URE & AL HYS, VOL 5, OCTOBER = () = () provided the ollowing conditions are satisied in the limit o extreme compression (. x = (3) x 3 x = (4) = (5) In addition to the condition ( at x. The unction ( expressed by Eqs (3), (8), (), () (4) based on dierent EOS satisy these conditions. The quantities in Eqs () to () become individually zero, but their ratios remain inite at extreme compression 3,4. To demonstrate this we divide Eq. (7) by Eq. (6) to get : x x + x x + = (6) x Now taking the limit x, using Eq. () evaluating x x with the help o Eqs (3), (8), (), () (4), we ind: = 3 = = 3 (7) (8) () These equations are the results based on dierent EOS. The Brich-Murnaghan third order as well as ourth order EOS (Eqs 3 8) yield the same result in the orm o Eq. (7). The oirier-tarantola logarithmic third order as well as ourth order EOS (Eqs ) both give the same result as given by Eq. (8). The generalized Rydberg-Vinet EOS based on Eq. (4) or ( yields Eq. (). The ratio o third order pressure derivative second pressure derivative o bulk modulus can be obtained by dividing Eq. (8) by Eq. (7), rearranging the terms in numerator (N) denominator (D) in the ollowing manner : 3 N x + 3 x + x = D x + x where N = + 4 + 6 D = + (3) (3) (3) In the limit o extreme compression or ininite pressure, we have: N D 3 = = + (33) (34) Eqs (33) (34) have been obtained by using = taking ( ) to be inite at extreme compression. In view o Eqs (33) (34), Eq. (3) becomes : N D 3 x ( + 3x ( + x = x ( + x ( x (35) Values o ( its derivatives are evaluated using the expressions or ( based on dierent EOS. The ollowing results are obtained using Eq. (35) with the help o Eqs (7-), (33) (34). = + 3 (36)

SINGH & DWIVEDI: RESSURE DERIVATIVES OF BUL MODULUS 737 = (37) = + (43) = + 3 (38) The Birch-Murnaghan third order ourth order EOS yield Eq. (36), the oirier-tarantola logarithmic third order ourth order EOS give Eq. (37), the generalized Rydberg-Vinet EOS yields Eq.(38). It should be mentioned that have the same units o Ga, so / is dimensionless. is also dimensionless. is multiplied by, by, so that both are dimensionless. In the extreme compression limit, Eq. (6) reduces to the ollowing: x + x + = x x (3) Using the calculus o indeterminates, we have the ollowing relationship in the limit x, 3 x + x x + 3 x + x = x x + x (4) Eq. (4) is valid under the conditions given by Eqs (3) to (5) in addition to the condition that ( tends to zero in the limit x. Eq. (4) gives the ollowing relationship with the help o Eqs (35) (3), N = D (4) Substituting Eq. (33) or N, Eq.(34) D in Eq.(4), we get: = (4) Eq. (4) is satisied not only by all such EOS which are represented by Eq. (), but also some EOS based on the expressions or. The eane EOS is written as 7 : which gives: = = (44) (45) Eqs (44) (45) satisy Eq. (4). The reciprocal -primed EOS is given below 8 : = + Eq.(46) gives: ( ) = ( + ) = (46) (47) (48) Eqs (47) (48) satisy Eq. (4). In act, all equations o state which satisy Eq. () also satisy Eq. (4). 3 Results Discussion It should be emphasized that the pressure derivatives o bulk modulus are o central importance or determining thermoelastic properties o materials at high pressures high temperatures. The ormulations presented here are related to dierent equations o state which have been used recently 5- or investigating properties o materials at high pressures. We may thus conclude that Eq. (4) has the status o an identity which can be used or determining the third-order Grüneisen parameter. Acknowledgement The authors are thankul to Dr Jai Shanker Dr B S Sharma or helpul discussion.

738 INDIAN J URE & AL HYS, VOL 5, OCTOBER Reerences Anderson O L, Equations o State o Solids or Geophysics Ceramic Science (Oxord University ress, Oxord), 5,. Stacey F D & Davis M, hys Earth lanet Inter, 4 (4) 37. 3 Birch F, J Geophys Res, (86) 44. 4 oirier J & Tarantola A, hys Earth lanet Inter, (8). 5 Rydberg R, Z hysik, 73 (3) 376. 6 Vinet, Rose J H, Ferrante J & Smith J R, J hys: Condens Matter, (87) 4. 7 eane A, Aust J hys, 7 (54) 33. 8 Stacey F D, Geophys J Int, 6 () 43. Stacey F D, hys Earth lanet Inter, (5) 8. Singh, Indian J ure & Appl hys, 48 () 43. Stacey F D, Rep rog hys, 34 (5) 68. Shanker J, Dulari & Singh, hysica B, 44 () 483. 3 Shanker J & Singh B, hysica B, 37 (5) 78. 4 Shanker J, Singh B & Jitendra, Condensed Matter hys, () 5. 5 Upadhyay A & Sharma B S, Indian J ure & Appl hys,47 () 36. 6 ushwah S S & Bhardwaj N, Indian J ure &Appl hys, 47 () 663. 7 umar M S & Sharma B S, Indian J ure & Appl hys, 48 () 8. 8 Upadhyay A & Sharma B S, Indian J ure & Appl hys, 4 () 3. ushwah S S & Tomar Y S, Indian J ure & Appl hys, 4 (). Singh B & Gajendra S, Indian J ure & Appl hys, 4 () 467.