Recent Advances in Time-Domain Maxwell s Equations in Metamaterials
|
|
- Diane Townsend
- 5 years ago
- Views:
Transcription
1 Recent Advances in Time-Domain Maxwell s Equations in Metamaterials Yunqing Huang 1, and Jicun Li, 1 Hunan Key Laboratory for Computation and Simulation in Science and Engineering, Xiangtan University, Cina Department of Matematical Sciences, University of Nevada, Las Vegas, Nevada , USA Abstract. In tis paper, we present tree second-order finite element metods solving Maxwell s equations wen metamaterials are involved. Te first metod is based on integral-differential equations transformed from te governing equations; te oter two metods solve te original governing differential equations. Numerical results are presented for te last two metods. 1 Introduction Te double negative (DNG metamaterials are artificially structured electromagnetic materials wit bot negative permittivity and permeability. Te successful construction of DNG metamaterials in triggered a new wave in te study and design of DNG metamaterials and exploration of applications in diverse areas suc as sub-wavelengt imaging and cloaking. Since, engineers and pysicists ave proposed some matematical models and initiated some numerical simulations for DNG metamaterials. Most metods are eiter based on te finite-difference timedomain (FDTD metods or finite element metods on frequency domain (mainly using commercial packages suc as HFSS and FEMLAB. Toug tere exist many work (e.g., [5] and references cited terein devoted to Maxwell s equations, to our best knowledge, most are for te simple media case. Tere are very few numerical analysis work wen metamaterials are involved except our recent effort in time-domain finite element metods [4,3,1]. In tis paper, we systematically present tree second-order metods: one solved as integral-differential equations; te oter two solved directly from te original differential equations eiter implicitly or explicitly. Te Governing Equations Te DNG metamaterials can be simulated using lossy Drude polarization and magnetization models. Te governing equations for modeling wave propagation in metamaterials are [4]: Supported by te NSFC for Distinguised Young Scolars (16516 and National Basic Researc Program of Cina under te grant 5CB3171. Supported by National Science Foundation grant DMS W. Zang et al. (Eds.: HPCA 9, LNCS 5938, pp , 1. c Springer-Verlag Berlin Heidelberg 1
2 Recent Advances in Time-Domain Maxwell s Equations in Metamaterials 49 E H ɛ = H J, t μ = E K, (1 t 1 J ɛ ωpe t + Γe 1 K J = E, ɛ ωpe μ ωpm t + Γm K = H, ( μ ωpm were ɛ and μ are te vacuum permittivity and permeability respectively, ω pe and ω pm are te electric and magnetic plasma frequencies respectively, Γ e and Γ m are te electric and magnetic damping frequencies respectively, E(x,tand H(x,t are te electric and magnetic fields respectively, and J(x,tandK(x,t are te induced electric and magnetic currents respectively. For simplicity, we assume tat te modeling domain be Ω (,T, were Ω is a bounded Lipscitz polyedral domain in R 3 wit connected boundary Ω. Furtermore, we assume tat te boundary of Ω is perfect conducting so tat n E = on Ω, (3 were n is te unit outward normal. Also we assume te initial conditions: E(x, = E (x, H(x, = H (x, J(x, = J (x, K(x, = K (x. Lemma 1. [1] Tere exists a unique solution for te system (1-(. Furtermore, te solution of (1-( satisfies te following stability estimate ɛ E(t + μ H(t + 1 K(t μ ωpm + 1 J(t ɛ ωpe ɛ E( + μ H( + 1 K( μ ωpm + 1 J( ɛ. (4 ωpe 3 Tree Different Numerical Metods To construct a finite element metod, we partition Ω by a family of regular tetraedral meses T wit maximum mes size. For simplicity, consider te lowest order Raviart-Tomas-Nédélec s spaces: for any K T, U = {u H(div; Ω u K = c K + d Kx}, V = {v H(curl; Ω v K = a K + b K x}, V = V {n v = on Ω}, were a K, b K, c K are constant vectors in R 3, and d K is a real constant. For any v H α (curl; Ω, 1 <α 1, it is known tat we can define te Nédélec interpolant Π v V on eac tetraedron K and ave ([5]: v Π v + (v Π v C α v α,curl. (5 Denoting by P u U te standard (L (Ω 3 -projection onto U,weave u P u C α u α u H α (Ω. (6 To define a fully discrete sceme, we divide te time interval (,TintoM uniform subintervals by points t k = k, were = T M. Moreover, we denote te k-t subinterval =[t k 1,t k ], and define u k = u(,k.
3 5 Y. Huang and J. Li 3.1 Treated as Integral-Differential Equations Solving ( wit initial currents J (x andk (x, we obtain J(x,t; E =e Γet J (x+ɛ ω pe t K(x,t; H =e Γmt K (x+μ ω pm e Γe(t s E(x,sds J N (x+j N(E, (7 t e Γm(t s H(x,sds K N (x+k N (H. Hence, te governing equations (1-( become: find (E, H suc tat E ɛ t H + J N(E = J N (x (x,t Ω (,T, (8 H μ t + E + K N (H = K N (x (x,t Ω (,T. (9 Now we can formulate a Crank-Nicolson mixed finite element sceme for (8-(9 as follows: for k =1,,,M, find E k V, Hk U suc tat ɛ (δ E k, φ (H k, φ +(J k, φ = ( J N, φ, (1 μ (δ H k, ψ +( E k, ψ +(K k, ψ = ( K N, ψ (11 for any φ V, ψ U and <t T, subject to te initial conditions E (x =Π E (x and H (x =P H (x. (1 Here we denote ũ = 1 I u(tdt for any u. Te J k k is defined as J k = 1 (1 + e Γe J k 1 J =, J k = e Γe J k 1 Similar formula will be used for K k. + ɛ ω pe 1 (e Γe E k 1 + E k, 1 k M, (13 + ɛ ω pe (e Γe E k 1 + E k, 1 k M. (14 Lemma. [1] For S = H 1 (curl; Ω or S =(H α (Ω 3 wit α, we ave (i δ u k S 1 (ii u k 1 t k 1 u t(t Sdt u(tdt S 3 t k 1 4 u H 1 (,T; S, t k 1 u tt(t Sdt u H (,T; S. be defined by (7 and (13, respec- Lemma 3. Let J k N J N (E(,t k and J k tively. Ten for any 1 k M, we ave J k J k N C[ E l E l + 4 ( E(t + E t (t + E t (t dt]. l=
4 Recent Advances in Time-Domain Maxwell s Equations in Metamaterials 51 Proof. Recall tat J N (E =ɛ ω pe t e Γe(t s E(sds. For clarity, we denote J k = J N (E(t k. By definitions of J and J k,weave J k J k = 1 (1 + e Γe (J k 1 J k 1 +ɛ ωpe [ 1 e Γ e(t k s E(sds 1 e t k 1 (e Γ E k 1 + E k ] +ɛ ωpe 1 [(Ek E k +e Γe (E k 1 E k 1 ] J k 1 J k 1 + ɛ ωpe t k 1 s (e Γe(t k s E(s ds +ɛ ωpe 4 ( Ek E k + E k 1 E k 1, (15 were in te last step we used te inequality 1 (uk + u k 1 1 u(tdt t k 1 for u(s =e Γe(tk s E(s. From [, (Eq. (7], for k, we ave l= t k 1 u tt (t dt k 1 1 J k 1 J k 1 ɛ ωpe[ E l E l + C ( E + E s + E ss ds]. (16 Wen k =1,weaveJ = J =. Hence (16 olds true for k =1. Combining (15-(16, we obtain (for any k 1 J k J k C[ E l E l + ( E + E t + E tt dt], (17 l= were we absorbed te dependence of ɛ,ω pe,γ e into te generic constant C. Squaring bot sides of (17 and using Caucy-Scwarz inequality leads to J k J k C[ ( C[ E l E l + 4 ( l= ( E + E t + E tt dt ] E l E l + 4 ( E + E t + E tt dt] l= wic concludes te proof. Teorem 1. Let (E n, H n and (E n, H n be te solutions of (1-( and (1- (11 at time t = t n, respectively. Assume tat for 1 <α 1, E, E t L (,T; H α (curl; Ω, H, H t L (,T; H α (Ω, E tt, E tt, H tt, H tt L (,T;(L (Ω 3. Ten tere is a constant C = C(T,ɛ,μ,ω pe,ω pm,γ e,γ m, E, H, independent of bot te time step and te mes size, suc tat max 1 n M ( En E n + H n H n C( + α.
5 5 Y. Huang and J. Li Proof. Multiplying (8-(9 by test functions φ and ψ, integrating te resultants in time over, ten coosing φ = 1 φ, ψ = 1 ψ and subtracting (1 and (11, respectively, we obtain te error equations ɛ (δ ξ, k φ (η k, φ =ɛ (δ (Π E k E k, φ (P H k 1 H(sds, φ +(J k 1 J N (sds, φ, (18 μ (δ η, k ψ +( ξ k, ψ =μ (δ (P H k H k, ψ +( (Π E k 1 E(sds, ψ +(K k 1 K N (sds, ψ, (19 were we introduced te notation ξ k = Π E k E k and ηk = P H k H k. Coosing φ = (ξ k + ξk 1 =ξ k, ψ = (η k + ηk 1 =η k in (18-(19, and adding te resultants togeter, we ave ɛ ( ξ k ξ k 1 +μ ( η k η k 1 =ɛ (δ (Π E k E k, ξ k (P H k 1 H(sds, ξ k +(J k 1 J N (sds, ξ k +μ (δ (P H k H k, η k +( (Π E k 1 E(sds, η k +(K k 1 K N (sds, η k 6 = (Err i. i=1 Using te inequality ab δa + 1 4δ b and Lemma, we ave ( (Err 1 δ 1ɛ ξ k + ɛ δ (Π E k E k δ 1 δ 1ɛ ( ξ k + ξ k 1 + ɛ t(π E k E k dt. (1 δ 1 Similarly, by te projection property of P and integration by parts, we ave (Err = ( (H k 1 H(sds, ξ k δ ξ k + (H k 1 H(sds, δ δ ( ξ k + ξ k H tt(t dt. ( 8δ Similarly, by Lemmas and 3, we can obtain (Err 3 =(J k J k N + J k N 1 J N (sds, ξ k δ 3( ξ k + ξ k 1 + ( J k J k δ N + J k N 1 J N (sds 3 δ 3( ξ k + ξ k 1 + C E l E l δ 3 l=
6 Recent Advances in Time-Domain Maxwell s Equations in Metamaterials 53 + C5 tk ( E + E t + E tt δ dt + 4 J N 3 4δ 3 t ds δ 3( ξ k + ξ k 1 + C ( ξ δ l + Π E l E l 3 l= + C5 tk ( E + E t + E tt δ dt + 4 J N 3 4δ 3 t ds. (3 Similarly, we can easily obtain (Err 4 δ 4μ ( η k + η k 1 + μ t(p H k H k dt, δ 4 (Err 6 δ 6( η k + η k 1 + C ( η l + P H l H l δ 6 l= + C5 tk ( H + H t + H tt dt + 4 K N ds, δ 6 4δ 6 t (Err 5 =( (Π E k E k + E k 1 E(sds, η k δ 5 η k + ( (E k 1 E(sds + Π E k E k δ 5 δ 5( η k + η k ( E 4δ 5 t dt + Π E k E k δ. 5 Finally, summing up bot sides of ( from k =1uptoanyn M, usingte estimates of (Err i, te fact tat ξ = η =, and coosing δ i small enoug so tat ξ n and ηn can be controlled by te left and side, we ave ɛ ξ n + μ η n n 1 C 1 ( ξ k + η k + ɛ δ δ + C4 δ 3 + μ δ δ 5 + C4 δ 6 k=1 tn tn n k=1 tn tn n k=1 H tt(t dt + C δ 3 t(π E k E k dt n Π E l E l l=1 ( E + E t + E tt dt + 4 4δ 3 tn t(p H k H k dt + δ 5 Π E k E k ( E t dt + C δ 6 n P H k H k k=1 ( H + H t + H tt dt + 4 4δ 6 tn J N t ds K N t ds, (4 wereweusedtefactk T and absorbed T into te generic constant C. Using interpolation estimates (5 and (6 to (4 and te inequality J T N t C( E(s ds + E(t + E t(t, (5
7 54 Y. Huang and J. Li we obtain n 1 ɛ ξ n + μ η n C 1 ( ξ k + η k +C ( 4 + α, (6 k=1 wic concludes te proof by te standard tecnique [1]. 3. Treated as Differential Equations: Crank-Nicolson Metod Recently, we found tat we can solve (1-( directly and efficiently by constructing a Crank-Nicolson sceme[3]: for k =1,,,M, find E k V, J k V, H k, K k U suc tat ɛ (δ E k, φ (H k, φ +(J k, φ =, (7 μ (δ H k, ψ +( E k, ψ +(K k, ψ =, (8 1 (δ ɛ ωpe J k, φ + Γe (J k, φ ɛ ωpe =(E k, φ, (9 1 (δ μ ωpm K k, ψ + Γm (K k μ, ψ ωpm =(H k, ψ, (3 are true for any φ, φ V, and ψ, ψ U. In practical implementation, at eac time step we just need to first solve a smaller system ( ɛ =( ɛ ( μ =( μ + ɛω pe +Γ e (E k, φ (H k, φ ɛω pe (E k 1 +Γ e, φ +(H k 1, φ + μω pm +Γ m (H k, ψ +( E k, ψ μω pm 4 +Γ e (J k 1, φ, (31 (H k 1, ψ +Γ ( E k 1 4, ψ (K k 1, ψ m +Γ. (3 m for E k and H k,tenupdatej k and K k as follows: J k = ɛω pe (E k + E k 1 + Γe J k 1, +Γ e +Γ e (33 K k = μω pm (H k + H k 1 + Γm K k 1. +Γ m +Γ m (34 In [3], we proved te following optimal error estimate: max 1 n M ( En E n + H n H n + J n J n + K n K n C( + α.
8 Recent Advances in Time-Domain Maxwell s Equations in Metamaterials Treated as Differential Equations: Leap-Frog Metod Inspiring from te classic FDTD sceme, we can formulate a leap-frog sceme: Given initial approximations E, K, H 1, J 1, for k = 1,,, find E k V, J k+ 1 V, H k+ 1, K k U suc tat ɛ ( Ek E k 1, φ (H k 1, φ +(J k 1, φ =, (35 μ H k+ 1 H k 1 (, ψ +( E k, ψ +(K k, ψ =, (36 1 ( J k+ 1 J k 1, φ ɛ ωpe + Γe ( 1 ɛ ωpe (J k+ 1 + J k 1, φ =(E k, φ, (37 1 μ, ψ ωpm + Γm ( 1 μ ωpm (K k + K k 1, ψ =(H k 1, ψ, (38 for any φ V, ψ U, φ V, ψ U. Conditional stability and te following optimal error estimate were obtained in [1]: max ( E n E n + H n+ 1 n+ 1 H + J n+ 1 n+ J 1 + K n K n 1 n C( + α +C ( E E + H 1 H 1 + J 1 J 1 + K K. 4 Numerical Examples Here we sow some numerical results obtained by te Crank-Nicolson (CN sceme (7-(3 and Leap-Frog (LP sceme (35-(38. We implemented our algoritms using te lowest order rectangular edge element basis functions [1]. Our tests were carried out using MATLAB 7. running on Dell Latitude D63 laptop wit 1GB memory and. GHz CPU. In order to test our algoritm wit an analytical solution, we add source terms to te original governing equations, i.e., te LP sceme (35-(38 becomes: (E k, φ =(E k 1, φ + ɛ [(H k 1, φ (J k 1, φ +(f k 1, φ ], (39 (H k+ 1, ψ =(H k 1, ψ μ [( E k, ψ +(K k, ψ (g k, ψ ], (4 J k+ 1 = ɛω pe E k + Γe J k 1, (41 +Γ e +Γ e K k = μω pm H k 1 + Γm K k 1, (4 +Γ m +Γ m were f k 1 and g k are added source terms. Wile for CN sceme, we add f k and g k to te rigt and side of (7 and (8, respectively. To compare te performance between CN and LP scemes, we used te exact solution developed for te -D transverse electrical model in [1]. Troug many
9 56 Y. Huang and J. Li Table 1. L errors obtained by te CN sceme wit =1 8 Time steps nt =1 nt = 1 DOF Mes sizes E x H z E x H z Time (sec e e e e e e e e e e e e e e e e e e e e Table. L errors obtained by te CN sceme wit =1 1 Time steps nt =1 nt = 1 DOF Mes sizes E x H z E x H z Time (sec e e e e e e e e e e e e e e e e e e e e Fig. 1. Numerical solution on mes: (Left electric field; (Rigt magnetic field tests, we found tat even for CN sceme good results can only be acieved until =1 8 or smaller. Our numerical results obtained wit CN sceme using =1 8 and =1 1 are presented in Table 1 and Table, respectively. Also presented are te total number of degree of freedom (DOF and computational time. Tables 1- sow clearly O( convergence for E, wic is superconvergent for rectangular meses [1]. Te accuracies for H are really excellent, but it sows O( as mes gets finer. We tink tis is due to te larger sampling errors coming wit finer meses since H are approximated by piecewise constants.
10 Recent Advances in Time-Domain Maxwell s Equations in Metamaterials 57 Comparing Tables 1- to Tables -3 of [1] (results obtained by LP sceme, we see tat wit te same mes and time step CN sceme acieves similar accuracy for E as LP sceme and as muc more accurate results for H tan LP sceme. However, CN sceme takes muc longer computational time tan LP sceme. 5 Concluding Remarks In tis paper, we presented tree second order metods for solving te metamaterial Maxwell s equations in time-domain. Considering te accuracy and computational cost, it seems tat te LP sceme is te best coice of te tree metods. More rigorous comparisons will be carried out in our future work. For practical applications, we ave to work on ig performance computers, since our laptop runs out of memory for 3 3 mes. References 1. Li, J.: Numerical convergence and pysical fidelity analysis for Maxwell s equations in metamaterials. Comput. Metods Appl. Mec. Engrg. 198, (9. Li, J., Cen, Y.: Analysis of a time-domain finite element metod for 3-D Maxwell s equations in dispersive media. Comput. Metods Appl. Mec. Engrg. 195, 4 49 (6 3. Li, J., Cen, Y., Elander, V.: Matematical and numerical study of wave propagation in negative-index materials. Comput. Metods Appl. Mec. Engrg. 197, (8 4. Li, J., Wood, A.: Finite element analysis for wave propagation in double negative metamaterials. J. Sci. Comput. 3, (7 5. Monk, P.: Finite Element Metods for Maxwell s Equations. Oxford University Press, Oxford (3
1. Introduction. We consider the model problem: seeking an unknown function u satisfying
A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced
More informationNumerical Study of the Plasma-Lorentz Model in Metamaterials
J Sci Comput 3 54: 44 DOI.7/s95--968-5 Numerical Study of te Plasma-Lorentz Model in Metamaterials Jicun Li Yunqing Huang Wei Yang Received: 7 February / Revised: 3 May / Accepted: 3 May / Publised online:
More informationNumerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems
Applied Matematics, 06, 7, 74-8 ttp://wwwscirporg/journal/am ISSN Online: 5-7393 ISSN Print: 5-7385 Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for
More informationAnalysis of A Continuous Finite Element Method for H(curl, div)-elliptic Interface Problem
Analysis of A Continuous inite Element Metod for Hcurl, div)-elliptic Interface Problem Huoyuan Duan, Ping Lin, and Roger C. E. Tan Abstract In tis paper, we develop a continuous finite element metod for
More informationA finite element approximation for the quasi-static Maxwell Landau Lifshitz Gilbert equations
ANZIAM J. 54 (CTAC2012) pp.c681 C698, 2013 C681 A finite element approximation for te quasi-static Maxwell Landau Lifsitz Gilbert equations Kim-Ngan Le 1 T. Tran 2 (Received 31 October 2012; revised 29
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationA Weak Galerkin Method with an Over-Relaxed Stabilization for Low Regularity Elliptic Problems
J Sci Comput (07 7:95 8 DOI 0.007/s095-06-096-4 A Weak Galerkin Metod wit an Over-Relaxed Stabilization for Low Regularity Elliptic Problems Lunji Song, Kaifang Liu San Zao Received: April 06 / Revised:
More informationAN ANALYSIS OF NEW FINITE ELEMENT SPACES FOR MAXWELL S EQUATIONS
Journal of Matematical Sciences: Advances and Applications Volume 5 8 Pages -9 Available at ttp://scientificadvances.co.in DOI: ttp://d.doi.org/.864/jmsaa_7975 AN ANALYSIS OF NEW FINITE ELEMENT SPACES
More informationCONVERGENCE ANALYSIS OF YEE SCHEMES FOR MAXWELL S EQUATIONS IN DEBYE AND LORENTZ DISPERSIVE MEDIA
INTRNATIONAL JOURNAL OF NUMRICAL ANALYSIS AND MODLING Volume XX Number 0 ages 45 c 03 Institute for Scientific Computing and Information CONVRGNC ANALYSIS OF Y SCHMS FOR MAXWLL S QUATIONS IN DBY AND LORNTZ
More informationCONVERGENCE ANALYSIS OF YEE SCHEMES FOR MAXWELL S EQUATIONS IN DEBYE AND LORENTZ DISPERSIVE MEDIA
INTRNATIONAL JOURNAL OF NUMRICAL ANALYSIS AND MODLING Volume Number 4 ages 657 687 c 04 Institute for Scientific Computing and Information CONVRGNC ANALYSIS OF Y SCHMS FOR MAXWLL S QUATIONS IN DBY AND
More informationNumerische Mathematik
Numer. Mat. (1999 82: 193 219 Numerisce Matematik c Springer-Verlag 1999 Electronic Edition Fully discrete finite element approaces for time-dependent Maxwell s equations P. Ciarlet, Jr 1, Jun Zou 2, 1
More informationarxiv: v1 [math.na] 12 Mar 2018
ON PRESSURE ESTIMATES FOR THE NAVIER-STOKES EQUATIONS J A FIORDILINO arxiv:180304366v1 [matna 12 Mar 2018 Abstract Tis paper presents a simple, general tecnique to prove finite element metod (FEM) pressure
More informationMIXED DISCONTINUOUS GALERKIN APPROXIMATION OF THE MAXWELL OPERATOR. SIAM J. Numer. Anal., Vol. 42 (2004), pp
MIXED DISCONTINUOUS GALERIN APPROXIMATION OF THE MAXWELL OPERATOR PAUL HOUSTON, ILARIA PERUGIA, AND DOMINI SCHÖTZAU SIAM J. Numer. Anal., Vol. 4 (004), pp. 434 459 Abstract. We introduce and analyze a
More informationSuperconvergence of energy-conserving discontinuous Galerkin methods for. linear hyperbolic equations. Abstract
Superconvergence of energy-conserving discontinuous Galerkin metods for linear yperbolic equations Yong Liu, Ci-Wang Su and Mengping Zang 3 Abstract In tis paper, we study superconvergence properties of
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More information4.2 - Richardson Extrapolation
. - Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Definition Let x n n converge to a number x. Suppose tat n n is a sequence
More informationUniversity Mathematics 2
University Matematics 2 1 Differentiability In tis section, we discuss te differentiability of functions. Definition 1.1 Differentiable function). Let f) be a function. We say tat f is differentiable at
More informationPreconditioning in H(div) and Applications
1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition
More informationarxiv: v1 [math.na] 17 Jul 2014
Div First-Order System LL* FOSLL* for Second-Order Elliptic Partial Differential Equations Ziqiang Cai Rob Falgout Sun Zang arxiv:1407.4558v1 [mat.na] 17 Jul 2014 February 13, 2018 Abstract. Te first-order
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationDownloaded 11/15/17 to Redistribution subject to SIAM license or copyright; see
SIAM J. NUMER. ANAL. Vol. 55, No. 6, pp. 2787 2810 c 2017 Society for Industrial and Applied Matematics EDGE ELEMENT METHOD FOR OPTIMAL CONTROL OF STATIONARY MAXWELL SYSTEM WITH GAUSS LAW IRWIN YOUSEPT
More informationDifferentiation in higher dimensions
Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends
More informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More informationDiscontinuous Galerkin Methods for Relativistic Vlasov-Maxwell System
Discontinuous Galerkin Metods for Relativistic Vlasov-Maxwell System He Yang and Fengyan Li December 1, 16 Abstract e relativistic Vlasov-Maxwell (RVM) system is a kinetic model tat describes te dynamics
More information= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)
Paper 1: Pure Matematics 1 Mark Sceme 1(a) (i) (ii) d d y 3 1x 4x x M1 A1 d y dx 1.1b 1.1b 36x 48x A1ft 1.1b Substitutes x = into teir dx (3) 3 1 4 Sows d y 0 and states ''ence tere is a stationary point''
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationA THIRD ORDER LINEARIZED BDF SCHEME FOR MAXWELL S EQUATIONS WITH NONLINEAR CONDUCTIVITY USING FINITE ELEMENT METHOD
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 4, Number 4-5, Pages 5 5 c 7 Institute for Scientific Computing and Information A THIRD ORDER LINEARIZED BDF SCHEME FOR MAXWELL S EQUATIONS
More informationMath Spring 2013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, (1/z) 2 (1/z 1) 2 = lim
Mat 311 - Spring 013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, 013 Question 1. [p 56, #10 (a)] 4z Use te teorem of Sec. 17 to sow tat z (z 1) = 4. We ave z 4z (z 1) = z 0 4 (1/z) (1/z
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More informationERROR ESTIMATES FOR A FULLY DISCRETIZED SCHEME TO A CAHN-HILLIARD PHASE-FIELD MODEL FOR TWO-PHASE INCOMPRESSIBLE FLOWS
ERROR ESTIMATES FOR A FULLY DISCRETIZED SCHEME TO A CAHN-HILLIARD PHASE-FIELD MODEL FOR TWO-PHASE INCOMPRESSIBLE FLOWS YONGYONG CAI, AND JIE SHEN Abstract. We carry out in tis paper a rigorous error analysis
More informationFourier Type Super Convergence Study on DDGIC and Symmetric DDG Methods
DOI 0.007/s095-07-048- Fourier Type Super Convergence Study on DDGIC and Symmetric DDG Metods Mengping Zang Jue Yan Received: 7 December 06 / Revised: 7 April 07 / Accepted: April 07 Springer Science+Business
More informationA Hybrid Mixed Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Hybrid Mixed Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationA SPLITTING LEAST-SQUARES MIXED FINITE ELEMENT METHOD FOR ELLIPTIC OPTIMAL CONTROL PROBLEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALSIS AND MODELING Volume 3, Number 4, Pages 6 626 c 26 Institute for Scientific Computing and Information A SPLITTING LEAST-SQUARES MIED FINITE ELEMENT METHOD FOR
More informationA SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION
A SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION GABRIEL R. BARRENECHEA, LEOPOLDO P. FRANCA 1 2, AND FRÉDÉRIC VALENTIN Abstract. Tis work introduces and analyzes novel stable
More informationarxiv: v1 [math.na] 27 Jan 2014
L 2 -ERROR ESTIMATES FOR FINITE ELEMENT APPROXIMATIONS OF BOUNDARY FLUXES MATS G. LARSON AND ANDRÉ MASSING arxiv:1401.6994v1 [mat.na] 27 Jan 2014 Abstract. We prove quasi-optimal a priori error estimates
More informationA Mixed-Hybrid-Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Mixed-Hybrid-Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationEfficient, unconditionally stable, and optimally accurate FE algorithms for approximate deconvolution models of fluid flow
Efficient, unconditionally stable, and optimally accurate FE algoritms for approximate deconvolution models of fluid flow Leo G. Rebolz Abstract Tis paper addresses an open question of ow to devise numerical
More informationResearch Article Error Analysis for a Noisy Lacunary Cubic Spline Interpolation and a Simple Noisy Cubic Spline Quasi Interpolation
Advances in Numerical Analysis Volume 204, Article ID 35394, 8 pages ttp://dx.doi.org/0.55/204/35394 Researc Article Error Analysis for a Noisy Lacunary Cubic Spline Interpolation and a Simple Noisy Cubic
More informationFinite Difference Methods Assignments
Finite Difference Metods Assignments Anders Söberg and Aay Saxena, Micael Tuné, and Maria Westermarck Revised: Jarmo Rantakokko June 6, 1999 Teknisk databeandling Assignment 1: A one-dimensional eat equation
More informationSolution. Solution. f (x) = (cos x)2 cos(2x) 2 sin(2x) 2 cos x ( sin x) (cos x) 4. f (π/4) = ( 2/2) ( 2/2) ( 2/2) ( 2/2) 4.
December 09, 20 Calculus PracticeTest s Name: (4 points) Find te absolute extrema of f(x) = x 3 0 on te interval [0, 4] Te derivative of f(x) is f (x) = 3x 2, wic is zero only at x = 0 Tus we only need
More informationGrad-div stabilization for the evolutionary Oseen problem with inf-sup stable finite elements
Noname manuscript No. will be inserted by te editor Grad-div stabilization for te evolutionary Oseen problem wit inf-sup stable finite elements Javier de Frutos Bosco García-Arcilla Volker Jon Julia Novo
More informationSMAI-JCM SMAI Journal of Computational Mathematics
SMAI-JCM SMAI Journal of Computational Matematics Compatible Maxwell solvers wit particles II: conforming and non-conforming 2D scemes wit a strong Faraday law Martin Campos Pinto & Eric Sonnendrücker
More informationPoisson Equation in Sobolev Spaces
Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on
More informationarxiv: v1 [math.na] 20 Nov 2018
An HDG Metod for Tangential Boundary Control of Stokes Equations I: Hig Regularity Wei Gong Weiwei Hu Mariano Mateos Jon R. Singler Yangwen Zang arxiv:1811.08522v1 [mat.na] 20 Nov 2018 November 22, 2018
More informationA Distributed Lagrange Multiplier Based Fictitious Domain Method for Maxwell s Equations
A Distributed Lagrange Multiplier Based Fictitious Domain Metod for Maxwell s Equations V. A. Bokil a and R. Glowinski b Center for Researc in Scientific Computation a Nort Carolina State University Box
More informationExercises for numerical differentiation. Øyvind Ryan
Exercises for numerical differentiation Øyvind Ryan February 25, 2013 1. Mark eac of te following statements as true or false. a. Wen we use te approximation f (a) (f (a +) f (a))/ on a computer, we can
More informationarxiv: v1 [math.na] 28 Apr 2017
THE SCOTT-VOGELIUS FINITE ELEMENTS REVISITED JOHNNY GUZMÁN AND L RIDGWAY SCOTT arxiv:170500020v1 [matna] 28 Apr 2017 Abstract We prove tat te Scott-Vogelius finite elements are inf-sup stable on sape-regular
More informationParameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationA UNIFORM INF SUP CONDITION WITH APPLICATIONS TO PRECONDITIONING
A UNIFORM INF SUP CONDIION WIH APPLICAIONS O PRECONDIIONING KEN ANDRE MARDAL, JOACHIM SCHÖBERL, AND RAGNAR WINHER Abstract. A uniform inf sup condition related to a parameter dependent Stokes problem is
More informationFunction Composition and Chain Rules
Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function
More informationLinearized Primal-Dual Methods for Linear Inverse Problems with Total Variation Regularization and Finite Element Discretization
Linearized Primal-Dual Metods for Linear Inverse Problems wit Total Variation Regularization and Finite Element Discretization WENYI TIAN XIAOMING YUAN September 2, 26 Abstract. Linear inverse problems
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More informationCONVERGENCE ANALYSIS OF FINITE ELEMENT SOLUTION OF ONE-DIMENSIONAL SINGULARLY PERTURBED DIFFERENTIAL EQUATIONS ON EQUIDISTRIBUTING MESHES
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume, Number, Pages 57 74 c 5 Institute for Scientific Computing and Information CONVERGENCE ANALYSIS OF FINITE ELEMENT SOLUTION OF ONE-DIMENSIONAL
More informationMore on generalized inverses of partitioned matrices with Banachiewicz-Schur forms
More on generalized inverses of partitioned matrices wit anaciewicz-scur forms Yongge Tian a,, Yosio Takane b a Cina Economics and Management cademy, Central University of Finance and Economics, eijing,
More information3.4 Worksheet: Proof of the Chain Rule NAME
Mat 1170 3.4 Workseet: Proof of te Cain Rule NAME Te Cain Rule So far we are able to differentiate all types of functions. For example: polynomials, rational, root, and trigonometric functions. We are
More informationAPPROXIMATION BY QUADRILATERAL FINITE ELEMENTS
MATHEMATICS OF COMPUTATION Volume 71, Number 239, Pages 909 922 S 0025-5718(02)01439-4 Article electronically publised on Marc 22, 2002 APPROXIMATION BY QUADRILATERAL FINITE ELEMENTS DOUGLAS N. ARNOLD,
More informationLIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT
LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as
More information(4.2) -Richardson Extrapolation
(.) -Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Suppose tat lim G 0 and lim F L. Te function F is said to converge to L as
More informationINTERIOR PENALTY DISCONTINUOUS GALERKIN METHODS WITH IMPLICIT TIME-INTEGRATION TECHNIQUES FOR NONLINEAR PARABOLIC EQUATIONS
INTERIOR PENALTY DISCONTINUOUS GALERKIN METHODS WITH IMPLICIT TIME-INTEGRATION TECHNIQUES FOR NONLINEAR PARABOLIC EQUATIONS LUNJI SONG,, GUNG-MIN GIE 3, AND MING-CHENG SHIUE 4 Scool of Matematics and Statistics,
More information3 Parabolic Differential Equations
3 Parabolic Differential Equations 3.1 Classical solutions We consider existence and uniqueness results for initial-boundary value problems for te linear eat equation in Q := Ω (, T ), were Ω is a bounded
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationAN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS
AN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS BERNARDO COCKBURN, JOHNNY GUZMÁN, SEE-CHEW SOON, AND HENRYK K. STOLARSKI Abstract. Te embedded discontinuous
More informationArbitrary order exactly divergence-free central discontinuous Galerkin methods for ideal MHD equations
Arbitrary order exactly divergence-free central discontinuous Galerkin metods for ideal MHD equations Fengyan Li, Liwei Xu Department of Matematical Sciences, Rensselaer Polytecnic Institute, Troy, NY
More informationOSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix
Opuscula Mat. 37, no. 6 (2017), 887 898 ttp://dx.doi.org/10.7494/opmat.2017.37.6.887 Opuscula Matematica OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS Sandra
More informationSmoothed projections in finite element exterior calculus
Smooted projections in finite element exterior calculus Ragnar Winter CMA, University of Oslo Norway based on joint work wit: Douglas N. Arnold, Minnesota, Ricard S. Falk, Rutgers, and Snorre H. Cristiansen,
More informationLEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS
SIAM J. NUMER. ANAL. c 998 Society for Industrial Applied Matematics Vol. 35, No., pp. 393 405, February 998 00 LEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS YANZHAO CAO
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationPart VIII, Chapter 39. Fluctuation-based stabilization Model problem
Part VIII, Capter 39 Fluctuation-based stabilization Tis capter presents a unified analysis of recent stabilization tecniques for te standard Galerkin approximation of first-order PDEs using H 1 - conforming
More informationROBUST MULTISCALE ITERATIVE SOLVERS FOR NONLINEAR FLOWS IN HIGHLY HETEROGENEOUS MEDIA
ROBUST MULTISCALE ITERATIVE SOLVERS FOR NONLINEAR FLOWS IN HIGHLY HETEROGENEOUS MEDIA Y. EFENDIEV, J. GALVIS, S. KI KANG, AND R.D. LAZAROV Abstract. In tis paper, we study robust iterative solvers for
More informationA Local Projection Stabilization/Continuous Galerkin Petrov Method for Incompressible Flow Problems
A Local Projection Stabilization/Continuous Galerkin Petrov Metod for Incompressible Flow Problems Naveed Amed, Volker Jon, Gunar Matties, Julia Novo Marc 22, 218 Abstract A local projection stabilization
More informationH(div) conforming and DG methods for incompressible Euler s equations
H(div) conforming and DG metods for incompressible Euler s equations Jonny Guzmán Filánder A. Sequeira Ci-Wang Su Abstract H(div) conforming and discontinuous Galerkin (DG) metods are designed for incompressible
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More information1. Introduction. Consider a semilinear parabolic equation in the form
A POSTERIORI ERROR ESTIMATION FOR PARABOLIC PROBLEMS USING ELLIPTIC RECONSTRUCTIONS. I: BACKWARD-EULER AND CRANK-NICOLSON METHODS NATALIA KOPTEVA AND TORSTEN LINSS Abstract. A semilinear second-order parabolic
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
More informationParabolic PDEs: time approximation Implicit Euler
Part IX, Capter 53 Parabolic PDEs: time approximation We are concerned in tis capter wit bot te time and te space approximation of te model problem (52.4). We adopt te metod of line introduced in 52.2.
More informationPhysically Based Modeling: Principles and Practice Implicit Methods for Differential Equations
Pysically Based Modeling: Principles and Practice Implicit Metods for Differential Equations David Baraff Robotics Institute Carnegie Mellon University Please note: Tis document is 997 by David Baraff
More informationc 2004 Society for Industrial and Applied Mathematics
SIAM J NUMER ANAL Vol 4, No, pp 86 84 c 004 Society for Industrial and Applied Matematics LEAST-SQUARES METHODS FOR LINEAR ELASTICITY ZHIQIANG CAI AND GERHARD STARKE Abstract Tis paper develops least-squares
More informationAnalysis of the grad-div stabilization for the time-dependent Navier Stokes equations with inf-sup stable finite elements
arxiv:161.517v3 [mat.na] 2 May 217 Analysis of te grad-div stabilization for te time-dependent Navier Stokes equations wit inf-sup stable finite elements Javier de Frutos Bosco García-Arcilla Volker Jon
More informationAnalysis of a second order discontinuous Galerkin finite element method for the Allen-Cahn equation
Graduate Teses and Dissertations Graduate College 06 Analysis of a second order discontinuous Galerin finite element metod for te Allen-Can equation Junzao Hu Iowa State University Follow tis and additional
More informationA PRIORI AND A POSTERIORI OF A LINEAR ELLIPTIC PROBLEM WITH DYNAMICAL BOUNDARY CONDITION.
A PRIORI AND A POSTERIORI OF A LINEAR ELLIPTIC PROBLEM WITH DYNAMICAL BOUNDARY CONDITION. Toufic El Arwadi, Séréna Dib, Toni Saya To cite tis version: Toufic El Arwadi, Séréna Dib, Toni Saya. A PRIORI
More informationarxiv: v1 [math.na] 9 Sep 2015
arxiv:509.02595v [mat.na] 9 Sep 205 An Expandable Local and Parallel Two-Grid Finite Element Sceme Yanren ou, GuangZi Du Abstract An expandable local and parallel two-grid finite element sceme based on
More informationA method of Lagrange Galerkin of second order in time. Une méthode de Lagrange Galerkin d ordre deux en temps
A metod of Lagrange Galerkin of second order in time Une métode de Lagrange Galerkin d ordre deux en temps Jocelyn Étienne a a DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, Great-Britain.
More informationGradient Descent etc.
1 Gradient Descent etc EE 13: Networked estimation and control Prof Kan) I DERIVATIVE Consider f : R R x fx) Te derivative is defined as d fx) = lim dx fx + ) fx) Te cain rule states tat if d d f gx) )
More informationHigh-Order Energy and Linear Momentum Conserving Methods for the Klein-Gordon Equation
matematics Article Hig-Order Energy and Linear Momentum Conserving Metods for te Klein-Gordon Equation He Yang Department of Matematics, Augusta University, Augusta, GA 39, USA; yang@augusta.edu; Tel.:
More informationCLEMSON U N I V E R S I T Y
A Fractional Step θ-metod for Convection-Diffusion Equations Jon Crispell December, 006 Advisors: Dr. Lea Jenkins and Dr. Vincent Ervin Fractional Step θ-metod Outline Crispell,Ervin,Jenkins Motivation
More informationERROR BOUNDS FOR FINITE-DIFFERENCE METHODS FOR RUDIN OSHER FATEMI IMAGE SMOOTHING
ERROR BOUNDS FOR FINITE-DIFFERENCE METHODS FOR RUDIN OSHER FATEMI IMAGE SMOOTHING JINGYUE WANG AND BRADLEY J. LUCIER Abstract. We bound te difference between te solution to te continuous Rudin Oser Fatemi
More informationFinite Element Methods for Linear Elasticity
Finite Element Metods for Linear Elasticity Ricard S. Falk Department of Matematics - Hill Center Rutgers, Te State University of New Jersey 110 Frelinguysen Rd., Piscataway, NJ 08854-8019 falk@mat.rutgers.edu
More informationRobotic manipulation project
Robotic manipulation project Bin Nguyen December 5, 2006 Abstract Tis is te draft report for Robotic Manipulation s class project. Te cosen project aims to understand and implement Kevin Egan s non-convex
More informationMANY scientific and engineering problems can be
A Domain Decomposition Metod using Elliptical Arc Artificial Boundary for Exterior Problems Yajun Cen, and Qikui Du Abstract In tis paper, a Diriclet-Neumann alternating metod using elliptical arc artificial
More informationComputing eigenvalues and eigenfunctions of Schrödinger equations using a model reduction approach
Computing eigenvalues and eigenfunctions of Scrödinger equations using a model reduction approac Suangping Li 1, Ziwen Zang 2 1 Program in Applied and Computational Matematics, Princeton University, New
More informationIsogeometric analysis and proper orthogonal decomposition for the acoustic wave equation
MATHICSE Matematics Institute of Computational Science and Engineering Scool of Basic Sciences - Section of Matematics MATHICSE Tecnical Report Nr. 02.206 January 206 Isogeometric analysis and proper ortogonal
More informationMixed Finite Element Methods for Incompressible Flow: Stationary Stokes Equations
Mixed Finite Element Metods for Incompressible Flow: Stationary Stoes Equations Ziqiang Cai, Carles Tong, 2 Panayot S. Vassilevsi, 2 Cunbo Wang Department of Matematics, Purdue University, West Lafayette,
More informationNotes on wavefunctions II: momentum wavefunctions
Notes on wavefunctions II: momentum wavefunctions and uncertainty Te state of a particle at any time is described by a wavefunction ψ(x). Tese wavefunction must cange wit time, since we know tat particles
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More informationJournal of Computational and Applied Mathematics
Journal of Computational and Applied Matematics 94 (6) 75 96 Contents lists available at ScienceDirect Journal of Computational and Applied Matematics journal omepage: www.elsevier.com/locate/cam Smootness-Increasing
More informationFEM solution of the ψ-ω equations with explicit viscous diffusion 1
FEM solution of te ψ-ω equations wit explicit viscous diffusion J.-L. Guermond and L. Quartapelle 3 Abstract. Tis paper describes a variational formulation for solving te D time-dependent incompressible
More informationKey words. Navier-Stokes, Unconditional stability, IMEX methods, second order convergence, Crank-Nicolson, BDF2
AN OPTIMALLY ACCURATE DISCRETE REULARIZATION FOR SECOND ORDER TIMESTEPPIN METHODS FOR NAVIER-STOKES EQUATIONS NAN JIAN, MUHAMMAD MOHEBUJJAMAN, LEO. REBHOLZ, AND CATALIN TRENCHEA Abstract. We propose a
More informationKey words. Finite element method; convection-diffusion-reaction; nonnegativity; boundedness
PRESERVING NONNEGATIVITY OF AN AFFINE FINITE ELEMENT APPROXIMATION FOR A CONVECTION-DIFFUSION-REACTION PROBLEM JAVIER RUIZ-RAMÍREZ Abstract. An affine finite element sceme approximation of a time dependent
More information