ON THE CONVERGENCE OF A DUAL-PRIMAL SUBSTRUCTURING METHOD. January 2000
|
|
- Eustacia Matilda Holt
- 6 years ago
- Views:
Transcription
1 ON THE CONVERGENCE OF A DUAL-PRIMAL SUBSTRUCTURING METHOD JAN MANDEL AND RADEK TEZAUR January 2000 Abstract In the Dual-Primal FETI method, introduced by Farhat et al [5], the domain is decomposed into non-overlapping subdomains, but the degrees of freedom on crosspoints remain common to all subdomains adjacent to the crosspoint The continuity of the remaining degrees of freedom on subdomain interfaces is enforced by Lagrange multipliers and all degrees of freedom are eliminated The resulting dual problem is solved by preconditioned conjugate gradients We give an algebraic bound on the condition number, assuming only a single inequality in discrete norms, and use the algebraic bound to show that the condition number is bounded by C(1 + log 2 (H/h)) for both second and fourth order elliptic selfadjoint problems discretized by conforming finite elements, as well as for a wide class of finite elements for the Reissner-Mindlin plate model 1 Introduction This article is concerned with convergence bounds for an iterative method for the parallel solution of symmetric, positive definite systems of linear equations that arise from elliptic boundary value problems discretized by finite elements The original Finite Element Tearing and Interconnecting method (FETI) was proposed by Farhat and Roux [9] The FETI method consists of decomposing the domain into non-overlapping subdomains, enforcing that the corresponding degrees of freedom on subdomain interfaces coincide by Lagrange multipliers, and eliminating all degrees of freedom, leaving a dual system for the Lagrange multipliers The dual system was solved by preconditioned conjugate gradients with a diagonal preconditioner Evaluation of the dual operator involves the solution of independent Neumann problems in all subdomains, and of a small system of equations for the nullspace component Farhat, Mandel, and Roux [8] recognized that this system for the nullspace components plays the role of a coarse problem that facilitates global exchange of information between the subdomains, causing the condition to be bounded as the number of subdomains increases They also replaced the diagonal preconditioner by a block preconditioner with the solution of independent Dirichlet problems in each subdomain and observed numerically that this Dirichlet preconditioner results in a very slow growth of the condition number with subdomain size Mandel and Tezaur [13] proved that the condition number grows at most as log 3 (H/h), where H is subdomain size and h is element size, both in 2D and 3D Tezaur [19] proved that a method by Park, Justino, and Felippa [16] is equivalent to the method of [9] with a special choice of the constraint matrices, and proved the log 2 (H/h) bound for this variant For further comparison, see Rixen et al [18] Klawonn and Widlund [10] have used preconditioned conjugate residuals to solve a saddle problem keeping both the original degrees of freedom and the Lagrange multipliers, and obtained the asymptotic bound log 2 (H/h) using an extension of the theory of [13] The saddle point approach has the advantage that approximate solvers can be used for both the Neumann and the Dirichlet subdomain problems, at the cost This work was supported by NSF grant ECS and by ONR grant N Department of Mathematics, University of Colorado at Denver, Denver, CO , USA, and Department of Aerospace Engineering Sciences, University of Colorado at Boulder, Boulder, CO , USA jmandel@coloradoedu Center for Aerospace Structures and Department of Aerospace Engineering Sciences, University of Colorado at Boulder, Boulder, CO , USA rtezaur@coloradoedu 1
2 Report Documentation Page Form Approved OMB No Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information Send comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington VA Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to a penalty for failing to comply with a collection of information if it does not display a currently valid OMB control number 1 REPORT DATE JAN REPORT TYPE 3 DATES COVERED - 4 TITLE AND SUBTITLE ON THE CONVERGENCE OF A DUAL-PRIMAL SUBSTRUCTURING METHOD 5a CONTRACT NUMBER 5b GRANT NUMBER 5c PROGRAM ELEMENT NUMBER 6 AUTHOR(S) 5d PROJECT NUMBER 5e TASK NUMBER 5f WORK UNIT NUMBER 7 PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) Office of Naval Research,875 North Randolph Streeet,Arlington,VA, PERFORMING ORGANIZATION REPORT NUMBER 9 SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10 SPONSOR/MONITOR S ACRONYM(S) 12 DISTRIBUTION/AVAILABILITY STATEMENT Approved for public release; distribution unlimited 13 SUPPLEMENTARY NOTES 11 SPONSOR/MONITOR S REPORT NUMBER(S) 14 ABSTRACT In the Dual-Primal FETI method, introduced by Farhat et al [5], the domain is decomposed into non-overlapping subdomains, but the degrees of freedom on crosspoints remain common to all subdomains adjacent to the crosspoint The continuity of the remaining degrees of freedom on subdomain interfaces is enforced by Lagrange multipliers and all degrees of freedom are eliminated The resulting dual problem is solved by preconditioned conjugate gradients We give an algebraic bound on the condition number, assuming only a single inequality in discrete norms and use the algebraic bound to show that the condition number is bounded by C(1+log2(H=h)) for both second and fourth order elliptic selfadjoint problems discretized by conforming nite elements as well as for a wide class of nite elements for the Reissner-Mindlin plate model 15 SUBJECT TERMS 16 SECURITY CLASSIFICATION OF: 17 LIMITATION OF ABSTRACT a REPORT unclassified b ABSTRACT unclassified c THIS PAGE unclassified 18 NUMBER OF PAGES 12 19a NAME OF RESPONSIBLE PERSON Standard Form 298 (Rev 8-98) Prescribed by ANSI Std Z39-18
3 of of solving a larger indefinite problem instead of a small positive definite problem In [10], Klawonn and Widlund have proposed new preconditioners, proved uniform bounds for modifications of the method for the case of coefficient jumps, which include an earlier algorithm of Rixen and Farhat [17], and provided further theoretical insights The original method of [9, 8] does not converge well for plate and shell problems, and the existence and the form of the coarse space depend on the singularity of the subdomain matrices Therefore, Mandel, Tezaur, and Farhat [14] and Farhat, Mandel, and Chen [7] proposed to project the Lagrange multipliers in each iteration on an auxiliary space An auxiliary space chosen so that the corresponding primal solutions are continuous on crosspoints made it possible to prove that the condition number does not grow faster than log 3 (H/h) for plate problems [14], and fast convergence was observed for plate [7] as well as shell problems [4] This method is now called FETI-2 For related results for symmetric positive definite problems, see [6, 18] and references therein The subject of this paper is the Dual-Primal FETI method (FETI-DP), introduced by Farhat et al [5] This method enforces the continuity of the primal solution at crosspoints directly by the formulation of the dual problem: the degrees of freedom on a crosspoint remain common to all subdomains sharing the crosspoint and the continuity of the remaining degrees of freedom on the interfaces in enforced by Lagrange multipliers The degrees of freedom are then eliminated and the resulting dual problem for the Lagrange multipliers is solved by preconditioned conjugate gradients with a Dirichlet preconditioner Evaluating the dual operator involves the solution of independent subdomain problems with nonsingular matrices and of a coarse problem based on subdomain corners The advantage of this method is a simpler formulation than the methods of [14, 7], there is no need to solve problems with singular matrices, and the method was observed to be significantly faster in practice for 2D problems However, the design of a good 3D variant of the method is an open problem [5] In this paper, we prove that the condition number of the FETI-DP method with the Dirichlet preconditioner does not grow faster than log 2 (H/h) for both second order and fourth order problems in 2D By spectral equivalence, the result for fourth order problems extends to a large class of Reissner-Mindlin elements for plate bending as in [11, 12, 14] After few initial definitions, which are substantially different, the analysis is related to the analysis developed in [13, 14] Just as the formulation of the present method is simpler, the analysis is simpler and more elegant than in [13, 14] The paper is organized as follows The notation and assumptions are introduced in Section 2 In Section 3, we review the algorithm from [5] Section 4 gives the algebraic condition number estimates; these estimates apply to partitioning any symmetric positive definite system, not necessarily originating from partial differential equations Finally, in Section 5, we prove the polylogarithmic condition number bounds for two model problems of the second and fourth order 2 Domain partitioning, notation, and assumptions We are concerned with iterative solution of symmetric, positive definite linear algebraic systems Partitioning of the system is motivated as follows Let Ω be a domain in R 2 decomposed into N s non-overlapping subdomains Ω 1, Ω 2,, Ω Ns, and each of the subdomains be a union of some of the elements Let u s be the vector of degrees of freedom for the subdomain Ω s corresponding to a conforming finite element discretization of a second order elliptic problem or a fourth order plate bending problem defined on Ω Let K s and f s be the local stiffness matrix and the load vector associated with the subdomain Ω s We denote the edges of the subdomains by Γ st = Ω s Ω t Corners 2
4 are endpoints of edges The subdomain vectors are partitioned as u s u s i (21) = u s r, u s c where u s i are the values of the degrees of freedom in the subdomain interior, us c the values of the degrees of freedom at the corners of the subdomain, and u s r are the remaining values of the degrees of freedom, ie those located on the edges of the subdomains between the corners The subdomain matrices are partitioned accordingly, K s = Ks ii Kir s Kic s (22) Kri s Krr s Krc s Kci s Kcr s Kcc s We use the block notation u = u 1 u Ns and K = diag(k s ) = K 1 K Ns The block vectors u i, u c and u r of all internal, all corner, and all remainder degrees of freedom, respectively, are then defined similarly, u i = u 1 i u Ns i, u r = u 1 r u Ns r, and u c = u 1 c u Ns c Vectors of values of degrees of freedom on the whole Ω s and the corresponding block vectors will be written as [ ] v 1 vr,c s v s r,c [ ] = r vr vc s, v r,c =, and also v r,c = v c vr,c Ns Let the global to local map be a 0 1 block matrix L = L 1 L Ns That is, for a global vector of degrees of freedom u g, L s u g is the vector of corresponding degrees of freedom on Ω s The map L is introduced so that we can state the problem to be solved, (23) L T KLu g = L T f, independently of the solution algorithm Note that ImL is the space of all vectors u that are continuous across the subdomain interfaces and that L is of full column rank We assume that each matrix K s is symmetric positive semidefinite and that K s is positive definite on the subspace of vectors that are zero at subdomain corners, {u s u s c = 0} We also assume that the global stiffness matrix L T KL is positive definite, or, equivalently, that K is positive definite on ImL These assumptions are satisfied in the intended finite element applications 3
5 3 Formulation of the algorithm In this section, we review the algorithm proposed in [5] in a form suitable for our purposes The degrees of freedom from both sides of each edge Γ st should coincide, (31) u s r Γ st u t r Γ st = 0 In (31), each pair of subdomains {s, t} is taken only once, with the order (s, t) chosen arbitrarily We write the constraints (31) as B r u r = 0, B r = [B 1 r,, B Ns r ] Note that it follows immediately from the definition of B r that (32) and, for any edge Γ st, (33) B r B T r = 2I (B T r B r u r ) s Γ st = ±(u s r Γ st u t r Γ st) Let B c be a matrix with 0, 1 entries implementing the global-to-local map on subdomain corners That is, the equation u c = B c u g c, B c = B 1 c B Ns c determines the common values of the degrees of freedom on subdomain corners from a global vector u g c From the construction, the space of all vectors of degrees of freedom continuous across the interfaces can be written as, ImL = {u B r u r = 0, u c ImB c } The problem (23) is reformulated as the equivalent constrained minimization problem 1 2 ut Ku u T f min, subject to B r u r = 0 and u c = B c u g c for some u g c, which is in turn equivalent to finding the stationary point of the Lagrangean L(u i, u r, u g c, λ) = 1 u s 2 vt Kv v T f + u T r Br T λ, v = [v s ], v s i = u s r Bcu s g c Eliminating u s i, us r, and u g c from the Euler equations, we obtain a dual system of the form, cf [5, Eq (14)], (34) F λ = g, and solve it using the preconditioned conjugate gradients method with the preconditioner (35) where (36) M = B r S rr B T r, S rr = diag(s s rr), S s rr = K s rr K s rik s ii 1 K s ir For details of the implementation and numerical results, see [5] 4
6 4 Algebraic bounds In this section, we prove bounds on the condition number of the iterative method defined by Eqs (34) and (35) Denote u = u T u and, for a symmetric positive semidefinite matrix A, denote the induced matrix seminorm u A = u T Au = A 1/2 u If A is known to be positive definite, we write u A instead or u A, because the seminorm is then known to be a norm From the minimization property of the Schur complement, we immediately obtain the following lemma, which characterizes the bilinear form associated with the preconditioner Lemma 41 It holds that u T r S rr u r = min{v t Kv v r = u r, v c = 0} The next lemma gives a more specific description of the matrix of the dual equation (34) Lemma 42 It holds that F = B S 1 r Br T, where the positive definite matrix S is defined by (41) u T r Su r = min{v t Kv v r = u r, v c ImB c } Proof Let L(u r, λ) = min L(u i, u r, u g c, λ) u i,u g c Then, (42) L(u r, λ) = 1 2 ut r Su r + u T r B T r λ u T r h r, with some h r Minimizing over u r, we get u r = S 1 (h r B T r λ) Substituting u r into (42) and taking the variation over λ gives (34) with F = B r S 1 B T r We can now characterize the norm induced by the dual matrix F Lemma 43 It holds that (43) Proof From Lemma 42, λ T F λ = max v r 0 vr T Br T λ 2 v r 2 S λ T F λ = S 1/2 Br T λ 2 w T = max S r 1/2 Br T λ 2 w r 0 w r 2 The substitution w r = S 1/2 v r yields (43) Let V = R p be the space of Lagrange multipliers In this space, define the norm and also the dual norm, µ V = B T r µ Srr (44) λ V = max µ 0 µ T λ µ V Since V = ImB r, substituting µ = B r w r, we can rewrite (44) as (45) wr T Br T λ λ V = max B rw r 0 Br T B r w r Srr 5
7 The main result of this section is the following theorem, which gives a bound on the minimal and maximal eigenvalues of the preconditioned operator M F Theorem 44 If there exists a constant c 1 such that for all w r, then Br T B r w r 2 S rr c 1 w r 2 S, λ max (MF ) λ min (MF ) c 1 4 Proof The proof is based on a comparison of (45) and (43) Using Lemma 43, the substitution v r = B T r B r w r, and the property (32), we find that λ T vr T Br T λ 2 vr T Br T λ 2 wr T Br T λ 2 F λ = max v r 0 v r 2 S max 0 v r ImBr T v Br r 2 S = 4 max B rw r 0 Br T B r w r 2 S Since, by definition, B T r B r w r 2 S B T r B r w r 2 S rr, we conclude, comparing with (45), that (46) λ T F λ 4 λ 2 V On the other hand, from the assumption, we obtain (47) λ T wr T Br T λ 2 F λ = max w r 0 w r 2 S c 1 max B rw r 0 w T r B T r λ B T r B r w r Srr = c 1 λ 2 V Trivially from the definition of M and the norm in V, we have (48) µ 2 V = µ T Mµ Using (46), (47), and (48) in [13, Lemma 31] completes the proof We now show how to verify the assumption of Theorem 44 from inequalities of a form that is more usual in substructuring and easier to estimate for boundary value problems Denote by E s,t the operator that extends the vector of values of degrees of freedom on Γ st, not including corners, by zero entries to a vector of values of degrees of freedom on the whole Ω s, and let E s be the set of all indices of neighbors Ω t of the domain Ω s, with a common edge Γ st Denote by S s the Schur complement on Ω s obtained by eliminating the interior degrees of freedom of Ω s, ie, u st S s u s = min{v st K s v s v s r,c = u s r,c} Our estimate is based on an a-priori bound of the error of approximating a vector of interface degrees of freedom that is continuous across the corners by a vector that is continuous also across the edges In the applications in Sec 5, the approximating vector will be chosen as the natural interpolation on the edges from the corners Theorem 45 Suppose there is a constant c 2 such that for every w r,c, w c ImB c, there exists u r such that B r u r = 0 and, for all s and all t E s, E s,t (w i r u i r) 2 S s c 2 w i r,c 2 S i, i = s, t (49) 6
8 Then, λ max (MF ) λ min (MF ) c 2n e, where n e is the maximum number of the edges of any subdomain Proof Let w r be given and define w c to be the optimal corner degrees of freedom from the definition of S, cf (41) Then (410) N s w r 2 S = wr,c s 2 S s s=1 Let u s r be as in the assumption of the theorem Then B r u r = 0, and, consequently, B T r B r w r = B T r B r (w r u r ) Extending B T r B r w r by zero values of all corner degrees of freedom, we get using the definition of S rr, cf, (36), that (411) N s Br T B r w r 2 S rr = vr,c s 2 S s, s=1 where Using the definition of E s,t, we have v s r = ( B T r B r (w r u r ) ) s, v s c = 0 v s r,c = t E s E s,t v s r Hence, from the triangle inequality, and then using the property (33) of B r and the triangle inequality again, it follows that v s r,c S s t E s E s,t v s r,c S s t E s ( E s,t (w s r u s r) S s + E s,t (w t r u t r) S s) Squaring the the last inequality, using the inequality (a + b) 2 2(a 2 + b 2 ), and the a-priori bound (49) yields vr,c s 2 S 2c ( s 2 w s r,c 2 S + s wt r,c 2 ) S t t E s By the summation over the subdomains and using (411) and (410), we can conclude that N s Br T B r w r 2 S rr = vr,c s 2 S 4c N s 2n s e wr,c s 2 S = 4c 2n s e w r 2 S s=1 s=1 It remains to use Theorem 44 7
9 5 Applications In this section, we verify the assumption of Theorem 45 for two model problems The Sobolev seminorm, denoted by u m,p,x, is the L p (X) norm of the generalized derivatives of of order m 0 of the function u The Sobolev norm is then defined by u m,p,x = [ u 0,p,X,, u m,p,x ] l p Sobolev norms for noninteger m are defined by interpolation Cf, eg, [2, 15] for details and references In particular, u 0,p,X = u 0,p,X = u Lp (X), and, on the boundary Γ of a domain in R 2, u Γ,2,Γ = u(x) u(y) 2 Γ x u 2 dx dy (51) where 51 A second order elliptic problem Consider the boundary value problem Au = g in Ω, u = 0 on Ω, Av = d i,j=1 x i ( α(x) v(x) ), x j with α(x) a measurable function such that 0 < α 0 α(x) α 1 ae in Ω We assume that model problem (51) is discretized using conforming P1 or Q1 elements and denote by Vh P 1 (Ω) the corresponding finite element space that satisfies the usual regularity and inverse properties [2] Let h denote the characteristic element size Assume that all functions in Vh P 1 (Ω) vanish on the boundary of Ω Assume that each subdomain is the union of some of the elements and denote the space of restrictions of functions from Vh P 1 P 1 (Ω) to subdomains by Vh (Ωs ) For every vector of degrees of freedom u s, denote by I P 1 u s the corresponding finite element function The trace of this function is determined by degrees of freedom on the boundary only and we often abuse the notation to define the trace from the boundary values only For simplicity, we assume that the subdomains Ω i, i = 1,, N s form a quasiregular triangulation of the domain Ω Denote the characteristic size of the subdomains by H Finally, let C denote a generic constant independent of h and H Theorem 51 For the model second order problem, it holds that ( λ max (MF ) λ min (MF ) C 1 + log H ) 2, h where the constant does not depend on H and h Proof We only need to verify the assumptions of Theorem 45 For a given w r,c, define u r,c by linear interpolation from w c ; that is, I P 1 u s r,c is linear on all edges of Ω s, and u c = w c Writing w r,c = (w r,c u r,c ) + u r,c, we get from [1, Lemma 35] that I P 1 E s,t (w s r u s r) 1 2,2, Ωs C (1 + log H h ) w i r,c S s Using the uniform equivalence of the seminorms I P 1 vr,c s 1 2,2, Ωs vs r,c S s, cf, eg, [1], and the uniform equivalence of the seminorms v 1 v 2,2, Ωs 1 2,2, Ωt for functions with support on the edge Γ st, which follows from the fact that the subdomains are shape regular, we get (49) with c 2 = C ( ) log H h Remark 52 The assumptions that the subdomains have straight edges can be relaxed to accommodate the case of only shape regular subdomains with edges that are not straight in a standard way by mapping from one or more reference subdomains 8
10 52 A fourth order problem Consider a biharmonic boundary value problem in a variational form: Find u H 2 0 (Ω) such that (52) a(u, v) = f(v), v H 2 0 (Ω), where a(u, v) = Ω 11 u 11 v + 12 u 12 v + 22 u 22 v, u, v H 2 0 (Ω) Let the model problem (52) be discretized by reduced HCT elements [2] We use the same assumptions on the decomposition as in Section 51; in particular, the subdomain edges are straight Let Vh HCT (Ω) be the finite element space of HCT elements satisfying the usual regularity and inverse properties and the essential boundary conditions Note that Vh HCT (Ω) C 1 (Ω) H0 2 (Ω) On each element, a function v in Vh HCT (Ω i ) is determined by the values v(a i ) and the values of its derivatives v x j (a i ), j = 1, 2, at the vertices of the element Denote by t, n, t, and n the tangential and normal directions, and the tangential and normal derivative, respectively The traces of functions from Vh HCT (Ω) on Ω are pairs of functions (u, v) such that u is piecewise cubic, n v is piecewise linear, and u and u are consistent, t u = t v We will abuse notation and write the trace functions simply as (u, u) or just u The space of HCT trace functions is denoted by Vh HCT ( Ω) Denote the finite element interpolation operator by I HCT As in Section 51, we abuse the notation by defining the trace of the interpolation from the boundary degrees of freedom and write, eg, I HCT u s r,c Vh HCT ( Ω s ) We adopt the convention that all functions are understood extended by zero outside of their stated domain; eg, u Γ is zero outside of Γ The following lemma gives a bound on the interpolation from subdomain corners to edges Lemma 53 For every w Vh HCT ( Ω s ), define u Vh HCT ( Ω s ) by HCT interpolation from the corners of Ω s to the edges; that is, u = v and u = w at the corners of Ω s, and, on each edge Γ st, u is a cubic polynomial and n u is linear Then ( (w u) Γ st) ( 1 C 1 + log H ) (53) w 2,2, Ωs 1 h,2, Ωs 2 Proof Denote v = w u Vh HCT ( Ω) Then v = 0 and v = 0 on the corners of Ω s, so from [12, Lemma 41], it follows that ( v Γ st) ( 2 1 v 2 2,2, Ωs 1 + C 1 + log H ) (54) v 2 2,2,Γst 0,,Γ h st Since u Γ st is from a space of dimension 6 and all norms on a finite dimensional space are equivalent, we have in the case when the length of Γ st is one that u 1 2,2,Γst u 1 2,2,Γst u 0,,Γ st Since u 1 and u 2,2,Γst 0,,Γst are invariant to stretching the edge, we get (55) u 1 2,2,Γst C u 0,,Γ st in the general case by scaling We will show at the end of the proof that (56) u 0,,Γ st 5 w 0,,Γ st 9
11 From the discrete Sobolev inequality as generalized by [12, Lemma 42], we have ( w 2 0,,Γ C 1 + log H ) ( w 2 st ) (57) h 2,2,Γst H w 2 0,2,Γ st Combining (54), (56), and (57), we obtain ( v 2 1 C 1 + log H ) ( w 2 2,2, Ωs ) (58) h 2,2, Ωs H w 2 0,2, Ω s First consider the case when Ω s Ω = Then (53) is invariant to adding a linear function to w because v is the error of linear interpolation of w and because this adds only a constant to w So, without loss of generality, let Ω s w = 0 Then by the Poincaré-Friedrichs inequality (59) 1 H y 2 0,2, Ω s C y 2 1 2,2, Ωs, proved first on a reference domain and then scaled [3, 20] Now (59) and (58) give (53) In the case when Ω s Ω, there are some essential boundary conditions on Ω s Since (53) has already been proved in the absence of essential boundary conditions on Ω s, it is sufficient to restrict (53) onto the subspace defined by the boundary conditions, noting that u satisfies the boundary conditions as well It remains to prove (56) The inequality is trivial for n u because the normal derivative is interpolated linearly between the corners Let u L be the linear function on the edge Γ st defined by the values of u on the corners Then, using the triangle inequality, (510) t u 0,,Γ st t (u u L ) 0,,Γ st + t u L 0,,Γ st By a simple computation, we see that the Hermite basis function φ, φ(0) = 0, φ (0) = 1, φ(1) = φ (1) = 0, φ a polynomial of order 3, attains the maximum of φ at 0 Mapping the interval (0, 1) on the edge Γ st and noting that u u L is zero at the endpoints x 1 and x 2 of Γ st, we get using the triangle inequality that (511) t (u u L ) 0,,Γ st t (u u L )(x 1 ) + t (u u L )(x 2 ) 2 t w 0,,Γ st + 2 t u L 0,,Γ st, because t u(x i ) = t w(x i ), i = 1, 2 From the mean value theorem and the fact that u L (x i ) = w(x i ), i = 1, 2, it follows that (512) t u L 0,,Γ st t w 0,,Γ st, hence (511) gives (513) t (u u L ) 0,,Γ st 4 t w 0,,Γ st Now (510), (512), and (513) give (56) Theorem 54 For the fourth order model problem, ( λ max (MF ) λ min (MF ) C 1 + log H ) 2, h where the constant does not depend on H and h 10
12 Proof The proof follows immediately from Theorem 45, with (49) being a consequence of Lemma 53 and the uniform equivalence of seminorms I HCT wr,c s 1 2,2, Ωs wr,c s S s, cf, [12] Remark 55 The result extends, by spectral equivalence, to DKT elements and a certain class of non-locking elements for the Reissner-Mindlin plate model as in [12] The result also extends to the case when the subdomain edges are not straight by considering the subdomains to be images of a reference domain REFERENCES [1] J H Bramble, J E Pasciak, and A H Schatz, The construction of preconditioners for elliptic problems by substructuring, I, Math Comp, 47 (1986), pp [2] P G Ciarlet, Basic error estimates for elliptic problems, in Handbook of Numerical Analysis, P Ciarlet and J L Lions, eds, vol II, North-Holland, Amsterdam, 1989, pp [3] M Dryja, A method of domain decomposition for 3-D finite element problems, in First International Symposium on Domain Decomposition Methods for Partial Differential Equations, R Glowinski, G H Golub, G A Meurant, and J Périaux, eds, Philadelphia, PA, 1988, SIAM [4] C Farhat, P-S Chen, J Mandel, and F-X Roux, The two-level FETI method - Part II: Extension to shell problems, parallel implementation and performance results, Comp Meth Appl Mech Engrg, 155 (1998), pp [5] C Farhat, M Lesoinne, P LeTallec, K Pierson, and D Rixen, FETI-DP: A dualprimal unified FETI method Part I: A faster alternative to the two-level FETI method, Tech Rep CU-CAS-99-15, Center for Aerospace Structures, University of Colorado at Boulder, August 1999 [6] C Farhat and J Mandel, Scalable substructuring by Lagrange multipliers in theory and practice, in Proceedings of the 9th International Conference on Domain Decomposition, Bergen, Norway, June 1996, DDMORG, 1998 wwwddmorg [7], The two-level FETI method for static and dynamic plate problems - Part I: An optimal iterative solver for biharmonic systems, Comp Meth Appl Mech Engrg, 155 (1998), pp [8] C Farhat, J Mandel, and F-X Roux, Optimal convergence properties of the FETI domain decomposition method, Comput Methods Appl Mech Engrg, 115 (1994), pp [9] C Farhat and F-X Roux, An unconventional domain decomposition method for an efficient parallel solution of large-scale finite element systems, SIAM J Sci Stat Comput, 13 (1992), pp [10] A Klawonn and O B Widlund, A domain decomposition method with Lagrange multipliers for linear elasticity, Tech Rep TR 780, Courant Institute, New York University, February 1999 [11] P Le Tallec, J Mandel, and M Vidrascu, Balancing domain decomposition for plates, Contemporary Mathematics, 180 (1994), pp Proceedings of the 7th International Symposium on Domain Decomposition Methods, Penn State, November 1993 [12], A Neumann-Neumann domain decomposition algorithm for solving plate and shell problems, SIAM J Numer Anal, 35 (1998), pp [13] J Mandel and R Tezaur, Convergence of a substructuring method with Lagrange multipliers, Numerische Mathematik, 73 (1996), pp [14] J Mandel, R Tezaur, and C Farhat, A scalable substructuring method by Lagrange multipliers for plate bending problems, SIAM J Numer Anal, 36 (1999), pp [15] J Nečas, Les méthodes directes en théorie des équations elliptiques, Academia, Prague, 1967 [16] K C Park, M R Justino, Jr, and C A Felippa, An algebraically partitioned FETI method for parallel structural analysis: Algorithm description, Int J Numer Meth Engrg, 40 (1997), pp [17] D Rixen and C Farhat, A simple and efficient extension of a class of substructure based preconditioners to heterogeneous structural mechanics problems, Int J Numer Meth Engrg, 46 (1999), pp [18] D J Rixen, C Farhat, R Tezaur, and J Mandel, Theoretical comparison of the FETI and algebraically partitioned FETI methods, and performance comparisons with a direct sparse solver, Intl Journal for Numerical Methods in Engineering, 46 (1999), pp [19] R Tezaur, Analysis of Lagrange Multiplier Based Domain Decomposition, PhD thesis, Uni- 11
13 versity of Colorado at Denver, [20] O B Widlund, Iterative substructuring methods: Algorithms and theory for elliptic problems in the plane, in First International Symposium on Domain Decomposition Methods for Partial Differential Equations, R Glowinski, G H Golub, G A Meurant, and J Périaux, eds, Philadelphia, PA, 1988, SIAM, pp
ON THE CONVERGENCE OF A DUAL-PRIMAL SUBSTRUCTURING METHOD. January 2000 Revised April 2000
ON THE CONVERGENCE OF A DUAL-PRIMAL SUBSTRUCTURING METHOD JAN MANDEL AND RADEK TEZAUR January 2000 Revised April 2000 Abstract In te Dual-Primal FETI metod, introduced by Farat et al [5], te domain is
More informationMultispace and Multilevel BDDC. Jan Mandel University of Colorado at Denver and Health Sciences Center
Multispace and Multilevel BDDC Jan Mandel University of Colorado at Denver and Health Sciences Center Based on joint work with Bedřich Sousedík, UCDHSC and Czech Technical University, and Clark R. Dohrmann,
More informationMultilevel and Adaptive Iterative Substructuring Methods. Jan Mandel University of Colorado Denver
Multilevel and Adaptive Iterative Substructuring Methods Jan Mandel University of Colorado Denver The multilevel BDDC method is joint work with Bedřich Sousedík, Czech Technical University, and Clark Dohrmann,
More informationShort title: Total FETI. Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ Ostrava, Czech Republic
Short title: Total FETI Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ-70833 Ostrava, Czech Republic mail: zdenek.dostal@vsb.cz fax +420 596 919 597 phone
More informationAdaptive Coarse Space Selection in BDDC and FETI-DP Iterative Substructuring Methods: Towards Fast and Robust Solvers
Adaptive Coarse Space Selection in BDDC and FETI-DP Iterative Substructuring Methods: Towards Fast and Robust Solvers Jan Mandel University of Colorado at Denver Bedřich Sousedík Czech Technical University
More informationA FETI-DP Method for Mortar Finite Element Discretization of a Fourth Order Problem
A FETI-DP Method for Mortar Finite Element Discretization of a Fourth Order Problem Leszek Marcinkowski 1 and Nina Dokeva 2 1 Department of Mathematics, Warsaw University, Banacha 2, 02 097 Warszawa, Poland,
More informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY 11, USA Dan Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
More informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY, USA. Dan_Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
More informationDiagonal Representation of Certain Matrices
Diagonal Representation of Certain Matrices Mark Tygert Research Report YALEU/DCS/RR-33 December 2, 2004 Abstract An explicit expression is provided for the characteristic polynomial of a matrix M of the
More informationReport Documentation Page
Report Documentation Page Form Approved OMB No. 0704-0188 Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,
More information20. A Dual-Primal FETI Method for solving Stokes/Navier-Stokes Equations
Fourteenth International Conference on Domain Decomposition Methods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 23 DDM.org 2. A Dual-Primal FEI Method for solving Stokes/Navier-Stokes
More informationParallel scalability of a FETI DP mortar method for problems with discontinuous coefficients
Parallel scalability of a FETI DP mortar method for problems with discontinuous coefficients Nina Dokeva and Wlodek Proskurowski University of Southern California, Department of Mathematics Los Angeles,
More informationParallel Scalability of a FETI DP Mortar Method for Problems with Discontinuous Coefficients
Parallel Scalability of a FETI DP Mortar Method for Problems with Discontinuous Coefficients Nina Dokeva and Wlodek Proskurowski Department of Mathematics, University of Southern California, Los Angeles,
More informationUse of Wijsman's Theorem for the Ratio of Maximal Invariant Densities in Signal Detection Applications
Use of Wijsman's Theorem for the Ratio of Maximal Invariant Densities in Signal Detection Applications Joseph R. Gabriel Naval Undersea Warfare Center Newport, Rl 02841 Steven M. Kay University of Rhode
More informationSelecting Constraints in Dual-Primal FETI Methods for Elasticity in Three Dimensions
Selecting Constraints in Dual-Primal FETI Methods for Elasticity in Three Dimensions Axel Klawonn 1 and Olof B. Widlund 2 1 Universität Duisburg-Essen, Campus Essen, Fachbereich Mathematik, (http://www.uni-essen.de/ingmath/axel.klawonn/)
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 31, pp. 384-402, 2008. Copyright 2008,. ISSN 1068-9613. ETNA ON THE EQUIVALENCE OF PRIMAL AND DUAL SUBSTRUCTURING PRECONDITIONERS BEDŘICH SOUSEDÍK
More informationThe All-floating BETI Method: Numerical Results
The All-floating BETI Method: Numerical Results Günther Of Institute of Computational Mathematics, Graz University of Technology, Steyrergasse 30, A-8010 Graz, Austria, of@tugraz.at Summary. The all-floating
More informationEstimation of Vertical Distributions of Water Vapor and Aerosols from Spaceborne Observations of Scattered Sunlight
Estimation of Vertical Distributions of Water Vapor and Aerosols from Spaceborne Observations of Scattered Sunlight Dale P. Winebrenner Polar Science Center/Applied Physics Laboratory University of Washington
More informationAn Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions
An Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions Leszek Marcinkowski Department of Mathematics, Warsaw University, Banacha
More informationAnalysis Comparison between CFD and FEA of an Idealized Concept V- Hull Floor Configuration in Two Dimensions. Dr. Bijan Khatib-Shahidi & Rob E.
Concept V- Hull Floor Configuration in Two Dimensions Dr. Bijan Khatib-Shahidi & Rob E. Smith 10 November 2010 : Dist A. Approved for public release Report Documentation Page Form Approved OMB No. 0704-0188
More informationCONVERGENCE ANALYSIS OF A BALANCING DOMAIN DECOMPOSITION METHOD FOR SOLVING A CLASS OF INDEFINITE LINEAR SYSTEMS
CONVERGENCE ANALYSIS OF A BALANCING DOMAIN DECOMPOSITION METHOD FOR SOLVING A CLASS OF INDEFINITE LINEAR SYSTEMS JING LI AND XUEMIN TU Abstract A variant of balancing domain decomposition method by constraints
More informationAttribution Concepts for Sub-meter Resolution Ground Physics Models
Attribution Concepts for Sub-meter Resolution Ground Physics Models 76 th MORS Symposium US Coast Guard Academy Approved for public release distribution. 2 Report Documentation Page Form Approved OMB No.
More informationAn additive average Schwarz method for the plate bending problem
J. Numer. Math., Vol. 10, No. 2, pp. 109 125 (2002) c VSP 2002 Prepared using jnm.sty [Version: 02.02.2002 v1.2] An additive average Schwarz method for the plate bending problem X. Feng and T. Rahman Abstract
More informationAn Efficient FETI Implementation on Distributed Shared Memory Machines with Independent Numbers of Subdomains and Processors
Contemporary Mathematics Volume 218, 1998 B 0-8218-0988-1-03024-7 An Efficient FETI Implementation on Distributed Shared Memory Machines with Independent Numbers of Subdomains and Processors Michel Lesoinne
More informationP. Kestener and A. Arneodo. Laboratoire de Physique Ecole Normale Supérieure de Lyon 46, allée d Italie Lyon cedex 07, FRANCE
A wavelet-based generalization of the multifractal formalism from scalar to vector valued d- dimensional random fields : from theoretical concepts to experimental applications P. Kestener and A. Arneodo
More informationSums of the Thue Morse sequence over arithmetic progressions
Sums of the Thue Morse sequence over arithmetic progressions T.W. Cusick 1, P. Stănică 2 1 State University of New York, Department of Mathematics Buffalo, NY 14260; Email: cusick@buffalo.edu 2 Naval Postgraduate
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationDevelopment and Application of Acoustic Metamaterials with Locally Resonant Microstructures
Development and Application of Acoustic Metamaterials with Locally Resonant Microstructures AFOSR grant #FA9550-10-1-0061 Program manager: Dr. Les Lee PI: C.T. Sun Purdue University West Lafayette, Indiana
More informationClosed-form and Numerical Reverberation and Propagation: Inclusion of Convergence Effects
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Closed-form and Numerical Reverberation and Propagation: Inclusion of Convergence Effects Chris Harrison Centre for Marine
More informationGRUPO DE GEOFÍSICA MATEMÁTICA Y COMPUTACIONAL MEMORIA Nº 8
GRUPO DE GEOFÍSICA MATEMÁTICA Y COMPUTACIONAL MEMORIA Nº 8 MÉXICO 2014 MEMORIAS GGMC INSTITUTO DE GEOFÍSICA UNAM MIEMBROS DEL GGMC Dr. Ismael Herrera Revilla iherrerarevilla@gmail.com Dr. Luis Miguel de
More informationSOME PRACTICAL ASPECTS OF PARALLEL ADAPTIVE BDDC METHOD
Conference Applications of Mathematics 2012 in honor of the 60th birthday of Michal Křížek. Institute of Mathematics AS CR, Prague 2012 SOME PRACTICAL ASPECTS OF PARALLEL ADAPTIVE BDDC METHOD Jakub Šístek1,2,
More informationREGENERATION OF SPENT ADSORBENTS USING ADVANCED OXIDATION (PREPRINT)
AL/EQ-TP-1993-0307 REGENERATION OF SPENT ADSORBENTS USING ADVANCED OXIDATION (PREPRINT) John T. Mourand, John C. Crittenden, David W. Hand, David L. Perram, Sawang Notthakun Department of Chemical Engineering
More informationHigh-Fidelity Computational Simulation of Nonlinear Fluid- Structure Interaction Problems
Aerodynamic Issues of Unmanned Air Vehicles Fluid-Structure Interaction High-Fidelity Computational Simulation of Nonlinear Fluid- Structure Interaction Problems Raymond E. Gordnier Computational Sciences
More informationXIAO-CHUAN CAI AND MAKSYMILIAN DRYJA. strongly elliptic equations discretized by the nite element methods.
Contemporary Mathematics Volume 00, 0000 Domain Decomposition Methods for Monotone Nonlinear Elliptic Problems XIAO-CHUAN CAI AND MAKSYMILIAN DRYJA Abstract. In this paper, we study several overlapping
More informationA GENERAL PROCEDURE TO SET UP THE DYADIC GREEN S FUNCTION OF MULTILAYER CONFORMAL STRUCTURES AND ITS APPLICATION TO MICROSTRIP ANTENNAS
A GENERAL PROCEDURE TO SET UP THE DYADIC GREEN S FUNCTION OF MULTILAYER CONFORMAL STRUCTURES AND ITS APPLICATION TO MICROSTRIP ANTENNAS Michael Thiel*, Truong Vu Bang Giang and Achim Dreher Institute of
More informationConvergence analysis of a balancing domain decomposition method for solving a class of indefinite linear systems
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2000; 00:1 6 [Version: 2002/09/18 v1.02] Convergence analysis of a balancing domain decomposition method for solving a class of indefinite
More informationREPORT DOCUMENTATION PAGE
REPORT DOCUMENTATION PAGE Form Approved OMB No. 0704-0188 The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,
More informationSystem Reliability Simulation and Optimization by Component Reliability Allocation
System Reliability Simulation and Optimization by Component Reliability Allocation Zissimos P. Mourelatos Professor and Head Mechanical Engineering Department Oakland University Rochester MI 48309 Report
More informationAn Iterative Domain Decomposition Method for the Solution of a Class of Indefinite Problems in Computational Structural Dynamics
An Iterative Domain Decomposition Method for the Solution of a Class of Indefinite Problems in Computational Structural Dynamics Charbel Farhat a, and Jing Li a a Department of Aerospace Engineering Sciences
More informationApplication of Preconditioned Coupled FETI/BETI Solvers to 2D Magnetic Field Problems
Application of Preconditioned Coupled FETI/BETI Solvers to 2D Magnetic Field Problems U. Langer A. Pohoaţǎ O. Steinbach 27th September 2004 Institute of Computational Mathematics Johannes Kepler University
More informationWavelet Spectral Finite Elements for Wave Propagation in Composite Plates
Wavelet Spectral Finite Elements for Wave Propagation in Composite Plates Award no: AOARD-0904022 Submitted to Dr Kumar JATA Program Manager, Thermal Sciences Directorate of Aerospace, Chemistry and Materials
More informationCrowd Behavior Modeling in COMBAT XXI
Crowd Behavior Modeling in COMBAT XXI Imre Balogh MOVES Research Associate Professor ilbalogh@nps.edu July 2010 831-656-7582 http://movesinstitute.org Report Documentation Page Form Approved OMB No. 0704-0188
More informationNumerische Mathematik
Numer. Math. (1996) 75: 59 77 Numerische Mathematik c Springer-Verlag 1996 Electronic Edition A preconditioner for the h-p version of the finite element method in two dimensions Benqi Guo 1, and Weiming
More informationDomain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions
Domain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions Ernst P. Stephan 1, Matthias Maischak 2, and Thanh Tran 3 1 Institut für Angewandte Mathematik, Leibniz
More informationarxiv: v1 [math.na] 28 Feb 2008
BDDC by a frontal solver and the stress computation in a hip joint replacement arxiv:0802.4295v1 [math.na] 28 Feb 2008 Jakub Šístek a, Jaroslav Novotný b, Jan Mandel c, Marta Čertíková a, Pavel Burda a
More informationA report (dated September 20, 2011) on. scientific research carried out under Grant: FA
A report (dated September 20, 2011) on scientific research carried out under Grant: FA2386-10-1-4150 First-principles determination of thermal properties in nano-structured hexagonal solids with doping
More informationUSMC Enlisted Endstrength Model
USMC Enlisted Endstrength Model MORS Workshop Personnel and National Security: A Quantitative Approach Working Group #2: Models for Managing Retention Molly F. McIntosh January 27, 2009 Report Documentation
More informationREPORT DOCUMENTATION PAGE. Theoretical Study on Nano-Catalyst Burn Rate. Yoshiyuki Kawazoe (Tohoku Univ) N/A AOARD UNIT APO AP
REPORT DOCUMENTATION PAGE Form Approved OMB No. 0704-0188 The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,
More informationFRACTAL CONCEPTS AND THE ANALYSIS OF ATMOSPHERIC PROCESSES
1 FRACTAL CONCEPTS AND THE ANALYSIS OF ATMOSPHERIC PROCESSES Robert L. Street Environmental Fluid Mechanics Laboratory Department of Civil Engineering Stanford University Stanford, CA 94305-4020 650-723-4969;
More informationReport Documentation Page
Report Documentation Page Form Approved OMB No. 0704-0188 Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,
More informationOptimizing Robotic Team Performance with Probabilistic Model Checking
Optimizing Robotic Team Performance with Probabilistic Model Checking Software Engineering Institute Carnegie Mellon University Pittsburgh, PA 15213 Sagar Chaki, Joseph Giampapa, David Kyle (presenting),
More informationAuxiliary space multigrid method for elliptic problems with highly varying coefficients
Auxiliary space multigrid method for elliptic problems with highly varying coefficients Johannes Kraus 1 and Maria Lymbery 2 1 Introduction The robust preconditioning of linear systems of algebraic equations
More informationTHE EULER FUNCTION OF FIBONACCI AND LUCAS NUMBERS AND FACTORIALS
Annales Univ. Sci. Budapest., Sect. Comp. 40 (2013) nn nnn THE EULER FUNCTION OF FIBONACCI AND LUCAS NUMBERS AND FACTORIALS Florian Luca (Morelia, Mexico) Pantelimon Stănică (Monterey, USA) Dedicated to
More informationuniform distribution theory
Uniform Distribution Theory 9 (2014), no.1, 21 25 uniform distribution theory ON THE FIRST DIGITS OF THE FIBONACCI NUMBERS AND THEIR EULER FUNCTION Florian Luca Pantelimon Stănică ABSTRACT. Here, we show
More informationCRS Report for Congress
CRS Report for Congress Received through the CRS Web Order Code RS21396 Updated May 26, 2006 Summary Iraq: Map Sources Hannah Fischer Information Research Specialist Knowledge Services Group This report
More informationNormic continued fractions in totally and tamely ramified extensions of local fields
Normic continued fractions in totally and tamely ramified extensions of local fields Pantelimon Stănică Naval Postgraduate School Applied Mathematics Department, Monterey, CA 93943 5216, USA; email: pstanica@nps.edu
More informationBabylonian resistor networks
IOP PUBLISHING Eur. J. Phys. 33 (2012) 531 537 EUROPEAN JOURNAL OF PHYSICS doi:10.1088/0143-0807/33/3/531 Babylonian resistor networks Carl E Mungan 1 and Trevor C Lipscombe 2 1 Physics Department, US
More informationfor three dimensional problems are often more complicated than the quite simple constructions that work well for problems in the plane; see [23] for a
DUAL-PRIMAL FETI METHODS FOR THREE-DIMENSIONAL ELLIPTIC PROBLEMS WITH HETEROGENEOUS COEFFICIENTS AEL KLAWONN Λ, OLOF B. WIDLUND y, AND MAKSYMILIAN DRYJA z Abstract. In this paper, certain iterative substructuring
More informationBroadband matched-field source localization in the East China Sea*
Broadband matched-field source localization in the East China Sea* Renhe Zhang Zhenglin Li Jin Yan Zhaohui Peng Fenghua Li National Laboratory of Acoustics, Institute of Acoustics, Chinese Academy of Sciences,
More information7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION REPORT NUMBER
REPORT DOCUMENTATION PAGE Form Approved OMB No. 0704-0188 The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,
More informationDetermining the Stratification of Exchange Flows in Sea Straits
Determining the Stratification of Exchange Flows in Sea Straits Lawrence J. Pratt Physical Oceanography Department, MS #21 Woods Hole Oceanographic Institution, Woods Hole, MA 02543 phone: (508) 289-2540
More informationMetrology Experiment for Engineering Students: Platinum Resistance Temperature Detector
Session 1359 Metrology Experiment for Engineering Students: Platinum Resistance Temperature Detector Svetlana Avramov-Zamurovic, Carl Wick, Robert DeMoyer United States Naval Academy Abstract This paper
More informationMultispace and Multilevel BDDC
Multispace and Multilevel BDDC Jan Mandel Bedřich Sousedík Clark R. Dohrmann February 11, 2018 arxiv:0712.3977v2 [math.na] 21 Jan 2008 Abstract BDDC method is the most advanced method from the Balancing
More informationExtending the theory for domain decomposition algorithms to less regular subdomains
Extending the theory for domain decomposition algorithms to less regular subdomains Olof Widlund Courant Institute of Mathematical Sciences New York University http://www.cs.nyu.edu/cs/faculty/widlund/
More informationNonoverlapping Domain Decomposition Methods with Simplified Coarse Spaces for Solving Three-dimensional Elliptic Problems
Nonoverlapping Domain Decomposition Methods with Simplified Coarse Spaces for Solving Three-dimensional Elliptic Problems Qiya Hu 1, Shi Shu 2 and Junxian Wang 3 Abstract In this paper we propose a substructuring
More informationGrant Number: N IP To compare obtained theoretical results with NPAL experimental data.
Coherence of Low-Frequency Sound Signals Propagating through a Fluctuating Ocean: Analysis and Theoretical Interpretation of 2004 NPAL Experimental Data Alexander G. Voronovich NOAA/ESRL, PSD4, 325 Broadway,
More informationSensitivity of West Florida Shelf Simulations to Initial and Boundary Conditions Provided by HYCOM Data-Assimilative Ocean Hindcasts
Sensitivity of West Florida Shelf Simulations to Initial and Boundary Conditions Provided by HYCOM Data-Assimilative Ocean Hindcasts George Halliwell, MPO/RSMAS, University of Miami Alexander Barth, University
More informationDomain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions
Domain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions Bernhard Hientzsch Courant Institute of Mathematical Sciences, New York University, 51 Mercer Street, New
More informationASM-BDDC Preconditioners with variable polynomial degree for CG- and DG-SEM
ASM-BDDC Preconditioners with variable polynomial degree for CG- and DG-SEM C. Canuto 1, L. F. Pavarino 2, and A. B. Pieri 3 1 Introduction Discontinuous Galerkin (DG) methods for partial differential
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 11, pp. 1-24, 2000. Copyright 2000,. ISSN 1068-9613. ETNA NEUMANN NEUMANN METHODS FOR VECTOR FIELD PROBLEMS ANDREA TOSELLI Abstract. In this paper,
More informationGroupe de Travail Méthodes Numériques
Groupe de Travail Méthodes Numériques Laboratoire Jacques-Louis Lions, Paris 6 Monday, May 11th, 2009 Domain Decomposition with Lagrange Multipliers: A continuous framework for a FETI-DP+Mortar method
More informationOn the Use of Inexact Subdomain Solvers for BDDC Algorithms
On the Use of Inexact Subdomain Solvers for BDDC Algorithms Jing Li a, and Olof B. Widlund b,1 a Department of Mathematical Sciences, Kent State University, Kent, OH, 44242-0001 b Courant Institute of
More informationOcean Acoustics Turbulence Study
Ocean Acoustics Turbulence Study PI John Oeschger Coastal Systems Station 6703 West Highway 98 Panama City, FL 32407 phone: (850) 230-7054 fax: (850) 234-4886 email: OeschgerJW@ncsc.navy.mil CO-PI Louis
More informationEstimation of Vertical Distributions of Water Vapor from Spaceborne Observations of Scattered Sunlight
Estimation of Vertical Distributions of Water Vapor from Spaceborne Observations of Scattered Sunlight Dale P. Winebrenner Applied Physics Laboratory, Box 355640 University of Washington Seattle, WA 98195
More informationNUMERICAL SOLUTIONS FOR OPTIMAL CONTROL PROBLEMS UNDER SPDE CONSTRAINTS
NUMERICAL SOLUTIONS FOR OPTIMAL CONTROL PROBLEMS UNDER SPDE CONSTRAINTS AFOSR grant number: FA9550-06-1-0234 Yanzhao Cao Department of Mathematics Florida A & M University Abstract The primary source of
More informationSW06 Shallow Water Acoustics Experiment Data Analysis
DISTRIBUTION STATEMENT A: Approved for public release; distribution is unlimited. SW06 Shallow Water Acoustics Experiment Data Analysis James F. Lynch MS #12, Woods Hole Oceanographic Institution Woods
More informationReport Documentation Page
Inhibition of blood cholinesterase activity is a poor predictor of acetylcholinesterase inhibition in brain regions of guinea pigs exposed to repeated doses of low levels of soman. Sally M. Anderson Report
More informationDIRECTIONAL WAVE SPECTRA USING NORMAL SPREADING FUNCTION
CETN-I-6 3185 DIRECTIONAL WAVE SPECTRA USING NORMAL SPREADING FUNCTION PURPOSE : To present a parameterized model of a directional spectrum of the sea surface using an energy spectrum and a value for the
More information33 RASHO: A Restricted Additive Schwarz Preconditioner with Harmonic Overlap
Thirteenth International Conference on Domain Decomposition ethods Editors: N. Debit,.Garbey, R. Hoppe, J. Périaux, D. Keyes, Y. Kuznetsov c 001 DD.org 33 RASHO: A Restricted Additive Schwarz Preconditioner
More informationQuantitation and Ratio Determination of Uranium Isotopes in Water and Soil Using Inductively Coupled Plasma Mass Spectrometry (ICP-MS)
Quantitation and Ratio Determination of Uranium Isotopes in Water and Soil Using Inductively Coupled Plasma Mass Spectrometry (ICP-MS) D.N. Kurk, T.E. Beegle, S.C. Spence and R.J. Swatski Report Documentation
More informationReal-Time Environmental Information Network and Analysis System (REINAS)
Calhoun: The NPS Institutional Archive Faculty and Researcher Publications Faculty and Researcher Publications 1998-09 Real-Time Environmental Information Network and Analysis System (REINAS) Nuss, Wendell
More informationParametric Models of NIR Transmission and Reflectivity Spectra for Dyed Fabrics
Naval Research Laboratory Washington, DC 20375-5320 NRL/MR/5708--15-9629 Parametric Models of NIR Transmission and Reflectivity Spectra for Dyed Fabrics D. Aiken S. Ramsey T. Mayo Signature Technology
More informationPIPS 3.0. Pamela G. Posey NRL Code 7322 Stennis Space Center, MS Phone: Fax:
PIPS 3.0 Ruth H. Preller Naval Research Laboratory, Code 7322 Stennis Space Center, MS 39529 phone: (228) 688-5444 fax: (228)688-4759 email: preller@nrlssc.navy.mil Pamela G. Posey NRL Code 7322 Stennis
More informationAnalysis of Infrared Measurements of Microbreaking and Whitecaps
Analysis of Infrared Measurements of Microbreaking and Whitecaps Andrew T. Jessup Applied Physics Laboratory, University of Washington 1013 NE 40th St. Seattle, WA 98105-6698 phone (206) 685-2609 fax (206)
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES. A virtual overlapping Schwarz method for scalar elliptic problems in two dimensions
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES A virtual overlapping Schwarz method for scalar elliptic problems in two dimensions Juan Gabriel Calvo Preprint No. 25-2017 PRAHA 2017 A VIRTUAL
More informationOn Applying Point-Interval Logic to Criminal Forensics
On Applying Point-Interval Logic to Criminal Forensics (Student Paper) Mashhood Ishaque Abbas K. Zaidi Alexander H. Levis Command and Control Research and Technology Symposium 1 Report Documentation Page
More informationSuper-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models PI Isaac Ginis Graduate School of Oceanography
More informationCalculation of the Viscous Drag on Submerged Bodies From the Far Field Behavior of the Flow
Calculation of the Viscous Drag on Submerged Bodies From the Far Field Behavior of the Flow Robert T. Hudspeth 202 Apperson Hall Department of Civil Engineering Oregon State University Corvallis, OR 97331-2302
More informationStudy of Electromagnetic Scattering From Material Object Doped Randomely WithThin Metallic Wires Using Finite Element Method
Study of Electromagnetic Scattering From Material Object Doped Randomely WithThin Metallic Wires Using Finite Element Method Manohar D. Deshpande, NASA Langley Research Center, Hmapton, VA., USA Abstract
More informationREPORT DOCUMENTATION PAGE
REPORT DOCUMENTATION PAGE Form Approved OMB No. 0704-0188 Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,
More informationDomain Decomposition solvers (FETI)
Domain Decomposition solvers (FETI) a random walk in history and some current trends Daniel J. Rixen Technische Universität München Institute of Applied Mechanics www.amm.mw.tum.de rixen@tum.de 8-10 October
More informationImprovements in Modeling Radiant Emission from the Interaction Between Spacecraft Emanations and the Residual Atmosphere in LEO
Improvements in Modeling Radiant Emission from the Interaction Between Spacecraft Emanations and the Residual Atmosphere in LEO William L. Dimpfl Space Science Applications Laboratory The Aerospace Corporation
More informationUnderstanding Near-Surface and In-cloud Turbulent Fluxes in the Coastal Stratocumulus-topped Boundary Layers
Understanding Near-Surface and In-cloud Turbulent Fluxes in the Coastal Stratocumulus-topped Boundary Layers Qing Wang Meteorology Department, Naval Postgraduate School Monterey, CA 93943 Phone: (831)
More informationAn Invariance Property of the Generalized Likelihood Ratio Test
352 IEEE SIGNAL PROCESSING LETTERS, VOL. 10, NO. 12, DECEMBER 2003 An Invariance Property of the Generalized Likelihood Ratio Test Steven M. Kay, Fellow, IEEE, and Joseph R. Gabriel, Member, IEEE Abstract
More informationA Neumann-Dirichlet Preconditioner for FETI-DP 2 Method for Mortar Discretization of a Fourth Order 3 Problems in 2D 4 UNCORRECTED PROOF
1 A Neumann-Dirichlet Preconditioner for FETI-DP 2 Method for Mortar Discretization of a Fourth Order 3 Problems in 2D 4 Leszek Marcinkowski * 5 Faculty of Mathematics, University of Warsaw, Banacha 2,
More informationA Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
More informationCoastal Engineering Technical Note
Coastal Engineering Technical Note CETN I-48 (12/91) Evaluation and Application of the Wave Information Study for the Gulf of Mexico INTRODUCTION The Wave Information Study (WIS) for the Gulf of Mexico
More informationComputer Simulation of Sand Ripple Growth and Migration.
Computer Simulation of Sand Ripple Growth and Migration. Douglas J. Wilson OGI School of Environmental Science and Engineering at the Oregon Health and Sciences University 20000 N.W. Walker Road, Beaverton,
More information18. Balancing Neumann-Neumann for (In)Compressible Linear Elasticity and (Generalized) Stokes Parallel Implementation
Fourteenth nternational Conference on Domain Decomposition Methods Editors: smael Herrera, David E Keyes, Olof B Widlund, Robert Yates c 23 DDMorg 18 Balancing Neumann-Neumann for (n)compressible Linear
More informationEvolution of Tropical Cyclone Characteristics
Evolution of Tropical Cyclone Characteristics Patrick A. Harr Department of Meteorology Naval Postgraduate School Monterey, CA 93943-5114 Telephone : (831) 656-3787 FAX: (831) 656-3061 email: paharr@nps.navy.mil
More information