Two numerical methods for solving a backward heat conduction problem
|
|
- Simon Norton
- 5 years ago
- Views:
Transcription
1 pplied Mathematics and Computation 79 (6) wwwelseviercom/locate/amc Two numerical methods for solving a backward heat conduction problem Xiang-Tuan Xiong, Chu-Li Fu *, hi Qian epartment of Mathematics, Lanzhou University, Lanzhou 73, People s Republic of China bstract We introduce a central difference method and a quasi-reversibility method for solving a backward heat conduction problem (BHCP) numerically For these two numerical methods, we give the stability analysis Meanwhile, we investigate the roles of regularization parameters in these two methods Numerical results show that our algorithm is effective Ó 5 Elsevier Inc ll rights reserved Keywords: Inverse problem; Backward heat conduction; Central difference; Quasi-reversibility; Regularization Introduction The backward heat conduction problem (BHCP) is also referred to as final boundary value problem In general, no solution which satisfies the heat conduction equation with final data and the boundary conditions exists Even if a solution exists, it will not be continuously dependent on the final data The BHCP is a typical example of an ill-posed problem which is unstable by numerical methods and requires special regularization methods In the context of approximation method for this problem, many approaches have been investigated Such authors as Lattes and Lions [], Showalter [], mes et al [3], Miller [4] have approximated the BHCP by quasi-reversibility methods Schröter and Tautenhahn [5] established an optimal error estimate for a special BHCP Mera and Jourhmane used many numerical methods with regularization techniques to approximate the problem in [6 8], etc mollification method has been studied by Hào in[9] Kirkup and Wadsworth used an operator-splitting method in [] difference approximation method for solving sideways heat equation was provided by Eldénin[] So far in the literature (cf [3,] and the references therein), most of the authors used the eigenfunctions and eigenvalues to reconstruct the solution of the BHCP by many quasi-reversibility methods numerically, However, the eigenfunctions and eigenvalues are in general not available and the labor needed to compute these and the corresponding Fourier coefficients is very onerous In this paper, we use two regularization methods to solve the BHCP in one-dimensional setting numerically, but these methods can be generalized to two-dimensional case * Corresponding author addresses: xiongxt4@stlzueducn (X-T Xiong), fuchuli@lzueducn (C-L Fu) 96-33/$ - see front matter Ó 5 Elsevier Inc ll rights reserved doi:6/jamc54
2 The paper is organized as follows In the forthcoming section, we will present the mathematical problem on a BHCP; in Section 3, we review two regularization methods, one is central difference regularization method, the other is a special quasi-reversibility regularization method; in Section 4, some finite difference schemes are constructed for the inverse problem and the numerical stability analysis is provided; in Section 5, numerical examples are tested to verify the effect of the numerical schemes Mathematical problem The direct problem X-T Xiong et al / pplied Mathematics and Computation 79 (6) We consider the following heat equation: u t ðx; tþ u xx ðx; tþ; p < x < p; < t < T ; uðp; tþ sðtþ; < t < T ; uð p; tþ lðtþ; < t < T ; uðx; Þ f ðxþ; p < x < p ð:þ Solving the equation with given s(t); l(t) and f(x) is called a direct problem From the theory of heat equation, we can see that for s(t); l(t) and f(x) in some function space there exists a unique solution [3] The inverse problem Consider the following problem: u t ðx; tþ u xx ðx; tþ; p < x < p; < t < T ; uðp; tþ sðtþ; < t < T ; uð p; tþ lðtþ; < t < T ; uðx; T ÞgðxÞ; p < x < p The inverse problem is to determine the value of u(x, t) for 6 t < T from the data s(t); l(t) and g(x) If the solution exists, then the problem has a unique solution (cf [4, p 64]) The data g(x) are based on (physical) observations and are not known with complete accuracy, due to the ill-posedness of the BHCP, a small error in the data g(x) can cause an arbitrarily large error in the solution u(x, t) Now we want to reconstruct the temperature distribution u(x, t) for 6 t < T by two different regularization methods 3 Central difference regularization By using the central difference with step length h to approximate the second derivative u xx, we can get the following problem: u t ðx; tþ uðxþh;tþ uðx;tþþuðx h;tþ ; p < x < p; < t < T ; h uðp; tþ sðtþ; < t < T ; ð3:þ uð p; tþ lðtþ; < t < T ; uðx; T ÞgðxÞ; p < x < p Let ~t T t; wðx;~tþ uðx; tþ, we have w ~t w xx ; p < x < p; < ~t < T ð3:þ If w xx at x i is replaced by a second-order central difference, then (3) becomes ð:þ w ~t ðx i ;~tþ h ½wðx i þ h;~tþ wðx i ;~tþþwðx i h;~tþš ð3:3þ
3 37 X-T Xiong et al / pplied Mathematics and Computation 79 (6) Let x i = p +(i )h, for i =,,n + h p ; w n i w i ð~tþ wð p þði Þh;~tÞ Then by the boundary condition of (3) there hold w ð~tþ lð~tþ and w nþ ð~tþ sð~tþ Meanwhile, Eq (33) with initial boundary conditions can be discretized as w ð~tþ w ð~tþ w ð~tþ h h h h h w i ð~tþ B B w i ð~tþ B C B h w n ð~tþ h fflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl} w ~t n ð~tþ w nþ ð~tþ ðn Þðn Þ fflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflffl} n w ðþ gðx Þ w i ðþ gðx i Þ ð3:5þ B C w n ðþ gðx n Þ This is an ordinary differential system and there are many numerical methods such as Euler method, Runge Kutta method for the system But we find that the eigenvalues of are k k 4 sin kp P, k (the positive integer set) Therefore, the numerical method for (34) and (35) is unstable However, via the variable h ðn Þ transformation [5] v i ð~tþ e a~t w i ð~tþ; i ; ; n þ ; ð3:6þ where a > is a contraction factor to be determined such that (33) in terms of v i ð~tþ is changed to ov i ð~tþ o~t h a v i ð~tþ h ½v iþð~tþþv i ð~tþš; ð3:7þ similarly, we have v ð~tþ v ð~tþ e a~t lð~tþ a h h a h v i ð~tþ h h B B C C i ð~tþ B h B a h v n ð~tþ h fflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl} v n ð~tþ e a~t sð~tþ ~t B ðn Þðn Þ v ðþ gðx Þ v i ðþ gðx i Þ ð3:9þ B C v n ðþ gðx n Þ We choose a b=h P h such that the eigenvalues of matrix B are negative, thus the numerical methods such as Runge Kutta method for (38) and (39) are stable fter v i ð~tþ is obtained numerically, we can obtain w i ð~tþ e a~t v i ð~tþ, furthermore wðih;~tþ uðih; T ~tþ u i ðtþ In the central difference regularization method, the space step length h plays the role of regularization parameter ccording to the general regularization theory, the h should be neither too small nor too large ð3:4þ
4 X-T Xiong et al / pplied Mathematics and Computation 79 (6) By our experience, the parameter h is selected as h ln =, where d is the noise level of the data g(x) d We will see this choice rule is very effective in the subsequent numerical tests 4 Quasi-reversibility regularization method The initial boundary value problem () is replaced by the following problem: u t ðx; tþ u xx ðx; tþþu xxt ðx; tþ; p < x < p; < t < T ; uðp; tþ sðtþ; < t < T ; uð p; tþ lðtþ; < t < T ; uðx; T ÞgðxÞ; p < x < p; ð4:þ where is a small positive parameter For sufficiently small the solution of (4) approximates the solution (if it exists) of () in some sense This is one of well-known quasi-reversibility methods [] For the above mentioned problem, Ewing [] has presented a choice rule of the regularization parameter, ie, ln, d where d denotes the noise level of the data g(x), and the error estimate between the approximate solution and the exact solution is given in L (R)-norm Similarly, we take ~t T t, then problem (4) becomes u ~t ðx;~tþþu xx ðx;~tþ u xx~t ðx;~tþ ; p < x < p; < ~t < T ; uðp;~tþ sð~tþ; < ~t < T ; ð4:þ uð p;~tþ lð~tþ; < ~t < T ; uðx; Þ gðxþ; p < x < p The problem has a unique solution if a solution exists Now we prove it for two-dimensional case Theorem There exists a unique solution (if it exists) for the problem: u ~t ðx;~tþþuðx;~tþ u ~t ðx;~tþ ; ð; T Þ; uðx;~tþ hðx;~tþ; on o ð; T Þ; uðx; Þ gðxþ; in ; where is a bounded subset in R, is the Laplace operator, > ð4:3þ Proof We only need to prove the following problem has the zero solution: w ~t ðx;~tþþwðx;~tþ w ~t ðx;~tþ ; ð; T Þ; wðx;~tþ ; on o ð; T Þ; wðx; Þ ; in Set uð~tþ w ð~tþþjrwð~tþj dx; then duð~tþ ww ~t þ rw rw ~t dx d~t ue to the Green s second formula, we have duð~tþ ww ~t ww ~t dx wðw ~t w ~t Þdx d~t 6 w ð~tþþjrwð~tþj dx uð ~tþ wð wþdx jrwj dx ð4:4þ
5 374 X-T Xiong et al / pplied Mathematics and Computation 79 (6) Therefore, we have uð~tþ 6 uðþe ~t ð4:5þ since u() =, there holds uð~tþ 6 Hence w = h Now we construct the finite difference schemes for solving problem (4), let x i = p +(i )h, i =,,n +,~t j ðj Þs, j =,,m +, where h p ; s T Let n m uj i uðx i ;~t j Þ represent the value of the numerical solution of (4) at the mesh point ðx i ; ~t j Þ, then Eq (4) is discretized as ru jþ iþ þ r þ u jþ i ru jþ i s h þ r u j iþ þ h þ r þ u j i s h þ r u j i ; ð4:6þ where r, i =,,n; j =,,m h s Now we discuss the stability of difference schemes (46) by verifying the Von Neumann condition The propagation factor can be found þ 4 þ r sin rh s h Gðr; sþ þ 4r s sin rh It is easy to verify the fact that the Von Neumann condition jgðr; sþj 6 þ cs holds with c Hence, the numerical algorithm (46) is stable 5 Numerical examples For convenience, we take s(t) = l(t) = in() Example We consider the following direct problem: u t ðx; tþ u xx ðx; tþ; p < x < p; < t < ; uðp; tþ ; < t < ; uð p; tþ ; < t < with the initial condition: p þ x; p 6 x 6 ; uðx; Þ p x; 6 x 6 p Let x = p, x i = p +(i )h, i =,,n; x n + =, t j =(j )s, j =,,m +, the space step length h p and the time step length s, then we solve this problem by an explicit difference scheme in n m the following form: with u jþ i ru j iþ þð rþuj i þ ru j i ; j ; ; m; i ; ; n uðx nþ ; t j Þu j nþ ; j ; ; m þ ; ð5:þ uðx ; t j Þu j ; j ; ; m þ uðx i ; t Þ p þ x i; p 6 x i 6 ; p x i ; 6 x i 6 p; where r s, and it requires r < for numerical stability reasons h The numerical result for g(x) = u(x, T = ) is shown in Fig, where n =, m = 5 Now we solve the inverse problem by the g(x) generated numerically by the direct problem via the two different regularization methods We choose to restore the solution f(x) att = We denote the numerical result ð5:þ
6 X-T Xiong et al / pplied Mathematics and Computation 79 (6) g(x) Fig g(x) computed by (5) of the inverse problem as f % (x) If we introduce random noises e to data g(x), ie, g e (x i )=g(x i )+erand(i), where rand(i) is a random number between [, ], then the total noise d can be measured in the sense of root mean square (RMS) error according to vffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi u Xn d : t ðg n e ðx i Þ gðx i ÞÞ ð5:3þ i In the computation, we choose the regularization parameters h ln for the central difference regularization method and ln d for the quasi-reversibility method, respectively The numerical results are d shown in Fig (quasi-reversibility), Fig 3(central difference) with d = Example Consider the problem u t u xx ; p < x < p; < t < ; uð p; tþ uðp; tþ ; < t < ; ð5:4þ uðx; Þ e sin x; p < x < p 35 3 numerical solution exact solution Fig f % and f, m = 5, n =, = 45
7 376 X-T Xiong et al / pplied Mathematics and Computation 79 (6) numerical solution exact solution Fig 3 f % and f, b =,m = 5, n =4,h = p/4 8 numerical solution exact solution Fig 4 f % and f, m =, n = 3, = 9 8 numerical solution exact solution Fig 5 f % and f, b =3,m = 5, n =5,h = p/5
8 The exact solution is uðx; tþ e t sin x ð5:5þ The numerical result for the inverse problem is also denoted as f % (x) See Fig 4 (quasi-reversibility) and Fig 5 (central difference), where d = From Examples and, we conclude that the choice rules of the regularization parameters h and are very effective By our numerical experiments, we can see that the accuracy of the numerical results increases with the decreasing T; at the same time, the numerical solutions of the quasi-reversibility method depend on the parameter continuously This accords with the theoretical result (cf [6]) Here, we will not give the numerical results 6 Conclusions In this paper, we discussed two regularization methods for the one-dimensional backward heat conduction problem We presented two algorithms for the inverse problem Numerical results show that these method are effective with two appropriately chosen regularization parameters The algorithm for two-dimensional case is to be considered cknowledgements The work is supported by the NNSF of China (Nos 75 and 5779) and the NSF of Gansu Province of China (No 3S5-5-5) References X-T Xiong et al / pplied Mathematics and Computation 79 (6) [] R Lattes, JL Lions, Methode de Quasi-Reversibility et pplications, unod, Paris, 967 (English Transl: R Bellman, Elsevier, New York, 969) [] RE Showalter, The final value problem for evolution equations, J Math nal ppl 47 (974) [3] K mes, WC Gordon, JF Epperson, SF Oppenhermer, comparison of regularizations for an ill-posed problem, Math Comput 67 (998) [4] K Miller, Stabilized quasireversibility and other nearly best possible methods for non-well-posed problems, in: Symposium on Non- Well-Posed Problems and Logarithmic Convexity, Lecture Notes in Mathematics, vol 36, Springer-Verlag, Berlin, 973, pp 6 76 [5] T Schröter, U Tautenhahn, On optimal regularization methods for the backward heat equation, nal nw 5 (996) [6] NS Mera, L Elliott, B Ingham, Lesnic, n iterative boundary element method for solving the one dimensional backward heat conduction problem, Int J Heat Mass Transf 44 () [7] M Jourhmane, NS Mera, n iterative algorithm for the backward heat conduction problem based on variable relaxtion factors, Inverse Probl Eng () [8] NS Mera, The method of fundamental solutions for the backward heat conduction problem, Inverse Probl Sci Eng 3 () (5) [9] N Hào, mollification method for ill-posed problems, Numer Math 68 (994) [] SM Kirkup, M Wadsworth, Solution of inverse diffusion problems by operator-splitting methods, ppl Math Model 6 () () 3 8 [] L Eldén, Numerical solution of the sideways heat equation by difference approximation in time, Inverse Probl (995) [] RE Ewing, The approximation of certain parabolic equations backward in time by Sobolev equations, SIM J Math nal 6 (975) [3] JR Cannon, The One-imensional Heat Equation, ddison-wesley Publishing Company, California, 984 [4] LC Evans, Partial ifferential Equations, merican Mathematical Society, Providence, Rhode Island, 998 [5] CS Liu, Group preserving scheme for backward heat conduction problems, Int J Heat Mass Transf 47 (4) [6] K mes, LE Payne, Continuous dependence on modeling for some well-posed perturbations of the backward heat equation, J Inequal ppl 3 (999) 5 64
Initial inverse problem in heat equation with Bessel operator
International Journal of Heat and Mass Transfer 45 (22) 2959 2965 www.elsevier.com/locate/ijhmt Initial inverse problem in heat equation with Bessel operator Khalid Masood *, Salim Messaoudi, F.D. Zaman
More informationWavelets and regularization of the Cauchy problem for the Laplace equation
J. Math. Anal. Appl. 338 008440 1447 www.elsevier.com/locate/jmaa Wavelets and regularization of the Cauchy problem for the Laplace equation Chun-Yu Qiu, Chu-Li Fu School of Mathematics and Statistics,
More informationNumerical solution of hyperbolic heat conduction in thin surface layers
International Journal of Heat and Mass Transfer 50 (007) 9 www.elsevier.com/locate/ijhmt Numerical solution of hyperbolic heat conduction in thin surface layers Tzer-Ming Chen * Department of Vehicle Engineering,
More informationA WAVELET METHOD FOR SOLVING BACKWARD HEAT CONDUCTION PROBLEMS
Electronic Journal of Differential Equations, Vol. 7 (7), No. 9, pp. 9. ISSN: 7-669. URL: http://ejde.math.tstate.edu or http://ejde.math.unt.edu A WAVELET METHOD FOR SOLVING BACKWARD HEAT CONDUCTION PROBLEMS
More informationQuasireversibility Methods for Non-Well-Posed Problems
Kennesaw State University DigitalCommons@Kennesaw State University Faculty Publications 11-29-1994 Quasireversibility Methods for Non-Well-Posed Problems Gordon W. Clark Kennesaw State University Seth
More informationHaar wavelet method for nonlinear integro-differential equations
Applied Mathematics and Computation 176 (6) 34 333 www.elsevier.com/locate/amc Haar wavelet method for nonlinear integro-differential equations Ülo Lepik Institute of Applied Mathematics, University of
More informationA Fictitious Time Integration Method for Multi-Dimensional Backward Heat Conduction Problems
Copyright 2010 Tech Science Press CMC, vol.19, no.3, pp.285-314, 2010 A Fictitious Time Integration Method for Multi-Dimensional Backward Heat Conduction Problems Chih-Wen Chang 1 Abstract: In this article,
More informationExistence of positive periodic solutions for a periodic logistic equation
Applied Mathematics and Computation 139 (23) 311 321 www.elsevier.com/locate/amc Existence of positive periodic solutions for a periodic logistic equation Guihong Fan, Yongkun Li * Department of Mathematics,
More informationPast Cone Dynamics and Backward Group Preserving Schemes for Backward Heat Conduction Problems
Copyright c 2006 Tech Science Press CMES, vol12, no1, pp67-81, 2006 Past Cone Dynamics and Backward Group Preserving Schemes for Backward Heat Conduction Problems C-S Liu 1,C-WChang 2, J-R Chang 2 Abstract:
More informationThis article was published in an Elsevier journal. The attached copy is furnished to the author for non-commercial research and education use, including for instruction at the author s institution, sharing
More informationAn Energy Regularization Method for the Backward Diffusion Problem and its Applications to Image Deblurring
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 4, No. 1, pp. 177-194 Commun. Comput. Phys. July 8 An Energy Regularization Method for the Backward Diffusion Problem and its Applications to Image Deblurring
More informationA new modified Halley method without second derivatives for nonlinear equation
Applied Mathematics and Computation 189 (2007) 1268 1273 www.elsevier.com/locate/amc A new modified Halley method without second derivatives for nonlinear equation Muhammad Aslam Noor *, Waseem Asghar
More informationShooting Methods for Numerical Solution of Stochastic Boundary-Value Problems
STOCHASTIC ANALYSIS AND APPLICATIONS Vol. 22, No. 5, pp. 1295 1314, 24 Shooting Methods for Numerical Solution of Stochastic Boundary-Value Problems Armando Arciniega and Edward Allen* Department of Mathematics
More informationMATH 251 Final Examination August 14, 2015 FORM A. Name: Student Number: Section:
MATH 251 Final Examination August 14, 2015 FORM A Name: Student Number: Section: This exam has 11 questions for a total of 150 points. Show all your work! In order to obtain full credit for partial credit
More informationChebyshev pseudospectral method for computing numerical solution of convection diffusion equation
Available online at wwwsciencedirectcom Applied Mathematics and Computation (8) 537 546 wwwelseviercom/locate/amc Chebyshev pseudospectral method for computing numerical solution of convection diffusion
More informationDiffusion on the half-line. The Dirichlet problem
Diffusion on the half-line The Dirichlet problem Consider the initial boundary value problem (IBVP) on the half line (, ): v t kv xx = v(x, ) = φ(x) v(, t) =. The solution will be obtained by the reflection
More informationTwo-parameter regularization method for determining the heat source
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 8 (017), pp. 3937-3950 Research India Publications http://www.ripublication.com Two-parameter regularization method for
More informationStabilization of a solid propellant rocket instability by state feedback
INTERNATIONAL JOURNAL OF ROBUST AND NONLINEAR CONTROL Int. J. Robust Nonlinear Control 23; 13:483 495 (DOI: 1.12/rnc.732) Stabilization of a solid propellant rocket instability by state feedback Dejan
More informationPreservation of local dynamics when applying central difference methods: application to SIR model
Journal of Difference Equations and Applications, Vol., No. 4, April 2007, 40 Preservation of local dynamics when applying central difference methods application to SIR model LIH-ING W. ROEGER* and ROGER
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationNewton-homotopy analysis method for nonlinear equations
Applied Mathematics and Computation 188 (2007) 1794 1800 www.elsevier.com/locate/amc Newton-homotopy analysis method for nonlinear equations S. Abbasbandy a, *, Y. Tan b, S.J. Liao b a Department of Mathematics,
More informationDetermination of a source term and boundary heat flux in an inverse heat equation
ISSN 1746-7659, England, UK Journal of Information and Computing Science Vol. 8, No., 013, pp. 3-114 Determination of a source term and boundary heat flux in an inverse heat equation A.M. Shahrezaee 1
More informationA null space method for solving system of equations q
Applied Mathematics and Computation 149 (004) 15 6 www.elsevier.com/locate/amc A null space method for solving system of equations q Pu-yan Nie 1 Department of Mathematics, Jinan University, Guangzhou
More informationChaos synchronization of nonlinear Bloch equations
Chaos, Solitons and Fractal7 (26) 357 361 www.elsevier.com/locate/chaos Chaos synchronization of nonlinear Bloch equations Ju H. Park * Robust Control and Nonlinear Dynamics Laboratory, Department of Electrical
More informationA two level finite difference scheme for one dimensional PennesÕ bioheat equation
Applied Mathematics and Computation 171 (5) 3 331 wwwelseviercom/locate/amc A two level finite difference scheme for one dimensional PennesÕ bioheat equation Jennifer J Zhao a,1, Jun Zhang b, *,, Ning
More informationMATH 251 Final Examination May 3, 2017 FORM A. Name: Student Number: Section:
MATH 5 Final Examination May 3, 07 FORM A Name: Student Number: Section: This exam has 6 questions for a total of 50 points. In order to obtain full credit for partial credit problems, all work must be
More informationPARTIAL DIFFERENTIAL EQUATIONS. Lecturer: D.M.A. Stuart MT 2007
PARTIAL DIFFERENTIAL EQUATIONS Lecturer: D.M.A. Stuart MT 2007 In addition to the sets of lecture notes written by previous lecturers ([1, 2]) the books [4, 7] are very good for the PDE topics in the course.
More informationTwo new regularization methods for solving sideways heat equation
Nguyen and uu Journal of Inequalities and Applications 015) 015:65 DOI 10.1186/s13660-015-0564-0 R E S E A R C H Open Access Two new regularization methods for solving sideways heat equation Huy Tuan Nguyen
More informationChapter 10: Partial Differential Equations
1.1: Introduction Chapter 1: Partial Differential Equations Definition: A differential equations whose dependent variable varies with respect to more than one independent variable is called a partial differential
More informationAn inverse problem for Helmholtz equation
Inverse Problems in Science and Engineering Vol. 9, No. 6, September, 839 854 An inverse problem for Helmholtz equation M. Tadi a *, A.K. Nandakumaran b and S.S. Sritharan c a Department of Mechanical
More informationMATH 251 Final Examination December 19, 2012 FORM A. Name: Student Number: Section:
MATH 251 Final Examination December 19, 2012 FORM A Name: Student Number: Section: This exam has 17 questions for a total of 150 points. In order to obtain full credit for partial credit problems, all
More informationArbitrary order Krylov deferred correction methods for differential algebraic equations
Journal of Computational Physics 221 (2007) 739 760 www.elsevier.com/locate/jcp Arbitrary order Krylov deferred correction methods for differential algebraic equations Jingfang Huang *,1, Jun Jia 1, Michael
More informationDynamic scaling in a simple one-dimensional model of dislocation activity
Philosophical Magazine, 11 August 24 Vol. 84, No. 23, 2445 2454 Dynamic scaling in a simple one-dimensional model of dislocation activity Jack Deslippey, R. Tedstromy, Murray S. Daw*y, D. Chrzanz}, T.
More informationMathematics (JUN11MPC301) General Certificate of Education Advanced Level Examination June Unit Pure Core TOTAL
Centre Number Candidate Number For Examiner s Use Surname Other Names Candidate Signature Examiner s Initials Mathematics Unit Pure Core 3 Monday 13 June 2011 General Certificate of Education Advanced
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationJournal of Computational Physics
Journal of Computational Physics 9 () 759 763 Contents lists available at ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp Short Note A comment on the computation
More information75RI ARTICLE IN PRESS. A. Romeo a,, G. Finocchio a, M. Carpentieri a, L. Torres b, G. Consolo a, B. Azzerboni a
7RI Physica B 43 (28) 464 468 www.elsevier.com/locate/physb A numerical solution of the magnetization reversal modeling in a permalloy thin film using fifth order Runge Kutta method with adaptive step
More informationNumerical comparison of methods for solving linear differential equations of fractional order
Chaos, Solitons and Fractals 31 (27) 1248 1255 www.elsevier.com/locate/chaos Numerical comparison of methods for solving linear differential equations of fractional order Shaher Momani a, *, Zaid Odibat
More informationCOMPUTATIONAL METHODS AND ALGORITHMS Vol. II - Numerical Algorithms for Inverse and Ill-Posed Problems - A.М. Denisov
NUMERICAL ALGORITHMS FOR INVERSE AND ILL-POSED PROBLEMS Faculty of Computational Mathematics and Cybernetics, Moscow State University, Moscow, Russia Keywords: Direct problem, Inverse problem, Well-posed
More informationAn interior point type QP-free algorithm with superlinear convergence for inequality constrained optimization
Applied Mathematical Modelling 31 (2007) 1201 1212 www.elsevier.com/locate/apm An interior point type QP-free algorithm with superlinear convergence for inequality constrained optimization Zhibin Zhu *
More informationSequential dynamical systems over words
Applied Mathematics and Computation 174 (2006) 500 510 www.elsevier.com/locate/amc Sequential dynamical systems over words Luis David Garcia a, Abdul Salam Jarrah b, *, Reinhard Laubenbacher b a Department
More informationA discontinuous Galerkin finite element method for directly solving the Hamilton Jacobi equations
Journal of Computational Physics 3 (7) 39 415 www.elsevier.com/locate/jcp A discontinuous Galerkin finite element method for directly solving the Hamilton Jacobi equations Yingda Cheng, Chi-Wang Shu *
More informationNonlinear transient heat conduction problems for a class of inhomogeneous anisotropic materials by BEM
Engineering Analysis with Boundary Elements 32 (28) 154 16 www.elsevier.com/locate/enganabound Nonlinear transient heat conduction problems for a class of inhomogeneous anisotropic materials by BEM M.I.
More informationA multiple Riccati equations rational expansion method and novel solutions of the Broer Kaup Kupershmidt system
Chaos, Solitons and Fractals 30 (006) 197 03 www.elsevier.com/locate/chaos A multiple Riccati equations rational expansion method and novel solutions of the Broer Kaup Kupershmidt system Qi Wang a,c, *,
More informationA numerical analysis of a model of growth tumor
Applied Mathematics and Computation 67 (2005) 345 354 www.elsevier.com/locate/amc A numerical analysis of a model of growth tumor Andrés Barrea a, *, Cristina Turner b a CIEM, Universidad Nacional de Córdoba,
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat xxx ) xxx xxx Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: wwwelseviercom/locate/cnsns A note on the use of Adomian
More informationLocally divergence-free discontinuous Galerkin methods for the Maxwell equations
Journal of Computational Physics 194 (4) 588 61 www.elsevier.com/locate/jcp Locally divergence-free discontinuous Galerkin methods for the Maxwell equations Bernardo Cockburn a,1, Fengyan Li b,, Chi-Wang
More informationPARAREAL TIME DISCRETIZATION FOR PARABOLIC CONTROL PROBLEM
PARAREAL TIME DISCRETIZATION FOR PARABOLIC CONTROL PROBLEM Daoud S Daoud Dept. of Mathematics Eastern Mediterranean University Famagusta-North Cyprus via Mersin 10-Turkey. e mail: daoud.daoud@emu.edu.tr
More informationApplied Mathematics and Computation
Applied Mathematics and Computation 245 (2014) 86 107 Contents lists available at ScienceDirect Applied Mathematics and Computation journal homepage: www.elsevier.com/locate/amc An analysis of a new family
More informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationTutorial 2. Introduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2
More information1 Basic tools. 1.1 The Fourier transform. Synopsis. Evenness and oddness. The Fourier transform of f(x) is defined as Z 1
1 Basic tools 1.1 The Fourier transform 1 Synopsis The Fourier transform of f(x) is defined as FðsÞ ¼ f ðxþe i2pxs dx The inverse Fourier transform is given by f ðxþ ¼ FðsÞe þi2pxs ds Evenness and oddness
More informationMATH 251 Final Examination December 16, 2015 FORM A. Name: Student Number: Section:
MATH 5 Final Examination December 6, 5 FORM A Name: Student Number: Section: This exam has 7 questions for a total of 5 points. In order to obtain full credit for partial credit problems, all work must
More informationMATH 614 Dynamical Systems and Chaos Lecture 3: Classification of fixed points.
MATH 614 Dynamical Systems and Chaos Lecture 3: Classification of fixed points. Periodic points Definition. A point x X is called a fixed point of a map f : X X if f(x) = x. A point x X is called a periodic
More informationBoundedness of solutions to a retarded Liénard equation
Electronic Journal of Qualitative Theory of Differential Equations 21, No. 24, 1-9; http://www.math.u-szeged.hu/ejqtde/ Boundedness of solutions to a retarded Liénard equation Wei Long, Hong-Xia Zhang
More informationApplied Mathematics and Computation 169 (2005)
Applied Mathematics and Computation 169 (2005) 485 499 www.elsevier.com/locate/amc Solving Laplacian problems with boundary singularities: a comparison of a singular function boundary integral method with
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationNON-HYPERELLIPTIC RIEMANN SURFACES OF GENUS FIVE ALL OF WHOSE WEIERSTRASS POINTS HAVE MAXIMAL WEIGHT1. Ryutaro Horiuchi. Abstract
R HORIUCHI KODAI MATH J 0 (2007), 79 NON-HYPERELLIPTIC RIEMANN SURFACES OF GENUS FIVE ALL OF WHOSE WEIERSTRASS POINTS HAVE MAXIMAL WEIGHT1 Ryutaro Horiuchi Abstract It is proved that any non-hyperelliptic
More informationNumerical solution of a Fredholm integro-differential equation modelling _ h-neural networks
Available online at www.sciencedirect.com Applied Mathematics and Computation 195 (2008) 523 536 www.elsevier.com/locate/amc Numerical solution of a Fredholm integro-differential equation modelling _ h-neural
More informationApplied Mathematics and Computation
Applied Mathematics and Computation 8 (0) 584 599 Contents lists available at SciVerse ScienceDirect Applied Mathematics and Computation journal homepage: www.elsevier.com/locate/amc Basin attractors for
More informationLecture Introduction
Lecture 1 1.1 Introduction The theory of Partial Differential Equations (PDEs) is central to mathematics, both pure and applied. The main difference between the theory of PDEs and the theory of Ordinary
More informationQuantum lattice representations for vector solitons in external potentials
Physica A 362 (2006) 215 221 www.elsevier.com/locate/physa Quantum lattice representations for vector solitons in external potentials George Vahala a,, Linda Vahala b, Jeffrey Yepez c a Department of Physics,
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat xxx (2009) xxx xxx Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns A one-step optimal homotopy
More informationJournal of Computational Physics
Journal of Computational Physics 9 (010) 3019 3045 Contents lists available at ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp BDF-like methods for nonlinear
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationConvergence rate estimates for the gradient differential inclusion
Convergence rate estimates for the gradient differential inclusion Osman Güler November 23 Abstract Let f : H R { } be a proper, lower semi continuous, convex function in a Hilbert space H. The gradient
More informationStochastic approximations of perturbed Fredholm Volterra integro-differential equation arising in mathematical neurosciences
Applied athematics and Computation 6 (27) 73 2 www.elsevier.com/locate/amc Stochastic approximations of perturbed Fredholm Volterra integro-differential equation arising in mathematical neurosciences.
More informationEvolution equations with spectral methods: the case of the wave equation
Evolution equations with spectral methods: the case of the wave equation Jerome.Novak@obspm.fr Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris, France in collaboration with
More informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
More informationTHE POINT SPECTRUM OF FROBENIUS-PERRON AND KOOPMAN OPERATORS
PROCEEDINGS OF THE MERICN MTHEMTICL SOCIETY Volume 126, Number 5, May 1998, Pages 1355 1361 S 0002-9939(98)04188-4 THE POINT SPECTRUM OF FROBENIUS-PERRON ND KOOPMN OPERTORS J. DING (Communicated by Palle
More informationAsymptotics of minimax stochastic programs $
Statistics & Probability Letters 78 (2008) 150 157 www.elsevier.com/locate/stapro Asymptotics of minimax stochastic programs $ Alexander Shapiro School of Industrial and Systems Engineering, Georgia Institute
More informationJournal of Computational Physics
Journal of Computational Physics 229 (2010 7747 7764 Contents lists available at ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp A unified gas-kinetic scheme
More informationDETERMINATION OF AN UNKNOWN SOURCE TERM IN A SPACE-TIME FRACTIONAL DIFFUSION EQUATION
Journal of Fractional Calculus and Applications, Vol. 6(1) Jan. 2015, pp. 83-90. ISSN: 2090-5858. http://fcag-egypt.com/journals/jfca/ DETERMINATION OF AN UNKNOWN SOURCE TERM IN A SPACE-TIME FRACTIONAL
More informationSOURCE IDENTIFICATION FOR THE HEAT EQUATION WITH MEMORY. Sergei Avdonin, Gulden Murzabekova, and Karlygash Nurtazina. IMA Preprint Series #2459
SOURCE IDENTIFICATION FOR THE HEAT EQUATION WITH MEMORY By Sergei Avdonin, Gulden Murzabekova, and Karlygash Nurtazina IMA Preprint Series #2459 (November 215) INSTITUTE FOR MATHEMATICS AND ITS APPLICATIONS
More informationNon-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3
Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Tommaso Ruggeri Department of Mathematics and Research Center of Applied Mathematics University of Bologna January 21, 2017 ommaso
More informationAn inverse vibration analysis of a tower subjected to wind drags on a shaking ground
Applied Mathematical Modelling 26 (2002) 517 528 www.elsevier.com/locate/apm An inverse vibration analysis of a tower subjected to wind drags on a shaking ground C.C. Kang a, C.Y. Lo b, * a Department
More informationA MODIFIED QUASI-BOUNDARY VALUE METHOD FOR A CLASS OF ABSTRACT PARABOLIC ILL-POSED PROBLEMS
A MODIFIED QUASI-BOUNDARY VALUE METHOD FOR A CLASS OF ABSTRACT PARABOLIC ILL-POSED PROBLEMS M. DENCHE AND S. DJEZZAR Received 14 October 24; Accepted 9 August 25 We study a final value problem for first-order
More informationMATH 251 Final Examination December 16, 2014 FORM A. Name: Student Number: Section:
MATH 2 Final Examination December 6, 204 FORM A Name: Student Number: Section: This exam has 7 questions for a total of 0 points. In order to obtain full credit for partial credit problems, all work must
More informationExamination paper for TMA4130 Matematikk 4N: SOLUTION
Department of Mathematical Sciences Examination paper for TMA4 Matematikk 4N: SOLUTION Academic contact during examination: Morten Nome Phone: 9 84 97 8 Examination date: December 7 Examination time (from
More informationSEMIGROUP APPROACH FOR PARTIAL DIFFERENTIAL EQUATIONS OF EVOLUTION
SEMIGROUP APPROACH FOR PARTIAL DIFFERENTIAL EQUATIONS OF EVOLUTION Istanbul Kemerburgaz University Istanbul Analysis Seminars 24 October 2014 Sabanc University Karaköy Communication Center 1 2 3 4 5 u(x,
More informationThe Heat Equation John K. Hunter February 15, The heat equation on a circle
The Heat Equation John K. Hunter February 15, 007 The heat equation on a circle We consider the diffusion of heat in an insulated circular ring. We let t [0, ) denote time and x T a spatial coordinate
More informationFINAL EXAM, MATH 353 SUMMER I 2015
FINAL EXAM, MATH 353 SUMMER I 25 9:am-2:pm, Thursday, June 25 I have neither given nor received any unauthorized help on this exam and I have conducted myself within the guidelines of the Duke Community
More informationMAXIMUM PRINCIPLES FOR THE RELATIVISTIC HEAT EQUATION
MAXIMUM PRINCIPLES FOR THE RELATIVISTIC HEAT EQUATION EVAN MILLER AND ARI STERN arxiv:1507.05030v1 [math.ap] 17 Jul 2015 Abstract. The classical heat equation is incompatible with relativity, since the
More informationApplied Mathematical Modelling
Applied Mathematical Modelling 33 (29) 2293 2297 Contents lists available at ScienceDirect Applied Mathematical Modelling journal homepage: www.elsevier.com/locate/apm Stability criteria for a nonlinear
More informationMATH 241 Practice Second Midterm Exam - Fall 2012
MATH 41 Practice Second Midterm Exam - Fall 1 1. Let f(x = { 1 x for x 1 for 1 x (a Compute the Fourier sine series of f(x. The Fourier sine series is b n sin where b n = f(x sin dx = 1 = (1 x cos = 4
More informationNon-linear dynamic analysis of beams with variable stiffness
Journal of Sound and Vibration 27 (24) 847 863 JOURNAL OF SOUND AND VIBRATION www.elsevier.com/locate/jsvi Non-linear dynamic analysis of beams with variable stiffness J.T. Katsikadelis*, G.C. Tsiatas
More informationRounding Error in Numerical Solution of Stochastic Differential Equations
STOCHASTIC ANALYSIS AND APPLICATIONS Vol. 21, No. 2, pp. 281 300, 2003 Rounding Error in Numerical Solution of Stochastic Differential Equations Armando Arciniega* and Edward Allen Department of Mathematics
More informationPierre Lissy. 19 mars 2015
Explicit lower bounds for the cost of fast controls for some -D parabolic or dispersive equations, and a new lower bound concerning the uniform controllability of the -D transport-diffusion equation Pierre
More informationMATH 251 Final Examination May 4, 2015 FORM A. Name: Student Number: Section:
MATH 251 Final Examination May 4, 2015 FORM A Name: Student Number: Section: This exam has 16 questions for a total of 150 points. In order to obtain full credit for partial credit problems, all work must
More informationDeterminants and polynomial root structure
International Journal of Mathematical Education in Science and Technology, Vol 36, No 5, 2005, 469 481 Determinants and polynomial root structure L G DE PILLIS Department of Mathematics, Harvey Mudd College,
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationImproved Lebesgue constants on the triangle
Journal of Computational Physics 207 (2005) 625 638 www.elsevier.com/locate/jcp Improved Lebesgue constants on the triangle Wilhelm Heinrichs Universität Duisburg-Essen, Ingenieurmathematik (FB 10), Universitätsstrasse
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat 7 (0) 3499 3507 Contents lists available at SciVerse ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns Application of the
More informationDipartimento di Matematica, Politecnico di Torino, Torino 10129, Italy PLEASE SCROLL DOWN FOR ARTICLE
This article was downloaded by: [Camporesi, Roberto] On: 19 May 211 Access details: Access Details: [subscription number 936142537] Publisher Taylor & Francis Informa Ltd Registered in England and Wales
More informationMAIN ARTICLES. In present paper we consider the Neumann problem for the operator equation
Volume 14, 2010 1 MAIN ARTICLES THE NEUMANN PROBLEM FOR A DEGENERATE DIFFERENTIAL OPERATOR EQUATION Liparit Tepoyan Yerevan State University, Faculty of mathematics and mechanics Abstract. We consider
More informationAnalysis of a Chemotaxis Model for Multi-Species Host-Parasitoid Interactions
Applied Mathematical Sciences, Vol. 2, 2008, no. 25, 1239-1252 Analysis of a Chemotaxis Model for Multi-Species Host-Parasitoid Interactions Xifeng Tang and Youshan Tao 1 Department of Applied Mathematics
More informationand finally, any second order divergence form elliptic operator
Supporting Information: Mathematical proofs Preliminaries Let be an arbitrary bounded open set in R n and let L be any elliptic differential operator associated to a symmetric positive bilinear form B
More informationNontrivial Solutions for Boundary Value Problems of Nonlinear Differential Equation
Advances in Dynamical Systems and Applications ISSN 973-532, Volume 6, Number 2, pp. 24 254 (2 http://campus.mst.edu/adsa Nontrivial Solutions for Boundary Value Problems of Nonlinear Differential Equation
More informationARTICLE IN PRESS Mathematical and Computer Modelling ( )
Mathematical and Computer Modelling Contents lists available at ScienceDirect Mathematical and Computer Modelling ournal homepage: wwwelseviercom/locate/mcm Total variation diminishing nonstandard finite
More informationEstimation of carbon transfer coefficients using Duke Forest free-air CO 2 enrichment data
Applied Mathematics and Computation 1 () 11 1 wwwelseviercom/locate/amc Estimation of carbon transfer coefficients using Duke Forest free-air CO enrichment data Luther White a, *, Yiqi Luo b a Department
More information