FINITE ELEMENT DUAL SINGULAR FUNCTION METHODS FOR HELMHOLTZ AND HEAT EQUATIONS
|
|
- Sherman Davis
- 5 years ago
- Views:
Transcription
1 J. KSIAM Vol.22, No.2, , 2018 ttp://dx.doi.org/ /jksiam FINITE ELEMENT DUAL SINGULAR FUNCTION METHODS FOR HELMHOLTZ AND HEAT EQUATIONS DEOK-KYU JANG AND JAE-HONG PYO DEPARTMENT OF MATHEMATICS, KANGWON NATIONAL UNIVERSITY, CHUNCHEON 24341, KOREA address: ABSTRACT. Te dual singular function metod(dsfm is a numerical algoritm to get optimal solution including corner singularities for Poisson and Helmoltz equations. In tis paper, we apply DSFM to solve eat equation wic is a time dependent problem. Since te DSFM for eat equation is based on DSFM for Helmoltz equation, it also need to use Serman- Morrison formula. Tis formula requires linear solver n + 1 times for elliptic problems on a domain including n reentrant corners. However, te DSFM for eat equation needs to pay only linear solver once per eac time iteration to standard numerical metod and perform optimal numerical accuracy for corner singularity problems. Because te Serman-Morrison formula is rater complicated to apply computation, we introduce a simplified formula by reanalyzing te Serman-Morrison metod. 1. INTRODUCTION In order to deal wit corner singularity solution, let te computational domain Ω be an open and bounded concave polygon in R 2. We assume tat Ω as one reentrant corner for simple explanation. It can be readily applied on multiple reentrant corner domain and is performed in 4. Te goal of tis paper is to construct te dual singular function metod to solve eat equation: u t µ u = f in Ω, (1.1 u = g on Ω, wit given initial value u(t = 0, x = (x in Ω and termal diffusivity µ > 0. Altoug te given function f is very smoot, te solution of (1.1 ave singular beavior near te reentrant corner. Tis corner singularity make lose accuracy of numerical solution trougout te wole domain. To overcome tis difficulty, tere are several primary ways. One is locally adaptive mes refinement or moving mes tecnique. Oter ways are so called Singular Function Metod (SFM (see, e.g., [9] and Dual Singular Function Metod (DSFM constructed in Received by te editors Marc ; Accepted June ; Publised online June Matematics Subject Classification. 65N12, 65M12, 65M60, 35A20. Key words and prases. Dual singular function metod, corner singularity, eat equation. Corresponding autor. 101
2 102 D.-K. JANG AND J.-H. PYO [1, 2, 7, 8]. DSFM for Poisson equation as been introduced in [3, 4] and ten it is extended to solve Helmoltz equations in [10]: µ u + ku = f in Ω, u = g on Ω. Because it is known in [5] tat te singular function of eat equation (1.1 can be written by a linear combination of tat of Helmoltz equation (1.2, we can construct DSFM to solve eat equation (1.1 by combining a time discrete sceme and results in [10]. DSFM as to use Serman-Morrison formula wic requests linear solver 2 times for any elliptic problem, but DSFM for eat equation does not need linear solvers per time iteration cost to standard numerical algoritm. Because te Serman-Morrison formula is rater complicated to apply real computation, we introduce a simplified version by reanalyzing te Serman-Morrison formula. Tis paper organized as follows. In section 2, we summarize te results of DSFM for Helmoltz equation in [10], and we establis te relationsip between regular part of te solution and discrete singular functions by analyzing Serman-Morrison formulation (see Teorem 2.4. We will construct DSFM to solve eat equation (1.1 in section 3 and present several numerical results in section DUAL SINGULAR FUNCTION METHOD TO SOLVE HELMHOLTZ EQUATION (1.2 We denote Ω = Γ in Γ out, were Γ in is part of boundary including reentrant corner and Γ out = Ω Γ in. Let ω be te internal angle of Γ in satisfying π < ω < 2π. Te singular and te dual singular functions of Poisson equation are summarized in [3, 4] as were s P (r, θ := r π/ω Θ(θ and s P d (r, θ := r π/ω Θ(θ, (2.1 Θ(θ = c cos(αθ + d sin(αθ. (2.2 We note tat c and d will be determined to old te boundary conditions of u near te corner on Γ in. We introduce te following lemma in [3]. Lemma 2.1. Te functions s P and s P d is armonic functions. Te singular function s H H 1 (Ω and te dual singular function s Hd L 2 (Ω of Helmoltz equation (1.2 ave to satisfy µ s H + ks H = 0 and µ s Hd + ks Hd = 0. (2.3 Te functions are written by power series s H (r, θ = r (1 π/ω + and s Hd (r, θ = r π/ω (1 + k i ν i 4 i i! i m=1 (m + αr2i k i ν i 4 i i! i m=1 (m αr2i Θ(θ (2.4 Θ(θ, (2.5
3 A FEM-DSFM FOR HELMHOLTZ AND HEAT EQUATIONS 103 were Θ(θ is defined in (2.2. We note ere s Hd / H 1 (Ω. Because r β is in H 2 (Ω if β 1, te singular parts of (2.4 and (2.5 are only S(r, θ := r π/ω Θ(θ and S d (r, θ := r π/ω Θ(θ. (2.6 So we can define singular and dual singular functions as in (2.1. We note ere tat S(r, θ and S d (r, θ are not satisfying (2.3 and it is an obstacle to construct DSFM for Helmoltz equation. We need to use a cut-off function to be 0 boundary condition, because S(r, θ and S d (r, θ of (2.6 is not 0 on Γ out. To do tis, we use notation B(r 1 ; r 2 = {(r, θ : r 1 < r < r 2 and 0 < θ < ω} Ω. We now define smoot cut-off function η ρ as 1 in B(0; 1 2 ρ, η ρ (r = p(r in B( 1 2ρ; ρ, 0 in Ω\B(0; ρ, (2.7 wit p(r is a very smoot function and ρ is a parameter determining te range of p(x. Ten te solution u of (1.2 can be rewritten by u = w + αη ρ S, (2.8 were α is called te stress intensity factor and w H 2 (Ω is regular part of solution. Te main idea of DSFM is to find α and w instead of computing u / H 2 (Ω. We start to construct DSFM wit rewriting te Helmoltz equation by inserting (2.8 into (1.2 µ w + kw + α( µ η ρ S + kη ρ S = f in Ω, w = g on Ω. (2.9 Since η ρ equals 0 identically in Γ out and S = 0 on Γ in, te function η ρ u S becomes 0 on Ω. So u and w ave same boundary values. To derive DSFM for eat equation, we test (2.9 by v H 1 0 to get w, v + kw, v + α µ η ρ S + ks, v = f, v. (2.10 Te variational formulation (2.10 as 2 unknowns (w, α, so we need one more equation and te equation is a single equation, because α is a real number. linearly independent to above system. not disappeared α. computable near te singularity corner. So we test wit te dual singular function S d (r, θ / H 1 (Ω to obtain µ w + kw, η 2ρ S d + α µ η ρ S + kη ρ S, η 2ρ S d = f, η 2ρ S d. (2.11 We ave a crucial lemma wic is proved in [10]. Lemma 2.2. Dual Singular function S d as te following properties.
4 104 D.-K. JANG AND J.-H. PYO 1. S d L 2 and S d / H ϕ, η 2ρ S d = ϕ, η 2ρ S d, for all ϕ H 1 0. In conjunction wit Lemma 2.2, (2.11 can be re written by α = 1 (β f w, µ η 2ρ S d + ks d, (2.12 were = µ η ρ S + ks, η 2ρ S d and β f = f, η 2ρ S d. DSFM is an algoritm to solve te system (2.11 and (2.12. In order to construct fully discrete algoritm of te system, we introduce finite element space V := {v H 1 0 (Ω : v K P(K K T}, were P(K is a polynomial function space degree p. In ligt of inserting α in (2.12 into (2.11, te discrete system becomes w V and α R satisfying, for all v V, w, v 1 w, µ η 2ρ S d + kη 2ρ S d µ η ρ S + kη ρ S, v Te matrix form of (2.13 is + kw, v = f, v β f µ η ρ S + kη ρ S, v. (2.13 (A + ab T w = F, (2.14 were te matrix A is te discrete linear operator of te Hemoltz equation. Also a, b T and 1 F in are generated by µ η ρ S + kη ρ S, v, w, µ η 2ρ S d + kη 2ρ S d and f, v β f µ η ρ S + kη ρ S, v, respectively. Te linear system (2.14 can be solved by te Serman- Morrison formulation in [6]: (A + ab T 1 = A 1 A 1 ab T A 1 a + b T A 1 a. (2.15 Te well posedness of (2.13 as been proved and obtain te following finite element approximation in [10]. Teorem 2.3. Let (w, α be te solution of (2.9 and w H 2 (Ω H 1 0 (Ω. And let (w, α be te solution of te system (2.12 and (2.13. Ten tere exists a positive constant 0 suc tat 0, for all, satisfying w w 1 C f 0, w w 0 C 1+α f 0, α α C 1+α f 0. Te Serman-Morrison formulation (2.15 is crucial to solve (2.14, but (2.15 requires linear solver 2 times and is rater complicated to apply to real computation, especially multiple singularities problems and time evolution problems like eat equation. So we reanalyze te Serman-Morrison formulation (2.15 wit solutions u and z of u, v + ku, v = f, v, v V, (2.16
5 and A FEM-DSFM FOR HELMHOLTZ AND HEAT EQUATIONS 105 z, v + kz, v = µ η ρ S + kη ρ S, v, v V, (2.17 receptively. Ten we arrive at te following result. Teorem 2.4. If we define ζ v = v, µ η 2ρ S d + kη 2ρ S d, for any function v V, ten w wic is te solution of (2.13 is ( βf ζ u w = u z. (2.18 ζ z Proof. In ligt of definition of ζ w, (2.13 can be rewritten by, for all v V, w, v + kw, v ζ w µ η ρ S + kη ρ S, v = f, v β f µ η ρ S + kη ρ S, v. Since b T A 1 F R, te Serman-Morrison formulation (2.15 yields From (2.16 and (2.17, w = A 1 F A 1 ab T A 1 a + b T A 1 a F = A 1 F bt A 1 F 1 + b T A 1 a A 1 a. A 1 F = u β f z. Also, A 1 a = z and b T v = ζ v /, for v V implies tat Terefore, (2.19 becomes b T A 1 a = ζ z / and b T A 1 F = ζ u w = u β f z ζ u ( βf = u te proof is complete. + β f ζ βs 2 z z 1 ζz β sζ u β f ζ z z ( ζ z ( βf ( ζ z ζ u + β f ζ z = u ( ζ z We so construct te following DSFM using (2.4. z = u + β f βs 2 ζ z. ( βf ζ u ζ z z. (2.19 Algoritm 1 (FE-DSFM for Helmoltz equation. Let w be discrete smoot part and let α be discrete stress intensity factors in (2.13.
6 106 D.-K. JANG AND J.-H. PYO Step 1: Find (u, z as te solution of (2.16 and (2.17. Step 2: Compute w by (2.18. Step 3: Compute α by second part of (2.13. We now apply Algoritm 1 to more complicated problem including several reentrant corners. Te solution u can be expressed by n u = w + α i η ρi S i, and te Helmoltz equation becomes n µ w + kw + α i ( µ η ρi S i + kη ρi S i = f in Ω, w = g on Ω. (2.20 We can coose cut off functions η ρi and η 2ρj to ave mutually disjoint support to old µ η ρi S i + ks i, η 2ρj S dj = 0 if i j. In ligt of testing (2.20 by η 2ρi S di, we readily obtain α i = 1 (β fi w, µ η 2ρi S di + ks di. (2.21 i Te finite element weak form of system (2.20 and (2.21 is to find w V and α i R, i = 1 n, satisfying, for all v V, n w, v + kw, v + α i µ η ρi S i + kη ρi S i, v = f, v, α i = 1 (2.22 (β fi w, µ η 2ρi S di + kη 2ρi S di, i and its matrix form becomes n (A + a i b T i w = F. Let ζ i,v = v, µ η 2ρi S di + ks di and z i be te solution of z i, v + kz i, v = µ η ρi S i + kη ρi S i, v, v V. (2.23 In conjunction wit Teorem 2.4, w would be n ( βfi ζ u w = u i ζ zi z i. (2.24 Finally we arrive at DSFM to solve Helmoltz equations (1.2 on domain Ω including several reentrant corners. Algoritm 2 (Multiple singularities version of Algoritm 1. Let w be discrete smoot part and let α i be discrete stress intensity factors of α i Step 1: Find u and z i, i = 1... n, as te solution of (2.16 and (2.23.
7 A FEM-DSFM FOR HELMHOLTZ AND HEAT EQUATIONS 107 Step 2: Compute w by (2.24. Step 3: Compute α by second part of ( DUAL SINGULAR FUNCTION METHOD FOR THE HEAT EQUATION In tis section, we extend Algoritm 1 to solve te eat equation (1.1 using te backward Euler time discrete formula (BDF1 and te second order backward Euler time discrete formula (BDF2. We first consider BDF1 for eat equation: u u n µ u = f. (3.1 τ We already know in [5] tat te form of singular solution is same to tat of Helmoltz equation. It means tat te solution of (1.1 can be expressed by and (3.1 can be rewritten by u(x, t = w(x, t + α(tη ρ S, (3.2 (w µτ w + α (η ρ S µτ η ρ S = τf + w n + α n η ρ S. (3.3 Since (3.3 and (2.9 ave similar expression, te same strategy to Algoritm 1 leads us te solver of (3.3: find w V and α R satisfying, for all v V, were µτ w, v + w, v + α η ρ S µτ η ρ S, v = τ f, v + w n, v + α n η ρs, v, α = 1 ( τβ β f + w n + α n η ρ S, η 2ρ S d ζ w, s = η ρ S µτ η ρ S, η 2ρ S d, β f = f(t, η 2ρ S d, (3.4 ζ w = w, η 2ρ S d µτ η 2ρ S d. Te matrix form of te system becomes (2.14 and it can be solved by (2.15. To construct algoritm for eat equations, we denote z be te solution of µτ z, v + z, v = µτ η ρ S + η ρ S, v, v V. (3.5 Ten we arrive at DSFM wit BDF1 sceme to solve eat equation. Algoritm 3 ( FE-DSFM wit backward Euler for Heat equation. Let z be te solution of (3.5 and set û 0 = u0. Repeat for 1 n N: Step 1: Find u as te solution of µτ u, v + u, v = τ f, v + û n, v, v V. Step 2: Compute w w = u ( τβf + û n, η 2ρS d ζ u z. ζ z
8 108 D.-K. JANG AND J.-H. PYO Step 3: Compute α by second part of (3.4. Step 4: Update û = w + α η ρ S. Remark 1. Te Serman-Morrison formulation (2.15 requires to solve linear system 2 times for Poisson and Helmoltz equations, but Algoritm 3 need to apply linear solver only 1 at eac iteration, so tis algoritm use same linear solvers per eac time iteration to standard metod. And Algoritm 3 as optimal accuracy beavior for corner singularities problem, in contrast standard numerical metod. We now construct DSFM for eat equation wit backward Euler formula (BDF2 time discretization: 3u 4u n + u n 1 µ u = f. 2τ By te same manner wit te case BDF1, (3.2 leads us 3w 2µτ w + α (3η ρ S 2µτ η ρ S and testing (3.6 wit v V and wit η 2ρi S di yield were = 2τf + 4(w n + α n η ρ S (w n 1 + α n 1 η ρ S, 2µτ w, v + 3 w, v + α 3η ρ S 2µτ η ρ S, v = 2τ f, v + 4( w n, v + α n η ρs, v ( w n 1, v + α n 1 η ρ S, v, v V, α = 1 (F β ζ w, s F = 2τ f, η 2ρ S d + 4( w n, η 2ρ S d + α n η ρ S, η 2ρ S d ( w n 1, η 2ρ S d + α n 1 η ρ S, η 2ρ S d, ζ w = w, 3η 2ρ S d 2µτ η 2ρ S d and βs = 3η ρ S 2µτ η ρ S, η 2ρ S d. Since tis system also constructs a linear system of te form we can find w by Teorem 2.4, tat is were z is te solution of w = u (A + ab T w = F, ( F ζ u ζ z (3.6 (3.7 z, (3.8 2µτ z, v + 3 z, v = 2µτ η ρ S + 3η ρ S, v, v V.
9 If we denote û n = wn +αn η ρs and û n 1 of, for all v V, A FEM-DSFM FOR HELMHOLTZ AND HEAT EQUATIONS 109 2µτ u, v + 3 u, v Finally, we arrive at BDF2 DSFM to solve eat equation (1.1. = w n 1 +α n 1 η ρ S, ten u becomes te solution = 2τ f, v + 4 û n, v û n 1, v. (3.9 Algoritm 4 ( FE-DSFM wit BDF2 for Heat equation. Let w be discrete smoot part and α be discrete stress intensity factors of α(t of n + 1 step. Initially given u 0, set w 1 and α 1 using te backward FE-DSFM. Repeat for 1 n N Step 1: Find u as te solution of (3.9. Step 2: Compute w by (3.8. Step 3: Compute α by second part of ( NUMERICAL TEST In tis section, we document te computational performance of eac algoritm witin a polygonal domain wit reentrant corners. We use cut-off function in (2.7 wit p(r = 15 { 8 ( 4r ρ ( 4r 3 3 ρ 3 1 ( 4r 5 } 5 ρ 3. All computations are performed wit te conforming P 1 finite element and 6 points quadrature rule. Example 4.1. We consider Helmoltz equations (1.2 on te Γ sape computational domain ([ 1, 1] [ 1, 1] ([0, 1] [ 1, 0]. We coose te smoot part of te solution as { sin(2πx(1/2y 2 + y(y 2 1, (y < 0 w = and set te exact solution sin(2πx( 1/2y 2 + y(y 2 1, (y 0, u = w + αη ρ S wit α = 3.0 and ρ = 0.4. Also, te singular function S is given by S(r, θ = r π/ω sin( π ω θ wit ω = 3 2π. Te forcing term f is determined to become k = 1.0 and ν = 1.0. Table 1 is te error decay of te standard finite element metod for Example 4.1. Te losing convergence orders is natural beavior because of u / H 2 (Ω. Table 2 is te error decay of te DSFM wit Algoritm 1 and displays optimal error decay in contrast wit results in Table 1 Te next example is a multiple corner singularities case.
10 110 D.-K. JANG AND J.-H. PYO TABLE 1. Error table for Example 4.1 wit Standard FEM. u u 0 Iu u u u 1 Errors Order Errors Order Errors Order 1/ e e e-00 1/ e e e / e e e / e e e / e e e / e e e TABLE 2. Error table for Example 4.1 wit Algoritm 1. w w 0 Iw w w w 1 α α Errors Order Errors Order Errors Order Errors Order 1/ e e e e-03 1/ e e e e / e e e e / e e e e / e e e e / e e e e Example 4.2. Let te computational domain be ([ 2, 2] [ 2, 2] ([ 1, 1] [ 1.1] including 4 reentrant corners. Let smoot part of te solution be given by { sin(πx(1/2y 2 + y(y + 1(y + 2, (y < 0 w = sin(πx( 1/2y 2 + y(y 1(y 2, (y 0. In tis experiment, we coose k = 1.0, ν = 1.0 and set te forcing term f to be te exact solution u = w + α 1 η ρ1 S 1 + α 2 η ρ2 S 2 + α 3 η ρ3 S 3 + α 4 η ρ4 S 4, were ρ 1 = ρ 2 = ρ 3 = ρ 4 = 0.4 and stress intensity factors α 1 = 1.0, α 2 = 2.0, α 3 = 3.0 and α 4 = 4.0. Table 3 is te mes analysis of Algoritm 2 wic is based on (2.24. Tis results also display optimal error decay rates. We now perform numerical simulations for te eat equations. Example 4.3. We consider eat equation wit Γ sape computational domain ([ 1, 1] [ 1, 1] ([0, 1] [ 1, 0]. We use t as time variable and we perform Algoritms 3 and
11 A FEM-DSFM FOR HELMHOLTZ AND HEAT EQUATIONS FIGURE 1. Domain, mes and η ρi S i, (i = 1 4 for Example 4.2 TABLE 3. Error table for Example 4.2 wit Algoritm 2. w w 0 Iw w w w 1 4 α i α i Errors Order Errors Order Errors Order Errors Order 1/ e e e e-02 1/ e e e e / e e e e / e e e e / e e e e / e e e e wit smoot part of solution { sin(t sin(2πx(1/2y 2 + y(y 2 1, (y < 0 w = sin(t sin(2πx( 1/2y 2 + y(y 2 1, (y 0, and exact solution u = w + exp(tη ρ S in [0, 1] Ω, were ρ = 0.4 and forcing term is also determined by ν = 1.0. In tis example, te stress intensity factor depends on variable t and Tables 4 6 are mes analysis at t = 1. Since backward Euler metod is only 1st order sceme for time, we set τ = 16 2 to ceck te optimal convergence rates and Table 4 sows te optimal results of Algoritm 3. We apply te BDF2 time discretization sceme, Algoritm 4, for te same example. Since it is second order metod, we set τ = 1 2, and so tis computational cost is muc more effective tan Algoritm 3. Table 5 is results of te standard finite element metod wit BDF2 time discretization for Example 4.3. Since u / L 2 ([0, 1] : H 2, te errors do not decay wit optimal order.
12 112 D.-K. JANG AND J.-H. PYO TABLE 4. Error table for Example 4.3 wit Algoritm 3. w w 0 Iw w w w 1 α α τ Errors Order Errors Order Errors Order Errors Order 1/8 1/ e e e e-03 1/16 1/ e e e e /32 1/ e e e e /64 1/ e e e e /128 1/ e e e e /256 1/ e e e e TABLE 5. Error table for Example 4.3 wit Standard FEM wit BDF2 time discretization. u u 0 Iu u u u 1 τ Errors Order Errors Order Errors Order 1/8 1/ e e e-00 1/16 1/ e e e /32 1/ e e e /64 1/ e e e /128 1/ e e e /256 1/ e e e TABLE 6. Error table for Example 4.3 wit Algoritm 4. w w 0 Iw w w w 1 α α τ Errors Order Errors Order Errors Order Errors Order 1/8 1/ e e e e-03 1/16 1/ e e e e /32 1/ e e e e /64 1/ e e e e /128 1/ e e e e /256 1/ e e e e Table 6 is error table for Algoritm 4 for Example 4.3 and te errors converge to 0 wit optimal rate. So we can conclude tat DSFM for eat equation as optimal convergence rate for corner singularity problems. Finally, Figure 2 is error evolutions of smoot part of u on te L 2 space and errors of stress intensity factors. It displays smoot error propagation on time.
13 A FEM-DSFM FOR HELMHOLTZ AND HEAT EQUATIONS 113 Error 1.5 x 10 3 w w wit =1/32 τ =1/64 0 α α wit =1/32 τ =1/64 w w wit =1/64 τ =1/ α α wit =1/64 τ =1/ Time FIGURE 2. Error evolution of Example 4.3 wit Algoritm 4. ACKNOWLEDGMENT Tis study was supported by 2015 Researc Grant from Kangwon National University (No REFERENCES [1] H. Blum and M. Dobrowolski, On finite element metods for elliptic equations on domains wit corners, Computing 28 (1982, [2] M. Bourlard, M. Dauge, M.-S. Lubuma, and S. Nicaise, Coeficients of te singularities for elliptic boundary value problems on domains wit conical points III. Finite element metods on polygonal domains, SIAM Numer. Anal. 29 (1992, [3] Z. Cai and S. Kim, A finite element metod using singular functions for te Poisson equation: Corner singularities, SIAM J. Numer. Anal. 39 (2001, [4] Z. Cai, S. C. Kim and B. C. Sin, Solution metods for te Poisson equation: corner singularities, SIAM J. SCI. COMPut. 23 (2001, [5] H.J. Coi and J.R. Kweon Te fourier finite element metod for te corner singularity expansion of te eat equation Computers Mat. Applic. 69 (2015, [6] Gene H. Colub and Carles F. Van Loan Matrix Computations, Te Jon Hopkins University Press (1996. [7] M. Djaoua, Equations Intégrales pour un Probleme Singulier dans le Plan, Tese de Troisieme Cycle, Universite Pierre et Marie Curie, Paris, [8] M. Dobrowolski, Numerical approximation of elliptic interface and corner problems, Habilitation-Scrift, Bonn, [9] G. J. Fix, S. Gulati and G. I. Wakoff, On te use of singular functions wit finite elements approximations, J. Comput. Py. 13 (1973, [10] S.C. Kim, J.S. Lee and J.H. Pyo A finite element metod using singular functions for Helmoltz equations: Part I, J. KSIAM 12 (2008,
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
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 informationA Finite Element Method Using Singular Functions: Interface Problems
A Finite Element Method Using Singular Functions: Interface Problems Seokchan Kim Zhiqiang Cai Jae-Hong Pyo Sooryoun Kong Abstract The solution of the interface problem is only in H 1+α (Ω) with α > 0
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 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 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 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 informationDifferent Approaches to a Posteriori Error Analysis of the Discontinuous Galerkin Method
WDS'10 Proceedings of Contributed Papers, Part I, 151 156, 2010. ISBN 978-80-7378-139-2 MATFYZPRESS Different Approaces to a Posteriori Error Analysis of te Discontinuous Galerkin Metod I. Šebestová Carles
More informationImplicit-explicit variational integration of highly oscillatory problems
Implicit-explicit variational integration of igly oscillatory problems Ari Stern Structured Integrators Worksop April 9, 9 Stern, A., and E. Grinspun. Multiscale Model. Simul., to appear. arxiv:88.39 [mat.na].
More information1. 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 informationOn convexity of polynomial paths and generalized majorizations
On convexity of polynomial pats and generalized majorizations Marija Dodig Centro de Estruturas Lineares e Combinatórias, CELC, Universidade de Lisboa, Av. Prof. Gama Pinto 2, 1649-003 Lisboa, Portugal
More informationDecay of solutions of wave equations with memory
Proceedings of te 14t International Conference on Computational and Matematical Metods in Science and Engineering, CMMSE 14 3 7July, 14. Decay of solutions of wave equations wit memory J. A. Ferreira 1,
More informationContinuity and Differentiability of the Trigonometric Functions
[Te basis for te following work will be te definition of te trigonometric functions as ratios of te sides of a triangle inscribed in a circle; in particular, te sine of an angle will be defined to be te
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 34, pp. 14-19, 2008. Copyrigt 2008,. ISSN 1068-9613. ETNA A NOTE ON NUMERICALLY CONSISTENT INITIAL VALUES FOR HIGH INDEX DIFFERENTIAL-ALGEBRAIC EQUATIONS
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 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 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 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 informationStability properties of a family of chock capturing methods for hyperbolic conservation laws
Proceedings of te 3rd IASME/WSEAS Int. Conf. on FLUID DYNAMICS & AERODYNAMICS, Corfu, Greece, August 0-, 005 (pp48-5) Stability properties of a family of cock capturing metods for yperbolic conservation
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 informationParametric Spline Method for Solving Bratu s Problem
ISSN 749-3889 print, 749-3897 online International Journal of Nonlinear Science Vol4202 No,pp3-0 Parametric Spline Metod for Solving Bratu s Problem M Zarebnia, Z Sarvari 2,2 Department of Matematics,
More informationOverlapping domain decomposition methods for elliptic quasi-variational inequalities related to impulse control problem with mixed boundary conditions
Proc. Indian Acad. Sci. (Mat. Sci.) Vol. 121, No. 4, November 2011, pp. 481 493. c Indian Academy of Sciences Overlapping domain decomposition metods for elliptic quasi-variational inequalities related
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 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 informationBlanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS
Opuscula Matematica Vol. 26 No. 3 26 Blanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS Abstract. In tis work a new numerical metod is constructed for time-integrating
More informationLecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator.
Lecture XVII Abstract We introduce te concept of directional derivative of a scalar function and discuss its relation wit te gradient operator. Directional derivative and gradient Te directional derivative
More informationA h u h = f h. 4.1 The CoarseGrid SystemandtheResidual Equation
Capter Grid Transfer Remark. Contents of tis capter. Consider a grid wit grid size and te corresponding linear system of equations A u = f. Te summary given in Section 3. leads to te idea tat tere migt
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 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 informationThe Laplace equation, cylindrically or spherically symmetric case
Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,
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 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 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 informationCombining functions: algebraic methods
Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)
More informationarxiv: v1 [math.ap] 4 Aug 2017
New Two Step Laplace Adam-Basfort Metod for Integer an Non integer Order Partial Differential Equations arxiv:178.1417v1 [mat.ap] 4 Aug 17 Abstract Rodrigue Gnitcogna*, Abdon Atangana** *Department of
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 informationBOUNDARY ELEMENT METHODS FOR POTENTIAL PROBLEMS WITH NONLINEAR BOUNDARY CONDITIONS
MATHEMATICS OF COMPUTATION Volume 70, Number 235, Pages 1031 1042 S 0025-5718(00)01266-7 Article electronically publised on June 12, 2000 BOUNDARY ELEMENT METHODS FOR POTENTIAL PROBLEMS WITH NONLINEAR
More informationNew Fourth Order Quartic Spline Method for Solving Second Order Boundary Value Problems
MATEMATIKA, 2015, Volume 31, Number 2, 149 157 c UTM Centre for Industrial Applied Matematics New Fourt Order Quartic Spline Metod for Solving Second Order Boundary Value Problems 1 Osama Ala yed, 2 Te
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 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 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 informationError estimates for a semi-implicit fully discrete finite element scheme for the mean curvature flow of graphs
Interfaces and Free Boundaries 2, 2000 34 359 Error estimates for a semi-implicit fully discrete finite element sceme for te mean curvature flow of graps KLAUS DECKELNICK Scool of Matematical Sciences,
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 informationDELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science
DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Matematics and Computer Science. ANSWERS OF THE TEST NUMERICAL METHODS FOR DIFFERENTIAL EQUATIONS (WI3097 TU) Tuesday January 9 008, 9:00-:00
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 informationParallel algorithm for convection diffusion system based on least squares procedure
Zang et al. SpringerPlus (26 5:69 DOI.86/s464-6-3333-8 RESEARCH Parallel algoritm for convection diffusion system based on least squares procedure Open Access Jiansong Zang *, Hui Guo 2, Hongfei Fu 3 and
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More informationCONVERGENCE OF AN IMPLICIT FINITE ELEMENT METHOD FOR THE LANDAU-LIFSHITZ-GILBERT EQUATION
CONVERGENCE OF AN IMPLICIT FINITE ELEMENT METHOD FOR THE LANDAU-LIFSHITZ-GILBERT EQUATION SÖREN BARTELS AND ANDREAS PROHL Abstract. Te Landau-Lifsitz-Gilbert equation describes dynamics of ferromagnetism,
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 informationVolume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households
Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of
More informationA proof in the finite-difference spirit of the superconvergence of the gradient for the Shortley-Weller method.
A proof in te finite-difference spirit of te superconvergence of te gradient for te Sortley-Weller metod. L. Weynans 1 1 Team Mempis, INRIA Bordeaux-Sud-Ouest & CNRS UMR 551, Université de Bordeaux, France
More informationPOLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY
APPLICATIONES MATHEMATICAE 36, (29), pp. 2 Zbigniew Ciesielski (Sopot) Ryszard Zieliński (Warszawa) POLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY Abstract. Dvoretzky
More informationarxiv: v3 [math.na] 15 Dec 2009
THE NAVIER-STOKES-VOIGHT MODEL FOR IMAGE INPAINTING M.A. EBRAHIMI, MICHAEL HOLST, AND EVELYN LUNASIN arxiv:91.4548v3 [mat.na] 15 Dec 9 ABSTRACT. In tis paper we investigate te use of te D Navier-Stokes-Voigt
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 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 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 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 informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
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 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 informationch (for some fixed positive number c) reaching c
GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan
More informationBOUNDARY REGULARITY FOR SOLUTIONS TO THE LINEARIZED MONGE-AMPÈRE EQUATIONS
BOUNDARY REGULARITY FOR SOLUTIONS TO THE LINEARIZED MONGE-AMPÈRE EQUATIONS N. Q. LE AND O. SAVIN Abstract. We obtain boundary Hölder gradient estimates and regularity for solutions to te linearized Monge-Ampère
More informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
More informationFinite Difference Method
Capter 8 Finite Difference Metod 81 2nd order linear pde in two variables General 2nd order linear pde in two variables is given in te following form: L[u] = Au xx +2Bu xy +Cu yy +Du x +Eu y +Fu = G According
More informationConvexity and Smoothness
Capter 4 Convexity and Smootness 4.1 Strict Convexity, Smootness, and Gateaux Differentiablity Definition 4.1.1. Let X be a Banac space wit a norm denoted by. A map f : X \{0} X \{0}, f f x is called a
More informationarxiv: v1 [math.dg] 4 Feb 2015
CENTROID OF TRIANGLES ASSOCIATED WITH A CURVE arxiv:1502.01205v1 [mat.dg] 4 Feb 2015 Dong-Soo Kim and Dong Seo Kim Abstract. Arcimedes sowed tat te area between a parabola and any cord AB on te parabola
More informationDepartment of Mathematical Sciences University of South Carolina Aiken Aiken, SC 29801
RESEARCH SUMMARY AND PERSPECTIVES KOFFI B. FADIMBA Department of Matematical Sciences University of Sout Carolina Aiken Aiken, SC 29801 Email: KoffiF@usca.edu 1. Introduction My researc program as focused
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 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 informationSmoothness of solutions with respect to multi-strip integral boundary conditions for nth order ordinary differential equations
396 Nonlinear Analysis: Modelling and Control, 2014, Vol. 19, No. 3, 396 412 ttp://dx.doi.org/10.15388/na.2014.3.6 Smootness of solutions wit respect to multi-strip integral boundary conditions for nt
More informationClick here to see an animation of the derivative
Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,
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 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 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 informationarxiv: v1 [math.na] 3 Nov 2011
arxiv:.983v [mat.na] 3 Nov 2 A Finite Difference Gost-cell Multigrid approac for Poisson Equation wit mixed Boundary Conditions in Arbitrary Domain Armando Coco, Giovanni Russo November 7, 2 Abstract In
More informationarxiv: v1 [math.na] 18 Jul 2015
ENERGY NORM ERROR ESTIMATES FOR FINITE ELEMENT DISCRETIZATION OF PARABOLIC PROBLEMS HERBERT EGGER arxiv:150705183v1 [matna] 18 Jul 2015 Department of Matematics, TU Darmstadt, Germany Abstract We consider
More informationThe Verlet Algorithm for Molecular Dynamics Simulations
Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical
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 informationNumerical Solution of One Dimensional Nonlinear Longitudinal Oscillations in a Class of Generalized Functions
Proc. of te 8t WSEAS Int. Conf. on Matematical Metods and Computational Tecniques in Electrical Engineering, Bucarest, October 16-17, 2006 219 Numerical Solution of One Dimensional Nonlinear Longitudinal
More informationLIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION
LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LAURA EVANS.. Introduction Not all differential equations can be explicitly solved for y. Tis can be problematic if we need to know te value of y
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 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 informationInfluence of the Stepsize on Hyers Ulam Stability of First-Order Homogeneous Linear Difference Equations
International Journal of Difference Equations ISSN 0973-6069, Volume 12, Number 2, pp. 281 302 (2017) ttp://campus.mst.edu/ijde Influence of te Stepsize on Hyers Ulam Stability of First-Order Homogeneous
More informationFINITE ELEMENT APPROXIMATIONS AND THE DIRICHLET PROBLEM FOR SURFACES OF PRESCRIBED MEAN CURVATURE
FINITE ELEMENT APPROXIMATIONS AND THE DIRICHLET PROBLEM FOR SURFACES OF PRESCRIBED MEAN CURVATURE GERHARD DZIUK AND JOHN E. HUTCHINSON Abstract. We give a finite element procedure for te Diriclet Problem
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 informationSome Applications of Fractional Step Runge-Kutta Methods
Some Applications of Fractional Step Runge-Kutta Metods JORGE, J.C., PORTERO, L. Dpto. de Matematica e Informatica Universidad Publica de Navarra Campus Arrosadia, s/n 3006 Pamplona Navarra SPAIN Abstract:
More informationA First-Order System Approach for Diffusion Equation. I. Second-Order Residual-Distribution Schemes
A First-Order System Approac for Diffusion Equation. I. Second-Order Residual-Distribution Scemes Hiroaki Nisikawa W. M. Keck Foundation Laboratory for Computational Fluid Dynamics, Department of Aerospace
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 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 information158 Calculus and Structures
58 Calculus and Structures CHAPTER PROPERTIES OF DERIVATIVES AND DIFFERENTIATION BY THE EASY WAY. Calculus and Structures 59 Copyrigt Capter PROPERTIES OF DERIVATIVES. INTRODUCTION In te last capter you
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 information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
More information5.74 Introductory Quantum Mechanics II
MIT OpenCourseWare ttp://ocw.mit.edu 5.74 Introductory Quantum Mecanics II Spring 9 For information about citing tese materials or our Terms of Use, visit: ttp://ocw.mit.edu/terms. Andrei Tokmakoff, MIT
More informationAPPROXIMATION OF CRYSTALLINE DENDRITE GROWTH IN TWO SPACE DIMENSIONS. Introduction
Acta Mat. Univ. Comenianae Vol. LXVII, 1(1998), pp. 57 68 57 APPROXIMATION OF CRYSTALLINE DENDRITE GROWTH IN TWO SPACE DIMENSIONS A. SCHMIDT Abstract. Te pase transition between solid and liquid in an
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 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 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 Difference-Analytical Method for Solving Laplace s Boundary Value Problem with Singularities
5 10 July 2004, Antalya, urkey Dynamical Systems and Applications, Proceedings, pp. 339 360 A Difference-Analytical Metod for Solving Laplace s Boundary Value Problem wit Singularities A. A. Dosiyev and
More informationImproved Rotated Finite Difference Method for Solving Fractional Elliptic Partial Differential Equations
American Scientific Researc Journal for Engineering, Tecnolog, and Sciences (ASRJETS) ISSN (Print) 33-44, ISSN (Online) 33-44 Global Societ of Scientific Researc and Researcers ttp://asrjetsjournal.org/
More information