ICES REPORT Discontinuous and enriched Galerkin methods for phase-field fracture propagation in elasticity
|
|
- Marion Hines
- 5 years ago
- Views:
Transcription
1 ICES REPORT Fbruary 2016 Discontinuous and nrichd Galrkin mthods for phas-fild fractur propagation in lasticity by Prashant Mital, Thomas Wick, Mary F. Whlr, Grgina Pnchva Th Institut for Computational Enginring and Scincs Th Univrsity of Txas at Austin Austin, Txas Rfrnc: Prashant Mital, Thomas Wick, Mary F. Whlr, Grgina Pnchva, "Discontinuous and nrichd Galrkin mthods for phas-fild fractur propagation in lasticity," ICES REPORT 16-05, Th Institut for Computational Enginring and Scincs, Th Univrsity of Txas at Austin, Fbruary 2016.
2 Accptd for publication in Numrical Mathmatics and Advancd Applications, ENUMATH-2015 Procdings, Springr, 2016 Discontinuous and nrichd Galrkin mthods for phas-fild fractur propagation in lasticity Prashant Mital Thomas Wick Mary F. Whlr Grgina Pnchva In this work, w introduc discontinuous Galrkin and nrichd Galrkin formulations for th spatial discrtization of phas-fild fractur propagation. Th nonlinar coupld systm is formulatd in trms of th Eulr-Lagrang quations, which ar subjct to a crack irrvrsibility condition. Th rsulting variational inquality is solvd in a quasi-monolithic way in which th irrvrsibility condition is incorporatd with th hlp of an augmntd Lagrangian tchniqu. Th rlaxd nonlinar systm is tratd with Nwton s mthod. Numrical rsults complt th prsnt study. Kywords: Phas fild fractur; Discontinuous Galrkin; SIPG; NIPG; IIPG; augmntd Lagrangian mthod 1 Introduction Fractur propagation in lasticity, plasticity, and porous mdia is currntly on of th major rsarch topics in mchanical, nrgy, and nvironmntal nginring. In this papr, w concntrat spcifically on fractur propagation in lasticity. W considr a variational approach for brittl fractur introducd in [6], which has bn latr formulatd in trms of a thrmodynamically-consistnt phas fild tchniqu [8]. In fact, variational and phas fild formulations for fractur ar activ rsarch aras as attstd in rcnt yars,.g., [3, 9, 2, 4, 1, 10]. Our motivations for mploying a phas fild modl ar that fractur nuclation, propagation, kinking, and curvilinar paths ar automatically includd in th modl; post-procssing of strss intnsity factors and rmshing rsolving th crack path ar avoidd. Furthrmor, th undrlying quations ar basd on continuum mchanics principls that can b tratd with adaptiv Galrkin finit lmnts. In this work, w xtnd xisting Galrkin formulations for phas-fild fractur with rgard to two major aspcts: Spatial discrtization of th displacmnt fild with discontinuous Galrkin (DG) finit lmnts rsulting in NIPG [12] and IIPG mthods [5] and an nrichd Galrkin (EG) formulation [13]; Formulation of a quasi-monolithic augmntd Lagrangian itration for th nonlinar coupld displacmntphas-fild systm. Ths framworks ar formulatd in Sctions 2-4 and ar substantiatd with numrical tsts in Sction 5. Cntr for Subsurfac Modling, Th Institut for Computational Enginring and Scincs, Th Univrsity of Txas at Austin, Austin, Txas 78712, USA RICAM, Austrian Acadmy of Scincs Altnbrgr Str Linz, Austria Cntr for Subsurfac Modling, Th Institut for Computational Enginring and Scincs, Th Univrsity of Txas at Austin, Austin, Txas 78712, USA Cntr for Subsurfac Modling, Th Institut for Computational Enginring and Scincs, Th Univrsity of Txas at Austin, Austin, Txas 78712, USA 1
3 Mital, Wick, Whlr, Pnchva 2 Th phas-fild fractur modl W limit our attntion to 2-dimnsional problms and lt Ω R 2, b a smooth, opn, connctd and boundd st. W dnot th L 2 scalar product with (, ), and assum that th crack C is a 1-dimnsional st, not ncssarily connctd, containd in Ω. Using th variational/phas-fild approach to fractur [6, 3], th crack C is rprsntd using a continuous phas fild variabl ϕ : Ω [0, 1]. This valu of th phas fild variabl intrpolats btwn th brokn (ϕ = 0) and unbrokn (ϕ = 1) stats of th matrial. Th diffusiv transition zon btwn ths two stats is controlld by a rgularization paramtr ε > 0. Imposing a crack irrvrsibility condition ϕ ϕ n 1 (whr ϕ n 1 := ϕ(t n 1 ) dnots th prvious tim stp solution), and furthr ingrdints for a thrmodynamically consistnt phas-fild framwork [8] rsult in th following Eulr-Lagrang formulation: Formulation 1: Find vctor-valud displacmnts and a scalar-valud phas-fild variabl (u, ϕ) {ū+v } W such that ( ((1 κ)ϕ 2 + κ ) σ(u), (w)) = 0 w V, (1) as wll as, (1 κ)(ϕσ(u) : (u), ψ ϕ) + G c ( 1 ε (1 ϕ, ψ ϕ) + ε( ϕ, ψ ϕ) ) 0 ψ W in L (Ω), (2) whr V := H 1 0 (Ω), W in := {w H 1 (Ω) w ϕ n 1 1 a.. on Ω} and W := H 1 (Ω). Furthrmor, σ = σ(u) = 2µ s +λ s tr()i is th strss tnsor with µ s, λ s > 0, and (u) = 0.5( u+ u T ) is th linarizd strain tnsor. Th critical nrgy rlas rat is G c > 0. Th domain is subjct to boundary conditions, and w assum Γ D, with th possibly non-homognous and tim-dpndnt Dirichlt boundary conditions ū. Morovr, κ is a rgularization paramtr for th lastic nrgy boundd blow by 0, such that κ ε, s.g., [3]. To trat crack irrvrsibility, w us th augmntd-lagrangian formulation dscribd in [14]. To apply this mthod, w bgin by approximating th tim drivativ t ϕ using th backward diffrnc t ϕ t ϕ = ϕ ϕn 1 t 1 t ((λ + γ(ϕ ϕ n 1 )) + ), t = t n t n 1. Hr, λ and γ ar a pnalization function and paramtr, rspctivly, and ϕ n 1 is th phas fild solution at th prvious tim stp. Morovr, (x) + := max{0, x}. 3 A quasi-monolithic incrmntal formulation W choos a quasi-monolithic approach [7] as this rducs algorithmic complxity and has bn dmonstratd to b numrically robust and fficint whn ϕ is rplacd by th xtrapolation ϕ in th first trm of Formulation 2. Th rason for choosing ϕ is th nd to circumvnt th non-convxity of th undrlying nrgy functional. 2
4 DG and EG for phas-fild fractur Algorithm 1 Solution algorithm For ach tim t n : Lt m = 0; choos initial λ m L 2 (Ω), γ > 0. rpat Lt k = 0; choos initial Ũk V h W h. rpat Find δu k solving A (U k )(δu k, Ψ) = A(U k )(Ψ) Updat Ũk+1 Ũk + δu k Updat k k + 1 until Stopping critrion U k U k 1 TOL 2 is satisfid. St U m+1 = (u m+1, ϕ m+1 ) = Ũk Updat λ m+1 = min(0, λ m + γϕ m+1 ) + (λ m + γ(ϕ m+1 ϕ n 1 )) + Updat m m + 1 until Stopping critrion max( u m 1 u m, λ m 1 λ m ) TOL i is satisfid. (u n, ϕ n ) = (u m+1, ϕ m+1 ) Formulation 2: Find U := {u, ϕ} {ū + V } W, whr V := H 1 0 (Ω) and W := H1 (Ω), such that A(U)(Ψ) = 0 Ψ := {w, ψ} V W, (3) whr A( )( ) is th following smi-linar form ( ((1 A(U)(Ψ) = κ) ϕ 2 + κ ) ) σ(u), (w) + (1 κ)(ϕσ(u) : (u), ψ) ( + G c 1 ) ε (1 ϕ, ψ) + ε( ϕ, ψ) + 1 ( ) (λ + γ(ϕ ϕ n 1 )) +, ψ. (4) t Hr ϕ is a linar xtrapolation of tim-laggd ϕ, i.. ϕ ϕ := ϕ(ϕ n 1, ϕ n 2 ), with ϕ n 1, ϕ n 2 dnoting solutions to prvious tim stps. Solving th nonlinar variational problm (4) is prformd with Nwton s mthod and lin sarch backtracking. Th rsulting solution algorithm is outlind in Algorithm 1. 4 Spatial discrtization with DG and EG In this sction w stablish ky notations for DG and EG followd by th mathmatical statmnt of th discrt variational forms. On a conforming subdivision E h of a polygonal domain Ω subdividd into lmnts E w dfin th discontinuous finit lmnt subspac to b D k (E h ) = {v L 2 (Ω) : E E h, v E P k (E)}, (5) whr P k (E) dnots th spac of picwis polynomials of total dgr lss than or qual to k on E. W also dfin th spac of CG approximating polynomials nrichd with discontinuous picwis constants Hr D C k (E h) is th CG approximating spac dfind as D C0 k (E h) := D C k (E h) D 0 (E h ). (6) D C k (E h) = {v C(Ω) : E E h, v E P C k (E), v Γ D = 0}, (7) 3
5 Mital, Wick, Whlr, Pnchva whr P C k (E) dnots th spac of continuous picwis polynomials of total dgr lss than or qual to k on E. (a) D 1 C (E h ) (b) D 1 (E h ) (c) D C0 1 (E h ) D 2 C (E h ) D 2 (E h ) D C0 2 (E h ) Figur 1: Support points for bilinar and biquadratic basis functions. (a) CG: all support points in rd. (b) DG: support points for lft lmnt in rd, for right lmnt in blu. Th common dg has two sts of support points - on from ach lmnt. (c) EG: support points from CG approximating spac in rd, picwis constants in blu. Only th picwis constant dgr of frdom is discontinuous across th common dg. In ordr to dscrib th vctor-valud displacmnts, w considr th spacs of vctor functions that gnraliz (5) and (6): D k (E h ) = (D k (E h )) d, D C0 k (E h) = (Dk C0(E h)) d, whr d is th numbr of spatial dimnsions. W not that th functions in D k (E h ) and D C0 k (E h) ar discontinuous along th dgs (or facs) of th msh. Now, considr two nighboring lmnts E1 and E 2 that shar a common sid. Naturally thn, thr ar two tracs of w D k (E h ) along. W considr n to b th normal vctor associatd with to b orintd from E1 to E2 and dfin: {w} = 1 2 (w E1 ) (w E2 ), [w] = (w E1 ) (w E2 ) = E 1 E 2. W xtnd this dfinition to lmnts on th boundary Ω as: {w} = [w] = (w E 1 ) = E1 Ω. Furthr, w dnot by th lngth of an dg in d = 2. W now stat th quations corrsponding to a discontinuous spatial discrtization dirctly from inspction of th monolithic formulation (4) and th DG-schm for pur linar lasticity,.g., [11]. W pursu a discontinuous rprsntation of th displacmnt variabl u only, rcognizing that th rgularization in th cas of th phas fild variabl ϕ nforcs its continuity across th crack. W augmnt (4) with th jump and pnalization trms to dfin th discrt incrmntal smi-linar form A(U h )(Ψ h ) = ( (1 κ) ϕ 2 + κ ) ) σ(u), (w) Γ h Γ D +η Γ D E E h E { ((1 κ) ϕ 2 + κ ) } { σ(u) n [w] +η ((1 κ) ϕ 2 +κ ) } σ(w) n [u] Γ h (((1 κ) ϕ 2 + κ ) σ(w) n )(u g D ) + Γ h δ β [u][w] + +(1 κ)(ϕσ(u) : (u), ψ) +G c ( 1 ε (1 ϕ, ψ) +ε( ϕ, ψ) ) + 1 t Γ D δ β (u g D )w ) ( (λ + γ(ϕ ϕ n 1 )) +, ψ whr η = 1 (NIPG) or η = 0 (IIPG); and δ > 0 and β > 0 (hr β = 1) ar th DG-pnalization and, (8) 4
6 DG and EG for phas-fild fractur suprpnalization paramtrs, rspctivly. Th DG-CG variational problm rads: Find U h := {u, ϕ} {ū h + Vh DG } Wh CG such that A(U h )(Ψ h ) = 0 Ψ h := {w, ψ} V DG h W CG h. (9) Th EG-CG variational problm rads: Find U h := {u, ϕ} {ū h + Vh EG } Wh CG such that A(U h )(Ψ h ) = 0 Ψ h := {w, ψ} V EG h W CG h. (10) Th tst and trial spacs ar Vh DG := [D k (E h )] 2, Vh EG := [Dk C0(E h)] 2, Wh CG := Dk C(E h). Our formulation nforcs both Dirichlt and Numann boundary conditions wakly. Th us of homognous Numann boundary conditions for u and ϕ rsults in a formulation xclusivly dpndnt on Γ D. Th dirctional drivativ of (8) ndd for th Nwton itrations is computd analytically. 5 A numrical tst: singl dg notchd tnsion Th singl dg notchd tnsion tst is a widly usd xprimntal mthodology usd to charactriz th fractur toughnss of various matrials in plan-strain. ū = (0, ū y(t)) Figur 2: Schmatic of th singl-dg-notchd tnsion tst (lft), th final phas-fild crack pattrn at T = s (middl), and comparison of load-displacmnt curvs from our monolithic schm with CG, DG-IIPG and EG-IIPG against rsults rportd by Mih t al. [9] and Histr t al. [7]. W considr a squar plat with a horizontal notch placd at half-hight, running from th right outr surfac to th cntr of th spcimn. Th plat is subjct to zro displacmnt boundary conditions on th bottom surfac, and tim-dpndnt displacmnt on th top surfac. Th lft and right surfacs ar considrd to b tractionfr. Th problm stup is shown in Figur 2. Th matrial paramtrs ar chosn as λ = kN/mm 2, µ = kN/mm 2 and G c = kn/mm. Th displacmnt boundary condition on th top surfac is takn to b ū y (t) = tᾱ with ᾱ = 1 mm/s. Th xpctd rspons of this tst is th build-up of th strss concntration in th vicinity of th crack-tip, followd by unstabl, catastrophic crack growth. Our first objctiv is to study h-convrgnc for fixd ε. W choos t = s for th first 50 loading stps, aftr which t = s. This adaptivity in th tim stp is ncssary to captur th rapid movmnt of th crack tip. W choos ε = mm, κ = , and run our cod for h 1 = mm, h 2 = mm, and h 3 = mm. W valuat th surfac load vctor on th top surfac of Ω as τ = (F x, F y ) = Ω top σ(u)nds. In this xampl, w ar particularly intrstd in F y. Our findings for th surfac load volution with varying h ar shown in Figur 3 for th IIPG flavors of EG and DG. It is obsrvd 5
7 Mital, Wick, Whlr, Pnchva that our approach is stabl with spatial msh rfinmnt, and that our solution convrgs as w us finr mshs. Comparison to litratur valus ar displayd in Figur 2 at right. Figur 3: W fix κ and ε on th coarsst msh lvl and vary h. Lft: CG. Middl: DG-IIPG. Right: EG-IIPG. Spatial convrgnc is obsrvd for all schms. Rsults obtaind from DG-NIPG and EG-NIPG ar vry similar and ar thrfor not prsntd hr. Th SIPG mthod (η = 1) yilds unsatisfactory findings, which ar not shown in this work. With th rsults of our schm duly validatd, w procd to study th rlativ fficincy of th schms by comparing th numbr of Nwton itrations takn by ach of thm to convrg. W first invstigat th variation in th numbr of Nwton stps takn with th pnalization paramtr δ. Not that whn w multiply Equations (8) throughout by t, our ffctiv pnalization of th jump bcoms δ t. This is an important dtail that cannot b ovrlookd whil using DG/EG for th phas fild quations bcaus for instanc, using δ = 10 5 with t = s givs an ffctiv pnalization of δ t = 1 which is not sufficintly larg and producs spurious rsults. In th cas of an adaptiv tim stp siz (w usually tak t = s for th first 50 stps, and a smallr tim stp thraftr), th product δ t is rportd for th smallr tim stp. W vary th valus of th ffctiv pnalization and plot th cumulativ numbr of Nwton stps as a function of tim for h = mm, ε = 2h[mm], and κ = Th rsults of this study with IIPG and NIPG ar shown in Figur 4. In Figur 5, w obsrv that DG and EG schms tak much fwr Nwton itrations to convrg than CG spcially aftr th onst of crack growth (approximatly t = s). Figur 4: Nwton convrgnc prformanc with h = mm, t = {10 4, t < 0.005; 10 5, t 0.005}, ε = 2h[mm], and κ = and varying pnalty δ. Lft: DG-NIPG. Right: EG-NIPG. Lft: DG-IIPG. Right: EG-IIPG. Convrgnc is fastr for highr valus of pnalization. For a bttr comparison of th fficincy, w run a tst with th sam physical paramtrs as abov, but with 6
8 DG and EG for phas-fild fractur Figur 5: Singl dg notchd tnsion tst rsults using CG, DG-IIPG and EG-IIPG. Lft: load vs. displacmnt curv. Right: Nwton convrgnc prformanc for constant t = 10 5 s. a uniform tim stp of t = 10 5 s throughout. Th motivation is to supprss th ffct of adaptiv tim stpping on th Nwton prformanc and to giv an unbiasd comparison. Sinc th computational burdn with such a simulation is significant, w only considr th IIPG cas with δ t = Ths rsults ar shown in Figur 5. As w can s, DG and EG tak roughly th sam numbr of Nwton itrations (1800) whil CG taks significantly mor (2520). W also obsrv that th load-displacmnt curvs for all thr mthods ar in rasonabl agrmnt. Hnc, w can conclusivly stat that th Nwton mthod convrgs in fwr itrations for th DG and EG schms than for th CG schm. Furthrmor by inspcting Figur 1, w s that EG has significantly fwr dgrs of frdom than DG. Rfrncs [1] M. AMBATI, T. GERASIMOV, AND L. DE LORENZIS, A rviw on phas-fild modls of brittl fractur and a nw fast hybrid formulation, Computational Mchanics 55:2 (2015), [2] M. J. BORDEN, C. V. VERHOOSEL, M. A. SCOTT, T. J. HUGHES, AND C. M. LANDIS, A phas-fild dscription of dynamic brittl fractur, Computr Mthods in Applid Mchanics and Enginring 217 (2012), [3] B. BOURDIN, G. A. FRANCFORT, AND J. J. MARIGO, Th variational approach to fractur, Journal of lasticity 91:1-3 (2008), [4] S. BURKE, C. ORTNER, AND E. SÜLI, An adaptiv finit lmnt approximation of a variational modl of brittl fractur, SIAM Journal on Numrical Analysis 48:3 (2010), [5] C. DAWSON, S. SUN, AND M. F. WHEELER, Compatibl algorithms for coupld flow and transport, Computr Mthods in Applid Mchanics and Enginring 193:23 (2004), [6] G. A. FRANCFORT AND J.-J. MARIGO, Rvisiting brittl fractur as an nrgy minimization problm, Journal of th Mchanics and Physics of Solids 46:8 (1998), [7] T. HEISTER, M. F. WHEELER, AND T. WICK, A primal-dual activ st mthod and prdictor-corrctor msh adaptivity for computing fractur propagation using a phas-fild approach, Computr Mthods in Applid Mchanics and Enginring 290 (2015),
9 Mital, Wick, Whlr, Pnchva [8] C. MIEHE, F. WELSCHINGER, AND M. HOFACKER, Thrmodynamically consistnt phas-fild modls of fractur: Variational principls and multi-fild f implmntations, Intrnational Journal for Numrical Mthods in Enginring 83:10 (2010), [9] C. MIEHE, M. HOFACKER, AND F. WELSCHINGER, A phas fild modl for rat-indpndnt crack propagation: Robust algorithmic implmntation basd on oprator splits, Computr Mthods in Applid Mchanics and Enginring 199:45 (2010), [10] A. MIKELIĆ, M. F. WHEELER, AND T. WICK, A quasi-static phas-fild approach to prssurizd fracturs, Nonlinarity 28:5 (2015), [11] B. RIVIÈRE, Discontinuous galrkin mthods for solving lliptic and parabolic quations: thory and implmntation, SIAM, [12] B. RIVIÈRE, M. F. WHEELER, K. BANAŚ, ET AL., Discontinuous galrkin mthod applid to a singl phas flow in porous mdia, Computational Goscincs 4:4 (2000), [13] S. SUN AND J. LIU, A locally consrvativ finit lmnt mthod basd on picwis constant nrichmnt of th continuous galrkin mthod, SIAM Journal on Scintific Computing 31:4 (2009), [14] M. F. WHEELER, T. WICK, AND W. WOLLNER, An augmntd-lagrangian mthod for th phasfild approach for prssurizd fracturs, Computr Mthods in Applid Mchanics and Enginring 271 (2014),
Symmetric Interior Penalty Galerkin Method for Elliptic Problems
Symmtric Intrior Pnalty Galrkin Mthod for Elliptic Problms Ykatrina Epshtyn and Béatric Rivièr Dpartmnt of Mathmatics, Univrsity of Pittsburgh, 3 Thackray, Pittsburgh, PA 56, U.S.A. Abstract This papr
More informationAS 5850 Finite Element Analysis
AS 5850 Finit Elmnt Analysis Two-Dimnsional Linar Elasticity Instructor Prof. IIT Madras Equations of Plan Elasticity - 1 displacmnt fild strain- displacmnt rlations (infinitsimal strain) in matrix form
More informationDynamic Modelling of Hoisting Steel Wire Rope. Da-zhi CAO, Wen-zheng DU, Bao-zhu MA *
17 nd Intrnational Confrnc on Mchanical Control and Automation (ICMCA 17) ISBN: 978-1-6595-46-8 Dynamic Modlling of Hoisting Stl Wir Rop Da-zhi CAO, Wn-zhng DU, Bao-zhu MA * and Su-bing LIU Xi an High
More information2008 AP Calculus BC Multiple Choice Exam
008 AP Multipl Choic Eam Nam 008 AP Calculus BC Multipl Choic Eam Sction No Calculator Activ AP Calculus 008 BC Multipl Choic. At tim t 0, a particl moving in th -plan is th acclration vctor of th particl
More informationEinstein Equations for Tetrad Fields
Apiron, Vol 13, No, Octobr 006 6 Einstin Equations for Ttrad Filds Ali Rıza ŞAHİN, R T L Istanbul (Turky) Evry mtric tnsor can b xprssd by th innr product of ttrad filds W prov that Einstin quations for
More informationA Weakly Over-Penalized Non-Symmetric Interior Penalty Method
Europan Socity of Computational Mthods in Scincs and Enginring (ESCMSE Journal of Numrical Analysis, Industrial and Applid Mathmatics (JNAIAM vol.?, no.?, 2007, pp.?-? ISSN 1790 8140 A Wakly Ovr-Pnalizd
More informationHomotopy perturbation technique
Comput. Mthods Appl. Mch. Engrg. 178 (1999) 257±262 www.lsvir.com/locat/cma Homotopy prturbation tchniqu Ji-Huan H 1 Shanghai Univrsity, Shanghai Institut of Applid Mathmatics and Mchanics, Shanghai 272,
More informationLarge Scale Topology Optimization Using Preconditioned Krylov Subspace Recycling and Continuous Approximation of Material Distribution
Larg Scal Topology Optimization Using Prconditiond Krylov Subspac Rcycling and Continuous Approximation of Matrial Distribution Eric d Sturlr*, Chau L**, Shun Wang***, Glaucio Paulino** * Dpartmnt of Mathmatics,
More informationDivision of Mechanics Lund University MULTIBODY DYNAMICS. Examination Name (write in block letters):.
Division of Mchanics Lund Univrsity MULTIBODY DYNMICS Examination 7033 Nam (writ in block lttrs):. Id.-numbr: Writtn xamination with fiv tasks. Plas chck that all tasks ar includd. clan copy of th solutions
More information(Upside-Down o Direct Rotation) β - Numbers
Amrican Journal of Mathmatics and Statistics 014, 4(): 58-64 DOI: 10593/jajms0140400 (Upsid-Down o Dirct Rotation) β - Numbrs Ammar Sddiq Mahmood 1, Shukriyah Sabir Ali,* 1 Dpartmnt of Mathmatics, Collg
More informationDifference -Analytical Method of The One-Dimensional Convection-Diffusion Equation
Diffrnc -Analytical Mthod of Th On-Dimnsional Convction-Diffusion Equation Dalabav Umurdin Dpartmnt mathmatic modlling, Univrsity of orld Economy and Diplomacy, Uzbistan Abstract. An analytical diffrncing
More informationA LOCALLY CONSERVATIVE LDG METHOD FOR THE INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
A LOCALLY CONSERVATIVE LDG METHOD FOR THE INCOMPRESSIBLE NAVIER-STOES EQUATIONS BERNARDO COCBURN, GUIDO ANSCHAT, AND DOMINI SCHÖTZAU Abstract. In this papr, a nw local discontinuous Galrkin mthod for th
More informationMAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design
MAE4700/5700 Finit Elmnt Analysis for Mchanical and Arospac Dsign Cornll Univrsity, Fall 2009 Nicholas Zabaras Matrials Procss Dsign and Control Laboratory Sibly School of Mchanical and Arospac Enginring
More information1 Isoparametric Concept
UNIVERSITY OF CALIFORNIA BERKELEY Dpartmnt of Civil Enginring Spring 06 Structural Enginring, Mchanics and Matrials Profssor: S. Govindj Nots on D isoparamtric lmnts Isoparamtric Concpt Th isoparamtric
More informationA Sub-Optimal Log-Domain Decoding Algorithm for Non-Binary LDPC Codes
Procdings of th 9th WSEAS Intrnational Confrnc on APPLICATIONS of COMPUTER ENGINEERING A Sub-Optimal Log-Domain Dcoding Algorithm for Non-Binary LDPC Cods CHIRAG DADLANI and RANJAN BOSE Dpartmnt of Elctrical
More informationSearch sequence databases 3 10/25/2016
Sarch squnc databass 3 10/25/2016 Etrm valu distribution Ø Suppos X is a random variabl with probability dnsity function p(, w sampl a larg numbr S of indpndnt valus of X from this distribution for an
More informationAnother view for a posteriori error estimates for variational inequalities of the second kind
Accptd by Applid Numrical Mathmatics in July 2013 Anothr viw for a postriori rror stimats for variational inqualitis of th scond kind Fi Wang 1 and Wimin Han 2 Abstract. In this papr, w giv anothr viw
More informationNumerical considerations regarding the simulation of an aircraft in the approaching phase for landing
INCAS BULLETIN, Volum, Numbr 1/ 1 Numrical considrations rgarding th simulation of an aircraft in th approaching phas for landing Ionl Cristinl IORGA ionliorga@yahoo.com Univrsity of Craiova, Alxandru
More informationElements of Statistical Thermodynamics
24 Elmnts of Statistical Thrmodynamics Statistical thrmodynamics is a branch of knowldg that has its own postulats and tchniqus. W do not attmpt to giv hr vn an introduction to th fild. In this chaptr,
More informationSupplementary Materials
6 Supplmntary Matrials APPENDIX A PHYSICAL INTERPRETATION OF FUEL-RATE-SPEED FUNCTION A truck running on a road with grad/slop θ positiv if moving up and ngativ if moving down facs thr rsistancs: arodynamic
More informationAnalysis of a Discontinuous Finite Element Method for the Coupled Stokes and Darcy Problems
Journal of Scintific Computing, Volums and 3, Jun 005 ( 005) DOI: 0.007/s095-004-447-3 Analysis of a Discontinuous Finit Elmnt Mthod for th Coupld Stoks and Darcy Problms Béatric Rivièr Rcivd July 8, 003;
More informationReparameterization and Adaptive Quadrature for the Isogeometric Discontinuous Galerkin Method
Rparamtrization and Adaptiv Quadratur for th Isogomtric Discontinuous Galrkin Mthod Agns Silr, Brt Jüttlr 2 Doctoral Program Computational Mathmatics 2 Institut of Applid Gomtry Johanns Kplr Univrsity
More informationDiscontinuous Galerkin and Mimetic Finite Difference Methods for Coupled Stokes-Darcy Flows on Polygonal and Polyhedral Grids
Discontinuous Galrkin and Mimtic Finit Diffrnc Mthods for Coupld Stoks-Darcy Flows on Polygonal and Polyhdral Grids Konstantin Lipnikov Danail Vassilv Ivan Yotov Abstract W study locally mass consrvativ
More informationA Propagating Wave Packet Group Velocity Dispersion
Lctur 8 Phys 375 A Propagating Wav Packt Group Vlocity Disprsion Ovrviw and Motivation: In th last lctur w lookd at a localizd solution t) to th 1D fr-particl Schrödingr quation (SE) that corrsponds to
More informationIntroduction to Arithmetic Geometry Fall 2013 Lecture #20 11/14/2013
18.782 Introduction to Arithmtic Gomtry Fall 2013 Lctur #20 11/14/2013 20.1 Dgr thorm for morphisms of curvs Lt us rstat th thorm givn at th nd of th last lctur, which w will now prov. Thorm 20.1. Lt φ:
More informationThat is, we start with a general matrix: And end with a simpler matrix:
DIAGON ALIZATION OF THE STR ESS TEN SOR INTRO DUCTIO N By th us of Cauchy s thorm w ar abl to rduc th numbr of strss componnts in th strss tnsor to only nin valus. An additional simplification of th strss
More informationSymmetric centrosymmetric matrix vector multiplication
Linar Algbra and its Applications 320 (2000) 193 198 www.lsvir.com/locat/laa Symmtric cntrosymmtric matrix vctor multiplication A. Mlman 1 Dpartmnt of Mathmatics, Univrsity of San Francisco, San Francisco,
More informationME469A Numerical Methods for Fluid Mechanics
ME469A Numrical Mthods for Fluid Mchanics Handout #5 Gianluca Iaccarino Finit Volum Mthods Last tim w introducd th FV mthod as a discrtization tchniqu applid to th intgral form of th govrning quations
More informationThe Matrix Exponential
Th Matrix Exponntial (with xrciss) by Dan Klain Vrsion 28928 Corrctions and commnts ar wlcom Th Matrix Exponntial For ach n n complx matrix A, dfin th xponntial of A to b th matrix () A A k I + A + k!
More informationHigher order derivatives
Robrto s Nots on Diffrntial Calculus Chaptr 4: Basic diffrntiation ruls Sction 7 Highr ordr drivativs What you nd to know alrady: Basic diffrntiation ruls. What you can larn hr: How to rpat th procss of
More informationCramér-Rao Inequality: Let f(x; θ) be a probability density function with continuous parameter
WHEN THE CRAMÉR-RAO INEQUALITY PROVIDES NO INFORMATION STEVEN J. MILLER Abstract. W invstigat a on-paramtr family of probability dnsitis (rlatd to th Parto distribution, which dscribs many natural phnomna)
More informationNEW APPLICATIONS OF THE ABEL-LIOUVILLE FORMULA
NE APPLICATIONS OF THE ABEL-LIOUVILLE FORMULA Mirca I CÎRNU Ph Dp o Mathmatics III Faculty o Applid Scincs Univrsity Polithnica o Bucharst Cirnumirca @yahoocom Abstract In a rcnt papr [] 5 th indinit intgrals
More informationConstruction of asymmetric orthogonal arrays of strength three via a replacement method
isid/ms/26/2 Fbruary, 26 http://www.isid.ac.in/ statmath/indx.php?modul=prprint Construction of asymmtric orthogonal arrays of strngth thr via a rplacmnt mthod Tian-fang Zhang, Qiaoling Dng and Alok Dy
More informationFinite Strain Elastic-Viscoplastic Model
Finit Strain Elastic-Viscoplastic Modl Pinksh Malhotra Mchanics of Solids,Brown Univrsity Introduction Th main goal of th projct is to modl finit strain rat-dpndnt plasticity using a modl compatibl for
More informationBasic Polyhedral theory
Basic Polyhdral thory Th st P = { A b} is calld a polyhdron. Lmma 1. Eithr th systm A = b, b 0, 0 has a solution or thr is a vctorπ such that π A 0, πb < 0 Thr cass, if solution in top row dos not ist
More informationFourier Transforms and the Wave Equation. Key Mathematics: More Fourier transform theory, especially as applied to solving the wave equation.
Lur 7 Fourir Transforms and th Wav Euation Ovrviw and Motivation: W first discuss a fw faturs of th Fourir transform (FT), and thn w solv th initial-valu problm for th wav uation using th Fourir transform
More informationMachine Detector Interface Workshop: ILC-SLAC, January 6-8, 2005.
Intrnational Linar Collidr Machin Dtctor Intrfac Workshop: ILCSLAC, January 68, 2005. Prsntd by Brtt Parkr, BNLSMD Mssag: Tools ar now availabl to optimiz IR layout with compact suprconducting quadrupols
More informationMiddle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis
Middl East Tchnical Univrsity Dpartmnt of Mchanical Enginring ME 43 Introduction to Finit Elmnt Analysis Chaptr 3 Computr Implmntation of D FEM Ths nots ar prpard by Dr. Cünyt Srt http://www.m.mtu.du.tr/popl/cunyt
More informationLINEAR DELAY DIFFERENTIAL EQUATION WITH A POSITIVE AND A NEGATIVE TERM
Elctronic Journal of Diffrntial Equations, Vol. 2003(2003), No. 92, pp. 1 6. ISSN: 1072-6691. URL: http://jd.math.swt.du or http://jd.math.unt.du ftp jd.math.swt.du (login: ftp) LINEAR DELAY DIFFERENTIAL
More informationAn Efficiency Substructure Method for Nonlinear SSI Analysis of Large-scale Concrete Structures in Time Domain on the ANSYS Platform
An Efficincy Substructur Mthod for Nonlinar SSI Analysis of Larg-scal Concrt Structurs in Tim Domain on th ANSYS Platform J. B. Li, X. Q. Yin, G. Lin School of Civil and Hydraulic Enginring, Dalian Univrsity
More informationThe Matrix Exponential
Th Matrix Exponntial (with xrciss) by D. Klain Vrsion 207.0.05 Corrctions and commnts ar wlcom. Th Matrix Exponntial For ach n n complx matrix A, dfin th xponntial of A to b th matrix A A k I + A + k!
More informationUNIFIED A POSTERIORI ERROR ESTIMATOR FOR FINITE ELEMENT METHODS FOR THE STOKES EQUATIONS
UNIFIED A POSTERIORI ERROR ESTIMATOR FOR FINITE ELEMENT METHODS FOR THE STOKES EQUATIONS JUNPING WANG, YANQIU WANG, AND XIU YE Abstract. This papr is concrnd with rsidual typ a postriori rror stimators
More informationVII. Quantum Entanglement
VII. Quantum Entanglmnt Quantum ntanglmnt is a uniqu stat of quantum suprposition. It has bn studid mainly from a scintific intrst as an vidnc of quantum mchanics. Rcntly, it is also bing studid as a basic
More informationMA 262, Spring 2018, Final exam Version 01 (Green)
MA 262, Spring 218, Final xam Vrsion 1 (Grn) INSTRUCTIONS 1. Switch off your phon upon ntring th xam room. 2. Do not opn th xam booklt until you ar instructd to do so. 3. Bfor you opn th booklt, fill in
More informationCOMPUTER GENERATED HOLOGRAMS Optical Sciences 627 W.J. Dallas (Monday, April 04, 2005, 8:35 AM) PART I: CHAPTER TWO COMB MATH.
C:\Dallas\0_Courss\03A_OpSci_67\0 Cgh_Book\0_athmaticalPrliminaris\0_0 Combath.doc of 8 COPUTER GENERATED HOLOGRAS Optical Scincs 67 W.J. Dallas (onday, April 04, 005, 8:35 A) PART I: CHAPTER TWO COB ATH
More informationInference Methods for Stochastic Volatility Models
Intrnational Mathmatical Forum, Vol 8, 03, no 8, 369-375 Infrnc Mthods for Stochastic Volatility Modls Maddalna Cavicchioli Cá Foscari Univrsity of Vnic Advancd School of Economics Cannargio 3, Vnic, Italy
More informationDynamic response of a finite length euler-bernoulli beam on linear and nonlinear viscoelastic foundations to a concentrated moving force
Journal of Mchanical Scinc and Tchnology 2 (1) (21) 1957~1961 www.springrlink.com/contnt/1738-9x DOI 1.17/s1226-1-7-x Dynamic rspons of a finit lngth ulr-brnoulli bam on linar and nonlinar viscolastic
More informationPipe flow friction, small vs. big pipes
Friction actor (t/0 t o pip) Friction small vs larg pips J. Chaurtt May 016 It is an intrsting act that riction is highr in small pips than largr pips or th sam vlocity o low and th sam lngth. Friction
More informationTOPOLOGY DESIGN OF STRUCTURE LOADED BY EARTHQUAKE. Vienna University of Technology
Bluchr Mchanical Enginring Procdings May 2014, vol. 1, num. 1 www.procdings.bluchr.com.br/vnto/10wccm TOPOLOGY DESIG OF STRUCTURE LOADED BY EARTHQUAKE P. Rosko 1 1 Cntr of Mchanics and Structural Dynamics,
More informationThe Equitable Dominating Graph
Intrnational Journal of Enginring Rsarch and Tchnology. ISSN 0974-3154 Volum 8, Numbr 1 (015), pp. 35-4 Intrnational Rsarch Publication Hous http://www.irphous.com Th Equitabl Dominating Graph P.N. Vinay
More informationA New Approach to the Fatigue Life Prediction for Notched Components Under Multiaxial Cyclic Loading. Zhi-qiang TAO and De-guang SHANG *
2017 2nd Intrnational Conrnc on Applid Mchanics, Elctronics and Mchatronics Enginring (AMEME 2017) ISBN: 978-1-60595-497-4 A Nw Approach to th Fatigu Li Prdiction or Notchd Componnts Undr Multiaxial Cyclic
More informationu 3 = u 3 (x 1, x 2, x 3 )
Lctur 23: Curvilinar Coordinats (RHB 8.0 It is oftn convnint to work with variabls othr than th Cartsian coordinats x i ( = x, y, z. For xampl in Lctur 5 w mt sphrical polar and cylindrical polar coordinats.
More informationFinite Element Model of a Ferroelectric
Excrpt from th Procdings of th COMSOL Confrnc 200 Paris Finit Elmnt Modl of a Frrolctric A. Lópz, A. D Andrés and P. Ramos * GRIFO. Dpartamnto d Elctrónica, Univrsidad d Alcalá. Alcalá d Hnars. Madrid,
More informationAddition of angular momentum
Addition of angular momntum April, 0 Oftn w nd to combin diffrnt sourcs of angular momntum to charactriz th total angular momntum of a systm, or to divid th total angular momntum into parts to valuat th
More informationarxiv: v1 [physics.comp-ph] 30 Jun 2016
On anisotropy function in crystal growth simulations using Lattic Boltzmann quation AMINA YOUNSI a,1, ALAIN CARTALADE a, a Dn DM2S, STMF, LMSF, CEA, Univrsité d Paris-Saclay, F-91191, Gif-sur-Yvtt, Franc.
More informationQuasi-Classical States of the Simple Harmonic Oscillator
Quasi-Classical Stats of th Simpl Harmonic Oscillator (Draft Vrsion) Introduction: Why Look for Eignstats of th Annihilation Oprator? Excpt for th ground stat, th corrspondnc btwn th quantum nrgy ignstats
More informationFinite element discretization of Laplace and Poisson equations
Finit lmnt discrtization of Laplac and Poisson quations Yashwanth Tummala Tutor: Prof S.Mittal 1 Outlin Finit Elmnt Mthod for 1D Introduction to Poisson s and Laplac s Equations Finit Elmnt Mthod for 2D-Discrtization
More information22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches.
Subjct Chmistry Papr No and Titl Modul No and Titl Modul Tag 8/ Physical Spctroscopy / Brakdown of th Born-Oppnhimr approximation. Slction ruls for rotational-vibrational transitions. P, R branchs. CHE_P8_M
More informationA Simple Formula for the Hilbert Metric with Respect to a Sub-Gaussian Cone
mathmatics Articl A Simpl Formula for th Hilbrt Mtric with Rspct to a Sub-Gaussian Con Stéphan Chrétin 1, * and Juan-Pablo Ortga 2 1 National Physical Laboratory, Hampton Road, Tddinton TW11 0LW, UK 2
More informationNONCONFORMING FINITE ELEMENTS FOR REISSNER-MINDLIN PLATES
NONCONFORMING FINITE ELEMENTS FOR REISSNER-MINDLIN PLATES C. CHINOSI Dipartimnto di Scinz Tcnologi Avanzat, Univrsità dl Pimont Orintal, Via Bllini 5/G, 5 Alssandria, Italy E-mail: chinosi@mfn.unipmn.it
More informationCHAPTER 1. Introductory Concepts Elements of Vector Analysis Newton s Laws Units The basis of Newtonian Mechanics D Alembert s Principle
CHPTER 1 Introductory Concpts Elmnts of Vctor nalysis Nwton s Laws Units Th basis of Nwtonian Mchanics D lmbrt s Principl 1 Scinc of Mchanics: It is concrnd with th motion of matrial bodis. odis hav diffrnt
More informationTwo Products Manufacturer s Production Decisions with Carbon Constraint
Managmnt Scinc and Enginring Vol 7 No 3 pp 3-34 DOI:3968/jms9335X374 ISSN 93-34 [Print] ISSN 93-35X [Onlin] wwwcscanadant wwwcscanadaorg Two Products Manufacturr s Production Dcisions with Carbon Constraint
More informationDevelopments in Geomathematics 5
~~": ~ L " r. :.. ~,.!.-. r,:... I Rprintd from I Dvlopmnts in Gomathmatics 5 Procdings of th Intrnational Symposium on Variational Mthods in ' Goscincs hld at th Univrsity of Oklahoma, Norman, Oklahoma,
More information3 Finite Element Parametric Geometry
3 Finit Elmnt Paramtric Gomtry 3. Introduction Th intgral of a matrix is th matrix containing th intgral of ach and vry on of its original componnts. Practical finit lmnt analysis rquirs intgrating matrics,
More informationOn spanning trees and cycles of multicolored point sets with few intersections
On spanning trs and cycls of multicolord point sts with fw intrsctions M. Kano, C. Mrino, and J. Urrutia April, 00 Abstract Lt P 1,..., P k b a collction of disjoint point sts in R in gnral position. W
More informationCOHORT MBA. Exponential function. MATH review (part2) by Lucian Mitroiu. The LOG and EXP functions. Properties: e e. lim.
MTH rviw part b Lucian Mitroiu Th LOG and EXP functions Th ponntial function p : R, dfind as Proprtis: lim > lim p Eponntial function Y 8 6 - -8-6 - - X Th natural logarithm function ln in US- log: function
More informationAddition of angular momentum
Addition of angular momntum April, 07 Oftn w nd to combin diffrnt sourcs of angular momntum to charactriz th total angular momntum of a systm, or to divid th total angular momntum into parts to valuat
More informationEstimation of apparent fraction defective: A mathematical approach
Availabl onlin at www.plagiarsarchlibrary.com Plagia Rsarch Library Advancs in Applid Scinc Rsarch, 011, (): 84-89 ISSN: 0976-8610 CODEN (USA): AASRFC Estimation of apparnt fraction dfctiv: A mathmatical
More informationDirect Discontinuous Galerkin Method and Its Variations for Second Order Elliptic Equations
DOI 10.1007/s10915-016-0264-z Dirct Discontinuous Galrkin Mthod and Its Variations for Scond Ordr Elliptic Equations Hongying Huang 1,2 Zhng Chn 3 Jin Li 4 Ju Yan 3 Rcivd: 15 Dcmbr 2015 / Rvisd: 10 Jun
More information2.3 Matrix Formulation
23 Matrix Formulation 43 A mor complicatd xampl ariss for a nonlinar systm of diffrntial quations Considr th following xampl Exampl 23 x y + x( x 2 y 2 y x + y( x 2 y 2 (233 Transforming to polar coordinats,
More informationRotor Stationary Control Analysis Based on Coupling KdV Equation Finite Steady Analysis Liu Dalong1,a, Xu Lijuan2,a
204 Intrnational Confrnc on Computr Scinc and Elctronic Tchnology (ICCSET 204) Rotor Stationary Control Analysis Basd on Coupling KdV Equation Finit Stady Analysis Liu Dalong,a, Xu Lijuan2,a Dpartmnt of
More informationVSMN30 FINITA ELEMENTMETODEN - DUGGA
VSMN3 FINITA ELEMENTMETODEN - DUGGA 1-11-6 kl. 8.-1. Maximum points: 4, Rquird points to pass: Assistanc: CALFEM manual and calculator Problm 1 ( 8p ) 8 7 6 5 y 4 1. m x 1 3 1. m Th isotropic two-dimnsional
More informationBifurcation Theory. , a stationary point, depends on the value of α. At certain values
Dnamic Macroconomic Thor Prof. Thomas Lux Bifurcation Thor Bifurcation: qualitativ chang in th natur of th solution occurs if a paramtr passs through a critical point bifurcation or branch valu. Local
More informationHigher-Order Discrete Calculus Methods
Highr-Ordr Discrt Calculus Mthods J. Blair Prot V. Subramanian Ralistic Practical, Cost-ctiv, Physically Accurat Paralll, Moving Msh, Complx Gomtry, Slid 1 Contxt Discrt Calculus Mthods Finit Dirnc Mimtic
More informationSimulated Analysis of Tooth Profile Error of Cycloid Steel Ball Planetary Transmission
07 4th Intrnational Matrials, Machinry and Civil Enginring Confrnc(MATMCE 07) Simulatd Analysis of Tooth Profil Error of Cycloid Stl Ball Plantary Transmission Ruixu Hu,a, Yuquan Zhang,b,*, Zhanliang Zhao,c,
More informationFinite Element Models for Steady Flows of Viscous Incompressible Fluids
Finit Elmnt Modls for Stad Flows of Viscous Incomprssibl Fluids Rad: Chaptr 10 JN Rdd CONTENTS Govrning Equations of Flows of Incomprssibl Fluids Mid (Vlocit-Prssur) Finit Elmnt Modl Pnalt Function Mthod
More informationCoupled Pendulums. Two normal modes.
Tim Dpndnt Two Stat Problm Coupld Pndulums Wak spring Two normal mods. No friction. No air rsistanc. Prfct Spring Start Swinging Som tim latr - swings with full amplitud. stationary M +n L M +m Elctron
More informationAnalysis of potential flow around two-dimensional body by finite element method
Vol. 7(2), pp. 9-22, May, 2015 DOI: 10.5897/JMER2014.0342 rticl Numbr: 20E80053033 ISSN 2141 2383 Copyright 2015 uthor(s) rtain th copyright of this articl http://www.acadmicjournals.org/jmer Journal of
More informationTHE ANALYSIS OF THE ITERATIONS PROCESS IN THE ELASTO-PLASTIC STRESS MODEL
Plas cit this articl as: Joanna Wróbl, Adam Kulawik, Th analysis of th itrations procss in th lasto-plastic strss modl, Scintific Rsarch of th Institut of Mathmatics and Computr Scinc, 2012, Volum 11,
More informationElement connectivity parameterization method for the stress based topology optimization for geometrically nonlinear structure
0 th World Congrss on Structural and Multidisciplinary Optimization May 9-4, 03, Orlando, Florida, USA Elmnt connctivity paramtrization mthod for th strss basd topology optimization for gomtrically nonlinar
More informationInternational Journal of Scientific & Engineering Research, Volume 6, Issue 7, July ISSN
Intrnational Journal of Scintific & Enginring Rsarch, Volum 6, Issu 7, July-25 64 ISSN 2229-558 HARATERISTIS OF EDGE UTSET MATRIX OF PETERSON GRAPH WITH ALGEBRAI GRAPH THEORY Dr. G. Nirmala M. Murugan
More informationDIFFERENTIAL EQUATION
MD DIFFERENTIAL EQUATION Sllabus : Ordinar diffrntial quations, thir ordr and dgr. Formation of diffrntial quations. Solution of diffrntial quations b th mthod of sparation of variabls, solution of homognous
More informationApproximation and Inapproximation for The Influence Maximization Problem in Social Networks under Deterministic Linear Threshold Model
20 3st Intrnational Confrnc on Distributd Computing Systms Workshops Approximation and Inapproximation for Th Influnc Maximization Problm in Social Ntworks undr Dtrministic Linar Thrshold Modl Zaixin Lu,
More informationThe van der Waals interaction 1 D. E. Soper 2 University of Oregon 20 April 2012
Th van dr Waals intraction D. E. Sopr 2 Univrsity of Orgon 20 pril 202 Th van dr Waals intraction is discussd in Chaptr 5 of J. J. Sakurai, Modrn Quantum Mchanics. Hr I tak a look at it in a littl mor
More informationDirect Approach for Discrete Systems One-Dimensional Elements
CONTINUUM & FINITE ELEMENT METHOD Dirct Approach or Discrt Systms On-Dimnsional Elmnts Pro. Song Jin Par Mchanical Enginring, POSTECH Dirct Approach or Discrt Systms Dirct approach has th ollowing aturs:
More informationA Family of Discontinuous Galerkin Finite Elements for the Reissner Mindlin plate
A Family of Discontinuous Galrkin Finit Elmnts for th Rissnr Mindlin plat Douglas N. Arnold Franco Brzzi L. Donatlla Marini Sptmbr 28, 2003 Abstract W dvlop a family of locking-fr lmnts for th Rissnr Mindlin
More informationApplication of Vague Soft Sets in students evaluation
Availabl onlin at www.plagiarsarchlibrary.com Advancs in Applid Scinc Rsarch, 0, (6):48-43 ISSN: 0976-860 CODEN (USA): AASRFC Application of Vagu Soft Sts in studnts valuation B. Chtia*and P. K. Das Dpartmnt
More information3-D SQCE Model and Its Application in Fracture Mechanics *
3-D SQCE Modl and Its Application in Fractur Mchanics * Zhichao Wang Sr. ad Enginr Applid Mchanics Dpt., Emrson Climat Tchnology, USA Tribikram Kundu - Profssor Enginring Mchanics Dpt.,Th Univrsity of
More informationDerangements and Applications
2 3 47 6 23 Journal of Intgr Squncs, Vol. 6 (2003), Articl 03..2 Drangmnts and Applications Mhdi Hassani Dpartmnt of Mathmatics Institut for Advancd Studis in Basic Scincs Zanjan, Iran mhassani@iasbs.ac.ir
More informationSliding Mode Flow Rate Observer Design
Sliding Mod Flow Rat Obsrvr Dsign Song Liu and Bin Yao School of Mchanical Enginring, Purdu Univrsity, Wst Lafaytt, IN797, USA liu(byao)@purdudu Abstract Dynamic flow rat information is ndd in a lot of
More informationNumerische Mathematik
Numr. Math. 04 6:3 360 DOI 0.007/s00-03-0563-3 Numrisch Mathmatik Discontinuous Galrkin and mimtic finit diffrnc mthods for coupld Stoks Darcy flows on polygonal and polyhdral grids Konstantin Lipnikov
More informationAn adaptive Strategy for the Multi-scale Analysis of Plate and Shell Structures with Elasto-plastic Material Behaviour
TECHNISCHE MECHANIK, 36, 1-2, (2016), 142 154 submittd: Sptmbr 7, 2015 An adaptiv Stratgy for th Multi-scal Analysis of Plat and Shll Structurs with Elasto-plastic Matrial Bhaviour W Wagnr, F Gruttmann
More informationScattering States of l-wave Schrödinger Equation with Modified Rosen Morse Potential
Commun. Thor. Phys. 66 06 96 00 Vol. 66, No., August, 06 Scattring Stats of l-wav Schrödingr Equation with Modifid Rosn Mors Potntial Wn-Li Chn í,, Yan-Wi Shi á, and Gao-Fng Wi Ôô, Gnral Education Cntr,
More informationData Assimilation 1. Alan O Neill National Centre for Earth Observation UK
Data Assimilation 1 Alan O Nill National Cntr for Earth Obsrvation UK Plan Motivation & basic idas Univariat (scalar) data assimilation Multivariat (vctor) data assimilation 3d-Variational Mthod (& optimal
More informationc 2017 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 55, No. 4, pp. 1719 1739 c 017 Socity for Industrial and Applid Mathmatics ON HANGING NODE CONSTRAINTS FOR NONCONFORMING FINITE ELEMENTS USING THE DOUGLAS SANTOS SHEEN YE ELEMENT
More informationL ubomír Baňas 1 and Robert Nürnberg Introduction A POSTERIORI ESTIMATES FOR THE CAHN HILLIARD EQUATION WITH OBSTACLE FREE ENERGY
ESAIM: MAN 43 (009) 1003 106 DOI: 10.1051/man/009015 ESAIM: Mathmatical Modlling and Numrical Analysis www.saim-man.org A POSERIORI ESIMAES FOR HE CAHN HILLIARD EQUAION WIH OBSACLE FREE ENERGY L ubomír
More informationSquare of Hamilton cycle in a random graph
Squar of Hamilton cycl in a random graph Andrzj Dudk Alan Friz Jun 28, 2016 Abstract W show that p = n is a sharp thrshold for th random graph G n,p to contain th squar of a Hamilton cycl. This improvs
More informationObserver Bias and Reliability By Xunchi Pu
Obsrvr Bias and Rliability By Xunchi Pu Introduction Clarly all masurmnts or obsrvations nd to b mad as accuratly as possibl and invstigators nd to pay carful attntion to chcking th rliability of thir
More informationSymmetric Interior penalty discontinuous Galerkin methods for elliptic problems in polygons
Symmtric Intrior pnalty discontinuous Galrkin mthods for lliptic problms in polygons F. Müllr and D. Schötzau and Ch. Schwab Rsarch Rport No. 2017-15 March 2017 Sminar für Angwandt Mathmatik Eidgnössisch
More informationNumerical Analysis of Transient Responses for Elastic Structures Connected to a Viscoelastic Shock Absorber Using FEM with a Nonlinear Complex Spring
Numrical Analysis of Transint Rsponss for Elastic Structurs Connctd to a Viscolastic Shock Absorbr Using FEM with a Nonlinar Complx Spring Takao Yamaguchi, Yusaku Fujii, Toru Fukushima, Akihiro Takita,
More information