NORMALE. A modied structured central scheme for. 2D hyperbolic conservation laws. Theodoros KATSAOUNIS. Doron LEVY
|
|
- Lambert Carter
- 6 years ago
- Views:
Transcription
1 E COLE NORMALE SUPERIEURE A moded structured central scheme for 2D hyperbolc conservaton laws Theodoros KATSAOUNIS Doron LEVY LMENS Département de Mathématques et Informatque CNRS URA 762
2 A moded structured central scheme for 2D hyperbolc conservaton laws Theodoros KATSAOUNIS Doron LEVY LMENS June 1998 Laboratore de Mathematques de l'ecole Normale Supereure 45 rue d'ulm PARIS Cedex 05 Tel : (33)(1) Adresse electronque : katsaoun@dm.ens.fr, dlevy@dm.ens.fr
3 A Moded Structured Central Scheme for 2D Hyperbolc Conservaton Laws Theodoros Katsaouns y Doron Levy y Abstract We present a new central scheme for approxmatng solutons of two-dmensonal systems of hyperbolc conservaton laws. Ths method s based on a modcaton of the staggered grd proposed n [5] whch prevents the crossngs of dscontnutes n the normal drecton, whle retanng the smplcty of the central framework. Our method satses a local maxmum prncple whch s based on a more compact stencl. Unlke the prevous method, t enables a natural extenson to adaptve methods on structured grds. Key words. Hyperbolc conservaton laws, central derence schemes, non-oscllatory schemes, adaptve methods. 1 Introducton In [5] Jang and Tadmor presented a second-order two-dmensonal central scheme for approxmatng solutons of systems of hyperbolc conservaton laws whch extends the one-dmensonal Nessyahu-Tadmor (NT) scheme [7]. Smlar approach was taken by Armnjon, et al. n [1], [3]. Followng the central framework whose prototype s the Lax and Fredrchs scheme [4], a Godunov-type scheme was constructed. Frst, a pecewse-lnear MUSCL-type [6] nterpolant was reconstructed from the gven cell-averages. Spurous oscllatons n the reconstructon were avoded by mplementng a non-lnear lmtng mechansm [8]. Ths nterpolant was then evolved exactly n tme and nally projected on ts staggered cell-averages. Due to the staggerng, there was no need to solve two-dmensonal Remann problems. Unfortunately, the staggerng was not sucent to elmnate the dscontnutes from the problem, and one was actually left wth one-dmensonal Remann problems n the normal drecton. In the method proposed n [5], nstead of explctly solvng these 1D Remann problems, the values around the dscontnutes were averaged. The dsspatve treatment of those dscontnutes n the normal drecton resulted n several numercal consequences such as smearng of the dscontnutes as evdent n the numercal results presented n [5]. In ths work we present a new central scheme whch was desgned to avod the crossngs of dscontnutes n the normal drecton. Our goal s obtaned by exchangng the orgnal staggered mesh wth an alternatve rotated and stretched mesh. All that follows, s a drect mplementaton of the prevously desgned methods wth our new meshes. Ths new structured mesh can be vewed as a degenerate verson of the unstructured mesh used n [2]. Along wth the man advantage of our new method, several byproducts are n hand. Frst, the method n [5] when vewed n every two tme steps conssted of a 9-ponts stencl. Our method, however, s based on a more compact 5-ponts stencl. Moreover, unlke the method n [5], our y Departement de Mathematques et d'informatque, Ecole Normale Supereure, 45 rue d'ulm, Pars Cedex 05, France; dlevy@dm.ens.fr, katsaoun@dm.ens.fr 1
4 new method can be easly extended to adaptve central schemes whch are based on structured staggered grds. Such an extenson seems to be hghly non-trval n the prevous framework. Ths paper s organzed as follows: we start n x2 by presentng our new central method. The smplcty compared wth an adaptve unstructured framework s emphaszed by the explct formulaton of the method outlned below. We then proceed n x3 to formulate and prove a maxmum prncple on the scheme. We end n x4 wth several numercal examples. Acknowledgment: Research was supported by Hyperbolc systems of conservaton laws TMR grant #ERBFMRXCT We would lke to thank Prof. E. Tadmor for hs remarks. 2 The 2D Method We consder the two-dmensonal system of conservaton laws v t + f(v) x + g(v) y = 0; (2.1) augmented wth the ntal data, v 0 (x; y) = v(x; y; t = 0). To approxmate solutons of (2.1), we rst ntroduce a unform rectangular mesh n the (x; y) plane, wth spacngs taken as x; y. On top of ths mesh, we then buld a staggered mesh, whose cells are of the shape of damonds, consult Fgure 2.1 (a). x y (a) (b) Fgure 2.1: (a) The Staggered Mesh. Sold lnes - rectangular grd. Dotted lnes - damond-shaped grd. (b) A Structured Adaptve Grd. By and ~ we denote the rectangular S cells and the damond cells, respectvely. The resultng meshes are therefore abbrevated as = and ~ S = ~. For smplcty of notatons we use a xed tme-step t, and denote the dscrete tme by t n = nt. By u n we denote an approxmaton to the cell-average v n n cell and at tme t n. In the rst phase of the staggerng we assume that u n ~ are gven and we wsh to compute u n+1. In the second phase, the roles of ~ and are exchanged. The reconstructon of our Godunov-type second-order method starts by reconstructng a pecewselnear nterpolant from the gven cell-averages u(x; y; t n ) = X ~ P (x; y; t n ) ~ : (2.2) Here, ~ denoted the characterstc functon of the cell ~ whle P s a lnear polynomal ; P (x; y; t n ) = u n ~ + u 0 ~ x? xc x + u 8 ~ y? yc y 2
5 where the center of ~ s denoted by (x c ; y c ), and u 0 ~, u 8 ~, are the dscrete slopes n the x- and y-drectons respectvely, u 0 ~ x u x (x c ; y c ; t n ) + O(x) 2 ; u 8 ~ y u y (x c ; y c ; t n ) + O(y) 2. The reconstructon of the slopes utlzes nonlnear lmters descrbed n the remarks below. N N W E NW NE (a) S (b) W SW SE E (c) S Fgure 2.2: Reconstructon: (a) squares from damonds, (b)-(c) damonds from squares. An exact evoluton of (2.2) followed by a projecton on ts staggered cell-averages results wth (consult Fgure 2.2 (a)) { u n+1 = u n? 1 j j Z t n+1 t n Z Z f x + The rst term of the RHS of (2.3), I 1, equals I 1 = u n = un N + u n E + u n S + u n W Z Z g y dxdyd := I 1 + I 2 : (2.3) x(u 0 W? u0 E) + y(u 8 S? u8 N) The ntegrals over the uxes n the second term of the RHS of (2.3), I 2, are approxmated by a second-order md-pont quadrature rule, an approxmaton whch s vald as long as t s done n a smooth regon. The smoothness requrement results wth a CFL condton on the stablty of the method whch due to geometrc consderaton equals 1=4 (compared wth 1=2 n [5]). Ths quadrature can be then explctly wrtten as I 2 =? f(u n+1=2 E )? f(u n+1=2 W )? g(u n+1=2 N )? g(u n+1=2 S ) : ; (2.4) where = t=x and = t=y are the usual xed mesh ratos. The values at tme t n+1=2 requred n (2.4) are predcted usng a rst-order Taylor expanson (whch s sucent for overall second-order accuracy due to the predctor-corrector structure of the method). Hence, for example, = u n E? f(u 2 E) 0? g(u 2 E) 8, and analogously for the other ntermedate values. In the second phase of the staggerng we assume that the values, u n are known and we wsh to compute u n+1 ~. In order to smplfy the notatons we agan advance n tme from tme to tn t n+1. The computaton s analogous to the rst phase, only ths tme due to the lack of symmetry between the two phases of the staggerng, the resultng formulas are derent. Here, there are two possbltes whch are schematcally drawn P n Fgure 2.2 (b)-(c). For both cases an exact evoluton n tme of our reconstructon, u(x; y; t n ) = P (x; y; t n ), yelds (compare wth (2.3)), u n+1=2 E u n+1 ~ = u n ~? 1 j~ j Z t n+1 t n Z Z ~ f x + Z Z ~ g y dxdyd := I 1 + I 2 : (2.5) 3
6 We start wth the case descrbed n Fgure (b). Here, I 1 = u n ~ = un W +un E + x 2 6 (u0 W? u0 E), and X I 2 =?p f(u c ~ ) ~n x + g(u n+1=2 c ~ ) ~n y : (2.6) j u c j ~ = u c j ~? t j ~ The md-values, u c ~, requred n (2.6) are once agan predcted by Taylor expanson j hf u (u n@cj ) u 0 + g ~ u (u n@cj ) u 8@cj ; ~ ~ c ~ j represents the center of the edges of the damond cells, j = fne,se,nw,swg (see Fgure 2.2 (b)). Here, Du NE = Du SE = Du E and Du NW = Du SW = Du W, wth D denotng the dscrete dervatve ether n the x- or n the y-drecton. The pont-values, u c ~, are also j computed by a Taylor expanson, e.g., u n NE = u n E? x 4 u0 E + y 4 E. u8 Analogous computatons hold for the last case descrbed n Fgure 2.2 (c). For the sake of brevty we lst only the rst term on the RHS of (2.5), I 1, whch n ths case equals I 1 = u n ~ = u n N +un S 2 + y 6 (u8 S? u8 N). Remarks: 1. Reconstructon of the Dervatves. A reconstructon of the dervatves wthout creatng spurous oscllatons requres non-lnear bult-n lmters. One can use, for example, for the x-dervatve n a rectangular cell (j; k), a lmtng on the rght/centered/left dervatves u 0 = MM (u n? j;k j+1;k j;k); 1 un 2 (un? j+1;k j?1;k); (u un n? j;k j?1;k) ; (2.7) un wth 1 2, and MMfu1; u2; : : :g = 8 < : mn fu g; f v > 0; 8; max fu g; f v < 0; 8; 0; otherwse: The choce of = 1 agrees wth the classcal Mn-Mod lmter (see [8] and [5] for more detals). Analogous expresson holds for the y-dervatve. The same routne s repeated for the dervatves n the damond cells. Only ths tme, snce the cells are not algned wth the axes, one has to lmt the dervatves after projectng them on the x and y drectons. For systems, the dervatves are computed component-wse,.e., f 0 j;k are gven by (2.7), consult [7, 5]. = f u(u j;k )u 0 j;k, where u0 j;k 2. Adaptve Mesh. A possble extenson of the method to adaptve method on unstructured meshes s demonstrated n Fgure 2.1 (b). An equvalent extenson wth the prevous method n [5] seems to be mpossble, at least wthout dealng wth complcated cases at the boundares (corners, etc...) and wth uneven dvsons of the cells. Our method formulated on the new mesh requres no specal corner treatment. We consder ths to be the great advantage of our method over the other avalable structured 2D methods. Moreover, snce there s no upwndng nvolved, none of the reected-waves problems whch are typcal to upwndng methods on adaptve structured meshes should appear. 4
7 3. Ecent Implementaton. We note that the smplcty of the scheme can be drectly projected onto ts mplementaton. It s unnecessary to use the standard methods of the unstructured framework n order to mplement our method. The smplest data structure to store the values n the damonds would be to dvde them nto trangles and to store the values n the trangles n a two-dmensonal array, such that every pont n the array corresponds to the rectangle whch these trangles belong to. In fact, the values of only two trangles out of four n each rectangle should be stored as the values of the other two can be retreved from the neghborng cells. 3 A Maxmum Prncple for Scalar Approxmatons An equvalent maxmum prncple to Theorem 1 n [5] mples to our new scheme. Snce our scheme s non-symmetrc between the two phases of the staggerng, t s natural to formulate the theorem n a non-staggered verson by consderng two jont tme-steps. Ths results wth a local bound on the cell-averages based on values taken from a 5-pont stencl. An equvalent two tme-steps formulaton of Theorem 1 n [5] would be based on a 9-pont stencl. Theorem 3.1 Consder the two-dmensonal scalar scheme (2.3),(2.5). Assume that the dscrete slopes satsfy the lmter property (2.7). Then for any 1 < 2 there exsts a sucently small CFL number C, such that f the followng CFL condton s fullled max( max u jf u (u)j; max jg u (u)j) C ; u then the followng local maxmum prncple holds mnfu n j?1;k; u n j;k; u n j+1;k; u n j;k?1; u n j;k+1 g un+2 j;k maxfun j?1;k; u n j;k; u n j+1;k; u n j;k?1; u n j;k+1 g: The proof of Theorem 3.1 s analogous to the proof of [5, Theorem 1] and we omt t for brevty. The key observaton for the proof s that every new staggered cell average can be wrtten as a convex combnaton of sums and derences of the cell-averages n the supportng cells. 4 Numercal Examples In Table 4.1 we present the L 1 and L 1 errors and convergence rate estmates for the lnear oblque advecton v t +v x +v y = 0 subject to v 0 = sn((x+y)). Equal spacngs were used, x = y = 1=N. The CFL was taken as 0:2 and the tme T = 0:5. We used the MM lmter wth = 1. These results are ndeed comparable wth those presented n [5, Table 4.1]. N L 1 error L 1 order L 1 error L 1 order Table 4.1: Lnear Oblque Advecton. L 1 and L 1 errors and convergence rates. 5
8 y y y (a) x (b) x (c) x Fgure 4.1: Perodc 2D-Burgers. N = 160; (a) - JT, CFL=0.4; (b) - JT, CFL=0.2; (c) - The new method, CFL=0.2. We end by presentng an example demonstratng the non-oscllatory behavor of our method. In Fgure 4.1 we show the results obtaned for the perodc two-dmensonal Burgers' equaton, v t + vv x + vv y = 0 n [?1; 1] [?1; 1], for tme T = 0:5, subject to the ntal condtons, v 0 (x; y) =?1:0; x < 0; y < 0;?0:2; x > 0; y < 0; v 0 (x; y) = 0:8; x < 0; y > 0; 0:5; x > 0; y > 0: The label JT refers to the method of [5]. The contour plots n Fgure 4.1 are zoomed nto [0; 1] [?1; 0]. Clearly, our method handles better 'dagonal' waves compared wth the method of [5], whle less smearng the dscontnutes (compare the results of both methods for the same CFL). References [1] Armnjon P., Stanescu D., Vallon M.-C., A Two-Dmensonal Fnte Volume Extenson of the Lax-Fredrchs and Nessyahu-Tadmor Schemes for Compressble Flow, Proc. 6th. Int. Symp. on CFD, Lake Tahoe, 1995, M. Hafez and K. Oshma, edtors, Vol. IV, pp [2] Armnjon P., Vallon M.-C. and Madrane A., A Fnte Volume Extenson of the Lax-Fredrchs and Nessyahu-Tadmor Schemes for Conservaton Laws on Unstructured Grds, IJCFD, 9, (1997), pp
9 [3] Armnjon P., Vallon M.-C., Generalsaton du Schema de Nessyahu-Tadmor pour Une Equaton Hyperbolque a Deux Dmensons D'espace, C.R. Acad. Sc. Pars, t. 320, sere I. (1995), pp [4] Fredrchs K. O., Lax P. D., Systems of Conservaton Equatons wth a Convex Extenson, Proc. Nat. Acad. Sc., 68, (1971), pp [5] Jang G.-S., Tadmor E., Nonoscllatory Central Schemes for Multdmensonal Hyperbolc Conservaton Laws, SIAM J. Sc. Comp., to appear. [6] van Leer B., Towards the Ultmate Conservatve Derence Scheme, V. A Second-Order Sequel to Godunov's Method, JCP, 32, (1979), pp [7] Nessyahu H., Tadmor E., Non-oscllatory Central Derencng for Hyperbolc Conservaton Laws, JCP, 87, no. 2 (1990), pp [8] Sweby P. K., Hgh Resoluton Schemes Usng Flux Lmters for Hyperbolc Conservaton Laws, SINUM, 21, no. 5 (1984), pp
Very simple computational domains can be discretized using boundary-fitted structured meshes (also called grids)
Structured meshes Very smple computatonal domans can be dscretzed usng boundary-ftted structured meshes (also called grds) The grd lnes of a Cartesan mesh are parallel to one another Structured meshes
More informationA HIGH-ORDER SPECTRAL (FINITE) VOLUME METHOD FOR CONSERVATION LAWS ON UNSTRUCTURED GRIDS
AIAA-00-058 A HIGH-ORDER SPECTRAL (FIITE) VOLUME METHOD FOR COSERVATIO LAWS O USTRUCTURED GRIDS Z.J. Wang Department of Mechancal Engneerng Mchgan State Unversty, East Lansng, MI 88 Yen Lu * MS T7B-, ASA
More informationHermite Splines in Lie Groups as Products of Geodesics
Hermte Splnes n Le Groups as Products of Geodescs Ethan Eade Updated May 28, 2017 1 Introducton 1.1 Goal Ths document defnes a curve n the Le group G parametrzed by tme and by structural parameters n the
More informationMathematics 256 a course in differential equations for engineering students
Mathematcs 56 a course n dfferental equatons for engneerng students Chapter 5. More effcent methods of numercal soluton Euler s method s qute neffcent. Because the error s essentally proportonal to the
More informationGSLM Operations Research II Fall 13/14
GSLM 58 Operatons Research II Fall /4 6. Separable Programmng Consder a general NLP mn f(x) s.t. g j (x) b j j =. m. Defnton 6.. The NLP s a separable program f ts objectve functon and all constrants are
More informationAn Accurate Evaluation of Integrals in Convex and Non convex Polygonal Domain by Twelve Node Quadrilateral Finite Element Method
Internatonal Journal of Computatonal and Appled Mathematcs. ISSN 89-4966 Volume, Number (07), pp. 33-4 Research Inda Publcatons http://www.rpublcaton.com An Accurate Evaluaton of Integrals n Convex and
More information2x x l. Module 3: Element Properties Lecture 4: Lagrange and Serendipity Elements
Module 3: Element Propertes Lecture : Lagrange and Serendpty Elements 5 In last lecture note, the nterpolaton functons are derved on the bass of assumed polynomal from Pascal s trangle for the fled varable.
More informationSupport Vector Machines
/9/207 MIST.6060 Busness Intellgence and Data Mnng What are Support Vector Machnes? Support Vector Machnes Support Vector Machnes (SVMs) are supervsed learnng technques that analyze data and recognze patterns.
More informationR s s f. m y s. SPH3UW Unit 7.3 Spherical Concave Mirrors Page 1 of 12. Notes
SPH3UW Unt 7.3 Sphercal Concave Mrrors Page 1 of 1 Notes Physcs Tool box Concave Mrror If the reflectng surface takes place on the nner surface of the sphercal shape so that the centre of the mrror bulges
More informationChapter 6 Programmng the fnte element method Inow turn to the man subject of ths book: The mplementaton of the fnte element algorthm n computer programs. In order to make my dscusson as straghtforward
More informationS1 Note. Basis functions.
S1 Note. Bass functons. Contents Types of bass functons...1 The Fourer bass...2 B-splne bass...3 Power and type I error rates wth dfferent numbers of bass functons...4 Table S1. Smulaton results of type
More informationParallel Numerics. 1 Preconditioning & Iterative Solvers (From 2016)
Technsche Unverstät München WSe 6/7 Insttut für Informatk Prof. Dr. Thomas Huckle Dpl.-Math. Benjamn Uekermann Parallel Numercs Exercse : Prevous Exam Questons Precondtonng & Iteratve Solvers (From 6)
More informationAnalysis of Continuous Beams in General
Analyss of Contnuous Beams n General Contnuous beams consdered here are prsmatc, rgdly connected to each beam segment and supported at varous ponts along the beam. onts are selected at ponts of support,
More informationWavefront Reconstructor
A Dstrbuted Smplex B-Splne Based Wavefront Reconstructor Coen de Vsser and Mchel Verhaegen 14-12-201212 2012 Delft Unversty of Technology Contents Introducton Wavefront reconstructon usng Smplex B-Splnes
More informationModule 6: FEM for Plates and Shells Lecture 6: Finite Element Analysis of Shell
Module 6: FEM for Plates and Shells Lecture 6: Fnte Element Analyss of Shell 3 6.6. Introducton A shell s a curved surface, whch by vrtue of ther shape can wthstand both membrane and bendng forces. A shell
More informationS.P.H. : A SOLUTION TO AVOID USING EROSION CRITERION?
S.P.H. : A SOLUTION TO AVOID USING EROSION CRITERION? Célne GALLET ENSICA 1 place Emle Bloun 31056 TOULOUSE CEDEX e-mal :cgallet@ensca.fr Jean Luc LACOME DYNALIS Immeuble AEROPOLE - Bat 1 5, Avenue Albert
More informationQuasi-ENO Schemes for Unstructured Meshes. Carl F. Ollivier-Gooch. Mathematics and Computer Science Division. Argonne National Laboratory
Quas-ENO Schemes for Unstructured Meshes Based on Unlmted Data-Dependent Least-Squares Reconstructon Carl F. Ollver-Gooch Mathematcs and Computer Scence Dvson Argonne Natonal Laboratory Subject classcaton:
More informationHigh-Boost Mesh Filtering for 3-D Shape Enhancement
Hgh-Boost Mesh Flterng for 3-D Shape Enhancement Hrokazu Yagou Λ Alexander Belyaev y Damng We z Λ y z ; ; Shape Modelng Laboratory, Unversty of Azu, Azu-Wakamatsu 965-8580 Japan y Computer Graphcs Group,
More informationOrder of Accuracy Study of Unstructured Grid Finite Volume Upwind Schemes
João Luz F. Azevedo et al. João Luz F. Azevedo joaoluz.azevedo@gmal.com Comando-Geral de Tecnologa Aeroespacal Insttuto de Aeronáutca e Espaço IAE 12228-903 São José dos Campos, SP, Brazl Luís F. Fguera
More informationSolitary and Traveling Wave Solutions to a Model. of Long Range Diffusion Involving Flux with. Stability Analysis
Internatonal Mathematcal Forum, Vol. 6,, no. 7, 8 Soltary and Travelng Wave Solutons to a Model of Long Range ffuson Involvng Flux wth Stablty Analyss Manar A. Al-Qudah Math epartment, Rabgh Faculty of
More informationCMPS 10 Introduction to Computer Science Lecture Notes
CPS 0 Introducton to Computer Scence Lecture Notes Chapter : Algorthm Desgn How should we present algorthms? Natural languages lke Englsh, Spansh, or French whch are rch n nterpretaton and meanng are not
More informationFor instance, ; the five basic number-sets are increasingly more n A B & B A A = B (1)
Secton 1.2 Subsets and the Boolean operatons on sets If every element of the set A s an element of the set B, we say that A s a subset of B, or that A s contaned n B, or that B contans A, and we wrte A
More informationLecture 4: Principal components
/3/6 Lecture 4: Prncpal components 3..6 Multvarate lnear regresson MLR s optmal for the estmaton data...but poor for handlng collnear data Covarance matrx s not nvertble (large condton number) Robustness
More information6.854 Advanced Algorithms Petar Maymounkov Problem Set 11 (November 23, 2005) With: Benjamin Rossman, Oren Weimann, and Pouya Kheradpour
6.854 Advanced Algorthms Petar Maymounkov Problem Set 11 (November 23, 2005) Wth: Benjamn Rossman, Oren Wemann, and Pouya Kheradpour Problem 1. We reduce vertex cover to MAX-SAT wth weghts, such that the
More informationAn Optimal Algorithm for Prufer Codes *
J. Software Engneerng & Applcatons, 2009, 2: 111-115 do:10.4236/jsea.2009.22016 Publshed Onlne July 2009 (www.scrp.org/journal/jsea) An Optmal Algorthm for Prufer Codes * Xaodong Wang 1, 2, Le Wang 3,
More informationRadial Basis Functions
Radal Bass Functons Mesh Reconstructon Input: pont cloud Output: water-tght manfold mesh Explct Connectvty estmaton Implct Sgned dstance functon estmaton Image from: Reconstructon and Representaton of
More informationAMath 483/583 Lecture 21 May 13, Notes: Notes: Jacobi iteration. Notes: Jacobi with OpenMP coarse grain
AMath 483/583 Lecture 21 May 13, 2011 Today: OpenMP and MPI versons of Jacob teraton Gauss-Sedel and SOR teratve methods Next week: More MPI Debuggng and totalvew GPU computng Read: Class notes and references
More informationReading. 14. Subdivision curves. Recommended:
eadng ecommended: Stollntz, Deose, and Salesn. Wavelets for Computer Graphcs: heory and Applcatons, 996, secton 6.-6., A.5. 4. Subdvson curves Note: there s an error n Stollntz, et al., secton A.5. Equaton
More informationAn Application of the Dulmage-Mendelsohn Decomposition to Sparse Null Space Bases of Full Row Rank Matrices
Internatonal Mathematcal Forum, Vol 7, 2012, no 52, 2549-2554 An Applcaton of the Dulmage-Mendelsohn Decomposton to Sparse Null Space Bases of Full Row Rank Matrces Mostafa Khorramzadeh Department of Mathematcal
More informationOn the diameter of random planar graphs
On the dameter of random planar graphs Gullaume Chapuy, CNRS & LIAFA, Pars jont work wth Érc Fusy, Pars, Omer Gménez, ex-barcelona, Marc Noy, Barcelona. Probablty and Graphs, Eurandom, Endhoven, 04. Planar
More informationNAG Fortran Library Chapter Introduction. G10 Smoothing in Statistics
Introducton G10 NAG Fortran Lbrary Chapter Introducton G10 Smoothng n Statstcs Contents 1 Scope of the Chapter... 2 2 Background to the Problems... 2 2.1 Smoothng Methods... 2 2.2 Smoothng Splnes and Regresson
More information3D vector computer graphics
3D vector computer graphcs Paolo Varagnolo: freelance engneer Padova Aprl 2016 Prvate Practce ----------------------------------- 1. Introducton Vector 3D model representaton n computer graphcs requres
More informationInterpolation of the Irregular Curve Network of Ship Hull Form Using Subdivision Surfaces
7 Interpolaton of the Irregular Curve Network of Shp Hull Form Usng Subdvson Surfaces Kyu-Yeul Lee, Doo-Yeoun Cho and Tae-Wan Km Seoul Natonal Unversty, kylee@snu.ac.kr,whendus@snu.ac.kr,taewan}@snu.ac.kr
More informationLecture 5: Multilayer Perceptrons
Lecture 5: Multlayer Perceptrons Roger Grosse 1 Introducton So far, we ve only talked about lnear models: lnear regresson and lnear bnary classfers. We noted that there are functons that can t be represented
More informationFeature Reduction and Selection
Feature Reducton and Selecton Dr. Shuang LIANG School of Software Engneerng TongJ Unversty Fall, 2012 Today s Topcs Introducton Problems of Dmensonalty Feature Reducton Statstc methods Prncpal Components
More informationDiscontinuous Galerkin methods for flow and transport problems in porous media
T COMMUNICATIONS IN NUMERICA METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2; :1 6 [Verson: 2/3/22 v1.] Dscontnuous Galerkn methods for flow and transport problems n porous meda Béatrve Rvère and Mary
More informationApplication of Residual Distribution (RD) schemes to the geographical part of the Wave Action Equation
Applcaton of Resdual Dstrbuton (RD) schemes to the geographcal part of the Wave Acton Equaton Aron Roland Insttute for Hydraulc and Water Resources Engneerng, Technsche Unverstät Darmstadt, Germany Abstract
More informationRange images. Range image registration. Examples of sampling patterns. Range images and range surfaces
Range mages For many structured lght scanners, the range data forms a hghly regular pattern known as a range mage. he samplng pattern s determned by the specfc scanner. Range mage regstraton 1 Examples
More informationForm-factors Josef Pelikán CGG MFF UK Praha.
Form-factors 1996-2016 Josef Pelkán CGG MFF UK Praha pepca@cgg.mff.cun.cz http://cgg.mff.cun.cz/~pepca/ FormFactor 2016 Josef Pelkán, http://cgg.mff.cun.cz/~pepca 1 / 23 Form-factor F It ndcates the proporton
More informationSimulation Based Analysis of FAST TCP using OMNET++
Smulaton Based Analyss of FAST TCP usng OMNET++ Umar ul Hassan 04030038@lums.edu.pk Md Term Report CS678 Topcs n Internet Research Sprng, 2006 Introducton Internet traffc s doublng roughly every 3 months
More informationRECENT research on structured mesh flow solver for aerodynamic problems shows that for practical levels of
A Hgh-Order Accurate Unstructured GMRES Algorthm for Invscd Compressble Flows A. ejat * and C. Ollver-Gooch Department of Mechancal Engneerng, The Unversty of Brtsh Columba, 054-650 Appled Scence Lane,
More informationReview of approximation techniques
CHAPTER 2 Revew of appromaton technques 2. Introducton Optmzaton problems n engneerng desgn are characterzed by the followng assocated features: the objectve functon and constrants are mplct functons evaluated
More informationThe Codesign Challenge
ECE 4530 Codesgn Challenge Fall 2007 Hardware/Software Codesgn The Codesgn Challenge Objectves In the codesgn challenge, your task s to accelerate a gven software reference mplementaton as fast as possble.
More informationScan Conversion & Shading
Scan Converson & Shadng Thomas Funkhouser Prnceton Unversty C0S 426, Fall 1999 3D Renderng Ppelne (for drect llumnaton) 3D Prmtves 3D Modelng Coordnates Modelng Transformaton 3D World Coordnates Lghtng
More informationWhat are the camera parameters? Where are the light sources? What is the mapping from radiance to pixel color? Want to solve for 3D geometry
Today: Calbraton What are the camera parameters? Where are the lght sources? What s the mappng from radance to pel color? Why Calbrate? Want to solve for D geometry Alternatve approach Solve for D shape
More informationConditional Speculative Decimal Addition*
Condtonal Speculatve Decmal Addton Alvaro Vazquez and Elsardo Antelo Dep. of Electronc and Computer Engneerng Unv. of Santago de Compostela, Span Ths work was supported n part by Xunta de Galca under grant
More informationi v v 6 i 2 i 3 v + (1) (2) (3) (4) (5) Substituting (4) and (5) into (3) (6) = 2 (7) (5) and (6) (8) (4) and (6) ˆ
5V 6 v 6 î v v Ω î Ω v v 8Ω V î v 5 6Ω 5 Mesh : 6ˆ ˆ = Mesh : ˆ 8ˆ = Mesh : ˆ ˆ ˆ 8 0 = 5 Solvng ˆ ˆ ˆ from () = Solvng ˆ ˆ ˆ from () = 7 7 Substtutng () and (5) nto () (5) and (6) 9 ˆ = A 8 ˆ = A 0 ()
More informationModule Management Tool in Software Development Organizations
Journal of Computer Scence (5): 8-, 7 ISSN 59-66 7 Scence Publcatons Management Tool n Software Development Organzatons Ahmad A. Al-Rababah and Mohammad A. Al-Rababah Faculty of IT, Al-Ahlyyah Amman Unversty,
More informationScan Conversion & Shading
1 3D Renderng Ppelne (for drect llumnaton) 2 Scan Converson & Shadng Adam Fnkelsten Prnceton Unversty C0S 426, Fall 2001 3DPrmtves 3D Modelng Coordnates Modelng Transformaton 3D World Coordnates Lghtng
More informationComputer models of motion: Iterative calculations
Computer models o moton: Iteratve calculatons OBJECTIVES In ths actvty you wll learn how to: Create 3D box objects Update the poston o an object teratvely (repeatedly) to anmate ts moton Update the momentum
More informationEdge Detection in Noisy Images Using the Support Vector Machines
Edge Detecton n Nosy Images Usng the Support Vector Machnes Hlaro Gómez-Moreno, Saturnno Maldonado-Bascón, Francsco López-Ferreras Sgnal Theory and Communcatons Department. Unversty of Alcalá Crta. Madrd-Barcelona
More informationA Binarization Algorithm specialized on Document Images and Photos
A Bnarzaton Algorthm specalzed on Document mages and Photos Ergna Kavalleratou Dept. of nformaton and Communcaton Systems Engneerng Unversty of the Aegean kavalleratou@aegean.gr Abstract n ths paper, a
More informationProgramming in Fortran 90 : 2017/2018
Programmng n Fortran 90 : 2017/2018 Programmng n Fortran 90 : 2017/2018 Exercse 1 : Evaluaton of functon dependng on nput Wrte a program who evaluate the functon f (x,y) for any two user specfed values
More informationIn the planar case, one possibility to create a high quality. curve that interpolates a given set of points is to use a clothoid spline,
Dscrete Farng of Curves and Surfaces Based on Lnear Curvature Dstrbuton R. Schneder and L. Kobbelt Abstract. In the planar case, one possblty to create a hgh qualty curve that nterpolates a gven set of
More informationConstruction of second-order accurate monotone and stable residual distribution schemes for steady problems
Journal of Computatonal Physcs 195 (2004) 474 507 www.elsever.com/locate/cp Constructon of second-order accurate monotone and stable resdual dstrbuton schemes for steady problems Rem Abgrall *,1, Mohamed
More informationTopology Design using LS-TaSC Version 2 and LS-DYNA
Topology Desgn usng LS-TaSC Verson 2 and LS-DYNA Wllem Roux Lvermore Software Technology Corporaton, Lvermore, CA, USA Abstract Ths paper gves an overvew of LS-TaSC verson 2, a topology optmzaton tool
More information5 The Primal-Dual Method
5 The Prmal-Dual Method Orgnally desgned as a method for solvng lnear programs, where t reduces weghted optmzaton problems to smpler combnatoral ones, the prmal-dual method (PDM) has receved much attenton
More informationKinematics of pantograph masts
Abstract Spacecraft Mechansms Group, ISRO Satellte Centre, Arport Road, Bangalore 560 07, Emal:bpn@sac.ernet.n Flght Dynamcs Dvson, ISRO Satellte Centre, Arport Road, Bangalore 560 07 Emal:pandyan@sac.ernet.n
More informationAssignment # 2. Farrukh Jabeen Algorithms 510 Assignment #2 Due Date: June 15, 2009.
Farrukh Jabeen Algorthms 51 Assgnment #2 Due Date: June 15, 29. Assgnment # 2 Chapter 3 Dscrete Fourer Transforms Implement the FFT for the DFT. Descrbed n sectons 3.1 and 3.2. Delverables: 1. Concse descrpton
More informationTN348: Openlab Module - Colocalization
TN348: Openlab Module - Colocalzaton Topc The Colocalzaton module provdes the faclty to vsualze and quantfy colocalzaton between pars of mages. The Colocalzaton wndow contans a prevew of the two mages
More informationParallelism for Nested Loops with Non-uniform and Flow Dependences
Parallelsm for Nested Loops wth Non-unform and Flow Dependences Sam-Jn Jeong Dept. of Informaton & Communcaton Engneerng, Cheonan Unversty, 5, Anseo-dong, Cheonan, Chungnam, 330-80, Korea. seong@cheonan.ac.kr
More informationProblem Set 3 Solutions
Introducton to Algorthms October 4, 2002 Massachusetts Insttute of Technology 6046J/18410J Professors Erk Demane and Shaf Goldwasser Handout 14 Problem Set 3 Solutons (Exercses were not to be turned n,
More informationJournal of Computational Physics
Journal of Computatonal Physcs 230 (2011) 2540 2561 Contents lsts avalable at ScenceDrect Journal of Computatonal Physcs journal homepage: www.elsever.com/locate/jcp A grd based partcle method for solvng
More informationBLaC-Wavelets: A Multiresolution Analysis With Non-Nested Spaces. Georges-Pierre Bonneau Stefanie Hahmann Gregory M. Nielson z
BLaC-Wavelets: A Multresoluton Analyss Wth Non-Nested Spaces Georges-Perre Bonneau Stefane Hahmann Gregory M. Nelson z CNRS - Laboratore LMC Grenoble, France Arzona State Unversty Tempe, USA ABSTRACT In
More informationGeometric Error Estimation
Geometrc Error Estmaton Houman Borouchak Project-team GAMMA 3 UTT Troyes, France Emal: houman.borouchak@utt.fr Patrck Laug Project-team GAMMA 3 INRIA Pars - Rocquencourt, France Emal: patrck.laug@nra.fr
More informationOptimization Methods: Integer Programming Integer Linear Programming 1. Module 7 Lecture Notes 1. Integer Linear Programming
Optzaton Methods: Integer Prograng Integer Lnear Prograng Module Lecture Notes Integer Lnear Prograng Introducton In all the prevous lectures n lnear prograng dscussed so far, the desgn varables consdered
More informationSubspace clustering. Clustering. Fundamental to all clustering techniques is the choice of distance measure between data points;
Subspace clusterng Clusterng Fundamental to all clusterng technques s the choce of dstance measure between data ponts; D q ( ) ( ) 2 x x = x x, j k = 1 k jk Squared Eucldean dstance Assumpton: All features
More informationHarvard University CS 101 Fall 2005, Shimon Schocken. Assembler. Elements of Computing Systems 1 Assembler (Ch. 6)
Harvard Unversty CS 101 Fall 2005, Shmon Schocken Assembler Elements of Computng Systems 1 Assembler (Ch. 6) Why care about assemblers? Because Assemblers employ some nfty trcks Assemblers are the frst
More informationCE 221 Data Structures and Algorithms
CE 1 ata Structures and Algorthms Chapter 4: Trees BST Text: Read Wess, 4.3 Izmr Unversty of Economcs 1 The Search Tree AT Bnary Search Trees An mportant applcaton of bnary trees s n searchng. Let us assume
More informationDifferential formulation of discontinuous Galerkin and related methods for compressible Euler and Navier-Stokes equations
Graduate Theses and Dssertatons Graduate College 2011 Dfferental formulaton of dscontnuous Galerkn and related methods for compressble Euler and Naver-Stokes equatons Hayang Gao Iowa State Unversty Follow
More informationEcient Computation of the Most Probable Motion from Fuzzy. Moshe Ben-Ezra Shmuel Peleg Michael Werman. The Hebrew University of Jerusalem
Ecent Computaton of the Most Probable Moton from Fuzzy Correspondences Moshe Ben-Ezra Shmuel Peleg Mchael Werman Insttute of Computer Scence The Hebrew Unversty of Jerusalem 91904 Jerusalem, Israel Emal:
More informationLecture #15 Lecture Notes
Lecture #15 Lecture Notes The ocean water column s very much a 3-D spatal entt and we need to represent that structure n an economcal way to deal wth t n calculatons. We wll dscuss one way to do so, emprcal
More informationAVO Modeling of Monochromatic Spherical Waves: Comparison to Band-Limited Waves
AVO Modelng of Monochromatc Sphercal Waves: Comparson to Band-Lmted Waves Charles Ursenbach* Unversty of Calgary, Calgary, AB, Canada ursenbach@crewes.org and Arnm Haase Unversty of Calgary, Calgary, AB,
More informationHigh resolution 3D Tau-p transform by matching pursuit Weiping Cao* and Warren S. Ross, Shearwater GeoServices
Hgh resoluton 3D Tau-p transform by matchng pursut Wepng Cao* and Warren S. Ross, Shearwater GeoServces Summary The 3D Tau-p transform s of vtal sgnfcance for processng sesmc data acqured wth modern wde
More informationF Geometric Mean Graphs
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 2 (December 2015), pp. 937-952 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) F Geometrc Mean Graphs A.
More informationCHAPTER 2 DECOMPOSITION OF GRAPHS
CHAPTER DECOMPOSITION OF GRAPHS. INTRODUCTION A graph H s called a Supersubdvson of a graph G f H s obtaned from G by replacng every edge uv of G by a bpartte graph,m (m may vary for each edge by dentfyng
More informationTHE PULL-PUSH ALGORITHM REVISITED
THE PULL-PUSH ALGORITHM REVISITED Improvements, Computaton of Pont Denstes, and GPU Implementaton Martn Kraus Computer Graphcs & Vsualzaton Group, Technsche Unverstät München, Boltzmannstraße 3, 85748
More informationCircuit Analysis I (ENGR 2405) Chapter 3 Method of Analysis Nodal(KCL) and Mesh(KVL)
Crcut Analyss I (ENG 405) Chapter Method of Analyss Nodal(KCL) and Mesh(KVL) Nodal Analyss If nstead of focusng on the oltages of the crcut elements, one looks at the oltages at the nodes of the crcut,
More informationAmnon Shashua Shai Avidan Michael Werman. The Hebrew University, objects.
Trajectory Trangulaton over Conc Sectons Amnon Shashua Sha Avdan Mchael Werman Insttute of Computer Scence, The Hebrew Unversty, Jerusalem 91904, Israel e-mal: fshashua,avdan,wermang@cs.huj.ac.l Abstract
More informationNUMERICAL SOLVING OPTIMAL CONTROL PROBLEMS BY THE METHOD OF VARIATIONS
ARPN Journal of Engneerng and Appled Scences 006-017 Asan Research Publshng Network (ARPN). All rghts reserved. NUMERICAL SOLVING OPTIMAL CONTROL PROBLEMS BY THE METHOD OF VARIATIONS Igor Grgoryev, Svetlana
More informationCCT. Abdul-Jabbar: Design And Realization Of Circular Contourlet Transform ..VLSI. z 2. z 1.
Abdul-Jabbar: Desgn And Realzaton Of Crcular Contourlet Transform DESIGN AND REALIZATION OF CIRCULAR CONTOURLET TRANSFORM Dr Jassm M Abdul-Jabbar * Hala N Fathee ** * Ph D n Elect Eng, Dept of Computer
More informationImage Representation & Visualization Basic Imaging Algorithms Shape Representation and Analysis. outline
mage Vsualzaton mage Vsualzaton mage Representaton & Vsualzaton Basc magng Algorthms Shape Representaton and Analyss outlne mage Representaton & Vsualzaton Basc magng Algorthms Shape Representaton and
More informationan assocated logc allows the proof of safety and lveness propertes. The Unty model nvolves on the one hand a programmng language and, on the other han
UNITY as a Tool for Desgn and Valdaton of a Data Replcaton System Phlppe Quennec Gerard Padou CENA IRIT-ENSEEIHT y Nnth Internatonal Conference on Systems Engneerng Unversty of Nevada, Las Vegas { 14-16
More informationKiran Joy, International Journal of Advanced Engineering Technology E-ISSN
Kran oy, nternatonal ournal of Advanced Engneerng Technology E-SS 0976-3945 nt Adv Engg Tech/Vol. V/ssue /Aprl-une,04/9-95 Research Paper DETERMATO O RADATVE VEW ACTOR WTOUT COSDERG TE SADOWG EECT Kran
More informationTsinghua University at TAC 2009: Summarizing Multi-documents by Information Distance
Tsnghua Unversty at TAC 2009: Summarzng Mult-documents by Informaton Dstance Chong Long, Mnle Huang, Xaoyan Zhu State Key Laboratory of Intellgent Technology and Systems, Tsnghua Natonal Laboratory for
More informationAssembler. Shimon Schocken. Spring Elements of Computing Systems 1 Assembler (Ch. 6) Compiler. abstract interface.
IDC Herzlya Shmon Schocken Assembler Shmon Schocken Sprng 2005 Elements of Computng Systems 1 Assembler (Ch. 6) Where we are at: Human Thought Abstract desgn Chapters 9, 12 abstract nterface H.L. Language
More informationElectrical analysis of light-weight, triangular weave reflector antennas
Electrcal analyss of lght-weght, trangular weave reflector antennas Knud Pontoppdan TICRA Laederstraede 34 DK-121 Copenhagen K Denmark Emal: kp@tcra.com INTRODUCTION The new lght-weght reflector antenna
More informationAccounting for the Use of Different Length Scale Factors in x, y and z Directions
1 Accountng for the Use of Dfferent Length Scale Factors n x, y and z Drectons Taha Soch (taha.soch@kcl.ac.uk) Imagng Scences & Bomedcal Engneerng, Kng s College London, The Rayne Insttute, St Thomas Hosptal,
More informationPHANTOM VORTICITY IN EULER SOLUTIONS ON HIGHLY STRETCHED GRIDS
ICAS 000 CONGRESS PHANTOM VORTICITY IN EULER SOLUTIONS ON HIGHLY STRETCHED GRIDS S. A. Prnce, D. K. Ludlow, N. Qn Cranfeld College of Aeronautcs, Bedfordshre, UK Currently DERA Bedford, UK Keywords: Phantom
More informationPositivity-Preserving Well-Balanced Discontinuous Galerkin Methods for the Shallow Water Equations on Unstructured Triangular Meshes
Postvty-Preservng Well-Balanced Dscontnuous Galerkn Methods for the Shallow Water Equatons on Unstructured Trangular Meshes The MIT Faculty has made ths artcle openly avalable. Please share how ths access
More informationU.C. Berkeley CS294: Beyond Worst-Case Analysis Handout 5 Luca Trevisan September 7, 2017
U.C. Bereley CS294: Beyond Worst-Case Analyss Handout 5 Luca Trevsan September 7, 207 Scrbed by Haars Khan Last modfed 0/3/207 Lecture 5 In whch we study the SDP relaxaton of Max Cut n random graphs. Quc
More informationKFUPM. SE301: Numerical Methods Topic 8 Ordinary Differential Equations (ODEs) Lecture (Term 101) Section 04. Read
SE3: Numercal Metods Topc 8 Ordnar Dfferental Equatons ODEs Lecture 8-36 KFUPM Term Secton 4 Read 5.-5.4 6-7- C ISE3_Topc8L Outlne of Topc 8 Lesson : Introducton to ODEs Lesson : Talor seres metods Lesson
More informationsuch that is accepted of states in , where Finite Automata Lecture 2-1: Regular Languages be an FA. A string is the transition function,
* Lecture - Regular Languages S Lecture - Fnte Automata where A fnte automaton s a -tuple s a fnte set called the states s a fnte set called the alphabet s the transton functon s the ntal state s the set
More informationMath Homotopy Theory Additional notes
Math 527 - Homotopy Theory Addtonal notes Martn Frankland February 4, 2013 The category Top s not Cartesan closed. problem. In these notes, we explan how to remedy that 1 Compactly generated spaces Ths
More informationProper Choice of Data Used for the Estimation of Datum Transformation Parameters
Proper Choce of Data Used for the Estmaton of Datum Transformaton Parameters Hakan S. KUTOGLU, Turkey Key words: Coordnate systems; transformaton; estmaton, relablty. SUMMARY Advances n technologes and
More informationFair Triangle Mesh Generation with Discrete Elastica
Far Trangle Mesh Generaton wth Dscrete Elastca Shn Yoshzawa, and Alexander G. Belyaev, Computer Graphcs Group, Max-Planck-Insttut für Informatk, 66123 Saarbrücken, Germany Phone: [+49](681)9325-414 Fax:
More informationRepeater Insertion for Two-Terminal Nets in Three-Dimensional Integrated Circuits
Repeater Inserton for Two-Termnal Nets n Three-Dmensonal Integrated Crcuts Hu Xu, Vasls F. Pavlds, and Govann De Mchel LSI - EPFL, CH-5, Swtzerland, {hu.xu,vasleos.pavlds,govann.demchel}@epfl.ch Abstract.
More informationAn accelerated value/policy iteration scheme for the solution of DP equations
An accelerated value/polcy teraton scheme for the soluton of DP equatons Alessandro Alla 1, Maurzo Falcone 2, and Dante Kalse 3 1 SAPIENZA - Unversty of Rome, Ple. Aldo Moro 2, Rome, Italy alla@mat.unroma1.t
More informationJournal of Hydrology
Journal of Hydrology 59 (24) 23 32 Contents lsts avalable at ScenceDrect Journal of Hydrology journal homepage: www.elsever.com/locate/jhydrol Hgh resoluton numercal schemes for solvng knematc wave equaton
More informationPolyhedral Compilation Foundations
Polyhedral Complaton Foundatons Lous-Noël Pouchet pouchet@cse.oho-state.edu Dept. of Computer Scence and Engneerng, the Oho State Unversty Feb 8, 200 888., Class # Introducton: Polyhedral Complaton Foundatons
More information