Determining of temperature field in a L-shaped domain
|
|
- Hector Short
- 5 years ago
- Views:
Transcription
1 Available online at Advances in Applied Science Research, 0, (:-8 Determining of temperatre field in a L-shaped domain Oigo M. Zongo, Sié Kam, Kalifa Palm, and Alione Oedraogo ISSN: CODEN (USA: AASRFC Département de Physiqe, UFR-SEA, Université de Oagadogo, Oagadogo, Brkina Faso Institt de Recherches en Sciences Appliqées et Technologies (IRSAT, Oagadogo, Brkina Faso _ ABSTRACT This paper deals with solving Dirichlet s problem with nonhomogenos bondary conditions, sing the method of large singlar finite elements for Laplace s eqation in a L-shaped domain. This method is particlarly sitable for solving singlar problems since the analytical form of singlarities has been integrated in the approximate soltion. It is niqe in that the eqations are exactly verified except on those areas where internal segments are approached with great precision. Reslts are compared with those obtained throgh finite elements method by sing the COMSOL software. They dealt with soltion vales and those of its first derivatives. Both methods give reslts that align qite well everywhere, except near singlarities where significant differences exist. Keywords: Singlar points, least sqares, large elements, finite elements. _ INTRODUCTION The general theory of soltions to Laplace's eqation is known as potential theory. The soltions of Laplace's eqation are the harmonic fnctions, which are important in many fields of science, notably the fields of electromagnetism [] and flid dynamics [], becase they can be sed to accrately describe the behavior of electric, gravitational, and flid potentials. In the stdy of heat condction [, ], the Laplace eqation is the steadystate heat eqation. Nmerical solving of Laplace s eqation still remains very difficlt and sal methods of finite elements and finite differences give nsatisfactory reslts when they are sed in their standard form. These methods as evidenced by varios athors [-] can be significantly improved if they take into consideration the analytical form of the soltion near of the singlarities. We therefore apply the method of the large singlar finite elements [] to calclate the temperatre field in polygon L-shaped with Dirichlet non-homogenos bondary conditions. This gives good reslts all over the stdy domain while the finite element method gives good reslts only on areas located far from singlarities. This shows the power, efficiency and accracy of the method of large singlar finite elements for a limited nmber of conserved coefficients than the finite elements method. MATERIALS AND METHODS Either to determine the stationary temperatre field in a L-shaped polygon comprising eight segments of nit length kept at varios constant temperatres. In redced variables, the problem consists in solving Laplace s eqation with Dirichlet bondary conditions with mps. The L-shaped domain [] of side and bondary conditions are given in figre. The problem as posed is singlar de to discontinities in limit conditions on borders and the six smmits of the domain where there is discontinity in the external normal in the border of the domain. This is solved sing method of large singlar finite elements.
2 Oigo M. Zongo et al Adv. Appl. Sci. Res., 0, (:-8 Figre -L-shaped domain. Dirichlet bondary conditions. Figre - L-shaped domain. Smmary of local coordinates. We solve the Laplace s eqation on a polygonal domain sing MATLAB software. Bondaries Data can combine conditions of Dirichlet or Nemann. The method of large singlar finite elements comprises three steps []:
3 Oigo M. Zongo et al Adv. Appl. Sci. Res., 0, (:-8 Step : Decomposition of the domain. Althogh the domain L-shaped is geometrically symmetrical, it is not possible to replace the problem by an eqivalent problem in a redced domain becase conditions on limits are not symmetrical. The polygon mst therefore be taken as a whole. The first step of the method consists in dividing the domain into eight sb-domains identified in Figre, separated by ten sb-borders. Step : Solving axiliary problems The nmber of problems mst be eqal to the nmber of sb-domainsω. Therefore, to each sb-domain Ω is associated with an originσ, which is a singlarity, α the opening angle and a local system with polar coordinates, ( i i r (figre. i i i i First Axiliary Problem r = 0 ( r ( r,0 = 0 ( Ω π ( r, = 0. Second Axiliary Problem ( r = 0 ( r Ω ( r,0 = 0. ( r, π = 0 Third Axiliary Problem ( r, = 0 ( r, Ω ( r,0 = 0. π ( r, = 0. Forth Axiliary Problem ( r = 0 ( r Ω ( r,0 = 0. π ( r, = 0. Fifth Axiliary Problem ( Ω r, = 0 ( r, ( r,0 = 0. π ( r, = 0. Sixth Axiliary Problem ( r = 0 ( r Ω ( r,0 = 0. π ( r, = 0. Seventh Axiliary Problem ( r, = 0 ( r, Ω (-a (-b (-c (-a (-b (-c (-a (-b (-c (-a (-b (-c (-a (-b (-c (-a (-b (-c (-a
4 Oigo M. Zongo et al Adv. Appl. Sci. Res., 0, (:-8 ( r,0 = 0 (-b π ( r, = 0. (-c Eighth Axiliary Problem 8 ( r8 8 = 0 ( r8 8 Ω8 (8-a ( r 8,0 = 0. (8-b 8 8 ( r 8, π = 0 SOLUTIONS AND DISCUSSION General soltions of the eight axiliary problems are written taking into accont varios bondary conditions specified in figre []: i ( r = a ir sin i (9-a π i= π = ( r = ( a r sin k ( r = ( ak r sin π k = k m ( r = ( amr sin π m= ( r = ( π n= a n m r m m sin p ( r = + a pr sin π p= q ( r = aqr sin π q= = q p 8 8 ( r8 8 = ( a8 r8 sin 8 π (8-c (0-a (-a (-a (-a (-a (-a (-a Step: Connecting axiliary soltions The third step of the method consists in connecting these axiliary soltions by imposing the continity of the fnction and its normal derivatives all along the varios sb-borders Γ separating two adacent sb-domains Ωk and Ω. If conditions on interfaces Γ are met, therefore the exact soltion to the initial problem will be fond. l Actally, since we cannot solve an infinite system, we mst generally keep to approximate soltions. The approximation reslt on the one hand, the fact that we mst limit the development (9-a to (-a to a finite nmber of terms, and secondly, that 'we mst content orselves, with rare exceptions, a connection imperfect. Approximate soltions have been obtained by limiting the series sed in general soltions. The nmber of terms kept in each of the sms is chosen according to the principle proposed by Desclox and Tolley [] aiming at representing the approximate soltions sing fnctions of degrees as niform as possible. This is achieved by keeping more coefficients for sb-domains with larger openings. Approximate axiliary soltions are all of N degrees (N being the nmber of coefficients retained for the sqare sb-domains which opening angle is π /. Finally, the total nmber of nknown parameters a, which vale can be freely chosen will be: ( + + N = N
5 Oigo M. Zongo et al Adv. Appl. Sci. Res., 0, (:-8 This allows getting the following approximate soltions: N i ( r = + a ir sin i (9-b π i= N ( r = ( + a r sin π = ( k N k r = ( + ak r sin π k = ( m N m r = ( + amr sin π m= n= m N m ( r = ( + anr sin π ( p N p r = + + a pr sin π p= ( q N q r = + aqr sin π q= N 8 8 ( r8 8 = ( + a8 r8 sin 8 π = (0-b (-b (-b (-b (-b (-b (-b We are therefore connecting the soltions of axiliary problems according to least sqares method, i.e. we calclate coefficients I ( a mn = a that allow minimizing the fnction: ( ( i a a l ik ( i ( aik ( a l ds ni n < + + Γ i i These coefficients are soltion of the linear algebraic system positive definite sqare matrix of order N with N nknowns classically called normal eqations of Gass: I ( amn = O a m, k =,,..., 8; n, l =,.... N (8 The accracy of approximate soltions is directly linked to the qality of the connection of soltions to axiliary problems. It is therefore natral to characterize its precision by measring the imperfections of continity conditions. We will se the overall connecting error definite by (9: = ( + + ds η (9 Where k l k l k l S < Γ k l υ υ ds is the element which arch length is Γ, S its length and k i υ and ( υ l the normals to the sb-border separating both adacent sb-areas. If the overall error is nll, the approximate soltion got aligns with the exact soltion. Mode of convergence of the method of large singlar finite elements is exponential. Indeed, the base 0 logarithm of the overall error decreases linearly with N as shown in Figre. Moreover, this crve of the overall error allows s to say, for example, that by keeping 88 coefficients a the nmerical reslts are obtained, with an absolte accracy of less than 0.
6 Oigo M. Zongo et al Adv. Appl. Sci. Res., 0, (:-8 Figre - Evoltion of the overall error according to the nmber of coefficients a retained. In order to compare reslts obtained sing the method of large singlar finite elements with those provided by the method of conventional finite elements in its standard form, the problem was solved sing COMSOL software. The comparison of both methods is smmarized in few graphics shown in figres, and. In figre, approximations of the fnction and its derivatives were presented all along circles centered on singlarity (figre. The fll lines have been obtained sing the method of large singlar finite elements by keeping 88 coefficients a and the small circles are the reslts from the finite elements method for a range of, degrees of freedom. Figre is similar to the previos one, bt only data on fnction for the for circles are there presented. Finally, the bi-logarithmic chart in figre gives the maximal gaps between reslts obtained sing the method of large singlar finite elements and the finite elements method all along varios circles centered inσ according to their radis. Reslts obtained by both methods are in line with the whole domain of calclation. However, if we examine the nmerical vales obtained near the singlar points (discontinity points of the fnction, it appears that the method of large singlar finite elements is significantly well accrate enogh. If gaps between both methods are arond 0 over a circle which radis is.0, it increases to 0 for a circle which radis is0, then to 0 if the radis of the circle is 0. σ
7 Oigo M. Zongo et al Adv. Appl. Sci. Res., 0, (:-8 Figre Comparing the vales of (ble, (red and (black obtained throgh the method of large x y singlar finite elements (fll lines and finite elements method (circles. 8
8 Oigo M. Zongo et al Adv. Appl. Sci. Res., 0, (:-8 Figre - Comparing the vales of obtained sing the method of large singlar finite elements (fll lines and finite elements method (circles. 9
9 Oigo M. Zongo et al Adv. Appl. Sci. Res., 0, (:-8 Figre - Maximal local gap between the vales of and its derivatives calclated sing the method of large singlar finite elements and finite elements method all along the varios circles centered on the singlarity σ CONCLUSION Reslts obtained by the method of large singlar finite elements and finite elements method are in line all over the domain. However, by checking the nmerical vales obtained near the singlar points, this shows that the method of large singlar finite elements is well accrate enogh. The method of large singlar finite elements takes the existence of singlarities into accont, while trying to find analytical soltions near them; this allows therefore getting withot frther formlation, derivative vales. Moreover, the mode of convergence of this method is exponential. Acknowledgements The athors wold like to thank CUD, Belgim, for the grant (CIUF, P, Physics given to Mr. Zongo and for the provision of MATLAB Software. REFERENCES [] 0. Marice, A. Reineix, P. Hoffmann, B. Bernard and P. Poligen, Adv. in Appl. Sci. Res. 0, (, 9. [] D. S. Chahan and V. Kmar, Adv. in Appl. Sci. Res. 0, (,. [] R. Saravana, S. Sreekanth, S. Sreenadh and S. Hemadri Reddy. Adv. in Appl. Sci. Res. 0, (,. [] V. Sri Hari Bab and G.V. Ramana Reddy, Adv. in Appl. Sci. Res. 0, (, 8. [] R.E Barnhill and J. R. Whiteman, Singlarities de to re-entrant bondaries in elliptic problem. in L. Collatz (ed., I.S.M.N. 9, Birkhaser-Verlag, Basel, 9. [] A. F. Emery, Trans. ASME (C, 9, 9,. [] G. J. Fix, Math. Mech., 99, 8,. [8] Z.C. Li and T.T. L, Math. Compt. Model., 000,, 9. [9] H. Motz, Qart. Appl. Math., 9,,. [0] G. Strang and G.Fix, Analysis of the finite element method. Prentice Hall, 9. [] R. Wait and A. R. Mitchell, J. Compt. Phys. 9, 8,. [] Z. C. Li, T. T. L, H. Y. H and A. H. D. Cheng. Eng. Anal. Bond. Elem., 00, 9, 9, 80
10 Oigo M. Zongo et al Adv. Appl. Sci. Res., 0, (:-8 [] L. Fox, P. Henrici and C. B. Moler, SIAM, J. Nmer. Anal. 9,, 99. [] M. D. Tolley, PhD Thesis, Université Libre de Brxelles, (Brxelles, Belgim, 9. [] J. Desclox and M. D. Tolley, Compt. Method. Appl. Mech. Eng, 98, 9 (,. 8
Chapter 2 Difficulties associated with corners
Chapter Difficlties associated with corners This chapter is aimed at resolving the problems revealed in Chapter, which are cased b corners and/or discontinos bondar conditions. The first section introdces
More informationLecture Notes: Finite Element Analysis, J.E. Akin, Rice University
9. TRUSS ANALYSIS... 1 9.1 PLANAR TRUSS... 1 9. SPACE TRUSS... 11 9.3 SUMMARY... 1 9.4 EXERCISES... 15 9. Trss analysis 9.1 Planar trss: The differential eqation for the eqilibrim of an elastic bar (above)
More informationDiscontinuous Fluctuation Distribution for Time-Dependent Problems
Discontinos Flctation Distribtion for Time-Dependent Problems Matthew Hbbard School of Compting, University of Leeds, Leeds, LS2 9JT, UK meh@comp.leeds.ac.k Introdction For some years now, the flctation
More informationSecond-Order Wave Equation
Second-Order Wave Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 3 December 016 1 Introdction The classical wave eqation is a second-order
More informationFormal Methods for Deriving Element Equations
Formal Methods for Deriving Element Eqations And the importance of Shape Fnctions Formal Methods In previos lectres we obtained a bar element s stiffness eqations sing the Direct Method to obtain eact
More informationComputational Geosciences 2 (1998) 1, 23-36
A STUDY OF THE MODELLING ERROR IN TWO OPERATOR SPLITTING ALGORITHMS FOR POROUS MEDIA FLOW K. BRUSDAL, H. K. DAHLE, K. HVISTENDAHL KARLSEN, T. MANNSETH Comptational Geosciences 2 (998), 23-36 Abstract.
More informationA Computational Study with Finite Element Method and Finite Difference Method for 2D Elliptic Partial Differential Equations
Applied Mathematics, 05, 6, 04-4 Pblished Online November 05 in SciRes. http://www.scirp.org/jornal/am http://d.doi.org/0.46/am.05.685 A Comptational Stdy with Finite Element Method and Finite Difference
More informationCurves - Foundation of Free-form Surfaces
Crves - Fondation of Free-form Srfaces Why Not Simply Use a Point Matrix to Represent a Crve? Storage isse and limited resoltion Comptation and transformation Difficlties in calclating the intersections
More information08.06 Shooting Method for Ordinary Differential Equations
8.6 Shooting Method for Ordinary Differential Eqations After reading this chapter, yo shold be able to 1. learn the shooting method algorithm to solve bondary vale problems, and. apply shooting method
More informationPulses on a Struck String
8.03 at ESG Spplemental Notes Plses on a Strck String These notes investigate specific eamples of transverse motion on a stretched string in cases where the string is at some time ndisplaced, bt with a
More informationVectors in Rn un. This definition of norm is an extension of the Pythagorean Theorem. Consider the vector u = (5, 8) in R 2
MATH 307 Vectors in Rn Dr. Neal, WKU Matrices of dimension 1 n can be thoght of as coordinates, or ectors, in n- dimensional space R n. We can perform special calclations on these ectors. In particlar,
More informationNumerical verification of the existence of localization of the elastic energy for closely spaced rigid disks
Nmerical verification of the existence of localization of the elastic energy for closely spaced rigid disks S. I. Rakin Siberian State University of transport Rssia, 6349, Novosibirsk, Dsy Kovalchk street,
More information4 Exact laminar boundary layer solutions
4 Eact laminar bondary layer soltions 4.1 Bondary layer on a flat plate (Blasis 1908 In Sec. 3, we derived the bondary layer eqations for 2D incompressible flow of constant viscosity past a weakly crved
More informationTechnical Note. ODiSI-B Sensor Strain Gage Factor Uncertainty
Technical Note EN-FY160 Revision November 30, 016 ODiSI-B Sensor Strain Gage Factor Uncertainty Abstract Lna has pdated or strain sensor calibration tool to spport NIST-traceable measrements, to compte
More informationPhysicsAndMathsTutor.com
C Integration - By sbstittion PhysicsAndMathsTtor.com. Using the sbstittion cos +, or otherwise, show that e cos + sin d e(e ) (Total marks). (a) Using the sbstittion cos, or otherwise, find the eact vale
More informationStudy of the diffusion operator by the SPH method
IOSR Jornal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-684,p-ISSN: 2320-334X, Volme, Isse 5 Ver. I (Sep- Oct. 204), PP 96-0 Stdy of the diffsion operator by the SPH method Abdelabbar.Nait
More informationA Survey of the Implementation of Numerical Schemes for Linear Advection Equation
Advances in Pre Mathematics, 4, 4, 467-479 Pblished Online Agst 4 in SciRes. http://www.scirp.org/jornal/apm http://dx.doi.org/.436/apm.4.485 A Srvey of the Implementation of Nmerical Schemes for Linear
More informationA New Approach to Direct Sequential Simulation that Accounts for the Proportional Effect: Direct Lognormal Simulation
A ew Approach to Direct eqential imlation that Acconts for the Proportional ffect: Direct ognormal imlation John Manchk, Oy eangthong and Clayton Detsch Department of Civil & nvironmental ngineering University
More informationNonlinear parametric optimization using cylindrical algebraic decomposition
Proceedings of the 44th IEEE Conference on Decision and Control, and the Eropean Control Conference 2005 Seville, Spain, December 12-15, 2005 TC08.5 Nonlinear parametric optimization sing cylindrical algebraic
More informationOPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIELD OF A POLYHEDRAL BODY WITH LINEARLY INCREASING DENSITY 1
OPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIEL OF A POLYHERAL BOY WITH LINEARLY INCREASING ENSITY 1 V. POHÁNKA2 Abstract The formla for the comptation of the gravity field of a polyhedral body
More informationThe spreading residue harmonic balance method for nonlinear vibration of an electrostatically actuated microbeam
J.L. Pan W.Y. Zh Nonlinear Sci. Lett. Vol.8 No. pp.- September The spreading reside harmonic balance method for nonlinear vibration of an electrostatically actated microbeam J. L. Pan W. Y. Zh * College
More informationWorst-case analysis of the LPT algorithm for single processor scheduling with time restrictions
OR Spectrm 06 38:53 540 DOI 0.007/s009-06-043-5 REGULAR ARTICLE Worst-case analysis of the LPT algorithm for single processor schedling with time restrictions Oliver ran Fan Chng Ron Graham Received: Janary
More informationSimilarity Solution for MHD Flow of Non-Newtonian Fluids
P P P P IJISET - International Jornal of Innovative Science, Engineering & Technology, Vol. Isse 6, Jne 06 ISSN (Online) 48 7968 Impact Factor (05) - 4. Similarity Soltion for MHD Flow of Non-Newtonian
More informationReduction of over-determined systems of differential equations
Redction of oer-determined systems of differential eqations Maim Zaytse 1) 1, ) and Vyachesla Akkerman 1) Nclear Safety Institte, Rssian Academy of Sciences, Moscow, 115191 Rssia ) Department of Mechanical
More informationTransient Approach to Radiative Heat Transfer Free Convection Flow with Ramped Wall Temperature
Jornal of Applied Flid Mechanics, Vol. 5, No., pp. 9-1, 1. Available online at www.jafmonline.net, ISSN 175-57, EISSN 175-645. Transient Approach to Radiative Heat Transfer Free Convection Flow with Ramped
More informationHomotopy Perturbation Method for Solving Linear Boundary Value Problems
International Jornal of Crrent Engineering and Technolog E-ISSN 2277 4106, P-ISSN 2347 5161 2016 INPRESSCO, All Rights Reserved Available at http://inpressco.com/categor/ijcet Research Article Homotop
More informationInvestigation of Haar Wavelet Collocation Method to Solve Ninth Order Boundary Value Problems
Global Jornal of Pre and Applied Mathematics. ISSN 0973-1768 Volme 13, Nmber 5 (017), pp. 1415-148 Research India Pblications http://www.ripblication.com Investigation of Haar Wavelet Collocation Method
More informationKragujevac J. Sci. 34 (2012) UDC 532.5: :537.63
5 Kragjevac J. Sci. 34 () 5-. UDC 53.5: 536.4:537.63 UNSTEADY MHD FLOW AND HEAT TRANSFER BETWEEN PARALLEL POROUS PLATES WITH EXPONENTIAL DECAYING PRESSURE GRADIENT Hazem A. Attia and Mostafa A. M. Abdeen
More information10.2 Solving Quadratic Equations by Completing the Square
. Solving Qadratic Eqations b Completing the Sqare Consider the eqation ( ) We can see clearl that the soltions are However, What if the eqation was given to s in standard form, that is 6 How wold we go
More informationSetting The K Value And Polarization Mode Of The Delta Undulator
LCLS-TN-4- Setting The Vale And Polarization Mode Of The Delta Undlator Zachary Wolf, Heinz-Dieter Nhn SLAC September 4, 04 Abstract This note provides the details for setting the longitdinal positions
More informationElements of Coordinate System Transformations
B Elements of Coordinate System Transformations Coordinate system transformation is a powerfl tool for solving many geometrical and kinematic problems that pertain to the design of gear ctting tools and
More information5. The Bernoulli Equation
5. The Bernolli Eqation [This material relates predominantly to modles ELP034, ELP035] 5. Work and Energy 5. Bernolli s Eqation 5.3 An example of the se of Bernolli s eqation 5.4 Pressre head, velocity
More informationSubcritical bifurcation to innitely many rotating waves. Arnd Scheel. Freie Universitat Berlin. Arnimallee Berlin, Germany
Sbcritical bifrcation to innitely many rotating waves Arnd Scheel Institt fr Mathematik I Freie Universitat Berlin Arnimallee 2-6 14195 Berlin, Germany 1 Abstract We consider the eqation 00 + 1 r 0 k2
More informationOptimal Control of a Heterogeneous Two Server System with Consideration for Power and Performance
Optimal Control of a Heterogeneos Two Server System with Consideration for Power and Performance by Jiazheng Li A thesis presented to the University of Waterloo in flfilment of the thesis reqirement for
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT 02-01 The Krylov accelerated SIMPLER method for incompressible flow C. Vik and A. Saghir ISSN 1389-6520 Reports of the Department of Applied Mathematical Analysis
More informationDepartment of Industrial Engineering Statistical Quality Control presented by Dr. Eng. Abed Schokry
Department of Indstrial Engineering Statistical Qality Control presented by Dr. Eng. Abed Schokry Department of Indstrial Engineering Statistical Qality Control C and U Chart presented by Dr. Eng. Abed
More informationSuyeon Shin* and Woonjae Hwang**
JOURNAL OF THE CHUNGCHEONG MATHEMATICAL SOCIETY Volme 5, No. 3, Agst THE NUMERICAL SOLUTION OF SHALLOW WATER EQUATION BY MOVING MESH METHODS Syeon Sin* and Woonjae Hwang** Abstract. Tis paper presents
More informationEfficiency Increase and Input Power Decrease of Converted Prototype Pump Performance
International Jornal of Flid Machinery and Systems DOI: http://dx.doi.org/10.593/ijfms.016.9.3.05 Vol. 9, No. 3, Jly-September 016 ISSN (Online): 188-9554 Original Paper Efficiency Increase and Inpt Power
More informationSection 7.4: Integration of Rational Functions by Partial Fractions
Section 7.4: Integration of Rational Fnctions by Partial Fractions This is abot as complicated as it gets. The Method of Partial Fractions Ecept for a few very special cases, crrently we have no way to
More informationMEASUREMENT OF TURBULENCE STATISTICS USING HOT WIRE ANEMOMETRY
MEASUREMENT OF TURBULENCE STATISTICS USING HOT WIRE ANEMOMETRY Mrgan Thangadrai +, Atl Kmar Son *, Mritynjay Singh +, Sbhendra *, Vinoth Kmar ++, Ram Pyare Singh +, Pradip K Chatterjee + + Thermal Engineering,
More informationMODELLING OF TURBULENT ENERGY FLUX IN CANONICAL SHOCK-TURBULENCE INTERACTION
MODELLING OF TURBULENT ENERGY FLUX IN CANONICAL SHOCK-TURBULENCE INTERACTION Rssell Qadros, Krishnend Sinha Department of Aerospace Engineering Indian Institte of Technology Bombay Mmbai, India 476 Johan
More informationGeometric Image Manipulation. Lecture #4 Wednesday, January 24, 2018
Geometric Image Maniplation Lectre 4 Wednesda, Janar 4, 08 Programming Assignment Image Maniplation: Contet To start with the obvios, an image is a D arra of piels Piel locations represent points on the
More information1. INTRODUCTION. A solution for the dark matter mystery based on Euclidean relativity. Frédéric LASSIAILLE 2009 Page 1 14/05/2010. Frédéric LASSIAILLE
Frédéric LASSIAILLE 2009 Page 1 14/05/2010 Frédéric LASSIAILLE email: lmimi2003@hotmail.com http://lmi.chez-alice.fr/anglais A soltion for the dark matter mystery based on Eclidean relativity The stdy
More informationInertial Instability of Arbitrarily Meandering Currents Governed by the Eccentrically Cyclogeostrophic Equation
Jornal of Oceanography, Vol. 59, pp. 163 to 17, 3 Inertial Instability of Arbitrarily Meandering Crrents Governed by the Eccentrically Cyclogeostrophic Eqation HIDEO KAWAI* 131-81 Shibagahara, Kse, Joyo,
More informationChapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS
Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS 3. System Modeling Mathematical Modeling In designing control systems we mst be able to model engineered system dynamics. The model of a dynamic system
More informationDevelopment of Second Order Plus Time Delay (SOPTD) Model from Orthonormal Basis Filter (OBF) Model
Development of Second Order Pls Time Delay (SOPTD) Model from Orthonormal Basis Filter (OBF) Model Lemma D. Tfa*, M. Ramasamy*, Sachin C. Patwardhan **, M. Shhaimi* *Chemical Engineering Department, Universiti
More informationDirect linearization method for nonlinear PDE s and the related kernel RBFs
Direct linearization method for nonlinear PDE s and the related kernel BFs W. Chen Department of Informatics, Uniersity of Oslo, P.O.Box 1080, Blindern, 0316 Oslo, Norway Email: wenc@ifi.io.no Abstract
More informationInfluence of the Non-Linearity of the Aerodynamic Coefficients on the Skewness of the Buffeting Drag Force. Vincent Denoël *, 1), Hervé Degée 1)
Inflence of the Non-Linearity of the Aerodynamic oefficients on the Skewness of the Bffeting rag Force Vincent enoël *, 1), Hervé egée 1) 1) epartment of Material mechanics and Strctres, University of
More information3 2D Elastostatic Problems in Cartesian Coordinates
D lastostatic Problems in Cartesian Coordinates Two dimensional elastostatic problems are discssed in this Chapter, that is, static problems of either plane stress or plane strain. Cartesian coordinates
More informationCalculations involving a single random variable (SRV)
Calclations involving a single random variable (SRV) Example of Bearing Capacity q φ = 0 µ σ c c = 100kN/m = 50kN/m ndrained shear strength parameters What is the relationship between the Factor of Safety
More informationarxiv: v1 [physics.flu-dyn] 4 Sep 2013
THE THREE-DIMENSIONAL JUMP CONDITIONS FOR THE STOKES EQUATIONS WITH DISCONTINUOUS VISCOSITY, SINGULAR FORCES, AND AN INCOMPRESSIBLE INTERFACE PRERNA GERA AND DAVID SALAC arxiv:1309.1728v1 physics.fl-dyn]
More informationFOUNTAIN codes [3], [4] provide an efficient solution
Inactivation Decoding of LT and Raptor Codes: Analysis and Code Design Francisco Lázaro, Stdent Member, IEEE, Gianligi Liva, Senior Member, IEEE, Gerhard Bach, Fellow, IEEE arxiv:176.5814v1 [cs.it 19 Jn
More informationApproximate Solution for the System of Non-linear Volterra Integral Equations of the Second Kind by using Block-by-block Method
Astralian Jornal of Basic and Applied Sciences, (1): 114-14, 008 ISSN 1991-8178 Approximate Soltion for the System of Non-linear Volterra Integral Eqations of the Second Kind by sing Block-by-block Method
More informationOptimal Control, Statistics and Path Planning
PERGAMON Mathematical and Compter Modelling 33 (21) 237 253 www.elsevier.nl/locate/mcm Optimal Control, Statistics and Path Planning C. F. Martin and Shan Sn Department of Mathematics and Statistics Texas
More informationcalled the potential flow, and function φ is called the velocity potential.
J. Szantr Lectre No. 3 Potential flows 1 If the flid flow is irrotational, i.e. everwhere or almost everwhere in the field of flow there is rot 0 it means that there eists a scalar fnction ϕ,, z), sch
More informationKrauskopf, B., Lee, CM., & Osinga, HM. (2008). Codimension-one tangency bifurcations of global Poincaré maps of four-dimensional vector fields.
Kraskopf, B, Lee,, & Osinga, H (28) odimension-one tangency bifrcations of global Poincaré maps of for-dimensional vector fields Early version, also known as pre-print Link to pblication record in Explore
More informationFormules relatives aux probabilités qui dépendent de très grands nombers
Formles relatives ax probabilités qi dépendent de très grands nombers M. Poisson Comptes rends II (836) pp. 603-63 In the most important applications of the theory of probabilities, the chances of events
More informationSubstructure Finite Element Model Updating of a Space Frame Structure by Minimization of Modal Dynamic Residual
1 Sbstrctre Finite Element Model Updating of a Space Frame Strctre by Minimization of Modal Dynamic esidal Dapeng Zh, Xinn Dong, Yang Wang, Member, IEEE Abstract This research investigates a sbstrctre
More informationUNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL
8th International DAAAM Baltic Conference "INDUSTRIAL ENGINEERING - 19-1 April 01, Tallinn, Estonia UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL Põdra, P. & Laaneots, R. Abstract: Strength analysis is a
More informationPhysicsAndMathsTutor.com
. Two smooth niform spheres S and T have eqal radii. The mass of S is 0. kg and the mass of T is 0.6 kg. The spheres are moving on a smooth horizontal plane and collide obliqely. Immediately before the
More informationA Single Species in One Spatial Dimension
Lectre 6 A Single Species in One Spatial Dimension Reading: Material similar to that in this section of the corse appears in Sections 1. and 13.5 of James D. Mrray (), Mathematical Biology I: An introction,
More informationSimplified Identification Scheme for Structures on a Flexible Base
Simplified Identification Scheme for Strctres on a Flexible Base L.M. Star California State University, Long Beach G. Mylonais University of Patras, Greece J.P. Stewart University of California, Los Angeles
More informationIntegration of Basic Functions. Session 7 : 9/23 1
Integration o Basic Fnctions Session 7 : 9/3 Antiderivation Integration Deinition: Taking the antiderivative, or integral, o some nction F(), reslts in the nction () i ()F() Pt simply: i yo take the integral
More informationNumerical Study on Bouncing and Separation Collision Between Two Droplets Considering the Collision-Induced Breakup
Jornal of Mechanical Science and Technology (007) 585~59 Jornal of Mechanical Science and Technology Nmerical Stdy on Boncing and Separation Collision Between Two Droplets Considering the Collision-Indced
More informationSolving a System of Equations
Solving a System of Eqations Objectives Understand how to solve a system of eqations with: - Gass Elimination Method - LU Decomposition Method - Gass-Seidel Method - Jacobi Method A system of linear algebraic
More informationFEA Solution Procedure
EA Soltion Procedre (demonstrated with a -D bar element problem) EA Procedre for Static Analysis. Prepare the E model a. discretize (mesh) the strctre b. prescribe loads c. prescribe spports. Perform calclations
More informationTwo-media boundary layer on a flat plate
Two-media bondary layer on a flat plate Nikolay Ilyich Klyev, Asgat Gatyatovich Gimadiev, Yriy Alekseevich Krykov Samara State University, Samara,, Rssia Samara State Aerospace University named after academician
More informationDesign and Data Acquisition for Thermal Conductivity Matric Suction Sensors
68 TRANSPORTATION RSARCH RCORD 1432 Design and Data Acqisition for Thermal Condctivity Matric Sction Sensors J. K.-M. GAN, D. G. FRDLUND, A. XING, AND W.-X. LI The principles behind sing the thermal condctivity
More informationFixed points for discrete logarithms
Fixed points for discrete logarithms Mariana Levin 1, Carl Pomerance 2, and and K. Sondararajan 3 1 Gradate Grop in Science and Mathematics Edcation University of California Berkeley, CA 94720, USA levin@berkeley.ed
More informationBLOOM S TAXONOMY. Following Bloom s Taxonomy to Assess Students
BLOOM S TAXONOMY Topic Following Bloom s Taonomy to Assess Stdents Smmary A handot for stdents to eplain Bloom s taonomy that is sed for item writing and test constrction to test stdents to see if they
More informationFREQUENCY DOMAIN FLUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS
7 TH INTERNATIONAL CONGRESS O THE AERONAUTICAL SCIENCES REQUENCY DOMAIN LUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS Yingsong G, Zhichn Yang Northwestern Polytechnical University, Xi an, P. R. China,
More informationPrandl established a universal velocity profile for flow parallel to the bed given by
EM 0--00 (Part VI) (g) The nderlayers shold be at least three thicknesses of the W 50 stone, bt never less than 0.3 m (Ahrens 98b). The thickness can be calclated sing Eqation VI-5-9 with a coefficient
More informationSolving a Class of PDEs by a Local Reproducing Kernel Method with An Adaptive Residual Subsampling Technique
Copyright 25 Tech Science Press CMES, vol.8, no.6, pp.375-396, 25 Solving a Class of PDEs by a Local Reprodcing Kernel Method with An Adaptive Residal Sbsampling Techniqe H. Rafieayan Zadeh, M. Mohammadi,2
More informationAnalytical Investigation of Hyperbolic Equations via He s Methods
American J. of Engineering and Applied Sciences (4): 399-47, 8 ISSN 94-7 8 Science Pblications Analytical Investigation of Hyperbolic Eqations via He s Methods D.D. Ganji, M. Amini and A. Kolahdooz Department
More informationSTATIC, STAGNATION, AND DYNAMIC PRESSURES
STATIC, STAGNATION, AND DYNAMIC PRESSURES Bernolli eqation is g constant In this eqation is called static ressre, becase it is the ressre that wold be measred by an instrment that is static with resect
More informationAn optimizing reduced order FDS for the tropical Pacific Ocean reduced gravity model
An optimizing redced order FDS for the tropical Pacific Ocean redced gravity model Zhendong Lo a,, Jing Chen b,, Jiang Zh c,, Riwen Wang c,, and I. M. Navon d,, a School of Science, Beijing Jiaotong University,
More informationConditions for Approaching the Origin without Intersecting the x-axis in the Liénard Plane
Filomat 3:2 (27), 376 377 https://doi.org/.2298/fil7276a Pblished by Faclty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Conditions for Approaching
More informationVNVe 2017/ Final project
VNVe 2017/2018 - Final project Athor s name Febrary 21, 2018 Solve all examples and send yor final soltion (pdf file) and all sorce codes (LaTex, MATLAB,, ++, etc.) to e-mail address satek@fit.vtbr.cz
More informationFrequency Estimation, Multiple Stationary Nonsinusoidal Resonances With Trend 1
Freqency Estimation, Mltiple Stationary Nonsinsoidal Resonances With Trend 1 G. Larry Bretthorst Department of Chemistry, Washington University, St. Lois MO 6313 Abstract. In this paper, we address the
More informationStep-Size Bounds Analysis of the Generalized Multidelay Adaptive Filter
WCE 007 Jly - 4 007 London UK Step-Size onds Analysis of the Generalized Mltidelay Adaptive Filter Jnghsi Lee and Hs Chang Hang Abstract In this paper we analyze the bonds of the fixed common step-size
More informationFinite Difference Method of Modelling Groundwater Flow
Jornal of Water Resorce and Protection, 20, 3, 92-98 doi:0.4236/warp.20.33025 Pblished Online March 20 (http://www.scirp.org/ornal/warp) Finite Difference Method of Modelling Grondwater Flow Abstract Magns.
More informationApplying Laminar and Turbulent Flow and measuring Velocity Profile Using MATLAB
IOS Jornal of Mathematics (IOS-JM) e-issn: 78-578, p-issn: 319-765X. Volme 13, Isse 6 Ver. II (Nov. - Dec. 17), PP 5-59 www.iosrjornals.org Applying Laminar and Trblent Flow and measring Velocity Profile
More informationPrediction of Effective Asphalt Layer Temperature
TRANSPORTATION RESEARCH RECORD 1473 93 Prediction of Effective Asphalt Layer Temperatre EARL H. INGE, JR., AND Y. RICHARD KIM The most widely sed method for evalating deflection measrements for overlay
More informationComplex Variables. For ECON 397 Macroeconometrics Steve Cunningham
Comple Variables For ECON 397 Macroeconometrics Steve Cnningham Open Disks or Neighborhoods Deinition. The set o all points which satis the ineqalit
More informationEXERCISES WAVE EQUATION. In Problems 1 and 2 solve the heat equation (1) subject to the given conditions. Assume a rod of length L.
.4 WAVE EQUATION 445 EXERCISES.3 In Problems and solve the heat eqation () sbject to the given conditions. Assme a rod of length.. (, t), (, t) (, ),, > >. (, t), (, t) (, ) ( ) 3. Find the temperatre
More informationRadiation Effects on Heat and Mass Transfer over an Exponentially Accelerated Infinite Vertical Plate with Chemical Reaction
Radiation Effects on Heat and Mass Transfer over an Exponentially Accelerated Infinite Vertical Plate with Chemical Reaction A. Ahmed, M. N.Sarki, M. Ahmad Abstract In this paper the stdy of nsteady flow
More informationWEAR PREDICTION OF A TOTAL KNEE PROSTHESIS TIBIAL TRAY
APPLIED PHYSICS MEDICAL WEAR PREDICTION OF A TOTAL KNEE PROSTHESIS TIBIAL TRAY L. CÃPITANU, A. IAROVICI, J. ONIªORU Institte of Solid Mechanics, Romanian Academy, Constantin Mille 5, Bcharest Received
More informationComputational Fluid Dynamics Simulation and Wind Tunnel Testing on Microlight Model
Comptational Flid Dynamics Simlation and Wind Tnnel Testing on Microlight Model Iskandar Shah Bin Ishak Department of Aeronatics and Atomotive, Universiti Teknologi Malaysia T.M. Kit Universiti Teknologi
More informationFRTN10 Exercise 12. Synthesis by Convex Optimization
FRTN Exercise 2. 2. We want to design a controller C for the stable SISO process P as shown in Figre 2. sing the Yola parametrization and convex optimization. To do this, the control loop mst first be
More informationSareban: Evaluation of Three Common Algorithms for Structure Active Control
Engineering, Technology & Applied Science Research Vol. 7, No. 3, 2017, 1638-1646 1638 Evalation of Three Common Algorithms for Strctre Active Control Mohammad Sareban Department of Civil Engineering Shahrood
More informationTIME ACCURATE FAST THREE-STEP WAVELET-GALERKIN METHOD FOR PARTIAL DIFFERENTIAL EQUATIONS
International Jornal of Wavelets, Mltiresoltion and Information Processing Vol. 4, No. (26) 65 79 c World Scientific Pblishing Company TIME ACCURATE FAST THREE-STEP WAVELET-GALERKIN METHOD FOR PARTIAL
More informationMathematical Analysis of Nipah Virus Infections Using Optimal Control Theory
Jornal of Applied Mathematics and Physics, 06, 4, 099- Pblished Online Jne 06 in SciRes. http://www.scirp.org/jornal/jamp http://dx.doi.org/0.436/jamp.06.464 Mathematical Analysis of Nipah Virs nfections
More informationVariability sustained pattern formation in subexcitable media
Variability sstained pattern formation in sbexcitable media Erik Glatt, Martin Gassel, and Friedemann Kaiser Institte of Applied Physics, Darmstadt University of Technology, 64289 Darmstadt, Germany (Dated:
More informationThe Linear Quadratic Regulator
10 The Linear Qadratic Reglator 10.1 Problem formlation This chapter concerns optimal control of dynamical systems. Most of this development concerns linear models with a particlarly simple notion of optimality.
More informationBurgers Equation. A. Salih. Department of Aerospace Engineering Indian Institute of Space Science and Technology, Thiruvananthapuram 18 February 2016
Brgers Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 18 Febrary 216 1 The Brgers Eqation Brgers eqation is obtained as a reslt of
More informationMAXIMUM AND ANTI-MAXIMUM PRINCIPLES FOR THE P-LAPLACIAN WITH A NONLINEAR BOUNDARY CONDITION. 1. Introduction. ν = λ u p 2 u.
2005-Ojda International Conference on Nonlinear Analysis. Electronic Jornal of Differential Eqations, Conference 14, 2006, pp. 95 107. ISSN: 1072-6691. URL: http://ejde.math.txstate.ed or http://ejde.math.nt.ed
More informationA Regulator for Continuous Sedimentation in Ideal Clarifier-Thickener Units
A Reglator for Continos Sedimentation in Ideal Clarifier-Thickener Units STEFAN DIEHL Centre for Mathematical Sciences, Lnd University, P.O. Box, SE- Lnd, Sweden e-mail: diehl@maths.lth.se) Abstract. The
More informationShock wave structure for Generalized Burnett Equations
Shock wave strctre for Generalized Brnett Eqations A.V. Bobylev, M. Bisi, M.P. Cassinari, G. Spiga Dept. of Mathematics, Karlstad University, SE-65 88 Karlstad, Sweden, aleander.bobylev@ka.se Dip. di Matematica,
More informationarxiv: v1 [physics.flu-dyn] 11 Mar 2011
arxiv:1103.45v1 [physics.fl-dyn 11 Mar 011 Interaction of a magnetic dipole with a slowly moving electrically condcting plate Evgeny V. Votyakov Comptational Science Laboratory UCY-CompSci, Department
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 information