Meshless analysis of three-dimensional steady-state heat conduction problems
|
|
- Lizbeth Barnett
- 5 years ago
- Views:
Transcription
1 Cin. Pys. B ol. 19, No. 9 (010) Mesless analysis of tree-dimensional steady-state eat conduction problems Ceng Rong-Jun( 程荣军 ) a) and Ge Hong-Xia( 葛红霞 ) b) a) Ningbo Institute of ecnology, Zejiang University, Ningbo , Cina b) Faculty of Science, Ningbo University, Ningbo 31511, Cina (Received November 009; revised manuscript received 16 January 010) Steady-state eat conduction problems arisen in connection wit various pysical and engineering problems te functions satisfy a given partial differential equation and particular boundary conditions, ave attracted muc attention and researc recently. ese problems are independent of time and involve only space coordinates, as in Poisson s equation or te Laplace equation wit Diriclet, Neuman, or mixed conditions. Wen te problems are too complex, it is difficult to find an analytical solution, te only coice left is an approximate numerical solution. is paper deals wit te numerical solution of tree-dimensional steady-state eat conduction problems using te mesless reproducing kernel particle metod (RKPM). A variational metod is used to obtain te discrete equations. e essential boundary conditions are enforced by te penalty metod. e effectiveness of RKPM for tree-dimensional steady-state eat conduction problems is investigated by two numerical examples. Keywords: reproducing kernel particle metod, mesless metod, steady-state eat conduction problem PACC: 000, 060, Introduction Partial differential equations arise in connection wit various pysical and geometrical problems, te involved functions depend on two or more independent variables, usually on time t and on one or several space variables. [1] A steady-state eat conduction equation is one of te most important partial differential equations in engineering matematics, because it occurs in connection wit gravitational fields, electrostatics fields, incompressible fluid flow, and oter areas. [1] Matematically, te problem is a function wic satisfies a given partial differential equation and particular boundary conditions (boundary-value problems). Pysically speaking, te problem is independent of time and involves only space coordinates, just as initial-value problems are associated wit yperbolic partial differential equations of te elliptic type. [,3] Wen steady-state eat conduction problems are too complex, it is difficult to find an analytical solution for suc problems, so te only coice left is an approximate numerical solution. In tis paper, te reproducing kernel particle metod (RKPM) is used to obtain te numerical solution of a tree-dimensional steady-state eat conduction problem. In recent years mesless tecniques ave attracted te attention of researcers. In a mesless (mes-free) metod, a set of scattered nodes are used instead of mesing te domain of te problem. For many engineering problems, suc as large deformation and crack growt, it is necessary to deal wit extremely large deformations or fractures of te mes wit te remesing tecnique. e mesless metod is a new and interesting numerical tecnique wic can solve many engineering problems unsuited for conventional numerical metods wit a minimum of mesing or no mesing at all. [4] Some mesless metods ave been developed, suc as smoot particle ydrodynamics (SPH) metods, [5] radial basis function (RBF), [6] element free Galerkin (EFG) metod, [7 9] mesless local Petrov Galerkin (MLPG) metod, [10] RKPM, [11,1] finite point metod (FPM), [13,14] mesless metod based on local boundary integral equation, [15] boundary element-free metod (BEFM), [16,17] mesless metod wit complex variables, [18 1] and so on. e eat conduction problems arising in many Project supported by te Natural Science Foundation of Ningbo, Cina (Grant Nos. 009A and 009A610154) and te Natural Science Foundation of Zejiang Province, Cina (Grant No. Y ). Corresponding autor. cengrongjun76@16.com c 010 Cinese Pysical Society and IOP Publising Ltd ttp:// ttp://cpb.ipy.ac.cn
2 Cin. Pys. B ol. 19, No. 9 (010) kinds of engineering fields ave tus attracted muc attention and researc recently. It is difficult to find an analytical solution for suc problems, so it is important to find te numerical solution. Cen uses te RBF mesless metod to solve a simple eat conduction problem, [] Sing et al. ave done muc work on steady, unsteady, nonlinear eat transfer problems based on EFG metod, [3 7] Arefmanes et al. put forward a mesless metod for solving eat conduction problem based on MLPG, [8,9] Liu et al. use a mesless least-squares metod for solving te steadystate eat conduction problem, [30,31] Sladek et al. use a mesless local boundary integral equation metod (LBIE) for eat conduction analysis, [3,33] Ceng and Liew use RKPM for two-dimensional unsteady eat conduction problem and obtain te convergence analysis of tis metod. [34] e RKPM is greatly developed and widely applied in various fields. o te best knowledge of te autors, tis metod as not been used for a treedimensional (3D) steady-state eat conduction problem. In tis paper, RKPM is first used to obtain te numerical solution of a 3D steady-state eat conduction problem. In tis metod, te function over te solution domain needs only a set of nodes, it does not require element connectivity. e variational metod is used to obtain te corresponding discrete equations and te essential boundary conditions are enforced by a penalty metod. e numerical results are presented in comparison wit te exact results.. e RKPM [10] In RKPM, we let u (x) = Ω u(x ) w(x x )dx, (1) w(x x ) is te correction kernel function w(x x ) = c(x; x x )w(x x ), () c(x; x x ) is te correction function wic can be te combination of te polynomial basis. Let c(x; x x ) = m p i (x x )b i (x) i=1 = p (x x )b(x), (x Ω), (3) m is te number of terms in te basis, p i (x x ) are te monomial basis functions, and b i (x) are coefficients of te monomial basis functions. Generally, te basis can be cosen as follows: Linear basis: p = (1, x 1 x 1, x x ), (D) (4) p = (1, x 1 x 1, x x, x 3 x 3). (3D) (5) Quadratic basis: p = ( 1, x 1 x 1, (x 1 x 1) ), (1D) (6) p = ( 1, x 1 x 1, x x, (x 1 x 1), (x 1 x 1)(x x ), (x x ) ). (D) (7) Corresponding to te kernel approximation of Eq. (1), we can obtain te discretization approximation of Eq. (1) by te trapezoidal integral rule as follows: u (x) = w(x x I )u(x I ) I = c(x; x x I )w(x x I )u I I, (8) is related to te nodes x I, I =. (9) For a 3D problem, means te total volume of te problem domain. Equation (8) can be rewritten in te vector form u (x) = C(x)W (x) u, (10) w(x x I ) are te weigting functions wit property of compact support, x I are nodes in te domain of te node x, I is te domain measure wic u = (u 1, u,..., u n ), (11)
3 w(x x 1 ) w(x x )... 0 W (x) = = ten Let Cin. Pys. B ol. 19, No. 9 (010) w(x x n ) n, (1). (13) c I (x) = c(x; x x I ), (14) C(x) = (c 1 (x), c (x),..., c n (x)) = b (x)p, (15) p 1 (x x 1 ) p 1 (x x )... p 1 (x x n ) p (x x 1 ) p (x x )... p (x x n ) P =..,.... p m (x x 1 ) p m (x x )... p m (x x n ) (16) b (x) = (b 1 (x), b (x),..., b m (x)), (17) te coefficients b i (x) are determined by te reproducing conditions of te approximation function. In order to obtain te accuracy of n order, we obtain as M(x) = M(x) = M(x)b(x) = H, (18) p(x x I )p (x x I ) w(x x I ) I, (19) H = (1, 0,..., 0), (0) b(x) = M 1 (x)h. (1) e approximation u (x) can ten be expressed u (x) = Φ I (x)u I = Φ(x)u, () te sape function Φ(x) is given by Φ(x) = [Φ 1 (x), Φ (x),..., Φ n (x)] = C(x)W (x). (3) 3. e discretion and numerical implementation A tree-dimensional steady-state eat conduction equation in isotropic materials wit termal properties independent of temperature is governed by ( ) k x + y + z + Q = 0, (on ) (4a) wit boundary conditions: at surface x = 0(S 1 ), = s1, (4b) at surface x = L(S ), = s, (4c) at surface y = 0(S 3 ), at surface y = W (S 4 ), at surface z = 0(S 5 ), at surface z = H(S 6 ), k y n y = ( ), (4d) k y n y = ( ), (4e) k z n z = ( ), (4f) k z n z = ( ), (4g) is te temperature, k denotes te coefficient of termal conductivity and ρ represents te density of te material, c is specific eat of te material, Q denotes internal eat generation per unit volume. is convective eat transfer coefficient, denotes surrounding fluid temperature, is te 3D domain. L, W and H denote te lengt of domain in te direction of x, y and z, respectively. e weigted integral form of Eq. (4a) is obtained as follows [ ( ) ] w k x + y + z + Q d = 0. (5) e weak form of Eq. (5) is [ k w w Q ] ( d wk S x n x + x n y + ) x n z ds = 0. (6) e functional Π ( ) wic denotes te total potential energy can be written as
4 Π ( ) = 1 + Cin. Pys. B ol. 19, No. 9 (010) (k )d S 4 ( ) ds + Qd + S 5 ( ) S 3 ds 3 ( ) ( ) ds + ds. (7) Enforcing essential boundary conditions using a penalty metod, te modified functional Π ( ) can be obtained as Π ( ) = 1 ( ) (k )d Qd + S 3 ds ( ) ( ) ( ) + S 4 ds + S 5 ds + S 6 ds + α ( S1 ) ds + α ( S ) ds. (8) S 1 S aking variation, we ave δπ ( ) = (kδ )d δ Qd + δ ( )ds S 3 S 4 S 5 S 6 + δ α( S1 )ds + δ α( S )ds, (9) S 1 S α is te penalty factor. By Eq. (), we can obtain te approximation of te temperature function (x, t) = S 6 Φ I (x) I (t) = Φ(x), (30) Φ(x) = (Φ 1 (x), Φ (x),..., Φ n (x)), (31) = ( 1 (t), (t),..., n (t)), (3) x = (x, y, z). (33) Substituting Eq. (30) into Eq. (9) yields δπ ( ) = (Φ(x) )(kδ (Φ(x) )d [δ(φ(x) ) Qd + δ(φ(x) )(Φ(x) )ds + δ(φ(x) )α(φ(x) S1 )ds S 3 S 4 S 5 S 6 S 1 + δ(φ(x) )α(φ(x) S )ds. (34) S Since δπ ( ) = 0 and δ is arbitrary in Eq. (34), te following equations can be obtained K = F, (35) K = Φ (x)k Φ(x)d + Φ (x)αφ(x)ds + Φ (x)φ(x)ds, (36) S 1 S S 3 S 4 S 5 S 6 K IJ = ΦI (x)ρc Φ J (x)d + ΦI (x)αφ J (x)ds + ΦI (x)φ J (x)dω, (37) S 1 S S 3 S 4 S 5 S 6 F I = Φ I (x) Qd + αφ I (x) S1 ds + αφ I (x) S ds + Φ I (x) dγ. (38) S 1 S S 3 S 4 S 5 S
5 4. Numerical example o test te efficiency of te metod of te reproducing kernel particle approximation on a 3D steadystate eat conduction problem, a numerical analysis of te model sown in Fig. 1 as been carried out by using a linear basis function. e results obtained by te RKPM are also compared wit analytical metods. Cin. Pys. B ol. 19, No. 9 (010) is plotted in comparison wit numerical results wen x = 0.75 and y = 1.6. From Fig. 3, we can conclude tat numerical result obtained by RKPM metod is in good agreement wit te exact solution. Fig.. Nodes arrangement for example. Fig. 1. e 3D model (L = 1 m, W = m, H = 3 m, k = 1 W/mK, = 0 W/m K, = 0 C, S1 = 0 C, S = 0 C). Consider te following problem x + y + z = f(x, y, z), 0 x a, 0 y b, 0 z c, (39) wit boundary conditions (0, y, z) = (a, y, z) = 0, (x, 0, z) = (x, b, z) = 0, (x, y, 0) = (x, y, c) = 0, f(x, y, z) = π ( 1 a + 1 b + 1 c ) (40a) (40b) (40c) sin πx πy πz sin sin a b c. (41) e exact solution of Eq. (39) is (x, y, z) = sin πx πy πz sin sin a b c. (4) We use te RKPM to solve te above eat conduction problem wit penalty factor α = and scaling parameter d max = 1.13 and a = 1, b =, c = 3. e back-ground mes is adopted for integration in tis metod. Figure sows uniform nodes arrangement for example. e weigt function is cosen to be Gauss function and te bases are cosen to be linear. In Fig. 3, te exact solution (x, y, z) = sin πx a sin πy b sin πz c Fig. 3. Numerical solution and exact solution at location x = 0.75 and y = Conclusion e RKPM was used for solving te 3D steadystate eat conduction problem. e mesless property of RKPM is te most important advantage of tis sceme over te traditional mes dependent tecniques suc as te finite difference metod, finite element metod and boundary element metod. e mesless nature allows us to solve problems wit nonregular geometry. Our numerical results sow tat te tecnique is accurate and efficient. e RKPM can be extended for 3D unsteady-state eat conduction problems and nonlinear eat conduction problems
6 Cin. Pys. B ol. 19, No. 9 (010) References [1] Erwin K 1999 Advanced Engineering Matematics (8t ed.) (New York: Wiley and Sons) [] Selvadurai 1979 Partial Differential Equations in Mecanics (Berlin: Springe erlag) [3] Cristoper R J 00 Advanced Metods in Scientific (Ula: Universi ty of Uta, Spring Semester) [4] Belytscko, Krongauz Y and Organ D 1996 Comput. Met. Appl. Mec. Engng [5] Monagan J J 1988 Comput. Pys. Commun [6] Cen W 000 Mesfree Metod for Partial Differential Equations ol. 1 (Berlin: Springer erlag) p75 [7] Belytscko, Lu Y Y and Gu L 1994 Int. J. Numer. Met. Engng [8] Ceng R J and Ceng Y M 008 Acta Pys. Sin (in [9] Ceng R J and Ge H X 009 Cin. Pys. B [10] Atluri S N and Zu 1998 Comput. Mec. 117 [11] Liu W K, Jun S and Zang Y F 1995 Int. J. Numer. Met. Engng [1] Qin Y X and Ceng Y M 006 Acta Pys. Sin (in [13] Ceng R J and Ceng Y M 008 Appl. Numer. Mat [14] Ceng R J and Ceng Y M 007 Acta Pys. Sin (in [15] Dai B D and Ceng Y M 007 Acta Pys. Sin (in [16] Ceng Y M and Cen M J 003 Acta Mec. Sin (in [17] Ceng Y M and Peng M J 005 Science in Cina Ser. G Pysics, Mecanics & Astronomy [18] Ceng Y M and Li J H 005 Acta Pys. Sin (in [19] Cen L and Ceng Y M 008 Acta Pys. Sin (in [0] Cen L and Ceng Y M 008 Acta Pys. Sin (in [1] Peng M J and Ceng Y M 009 Engineering Analysis wit Boundary Elements 3 77 [] Cen C S 1998 Int. J. Numer. Metods Engng [3] Sing I 004 Numer. Heat ransfer Part A [4] Sing I 005 Numer. Heat ransfer Part A [5] Sing I and Prakas R 005 Heat and ecnology 1 73 [6] Sing I and Jain P K 005 Advanced in Engineering Software [7] Sing A and Sing I 007 Int. J. Heat Mass rans [8] Arefmanes A and Najafi M 005 J. Fluids Engng [9] Batra R C, Porfiri M and Spinello D 004 Int. J. Numer. Met. Engng [30] Liu Y, Zang X and Lu M W 005 singua Science and ecnology (in Englis) [31] Liu Y, Zang X and Lu M W 005 Numer. Heat ransfer [3] Sladek J, Sladek and Atluri S N 001 Comp. Modeling Engng. Sci. 43 [33] Sladek J, Sladek and Hellmic C 006 Comput. Mec [34] Ceng R J and Liew K M 009 Comput. Mec
Solving Poisson s equations by the Discrete Least Square meshless method
Boundary Elements and Oter Mes eduction Metods XXV 3 Solving Poisson s equations by te Discrete Least Square mesless metod H. Arzani & M. H. Afsar Said aaee University, Lavizan, eran, ran Department of
More informationAnalysis of the equal width wave equation with the mesh-free reproducing kernel particle Ritz method
Analysis of the equal width wave equation with the mesh-free reproducing kernel particle Ritz method Cheng Rong-Jun( 程荣军 ) a) and Ge Hong-Xia( 葛红霞 ) b) a) Ningbo Institute of Technology, Zhejiang University,
More informationMANY scientific and engineering problems can be
A Domain Decomposition Metod using Elliptical Arc Artificial Boundary for Exterior Problems Yajun Cen, and Qikui Du Abstract In tis paper, a Diriclet-Neumann alternating metod using elliptical arc artificial
More informationFourier Type Super Convergence Study on DDGIC and Symmetric DDG Methods
DOI 0.007/s095-07-048- Fourier Type Super Convergence Study on DDGIC and Symmetric DDG Metods Mengping Zang Jue Yan Received: 7 December 06 / Revised: 7 April 07 / Accepted: April 07 Springer Science+Business
More informationEffect of Numerical Integration on Meshless Methods
Effect of Numerical Integration on Mesless Metods Ivo Babuška Uday Banerjee Jon E. Osborn Qingui Zang Abstract In tis paper, we present te effect of numerical integration on mesless metods wit sape functions
More informationResearch Article Solution to Two-Dimensional Steady Inverse Heat Transfer Problems with Interior Heat Source Based on the Conjugate Gradient Method
Hindawi Matematical Problems in Engineering Volume 27, Article ID 286342, 9 pages ttps://doiorg/55/27/286342 Researc Article Solution to Two-Dimensional Steady Inverse Heat Transfer Problems wit Interior
More informationFinite Difference Methods Assignments
Finite Difference Metods Assignments Anders Söberg and Aay Saxena, Micael Tuné, and Maria Westermarck Revised: Jarmo Rantakokko June 6, 1999 Teknisk databeandling Assignment 1: A one-dimensional eat equation
More informationSuperconvergence of energy-conserving discontinuous Galerkin methods for. linear hyperbolic equations. Abstract
Superconvergence of energy-conserving discontinuous Galerkin metods for linear yperbolic equations Yong Liu, Ci-Wang Su and Mengping Zang 3 Abstract In tis paper, we study superconvergence properties of
More informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More informationThe Sinc-Collocation Method for Solving the Telegraph Equation
Te Sinc-Collocation Metod for Solving te Telegrap Equation E. Hesameddini *1, E. Asadolaifard Department of Matematics, Faculty of Basic Sciences, Siraz University of Tecnology, Siraz, Iran *1 esameddini@sutec.ac.ir;
More informationThe entransy dissipation minimization principle under given heat duty and heat transfer area conditions
Article Engineering Termopysics July 2011 Vol.56 No.19: 2071 2076 doi: 10.1007/s11434-010-4189-x SPECIAL TOPICS: Te entransy dissipation minimization principle under given eat duty and eat transfer area
More informationMathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative
Matematics 5 Workseet 11 Geometry, Tangency, and te Derivative Problem 1. Find te equation of a line wit slope m tat intersects te point (3, 9). Solution. Te equation for a line passing troug a point (x
More information1. Introduction. We consider the model problem: seeking an unknown function u satisfying
A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced
More informationInf sup testing of upwind methods
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Met. Engng 000; 48:745 760 Inf sup testing of upwind metods Klaus-Jurgen Bate 1; ;, Dena Hendriana 1, Franco Brezzi and Giancarlo
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
More informationAPPROXIMATION OF CRYSTALLINE DENDRITE GROWTH IN TWO SPACE DIMENSIONS. Introduction
Acta Mat. Univ. Comenianae Vol. LXVII, 1(1998), pp. 57 68 57 APPROXIMATION OF CRYSTALLINE DENDRITE GROWTH IN TWO SPACE DIMENSIONS A. SCHMIDT Abstract. Te pase transition between solid and liquid in an
More informationThe Laplace equation, cylindrically or spherically symmetric case
Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,
More informationNumerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems
Applied Matematics, 06, 7, 74-8 ttp://wwwscirporg/journal/am ISSN Online: 5-7393 ISSN Print: 5-7385 Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for
More informationBending analysis of a functionally graded piezoelectric cantilever beam
Science in Cina Series G: Pysics Mecanics & Astronomy 7 Science in Cina Press Springer-Verlag Bending analysis of a functionally graded pieoelectric cantilever beam YU Tao & ZHONG Zeng Scool of Aerospace
More informationParameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
More informationKhyber Pakhtunkhwa Agricultural University, Peshawar, Pakistan
ttp://www.lifesciencesite.com Septic B-Spline Collocation metod for numerical solution of te Equal Widt Wave (EW) equation 1 Fazal-i-Haq, 2 Inayat Ali sa and 3 Sakeel Amad 1 Department of Matematics, Statistics
More informationA method of Lagrange Galerkin of second order in time. Une méthode de Lagrange Galerkin d ordre deux en temps
A metod of Lagrange Galerkin of second order in time Une métode de Lagrange Galerkin d ordre deux en temps Jocelyn Étienne a a DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, Great-Britain.
More informationA New Fifth Order Finite Difference WENO Scheme for Hamilton-Jacobi Equations
A New Fift Order Finite Difference WENO Sceme for Hamilton-Jacobi Equations Jun Zu, 1 Jianxian Qiu 1 College of Science, Nanjing University of Aeronautics and Astronautics, Nanjing, Jiangsu 10016, People
More informationFEM solution of the ψ-ω equations with explicit viscous diffusion 1
FEM solution of te ψ-ω equations wit explicit viscous diffusion J.-L. Guermond and L. Quartapelle 3 Abstract. Tis paper describes a variational formulation for solving te D time-dependent incompressible
More informationResearch Article Cubic Spline Iterative Method for Poisson s Equation in Cylindrical Polar Coordinates
International Scolarly Researc Network ISRN Matematical Pysics Volume 202, Article ID 2456, pages doi:0.5402/202/2456 Researc Article Cubic Spline Iterative Metod for Poisson s Equation in Cylindrical
More informationDecay of solutions of wave equations with memory
Proceedings of te 14t International Conference on Computational and Matematical Metods in Science and Engineering, CMMSE 14 3 7July, 14. Decay of solutions of wave equations wit memory J. A. Ferreira 1,
More informationPreconditioning in H(div) and Applications
1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition
More informationFinding and Using Derivative The shortcuts
Calculus 1 Lia Vas Finding and Using Derivative Te sortcuts We ave seen tat te formula f f(x+) f(x) (x) = lim 0 is manageable for relatively simple functions like a linear or quadratic. For more complex
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationCONVERGENCE ANALYSIS OF FINITE ELEMENT SOLUTION OF ONE-DIMENSIONAL SINGULARLY PERTURBED DIFFERENTIAL EQUATIONS ON EQUIDISTRIBUTING MESHES
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume, Number, Pages 57 74 c 5 Institute for Scientific Computing and Information CONVERGENCE ANALYSIS OF FINITE ELEMENT SOLUTION OF ONE-DIMENSIONAL
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationN igerian Journal of M athematics and Applications V olume 23, (2014), 1 13
N igerian Journal of M atematics and Applications V olume 23, (24), 3 c N ig. J. M at. Appl. ttp : //www.kwsman.com CONSTRUCTION OF POLYNOMIAL BASIS AND ITS APPLICATION TO ORDINARY DIFFERENTIAL EQUATIONS
More informationResearch Article Smoothing Analysis of Distributive Red-Black Jacobi Relaxation for Solving 2D Stokes Flow by Multigrid Method
Matematical Problems in Engineering Volume 205, Article ID 57298, 7 pages ttp://dx.doi.org/0.55/205/57298 Researc Article Smooting Analysis of Distributive Red-Black Jacobi Relaxation for Solving 2D Stokes
More informationGENERALIZED DIFFERENTIAL TRANSFORM METHOD FOR SOLUTIONS OF NON-LINEAR PARTIAL DIFFERENTIAL EQUATIONS OF FRACTIONAL ORDER
December 7 Volume Issue JETIR (ISSN-39-56) GENERALIZED DIFFERENTIAL TRANSFORM METHOD FOR SOLTIONS OF NON-LINEAR PARTIAL DIFFERENTIAL EQATIONS OF FRACTIONAL ORDER Deepanjan Das Department of Matematics
More informationOn convergence of the immersed boundary method for elliptic interface problems
On convergence of te immersed boundary metod for elliptic interface problems Zilin Li January 26, 2012 Abstract Peskin s Immersed Boundary (IB) metod is one of te most popular numerical metods for many
More informationA C 2 -continuous control-volume technique based on Cartesian grids and two-node integrated-rbf elements for second-order elliptic problems
A C 2 -continuous control-volume tecnique based on Cartesian grids and two-node integrated-rbf elements for second-order elliptic problems D.-A. An-Vo 1, N. Mai-Duy 1 and T. Tran-Cong 1 Abstract: Tis paper
More informationFINITE ELEMENT DUAL SINGULAR FUNCTION METHODS FOR HELMHOLTZ AND HEAT EQUATIONS
J. KSIAM Vol.22, No.2, 101 113, 2018 ttp://dx.doi.org/10.12941/jksiam.2018.22.101 FINITE ELEMENT DUAL SINGULAR FUNCTION METHODS FOR HELMHOLTZ AND HEAT EQUATIONS DEOK-KYU JANG AND JAE-HONG PYO DEPARTMENT
More informationHigh-Order Extended Finite Element Methods for Solving Interface Problems
Hig-Order Extended Finite Element Metods for Solving Interface Problems arxiv:1604.06171v1 [mat.na] 21 Apr 2016 Fei Wang Yuanming Xiao Jincao Xu ey words. Elliptic interface problems, unfitted mes, extended
More informationch (for some fixed positive number c) reaching c
GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan
More informationSmoothness of solutions with respect to multi-strip integral boundary conditions for nth order ordinary differential equations
396 Nonlinear Analysis: Modelling and Control, 2014, Vol. 19, No. 3, 396 412 ttp://dx.doi.org/10.15388/na.2014.3.6 Smootness of solutions wit respect to multi-strip integral boundary conditions for nt
More informationExercise 19 - OLD EXAM, FDTD
Exercise 19 - OLD EXAM, FDTD A 1D wave propagation may be considered by te coupled differential equations u x + a v t v x + b u t a) 2 points: Derive te decoupled differential equation and give c in terms
More informationHeat Transfer and MHD Boundary Layer Flow over a Rotating Disk
J. Basic. Appl. Sci. Res., 4()37-33, 4 4, TextRoad Publication ISSN 9-434 Journal of Basic and Applied Scientific Researc www.textroad.com Heat Transfer and MHD Boundary Layer Flow over a Rotating Disk
More informationNew Fourth Order Quartic Spline Method for Solving Second Order Boundary Value Problems
MATEMATIKA, 2015, Volume 31, Number 2, 149 157 c UTM Centre for Industrial Applied Matematics New Fourt Order Quartic Spline Metod for Solving Second Order Boundary Value Problems 1 Osama Ala yed, 2 Te
More informationExplicit Hyperbolic Reconstructed Discontinuous Galerkin Methods for Time-Dependent Problems
AIAA AVIATION Forum June 25-29 218 Atlanta Georgia 218 Fluid Dynamics Conference 1.2514/6.218-427 Explicit Hyperbolic Reconstructed Discontinuous Galerkin Metods for Time-Dependent Problems Jialin Lou
More informationJian-Guo Liu 1 and Chi-Wang Shu 2
Journal of Computational Pysics 60, 577 596 (000) doi:0.006/jcp.000.6475, available online at ttp://www.idealibrary.com on Jian-Guo Liu and Ci-Wang Su Institute for Pysical Science and Tecnology and Department
More informationApplication of Quintic B-splines Collocation Method on Some Rosenau Type Nonlinear Higher Order Evolution Equations.
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.13(2012) No.2,pp.142-152 Application of Quintic B-splines Collocation Metod on Some Rosenau Type Nonlinear Higer
More informationReflection Symmetries of q-bernoulli Polynomials
Journal of Nonlinear Matematical Pysics Volume 1, Supplement 1 005, 41 4 Birtday Issue Reflection Symmetries of q-bernoulli Polynomials Boris A KUPERSHMIDT Te University of Tennessee Space Institute Tullaoma,
More informationWYSE Academic Challenge 2004 Sectional Mathematics Solution Set
WYSE Academic Callenge 00 Sectional Matematics Solution Set. Answer: B. Since te equation can be written in te form x + y, we ave a major 5 semi-axis of lengt 5 and minor semi-axis of lengt. Tis means
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
More informationImproved Rotated Finite Difference Method for Solving Fractional Elliptic Partial Differential Equations
American Scientific Researc Journal for Engineering, Tecnolog, and Sciences (ASRJETS) ISSN (Print) 33-44, ISSN (Online) 33-44 Global Societ of Scientific Researc and Researcers ttp://asrjetsjournal.org/
More informationArbitrary order exactly divergence-free central discontinuous Galerkin methods for ideal MHD equations
Arbitrary order exactly divergence-free central discontinuous Galerkin metods for ideal MHD equations Fengyan Li, Liwei Xu Department of Matematical Sciences, Rensselaer Polytecnic Institute, Troy, NY
More informationThermal Bending of Circular Plates for Non-axisymmetrical Problems
Copyrigt 2 Tec Science Press SL, vol.4, no.2, pp.5-2, 2 Termal Bending of Circular Plates for Non-axisymmetrical Problems Dong Zengzu Peng Weiong and Li Suncai Abstract: Due to te complexity of termal
More informationA First-Order System Approach for Diffusion Equation. I. Second-Order Residual-Distribution Schemes
A First-Order System Approac for Diffusion Equation. I. Second-Order Residual-Distribution Scemes Hiroaki Nisikawa W. M. Keck Foundation Laboratory for Computational Fluid Dynamics, Department of Aerospace
More informationAN ANALYSIS OF AMPLITUDE AND PERIOD OF ALTERNATING ICE LOADS ON CONICAL STRUCTURES
Ice in te Environment: Proceedings of te 1t IAHR International Symposium on Ice Dunedin, New Zealand, nd t December International Association of Hydraulic Engineering and Researc AN ANALYSIS OF AMPLITUDE
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationDIFFERENTIAL POLYNOMIALS GENERATED BY SECOND ORDER LINEAR DIFFERENTIAL EQUATIONS
Journal o Applied Analysis Vol. 14, No. 2 2008, pp. 259 271 DIFFERENTIAL POLYNOMIALS GENERATED BY SECOND ORDER LINEAR DIFFERENTIAL EQUATIONS B. BELAÏDI and A. EL FARISSI Received December 5, 2007 and,
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More informationCALCULATION OF COLLAPSE PRESSURE IN SHALE GAS FORMATION AND THE INFLUENCE OF FORMATION ANISOTROPY
CALCULATION OF COLLAPSE PRESSURE IN SHALE GAS FORMATION AND THE INFLUENCE OF FORMATION ANISOTROPY L.Hu, J.Deng, F.Deng, H.Lin, C.Yan, Y.Li, H.Liu, W.Cao (Cina University of Petroleum) Sale gas formations
More informationarxiv: v1 [math.na] 3 Nov 2011
arxiv:.983v [mat.na] 3 Nov 2 A Finite Difference Gost-cell Multigrid approac for Poisson Equation wit mixed Boundary Conditions in Arbitrary Domain Armando Coco, Giovanni Russo November 7, 2 Abstract In
More informationThe development of contact and noncontact technique to study the heat dissipation in metals under loading
Te development of contact and noncontact tecnique to study te eat dissipation in metals under loading More info about tis article: ttp://www.ndt.net/?id=73 Abstract * ICMM UB RAS, Ac. Koroleva Str., 63
More informationDIRECTLY SOLVING SECOND ORDER LINEAR BOUNDARY VALUE PROBLEMS OF ORDINARY DIFFERENTIAL EQUATIONS. Ra ft Abdelrahim 1, Z. Omar 2
International Journal of Pure and Applied Matematics Volume 6 No. 6, -9 ISSN: - (printed version); ISSN: -95 (on-line version) url: ttp://www.ijpam.eu doi:.7/ijpam.v6i. PAijpam.eu DIRECTLY SOLVING SECOND
More informationSeepage Analysis through Earth Dam Based on Finite Difference Method
J. Basic. Appl. Sci. Res., (11)111-1, 1 1, TetRoad Publication ISSN -44 Journal of Basic and Applied Scientific Researc www.tetroad.com Seepage Analysis troug Eart Dam Based on Finite Difference Metod
More informationarxiv: v1 [math.ap] 4 Aug 2017
New Two Step Laplace Adam-Basfort Metod for Integer an Non integer Order Partial Differential Equations arxiv:178.1417v1 [mat.ap] 4 Aug 17 Abstract Rodrigue Gnitcogna*, Abdon Atangana** *Department of
More informationEXTENSION OF A POSTPROCESSING TECHNIQUE FOR THE DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC EQUATIONS WITH APPLICATION TO AN AEROACOUSTIC PROBLEM
SIAM J. SCI. COMPUT. Vol. 26, No. 3, pp. 821 843 c 2005 Society for Industrial and Applied Matematics ETENSION OF A POSTPROCESSING TECHNIQUE FOR THE DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC EQUATIONS
More information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
More informationParametric Spline Method for Solving Bratu s Problem
ISSN 749-3889 print, 749-3897 online International Journal of Nonlinear Science Vol4202 No,pp3-0 Parametric Spline Metod for Solving Bratu s Problem M Zarebnia, Z Sarvari 2,2 Department of Matematics,
More informationLEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS
SIAM J. NUMER. ANAL. c 998 Society for Industrial Applied Matematics Vol. 35, No., pp. 393 405, February 998 00 LEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS YANZHAO CAO
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More informationComputers and Mathematics with Applications. A nonlinear weighted least-squares finite element method for Stokes equations
Computers Matematics wit Applications 59 () 5 4 Contents lists available at ScienceDirect Computers Matematics wit Applications journal omepage: www.elsevier.com/locate/camwa A nonlinear weigted least-squares
More informationarxiv: v1 [math.na] 9 Sep 2015
arxiv:509.02595v [mat.na] 9 Sep 205 An Expandable Local and Parallel Two-Grid Finite Element Sceme Yanren ou, GuangZi Du Abstract An expandable local and parallel two-grid finite element sceme based on
More informationJournal of Computational and Applied Mathematics
Journal of Computational and Applied Matematics 94 (6) 75 96 Contents lists available at ScienceDirect Journal of Computational and Applied Matematics journal omepage: www.elsevier.com/locate/cam Smootness-Increasing
More informationarxiv: v1 [math.na] 27 Jan 2014
L 2 -ERROR ESTIMATES FOR FINITE ELEMENT APPROXIMATIONS OF BOUNDARY FLUXES MATS G. LARSON AND ANDRÉ MASSING arxiv:1401.6994v1 [mat.na] 27 Jan 2014 Abstract. We prove quasi-optimal a priori error estimates
More informationNumerical study of the Benjamin-Bona-Mahony-Burgers equation
Bol. Soc. Paran. Mat. (3s.) v. 35 1 (017): 17 138. c SPM ISSN-175-1188 on line ISSN-0037871 in press SPM: www.spm.uem.br/bspm doi:10.569/bspm.v35i1.8804 Numerical study of te Benjamin-Bona-Maony-Burgers
More informationMultigrid Methods for Obstacle Problems
Multigrid Metods for Obstacle Problems by Cunxiao Wu A Researc Paper presented to te University of Waterloo in partial fulfillment of te requirement for te degree of Master of Matematics in Computational
More informationParallel algorithm for convection diffusion system based on least squares procedure
Zang et al. SpringerPlus (26 5:69 DOI.86/s464-6-3333-8 RESEARCH Parallel algoritm for convection diffusion system based on least squares procedure Open Access Jiansong Zang *, Hui Guo 2, Hongfei Fu 3 and
More informationA Hybrid Mixed Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Hybrid Mixed Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationCritical control in transcritical shallow-water flow over two obstacles
Lougboroug University Institutional Repository Critical control in transcritical sallow-water flow over two obstacles Tis item was submitted to Lougboroug University's Institutional Repository by te/an
More informationMIMO decorrelation for visible light communication based on angle optimization
MIMO decorrelation for visible ligt communication based on angle optimization Haiyong Zang, and Yijun Zu Citation: AIP Conference Proceedings 80, 09005 (07); View online: ttps://doi.org/0.03/.4977399 View
More informationLeast-squares collocation meshless method
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2001; 51:1089 1100 Least-squares collocation meshless method Xiong Zhang ;, Xiao-HuLiu, Kang-ZuSong and Ming-Wan Lu
More informationAnalytic study on cold field electron emission from metallic nanowall array *
Analytic study on cold field electron emission from metallic nanowall array * Xizou Qin, Weiliang Wang and Zibing Li a) State Key Laboratory of Optoelectronic aterials and Tecnologies, Scool of Pysics
More informationJournal of Applied Science and Agriculture. The Effects Of Corrugated Geometry On Flow And Heat Transfer In Corrugated Channel Using Nanofluid
Journal o Applied Science and Agriculture, 9() February 04, Pages: 408-47 AENSI Journals Journal o Applied Science and Agriculture ISSN 86-9 Journal ome page: www.aensiweb.com/jasa/index.tml Te Eects O
More informationDepartment of Mathematical Sciences University of South Carolina Aiken Aiken, SC 29801
RESEARCH SUMMARY AND PERSPECTIVES KOFFI B. FADIMBA Department of Matematical Sciences University of Sout Carolina Aiken Aiken, SC 29801 Email: KoffiF@usca.edu 1. Introduction My researc program as focused
More information4.2 - Richardson Extrapolation
. - Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Definition Let x n n converge to a number x. Suppose tat n n is a sequence
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More informationRecent Advances in Time-Domain Maxwell s Equations in Metamaterials
Recent Advances in Time-Domain Maxwell s Equations in Metamaterials Yunqing Huang 1, and Jicun Li, 1 Hunan Key Laboratory for Computation and Simulation in Science and Engineering, Xiangtan University,
More informationA Hybridizable Discontinuous Galerkin Method for the Compressible Euler and Navier-Stokes Equations
A Hybridizable Discontinuous Galerkin Metod for te Compressible Euler and Navier-Stokes Equations J. Peraire and N. C. Nguyen Massacusetts Institute of Tecnology, Cambridge, MA 02139, USA B. Cockburn University
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationA posteriori error estimation for unilateral contact with matching and non-matching meshes
Comput. Metods Appl. Mec. Engrg. 186 (2000) 65±83 www.elsevier.com/locate/cma A posteriori error estimation for unilateral contact wit matcing and non-matcing meses Patrice Coorevits a, *, Patrick Hild
More informationFlapwise bending vibration analysis of double tapered rotating Euler Bernoulli beam by using the differential transform method
Meccanica 2006) 41:661 670 DOI 10.1007/s11012-006-9012-z Flapwise bending vibration analysis of double tapered rotating Euler Bernoulli beam by using te differential transform metod Ozge Ozdemir Ozgumus
More informationWeak Solutions for Nonlinear Parabolic Equations with Variable Exponents
Communications in Matematics 25 217) 55 7 Copyrigt c 217 Te University of Ostrava 55 Weak Solutions for Nonlinear Parabolic Equations wit Variable Exponents Lingeswaran Sangerganes, Arumugam Gurusamy,
More informationNumerical analysis of a free piston problem
MATHEMATICAL COMMUNICATIONS 573 Mat. Commun., Vol. 15, No. 2, pp. 573-585 (2010) Numerical analysis of a free piston problem Boris Mua 1 and Zvonimir Tutek 1, 1 Department of Matematics, University of
More informationImplicit-explicit variational integration of highly oscillatory problems
Implicit-explicit variational integration of igly oscillatory problems Ari Stern Structured Integrators Worksop April 9, 9 Stern, A., and E. Grinspun. Multiscale Model. Simul., to appear. arxiv:88.39 [mat.na].
More informationDELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science
DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Matematics and Computer Science. ANSWERS OF THE TEST NUMERICAL METHODS FOR DIFFERENTIAL EQUATIONS (WI3097 TU) Tuesday January 9 008, 9:00-:00
More informationMoving Mesh Methods for Singular Problems on a Sphere Using Perturbed Harmonic Mappings
Moving Mes Metods for Singular Problems on a Spere Using Perturbed Harmonic Mappings Yana Di, Ruo Li, Tao Tang and Pingwen Zang Marc 4, 6 Abstract Tis work is concerned wit developing moving mes strategies
More informationA Polynomial Adaptive LCP Scheme for Viscous Compressible Flows
A Polynomial Adaptive LCP Sceme for Viscous Compressible Flows J.S. Cagnone, B.C. Vermeire, and S.. Nadaraja Department of Mecanical Engineering, McGill University, Montreal, Canada, H3A 2S6 Email: jean-sebastien.cagnone@mail.mcgill.ca
More informationThe Derivative as a Function
Section 2.2 Te Derivative as a Function 200 Kiryl Tsiscanka Te Derivative as a Function DEFINITION: Te derivative of a function f at a number a, denoted by f (a), is if tis limit exists. f (a) f(a + )
More informationA proof in the finite-difference spirit of the superconvergence of the gradient for the Shortley-Weller method.
A proof in te finite-difference spirit of te superconvergence of te gradient for te Sortley-Weller metod. L. Weynans 1 1 Team Mempis, INRIA Bordeaux-Sud-Ouest & CNRS UMR 551, Université de Bordeaux, France
More informationLarge eddy simulation of turbulent flow downstream of a backward-facing step
Available online at www.sciencedirect.com Procedia Engineering 31 (01) 16 International Conference on Advances in Computational Modeling and Simulation Large eddy simulation of turbulent flow downstream
More information