Extension of Smoothed Particle Hydrodynamics (SPH), Mathematical Background of Vortex Blob Method (VBM) and Moving Particle Semi-Implicit (MPS)
|
|
- Julius Hutchinson
- 6 years ago
- Views:
Transcription
1 Amercan Journal of Computatonal athematcs, 04, 5, Publshed Onlne December 04 n ScRes. Extenson of Smoothed Partcle Hydrodynamcs (SPH), athematcal Background of Vortex Blob ethod (VB) and ovng Partcle Sem-Implct (PS) Hrosh Isshk Insttute of athematcal Analyss (IA), Osaka, Japan Emal: sshk@dab.h-ho.ne.p Receved October 04; revsed 3 November 04; accepted 4 November 04 Copyrght 04 by author and Scentfc Research Publshng Inc. Ths work s lcensed under the Creatve Commons Attrbuton Internatonal Lcense (CC BY). Abstract SPH has a reasonable mathematcal background. Although VB and PS are smlar to SPH, ther mathematcal backgrounds seem fragle. VB has some problems n treatng the vscous dffuson of vortces but s known as a practcal method for calculatng vscous flows. The mathematcal background of PS s also not suffcent. Not wth standng, the numercal results seem reasonable n many cases. The problem common n both VB and PS s that the space dervatves necessary for calculatng vscous dffuson are not estmated reasonably, although the treatment of advecton s mathematcally correct. Ths paper dscusses a method to estmate the above mentoned problem of how to treat the space dervatves. The numercal results show the comparson among FD (Fnte Dfference ethod), SPH and PS n detal. In some cases, there are bg dfferences among them. An extenson of SPH s also gven. Keywords SPH, PS, VB, Concentraton of ass or Vortcty, Estmaton of Space Dervatves. Introducton The author has once shown how to obtan the space dervatve n an rregular mesh usng the movng least square method []. A smlar problem s dscussed from a dfferent vewpont n SPH. Gngold & onaghan [] and Lucy [3] have developed Smooth Partcle Hydrodynamcs method (SPH). In SPH, the contnuous mass dstrbuton s approxmated by the fnte number of partcles. Namely, the contnuous quanttes are represented by How to cte ths paper: Isshk, H. (04) Extenson of Smoothed Partcle Hydrodynamcs (SPH), athematcal Background of Vortex Blob ethod (VB) and ovng Partcle Sem-Implct (PS). Appled athematcs, 5,
2 H. Isshk the fnte number of the dscrete quanttes. At the same tme, the space dervatves of the dstrbuted quanttes are expressed algebracally by the dscrete quanttes. Then, they transform the contnuous system nto a dscrete system convenent for the numercal soluton of the ntal and boundary value problem. Vortex Blob ethod (VB) s one of the practcal numercal tools for hgh Reynolds number flows [4]. The abducton s approxmated reasonably. However, the dffuson of vortces can t be approxmated precsely. We have shown that ths problem can be solved, f we apply the deas developed by SPH. ovng Partcle Sem-mplct method (PS) uses a smlar dscretzaton of the ntal and boundary value problem as SPH and VB. However, the mathematcal background s not suffcent. Not wth standng, the numercal results seem reasonable n many cases [5]. Ths suggests us that we may fnd the mathematcal background [] [6]. We have shown an nterestng relatonshp between SPH and PS as far as the gradent operator s concerned. The numercal results show the comparson among FD, SPH and PS n detal. In some cases, there are bg dfferences among them. A classfcaton of numercal solutons s shown n Fgure. The weak solutons: FE, FV, IR and GIR are Fnte Element ethod, Fnte Volume ethod, Integral Representaton ethod and Generalzed Integral Representaton ethod, respectvely. The strong solutons: FD, SPH, PS, LB and Col are Fnte Dfference ethod, Smooth Partcle Hydrodynamcs, ovng Partcle Sem-Implct ethod, Voltex Blob ethod, Lattce Boltzmann ethod and Collocaton ethod, respectvely. The general characterstcs of the varous numercal solutons are shown n Table. Snce SPH, VB and PS belong to the strong solutons and are mesh-less methods, the low computatonal cost becomes very mportant n some cases.. Dscretzaton Used n SPH and the Extenson α 3 Let x = x eα = x e+ x e + x e3 = x+ y+ zk be a poston vector, where the suffx n Greek letter s used to refer to the component of the coordnates, and the summaton conventon s used for the Greek suffxes. We defne a functon η ( x ) as an ntegral of w( x x ) ξ( x ) ρ( x ) wth respect to a volume V, where w( x ), ρ x are a weght functon, an auxlary functon and the densty of a flud, respectvely: ξ ( x ) and The functon η ( x ) and = ( ) d η x x x ξ x ρ x V. () V w ξ x can be scalars, vectors and tensors. If we ntroduce an dscretzaton of the volume V and the densty ρ as V = V, = 0,,, N, () Fgure. Classfcaton of numercal solutons. Table. General characterstcs of weak and strong solutons ( : very good, : good, : not good). Week method Strong method Stablty Precson Computatonal cost (memory, tme) 45
3 H. Isshk where the suffx n Roman letter s used to dscrmnate the pont, and V. We have an approxmaton of Equaton (): d and Furthermore, f we assume for the weght w ρ ( x)d V =, = 0,,, N, (3) = w( ) ξ s the mass of the volume element η x x x x. (4) w ( x )d V = (5a) V w = 0 for x > σ. (5b) If σ tends to 0, w( x ) tends to δ ( x ), where δ ( x ) s Drac s delta functon, namely:, when 0. Usng Equaton (6), we have a followng approxmaton: From Equatons (4) and (7), we obtan w x δ x σ (6) = ( ) ( ) ( ) d V = η x x x ξ x ρ r ξ x ρ x. (7) V w = w( ) If we assume ξ ( x ) =, then, we have from Equaton (8) ξ x ρ x x x ξ x. (8) = ( ) ρ( ) = ( ) ρ x w x x x V w x x. (9) V d Substtutng x = x, we obtan w ρ x = x x. (0) Hence, we have w( 0) = w( 0d ) V ρ x ρ x ρ x. () Ths means that Equaton (9) s a plausble approxmaton. Substtutng Equaton (9) nto Equaton (8), we obtan w( x x) ξ ( x) w( x x) ρ ( x) ξ ξ x = = w x x x. () Settng x to x n Equaton (), we have usng Equaton (0) ξ( x ) = w x x x ρ ( x ) ξ. (3) Operatng x α on both sdes of Equaton (8), we obtan ( x) ( x) w( ) ξ ρ ρ( x) + ξ( x) = x x ξ( x ). (4) x x x α α α 46
4 H. Isshk If we assume We derve w d ( ) = w( ) x x x x. (5) w x x w x x α α x x x x = = w α α ( x x ) x d x x x x x, (6) ( x) ( x x ) ρ w x x = α = x α α w ( ) α x x x x x. (7) Substtutng Equatons (9), (6) and (7) nto Equaton (4), we obtan ξ ( x) x Settng x to x n Equaton (8), we have ξ x α ( x) α α α x x = ρ w ( x x ) ( ξ( x) ξ( x )) ( x) x x α α x x = ρ ρ x w ( x x ) ( ξ( x ) ξ( x )) ( x ) x x ( x) α α α w ( x x ) x x = where we assume w ( 0) = 0. Now, we derve formulas for second dervatves. Frstly, from Equaton (4), we have ( x x ) x x ξ x ξ x ρ x ξ x ρ x ρ x ρ( x) ξ( x) x x x x x x x x = w x x ( x ).. (8). (9), (0) α β β α α β α β α β ξ () From Equatons (6) and (7), we obtan ( x x ) w β x x = w α β α ( x x ) x x x x x β x x w ( x x ) x x x x x x x x = w + ( x) w( x x ) ρ = x x x x α β α β α α β β αβ α α β β δ, 3 x x x x x x ( x x )( x x ) x x x x = w x x + w ( x x ) α α β β αβ α α β β δ, 3 x x x x x x () (3) αβ αβ αβ where δ s Chronecker s delta: δ = and δ = 0 when α = β and α β, respectvely. Substtutng Equatons (9), (7), () and (3) nto Equaton (), we derve 47
5 H. Isshk ( x x ) ξ x w ξ x ρ x ξ x ρ x ρ x ρ( x) = ξ( x) ξ( x) x x x x x x x x x x α β α β β α α β α β ( x x )( x x ) x x x x = w x x + w ( x x ) x α α β β αβ α α β β δ ξ 3 x x x x x x x x x x + w ( x x ) ( ξ( x) ξ( x) ) w ( x x ) ρ ( x) x x x x β β α α x x x x + w w ρ x x x x x ( x x ) ( ξ( x ) ξ( x )) ( x x ) α α β β ( x x )( x x ) x x x x ξ x w x x + w ( x x ) α α β β αβ α α β β δ 3 x x x x x x ( α α )( β β ) α α β β x x x x αβ δ x x x x = w x x + w ( x x ) 3 ξ x ξ x x x x x x x ( ) ( β β )( α α ) ( α α )( β β x ) x x xk + x x x x k + w x xk w x x ( ξ( x) ξ( x) ) k. ρ x x x x x Settng x to k k ( x x ) x n Equatons (), (3) and (4), we have δ = w x x + w ( x x ) x x α α β β α α β β w x x x x αβ x x x x α β 3 x x x x x x ( α α )( β β ) α α β β x x x x αβ x x x x ρ x δ = w + w ( ) α β x x x x x x ( x) ξ α x x β ρ 3 x x x x x x ( x ) k (4), (5) ( α α )( β β ) α α β β x x x x αβ δ x x x x = w w x x + x x ( ξ( x) ξ( x )) 3 x x x x x x k k k, (6) ( β β )( α α ) ( α α )( β β x ) x x xk + x x x x k + w x x w x x ( ξ( x ) ξ( x )), ρ x x x x x where we assume w ( 0) s fnte. From Equaton (9), we obtan for scalar φ and vector a x x φ = φ φ ρ [ ] [ ] w ( x x ) ( ( x ) ( x )) ( x ) x x ( x ) ( ) x a x a x w ( x ) x x a = x x ρ From Equaton (7), we obtan for scalar φ (7), (8). (9) 48
6 H. Isshk φ ( x) ρ( x ) k d ( x x ) w ( x x ) ( φ( x ) φ( x )) ( α α )( α α ) α α α α x x x x d x x x x = w w x x + x x ( φ( x) φ( x )) 3 x x x x x x ( α α )( α α ) ( α α )( α α x ) x x xk + x x x x k + w x x w x x ( φ( x ) φ( x )) ρ x x x x x k k k = w + x x ( x ) x x xk + w x xk w x x ( φ( x) φ( x) ) k, ρ x x x x x k k where d s the number of the dmenson. 3. Applcaton to Vortex Blob ethod (VB) The basc equatons for an ncompressble vscous flow are gven by Naver-Stokes equaton [7]: or ρ u + t ρ = + µ where µ s the coeffcent of vscosty, and the contnuty equaton: The pressure p satsfy (30) u u p u, (3) du ρ = p + µ u, (3) dt u =0. (33) The vscous flow s determned by Equatons (3) or (3), (33) and (34). If we ntroduce vortcty ω, the basc equatons can be rewrtten as Equaton (36) can be wrtten as p = ρ u u. (34) ω = u, u= 0, (35) ω ρ + ρu ω = ρ ω u+ µ ω. (36) t dω = ω u + ν ω, (37) dt where ν = µ ρ s the knematc vscosty. If ω s determned, the velocty u s obtaned by an ntegral representaton [7]: ( ) S G( x ξ t) ω( ξ t) Vξ + u ( x ξ t) ε u x, t = G x, ξ, t u ξ, t n G x, ξ, t n u ξ, t ds V ξ ξ ξ,,, d,,, where ε =,, 0 when x V, x S, x V S, respectvely, and (, ) δ (, ) (38) xg x ξ = x ξ, (39a) 49
7 H. Isshk G ( x ξ) x ξ ( x ξ ) n D, = π ln n D ( 4π) x ξ n 3D ( x ξ) x ξ n D xg( x, ξ) = G( x, ξ) = G ( x, ξ) e = ( π)( x ξ) x ξ n D 3 ( 4π)( x ξ) x ξ n 3D The second term on the rght hand sde of Equaton (38) s usually called Bo-Savart law. Vortex Blob ethod (VB) uses Equatons (37) and (38) and dscretzes the vortex feld as an assembly of concentrated vortces. The poston of the concentrated vortces are determned by ( t) d x a, dt (, t) (39b) (39c) = ua, (40) where a s the Lagrangan coordnates. It becomes very mportant how to obtan ω u and ω. In ths problem, we use the vortcty ω corresponds to the densty ρ n Secton. Then, we consder nstead of Equaton () Equaton (3) s replaced by Then, we have nstead of Equatons (8) and (9) = ( ) d η x x x ξ x ω x V. (4) V w ω ( x )d V =Ω, = 0,,, N. (4) w ξ x ω x = x x ξ x Ω. (43) w( ) ω( ) V w V d ω x = x x x = x x Ω. (44) In all equatons n Secton, f we replace ρ and wth ω and Ω, respectvely, we obtan the dfferental formulas n Secton athematcal Background of ovng Partcle Sem-Implct ethod (PS) From Equatons (9), (8) and (30), we have ( x t) α β β u, x x α α w ( u (, t) u (, t) ) β = x x x x x ρ ( x ) x x [ ] x x p = ( p( x ) p ( x )) w ( x x ) ρ x x x = + d ( x, ) ρ ( x ) ( x x ) ( x x ) ( ( x, ) ( x, )) k, (45), (46) α α α u t w w u t u t x x ( x ) x x xk + w x xk w x x ( φ( x) φ( x) ) k. ρ x x x x x k If the weght wr satsfes wr d d ~ wr and δ ( δ ) ( d ) π wrr dr dr r + =, (48) 0 (47) 40
8 H. Isshk then, we obtan where c( d ) s a constant. Then, we have wr ~ c d r, (49) wr cd wr cd ~ ~ wr and wr 3. dr r r dr = r = r (50) d d Substtutng Equatons (49) and (50) nto Equatons (45), (46) and (47), we obtan ( x t) α β β u, x x α α w ( u (, t) u (, t) ) β x x x x x ρ ( x ) x x [ ] x x p ( p( x ) p ( x )) w( x x ) ρ x x x d 3 ( x, ) ρ ( x ) ( x x ) ( ( x, ) ( x, )) u t w u t u t α α α x x, (5), (5) ( x ) x x xk α α + w k w ( u u ) k. ρ ( ) x x x x x x x k x x x x Now, we consder a weght k wr wth a small parameter 0 < ε = r and constant α : wr r e α = r+ ε + ε re 0 ( re r), where r = x. The weght wr has an asymptotc form: Ths asymptotc form s smlar to e ( 0 r r ) re re wr α ~ re when r re. r r α ε ε e r α = = r r ε < + + e k r defned by Equaton (4) n Ref. [5]. If e wr satsfes Equaton (49). In ths case, Equaton (5) has smlar forms as gven by Equaton () n Ref. [5], f we assume =, = 0,,. The author has once dscussed ths problem from a dfferent angle [] and had the concluson that the dscrete dfferental operators used n PS, especally, Laplace operator don t have the strct background from the mathematcal vewpont. Those operators should be consdered to be a knd of expermental formulas that s derved numercally. If we apply them to rregular grds, the results vary rregularly. We should say the operator estmate the dervatves statstcally. Qute naturally, the results are not unque. It may gve a good estmate n one tme and a wrong result n the other tme. Hence, f we need to verfy the accuracy, we should obtan several numercal results of the same problem usng the varous dscretzatons of the regon and take the statstcal average. Although we do not deny ths knd of approach, we wsh to ensure the relablty. For the purpose, we need much dscusson. Recently, a new paper [6] dscusses the accuracy of Laplace operator n PS n detal. Hopefully, we wsh to prove mathematcally, f possble, that, f we decrease the grd sze zero, then, the numercal error approaches stably to zero. e r s bg, 5. Numercal Verfcaton of Dscretzed Dfferental Operators of SPH wth Gaussan Kernel Numercal calculatons were conducted to verfy the valdty of the dscretzed dfferental operators such as (53) (54) (55) 4
9 H. Isshk Equatons (8) and (30). For smplcty, one-dmensonal (D) cases were consdered. The frst and second de- φ x : rvates of a scalar functon were calculated usng Equatons (8) and (30), respectvely. φ x = sn π xl n 0 x L (56) The regon s dvded nto N ntervals. Frst, we consder the unform mesh. Let dx and x ( = 0,,, N ) be the length of the nterval and the mdpont of the nterval, respectvely: 0.5dx0 for = 0 x =, dv = dx and = ρdv = ρd x. x + 0.5( dx + dx) for =,, As the weght functon wr, we use Gaussan functon: (57), (58) The weght functon wr satsfes wr r = exp. (59) πγ γ 0 w x dx =. (60) The D verson of Equatons (8) and (30) are gven by dφ x x w x x x x dx = ρ x x ( φ φ) ( x ), (6) = ( )( ) d φ x w x x x x ρ φ φ dx ( x ) x x xk + w x xk w x x ( φ( x) φ( x) ) k. ρ x x x x x k k (6) 5.. Unform Densty and Regular esh The densty ρ and the length of element dx are gven by where the computatonal regon s defned as 0 < x< L Verfcaton of Gradent Operator From Equatons (6), (63) and (45), we have ρ = = 0,,, N, (63a) dx= LN= 0,,, N, (63b) dφ x x w ( x x ) ( φ( x) φ( x) ) dx =. (64) x x We used parameters: L = 4, N = 80 and γ = dx = The numercal results are shown n Fgure. Although there exst large errors at the boundares as shown n Fgure (a), the numercal results agree very well wth the exact except the neghborhood of the boundares. If we extend the regon beyond the boundares, the estmatons wthn the orgnal regon are mproved as shown n Fgure (b) Verfcaton of Laplace Operator From Equaton (6), we have 4
10 H. Isshk d φ dx ( x) Fgure. A comparson between calculated and exact values of dφ/dx ((a) Not usng extenson of regon; (b) Usng extenson of regon). ρ( x ) = w ( x x ) ( φ( x ) φ( x )) x x x x + w x x w x xk x x k. ρ k ( φ φ) ( x ) k x x x x k The numercal results are shown n Fgure 3. If we extend the regon beyond the boundares, the estmatons wthn the orgnal regon are mproved. 5.. Non-Unform Densty and Regular esh The densty ρ and the length of element dx are gven by (65) ρ = + x L = 0,,, N, (66a) d x= LN, = 0,,, N, (66b) respectvely. The other parameters are same as n Secton 5... The numercal results are shown n Fgure 4. The numercal result agrees well wth the exact result Unform Densty and Non-Unform esh The densty ρ and the length of element dx are gven by ρ = = 0,,, N, (67a) d (.05) 0,,,, = N L L x = α = N α (.05 ) = L, (67b) N N respectvely. The other parameters are same as n Secton 5... The numercal results are shown n Fgure 5. The numercal result agrees well wth the exact result Non-Unform Densty and Non-Unform esh The densty ρ and the length of element dx are gven by = 0 ρ = + x L = 0,,, N, (68a) 43
11 H. Isshk Fgure 3. A comparson between calculated and exact values of d φ/dx ((a) Not usng extenson of regon; (b) Usng extenson of regon). Fgure 4. A comparson between calculated and exact values ((a) Gradent operator dφ/dx; (b) Laplace operator d φ/dx ). Fgure 5. A comparson between calculated and exact values ((a) Gradent operator dφ/dx; (b) Laplace operator d φ/dx ). 44
12 H. Isshk = N L L dx = α (.05) = 0,,, N, α (.05) = L (68b) N N respectvely. The other parameter are same as n Secton 5... The numercal results are shown n Fgure 6. The numercal result agrees well wth the exact result. 6. Comparson between SPH and PS 6.. D Formulas for Dscrete Gradent and Laplacan Operator For convenence, we summarze the D dscrete dfferental operators used n SHP and PS as follows. From Equatons (6) and (6), we have for SHP = 0 x x w x x x x dφ ( ) dx = φ φ ρ ( x ) x x ( x) w ( x x )( ( x) ( x) ) φ ρ = φ φ k ( x ) x x xk + w x xk w x x ( φ( x) φ( x) ) k, ρ x x x x x From Refs. [5] and [8], we have for PS where k, (69a) (69b) = ρ ( x )dx. (69c) dφ φ φ = ( x x) w( x x ), (70a) dx n x x d φ = dx λ n λ = ( φ φ) w( x x ), (70b) n = w( x x ), (7a) ( ) ( x ) x x w x x w x. (7b) Fgure 6. A comparson between calculated and exact values ((a) Gradent operator dφ/dx; (b) Laplace operator d φ/dx ). 45
13 H. Isshk 6.. Effects of esh and Weght on Estmaton of Gradent and Laplacan of a Gven Functon We summarze the meshes and weghts used n the followng n Table and Table 3, respectvely. The parameters used n Sectons are gven below. L = 4, N = 40, ρ =, γ = dx for SPH0, h = 8dx for SPH and SPH, and dx for PS, N = 8 for SPH0, SPH and SPH, and for PS. end _ cor In Sectons , the computatonal results are shown n Fgures 7-. The results are summarzed Table 4 and Table 5 n Secton Regular esh As shown n Fgure 7, there are no bg errors n all methods SPH0, SPH, SPH and PS. The accuracy of SPH0 s very hgh. Table. Classfcaton of mesh. Name dx L Regular dx = N Algebrac L L dx α = N = + N 00, dx = L = 0 Geometrc dx L N α.05 =, = N = 0 dx = L L = +, N Random dx α ( 0.065drand ) L Snusodal dx = α 0.5sn 4π + N N, = N = 0 = N = 0 dx = L dx = L Table 3. Classfcaton of weght. Name Weght SPH0 Gauss w( x) SPH Lucy w( x ) SPH Opt.kernel PS w( x ) x = exp, x < πγ γ 3 5 x x + 3, x σ = 4 σ σ σ 0, x > σ 3 x x 6 + 6, 0 x < 0.5h h h 3 4 x w x =, 0.5 h x < h 3h h 0, h x h x, 0 x < h = 0, h x 46
14 H. Isshk Table 4. Summary on dφ/dx ( : good, : not good, : bad). Regular Algebrac Geometrc Random Snusodal SPH0 SPH SPH PS Table 5. Summary on d φ/dx ( : good, : not good, : bad). Regular Algebrac Geometrc Random Snusodal SPH0 SPH SPH PS 6... Algebrac esh As shown n Fgure 8, a bg error occurs n Geometrc esh As shown n Fgure 9, a bg error occurs n d φ dx of PS. d φ dx of PS Random esh As shown n Fgure 0, the errors are rather small n all methods SPH0, SPH, SPH and PS not only for dφ dx but also for d φ dx. The errors n PS are surprsngly small Snusodal esh As shown n Fgure, the errors are rather bg n all methods SPH0, SPH, SPH and PS not only for dφ dx but also for d φ dx except dφ dx n PS. The errors n SHP0 are smaller than those n the other methods Summary of Results From the above mentoned numercal results, the followng summares are obtaned Applcaton to Soluton of Intal Value Problem D Flud oton wthout Pressure and Vscosty The moton of vast number of partcles dstrbuted n space under the acton of the gravtatonal force may be treated as a flud moton wthout pressure and vscosty [9] [0]: ρ ρu + = 0, (7a) t x u u + u = Π t x x, (7b) Π = 4πGρ, (7c) x where ρ, u and Π s the densty, velocty and gravtatonal potental. G s the gravtatonal constant. The soluton of the problem defned above by Euleran method s gven as follows: ) At tme t, assume ρ, u and Π are gven. ) ρ t and u t are obtaned from Equatons (7a) and (7b), and Π s obtaned from Equaton (7c). 47
15 H. Isshk dφ dx d φ dx SPH0 SPH SPH PS Fgure 7. Regular mesh. 48
16 H. Isshk dφ dx d φ dx SPH0 SPH SPH PS Fgure 8. Algebrac mesh. 49
17 H. Isshk dφ dx d φ dx SPH0 SPH SPH PS Fgure 9. Geometrc mesh. 430
18 H. Isshk dφ dx d φ dx SPH0 SPH SPH PS Fgure 0. Random mesh. 43
19 H. Isshk dφ dx d φ dx SPH0 SPH SPH PS Fgure. Snusodal mesh. 43
20 H. Isshk 3) ρ, and u at tme t+ dt s calculated. 4) Repeat ths process. In the soluton by Lagrangan method, frst, the equatons n Euleran form are transformed nto Laglangan form: Dx u Dt =, (73a) Du u u + u = Π Dt t x x, (73b) D ρ ρ u u ρ + = ρ Dt t x x, (73c) Π = 4πGρ. (73d) x The followng procedures gve the soluton of the problem defned above by Lagrangan method: ) At tme t, assume ρ, u and Π are gven. ) Dx Dt, Dρ Dt and Du Dt s obtaned from Equatons (73a), (73b) and (73c), and Π s obtaned from Equaton (73d). 3) x, ρ and u of the materal pont at tme t+ dt s calculated. 4) Repeat ths process. If x s obtaned, then, dx and ρ s obtaned as shown below: and dx = ( x+ x) + ( x x ) = ( x+ x ) (74a) ρ =. (74b) dx Hence, from the theoretcal vewpont, ths problem can be solved wthout usng the gradent operator. However, we use the gradent operator to obtan ρ usng the contnuty equaton. ) Trapezodal Dstrbuton of the Intal Densty The ntal condtons are gven by The boundary condtons are specfed as The computatonal condtons are as follows: x = x0 = LN, (75a) 0 + ( x 0.35L) when 0.5L x < 0.35L L when 0.35L x < 0.65L ρ = 0 ( x 0.65L) when 0.65L x < 0.75L L otherwse, (75b) u = 0. (75c) ρ = ρ, ρ = ρ ; (76a) 0 N N u = u, u = u. (76b) 0 N N L = 4, N = 4, ν = 0.0, dt = , t = 30000dt, G = 0.005, γ = dx for SPH0, h = dx for PS, N _ = 8 for SPH0 and for PS. The numercal results are shown n Fgure. The FD (Fnte Dfference ethod) uses the central dffer- end end cor 433
21 H. Isshk ρ n x 0 ρ n x FD SPH0 PS Fgure. Trapezodal dstrbuton of the ntal densty. ences for the frst and second space dervatves, and the precson of FD s consdered hgh. The Euleran soluton s used for FD, and the Lagrangan soluton s used for SPH0 and PS. The dstrbuton pattern of PS s slghtly dfferent from those of FD and SPH0. ) Rectangular Dstrbuton of the Intal Densty wth G = 0.00 The ntal condtons are gven by x= x0 = ( + 0.5) LN, (77a) when x < 0.5L ρ = when 0.5L x 0.75L when 0.75 L x, (77b) 434
22 H. Isshk The boundary condtons are specfed as u = 0. (77c) ρ = ρ, ρ = ρ ; (78a) 0 N N u = u, u = u. (78b) 0 N N The computatonal condtons are as follows: L = 4, N = 4, ν = 0.0, dt = , t = 30000dt, G = 0.00, γ = dx for SPH0, h = dx for PS, N _ = 8 for SPH0 and for PS. The numercal results are shown n Fgure 3. The FD uses the central dfferences for the frst and second space dervatves, and the precson of FD s consdered hgh, f we neglect the spurous oscllaton. The dstrbuton pattern of PS s slghtly dfferent from those of FD and SPH0. end end cor ρ n x 0 ρ n x FD SPH0 PS Fgure 3. Rectangular dstrbuton of the ntal densty wth G =
23 H. Isshk 3) Rectangular Dstrbuton of the Intal Densty wth G = The ntal condtons are gven by The boundary condton are specfed as Computatonal condton x = x0 = LN, (79a) when x < 0.35L ρ = when 0.35L x < 0.75 L, when 0.35L x (79b) u = 0. (79c) ρ = ρ, ρ = ρ ; (80a) 0 N N u = u, u = u. (80b) 0 N N L = 4, N = 4, ν = 0.0, dt = , t = 30000dt, G = 0.005, γ = dx for SPH0, h = dx for PS, N _ = 8 for SPH0 and for PS. The numercal results are shown n Fgure 4. The FD uses the central dfferences for the frst and second space dervatves, and the precson of FD s consdered hgh, f we neglect the spurous oscllaton. The dstrbuton pattern of PS s dfferent from those of FD and SPH D Flud oton wthout Pressure but wth Vscosty In the present secton, the vscosty s ntroduced: end ρ ρu + = 0, t x end cor u + u u = + ν u, Π t x x x Π = 4π Gρ, (8c) x where ν s the knematc vscosty. Snce the dscrete Laplacan operator of SPH0 generates a bg error at the dscontnuty, the ntal densty dstrbuton was smoothened usng the fve pont runnng average, and, n the second example below, an small artfcal numercal vscosty μ = n case of N = 4 was added at every tme step: + dx (8a) (8b) ρ ρ + µ ρ ρ + ρ, (8a) u u + u u + u. (8b) µ + dx In the thrd examples below, a bg dfference has occurred between SPH and PS solutons. ) Exponental Dstrbuton of the Intal Densty The ntal condtons are gven by x = x0 = ( ) LN, (83a) The boundary condtons are specfed as x 0.5L ρ = exp, 0.L (83b) u = 0. (83c) ρ = ρ, ρ = ρ ; (84a) 0 N N 436
24 H. Isshk ρ n x 0 ρ n x FD SPH0 PS Fgure 4. Rectangular dstrbuton of the ntal densty wth G = The computatonal condtons are as follows: u = u, u = u. (84b) 0 N N L = 4, N = 4, ν = 0., dt = , t = 0000dt, γ = dx for SPH0, h = dx for PS, G = 0.005, N _ = 8 for SPH0 and for PS. The numercal results are shown n Fgure 5. The FD uses the central dfferences for the frst and second space dervatves, and the precson of FD s consdered hgh. In ths example, the FD, SPH and PS solutons becomes smlar. ) The Trapezodal Dstrbuton of the Intal Densty The ntal condtons are gven by x = x0 = ( + 0.5) LN, (85a) end end cor 437
25 H. Isshk ρ n x 0 ρ n x FD SPH0 PS Fgure 5. Exponental dstrbuton of the ntal densty. 0 + ( x 0.35L) when 0.5L x < 0.35L L when 0.35L x < 0.65L ρ = 0 ( x 0.65L) when 0.65L x < 0.75L L otherwse, (85b) u = 0, (85c) where the ntal densty dstrbuton was smoothened usng the fve pont runnng average: ( ρ + ρ + ρ + ρ+ + ρ+ ) ρ. (86) 5 438
26 H. Isshk The boundary condtons are specfed as The computatonal condtons are as follows: ρ = ρ, ρ = ρ ; (87a) 0 N N u = u, u = u. (87b) 0 N N L = 4, N = 4, ν = 0.0, dt = , t = 30000dt, G = 0.005, γ = dx for SPH0, h = dx for PS, N _ = 4 for SPH0 and for PS. The numercal results are shown n Fgure 6. The FD uses the central dfferences for the frst and second space dervatves, and the precson of FD s consdered hgh. In ths example, small dfference s observed n the FD, SPH and PS solutons. 3) The Flat-Slope-Flat Dstrbuton of the Intal Velocty end cor end ρ n x 0 ρ n x FD SPH0 PS Fgure 6. Trapezodal dstrbuton of the ntal densty. 439
27 H. Isshk The ntal condtons are gven by x = x0 = LN, (88a) ρ =, (88b) 0.5 when x < 0.5L u = ( x 0.5L) when 0.5L x 0.75L (88c) L 0.5 when 0.75L x The boundary condtons are specfed as ρ =, = ; (89a) The computatonal condtons are as follows: 0 ρn u = 0.5, = 0.5. (89b) 0 un L = 4, N = 4, ν = 0.0, dt = , t = 0000dt, G = 0, γ = dx for SPH0, h = d x for PS, N _ = 8 for SPH0 and for PS. end cor Snce we assume G = 0 n ths example, the equaton becomes Burgers equaton. a) Comparson of FD Soluton wth the Analytcal One The central dfferences were used for the frst and second space dervatves. As shown n Fgure 7. The FD soluton s a very good approxmaton of the analytcal soluton []. b) Comparson of FD, SPH0 and PS Solutons Fgure 8 and Fgure 9 show the comparsons of ρ and u among FD, SPH0 and PS solutons. SPH0 and PS solutons do not gve good approxmatons. Wth respect to the velocty u, the dfference between SPH0 and PS s bg. 7. odfed Gaussan Weght Gaussan-type weghts of fnte support wth C contnuty are gven as follows: end 4 x x x x x αexp exp = α + 4 γ h h h γ w( x) = x πx α exp cos, γ h (90a), (90b) Fgure 7. Comparson of FD soluton wth the analytcal one ((a) FD soluton; (b) Exact soluton). 440
28 H. Isshk ρ n x 0 ρ n x FD SPH0 PS Fgure 8. Comparson of densty ρ. where α satsfes h x x hx α exp 0 exp d d x αh 0 x x h h γ h γ = w( x) dx w( x) dx = h = = 0 h x πx hx πx α exp cos dx α h exp cos d x. 0 0 γ h γ The frst and the second dervatves of w( x ) gven by Equaton (90a) are gven by x x exp 4 4 w x = α + x+ + x, (9a) h γ h γ h γ h γ (9) 44
29 H. Isshk u n x 0 u n x FD SPH0 PS Fgure 9. Comparson of velocty u x exp w x = α x + x + x. (9b) h γ h γ h γ γ h γ h γ h γ Numercal examples are shown n Fgure 0. When h s small, the accuracy becomes low. As has already ponted out n Secton 5.., there exst large errors at the boundares. The numercal results agree well wth the exact ones on overall. In the examples, γ s equal to dx = L/N = 0. and the accuracy seems suffcent when h s bgger than or equal to 4dx. 8. Conclusons The author has once shown how to obtan the partal dervatve n an rregular mesh usng the movng least 44
30 H. Isshk h dφ dx d φ dx Fgure 0. Effect of support h ( φ = sn ( πx L) ; L = 4, N = 40, γ = 0., h = 0., 0.4, 0.8). square method []. A smlar problem s dscussed from a dfferent vewpont n SPH. Gngold & onaghan [] and Lucy [3] have developed Smooth Partcle Hydrodynamcs method (SPH). We have extended SPH theoretcally n the present paper. In Vortex Blob ethod (VB), the abducton s approxmated reasonably. However, the dffuson of vortces can t be approxmated precsely. We have shown n the present paper that ths problem can be solved, f we apply the deas developed by SPH. ovng Partcle Sem-mplct method (PS) uses a smlar dscretzaton of the ntal and boundary value problem as SPH and VB. However, the mathematcal background of PS s not suffcent. We have shown n the present paper that the mathematcal background of the dscrete gradent operator s strengthened by applyng the deas developed by SPH. However, that of the dscrete Laplacan operator could not gven. 443
31 H. Isshk The dscrete gradent and Laplacan operators of SPH nclude the frst, frst and second dervatves of the weght functon wth respect to the space coordnates, respectvely. On the other hand, those of PS nclude only the weght functon tself. Ths may be the bggest dfference between the dscrete dfferental operators n SPH and PH. In the present paper, solutons by FD (Fnte Dfference ethod), SPH and PS are compared numercally n detal. ) Effects of mesh on the dscrete gradent and Laplacan operators were studed. PS showed good results wth respect to the random mesh. ) The FD, SPH and PS were appled to the ntal value problems, and the effects of the dfference of the soluton method were studed. In some cases, the solutons of ntal value problem showed a dfference between SPH and PS. 3) The dscrete Laplacan operator of SPH s senstve to the spacal dscontnutes of the soluton functon. Hence, n the case of the ntal value problem, the dscontnuty n the ntal condton should be smoothened beforehand, and the small amount of the artfcal vscosty should be ntroduced. 4) The author has once studed a mathematcal background of PS theoretcally [] and dd some numercal calculatons. In very lmted cases, the dscrete Laplacan operator of PS can be obtaned mathematcally. However, the generalzaton to the general mesh was not obtaned. Hence, we are oblged to apply the Laplacan operator as a bold approxmaton n case of the general mesh. 5) Recently, Ng, Hwang and Sheu [6] dscussed the accuracy of the dscrete Laplacan operator of PS theoretcally and numercally. They clarfed an mportant aspect of the propertes of the operator. As one of the propertes, they ponted out n n concluson of Ref. [6] that PS gave generally a favorable result n case of the rregular mesh. We also had a smlar mpresson as expressed n () above. The support of Gaussan weght n SPH s nfnte. In the present paper, weghts of a Gaussan-type of fnte support wth C contnuty were also gven. Acknowledgements The author thanks sncerely Prof. G. Yagawa, Toyo Unv., Japan for hs support and encouragement. References [] Isshk, H. (0) Dscrete Dfferental Operators on Irregular Nodes (DDIN). Internatonal Journal for Numercal ethods n Engneerng, 88, [] Gngold, R.A. and onaghan, J.J. (977) Smoothed Partcle Hydrodynamcs: Theory and Applcaton to Non-Sphercal Stars. onthly Notces of the Royal Astronomcal Socety, 8, mp;whole_paper=yes&type=printer&fletype=.pdf [3] Lucy, L.B. (977) A Numercal Approach to the Testng of the Fsson Hypothess. The Astronomcal Journal, 8, [4] Chorn, A.J. (973) Numercal Study of Slghtly Vscous Flow. Journal of Flud echancs, 57, [5] Yokoyama,., Kubota, Y., Kkuch, K., Yagawa, G. and ochzuk, O. (04) Some Remarks on Surface Condtons of Sold Body Plungng nto Water wth Partcle ethod. Advanced odelng and Smulaton n Engneerng Scences,, [6] Ng, K.C., Hwang, Y.H. and Sheu, T.W.H. (04) On the Accuracy Assessment of Laplacan odels n PS. Computer Physcs Communcatons, 85, [7] Isshk, H., Nagata, S. and Ima, Y. (04) Soluton of Vscous Flow around a Crcular Cylnder by a New Integral Representaton ethod (NIR). The Assocaton for Japan Exchange and Teachng,, fle:///c:/users/l/downloads/ pb%0().pdf [8] Koshzuka, S. and Oka, Y. (996) ovng Partcle Sem-Implct ethod for Fragmentaton of Incompressble Flud. Nuclear Scence and Engneerng, 3, [9] oscardn, L. and Dolag, K. (0) Chapter 4: Cosmology wth Numercal Smulatons. In: atarrese, S., Colp,., 444
32 H. Isshk Gorn, V. and oschella, U., Eds., Dark atter and Dark Energy, Sprnger, Berln, [0] onaghan, J.J. (99) Smoothed Partcle Hydrodynamcs. Annual Revew of Astronomy and Astrophyscs, 30, [] Wkpeda, Burgers Equaton
33
Numerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationarxiv: v1 [physics.flu-dyn] 16 Sep 2013
Three-Dmensonal Smoothed Partcle Hydrodynamcs Method for Smulatng Free Surface Flows Rzal Dw Prayogo a,b, Chrstan Fredy Naa a a Faculty of Mathematcs and Natural Scences, Insttut Teknolog Bandung, Jl.
More informationModule 1 : The equation of continuity. Lecture 1: Equation of Continuity
1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More information(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate
Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1
More informationTensor Smooth Length for SPH Modelling of High Speed Impact
Tensor Smooth Length for SPH Modellng of Hgh Speed Impact Roman Cherepanov and Alexander Gerasmov Insttute of Appled mathematcs and mechancs, Tomsk State Unversty 634050, Lenna av. 36, Tomsk, Russa RCherepanov82@gmal.com,Ger@npmm.tsu.ru
More informationSTUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS
Blucher Mechancal Engneerng Proceedngs May 0, vol., num. www.proceedngs.blucher.com.br/evento/0wccm STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS Takahko Kurahash,
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationThe Tangential Force Distribution on Inner Cylinder of Power Law Fluid Flowing in Eccentric Annuli with the Inner Cylinder Reciprocating Axially
Open Journal of Flud Dynamcs, 2015, 5, 183-187 Publshed Onlne June 2015 n ScRes. http://www.scrp.org/journal/ojfd http://dx.do.org/10.4236/ojfd.2015.52020 The Tangental Force Dstrbuton on Inner Cylnder
More informationNUMERICAL MODEL FOR NON-DARCY FLOW THROUGH COARSE POROUS MEDIA USING THE MOVING PARTICLE SIMULATION METHOD
THERMAL SCIENCE: Year 2018, Vol. 22, No. 5, pp. 1955-1962 1955 NUMERICAL MODEL FOR NON-DARCY FLOW THROUGH COARSE POROUS MEDIA USING THE MOVING PARTICLE SIMULATION METHOD Introducton by Tomok IZUMI a* and
More informationThe Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method
Journal of Electromagnetc Analyss and Applcatons, 04, 6, 0-08 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More informationComputational Fluid Dynamics. Smoothed Particle Hydrodynamics. Simulations. Smoothing Kernels and Basis of SPH
Computatonal Flud Dynamcs If you want to learn a bt more of the math behnd flud dynamcs, read my prevous post about the Naver- Stokes equatons and Newtonan fluds. The equatons derved n the post are the
More informationAppendix B. The Finite Difference Scheme
140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton
More informationHigh resolution entropy stable scheme for shallow water equations
Internatonal Symposum on Computers & Informatcs (ISCI 05) Hgh resoluton entropy stable scheme for shallow water equatons Xaohan Cheng,a, Yufeng Ne,b, Department of Appled Mathematcs, Northwestern Polytechncal
More informationON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION
Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationIrregular vibrations in multi-mass discrete-continuous systems torsionally deformed
(2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected
More informationA new Approach for Solving Linear Ordinary Differential Equations
, ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of
More informationThe Finite Element Method
The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of
More informationThe Order Relation and Trace Inequalities for. Hermitian Operators
Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence
More informationTHE STURM-LIOUVILLE EIGENVALUE PROBLEM - A NUMERICAL SOLUTION USING THE CONTROL VOLUME METHOD
Journal of Appled Mathematcs and Computatonal Mechancs 06, 5(), 7-36 www.amcm.pcz.pl p-iss 99-9965 DOI: 0.75/jamcm.06..4 e-iss 353-0588 THE STURM-LIOUVILLE EIGEVALUE PROBLEM - A UMERICAL SOLUTIO USIG THE
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationFixed point method and its improvement for the system of Volterra-Fredholm integral equations of the second kind
MATEMATIKA, 217, Volume 33, Number 2, 191 26 c Penerbt UTM Press. All rghts reserved Fxed pont method and ts mprovement for the system of Volterra-Fredholm ntegral equatons of the second knd 1 Talaat I.
More informationQuantum Particle Motion in Physical Space
Adv. Studes Theor. Phys., Vol. 8, 014, no. 1, 7-34 HIKARI Ltd, www.-hkar.co http://dx.do.org/10.1988/astp.014.311136 Quantu Partcle Moton n Physcal Space A. Yu. Saarn Dept. of Physcs, Saara State Techncal
More informationApplication of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems
Mathematca Aeterna, Vol. 1, 011, no. 06, 405 415 Applcaton of B-Splne to Numercal Soluton of a System of Sngularly Perturbed Problems Yogesh Gupta Department of Mathematcs Unted College of Engneerng &
More informationA Solution of the Harry-Dym Equation Using Lattice-Boltzmannn and a Solitary Wave Methods
Appled Mathematcal Scences, Vol. 11, 2017, no. 52, 2579-2586 HIKARI Ltd, www.m-hkar.com https://do.org/10.12988/ams.2017.79280 A Soluton of the Harry-Dym Equaton Usng Lattce-Boltzmannn and a Soltary Wave
More informationA large scale tsunami run-up simulation and numerical evaluation of fluid force during tsunami by using a particle method
A large scale tsunam run-up smulaton and numercal evaluaton of flud force durng tsunam by usng a partcle method *Mtsuteru Asa 1), Shoch Tanabe 2) and Masaharu Isshk 3) 1), 2) Department of Cvl Engneerng,
More informationChapter 4 The Wave Equation
Chapter 4 The Wave Equaton Another classcal example of a hyperbolc PDE s a wave equaton. The wave equaton s a second-order lnear hyperbolc PDE that descrbes the propagaton of a varety of waves, such as
More informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationOptimal Control of Temperature in Fluid Flow
Kawahara Lab. 5 March. 27 Optmal Control of Temperature n Flud Flow Dasuke YAMAZAKI Department of Cvl Engneerng, Chuo Unversty Kasuga -3-27, Bunkyou-ku, Tokyo 2-855, Japan E-mal : d33422@educ.kc.chuo-u.ac.jp
More information6.3.4 Modified Euler s method of integration
6.3.4 Modfed Euler s method of ntegraton Before dscussng the applcaton of Euler s method for solvng the swng equatons, let us frst revew the basc Euler s method of numercal ntegraton. Let the general from
More informationResearch Article Green s Theorem for Sign Data
Internatonal Scholarly Research Network ISRN Appled Mathematcs Volume 2012, Artcle ID 539359, 10 pages do:10.5402/2012/539359 Research Artcle Green s Theorem for Sgn Data Lous M. Houston The Unversty of
More informationProcedia Computer Science
Avalable onlne at www.scencedrect.com Proceda Proceda Computer Computer Scence Scence 1 (01) 00 (009) 589 597 000 000 Proceda Computer Scence www.elsever.com/locate/proceda Internatonal Conference on Computatonal
More informationNumerical Solutions of a Generalized Nth Order Boundary Value Problems Using Power Series Approximation Method
Appled Mathematcs, 6, 7, 5-4 Publshed Onlne Jul 6 n ScRes. http://www.scrp.org/journal/am http://.do.org/.436/am.6.77 umercal Solutons of a Generalzed th Order Boundar Value Problems Usng Power Seres Approxmaton
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationCONTROLLED FLOW SIMULATION USING SPH METHOD
HERI COADA AIR FORCE ACADEMY ROMAIA ITERATIOAL COFERECE of SCIETIFIC PAPER AFASES 01 Brasov, 4-6 May 01 GEERAL M.R. STEFAIK ARMED FORCES ACADEMY SLOVAK REPUBLIC COTROLLED FLOW SIMULATIO USIG SPH METHOD
More informationA Solution of Porous Media Equation
Internatonal Mathematcal Forum, Vol. 11, 016, no. 15, 71-733 HIKARI Ltd, www.m-hkar.com http://dx.do.org/10.1988/mf.016.6669 A Soluton of Porous Meda Equaton F. Fonseca Unversdad Naconal de Colomba Grupo
More informationConsideration of 2D Unsteady Boundary Layer Over Oscillating Flat Plate
Proceedngs of the th WSEAS Internatonal Conference on Flud Mechancs and Aerodynamcs, Elounda, Greece, August -, (pp-) Consderaton of D Unsteady Boundary Layer Over Oscllatng Flat Plate N.M. NOURI, H.R.
More informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12
REVIEW Lecture 11: 2.29 Numercal Flud Mechancs Fall 2011 Lecture 12 End of (Lnear) Algebrac Systems Gradent Methods Krylov Subspace Methods Precondtonng of Ax=b FINITE DIFFERENCES Classfcaton of Partal
More informationVisco-Rubber Elastic Model for Pressure Sensitive Adhesive
Vsco-Rubber Elastc Model for Pressure Senstve Adhesve Kazuhsa Maeda, Shgenobu Okazawa, Koj Nshgch and Takash Iwamoto Abstract A materal model to descrbe large deformaton of pressure senstve adhesve (PSA
More informationOn an Extension of Stochastic Approximation EM Algorithm for Incomplete Data Problems. Vahid Tadayon 1
On an Extenson of Stochastc Approxmaton EM Algorthm for Incomplete Data Problems Vahd Tadayon Abstract: The Stochastc Approxmaton EM (SAEM algorthm, a varant stochastc approxmaton of EM, s a versatle tool
More informationUniqueness of Weak Solutions to the 3D Ginzburg- Landau Model for Superconductivity
Int. Journal of Math. Analyss, Vol. 6, 212, no. 22, 195-114 Unqueness of Weak Solutons to the 3D Gnzburg- Landau Model for Superconductvty Jshan Fan Department of Appled Mathematcs Nanjng Forestry Unversty
More information3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X
Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number
More informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georga Tech PHYS 624 Mathematcal Methods of Physcs I Instructor: Predrag Cvtanovć Fall semester 202 Homework Set #7 due October 30 202 == show all your work for maxmum credt == put labels ttle legends
More informationSignificance of Dirichlet Series Solution for a Boundary Value Problem
IOSR Journal of Engneerng (IOSRJEN) ISSN (e): 5-3 ISSN (p): 78-879 Vol. 6 Issue 6(June. 6) V PP 8-6 www.osrjen.org Sgnfcance of Drchlet Seres Soluton for a Boundary Value Problem Achala L. Nargund* and
More informationUsing T.O.M to Estimate Parameter of distributions that have not Single Exponential Family
IOSR Journal of Mathematcs IOSR-JM) ISSN: 2278-5728. Volume 3, Issue 3 Sep-Oct. 202), PP 44-48 www.osrjournals.org Usng T.O.M to Estmate Parameter of dstrbutons that have not Sngle Exponental Famly Jubran
More informationMore metrics on cartesian products
More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of
More informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
More informationResearch Article Cubic B-Spline Collocation Method for One-Dimensional Heat and Advection-Diffusion Equations
Appled Mathematcs Volume 22, Artcle ID 4587, 8 pages do:.55/22/4587 Research Artcle Cubc B-Splne Collocaton Method for One-Dmensonal Heat and Advecton-Dffuson Equatons Joan Goh, Ahmad Abd. Majd, and Ahmad
More informationConvexity preserving interpolation by splines of arbitrary degree
Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete
More informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More informationIntroduction to Computational Fluid Dynamics
Introducton to Computatonal Flud Dynamcs M. Zanub 1, T. Mahalakshm 2 1 (PG MATHS), Department of Mathematcs, St. Josephs College of Arts and Scence for Women-Hosur, Peryar Unversty 2 Assstance professor,
More information1-Dimensional Advection-Diffusion Finite Difference Model Due to a Flow under Propagating Solitary Wave
014 4th Internatonal Conference on Future nvronment and nergy IPCB vol.61 (014) (014) IACSIT Press, Sngapore I: 10.776/IPCB. 014. V61. 6 1-mensonal Advecton-ffuson Fnte fference Model ue to a Flow under
More informationAERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY
LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY The Bot-Savart Law The velocty nduced by the sngular vortex lne wth the crculaton can be determned by means of the Bot- Savart formula
More informationSystem in Weibull Distribution
Internatonal Matheatcal Foru 4 9 no. 9 94-95 Relablty Equvalence Factors of a Seres-Parallel Syste n Webull Dstrbuton M. A. El-Dacese Matheatcs Departent Faculty of Scence Tanta Unversty Tanta Egypt eldacese@yahoo.co
More informationKinematics of Fluids. Lecture 16. (Refer the text book CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlines) 17/02/2017
17/0/017 Lecture 16 (Refer the text boo CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlnes) Knematcs of Fluds Last class, we started dscussng about the nematcs of fluds. Recall the Lagrangan and Euleran
More informationSolution of the Navier-Stokes Equations
Numercal Flud Mechancs Fall 2011 Lecture 25 REVIEW Lecture 24: Soluton of the Naver-Stokes Equatons Dscretzaton of the convectve and vscous terms Dscretzaton of the pressure term Conservaton prncples Momentum
More informationA NUMERICAL COMPARISON OF LANGRANGE AND KANE S METHODS OF AN ARM SEGMENT
Internatonal Conference Mathematcal and Computatonal ology 0 Internatonal Journal of Modern Physcs: Conference Seres Vol. 9 0 68 75 World Scentfc Publshng Company DOI: 0.4/S009450059 A NUMERICAL COMPARISON
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 6
REVIEW of Lecture 5 2.29 Numercal Flud Mechancs Fall 2011 Lecture 6 Contnuum Hypothess and conservaton laws Macroscopc Propertes Materal covered n class: Dfferental forms of conservaton laws Materal Dervatve
More informationPART 8. Partial Differential Equations PDEs
he Islamc Unverst of Gaza Facult of Engneerng Cvl Engneerng Department Numercal Analss ECIV 3306 PAR 8 Partal Dfferental Equatons PDEs Chapter 9; Fnte Dfference: Ellptc Equatons Assocate Prof. Mazen Abualtaef
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationNew Method for Solving Poisson Equation. on Irregular Domains
Appled Mathematcal Scences Vol. 6 01 no. 8 369 380 New Method for Solvng Posson Equaton on Irregular Domans J. Izadan and N. Karamooz Department of Mathematcs Facult of Scences Mashhad BranchIslamc Azad
More informationDETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM
Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI
More informationComparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method
Appled Mathematcal Scences, Vol. 7, 0, no. 47, 07-0 HIARI Ltd, www.m-hkar.com Comparson of the Populaton Varance Estmators of -Parameter Exponental Dstrbuton Based on Multple Crtera Decson Makng Method
More informationGeneralized Linear Methods
Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set
More informationSolving Fractional Nonlinear Fredholm Integro-differential Equations via Hybrid of Rationalized Haar Functions
ISSN 746-7659 England UK Journal of Informaton and Computng Scence Vol. 9 No. 3 4 pp. 69-8 Solvng Fractonal Nonlnear Fredholm Integro-dfferental Equatons va Hybrd of Ratonalzed Haar Functons Yadollah Ordokhan
More informationAsymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation
Nonl. Analyss and Dfferental Equatons, ol., 4, no., 5 - HIKARI Ltd, www.m-har.com http://dx.do.org/.988/nade.4.456 Asymptotcs of the Soluton of a Boundary alue Problem for One-Characterstc Dfferental Equaton
More informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More informationStatistical Inference. 2.3 Summary Statistics Measures of Center and Spread. parameters ( population characteristics )
Ismor Fscher, 8//008 Stat 54 / -8.3 Summary Statstcs Measures of Center and Spread Dstrbuton of dscrete contnuous POPULATION Random Varable, numercal True center =??? True spread =???? parameters ( populaton
More informationMATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018
MATH 5630: Dscrete Tme-Space Model Hung Phan, UMass Lowell March, 08 Newton s Law of Coolng Consder the coolng of a well strred coffee so that the temperature does not depend on space Newton s law of collng
More informationResearch Article Relative Smooth Topological Spaces
Advances n Fuzzy Systems Volume 2009, Artcle ID 172917, 5 pages do:10.1155/2009/172917 Research Artcle Relatve Smooth Topologcal Spaces B. Ghazanfar Department of Mathematcs, Faculty of Scence, Lorestan
More informationSimulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests
Smulated of the Cramér-von Mses Goodness-of-Ft Tests Steele, M., Chaselng, J. and 3 Hurst, C. School of Mathematcal and Physcal Scences, James Cook Unversty, Australan School of Envronmental Studes, Grffth
More informationPerfect Fluid Cosmological Model in the Frame Work Lyra s Manifold
Prespacetme Journal December 06 Volume 7 Issue 6 pp. 095-099 Pund, A. M. & Avachar, G.., Perfect Flud Cosmologcal Model n the Frame Work Lyra s Manfold Perfect Flud Cosmologcal Model n the Frame Work Lyra
More informationResearch Article A Multilevel Finite Difference Scheme for One-Dimensional Burgers Equation Derived from the Lattice Boltzmann Method
Appled Mathematcs Volume 01, Artcle ID 9590, 13 pages do:10.1155/01/9590 Research Artcle A Multlevel Fnte Dfference Scheme for One-Dmensonal Burgers Equaton Derved from the Lattce Boltzmann Method Qaoe
More informationAdditional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty
Addtonal Codes usng Fnte Dfference Method Benamn Moll 1 HJB Equaton for Consumpton-Savng Problem Wthout Uncertanty Before consderng the case wth stochastc ncome n http://www.prnceton.edu/~moll/ HACTproect/HACT_Numercal_Appendx.pdf,
More informationSIMULATION OF SOUND WAVE PROPAGATION IN TURBULENT FLOWS USING A LATTICE-BOLTZMANN SCHEME. Abstract
SIMULATION OF SOUND WAVE PROPAGATION IN TURBULENT FLOWS USING A LATTICE-BOLTZMANN SCHEME PACS REFERENCE: 43.20.Mv Andreas Wlde Fraunhofer Insttut für Integrerte Schaltungen, Außenstelle EAS Zeunerstr.
More informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
More informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
More informationA new integrated-rbf-based domain-embedding scheme for solving fluid-flow problems
Home Search Collectons Journals About Contact us My IOPscence A new ntegrated-rbf-based doman-embeddng scheme for solvng flud-flow problems Ths artcle has been downloaded from IOPscence. Please scroll
More informationNumerical Solution of Ordinary Differential Equations
Numercal Methods (CENG 00) CHAPTER-VI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples
More informationEcon107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)
I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes
More informationOpen Systems: Chemical Potential and Partial Molar Quantities Chemical Potential
Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,
More informationThe equation of motion of a dynamical system is given by a set of differential equations. That is (1)
Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence
More informationComparative Studies of Law of Conservation of Energy. and Law Clusters of Conservation of Generalized Energy
Comparatve Studes of Law of Conservaton of Energy and Law Clusters of Conservaton of Generalzed Energy No.3 of Comparatve Physcs Seres Papers Fu Yuhua (CNOOC Research Insttute, E-mal:fuyh1945@sna.com)
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationJournal of Fluid Science and Technology
Journal of Flud Scence and Technology Numercal Smulaton of Incompressble Flows around a Fsh Model at Low Reynolds Number Usng Seamless Vrtual Boundary Method * Hdetosh NISHIDA ** and Kyohe TAJIRI ** **Department
More informationEEE 241: Linear Systems
EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they
More informationMathematical Preparations
1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the
More informationCubic Trigonometric B-Spline Applied to Linear Two-Point Boundary Value Problems of Order Two
World Academy of Scence Engneerng and echnology Internatonal Journal of Mathematcal and omputatonal Scences Vol: No:0 00 ubc rgonometrc B-Splne Appled to Lnear wo-pont Boundary Value Problems of Order
More informationAn Application of Fuzzy Hypotheses Testing in Radar Detection
Proceedngs of the th WSES Internatonal Conference on FUZZY SYSEMS n pplcaton of Fuy Hypotheses estng n Radar Detecton.K.ELSHERIF, F.M.BBDY, G.M.BDELHMID Department of Mathematcs Mltary echncal Collage
More informationPublication 2006/01. Transport Equations in Incompressible. Lars Davidson
Publcaton 2006/01 Transport Equatons n Incompressble URANS and LES Lars Davdson Dvson of Flud Dynamcs Department of Appled Mechancs Chalmers Unversty of Technology Göteborg, Sweden, May 2006 Transport
More informationHomework Assignment 3 Due in class, Thursday October 15
Homework Assgnment 3 Due n class, Thursday October 15 SDS 383C Statstcal Modelng I 1 Rdge regresson and Lasso 1. Get the Prostrate cancer data from http://statweb.stanford.edu/~tbs/elemstatlearn/ datasets/prostate.data.
More informationBasic concept of reactive flows. Basic concept of reactive flows Combustion Mixing and reaction in high viscous fluid Application of Chaos
Introducton to Toshhsa Ueda School of Scence for Open and Envronmental Systems Keo Unversty, Japan Combuston Mxng and reacton n hgh vscous flud Applcaton of Chaos Keo Unversty 1 Keo Unversty 2 What s reactve
More informationACTM State Calculus Competition Saturday April 30, 2011
ACTM State Calculus Competton Saturday Aprl 30, 2011 ACTM State Calculus Competton Sprng 2011 Page 1 Instructons: For questons 1 through 25, mark the best answer choce on the answer sheet provde Afterward
More informationWeek 9 Chapter 10 Section 1-5
Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,
More information