A Cartesian grid finite-difference method for 2D incompressible viscous flows in irregular geometries
|
|
- Lester Waters
- 6 years ago
- Views:
Transcription
1 Journal of Computatonal Physcs 204 (2005) A Cartesan grd fnte-dfference method for 2D ncompressble vscous flows n rregular geometres E. Sanmguel-Rojas a, J. Ortega-Casanova b, C. del Pno b, R. Fernandez-Fera b, * a Unversdad Poltécnca de Cartagena, E.T.S. Ingeneros Industrales, Cartagena, Murca, Span b Unversdad de Málaga, E.T.S. Ingeneros Industrales, Málaga, Span Receved 23 Aprl 2004 receved n revsed form 11 October 2004 accepted 12 October 2004 Avalable onlne 11 November 2004 Abstract A method for generatng a non-unform Cartesan grd for rregular two-dmensonal (2D) geometres such that all the boundary ponts are regular mesh ponts s gven. The resultng non-unform grd s used to dscretze the Naver Stokes equatons for 2D ncompressble vscous flows usng fnte-dfference approxmatons. To that end, fntedfference approxmatons of the dervatves on a non-unform mesh are gven. We test the method wth two dfferent examples: the shallow water flow on a lake wth rregular contour and the pressure drven flow through an rregular array of crcular cylnders. Ó 2004 Elsever Inc. All rghts reserved. 1. Introducton The numercal smulaton of flows wth rregular geometres s a problem of ncreasng nterest. In partcular, much effort have been dedcated n recent years to the use of Cartesan grds whch does not conform to the rregular boundares [1 6]. In relaton to the conventonal structured-grd approach wth curvlnear grds that conforms to the boundares, ths approach has the man advantage of ts smplcty, both n the grd generaton and n the governng equatons. In addton, the transformaton of the governng equatons to a curvlnear coordnate system that conforms to very complcated boundares s not an easy task, and usually affects to the stablty, convergence, and accuracy of the numercal solver. In ths paper we present a technque for Cartesan grd generaton that conforms to rregular twodmensonal (2D) boundares. It has the advantage of workng wth a Cartesan mesh n whch all the * Correspondng author. E-mal address: ramon.fernandez@uma.es. (R. Fernandez-Fera) /$ - see front matter Ó 2004 Elsever Inc. All rghts reserved. do: /j.jcp
2 E. Sanmguel-Rojas et al. / Journal of Computatonal Physcs 204 (2005) boundares nodes are regular nodes of the grd, thus avodng the usual complcated nterpolatons needed for the Cartesan cells cut by mmersed boundares, or the use of artfcal body forces, or other artfces such as a dstrbuton of vortcty sources, to mpose the boundary condtons (see, for example [1 6]). The prce one has to pay s that the Cartesan grd s non-unform. For ths reason we develop second-order accurate fnte-dfference approxmatons of the dervatves for non-equdstant grd ponts that substtutes the usual fnte-dfference approxmatons for unform grds (some of these expressons are already reported n the lterature see, e.g. [7]). Thus, we are led to dscrete equatons wth the same level of complexty than the Cartesan equatons dscretzed on a unform grd, but conformng to an rregular geometry. One dsadvantage n relaton to the nterpolaton technques gven n some of the references cted above s that the present method s not suted for movng nterfaces. However, second-order nterpolaton technques [4] are not easly appled to Neumann boundary condtons. The structure of the paper s the followng. Secton 2 ntroduces the grd generaton technque on a generc, rregular 2D doman. The expresson for the fnte-dfference approxmatons of the Cartesan dervatves on non-unform grds are gven n Secton 3. Sectons 4 valdates the method wth two examples qute dfferent to each other: the 2D shallow water, wnd-drven flow on a lake wth rregular contour, and the 2D ncompressble, pressure-drven flow around an rregular array of crcular cylnders. Some conclusons are drawn n the last secton. 2. Cartesan grd generaton that conforms to an rregular 2D doman Consder the 2D rregular doman of Fg. 1. Our objectve s to generate a Cartesan grd where all the boundary ponts are regular mesh ponts. Ths means that all the nteror ponts have to be collocated n relaton to a set of selected boundary ponts, and that the resultng Cartesan grd wll not be unform. In order to smplfy the storage of the grd ponts locaton n matrx form, what we propose here s a ray tracng technque. One starts at a gven boundary pont (marked wth a crcle n Fg. 1), and generates a set of boundary ponts by ÔCartesan reflectonsõ of the ray (squares n Fg. 1). In order to avod an nfnte regress, the process ends when the ray reaches a boundary perpendcular to t (.e., when t reaches a secton of the boundary parallel to one of the Cartesan axs), or when a boundary node s generated very close to a prevous one (ther separaton s less than a gven tolerance). It s mportant to detect frst the man boundary ponts or ponts of ntersecton between the several sectons of the boundary (crcles n Fg. 2). The frst rays wll start from these man ponts, dvdng the doman, and the boundary, n a number sectons. Each of these sectons s then dvded usng a number of ponts on each boundary secton whch depends on the desred precson. The resultng mesh (see Fg. 2(b)) concentrates the nodes wth the desred precson at the dfferent sectons of the boundary. Programng ths technque s relatvely easy and the storng n matrx form of the resultng grd ponts locatons s also straghtforward. The technque s a lttle more complcated n domans wth concave boundary sectons lke that depcted n Fg. 3. We have traced three rays startng from the three crcled ponts. These rays generate a seres of mesh ponts, both on the boundary and nsde the doman, before stopng at a boundary parallel to the x-axs. The last three nodes on that boundary have no correspondng ones on the lower boundary. Thus, n order to facltate the storage of the nodes n matrx form,.e., n order to have an structured grd, t s convenent to contnue these rays tll the lowest boundary (dashed lnes n Fg. 3(a)). Though we store all the nodes, we can dsregard the three trangled nodes n Fg. 3(b), beng the effectve node to the rght of node j that labelled wth j + 4. Once the grd s generated, one may create an ndrect access wth a new ndex, say jj, such that jj + 1 corresponds to j + 4. Thus, the fnal non-structured grd s stored wthn a structured grd form of n m nodes (Fg. 3(c)). The ndrect access through the ndexes (,jj) tell us whether a node of the structured grd s an actual node of the computatonal non-structured grd.
3 304 E. Sanmguel-Rojas et al. / Journal of Computatonal Physcs 204 (2005) (a) (b) Fg. 1. Generc 2D doman (a) and llustraton of grd generaton by ray tracng (b). In very rregular domans t may happen that ths ray tracng method may generate very small cells near some boundares, or even nsde the doman. To avod ths we use a lower lmt for the cell sze, n such a way that nodes that generate cells wth sze less than ths lmt are dscarded, n a smlar way to what s done n concave boundares, as descrbed above. 3. Fnte-dfference approxmaton on non-unform meshes In order to dscretze the flow equatons n the non-unform grd developed n the above secton, one has to use fnte-dfference approxmatons on a non-unform mesh. In ths secton we develop these fntedfference expressons for all the spatal dervatves appearng n the Naver Stokes equatons. Some of them have been prevously reported by Turkel [7]. In partcular, Turkel provdes the frst dervatve, and the centered form of the second dervatve wth frst-order truncaton error. All the expressons we gve below (some of them are gven n Appendx A) are second-order accurate, and we nclude forward and backwards expressons for the second-order dervatves, whch are needed at the boundares. It has been shown [,9] that second-order accuracy can be obtaned (on non-unform grds) even though local truncaton
4 E. Sanmguel-Rojas et al. / Journal of Computatonal Physcs 204 (2005) Fg. 2. Man ponts (crcles), and resultng non-unform grd (b). errors are of lower order. However, ths s vald only for non-unform grds n whch the varaton of the mesh sze s very small and for lnear equatons. Snce we want to apply the fnte-dfferences method to arbtrary non-unform meshes, and to the non-lnear Naver Stokes equatons, we need second-order truncaton errors to reach second-order accuracy. Ths s mportant n problems wth a long tme evoluton, such as the examples gven below, where second-order accuracy s needed at both nner and boundary nodes. Consder a 1D non-unform grd wth nx + 1 dscrete ponts (0 6 6 nx) located arbtrarly on the unt length (Fg. 4). If the value of a generc functon f(x) and ts dervatves are known at the pont th, x = Dx, Dx =1/nx, the values of f at the ponts ± 1 and ± 2 can be approxmated usng Taylor expansons. Indcatng wth a subscrpt the grd pont, and wth prmes the dervatves wth respect to x, the unknown values f ±1 and f ±2 can be wrtten as f þ1 ¼ f þ h þ1 f 0 þ h2 þ1 2 f 00 þ h3 þ1 6 f 000 þ Oðh 4 þ1þ ð1þ f 1 ¼ f þ h 1 f 0 þ h2 1 2 f 00 þ h3 1 6 f 000 þ Oðh 4 1Þ ð2þ
5 306 E. Sanmguel-Rojas et al. / Journal of Computatonal Physcs 204 (2005) (a) (b) (c) Fg. 3. (a) and (b) Fcttous grd ponts n a doman wth concave boundary sectons. (c) Fnal non-structured grd (squares). Fg. 4. Nonunform grd wth nx + 1 grd ponts dstrbuted arbtrarly.
6 E. Sanmguel-Rojas et al. / Journal of Computatonal Physcs 204 (2005) f þ2 ¼ f þ h þ2 f 0 þ h2 þ2 2 f 00 þ h3 þ2 6 f 000 þ Oðh 4 þ2þ ð3þ f 2 ¼ f þ h 2 f 0 þ h2 2 2 f 00 þ h3 2 6 f 000 þ Oðh 4 2Þ ð4þ where the last terms n these expressons ndcate the order of the truncaton error of the approxmaton, and h þ1 ¼ x þ1 x h 1 ¼ x 1 x h þ2 ¼ x þ2 x h 2 ¼ x 2 x : ð5þ These expansons can be used to approxmate the nth dervatve of f at the pont up to any order of the truncaton error, provded that they are convenently combned. Wth second-order accuracy, the fntedfference approxmaton n the centered form for the frst and second dervatves of f at the pont are (for hgher dervatves, or hgher order for the desred truncaton error, more ponts than those consdered n (1) (4) are needed): h m ¼ þ1 h þ1 h 1 þh 2 1 h f 0 n ¼ 1 ¼ m f 1 þ s f þ n f þ1 þ E h þ1 h 1 þh 2 þ1 ð6þ s ¼ ðm þ n Þ >: E ¼ 1 h 6 þ1h 1 f 000 f 00 ¼ m f 1 þ s f þ n f þ1 þ p f þ2 þ E f 00 ¼ p f 2 þ m f 1 þ s f þ n f þ1 þ E m ¼ n ¼ p ¼ 2ðh þ1 þh þ2 Þ ðh þ1 h þ2 h 1 h þ1 þh 2 1 h 1h þ2 Þh 1 2ðh 1 þh þ2 Þ ð h þ1 h 1 þh 1 h þ2 þh 2 þ1 h þ1h þ2 Þh þ1 2ðh þ1 þh 1 Þ ðh þ1 h þ2 h 1 h þ1 h 2 þ2 þh 1h þ2 Þh þ2 s ¼ ðm þ n þ p Þ >: E ¼ 1 m 12 h 4 1 þ n h 4 þ1 þ p h 4 þ2 f v m ¼ n ¼ p ¼ 2ðh þ1 þh 2 Þ ðh þ1 h 2 h 1 h þ1 þh 2 1 h 1h 2 Þh 1 2ðh 1 þh 2 Þ ð h þ1 h 1 þh 1 h 2 þh 2 þ1 h þ1h 2 Þh þ1 2ðh þ1 þh 1 Þ ðh þ1 h 2 h 1 h þ1 h 2 2 þh 1h 2 Þh 2 s ¼ ðm þ n þ p Þ >: E ¼ 1 p 12 h 4 2 þ m h 4 1 þ n h 4 þ1 f v where E s the truncaton error of each approxmaton. Ths truncaton error s a functon of the separaton between the local ponts around. Thus, n a non-unform grd, expressons (6) () wll be more accurate as the grd ponts becomes closer. Note that two dfferent expressons for the second dervatve are gven, (7) and (), nether of them fully centered on the pont. Ths s because one needs four grd ponts to have a second-order truncaton error n the second dervatve, so that two dfferent expressons result dependng on whether we use the pont + 2, or the pont 2. ð7þ ðþ
7 30 E. Sanmguel-Rojas et al. / Journal of Computatonal Physcs 204 (2005) If the grd were unform, h h +1, h 1 = h, h +2 =2h, h 2 = 2h, expressons (6) () are, obvously, the standard centered second-order fnte-dfference approxmaton for the frst and second dervatves: f 0 ¼ m f 1 þ s f þ n f þ1 þ E m ¼ 1 2h n ¼ 1 2h s ¼ 0 >: E ¼ 1 6 h2 f 000 ð9þ m ¼ 1 h 2 n ¼ 1 h 2 f 00 ¼ m f 1 þ s f þ n f þ1 þ p f 2 þ E p ¼ 0 ð10þ s ¼ 2 : h >: 2 E ¼ 1 12 h2 f v Now the grd ponts ± 2 do not appear n the approxmaton to the second dervatve because p =0. From (1) (4) one can obtan not only centered approxmatons, but also forward or backward ones. These expressons, whch are needed at the boundares, are gven n Appendx A wth second-order accuracy. We also gve there fnte-dfference expressons needed to apply Neumann boundary condtons wth second-order accuracy on a non-unform grd. 4. Results 4.1. Wnd-drven flow n a lake One of the man ntended applcatons of the present technque s the smulatons of 2D envronmental flows, such as the flows n shallow lakes and reservors, whch usually have very rregular geometres. For ths reason, as a frst example of a 2D flow n a complex doman we consder the wnddrven flow n Lake Belau (Northern Germany), for whch Podsetchne and Schernewsk [10] reported numercal and expermental results whch can be used to compare wth. The bathymetry of the lake s gven n Fg. 5. Snce the lake s not very deep, one may use the vertcally ntegrated equatons of contnuty and momentum on the horzontal plane (x, y), or shallow water approxmaton (see, for example, [11]): of þrv ¼ 0 ot ð11þ ov ot þr vv h þ ghrf m r 2 v 2rf r v v r 2 f f ^ v kw j W jþ gv j v j ¼ 0: ð12þ H H H C 2 2 H In these equatons, v =[u(x, y, t), v(x, y, t)] s the depth-averaged horzontal velocty, v Z f h v h dz wth v h (x,y,z,t) the local horzontal velocty, H(x,y,t) h(x,y) +f(x,y,t) s the total water depth, wth h the depth below the horzontal reference plane (z = 0) and f the water surface elevaton above z = 0, $ = o/ox + o/oy, g. 9.1 m/s 2 s the acceleraton due to gravty, f = fe z, wth f s 1 s the Corols parameter, m s the averaged, horzontal eddy vscosty, and C s the Chezy coeffcent. ð13þ
8 E. Sanmguel-Rojas et al. / Journal of Computatonal Physcs 204 (2005) Fg. 5. Dgtal contour and bathymetry (n meters) of Lake Belau such as they are used n the computatons (taken from Fg. 2 n [10]). The flow wll be drven by a wnd of velocty W =(W x,w y ), through the term kwjwj, where k = q a C W /q, wth q and q a the denstes of water and ar, respectvely, and C W the wnd drag coeffcent. The numercal values of the parameters that wll be used to solve these equatons are the followng [10]: m = 0.01 m 2 / s, C =40 m 1/2 /s, C W = 0.002, q =10 3 kg/m 3, q a = kg/m 3, and a spatally unform south-westerly wnd (a headng of 220 ) wth a speed jwj of 6 m/s. The boundary condtons at the contour of the lake are u = v =0. A mesh of 7517 grd ponts has been generated usng the technque of Secton 2 (see Fg. 6). The equatons have been dscretzed n ths non-unform grd usng the fnte-dfference approxmatons gven n Secton 3. In partcular, we have used ArakawaÕs grd of the type C [12], where the water elevaton f s evaluated at the grd ponts, whle the averaged velocty components are evaluated at the md ponts of ther respectve cell sdes (see Fg. 7). An explct, two-step, second-order accurate, predctor corrector scheme has been used to advance n tme. If one wrtes Eqs. (11) and (12) schematcally as of/ot = A(v), ov/ot = B(f, v), these two steps are gven by
9 310 E. Sanmguel-Rojas et al. / Journal of Computatonal Physcs 204 (2005) y (m) 250 y (m) x (m) x (m) Fg. 6. Two detals of the dscretzed geometry of the lake showng some of the grd ponts. Fg. 7. ArakawaÕs scheme used n the computatons, where f s evaluated at the grd ponts (trangles), u s evaluated at the squares, and v at the crcles.
10 predctor step: v ¼ v n þ Dt 2 Bðfn v n Þ f ¼ f n þ Dt 2 Aðv Þ E. Sanmguel-Rojas et al. / Journal of Computatonal Physcs 204 (2005) corrector step: v nþ1 ¼ v n þ DtBðf v Þ ð15þ f nþ1 ¼ f n þ DtAðv nþ1 Þ where the superscrpts denote the nstant of tme, and Dt s the tme step. The numercal computatons are started at t = 0 wth the flud at rest and f = 0. We use Dt = 1 s. The results for the velocty feld at t =3h are plotted n Fg.. These results compares very well wth those gven n Fg. 4(a) of Podsetchne and Schernewsk [10], who used a fnte-element method on a trangular mesh created wth a commercal 1000 ð14þ N W S Wnd E 600 y (m) Velocty scale 10 cm/s x (m) Fg.. Averaged velocty feld at t =3h.
11 312 E. Sanmguel-Rojas et al. / Journal of Computatonal Physcs 204 (2005) software package to solve the shallow water equatons. These authors also checked ther numercal results wth expermental measurements Pressure-drven flow through an array of crcular cylnders As a second example we have selected the 2D ncompressble flow around an rregular array of three crcular cylnders nsde a channel (see Fg. 9). In partcular we have consdered the pressure drven flow [13] orgnated by a gven pressure dfference set between the nlet and the outlet. The dmensonless equatons are rv ¼ 0 ð16þ ov ot þ v rv ¼ rpþ 1 Re r2 v ð17þ where v =(u,v) and p are the dmensonless velocty and pressure, respectvely. To non-dmensonalze these equatons we have used the dameter D of the cylnders as the length scale, and p a characterstc velocty based on the the pressure dfference Dp c between the nlet and the outlet, V c ¼ ffffffffffffffffffffffff Dp c =q, where q s the flud densty the Reynolds number s based on ths velocty [13] Re ¼ V cd m ¼ sffffffffffffff Dp c q D m where m s the knematc vscosty of the flud. The moton of the flud s set by the boundary condtons pðx ¼ 0 y ¼ 10 tþ ¼1 pðx ¼ 20 y ¼ 10 tþ ¼0: ð19þ The remanng boundary condtons are v = 0 on the cylnders, and at the channel walls, y = 0 and y = 10. The equatons are solved numercally wth the second-order (both n space and tme) projecton method descrbed n [13], usng a fnte-dfference scheme on a non-unform grd generated wth the method ð1þ y x Fg. 9. Geometry of the channel flow through three crcular cylnders.
12 E. Sanmguel-Rojas et al. / Journal of Computatonal Physcs 204 (2005) y Fg. 10. Mesh ponts concentrated around the cylnders. x Re q t Fg. 11. Reynolds number based on the ext flow rate as a functon of tme.
13 314 E. Sanmguel-Rojas et al. / Journal of Computatonal Physcs 204 (2005) descrbed above. In partcular we have used a grd of mesh ponts, and Dt = A detal of the grd near two of the cylnders s depcted n Fg. 10. To solve numercally the Posson equaton for the pressure we use an ADI based technque, and standard solvers for band matrces wth LU factorzaton from Blas and Lapack packages. The fact that we have now, n general, a non-structured grd does not affect to the effcency of these Posson solvers because what s suppled to them are the actual computatonal nodes and ther correspondng dscretzed equatons through the ndrect access mentoned n Secton 2. Results for Re = 50 are plotted n Fgs. 11 and 12. The flow evolves n tme from rest untl t reaches a quas-perodc state. Ths s shown n Fg. 11, where we plot as a functon of tme the Reynolds number based on the flow rate at the ext of the channel (x = L = 20), Re q q H Re Z H qðtþ ¼ uðx ¼ L y tþdy 0 where H = 10 s the wdth of the channel. Fnally, streamlnes and sobars at t = 106 are plotted n Fg. 12. ð20þ Fg. 12. Streamlnes (a), and sobars (b), at t = 106.
14 E. Sanmguel-Rojas et al. / Journal of Computatonal Physcs 204 (2005) y (a) x y (b) x Fg. 13. Flow nsde a channel wth a sngle crcular cylnder for Re q = 10. (a) Streamlnes (all the lengths are scaled wth the dameter of the cylnder). (b) Detal of the wake behnd the cylnder. The above results show that the current method has no dffculty n resolvng the complex flow pattern and recrculaton regons behnd the cylnders. However, snce no prevous results are reported for ths partcular flow, we consder next the case of a sngle cylnder nsde a channel (see Fg. 13). Ths flow, wth the dmensons gven n Fg. 13(a), has been studed expermentally and numercally [14]. It has been reported that the length l of the axsymmetrc wake behnd the cylnder for moderately small Reynolds number vares lnearly wth Re q. Thus, for nstance, for Re q = 10 (whch n our smulaton corresponds to Re. 17), the length of the wake behnd the cylnder reported s l (tmes de dameter of the cylnder) [14]. Fg. 13(b) shows a detal of the streamlnes behnd the cylnder obtaned numercally by our numercal code, showng that l. 0.21, n close agreement wth the prevous result. 5. Conclusons We have presented here a fnte-dfference method n a non-unform Cartesan grd whch allows us to solve 2D unsteady vscous flows n rregular geometres. It has all the advantages of solvng numercally the flow equatons n a Cartesan mesh (smplcty, accuracy, stablty), and has the mportant pecularty that all the boundares are ftted to the regular ponts of the mesh, so that complcated nterpolatons at the boundares are avoded. We have developed a smple method for generatng such a non-unform
15 316 E. Sanmguel-Rojas et al. / Journal of Computatonal Physcs 204 (2005) Cartesan mesh n complex geometres, and provded second-order fnte-dfference approxmatons n non-unform meshes, whch allow us to solve the flow equatons wth second-order accuracy n all the flow doman, ncludng the boundares. To check the valdty of the method n flows wth both rregular external boundares and complex mmersed boundares, we have selected two typcal examples where we thnk ths method wll be of most nterest: the 2D (shallow water) flow n a lake wth a very complcated external boundary, and a vscous flow wth a complex nternal boundary such as the flow through an rregular array of cylnders between parallel plates. These flows are solved wth second-order tme accuracy usng two dfferent numercal schemes, and the results show that the non-unform Cartesan mesh s able to accurately smulate these flows wth easy. Acknowledgement Ths research has been supported by the COPT of the Junta de Andalucía (Contract No. 07/ ). The numercal smulatons have been made n the computer facltes at the SAIT (U.P. Cartagena) and at the SCAI (U. Málaga). Appendx A. In ths appendx we gve some addtonal expressons for fnte-dfference approxmatons on a non-unform grd that we have used n our computatons. For example, to apply Drchlet boundary condtons we need forward or backward expressons for the frst and second dervatves. Wth second-order accuracy, these one sded approxmatons are: f 0 ¼ s f þ m f 1 þ p f 2 þ E f 0 ¼ s f þ n f þ1 þ p f þ2 þ E for the frst dervatve, and f 00 ¼ s f þ m f 1 þ p f 2 þ q f 3 þ E h m ¼ 2 h 1 ðh 1 h 2 Þ h p ¼ 1 h 2 ðh 1 h 2 Þ s ¼ ðm þ p Þ >: E ¼ 1 h 6 1h 2 f 000 h n ¼ þ2 h þ1 ðh þ1 h þ2 Þ h p ¼ þ1 h þ2 ðh þ1 h þ2 Þ s ¼ ðn þ p Þ >: E ¼ 1 h 6 þ1h þ2 f 000 m ¼ p ¼ q ¼ 2ðh 2 þh 3 Þ ð h 2 h 1 þh 2 h 3 þh 2 1 h 3h 1 Þh 1 2ðh 1 þh 3 Þ ð h 3 h 1 þh 2 h 1 h 2 2 þh 2h 3 Þh 2 2ðh 1 þh 2 Þ ð h 3 h 1 þh 2 h 1 þh 2 3 h 2h 3 Þh 3 s ¼ ðm þ p þ q Þ >: E ¼ 1 m 12 h 4 1 þ p h 4 2 þ q h 4 3 f v ða:1þ ða:2þ ða:3þ
16 E. Sanmguel-Rojas et al. / Journal of Computatonal Physcs 204 (2005) f 00 ¼ s f þ n f þ1 þ p f þ2 þ q f þ3 þ E n ¼ p ¼ q ¼ 2ðh þ2 þh þ3 Þ ð h þ2 h þ1 þh þ2 h þ3 þh 2 þ1 h þ3h þ1 Þh þ1 2ðh þ1 þh þ3 Þ ð h þ3 h þ1 þh þ2 h þ1 h 2 þ2 þh þ2h þ3 Þh þ2 2ðh þ1 þh þ2 Þ ð h þ3 h þ1 þh þ2 h þ1 þh 2 þ3 h þ2h þ3 Þh þ3 s ¼ ðn þ p þ q Þ >: E ¼ 1 n 12 h 4 þ1 þ p h 4 þ2 þ q h 4 þ3 f v ða:4þ for the second dervatve (h ±3 = x ±3 x ). To apply Neumann boundary condtons we need the second-order dervatve n terms of the the frstorder one. Wth second-order accuracy, these expressons (forward and backward) are: f 00 ¼ r f 0 þ s f þ n f þ1 þ p f þ2 þ E f 00 ¼ r f 0 þ s f þ n f 1 þ p f 2 þ E r ¼ 2ðh þ1þh þ2 Þ h þ1 h þ2 s ¼ 2ðh2 þ1 þh þ1h þ2 þh 2 þ2 Þ h 2 þ1 þh2 þ2 n ¼ 2h þ2 p ¼ >: ðh þ1 h þ2 Þh 2 þ1 2h þ1 ðh þ1 h þ2 Þh 2 þ2 E ¼ 1 h 12 þ1h þ2 f v r ¼ 2ðh 1þh 2 Þ h 1 h 2 s ¼ 2ðh2 1 þh 1h 2 þh 2 2 Þ h 2 1 þh2 2 n ¼ 2h 2 p ¼ >: ðh 1 h 2 Þh 2 1 2h 1 ðh 1 h 2 Þh 2 2 E ¼ 1 12 h 1h 2 f v : ða:5þ ða:6þ References [1] T. Ye, R. Mttal, H.S. Udaykumar, W. Shyy, An accurate cartesan grd method for vscous ncompressble flows wth complex mmersed boundares, J. Comput. Phys. 156 (1999) [2] D. Calhoun, R.J. LeVeque, An cartesan grd fnte-volume method for the advecton dffuson equaton n rregular geometres, J. Comput. Phys. 157 (2000) [3] E.A. Fadlun, R. Verzcco, P. Orland, J. Mohd-Yusof, Combned mmersed-boundary fnte-dfference methods for threedmensonal complex flow smulatons, J. Comput. Phys. 161 (2000) [4] F. Gbou, R. Fedkw, L.T. Cheng, M. Kang, A second-order-accurate symmetrc dscretzaton of the Posson equaton on rregular domans, J. Comput. Phys. 176 (2002) [5] D. Calhoun, A cartesan grd method for solvng the two-dmensonal streamfuncton vortcty equatons n rregular regons, J. Comput. Phys. 176 (2002) [6] A. Glmanov, F. Sotropoulos, E. Balaras, A general reconstructon algorthm for smulatng flows wth complex 3D mmersed boundares on cartesan grds, J. Comput. Phys. 191 (2003) [7] E. Turkel, Accuracy of schemes wth nonunform meshes for compressble flud-flows, J. Appl. Numer. Math. 2 (196) [] T.A. Manteuffel, A.B. Whte, The numercal soluton of second-order boundary value problems on nonunform meshes, Math. Comput. 47 (196)
17 31 E. Sanmguel-Rojas et al. / Journal of Computatonal Physcs 204 (2005) [9] H.O. Kress, T.A. Manteuffel, B. Swartz, B. Wendroff, A.B. Whte, Supra-convergent schemes on rregular grds, Math. Comput. 47 (196) [10] V. Podsetchne, G. Schernewsk, The nfluence of spatal wnd nhomogenety on flow patterns n a small lake, Wat. Res. 33 (1999) [11] T. Wejan, Shallow Water Hydrodynamcs, Elsever, Amsterdam, [12] A. Arakawa, V. Lamb, Computatonal desgn of the basc dynamcal processes of the UCLA general crculaton model, Methods Comput. Phys. 17 (1977) [13] R. Fernandez-Fera, E. Sanmguel-Rojas, An explct projecton method for solvng ncompressble flows drven by a pressure dfference, Comput. Fluds 33 (2004) [14] J.H. Chen, W.G. Prtchard, S.J. Tavener, Bfurcaton for flow past a cylnder between parallel planes, J. Flud Mech. 24 (1995)
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 informationModeling of Airfoil Trailing Edge Flap with Immersed Boundary Method
Downloaded from orbt.dtu.dk on: Sep 27, 2018 Modelng of Arfol Tralng Edge Flap wth Immersed Boundary Method Zhu, We Jun; Shen, Wen Zhong; Sørensen, Jens Nørkær Publshed n: ICOWEOE-2011 Publcaton date:
More informationImmersed Boundary Method for the Solution of 2D Inviscid Compressible Flow Using Finite Volume Approach on Moving Cartesian Grid
Journal of Appled Flud Mechancs, Vol. 4, No. 2, Specal Issue, pp. 27-36, 2011. Avalable onlne at www.jafmonlne.net, ISSN 1735-3572, EISSN 1735-3645. Immersed Boundary Method for the Soluton of 2D Invscd
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 informationVISCOELASTIC SIMULATION OF BI-LAYER COEXTRUSION IN A SQUARE DIE: AN ANALYSIS OF VISCOUS ENCAPSULATION
VISCOELASTIC SIMULATION OF BI-LAYER COEXTRUSION IN A SQUARE DIE: AN ANALYSIS OF VISCOUS ENCAPSULATION Mahesh Gupta Mchgan Technologcal Unversty Plastc Flow, LLC Houghton, MI 49931 Hancock, MI 49930 Abstract
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 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 informationStorm Surge and Tsunami Simulator in Oceans and Coastal Areas
Proc. Int. Conf. on Montorng, Predcton and Mtgaton of Water-Related Dsasters, Dsaster Preventon Research Insttute, Kyoto Unv. (2005) Storm Surge and Tsunam Smulator n Oceans and Coastal Areas Taro Kaknuma
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 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 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 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 informationApplications of DEC:
Applcatons of DEC: Flud Mechancs and Meshng Dscrete Dfferental Geometry Overvew Puttng DEC to good use Fluds, fluds, fluds geometrc nterpretaton of classcal models dscrete geometrc nterpretaton new geometry-based
More informationImprovement of Spatial Resolution Using BlockMatching Based Motion Estimation and Frame. Integration
Improvement of Spatal Resoluton Usng BlockMatchng Based Moton Estmaton and Frame Integraton Danya Suga and Takayuk Hamamoto Graduate School of Engneerng, Tokyo Unversty of Scence, 6-3-1, Nuku, Katsuska-ku,
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 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 informationComparison of h-, p- and hp-adaptation for convective heat transfer
Computatonal Methods and Expermental Measurements XIII 495 Comparson of h-, p- and hp-adaptaton for convectve heat transfer D. W. Pepper & X. Wang Nevada Center for Advanced Computatonal Methods, Unversty
More informationSimulation of a Ship with Partially Filled Tanks Rolling in Waves by Applying Moving Particle Semi-Implicit Method
Smulaton of a Shp wth Partally Flled Tanks Rollng n Waves by Applyng Movng Partcle Sem-Implct Method Jen-Shang Kouh Department of Engneerng Scence and Ocean Engneerng, Natonal Tawan Unversty, Tape, Tawan,
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 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 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 informationA Fast Visual Tracking Algorithm Based on Circle Pixels Matching
A Fast Vsual Trackng Algorthm Based on Crcle Pxels Matchng Zhqang Hou hou_zhq@sohu.com Chongzhao Han czhan@mal.xjtu.edu.cn Ln Zheng Abstract: A fast vsual trackng algorthm based on crcle pxels matchng
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 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 informationStructured Grid Generation Via Constraint on Displacement of Internal Nodes
Internatonal Journal of Basc & Appled Scences IJBAS-IJENS Vol: 11 No: 4 79 Structured Grd Generaton Va Constrant on Dsplacement of Internal Nodes Al Ashrafzadeh, Razeh Jalalabad Abstract Structured grd
More informationJournal of Computational Physics
Journal of Computatonal Physcs xxx (2011) xxx xxx Contents lsts avalable at ScenceDrect Journal of Computatonal Physcs journal homepage: www.elsever.com/locate/jcp A novel mesh regeneraton algorthm for
More informationREFRACTION. a. To study the refraction of light from plane surfaces. b. To determine the index of refraction for Acrylic and Water.
Purpose Theory REFRACTION a. To study the refracton of lght from plane surfaces. b. To determne the ndex of refracton for Acrylc and Water. When a ray of lght passes from one medum nto another one of dfferent
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 informationA MOVING MESH APPROACH FOR SIMULATION BUDGET ALLOCATION ON CONTINUOUS DOMAINS
Proceedngs of the Wnter Smulaton Conference M E Kuhl, N M Steger, F B Armstrong, and J A Jones, eds A MOVING MESH APPROACH FOR SIMULATION BUDGET ALLOCATION ON CONTINUOUS DOMAINS Mark W Brantley Chun-Hung
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 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 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 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 informationContent Based Image Retrieval Using 2-D Discrete Wavelet with Texture Feature with Different Classifiers
IOSR Journal of Electroncs and Communcaton Engneerng (IOSR-JECE) e-issn: 78-834,p- ISSN: 78-8735.Volume 9, Issue, Ver. IV (Mar - Apr. 04), PP 0-07 Content Based Image Retreval Usng -D Dscrete Wavelet wth
More informationFEDSM-ICNMM
Proceedngs of the of ASME 2010 3rd Jont US-European Fluds Engneerng Summer Meetng and 8th Internatonal Conference on Nanochannels, Mcrochannels, and Mnchannels FEDSM2010-ICNMM2010 FEDSM-ICNMM2010 August
More information3D numerical simulation of tsunami runup onto a complex beach
Proc. Int. Symp. on Fluval and Coastal Dsasters -Copng wth Extreme Events and Regonal Dversty-, Dsaster Preventon Research Insttute, Kyoto Unv. (2005), pp. 221-228 3D numercal smulaton of tsunam runup
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 informationA gradient smoothing method (GSM) for fluid dynamics problems
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluds 2008; 58:1101 1133 Publshed onlne 27 March 2008 n Wley InterScence (www.nterscence.wley.com)..1788 A gradent smoothng method
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 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 informationCAD and CAE Analysis for Siphon Jet Toilet
Avalable onlne at www.scencedrect.com Physcs Proceda 19 (2011) 472 476 Internatonal Conference on Optcs n Precson Engneerng and Nanotechnology 2011 CAD and CAE Analyss for Sphon Jet Tolet Yuhua Wang a
More informationBarycentric Coordinates. From: Mean Value Coordinates for Closed Triangular Meshes by Ju et al.
Barycentrc Coordnates From: Mean Value Coordnates for Closed Trangular Meshes by Ju et al. Motvaton Data nterpolaton from the vertces of a boundary polygon to ts nteror Boundary value problems Shadng Space
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 information(This is a sample cover image for this issue. The actual cover is not yet available at this time.)
(Ths s a sample cover mage for ths ssue. The actual cover s not yet avalable at ths tme.) Ths artcle appeared n a journal publshed by Elsever. The attached copy s furnshed to the author for nternal non-commercal
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 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 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 informationAP PHYSICS B 2008 SCORING GUIDELINES
AP PHYSICS B 2008 SCORING GUIDELINES General Notes About 2008 AP Physcs Scorng Gudelnes 1. The solutons contan the most common method of solvng the free-response questons and the allocaton of ponts for
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 informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
Ths artcle appeared n a journal publshed by Elsever. The attached copy s furnshed to the author for nternal non-commercal research and educaton use, ncludng for nstructon at the authors nsttuton and sharng
More informationThe Research of Ellipse Parameter Fitting Algorithm of Ultrasonic Imaging Logging in the Casing Hole
Appled Mathematcs, 04, 5, 37-3 Publshed Onlne May 04 n ScRes. http://www.scrp.org/journal/am http://dx.do.org/0.436/am.04.584 The Research of Ellpse Parameter Fttng Algorthm of Ultrasonc Imagng Loggng
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 informationUrbaWind, a Computational Fluid Dynamics tool to predict wind resource in urban area
UrbaWnd, a Computatonal Flud Dynamcs tool to predct wnd resource n urban area Karm FAHSSIS a, Gullaume DUPONT a, Perre LEYRONNAS a a Meteodyn, Nantes, France Presentng Author: Karm.fahsss@meteodyn.com,
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 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 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 informationAPPLICATION OF A COMPUTATIONALLY EFFICIENT GEOSTATISTICAL APPROACH TO CHARACTERIZING VARIABLY SPACED WATER-TABLE DATA
RFr"W/FZD JAN 2 4 1995 OST control # 1385 John J Q U ~ M Argonne Natonal Laboratory Argonne, L 60439 Tel: 708-252-5357, Fax: 708-252-3 611 APPLCATON OF A COMPUTATONALLY EFFCENT GEOSTATSTCAL APPROACH TO
More informationNUMERICAL SIMULATION OF WATER ENTRY OF WEDGES BASED ON THE CIP METHOD
Journal of Marne Scence and Technology, Vol., No., pp. -5 (5) DOI:.69/JMST--6- NUMERICAL SIMULATION OF WATER ENTRY OF WEDGES BASED ON THE METHOD Zhao-Yu We,, He-ye Xao, Zh-dong Wang, and Xu-Hua Sh Key
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 informationParallel matrix-vector multiplication
Appendx A Parallel matrx-vector multplcaton The reduced transton matrx of the three-dmensonal cage model for gel electrophoress, descrbed n secton 3.2, becomes excessvely large for polymer lengths more
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 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 informationFinite Element Analysis of Rubber Sealing Ring Resilience Behavior Qu Jia 1,a, Chen Geng 1,b and Yang Yuwei 2,c
Advanced Materals Research Onlne: 03-06-3 ISSN: 66-8985, Vol. 705, pp 40-44 do:0.408/www.scentfc.net/amr.705.40 03 Trans Tech Publcatons, Swtzerland Fnte Element Analyss of Rubber Sealng Rng Reslence Behavor
More informationA Fast Content-Based Multimedia Retrieval Technique Using Compressed Data
A Fast Content-Based Multmeda Retreval Technque Usng Compressed Data Borko Furht and Pornvt Saksobhavvat NSF Multmeda Laboratory Florda Atlantc Unversty, Boca Raton, Florda 3343 ABSTRACT In ths paper,
More informationLoad Balancing for Hex-Cell Interconnection Network
Int. J. Communcatons, Network and System Scences,,, - Publshed Onlne Aprl n ScRes. http://www.scrp.org/journal/jcns http://dx.do.org/./jcns.. Load Balancng for Hex-Cell Interconnecton Network Saher Manaseer,
More informationSPH and ALE formulations for sloshing tank analysis
Int. Jnl. of Multphyscs Volume 9 Number 3 2015 209 SPH and ALE formulatons for sloshng tank analyss Jngxao Xu 1, Jason Wang 1 and Mhamed Soul*, 2 1 LSTC, Lvermore Software Technology Corp. Lvermore CA
More informationType-2 Fuzzy Non-uniform Rational B-spline Model with Type-2 Fuzzy Data
Malaysan Journal of Mathematcal Scences 11(S) Aprl : 35 46 (2017) Specal Issue: The 2nd Internatonal Conference and Workshop on Mathematcal Analyss (ICWOMA 2016) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES
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 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 informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
Ths artcle appeared n a journal publshed by Elsever. The attached copy s furnshed to the author for nternal non-commercal research and educaton use, ncludng for nstructon at the authors nsttuton and sharng
More informationJournal of Computational Physics
Journal of Computatonal Physcs 233 (2013) 34 65 Contents lsts avalable at ScVerse ScenceDrect Journal of Computatonal Physcs journal homepage: www.elsever.com/locate/jcp Fnte element smulaton of dynamc
More informationSteps for Computing the Dissimilarity, Entropy, Herfindahl-Hirschman and. Accessibility (Gravity with Competition) Indices
Steps for Computng the Dssmlarty, Entropy, Herfndahl-Hrschman and Accessblty (Gravty wth Competton) Indces I. Dssmlarty Index Measurement: The followng formula can be used to measure the evenness between
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 informationComplex Numbers. Now we also saw that if a and b were both positive then ab = a b. For a second let s forget that restriction and do the following.
Complex Numbers The last topc n ths secton s not really related to most of what we ve done n ths chapter, although t s somewhat related to the radcals secton as we wll see. We also won t need the materal
More informationOn Some Entertaining Applications of the Concept of Set in Computer Science Course
On Some Entertanng Applcatons of the Concept of Set n Computer Scence Course Krasmr Yordzhev *, Hrstna Kostadnova ** * Assocate Professor Krasmr Yordzhev, Ph.D., Faculty of Mathematcs and Natural Scences,
More informationDynamic wetting property investigation of AFM tips in micro/nanoscale
Dynamc wettng property nvestgaton of AFM tps n mcro/nanoscale The wettng propertes of AFM probe tps are of concern n AFM tp related force measurement, fabrcaton, and manpulaton technques, such as dp-pen
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 informationON SOME ENTERTAINING APPLICATIONS OF THE CONCEPT OF SET IN COMPUTER SCIENCE COURSE
Yordzhev K., Kostadnova H. Інформаційні технології в освіті ON SOME ENTERTAINING APPLICATIONS OF THE CONCEPT OF SET IN COMPUTER SCIENCE COURSE Yordzhev K., Kostadnova H. Some aspects of programmng educaton
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 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 informationMultiblock method for database generation in finite element programs
Proc. of the 9th WSEAS Int. Conf. on Mathematcal Methods and Computatonal Technques n Electrcal Engneerng, Arcachon, October 13-15, 2007 53 Multblock method for database generaton n fnte element programs
More informationCluster Analysis of Electrical Behavior
Journal of Computer and Communcatons, 205, 3, 88-93 Publshed Onlne May 205 n ScRes. http://www.scrp.org/ournal/cc http://dx.do.org/0.4236/cc.205.350 Cluster Analyss of Electrcal Behavor Ln Lu Ln Lu, School
More informationMode III fracture mechanics analysis with Fourier series method
Mode III fracture mechancs analyss wth Fourer seres method F W M Kwok 3 L C Kang, and C R Steele Dvson of ppled Mechancs, Stanford Unversty, Stanford C 94305 numercal scheme based on the Fourer seres method
More information3D Virtual Eyeglass Frames Modeling from Multiple Camera Image Data Based on the GFFD Deformation Method
NICOGRAPH Internatonal 2012, pp. 114-119 3D Vrtual Eyeglass Frames Modelng from Multple Camera Image Data Based on the GFFD Deformaton Method Norak Tamura, Somsangouane Sngthemphone and Katsuhro Ktama
More informationSolutions to Programming Assignment Five Interpolation and Numerical Differentiation
College of Engneerng and Coputer Scence Mechancal Engneerng Departent Mechancal Engneerng 309 Nuercal Analyss of Engneerng Systes Sprng 04 Nuber: 537 Instructor: Larry Caretto Solutons to Prograng Assgnent
More informationTowards sibilant /s/ modelling: preliminary computational results
Acoustcs 8 Pars Towards sblant /s/ modellng: prelmnary computatonal results X. Grandchamp a, A. Van rtum a, X. Pelorson a, K. Nozak b and S. Shmoo b a Département Parole & Cognton, GIPSA-lab, 46, avenue
More informationSixth-Order Difference Scheme for Sigma Coordinate Ocean Models
064 JOURNAL OF PHYSICAL OCEANOGRAPHY VOLUME 7 Sxth-Order Dfference Scheme for Sgma Coordnate Ocean Models PETER C. CHU AND CHENWU FAN Department of Oceanography, Naval Postgraduate School, Monterey, Calforna
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 informationHelsinki University Of Technology, Systems Analysis Laboratory Mat Independent research projects in applied mathematics (3 cr)
Helsnk Unversty Of Technology, Systems Analyss Laboratory Mat-2.08 Independent research projects n appled mathematcs (3 cr) "! #$&% Antt Laukkanen 506 R ajlaukka@cc.hut.f 2 Introducton...3 2 Multattrbute
More informationIntroduction to Geometrical Optics - a 2D ray tracing Excel model for spherical mirrors - Part 2
Introducton to Geometrcal Optcs - a D ra tracng Ecel model for sphercal mrrors - Part b George ungu - Ths s a tutoral eplanng the creaton of an eact D ra tracng model for both sphercal concave and sphercal
More informationDetermining the Optimal Bandwidth Based on Multi-criterion Fusion
Proceedngs of 01 4th Internatonal Conference on Machne Learnng and Computng IPCSIT vol. 5 (01) (01) IACSIT Press, Sngapore Determnng the Optmal Bandwdth Based on Mult-crteron Fuson Ha-L Lang 1+, Xan-Mn
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 informationDesign for Reliability: Case Studies in Manufacturing Process Synthesis
Desgn for Relablty: Case Studes n Manufacturng Process Synthess Y. Lawrence Yao*, and Chao Lu Department of Mechancal Engneerng, Columba Unversty, Mudd Bldg., MC 473, New York, NY 7, USA * Correspondng
More informationExplicit Formulas and Efficient Algorithm for Moment Computation of Coupled RC Trees with Lumped and Distributed Elements
Explct Formulas and Effcent Algorthm for Moment Computaton of Coupled RC Trees wth Lumped and Dstrbuted Elements Qngan Yu and Ernest S.Kuh Electroncs Research Lab. Unv. of Calforna at Berkeley Berkeley
More informationOutline. Midterm Review. Declaring Variables. Main Variable Data Types. Symbolic Constants. Arithmetic Operators. Midterm Review March 24, 2014
Mdterm Revew March 4, 4 Mdterm Revew Larry Caretto Mechancal Engneerng 9 Numercal Analyss of Engneerng Systems March 4, 4 Outlne VBA and MATLAB codng Varable types Control structures (Loopng and Choce)
More informationTEST-05 TOPIC: OPTICS COMPLETE
Q. A boy s walkng under an nclned mrror at a constant velocty V m/s along the x-axs as shown n fgure. If the mrror s nclned at an angle wth the horzontal then what s the velocty of the mage? Y V sn + V
More informationNumerical simulation of flow past twin near-wall circular cylinders in tandem arrangement at low Reynolds number
Water Scence and Engneerng 2015, 8(4): 315e325 HOSTED BY Avalable onlne at www.scencedrect.com Water Scence and Engneerng journal homepage: http://www.waterjournal.cn Numercal smulaton of flow past twn
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 informationFlow over Broad Crested Weirs: Comparison of 2D and 3D Models
Journal of Cvl Engneerng and Archtecture 11 (2017) 769-779 do: 10.17265/1934-7359/2017.08.005 D DAVID PUBLISHING Flow over Broad Crested Wers: Comparson of 2D and 3D Models Shaymaa A. M. Al-Hashm 1, Huda
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 information