The solution of boundary element problems using symbolic computation A. Almeida, H.L. Pina Institute Superior Tecnico, Lisbon, Portugal
|
|
- Hortense Garrett
- 5 years ago
- Views:
Transcription
1 The solution of boundary element problems using symbolic computation A. Almeida, H.L. Pina Institute Superior Tecnico, Lisbon, Portugal Abstract The paper explores the potential of using Symbolic Computation in BEM problems. The work focuses on the following aspects: the development of analytic integration formulae to compute the element matrices with either regular or singular kernels, the study of how collocation node placement influences the results and the possibility of automatic FORTRAN code generation via Symbolic Computation Packages. 1 Introduction Let us consider the following boundary value problem (b.v.p.). C^(fi) (in order to be a solution in the classical sense) such that: Find u G v^(x) = o, Vxenc]R\ (i) Q is a open bounded subset of IR^, and F its Lipschitz-continuous boundary. The boundary conditions are taken as w(x) = iz(x), x 6 To (2) flu q(x)= = q(x), X l\ (3) where g(x) and %T(x) are given functions and n is the outer unit normal and r = TO u r\ (4) The previous b.v.p. can be solved by use of the Green function associated with the b.v.p. Let us introduce the following partial differential equation:
2 29 Boundary Element Technology W(x) + f(x-xi) = %'n,, (5) where u~ is the fundamental solution of Eq. (1) and 6 is the Dirac mass concentrated at point x. The fundamental solution for a two-dimensional space is given by (see [1]): «"(r) = ^log(^), (6) where r is the Euclidean distance in IR? from the field point x to the source point x;. Now we introduce the following Green's formula: Lemma 1 (Green's Formula) Let fl be a bounded open set and T its Lip shitz- continuous boundary, then the following formula holds: f (uv*v - v/n r an an Existence and uniqueness of the weak solution of Eq. (1) and Eq. (5) can easily be established. Both weak solutions are H* regular and so it is possible to apply Lemma 1 with v = u*. We obtain after some algebraic manipulation the following integral representation for u in fi, where z; H and u,u~ H*( l). Eq. (7) involves knowing the values of both u and du/dn on the boundary F. But by the problem definition (Eqs. (l)-(3)) we only know u in FQ and du/du in FI. In order to determine the value of u in FI and du/dn. in FQ instead of an integral representation in fi, we need now an integral representation on the boundary F. By a limiting process the interior point Xj can be taken to the boundary thus obtaining the integral representation for u on the boundary F. = - /r i^l> + /r where the fundamental solution u* has been substituted by Eq. (6) and C is a constant that is given by the fraction of solid angle subtended by F at x. For a regular boundary (that is piecewise smooth C* boundary) C = 1/2. We are now able to compute the values of u in any point of fi ( 1 is the closure of Q) by using Eq. (8) and/or Eq. (7).
3 Boundary Element Technology Problem discretization We can discretize the previous problem by replacing the unknown functions u and du/dn in the boundary (in TO and I\ respectively) by a finite element approximation on the boundary. Adopting a master element description of the discretized problem we will use the notation u and q for for the approximations of u and q = du/dn respectively in the master boundary finite element F. The elementary shape functions will be denoted by ^, where i = 1... N<, is the node number in F. As usual the master element coordinate is denoted by and 1 <(< 1. Thus we have: ( = <Mfc, (9) 1 = 1 «( = *4 (1) 1=1 where U{ and qi are the nodal values of u and g, respectively, in the master element. The discrete equation corresponding to Eq. (8) is thus: Z / %k, ir ^ = E / %^(xi j=l^3 OH y.' j=l ^J where x' is a node point and XQ is the collocation point; u^. and q^ are the restrictions of the approximation of u and q to the boundary element FJ; TV is the number of nodes. We can write Eq. (11) in matrix form: where ds. (12) (13) and N. i = E L (14) 3=1 *> The coordinate transformation between a generic boundary element Fj and the master element F is given by the parametric transformation: M X = t=l M (15) y =
4 292 Boundary Element Technology where (x;,2/{) is element Tj nodal coordinates and M is the number of nodes of transformation Tg. The basic assumptions that we will use throughout this paper are the following: Only constant and linear approximations of u and q in each element will be analyzed. Boundaries are C* piecewise smooth. Only symmetrical node placement w.r.t. = will be considered. In what concerns the continuity across elements we will deal only discontinuous approximations. As a direct consequence of the last assumption we will use different shape functions for geometry interpolation and u and q approximation. By Eq. (15) the coordinate transformation is: x = - -- f +, 2/2 + where Eq. (16) was obtained by adopting the Lagrange polynomials of first degree in the master element F for geometry interpolation, i.e.: (l-, (17) i(l +, (18) The subscript G indicates that these are shape functions for geometry interpolation. This assumption may seen, at first sight, strange. It was made in order to enable us to study the influence of the element collocation nodes position. We will return to this subject later on. 3 Formulae derivation using Maple We will enter now the specific subject of this paper which focus the use of a computer algebra system, in this case Maple [2, 3, 4], in the solution of problems by the Boundary Element Method (BEM). The first task for which we used Maple was the derivation of closed form expressions for the integrals that appear in Eqs. (13)-(14).
5 Boundary Element Technology 293 In the case of constant approximation of u and q in each element and by Eqs. (6), (13) we have the H and Q matrices entries: y = a, - E / ''fo where SH = - j- /' log(r)v'(*2-ai)* + (yj-»i)'#, (2) **" j=i-'-i r = (*' - ^ - 2i±Si)' + (,' - 2L^ - ^)>. (21) In what concerns linear approximation of u and q in each element Eqs. (16)-(18) still hold but the matrices H and Q entries are: ~ ~, \ > (22) (23) where r is given by Eq. (21). The shape functions are given by: ^ = o, v,,(( + 4-1), (24) 2 + &o + (%(, V"2 = ^ U + da ~ 1), (25) 2 + da 4- at where da and <f& are the distance from the element nodes (for function approximation) to the extreme nodes (for geometry interpolation ) as illustrated in Figure 1. For da. = db = the element nodes for geometry interpolation coincide with element nodes for function approximation. Only in this case continuity across elements is obtained. Note that for the constant element the shape function is ^i = 1 and the element node for function approximation is placed in the element center, corresponding to =. The integrals that appear in the expressions for computing %., and Qij are regular when the collocation point lies outside the domain of integration. Maple's int function can handle them quite easily. However if the collocation point x' lies inside the domain of integration the integrals for 7Y%j and Qij become singular.
6 294 Boundary Element Technology \ /2 da db node 1 node 2 Figure 1: Linear element and node location Obtaining closed form expressions for such integrals using Maple demands certain precautions in order to be able to compute such expressions successfully. We will compute the following integrals: Ik = (26) where fc =,1 and r is given by Eq. (21). Note that the nodal coordinates are given by: where a, b [,1] are such that: Xi = ax i + bxi Ui = &2/i + fy/2 (27) -6) (28) is the master element f length. In this case = 1 - (-1) = 2. Let us recall that d^ and d^ are the distance from the nearest element node to the left and right extremes ( = 1 and = 1) of the element respectively (see Eqs. (24) and (25)). Note that the integrals given by Eq. (26) are exactly the integrals contributing to Hij (Eqs. (2) and (21)). When x; belongs to the domain of integration, the integrals that appear in the expressions for computing Ti,ij are due to orthogonality between the distance r and the outer normal n. This very interesting result holds true only for straight elements. Therefore the computation of singular integrals reduces to computing closed form expressions for Eq. (26). Introducing the expression relative to Eq. (26) for k = in Maple and using the int function with -1 < < 1 we get a rather lengthy expression as
7 Boundary Element Technology 295 a result. The expression by itself is of little interest. Instead we will explain some details of the procedure followed in its deduction. In the expression that Maple returns as a result of computing IQ two terms involving the arctan function appear. Those terms have a denominator (after factorization) involving the term 1 (a-t-6). Note that according to our assumptions both element nodes are inside the element (i.e. in [ 1, 1]) and we must have a = 1 6. If a susbstitution of this result is made in the expression of the integral using Maple subs function in the arctan terms, a division by zero error will be issued. This is not what we expected. We would like Maple to give us the value of the arctan terms that will be %/2 or?r/2 according to the sign of the numerator of those terms. However as Maple has a fully recursive evaluation mode for this kind of expressions (only procedures and tables are evaluated by name), Maple will first try to evaluate the arctan function argument which gives a division by zero, and thus signaling an error. The solution is to "isolate" the arctan terms from the remainder of the expression and then extracting the numerator and denominator of the arctan terms argument handling them separately. By doing so we are able to sucessfully compute the wanted expressions. The arctan terms will be either 7T/2 or?r/2 according to the sign of the factor 6 which appears in the denominator of one term and 6 1 for the other. Because b belongs to [, 1] we have: 1 < b<, -1 < 6-1 <, concluding that the arctan terms cancel each other for b G (, 1) because they have opposite signs and obviously ir/2 -f Tr/2 =. Accordingly we can now obtain the formula for IQ with b (, 1). - log 6-1 ) + 21og Zi)2 + (7/2 - y,)2))(z2-zi)2 + (2/2-3/1)'. (29) Putting a = 6 = 1/2 we obtain the formula for the integral contributing to Ga for the constant element which can be found in any introductory text to the BEM (e.g see [1]). /o = -(log((z2-zi)' + (3/2-3/i)')-2(log(2)-l))\(z2-21)2 + (3/2 - For /i (putting k = 1 in Eq. (26)) following a similar procedure we have -log 6- l\)b(b- 1) -26+ l)\(z:-zi)2 + (3,2-3,1)2, (31)
8 296 Boundary Element Technology Combining Eq. (29) and Eq. (31) we can obtain closed form expressions for linear elements with shape functions as those given by Eqs. (24) and (25). These expressions are valid for an arbitrary node placement in the element. The relations between the parameters a, 6 and d^ 4 are established by Eqs. (28). This enables us to study node placement influence in BEM solution behaviour. 4 Comparison between Maple's and classical BEM code solutions The model problem that we have solved consists on the following: V*u = in, (32) with n = [,5] x [,5] and boundary conditions and The solution to this b.v.p. is given by: fjf(«,) = (*,5) =, (33) %(,y) = A, %(5,3/) = B. (34) N B-A u(x,y) = - z + A. (35) o In this case we have taken A 1 and B = 2 thus u(x,y) 2x + 1. This extremely simple problem enables us to thoroughly analyze the several aspects concerning the influence of numerical integration in BEM solution behavior. Let us now compare the Q and H matrices obtained by use of Maple's int function (analytical integration) with those obtained by a two point Gauss-Legendre rule. We have solved the previous b.v.p. using a 16 element mesh with both constant and linear elements. Below we present a few typical entries of the Q and H matrices. The 7i matrix for constant elements (analytical integration)
9 Boundary Element Technology 297 The H matrix for constant elements with Gauss-Legendre quadrature The Q matrix for constant elements using Eq. (3) (analytical integration) The Q matrix for constant elements with Gauss-Legendre quadrature For linear elements the H matrices the results are similar. The diagonal terms are for the Gauss-Legendre quadrature and for analytical integration. The Q matrix for linear elements using Eq. (31) (analytical integration) The Q matrix for linear elements with Gauss-Legendre quadrature A comparison between the entries of the Q matrix for numerical and analytical integration reveals sensible differences in the entries corresponding to the singular integrals. By using numerical integration the singular behaviour of the integrand cannot be captured whereas through the use of analytical integration the integrand singular behaviour is properly taken in to account.
10 298 Boundary Element Technology 5 Node placement influence in BEM solution behaviour When using analytical integration the node placement influence in solution behaviour can be analyzed. Although the problem presented here is extremely simple and consequently the node placement influence in the solution was less obvious than would be in a more elaborated problem, some variation in the solution was observed when the node placement is changed. Here we have considered only two positions for element node placement. This positions being = ±V3/3 (the Gauss-Legendre points) and = ±1/2. These positions are usually referred in the BEM literature as being the most appropriate for node placement, since they tend to contribute for a better performance of the BEM. The results we have obtained are presented in the table below. Nodes X y u f = W3/3 9= J %, u , f = ±1/2 <?* LI As we can observe, there are sizable diferences in the results particularly in the flux (derivatives) values. <7v 6 A brief comment on the Maple program The Maple program that was developed has basically seven different phases. The program layout is described in Figure 2. In all the phases full use of Maple's symbolic computation capability was made. Mesh generation was done by procedures that were developed specifically for the program. These procedures generate 3 type of curves: archs; parabolas and straight lines. The user specifies the number of points in the curve. The system of equations Tiu = gq was solved using Maple's built in function solve. This approach was preferred over a use of the linsolve function. In order to use the linsolve function the system should be in the form Ax = b. A somewhat complicated algorithm is needed for transforming Tiu = Qq in Ax = b. Using solve we can collect the unknown variables in to a set and give any set of equations (which is consistent) to solve that
11 Boundary Element Technology 299 Problem Data: V* = I Elementary Q Ji matrices I Mesh data reading] Q and H matrices assembly [Essential and natural B.C. imposition I Equation System solving I Computation of u and q in Figure 2: Program layout envolves those unknown variables. The disadvantage of following such an approach is that results become dependent on Maple's solving algorithm. However in the several examples that we have solved using this program the results obtained were always satisfactory. The computation of the diagonal terms of matrix Ji was done by considering a constant potential (or rigid body motions as in Mechanics). In such case the sum of all the off diagonal row elements is. Thus C is given by: N C = E%j. (36) The computation of the integrals that contribute to the Q and H matrices was done using the placeholder Int (for an explanation of the several Maple placeholders see [3]). Int is not an ordinary Maple function. Instead of trying to compute a value corresponding to the function invocation, a placeholder simply returns a string indicating the type of operation we would like to perform in a given expression. We can say that a placeholder "freezes" the evaluation of a Maple expression. The evaluation can be delayed until is strictly necessary. Maple's evaluator recognize placeholders and knows how to deal with such expressions in order to obtain (if possible) a numeric or algebraic expression. The placeholder concept enables Maple to simulate a stream, in the sense that when we use a placeholder an unevaluated expression and a way to evaluate it is being specified (albeit in an implicit fashion). In [3] a placeholder is designated also as a function inert form. In this paper we have prefered the first designation since, in our opinion, it describes more
12 3 Boundary Element Technology accurately how Maple deals with such objects. In the program that we have implemented the Q and *H matrices are assembled with their entries being placeholders for the int function. It is only after assembling the matrices that the numerical values of the entries are computed by using evalf. This approach avoids Maple from signaling divisions by zero while computing the integrals. Note that have we used int instead of Int and several division by zero errors would have been signaled when computing the integrals. A FORTRAN (C) output of Maple results can easily be obtained through the use of the f ortran (C) builtin function. This function invokes a parser to translate Maple algebraic expressions in to FORTRAN (C) code (see [3]). 7 Conclusions We have shown that the use of Symbolic Computation Packages can present advantages in the deduction of analytic (exact) integration formulae for computation of BEM matrices for the case of constant and linear shape functions. Uncertainties regarding numeric (approximate) integration are thus avoided and eventually more efficient computer programs are produced. Besides, the flexibility provided by these packages allow for easier code development and experimentation before one commits itself to a particular algorithm for production software as was exemplified with the collocation node localization study. Undoubtly we shall see in the near future mathematical software with fully integrated numerical and symbolic capabilities makeing knowledge of Symbolic Computation Packages indispensable. Although it is not yet possible to generate FORTRAN or C code directly from a symbolic manipulation language (only algebraic expressions are convertible to FORTRAN or C), compiled code is now obtainable (see [7]) making the symbolic approach to problem solving "almost" as efficient as a numerical approach coded in FOR- TRAN or C. Acknowledgements The present work was partially supported by JNICT, Junta Nacional de Investigate Cientifica, Programa CIENCIA, BM/2975/92-IB, and by Institute Superior Tecnico, Lisbon. References [1] C. A. Brebbia. The Boundary Element Method for Engineers. Pentech Press, 1978.
13 Boundary Element Technology 31 [2] B. W. Char et al. Maple Language Reference Manual Springer-Verlag, [3] B. W. Char et al. Maple V Library Reference Manual. Springer-Verlag, [4] B. W. Char et al. First Leaves: A Tutorial Introduction to Maple V. Springer-Verlag, [5] K.. Geddes. Numerical integration using symbolic analysis. Maple Technical Newsletter, (6):8-17, [6] J.-L. Lions and E. Magenes. Non-Homogeneous Boundary Value Problems and Applications, volume 1. Springer-Verlag, [7] Robert S. Tutor, editor. AXIOM User Guide: A System for doing Computer Mathematics. Numerical Algorithms Group, Inc., 1991.
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 36
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 36 In last class, we have derived element equations for two d elasticity problems
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationThe goal is the definition of points with numbers and primitives with equations or functions. The definition of points with numbers requires a
The goal is the definition of points with numbers and primitives with equations or functions. The definition of points with numbers requires a coordinate system and then the measuring of the point with
More informationMODELING MIXED BOUNDARY PROBLEMS WITH THE COMPLEX VARIABLE BOUNDARY ELEMENT METHOD (CVBEM) USING MATLAB AND MATHEMATICA
A. N. Johnson et al., Int. J. Comp. Meth. and Exp. Meas., Vol. 3, No. 3 (2015) 269 278 MODELING MIXED BOUNDARY PROBLEMS WITH THE COMPLEX VARIABLE BOUNDARY ELEMENT METHOD (CVBEM) USING MATLAB AND MATHEMATICA
More informationAMS527: Numerical Analysis II
AMS527: Numerical Analysis II A Brief Overview of Finite Element Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 1 / 25 Overview Basic concepts Mathematical
More informationLecture 3.2 Methods for Structured Mesh Generation
Lecture 3.2 Methods for Structured Mesh Generation 1 There are several methods to develop the structured meshes: Algebraic methods, Interpolation methods, and methods based on solving partial differential
More informationFinite Difference Calculus
Chapter 2 Finite Difference Calculus In this chapter we review the calculus of finite differences. The topic is classic and covered in many places. The Taylor series is fundamental to most analysis. A
More informationSpline Methods Draft. Tom Lyche and Knut Mørken
Spline Methods Draft Tom Lyche and Knut Mørken 24th May 2002 2 Contents 1 Splines and B-splines an introduction 3 1.1 Convex combinations and convex hulls..................... 3 1.1.1 Stable computations...........................
More informationMA 323 Geometric Modelling Course Notes: Day 21 Three Dimensional Bezier Curves, Projections and Rational Bezier Curves
MA 323 Geometric Modelling Course Notes: Day 21 Three Dimensional Bezier Curves, Projections and Rational Bezier Curves David L. Finn Over the next few days, we will be looking at extensions of Bezier
More informationParameterization. Michael S. Floater. November 10, 2011
Parameterization Michael S. Floater November 10, 2011 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to generate from point
More informationSpline Methods Draft. Tom Lyche and Knut Mørken
Spline Methods Draft Tom Lyche and Knut Mørken January 5, 2005 2 Contents 1 Splines and B-splines an Introduction 3 1.1 Convex combinations and convex hulls.................... 3 1.1.1 Stable computations...........................
More informationTIME 2014 Technology in Mathematics Education July 1 st -5 th 2014, Krems, Austria
TIME 2014 Technology in Mathematics Education July 1 st -5 th 2014, Krems, Austria Overview Introduction Using a 2D Plot Window in a CAS Perspective Plotting a circle and implicit differentiation Helping
More informationContents. I The Basic Framework for Stationary Problems 1
page v Preface xiii I The Basic Framework for Stationary Problems 1 1 Some model PDEs 3 1.1 Laplace s equation; elliptic BVPs... 3 1.1.1 Physical experiments modeled by Laplace s equation... 5 1.2 Other
More informationCS 450 Numerical Analysis. Chapter 7: Interpolation
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 informationQuasilinear First-Order PDEs
MODULE 2: FIRST-ORDER PARTIAL DIFFERENTIAL EQUATIONS 16 Lecture 3 Quasilinear First-Order PDEs A first order quasilinear PDE is of the form a(x, y, z) + b(x, y, z) x y = c(x, y, z). (1) Such equations
More informationSpline Methods Draft. Tom Lyche and Knut Mørken. Department of Informatics Centre of Mathematics for Applications University of Oslo
Spline Methods Draft Tom Lyche and Knut Mørken Department of Informatics Centre of Mathematics for Applications University of Oslo January 27, 2006 Contents 1 Splines and B-splines an Introduction 1 1.1
More informationDiscrete Cubic Interpolatory Splines
Publ RIMS, Kyoto Univ. 28 (1992), 825-832 Discrete Cubic Interpolatory Splines By Manjulata SHRIVASTAVA* Abstract In the present paper, existence, uniqueness and convergence properties of a discrete cubic
More informationTHE MORTAR FINITE ELEMENT METHOD IN 2D: IMPLEMENTATION IN MATLAB
THE MORTAR FINITE ELEMENT METHOD IN D: IMPLEMENTATION IN MATLAB J. Daněk, H. Kutáková Department of Mathematics, University of West Bohemia, Pilsen MECAS ESI s.r.o., Pilsen Abstract The paper is focused
More informationMath background. 2D Geometric Transformations. Implicit representations. Explicit representations. Read: CS 4620 Lecture 6
Math background 2D Geometric Transformations CS 4620 Lecture 6 Read: Chapter 2: Miscellaneous Math Chapter 5: Linear Algebra Notation for sets, functions, mappings Linear transformations Matrices Matrix-vector
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 24
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 24 So in today s class, we will look at quadrilateral elements; and we will
More informationVector Algebra Transformations. Lecture 4
Vector Algebra Transformations Lecture 4 Cornell CS4620 Fall 2008 Lecture 4 2008 Steve Marschner 1 Geometry A part of mathematics concerned with questions of size, shape, and relative positions of figures
More informationLecture 2.2 Cubic Splines
Lecture. Cubic Splines Cubic Spline The equation for a single parametric cubic spline segment is given by 4 i t Bit t t t i (..) where t and t are the parameter values at the beginning and end of the segment.
More informationFinite Element Method. Chapter 7. Practical considerations in FEM modeling
Finite Element Method Chapter 7 Practical considerations in FEM modeling Finite Element Modeling General Consideration The following are some of the difficult tasks (or decisions) that face the engineer
More informationCOMP30019 Graphics and Interaction Transformation geometry and homogeneous coordinates
COMP30019 Graphics and Interaction Transformation geometry and homogeneous coordinates Department of Computer Science and Software Engineering The Lecture outline Introduction Vectors and matrices Translation
More informationSolve Non-Linear Parabolic Partial Differential Equation by Spline Collocation Method
Solve Non-Linear Parabolic Partial Differential Equation by Spline Collocation Method P.B. Choksi 1 and A.K. Pathak 2 1 Research Scholar, Rai University,Ahemdabad. Email:pinalchoksey@gmail.com 2 H.O.D.
More informationChapter 7 Practical Considerations in Modeling. Chapter 7 Practical Considerations in Modeling
CIVL 7/8117 1/43 Chapter 7 Learning Objectives To present concepts that should be considered when modeling for a situation by the finite element method, such as aspect ratio, symmetry, natural subdivisions,
More informationA Connection between Network Coding and. Convolutional Codes
A Connection between Network Coding and 1 Convolutional Codes Christina Fragouli, Emina Soljanin christina.fragouli@epfl.ch, emina@lucent.com Abstract The min-cut, max-flow theorem states that a source
More information30. Constrained Optimization
30. Constrained Optimization The graph of z = f(x, y) is represented by a surface in R 3. Normally, x and y are chosen independently of one another so that one may roam over the entire surface of f (within
More informationSecond-order shape optimization of a steel bridge
Computer Aided Optimum Design of Structures 67 Second-order shape optimization of a steel bridge A.F.M. Azevedo, A. Adao da Fonseca Faculty of Engineering, University of Porto, Porto, Portugal Email: alvaro@fe.up.pt,
More informationThe DRM-MD integral equation method for the numerical solution of convection-diffusion
The DRM-MD integral equation method for the numerical solution of convection-diffusion equation V. Popov & H. Power Wessex Institute of Technology, Ashurst Lodge, Ashurst, Abstract This work presents a
More informationInterpolation and Splines
Interpolation and Splines Anna Gryboś October 23, 27 1 Problem setting Many of physical phenomenona are described by the functions that we don t know exactly. Often we can calculate or measure the values
More informationADAPTIVE APPROACH IN NONLINEAR CURVE DESIGN PROBLEM. Simo Virtanen Rakenteiden Mekaniikka, Vol. 30 Nro 1, 1997, s
ADAPTIVE APPROACH IN NONLINEAR CURVE DESIGN PROBLEM Simo Virtanen Rakenteiden Mekaniikka, Vol. 30 Nro 1, 1997, s. 14-24 ABSTRACT In recent years considerable interest has been shown in the development
More informationMaple Quick Start. Maplesoft, a division of Waterloo Maple Inc.
Maple Quick Start Maplesoft, a division of Waterloo Maple Inc. This tutorial is designed to help you become familiar with the Maple environment and teach you the few fundamental concepts and tools you
More informationCS231A Course Notes 4: Stereo Systems and Structure from Motion
CS231A Course Notes 4: Stereo Systems and Structure from Motion Kenji Hata and Silvio Savarese 1 Introduction In the previous notes, we covered how adding additional viewpoints of a scene can greatly enhance
More informationChapter 3. Quadric hypersurfaces. 3.1 Quadric hypersurfaces Denition.
Chapter 3 Quadric hypersurfaces 3.1 Quadric hypersurfaces. 3.1.1 Denition. Denition 1. In an n-dimensional ane space A; given an ane frame fo;! e i g: A quadric hypersurface in A is a set S consisting
More informationGraphics and Interaction Transformation geometry and homogeneous coordinates
433-324 Graphics and Interaction Transformation geometry and homogeneous coordinates Department of Computer Science and Software Engineering The Lecture outline Introduction Vectors and matrices Translation
More informationMATH 2400, Analytic Geometry and Calculus 3
MATH 2400, Analytic Geometry and Calculus 3 List of important Definitions and Theorems 1 Foundations Definition 1. By a function f one understands a mathematical object consisting of (i) a set X, called
More information2D/3D Geometric Transformations and Scene Graphs
2D/3D Geometric Transformations and Scene Graphs Week 4 Acknowledgement: The course slides are adapted from the slides prepared by Steve Marschner of Cornell University 1 A little quick math background
More informationDriven Cavity Example
BMAppendixI.qxd 11/14/12 6:55 PM Page I-1 I CFD Driven Cavity Example I.1 Problem One of the classic benchmarks in CFD is the driven cavity problem. Consider steady, incompressible, viscous flow in a square
More informationMatrices and Systems of Linear Equations
Chapter The variable x has now been eliminated from the first and third equations. Next, we eliminate x3 from the first and second equations and leave x3, with coefficient, in the third equation: System:
More informationA Hybrid Magnetic Field Solver Using a Combined Finite Element/Boundary Element Field Solver
A Hybrid Magnetic Field Solver Using a Combined Finite Element/Boundary Element Field Solver Abstract - The dominant method to solve magnetic field problems is the finite element method. It has been used
More informationNon-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, Politecnico di Milano, February 3, 2017, Lesson 1
Non-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, attilio.frangi@polimi.it Politecnico di Milano, February 3, 2017, Lesson 1 1 Politecnico di Milano, February 3, 2017, Lesson 1 2 Outline
More informationLecture notes on the simplex method September We will present an algorithm to solve linear programs of the form. maximize.
Cornell University, Fall 2017 CS 6820: Algorithms Lecture notes on the simplex method September 2017 1 The Simplex Method We will present an algorithm to solve linear programs of the form maximize subject
More informationFully discrete Finite Element Approximations of Semilinear Parabolic Equations in a Nonconvex Polygon
Fully discrete Finite Element Approximations of Semilinear Parabolic Equations in a Nonconvex Polygon Tamal Pramanick 1,a) 1 Department of Mathematics, Indian Institute of Technology Guwahati, Guwahati
More informationA Random Variable Shape Parameter Strategy for Radial Basis Function Approximation Methods
A Random Variable Shape Parameter Strategy for Radial Basis Function Approximation Methods Scott A. Sarra, Derek Sturgill Marshall University, Department of Mathematics, One John Marshall Drive, Huntington
More informationMathematics Curriculum
6 G R A D E Mathematics Curriculum GRADE 6 5 Table of Contents 1... 1 Topic A: Area of Triangles, Quadrilaterals, and Polygons (6.G.A.1)... 11 Lesson 1: The Area of Parallelograms Through Rectangle Facts...
More informationLet denote the number of partitions of with at most parts each less than or equal to. By comparing the definitions of and it is clear that ( ) ( )
Calculating exact values of without using recurrence relations This note describes an algorithm for calculating exact values of, the number of partitions of into distinct positive integers each less than
More informationA spectral boundary element method
Boundary Elements XXVII 165 A spectral boundary element method A. Calaon, R. Adey & J. Baynham Wessex Institute of Technology, Southampton, UK Abstract The Boundary Element Method (BEM) is not local and
More informationElectromagnetic Field Numerical Modeling using BEM2D
Volume 56, Number 5, 2015 197 Electromagnetic Field Numerical Modeling using BEM2D Adina Răcășan, Claudia Păcurar, Călin Munteanu, Vasile Țopa, Claudia Constantinescu, Lorand Szabo and Marius Dodea Faculty
More informationMaple as an Instructional Tool
Maple as an Instructional Tool Terence A. Weigel 1 Abstract Computer algebra software such as Maple is an important component of the engineer s toolkit, much as are Matlib, MathCAD and Excel. However,
More informationAlgorithms and Data Structures
Charles A. Wuethrich Bauhaus-University Weimar - CogVis/MMC June 22, 2017 1/51 Introduction Matrix based Transitive hull All shortest paths Gaussian elimination Random numbers Interpolation and Approximation
More informationAlmost Curvature Continuous Fitting of B-Spline Surfaces
Journal for Geometry and Graphics Volume 2 (1998), No. 1, 33 43 Almost Curvature Continuous Fitting of B-Spline Surfaces Márta Szilvási-Nagy Department of Geometry, Mathematical Institute, Technical University
More informationLimits and Their Properties. Copyright Cengage Learning. All rights reserved.
1 Limits and Their Properties Copyright Cengage Learning. All rights reserved. 1.1 A Preview of Calculus Copyright Cengage Learning. All rights reserved. What Is Calculus? 3 Calculus Calculus is the mathematics
More informationMAT175 Overview and Sample Problems
MAT175 Overview and Sample Problems The course begins with a quick review/overview of one-variable integration including the Fundamental Theorem of Calculus, u-substitutions, integration by parts, and
More informationTHE COMPUTER MODELLING OF GLUING FLAT IMAGES ALGORITHMS. Alekseí Yu. Chekunov. 1. Introduction
MATEMATIQKI VESNIK Corrected proof Available online 01.10.2016 originalni nauqni rad research paper THE COMPUTER MODELLING OF GLUING FLAT IMAGES ALGORITHMS Alekseí Yu. Chekunov Abstract. In this paper
More informationJournal of Engineering Research and Studies E-ISSN
Journal of Engineering Research and Studies E-ISS 0976-79 Research Article SPECTRAL SOLUTIO OF STEADY STATE CODUCTIO I ARBITRARY QUADRILATERAL DOMAIS Alavani Chitra R 1*, Joshi Pallavi A 1, S Pavitran
More informationTransactions on Modelling and Simulation vol 20, 1998 WIT Press, ISSN X
r-adaptive boundary element method Eisuke Kita, Kenichi Higuchi & Norio Kamiya Department of Mechano-Informatics and Systems, Nagoya University, Nagoya 464-01, Japan Email: kita@mech.nagoya-u.ac.jp Abstract
More informationParameterization of triangular meshes
Parameterization of triangular meshes Michael S. Floater November 10, 2009 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to
More informationSung-Eui Yoon ( 윤성의 )
CS480: Computer Graphics Curves and Surfaces Sung-Eui Yoon ( 윤성의 ) Course URL: http://jupiter.kaist.ac.kr/~sungeui/cg Today s Topics Surface representations Smooth curves Subdivision 2 Smooth Curves and
More information10. Cartesian Trajectory Planning for Robot Manipulators
V. Kumar 0. Cartesian rajectory Planning for obot Manipulators 0.. Introduction Given a starting end effector position and orientation and a goal position and orientation we want to generate a smooth trajectory
More informationLacunary Interpolation Using Quartic B-Spline
General Letters in Mathematic, Vol. 2, No. 3, June 2017, pp. 129-137 e-issn 2519-9277, p-issn 2519-9269 Available online at http:\\ www.refaad.com Lacunary Interpolation Using Quartic B-Spline 1 Karwan
More informationIsogeometric Collocation Method
Chair for Computational Analysis of Technical Systems Faculty of Mechanical Engineering, RWTH Aachen University Isogeometric Collocation Method Seminararbeit By Marko Blatzheim Supervisors: Dr. Stefanie
More information13.1. Functions of Several Variables. Introduction to Functions of Several Variables. Functions of Several Variables. Objectives. Example 1 Solution
13 Functions of Several Variables 13.1 Introduction to Functions of Several Variables Copyright Cengage Learning. All rights reserved. Copyright Cengage Learning. All rights reserved. Objectives Understand
More informationA mesh refinement technique for the boundary element method based on local error analysis
A mesh refinement technique for the boundary element method based on local error analysis J. J. Rodriguez
More informationCOMPUTER AIDED ENGINEERING DESIGN (BFF2612)
COMPUTER AIDED ENGINEERING DESIGN (BFF2612) BASIC MATHEMATICAL CONCEPTS IN CAED by Dr. Mohd Nizar Mhd Razali Faculty of Manufacturing Engineering mnizar@ump.edu.my COORDINATE SYSTEM y+ y+ z+ z+ x+ RIGHT
More informationCCSSM Curriculum Analysis Project Tool 1 Interpreting Functions in Grades 9-12
Tool 1: Standards for Mathematical ent: Interpreting Functions CCSSM Curriculum Analysis Project Tool 1 Interpreting Functions in Grades 9-12 Name of Reviewer School/District Date Name of Curriculum Materials:
More informationMathematica CalcCenter
Mathematica CalcCenter Basic features Wolfram Mathematica CalcCenter is based on Mathematica Professional and it is primarily designed for technical calculations. Information about this product can be
More informationChallenge Problem 5 - The Solution Dynamic Characteristics of a Truss Structure
Challenge Problem 5 - The Solution Dynamic Characteristics of a Truss Structure In the final year of his engineering degree course a student was introduced to finite element analysis and conducted an assessment
More informationGuidelines for proper use of Plate elements
Guidelines for proper use of Plate elements In structural analysis using finite element method, the analysis model is created by dividing the entire structure into finite elements. This procedure is known
More informationDOWNLOAD PDF BIG IDEAS MATH VERTICAL SHRINK OF A PARABOLA
Chapter 1 : BioMath: Transformation of Graphs Use the results in part (a) to identify the vertex of the parabola. c. Find a vertical line on your graph paper so that when you fold the paper, the left portion
More informationChapter 15 Introduction to Linear Programming
Chapter 15 Introduction to Linear Programming An Introduction to Optimization Spring, 2015 Wei-Ta Chu 1 Brief History of Linear Programming The goal of linear programming is to determine the values of
More informationA substructure based parallel dynamic solution of large systems on homogeneous PC clusters
CHALLENGE JOURNAL OF STRUCTURAL MECHANICS 1 (4) (2015) 156 160 A substructure based parallel dynamic solution of large systems on homogeneous PC clusters Semih Özmen, Tunç Bahçecioğlu, Özgür Kurç * Department
More informationChapter Seventeen. Gauss and Green. We shall do this by computing the surface integral over each of the six sides of B and adding the results.
Chapter Seventeen Gauss and Green 7 Gauss's Theorem Let B be the bo, or rectangular parallelepiped, given by B {(, y, z):, y y y, z z z } 0 0 0 ; and let S be the surface of B with the orientation that
More informationAn ε-uniform Initial Value Technique For Convection-Diffusion Singularly Perturbed Problems
An -Uniform Initial Value Technique For Convection-Diffusion Singularly Perturbed Problems RKBawa and Vinod Kumar Abstract In this paper, we have proposed an -uniform initial value technique for singularly
More informationGeneralized barycentric coordinates
Generalized barycentric coordinates Michael S. Floater August 20, 2012 In this lecture, we review the definitions and properties of barycentric coordinates on triangles, and study generalizations to convex,
More informationSection 17.7: Surface Integrals. 1 Objectives. 2 Assignments. 3 Maple Commands. 4 Lecture. 4.1 Riemann definition
ection 17.7: urface Integrals 1 Objectives 1. Compute surface integrals of function of three variables. Assignments 1. Read ection 17.7. Problems: 5,7,11,1 3. Challenge: 17,3 4. Read ection 17.4 3 Maple
More informationTHE COMPUTER MODELLING OF GLUING FLAT IMAGES ALGORITHMS. Alekseí Yu. Chekunov. 1. Introduction
MATEMATIČKI VESNIK MATEMATIQKI VESNIK 69, 1 (2017), 12 22 March 2017 research paper originalni nauqni rad THE COMPUTER MODELLING OF GLUING FLAT IMAGES ALGORITHMS Alekseí Yu. Chekunov Abstract. In this
More informationCalculus I Review Handout 1.3 Introduction to Calculus - Limits. by Kevin M. Chevalier
Calculus I Review Handout 1.3 Introduction to Calculus - Limits by Kevin M. Chevalier We are now going to dive into Calculus I as we take a look at the it process. While precalculus covered more static
More informationTransactions on Modelling and Simulation vol 20, 1998 WIT Press, ISSN X
Parallel indirect multipole BEM analysis of Stokes flow in a multiply connected domain M.S. Ingber*, A.A. Mammoli* & J.S. Warsa* "Department of Mechanical Engineering, University of New Mexico, Albuquerque,
More informationHands-On Standards Deluxe Grades: 7, 8 States: California Content Standards
Hands-On Standards Deluxe Grades: 7, 8 States: Hands-On Standards Deluxe Edition Kit, Grades 5-6: Algebra Summary: This resource guide meets 5 and 6 math curriculum standards by matching activities to
More information6 Mathematics Curriculum
New York State Common Core 6 Mathematics Curriculum GRADE GRADE 6 MODULE 5 Table of Contents 1 Area, Surface Area, and Volume Problems... 3 Topic A: Area of Triangles, Quadrilaterals, and Polygons (6.G.A.1)...
More informationLecture 25: Bezier Subdivision. And he took unto him all these, and divided them in the midst, and laid each piece one against another: Genesis 15:10
Lecture 25: Bezier Subdivision And he took unto him all these, and divided them in the midst, and laid each piece one against another: Genesis 15:10 1. Divide and Conquer If we are going to build useful
More informationFitting Uncertain Data with NURBS
Fitting Uncertain Data with NURBS Wolfgang Heidrich, Richard Bartels, George Labahn Abstract. Fitting of uncertain data, that is, fitting of data points that are subject to some error, has important applications
More informationHowever, m pq is just an approximation of M pq. As it was pointed out by Lin [2], more precise approximation can be obtained by exact integration of t
FAST CALCULATION OF GEOMETRIC MOMENTS OF BINARY IMAGES Jan Flusser Institute of Information Theory and Automation Academy of Sciences of the Czech Republic Pod vodarenskou vez 4, 82 08 Prague 8, Czech
More informationA TESSELLATION FOR ALGEBRAIC SURFACES IN CP 3
A TESSELLATION FOR ALGEBRAIC SURFACES IN CP 3 ANDREW J. HANSON AND JI-PING SHA In this paper we present a systematic and explicit algorithm for tessellating the algebraic surfaces (real 4-manifolds) F
More informationChapter 13. Boundary Value Problems for Partial Differential Equations* Linz 2002/ page
Chapter 13 Boundary Value Problems for Partial Differential Equations* E lliptic equations constitute the third category of partial differential equations. As a prototype, we take the Poisson equation
More informationPaul's Online Math Notes. Online Notes / Algebra (Notes) / Systems of Equations / Augmented Matricies
1 of 8 5/17/2011 5:58 PM Paul's Online Math Notes Home Class Notes Extras/Reviews Cheat Sheets & Tables Downloads Algebra Home Preliminaries Chapters Solving Equations and Inequalities Graphing and Functions
More informationChapter 6. Curves and Surfaces. 6.1 Graphs as Surfaces
Chapter 6 Curves and Surfaces In Chapter 2 a plane is defined as the zero set of a linear function in R 3. It is expected a surface is the zero set of a differentiable function in R n. To motivate, graphs
More informationMonday, 12 November 12. Matrices
Matrices Matrices Matrices are convenient way of storing multiple quantities or functions They are stored in a table like structure where each element will contain a numeric value that can be the result
More informationPreferred directions for resolving the non-uniqueness of Delaunay triangulations
Preferred directions for resolving the non-uniqueness of Delaunay triangulations Christopher Dyken and Michael S. Floater Abstract: This note proposes a simple rule to determine a unique triangulation
More informationADVANCES AT TELLES METHOD APPLIED IN SCIENTIFIC VISUALIZATION
Copyright 009 by ABCM January 04-08, 010, Foz do Iguaçu, PR, Brazil ADVANCES AT TELLES METHOD APPLIED IN SCIENTIFIC VISUALIZATION Carlos Andrés Reyna Vera-Tudela, candres@ufrrj.br Universidade Federal
More informationVisual Recognition: Image Formation
Visual Recognition: Image Formation Raquel Urtasun TTI Chicago Jan 5, 2012 Raquel Urtasun (TTI-C) Visual Recognition Jan 5, 2012 1 / 61 Today s lecture... Fundamentals of image formation You should know
More informationCurve Representation ME761A Instructor in Charge Prof. J. Ramkumar Department of Mechanical Engineering, IIT Kanpur
Curve Representation ME761A Instructor in Charge Prof. J. Ramkumar Department of Mechanical Engineering, IIT Kanpur Email: jrkumar@iitk.ac.in Curve representation 1. Wireframe models There are three types
More informationGeometric and Solid Modeling. Problems
Geometric and Solid Modeling Problems Define a Solid Define Representation Schemes Devise Data Structures Construct Solids Page 1 Mathematical Models Points Curves Surfaces Solids A shape is a set of Points
More informationAdvanced Operations Research Techniques IE316. Quiz 1 Review. Dr. Ted Ralphs
Advanced Operations Research Techniques IE316 Quiz 1 Review Dr. Ted Ralphs IE316 Quiz 1 Review 1 Reading for The Quiz Material covered in detail in lecture. 1.1, 1.4, 2.1-2.6, 3.1-3.3, 3.5 Background material
More informationarxiv: v1 [math.na] 20 Sep 2016
arxiv:1609.06236v1 [math.na] 20 Sep 2016 A Local Mesh Modification Strategy for Interface Problems with Application to Shape and Topology Optimization P. Gangl 1,2 and U. Langer 3 1 Doctoral Program Comp.
More informationA new Eulerian computational method for the propagation of short acoustic and electromagnetic pulses
A new Eulerian computational method for the propagation of short acoustic and electromagnetic pulses J. Steinhoff, M. Fan & L. Wang. Abstract A new method is described to compute short acoustic or electromagnetic
More informationVW 1LQH :HHNV 7KH VWXGHQW LV H[SHFWHG WR
PreAP Pre Calculus solve problems from physical situations using trigonometry, including the use of Law of Sines, Law of Cosines, and area formulas and incorporate radian measure where needed.[3e] What
More informationNumerical Integration
Numerical Integration Numerical Integration is the process of computing the value of a definite integral, when the values of the integrand function, are given at some tabular points. As in the case of
More informationGeometric transformations assign a point to a point, so it is a point valued function of points. Geometric transformation may destroy the equation
Geometric transformations assign a point to a point, so it is a point valued function of points. Geometric transformation may destroy the equation and the type of an object. Even simple scaling turns a
More information