Numerical Solution of two dimensional coupled viscous Burgers Equation using the Modified Cubic B-Spline Differential Quadrature Method
|
|
- Tyler Moody
- 5 years ago
- Views:
Transcription
1 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 athematcs & Statstcs, DDU Gorakhpur Unversty, Gorakhpur-73009, Inda ISRO Telemetry, Trackng and Command etwork (ISTRAC), Bangalore , Inda 3 ISRO Satellte Center (ISAC), Bangalore , Inda ABSTRACT In ths paper, a numercal soluton of the two dmensonal nonlnear coupled vscous Burgers equaton s dscussed wth the approprate ntal and boundary condtons usng the modfed cubc B-splne dfferental quadrature method (CB-DQ). In ths method, the weghtng coeffcents are computed usng the modfed cubc B-splne as a bass functon n the dfferental quadrature method. Thus, the coupled Burgers equatons are reduced nto a system of ordnary dfferental equatons (ODEs). An optmal fve stage and fourth-order strong stablty preservng Runge Kutta (SSP-RK54) scheme s appled to solve the resultng system of ODEs. The accuracy of the scheme s llustrated va two numercal examples. Computed results are compared wth the exact solutons and other results avalable n the lterature. umercal results show that the CB-DQ s effcent and relable scheme for solvng the two dmensonal coupled Burgers equaton. Keywords: Burgers equaton; Reyonlds number; odfed cubc B-splne functon; CB- DQ; SSP-RK Introducton Consder the two dmensonal nonlnear unsteady coupled vscous Burgers equatons: u 1 u u v u u u, (1.1) t x y Re x y u v t x y Re x y v v v 1 v v wth the ntal condtons:, (1.)
2 ,,0 1, ;,,,0, ;, u x y x y x y, v x y x y x y (1.3) and Drchlet boundary condtons: u x, y, t x, y, t; x, y vx, y, t x, y, t; x, y where, :, u x, y, t and,,, t 0 (1.4) x y a x b c x d s the computatonal doman and s ts boundary, v x y t are the velocty components to be determned, 1,, and are u u known functons, s unsteady term, u t s the nonlnear convecton term, 1 u u x Re x y s the dffuson term, and Re s the Reynolds number. Coupled vscous Burgers equaton s a more approprate form of the aver-stokes equatons havng the exact solutons. It has the same convectve and dffuson form as the ncompressble aver-stokes equatons, and s a smple model for understandng of varous physcal flows and problems, such as hydrodynamc turbulence, shock wave theory, wave processes n thermo-elastc medum, vortcty transport, dsperson n porous meda 1-3. umercal soluton of Burgers equaton s a natural and frst step towards developng methods for the computaton of complex flows. Thus, t has become customary to test new approaches n computatonal flud dynamcs by mplementng novel and new approaches to the Burgers equaton yeldng n varous fnte-dfferences, fnte volume, fnte-element and boundary element methods etc. Analytc soluton of two dmensonal coupled Burgers equatons was frst gven by Fletcher 4 usng the Hopf-Cole transformaton. The numercal soluton of coupled Burgers equatons are numercally solved by many researchers In recent years, varous researchers 16- proposed varant of dfferental quadrature method for the numercal soluton 3, 4 of the one and two dmensonal lnear/nonlnear dfferental equatons. Korkmaz & Dag proposed cubc B-splne and snc dfferental quadrature methods. Arora & Sngh 5 proposed the modfed cubc B-splne dfferental quadrature method (CB-DQ) and appled on one dmensonal Burgers equaton to checked ts effcency and accuracy. They found that CB- DQ s very powerful and effcent scheme as compared to other exstng numercal methods. Recently, an extenton of the CB-DQ s proposed by Jwar & Yuan 6 to show the computatonal modelng of two-dmensonal reacton dffuson Brusselator model wth approprate ntal and eumann boundary condtons. Researchers 7-9 developed an optmal
3 strong stablty preservng (SSP) hgh order tme dscretzaton schemes. Strong stablty propertes of SSP methods s preserved n any norm, sem norm or convex functonal of the spatal dscretzaton coupled wth the frst order Euler tme steppng. A descrpton of the optmal explct and mplct SSP Runge-Kutta and multstep methods s also dscussed by the authors. In ths paper, we study the numercal smulaton of the two-dmensonal unsteady nonlnear coupled vscous Burgers equatons for dfferent Reynolds number. The effcacy and accuracy of the method s confrmed by takng two test problem wth sutable ntal and boundary condtons. Ths study shows that the CB-DQ results are acceptable and n good agreement wth the exact solutons and earler results avalable n the lterature. Rest of the artcle s prepared as: In Secton, the modfed cubc B-splne dfferental quadrature method s ntroduced. In Secton 3, the mplementaton procedure for the problem (1.1) (1.) together wth the ntal condtons (1.3) and boundary condtons (1.4) s llustrated; In Secton 4, two test problems are gven to establsh the applcablty and accuracy of the method, whle the Secton 5 concludes our study.. odfed cubc B-splne dfferental quadrature method In 197, Bellman et al. 16 ntroduced dfferental quadrature method (DQ). Ths method approxmates the spatal dervatves of a functon usng the weghted sum of the functonal values at the certan dscrete ponts. In DQ, the weghtng coeffcents are determned usng several knds of test functons such as splne functon 3, snc functon 4, Lagrange nterpolaton polynomals, Legendre polynomals 17- etc. Ths secton revsts the CB- DQ 5-6 n order to complete our problem n two dmenson. It s assumed that the and grd ponts: a x,... 1 x x b and c y,... 1 y y d are unformly dstrbuted wth the spatal step sze x x 1 x and y y 1 y n x and y drectons, respectvely. (keepng The frst and second order spatal partal dervatves of u x, y, t wth respect to x y as fxed) and wth respect to y (keepng x as fxed), approxmated at x and y, respectvely, are defned as:,, 1 k k u x y t x k1 w u x, y, t, 1,,..., (.1)
4 ,, k k u x y t x k1,, 1 k u x y t y k1 w u x, y, t, 1,,..., w u x, y, t, 1,,..., k (.) (.3),, k u x y t y k1 w u x, y, t, 1,,...,. k (.4) In the same way, the frst and second order spatal partal dervatves of vx, y, t wth respect to x and wth respect to y are approxmated as:,, 1 k k v x y t x k1 w v x, y, t, 1,,...,,, k k v x y t x k1,, 1 k v x y t y k1,, k v x y t y k1 w v x, y, t, 1,,..., w v x, y, t, 1,,..., k w v x, y, t, 1,,...,, k (.5) (.6) (.7) (.8) where w r and r w, r 1, are the weghtng coeffcents of the rth-order spatal partal dervatves wth respect to x and y. The cubc B-splne bass functons at the knots are defned as: m x x 3 xm, x xm, xm 1 x 3 3 xm x xm 1 x xm1 xm 3,, 1 4,, 0, otherwse, 3 h xm x x xm 1 xm (.9) where,,...,, s chosen n such a way that t forms a bass over the doman x, y : a x b; c y d. The values of cubc B-splnes and ts dervatves at the nodal ponts are depcted n Table 1.
5 Table 1: Coeffcents of the cubc B-splne m and ts dervatves at the node x m. m x xm xm 1 x m xm 1 xm m x 0 3/h 0 3/h 0 m x 0 6/h 1/ h 6/h 0 ow, to get a dagonally domnant system of the lnear equatons, the cubc B-splne bass functons are modfed as 5 : x x x x x x m 1 x 1 x 1 x x x x m x x for m 3,..., -, 1 where forms a bass over the doman,,..., 1 m x, y : a x b; c y d. (.10) In Eq. (.1), substtutng the values of ( x ), m 1,,...,, we get a system of lnear equatons: (1) x w x, 1,,...,. (.11) m k m k k1 Wth the help of Eq. (.9), (.10) and Table 1, Eq. (.11) reduces nto a trdagonal system of equatons: 1 Aw R, for 1,,...,, where w w, w,..., w 1 T (.1) s the weghtng coeffcent vector correspondng to x, 1,,,,......, 1,,, T R, and the coeffcent matrx A s gven by: 1,1 1,,1,,3 3, 3,3 3,4 A, 3,, 1 1, 1, 1 1,, 1,
6 Here, we pont out that the coeffcent matrx A s nvertble. The trdagonal system of lnear equatons (.1) s solved for each usng Thomas algorthm, whch gves the weghtng coeffcents w 1 k of the frst order partal dervatve. In a smlar way, the weghtng coeffcents w,1, are determned. Weghtng coeffcents w,1,, can be computed as : r1 r 1 r1 w w r w w, for and 1,,3,..., ; r,3,..., 1 x x r r w w, for, 1, (.13) where r 1 w and r w are the weghtng coeffcents of the th r 1 and th r order partal dervatves wth respect to x. Smlarly, the weghtng coeffcents w 1 k of the frst order partal dervatves wth respect to y usng the modfed cubc B-Splne functons n Eq. (11) s obtaned. Weghtng coeffcents w,1, for the second dervatves can be computed from the formule: r1 r 1 r1 w w r w w, for and 1,,3,..., ; r,3,..., 1 x x r r w w, for, 1, (.14) where r 1 w and r w are the weghtng coeffcents of the th r 1 and th r order partal dervatves wth respect to y. 3. CB-DQ for two-dmensonal coupled Burgers equaton On substtutng the approxmate values of the spatal dervatves computed by the CB- DQ, Eq. (1.1) can be wrtten as: u x, y, t t (1) (1),,,, k k k k u x y w u x y v x y w u x y k1 k1 1 () (),,,,, 0, 1,,...,, 1,,...,. k k k k Re w u x y w u x y x y R t k1 k1 (3.1) Smlarly, Eq. () can be wrtten as:
7 u x, y, t t (1) (1),,,, k k k k u x y w v x y v x y w v x y k1 k1 1 () (),,,,, 0, 1,,...,, 1,,...,. k k k k Re w v x y w v x y x y R t k1 k1 (3.) Eq. (3.1) and Eq. (3.) reduce nto a system of ordnary dfferental equatons: du x, y, t dt dv x, y, t dt 1 F u x, y, t, 1,,..., and 1,,...,. (3.3) F u x, y, t, 1,,..., and 1,,...,. (3.4) Eqs. (3.3) and (3.4) together wth the ntal condtons (1.3) and Drchlet boundary condtons (1.4) are solved by SSP-RK54 scheme. 4. Results and dscusson Here we consder two test problems of two dmensonal coupled Burgers equaton as gven n the ntroducton part to provde the CB-DQ numercal solutons. The accuracy and consstency of the scheme s measured n terms of error norms L and L, defned as: n n exact computed : exact computed,, 0 0 exact computed uexact ucomputed u, u,, L u u u u L := max where uexact and u computed represent the exact and computed solutons, respectvely., (4.1) problems. umercal solutons to Eqs. (1) and () wll be tested for the followng two test 4.1. Problem: The analytcal solutons of Eqs. (1.1) and (1.) can be generated as 4 : 3 1 u x, y, t 4 4 1exp4x 4y tre/ 3 x,y, 3 1 v x, y, t exp4x 4y tre/ 3 (4.)
8 The square doman x, y: 0 x 1,0 y 1 and the ntal and boundary condtons for,, s consdered as the computatonal doman, u x y t and vx, y, tare taken from the analytcal solutons Eq. (4.). 4.. Problem: In ths problem, we take the computatonal doman x, y :0 x 0.5,0 y 0.5 wth the ntal condtons:,,0 sn cos,,0 u x y x y x,y, v x y x y (4.3) and Drchet boundary condtons: u u v v u u v v 0, y, t cos y, 0.5, y, t y 0, yt, y, 0.5, yt, 0.5 y, 1 cos, x,0, t x x,0.5, t sn x, x,0, t x, x,0.5, t x0.5, 1 sn, 0 y 0.5, t 0, 0 x 0.5, t 0, (4.4) (4.5) For the test problem 4.1, we have taken a grd sze 0 0 wth tme step t and Re 100. Computed and exact values of the results gven by Srvastava et al u and v are shown n Tables and 3 along wth and Bahadr 9 at some typcal grd pont. The tabulated results show that the proposed scheme produces better result than Bahadr 9. Tables 4 and 5 show the errors L and L, and also the rate of convergence of u and v components, respectvely, at Re 100, t 1.0 for t From Tables 4 and 5, t can be observed that the CB-DQ performs better than Srvastava et al. 15 and gves more than quadratc rate of convergence. The CB-DQ computed solutons of u and v for Re 100 att 0.5 are depcted n Fg. 1 whle Fg. shows analytcal solutons of u and v, respectvely. For the test problem 4., numercal computatons are carred out wth the parameters: Re 50, 100, and 0 0 grd, t 0.65 for t n order to compare the computed results wth those gven by Jan & Holla 5, Bahadr 9, and Srvastava et al Table 6 shows
9 the comparsons of numercal results obtaned usng the CB-DQ scheme at t 0.65, wth the methods of Jan & Holla 5, Bahadr 9, and Srvastava et al From Table 6, t can be notced that the computed CB-DQ results are n good agreement wth Jan & Holla 5, Bahadr 9, and Srvastava et al Fg. 3 depcts the CB-DQ computed u and v solutons correspondng to Re =50, 100 and 500 at t Table : Comparson of CB-DQ and exact solutons of u for Re 100, 0 0 t grd and Grd xy, CB- DQ t 0.5 t.0 Exact I-LFD 14 Expo- Bahadr 9 CB- FD 15 DQ Exact I- LFD 14 Expo- Bahadr 9 FD 15 (0.1,0.1) (0.5,0.1) (0.9,0.1) (0.3, 0.3) (0.7, 0.3) (0.1, 0.5) (0.5, 0.5) (0.9, 0.5) (0.3, 0.7) (0.7, 0.7) (0.1, 0.9) (0.5, 0.9) (0.9, 0.9) Grd xy, Table 3: Comparson of CB-DQ and exact solutons of v for Re 100, 0 0 grd and t t 0.5 t.0 CB-DQ Exact I- Expo- Bahadr 9 CB-DQ Exact I- Expo- Bahadr 9 LFD 14 FD 15 LFD 14 FD 15 (0.1,0.1) (0.5,0.1) (0.9,0.1) (0.3, 0.3) (0.7, 0.3) (0.1, 0.5) (0.5, 0.5) (0.9, 0.5)
10 (0.3, 0.7) (0.7, 0.7) (0.1, 0.9) (0.5, 0.9) (0.9, 0.9) Table 4. L, L errors and rate of convergence for the u -component for Re 100, t at t 1.0. Grd L Expo-FD 15 CB-DQ Expo-FD 15 CB-DQ L e e e e-003 Rate e e e e e e e e e e e-003.0e Rate e e Table 5. L, L errors and rate of convergence for the v -component for Re 100, t at t 1.0. Grd L Expo-FD 15 CB-DQ Expo-FD 15 CB-DQ L Rate e e e e e e e e e e e e e e e-003.0e Rate e e Table 6: Comparson of computed results of u and v for Re 50, grd sze 0 0 and t at t Grd Computed values of u Computed values of v xy, CB- I- Expo- DQ LFD 14 FD 15 Bahadr 9 Jan CB- I- Expoand DQ LFD 14 FD 15 Bahadr 9 Jan and Holla 5 Holla 5 (0.1, 0.1) (0.3, 0.1) (0., 0.) (0.4, 0.)
11 (0.1, 0.3) (0.3, 0.3) (0., 0.4) (0.4, 0.4) Fg. 1: umercal soluton at t 0.5 wth t , Re 100 and grd sze 0 0. for the test problem 4.1. Fg.. Exact soluton at t 0.5 wth t , Re 100 and grd sze 0 0. for the test problem 4.1.
12 Fg. 3. umercal solutons at t 0.65 wth t and grd sze problem for the 5. Conclusons A modfed cubc B-splne dfferental quadrature method s presented for the numercal solutons of two dmensonal nonlnear coupled vscous Burgers equatons. The computed
13 results show that the soluton obtaned by ths scheme s hghly accurate and very close to the exact solutons. We also notce that the scheme has more than quadratc rate of convergance. The obtaned results show that the CB-DQ s a promsng numercal scheme for solvng the hgher dmensonal nonlnear physcal problems governed by partal dfferental equatons. References 1 Cole J. D., On a quaslnear parabolc equatons occurrng n aerodynamcs, Quart Appl ath 9, 5 (1951). Espov S. E., Coupled Burgers equatons: a model of poly-dspersve sedmentaton, Phys Rev 5, 3711 (1995). 3 J. D. Logan, An ntroducton to nonlnear partal dfferental equatons, Wly-Interscence, ew York (1994). 4 C. A. J. Fletcher, Generatng exact solutons of the two-dmensonal Burgers equaton, Int. J. umer. eth. Fluds 3, 13 (1983). 5 P. C. Jan, D.. Holla, umercal soluton of coupled Burgers equatons, Int. J. umer. eth.eng. 1, 13 (1978). 6 C. A. J. Fletcher, A comparson of fnte element and fnte dfference of the one- and twodmensonal Burgers equatons, J. Comput. Phys. 51, 159 (1983). 7 F.W. Wubs, E.D. de Goede, An explct mplct method for a class of tme-dependent partal dfferental equatons, Appl. umer. ath. 9,157 (199). 8 O. Goyon, ultlevel schemes for solvng unsteady equatons, Int. J. umer. eth. Fluds, 937 (1996). 9 Bahadr A. R., A fully mplct fnte-dfference scheme for two-dmensonal Burgers equaton, Appled athematcs and Computaton 137, 131 (003). 10 V. K. Srvastava,. Tamsr, U. Bhardwa, Y. Sanyasrau, Crank-colson scheme for numercal solutons of two dmensonal coupled Burgers equatons, IJSER (5), 44 (011). 11. Tamsr, V. K. Srvastava, A sem-mplct fnte-dfference approach for twodmensonal coupled Burgers equatons, IJSER (6), 46 (011). 1 V. K. Srvastava,. Tamsr, Crank-colson sem-mplct approach for numercal solutons of two-dmensonal coupled nonlnear Burgers equatons, Int. J. Appl. ech. Eng. 17 (), 571 (01).
14 13 V. K. Srvastava,. K. Awasth,. Tamsr, A fully mplct Fnte-dfference soluton to one-dmensonal Coupled onlnear Burgers equatons, Int. J. ath. Sc. 7(4), 3 (013). 14 V. K. Srvastava,. K. Awasth, S. Sngh, An mplct logarthm fnte dfference technque for two dmensonal coupled vscous Burgers equaton, AIP Advances 3, 1105 (013). 15 V. K. Srvastava, S. Sngh,. K. Awasth, umercal solutons of coupled Burgers equatons by an mplct fnte-dfference scheme, AIP Advances 3, (013). 16 R. Bellman, B. G. Kashef, J. Cast, Dfferental quadrature: a technque for the rapd soluton of nonlnear dfferental equatons, J. Comput. Phy. 10, 40-5, (197). 17 C. Shu, B. E. Rchards, Applcaton of generalzed dfferental quadrature to solve two dmensonal ncompressble naver-stokes equatons, Int. J. umer. eth. Fluds, 15, , (199). 18 J. R. Quan, C.T. Chang, ew nsghts n solvng dstrbuted system equatons by the quadrature methods-i, Comput. Chem. Eng. 13, , (1989). 19 J. R. Quan, C.T. Chang, ew nsghts n solvng dstrbuted system equatons by the quadrature methods-ii, Comput. Chem. Eng. 13, , (1989). 0 C. Shu, Y.T. Chew, Fourer expanson-based dfferental quadrature and ts applcaton to Helmholtz egenvalue problems, Commun. umer. ethods Eng. 13 (8), , (1997). 1 C. Shu, H. Xue, Explct computaton of weghtng coeffcents n the harmonc dfferental quadrature, J. Sound Vb. 04 (3), , (1997). C. Shu, Dfferental Quadrature and ts Applcaton n Engneerng, Athenaeum Press Ltd., Great Brtan, (000). 3 A. Korkmaz, I. Dag, Cubc B-splne dfferental quadrature methods and stablty for Burgers equaton, Eng. Comput. Int. J. Comput. Aded Eng. Software, 30 (3), , (013). 4 A. Korkmaz, I. Dag, Shock wave smulatons usng snc dfferental quadrature method, Eng. Comput. Int. J. Comput. Aded Eng. Software, 8(6), , (011). 5 G. Arora, B. K. Sngh, umercal soluton of Burgers equaton wth modfed cubc B- splne dfferental quadrature method, Appled ath. Comput., 4 (1), , (013). 6 R. Jwar, J. Yuan, A computatonal modelng of two dmensonal reacton dffuson Brusselator system arsng n chemcal processes, J. ath. Chem., 5 (6), , (014). 7 S. Gottleb, C. W. Shu, E. Tadmor, Strong Stablty-Preservng Hgh-Order Tme Dscretzaton ethods, SIA REVIEW, 43(1), 89-11, (001).
15 8 J. R. Spter, S. J. Ruuth, A new class of optmal hgh-order strong stablty-preservng tmesteppng schemes, SIA J. umer. Analyss. 40 (), , (00). 9 S. Gottleb, D. I. Ketcheson, C. W. Shu, Hgh Order Strong Stablty Preservng Tme Dscretzatons, J. Sc. Comput., 38, 51-89, (009).
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 informationNON-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 informationCHAPTER 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 informationAppendix 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 informationPOLYNOMIAL 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 informationLecture 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 informationHigh 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 informationTHE 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 informationNumerical 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 informationResearch 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(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate
Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1
More informationA 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 informationSolution 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 informationNumerical Simulation of One-Dimensional Wave Equation by Non-Polynomial Quintic Spline
IOSR Journal of Matematcs (IOSR-JM) e-issn: 78-578, p-issn: 319-765X. Volume 14, Issue 6 Ver. I (Nov - Dec 018), PP 6-30 www.osrournals.org Numercal Smulaton of One-Dmensonal Wave Equaton by Non-Polynomal
More informationConvexity 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 informationOne-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 informationA 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 informationLecture 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 informationIrregular 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 informationNew 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 informationME 501A Seminar in Engineering Analysis Page 1
umercal Solutons of oundary-value Problems n Os ovember 7, 7 umercal Solutons of oundary- Value Problems n Os Larry aretto Mechancal ngneerng 5 Semnar n ngneerng nalyss ovember 7, 7 Outlne Revew stff equaton
More informationALGORITHM 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 informationpage 2 2 dscretzaton mantans ths stablty under a sutable restrcton on the tme step. SSP tme dscretzaton methods were frst developed by Shu n [20] and
page 1 A Survey of Strong Stablty Preservng Hgh Order Tme Dscretzatons Ch-Wang Shu Λ 1 Introducton Numercal soluton for ordnary dfferental equatons (ODEs) s an establshed research area. There are many
More informationNumerical 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 informationFormal 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 informationOptimal Control of Temperature in Fluid Flow
Kawahara Lab. 5 March. 27 Optmal Control of Temperature n Flud Flow Dasuke YAMAZAKI Department of Cvl Engneerng, Chuo Unversty Kasuga -3-27, Bunkyou-ku, Tokyo 2-855, Japan E-mal : d33422@educ.kc.chuo-u.ac.jp
More informationSeptic 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 informationChapter 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 informationMaejo 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[7] R.S. Varga, Matrix Iterative Analysis, Prentice-Hall, Englewood Clis, New Jersey, (1962).
[7] R.S. Varga, Matrx Iteratve Analyss, Prentce-Hall, Englewood ls, New Jersey, (962). [8] J. Zhang, Multgrd soluton of the convecton-duson equaton wth large Reynolds number, n Prelmnary Proceedngs of
More informationHigher Order Wall Boundary Conditions for Incompressible Flow Simulations
THE 5 TH ASIAN COMPUTAITIONAL FLUID DYNAMICS BUSAN KOREA OCTOBER 7-30 003 Hgher Order Wall Boundary Condtons for Incompressble Flow Smulatons Hdetosh Nshda. Department of Mechancal and System Engneerng
More informationCubic 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 informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationTensor Smooth Length for SPH Modelling of High Speed Impact
Tensor Smooth Length for SPH Modellng of Hgh Speed Impact Roman Cherepanov and Alexander Gerasmov Insttute of Appled mathematcs and mechancs, Tomsk State Unversty 634050, Lenna av. 36, Tomsk, Russa RCherepanov82@gmal.com,Ger@npmm.tsu.ru
More informationA Solution of the Harry-Dym Equation Using Lattice-Boltzmannn and a Solitary Wave Methods
Appled Mathematcal Scences, Vol. 11, 2017, no. 52, 2579-2586 HIKARI Ltd, www.m-hkar.com https://do.org/10.12988/ams.2017.79280 A Soluton of the Harry-Dym Equaton Usng Lattce-Boltzmannn and a Soltary Wave
More informationChapter 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 informationLecture 21: Numerical methods for pricing American type derivatives
Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)
More information6.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 informationLecture 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 informationTransfer 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 informationProcedia 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 informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More informationSTUDY 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 informationConsistency & 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 informationLecture 5.8 Flux Vector Splitting
Lecture 5.8 Flux Vector Splttng 1 Flux Vector Splttng The vector E n (5.7.) can be rewrtten as E = AU (5.8.1) (wth A as gven n (5.7.4) or (5.7.6) ) whenever, the equaton of state s of the separable form
More informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More informationNumerical Solution of Ordinary Differential Equations
Numercal Methods (CENG 00) CHAPTER-VI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples
More informationDETERMINATION 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 informationChapter 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 informationSignificance of Dirichlet Series Solution for a Boundary Value Problem
IOSR Journal of Engneerng (IOSRJEN) ISSN (e): 5-3 ISSN (p): 78-879 Vol. 6 Issue 6(June. 6) V PP 8-6 www.osrjen.org Sgnfcance of Drchlet Seres Soluton for a Boundary Value Problem Achala L. Nargund* and
More informationSolving 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 information2.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 informationFTCS Solution to the Heat Equation
FTCS Soluton to the Heat Equaton ME 448/548 Notes Gerald Recktenwald Portland State Unversty Department of Mechancal Engneerng gerry@pdx.edu ME 448/548: FTCS Soluton to the Heat Equaton Overvew 1. Use
More informationNUMERICAL 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 informationResearch Article A Multilevel Finite Difference Scheme for One-Dimensional Burgers Equation Derived from the Lattice Boltzmann Method
Appled Mathematcs Volume 01, Artcle ID 9590, 13 pages do:10.1155/01/9590 Research Artcle A Multlevel Fnte Dfference Scheme for One-Dmensonal Burgers Equaton Derved from the Lattce Boltzmann Method Qaoe
More informationAdditional 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 information1-Dimensional Advection-Diffusion Finite Difference Model Due to a Flow under Propagating Solitary Wave
014 4th Internatonal Conference on Future nvronment and nergy IPCB vol.61 (014) (014) IACSIT Press, Sngapore I: 10.776/IPCB. 014. V61. 6 1-mensonal Advecton-ffuson Fnte fference Model ue to a Flow under
More informationModified Mass Matrices and Positivity Preservation for Hyperbolic and Parabolic PDEs
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2000; 00:6 Prepared usng cnmauth.cls [Verson: 2000/03/22 v.0] Modfed Mass Matrces and Postvty Preservaton for Hyperbolc and
More informationNumerical Solutions of a Generalized Nth Order Boundary Value Problems Using Power Series Approximation Method
Appled Mathematcs, 6, 7, 5-4 Publshed Onlne Jul 6 n ScRes. http://www.scrp.org/journal/am http://.do.org/.436/am.6.77 umercal Solutons of a Generalzed th Order Boundar Value Problems Usng Power Seres Approxmaton
More informationCOMPARISON 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 informationA NUMERICAL COMPARISON OF LANGRANGE AND KANE S METHODS OF AN ARM SEGMENT
Internatonal Conference Mathematcal and Computatonal ology 0 Internatonal Journal of Modern Physcs: Conference Seres Vol. 9 0 68 75 World Scentfc Publshng Company DOI: 0.4/S009450059 A NUMERICAL COMPARISON
More informationConstruction of Serendipity Shape Functions by Geometrical Probability
J. Basc. Appl. Sc. Res., ()56-56, 0 0, TextRoad Publcaton ISS 00-0 Journal of Basc and Appled Scentfc Research www.textroad.com Constructon of Serendpty Shape Functons by Geometrcal Probablty Kamal Al-Dawoud
More informationDifference 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 informationTHE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD
Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS
More informationThe 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 informationCOMPOSITE 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 informationGlobal 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 informationConsideration of 2D Unsteady Boundary Layer Over Oscillating Flat Plate
Proceedngs of the th WSEAS Internatonal Conference on Flud Mechancs and Aerodynamcs, Elounda, Greece, August -, (pp-) Consderaton of D Unsteady Boundary Layer Over Oscllatng Flat Plate N.M. NOURI, H.R.
More information2.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 informationSuppose 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 informationHaar 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 informationSolving 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 informationReport 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 informationLeast squares cubic splines without B-splines S.K. Lucas
Least squares cubc splnes wthout B-splnes S.K. Lucas School of Mathematcs and Statstcs, Unversty of South Australa, Mawson Lakes SA 595 e-mal: stephen.lucas@unsa.edu.au Submtted to the Gazette of the Australan
More informationSolution of the Navier-Stokes Equations
Numercal Flud Mechancs Fall 2011 Lecture 25 REVIEW Lecture 24: Soluton of the Naver-Stokes Equatons Dscretzaton of the convectve and vscous terms Dscretzaton of the pressure term Conservaton prncples Momentum
More informationLecture 10 Support Vector Machines II
Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed
More informationPART 8. Partial Differential Equations PDEs
he Islamc Unverst of Gaza Facult of Engneerng Cvl Engneerng Department Numercal Analss ECIV 3306 PAR 8 Partal Dfferental Equatons PDEs Chapter 9; Fnte Dfference: Ellptc Equatons Assocate Prof. Mazen Abualtaef
More informationA PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.
Proceedngs of the th Brazlan Congress of Thermal Scences and Engneerng -- ENCIT 006 Braz. Soc. of Mechancal Scences and Engneerng -- ABCM, Curtba, Brazl,- Dec. 5-8, 006 A PROCEDURE FOR SIMULATING THE NONLINEAR
More information1 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 informationPHYS 705: Classical Mechanics. Canonical Transformation II
1 PHYS 705: Classcal Mechancs Canoncal Transformaton II Example: Harmonc Oscllator f ( x) x m 0 x U( x) x mx x LT U m Defne or L p p mx x x m mx x H px L px p m p x m m H p 1 x m p m 1 m H x p m x m m
More informationSalmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2
Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to
More informationA 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 informationHongyi 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 informationThe equation of motion of a dynamical system is given by a set of differential equations. That is (1)
Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence
More informationDiscretization. Consistency. Exact Solution Convergence x, t --> 0. Figure 5.1: Relation between consistency, stability, and convergence.
Chapter 5 Theory The numercal smulaton of PDE s requres careful consderatons of propertes of the approxmate soluton. A necessary condton for the scheme used to model a physcal problem s the consstency
More informationFinite 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 information2 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 informationThe 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 informationNote 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 informationCHAPTER - 3 Orthogonal Collocation of Finite Elements using Hermite Basis
CHAPTER - 3 Orthogonal Collocaton of Fnte Elements usng Hermte Bass Contents of ths chapter are publshed n: 1. Computers and Chemcal Engneerng, 58, 03-10, 013.. Appled Mathematcal Scences, 7(34), 1663-1673,
More informationn α 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 informationThe Discretization Process
FMIA F Moukalled L Mangan M Darwsh An Advanced Introducton wth OpenFOAM and Matlab Ths textbook explores both the theoretcal foundaton of the Fnte Volume Method (FVM) and ts applcatons n Computatonal Flud
More informationFinding The Rightmost Eigenvalues of Large Sparse Non-Symmetric Parameterized Eigenvalue Problem
Fndng he Rghtmost Egenvalues of Large Sparse Non-Symmetrc Parameterzed Egenvalue Problem Mnghao Wu AMSC Program mwu@math.umd.edu Advsor: Professor Howard Elman Department of Computer Scences elman@cs.umd.edu
More informationVisco-Rubber Elastic Model for Pressure Sensitive Adhesive
Vsco-Rubber Elastc Model for Pressure Senstve Adhesve Kazuhsa Maeda, Shgenobu Okazawa, Koj Nshgch and Takash Iwamoto Abstract A materal model to descrbe large deformaton of pressure senstve adhesve (PSA
More informationComparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method
Appled Mathematcal Scences, Vol. 7, 0, no. 47, 07-0 HIARI Ltd, www.m-hkar.com Comparson of the Populaton Varance Estmators of -Parameter Exponental Dstrbuton Based on Multple Crtera Decson Makng Method
More informationThe Quadratic Trigonometric Bézier Curve with Single Shape Parameter
J. Basc. Appl. Sc. Res., (3541-546, 01 01, TextRoad Publcaton ISSN 090-4304 Journal of Basc and Appled Scentfc Research www.textroad.com The Quadratc Trgonometrc Bézer Curve wth Sngle Shape Parameter Uzma
More informationPrinciples of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Heat Transfer
Prncples of Food and Boprocess Engneerng (FS 31) Solutons to Example Problems on Heat Transfer 1. We start wth Fourer s law of heat conducton: Q = k A ( T/ x) Rearrangng, we get: Q/A = k ( T/ x) Here,
More informationMMA 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 informationA 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 informationCOEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN
Int. J. Chem. Sc.: (4), 04, 645654 ISSN 097768X www.sadgurupublcatons.com COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN R. GOVINDARASU a, R. PARTHIBAN a and P. K. BHABA b* a Department
More information