Two Dimensional Consolidations for Clay Soil of Non-Homogeneous and Anisotropic Permeability

Size: px
Start display at page:

Download "Two Dimensional Consolidations for Clay Soil of Non-Homogeneous and Anisotropic Permeability"

Transcription

1 Two Diensional Consolidations for Clay Soil of Non-Hoogeneous and Anisotropic Pereability Ressol R. Shakir, Muhaed Majeed Thiqar University, College of Engineering, Thiqar, Iraq University of Technology, Building and Construction Departent, Baghdad, Iraq Abstract: In this research a coputer progra was odified to solve two-diensional consolidation probles. The odified coputer progra results showed a very good agreeent copared with the analytical solution of two-diensional consolidation. The nonhoogeneity in one and two directions with the an-isotropic effect of pereability was included in the progra. It was found that the increasing of anisotropic ratio of pereability increases the average degree of consolidation while the increase in non-hoogeneous pereability ratio decreases the average degree of consolidation according to the drainage face taken there. This effect should be taken in to account when dealing with the analysis of consolidation of copressible soil and average degree of consolidation. Keywords: Finite Eleent, Consolidation, pereability. Introduction The conventional Terzaghi theory [] is derived based on assuptions such as the soil fluid flow is in one direction only and also that the soil is hoogenous and isotropic. These assuptions are ost often-violated in real life situations. They are functioned to siulate a special case of subjecting soil layer to unifor load such as fill load. The hoogeneous and isotropic soil eans that every soil paraeters are constant with any direction and so the pereability of soil is assued constant along the depth of soil layer. Two and three-diensional consolidation behavior that governed by two and three diffusion equations were presented by any investigators to overcoe the defect of assuing fluid flow in one direction [, 3, 4]. Rendulic- Terzaghi equation is considered as extending to Terzaghi's equation. For siple geoetry and one diensional consolidation, analytical solution can be used, otherwise nuerical analysis such finite difference and finite eleent ethod are utilized. Using two diensional consolidation analyses is ore realistic than one diensional consolidation. A considerable aount of research is achieved to include the effect of non hoogeneous and isotropic soil paraeters into account when analyzing consolidation behavior. Analytical solution for consolidation in two layer soil colun was provided by Zhu and Yin (005)[5]. A finite difference approach was used to solve the consolidation settleent for a layer with variable copressibility, pereability and coefficient of consolidation [6]. Xie et al. presented a fully explicit analytical solution for one diensional consolidation of two layered soils with partially drained boundaries. Airi and Esaelly [8] investigated the effect of variable pereability and copressibility on the consolidation behavior. In the present work, the pereability is considered variable in one or two different direction and is included in the finite eleent code of two diensional consolidation settleent proble. The study contribute in investigating the effect of an isotropic pereability and non hoogeneous pereability on the average degree of consolidation. Geverning Equations The consolidation governing equation for excess pore water pressure u is c vx u u u + cvy = () x y t Where c vx, c vy are the coefficient of consolidation in x, and y k direction and Cv = γ w v It is practically accepted to assue C v as a constant through the consolidation progressing since k and v both decrease through the process. Where K is the pereability and v is the coefficient of volue change. According to Griffiths and Sith [9], the discritization of equation using Galerkin ethod leads to the following equation du []{} u + [ ] = {} 0 k () Where [ K ] and [ ] are the fluid conductivity and ass atrices respectively. The above first order tie dependent equation can be solved by finite difference as follow using eleent assebly ethod du [ K]{} U 0 + [ M ] = {} 0 (3) 0 0 du [ K]{} U + [ M ] = {} 0 (4) Paper ID: J03370 of 6

2 A third equation advances the solution fro 0 to using a weighted average of the gradient at the beginning and end of tie interval, thus, du du {} U = {} U + Δt ( -θ) + (5) ([ M ] + [ ]){ } ([ ] ( ) [ K ]){ U } 0 (6) 0 θ 0 θ Δt K U = M - θ Δt The above equation gives the distribution of excess pore water pressure [0] 3. Verification of the Coputer Progra To check the accuracy of the results of finite eleent progra, a proble which represents the dissipation of excess pore water pressure fro a square of soil [9]. The square of soil is previous fro four directions, and because of syetry, it can be used as a part of square of soil (Figure ). The load is unifor and the coefficient of consolidation is constant with tie and equal in agnitude. The pressure at the center of the esh is calculated at tie interval 0., and then the results are plotted in Figure, where the results are copared with series solution values obtained by []. The crude finite eleent idealization gives excellent results. The solution of the finite eleent progra is copared with analytical solution giving good agreeent. Figure : Coparisons of finite eleent results with series results, After Sith and Griffiths Variation of Pereability The finite eleent progra was used to study the effect of the variable pereability on the average degree of consolidation for the two-diensional consolidation proble. The variation of pereability can be discussed as follows:. The pereability in x direction is different fro that in y direction (anisotropic pereability).. The pereability is variable with depth in one direction (non- hoogeneous pereability in one direction) 3. The pereability is variable with depth in each direction (non-hoogeneous pereability in two direction) The proble will depend on the type of variation of pereability (if it is stratified, linear or nonlinear variation) and if the type of variation is the sae in two directions or different. The progra can be odified to include all these types of variation of pereability. In this study the variation of pereability is used as linear variation with depth. The sketch for the pereability variation with depth has been introduced to siplify obtaining the average degree of consolidation fro these variables (Figure 3). The variables considered in this study are as follows: Figure : Grid of the verification proble. k x is the pereability in x direction. k y is the pereability in y direction 3. k x /k y is the anisotropic pereability ratio (-0) 4. k x /k x is the non-hoogeneous pereability ratio in the x direction (-8) 5. k y /k y is the non-hoogeneous pereability ratio in the y direction (-8) 6. H h /H v is the ratio of drainage path in x direction to drainage path in y direction. 7. X is the distance fro the centerline of the esh to the right side in x direction. To study the effect of the factor (X/H h ) and the anisotropic pereability ratio (k x /k y ) on the average degree of consolidation, two cases are studied. For the first case the Paper ID: J03370 of 6

3 pereability of soil is considered isotropic and hoogeneous and the excess pore pressures are recorded at different location fro the centerline of the esh, X/H h. The second case considered the pereability of soil as anisotropic paraeter and the excess pore pressures are recorded at two locations fro the centerline of the esh X/H h = 0 and X/H h = 0.4. Figure (4) shows the effects of factor (X/H h ) on the average degree of consolidation for values (0, 0., 0.4, 0.6, 0.8). It shows that the average degree of consolidation increases with the increasing of (X/H h ). The increase of (X/H h ) eans the region that the average degree of consolidation is calculated close to the drainage face. Figure (5) shows the effect of anisotropic ratio of pereability on the average degree of consolidation for different ratio of X/H h (0, and 0.4), and different ratio of anisotropic ratio of pereability (, 4, 6, 0). The average degree of consolidation increases at anisotropic ratio of pereability greater than one. Also increase of X/H h will increase the average degree of consolidation. K NY =k y -(NY(k y -k y )/NYE) K NX =k x -(NX(k x -k x )/NXE) Figure 3: Linear variation of pereability 5. Effect of Non Hoogeneous Pereability Ratio The results of the finite eleent progra for all cases studied are presented as average degree of consolidation, U av. It can be coputed at any location in the esh by the following forula (U av = -u e /u o ) where u e is the excess pore water pressure and u o is the initial pore water pressure. For the case of non-hoogeneous pereability ratio in y direction k y /k y =8 and anisotropic pereability ratio k x /k y >, when the anisotropic pereability ratio increases fro 4 to 0, keeping constant value of non-hoogeneous pereability ratio, the average degree of consolidation increases fro 0.7 to 0.9 (Table ). The increase in U av is due to the increase in anisotropic pereability ratio. Effect of the factor (X/H h ) at which the average degree of consolidation was calculated are also presented in Table. Two cases are discussed the first is for U av coputed at X/H h = 0 and the second is for U av coputed at X/H h = 0.4. The U av calculated at X/H h =0.4 is greater than U av calculated at X/H h =0 because this region is close to the right face of the esh which is idealized as a drainage face in the presented proble. Figure 4: average degree of consolidation at different location X/H h for hoogeneous and isotropic soil. Table : The effect of the anisotropic ratio of pereability and distance ratio X/H h on the average degree of consolidation H v /H h K x /K y X/H h Tv K y /K y U av For the case of non-hoogeneous pereability in two directions together with the anisotropic pereability the U av increases due to anisotropic of pereability of soil. With high pereability ratio the U av decreases because the high pereability is concentrated close to the ipervious face. The high ratio of non-hoogeneous pereability decreases the average degree of consolidation (Table ) because the high pereability in x direction and y direction are close to the ipervious side of the esh. Table : show the effect of the non-hoogeneity in two direction and the distance ratio X/H h on the average degree of consolidation H v /H h K x /K y X/H h Tv K y /K y, K x /K x U av , Table also shows the result of average degree of consolidation for the case of two diensional consolidation for soil with anisotropic pereability and non hoogeneous pereability in x and y direction. The average degree of consolidation is calculated at two location represented by X/H h. It is shown that the average degree of consolidation increases as the location closes to the face of drainage. Paper ID: J of 6

4 non-hoogeneous pereability in x-direction eans increase in pereability in the region that is far fro the drainage face of the soil while the region close to the drainage face still has the sae value of pereability. Figure 9a and Figure0a show the U av versus tie factor, Tv, for the sae previous case of soil with non-hoogeneous pereability in y direction and hoogenous pereability in x direction but for anisotropic pereability k x /k y = 4 and k x /k y = 0. The results of U av are approxiately the sae. The anisotropic pereability has no effect on the U av. Figure 5: average degree of consolidation at X/H h =0 for different anisotropic pereability ratio Figure 7: Grid of the studied proble Figure 6: average degree of consolidation at X/H h =0.4for different anisotropic pereability ratio 6. Effect of Aniostropic and Hoogeneous Pereability The study also considered rectangular esh of diensions H v and H h = H v (Figure 7). Three cases are studied, the first case is for soil with isotropic pereability and nonhoogeneous pereability in y or/and in x direction where the U av is calculated at X/H h = 0.4. The effect of anisotropic pereability on the U av are studied in the other two cases by using anisotropic pereability ratio of k x /k y = 4 and 0. The non-hoogeneous and anisotropic pereability reflects the nature of assued soils which approxiately represent real soil. Figure 8a shows the result of U av for case of soil with isotropic and non-hoogeneous pereability in y direction and hoogeneous pereability in x direction. Figure 8b shows the results of U av for the case of soil with isotropic and non-hoogeneous pereability in x and y direction. Fro the two Figures, it is found the effect of nonhoogeneous pereability ratio in x-direction has slightly effect on U av. It is attributed to that the increase of ratio of Figure 8a: Average degree of consolidation, U av, versus tie factor, Tv, for soil with isotropic pereability, k x /k y =, hoogeneous pereability in x direction and nonhoogeneous pereability in x direction This ay be attributed to that the high pereability is close to the ipervious face. Figure 9b and Figure 0b show the result of U av against Tv for soil with non-hoogeneous pereability in both x and y direction and with sae ratio of anisotropic pereability. It is obvious that the U av increases as non-hoogeneous pereability increases fro k x /k x = to k x /k x =8. Paper ID: J of 6

5 Figure 8b: Average degree of consolidation, U av, versus tie factor, Tv, for soil with isotropic pereability, k x /k y = 0, non-hoogeneous pereability in y direction and nonhoogeneous pereability in x direction Figure 9b: Average degree of consolidation, U av, versus tie factor, Tv, for soil with anisotropic pereability, k x /k y = 4, non-hoogeneous pereability in y direction and nonhoogeneous pereability in x direction Figure 9a: Average degree of consolidation, U av, versus tie factor, Tv, for soil with anisotropic pereability, k x /k y = 4, non-hoogeneous pereability in y direction and hoogeneous pereability in x direction Figure 0a: Average degree of consolidation, U av, versus tie factor, Tv, for soil with anisotropic pereability, k x /k y =0, non-hoogeneous pereability in y direction and hoogeneous pereability in x direction Paper ID: J of 6

6 Figure 0b: Average degree of consolidation, U av, versus tie factor, Tv, for soil with anisotropic pereability k x /k y =0, non-hoogeneous pereability in y direction and nonhoogeneous pereability in x direction 7. Conclusion In this study a coputer progra which depends on the uncoupled solution using finite eleent ethod has been used to solve proble of two- diensional consolidation settleent. This coputer progra is a good tool to study the non-hoogeneous and anisotropic pereability effects on the average degree of consolidation. According to the study, the following conclusions are drawn:. It is concluded that the average degree of consolidation decreases with the decrease of pereability with depth "non-hoogeneous pereability in y direction.. The soil layer of low pereability which located close to the drainage face reduces the average degree of consolidation. 3. The average degree of consolidation increases when the anisotropic pereability ratio increases ore than one 4. Including real variation of pereability in the soil to copute the average degree of consolidation has to be depended on in design. References [] Terzaghi, K.,(943) Theoretical Soil Mechanics, Wiley Publications, New York. [] Rendulic L.,(936), Porenziffer and Porenwasserdruck in Tonen entioned by N. Foldin and B. Brows, In Soft clay in engineering, 98, Elsevier Scientific publishing Copany. [3] Carillo, (94), Siple Two and Three Diensional Cases in the Theory of Consolidation of Soils, J. Math. And Phys., No., P. -5. [4] Poskitt, T.,J.(970) Settelent Charts for An Isotropic Soil, Geotechnique, Vol.0, No.3, pp [5] Zhu, G. and Yin, J-H, 005. Solution charts for the consolidation of double soil layers. Canadian Geotechnical Journal, 4: [6] Abbasi, N., Rahii, H., and Javadi, A. A. And Fakher, A. (007) Finite difference approach for consolidation with variable copressibility and pereability. Coputer and Geotechnics. 34(), (007), pp [7] Kang-He Xie, Xin-Yu Xie, Xiang Gao, Theory of one diensional consolidation of two-layered soil with partially drained boundaries, Coputers and Geotechnics, Volue 4, Issue 4, June 999, Pages 65-78, ISSN X. [8] Airi S. N and Eaaely A. (04) A realistic theory of soils consolidation" Journal of science and geotechnical engineering" Vol 4, no, [9] Sith, I. M. and Griffiths, D.V. (998) Prograing the Finite Eleent Method third edition, Jon Wiley & Sons, England. [0] J Huang, DV Griffiths (00) "One-diensional consolidation theories for layered soil and coupled and uncoupled solutions by the finite-eleent ethod" Geotechnique, Volue 60. Issue 9 Pages [] Carslaw, H.S., and Jaeger, J.C.(959) Conduction of Heat in Solids, nd edn. Clarendon Press, Oxford Paper ID: J of 6

Stress Increment Solution Charts for Soil Consolidation Analysis

Stress Increment Solution Charts for Soil Consolidation Analysis Global Journal of Researches in Engineering General Engineering Volue 1 Issue 4 Version Year 01 Type: Double Blind Peer Reviewed International Research Journal Publisher: Global Journals Inc. (USA) Online

More information

On the characterization of non-linear diffusion equations. An application in soil mechanics

On the characterization of non-linear diffusion equations. An application in soil mechanics On the characterization of non-linear diffusion equations. An application in soil echanics GARCÍA-ROS, G., ALHAMA, I., CÁNOVAS, M *. Civil Engineering Departent Universidad Politécnica de Cartagena Paseo

More information

SD Numerical Simulation Technique for Hydrodynamic Flow Gas- Solids Mixing Mantilla Núñez, Irla* 1 S. De Vicente C. 2

SD Numerical Simulation Technique for Hydrodynamic Flow Gas- Solids Mixing Mantilla Núñez, Irla* 1 S. De Vicente C. 2 SD Nuerical Siulation Technique for Hydrodynaic Flow Gas- Solids Mixing Mantilla Núñez, Irla* S. De Vicente C. National University of Engineering, Lia, Perú, Polytechnic University of Madrid. *Lia 5, irlan@uni.edu.pe.

More information

Projectile Motion with Air Resistance (Numerical Modeling, Euler s Method)

Projectile Motion with Air Resistance (Numerical Modeling, Euler s Method) Projectile Motion with Air Resistance (Nuerical Modeling, Euler s Method) Theory Euler s ethod is a siple way to approxiate the solution of ordinary differential equations (ode s) nuerically. Specifically,

More information

Spine Fin Efficiency A Three Sided Pyramidal Fin of Equilateral Triangular Cross-Sectional Area

Spine Fin Efficiency A Three Sided Pyramidal Fin of Equilateral Triangular Cross-Sectional Area Proceedings of the 006 WSEAS/IASME International Conference on Heat and Mass Transfer, Miai, Florida, USA, January 18-0, 006 (pp13-18) Spine Fin Efficiency A Three Sided Pyraidal Fin of Equilateral Triangular

More information

Extension of CSRSM for the Parametric Study of the Face Stability of Pressurized Tunnels

Extension of CSRSM for the Parametric Study of the Face Stability of Pressurized Tunnels Extension of CSRSM for the Paraetric Study of the Face Stability of Pressurized Tunnels Guilhe Mollon 1, Daniel Dias 2, and Abdul-Haid Soubra 3, M.ASCE 1 LGCIE, INSA Lyon, Université de Lyon, Doaine scientifique

More information

REDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION

REDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION ISSN 139 14X INFORMATION TECHNOLOGY AND CONTROL, 008, Vol.37, No.3 REDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION Riantas Barauskas, Vidantas Riavičius Departent of Syste Analysis, Kaunas

More information

AN APPLICATION OF CUBIC B-SPLINE FINITE ELEMENT METHOD FOR THE BURGERS EQUATION

AN APPLICATION OF CUBIC B-SPLINE FINITE ELEMENT METHOD FOR THE BURGERS EQUATION Aksan, E..: An Applıcatıon of Cubıc B-Splıne Fınıte Eleent Method for... THERMAL SCIECE: Year 8, Vol., Suppl., pp. S95-S S95 A APPLICATIO OF CBIC B-SPLIE FIITE ELEMET METHOD FOR THE BRGERS EQATIO by Eine

More information

A GENERAL FORM FOR THE ELECTRIC FIELD LINES EQUATION CONCERNING AN AXIALLY SYMMETRIC CONTINUOUS CHARGE DISTRIBUTION

A GENERAL FORM FOR THE ELECTRIC FIELD LINES EQUATION CONCERNING AN AXIALLY SYMMETRIC CONTINUOUS CHARGE DISTRIBUTION A GENEAL FOM FO THE ELECTIC FIELD LINES EQUATION CONCENING AN AXIALLY SYMMETIC CONTINUOUS CHAGE DISTIBUTION BY MUGU B. ăuţ Abstract..By using an unexpected approach it results a general for for the electric

More information

Simulation of Geomechanical Behavior during SAGD Process using COMSOL Multiphysics

Simulation of Geomechanical Behavior during SAGD Process using COMSOL Multiphysics Siulation of Geoechanical Behavior during SAGD Process using COMSOL Multiphysics X. Gong 1, R.Wan *2 Departent of Civil Engineering, University of Calgary, Calgary, Alberta, Canada. *Corresponding author:

More information

A BLOCK MONOTONE DOMAIN DECOMPOSITION ALGORITHM FOR A NONLINEAR SINGULARLY PERTURBED PARABOLIC PROBLEM

A BLOCK MONOTONE DOMAIN DECOMPOSITION ALGORITHM FOR A NONLINEAR SINGULARLY PERTURBED PARABOLIC PROBLEM INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volue 3, Nuber 2, Pages 211 231 c 2006 Institute for Scientific Coputing and Inforation A BLOCK MONOTONE DOMAIN DECOMPOSITION ALGORITHM FOR A NONLINEAR

More information

HORIZONTAL MOTION WITH RESISTANCE

HORIZONTAL MOTION WITH RESISTANCE DOING PHYSICS WITH MATLAB MECHANICS HORIZONTAL MOTION WITH RESISTANCE Ian Cooper School of Physics, Uniersity of Sydney ian.cooper@sydney.edu.au DOWNLOAD DIRECTORY FOR MATLAB SCRIPTS ec_fr_b. This script

More information

Dual porosity DRM formulation for flow and transport through fractured porous media

Dual porosity DRM formulation for flow and transport through fractured porous media Boundary Eleents XXVII 407 Dual porosity DRM orulation or lo and transport through ractured porous edia T. aardzioska & V. Popov Wessex Institute o Technology, UK Abstract The ain objective o this ork

More information

IN A SENSE, every material is a composite, even if the

IN A SENSE, every material is a composite, even if the IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES, VOL. 47, NO. 11, NOVEMBER 1999 2075 Magnetis fro Conductors and Enhanced Nonlinear Phenoena J. B. Pendry, A. J. Holden, D. J. Robbins, and W. J. Stewart,

More information

Kinetic Theory of Gases: Elementary Ideas

Kinetic Theory of Gases: Elementary Ideas Kinetic Theory of Gases: Eleentary Ideas 17th February 2010 1 Kinetic Theory: A Discussion Based on a Siplified iew of the Motion of Gases 1.1 Pressure: Consul Engel and Reid Ch. 33.1) for a discussion

More information

SIMULATION OF SOIL TEMPERATURE VARIATION FOR GEOTHERMAL APPLICATIONS

SIMULATION OF SOIL TEMPERATURE VARIATION FOR GEOTHERMAL APPLICATIONS International Journal of Mechanical Engineering and Technology (IJMET) Volue 8, Issue 2, February 2017, pp. 167 175 Article ID: IJMET_08_02_020 Available online at http://www.iaee.co/ijmet/issues.asp?jtype=ijmet&vtype=8&itype=2

More information

Keywords: Estimator, Bias, Mean-squared error, normality, generalized Pareto distribution

Keywords: Estimator, Bias, Mean-squared error, normality, generalized Pareto distribution Testing approxiate norality of an estiator using the estiated MSE and bias with an application to the shape paraeter of the generalized Pareto distribution J. Martin van Zyl Abstract In this work the norality

More information

Kinetic Theory of Gases: Elementary Ideas

Kinetic Theory of Gases: Elementary Ideas Kinetic Theory of Gases: Eleentary Ideas 9th February 011 1 Kinetic Theory: A Discussion Based on a Siplified iew of the Motion of Gases 1.1 Pressure: Consul Engel and Reid Ch. 33.1) for a discussion of

More information

POD-DEIM MODEL ORDER REDUCTION FOR THE MONODOMAIN REACTION-DIFFUSION EQUATION IN NEURO-MUSCULAR SYSTEM

POD-DEIM MODEL ORDER REDUCTION FOR THE MONODOMAIN REACTION-DIFFUSION EQUATION IN NEURO-MUSCULAR SYSTEM 6th European Conference on Coputational Mechanics (ECCM 6) 7th European Conference on Coputational Fluid Dynaics (ECFD 7) 1115 June 2018, Glasgow, UK POD-DEIM MODEL ORDER REDUCTION FOR THE MONODOMAIN REACTION-DIFFUSION

More information

TABLE FOR UPPER PERCENTAGE POINTS OF THE LARGEST ROOT OF A DETERMINANTAL EQUATION WITH FIVE ROOTS. By William W. Chen

TABLE FOR UPPER PERCENTAGE POINTS OF THE LARGEST ROOT OF A DETERMINANTAL EQUATION WITH FIVE ROOTS. By William W. Chen TABLE FOR UPPER PERCENTAGE POINTS OF THE LARGEST ROOT OF A DETERMINANTAL EQUATION WITH FIVE ROOTS By Willia W. Chen The distribution of the non-null characteristic roots of a atri derived fro saple observations

More information

Numerical Studies of a Nonlinear Heat Equation with Square Root Reaction Term

Numerical Studies of a Nonlinear Heat Equation with Square Root Reaction Term Nuerical Studies of a Nonlinear Heat Equation with Square Root Reaction Ter Ron Bucire, 1 Karl McMurtry, 1 Ronald E. Micens 2 1 Matheatics Departent, Occidental College, Los Angeles, California 90041 2

More information

Transactions on Engineering Sciences vol 3, 1993 WIT Press, ISSN

Transactions on Engineering Sciences vol 3, 1993 WIT Press,  ISSN A 3D Monte Carlo seiconductor device siulator for subicron silicon transistors* MOS K. Tarnay,^ p Masszi,& A. Poppe,* R Verhas," T. Kocsis," Zs. Kohari* " Technical University of Budapest, Departent of

More information

Generalized AOR Method for Solving System of Linear Equations. Davod Khojasteh Salkuyeh. Department of Mathematics, University of Mohaghegh Ardabili,

Generalized AOR Method for Solving System of Linear Equations. Davod Khojasteh Salkuyeh. Department of Mathematics, University of Mohaghegh Ardabili, Australian Journal of Basic and Applied Sciences, 5(3): 35-358, 20 ISSN 99-878 Generalized AOR Method for Solving Syste of Linear Equations Davod Khojasteh Salkuyeh Departent of Matheatics, University

More information

Ph 20.3 Numerical Solution of Ordinary Differential Equations

Ph 20.3 Numerical Solution of Ordinary Differential Equations Ph 20.3 Nuerical Solution of Ordinary Differential Equations Due: Week 5 -v20170314- This Assignent So far, your assignents have tried to failiarize you with the hardware and software in the Physics Coputing

More information

Block designs and statistics

Block designs and statistics Bloc designs and statistics Notes for Math 447 May 3, 2011 The ain paraeters of a bloc design are nuber of varieties v, bloc size, nuber of blocs b. A design is built on a set of v eleents. Each eleent

More information

Introduction to Robotics (CS223A) (Winter 2006/2007) Homework #5 solutions

Introduction to Robotics (CS223A) (Winter 2006/2007) Homework #5 solutions Introduction to Robotics (CS3A) Handout (Winter 6/7) Hoework #5 solutions. (a) Derive a forula that transfors an inertia tensor given in soe frae {C} into a new frae {A}. The frae {A} can differ fro frae

More information

An Approximate Model for the Theoretical Prediction of the Velocity Increase in the Intermediate Ballistics Period

An Approximate Model for the Theoretical Prediction of the Velocity Increase in the Intermediate Ballistics Period An Approxiate Model for the Theoretical Prediction of the Velocity... 77 Central European Journal of Energetic Materials, 205, 2(), 77-88 ISSN 2353-843 An Approxiate Model for the Theoretical Prediction

More information

DESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS *

DESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS * IJST, Transactions of Mechanical Engineering, Vol. 39, No. M1, pp 89-100 Printed in The Islaic Republic of Iran, 2015 Shira University DESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS

More information

Numerical Solution of the MRLW Equation Using Finite Difference Method. 1 Introduction

Numerical Solution of the MRLW Equation Using Finite Difference Method. 1 Introduction ISSN 1749-3889 print, 1749-3897 online International Journal of Nonlinear Science Vol.1401 No.3,pp.355-361 Nuerical Solution of the MRLW Equation Using Finite Difference Method Pınar Keskin, Dursun Irk

More information

Uniaxial compressive stress strain model for clay brick masonry

Uniaxial compressive stress strain model for clay brick masonry Uniaxial copressive stress strain odel for clay brick asonry Heant B. Kaushik, Durgesh C. Rai* and Sudhir K. Jain Departent of Civil Engineering, Indian Institute of Technology Kanpur, Kanpur 208 016,

More information

3D acoustic wave modeling with a time-space domain dispersion-relation-based Finite-difference scheme

3D acoustic wave modeling with a time-space domain dispersion-relation-based Finite-difference scheme P-8 3D acoustic wave odeling with a tie-space doain dispersion-relation-based Finite-difference schee Yang Liu * and rinal K. Sen State Key Laboratory of Petroleu Resource and Prospecting (China University

More information

Quantum Chemistry Exam 2 Take-home Solutions

Quantum Chemistry Exam 2 Take-home Solutions Cheistry 60 Fall 07 Dr Jean M Standard Nae KEY Quantu Cheistry Exa Take-hoe Solutions 5) (0 points) In this proble, the nonlinear variation ethod will be used to deterine an approxiate solution for the

More information

COULD A VARIABLE MASS OSCILLATOR EXHIBIT THE LATERAL INSTABILITY?

COULD A VARIABLE MASS OSCILLATOR EXHIBIT THE LATERAL INSTABILITY? Kragujevac J. Sci. 3 (8) 3-44. UDC 53.35 3 COULD A VARIABLE MASS OSCILLATOR EXHIBIT THE LATERAL INSTABILITY? Nebojša Danilović, Milan Kovačević and Vukota Babović Institute of Physics, Faculty of Science,

More information

Ştefan ŞTEFĂNESCU * is the minimum global value for the function h (x)

Ştefan ŞTEFĂNESCU * is the minimum global value for the function h (x) 7Applying Nelder Mead s Optiization Algorith APPLYING NELDER MEAD S OPTIMIZATION ALGORITHM FOR MULTIPLE GLOBAL MINIMA Abstract Ştefan ŞTEFĂNESCU * The iterative deterinistic optiization ethod could not

More information

SIMULATION OF THE HEATING STEP WITHIN THE THERMOFORMING PROCESS USING THE FINITE DIFFERENCE METHOD

SIMULATION OF THE HEATING STEP WITHIN THE THERMOFORMING PROCESS USING THE FINITE DIFFERENCE METHOD SIMULATION OF THE HEATING STEP WITHIN THE THERMOFORMING PROCESS USING THE FINITE DIFFERENCE METHOD A. Fertschej 1 *, G.R. Langecker 1 University of Leoben artur.fertschej@u-leoben.at; Franz-Josef Strasse

More information

e-companion ONLY AVAILABLE IN ELECTRONIC FORM

e-companion ONLY AVAILABLE IN ELECTRONIC FORM OPERATIONS RESEARCH doi 10.1287/opre.1070.0427ec pp. ec1 ec5 e-copanion ONLY AVAILABLE IN ELECTRONIC FORM infors 07 INFORMS Electronic Copanion A Learning Approach for Interactive Marketing to a Custoer

More information

MODIFIED SERIES RESISTANCE MODEL - DETERMINATION OF MEAN CONCENTRATION BY INTEGRAL TRANSFORMATION

MODIFIED SERIES RESISTANCE MODEL - DETERMINATION OF MEAN CONCENTRATION BY INTEGRAL TRANSFORMATION MODIFIED SERIES RESISTANCE MODEL - DETERMINATION OF MEAN CONCENTRATION BY INTEGRAL TRANSFORMATION A. L. Venezuela a, R. F. Cantão a, R. S. Ongaratto b, and R. N. Haneda c a Universidade Federal de São

More information

Ocean 420 Physical Processes in the Ocean Project 1: Hydrostatic Balance, Advection and Diffusion Answers

Ocean 420 Physical Processes in the Ocean Project 1: Hydrostatic Balance, Advection and Diffusion Answers Ocean 40 Physical Processes in the Ocean Project 1: Hydrostatic Balance, Advection and Diffusion Answers 1. Hydrostatic Balance a) Set all of the levels on one of the coluns to the lowest possible density.

More information

INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY

INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY IJESRT INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY Software for Design of Shallow Foundation using Matlab Vaishali atel Departent of Civil Engineering,, GTU / Bhailalabhai & Bhikhbhai

More information

Motion of Charges in Uniform E

Motion of Charges in Uniform E Motion of Charges in Unifor E and Fields Assue an ionized gas is acted upon by a unifor (but possibly tie-dependent) electric field E, and a unifor, steady agnetic field. These fields are assued to be

More information

Accuracy of the Scaling Law for Experimental Natural Frequencies of Rectangular Thin Plates

Accuracy of the Scaling Law for Experimental Natural Frequencies of Rectangular Thin Plates The 9th Conference of Mechanical Engineering Network of Thailand 9- October 005, Phuket, Thailand Accuracy of the caling Law for Experiental Natural Frequencies of Rectangular Thin Plates Anawat Na songkhla

More information

Explicit Approximate Solution for Finding the. Natural Frequency of the Motion of Pendulum. by Using the HAM

Explicit Approximate Solution for Finding the. Natural Frequency of the Motion of Pendulum. by Using the HAM Applied Matheatical Sciences Vol. 3 9 no. 1 13-13 Explicit Approxiate Solution for Finding the Natural Frequency of the Motion of Pendulu by Using the HAM Ahad Doosthoseini * Mechanical Engineering Departent

More information

Monitoring and system identification of suspension bridges: An alternative approach

Monitoring and system identification of suspension bridges: An alternative approach Monitoring and syste identification of suspension bridges: An alternative approach Erdal Şafak Boğaziçi University, Kandilli Observatory and Earthquake Reseach Institute, Istanbul, Turkey Abstract This

More information

A Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair

A Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair Proceedings of the 6th SEAS International Conference on Siulation, Modelling and Optiization, Lisbon, Portugal, Septeber -4, 006 0 A Siplified Analytical Approach for Efficiency Evaluation of the eaving

More information

2.141 Modeling and Simulation of Dynamic Systems Assignment #2

2.141 Modeling and Simulation of Dynamic Systems Assignment #2 2.141 Modeling and Siulation of Dynaic Systes Assignent #2 Out: Wednesday Septeber 20, 2006 Due: Wednesday October 4, 2006 Proble 1 The sketch shows a highly siplified diagra of a dry-dock used in ship

More information

Optimum Design of Assembled Cavity Dies for Precision Forging Process

Optimum Design of Assembled Cavity Dies for Precision Forging Process International Syposiu on Material, Energy and Environent Enginee (ISM3E 2015) Optiu Design of Assebled Cavity Dies for Precision Forging Process Jun-song Jin1 and Xin-yun Wang1,* 1 State Key Laboratory

More information

BUCKLING OF WING SPARS UNDER COMBINED LOADING

BUCKLING OF WING SPARS UNDER COMBINED LOADING 5 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES BUCLING OF WING SPARS UNDER COMBINED LOADING David ennedy*, Dharesh C. Patel*, Carol A. Featherston* *Cardiff School of Engineering, Cardiff University,

More information

Anisotropic reference media and the possible linearized approximations for phase velocities of qs waves in weakly anisotropic media

Anisotropic reference media and the possible linearized approximations for phase velocities of qs waves in weakly anisotropic media INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS D: APPLIED PHYSICS J. Phys. D: Appl. Phys. 5 00 007 04 PII: S00-770867-6 Anisotropic reference edia and the possible linearized approxiations for phase

More information

Numerical simulations of isotropic and die compaction of powder by the discrete element method

Numerical simulations of isotropic and die compaction of powder by the discrete element method Nuerical siulations of isotropic and die copaction of powder by the discrete eleent ethod J-F. Jerier, B. Harthong, B. Chareyre, D. Ibault, F-V. Donzé & P. Doréus Laboratoire Sols, Solides, Structures,

More information

Part I: How Dense Is It? Fundamental Question: What is matter, and how do we identify it?

Part I: How Dense Is It? Fundamental Question: What is matter, and how do we identify it? Part I: How Dense Is It? Fundaental Question: What is atter, and how do we identify it? 1. What is the definition of atter? 2. What do you think the ter ass per unit volue eans? 3. Do you think that a

More information

of daping test interfered by the noise of torsional vibration test syste was lager for low signal-to-noise ratio. Therefore, only a test ethod studyin

of daping test interfered by the noise of torsional vibration test syste was lager for low signal-to-noise ratio. Therefore, only a test ethod studyin Threshold shear strain for dynaic Raberg-Osgood odel in soils based on theroechanics Xiaoxia Guo,Xiang Sun The State Key Laboratory of Structural Analysis for Industrial Equipent, Dalian University of

More information

THE EFFECT OF SOLID PARTICLE SIZE UPON TIME AND SEDIMENTATION RATE

THE EFFECT OF SOLID PARTICLE SIZE UPON TIME AND SEDIMENTATION RATE Bulletin of the Transilvania University of Braşov Series II: Forestry Wood Industry Agricultural Food Engineering Vol. 5 (54) No. 1-1 THE EFFECT OF SOLID PARTICLE SIZE UPON TIME AND SEDIMENTATION RATE

More information

EFFECT OF MATERIAL PROPERTIES ON VIBRATIONS OF NONSYMMETRICAL AXIALLY LOADED THIN-WALLED EULER-BERNOULLI BEAMS

EFFECT OF MATERIAL PROPERTIES ON VIBRATIONS OF NONSYMMETRICAL AXIALLY LOADED THIN-WALLED EULER-BERNOULLI BEAMS Matheatical and Coputational Applications, Vol. 5, No., pp. 96-07, 00. Association for Scientific Research EFFECT OF MATERIAL PROPERTIES ON VIBRATIONS OF NONSYMMETRICAL AXIALLY LOADED THIN-WALLED EULER-BERNOULLI

More information

Direct numerical simulation of matrix diffusion across fracture/matrix interface

Direct numerical simulation of matrix diffusion across fracture/matrix interface Water Science and Engineering, 2013, 6(4): 365-379 doi:10.3882/j.issn.1674-2370.2013.04.001 http://www.waterjournal.cn e-ail: wse2008@vip.163.co Direct nuerical siulation of atrix diffusion across fracture/atrix

More information

CONTINUOUS THERMODYNAMICS FINITE DIFFUSION MODEL FOR MULTICOMPONENT FUEL SPRAY EVAPORATION

CONTINUOUS THERMODYNAMICS FINITE DIFFUSION MODEL FOR MULTICOMPONENT FUEL SPRAY EVAPORATION CONTINUOUS THERMODYNAMICS FINITE DIFFUSION MODEL FOR MULTICOMPONENT FUEL SPRAY EVAPORATION Dongyao Wang 1 and Chia-fon F. Lee Departent of Theoretical and Applied Mechanics 1 Departent of Mechanical and

More information

Chapter 4: Hypothesis of Diffusion-Limited Growth

Chapter 4: Hypothesis of Diffusion-Limited Growth Suary This section derives a useful equation to predict quantu dot size evolution under typical organoetallic synthesis conditions that are used to achieve narrow size distributions. Assuing diffusion-controlled

More information

IDENTIFICATION OF STABILITY OF CONCRETE TUNNEL LINING USING COUPLED MODELING

IDENTIFICATION OF STABILITY OF CONCRETE TUNNEL LINING USING COUPLED MODELING IDENTIFICATION OF STABILITY OF CONCRETE TUNNEL LINING USING COUPLED MODELING Kaila Weiglová, Technical University in Brno, Institute of Geoechanics, Brno, Czech Republic Petr Procházka*, Czech Association

More information

Measuring Temperature with a Silicon Diode

Measuring Temperature with a Silicon Diode Measuring Teperature with a Silicon Diode Due to the high sensitivity, nearly linear response, and easy availability, we will use a 1N4148 diode for the teperature transducer in our easureents 10 Analysis

More information

RECOVERY OF A DENSITY FROM THE EIGENVALUES OF A NONHOMOGENEOUS MEMBRANE

RECOVERY OF A DENSITY FROM THE EIGENVALUES OF A NONHOMOGENEOUS MEMBRANE Proceedings of ICIPE rd International Conference on Inverse Probles in Engineering: Theory and Practice June -8, 999, Port Ludlow, Washington, USA : RECOVERY OF A DENSITY FROM THE EIGENVALUES OF A NONHOMOGENEOUS

More information

THM Processes in a Fluid-Saturated Poroelastic Geomaterial: Comparison of Analytical Results and Computational Estimates

THM Processes in a Fluid-Saturated Poroelastic Geomaterial: Comparison of Analytical Results and Computational Estimates ROCKENG09: Proceedings of the rd CANUS Rock Mechanics Syposiu, Toronto, May 009 (Ed: M.Diederichs and G.Grasselli) THM Processes in a Fluid-Saturated Poroelastic Geoaterial: Coparison of Analytical Results

More information

Net Exchange Reformulation of Radiative Transfer in the CO m Band on Mars

Net Exchange Reformulation of Radiative Transfer in the CO m Band on Mars SEPTEMBER 2005 D U F R E S N E E T A L. 3303 Net Exchange Reforulation of Radiative Transfer in the CO 2 15- Band on Mars JEAN-LOUIS DUFRESNE Laboratoire de Météorologie Dynaique, Institut Pierre Sion

More information

Settlement Estimations of Peat during Centrifuge Experiments

Settlement Estimations of Peat during Centrifuge Experiments Settleent Estiations of Peat during Centrifuge Experients Riccardo Cappa 1, Sauel Yniesta 2, Anne Lenitzer 3, Scott J. Brandenberg 4, Ali Shafiee 5 1 PhD Candidate, S.M.ASCE, University of California at

More information

INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 1, 2011

INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 1, 2011 INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volue, No 1, 11 Copyright 1 All rights reserved Integrated Publishing services Research article ISSN 976 499 Coparison of various shear deforation

More information

Physically Based Modeling CS Notes Spring 1997 Particle Collision and Contact

Physically Based Modeling CS Notes Spring 1997 Particle Collision and Contact Physically Based Modeling CS 15-863 Notes Spring 1997 Particle Collision and Contact 1 Collisions with Springs Suppose we wanted to ipleent a particle siulator with a floor : a solid horizontal plane which

More information

ANALYSIS OF THE EFFECT OF THE CHEMICAL SPECIES CONCENTRATIONS ON THE RADIATION HEAT TRANSFER IN PARTICIPATING GASES USING A MONTE CARLO METHODOLOGY

ANALYSIS OF THE EFFECT OF THE CHEMICAL SPECIES CONCENTRATIONS ON THE RADIATION HEAT TRANSFER IN PARTICIPATING GASES USING A MONTE CARLO METHODOLOGY Proceedings of the 11 th Brazilian Congress of Theral Sciences and Engineering -- ECIT 2006 Braz. Soc. of Mechanical Sciences and Engineering -- ABCM, Curitiba, Brazil,- Dec. 5-8, 2006 Paper CIT06-0366

More information

Estimation of the Mean of the Exponential Distribution Using Maximum Ranked Set Sampling with Unequal Samples

Estimation of the Mean of the Exponential Distribution Using Maximum Ranked Set Sampling with Unequal Samples Open Journal of Statistics, 4, 4, 64-649 Published Online Septeber 4 in SciRes http//wwwscirporg/ournal/os http//ddoiorg/436/os4486 Estiation of the Mean of the Eponential Distribution Using Maiu Ranked

More information

Sharp Time Data Tradeoffs for Linear Inverse Problems

Sharp Time Data Tradeoffs for Linear Inverse Problems Sharp Tie Data Tradeoffs for Linear Inverse Probles Saet Oyak Benjain Recht Mahdi Soltanolkotabi January 016 Abstract In this paper we characterize sharp tie-data tradeoffs for optiization probles used

More information

DISTRIBUTION OF THE HYDRAULIC PARAMETERS AT RIVER BENDS

DISTRIBUTION OF THE HYDRAULIC PARAMETERS AT RIVER BENDS DISTRIBUTION OF THE HYDRAULIC PARAMETERS AT RIVER BENDS Isa Issa Oran *, Riyad Hassan Al-Anbari ** and Walaa Khalil Ali *** * Assist. Professor, Foundation of Technical Education ** Assist. Professor,

More information

Ufuk Demirci* and Feza Kerestecioglu**

Ufuk Demirci* and Feza Kerestecioglu** 1 INDIRECT ADAPTIVE CONTROL OF MISSILES Ufuk Deirci* and Feza Kerestecioglu** *Turkish Navy Guided Missile Test Station, Beykoz, Istanbul, TURKEY **Departent of Electrical and Electronics Engineering,

More information

1 Analysis of heat transfer in a single-phase transformer

1 Analysis of heat transfer in a single-phase transformer Assignent -7 Analysis of heat transr in a single-phase transforer The goal of the first assignent is to study the ipleentation of equivalent circuit ethod (ECM) and finite eleent ethod (FEM) for an electroagnetic

More information

Figure 1: Equivalent electric (RC) circuit of a neurons membrane

Figure 1: Equivalent electric (RC) circuit of a neurons membrane Exercise: Leaky integrate and fire odel of neural spike generation This exercise investigates a siplified odel of how neurons spike in response to current inputs, one of the ost fundaental properties of

More information

REPRODUCING KERNEL PARTICLE METHOD FOR CONDUCTION RADIATION INTERACTION IN 3D PARTICIPATING MEDIA

REPRODUCING KERNEL PARTICLE METHOD FOR CONDUCTION RADIATION INTERACTION IN 3D PARTICIPATING MEDIA Proceedings of the Asian Conference on Theral Sciences 2017, 1st ACTS March 26-30, 2017, Jeu Island, Korea ACTS-P00337 REPRODUCING KERNEL PARTICLE METHOD FOR CONDUCTION RADIATION INTERACTION IN 3D PARTICIPATING

More information

A note on the multiplication of sparse matrices

A note on the multiplication of sparse matrices Cent. Eur. J. Cop. Sci. 41) 2014 1-11 DOI: 10.2478/s13537-014-0201-x Central European Journal of Coputer Science A note on the ultiplication of sparse atrices Research Article Keivan Borna 12, Sohrab Aboozarkhani

More information

ESTIMATE OF THE TRUNCATION ERROR OF FINITE VOLUME DISCRETISATION OF THE NAVIER-STOKES EQUATIONS ON COLOCATED GRIDS.

ESTIMATE OF THE TRUNCATION ERROR OF FINITE VOLUME DISCRETISATION OF THE NAVIER-STOKES EQUATIONS ON COLOCATED GRIDS. ESTIMATE OF THE TRUNCATION ERROR OF FINITE VOLUME DISCRETISATION OF THE NAVIER-STOKES EQUATIONS ON COLOCATED GRIDS. Alexandros Syrakos, Apostolos Goulas Departent of Mechanical Engineering, Aristotle University

More information

Answers to assigned problems from Chapter 1

Answers to assigned problems from Chapter 1 Answers to assigned probles fro Chapter 1 1.7. a. A colun of ercury 1 in cross-sectional area and 0.001 in height has a volue of 0.001 and a ass of 0.001 1 595.1 kg. Then 1 Hg 0.001 1 595.1 kg 9.806 65

More information

2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all

2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all Lecture 6 Introduction to kinetic theory of plasa waves Introduction to kinetic theory So far we have been odeling plasa dynaics using fluid equations. The assuption has been that the pressure can be either

More information

Reed-Muller Codes. m r inductive definition. Later, we shall explain how to construct Reed-Muller codes using the Kronecker product.

Reed-Muller Codes. m r inductive definition. Later, we shall explain how to construct Reed-Muller codes using the Kronecker product. Coding Theory Massoud Malek Reed-Muller Codes An iportant class of linear block codes rich in algebraic and geoetric structure is the class of Reed-Muller codes, which includes the Extended Haing code.

More information

Bootstrapping Dependent Data

Bootstrapping Dependent Data Bootstrapping Dependent Data One of the key issues confronting bootstrap resapling approxiations is how to deal with dependent data. Consider a sequence fx t g n t= of dependent rando variables. Clearly

More information

Secondary compression of clays as a local dehydration process

Secondary compression of clays as a local dehydration process Navarro, V. & Alonso, E. E. (21). Géotechnique 51, No. 1, 859 869 Secondary copression of clays as a local dehydration process V. NAVA R R O and E. E. ALONSO{ Long-ter voluetric deforations of clays are

More information

Modelling of damage in composite materials using interface elements

Modelling of damage in composite materials using interface elements 5 th European LS-DYNA Users Conference Coposites Modelling of daage in coposite aterials using interface eleents Authors: W.G. Jiang, Departent of Aerospace Engineering, University of Bristol S.R. Hallett,

More information

Seismic Analysis of Structures by TK Dutta, Civil Department, IIT Delhi, New Delhi.

Seismic Analysis of Structures by TK Dutta, Civil Department, IIT Delhi, New Delhi. Seisic Analysis of Structures by K Dutta, Civil Departent, II Delhi, New Delhi. Module 5: Response Spectru Method of Analysis Exercise Probles : 5.8. or the stick odel of a building shear frae shown in

More information

Comput. Methods Appl. Mech. Engrg.

Comput. Methods Appl. Mech. Engrg. Coput. Methods Appl. Mech. Engrg. 200 (2011) 2768 2781 Contents lists available at ScienceDirect Coput. Methods Appl. Mech. Engrg. journal hoepage: www.elsevier.co/locate/ca A posteriori error analysis

More information

Analysis of Impulsive Natural Phenomena through Finite Difference Methods A MATLAB Computational Project-Based Learning

Analysis of Impulsive Natural Phenomena through Finite Difference Methods A MATLAB Computational Project-Based Learning Analysis of Ipulsive Natural Phenoena through Finite Difference Methods A MATLAB Coputational Project-Based Learning Nicholas Kuia, Christopher Chariah, Mechatronics Engineering, Vaughn College of Aeronautics

More information

Data-Driven Imaging in Anisotropic Media

Data-Driven Imaging in Anisotropic Media 18 th World Conference on Non destructive Testing, 16- April 1, Durban, South Africa Data-Driven Iaging in Anisotropic Media Arno VOLKER 1 and Alan HUNTER 1 TNO Stieltjesweg 1, 6 AD, Delft, The Netherlands

More information

Intelligent Systems: Reasoning and Recognition. Perceptrons and Support Vector Machines

Intelligent Systems: Reasoning and Recognition. Perceptrons and Support Vector Machines Intelligent Systes: Reasoning and Recognition Jaes L. Crowley osig 1 Winter Seester 2018 Lesson 6 27 February 2018 Outline Perceptrons and Support Vector achines Notation...2 Linear odels...3 Lines, Planes

More information

Ch 12: Variations on Backpropagation

Ch 12: Variations on Backpropagation Ch 2: Variations on Backpropagation The basic backpropagation algorith is too slow for ost practical applications. It ay take days or weeks of coputer tie. We deonstrate why the backpropagation algorith

More information

A Numerical Investigation of Turbulent Magnetic Nanofluid Flow inside Square Straight Channel

A Numerical Investigation of Turbulent Magnetic Nanofluid Flow inside Square Straight Channel A Nuerical Investigation of Turbulent Magnetic Nanofluid Flow inside Square Straight Channel M. R. Abdulwahab Technical College of Mosul, Mosul, Iraq ohaedalsafar2009@yahoo.co Abstract A nuerical study

More information

PHY307F/407F - Computational Physics Background Material for Expt. 3 - Heat Equation David Harrison

PHY307F/407F - Computational Physics Background Material for Expt. 3 - Heat Equation David Harrison INTRODUCTION PHY37F/47F - Coputational Physics Background Material for Expt 3 - Heat Equation David Harrison In the Pendulu Experient, we studied the Runge-Kutta algorith for solving ordinary differential

More information

Numerical issues in the implementation of high order polynomial multidomain penalty spectral Galerkin methods for hyperbolic conservation laws

Numerical issues in the implementation of high order polynomial multidomain penalty spectral Galerkin methods for hyperbolic conservation laws Nuerical issues in the ipleentation of high order polynoial ultidoain penalty spectral Galerkin ethods for hyperbolic conservation laws Sigal Gottlieb 1 and Jae-Hun Jung 1, 1 Departent of Matheatics, University

More information

12 Towards hydrodynamic equations J Nonlinear Dynamics II: Continuum Systems Lecture 12 Spring 2015

12 Towards hydrodynamic equations J Nonlinear Dynamics II: Continuum Systems Lecture 12 Spring 2015 18.354J Nonlinear Dynaics II: Continuu Systes Lecture 12 Spring 2015 12 Towards hydrodynaic equations The previous classes focussed on the continuu description of static (tie-independent) elastic systes.

More information

An Extension to the Tactical Planning Model for a Job Shop: Continuous-Time Control

An Extension to the Tactical Planning Model for a Job Shop: Continuous-Time Control An Extension to the Tactical Planning Model for a Job Shop: Continuous-Tie Control Chee Chong. Teo, Rohit Bhatnagar, and Stephen C. Graves Singapore-MIT Alliance, Nanyang Technological Univ., and Massachusetts

More information

Numerical Solution of the Extend Graetz Problem for a Bingham Plastic Fluid in Laminar Tube Flow

Numerical Solution of the Extend Graetz Problem for a Bingham Plastic Fluid in Laminar Tube Flow Applied Matheatics 13, 3(1): 7-37 DOI: 1.593/j.a.1331.4 erical Solution of the Extend Graetz Proble for a Bingha Plastic Fluid in Lainar Tube Flow M. E. Sayed-Ahe d 1,, H. K. Saleh, Ibrahi Hady 1 Engineering

More information

Jordan Journal of Physics

Jordan Journal of Physics Volue 5, Nuber 3, 212. pp. 113-118 ARTILE Jordan Journal of Physics Networks of Identical apacitors with a Substitutional apacitor Departent of Physics, Al-Hussein Bin Talal University, Ma an, 2, 71111,

More information

Comparison of Stability of Selected Numerical Methods for Solving Stiff Semi- Linear Differential Equations

Comparison of Stability of Selected Numerical Methods for Solving Stiff Semi- Linear Differential Equations International Journal of Applied Science and Technology Vol. 7, No. 3, Septeber 217 Coparison of Stability of Selected Nuerical Methods for Solving Stiff Sei- Linear Differential Equations Kwaku Darkwah

More information

Modeling Chemical Reactions with Single Reactant Specie

Modeling Chemical Reactions with Single Reactant Specie Modeling Cheical Reactions with Single Reactant Specie Abhyudai Singh and João edro Hespanha Abstract A procedure for constructing approxiate stochastic odels for cheical reactions involving a single reactant

More information

Dr. Naser Abu-Zaid; Lecture notes in electromagnetic theory 1; Referenced to Engineering electromagnetics by Hayt, 8 th edition 2012; Text Book

Dr. Naser Abu-Zaid; Lecture notes in electromagnetic theory 1; Referenced to Engineering electromagnetics by Hayt, 8 th edition 2012; Text Book Text Book Dr. Naser Abu-Zaid Page 1 9/4/2012 Course syllabus Electroagnetic 2 (63374) Seester Language Copulsory / Elective Prerequisites Course Contents Course Objectives Learning Outcoes and Copetences

More information

An Application of HAM for MHD Heat Source Problem. with Variable Fluid Properties

An Application of HAM for MHD Heat Source Problem. with Variable Fluid Properties Advances in Theoretical and Applied Mechanics, Vol. 7, 14, no., 79-89 HIKARI Ltd, www.-hiari.co http://dx.doi.org/1.1988/ata.14.4814 An Application of HAM for MHD Heat Source Proble with Variable Fluid

More information

Handwriting Detection Model Based on Four-Dimensional Vector Space Model

Handwriting Detection Model Based on Four-Dimensional Vector Space Model Journal of Matheatics Research; Vol. 10, No. 4; August 2018 ISSN 1916-9795 E-ISSN 1916-9809 Published by Canadian Center of Science and Education Handwriting Detection Model Based on Four-Diensional Vector

More information

An improved self-adaptive harmony search algorithm for joint replenishment problems

An improved self-adaptive harmony search algorithm for joint replenishment problems An iproved self-adaptive harony search algorith for joint replenishent probles Lin Wang School of Manageent, Huazhong University of Science & Technology zhoulearner@gail.co Xiaojian Zhou School of Manageent,

More information

Bayesian Approach for Fatigue Life Prediction from Field Inspection

Bayesian Approach for Fatigue Life Prediction from Field Inspection Bayesian Approach for Fatigue Life Prediction fro Field Inspection Dawn An and Jooho Choi School of Aerospace & Mechanical Engineering, Korea Aerospace University, Goyang, Seoul, Korea Srira Pattabhiraan

More information