Transactions on Modelling and Simulation vol 8, 1994 WIT Press, ISSN X
|
|
- Kerrie Stevens
- 5 years ago
- Views:
Transcription
1 Boundary element method for an improperly posed problem in unsteady heat conduction D. Lesnic, L. Elliott & D.B. Ingham Department of Applied Mathematical Studies, University of Leeds, Leeds LS2 9JT, UK Abstract In this one-dimensional, unsteady, inverse heat conduction problem (IHCP) the boundary conditions are enforced at time t = 0, at x = d (0 < d 3 1) and x = 1 and the solution is sought using the Boundary Element Method (BEM). Because this is an improperly posed problem the direct method cannot be applied and special corrective procedures, such as least squares and minimal energy methods, have been introduced. From all the numerical results performed it has been found that both the least squares and the minimal energy methods give good agreement with the analytical solutions provided that the known boundary condition is extended over a suitable additional time interval such that accurate results can be obtained over the complete earlier time domain. However, when a sufficient amount of noise is included in the data then the least squares method produced substantially less accurate and oscillatory results. In these cases, the minimal energy technique overcomes these problems by minimizing the kinetic energy functional. 1 Introduction Inverse heat conduction problems (IHCP) arise in many heat transfer situations when experimental difficulties are uncounted in measuring or producing the appropriate boundary conditions. In practice it is often difficult to measure the surface temperature or the heat flux at the surface of a solid as a function of time. Sometimes it is easier to record the temperature at one or more interior locations within the solid. This is only one of the many examples of inverse heat conduction problems discussed in the excellent review by Beck et al. [11. In this study the solution of the one-dimensional, unsteady linear heat conduction equation in a slab geometry with constant physical properties, namely
2 78 Boundary Element Technology is analysed. The governing eqn (1) has to be solved subject to the boundary conditions T(x,0) = T (x) for x [0,1] (2) o a q(l,t) + p T(l,t) = f(t) for t e (0,oo) (3) where T is the temperature, q = dt/dx is the heat flux, T and f are prescribed functions and a and /3 are known constants which are not simultaneously zero. In addition, temperature readings are provided at an arbitrary space location x = d (0,1], namely T(x,t) = g(t) for x = d (4) where g is a known function. Based on the formulation in eqns (l)-(4) the IHCP considered here requires the determination only of the temperature T(0,t) and the heat flux q(0,t) on the remaining part of the boundary, i.e. x = 0. In practice only a finite time domain is analysed such that the time t (0,t ], where t is an arbitrary specified final time of interest. 2 The boundary element method The BEM for the numerical solution of this class of improperly posed problems is applied since this choice of method does not require any domain discretisation. In addition, the existence of a fundamental solution for the governing differential equation enables the problem to be reformulated in an integral representation form involving only the boundary integrals. The one-dimensional fundamental solution of the heat conduction eqn (1) is of the form, see Brebbia et al. [2], ^\2 1/7 exp - 4- I J^L_ t T 1 H(t-%) (5) [47l(t-T)]^ I 4lt TJ where H is the Heaviside function and this enables the differential eqn (1) to be transformed into the following boundary integral equation i 7)(x,t) T(x,t) = f T'(,T) F(x,t;g,%) dx - f T(,T) F' dt Jo o ^0 0! T(y,0) F(x,t;y,0) dy (6) r. 'o where (x,t) [0,l]x[0,oo), primes denote differentiation with respect to the outward normal, TJ = 1 for (x,t) <E (0,l)x(0,oo), and 17 = 0.5 for x = 0, 1 or t = 0. In eqn (6), T' = q on = 1 and T' = -q on = 0. For simplicity, in the formulation of the BEM applied here constant boundary elements, i.e. the temperature and the heat flux are assumed constant over each element, are used when the integral eqn (6) is discretised. The time interval [0,t ] is divided into N elements on each boundary x = 0 and x = 1 and the space interval [0,1] into N elements. If one takes x on the boundary, at x = 0 and x = 1, then the integral eqn (6) results in a system of linear algebraic equations, namely, (1)
3 Boundary Element Technology 79 2N 2 Mr 1 y G(x,t; ;t,t ) T' -y E(x,t ;,t.,t.) + 7). 5 I T. L I' l'^' j-l' J J L[ I. j j-ij i uj j N 0 H(x,t;y,y ) T = 0, i = 1,2N (7) j=i ' l J-' J ^ where 77 = 0.5, 6 is the Kronecker symbol, t _^ and t. are the endpoints of an element on the boundaries x = 0 and x = 1 and t^ its midpoint, y t = 0, T, T' =-q and y, are the endpoints of an element on the boundary for j = T7N and T' = q. for j = N+1,2N are the values of the temperature T and heat flux ±q at the midpoint of the element [t,t ], T takes the value of T at the midpoint of the j-l j Oj 0 element [y._,y.], x. = 0 for i < N, x = 1 for i > N,?. = 0 for j < N,, = 1 for j > N and the influence matrices are defined as ^b ~(x.t:.t) H(t-x) dx (8) G(x,t;g;a,b) = [ r^ E(x,t; ;a,b) = F' (x,t;?,x) H(t-x) dx (9) H(x,t;c,d) = [ F(x,t;y,0) H(t) dy (10) The discretised form of the boundary condition (3) is a q. + 3 T. = f., i = N + 1.2N (11) where f = f(t ), and in order to complete the system of equations one must add further information represented by the internal measurements (4). The time interval [0,t ] on the space location x = d is further discretised into N elements with t' denoting the midpoint of the i-th element on this interval. Applying the discretised form of equation (4) into the boundary integral (6} results in another N linear equations, namely 2N 2N G(d,t';g;t,t)T' -y E(d,t';^,t,t)T. i J J-l J J L i J J-l J J J=l J=l N 0 V H(d,t'.;y.,y.) T. = 7)(d,t'. ) g., i = 1,N (12) Lt i j-l j oj i i * j-i where g = g(t'). At this stage one recalls that the application of BEM for the IHCP produced an algebraic system (7), (11) and (12) of (3N + NM linear equations with 4N unknowns. In particular, eqns (7) and (11) enable some of the unknowns, preferably the discretised temperatures T. for j = N+1,2N, i.e. on the boundary x = 1, and the heat fluxes q. for j = 1,2N, i.e. on the boundaries x = 0 and x = 1, to be expressed
4 80 Boundary Element Technology as a function of the remaining discretised temperatures T for j = i,n, i.e. on the boundary x = 0. These substitutions, when introduced into eqns (12), lead to a lower order system of N linear equations with N unknowns, which in a generic form can be written in the form A T = b (13) where A = (A ) for i = 1,N, j = 1,N is a matrix depending on the geometrical matrices introduced in expressions (8)-(10), b = (by) for i = IN is a known vector and T = (T ) for j = 1,N is the unknown T - j vector of the temperature on the boundary x = 0. So far, the application of the BEM and boundary conditions (2)-(4) has reduced the IHCP to the linear system of algebraic eqns (13). A necessary condition for a solution is that N ^ N. Taking N = N, then the system of eqns (13) contains N linear equations with N unknowns and, if the matrix A has an invers, then one simply has T = A~* b. However, because the IHCP (l)-(4) is improperly posed the system of eqns (13) is ill-conditioned, i.e. the matrix A is singular, and hence no direct solution is possible. Since the system of eqns (13) is ill-conditioned more information is necessary and this is achieved from the known data (4) by taking N > N and the system becomes overdetermined. Using the least-squares method, based on the minimization of the Euclidian norm inf A T - b I ^, the solution of this system is given by, see Pasquetti & Le Niliot [3] T = (A^ A)~* A^ b (14) However, it has been found that some of the results are not sufficiently accurate nor stable and convergent. Therefore we replace the exact internal condition (4) with the inequality T(d,t) - g(t) * c for t e (0,oo) (15) where e > 0 is a preassigned small quantity and then the system of eqns (11) becomes a system of N linear inequalities with N unknowns, namely, A T - b * c (16) where the value of c > 0 should be chosen as small as possible. The system of inequalities (26) will form the set of constraints which should be satisfied by an energy functional introduced in the next section. 3 Minimal Energy Method In order to improve the stability of the solution with respect to large experimental errors, the minimal energy technique has been considered. This method was developed by Han et al. [4] for the backward heat conduction problem and it has been extended to the present IHCP in this paper. A physical argument is now adopted, in an analogy with the
5 Boundary Element Technology 81 steady case, and consequently one minimizes the kinetic energy - f' f' I %^1'dxdT (17) ^ ^ In order to put the relation (17) into a boundary integral form the energy equation is derived by multiplying the governing eqn (1) by u(x,t) and integrating over the domain [0,l]x[0,t^], namely, _ r'f T* 9T(x,%),2 ^ ^ ^ 1 r ^(^ )dx - 1 [ f (x,0)dx (18) Jo Jo ^^ ^o * ^o where the left hand side denotes the absorbed (or released) thermal energy through both ends x = 0 and x = 1 and the kinetic energy KE(T), whilst the right hand side denotes the inner energy of the system at times t = 0 and t = t. Introducing the identity (18) into relation (17) after discretisation yields N 2N, 0 _. _ KE(T) = U - t^j T. r + I[ <y, - y^ - T <y,t,» (19) where y. is the midpoint of the i-th element on t = t^. the boundary 2N integral (6) the terms T(y.,t ), can be expressed as 2N G(y,t;^;t,t)T' -Y E(y,t;g.,t,t)T. i f j j-l j j L i f J j-l J J According to N 0 _ + V L H(y,t ;y,y ) T = n(y,t ) T(y.,t ), i = 1,N I f j-l j oj if 11 u (20) The process of substitution mentioned in section 2, when introduced into expression (20), results in a quadratic form for the kinetic energy as a function of the unknowns temperatures T = (T.) for j = Ij\J which can be written in the following generic form KE(T) = T'" Q T + q'" T + q^ (21) where Q = (QJ for i, j = T^N, q^ = (q^) for i = TIN and q^ are known matrix, vector and constant, respectively, depending on the geometry and the known boundary conditions (2) and (3). The minimal energy method now reduces to finding the temperature vector T which minimizes the kinetic energy (21), subject to the linear constraints of the type (16). This problem is solved using the NAG routine E04UCF based on the minimization of an arbitrary smooth function subject to certain constraints which may include simple bounds of variables and both linear and nonlinear constraints. 4 Numerical results and discussion As expected, the IHCP is ill-posed and therefore the direct method was found to be inappropriate since the system of eqns (13) is ill-conditioned. For illustrating various features of the numerical
6 82 Boundary Element Technology schemes a simple the test function is considered, namely, T(x,t) = 2t + x* (22) The function T satisfies equation (1), together with the general boundary condition at x = 1, namely a = = 1, and the internal condition is taken on the boundary, i.e. d = 1. In all the numerical calculations presented herein the mesh size was fixed at N = N = 40, N =80 and the small control quantity c was chosen to be O(10 ), so as to ensure that the results were accurate. In the many examples performed the numerical results for the surface temperature were found to agree well with the analytical solution and therefore, in this study the results are presented only for the surface heat flux which is more difficult to accurately calculate. The effect of the choice of the value of t is illustrated by solving the problem in the domain [0,l]x[0,l] and in the domain [0,l]x[0,2] with the numerical results for the heat flux on the boundary x = 0 shown in tables la and Ib, respectively. From table la it can be seen that the numerical results for the surface heat flux obtained by using the minimal energy methods are very similar to the analytical solution up to a value of approximately 0.9 t. However, inaccurate and oscillatory results in predicting the surface heat flux are strongly present when using the least squares method. However, from table Ib, where t = 2, it can be observed that the numerical results for the surface heat flux obtained by using both the least squares and minimal energy methods are sufficiently close to the analytical solution for t ^ t /2 = 1. Table la. Surface heat flux for Table Ib. Surface heat flux for t Leas t squares Mini ma 1 energy Eqn. (22) t Least squares Mini ma 1 energy Eqn. (22) Hence, it is concluded that if [0,t ] is the time interval of f interest, then by solving the problem on an extended time interval results which are in good agreement with the analytical solution over
7 Boundary Element Technology 83 the earlier time interval are obtained. At this stage, because IHCP is an improperly posed problem it is necessary to investigate the stability of the numerical schemes under discussion. Therefore some noise is introduced into both the internal and boundary data through a small perturbation function f(x,t). Here we take T (x,t) = T(x,t) + f (x,t) (23) i where the perturbation function T(x,t) is given by i i f(x,t) = %(?) /2 cos( 2?^ t - y x - ^ ) exp( -? x ) (24) and? = k^ T^, %(y) = 0 if k = 0, %(y) = 1 if k * 0 and k_ is a constant which has to be prescribed. Clearly, the values T(x,0), T(l,t), q(l,t), T(d,t) can be made as small as one wishes by increasing the value of y. However, at the surface x = 0, T(0,t) and q(0,t) are of order of unity and thus the inequality (15) contains solutions which are discontinuously dependent on the initial data and consequently the problem is ill-posed. Tables 2a and 2b show the numerical results for the surface heat flux obtained by solving the problem with t = 2 and for k = 0 and 3, respectively, and the internal condition (15) is imposed at d = It is observed that the results for the surface heat flux differ slightly, but are substantially different from the large values of the heat flux given by the perturbed function T. Consequently, the least squares and the minimal energy methods are able to recognize the correct solution when a small amount of noise is included into the input data. Table 2a. Surface heat flux for Table 2b. Surface heat flux for k = 0. k = 3. t Least squares Mini ma i energy Eqn. (22) t Least squares Mini mal energy Eqn. (23) Figure 1 shows the numerical results for the surface heat flux obtained by using the numerical methods when a large amount of noise is included into the data in comparison with the unperturbed solution (22). This large noise is simulated by taking 7 = 2n, i.e.
8 84 Boundary Element Technology k = V 2 in expression (24) which defines the perturbation function T. Some of the numerical results for the surface heat flux exhibit oscillatory behaviour when using the least squares method but these large oscillations are substantially reduced when using the minimal energy technique q<o,t) I I I 1 I Figure 1: The surface heat flux when there is a large amount of noise, y = 2n in the input data, where analytical solution, a a a the least squares method, «the minimal energy technique. REFERENCES [11 Beck, J.V., Blackwell, B. & Clair JR., C.R.St. Inverse Heat Conduction, J.Wiley-Intersc.Publ., New York, [21 Brebbia, C.A., Telles, J.C.F. & Wrobel, L.C. Boundary Element Techniques, Springer-Verlag, New York, [31 Pasquetti, R. & Le Niliot, C. Boundary element approach for inverse heat conduction problems: application to a bidimensional transient numerical experiment, Num.Heat Transfer, 1991, 20, [41 Han, H., Ingham, D.B. & Yuan, Y. The boundary-element method for the solution of the backward heat conduction equation, submitted to IMA J.Appl.Math., 1993.
Conjugate gradient-boundary element method for a Cauchy problem in the Lam6 system
Conjugate gradient-boundary element method for a Cauchy problem in the Lam6 system L. Marin, D.N. Hho & D. Lesnic I Department of Applied Mathematics, University of Leeds, UK 'Hanoi Institute of Mathematics,
More informationThe solution of a nonlinear inverse problem in heat transfer
IMA Journal of Applied Mathematics (1993) 50, 113-132 The solution of a nonlinear inverse problem in heat transfer D. B. INGHAM AND Y. YUAN Department of Applied Mathematical Studies, University of Leeds,
More informationDetermination of a space-dependent force function in the one-dimensional wave equation
S.O. Hussein and D. Lesnic / Electronic Journal of Boundary Elements, Vol. 12, No. 1, pp. 1-26 (214) Determination of a space-dependent force function in the one-dimensional wave equation S.O. Hussein
More informationThe detection of subsurface inclusions using internal measurements and genetic algorithms
The detection of subsurface inclusions using internal measurements and genetic algorithms N. S. Meral, L. Elliott2 & D, B, Ingham2 Centre for Computational Fluid Dynamics, Energy and Resources Research
More information3rd Int. Conference on Inverse Problems in Engineering AN INVERSE PROBLEM FOR SLOW VISCOUS INCOMPRESSIBLE FLOWS. University of Leeds
Inverse Problems in Engineering: Theory and Practice 3rd Int. Conference on Inverse Problems in Engineering June 13-18, 1999, Port Ludlow, WA, USA ME06 AN INVERSE PROBLEM FOR SLOW VISCOUS INCOMPRESSIBLE
More informationSolving the complete-electrode direct model of ERT using the boundary element method and the method of fundamental solutions
Solving the complete-electrode direct model of ERT using the boundary element method and the method of fundamental solutions T. E. Dyhoum 1,2, D. Lesnic 1, and R. G. Aykroyd 2 1 Department of Applied Mathematics,
More informationControllability of the linear 1D wave equation with inner moving for
Controllability of the linear D wave equation with inner moving forces ARNAUD MÜNCH Université Blaise Pascal - Clermont-Ferrand - France Toulouse, May 7, 4 joint work with CARLOS CASTRO (Madrid) and NICOLAE
More informationAN EVOLUTIONARY ALGORITHM TO ESTIMATE UNKNOWN HEAT FLUX IN A ONE- DIMENSIONAL INVERSE HEAT CONDUCTION PROBLEM
Proceedings of the 5 th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15 th July 005. AN EVOLUTIONARY ALGORITHM TO ESTIMATE UNKNOWN HEAT FLUX IN A
More informationTHE WAVE EQUATION. d = 1: D Alembert s formula We begin with the initial value problem in 1 space dimension { u = utt u xx = 0, in (0, ) R, (2)
THE WAVE EQUATION () The free wave equation takes the form u := ( t x )u = 0, u : R t R d x R In the literature, the operator := t x is called the D Alembertian on R +d. Later we shall also consider the
More information2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1
Lectures - Week 11 General First Order ODEs & Numerical Methods for IVPs In general, nonlinear problems are much more difficult to solve than linear ones. Unfortunately many phenomena exhibit nonlinear
More information(1 + 2y)y = x. ( x. The right-hand side is a standard integral, so in the end we have the implicit solution. y(x) + y 2 (x) = x2 2 +C.
Midterm 1 33B-1 015 October 1 Find the exact solution of the initial value problem. Indicate the interval of existence. y = x, y( 1) = 0. 1 + y Solution. We observe that the equation is separable, and
More informationThe Finite Difference Method
Chapter 5. The Finite Difference Method This chapter derives the finite difference equations that are used in the conduction analyses in the next chapter and the techniques that are used to overcome computational
More information3rd Int. Conference on Inverse Problems in Engineering. Daniel Lesnic. University of Leeds. Leeds, West Yorkshire LS2 9JT
Inverse Problems in Engineering: Theory and Practice 3rd Int. Conference on Inverse Problems in Engineering June 13-18, 1999, Port Ludlow, WA, USA APP RETRIEVING THE FLEXURAL RIGIDITY OF A BEAM FROM DEFLECTION
More informationOrdinary Differential Equation Theory
Part I Ordinary Differential Equation Theory 1 Introductory Theory An n th order ODE for y = y(t) has the form Usually it can be written F (t, y, y,.., y (n) ) = y (n) = f(t, y, y,.., y (n 1) ) (Implicit
More informationDegenerate Perturbation Theory. 1 General framework and strategy
Physics G6037 Professor Christ 12/22/2015 Degenerate Perturbation Theory The treatment of degenerate perturbation theory presented in class is written out here in detail. The appendix presents the underlying
More informationValue Problems Introduction
Chapter 17. 17.0 Introduction Two Point Boundary Value Problems When ordinary differential equations are to satisfy conditions at more than one of the independent variable, the resulting problem is called
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationInterpolation. Chapter Interpolation. 7.2 Existence, Uniqueness and conditioning
76 Chapter 7 Interpolation 7.1 Interpolation Definition 7.1.1. Interpolation of a given function f defined on an interval [a,b] by a polynomial p: Given a set of specified points {(t i,y i } n with {t
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
More informationwe get y 2 5y = x + e x + C: From the initial condition y(0) = 1, we get 1 5 = 0+1+C; so that C = 5. Completing the square to solve y 2 5y = x + e x 5
Math 24 Final Exam Solution 17 December 1999 1. Find the general solution to the differential equation ty 0 +2y = sin t. Solution: Rewriting the equation in the form (for t 6= 0),we find that y 0 + 2 t
More informationModels for Control and Verification
Outline Models for Control and Verification Ian Mitchell Department of Computer Science The University of British Columbia Classes of models Well-posed models Difference Equations Nonlinear Ordinary Differential
More informationLecture 13 - Wednesday April 29th
Lecture 13 - Wednesday April 29th jacques@ucsdedu Key words: Systems of equations, Implicit differentiation Know how to do implicit differentiation, how to use implicit and inverse function theorems 131
More informationAlgebraic Properties of Solutions of Linear Systems
Algebraic Properties of Solutions of Linear Systems In this chapter we will consider simultaneous first-order differential equations in several variables, that is, equations of the form f 1t,,,x n d f
More informationAPPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS
LECTURE 10 APPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS Ordinary Differential Equations Initial Value Problems For Initial Value problems (IVP s), conditions are specified
More informationChap. 1. Some Differential Geometric Tools
Chap. 1. Some Differential Geometric Tools 1. Manifold, Diffeomorphism 1.1. The Implicit Function Theorem ϕ : U R n R n p (0 p < n), of class C k (k 1) x 0 U such that ϕ(x 0 ) = 0 rank Dϕ(x) = n p x U
More informationNUMERICAL METHOD FOR THE MIXED VOLTERRA-FREDHOLM INTEGRAL EQUATIONS USING HYBRID LEGENDRE FUNCTIONS
Conference Applications of Mathematics 215, in honor of the birthday anniversaries of Ivo Babuška (9), Milan Práger (85), and Emil Vitásek (85) J. Brandts, S. Korotov, M. Křížek, K. Segeth, J. Šístek,
More informationMathematical Methods - Lecture 9
Mathematical Methods - Lecture 9 Yuliya Tarabalka Inria Sophia-Antipolis Méditerranée, Titane team, http://www-sop.inria.fr/members/yuliya.tarabalka/ Tel.: +33 (0)4 92 38 77 09 email: yuliya.tarabalka@inria.fr
More informationOctober 25, 2013 INNER PRODUCT SPACES
October 25, 2013 INNER PRODUCT SPACES RODICA D. COSTIN Contents 1. Inner product 2 1.1. Inner product 2 1.2. Inner product spaces 4 2. Orthogonal bases 5 2.1. Existence of an orthogonal basis 7 2.2. Orthogonal
More informationLecture Notes for Math 251: ODE and PDE. Lecture 7: 2.4 Differences Between Linear and Nonlinear Equations
Lecture Notes for Math 51: ODE and PDE. Lecture 7:.4 Differences Between Linear and Nonlinear Equations Shawn D. Ryan Spring 01 1 Existence and Uniqueness Last Time: We developed 1st Order ODE models for
More informationConvergence and Error Bound Analysis for the Space-Time CESE Method
Convergence and Error Bound Analysis for the Space-Time CESE Method Daoqi Yang, 1 Shengtao Yu, Jennifer Zhao 3 1 Department of Mathematics Wayne State University Detroit, MI 480 Department of Mechanics
More informationResearch Article Numerical Solution of the Inverse Problem of Determining an Unknown Source Term in a Heat Equation
Applied Mathematics Volume 22, Article ID 39876, 9 pages doi:.55/22/39876 Research Article Numerical Solution of the Inverse Problem of Determining an Unknown Source Term in a Heat Equation Xiuming Li
More informationORDINARY DIFFERENTIAL EQUATIONS
ORDINARY DIFFERENTIAL EQUATIONS GABRIEL NAGY Mathematics Department, Michigan State University, East Lansing, MI, 48824. JANUARY 3, 25 Summary. This is an introduction to ordinary differential equations.
More informationREVIEW OF DIFFERENTIAL CALCULUS
REVIEW OF DIFFERENTIAL CALCULUS DONU ARAPURA 1. Limits and continuity To simplify the statements, we will often stick to two variables, but everything holds with any number of variables. Let f(x, y) be
More informationCHAPTER 10: Numerical Methods for DAEs
CHAPTER 10: Numerical Methods for DAEs Numerical approaches for the solution of DAEs divide roughly into two classes: 1. direct discretization 2. reformulation (index reduction) plus discretization Direct
More informationDeterministic Dynamic Programming
Deterministic Dynamic Programming 1 Value Function Consider the following optimal control problem in Mayer s form: V (t 0, x 0 ) = inf u U J(t 1, x(t 1 )) (1) subject to ẋ(t) = f(t, x(t), u(t)), x(t 0
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 7 Interpolation Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted
More informationLecture 10: Singular Perturbations and Averaging 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.243j (Fall 2003): DYNAMICS OF NONLINEAR SYSTEMS by A. Megretski Lecture 10: Singular Perturbations and
More informationSolutions to Homework 1, Introduction to Differential Equations, 3450: , Dr. Montero, Spring y(x) = ce 2x + e x
Solutions to Homewor 1, Introduction to Differential Equations, 3450:335-003, Dr. Montero, Spring 2009 problem 2. The problem says that the function yx = ce 2x + e x solves the ODE y + 2y = e x, and ass
More informationTransactions on Modelling and Simulation vol 12, 1996 WIT Press, ISSN X
Simplifying integration for logarithmic singularities R.N.L. Smith Department ofapplied Mathematics & OR, Cranfield University, RMCS, Shrivenham, Swindon, Wiltshire SN6 SLA, UK Introduction Any implementation
More informationSolving Nonlinear Two-Dimensional Volterra Integral Equations of the First-kind Using the Bivariate Shifted Legendre Functions
International Journal of Mathematical Modelling & Computations Vol. 5, No. 3, Summer 215, 219-23 Solving Nonlinear Two-Dimensional Volterra Integral Equations of the First-kind Using the Bivariate Shifted
More informationTHE INVERSE FUNCTION THEOREM
THE INVERSE FUNCTION THEOREM W. PATRICK HOOPER The implicit function theorem is the following result: Theorem 1. Let f be a C 1 function from a neighborhood of a point a R n into R n. Suppose A = Df(a)
More informationSection 2.4 Linear Equations
Section 2.4 Linear Equations Key Terms: Linear equation Homogeneous linear equation Nonhomogeneous (inhomogeneous) linear equation Integrating factor General solution Variation of parameters A first-order
More informationIntroduction to Initial Value Problems
Chapter 2 Introduction to Initial Value Problems The purpose of this chapter is to study the simplest numerical methods for approximating the solution to a first order initial value problem (IVP). Because
More informationStability investigation for the inverse heat conduction problems Krzysztof Grysa & Artur Maciqg Kielce University of Technology, AL 1000-leciaPP 7,
Stability investigation for the inverse heat conduction problems Krzysztof Grysa & Artur Maciqg Kielce University of Technology, AL 1000-leciaPP 7, Abstract Discretisation with respect to time changes
More informationUniformly Uniformly-ergodic Markov chains and BSDEs
Uniformly Uniformly-ergodic Markov chains and BSDEs Samuel N. Cohen Mathematical Institute, University of Oxford (Based on joint work with Ying Hu, Robert Elliott, Lukas Szpruch) Centre Henri Lebesgue,
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationL2 gains and system approximation quality 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 24: MODEL REDUCTION L2 gains and system approximation quality 1 This lecture discusses the utility
More informationA Gauss Lobatto quadrature method for solving optimal control problems
ANZIAM J. 47 (EMAC2005) pp.c101 C115, 2006 C101 A Gauss Lobatto quadrature method for solving optimal control problems P. Williams (Received 29 August 2005; revised 13 July 2006) Abstract This paper proposes
More informationON WEAK SOLUTION OF A HYPERBOLIC DIFFERENTIAL INCLUSION WITH NONMONOTONE DISCONTINUOUS NONLINEAR TERM
Internat. J. Math. & Math. Sci. Vol. 22, No. 3 (999 587 595 S 6-72 9922587-2 Electronic Publishing House ON WEAK SOLUTION OF A HYPERBOLIC DIFFERENTIAL INCLUSION WITH NONMONOTONE DISCONTINUOUS NONLINEAR
More informationVariables and Expressions
Variables and Expressions A variable is a letter that represents a value that can change. A constant is a value that does not change. A numerical expression contains only constants and operations. An algebraic
More informationBMIR Lecture Series on Probability and Statistics Fall 2015 Discrete RVs
Lecture #7 BMIR Lecture Series on Probability and Statistics Fall 2015 Department of Biomedical Engineering and Environmental Sciences National Tsing Hua University 7.1 Function of Single Variable Theorem
More informationFirst Order Differential Equations Lecture 3
First Order Differential Equations Lecture 3 Dibyajyoti Deb 3.1. Outline of Lecture Differences Between Linear and Nonlinear Equations Exact Equations and Integrating Factors 3.. Differences between Linear
More informationMath 3313: Differential Equations Laplace transforms
Math 3313: Differential Equations Laplace transforms Thomas W. Carr Department of Mathematics Southern Methodist University Dallas, TX Outline Introduction Inverse Laplace transform Solving ODEs with Laplace
More informationCONTROLLABILITY OF NONLINEAR SYSTEMS WITH DELAYS IN BOTH STATE AND CONTROL VARIABLES
KYBERNETIKA VOLUME 22 (1986), NUMBER 4 CONTROLLABILITY OF NONLINEAR SYSTEMS WITH DELAYS IN BOTH STATE AND CONTROL VARIABLES K. BALACHANDRAN Relative controllability of nonlinear systems with distributed
More informationAMATH 351 Mar 15, 2013 FINAL REVIEW. Instructor: Jiri Najemnik
AMATH 351 Mar 15, 013 FINAL REVIEW Instructor: Jiri Najemni ABOUT GRADES Scores I have so far will be posted on the website today sorted by the student number HW4 & Exam will be added early next wee Let
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 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 informationLaplace Transforms Chapter 3
Laplace Transforms Important analytical method for solving linear ordinary differential equations. - Application to nonlinear ODEs? Must linearize first. Laplace transforms play a key role in important
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 informationNon-linear wave equations. Hans Ringström. Department of Mathematics, KTH, Stockholm, Sweden
Non-linear wave equations Hans Ringström Department of Mathematics, KTH, 144 Stockholm, Sweden Contents Chapter 1. Introduction 5 Chapter 2. Local existence and uniqueness for ODE:s 9 1. Background material
More informationChapter 5 Types of Governing Equations. Chapter 5: Governing Equations
Chapter 5 Types of Governing Equations Types of Governing Equations (1) Physical Classification-1 Equilibrium problems: (1) They are problems in which a solution of a given PDE is desired in a closed domain
More informationThe Model and Preliminaries Problems and Solutions
Chapter The Model and Preliminaries Problems and Solutions Facts that are recalled in the problems The wave equation u = c u + G(x,t), { u(x,) = u (x), u (x,) = v (x), u(x,t) = f (x,t) x Γ, u(x,t) = x
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 informationSection 2.1 (First Order) Linear DEs; Method of Integrating Factors. General first order linear DEs Standard Form; y'(t) + p(t) y = g(t)
Section 2.1 (First Order) Linear DEs; Method of Integrating Factors Key Terms/Ideas: General first order linear DEs Standard Form; y'(t) + p(t) y = g(t) Integrating factor; a function μ(t) that transforms
More informationAbsorbing Boundary Conditions for Nonlinear Wave Equations
Absorbing Boundary Conditions for Nonlinear Wave Equations R. B. Gibson December, 2010 1 1 Linear Wave Equations 1.1 Boundary Condition 1 To begin, we wish to solve where with the boundary condition u
More informationYURI LEVIN, MIKHAIL NEDIAK, AND ADI BEN-ISRAEL
Journal of Comput. & Applied Mathematics 139(2001), 197 213 DIRECT APPROACH TO CALCULUS OF VARIATIONS VIA NEWTON-RAPHSON METHOD YURI LEVIN, MIKHAIL NEDIAK, AND ADI BEN-ISRAEL Abstract. Consider m functions
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationDiffusion Processes. Lectures INF2320 p. 1/72
Diffusion Processes Lectures INF2320 p. 1/72 Lectures INF2320 p. 2/72 Diffusion processes Examples of diffusion processes Heat conduction Heat moves from hot to cold places Diffusive (molecular) transport
More informationSimulation of an Inverse Heat Conduction Boundary Estimation Problem Based on State Space Model
Simulation of an Inverse Heat Conduction Boundary Estimation Problem Based on State Space Model Purna Chandra Mishra # # School of Mechanical Engineering, KIIT University, Patia, Bhubaneswar- 7, Odisha,
More informationSec. 1.1: Basics of Vectors
Sec. 1.1: Basics of Vectors Notation for Euclidean space R n : all points (x 1, x 2,..., x n ) in n-dimensional space. Examples: 1. R 1 : all points on the real number line. 2. R 2 : all points (x 1, x
More informationNon-linear Scalar Equations
Non-linear Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro August 24, 2014 1 / 44 Overview Here
More information4.7 Dispersion in an oscillatory shear flow
Lecture notes in Fluid ynamics.63j/.0j) by Chiang C. Mei, MIT, Spring, 007 4-6dispersion.tex March, 007 [Refs]:. Aris:. Fung, Y. C. Biomechanics 4.7 ispersion in an oscillatory shear flow Relevant to the
More information2. Time-Domain Analysis of Continuous- Time Signals and Systems
2. Time-Domain Analysis of Continuous- Time Signals and Systems 2.1. Continuous-Time Impulse Function (1.4.2) 2.2. Convolution Integral (2.2) 2.3. Continuous-Time Impulse Response (2.2) 2.4. Classification
More informationSymplectic integration with Runge-Kutta methods, AARMS summer school 2015
Symplectic integration with Runge-Kutta methods, AARMS summer school 2015 Elena Celledoni July 13, 2015 1 Hamiltonian systems and their properties We consider a Hamiltonian system in the form ẏ = J H(y)
More informationN.G.Bean, D.A.Green and P.G.Taylor. University of Adelaide. Adelaide. Abstract. process of an MMPP/M/1 queue is not a MAP unless the queue is a
WHEN IS A MAP POISSON N.G.Bean, D.A.Green and P.G.Taylor Department of Applied Mathematics University of Adelaide Adelaide 55 Abstract In a recent paper, Olivier and Walrand (994) claimed that the departure
More informationelastoplastic contact problems D. Martin and M.H. Aliabadi Wessex Institute of Technology, Ashurst Lodge, Ashurst, Southampton, SO40 7AA, UK
Non-conforming BEM elastoplastic contact problems D. Martin and M.H. Aliabadi discretisation in Wessex Institute of Technology, Ashurst Lodge, Ashurst, Southampton, SO40 7AA, UK Abstract In this paper,
More informationSpectrum and Exact Controllability of a Hybrid System of Elasticity.
Spectrum and Exact Controllability of a Hybrid System of Elasticity. D. Mercier, January 16, 28 Abstract We consider the exact controllability of a hybrid system consisting of an elastic beam, clamped
More informationDepartment of Structural, Faculty of Civil Engineering, Architecture and Urban Design, State University of Campinas, Brazil
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm A SIMPLIFIED FORMULATION FOR STRESS AND TRACTION BOUNDARY IN- TEGRAL EQUATIONS USING THE
More informationy + p(t)y = g(t) for each t I, and that also satisfies the initial condition y(t 0 ) = y 0 where y 0 is an arbitrary prescribed initial value.
p1 Differences Between Linear and Nonlinear Equation Theorem 1: If the function p and g are continuous on an open interval I : α < t < β containing the point t = t 0, then there exists a unique function
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationOld Math 330 Exams. David M. McClendon. Department of Mathematics Ferris State University
Old Math 330 Exams David M. McClendon Department of Mathematics Ferris State University Last updated to include exams from Fall 07 Contents Contents General information about these exams 3 Exams from Fall
More informationMA8109 Stochastic Processes in Systems Theory Autumn 2013
Norwegian University of Science and Technology Department of Mathematical Sciences MA819 Stochastic Processes in Systems Theory Autumn 213 1 MA819 Exam 23, problem 3b This is a linear equation of the form
More informationTransactions on Modelling and Simulation vol 8, 1994 WIT Press, ISSN X
A BEM model for convection- diffusion- migration problems in binary electrolytes Z.H. Qiu, L.C. Wrobel & H. Power Wessex Institute of Technology, University of Portsmouth, Ashurst Lodge, Ashurst, 7A.4,
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationSome recent results on controllability of coupled parabolic systems: Towards a Kalman condition
Some recent results on controllability of coupled parabolic systems: Towards a Kalman condition F. Ammar Khodja Clermont-Ferrand, June 2011 GOAL: 1 Show the important differences between scalar and non
More informationMaxima and Minima. (a, b) of R if
Maxima and Minima Definition Let R be any region on the xy-plane, a function f (x, y) attains its absolute or global, maximum value M on R at the point (a, b) of R if (i) f (x, y) M for all points (x,
More informationDifferential Games II. Marc Quincampoix Université de Bretagne Occidentale ( Brest-France) SADCO, London, September 2011
Differential Games II Marc Quincampoix Université de Bretagne Occidentale ( Brest-France) SADCO, London, September 2011 Contents 1. I Introduction: A Pursuit Game and Isaacs Theory 2. II Strategies 3.
More informationThe formulas for derivatives are particularly useful because they reduce ODEs to algebraic expressions. Consider the following ODE d 2 dx + p d
Solving ODEs using Fourier Transforms The formulas for derivatives are particularly useful because they reduce ODEs to algebraic expressions. Consider the following ODE d 2 dx + p d 2 dx + q f (x) R(x)
More informationFirst-Order Ordinary Differntial Equations: Classification and Linear Equations. David Levermore Department of Mathematics University of Maryland
First-Order Ordinary Differntial Equations: Classification and Linear Equations David Levermore Department of Mathematics University of Maryland 1 February 2009 These notes cover some of the material that
More informationThe method of lines (MOL) for the diffusion equation
Chapter 1 The method of lines (MOL) for the diffusion equation The method of lines refers to an approximation of one or more partial differential equations with ordinary differential equations in just
More informationChapter 2: First Order DE 2.4 Linear vs. Nonlinear DEs
Chapter 2: First Order DE 2.4 Linear vs. Nonlinear DEs First Order DE 2.4 Linear vs. Nonlinear DE We recall the general form of the First Oreder DEs (FODE): dy = f(t, y) (1) dt where f(t, y) is a function
More informationECON 582: An Introduction to the Theory of Optimal Control (Chapter 7, Acemoglu) Instructor: Dmytro Hryshko
ECON 582: An Introduction to the Theory of Optimal Control (Chapter 7, Acemoglu) Instructor: Dmytro Hryshko Continuous-time optimization involves maximization wrt to an innite dimensional object, an entire
More informationCHAPTER 5: Linear Multistep Methods
CHAPTER 5: Linear Multistep Methods Multistep: use information from many steps Higher order possible with fewer function evaluations than with RK. Convenient error estimates. Changing stepsize or order
More informationHOMEWORK 7 SOLUTIONS
HOMEWORK 7 SOLUTIONS MA11: ADVANCED CALCULUS, HILARY 17 (1) Using the method of Lagrange multipliers, find the largest and smallest values of the function f(x, y) xy on the ellipse x + y 1. Solution: The
More information10. Linear Systems of ODEs, Matrix multiplication, superposition principle (parts of sections )
c Dr. Igor Zelenko, Fall 2017 1 10. Linear Systems of ODEs, Matrix multiplication, superposition principle (parts of sections 7.2-7.4) 1. When each of the functions F 1, F 2,..., F n in right-hand side
More informationIntroduction to Algebraic and Geometric Topology Week 14
Introduction to Algebraic and Geometric Topology Week 14 Domingo Toledo University of Utah Fall 2016 Computations in coordinates I Recall smooth surface S = {f (x, y, z) =0} R 3, I rf 6= 0 on S, I Chart
More informationODEs Cathal Ormond 1
ODEs Cathal Ormond 2 1. Separable ODEs Contents 2. First Order ODEs 3. Linear ODEs 4. 5. 6. Chapter 1 Separable ODEs 1.1 Definition: An ODE An Ordinary Differential Equation (an ODE) is an equation whose
More informationReview for the Final Exam
Math 171 Review for the Final Exam 1 Find the limits (4 points each) (a) lim 4x 2 3; x x (b) lim ( x 2 x x 1 )x ; (c) lim( 1 1 ); x 1 ln x x 1 sin (x 2) (d) lim x 2 x 2 4 Solutions (a) The limit lim 4x
More informationNonlinear Optimization
Nonlinear Optimization (Com S 477/577 Notes) Yan-Bin Jia Nov 7, 2017 1 Introduction Given a single function f that depends on one or more independent variable, we want to find the values of those variables
More informationOrdinary Differential Equations
Chapter 7 Ordinary Differential Equations Differential equations are an extremely important tool in various science and engineering disciplines. Laws of nature are most often expressed as different equations.
More information