A MANIFOLD-BASED EXPONENTIALLY CONVERGENT ALGORITHM FOR SOLVING NON-LINEAR PARTIAL DIFFERENTIAL EQUATIONS
|
|
- Hugh Wilson
- 5 years ago
- Views:
Transcription
1 Journal of Marine Science and Technology Vol. No. pp. -9 () DOI:.69/JMST--- A MANIOLD-ASED EXPONENTIALLY CONVERGENT ALGORITHM OR SOLVING NON-LINEAR PARTIAL DIERENTIAL EQUATIONS Chein-Shan Liu Key words: non-linear algebraic equations (NAEs) non-linear partial differential equations manifold-based eponentially convergent algorithm (MECA) fictitious time integration method (TIM). ASTRACT or solving a non-linear system of algebraic equations of the type: i ( j ) i j n a Newton-like algorithm is still the most popular one; however it had some drawbacks as being locally convergent sensitive to initial guess and time consumption in finding the inversion of the Jacobian matri i / j. ased-on a manifold defined in the space of ( i t) we can derive a system of non-linear Ordinary Differential Equations (ODEs) in terms of the fictitious time-like variable t and the residual error is eponentially decreased to zero along the path of (t) by solving the resultant ODEs. We apply it to solve D non-linear PDEs and the vector-form of the matri-type non-linear algebraic equations (NAEs) is derived. Several numerical eamples of non-linear PDEs show the efficiency and accuracy of the present algorithm. A scalar equation is derived to find the adjustive fictitious time stepsize such that the irregular bursts appeared in the residual error curve can be overcome. We propose a future direction to construct a really eponentially convergent algorithm according to a manifold setting. I. INTRODUCTION Paper submitted /3/; accepted //. Author for correspondence: Chein-Shan Liu ( liucs@ntu.edu.tw). Department of Civil Engineering National Taiwan University Taipei Taiwan R.O.C. In many practical non-linear engineering problems governed by partial differential equations the methods such as the finite element method boundary element method finite volume method the meshless method etc. eventually led to a system of non-linear algebraic equations (NAEs). Many numerical methods used in the computational mechanics such as demonstrated by Zhu et al. [3] Atluri and Zhu [6] Atluri [] Atluri and Shen [5] and Atluri et al. [3] were always led to the solution of a system of linear algebraic equations for a linear problem and of an NAEs system for a non-linear problem. Over the past forty years two important contributions have been made towards the numerical solutions of NAEs. One of the methods has been called the predictor-corrector or pseudo-arclength continuation method. This method has its historical roots in the embedding and incremental loading methods which have been successfully used for several decades by engineers to improve the convergence properties when an adequate starting value for an iterative method is not available. Another is the so-called simplical or piecewise linear method. The monographs by Allgower and Georg [] and Deuflhard [] are devoted to the continuation methods for solving NAEs. Liu and Atluri [3] have employed the technique of ictitious Time Integration Method (TIM) to solve a large system of NAEs and showed that high performance can be achieved by using the TIM. More recently Liu [7] has used the TIM technique to solve the non-linear complementarity problems. Then Liu [8 9] has used the TIM to solve the boundary value problems of elliptic type partial differential equations. Liu and Atluri [] also employed this technique of TIM to solve mied-complementarity problems and optimization problems. Liu and Atluri [5] using the technique of TIM solved the inverse Sturm-Liouville problem under specified eigenvalues. Liu [] has utilized a simple TIM to compute both the forward and backward in time problems of urgers equation and has found that the TIM is robust against the disturbance of noise. Liu and Chang [8] have used the TIM to solve ill-posed linear problems. or its numerical implementation being quite simple the TIM was also used in other problems like as Liu [ ] Liu and Atluri [6] Chang and Liu [8] Chi et al. [] and Tsai et al. [3]. In spite of its success the TIM is local convergence and needs to determine viscous damping coefficients for different equations in one problem. Liu and Atluri [7] were the first to
2 Journal of Marine Science and Technology Vol. No. () make a breakthrough by deriving a residual-norm based algorithm which can circumvent the drawbacks of TIM and Newton s algorithm; also refer Atluri et al. []. In this paper we introduce a novel continuation method of Manifold-ased Eponentially Convergent Algorithm (MECA) which can be easily implemented to solve nonlinear partial differential equations after some suitable discretizations by using the finite difference method or epansion by the radial basis function []. To remedy the shortcoming of vector homotopy method as initiated by Davidenko [] Liu et al. [9] have proposed a scalar homotopy method and Ku et al. [5] combined this idea with the eponentially decayed scalar homotopy function developing a Manifold- ased Eponentially Convergent Algorithm (MECA). In this paper we will point out the limitation of MECA and eplain why it may stuck and after that the algorithm cannot proceed to find solution. II. THEORETICAL ASIS or the following non-linear algebraic equations (NAEs) in a vector form: ( ) () we will propose a numerical algorithm to solve it based on a space-time manifold.. Motivation To motivate the present approach let us consider an uncoupled algebraic system: y. () The above two equations can be combined into a single one: [( ) ( ) ]. + y (3) We can see that it is a circle in the spatial-plane ( y) with a center () but with a zero radius. Now we introduce an etra variable t as a fictitious timelike variable and let and y be functions of t. We consider the following equation defined in the space-time domain by αt αt e e [( ) + ( y ) ] C () where α > and C is determined by the initial values of () and y() y by C [( ) + ( y ) ]. Eq. () can be written as ( ) ( y ) Ce αt + (5) which is a manifold in the space-time domain ( y t) at each cross-section of a fied t it being a circle with center () and with a radius Ce αt. We can construct the following ODEs: α( ) y α( y ) (6) such that the path of ((t) y(t)) generated from the above ODEs is located on the manifold defined by Eq. (5). Solving Eq. (6) we have αt t () [ ] e + yt () [ y ] e +. (7) αt Inserting them into Eq. () we indeed can prove that ((t) y(t)) is located on the manifold and ((t) y(t)) fast tends to the solution ().. A Setting ased on Manifold rom the above idea of space-time manifold for Eq. () we can consider h( t): Q( t) ( ) C (8) where Q(t) > is a given function of t monotonically increasing with t and with Q() and C is determined by the initial condition () with C ( ). (9) Usually C > and C when the initial value is just the root of Eq. (). However it is rare of this lucky case. When C > and Q > the manifold defined by Eq. (8) is continuous and thus the differential operation being taken on the manifold makes sense. or the requirement of consistency condition by taking the differential of Eq. (8) with respect to t and considering (t) we have () ( ) ()( ) Qt Qt + () where is the Jacobian matri with its ij-component given by ij i / j. Eq. () cannot uniquely determine the governing equation of ; however we suppose that h λ λqt () ()
3 C.-S. Liu: A Manifold-ased Algorithm 3 where λ is to be determined. Inserting Eq. () into Eq. () we can solve for λ by λ Qt () Q ( t). () Thus we obtain an evolution equation for defined by the following ODEs: III. NUMERICAL METHODS. Adjusting the ictitious Time Step In order to keep on the manifold defined by Eq. (8) we can consider the evolution of along the path (t) by qt () A (9) qt () (3) where A:. () where Qt () qt ():. Qt ( ) () There are many ways to choose a suitable function of Q(t); however we can let Hence we have Qt () ν qt () < m. Qt ( ) ( + t) m ν Qt + m m () ep [( t) ]. Therefore we can derive the following equation:. (5) (6) ν (7) ( + t) m It is not difficult to prove that the orbit (t) with () solved from Eq. (7) is located on the space-time manifold defined by Eq. (8) with Q(t) given by Eq. (6). Hence we have an eponentially convergent property in solving the NAEs in Eq. (): Suppose that we simply use the Euler forward scheme to integrate Eq. (9) which yields ( t+ t) () t q() t t A. () Taking the square-norms of both the sides and using Eq. (8) we can obtain C C C ( A) qt ( ) t Qt ( + t) Qt () Qt () C ( qt ( ) t) A. + Qt () Thus we have the following scalar equation to solve t: where Qt () f( t): a( t) b t + Qt ( + t) A : ( ) a q t () (3) () b: q( t). (5) C ( ). (8) Qt () When t is increasing the above equation enforces the residual error ( ) tending to zero eponentially and meanwhile the solution of Eq. () is obtained approimately. However it is a great challenge by developing a suitable numerical integrator for Eq. (7) such that the orbit of can really retain on the manifold. Obviously t is a root of Eq. (3); however we prefer to search another one with t >. On the other hand we note that ν in Eq. (6) cannot be too large say ν > ; otherwise the term Q(t)/Q(t + t) in Eq. (3) can be very near to zero and thus makes Eq. (3) having no real solution because of ( b) a q A <
4 Journal of Marine Science and Technology Vol. No. () which is due to ( A) < A. (6) We can apply the TIM to find the solution of Eq. (3): µ f( ). (7) ( + t) γ variables K and the vectorial algebraic equations K by Do i n i j i j Do j n K n ( i ) + j K K u. (3) When t is solved we can use the following iteration to calculate (t + t): ( t+ t) () t q() t t. (8) The above algorithm is indeed a very powerful numerical method with eponentially convergent speed and the orbit of (t) being retained on the manifold. Thus we may call the present algorithm a Manifold-ased Eponentially- Convergent Algorithm (MECA).. A Matri Type NAEs Sometimes we may encounter the NAEs generated from the discretization of PDE with a matri type. In this situation it is not so straightforward to write its counterpart as being a vector-form NAEs. Let us consider uy ( ) yuu ( u )( y ) Ω (9) y uy ( ) Hy ( )( y ) Γ (3) where is the Laplacian operator Γ is the boundary of a problem domain Ω : [a a ] [b b ] and and H are given functions. y a standard finite difference applied to Eq. (9) one has a system of NAEs of matri type: i j [ u i+ j ui j+ ui j] + [ u i j+ ui j+ ui j ] ( ) ( y) u u u u y i n j n. i+ j i j i j+ i j i yj ui j (3) Here we divide the rectangle Ω by a uniform grid with (a a )/(n + ) and y (b b )/(n + ) being the uniform spatial grid lengths in the - and y-direction and u i j : u( i y j ) be a numerical value of u at the grid point ( i y j ) Ω. Let K n (i ) + j and with i running from to n and j running from to n we can respectively set the vectorial At the same time the components of the Jacobian matri are constructed by Do i n Do j n K n ( i ) + j L n ( i ) + j L n ( i ) + j L n ( i ) + j 3 L n ( i ) + j+ L5 ni+ j ( ) ui+ j ui j ui j+ ui j + u i yj ui j y KL ( y) ui+ j ui j ui j+ ui j + u y i yj ui j y y KL ( ) ( y) KL 3 u u u u y i+ j i j i j+ i j u i yj ui j ( y) ui+ j ui j ui j+ ui j u y i yj ui j y y KL ( ) KL 5 ui+ j ui j ui j+ ui j u. i yj ui j y (33)
5 C.-S. Liu: A Manifold-ased Algorithm 5 In above u u and uy denote respectively the partial differentials of the function ( y u u u y ) with respect to u u and u y. When the above quantities are available we can apply the vector-form MECA to solve the non-linear PDE in Eq. (9). IV. NUMERICAL TESTS O MECA In this section we apply the new method of MECA to solve some non-linear PDEs.. Eample We consider a non-linear heat conduction equation: 8 6 α u α u + α u + u + h t (3) t ( ) '( ) ( ) ( ) ( 3) ( ) 7( 3) t t h t e ( 3) e (35).7.6 with a closed-form solution u( t) ( 3) e t. y applying the MECA to solve the above equation in the domain of and t we fi n n 3 ν 5 and m. We apply the group-preserving scheme (GPS) developed by Liu [6] to integrate the resultant ODEs with a fictitious time stepsize t.5. In ig. we show the residual errors with respect to the number of steps up to. The absolute errors of numerical solution are plotted in ig. which reveal an accurate numerical result with the maimum error being Eample One famous mesh-less numerical method to solve the non-linear PDE of elliptic type is the radial basis function (R) method which epands the solution u by Error of u ig.. or eample : showing the convergence speed of residual errors and displaying the numerical error. t n uy ( ) aφ (36) k where a k are the epansion coefficients to be determined and φ k is a set of Rs for eample N 3/ φk ( rk + c ) N N φk rk ln rk N r k φk ep a r φk + a N 3/ k ( rk c ) ep N k k (37) where the radius function r k is given by r k ( k) + ( y yk) while ( k y k ) k n are called source points. The constants a and c are shape parameters. In the below we take the first set of φ k as trial functions which is known as a multi-quadric R [9 3] with N. In this eample we apply the multi-quadric radial basis function to solve the following non-linear PDE: 3 u u ( + y + a ) (38) where a was fied. The domain is an irregular domain with ρ( θ) (sin θ) ep(sin θ) + (cos θ) ep(cos θ). (39) The analytic solution is given by uy ( ) + y a which is singular on the circle with a radius a. ()
6 6 Journal of Marine Science and Technology Vol. No. (). E Residual error and stepsize E- E- Residual error by solving Δt Δt Residual error with a fied stepsize Error of u ig.. or eample : showing the convergence speed of residual errors and displaying the numerical error. Inserting Eq. (36) into Eq. (38) and placing some field points inside the domain to satisfy the governing equation and some points on the boundary to satisfy the boundary condition we can derive n NAEs to determine the n coefficients a k. The source points ( k y k ) k n are uniformly distributed on a contour given by R + ρ(θ k ) where θ k kπ/n. Under the following parameters R.5 c.5 ν 5 m. and t. in ig. we show the residual errors up to steps. It can be seen that in the residual-error curve there is a little oscillation and it decays very fast at the first few steps. The absolute error of numerical solution is plotted in ig. which is accurate with the maimum error being Now we test the effect by adjusting the stepsize. Under the following parameters R.5 c. ν 5 and m. in ig. 3 we compare the residual errors obtained by fiing a stepsize with t. which ehibits many irregular bursts as shown by the dashed line and that obtained by adjusting the stepsize as shown by a solid line which is smooth and decreased. The adaptive stepsize is also shown by the dasheddotted line up to 5 steps. - - y E Number of steps ig. 3. or eample comparing the residual errors obtained by a fied stepsize and an adaptive stepsize. 3. Eample 3 In this eample we apply the MECA to solve the following boundary value problem of non-linear elliptic equation [ ]: 3 uy ( ) + ω uy ( ) + εu( y ) py ( ). () While the eact solution is uy ( ) ( + y) + 3( y+ y) () 6 the eact p can be obtained by inserting the above u into Eq. (). y introducing a finite difference discretization of u at the grid points we can obtain i j ( u i+ j ui j+ ui j) + ( u i j+ ui j+ ui j ) ( ) ( y) + ω u + εu p. 3 i j i j i j (3) The boundary conditions can be obtained from the eact solution in Eq. (). Under the following parameters n n 3 t.5 ν 5 ω and ε. we compute the solution of the above system and compare them with the eact solution. In ig. we show the residual errors up to steps. The absolute errors of numerical solution are plotted in ig. of which the maimum error is smaller than.66 3.
7 C.-S. Liu: A Manifold-ased Algorithm 7.3 Stepsize Error of u. 6 8 ig. 5. or eample : displaying time stepsize and showing the convergence speed of residual errors. However the stepsize adjusting technique is computational time consuming ig.. or eample 3: showing the convergence speed of residual errors and displaying the numerical error..8 V. LIMITATION efore Eq. () we have mentioned that the governing equation of is not unique. Indeed it is not a necessary condition but a sufficient condition to satisfy Eq. (). Thus the following algorithm:.. ( t+ t) () t q() t t.6 y.8 () based-on Eq. (3) has some limitations. Obviously from the residual error curves as shown in igs. 3 and there is a common pattern that the curves decrease very fast at the beginning and then they all tend to be flattened with a very slow convergence without having an eponential convergence no longer. Another feature as obviously shown in igs. and is that the curves ehibit irregular bursts. To solve the latter problem we can employ the technology in Section III. by adjusting the time stepsize. As shown in ig. 5 we plot the adjustive time stepsize and residual error for eample. In the curve of residual error the bursts disappear. However it still has the first problem. More difficultly the time stepsize adjusting technique spends a lot of computational time to find t. In order to find the reason for a slow convergence of the residual error curve in its later stage we give the following analysis. Inserting Eqs. () and (5) into Eq. (3) we can derive where Qt () Qt ( + t) ( ) ( ) + a q t q t A a : in view of Eq. (6). rom Eq. () and the approimation of we have Q () t t Q( t+ t) Q() t (5) (6) qt () t [ Rt () ] (7)
8 8 Journal of Marine Science and Technology Vol. No. () a R (c) When R(t) tends to i.e. a tends to the dynamical force on the right-hand side is lost for the above equation. In ig. 6 we show a R and the residual error with respect to the number of steps. It can be seen that the above dynamic equation loses its pushing force and remains a large residual error due to a tending to and thus R tending to. VI. CONCLUSIONS A manifold-based eponentially convergent algorithm (MECA) was established in this paper. When we apply it to solve some D non-linear PDEs the vector-form was derived. Several numerical eamples of non-linear PDEs showed the efficiency and accuracy of MECA. The irregular bursts appeared in the residual error curve can be overcome by solving a scalar equation to find the adjustable fictitious time stepsize. However the governing Eq. (3) has faced a bottle-neck by only considering with the gradient vector as a unique source of dynamic force. An open problem is that how to construct a really eponentially convergent algorithm according to the setting of a manifold defined in Eq. (8). 6 8 ig. 6. or eample : displaying a showing R and (c) the convergence speed of residual errors. where the ratio Rt () is defined by Qt ( + t) Rt (). (8) Qt () As a requirement of Qt () > we need R(t) >. Thus through some manipulations Eq. (5) becomes ar() t (a + ) R() t + ( a + 8) Rt (). (9) 3 It is interesting that the above equation can be written as [ Rt ( ) ] [ art ( ) ]. (5) ecause R is a double roots and it is not the desired one we take Rt (). a A y using Eq. (7) Eq. () can now be written as ( t+ t) () t [ R() t ]. (5) (5) ACKNOWLEDGMENTS Taiwan s National Science Council project NSC-99-- E--7-MY3 granted to the author is highly appreciated. REERENCES. Allgower E. L. and Georg K. Numerical Continuation Methods: an Introduction Springer New York (99).. Atluri S. N. Methods of Computer Modeling in Engineering and Sciences Tech Science Press (). 3. Atluri S. N. Liu H. T. and Han Z. D. Meshless local Petrov-Galerkin (MLPG) mied collocation method for elasticity problems CMES: Computer Modeling in Engineering & Sciences Vol. pp. -5 (6).. Atluri S. N. Liu C.-S. and Kuo C. L. A modified Newton method for solving non-linear algebraic equations Journal of Marine Science and Technology Vol. 7 pp (9). 5. Atluri S. N. and Shen S. The meshless local Petrov-Galerkin (MLPG) method: a simple & less-costly alternative to the finite and boundary element methods CMES: Computer Modeling in Engineering & Sciences Vol. 3 pp. -5 (). 6. Atluri S. N. and Zhu T. L. A new meshless local Petrov-Galerkin (MLPG) approach in computational mechanics Computational Mechanics Vol. pp. 7-7 (998). 7. Atluri S. N. and Zhu T. L. A new meshless local Petrov-Galerkin (MLPG) approach to nonlinear problems in computer modeling and simulation Computational Modeling and Simulation in Engineering Vol. 3 pp (998). 8. Chang C. W. and Liu C.-S. A fictitious time integration method for backward advection-dispersion equation CMES: Computer Modeling in Engineering & Sciences Vol. 5 pp (9). 9. Cheng A. H. D. Golberg M. A. Kansa E. J. and Zammito G. Eponential convergence and H-c multiquadric collocation method for partial differential equations Numerical Methods for Partial Differential Equations Vol. 9 pp (3).. Chi C. C. Yeih W. and Liu C.-S. A novel method for solving the Cauchy problem of Laplace equation using the fictitious time integration
9 C.-S. Liu: A Manifold-ased Algorithm 9 method CMES: Computer Modeling in Engineering & Sciences Vol. 7 pp (9).. Davidenko D. On a new method of numerically integrating a system of nonlinear equations Doklady Akademii Nauk SSSR Vol. 88 pp. 6-6 (953).. Deuflhard P. Newton Methods for Nonlinear Problems: Affine Invariance and Adaptive Algorithms Springer New York (). 3. Golberg M. A. Chen C. S. and Karur S. R. Improved multiquadric approimation for partial differential equations Engineering Analysis with oundary Elements Vol. 8 pp. 9-7 (996).. Hu H. Y. Li Z. C. and Cheng A. H. D. Radial basis collocation methods for elliptic boundary value problems Computers and Mathematics with Applications Vol. 5 pp (5). 5. Ku C. Y. Yeih W. and Liu C.-S. Solving non-linear algebraic equations by a scalar Newton-homotopy continuation method International Journal of Nonlinear Sciences & Numerical Simulation Vol. pp (). 6. Liu C.-S. Cone of non-linear dynamical system and group preserving schemes International Journal of Non-Linear Mechanics Vol. 36 pp (). 7. Liu C.-S. A time-marching algorithm for solving non-linear obstacle problems with the aid of an NCP-function CMC: Computers Materials & Continua Vol. 8 pp (8). 8. Liu C.-S. A fictitious time integration method for two-dimensional quasilinear elliptic boundary value problems CMES: Computer Modeling in Engineering & Sciences Vol. 33 pp (8). 9. Liu C.-S. A fictitious time integration method for a quasilinear elliptic boundary value problem defined in an arbitrary plane domain CMC: Computers Materials & Continua Vol. pp. 5-3 (9).. Liu C.-S. A fictitious time integration method for the urgers equation CMC: Computers Materials & Continua Vol. 9 pp. 9-5 (9).. Liu C.-S. A fictitious time integration method for solving delay ordinary differential equations CMC: Computers Materials & Continua Vol. pp (9).. Liu C.-S. A fictitious time integration method for solving m-point boundary value problems CMES: Computer Modeling in Engineering & Sciences Vol. 39 pp. 5-5 (9). 3. Liu C.-S. and Atluri S. N. A novel time integration method for solving a large system of non-linear algebraic equations CMES: Computer Modeling in Engineering & Sciences Vol. 3 pp (8).. Liu C.-S. and Atluri S. N. A fictitious time integration method (TIM) for solving mied complementarity problems with applications to non-linear optimization CMES: Computer Modeling in Engineering & Sciences Vol. 3 pp (8). 5. Liu C.-S. and Atluri S. N. A novel fictitious time integration method for solving the discretized inverse Sturm-Liouville problems for specified eigenvalues CMES: Computer Modeling in Engineering & Sciences Vol. 36 pp (8). 6. Liu C.-S. and Atluri S. N. A fictitious time integration method for the numerical solution of the redholm integral equation and for numerical differentiation of noisy data and its relation to the filter theory CMES: Computer Modeling in Engineering & Sciences Vol. pp. 3-6 (9). 7. Liu C.-S. and Atluri S. N. Simple residual-norm based algorithms for the solution of a large system of non-linear algebraic equations which converge faster than the Newton s method CMES: Computer Modeling in Engineering & Sciences Vol. 7 pp (). 8. Liu C.-S. and Chang C. W. Novel methods for solving severely ill-posed linear equations system Journal of Marine Science and Technology Vol. 7 pp. 6-7 (9). 9. Liu C.-S. Yeih W. Kuo C. L. and Atluri S. N. A scalar homotopy method for solving an over/under-determined system of non-linear algebraic equations CMES: Computer Modeling in Engineering & Sciences Vol. 53 pp. 7-7 (9). 3. Tsai C. C. Liu C.-S. and Yeih W. ictitious time integration method of fundamental solutions with Chebyshev polynomials for solving Poisson-type nonlinear PDEs CMES: Computer Modeling in Engineering & Sciences Vol. 56 pp. 3-5 (). 3. Zhu T. Zhang J. and Atluri S. N. A meshless local boundary integral equation (LIE) method for solving nonlinear problems Computational Mechanics Vol. pp (998). 3. Zhu T. Zhang J. and Atluri S. N. A meshless numerical method based on the local boundary integral equation (LIE) to solve linear and nonlinear boundary value problems Engineering Analysis with oundary Elements Vol. 3 pp (999).
A RESIDUAL-NORM BASED ALGORITHM FOR SOLVING DETERMINATE/INDETERMINATE SYSTEMS OF NON-LINEAR ALGEBRAIC EQUATIONS
Journal of Marine cience and echnology, Vol., No., pp. 8-97 () 8 DOI:.69/JM--9- A REIDUAL-NORM BAED ALGORIHM FOR OLVING DEERMINAE/INDEERMINAE EM OF NON-LINEAR ALGEBRAIC EQUAION ung-wei Chen Key words:
More information(P = R F(F R)/ F 2 ) in Iteratively Solving the Nonlinear System of Algebraic Equations F(x) = 0
Copyright 2011 Tech Science Press CMES, vol.81, no.2, pp.195-227, 2011 A Further Study on Using ẋ = λ[αr + βp] (P = F R(F R)/ R 2 ) and ẋ = λ[αf + βp ] (P = R F(F R)/ F 2 ) in Iteratively Solving the Nonlinear
More informationA Novel Time Integration Method for Solving A Large System of Non-Linear Algebraic Equations
Copyright c 2008 Tech Science Press CMES, vol.31, no.2, pp.71-83, 2008 A Novel Time Integration Method for Solving A Large System of Non-Linear Algebraic Equations Chein-Shan Liu 1 and Satya N. Atluri
More informationIterative Solution of a System of Nonlinear Algebraic
Copyright 2011 Tech Science Press CMES, vol.81, no.4, pp.335-362, 2011 Iterative Solution of a System of Nonlinear Algebraic Equations F(x) = 0, Using ẋ = λ[αr + βp] or ẋ = λ[αf + βp ] R is a Normal to
More informationSolution of Post-Buckling & Limit Load Problems, Without Inverting the Tangent Stiffness Matrix & Without Using Arc-Length Methods
Copyright 214 Tech Science Press CMES, vol.98, no.6, pp.543-563, 214 Solution of Post-Buckling & Limit Load Problems, Without Inverting the Tangent Stiffness Matrix & Without Using Arc-Length Methods T.A.
More informationCopyright 2009 Tech Science Press CMES, vol.41, no.3, pp , 2009
Copyright 29 Tech Science Press CMES, vol.41, no.3, pp.243-261, 29 A Fictitious Time Integration Method for the Numerical Solution of the Fredholm Integral Equation and for Numerical Differentiation of
More informationA Globally Optimal Iterative Algorithm to Solve an Ill-Posed Linear System
Copyright 2012 Tech Science Press CMES, vol.84, no.4, pp.383-403, 2012 A Globally Optimal Iterative Algorithm to Solve an Ill-Posed Linear System Chein-Shan Liu 1 Abstract: An iterative algorithm based
More informationLecture 4.2 Finite Difference Approximation
Lecture 4. Finite Difference Approimation 1 Discretization As stated in Lecture 1.0, there are three steps in numerically solving the differential equations. They are: 1. Discretization of the domain by
More informationA Two-Side Equilibration Method to Reduce the Condition Number of an Ill-Posed Linear System
Copyright 2013 Tech Science Press CMES, vol.91, no.1, pp.17-42, 2013 A Two-Side Equilibration Method to Reduce the Condition Number of an Ill-Posed Linear System Chein-Shan Liu 1 Abstract: In the present
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 4, pp. 373-38, 23. Copyright 23,. ISSN 68-963. ETNA A NOTE ON THE RELATION BETWEEN THE NEWTON HOMOTOPY METHOD AND THE DAMPED NEWTON METHOD XUPING ZHANG
More informationThe Jordan Structure of Residual Dynamics Used to Solve Linear Inverse Problems
Copyright 2012 Tech Science Press CMES, vol.88, no.1, pp.29-47, 2012 The Jordan Structure of Residual Dynamics Used to Solve Linear Inverse Problems Chein-Shan Liu 1, Su-Ying Zhang 2 and Satya N. Atluri
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 informationModifications of Steepest Descent Method and Conjugate Gradient Method Against Noise for Ill-posed Linear Systems
Available online at www.ispacs.com/cna Volume 2012, Year 2012 Article ID cna-00115, 24 pages doi: 10.5899/2012/cna-00115 Research Article Modifications of Steepest Descent Method and Conjugate Gradient
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 informationMeshless Local Petrov-Galerkin Mixed Collocation Method for Solving Cauchy Inverse Problems of Steady-State Heat Transfer
Copyright 2014 Tech Science Press CMES, vol.97, no.6, pp.509-533, 2014 Meshless Local Petrov-Galerkin Mixed Collocation Method for Solving Cauchy Inverse Problems of Steady-State Heat Transfer Tao Zhang
More informationNEW INVESTIGATIONS INTO THE MFS FOR 2D INVERSE LAPLACE PROBLEMS. Received June 2015; revised October 2015
International Journal of Innovative Computing, Information and Control ICIC International c 2016 ISSN 1349-4198 Volume 12, Number 1, February 2016 pp. 55 66 NEW INVESTIGATIONS INTO THE MFS FOR 2D INVERSE
More informationCopyright 2010 Tech Science Press CMES, vol.60, no.3, pp , 2010
Copyright 2010 Tech Science Press CMES, vol.60, no.3, pp.279-308, 2010 Novel Algorithms Based on the Conjugate Gradient Method for Inverting Ill-Conditioned Matrices, and a New Regularization Method to
More informationCopyright 2013 Tech Science Press CMES, vol.94, no.1, pp.1-28, 2013
Copyright 2013 Tech Science Press CMES, vol.94, no.1, pp.1-28, 2013 Application of the MLPG Mixed Collocation Method for Solving Inverse Problems of Linear Isotropic/Anisotropic Elasticity with Simply/Multiply-Connected
More informationSolving Elastic Problems with Local Boundary Integral Equations (LBIE) and Radial Basis Functions (RBF) Cells
Copyright 2010 Tech Science Press CMES, vol.57, no.2, pp.109-135, 2010 Solving Elastic Problems with Local Boundary Integral Equations (LBIE) and Radial Basis Functions (RBF) Cells E. J. Sellountos 1,
More informationStatistical Geometry Processing Winter Semester 2011/2012
Statistical Geometry Processing Winter Semester 2011/2012 Linear Algebra, Function Spaces & Inverse Problems Vector and Function Spaces 3 Vectors vectors are arrows in space classically: 2 or 3 dim. Euclidian
More informationChapter 6. Nonlinear Equations. 6.1 The Problem of Nonlinear Root-finding. 6.2 Rate of Convergence
Chapter 6 Nonlinear Equations 6. The Problem of Nonlinear Root-finding In this module we consider the problem of using numerical techniques to find the roots of nonlinear equations, f () =. Initially we
More informationExact and Approximate Numbers:
Eact and Approimate Numbers: The numbers that arise in technical applications are better described as eact numbers because there is not the sort of uncertainty in their values that was described above.
More informationarxiv:gr-qc/ v1 6 Sep 2006
Introduction to spectral methods Philippe Grandclément Laboratoire Univers et ses Théories, Observatoire de Paris, 5 place J. Janssen, 995 Meudon Cede, France This proceeding is intended to be a first
More informationA Harmonic Balance Approach for Large-Scale Problems in Nonlinear Structural Dynamics
A Harmonic Balance Approach for Large-Scale Problems in Nonlinear Structural Dynamics Allen R, PhD Candidate Peter J Attar, Assistant Professor University of Oklahoma Aerospace and Mechanical Engineering
More informationNumerical Integration (Quadrature) Another application for our interpolation tools!
Numerical Integration (Quadrature) Another application for our interpolation tools! Integration: Area under a curve Curve = data or function Integrating data Finite number of data points spacing specified
More informationEfficient solution of stationary Euler flows with critical points and shocks
Efficient solution of stationary Euler flows with critical points and shocks Hans De Sterck Department of Applied Mathematics University of Waterloo 1. Introduction consider stationary solutions of hyperbolic
More informationThe Method of Fundamental Solutions for One-Dimensional Wave Equations
Copyright 9 Tech Science Press CMC, vol., no.3, pp.85-8, 9 The Method of Fundamental Solutions for One-Dimensional Wave Equations Gu, M. H., Young, D. L., and Fan, C. M. Abstract: A meshless numerical
More informationNational Taiwan University
National Taiwan University Meshless Methods for Scientific Computing (Advisor: C.S. Chen, D.L. oung) Final Project Department: Mechanical Engineering Student: Kai-Nung Cheng SID: D9956 Date: Jan. 8 The
More informationNumerical solution of Maxwell equations using local weak form meshless techniques
Journal of mathematics and computer science 13 2014), 168-185 Numerical solution of Maxwell equations using local weak form meshless techniques S. Sarabadan 1, M. Shahrezaee 1, J.A. Rad 2, K. Parand 2,*
More informationLong Time Dynamics of Forced Oscillations of the Korteweg-de Vries Equation Using Homotopy Perturbation Method
Studies in Nonlinear Sciences 1 (3): 57-6, 1 ISSN 1-391 IDOSI Publications, 1 Long Time Dynamics of Forced Oscillations of the Korteweg-de Vries Equation Using Homotopy Perturbation Method 1 Rahmat Ali
More informationA high order adaptive finite element method for solving nonlinear hyperbolic conservation laws
A high order adaptive finite element method for solving nonlinear hyperbolic conservation laws Zhengfu Xu, Jinchao Xu and Chi-Wang Shu 0th April 010 Abstract In this note, we apply the h-adaptive streamline
More informationRegularized Meshless Method for Solving Acoustic Eigenproblem with Multiply-Connected Domain
Copyright c 2006 Tech Science Press CMES, vol16, no1, pp27-39, 2006 Regularized Meshless Method for Solving Acoustic Eigenproblem with Multiply-Connected Domain KH Chen 1,JTChen 2 and JH Kao 3 Abstract:
More informationTheory, Solution Techniques and Applications of Singular Boundary Value Problems
Theory, Solution Techniques and Applications of Singular Boundary Value Problems W. Auzinger O. Koch E. Weinmüller Vienna University of Technology, Austria Problem Class z (t) = M(t) z(t) + f(t, z(t)),
More information3.1: 1, 3, 5, 9, 10, 12, 14, 18
3.:, 3, 5, 9,,, 4, 8 ) We want to solve d d c() d = f() with c() = c = constant and f() = for different boundary conditions to get w() and u(). dw d = dw d d = ( )d w() w() = w() = w() ( ) c d d = u()
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationSOLVING INHOMOGENEOUS PROBLEMS BY SINGULAR BOUNDARY METHOD
8 Journal of Marine Science and Technology Vol. o. pp. 8-4 03) DOI: 0.69/JMST-0-0704- SOLVIG IHOMOGEEOUS PROBLEMS BY SIGULAR BOUDARY METHOD Xing Wei Wen Chen and Zhuo-Jia Fu Key words: singular boundary
More informationComputation Time Assessment of a Galerkin Finite Volume Method (GFVM) for Solving Time Solid Mechanics Problems under Dynamic Loads
Proceedings of the International Conference on Civil, Structural and Transportation Engineering Ottawa, Ontario, Canada, May 4 5, 215 Paper o. 31 Computation Time Assessment of a Galerkin Finite Volume
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationA numerical study of a technique for shifting eigenvalues of radial basis function differentiation matrices
A numerical study of a technique for shifting eigenvalues of radial basis function differentiation matrices Scott A. Sarra Marshall University and Alfa R.H. Heryudono University of Massachusetts Dartmouth
More informationElastic Buckling Behavior of Beams ELASTIC BUCKLING OF BEAMS
Elastic Buckling Behavior of Beams CE579 - Structural Stability and Design ELASTIC BUCKLING OF BEAMS Going back to the original three second-order differential equations: 3 Therefore, z z E I v P v P 0
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. II - Solutions to the Cauchy problem. Alberto Bressan
Hyperbolic Systems of Conservation Laws in One Space Dimension II - Solutions to the Cauchy problem Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 Global
More informationPerformance of Multiquadric Collocation Method in Solving Lid-driven Cavity Flow Problem with Low Reynolds Number
Copyright c 2006 Tech Science Press CMES, vol.15, no.3, pp.137-146, 2006 Performance of Multiquadric Collocation Method in Solving Lid-driven Cavity Flow Problem with Low Reynolds umber S. Chantasiriwan
More informationNonconstant Coefficients
Chapter 7 Nonconstant Coefficients We return to second-order linear ODEs, but with nonconstant coefficients. That is, we consider (7.1) y + p(t)y + q(t)y = 0, with not both p(t) and q(t) constant. The
More informationIndex. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2
Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604
More informationA truly meshless Galerkin method based on a moving least squares quadrature
A truly meshless Galerkin method based on a moving least squares quadrature Marc Duflot, Hung Nguyen-Dang Abstract A new body integration technique is presented and applied to the evaluation of the stiffness
More informationName of the Student: Unit I (Solution of Equations and Eigenvalue Problems)
Engineering Mathematics 8 SUBJECT NAME : Numerical Methods SUBJECT CODE : MA6459 MATERIAL NAME : University Questions REGULATION : R3 UPDATED ON : November 7 (Upto N/D 7 Q.P) (Scan the above Q.R code for
More informationSolving the One Dimensional Advection Diffusion Equation Using Mixed Discrete Least Squares Meshless Method
Proceedings of the International Conference on Civil, Structural and Transportation Engineering Ottawa, Ontario, Canada, May 4 5, 2015 Paper No. 292 Solving the One Dimensional Advection Diffusion Equation
More informationA Fictitious Time Integration Method (FTIM) for Solving Mixed Complementarity Problems with Applications to Non-Linear Optimization
Copyright 2008 Tech Science Press CMES, vol.34, no.2, pp.155-178, 2008 A Fictitious Time Integration Method (FTIM) for Solving Mixed Complementarity Problems with Applications to Non-Linear Optimization
More informationExplicit Jump Immersed Interface Method: Documentation for 2D Poisson Code
Eplicit Jump Immersed Interface Method: Documentation for 2D Poisson Code V. Rutka A. Wiegmann November 25, 2005 Abstract The Eplicit Jump Immersed Interface method is a powerful tool to solve elliptic
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
More informationFractional Spectral and Spectral Element Methods
Fractional Calculus, Probability and Non-local Operators: Applications and Recent Developments Nov. 6th - 8th 2013, BCAM, Bilbao, Spain Fractional Spectral and Spectral Element Methods (Based on PhD thesis
More informationPhysics 202 Laboratory 5. Linear Algebra 1. Laboratory 5. Physics 202 Laboratory
Physics 202 Laboratory 5 Linear Algebra Laboratory 5 Physics 202 Laboratory We close our whirlwind tour of numerical methods by advertising some elements of (numerical) linear algebra. There are three
More informationHigh-order solution of elliptic partial differential equations in domains containing conical singularities
High-order solution of elliptic partial differential equations in domains containing conical singularities Thesis by Zhiyi Li In Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationNumerical Analysis Comprehensive Exam Questions
Numerical Analysis Comprehensive Exam Questions 1. Let f(x) = (x α) m g(x) where m is an integer and g(x) C (R), g(α). Write down the Newton s method for finding the root α of f(x), and study the order
More informationNonlinear Diffusion. Journal Club Presentation. Xiaowei Zhou
1 / 41 Journal Club Presentation Xiaowei Zhou Department of Electronic and Computer Engineering The Hong Kong University of Science and Technology 2009-12-11 2 / 41 Outline 1 Motivation Diffusion process
More informationQuantum Dynamics. March 10, 2017
Quantum Dynamics March 0, 07 As in classical mechanics, time is a parameter in quantum mechanics. It is distinct from space in the sense that, while we have Hermitian operators, X, for position and therefore
More informationMcGill University Department of Mathematics and Statistics. Ph.D. preliminary examination, PART A. PURE AND APPLIED MATHEMATICS Paper BETA
McGill University Department of Mathematics and Statistics Ph.D. preliminary examination, PART A PURE AND APPLIED MATHEMATICS Paper BETA 17 August, 2018 1:00 p.m. - 5:00 p.m. INSTRUCTIONS: (i) This paper
More informationIndex. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems
Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,
More informationCS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations
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 informationHIGH ACCURACY NUMERICAL METHODS FOR THE SOLUTION OF NON-LINEAR BOUNDARY VALUE PROBLEMS
ABSTRACT Of The Thesis Entitled HIGH ACCURACY NUMERICAL METHODS FOR THE SOLUTION OF NON-LINEAR BOUNDARY VALUE PROBLEMS Submitted To The University of Delhi In Partial Fulfillment For The Award of The Degree
More informationInverse Gravimetry Problem
Inverse Gravimetry Problem Victor Isakov September 21, 2010 Department of M athematics and Statistics W ichita State U niversity W ichita, KS 67260 0033, U.S.A. e mail : victor.isakov@wichita.edu 1 Formulation.
More informationUnconstrained Multivariate Optimization
Unconstrained Multivariate Optimization Multivariate optimization means optimization of a scalar function of a several variables: and has the general form: y = () min ( ) where () is a nonlinear scalar-valued
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 informationMOL 410/510: Introduction to Biological Dynamics Fall 2012 Problem Set #4, Nonlinear Dynamical Systems (due 10/19/2012) 6 MUST DO Questions, 1
MOL 410/510: Introduction to Biological Dynamics Fall 2012 Problem Set #4, Nonlinear Dynamical Systems (due 10/19/2012) 6 MUST DO Questions, 1 OPTIONAL question 1. Below, several phase portraits are shown.
More informationPoint interpolation method based on local residual formulation using radial basis functions
Structural Engineering and Mechanics, Vol. 14, No. 6 (2002) 713-732 713 Point interpolation method based on local residual formulation using radial basis functions G.R. Liu, L. Yan, J.G. Wang and Y.T.
More informationLeast-Squares Finite Element Methods
Pavel В. Bochev Max D. Gunzburger Least-Squares Finite Element Methods Spri ringer Contents Part I Survey of Variational Principles and Associated Finite Element Methods 1 Classical Variational Methods
More informationMethod of Fundamental Solutions for Scattering Problems of Electromagnetic Waves
Copyright c 2005 Tech Science Press CMES, vol.7, no.2, pp.223-232, 2005 Method of Fundamental Solutions for Scattering Problems of Electromagnetic Waves D.L. Young 1,2 and J.W. Ruan 2 Abstract: The applications
More informationA Lie-Group Shooting Method Estimating Nonlinear Restoring Forces in Mechanical Systems
Copyright 2008 Tech Science Press CMES, vol.35, no.2, pp.157-180, 2008 A Lie-Group Shooting Method Estimating Nonlinear Restoring Forces in Mechanical Systems Chein-Shan Liu 1 Abstract: For an inverse
More informationNewton Method with Adaptive Step-Size for Under-Determined Systems of Equations
Newton Method with Adaptive Step-Size for Under-Determined Systems of Equations Boris T. Polyak Andrey A. Tremba V.A. Trapeznikov Institute of Control Sciences RAS, Moscow, Russia Profsoyuznaya, 65, 117997
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationHIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS
HIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS JASON ALBRIGHT, YEKATERINA EPSHTEYN, AND QING XIA Abstract. Highly-accurate numerical methods that can efficiently
More informationOrdinary differential equations - Initial value problems
Education has produced a vast population able to read but unable to distinguish what is worth reading. G.M. TREVELYAN Chapter 6 Ordinary differential equations - Initial value problems In this chapter
More informationRadial point interpolation based finite difference method for mechanics problems
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 26; 68:728 754 Published online 11 April 26 in Wiley InterScience (www.interscience.wiley.com). DOI: 1.12/nme.1733
More informationRank-Deficient Nonlinear Least Squares Problems and Subset Selection
Rank-Deficient Nonlinear Least Squares Problems and Subset Selection Ilse Ipsen North Carolina State University Raleigh, NC, USA Joint work with: Tim Kelley and Scott Pope Motivating Application Overview
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More informationWavelet solution of a class of two-dimensional nonlinear boundary value problems
Copyright 2013 Tech Science Press CMES, vol.92, no.5, pp.493-505, 2013 Wavelet solution of a class of two-dimensional nonlinear boundary value problems Xiaojing Liu 1, Jizeng Wang 1,2, Youhe Zhou 1,2 Abstract:
More informationGeometric Modeling Summer Semester 2010 Mathematical Tools (1)
Geometric Modeling Summer Semester 2010 Mathematical Tools (1) Recap: Linear Algebra Today... Topics: Mathematical Background Linear algebra Analysis & differential geometry Numerical techniques Geometric
More informationAffine covariant Semi-smooth Newton in function space
Affine covariant Semi-smooth Newton in function space Anton Schiela March 14, 2018 These are lecture notes of my talks given for the Winter School Modern Methods in Nonsmooth Optimization that was held
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 informationMeshfree Approximation Methods with MATLAB
Interdisciplinary Mathematical Sc Meshfree Approximation Methods with MATLAB Gregory E. Fasshauer Illinois Institute of Technology, USA Y f? World Scientific NEW JERSEY LONDON SINGAPORE BEIJING SHANGHAI
More informationDiscrete Projection Methods for Integral Equations
SUB Gttttingen 7 208 427 244 98 A 5141 Discrete Projection Methods for Integral Equations M.A. Golberg & C.S. Chen TM Computational Mechanics Publications Southampton UK and Boston USA Contents Sources
More informationOrdinary Differential Equations
CHAPTER 8 Ordinary Differential Equations 8.1. Introduction My section 8.1 will cover the material in sections 8.1 and 8.2 in the book. Read the book sections on your own. I don t like the order of things
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 34: Improving the Condition Number of the Interpolation Matrix Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
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 informationState Space and Hidden Markov Models
State Space and Hidden Markov Models Kunsch H.R. State Space and Hidden Markov Models. ETH- Zurich Zurich; Aliaksandr Hubin Oslo 2014 Contents 1. Introduction 2. Markov Chains 3. Hidden Markov and State
More informationChapter 9. Derivatives. Josef Leydold Mathematical Methods WS 2018/19 9 Derivatives 1 / 51. f x. (x 0, f (x 0 ))
Chapter 9 Derivatives Josef Leydold Mathematical Methods WS 208/9 9 Derivatives / 5 Difference Quotient Let f : R R be some function. The the ratio f = f ( 0 + ) f ( 0 ) = f ( 0) 0 is called difference
More informationSpotlight on Laplace s Equation
16 Spotlight on Laplace s Equation Reference: Sections 1.1,1.2, and 1.5. Laplace s equation is the undriven, linear, second-order PDE 2 u = (1) We defined diffusivity on page 587. where 2 is the Laplacian
More informationNumerical Mathematics
Alfio Quarteroni Riccardo Sacco Fausto Saleri Numerical Mathematics Second Edition With 135 Figures and 45 Tables 421 Springer Contents Part I Getting Started 1 Foundations of Matrix Analysis 3 1.1 Vector
More informationNONLINEAR DIFFUSION PDES
NONLINEAR DIFFUSION PDES Erkut Erdem Hacettepe University March 5 th, 0 CONTENTS Perona-Malik Type Nonlinear Diffusion Edge Enhancing Diffusion 5 References 7 PERONA-MALIK TYPE NONLINEAR DIFFUSION The
More informationME Computational Fluid Mechanics Lecture 5
ME - 733 Computational Fluid Mechanics Lecture 5 Dr./ Ahmed Nagib Elmekawy Dec. 20, 2018 Elliptic PDEs: Finite Difference Formulation Using central difference formulation, the so called five-point formula
More informationis represented by a convolution of a gluon density in the collinear limit and a universal
. On the k? dependent gluon density in hadrons and in the photon J. Blumlein DESY{Zeuthen, Platanenallee 6, D{15735 Zeuthen, Germany Abstract The k? dependent gluon distribution for protons, pions and
More information17th European Conference on Fracture 2-5 September,2008, Brno, Czech Republic. Thermal Fracture of a FGM/Homogeneous Bimaterial with Defects
-5 September,8, Brno, Czech Republic Thermal Fracture of a FGM/Homogeneous Bimaterial with Defects Vera Petrova, a, Siegfried Schmauder,b Voronezh State University, University Sq., Voronezh 3946, Russia
More informationNonlinear Oscillations and Chaos
CHAPTER 4 Nonlinear Oscillations and Chaos 4-. l l = l + d s d d l l = l + d m θ m (a) (b) (c) The unetended length of each spring is, as shown in (a). In order to attach the mass m, each spring must be
More informationHIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS
HIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS JASON ALBRIGHT, YEKATERINA EPSHTEYN, AND QING XIA Abstract. Highly-accurate numerical methods that can efficiently
More informationNonlinear Evolution of a Vortex Ring
Nonlinear Evolution of a Vortex Ring Yuji Hattori Kyushu Institute of Technology, JAPAN Yasuhide Fukumoto Kyushu University, JAPAN EUROMECH Colloquium 491 Vortex dynamics from quantum to geophysical scales
More informationM445: Heat equation with sources
M5: Heat equation with sources David Gurarie I. On Fourier and Newton s cooling laws The Newton s law claims the temperature rate to be proportional to the di erence: d dt T = (T T ) () The Fourier law
More informationStochastic Collocation Methods for Polynomial Chaos: Analysis and Applications
Stochastic Collocation Methods for Polynomial Chaos: Analysis and Applications Dongbin Xiu Department of Mathematics, Purdue University Support: AFOSR FA955-8-1-353 (Computational Math) SF CAREER DMS-64535
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationH. L. Atkins* NASA Langley Research Center. Hampton, VA either limiters or added dissipation when applied to
Local Analysis of Shock Capturing Using Discontinuous Galerkin Methodology H. L. Atkins* NASA Langley Research Center Hampton, A 68- Abstract The compact form of the discontinuous Galerkin method allows
More information