THE ACCURATE ELEMENT METHOD: A NOVEL METHOD FOR INTEGRATING ORDINARY DIFFERENTIAL EQUATIONS
|
|
- Arron Hicks
- 6 years ago
- Views:
Transcription
1 AAS THE ACCURATE EEMENT METHOD: A NOVE METHOD FOR INTEGRATING ORDINARY DIFFERENTIA EQUATIONS Paul Cizmas and Maty Blumenfeld This paper presents the development of a new methodology for numerically solving ordinary differential equations. This methodology, which we call the accurate element method (AEM), breaks the connection between the solution accuracy and the number of unknowns. The differential equations are discretized by using finite elements, like in the finite element method (FEM). A key feature of the AEM is the methodology developed for eliminating unknowns inside the element by using the relations provided by the governing equations. The results of the AEM applied to a second-order ordinary differential equation (ODE) show that the degree of the local approximation function does not affect the number of unknowns used to discretize the differential equation. For a second-order ODE it is shown that the discretized solution has the same number of unknowns, whether the local approximation (or interpolation) function is a third-degree or a nineteen-degree polynomial. The implication of this result is that higher-degree interpolation functions can be used without increasing the number of unknowns. For a specified accuracy, the usage of elements with high-degree interpolation functions leads to a reduction in the necessary number of elements. Numerical results using the AEM, the shooting method and the relaxation method are presented and compared for several ODE boundary value problems. INTRODUCTION Numerical solution of differential equations is of paramount importance for most of the fields of engineering and science because it provides a way to investigate the governing equations of various physical phenomena. As the field of numerical methods for solving differential equations matured, some of the algorithms were preferred for various reasons and are currently part of the classical methodologies. If one restricts this discussion to ordinary differential equations, it is easy to identify the most popular methods of the moment: () Runge-Kutta methods, () Richardson extrapolation and its particular implementation as Bulirsch-Stoer method, and (3) predictor-corrector methods, 4 p.70. The most general and widely used methods for partial differential equations include the finite difference method, finite volume method and finite element method (FEM). The common Assistant Professor, Department of Aerospace Engineering, Texas A&M University, College Station, TX cizmas@tamu.edu. Professor Emeritus, Chair of Strength of Materials, Politehnica University Bucharest, Romania, E- mail:mblum@itcnet.ro.
2 goal of these methods is to approximate as well as possible the solution of the governing equations describing various physical phenomena. Several options are currently used to improve the accuracy of these methods: () increasing the degree of the polynomials that locally approximate the solution, that is, p-type refinement; () decreasing the size of the elements which discretize the computational domain in order to reduce the variation of the unknowns over the element, that is, h-type refinement; (3) using anisotropic elements, that is, having different orders of expansion for each spatial direction in an element; (4) using non-conforming elements with either h-type or p-type non-conformities; (5) using adaptive refinement strategies. High-order finite element methods 9, 6 and spectral element methods 7, 0, 3, 7 are commonly used to produce high accuracy solutions. Spectral element methods combine the high accuracy of spectral methods 6, and the flexibility of the finite element method. The Galerkin spectral element method, a subset of the spectral element method, is similar to the high-order finite element method in the sense that both are weighted residual techniques. A basic difference between the two methods is the fact that the weighting functions used in the Galerkin spectral element method are either Chebyshev or egendre polynomials. 5, 6 Anisotropic and non-conforming elements are also used to improve the solution accuracy of certain problems. The non-conforming elements can be either h-type or p-type non-conforming. Various special formulations have been developed for anisotropic and non-conforming elements, such as the constrained approximation method, 8,, 5 the mortar element method,, the FETI method, 9 the transition element method 3, and the iterative patching method. 0, 5 Another option for producing high accuracy solutions is to use adaptive refinement strategies. Adaptive methods have the ability to automatically adjust themselves to improve the solution of a problem. The adaptive methods are mainly used for: () problems for which a properly resolved initial mesh is not available; and () unsteady problems for which the initial mesh is not good enough at some instants of the solution. The local error estimation and the adaptive strategy are two important aspects of the adaptive methods. 7, 4, 8 For all the methods currently used to solve differential equations, the outcome of increasing solution accuracy is an increase in the number of unknowns. To minimize the number of unknowns, each method tries to optimize the relation between the degree of the approximation function and the size of the element, such as in the h p adaptive finite element or spectral element methods. This paper presents a new methodology for the numerical solution of differential equations which breaks the connection between the solution accuracy and the number of unknowns. Consequently, the number of unknowns corresponding to high-accuracy solutions obtained using high-order approximation polynomials will be equal to the number of unknowns corresponding to low-order approximation polynomials. This new methodology, which we call the accurate element method (AEM), will be applied herein to numerically solve linear and nonlinear ODEs. The results shown in this paper indicate that the degree of the local approximation function does not affect the number of unknowns used to discretize the differential equation. THE ACCURATE EEMENT METHOD (AEM) The AEM will be presented by comparing it to the FEM. The strategies of both the AEM and the FEM can be summarized by the following three steps: () discretization of the domain D on which the governing equations must be integrated; () local approximation of the solution of the governing equations, in which the information given by each element is concentrated at the nodes of the element; and (3) reconstruction of the domain D, obtained by writing nodal equations that bring together the information given by the elements adjacent to each node. The solution of the
3 system of equations obtained in this way represents the nodal unknowns. The first and third steps of the AEM are identical to those of the FEM. The two methods differ only in how the second step of their common strategy is accomplished. Note that the second step is crucial in obtaining an accurate solution of the governing equations. As will be shown herein, the AEM uses a much smaller number of elements than the FEM to obtain a similar accuracy. In step two, the local approximation of the solution in the FEM is done by using interpolation functions (IF). Usually these interpolation functions are low degree polynomials that lead to poor approximations. Consequently, the number of elements used in the FEM must be large in order to obtain a result with reasonable accuracy. The local approximation in the AEM will be done by high-degree polynomial functions. Because the functions used for local approximation are highdegree interpolation functions (HIF), only a small number of elements is needed to obtain an accurate solution in the AEM. The example presented in the next section shows that the accuracy obtained in the AEM by using two elements with a seventh-degree polynomial HIF is better than the accuracy obtained in the FEM by using 00 elements with a linear IF. The degree of the polynomials used in the AEM is usually higher than seven, since, as will be shown herein, the computational effort in the AEM is not affected by the degree of the polynomial. This is an important difference between the AEM and the FEM. High-degree polynomials may be used for the IF in the FEM, but this increases the number of unknowns of the discretized system of equations. In the AEM, the degree of the interpolation function used in the local approximation does not affect the number of unknowns of the discretized system of equations. This is an important advantage of the AEM compared to the finite element, finite volume or finite difference methods. To illustrate how the AEM can be implemented to solve ODEs, three cases will be presented: () a linear ODE with constants coefficients; () a linear ODE with variable coefficients; (3) a nonlinear ODE. All cases shown herein are second-order ODEs. AEM Applied to inear ODEs To introduce the methodology of the AEM, let us start by presenting a simple but yet meaningful case, the solution of a linear, second-order ordinary differential equation: d φ + cφ(x) + d(x) = 0, () where c is a constant and d(x) is a function of the independent variable x. et us consider the following boundary conditions φ := φ(x = 0) = 0 () φ := φ(x = ) = 0. et us consider a HIF that is a seventh-degree polynomial: φ(x) = C 0 + C x + C x + C 3 x C 7 x 7 (3) and let us write φ using a natural coordinate η = x/, where is the length of the element. Equation (3) can be written in vectorial form using the natural coordinate η: φ(η) = ηc, (4) 3
4 where η = η η... η 7 C = C 0 C C... C 7 7 T. The first goal is to write the local interpolation function φ(x) as a function of the values of φ(x) and its derivatives at nodal locations. The values of the function φ and its derivatives at both ends of the element can be written as: or φ () φ () 3 φ (3) φ () φ () 3 φ (3) = C (5) B = AC. (6) The notation i := φ i and φ (k) i B matrix can be rearranged as: 0 B = 0 + := dk φ i k, i, and k has been used in equation (5). The 0 0 φ () φ () φ () φ () and the inverse of the A matrix needs to be computed only once. The vector of coefficients C can be written using equations (6) and (7) as: C = A B = F 80 φ () + F 8 φ () + φ () F 8 φ () + 3 φ (3) F 83 φ (3) (8) φ (3) φ (3) (7) where F 80 = , F 8 = , F 8 = , 4
5 F 83 = By substituting the C vector from equation (8) into equation (4), the local approximation function φ(x) is written as a function of the nodal values of the function and its derivatives. et us derive a relation similar to but more general than the stiffness matrix relation. This and the variable derivatives φ (k) i, k. relation should provide a link between the variables i This relation, which is a generalization of the stiffness matrix relation, will be called the complete transfer relation (CTR). 3 The derivation of the CTR can be done based on the physical phenomenon that is described by the governing equations or by using the Galerkin method. The latter approach will be used herein, because of its generality. et us use the following weighting function w(x) = x φ + x φ = N φ φ. (9) Using the Galerkin projection one obtains 0 N T d φ + c 0 N T φ(x) + After integrating by parts the first term, equation (0) becomes + 0 φ () 0 φ () c x 0 x = 0 x x 0 N T d(x) = 0. (0) φ(x) = d(x). () et us now replace the function φ of the third term of the left-hand side by the expression given in equation (4), where the vector of coefficients C is written as in equation (8). Equation () becomes where c 3 G0F 8 G0F 0 8 = 36 φ () φ () G0F 8 = c 4 G0F 3 8 φ () φ () φ (3) φ (3), G0F 8 = 5 cg0f 0 8, G0F 3 8 = 50 = x x,. c G0F 8 φ () φ () d(x) () Note that the third term of left-hand side of equation () is responsible for introducing errors in FEM, because the function φ in the integral is approximated by a low degree polynomial. This 5
6 term is responsible for introducing errors in the FEM. Unlike the FEM, the AEM uses a high degree polynomial to approximate φ(x). The result of this approach is that equation () has 8 unknowns,..., φ(3). There are, however, only 4 available relations: equations () and boundary conditions (8). Consequently, 4 additional relations are necessary. Usually, when more information is needed, this information is obtained either from some relations with other elements or by accepting a reasonable approximation. Neither of these options are used in the AEM. Instead, the necessary information is obtained accurately from inside the, φ(0) element, that is, from the governing equations themselves. derivatives are φ () φ () = c d (0) d (0) Consequently, the second and third (3) and φ (3) φ (3) = c φ () φ () d () d (), (4) where d (0) := d(x = 0), d (0) := d(x = ), d () := d d(x) x=0, d () := d d(x) x=. Consequently, the equation () becomes A φ () + B φ () = C, (5) where A = c G0F c 3 G0F 8 B = 0 0 c G0F 8 + c 4 G0F 3 8 (6) C = x 0 x d(x) c 3 d (0) G0F 8 d (0) c 4 d () G0F 38 d () The unknowns φ (), φ(), φ(3) and φ (3), which we call the apparent unknowns,4 have been eliminated by using the governing equations. To illustrate the methodology, let us consider that in equation () c = 4 and d(x) = 4x 3 8x 8x 8, that is d φ + 4φ 4x3 8x 8x 8 = 0 (7) with the boundary conditions φ := φ(x = 0) = 0 (8) φ := φ(x = ) = 0. 6
7 The analytical solution of equation (7) with the specified boundary conditions is: φ(x) = (cos() 7)/ sin() sin(x) cos(x) + x 3 + x + 3x +. (9) The values of the A, B and C matrices corresponding to the seventh- and nineteen-degree interpolation function are obtained from equations (6): HIF8 A = B = C = / / HIF0 A = B = C = For comparison, let us show the A, B and C matrices obtained using a linear interpolation function (IF) A = /3 5/3 5/3 /3 B = 0 0 C = 3.6/3.8. Table shows the variation of φ () as a function of the degree of the polynomial in the HIF. This solution is based on one element. For comparison, the value of φ () calculated using the FEM with an IF is also included. The exact value is obtained from the solution (9) as φ () = The variation of the error shows that an accurate solution can be obtained with only one AEM element with HIF. φ () Error % IF HIF HIF e-3 HIF e-7 HIF e- HIF e-4 exact Table Solution based on one element. The variation of φ () as a function of the number of elements is shown in Table. Results are shown for the FEM with IF and the AEM with HIF4 and HIF8. The advantage of using the AEM 7
8 and high-degree polynomials is obvious: two AEM elements with HIF8 produce more accurate results than 00 FEM elements with IF. Four AEM elements with HIF8 produce a solution with an error of the order of 0 8 %. IF HIF4 HIF8 Number Elem. φ () Error % φ () Error % Error % e e e e- -.84e e-3 * e-4 * e-5 * e e-7 * e-3 * * * Table Variation of the value of φ () as a function of the number of elements and degree of the HIF/IF. AEM Applied to Nonlinear ODEs This section presents the methodology for solving nonlinear ODEs using the AEM. In order to solve nonlinear ODEs using the AEM, a methodology for solving linear ODEs with variable coefficients must be developed first. For this reason, the first part of this section is devoted to solving linear ODEs with variable coefficients. The second part of this section presents the methodology for solving nonlinear ODEs using the AEM. et us consider the following second-order ODE with variable coefficients a(x) d φ + b(x)dφ + c(x)φ(x) + d(x) = 0. (0) The CTR corresponding to equation (0) will be derived using the Galerkin method. reason, it is convenient to rewrite equation (0) in a transformed form For this d a(x) dφ + d ( b(x) da(x) ) ( φ + c(x) db(x) ) + d a(x) φ + d(x) = 0. () et us use for the Galerkin projection the weighting function w(x) from equation (9). The result of the Galerkin projection applied to equation () is: T d N a(x) dφ ( T d + N b(x) da(x) ) φ + T N (c(x) db(x) ) + d a(x) T φ + N d(x) = 0. () For brevity, only the derivation of the first term of equation () will be presented herein. The same methodology is applied for the derivation of the other terms in equation (). The first term of equation () will be calculated using integration by parts: T = T d N a(x) dφ d T N = a(x)dφ + T dφ N a(x). (3) The first term of the right-hand side is denoted T a and the second term is denoted T b. Taking 8
9 into account the expression of N from equation (9), the term T a becomes d T N T a = a(x)dφ = a(x) dφ. Integrating by parts again, T a becomes T a = { } da(x) φ + (a()φ a(0)φ ). Using the notations introduced in equation (5) and a (k) := dk a(0) and a (k) k := dk a(), the term T k becomes T = da(x) φ a (0) a (0) a (0) a (0) a (0) 0 φ () 0 a (0) φ (). A similar approach is used for the other terms in equation (), such that the general form of the CTR is a (0) / + b(0) a () a (0) a (0) / a(0) / b(0) a () a (0) + 0 φ () 0 a (0) + ( b(x) da(x) ) φ ( c(x) db(x) + d a(x) ) N T φ = φ () N T.(4) Having derived the CTR for linear ODEs with variable coefficients, the next step is to produce an algorithm for solving nonlinear ODEs. To present the methodology for solving nonlinear ODEs, let us choose the following second-order, nonlinear ODE φ () + (3 + 4φ () 4) + (8 40x + 85x 6x 3 + x 4 + x 5 3x 6 ) = 0 (5) with the boundary conditions := φ(x = ) = 0 and φ () := dφ (x = ) = 0. For simplicity, let us assume =. The nonlinearity is due to the 3 and 4φ () terms. In order to eliminate the nonlinearity, the terms in the parentheses (3 + 4φ () 4) will be approximated using HIF. To simplify the algebra, a third-degree polynomial HIF will be used in this example. The derivation of the vector of coefficients C is similar to the derivation presented in section, where C is given by equation (8). For a third-degree polynomial HIF, the expression of the vector of coefficients C is C = F 40 + F 4 φ () φ (), (6) where F 40 = and F 4 = 0. 9
10 Using equation (4), the HIF turns out to be the well known Hermite function φ(x) = 3 3x + x x x x x + x 3 φ () + x + x 3 φ (). (7) Replacing the terms in the parenthesis (3 + 4φ () 4) by the HIF from equation (7), equation (5) becomes where φ () + (K 0 + K x + K x + K 3 x 3 ) + (8 40x + 85x 6x 3 + x 4 + x 5 3x 6 ) = 0, (8) K 0 = 7 4, K = 4 3φ (), K = 5 + 6φ (), K 3 = 6 + 3φ (). (9) As a result, the nonlinear ordinary differential equation (5) has been replaced by a linear ODE with variable coefficients (8). The CTR of equation (8) is derived using the general form obtained in equation (4). Having b(x) = 0 and =, equation (4) becomes = 0 ( K0 + K x + K x + K 3 x 3) T N φ(x) = φ () φ () = 0 + d(x) N T. (30) After substituting the HIF, φ(x), from equation (7), the K i coefficients from equation (9) and d(x) from equation (5), the CTR equation (30) becomes ( + 0 φ () ) + /5 0 7/ /0 3/0 φ() + ( ) 3/840 φ () /0 + 49/0 /0 + 7/0 φ () 7/5 357/840 = 0, (3) 333/840 or using the notation u := and v := φ () { 336u 356uv 3v u + 48v = u + 9uv + 7v + 94u + 39v = 999 (3) Two sets of real solutions are obtained by solving the nonlinear system of equations (3) u =, v = 5 and u = 5.8, v = Only the first set of solutions leads to an acceptable solution which is φ(x) = (x 3 3x + ) + (x 3 x + x)φ () = (x 3 3x + )( ) + (x 3 x + x)(5) = x 3 4x + 5x. (33) The second set of solutions (u := = 5.8, v := φ () = 60.3) produces the solution φ(x) = 9.57(x ) (x 0.799) which does not satisfy the governing equation (5). 0
11 AEM vs. Existing ODE Solvers To compare the AEM against existing ODE solvers, let us consider the second-order ODE with variable coefficients: (+x+3x +3x 3 +x 4 ) dφ +( 5x+x +3x 3 +4x 4 ) dφ +(+x+x +3x 3 +4x 4 )φ+f(x) = 0, where f(x) has been chosen to be f(x) = x + 8x + 66x 3 59x 4 8x 5 98x 6 75x 7 44x 8 053x 9 056x 0 447x 60x 6x 3, such that the solution of the ODE is the ninth-degree polynomial φ(x) = 3 + 5x 4x + x 3 6x 4 + x 5 + x 6 + x 7 + x 8 + 4x 9. Various integration lengths have chosen with various two-point boundary conditions, as shown in Table 3. The AEM has been compared against the shooting method and the relaxation method. All the methods were implemented in FORTRAN. The shooting method used the Runge-Kutta method implemented in the odeint routine. 4 One hundred years after they have been introduced, 8 the Runge-Kutta methods are still the workhorse for solving ODEs, 4 p There were several improvements of the Runge-Kutta method, such as Runge-Kutta-Gill and (3/8) Runge-Kutta, but these modifications did not fundamentally changed the basic idea of the method. Other method for solving ODEs, such as the Bulirsch-Stoer method and the predictor-corrector method, can be faster than the Runge-Kutta in certain circumstances, but the differences are rather limited. Given these facts and due to its robustness, the fourth-order Runge-Kutta method has been selected as one of the methods to be compared against the AEM. The relaxation method used the solvde routine. 4 The computations were done on an SGI Indigo computer, with an R8000 processor. The results shown in Fig. indicate that for the same computational time, the accuracy of the AEM results with a 3th-degree polynomial (HIF4) was approximately five orders of magnitude better than the accuracy of the shooting and relaxation methods. In addition, for the same error level, the AEM with HIF4 interpolation functions produced solutions at least two orders of magnitude faster than the shooting and relaxation methods. The relaxation method diverged when the integration was equal or exceeded x =. Domain Boundary Values Case x left x right φ left φ right Table 3 Integration cases. CONCUSIONS A new method, which we call the accurate element method, has been presented herein and applied to linear and nonlinear ODEs. A key feature of the AEM was the methodology developed for eliminating unknowns inside the element by using the relations provided by the governing equations. The results of the AEM applied to a second-order ODE showed that the order of the local approximation function did not affect the number of unknowns used to discretize the differential equation. This is a breakthrough in the paradigm used to approach the numerical solution of differential equations. For a second-order ODE it was shown that the discretized solution has the same number of unknowns, whether the local approximation (or interpolation) function is a
12 0<=x<= 0<=x<= log0(error) % AEM, CF0 AEM, CF4 Relax, eps=5e 6 Relax, eps=5e 7 Relax, eps=5e 9 Shoot, eps=e, h=e 8 Shoot, eps=e, h=e 6 Shoot, eps=e 8, h=e 6 log0(error) % AEM, CF0 AEM, CF4 Shoot, eps=(e 8,e 4), h=e log0(time) sec (i) log0(time) sec (ii) 0<=x<=3 0<=x<= log0(error) % log0(error) % AEM, CF0 AEM, CF4 Shoot, eps=(e,e 8), h=e 8 0 AEM, CF0 AEM, CF4 Shoot, eps=(e,e 8), h=e log0(time) sec (iii) log0(time) sec (iv) Figure AEM vs. Relaxation and Shooting Methods
13 third-degree or a nineteen-degree polynomial. The implication of this result is that higher-degree interpolation functions can be used without increasing the number of unknowns. The AEM has been used to solve linear ODEs with constant and variable coefficients and nonlinear ODEs. For a linear ODE with variable coefficients, the AEM has been compared against the relaxation method and the shooting method. The accuracy of the AEM results with a 3thdegree polynomial (HIF4) was approximately five orders of magnitude better than the accuracy of the shooting and relaxation methods. In addition, for the same error level, the AEM with HIF4 interpolation functions produced solutions at least two orders of magnitude faster than traditional ODE solvers. REFERENCES. C. Bernardi, Y. Maday, C. Mavriplis, and A. T. Patera. The mortar element method applied to spectral discretizations. In T. J. Chung and G. R. Karr, editors, Finite Element Analysis in Fluids, The Seventh Int. Conf. on Finite Element Methods in Flow Problems. UAH Press, C. Bernardi, Y. Maday, and A. T. Patera. A new nonconforming approach to domain decomposition: The mortar element method. In H. Brezis and J.. ions, editors, Nonlinear Partial Differential Equations and Their Applications, number 99 in Pitman Research Notes in Mathematics. Wiley, M. Blumenfeld. A New Method for Accurate Solving of Ordinary Differential Equations. Editura Tehnica, Bucharest, M. Blumenfeld and P. G. A. Cizmas. The accurate element method: A new paradigm for numerical solution of ordinary differential equations. submitted for publication to the Journal of the Romanian Academy, January J. P. Boyd. Chebyshev and Fourier Spectral Methods. Springer, New York, C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang. Spectral Methods in Fluid Dynamics. Springer, New York, G. F. Carey and A. I. Pehlivanov. Computer methods in applied mechanics and engineering. Computer Methods in Applied Mechanics and Engineering, 50(-4):5 3, Dec R. Y. Chang and C. H. Hsu. A variable-order spectral element method for incompressible viscous flow simulation. International Journal for Numerical Methods in Engineering, 39(7): , Demkowicz, J. T. Oden, W. Rachowicz, and O. Hardy. Toward a universal h p adaptive finite element strategy: Part : Constrained approximation and data structure. Computer Methods in Applied Mechanics and Engineering, 77(-):79, D. Funaro, A. Quarteroni, and P. Zanoli. An iterative procedure with interface relaxation for domain decomposition methods. SIAM Journal on Numerical Analysis, 5(6):3 36, Dec S. Gill. A process for the step-by-step integration of differential equations in an automatic digital computing machine. Proc. Cambridge Philos. Soc., 95.. D. Gottlieb and S. A. Orszag. Numerical Analysis of Spectral Methods: Theory and Applications. SIAM-CMBS, Philadelphia,
14 3. A. K. Gupta. A finite element method for transition from a fine to a coarse grid. Int. J. Numer. Meth. Eng., :395 4, R. D. Henderson. Dynamic refinement algorithms for spectral element methods. Computer Methods in Applied Mechanics and Engineering, 75(3-4):395 4, R. D. Henderson and G. E. Karniadakis. Hybrid spectral element-low order methods for incompressible flows. J. Sci. Comp., 6():79, G. E. Karniadakis and S. J. Sherwin. Spectral/hp Element Methods for CFD. Oxford University Press, Oxford, K. Z. Korczak and A. T. Patera. An isotropic spectral element method for solution of the navier-stokes equations in complex geometry. Journal of Computational Physics, 6():36 38, Feb W. Kutta. Beitrag zur naherungsweisen integration totaler differentialgleichunge. Zeit. Math. Physik, 46: , C. acour and Y. Maday. Two different approaches for matching nonconforming grids: The mortar element method and the feti method. BIT, 37(3):70 738, Y. Maday and A. T. Patera. Spectral element methods for the incompressible navier-stokes equations. In A. K. Noor and J. T. Oden, editors, State of the Art Surveys in Computational Mechanics, chapter D. J. Morton, J. M. Tyler, and J. R. Dorroh. A new 3D finite element for adaptive h- refinement in -irregular meshes. International Journal for Numerical Methods in Engineering, 38(3): , Dec T. J. Oden and. Demkowicz. Advances in adaptive improvement: A survey of adaptive finite element methods in computational mechanics. In A. K. Noor and J. T. Oden, editors, State of the Art Surveys in Computational Mechanics, chapter A. T. Patera. A spectral method for fluid dynamics: aminar flow in a channel expansion. Journal of Computational Physics, 54(3): , W. H. Press, W. T. Vetterling, S. A. Teukolsky, and B. P. Flannery. Numerical Recipes in FORTRAN - The Art of Scientific Computing. Cambridge, second edition, W Rachowicz, J. T. Oden, and. Demkowicz. Toward a universal h p adaptive finite element strategy: Part iii: Design of h p meshes. Computer Methods in Applied Mechanics and Engineering, 77(-):8, J. N. Reddy. An Introduction to the Finite Element Method. McGraw-Hill, New York, second edition, E. M. Ronquist and A. T. Patera. Spectral element methods for unsteady navier-stokes equations. In Proceedings of the Seventh GAMM Conference on Numerical Methods in Fluid Mechanics, page 38. Viewig, W. Wang. Special quadratic quadrilateral finite elements for local refinement with irregular nodes. Computer Methods in Applied Mechanics and Engineering, 8(-):09 34, Feb
THE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS
THE PUBISHING HOUSE PROCEEDINGS O THE ROMANIAN ACADEMY, Series A, O THE ROMANIAN ACADEMY Volume, Number /, pp. 6 THE ACCURATE EEMENT METHOD: A NEW PARADIGM OR NUMERICA SOUTION O ORDINARY DIERENTIA EQUATIONS
More informationEigenvalues of Trusses and Beams Using the Accurate Element Method
Eigenvalues of russes and Beams Using the Accurate Element Method Maty Blumenfeld Department of Strength of Materials Universitatea Politehnica Bucharest, Romania Paul Cizmas Department of Aerospace Engineering
More informationModified Fourier Spectral Method for Non-periodic CFD
Modified Fourier Spectral Method for Non-periodic CFD Huankun Fu Chaoqun Liu Technical Report 2011-14 http://www.uta.edu/math/preprint/ Modified Fourier Spectral Method for Non-periodic CFD Huankun Fu
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More informationA FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS
Proceedings of ALGORITMY 2005 pp. 222 229 A FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS ELENA BRAVERMAN, MOSHE ISRAELI, AND ALEXANDER SHERMAN Abstract. Based on a fast subtractional
More informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
More informationOrdinary Differential Equations. Monday, October 10, 11
Ordinary Differential Equations Monday, October 10, 11 Problems involving ODEs can always be reduced to a set of first order differential equations. For example, By introducing a new variable z, this can
More informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY 11, USA Dan Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
More informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
More informationA Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning
Int. J. Contemp. Math. Sciences, Vol. 2, 2007, no. 22, 1097-1106 A Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning M. T. Darvishi a,, S.
More informationStability of the Parareal Algorithm
Stability of the Parareal Algorithm Gunnar Andreas Staff and Einar M. Rønquist Norwegian University of Science and Technology Department of Mathematical Sciences Summary. We discuss the stability of the
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationAn Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions
An Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions Leszek Marcinkowski Department of Mathematics, Warsaw University, Banacha
More informationNewton s Method and Efficient, Robust Variants
Newton s Method and Efficient, Robust Variants Philipp Birken University of Kassel (SFB/TRR 30) Soon: University of Lund October 7th 2013 Efficient solution of large systems of non-linear PDEs in science
More informationSpectral Collocation Method for Parabolic Partial Differential Equations with Neumann Boundary Conditions
Applied Mathematical Sciences, Vol. 1, 2007, no. 5, 211-218 Spectral Collocation Method for Parabolic Partial Differential Equations with Neumann Boundary Conditions M. Javidi a and A. Golbabai b a Department
More informationLeast-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations
Least-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations by Wilhelm Heinrichs Universität Duisburg Essen, Ingenieurmathematik Universitätsstr.
More informationMIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING
144 MIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING J. N. Reddy* and Chen-Shyh-Tsay School of Aerospace, Mechanical and Nuclear Engineering, University of Oklahoma, Norman, Oklahoma The paper describes
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationThe need for de-aliasing in a Chebyshev pseudo-spectral method
The need for de-aliasing in a Chebyshev pseudo-spectral method Markus Uhlmann Potsdam Institut für Klimafolgenforschung, D-442 Potsdam uhlmann@pik-potsdam.de (March 2) Abstract In the present report, we
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More informationSpectral Methods for Reaction Diffusion Systems
WDS'13 Proceedings of Contributed Papers, Part I, 97 101, 2013. ISBN 978-80-7378-250-4 MATFYZPRESS Spectral Methods for Reaction Diffusion Systems V. Rybář Institute of Mathematics of the Academy of Sciences
More informationMATHEMATICAL METHODS INTERPOLATION
MATHEMATICAL METHODS INTERPOLATION I YEAR BTech By Mr Y Prabhaker Reddy Asst Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad SYLLABUS OF MATHEMATICAL METHODS (as per JNTU
More informationarxiv: v1 [math.na] 1 May 2013
arxiv:3050089v [mathna] May 03 Approximation Properties of a Gradient Recovery Operator Using a Biorthogonal System Bishnu P Lamichhane and Adam McNeilly May, 03 Abstract A gradient recovery operator based
More informationNUMERICAL METHODS FOR ENGINEERING APPLICATION
NUMERICAL METHODS FOR ENGINEERING APPLICATION Second Edition JOEL H. FERZIGER A Wiley-Interscience Publication JOHN WILEY & SONS, INC. New York / Chichester / Weinheim / Brisbane / Singapore / Toronto
More informationA Posteriori Error Estimator for a Non-Standard Finite Difference Scheme Applied to BVPs on Infinite Intervals
arxiv:1503.00037v2 [math.na] 19 Mar 2015 A Posteriori Error Estimator for a Non-Standard Finite Difference Scheme Applied to BVPs on Infinite Intervals Riccardo Fazio and Alessandra Jannelli Department
More informationRATIONAL CHEBYSHEV COLLOCATION METHOD FOR SOLVING NONLINEAR ORDINARY DIFFERENTIAL EQUATIONS OF LANE-EMDEN TYPE
INTERNATIONAL JOURNAL OF INFORMATON AND SYSTEMS SCIENCES Volume 6, Number 1, Pages 72 83 c 2010 Institute for Scientific Computing and Information RATIONAL CHEBYSHEV COLLOCATION METHOD FOR SOLVING NONLINEAR
More informationNumerical solution of linear time delay systems using Chebyshev-tau spectral method
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 93-9466 Applications and Applied Mathematics: An International Journal (AAM) Vol., Issue (June 07), pp. 445 469 Numerical solution of linear time
More informationFinite Difference and Finite Element Methods
Finite Difference and Finite Element Methods Georgy Gimel farb COMPSCI 369 Computational Science 1 / 39 1 Finite Differences Difference Equations 3 Finite Difference Methods: Euler FDMs 4 Finite Element
More informationNUMERICAL MATHEMATICS AND COMPUTING
NUMERICAL MATHEMATICS AND COMPUTING Fourth Edition Ward Cheney David Kincaid The University of Texas at Austin 9 Brooks/Cole Publishing Company I(T)P An International Thomson Publishing Company Pacific
More informationPSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS Jie Shen Department of Mathematics, Penn State University University Par, PA 1680, USA Abstract. We present in this
More informationIntroduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
More informationPARTIAL DIFFERENTIAL EQUATIONS
MATHEMATICAL METHODS PARTIAL DIFFERENTIAL EQUATIONS I YEAR B.Tech By Mr. Y. Prabhaker Reddy Asst. Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad. SYLLABUS OF MATHEMATICAL
More informationEquations Introduction
Chapter 9 90 Introduction Partial Differential Equations The numerical treatment of partial differential equations is, by itself, a vast subject Partial differential equations are at the heart of many,
More informationNumAn2014 Conference Proceedings
OpenAccess Proceedings of the 6th International Conference on Numerical Analysis, pp 198-03 Contents lists available at AMCL s Digital Library. NumAn014 Conference Proceedings Digital Library Triton :
More informationAn integral-collocation-based fictitious domain technique for solving elliptic problems
An integral-collocation-based fictitious domain technique for solving elliptic problems N. Mai-Duy,, H. See and T. Tran-Cong Computational Engineering and Science Research Centre (CESRC) The University
More informationSolving Poisson s Equations Using Buffered Fourier Spectral Method
Solving Poisson s Equations Using Buffered Fourier Spectral Method Yinlin Dong Hassan Abd Salman Al-Dujaly Chaoqun Liu Technical Report 2015-12 http://www.uta.edu/math/preprint/ Solving Poisson s Equations
More informationUtilising high-order direct numerical simulation for transient aeronautics problems
Utilising high-order direct numerical simulation for transient aeronautics problems D. Moxey, J.-E. Lombard, J. Peiró, S. Sherwin Department of Aeronautics, Imperial College London! WCCM 2014, Barcelona,
More informationIntegration, differentiation, and root finding. Phys 420/580 Lecture 7
Integration, differentiation, and root finding Phys 420/580 Lecture 7 Numerical integration Compute an approximation to the definite integral I = b Find area under the curve in the interval Trapezoid Rule:
More informationApplied Numerical Analysis
Applied Numerical Analysis Using MATLAB Second Edition Laurene V. Fausett Texas A&M University-Commerce PEARSON Prentice Hall Upper Saddle River, NJ 07458 Contents Preface xi 1 Foundations 1 1.1 Introductory
More informationDe-aliasing on non-uniform grids: algorithms and applications
Journal of Computational Physics 191 (2003) 249 264 www.elsevier.com/locate/jcp De-aliasing on non-uniform grids: algorithms and applications Robert M. Kirby a, *, George Em Karniadakis b a School of Computing,
More information9.6 Predictor-Corrector Methods
SEC. 9.6 PREDICTOR-CORRECTOR METHODS 505 Adams-Bashforth-Moulton Method 9.6 Predictor-Corrector Methods The methods of Euler, Heun, Taylor, and Runge-Kutta are called single-step methods because they use
More informationSolution of Parabolic Partial Differential Equations
Solution of Parabolic Partial Differential Equations in Complex Geometries by a Modified Fourier Collocation Method Knut S. Eckhoff* Carl Erik Wasberg* Abstract The heat equation with Dirichlet boundary
More informationMultiscale method and pseudospectral simulations for linear viscoelastic incompressible flows
Interaction and Multiscale Mechanics, Vol. 5, No. 1 (2012) 27-40 27 Multiscale method and pseudospectral simulations for linear viscoelastic incompressible flows Ling Zhang and Jie Ouyang* Department of
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 5 Chapter 21 Numerical Differentiation PowerPoints organized by Dr. Michael R. Gustafson II, Duke University 1 All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationA posteriori error estimates applied to flow in a channel with corners
Mathematics and Computers in Simulation 61 (2003) 375 383 A posteriori error estimates applied to flow in a channel with corners Pavel Burda a,, Jaroslav Novotný b, Bedřich Sousedík a a Department of Mathematics,
More informationBASIC EXAM ADVANCED CALCULUS/LINEAR ALGEBRA
1 BASIC EXAM ADVANCED CALCULUS/LINEAR ALGEBRA This part of the Basic Exam covers topics at the undergraduate level, most of which might be encountered in courses here such as Math 233, 235, 425, 523, 545.
More informationNumerical investigation on vortex-induced motion of a pivoted cylindrical body in uniform flow
Fluid Structure Interaction VII 147 Numerical investigation on vortex-induced motion of a pivoted cylindrical body in uniform flow H. G. Sung 1, H. Baek 2, S. Hong 1 & J.-S. Choi 1 1 Maritime and Ocean
More informationNew Diagonal-Norm Summation-by-Parts Operators for the First Derivative with Increased Order of Accuracy
AIAA Aviation -6 June 5, Dallas, TX nd AIAA Computational Fluid Dynamics Conference AIAA 5-94 New Diagonal-Norm Summation-by-Parts Operators for the First Derivative with Increased Order of Accuracy David
More informationParallel domain decomposition meshless modeling of dilute convection-diffusion of species
Boundary Elements XXVII 79 Parallel domain decomposition meshless modeling of dilute convection-diffusion of species E. Divo 1,. Kassab 2 & Z. El Zahab 2 1 Department of Engineering Technology, University
More informationACCURATE NUMERICAL METHOD FOR BLASIUS PROBLEM FOR FLOW PAST A FLAT PLATE WITH MASS TRANSFER
FDS-2000 Conference, pp. 1 10 M. Sapagovas et al. (Eds) c 2000 MII ACCURATE UMERICAL METHOD FOR BLASIUS PROBLEM FOR FLOW PAST A FLAT PLATE WITH MASS TRASFER B. GAHA 1, J. J. H. MILLER 2, and G. I. SHISHKI
More informationNumerical Methods in Physics and Astrophysics
Kostas Kokkotas 2 October 20, 2014 2 http://www.tat.physik.uni-tuebingen.de/ kokkotas Kostas Kokkotas 3 TOPICS 1. Solving nonlinear equations 2. Solving linear systems of equations 3. Interpolation, approximation
More informationFour Point Gauss Quadrature Runge Kuta Method Of Order 8 For Ordinary Differential Equations
International journal of scientific and technical research in engineering (IJSTRE) www.ijstre.com Volume Issue ǁ July 206. Four Point Gauss Quadrature Runge Kuta Method Of Order 8 For Ordinary Differential
More informationFourier analysis for discontinuous Galerkin and related methods. Abstract
Fourier analysis for discontinuous Galerkin and related methods Mengping Zhang and Chi-Wang Shu Abstract In this paper we review a series of recent work on using a Fourier analysis technique to study the
More informationA THEORETICAL INTRODUCTION TO NUMERICAL ANALYSIS
A THEORETICAL INTRODUCTION TO NUMERICAL ANALYSIS Victor S. Ryaben'kii Semyon V. Tsynkov Chapman &. Hall/CRC Taylor & Francis Group Boca Raton London New York Chapman & Hall/CRC is an imprint of the Taylor
More informationComparison of Numerical Ordinary Differential Equation Solvers
Adrienne Criss Due: October, 008 Comparison of Numerical Ordinary Differential Equation Solvers Many interesting physical systems can be modeled by ordinary differential equations (ODEs). Since it is always
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 informationNumerical Methods in Physics and Astrophysics
Kostas Kokkotas 2 October 17, 2017 2 http://www.tat.physik.uni-tuebingen.de/ kokkotas Kostas Kokkotas 3 TOPICS 1. Solving nonlinear equations 2. Solving linear systems of equations 3. Interpolation, approximation
More informationA Non-Intrusive Polynomial Chaos Method For Uncertainty Propagation in CFD Simulations
An Extended Abstract submitted for the 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada January 26 Preferred Session Topic: Uncertainty quantification and stochastic methods for CFD A Non-Intrusive
More informationThird Edition. William H. Press. Raymer Chair in Computer Sciences and Integrative Biology The University of Texas at Austin. Saul A.
NUMERICAL RECIPES The Art of Scientific Computing Third Edition William H. Press Raymer Chair in Computer Sciences and Integrative Biology The University of Texas at Austin Saul A. Teukolsky Hans A. Bethe
More informationnek5000 massively parallel spectral element simulations
nek5000 massively parallel spectral element simulations PRACE Scientific Seminar HPC Boosts Science, 22th February 2011 P. Schlatter & D. S. Henningson Linné Flow Centre, KTH Mechanics Fluid flows Tornado,
More informationResearch Article A Note on the Solutions of the Van der Pol and Duffing Equations Using a Linearisation Method
Mathematical Problems in Engineering Volume 1, Article ID 693453, 1 pages doi:11155/1/693453 Research Article A Note on the Solutions of the Van der Pol and Duffing Equations Using a Linearisation Method
More informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
More informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY, USA. Dan_Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
More information16.7 Multistep, Multivalue, and Predictor-Corrector Methods
740 Chapter 16. Integration of Ordinary Differential Equations 16.7 Multistep, Multivalue, and Predictor-Corrector Methods The terms multistepand multivaluedescribe two different ways of implementing essentially
More informationVolterra integral equations solved in Fredholm form using Walsh functions
ANZIAM J. 45 (E) ppc269 C282, 24 C269 Volterra integral equations solved in Fredholm form using Walsh functions W. F. Blyth R. L. May P. Widyaningsih (Received 8 August 23; revised 6 Jan 24) Abstract Recently
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationhp-adaptive least squares spectral element method for hyperbolic partial differential equations
Journal of Computational and Applied Mathematics 215 (2008) 409 418 www.elsevier.com/locate/cam hp-adaptive least squares spectral element method for hyperbolic partial differential equations Árpád Galvão
More informationIntegration of Ordinary Differential Equations
Integration of Ordinary Differential Equations Com S 477/577 Nov 7, 00 1 Introduction The solution of differential equations is an important problem that arises in a host of areas. Many differential equations
More informationResearch Article Solution of the Porous Media Equation by a Compact Finite Difference Method
Hindawi Publishing Corporation Mathematical Problems in Engineering Volume 2009, Article ID 9254, 3 pages doi:0.55/2009/9254 Research Article Solution of the Porous Media Equation by a Compact Finite Difference
More informationThe JHU Turbulence Databases (JHTDB)
The JHU Turbulence Databases (JHTDB) TURBULENT CHANNEL FLOW DATA SET Data provenance: J. Graham 1, M. Lee 2, N. Malaya 2, R.D. Moser 2, G. Eyink 1 & C. Meneveau 1 Database ingest and Web Services: K. Kanov
More informationNUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING
NUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING C. Pozrikidis University of California, San Diego New York Oxford OXFORD UNIVERSITY PRESS 1998 CONTENTS Preface ix Pseudocode Language Commands xi 1 Numerical
More information, /'/ NI'9 ; / Proceedings of the Third International. High Order Methods. Conference on Spectral and. Houston, Texas
; /, /'/ NI'9 - Edited by Andrew V. Ilin and L. Ridgway Scott Proceedings of the Third International Conference on Spectral and High Order Methods Houston, Texas 5-9 June 1995 Houston Journal of Mathematics
More informationResearch Article Evaluation of the Capability of the Multigrid Method in Speeding Up the Convergence of Iterative Methods
International Scholarly Research Network ISRN Computational Mathematics Volume 212, Article ID 172687, 5 pages doi:1.542/212/172687 Research Article Evaluation of the Capability of the Multigrid Method
More informationNumerical Methods for Engineers. and Scientists. Applications using MATLAB. An Introduction with. Vish- Subramaniam. Third Edition. Amos Gilat.
Numerical Methods for Engineers An Introduction with and Scientists Applications using MATLAB Third Edition Amos Gilat Vish- Subramaniam Department of Mechanical Engineering The Ohio State University Wiley
More informationOptimal multilevel preconditioning of strongly anisotropic problems.part II: non-conforming FEM. p. 1/36
Optimal multilevel preconditioning of strongly anisotropic problems. Part II: non-conforming FEM. Svetozar Margenov margenov@parallel.bas.bg Institute for Parallel Processing, Bulgarian Academy of Sciences,
More information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
More informationSpectral method for the unsteady incompressible Navier Stokes equations in gauge formulation
Report no. 04/09 Spectral method for the unsteady incompressible Navier Stokes equations in gauge formulation T. W. Tee I. J. Sobey A spectral method which uses a gauge method, as opposed to a projection
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 informationOn A Comparison of Numerical Solution Methods for General Transport Equation on Cylindrical Coordinates
Appl. Math. Inf. Sci. 11 No. 2 433-439 (2017) 433 Applied Mathematics & Information Sciences An International Journal http://dx.doi.org/10.18576/amis/110211 On A Comparison of Numerical Solution Methods
More informationPoisson Equation in 2D
A Parallel Strategy Department of Mathematics and Statistics McMaster University March 31, 2010 Outline Introduction 1 Introduction Motivation Discretization Iterative Methods 2 Additive Schwarz Method
More informationChapter 1 Direct Modeling for Computational Fluid Dynamics
Chapter 1 Direct Modeling for Computational Fluid Dynamics Computational fluid dynamics (CFD) is a scientific discipline, which aims to capture fluid motion in a discretized space. The description of the
More informationarxiv: v4 [physics.comp-ph] 21 Jan 2019
A spectral/hp element MHD solver Alexander V. Proskurin,Anatoly M. Sagalakov 2 Altai State Technical University, 65638, Russian Federation, Barnaul, Lenin prospect,46, k2@list.ru 2 Altai State University,
More informationAn Overly Simplified and Brief Review of Differential Equation Solution Methods. 1. Some Common Exact Solution Methods for Differential Equations
An Overly Simplified and Brief Review of Differential Equation Solution Methods We will be dealing with initial or boundary value problems. A typical initial value problem has the form y y 0 y(0) 1 A typical
More informationOn the Non-linear Stability of Flux Reconstruction Schemes
DOI 10.1007/s10915-011-9490-6 TECHNICA NOTE On the Non-linear Stability of Flux econstruction Schemes A. Jameson P.E. Vincent P. Castonguay eceived: 9 December 010 / evised: 17 March 011 / Accepted: 14
More informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
More informationAN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION. Abstract
nd Workshop on CFD Uncertainty Analysis - Lisbon, 19th and 0th October 006 AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION Alfredo Iranzo 1, Jesús Valle, Ignacio Trejo 3, Jerónimo
More informationContents. Preface to the Third Edition (2007) Preface to the Second Edition (1992) Preface to the First Edition (1985) License and Legal Information
Contents Preface to the Third Edition (2007) Preface to the Second Edition (1992) Preface to the First Edition (1985) License and Legal Information xi xiv xvii xix 1 Preliminaries 1 1.0 Introduction.............................
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationGauge finite element method for incompressible flows
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2000; 34: 701 710 Gauge finite element method for incompressible flows Weinan E a, *,1 and Jian-Guo Liu b,2 a Courant Institute
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationTermination criteria for inexact fixed point methods
Termination criteria for inexact fixed point methods Philipp Birken 1 October 1, 2013 1 Institute of Mathematics, University of Kassel, Heinrich-Plett-Str. 40, D-34132 Kassel, Germany Department of Mathematics/Computer
More informationMEEN 673 ********************************* NONLINEAR FINITE ELEMENT ANALYSIS ********************************* J. N. Reddy
MEEN 673 ********************************* NONLINEAR FINITE ELEMENT ANALYSIS ********************************* J. N. Reddy e-mail: jnreddy@tamu.edu Texas A&M University College Station, TX 77843-3123 USA
More informationSpectral Methods and Inverse Problems
Spectral Methods and Inverse Problems Omid Khanmohamadi Department of Mathematics Florida State University Outline Outline 1 Fourier Spectral Methods Fourier Transforms Trigonometric Polynomial Interpolants
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationUNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES
UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the
More informationAN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION
AN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION Fareed Hussain Mangi*, Umair Ali Khan**, Intesab Hussain Sadhayo**, Rameez Akbar Talani***, Asif Ali Memon* ABSTRACT High order
More informationA LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES
A LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES P. HANSBO Department of Applied Mechanics, Chalmers University of Technology, S-4 96 Göteborg, Sweden E-mail: hansbo@solid.chalmers.se
More informationComputational Methods for Optoelectronics
IPT 526 Computational Methods for Optoelectronics Ray-Kuang Lee Department of Electrical Engineering and Institute of Photonics Technologies National Tsing-Hua University, Hsinchu, Taiwan e-mail: rklee@ee.nthu.edu.tw
More information