Numerical Simulation of One-Dimensional Wave Equation by Non-Polynomial Quintic Spline

Size: px
Start display at page:

Download "Numerical Simulation of One-Dimensional Wave Equation by Non-Polynomial Quintic Spline"

Transcription

1 IOSR Journal of Matematcs (IOSR-JM) e-issn: , p-issn: X. Volume 14, Issue 6 Ver. I (Nov - Dec 018), PP Numercal Smulaton of One-Dmensonal Wave Equaton by Non-Polynomal Quntc Splne Md. Jon Alam 1, Laek Sazzad Andalla 1 Department of Matematcs, IUBAT-Internatonal Unversty of Busness Agrculture and Tecnology, Uttara, Daka, Banglades. Department of Matematcs, Jaangrnagar Unversty, Savar, Daka, Banglades. Correspondng Autor: Md. Jon Alam Abstract: Ts paper present a numercal algortm for te lnear one-dmensonal wave equaton. In ts metod a fnte dfference approac ad been used to descrbe te tme dervatve wle quntc splne s used as an nterpolaton functon n te space dmenson. We dscuss te accuracy of te metod by expandng te equaton based on Taylor seres and mnmzng te error. Te proposed metod was egt-order accuracy n space and fourt-order accuracy n tme varables. From te computatonal pont of vew, te soluton obtaned by ts metod s n excellent agreement wt tose obtaned by prevous works and also t s effcent to use. Numercal examples are gven to sow te applcablty and effcency of te metod Date of Submsson: Date of acceptance: I. Introducton One dmensonal wave equaton arse n many pyscal and engneerng applcatons suc as contnuum pyscs, mxed models of transonc flows, flud dynamcs, etc. Many autors ave studed te numercal solutons of lnear yperbolc wave equatons by usng varous tecnques. Consder one-dmensonal wave equaton of te form u t = u c, 0 x l, t 0 (1) x wt ntal condton u(x, 0) u x, 0 = f x, = f t 3 t, 0 x l () and boundary condtons u 0, t = p t, t 0 and u l, t = q t, t 0 (3) were c and l are postve fnte real constants and f 1 x, f x, f 3 x, p 1 t, p t, q 1 (t) and q (t) are real contnues functons. Some penomena, wc arse n many felds of scence suc as sold state pyscs, plasma pyscs, flud dynamcs, matematcal bology and cemcal knetcs, can be modeled by partal dfferental equatons. Te wave equaton s of prmary mportance n many pyscal systems suc as electro termal analogy, sgnal formaton, dranng flm, water transfer n sols, mecancs and pyscs, elastcty and etc. Tere are several numercal scemes tat ave been developed for te soluton of wave equaton. In ts paper te formaton of non-polynomal quntc splne as been developed and te consstency relaton obtaned s useful to dscretze wave equaton. We present dscretzaton of te equaton by a fnte dfference approxmaton to obtan te formulaton of proposed metod. Truncaton error and stablty analyss are dscussed. In ts study, we approxmate te functons based on Taylor seres to mnmze te error term and to obtan te class of metods. Numercal experments are conducted to demonstrate te vablty and te effcency of te proposed metod computatonally. II. Formulaton of Tenson Quntc Splne Functon In recent years, many scolars ave used non-polynomal splne for solvng dfferental equatons. Te splne functon s a pecewse polynomal or non-polynomal of degree n satsfyng te contnuty of (n 1)t dervatve. Tenson quntc splne s a non-polynomal functon tat as sx parameters to be determned ence t can satsfy te condtons of two end ponts of te nterval and contnuty of frst, second, trd and fourt dervatves. We ntroduce te set of grd ponts n te nterval [0, l] n space drecton. x =, = l, = 0,1,,, n + 1 n + 1 For eac segment quntc splne p (x) s defned by DOI: / Page

2 p x = a + b x x + c (x x ) + d (x x ) 3 + e e w x x e w x x + f e w x x e w x x, = 0,1,,, n (4) Were a, b, c, d, e and f are te unknown coeffcents to be determned also ω s free parameter. If ω 0 ten p (x) reduces to quntc splne n te nterval [0, l]. To derve te unknown coeffcents, we defne p x = u, p x = u, p x = M, p x = M 4, p x = S, 4 p x = S (5) From Eq. (4) and Eq.(5) we can determne te unknown coeffcents a = u s ω 4, b = u u + S ω 4 + 3ω + S 6ω 1 ω 4 3 M 6 M c = M S ω, d = 1 M 6 M + S ω S ω S e = ω 4 (e θ e θ ) S e θ + e θ ω 4 e θ e θ, f = S ω 4 Were θ = ω and = 0, 1, n Fnally usng te contnuty of frst dervatve at te support ponts for =, 3,.., n 1, we ave u u + u 1 6 M 1 3 M 6 M = S 1 ω 4 6ω eθ + e θ ω 3 e θ e θ + eθ e θ ω 3 (6) and from contnuty of trd dervatves, we ave M M + M 1 = S 1 ω eθ + e θ ω e θ e θ + eθ e θ +S ω + eθ + e θ ω e θ e θ + S 1 ω + ω e θ e θ (7) From Eq. (6) and Eq. (7), after elmnatng S we ave te followng useful relaton for =,3,, n. u + + u 6u + u 1 + u = 0 αm + + βm + γm + βm 1 + αm (8) Were α = θ 4 θ 3 e θ e θ 1 3θ e θ e θ, β = 4θ eθ + e θ θ 3 e θ e θ + eθ + e θ 4 3θ e θ e θ, γ = 6θ 4 + 4eθ + 4e θ 3θ e θ e θ 4 + 4eθ + 4e θ θ 3 e θ e θ Wen ω 0so θ 0, ten (α, β, γ) (1, 6, 66) and te relaton defned by Eq.(8) reduce nto ordnary quntc splne u + + u 6u + u 1 + u = 0 M M + 6M 1 + M (9) III. Numercal Tecnque By usng Eq.(9) for ( + 1)t, ()t and ( 1)t tme level we ave u u +1 6u +1 + u u = 0 αm βm +1 + γm +1 + βm αm ω (10) u + + u 6u + u 1 + u = 0 αm + + βm + γm + βm 1 + αm (11) u u 1 6u 1 + u u 0 αm βm 1 + γm 1 + βm αm tat we wll use tese equatons to dscretze wave equaton. = DOI: / Page (1)

3 Wave Equaton Fnte dfference approxmaton for second order tme dervatve s u tt = u +1 u 1 + u k = u tt + O k As we consder te space dervatve s approxmated by non-polynomal tenson splne (13) u xx = p x, t = M (14) By usng Eq.(13) and Eq.(14), we can develop a new approxmaton for te soluton of wave equaton, so tat te Eq.(1) can be replaced by ηm η M + η M +1 = u +1 u 1 + u c k were 0 η 1 s a free varable (15) Agan we multply Eq.(11) by 1 η and add ts to Eq.(10) and Eq.(1) multpled by η and elmnate M, ten we obtan te followng relaton for wave Eq(1) η u +, 1 + u, η u +, + u, + η u +, u, η u + 1, 1 + u 1, η u + 1, + u 1, 6η u, η u, 6η u, c α u +, + 1 u +, + u +, 1 + u, + 1 u, k + u, 1 40c β u + 1, + 1 u + 1, + u + 1, 1 + u 1, + 1 u 1, k + u 1, 1 40c γ u, + 1 u, + u, 1 = 0, k = 1,,3 =,.. N (16) Error Estmate n Splne Approxmaton To estmate te error for wave equaton we expand Eq.(16) n Taylor seres about u(x, t ) and ten we fnd te optmal value for α, β and γ. Error Estmate for Wave Equaton For wave equaton, we expand Eq.(16) n Taylor seres and replace te dervatves nvolvng t by te relaton + u x t = + u c (17) x+ ten we derve te local truncaton error. Te prncpal part of te local truncaton error of te proposed metod for wave equaton s T = α 10 γ 0 β D x,x U 0, 0 + α 5 β α 10 γ 40 β η k c D x,x,x,x U 0, 0 + α 15 β α 60 β η 4 k c D x,x,x,x,x,x U 0, 0 + α 3600 γ 700 β η k 4 c D x,x,x,x,x,x U 0, 0 + α 5 β D x,x,x,x,x,x,x,x U 0, 0 + (18) If we coose α = 7, β = 76 and γ = 67 n Eq.(18) we obtan a new sceme of order O( k 4 ), 6 3 furtermore by coosng η = 1 we can optmze our sceme, too. 1 Stablty Analyss In ts secton, we dscuss stablty of te proposed metod for te numercal soluton of wave equaton. We assume tat te soluton of Eq. (16) at te grd pont (l, k) s u l = ξ e lθ (19) Were = 1, θ s a real number and ξ s a complex number. DOI: / Page

4 By substtutng Eq.(19) n Eq.(16) we obtan a quadratc equaton as follow Qξ + ξ + φ = 0 (0) For 1D wave equaton we ave Q = cos θ η α 10ck + cos θ 4η β 10ck 6η γ 0ck = cos θ 4η α 5ck + cos θ 4 8η β 5ck 6 + 1η + γ 10ck φ = cos θ η 10ck + cos θ 4η β 10ck 6η γ 0ck Tus we ave Q + φ ξ + Q + + φ = 0. In order to ξ < 1, we must ave < 0 and Q + φ > 0. Obvously Q + φ > 0 for eac η, f η > 1 + equaton. 7 10ck ten < 0, terefore our sceme wll be stable for wave Numercal Example We appled te presented metod to te followng wave equaton. For ts purpose we two examples for wave equaton We appled proposed metod wt α, β, γ = 1, 6, 66 (metod I) wc s ordnary quntc splne of wt order O 6 + k 4 and f we select α, β, γ = 7 6, 76 3, 67 we obtan a new metod wc s of order O 8 + 4k4(metod II). Example 1: We consder Eq.(1) wt c = 1, f x = 0, f 3 x = π cos xπ, p t = sn πt, q x = sn πt. Te exact soluton for te problem s u x, t = π cos πx sn(πt) Absolute Error for Example 1. x t Metod I Metod II Ts problem s solved by dfferent values of step sze n te x-drecton and te tme step sze t = Te soluton by proposed metod are compared wt te exact soluton at te grd ponts and te maxmum absolute errors are tabulated n above table. Te space-tme of te estmated soluton s gven n te fgure. Te maxmum absolute error of ts example by metod II s and by metod I s DOI: / Page

5 Example. Now we consder Eq.(1) wt c = 1, f x = cos πx, f 3 x = 0, p t = cos πt, q x = cos πt. Te exact soluton for te problem s Absolute Error for Example. u x, t = 1 cos π x + t + 1 cos π x t x t Metod I Metod II Ts problem s solved by dfferent values of step sze n te x-drecton and te tme step sze t = Te soluton by proposed metod are compared wt te exact soluton at te grd ponts and te maxmum absolute errors are tabulated n above table. Te space-tme of te estmated soluton s gven n te fgure. Te maxmum absolute error of ts example by metod II s and by metod I s Reference [1]. Saadatmand and M. Degan, Numercal Soluton of te One-Dmensonal Wave Equaton wt An ntegral Condton, Numercal Metods for Partal Dfferental Equatons, Vol, 3, no., pp. 8-9, 007. []. J. Rasdna, R. Jallan, and V. Kazem, Splne metods for te solutons of yperbolc equatons, Appled Matematcs and Computaton, vol. 190, no. 1, pp , 007. [3]. R. K. Moanty and V. Gopal, Hg accuracy cubc splne fnte dfference approxmaton for te soluton of one-space dmensonal non- lnear wave equatons, Appled Matematcs and Computaton, vol. 18, no. 8, pp , 011. [4]. Ryo, Fast decay of solutons for lnear wave equatons wt dsspaton localzed near nfnty n an exteror doman, J. Dfferental Equatons, 188 () (003) [5]. Rasdna J. and Mosenyzadea, Numercal Soluton of One-Dmensonal Heat and Wave Equaton by Non-Polynomal Quntc Splne, Internatonal Journal of Matematcal Modellng & Computatons Vol. 05, No. 04, ( 015), [6]. Won Y. Yang, Wenwu Cao, Tae-Sang Cung, Jon Morrs, Appled Numercal metods usng MATLAB, Wley student Edton, 013. [7]. M. Degan, On te Soluton of An ntal-boundary Value Problem ta Combnes Neumann and Integral Condton for te Wave Equaton, Numercal Metods for Partal Dfferental Equatons, Vol. 1, no. 1, pp:4-40 [8]. P. Karageorgs and Walter A. Strauss, Instablty of steady states for nonlnear wave and eat equatons, J. Dfferental Equatons, 41 (1) (007) Md. Jon Alam" Numercal Smulaton of One-Dmensonal Wave Equaton by Non-Polynomal Quntc Splne." IOSR Journal of Matematcs (IOSR-JM) 14.6 (018): DOI: / Page

Numerical Solution of One-Dimensional Heat and Wave Equation by Non-Polynomial Quintic Spline

Numerical Solution of One-Dimensional Heat and Wave Equation by Non-Polynomial Quintic Spline Internatonal Journal of Mathematcal Modellng & Computatons Vol. 05, No. 04, Fall 2015, 291-305 Numercal Soluton of One-Dmensonal Heat and Wave Equaton by Non-Polynomal Quntc Splne J. Rashdna a, and M.

More information

Solution for singularly perturbed problems via cubic spline in tension

Solution for singularly perturbed problems via cubic spline in tension ISSN 76-769 England UK Journal of Informaton and Computng Scence Vol. No. 06 pp.6-69 Soluton for sngularly perturbed problems va cubc splne n tenson K. Aruna A. S. V. Rav Kant Flud Dynamcs Dvson Scool

More information

A Spline based computational simulations for solving selfadjoint singularly perturbed two-point boundary value problems

A Spline based computational simulations for solving selfadjoint singularly perturbed two-point boundary value problems ISSN 746-769 England UK Journal of Informaton and Computng Scence Vol. 7 No. 4 pp. 33-34 A Splne based computatonal smulatons for solvng selfadjont sngularly perturbed two-pont boundary value problems

More information

Solving Singularly Perturbed Differential Difference Equations via Fitted Method

Solving Singularly Perturbed Differential Difference Equations via Fitted Method Avalable at ttp://pvamu.edu/aam Appl. Appl. Mat. ISSN: 193-9466 Vol. 8, Issue 1 (June 013), pp. 318-33 Applcatons and Appled Matematcs: An Internatonal Journal (AAM) Solvng Sngularly Perturbed Dfferental

More information

The Finite Element Method: A Short Introduction

The Finite Element Method: A Short Introduction Te Fnte Element Metod: A Sort ntroducton Wat s FEM? Te Fnte Element Metod (FEM) ntroduced by engneers n late 50 s and 60 s s a numercal tecnque for solvng problems wc are descrbed by Ordnary Dfferental

More information

TR/28. OCTOBER CUBIC SPLINE INTERPOLATION OF HARMONIC FUNCTIONS BY N. PAPAMICHAEL and J.R. WHITEMAN.

TR/28. OCTOBER CUBIC SPLINE INTERPOLATION OF HARMONIC FUNCTIONS BY N. PAPAMICHAEL and J.R. WHITEMAN. TR/8. OCTOBER 97. CUBIC SPLINE INTERPOLATION OF HARMONIC FUNCTIONS BY N. PAPAMICHAEL and J.R. WHITEMAN. W960748 ABSTRACT It s sown tat for te two dmensonal Laplace equaton a unvarate cubc splne approxmaton

More information

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES COMPUTATIONAL FLUID DYNAMICS: FDM: Appromaton of Second Order Dervatves Lecture APPROXIMATION OF SECOMD ORDER DERIVATIVES. APPROXIMATION OF SECOND ORDER DERIVATIVES Second order dervatves appear n dffusve

More information

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems Mathematca Aeterna, Vol. 1, 011, no. 06, 405 415 Applcaton of B-Splne to Numercal Soluton of a System of Sngularly Perturbed Problems Yogesh Gupta Department of Mathematcs Unted College of Engneerng &

More information

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc

More information

A New Recursive Method for Solving State Equations Using Taylor Series

A New Recursive Method for Solving State Equations Using Taylor Series I J E E E C Internatonal Journal of Electrcal, Electroncs ISSN No. (Onlne) : 77-66 and Computer Engneerng 1(): -7(01) Specal Edton for Best Papers of Mcael Faraday IET Inda Summt-01, MFIIS-1 A New Recursve

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

TR/95 February Splines G. H. BEHFOROOZ* & N. PAPAMICHAEL

TR/95 February Splines G. H. BEHFOROOZ* & N. PAPAMICHAEL TR/9 February 980 End Condtons for Interpolatory Quntc Splnes by G. H. BEHFOROOZ* & N. PAPAMICHAEL *Present address: Dept of Matematcs Unversty of Tabrz Tabrz Iran. W9609 A B S T R A C T Accurate end condtons

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

Multivariate Ratio Estimator of the Population Total under Stratified Random Sampling

Multivariate Ratio Estimator of the Population Total under Stratified Random Sampling Open Journal of Statstcs, 0,, 300-304 ttp://dx.do.org/0.436/ojs.0.3036 Publsed Onlne July 0 (ttp://www.scrp.org/journal/ojs) Multvarate Rato Estmator of te Populaton Total under Stratfed Random Samplng

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information

Multigrid Methods and Applications in CFD

Multigrid Methods and Applications in CFD Multgrd Metods and Applcatons n CFD Mcael Wurst 0 May 009 Contents Introducton Typcal desgn of CFD solvers 3 Basc metods and ter propertes for solvng lnear systems of equatons 4 Geometrc Multgrd 3 5 Algebrac

More information

Positivity Preserving Interpolation by Using

Positivity Preserving Interpolation by Using Appled Matematcal Scences, Vol. 8, 04, no. 4, 053-065 HIKARI Ltd, www.m-kar.com ttp://dx.do.org/0.988/ams.04.38449 Postvty Preservng Interpolaton by Usng GC Ratonal Cubc Splne Samsul Arffn Abdul Karm Department

More information

1 Introduction We consider a class of singularly perturbed two point singular boundary value problems of the form: k x with boundary conditions

1 Introduction We consider a class of singularly perturbed two point singular boundary value problems of the form: k x with boundary conditions Lakshm Sreesha Ch. Non Standard Fnte Dfference Method for Sngularly Perturbed Sngular wo Pont Boundary Value Problem usng Non Polynomal Splne LAKSHMI SIREESHA CH Department of Mathematcs Unversty College

More information

New Method for Solving Poisson Equation. on Irregular Domains

New Method for Solving Poisson Equation. on Irregular Domains Appled Mathematcal Scences Vol. 6 01 no. 8 369 380 New Method for Solvng Posson Equaton on Irregular Domans J. Izadan and N. Karamooz Department of Mathematcs Facult of Scences Mashhad BranchIslamc Azad

More information

Stanford University CS254: Computational Complexity Notes 7 Luca Trevisan January 29, Notes for Lecture 7

Stanford University CS254: Computational Complexity Notes 7 Luca Trevisan January 29, Notes for Lecture 7 Stanford Unversty CS54: Computatonal Complexty Notes 7 Luca Trevsan January 9, 014 Notes for Lecture 7 1 Approxmate Countng wt an N oracle We complete te proof of te followng result: Teorem 1 For every

More information

Convexity preserving interpolation by splines of arbitrary degree

Convexity preserving interpolation by splines of arbitrary degree Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete

More information

Initial-value Technique For Singularly Perturbed Two Point Boundary Value Problems Via Cubic Spline

Initial-value Technique For Singularly Perturbed Two Point Boundary Value Problems Via Cubic Spline Unversty of Central Florda Electronc Teses and Dssertatons Masters Tess (Open Access) Intal-value Tecnque For Sngularly Perturbed Two Pont Boundary Value Problems Va Cubc Splne 010 Lus G. Negron Unversty

More information

Maejo International Journal of Science and Technology

Maejo International Journal of Science and Technology Maejo Int. J. Sc. Technol. () - Full Paper Maejo Internatonal Journal of Scence and Technology ISSN - Avalable onlne at www.mjst.mju.ac.th Fourth-order method for sngularly perturbed sngular boundary value

More information

5 The Laplace Equation in a convex polygon

5 The Laplace Equation in a convex polygon 5 Te Laplace Equaton n a convex polygon Te most mportant ellptc PDEs are te Laplace, te modfed Helmoltz and te Helmoltz equatons. Te Laplace equaton s u xx + u yy =. (5.) Te real and magnary parts of an

More information

The finite element method explicit scheme for a solution of one problem of surface and ground water combined movement

The finite element method explicit scheme for a solution of one problem of surface and ground water combined movement IOP Conference Seres: Materals Scence and Engneerng PAPER OPEN ACCESS e fnte element metod explct sceme for a soluton of one problem of surface and ground water combned movement o cte ts artcle: L L Glazyrna

More information

A new Approach for Solving Linear Ordinary Differential Equations

A new Approach for Solving Linear Ordinary Differential Equations , ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of

More information

Cubic Trigonometric B-Spline Applied to Linear Two-Point Boundary Value Problems of Order Two

Cubic Trigonometric B-Spline Applied to Linear Two-Point Boundary Value Problems of Order Two World Academy of Scence Engneerng and echnology Internatonal Journal of Mathematcal and omputatonal Scences Vol: No:0 00 ubc rgonometrc B-Splne Appled to Lnear wo-pont Boundary Value Problems of Order

More information

On Pfaff s solution of the Pfaff problem

On Pfaff s solution of the Pfaff problem Zur Pfaff scen Lösung des Pfaff scen Probles Mat. Ann. 7 (880) 53-530. On Pfaff s soluton of te Pfaff proble By A. MAYER n Lepzg Translated by D. H. Delpenc Te way tat Pfaff adopted for te ntegraton of

More information

FINITE DIFFERENCE SOLUTION OF MIXED BOUNDARY-VALUE ELASTIC PROBLEMS

FINITE DIFFERENCE SOLUTION OF MIXED BOUNDARY-VALUE ELASTIC PROBLEMS 4 t Internatonal Conference on Mecancal Engneerng, December 6-8, 1, Daa, Banglades/pp. V 171-175 FINITE DIFFERENCE SOLUTION OF MIXED BOUNDARY-VALUE ELASTIC PROBLEMS S. Reaz Amed, Noor Al Quddus and M.

More information

Research Article Cubic B-Spline Collocation Method for One-Dimensional Heat and Advection-Diffusion Equations

Research Article Cubic B-Spline Collocation Method for One-Dimensional Heat and Advection-Diffusion Equations Appled Mathematcs Volume 22, Artcle ID 4587, 8 pages do:.55/22/4587 Research Artcle Cubc B-Splne Collocaton Method for One-Dmensonal Heat and Advecton-Dffuson Equatons Joan Goh, Ahmad Abd. Majd, and Ahmad

More information

2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12

2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12 REVIEW Lecture 11: 2.29 Numercal Flud Mechancs Fall 2011 Lecture 12 End of (Lnear) Algebrac Systems Gradent Methods Krylov Subspace Methods Precondtonng of Ax=b FINITE DIFFERENCES Classfcaton of Partal

More information

ORDINARY DIFFERENTIAL EQUATIONS EULER S METHOD

ORDINARY DIFFERENTIAL EQUATIONS EULER S METHOD Numercal Analss or Engneers German Jordanan Unverst ORDINARY DIFFERENTIAL EQUATIONS We wll eplore several metods o solvng rst order ordnar derental equatons (ODEs and we wll sow ow tese metods can be appled

More information

THE STURM-LIOUVILLE EIGENVALUE PROBLEM - A NUMERICAL SOLUTION USING THE CONTROL VOLUME METHOD

THE STURM-LIOUVILLE EIGENVALUE PROBLEM - A NUMERICAL SOLUTION USING THE CONTROL VOLUME METHOD Journal of Appled Mathematcs and Computatonal Mechancs 06, 5(), 7-36 www.amcm.pcz.pl p-iss 99-9965 DOI: 0.75/jamcm.06..4 e-iss 353-0588 THE STURM-LIOUVILLE EIGEVALUE PROBLEM - A UMERICAL SOLUTIO USIG THE

More information

Numerical Solution of Singular Perturbation Problems Via Deviating Argument and Exponential Fitting

Numerical Solution of Singular Perturbation Problems Via Deviating Argument and Exponential Fitting Amercan Journal of Computatonal and Appled Matematcs 0, (): 49-54 DOI: 0.593/j.ajcam.000.09 umercal Soluton of Sngular Perturbaton Problems Va Devatng Argument and Eponental Fttng GBSL. Soujanya, Y.. Reddy,

More information

Chapter 4 The Wave Equation

Chapter 4 The Wave Equation Chapter 4 The Wave Equaton Another classcal example of a hyperbolc PDE s a wave equaton. The wave equaton s a second-order lnear hyperbolc PDE that descrbes the propagaton of a varety of waves, such as

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family IOSR Journal of Mathematcs IOSR-JM) ISSN: 2278-5728. Volume 3, Issue 3 Sep-Oct. 202), PP 44-48 www.osrjournals.org Usng T.O.M to Estmate Parameter of dstrbutons that have not Sngle Exponental Famly Jubran

More information

High resolution entropy stable scheme for shallow water equations

High resolution entropy stable scheme for shallow water equations Internatonal Symposum on Computers & Informatcs (ISCI 05) Hgh resoluton entropy stable scheme for shallow water equatons Xaohan Cheng,a, Yufeng Ne,b, Department of Appled Mathematcs, Northwestern Polytechncal

More information

Numerical Differentiation

Numerical Differentiation Part 5 Capter 19 Numercal Derentaton PowerPonts organzed by Dr. Mcael R. Gustason II, Duke Unversty Revsed by Pro. Jang, CAU All mages copyrgt Te McGraw-Hll Companes, Inc. Permsson requred or reproducton

More information

Polynomial Regression Models

Polynomial Regression Models LINEAR REGRESSION ANALYSIS MODULE XII Lecture - 6 Polynomal Regresson Models Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur Test of sgnfcance To test the sgnfcance

More information

Numerical Solution of two dimensional coupled viscous Burgers Equation using the Modified Cubic B-Spline Differential Quadrature Method

Numerical Solution of two dimensional coupled viscous Burgers Equation using the Modified Cubic B-Spline Differential Quadrature Method umercal Soluton of two dmensonal coupled vscous Burgers Equaton usng the odfed Cubc B-Splne Dfferental Quadrature ethod H. S. Shukla 1, ohammad Tamsr 1*, Vneet K. Srvastava, Ja Kumar 3 1 Department of

More information

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION

More information

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

More information

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI

More information

A Linearized Finite Difference Scheme with Non-uniform Meshes for Semi-linear Parabolic Equation

A Linearized Finite Difference Scheme with Non-uniform Meshes for Semi-linear Parabolic Equation Jun Zou A Lnearzed Fnte Dfference Sceme wt Non-unform Meses for Sem-lnear Parabolc Equaton Jun Zou Yangtze Normal Unversty Scool of Matematcs and Computer Scence Congqng 40800 Cna flzjzklm@6.com Abstract:

More information

TR/01/89 February An O(h 6 ) cubic spline interpolating procedure for harmonic functions. N. Papamichael and Maria Joana Soares*

TR/01/89 February An O(h 6 ) cubic spline interpolating procedure for harmonic functions. N. Papamichael and Maria Joana Soares* TR/0/89 February 989 An O( 6 cubc splne nterpolatng procedure for armonc functons N. Papamcael Mara Joana Soares* *Área de Matematca, Unversdade do Mno, 4700 Braga, Portugal. z 6393 ABSTRACT An O( 6 metod

More information

On a nonlinear compactness lemma in L p (0, T ; B).

On a nonlinear compactness lemma in L p (0, T ; B). On a nonlnear compactness lemma n L p (, T ; B). Emmanuel Matre Laboratore de Matématques et Applcatons Unversté de Haute-Alsace 4, rue des Frères Lumère 6893 Mulouse E.Matre@ua.fr 3t February 22 Abstract

More information

Appendix B. The Finite Difference Scheme

Appendix B. The Finite Difference Scheme 140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton

More information

ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION

ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION ANALELE ŞTIINŢIFICE ALE UNIVERSITĂŢII AL.I. CUZA DIN IAŞI (S.N.) MATEMATICĂ, Tomul LIX, 013, f.1 DOI: 10.478/v10157-01-00-y ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION BY ION

More information

Report on Image warping

Report on Image warping Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.

More information

Consistency & Convergence

Consistency & Convergence /9/007 CHE 374 Computatonal Methods n Engneerng Ordnary Dfferental Equatons Consstency, Convergence, Stablty, Stffness and Adaptve and Implct Methods ODE s n MATLAB, etc Consstency & Convergence Consstency

More information

Lecture 2: Numerical Methods for Differentiations and Integrations

Lecture 2: Numerical Methods for Differentiations and Integrations Numercal Smulaton of Space Plasmas (I [AP-4036] Lecture 2 by Lng-Hsao Lyu March, 2018 Lecture 2: Numercal Methods for Dfferentatons and Integratons As we have dscussed n Lecture 1 that numercal smulaton

More information

POLYNOMIAL BASED DIFFERENTIAL QUADRATURE FOR NUMERICAL SOLUTIONS OF KURAMOTO-SIVASHINSKY EQUATION

POLYNOMIAL BASED DIFFERENTIAL QUADRATURE FOR NUMERICAL SOLUTIONS OF KURAMOTO-SIVASHINSKY EQUATION POLYOMIAL BASED DIFFERETIAL QUADRATURE FOR UMERICAL SOLUTIOS OF KURAMOTO-SIVASHISKY EQUATIO Gülsemay YİĞİT 1 and Mustafa BAYRAM *, 1 School of Engneerng and atural Scences, Altınbaş Unversty, Istanbul,

More information

Chapter 3 Differentiation and Integration

Chapter 3 Differentiation and Integration MEE07 Computer Modelng Technques n Engneerng Chapter Derentaton and Integraton Reerence: An Introducton to Numercal Computatons, nd edton, S. yakowtz and F. zdarovsky, Mawell/Macmllan, 990. Derentaton

More information

Pathan Mahabub Basha *, Vembu Shanthi

Pathan Mahabub Basha *, Vembu Shanthi Amercan Journal of umercal Analyss 5 Vol 3 o 39-48 Avalable onlne at ttp://pubsscepubcom/ajna/3// Scence Educaton Publsng DOI:69/ajna-3-- A Unformly Convergent Sceme for A System of wo Coupled Sngularly

More information

6.3.4 Modified Euler s method of integration

6.3.4 Modified Euler s method of integration 6.3.4 Modfed Euler s method of ntegraton Before dscussng the applcaton of Euler s method for solvng the swng equatons, let us frst revew the basc Euler s method of numercal ntegraton. Let the general from

More information

2.29 Numerical Fluid Mechanics

2.29 Numerical Fluid Mechanics REVIEW Lecture 10: Sprng 2015 Lecture 11 Classfcaton of Partal Dfferental Equatons PDEs) and eamples wth fnte dfference dscretzatons Parabolc PDEs Ellptc PDEs Hyperbolc PDEs Error Types and Dscretzaton

More information

The Finite Element Method

The Finite Element Method The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of

More information

Hongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)

Hongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k) ISSN 1749-3889 (prnt), 1749-3897 (onlne) Internatonal Journal of Nonlnear Scence Vol.17(2014) No.2,pp.188-192 Modfed Block Jacob-Davdson Method for Solvng Large Sparse Egenproblems Hongy Mao, College of

More information

Suppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl

Suppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl RECURSIVE SPLINE INTERPOLATION METHOD FOR REAL TIME ENGINE CONTROL APPLICATIONS A. Stotsky Volvo Car Corporaton Engne Desgn and Development Dept. 97542, HA1N, SE- 405 31 Gothenburg Sweden. Emal: astotsky@volvocars.com

More information

Fixed point method and its improvement for the system of Volterra-Fredholm integral equations of the second kind

Fixed point method and its improvement for the system of Volterra-Fredholm integral equations of the second kind MATEMATIKA, 217, Volume 33, Number 2, 191 26 c Penerbt UTM Press. All rghts reserved Fxed pont method and ts mprovement for the system of Volterra-Fredholm ntegral equatons of the second knd 1 Talaat I.

More information

Additional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty

Additional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty Addtonal Codes usng Fnte Dfference Method Benamn Moll 1 HJB Equaton for Consumpton-Savng Problem Wthout Uncertanty Before consderng the case wth stochastc ncome n http://www.prnceton.edu/~moll/ HACTproect/HACT_Numercal_Appendx.pdf,

More information

Applied Mathematics Letters

Applied Mathematics Letters Appled Matheatcs Letters 2 (2) 46 5 Contents lsts avalable at ScenceDrect Appled Matheatcs Letters journal hoepage: wwwelseverco/locate/al Calculaton of coeffcents of a cardnal B-splne Gradr V Mlovanovć

More information

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS Blucher Mechancal Engneerng Proceedngs May 0, vol., num. www.proceedngs.blucher.com.br/evento/0wccm STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS Takahko Kurahash,

More information

A Hybrid Variational Iteration Method for Blasius Equation

A Hybrid Variational Iteration Method for Blasius Equation Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method

More information

2 Finite difference basics

2 Finite difference basics Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T

More information

Note 10. Modeling and Simulation of Dynamic Systems

Note 10. Modeling and Simulation of Dynamic Systems Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada

More information

The Study of Teaching-learning-based Optimization Algorithm

The Study of Teaching-learning-based Optimization Algorithm Advanced Scence and Technology Letters Vol. (AST 06), pp.05- http://dx.do.org/0.57/astl.06. The Study of Teachng-learnng-based Optmzaton Algorthm u Sun, Yan fu, Lele Kong, Haolang Q,, Helongang Insttute

More information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

More information

A Computational Method for Solving Two Point Boundary Value Problems of Order Four

A Computational Method for Solving Two Point Boundary Value Problems of Order Four Yoge Gupta et al, Int. J. Comp. Tec. Appl., Vol (5), - ISSN:9-09 A Computatonal Metod for Solvng Two Pont Boundary Value Problem of Order Four Yoge Gupta Department of Matematc Unted College of Engg and

More information

Solving Fractional Nonlinear Fredholm Integro-differential Equations via Hybrid of Rationalized Haar Functions

Solving Fractional Nonlinear Fredholm Integro-differential Equations via Hybrid of Rationalized Haar Functions ISSN 746-7659 England UK Journal of Informaton and Computng Scence Vol. 9 No. 3 4 pp. 69-8 Solvng Fractonal Nonlnear Fredholm Integro-dfferental Equatons va Hybrd of Ratonalzed Haar Functons Yadollah Ordokhan

More information

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed (2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected

More information

Septic B-Spline Collocation Method for the Numerical Solution of the Modified Equal Width Wave Equation

Septic B-Spline Collocation Method for the Numerical Solution of the Modified Equal Width Wave Equation Appled Mathematcs 79-749 do:.46/am..698 Publshed Onlne June (http://www.scrp.org/ournal/am) Septc B-Splne Collocaton Method for the umercal Soluton of the Modfed Equal Wdth Wave Equaton Abstract Turab

More information

Haar wavelet collocation method to solve problems arising in induction motor

Haar wavelet collocation method to solve problems arising in induction motor ISSN 746-7659, England, UK Journal of Informaton and Computng Scence Vol., No., 07, pp.096-06 Haar wavelet collocaton method to solve problems arsng n nducton motor A. Padmanabha Reddy *, C. Sateesha,

More information

COMPUTATIONAL METHODS AND ALGORITHMS Vol. II - Finite Element Method - Jacques-Hervé SAIAC

COMPUTATIONAL METHODS AND ALGORITHMS Vol. II - Finite Element Method - Jacques-Hervé SAIAC COMPUTATIONAL METHODS AND ALGORITHMS Vol. II - Fnte Element Method - Jacques-Hervé SAIAC FINITE ELEMENT METHOD Jacques-Hervé SAIAC Départment de Mathématques, Conservatore Natonal des Arts et Méters, Pars,

More information

A MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS

A MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS Journal of Mathematcs and Statstcs 9 (1): 4-8, 1 ISSN 1549-644 1 Scence Publcatons do:1.844/jmssp.1.4.8 Publshed Onlne 9 (1) 1 (http://www.thescpub.com/jmss.toc) A MODIFIED METHOD FOR SOLVING SYSTEM OF

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

More information

Adaptive Kernel Estimation of the Conditional Quantiles

Adaptive Kernel Estimation of the Conditional Quantiles Internatonal Journal of Statstcs and Probablty; Vol. 5, No. ; 206 ISSN 927-7032 E-ISSN 927-7040 Publsed by Canadan Center of Scence and Educaton Adaptve Kernel Estmaton of te Condtonal Quantles Rad B.

More information

EEE 241: Linear Systems

EEE 241: Linear Systems EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Procedia Computer Science

Procedia Computer Science Avalable onlne at www.scencedrect.com Proceda Proceda Computer Computer Scence Scence 1 (01) 00 (009) 589 597 000 000 Proceda Computer Scence www.elsever.com/locate/proceda Internatonal Conference on Computatonal

More information

Formal solvers of the RT equation

Formal solvers of the RT equation Formal solvers of the RT equaton Formal RT solvers Runge- Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan

More information

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method Journal of Electromagnetc Analyss and Applcatons, 04, 6, 0-08 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse

More information

Research Article Numerov s Method for a Class of Nonlinear Multipoint Boundary Value Problems

Research Article Numerov s Method for a Class of Nonlinear Multipoint Boundary Value Problems Hndaw Publsng Corporaton Matematcal Problems n Engneerng Volume 2012, Artcle ID 316852, 29 pages do:10.1155/2012/316852 Researc Artcle Numerov s Metod for a Class of Nonlnear Multpont Boundary Value Problems

More information

Global Sensitivity. Tuesday 20 th February, 2018

Global Sensitivity. Tuesday 20 th February, 2018 Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

Finite Element Modelling of truss/cable structures

Finite Element Modelling of truss/cable structures Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures

More information

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(

More information

Adjoint Methods of Sensitivity Analysis for Lyapunov Equation. Boping Wang 1, Kun Yan 2. University of Technology, Dalian , P. R.

Adjoint Methods of Sensitivity Analysis for Lyapunov Equation. Boping Wang 1, Kun Yan 2. University of Technology, Dalian , P. R. th World Congress on Structural and Multdscplnary Optmsaton 7 th - th, June 5, Sydney Australa Adjont Methods of Senstvty Analyss for Lyapunov Equaton Bopng Wang, Kun Yan Department of Mechancal and Aerospace

More information

A Local Variational Problem of Second Order for a Class of Optimal Control Problems with Nonsmooth Objective Function

A Local Variational Problem of Second Order for a Class of Optimal Control Problems with Nonsmooth Objective Function A Local Varatonal Problem of Second Order for a Class of Optmal Control Problems wth Nonsmooth Objectve Functon Alexander P. Afanasev Insttute for Informaton Transmsson Problems, Russan Academy of Scences,

More information

Application of Pontryagin s Maximum Principles and Runge- Kutta Methods in Optimal Control Problems

Application of Pontryagin s Maximum Principles and Runge- Kutta Methods in Optimal Control Problems IOSR Journal of Matematcs (IOSR-JM) e-issn: 78-578, p-issn: 39-765X. Volume, Issue 5 Ver. II (Sep. - Oct. 5), PP 43-63 www.osrjournals.org Applcaton of Pontryagn s Maxmum Prncples and Runge- Kutta Metods

More information

Nice plotting of proteins II

Nice plotting of proteins II Nce plottng of protens II Fnal remark regardng effcency: It s possble to wrte the Newton representaton n a way that can be computed effcently, usng smlar bracketng that we made for the frst representaton

More information

ACTM State Calculus Competition Saturday April 30, 2011

ACTM State Calculus Competition Saturday April 30, 2011 ACTM State Calculus Competton Saturday Aprl 30, 2011 ACTM State Calculus Competton Sprng 2011 Page 1 Instructons: For questons 1 through 25, mark the best answer choce on the answer sheet provde Afterward

More information

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,

More information

, rst we solve te PDE's L ad L ad n g g (x) = ; = ; ; ; n () (x) = () Ten, we nd te uncton (x), te lnearzng eedbac and coordnates transormaton are gve

, rst we solve te PDE's L ad L ad n g g (x) = ; = ; ; ; n () (x) = () Ten, we nd te uncton (x), te lnearzng eedbac and coordnates transormaton are gve Freedom n Coordnates Transormaton or Exact Lnearzaton and ts Applcaton to Transent Beavor Improvement Kenj Fujmoto and Tosaru Suge Dvson o Appled Systems Scence, Kyoto Unversty, Uj, Kyoto, Japan suge@robotuassyoto-uacjp

More information

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0 MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector

More information

Chapter Newton s Method

Chapter Newton s Method Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve

More information

Some modelling aspects for the Matlab implementation of MMA

Some modelling aspects for the Matlab implementation of MMA Some modellng aspects for the Matlab mplementaton of MMA Krster Svanberg krlle@math.kth.se Optmzaton and Systems Theory Department of Mathematcs KTH, SE 10044 Stockholm September 2004 1. Consdered optmzaton

More information