BASE VECTORS FOR SOLVING OF PARTIAL DIFFERENTIAL EQUATIONS
|
|
- Jewel Nichols
- 6 years ago
- Views:
Transcription
1 BASE VECTORS FOR SOLVING OF PARTIAL DIFFERENTIAL EQUATIONS J. Roubal, V. Havlena Department of Control Engineering, Facult of Electrical Engineering, Czech Technical Universit in Prague Abstract The distributed parameters sstems can be described b linear two-dimensional (dependent on two spatial directions) parabolic partial differential equations. Using the finite difference method a distributed parameters sstem can be transformed to a linear discrete state space model. The controller design based on this description is complicated because of the large dimension of the model. Therefore, a model reduction method has to be used. We transform the state space model to the balanced realization of the sstem and show that the state vector of the model can be epressed as the series of columns of the transformation matri. These columns can be imaged as base vectors of the state space. 1 Introduction There are man industrial processes that have distributed parameters behaviour. Consequentl, these processes cannot be modelled b lumped inputs and/or lumped outputs models for correct representation. This paper deals with two-dimensional dnamic processes (sstems with parameters dependent on two spatial directions) that can be described b lumped inputs and distributed output models. These models can be mathematicall described b partial differential equations (PDE) [5]. Unlike ordinar differential equations, the PDEs contain, in addition, derivatives with respect to spatial directions. Consequentl, the partial differential equations lead to more accurate models but their compleit is larger. The dnamic behaviour of the distributed parameters sstem, which is described b the PDE, can be approimatel described b a finite-dimensional model, for eample, b using the finite difference method [1]. Then the ordinar differential equation model with large dimension is obtained and can be used for a finite-dimensional controller design. Unfortunatel, for online solving of an optimization problem, e.g. the model predictive control approach [6, 2], the large model dimension introduces a problem for the control design. Therefore, a model reduction method has to be used. Partial differential equations are usuall solved b base functions. This approach is based on that the solution of the PDEs can be epressed as series of the base functions. This paper presents similar methodolog. The discrete state space description of a distributed parameters sstem is transformed to the balanced realization [8]. Then the sstem state can be epressed as the product of the transformation matri and the balanced state. We can imagine the columns of the transformation matri as base vectors to which the sstem state can be epanded. Then using onl the first few base vectors corresponds to the truncation of states with smaller influence on the input/output behaviour of the sstem. The paper is organized as follows. In section 2, the distributed parameters model for the finite controller design is developed. In section 3, the balanced truncation method is described. In section 4, the relationship between the balanced realization and the sstem state epansion to an orthonormal base vectors is developed. In section 5, this approach is demonstrated on a heat transfer process.
2 2 Distributed Parameter Process Description In this section, the model of a heat transfer process described b a linear two-dimensional parabolic PDE [5] is developed for a finite-dimensional controller design. At first, the stationar PDE is transformed to a linear equation sstem using the finite difference approimation [1]. Then the implicit scheme [1] and this equation sstem are used for the transformation of the evolutionar PDE to a linear dnamic discrete sstem. 2.1 Stationar Partial Differential Equation For the surface thermal conductivit λ [W/K] independent on the temperature [K] and a surface heat source f [W/m 2 ], the heat transfer process in the stationar case can be described b a parabolic PDE ( 2 ) (, ) λ (, ) 2 = λ (, ) = f(, ), (1) where is the Laplace operator. The unknown temperature must satisf equation (1) on an open set Ω = (, L 1 ) (, L 2 ) and boundar condition on Ω ( Ω means the boundar of set Ω). In this paper, the boundar condition which specifies the temperature gradient on the boundar Ω is described b the following statement (, ) λ = α ( (, ) s (, ) ), (2) n where n is the unit normal vector, α [W/(mK)] is an eternal heat transfer coefficient and s is the surrounding temperature. Note that equation (2) is known as Newton (or the third kind) boundar condition [5]. i=,j= δ i=,j=n2 δ i,j-1 i-1,j i,j Ω i,j+1 i+1,j i=n1,j= Ω i=n1,j=n2 Figure 1: The mesh on the set Ω For the transformation of PDE (1) with boundar condition (2) to the finite dimensional model, the set Ω is covered b an imaginar mesh so that the values of the mesh points satisf i,j = (i δ, j δ) on the closed set Ω, F i,j = f(i δ, j δ) on the open set Ω and F i,j = s (i δ, j δ) on the set Ω, where δ and δ are the grid sizes of the imaginar mesh and i, j are row and column indices, respectivel (see Figure 1). Matri is the matri of the temperature values in the mesh points and matri F represents the heat source f and the surrounding temperature s.
3 When we define vectors ( :, ) ( :, 1) θ =., f = ( :, N 2 1) ( :, N 2 ) with indices l = j(n 1 + 1) + i, F ( :, ) F ( :, 1). F ( :, N 2 1) F ( :, N 2 ), (3) where (:, ) means the zero column of the matri, (:,1) the first column and so on, partial differential equation (1) with Newton boundar condition equations (2) can be written in compact form P θ = f. (4) More details were described in [7]. 2.2 Evolutionar Partial Differential Equation In the non stationar case, PDE (1) can be written as ρ c d(,, t) dt λ (,, t) = f(,, t), (5) where ρ [kg/m 2 ] is the surface densit of the medium and c [Ws/kgK] is its thermal capacit. In this case, the unknown temperature profile (,, t), dependent on time t, must satisf, for an initial condition (,, t ) = init (, ), equation (5) on the open set Ω and boundar condition (2) on Ω for all time horizon t t, t end. Using equation (4) and the implicit discretization scheme [1] with a sampling period δt, evolutionar PDE (5) with Newton boundar condition (2) can be approimated as θ(k + 1) = Mθ(k)+Nf(k), θ(k ) = θ init, (6) ( M = I + δt ) 1 ( P, N = I + δt ) δt P 1, ρ c ρ c ρ c where I is the identit matri with the corresponding dimension. 3 Model Reduction Method The accurac of model (6) increases with decreasing grid sizes δ and δ. Unfortunatel, for the advanced controller design, for eample the predictive controller, a low dimension model is needed. In this section, one reduction method is shortl described. 3.1 Model Reduction b the Balanced Truncation Method There are infinitel man different state space realizations for a given transfer function. But some realizations are more useful in control design. One of these realizations is the balanced realization which gives balanced Gramians for controllabilit W c and observabilit W o [8]. In addition, these Gramians are equal to the diagonal matri Σ W c = W o = Σ = diag (σ 1, σ 2,..., σ n ).
4 Note that the decreasingl ordered numbers, σ 1 σ 2... σ n, are called the Hankel singular values of the sstem. We suppose σ r σ r+1 for some r 1; n). Then the balanced realization implies that the states corresponding to the singular values of σ r+1,..., σ n are less controllable and observable than the states corresponding to σ 1,..., σ r. The states corresponding to the singular values of σ r+1,..., σ n have smaller influence on the input/output behaviour of the sstem. Therefore, truncating those less controllable and observable states will not lose much information about the sstem input/output behaviour and the dimension of the model can be significantl reduced. 3.2 Transformation to the Balanced Realization Model (6) can be epressed as (k+1) = A(k) + Bu(k) + Ez(k), (k )=, (k) = C(k) + Du(k), (7) where = θ is a state vector of the model (temperature profile), is its output vector (temperature measured in several points on the set Ω), u is its input vector (manipulated variable), z represents the surrounding temperature profile (measurable disturbance) and A, B, C, D, E are state matrices. Gramians for controllabilit W c and observabilit W o [4] of model (7) can be epressed b W c = W o = A i [ B, E ] [ B, E ] T ( A T ) i, (8) i= ( A T ) i C T C A i. (9) i= Model (7) can be transformed, b a matri Q to another realization = Qv, (1) v(k+1) = Av(k) + Bu(k) + Ez(k), v(k )=v, (k) = Cv(k) + Du(k), (11) where the matrices satisf A = Q 1 AQ, B = Q 1 B, E = Q 1 E, C = CQ, D = D. Note that Gramians of the transformed realization (11) are equal to W c = Q 1 W c (Q 1 ) T, W o = Q T W o Q, (12) = W c W o = Q 1 W c W o Q. (13) Now we have to find transformation matri Q, such that realization (11) is the balanced realization of the sstem. Since Gramians are smmetric positive definite matrices, we can factor Gramian for observabilit as W o = P T P, where matri P is the Cholesk factor of matri W o. Then equation (13) can be written as Q 1 W c P T P Q = Σ 2. Using elementar rearrangement, the equation can be obtained as Q 1 P 1 P W c P T P Q= ( P Q ) 1 P Wc P T( P Q ) =Σ 2. (14)
5 Equation (14) means that P W c P T is similar to Σ 2 and is positive definite. Therefore, there eists an orthonormal transformation matri U (U T U = I), such that P W c P T = UΣ 2 U T. (15) From equations (12)a and (15), we can derive the formula for the transformation matri Q = P 1 UΣ 1/2. (16) Then model (11) is the balanced realization of model (7) and its Gramians for controllabilit and observabilit are W c = W o = Σ [8]. So if model (7) is a minimal realization and matri A is stable, the transformation matri Q can be obtained through the following procedure: compute the controllabilit and observabilit Gramians of model (7) W c >, W o >, find matri P such that W o = P T P, diagonalize P W c P T = UΣ 2 U T b using the singular value decomposition [3], let the Q = P 1 UΣ 1/2. 4 Base Vectors Considering the Hankel singular values are decreasingl ordered, the first state of state vector v has the biggest influence on the input/output behaviour of the sstem, the second state has less influence and so on. As it has been said in section 3, we can truncate states corresponding to small singular values σ i and the model keeps enough information about the input/output behaviour of the sstem. In this section, we show that the sstem state vector can be epressed in an orthonormal base. Then using onl the first r orthonormal base vectors (for epression of the sstem state vector) corresponds to the truncation of the states with smaller influence on the input/output behaviour of the sstem. 4.1 Decomposition of the Sstem State Vector to the Balanced Base Vectors Equation (1) can be rewritten as = Qv = n v i q i, (17) i=1 where v i are the elements of state vector v in balanced realization (11) and vectors q i are the columns of matri Q = [ ] q 1, q 2,..., q n. Since the transformation matri Q is not singular, the vectors q i are linearl independent. So equation (17) means that state vector can be epressed as the series of base vectors q i that are multiplied b weighted coefficients v i. Note that if we want to truncate states corresponding to the singular values of σ r+1,..., σ n, we epress state vector onl b the first r base vectors q i = r v i q i. i=1
6 4.2 Orthonormal Base Vectors The base vectors q i can be transformed to an orthonormal base b using the Gram-Schmidt orthonormalization. We decompose the matri Q such that Q = OR, where R is an upper triangular matri and O is an orthonormal matri (O T O = I). Then equation (17) can be rewritten as n = Qv = ORv = Oṽ = ṽ i o i, (18) where ṽ = Rv and vectors o i are the columns of matri O = [ o 1, o 2,..., o n ]. So equation (18) means that state vector can be epressed b the orthonormal base vectors o i. i=1 5 DEMONSTRATION EXAMPLE Consider a heat transfer process in a furnace where L 1 = L 2 =.9 m described b equation (5) with constants λ = 51 W/K, ρ = 25 kg/m 2, c = 1259 Ws/(kg K) and α = 1.14 W/(mK). The grid sizes are δ = δ =.2 m and the sampling period is δt = 3 s. The sstem has five lumped inputs (manipulated variables, see Figure 2) and the surface temperature is measured in 64 points which are uniforml distributed over the area Ω, see Figure 3b. Source distribution Figure 2: The heat source distribution f(, ) Stead state temperature profile Sstem output end Figure 3: (a) The stead-state temperature distribution (, ); (b) The sstem output (the temperature in 64 measurement points)
7 Hankel Singular Sstem Values σ > : (129) σ >.1: (18) σ >.1: (13) σ > 1: (9) σ > 1: (5) σ > 1: (1) σ Figure 4: The Hankel singular values of the sstem Figure 3a presents the stead-state temperature distribution (sstem state) for the unit step as input signal (see Figure 2) and the surrounding temperature s = 34 K. Figure 3b shows the sstem output temperature in 64 measurement points. Figure 4 shows the Hankel singular values of the sstem. Figure 5a presents the sstem output based on the first sstem state (in balanced realization), Figure 5b shows the base function corresponding to the first base vector (17) and Figure 5c shows the orthonormal base function corresponding to the first orthonormal base vector (18). Figure 6a presents the sstem output based on the first two sstem states and Figures 6b and 6c show the base function corresponding to the second base vector and the orthonormal base function corresponding to the second orthonormal base vector, respectivel. From Figures 5-14, it follows that the sstem state vector can be epressed in the series of the base functions (the columns of the transformation matri) or the orthonormal base function, respectivel. This methodolog is similar to the analtical solution of PDEs based on the nondiscrete base functions that are usuall nonlinear even if the PDEs are linear. Sstem output Base vector q 1 Orthonormal base vector o Figure 5: (a) Stead-state sstem output based on the first state; function (q 1 ); (c) The first orthonormal base function (o 1 ) (b) The first base Sstem output Base vector q 2 Orthonormal base vector o Figure 6: (a) Stead-state sstem output based on the first two states; (b) The second base function (q 2 ; (c) The second orthonormal base function (o 2 )
8 Sstem output Base vector q 3 Orthonormal base vector o Figure 7: (a) Stead-state sstem output based on the first three states; (b) The third base function (q 3 ); (c) The third orthonormal base function (o 3 ) Sstem output Base vector q 4 Orthonormal base vector o Figure 8: (a) Stead-state sstem output based on the first four states; (b) The fourth base function (q 4 ); (c) The fourth orthonormal base function (o 4 ) Sstem output Base vector q 5 Orthonormal base vector o Figure 9: (a) Stead-state sstem output based on the first five states; (b) The fifth base function (q 5 ); (c) The fifth orthonormal base function (o 5 ) Sstem output Base vector q 6 Orthonormal base vector o Figure 1: (a) Stead-state sstem output based on the first si states; (b) The sith base function (q 6 ); (c) The sith orthonormal base function (o 6 )
9 Sstem output Base vector q 7 Orthonormal base vector o Figure 11: (a) Stead-state sstem output based on the first seven states; (b) The seventh base function (q 7 ); (c) The seventh orthonormal base function (o 7 ) Sstem output Base vector q 8 Orthonormal base vector o Figure 12: (a) Stead-state sstem output based on the first eight states; (b) The eighth base function (q 8 ); (c) The eighth orthonormal base function (o 8 ) Sstem output Base vector q 9 Orthonormal base vector o Figure 13: (a) Stead-state sstem output based on the first nine states; (b) The ninth base function (q 9 ); (c) The ninth orthonormal base function (o 9 ) Sstem output Base vector q 1 Orthonormal base vector o Figure 14: (a) Stead-state sstem output based on the first five states; (b) The tenth base function (q 1 ); (c) The tenth orthonormal base function (o 1 )
10 6 CONCLUSION The state space model of the distributed parameters sstem described b the linear two-dimensional parabolic partial differential equation and the model reduction b the balanced truncation method are described. The connection of the model reduction b the balanced truncation method and the solution of the partial differential equations based on the base functions is developed. The result is demonstrated on the heat transfer process eample. You can observe from the figures that the base functions q i and the orthonormal base functions o i, respectivel, are similar to the Fourier decomposition of a signal (the first harmonic, the second harmonic, etc.) Acknowledgements This work was partl supported b the Ministr of Industr and Trade of the Czech Republic under project 1H-PK/22. References [1] Babuška, I., Práger, M., and Vitásek, E. Numerical Processes in Differential Equations. London New York Sdne, Interscience Publishers, [2] Havlena, V. and Findejs, J. Application of Model Predictive Control to Advanced Combustion Control. Control Engineering Practice, 13 (No.6):671 68, 25. [3] Horn, R. A. and Johnson, Ch. R. Matri Analsis. Cambridge Universit Press, United Kingdom, [4] Kailath, T. Linear Sstems. Prentice Hall, New Jerse, 198. [5] Long, C. A. Essential Heat Transfer. Longman Group, England, [6] Roubal, J. and Havlena, V. Range Control MPC Approach for Two-Dimensional Sstem. In Proceedings of IFAC 16th World Congress, Prague, Czech Republic, 25. [7] Roubal, J., Trnka, P., and Havlena, V. Simulation of Two-Dimensional Distributed Parameter Sstems. In Matlab Conference, Prague, Czech Republic, [8] Zhou, K., Dole, J. C., and Glover, K. Robust and Optimal Control. Prentice Hall, New Jerse, Jirka Roubal Department of Control Engineering, Facult of Electrical Engineering, Czech Technical Universit in Prague, Karlovo náměstí 13, Prague 2, Czech Republic, roubalj@control.felk.cvut.cz Vladimír Havlena Department of Control Engineering, Facult of Electrical Engineering, Czech Technical Universit in Prague, Karlovo náměstí 13, Prague 2, Czech Republic, havlena@control.felk.cvut.cz
RANGE CONTROL MPC APPROACH FOR TWO-DIMENSIONAL SYSTEM 1
RANGE CONTROL MPC APPROACH FOR TWO-DIMENSIONAL SYSTEM Jirka Roubal Vladimír Havlena Department of Control Engineering, Facult of Electrical Engineering, Czech Technical Universit in Prague Karlovo náměstí
More informationRECURSIVE SUBSPACE IDENTIFICATION IN THE LEAST SQUARES FRAMEWORK
RECURSIVE SUBSPACE IDENTIFICATION IN THE LEAST SQUARES FRAMEWORK TRNKA PAVEL AND HAVLENA VLADIMÍR Dept of Control Engineering, Czech Technical University, Technická 2, 166 27 Praha, Czech Republic mail:
More informationFACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING. Lectures AB = BA = I,
FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING Lectures MODULE 7 MATRICES II Inverse of a matri Sstems of linear equations Solution of sets of linear equations elimination methods 4
More informationIntroduction to Differential Equations. National Chiao Tung University Chun-Jen Tsai 9/14/2011
Introduction to Differential Equations National Chiao Tung Universit Chun-Jen Tsai 9/14/011 Differential Equations Definition: An equation containing the derivatives of one or more dependent variables,
More informationEigenvectors and Eigenvalues 1
Ma 2015 page 1 Eigenvectors and Eigenvalues 1 In this handout, we will eplore eigenvectors and eigenvalues. We will begin with an eploration, then provide some direct eplanation and worked eamples, and
More information1.1 The Equations of Motion
1.1 The Equations of Motion In Book I, balance of forces and moments acting on an component was enforced in order to ensure that the component was in equilibrium. Here, allowance is made for stresses which
More informationDEVELOPMENT AND USE OF OPERATIONAL MODAL ANALYSIS
Università degli Studi di Napoli Federico II DEVELOPMENT AND USE OF OPERATIONAL MODAL ANALYSIS Examples of Civil, Aeronautical and Acoustical Applications Francesco Marulo Tiziano Polito francesco.marulo@unina.it
More informationTowards Higher-Order Schemes for Compressible Flow
WDS'6 Proceedings of Contributed Papers, Part I, 5 4, 6. ISBN 8-867-84- MATFYZPRESS Towards Higher-Order Schemes for Compressible Flow K. Findejs Charles Universit, Facult of Mathematics and Phsics, Prague,
More informationChapter Adequacy of Solutions
Chapter 04.09 dequac of Solutions fter reading this chapter, ou should be able to: 1. know the difference between ill-conditioned and well-conditioned sstems of equations,. define the norm of a matri,
More informationSimultaneous Orthogonal Rotations Angle
ELEKTROTEHNIŠKI VESTNIK 8(1-2): -11, 2011 ENGLISH EDITION Simultaneous Orthogonal Rotations Angle Sašo Tomažič 1, Sara Stančin 2 Facult of Electrical Engineering, Universit of Ljubljana 1 E-mail: saso.tomaic@fe.uni-lj.si
More informationSection B. Ordinary Differential Equations & its Applications Maths II
Section B Ordinar Differential Equations & its Applications Maths II Basic Concepts and Ideas: A differential equation (D.E.) is an equation involving an unknown function (or dependent variable) of one
More informationInitial Value Problems for. Ordinary Differential Equations
Initial Value Problems for Ordinar Differential Equations INTRODUCTION Equations which are composed of an unnown function and its derivatives are called differential equations. It becomes an initial value
More informationUnit 12 Study Notes 1 Systems of Equations
You should learn to: Unit Stud Notes Sstems of Equations. Solve sstems of equations b substitution.. Solve sstems of equations b graphing (calculator). 3. Solve sstems of equations b elimination. 4. Solve
More informationLocal Phase Portrait of Nonlinear Systems Near Equilibria
Local Phase Portrait of Nonlinear Sstems Near Equilibria [1] Consider 1 = 6 1 1 3 1, = 3 1. ( ) (a) Find all equilibrium solutions of the sstem ( ). (b) For each equilibrium point, give the linear approimating
More informationIntroduction to Differential Equations
Math0 Lecture # Introduction to Differential Equations Basic definitions Definition : (What is a DE?) A differential equation (DE) is an equation that involves some of the derivatives (or differentials)
More informationGet Solution of These Packages & Learn by Video Tutorials on Matrices
FEE Download Stud Package from website: wwwtekoclassescom & wwwmathsbsuhagcom Get Solution of These Packages & Learn b Video Tutorials on wwwmathsbsuhagcom Matrices An rectangular arrangement of numbers
More informationAll parabolas through three non-collinear points
ALL PARABOLAS THROUGH THREE NON-COLLINEAR POINTS 03 All parabolas through three non-collinear points STANLEY R. HUDDY and MICHAEL A. JONES If no two of three non-collinear points share the same -coordinate,
More informationAssignment #9: Orthogonal Projections, Gram-Schmidt, and Least Squares. Name:
Assignment 9: Orthogonal Projections, Gram-Schmidt, and Least Squares Due date: Friday, April 0, 08 (:pm) Name: Section Number Assignment 9: Orthogonal Projections, Gram-Schmidt, and Least Squares Due
More informationc 1999 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol., No. 6, pp. 978 994 c 999 Societ for Industrial and Applied Mathematics A STUDY OF MONITOR FUNCTIONS FOR TWO-DIMENSIONAL ADAPTIVE MESH GENERATION WEIMING CAO, WEIZHANG HUANG,
More informationMMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM
B Course Content: A INTRODUCTION AND OVERVIEW Numerical method and Computer-Aided Engineering; Phsical problems; Mathematical models; Finite element method;. B Elements and nodes, natural coordinates,
More informationNONLINEAR DYNAMICS AND CHAOS. Numerical integration. Stability analysis
LECTURE 3: FLOWS NONLINEAR DYNAMICS AND CHAOS Patrick E McSharr Sstems Analsis, Modelling & Prediction Group www.eng.o.ac.uk/samp patrick@mcsharr.net Tel: +44 83 74 Numerical integration Stabilit analsis
More informationDemonstrate solution methods for systems of linear equations. Show that a system of equations can be represented in matrix-vector form.
Chapter Linear lgebra Objective Demonstrate solution methods for sstems of linear equations. Show that a sstem of equations can be represented in matri-vector form. 4 Flowrates in kmol/hr Figure.: Two
More informationLecture 4.2 Finite Difference Approximation
Lecture 4. Finite Difference Approimation 1 Discretization As stated in Lecture 1.0, there are three steps in numerically solving the differential equations. They are: 1. Discretization of the domain by
More informationOBJECTIVE 5 SOLVING SYSTEMS 5/19/2016 SOLVING SYSTEMS OF TWO EQUATIONS BY SUBSTITUTION:
/9/ OBJECTIVE Sstems & Matrices SOLVING SYSTEMS OF TWO EQUATIONS BY SUBSTITUTION:. Solve one of the equations for one of the variables in terms of the other.. Substitute this epression into the nd equation,
More informationHandout for Adequacy of Solutions Chapter SET ONE The solution to Make a small change in the right hand side vector of the equations
Handout for dequac of Solutions Chapter 04.07 SET ONE The solution to 7.999 4 3.999 Make a small change in the right hand side vector of the equations 7.998 4.00 3.999 4.000 3.999 Make a small change in
More informationTransformation of kinematical quantities from rotating into static coordinate system
Transformation of kinematical quantities from rotating into static coordinate sstem Dimitar G Stoanov Facult of Engineering and Pedagog in Sliven, Technical Universit of Sofia 59, Bourgasko Shaussee Blvd,
More informationSolving Variable-Coefficient Fourth-Order ODEs with Polynomial Nonlinearity by Symmetric Homotopy Method
Applied and Computational Mathematics 218; 7(2): 58-7 http://www.sciencepublishinggroup.com/j/acm doi: 1.11648/j.acm.21872.14 ISSN: 2328-565 (Print); ISSN: 2328-5613 (Online) Solving Variable-Coefficient
More informationTHE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY. V. A. Sharifulin.
THE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY 1. Introduction V. A. Sharifulin Perm State Technical Universit, Perm, Russia e-mail: sharifulin@perm.ru Water
More informationEVALUATION OF THERMAL TRANSPORT PROPERTIES USING A MICRO-CRACKING MODEL FOR WOVEN COMPOSITE LAMINATES
EVALUATION OF THERMAL TRANSPORT PROPERTIES USING A MICRO-CRACKING MODEL FOR WOVEN COMPOSITE LAMINATES C. Luo and P. E. DesJardin* Department of Mechanical and Aerospace Engineering Universit at Buffalo,
More informationFully Discrete Energy Stable High Order Finite Difference Methods for Hyperbolic Problems in Deforming Domains: An Initial Investigation
Full Discrete Energ Stable High Order Finite Difference Methods for Hperbolic Problems in Deforming Domains: An Initial Investigation Samira Nikkar and Jan Nordström Abstract A time-dependent coordinate
More informationEstimation Of Linearised Fluid Film Coefficients In A Rotor Bearing System Subjected To Random Excitation
Estimation Of Linearised Fluid Film Coefficients In A Rotor Bearing Sstem Subjected To Random Ecitation Arshad. Khan and Ahmad A. Khan Department of Mechanical Engineering Z.. College of Engineering &
More information15. Eigenvalues, Eigenvectors
5 Eigenvalues, Eigenvectors Matri of a Linear Transformation Consider a linear ( transformation ) L : a b R 2 R 2 Suppose we know that L and L Then c d because of linearit, we can determine what L does
More informationCONSERVATION LAWS AND CONSERVED QUANTITIES FOR LAMINAR RADIAL JETS WITH SWIRL
Mathematical and Computational Applications,Vol. 15, No. 4, pp. 742-761, 21. c Association for Scientific Research CONSERVATION LAWS AND CONSERVED QUANTITIES FOR LAMINAR RADIAL JETS WITH SWIRL R. Naz 1,
More informationKINEMATIC RELATIONS IN DEFORMATION OF SOLIDS
Chapter 8 KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS Figure 8.1: 195 196 CHAPTER 8. KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS 8.1 Motivation In Chapter 3, the conservation of linear momentum for a
More informationWavelet-Based Linearization Method in Nonlinear Systems
Wavelet-Based Linearization Method in Nonlinear Sstems Xiaomeng Ma A Thesis In The Department of Mechanical and Industrial Engineering Presented in Partial Fulfillment of the Requirements for the Degree
More informationINTRODUCTION TO DIFFERENTIAL EQUATIONS
INTRODUCTION TO DIFFERENTIAL EQUATIONS. Definitions and Terminolog. Initial-Value Problems.3 Differential Equations as Mathematical Models CHAPTER IN REVIEW The words differential and equations certainl
More informationA comparison of estimation accuracy by the use of KF, EKF & UKF filters
Computational Methods and Eperimental Measurements XIII 779 A comparison of estimation accurac b the use of KF EKF & UKF filters S. Konatowski & A. T. Pieniężn Department of Electronics Militar Universit
More informationChapter 3. Theory of measurement
Chapter. Introduction An energetic He + -ion beam is incident on thermal sodium atoms. Figure. shows the configuration in which the interaction one is determined b the crossing of the laser-, sodium- and
More informationInternational Journal of Applied Mathematics and Physics, 3(2), July-December 2011, pp Global Research Publications, India
International Journal of Applied Mathematics and Phsics, 3(), Jul-December 0, pp. 55-67 Global Research Publications, India Effects of Chemical Reaction with Heat and Mass Transfer on Peristaltic Flow
More informationIn this chapter a student has to learn the Concept of adjoint of a matrix. Inverse of a matrix. Rank of a matrix and methods finding these.
MATRICES UNIT STRUCTURE.0 Objectives. Introduction. Definitions. Illustrative eamples.4 Rank of matri.5 Canonical form or Normal form.6 Normal form PAQ.7 Let Us Sum Up.8 Unit End Eercise.0 OBJECTIVES In
More informationLATERAL BUCKLING ANALYSIS OF ANGLED FRAMES WITH THIN-WALLED I-BEAMS
Journal of arine Science and J.-D. Technolog, Yau: ateral Vol. Buckling 17, No. Analsis 1, pp. 9-33 of Angled (009) Frames with Thin-Walled I-Beams 9 ATERA BUCKING ANAYSIS OF ANGED FRAES WITH THIN-WAED
More informationKrylov Integration Factor Method on Sparse Grids for High Spatial Dimension Convection Diffusion Equations
DOI.7/s95-6-6-7 Krlov Integration Factor Method on Sparse Grids for High Spatial Dimension Convection Diffusion Equations Dong Lu Yong-Tao Zhang Received: September 5 / Revised: 9 March 6 / Accepted: 8
More information12.1 Systems of Linear equations: Substitution and Elimination
. Sstems of Linear equations: Substitution and Elimination Sstems of two linear equations in two variables A sstem of equations is a collection of two or more equations. A solution of a sstem in two variables
More informationExact and Approximate Numbers:
Eact and Approimate Numbers: The numbers that arise in technical applications are better described as eact numbers because there is not the sort of uncertainty in their values that was described above.
More informationApplications of Gauss-Radau and Gauss-Lobatto Numerical Integrations Over a Four Node Quadrilateral Finite Element
Avaiable online at www.banglaol.info angladesh J. Sci. Ind. Res. (), 77-86, 008 ANGLADESH JOURNAL OF SCIENTIFIC AND INDUSTRIAL RESEARCH CSIR E-mail: bsir07gmail.com Abstract Applications of Gauss-Radau
More informationClosed form expressions for the gravitational inner multipole moments of homogeneous elementary solids
Closed form epressions for the gravitational inner multipole moments of homogeneous elementar solids Julian Stirling 1,2, and Stephan Schlamminger 1 1 National Institute of Standards and Technolog, 1 Bureau
More informationA Low-Complexity Approach to Computation of the Discrete Fractional Fourier Transform
Circuits Sst Signal Process (1) 36:411 4144 DOI 11/s34-1-53-z A Low-Compleit Approach to Computation of the Discrete Fractional Fourier Transform Dorota Majorkowska-Mech 1 Aleksandr Cariow 1 Received:
More informationSlovak University of Technology in Bratislava Institute of Information Engineering, Automation, and Mathematics PROCEEDINGS
Slovak Universit of Technolog in Bratislava Institute of Information Engineering, Automation, and Mathematics PROCEEDINGS of the 8 th International Conference on Process Control Hotel Titris, Tatranská
More informationCOMPACT IMPLICIT INTEGRATION FACTOR METHODS FOR A FAMILY OF SEMILINEAR FOURTH-ORDER PARABOLIC EQUATIONS. Lili Ju. Xinfeng Liu.
DISCRETE AND CONTINUOUS doi:13934/dcdsb214191667 DYNAMICAL SYSTEMS SERIES B Volume 19, Number 6, August 214 pp 1667 1687 COMPACT IMPLICIT INTEGRATION FACTOR METHODS FOR A FAMILY OF SEMILINEAR FOURTH-ORDER
More informationModule 1 : The equation of continuity. Lecture 4: Fourier s Law of Heat Conduction
1 Module 1 : The equation of continuit Lecture 4: Fourier s Law of Heat Conduction NPTEL, IIT Kharagpur, Prof. Saikat Chakrabort, Department of Chemical Engineering Fourier s Law of Heat Conduction According
More informationVibrational Power Flow Considerations Arising From Multi-Dimensional Isolators. Abstract
Vibrational Power Flow Considerations Arising From Multi-Dimensional Isolators Rajendra Singh and Seungbo Kim The Ohio State Universit Columbus, OH 4321-117, USA Abstract Much of the vibration isolation
More informationNUMERICAL SIMULATION OF DISLOCATION DYNAMICS BY MEANS OF PARAMETRIC APPROACH
Proceedings of Czech Japanese Seminar in Applied Mathematics August -7,, Czech Technical Universit in Prague http://geraldine.fjfi.cvut.cz pp. 1 13 NUMERICAL SIMULATION OF DISLOCATION DYNAMICS BY MEANS
More information4 Strain true strain engineering strain plane strain strain transformation formulae
4 Strain The concept of strain is introduced in this Chapter. The approimation to the true strain of the engineering strain is discussed. The practical case of two dimensional plane strain is discussed,
More informationm x n matrix with m rows and n columns is called an array of m.n real numbers
LINEAR ALGEBRA Matrices Linear Algebra Definitions m n matri with m rows and n columns is called an arra of mn real numbers The entr a a an A = a a an = ( a ij ) am am amn a ij denotes the element in the
More informationModule 2 Stresses in machine elements. Version 2 ME, IIT Kharagpur
Module Stresses in machine elements Version M, IIT Kharagpur Lesson 3 Strain analsis Version M, IIT Kharagpur Instructional Objectives At the end of this lesson, the student should learn Normal and shear
More informationAn Optimization Method for Numerically Solving Three-point BVPs of Linear Second-order ODEs with Variable Coefficients
An Optimization Method for Numericall Solving Three-point BVPs of Linear Second-order ODEs with Variable Coefficients Liaocheng Universit School of Mathematics Sciences 252059, Liaocheng P.R. CHINA hougbb@26.com;
More informationMATHEMATICAL FUNDAMENTALS I. Michele Fitzpatrick
MTHEMTICL FUNDMENTLS I Michele Fitpatrick OVERVIEW Vectors and arras Matrices Linear algebra Del( operator Tensors DEFINITIONS vector is a single row or column of numbers. n arra is a collection of vectors
More informationFixed Point Theorem and Sequences in One or Two Dimensions
Fied Point Theorem and Sequences in One or Two Dimensions Dr. Wei-Chi Yang Let us consider a recursive sequence of n+ = n + sin n and the initial value can be an real number. Then we would like to ask
More informationc 2016 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 38, No. 2, pp. A765 A788 c 206 Societ for Industrial and Applied Mathematics ROOTS OF BIVARIATE POLYNOMIAL SYSTEMS VIA DETERMINANTAL REPRESENTATIONS BOR PLESTENJAK AND MICHIEL
More informationEFFECT OF DAMPING AND THERMAL GRADIENT ON VIBRATIONS OF ORTHOTROPIC RECTANGULAR PLATE OF VARIABLE THICKNESS
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 193-9466 Vol. 1, Issue 1 (June 17), pp. 1-16 Applications and Applied Mathematics: An International Journal (AAM) EFFECT OF DAMPING AND THERMAL
More information16.5. Maclaurin and Taylor Series. Introduction. Prerequisites. Learning Outcomes
Maclaurin and Talor Series 6.5 Introduction In this Section we eamine how functions ma be epressed in terms of power series. This is an etremel useful wa of epressing a function since (as we shall see)
More informationAnswers Investigation 3
Answers Investigation Applications. a., b. s = n c. The numbers seem to be increasing b a greater amount each time. The square number increases b consecutive odd integers:,, 7,, c X X=. a.,,, b., X 7 X=
More informationDynamics of multiple pendula without gravity
Chaotic Modeling and Simulation (CMSIM) 1: 57 67, 014 Dnamics of multiple pendula without gravit Wojciech Szumiński Institute of Phsics, Universit of Zielona Góra, Poland (E-mail: uz88szuminski@gmail.com)
More informationAN EFFICIENT MOVING MESH METHOD FOR A MODEL OF TURBULENT FLOW IN CIRCULAR TUBES *
Journal of Computational Mathematics, Vol.27, No.2-3, 29, 388 399. AN EFFICIENT MOVING MESH METHOD FOR A MODEL OF TURBULENT FLOW IN CIRCULAR TUBES * Yin Yang Hunan Ke Laborator for Computation and Simulation
More informationAnalytic Solutions of Partial Differential Equations
ii Analtic Solutions of Partial Differential Equations 15 credits Taught Semester 1, Year running 3/4 MATH3414 School of Mathematics, Universit of Lee Pre-requisites MATH36 or MATH4 or equivalent. Co-requisites
More informationStability Analysis of Laminated Composite Thin-Walled Beam Structures
Paper 224 Civil-Comp Press, 2012 Proceedings of the Eleventh International Conference on Computational Structures echnolog, B.H.V. opping, (Editor), Civil-Comp Press, Stirlingshire, Scotland Stabilit nalsis
More informationR S T. PreCalculus AB Final Exam SHOW YOUR WORK May 20, Name: 1. Find the area of this triangle. 2. Find the area of this trapezoid.
1. Find the area of this triangle. 138 ft 6 18 ft. Find the area of this trapezoid. 10 ft 8 ft 57 11 ft 3. Find the area of this trapezoid. 10 ft 8 ft 59 1 ft [A] 88 ft [B] 176 ft [C] 75.3 ft [D] 8.9 ft.
More informationEntropy ISSN
Entrop 5, 7[4], 3 37 3 Entrop ISSN 99-43 www.mdpi.org/entrop/ Full Paper Second law analsis of laminar flow in a channel filled with saturated porous media: a numerical solution Kamel Hooman, Arash Ejlali
More informationCONSERVATION OF ANGULAR MOMENTUM FOR A CONTINUUM
Chapter 4 CONSERVATION OF ANGULAR MOMENTUM FOR A CONTINUUM Figure 4.1: 4.1 Conservation of Angular Momentum Angular momentum is defined as the moment of the linear momentum about some spatial reference
More informationIdentification of Nonlinear Dynamic Systems with Multiple Inputs and Single Output using discrete-time Volterra Type Equations
Identification of Nonlinear Dnamic Sstems with Multiple Inputs and Single Output using discrete-time Volterra Tpe Equations Thomas Treichl, Stefan Hofmann, Dierk Schröder Institute for Electrical Drive
More information1 Differential Equations for Solid Mechanics
1 Differential Eqations for Solid Mechanics Simple problems involving homogeneos stress states have been considered so far, wherein the stress is the same throghot the component nder std. An eception to
More informationGROUP CLASSIFICATION OF LINEAR SECOND-ORDER DELAY ORDINARY DIFFERENTIAL EQUATION
Proceedings of the nd IMT-GT Regional onference on Mathematics Statistics and Applications Universiti Sains Malasia GROUP LASSIFIATION OF LINEAR SEOND-ORDER DELAY ORDINARY DIFFERENTIAL EQUATION Prapart
More informationConservation of Linear Momentum for a Differential Control Volume
Conservation of Linear Momentum for a Differential Control Volume When we applied the rate-form of the conservation of mass equation to a differential control volume (open sstem in Cartesian coordinates,
More informationOrdinary Differential Equations n
Numerical Analsis MTH63 Ordinar Differential Equations Introduction Talor Series Euler Method Runge-Kutta Method Predictor Corrector Method Introduction Man problems in science and engineering when formulated
More informationLecture 5. Equations of Lines and Planes. Dan Nichols MATH 233, Spring 2018 University of Massachusetts.
Lecture 5 Equations of Lines and Planes Dan Nichols nichols@math.umass.edu MATH 233, Spring 2018 Universit of Massachusetts Februar 6, 2018 (2) Upcoming midterm eam First midterm: Wednesda Feb. 21, 7:00-9:00
More informationVibration of Plate on Foundation with Four Edges Free by Finite Cosine Integral Transform Method
854 Vibration of Plate on Foundation with Four Edges Free b Finite Cosine Integral Transform Method Abstract The analtical solutions for the natural frequencies and mode shapes of the rectangular plate
More informationConsider a slender rod, fixed at one end and stretched, as illustrated in Fig ; the original position of the rod is shown dotted.
4.1 Strain If an object is placed on a table and then the table is moved, each material particle moves in space. The particles undergo a displacement. The particles have moved in space as a rigid bod.
More informationMATRIX TRANSFORMATIONS
CHAPTER 5. MATRIX TRANSFORMATIONS INSTITIÚID TEICNEOLAÍOCHTA CHEATHARLACH INSTITUTE OF TECHNOLOGY CARLOW MATRIX TRANSFORMATIONS Matri Transformations Definition Let A and B be sets. A function f : A B
More informationInternational Journal of Mathematical Engineering and Science ISSN : Volume 1 Issue 2
ISSN : 77-698 Volume Issue The Mathematical Structure and Analsis of an MHD Flow in Porous Media () () Naji Qatanani and Mai Musmar () Department of Mathematics, An-Najah National Universit Nablus-Palestine
More informationPrepared by: M. S. KumarSwamy, TGT(Maths) Page
Prepared b: M. S. KumarSwam, TGT(Maths) Page - 77 - CHAPTER 4: DETERMINANTS QUICK REVISION (Important Concepts & Formulae) Determinant a b If A = c d, then determinant of A is written as A = a b = det
More informationAdditional Topics in Differential Equations
6 Additional Topics in Differential Equations 6. Eact First-Order Equations 6. Second-Order Homogeneous Linear Equations 6.3 Second-Order Nonhomogeneous Linear Equations 6.4 Series Solutions of Differential
More informationTwo-Dimensional Analysis of the Power Transfer between Crossed Laser Beams
Two-Dimensional Analsis of the Power Transfer between Crossed Laser Beams The indirect-drive approach to inertial confinement fusion involves laser beams that cross as the enter the hohlraum. Ionacoustic
More informationCurriculum: MechEng, UC Berkeley
Curriculum: MechEng, UC Berkele ME 132, 3-units, required, 2 nd or 3 rd ear students Prerequisites: Programming and Scientific computing (E7); Linear Algebra and ODEs (Math 54) Emphasis: mathematics of
More informationFréchet Sensitivity Analysis for Partial Differential Equations with Distributed Parameters
211 American Control Conference on O'Farrell Street, San Francisco, CA, USA June 29 - July 1, 211 Fréchet Sensitivity Analysis for Partial Differential Equations with Distributed Parameters Jeff Borggaard*
More informationCramér-Rao Bounds for Estimation of Linear System Noise Covariances
Journal of Mechanical Engineering and Automation (): 6- DOI: 593/jjmea Cramér-Rao Bounds for Estimation of Linear System oise Covariances Peter Matiso * Vladimír Havlena Czech echnical University in Prague
More informationChapter 6 Controllability and Obervability
Chapter 6 Controllability and Obervability Controllability: whether or not the state-space equation can be controlled from input. Observability: whether or not the initial state can be observed from output.
More informationIn applications, we encounter many constrained optimization problems. Examples Basis pursuit: exact sparse recovery problem
1 Conve Analsis Main references: Vandenberghe UCLA): EECS236C - Optimiation methods for large scale sstems, http://www.seas.ucla.edu/ vandenbe/ee236c.html Parikh and Bod, Proimal algorithms, slides and
More informationarxiv: v2 [physics.comp-ph] 24 Sep 2015
Minimal subspace rotation on the Stiefel manifold for stabilization and enhancement of projection-based reduced order models for the compressible Navier-Stokes equations Maciej Balajewicz a,1,, Irina Tezaur
More informationExponential time differencing methods for nonlinear PDEs
XXI Congreso de Ecuaciones Diferenciales Aplicaciones XI Congreso de Matemática Aplicada Ciudad Real, 1- septiembre 9 (pp. 1 8) Exponential time differencing methods for nonlinear PDEs F. de la Hoz 1,
More informationNew Scientific Contribution on the 2-D Subdomain Technique in Cartesian Coordinates: Taking into Account of Iron Parts
Mathematical and Computational pplications rticle New Scientific Contribution on the 2-D Subdomain Technique in Cartesian Coordinates: Taking into ccount of Iron Parts Frédéric Dubas 1, * and Kamel Boughrara
More informationAPPENDIX D Rotation and the General Second-Degree Equation
APPENDIX D Rotation and the General Second-Degree Equation Rotation of Aes Invariants Under Rotation After rotation of the - and -aes counterclockwise through an angle, the rotated aes are denoted as the
More informationRealization of Quantum Hadamard Gate by Applying Optimal Control Fields to a Spin Qubit
00 nd International Conference on Mechanical and Electronics Engineering (ICMEE 00) Realiation of Quantum Hadamard Gate b Appling Optimal Control Fields to a Spin Qubit Sh Mohammad Nejad Nanoptronics Research
More informationTime-Frequency Analysis: Fourier Transforms and Wavelets
Chapter 4 Time-Frequenc Analsis: Fourier Transforms and Wavelets 4. Basics of Fourier Series 4.. Introduction Joseph Fourier (768-83) who gave his name to Fourier series, was not the first to use Fourier
More informationAn improved ADI-DQM based on Bernstein polynomial for solving two-dimensional convection-diffusion equations
Mathematical Theor and Modeling ISSN -50 (Paper) ISSN 5-05(Online) Vol 3, No.1, 013 An improved ADI-DQM based on Bernstein polnomial for solving two-dimensional convection-diffusion equations A.S.J. Al-
More informationPOLYNOMIAL EXPANSION OF THE PROBABILITY DENSITY FUNCTION ABOUT GAUSSIAN MIXTURES
POLYNOMIAL EXPANSION OF THE PROBABILITY DENSITY FUNCTION ABOUT GAUSSIAN MIXTURES Charles C. Cavalcante, J. C. M. Mota and J. M. T. Romano GTEL/DETI/CT/UFC C.P. 65, Campus do Pici, CEP: 6.455-76 Fortaleza-CE,
More informationPart D. Complex Analysis
Part D. Comple Analsis Chapter 3. Comple Numbers and Functions. Man engineering problems ma be treated and solved b using comple numbers and comple functions. First, comple numbers and the comple plane
More informationMathematical aspects of mechanical systems eigentones
Seminar: Vibrations and Structure-Borne Sound in Civil Engineering Theor and Applications Mathematical aspects of mechanical sstems eigentones Andre Kuzmin April st 6 Abstract Computational methods of
More informationInterspecific Segregation and Phase Transition in a Lattice Ecosystem with Intraspecific Competition
Interspecific Segregation and Phase Transition in a Lattice Ecosstem with Intraspecific Competition K. Tainaka a, M. Kushida a, Y. Ito a and J. Yoshimura a,b,c a Department of Sstems Engineering, Shizuoka
More informationPHYSICAL REVIEW E 78, Lattice Boltzmann method for simulation of compressible flows on standard lattices
PHYSICAL REVIEW E 7, 016704 00 Lattice Boltzmann method for simulation of compressible flows on standard lattices Nikolaos I. Prasianakis 1, * and Ilia V. Karlin 1,, 1 Aerothermochemistr and Combustion
More informationPeriodic Structures in FDTD
EE 5303 Electromagnetic Analsis Using Finite Difference Time Domain Lecture #19 Periodic Structures in FDTD Lecture 19 These notes ma contain coprighted material obtained under fair use rules. Distribution
More information