Preconditioning techniques for Newton s method for the incompressible Navier Stokes equations
|
|
- Felicity Harvey
- 5 years ago
- Views:
Transcription
1 Preconditioning techniques for Newton s method for the incomressible Navier Stokes equations H. C. ELMAN 1, D. LOGHIN 2 and A. J. WATHEN 3 1 Deartment of Comuter Science, University of Maryland, College Park, MD 20742, USA, elman@cs.umd.edu 2 Parallel Algorithms Project, CERFACS, 42 ave G. Coriolis, Toulouse, 31057, France, loghin@cerfacs.fr 3 Oxford University Comuting Laboratory, Parks Road, Oxford, OX1 3QD, wathen@comlab.ox.ac.uk CERFACS Technical Reort TR/PA/03/31 Abstract Newton s method for the incomressible Navier Stokes equations gives rise to large sarse non-symmetric indefinite matrices with a so-called saddle-oint structure for which Schur comlement reconditioners have roven to be effective when couled with iterative methods of Krylov tye. In this work we investigate the erformance of two reconditioning techniques introduced originally for the Picard method for which both roved significantly suerior to other aroaches such as the Uzawa method. The first is a block reconditioner which is based on the algebraic structure of the system matrix [4]. The other aroach uses also a block reconditioner which is derived by considering the underlying artial differential oerator matrix [9]. Analysis and numerical comarison of the methods are resented. Keywords: Block reconditioners, Navier Stokes, Newton s method. 1 Problem descrition Mixed finite element discretizations of Navier Stokes equations give rise to nonlinear systems with the following saddle-oint structure (see, e.g.,[7], [14]) ( ( ) ( u F (u) B t u K(u) = = f (1) ) B 0 ) where K(u) R n n, F (u) R n1 n1 is nonsymmetric and ossibly indefinite and B t R n1 n2 has nontrivial kernel sanned by the constant vector. 1
2 Since K is large and sarse, a cometitive class of solution methods is that of iterative solvers of Krylov tye combined with suitable reconditioning. In this work we are interested in the case where a Newton linearization is emloyed to solve (1); this leads to an iteration of the form [14, 10.3] ( ) ( ) ( ) J(u k δu k+1 F (u ) = k ) + M(u k ) B t δu k+1 B 0 δ k+1 = r k, (2) δ k+1 where M(u) = A u (u) u, u k+1 = u k + δu k+1, r k = f K(u k )u k. On the other hand, the Jacobian matrix J(u k ) is related to the system matrix arising from a Picard linearization of (1), ( ) ( ) ( ) K(u k δu k+1 F (u ) = k ) B t δu k+1 B 0 δ k+1 = r k, (3) δ k+1 for which efficient reconditioners are available. Our aim is to adat and analyze the reconditioned iterative solvers devised for the Picard iteration for roblems where Newton s method is emloyed. Henceforth, we will use the notation F (ρ, u k ) = F (u k ) + ρm(u k ); (4) where ρ {0, 1} to denote the matrices arising from the above two discretizations. To comlete the icture, we note here that for flow roblems in R d (d = 2, 3) there holds F (0, u k ) = I d F (u k ), where F (u k ) is nonsymmetric and ositive-definite and comes from the discretization of an advection-diffusion oerator F (u k ) = νl + N(u k ), where ν is the diffusion arameter and L is a discrete Lalacian. The matrix M(u k ) is also a d d block-matrix and reresents the discretization of a zero-order oerator. Hence it is well-conditioned; however, deending on the roblem, it can be indefinite which leads to the ossibility of F (1, u k ) being indefinite. Finally, we note that the arameter ν is very imortant in ractical alications. Its recirocal, ν, is the Reynolds number and given the size of the roblem n, there are arameters ν for which the nonlinear roblem (1) does not have a unique solution. These critical values of ν where non-uniqueness sets in can be detected in various ways, e.g., continuation methods. For this reason, otimality of our solvers will be measured in terms of deendence of the number of iterations on two arameters: n and ν. Whether Picard, Newton or variations thereof, linearizations of (1) reserve the saddle-oint structure of the system matrix for which Schur comlement reconditioners erform otimally. More recisely, it is known that the key to efficient reconditioning is the Schur comlement S = BF (ρ, u k ) B t ; (5) 2
3 in articular, a block (right) reconditioner of the form ( ) F (ρ, u P = k ) B t 0 S leads to convergence of any Krylov subsace iterative solver in 2 iterations [13], [8]. The difficulty with this aroach obviously relates to alying the action of the inverse of the Schur comlement (5) as art of the reconditioning technique. It is usually the case that the Schur comlement is not available, due to its dense structure and that suitable aroximations (or factorizations) are required if we are to stand any chance to make this aroach cometitive. Two successful Schur comlement aroximations have been devised, analyzed and tested for the case of the Picard iteration (ρ = 0). We consider these aroaches below. The first attemt to aroximate the Schur comlement corresonding to the Picard iteration with a nonsymmetric matrix was resented by Elman in ([4]) who suggested the aroximation S (BB t ) (BF (u k )B t )(BB t ). It is shown in ([4]) that the erformance of the block reconditioner (6) with this choice of aroximation has a mild deendence on both ν and n for the Picard iteration. We refer to ([4]) for the derivation and details; we note however that this choice of aroximation is algebraic, i.e., it can be generalized to any roblem of saddle-oint tye and is indeendent of the underlying differential oerator corresonding to F (ρ, u k ). For this reason, we will take our first aroximation of the Schur comlement of the Newton method to be S 1 = (BB t ) (BF (1, u k )B t )(BB t ). The other aroach we consider here is derived in ([9]). In this aroach, the Schur comlement is seen as the discretization of the fundamental solution for the differential oerator which gives rise to K(u k ), the Picard linearization matrix. The resulting aroximation of the Schur comlement is also defined as a factorization of sarse matrices S 2 = M F A where M, F, A are rojections onto the ressure sace of the identity, advection-diffusion and Lalace oerators, resectively (see ([9]) for details. The cost of assembling the reconditioner in this fashion is justified by the imroved erformance which becomes linear in the size of the roblem, although it still exhibits a mild deendence on ν. Unlike the first choice above, this reconditioner deends strongly on what the underlying oerator is: if F (ρ, u k ) changes, the new underlying oerator will lead to a new fundamental solution which will have to be comuted (if it exists) and discretized. Since F (1, u k ) corresonds to an oerator with non-constant coefficients, the fundamental solution is not available to our knowledge. For this reason we have not been able to adat this reconditioner to Newton s method. We therefore analyze and test the above Schur comlement aroximation S 2 (derived for the Picard iteration) on Newton s method. (6) 3
4 2 Eigenvalue analysis In this section we resent an analysis of the above reconditioning techniques for the case where the system matrix results from a mixed finite element discretization of the Navier Stokes equations. These discretizations have the advantage of satisfying certain stability conditions which arise in the formulation of the roblem. More recisely, we assume that the following conditions (known as the Babuška-Brezzi conditions) are satisfied by the Jacobian matrix J(u k ) for all u k and for all ν max max w R n \{0} v R n \{0} min max w R n \{0} v R n \{0} w t J(u k )v v H w H Γ (7a) w t J(u k )v v H w H γ, (7b) where Γ, γ > 0. These conditions guarantee existence and uniqueness of a finite element aroximation ([6, Cha. IV]). Here H is a symmetric ositive-definite matrix which defines the discrete norm of the roblem; this norm deends on the choice of function saces where the solution is sought. For the case of the standard mixed finite-element aroximations of the Navier Stokes equations the norm-matrix H takes the form ( ) A 0 H =, 0 M where A = I d L is a d-dimensional Lalacian matrix and M is a Grammian matrix (mass matrix) assembled on the ressure sace. Let us consider a general block-triangular reconditioner of the form (6) with F (ρ, u k ) and S relaced by the resective aroximations F, Ŝ ( ) F B t P =. (8) 0 Ŝ The reconditioned matrix is then ( F JP ˆF (I F = ˆF )B t Ŝ B ˆF SŜ and it aears that the derivation of eigenvalue bounds for the above system is not a simle task. For this reason we make recourse to the general analysis of block-reconditioners for saddle-oint roblems resented in [11]. We recall here the results of interest. Definition 2.1 H-norm equivalence Non-singular matrices M, N R n n are said to be H-norm equivalent if there exist constants α, β indeendent of n such that for all x R n \ {0} ). We write α Mx H Nx H β. M H N. 4
5 Thus, H-norm equivalence simly means that the H-singular values of MN are bounded. An immediate consequence is that the eigenvalues of the reconditioned system are bounded indeendently of the size of the roblem α < λ(mn ) < β; moreover, H-norm equivalence is an equivalence relation on R n n. We will need the following results roved in [11]. Lemma 2.2 Let (7) hold. Then Moreover, if for all v R n1 n1 H H J. v t F (1, u k )v η v 2 A (9) then the Schur comlement S = BF (1, u k ) B t satisfies S M M. (10) Lemma 2.3 Let (7) hold with H defined as above and let P be defined in (8). Then P H J if F A A, Ŝ M M. The above results simlify our task considerably; in articular, under standard assumtions on the matrix J, it is sufficient to erform our analysis on the equivalence of the diagonal blocks of H and P. Let us consider first a norm-equivalent aroximation F. Since F (0, u k ) is the discretization of a (ositive-definite) advection-diffusion oerator it satisfies (see, e.g., [15]) v t F (0, u k )v η v 2 A v R n1 n1. (11) This bound imlies that F (0, u k ) satisfies a min-max bound of tye (7b) with resect to the A-norm; since F (0, u k ) also satisfies a max-max bound with resect to the same norm, we can aly Lemma 2.2 to deduce that F (0, u k ) A A. It can be shown that there exists a value ν 0 such that a bound of tye (9) holds also for F (1, u k ) [6,. 300] and thus F (1, u k ) A A. We remark here that the value of ν 0 is not known a riori and that it deends on alication. We will investigate this issue in the numerics section. Remark 2.1 In general, even the weaker bound min max w R n 1 \{0} v R n 1 \{0} w t F (1, u k )v v A w A η (12) is not imlied by the stability conditions (7). In fact, the Babuška-Brezzi conditions imly that the above min-max bound only holds over kerb [1]. 5
6 Since the inverse of the matrix F (1, u k ) can be quite exensive to comute in ractice we consider below two block aroximations. For this urose, let us look at the structure of F (1, u k ) for two-dimensional roblems (cf. (4)) ( ) ( ) F (1, u k F (u ) = k ) 0 M11 (u 0 F (u k + k ) M 12 (u k ) ) M 21 (u k ) M 22 (u k. ) Simle aroximations are given by ( F T (u k F (u ) = k ) + M 11 (u k ) M 12 (u k ) 0 F (u k ) + M 22 (u k ) ). (13) and F D (u k ) = ( F (u k ) + M 11 (u k ) 0 0 F (u k ) + M 22 (u k ) ). (14) Given their structure, F T (u k ), F D (u k ) satisfy (7a) if F (1, u k ) does. If (9) holds for F (1, u k ) then F D (u k ) satisfies (9) also and thus F D (u k ) A A. (15) If moreover M 12 (u k ) < η, then one can show that v t F T (u k )v 1 2 η v 2 A ; thus F T (u k ) A A. (16) We note that for a quasi-uniform discretization (see below) M 12 (u k ) = O(h 2 ) so that the above restriction is satisfied if η > ch 2. However, η = O(ν) and ν > ν 0 is a stronger restriction on M 12 (u k ). Henceforth, we restrict our attention to the regime ν > ν 0, for which F (1, u k ), F T (u k ), F D (u k ) are all norm equivalent to A. This equivalence is useful in the analysis we resent below. In the following we assume a quasi-uniform discretization of our hysical domain, i.e., we assume that we have a subdivision of our comutational domain in which all the simlices do not exceed a diameter h. Moreover, the following bounds can be shown to hold for all v R n1, q R n2 \ kerb t and Ω R 2 c1h 2 v 2 v t Av c2 v 2, (17a) c 3 h 4 q 2 q t BB t q c 4 h 2 q 2, c 5 h 2 q 2 q t M q c 6 h 2 q 2. (17b) (17c) The following results rovide eigenvalue bounds for the system matrix reconditioned in the fashion described in the revious section. Theorem 2.4 Let (7) hold and let ν > ν 0 be a viscosity arameter such that F (1, u k ) satisfies (9). Let ( ) F (1, u P i = k ) B t (18) 0 S i 6
7 where S i, i = 1, 2 are defined as above. Then there exist constants C i, i = 1,..., 4 indeendent of n such that C 1 λ(jp 1 ) C 2 h 2 C 3 λ(jp 2 ) C 4. Proof We first rove the bounds for P 2. Since by hyothesis (9) holds, by the above discussion F (1, u k ) A A; in articular, (7a) imlies F (1, u k )A A = F (1, u k )v max A v n 1\0 Av A Γ. (19) According to Lemma 2.3 we only have to check that S 2 M M. This is equivalent to showing that there exist constants α, β such that α AF x M x M β. On the other hand the above inequality holds with resect to the l 2- norm ([12], see also [10] ) and the required result follows from the sectral equivalence between the mass matrix M and the identity [16]. The bound for P 1 follows in a similar way, the only difference being that we have to show that β is of order h 2 in this case. We have M S 1 M = M 1/2 S 1 M 1/2 c 6h 2 (BB t ) BF (1, u k )B t (BB t ) c 6h 2 A /2 F (1, u k )A /2 (BB t ) BAB t (BB t ) c 6Γh 2 (BB t ) BAB t (BB t ) c 6Γh 2 A (BB t ) c 6c 2c 3 Γh 2 where we used (19) and the fact that A /2 F (1, u k )A /2 = F (1, u k )A A. For the lower bound we note that S 1M M = M /2 S 1M /2 1 σ min(m 1/2 S 1 M 1/2 ). Let S 0 = (BB t )(BAB t ) (BB t ) and let n 2 = n 2 \ kerb t. We have σ min(m 1/2 S 1 M 1/2 ) min q n 2 q t M 1/2 S 1 M 1/2 q t q r t S 1 = min r r t S 0 r r n 2 r t S 0 r r t (BA B t ) r z t F (1, u k )z min z n 2 z t Az η 1 c q min r n 2 r t (BA B t ) r r t M r r t S 0 r r t (BA B t ) r r t (BA B t ) r r t M r 7
8 where we used Lemma (2.2), the fact that F (1, u k ) satisfies (9) and the result roved in ([4]) that the smallest eigenvalue of is greater or equal to 1. (BA B t )x = λs 0x x n 2 Let Pi T, P i D denote reconditioners of tye (18) with F (1, u k ) relaced by F T (u k ), F D (u k ) resectively. The following result shows that the bounds of Theorem 2.4 hold also for these aroximations. Corollary 2.5 Let the conditions of Theorem 2.4 hold and let Pi denote either of Pi T, P i D. Then there exist constants Ci, i = 1,..., 4 indeendent of n such that C1 λ(jp1 ) C2 h 2 C3 λ(jp 2 ) C4. Proof The roof follows similarly by using the norm-equivalences (15), (16). Remark 2.2 The above analysis relies on relation (9). Under this assumtion, the analysis for the reconditioned matrix JP2 is equivalent to that erformed in [10] for reconditioning the Picard iteration. However, the general analysis in [11] allowed us to simlify that analysis and also to rovide for the first time analytic bounds for the eigenvalues of JP1. Remark 2.3 The above result exhibits only the deendence of the eigenvalues on the mesh arameter. Since by hyothesis, the analysis is restricted to a sufficiently large ν, we will not ursue the deendence on ν here. For the Picard iteration, analyses and exeriments on ν deendence are contained in [4], [5], [10]. Numerical exeriments indicate that the bounds on the eigenvalues are tight. However, it is well-known that the eigenvalues alone do not always describe convergence of nonsymmetric iterative solvers. This is also the case here: although the largest modulus eigenvalues of the reconditioned systems behave as described theoretically, they actually belong to a small grou of outliers. We found exerimentally that most of the eigenvalues are clustered in a region of the lane close to 1, with a small number of eigenvalues migrating from the cluster as the arameters ν, h are reduced. A general qualitative model for convergence of GMRES when clustering occurs can be found in [2]. We recall here their result for one cluster, since it aears to be ertinent to the distribution of eigenvalues resulting from our reconditioning techniques. Given a cluster of eigenvalues C = {λ : λ c < ρ c } to which there corresond M = n C outliers, whose maximum distance δ away from C is λ j z δ = max max z c =ρ c 1 j M λ j 8
9 h(n) maxre(λ) maxim(λ) n o (2, 1) n o (3, 2) n o (4, 3) 1/16(289) /32(1089) /64(4096) Table 1: Sectral information for JP 1 ; n o (c, r) is the ercentage of eigenvalues outside the circle of radius r centered at c. the GMRES residual at ste k + d satisfies r k+d Cρ k r 0 where d is the degree of the minimum olynomial associated with the outliers and C ρδ d is indeendent of k. A secific numerical study of this henomenon for reconditioning with P 2 is contained in ([5]). In this case the location of the eigenvalues does not deend on the mesh-size and the migration of eigenvalues from the cluster occurs with reducing ν. This leads to a delay in convergence deendent on the number of eigenvalues outside the cluster. A very similar behaviour is observed for fixed ν in the case of reconditioning with P 1. As described by Thm 2.4, the maximum modulus eigenvalue of the reconditioned system is bounded from above by a factor of order O(h 2 ). This is in fact a descritive bound. Table 1 exhibits this deendence for the case ν = 1/10 for the driven cavity test roblem (see next section). However, convergence is not entirely described by the behaviour of the largest modulus eigenvalues indeed, the sectrum is mostly clustered in a small region of the comlex lane with only a handful of outliers exhibiting the O(h 2 ) behaviour. Thus, there aears to occur a clustering similar to the case described above where we recondition with P 2 and where convergence is described through a delay deendent on the (increasing) number of outliers. Table 1 resents also the ercentage of (outlying) eigenvalues n o (c, r) outside the circle of radius r centered at c. We see indeed that the number of eigenvalues outside each disc close to the unity, though in small roortion, is increasing with n. This aears to match the behaviour exhibited in the tests below, where a slight deendence on the size of the roblem is noticed, though never a linear one. 3 Numerical exeriments In this section we resent numerical exeriments obtained for a standard test roblem, the driven cavity flow. This involves solving the Navier Stokes equations inside the unit square with Dirichlet boundary conditions equal to zero everywhere, excet for the boundary y = 1 where the horizontal velocity has a ositive rescribed rofile. Although our analysis considered the exact Newton s method, in ractice it is of considerable advantage to use an inexact version for reasons of comutational efficiency. The algorithm we use is a Newton-GMRES method which involves solving (2) using GMRES with reconditioning and 9
10 the following stoing criterion suggested in [3]: writing (1) as F(w) = 0, where w = (u, ), at each Newton ste i we sto after k iterations of GMRES if the residual r k satisfies r k / F(w i ) c F(w i ) q, (20) where the choice c = 10 2, q = 1/4 does not affect the number of nonlinear iterations in our tests. Moreover, the order of convergence is only marginally affected, as we show below. We reort below the erformance of reconditioners P i, Pi T, P i D, i = 1, 2. The tests were erformed on discretizations resulting from three successively refined meshes and a range of diffusion arameters ν. The value ν 0 in Thm 2.4 deends on the roblem; for our examle we found exerimentally that ν 0 1/80. As the results below demonstrate, the regime ν ν 0 is also the regime where our reconditioners erform best. The finite element discretization is the so-called Q2 Q1 discretization, using a regular mesh of rectangles with finite element bases of quadratic and linear iecewise olynomials for the velocity and ressure, resectively. The initial guess for the Newton iteration was the zero vector, excet for the range ν 1/640 for which three stes of the Picard iteration were necessary to rovide a initial iterate which ensured convergence of Newton s method. We note that this will lead to an aarent imrovement of erformance in this regime, although not to a qualitative change. We stoed the Newton (outer) iteration when F(w i ) / F(w 0 ) The comutational cost of our reconditioners is dominated by the inversion of the (1,1)-block F (1, u k ). For the Q2 Q1 discretization the cost of storing the Schur comlement aroximations S i is aroximately 25n 2 (BB t ) +45n 2 (BF B t ) for S 1 and 3 9n 2 (A, F, M ) for S 2. The cost of alying these reconditioners will of course deend on the choice of imlementation. In the exeriments below we imlemented them exactly; however, iterative techniques are certainly ossible [4], [9] in which case the choice S 2 becomes clearly more cometitive. The results for the exact case, when reconditioners P 1, P 2 are emloyed are shown in Tables 2, 3. We see indeed that while the erformance of P 2 is mesh-indeendent, the erformance of P 1 exhibits a sub-linear deendence on the size of the roblem. This we believe to be exlained by the ν = 1/10 1/20 1/40 1/80 1/160 1/320 1/640 1/1280 n = 2, , , Table 2: Average number of GMRES iterations er Newton ste for P 1. ν = 1/10 1/20 1/40 1/80 1/160 1/320 1/640 1/1280 n = 2, , , Table 3: Average number of GMRES iterations er Newton ste for P 2. 10
11 ν = 1/10 1/20 1/40 1/80 1/160 1/320 1/640 1/1280 n = 2, , , Table 4: Average number of GMRES iterations er Newton ste for P T 1. ν = 1/10 1/20 1/40 1/80 1/160 1/320 1/640 1/1280 n = 2, , , Table 5: Average number of GMRES iterations er Newton ste for P T 2. distribution of eigenvalues which was exhibited in Table 1, which showed an increasing number of outliers with decreasing mesh-size. As for the viscosity arameter ν, reconditioner P 1 aears to exhibit a milder deendence than P 2, which makes it comarable to P 1 for high-reynolds number flows, desite its mesh-deendence. The same behaviour is exhibited by the reconditioners Pi T, P i D, i = 1, 2. In articular, the erformance remains mesh-indeendent in the case of P2 T, P2 D as shown in Tables 5, 7; as a consequence, the number of extra iterations comared to the exact case dislayed in Table 3 is indeendent of the mesh. This is not the case for P1 T, P 1 D as shown in Tables 4, 6, for which this difference (cf. Table 2) grows with the mesh arameter, although the number of iterations still grows sub-linearly with resect to the size of the roblem. From a ractical oint of view it is imortant to study the deterioration in erformance when the exact reconditioners P i are relaced with Pi T, P i D, as solving systems with F can be quite costly. As the above numerics show, this deterioration is negligible for large ν and quite accetable for smaller ν. To exemlify this we relaced the exact solution of systems ν = 1/10 1/20 1/40 1/80 1/160 1/320 1/640 1/1280 n = 2, , , Table 6: Average number of GMRES iterations er Newton ste for P D 1. ν = 1/10 1/20 1/40 1/80 1/160 1/320 1/640 1/1280 n = 2, , , Table 7: Average number of GMRES iterations er Newton ste for P D 2. 11
12 ν = 1/10 1/20 1/40 1/80 1/160 1/320 1/640 1/1280 F T F D Table 8: Average number of iterations er outer-gmres ste for inner- GMRES solver using reconditioners F T, F D. involving F (1, u k ) with a GMRES iteration with reconditioners F T, F D. The stoing criterion was the reduction of the relative 2-norm residual by a factor of 10 6 which insured that this inner iteration did not affect the outer iteration count dislayed in Tables 2, 3. The number of iterations was always mesh-indeendent. We resent these results in Table 8. We see indeed that the modified reconditioners Pi T, P i D (which require just one solve with F T and F D resectively) outerform by far P i when the inner-gmres iteration is used to aroximate F. Thus, unless a fast solution algorithm is available for F (1, u k ), the ideal erformance resented in Tables 2, 3 may rove too exensive to achieve. We end this section with a note on the convergence rate of Newton s method observed for our inexact Newton-GMRES algorithm. The choice of this algorithm and in articular of stoing criterion (20) was justified by both ractical and theoretical considerations [3]. In Table 9 we dislay the average convergence rate observed in our tests. More recisely, we tabulated the mean of the ratio q k of logs of consecutive residuals q k = log F(wk ) log F(w k ), for the cases q = 1/4 and q = 1 in (20) and for the exact case (q = ). The choice of reconditioner does not affect the convergence rate dislayed below. We see indeed that the choice of stoing criterion (i.e., q in (20)) does not affect significantly the rate of convergence of the method excet ossibly for the low Reynolds number (large ν) regime. In rincile, one could recover this convergence rate by a stricter criterion (q = 1, say), at the cost of more iterations. Given our choice of tolerance for the outer nonlinear iteration (10 6 ), q = 1/4 was sufficient to kee the same number of Newton stes as in the exact case. A stricter tolerance may require a larger value of q. However, for high Reynolds number flows a more economical stoing criterion is sufficient due to the sub-otimal convergence of Newton s method. ν = 1/10 1/20 1/40 1/80 1/160 1/320 q = 1/ q = q = Table 9: Convergence rate of Newton s method: exact and inexact with q = 1/4 and q = 1; n = 9,
13 4 Summary The urose of this work was to analyze and comare block reconditioners for the system matrix arising from a Newton iteration for the Navier Stokes equations. The reconditioners are saddle-oint reconditioners using an algebraic and oerator based aroximation of the Schur comlement. They were introduced originally for the case of the Picard iteration. In each case we give analytic bounds for the sectrum of the reconditioned system under standard stability assumtions for the finite element discretization. The bounds require a further assumtion of ositivedefiniteness of the (1, 1)-block of the matrix, which holds for large viscosity arameters. The numerical exeriments however show no qualitative change even for smaller values of viscosity. In each case the imlementation is modular we only require the action of inverses of sub-blocks which have to be constructed algebraically or assembled together with the system matrix. The erformance is very similar to the case of the Picard iteration which has been recently reorted. The oerator-based reconditioner aears to be more robust and less costly than the algebraic reconditioner and esecially the indeendence of the size of the roblem is an attractive feature. However, otimality with resect to the Reynolds number is still to be achieved. References [1] F. Brezzi and M. Fortin, Mixed and hybrid finite element methods, Sringer-Verlag, New York, [2] S. L. Cambell, I. C. F. Isen, C. T. Kelley, and C. D. Meyer, GMRES and the minimal olynomial, BIT, 36 (1996), [3] R. S. Dembo, S. C. Eisenstat, and T. Steihaug, Inexact Newton methods, SIAM J. Numer. Anal., 19 (1982), [4] H. C. Elman, Preconditioning for the steady-state Navier-Stokes equations with low viscosity, SIAM J. Sci. Com., 20 (1999), [5] H. C. Elman, D. J. Silvester, and A. J. Wathen, Performance and analysis of saddle oint reconditioners for the discrete steadystate Navier-Stokes equations, Numer. Math., 90 (2001), [6] V. Girault and P. A. Raviart, Finite Element Methods for Navier Stokes Equations, Sringer-Verlag, [7] P. M. Gresho and R. L. Sani, Incomressible flow and the finite element method: isothermal laminar flow, Wiley, Chichester, [8] I. C. F. Isen, A note on reconditioning non-symmetric matrices, SIAM J. Sci. Comut., 23(2001), [9] D. Kay, D. Loghin, and A. J. Wathen, A reconditioner for the Steady-State Navier-Stokes equations, SIAM J. Sci. Comut., 24 (2002),
14 [10] D. Loghin, Analysis of reconditioned Picard iterations for the Navier-Stokes equations. Submitted to Numer. Math, [11] D. Loghin and A. J. Wathen, Analysis of block reconditioners for saddle-oint roblems, Tech. Re. 13, Oxford University Comuting Laboratory, Submitted to SISC. [12] T. A. Manteuffel and S. V. Parter, Preconditioning and boundary conditions, SIAM J. Num. Anal., 27 (1989), [13] M. F. Murhy, G. H. Golub, and A. J. Wathen, A note on reconditioning for indefinite linear systems, SIAM J. Sci. Com., 21 (2000), [14] A. Quarteroni and A. Valli, Numerical Aroximation of Partial Differential Equations, Sringer-Verlag, [15] H. G. Roos, M. Stynes, and L. Tobiska, Numerical Methods for Singularly Perturbed Differential Equations, Sringer Series in Comutational Mathematics, Sringer, [16] A. J. Wathen, Realistic eigenvalue bounds for the Galerkin mass matrix, IMA J. Numer. Anal., (1987),
PRECONDITIONING TECHNIQUES FOR NEWTON S METHOD FOR THE INCOMPRESSIBLE NAVIER STOKES EQUATIONS
BIT Numerical Mathematics 43: 961 974, 2003. 2003 Kluwer Academic Publishers. Printed in the Netherlands. 961 PRECONDITIONING TECHNIQUES FOR NEWTON S METHOD FOR THE INCOMPRESSIBLE NAVIER STOKES EQUATIONS
More informationON THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS AND COMMUTATOR ARGUMENTS FOR SOLVING STOKES CONTROL PROBLEMS
Electronic Transactions on Numerical Analysis. Volume 44,. 53 72, 25. Coyright c 25,. ISSN 68 963. ETNA ON THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS AND COMMUTATOR ARGUMENTS FOR SOLVING STOKES
More informationON THE LEAST SIGNIFICANT p ADIC DIGITS OF CERTAIN LUCAS NUMBERS
#A13 INTEGERS 14 (014) ON THE LEAST SIGNIFICANT ADIC DIGITS OF CERTAIN LUCAS NUMBERS Tamás Lengyel Deartment of Mathematics, Occidental College, Los Angeles, California lengyel@oxy.edu Received: 6/13/13,
More informationPaper C Exact Volume Balance Versus Exact Mass Balance in Compositional Reservoir Simulation
Paer C Exact Volume Balance Versus Exact Mass Balance in Comositional Reservoir Simulation Submitted to Comutational Geosciences, December 2005. Exact Volume Balance Versus Exact Mass Balance in Comositional
More information4. Score normalization technical details We now discuss the technical details of the score normalization method.
SMT SCORING SYSTEM This document describes the scoring system for the Stanford Math Tournament We begin by giving an overview of the changes to scoring and a non-technical descrition of the scoring rules
More informationHENSEL S LEMMA KEITH CONRAD
HENSEL S LEMMA KEITH CONRAD 1. Introduction In the -adic integers, congruences are aroximations: for a and b in Z, a b mod n is the same as a b 1/ n. Turning information modulo one ower of into similar
More informationLocation of solutions for quasi-linear elliptic equations with general gradient dependence
Electronic Journal of Qualitative Theory of Differential Equations 217, No. 87, 1 1; htts://doi.org/1.14232/ejqtde.217.1.87 www.math.u-szeged.hu/ejqtde/ Location of solutions for quasi-linear ellitic equations
More informationFeedback-error control
Chater 4 Feedback-error control 4.1 Introduction This chater exlains the feedback-error (FBE) control scheme originally described by Kawato [, 87, 8]. FBE is a widely used neural network based controller
More information216 S. Chandrasearan and I.C.F. Isen Our results dier from those of Sun [14] in two asects: we assume that comuted eigenvalues or singular values are
Numer. Math. 68: 215{223 (1994) Numerische Mathemati c Sringer-Verlag 1994 Electronic Edition Bacward errors for eigenvalue and singular value decomositions? S. Chandrasearan??, I.C.F. Isen??? Deartment
More information#A64 INTEGERS 18 (2018) APPLYING MODULAR ARITHMETIC TO DIOPHANTINE EQUATIONS
#A64 INTEGERS 18 (2018) APPLYING MODULAR ARITHMETIC TO DIOPHANTINE EQUATIONS Ramy F. Taki ElDin Physics and Engineering Mathematics Deartment, Faculty of Engineering, Ain Shams University, Cairo, Egyt
More informationAn Analysis of Reliable Classifiers through ROC Isometrics
An Analysis of Reliable Classifiers through ROC Isometrics Stijn Vanderlooy s.vanderlooy@cs.unimaas.nl Ida G. Srinkhuizen-Kuyer kuyer@cs.unimaas.nl Evgueni N. Smirnov smirnov@cs.unimaas.nl MICC-IKAT, Universiteit
More informationProbability Estimates for Multi-class Classification by Pairwise Coupling
Probability Estimates for Multi-class Classification by Pairwise Couling Ting-Fan Wu Chih-Jen Lin Deartment of Comuter Science National Taiwan University Taiei 06, Taiwan Ruby C. Weng Deartment of Statistics
More informationA note on variational representation for singular values of matrix
Alied Mathematics and Comutation 43 (2003) 559 563 www.elsevier.com/locate/amc A note on variational reresentation for singular values of matrix Zhi-Hao Cao *, Li-Hong Feng Deartment of Mathematics and
More information#A37 INTEGERS 15 (2015) NOTE ON A RESULT OF CHUNG ON WEIL TYPE SUMS
#A37 INTEGERS 15 (2015) NOTE ON A RESULT OF CHUNG ON WEIL TYPE SUMS Norbert Hegyvári ELTE TTK, Eötvös University, Institute of Mathematics, Budaest, Hungary hegyvari@elte.hu François Hennecart Université
More informationIntrinsic Approximation on Cantor-like Sets, a Problem of Mahler
Intrinsic Aroximation on Cantor-like Sets, a Problem of Mahler Ryan Broderick, Lior Fishman, Asaf Reich and Barak Weiss July 200 Abstract In 984, Kurt Mahler osed the following fundamental question: How
More informationAdaptive preconditioners for nonlinear systems of equations
Adaptive preconditioners for nonlinear systems of equations L. Loghin D. Ruiz A. Touhami CERFACS Technical Report TR/PA/04/42 Also ENSEEIHT-IRIT Technical Report RT/TLSE/04/02 Abstract The use of preconditioned
More informationLinear diophantine equations for discrete tomography
Journal of X-Ray Science and Technology 10 001 59 66 59 IOS Press Linear diohantine euations for discrete tomograhy Yangbo Ye a,gewang b and Jiehua Zhu a a Deartment of Mathematics, The University of Iowa,
More informationUnderstanding DPMFoam/MPPICFoam
Understanding DPMFoam/MPPICFoam Jeroen Hofman March 18, 2015 In this document I intend to clarify the flow solver and at a later stage, the article-fluid and article-article interaction forces as imlemented
More informationRecursive Estimation of the Preisach Density function for a Smart Actuator
Recursive Estimation of the Preisach Density function for a Smart Actuator Ram V. Iyer Deartment of Mathematics and Statistics, Texas Tech University, Lubbock, TX 7949-142. ABSTRACT The Preisach oerator
More informationImproving AOR Method for a Class of Two-by-Two Linear Systems
Alied Mathematics 2 2 236-24 doi:4236/am22226 Published Online February 2 (htt://scirporg/journal/am) Imroving AOR Method for a Class of To-by-To Linear Systems Abstract Cuixia Li Shiliang Wu 2 College
More informationPositive decomposition of transfer functions with multiple poles
Positive decomosition of transfer functions with multile oles Béla Nagy 1, Máté Matolcsi 2, and Márta Szilvási 1 Deartment of Analysis, Technical University of Budaest (BME), H-1111, Budaest, Egry J. u.
More informationLORENZO BRANDOLESE AND MARIA E. SCHONBEK
LARGE TIME DECAY AND GROWTH FOR SOLUTIONS OF A VISCOUS BOUSSINESQ SYSTEM LORENZO BRANDOLESE AND MARIA E. SCHONBEK Abstract. In this aer we analyze the decay and the growth for large time of weak and strong
More informationFast solvers for steady incompressible flow
ICFD 25 p.1/21 Fast solvers for steady incompressible flow Andy Wathen Oxford University wathen@comlab.ox.ac.uk http://web.comlab.ox.ac.uk/~wathen/ Joint work with: Howard Elman (University of Maryland,
More informationSUPER-GEOMETRIC CONVERGENCE OF A SPECTRAL ELEMENT METHOD FOR EIGENVALUE PROBLEMS WITH JUMP COEFFICIENTS *
Journal of Comutational Mathematics Vol.8, No.,, 48 48. htt://www.global-sci.org/jcm doi:.48/jcm.9.-m6 SUPER-GEOMETRIC CONVERGENCE OF A SPECTRAL ELEMENT METHOD FOR EIGENVALUE PROBLEMS WITH JUMP COEFFICIENTS
More informationElementary Analysis in Q p
Elementary Analysis in Q Hannah Hutter, May Szedlák, Phili Wirth November 17, 2011 This reort follows very closely the book of Svetlana Katok 1. 1 Sequences and Series In this section we will see some
More informationCFL Conditions for Runge-Kutta Discontinuous Galerkin Methods on Triangular Grids
CFL Conditions for Runge-Kutta Discontinuous Galerkin Methods on Triangular Grids T. Toulorge a,, W. Desmet a a K.U. Leuven, Det. of Mechanical Engineering, Celestijnenlaan 3, B-31 Heverlee, Belgium Abstract
More informationNotes on duality in second order and -order cone otimization E. D. Andersen Λ, C. Roos y, and T. Terlaky z Aril 6, 000 Abstract Recently, the so-calle
McMaster University Advanced Otimization Laboratory Title: Notes on duality in second order and -order cone otimization Author: Erling D. Andersen, Cornelis Roos and Tamás Terlaky AdvOl-Reort No. 000/8
More informationPositivity, local smoothing and Harnack inequalities for very fast diffusion equations
Positivity, local smoothing and Harnack inequalities for very fast diffusion equations Dedicated to Luis Caffarelli for his ucoming 60 th birthday Matteo Bonforte a, b and Juan Luis Vázquez a, c Abstract
More informationUncorrelated Multilinear Principal Component Analysis for Unsupervised Multilinear Subspace Learning
TNN-2009-P-1186.R2 1 Uncorrelated Multilinear Princial Comonent Analysis for Unsuervised Multilinear Subsace Learning Haiing Lu, K. N. Plataniotis and A. N. Venetsanooulos The Edward S. Rogers Sr. Deartment
More informationGOOD MODELS FOR CUBIC SURFACES. 1. Introduction
GOOD MODELS FOR CUBIC SURFACES ANDREAS-STEPHAN ELSENHANS Abstract. This article describes an algorithm for finding a model of a hyersurface with small coefficients. It is shown that the aroach works in
More informationSome results of convex programming complexity
2012c12 $ Ê Æ Æ 116ò 14Ï Dec., 2012 Oerations Research Transactions Vol.16 No.4 Some results of convex rogramming comlexity LOU Ye 1,2 GAO Yuetian 1 Abstract Recently a number of aers were written that
More informationA Bound on the Error of Cross Validation Using the Approximation and Estimation Rates, with Consequences for the Training-Test Split
A Bound on the Error of Cross Validation Using the Aroximation and Estimation Rates, with Consequences for the Training-Test Slit Michael Kearns AT&T Bell Laboratories Murray Hill, NJ 7974 mkearns@research.att.com
More informationEfficient algorithms for the smallest enclosing ball problem
Efficient algorithms for the smallest enclosing ball roblem Guanglu Zhou, Kim-Chuan Toh, Jie Sun November 27, 2002; Revised August 4, 2003 Abstract. Consider the roblem of comuting the smallest enclosing
More informationApproximating min-max k-clustering
Aroximating min-max k-clustering Asaf Levin July 24, 2007 Abstract We consider the roblems of set artitioning into k clusters with minimum total cost and minimum of the maximum cost of a cluster. The cost
More informationUniform Law on the Unit Sphere of a Banach Space
Uniform Law on the Unit Shere of a Banach Sace by Bernard Beauzamy Société de Calcul Mathématique SA Faubourg Saint Honoré 75008 Paris France Setember 008 Abstract We investigate the construction of a
More informationReal Analysis 1 Fall Homework 3. a n.
eal Analysis Fall 06 Homework 3. Let and consider the measure sace N, P, µ, where µ is counting measure. That is, if N, then µ equals the number of elements in if is finite; µ = otherwise. One usually
More informationc Copyright by Helen J. Elwood December, 2011
c Coyright by Helen J. Elwood December, 2011 CONSTRUCTING COMPLEX EQUIANGULAR PARSEVAL FRAMES A Dissertation Presented to the Faculty of the Deartment of Mathematics University of Houston In Partial Fulfillment
More informationMultiplicity of weak solutions for a class of nonuniformly elliptic equations of p-laplacian type
Nonlinear Analysis 7 29 536 546 www.elsevier.com/locate/na Multilicity of weak solutions for a class of nonuniformly ellitic equations of -Lalacian tye Hoang Quoc Toan, Quô c-anh Ngô Deartment of Mathematics,
More informationCombining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)
Combining Logistic Regression with Kriging for Maing the Risk of Occurrence of Unexloded Ordnance (UXO) H. Saito (), P. Goovaerts (), S. A. McKenna (2) Environmental and Water Resources Engineering, Deartment
More informationSYMPLECTIC STRUCTURES: AT THE INTERFACE OF ANALYSIS, GEOMETRY, AND TOPOLOGY
SYMPLECTIC STRUCTURES: AT THE INTERFACE OF ANALYSIS, GEOMETRY, AND TOPOLOGY FEDERICA PASQUOTTO 1. Descrition of the roosed research 1.1. Introduction. Symlectic structures made their first aearance in
More informationQuantitative estimates of propagation of chaos for stochastic systems with W 1, kernels
oname manuscrit o. will be inserted by the editor) Quantitative estimates of roagation of chaos for stochastic systems with W, kernels Pierre-Emmanuel Jabin Zhenfu Wang Received: date / Acceted: date Abstract
More informationEstimation of the large covariance matrix with two-step monotone missing data
Estimation of the large covariance matrix with two-ste monotone missing data Masashi Hyodo, Nobumichi Shutoh 2, Takashi Seo, and Tatjana Pavlenko 3 Deartment of Mathematical Information Science, Tokyo
More informationarxiv: v2 [stat.me] 3 Nov 2014
onarametric Stein-tye Shrinkage Covariance Matrix Estimators in High-Dimensional Settings Anestis Touloumis Cancer Research UK Cambridge Institute University of Cambridge Cambridge CB2 0RE, U.K. Anestis.Touloumis@cruk.cam.ac.uk
More informationPositive Definite Uncertain Homogeneous Matrix Polynomials: Analysis and Application
BULGARIA ACADEMY OF SCIECES CYBEREICS AD IFORMAIO ECHOLOGIES Volume 9 o 3 Sofia 009 Positive Definite Uncertain Homogeneous Matrix Polynomials: Analysis and Alication Svetoslav Savov Institute of Information
More informationApplications to stochastic PDE
15 Alications to stochastic PE In this final lecture we resent some alications of the theory develoed in this course to stochastic artial differential equations. We concentrate on two secific examles:
More informationNONLINEAR OPTIMIZATION WITH CONVEX CONSTRAINTS. The Goldstein-Levitin-Polyak algorithm
- (23) NLP - NONLINEAR OPTIMIZATION WITH CONVEX CONSTRAINTS The Goldstein-Levitin-Polya algorithm We consider an algorithm for solving the otimization roblem under convex constraints. Although the convexity
More informationOn the Properties for Iteration of a Compact Operator with Unstructured Perturbation
On the Proerties for Iteration of a Comact Oerator with Unstructured Perturbation Jarmo Malinen Helsinki University of Technology Institute of Mathematics Research Reorts A360 (1996) 1 2 Jarmo Malinen:
More informationA Review of Preconditioning Techniques for Steady Incompressible Flow
Zeist 2009 p. 1/43 A Review of Preconditioning Techniques for Steady Incompressible Flow David Silvester School of Mathematics University of Manchester Zeist 2009 p. 2/43 PDEs Review : 1984 2005 Update
More informationAn Estimate For Heilbronn s Exponential Sum
An Estimate For Heilbronn s Exonential Sum D.R. Heath-Brown Magdalen College, Oxford For Heini Halberstam, on his retirement Let be a rime, and set e(x) = ex(2πix). Heilbronn s exonential sum is defined
More information5. PRESSURE AND VELOCITY SPRING Each component of momentum satisfies its own scalar-transport equation. For one cell:
5. PRESSURE AND VELOCITY SPRING 2019 5.1 The momentum equation 5.2 Pressure-velocity couling 5.3 Pressure-correction methods Summary References Examles 5.1 The Momentum Equation Each comonent of momentum
More informationk- price auctions and Combination-auctions
k- rice auctions and Combination-auctions Martin Mihelich Yan Shu Walnut Algorithms March 6, 219 arxiv:181.3494v3 [q-fin.mf] 5 Mar 219 Abstract We rovide for the first time an exact analytical solution
More informationE-companion to A risk- and ambiguity-averse extension of the max-min newsvendor order formula
e-comanion to Han Du and Zuluaga: Etension of Scarf s ma-min order formula ec E-comanion to A risk- and ambiguity-averse etension of the ma-min newsvendor order formula Qiaoming Han School of Mathematics
More informationarxiv:cond-mat/ v2 25 Sep 2002
Energy fluctuations at the multicritical oint in two-dimensional sin glasses arxiv:cond-mat/0207694 v2 25 Se 2002 1. Introduction Hidetoshi Nishimori, Cyril Falvo and Yukiyasu Ozeki Deartment of Physics,
More informationApplied Mathematics and Computation
Alied Mathematics and Comutation 217 (2010) 1887 1895 Contents lists available at ScienceDirect Alied Mathematics and Comutation journal homeage: www.elsevier.com/locate/amc Derivative free two-oint methods
More informationBoundary problems for fractional Laplacians and other mu-transmission operators
Boundary roblems for fractional Lalacians and other mu-transmission oerators Gerd Grubb Coenhagen University Geometry and Analysis Seminar June 20, 2014 Introduction Consider P a equal to ( ) a or to A
More informationarxiv: v1 [math.na] 5 Dec 2018
arxiv:8.0957v [math.na] 5 Dec 08 Residual-tye a osteriori estimators for a singularly erturbed reaction-diffusion variational inequality - reliability, efficiency and robustness Mirjam Walloth Fachbereich
More informationCFD AS A DESIGN TOOL FOR FLUID POWER COMPONENTS
CFD AS A DESIGN TOOL FOR FLUID POWER COMPONENTS M. BORGHI - M. MILANI Diartimento di Scienze dell Ingegneria Università degli Studi di Modena Via Cami, 213/b 41100 Modena E-mail: borghi@omero.dsi.unimo.it
More informationThe Fekete Szegő theorem with splitting conditions: Part I
ACTA ARITHMETICA XCIII.2 (2000) The Fekete Szegő theorem with slitting conditions: Part I by Robert Rumely (Athens, GA) A classical theorem of Fekete and Szegő [4] says that if E is a comact set in the
More informationFE FORMULATIONS FOR PLASTICITY
G These slides are designed based on the book: Finite Elements in Plasticity Theory and Practice, D.R.J. Owen and E. Hinton, 1970, Pineridge Press Ltd., Swansea, UK. 1 Course Content: A INTRODUCTION AND
More informationLilian Markenzon 1, Nair Maria Maia de Abreu 2* and Luciana Lee 3
Pesquisa Oeracional (2013) 33(1): 123-132 2013 Brazilian Oerations Research Society Printed version ISSN 0101-7438 / Online version ISSN 1678-5142 www.scielo.br/oe SOME RESULTS ABOUT THE CONNECTIVITY OF
More informationAI*IA 2003 Fusion of Multiple Pattern Classifiers PART III
AI*IA 23 Fusion of Multile Pattern Classifiers PART III AI*IA 23 Tutorial on Fusion of Multile Pattern Classifiers by F. Roli 49 Methods for fusing multile classifiers Methods for fusing multile classifiers
More informationMATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK
Comuter Modelling and ew Technologies, 5, Vol.9, o., 3-39 Transort and Telecommunication Institute, Lomonosov, LV-9, Riga, Latvia MATHEMATICAL MODELLIG OF THE WIRELESS COMMUICATIO ETWORK M. KOPEETSK Deartment
More informationMODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL
Technical Sciences and Alied Mathematics MODELING THE RELIABILITY OF CISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL Cezar VASILESCU Regional Deartment of Defense Resources Management
More informationCERIAS Tech Report The period of the Bell numbers modulo a prime by Peter Montgomery, Sangil Nahm, Samuel Wagstaff Jr Center for Education
CERIAS Tech Reort 2010-01 The eriod of the Bell numbers modulo a rime by Peter Montgomery, Sangil Nahm, Samuel Wagstaff Jr Center for Education and Research Information Assurance and Security Purdue University,
More informationFast Iterative Solution of Saddle Point Problems
Michele Benzi Department of Mathematics and Computer Science Emory University Atlanta, GA Acknowledgments NSF (Computational Mathematics) Maxim Olshanskii (Mech-Math, Moscow State U.) Zhen Wang (PhD student,
More informationOn the Chvatál-Complexity of Knapsack Problems
R u t c o r Research R e o r t On the Chvatál-Comlexity of Knasack Problems Gergely Kovács a Béla Vizvári b RRR 5-08, October 008 RUTCOR Rutgers Center for Oerations Research Rutgers University 640 Bartholomew
More informationAn Investigation on the Numerical Ill-conditioning of Hybrid State Estimators
An Investigation on the Numerical Ill-conditioning of Hybrid State Estimators S. K. Mallik, Student Member, IEEE, S. Chakrabarti, Senior Member, IEEE, S. N. Singh, Senior Member, IEEE Deartment of Electrical
More informationAdditive results for the generalized Drazin inverse in a Banach algebra
Additive results for the generalized Drazin inverse in a Banach algebra Dragana S. Cvetković-Ilić Dragan S. Djordjević and Yimin Wei* Abstract In this aer we investigate additive roerties of the generalized
More informationResearch Article Comparison of HPM and PEM for the Flow of a Non-newtonian Fluid between Heated Parallel Plates
Research Journal of Alied Sciences, Engineering and Technology 7(): 46-434, 4 DOI:.96/rjaset.7.793 ISSN: 4-7459; e-issn: 4-7467 4 Maxwell Scientific Publication Cor. Submitted: November, 3 Acceted: January
More informationA SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE
THE 19 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE K.W. Gan*, M.R. Wisnom, S.R. Hallett, G. Allegri Advanced Comosites
More informationSpecialized and hybrid Newton schemes for the matrix pth root
Secialized and hybrid Newton schemes for the matrix th root Braulio De Abreu Marlliny Monsalve Marcos Raydan June 15, 2006 Abstract We discuss different variants of Newton s method for comuting the th
More informationOn the Toppling of a Sand Pile
Discrete Mathematics and Theoretical Comuter Science Proceedings AA (DM-CCG), 2001, 275 286 On the Toling of a Sand Pile Jean-Christohe Novelli 1 and Dominique Rossin 2 1 CNRS, LIFL, Bâtiment M3, Université
More informationON POLYNOMIAL SELECTION FOR THE GENERAL NUMBER FIELD SIEVE
MATHEMATICS OF COMPUTATIO Volume 75, umber 256, October 26, Pages 237 247 S 25-5718(6)187-9 Article electronically ublished on June 28, 26 O POLYOMIAL SELECTIO FOR THE GEERAL UMBER FIELD SIEVE THORSTE
More informationSCHUR S LEMMA AND BEST CONSTANTS IN WEIGHTED NORM INEQUALITIES. Gord Sinnamon The University of Western Ontario. December 27, 2003
SCHUR S LEMMA AND BEST CONSTANTS IN WEIGHTED NORM INEQUALITIES Gord Sinnamon The University of Western Ontario December 27, 23 Abstract. Strong forms of Schur s Lemma and its converse are roved for mas
More informationRadial Basis Function Networks: Algorithms
Radial Basis Function Networks: Algorithms Introduction to Neural Networks : Lecture 13 John A. Bullinaria, 2004 1. The RBF Maing 2. The RBF Network Architecture 3. Comutational Power of RBF Networks 4.
More informationGENERALIZED NORMS INEQUALITIES FOR ABSOLUTE VALUE OPERATORS
International Journal of Analysis Alications ISSN 9-8639 Volume 5, Number (04), -9 htt://www.etamaths.com GENERALIZED NORMS INEQUALITIES FOR ABSOLUTE VALUE OPERATORS ILYAS ALI, HU YANG, ABDUL SHAKOOR Abstract.
More informationDETC2003/DAC AN EFFICIENT ALGORITHM FOR CONSTRUCTING OPTIMAL DESIGN OF COMPUTER EXPERIMENTS
Proceedings of DETC 03 ASME 003 Design Engineering Technical Conferences and Comuters and Information in Engineering Conference Chicago, Illinois USA, Setember -6, 003 DETC003/DAC-48760 AN EFFICIENT ALGORITHM
More informationSharp gradient estimate and spectral rigidity for p-laplacian
Shar gradient estimate and sectral rigidity for -Lalacian Chiung-Jue Anna Sung and Jiaing Wang To aear in ath. Research Letters. Abstract We derive a shar gradient estimate for ositive eigenfunctions of
More informationHomogeneous and Inhomogeneous Model for Flow and Heat Transfer in Porous Materials as High Temperature Solar Air Receivers
Excert from the roceedings of the COMSOL Conference 1 aris Homogeneous and Inhomogeneous Model for Flow and Heat ransfer in orous Materials as High emerature Solar Air Receivers Olena Smirnova 1 *, homas
More informationBrownian Motion and Random Prime Factorization
Brownian Motion and Random Prime Factorization Kendrick Tang June 4, 202 Contents Introduction 2 2 Brownian Motion 2 2. Develoing Brownian Motion.................... 2 2.. Measure Saces and Borel Sigma-Algebras.........
More informationSums of independent random variables
3 Sums of indeendent random variables This lecture collects a number of estimates for sums of indeendent random variables with values in a Banach sace E. We concentrate on sums of the form N γ nx n, where
More informationUniformly best wavenumber approximations by spatial central difference operators: An initial investigation
Uniformly best wavenumber aroximations by satial central difference oerators: An initial investigation Vitor Linders and Jan Nordström Abstract A characterisation theorem for best uniform wavenumber aroximations
More informationImproved Bounds on Bell Numbers and on Moments of Sums of Random Variables
Imroved Bounds on Bell Numbers and on Moments of Sums of Random Variables Daniel Berend Tamir Tassa Abstract We rovide bounds for moments of sums of sequences of indeendent random variables. Concentrating
More informationPOINTS ON CONICS MODULO p
POINTS ON CONICS MODULO TEAM 2: JONGMIN BAEK, ANAND DEOPURKAR, AND KATHERINE REDFIELD Abstract. We comute the number of integer oints on conics modulo, where is an odd rime. We extend our results to conics
More informationDetermining Momentum and Energy Corrections for g1c Using Kinematic Fitting
CLAS-NOTE 4-17 Determining Momentum and Energy Corrections for g1c Using Kinematic Fitting Mike Williams, Doug Alegate and Curtis A. Meyer Carnegie Mellon University June 7, 24 Abstract We have used the
More informationExistence of solutions to a superlinear p-laplacian equation
Electronic Journal of Differential Equations, Vol. 2001(2001), No. 66,. 1 6. ISSN: 1072-6691. URL: htt://ejde.math.swt.edu or htt://ejde.math.unt.edu ft ejde.math.swt.edu (login: ft) Existence of solutions
More informationEigenanalysis of Finite Element 3D Flow Models by Parallel Jacobi Davidson
Eigenanalysis of Finite Element 3D Flow Models by Parallel Jacobi Davidson Luca Bergamaschi 1, Angeles Martinez 1, Giorgio Pini 1, and Flavio Sartoretto 2 1 Diartimento di Metodi e Modelli Matematici er
More informationNUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS
NUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS Tariq D. Aslam and John B. Bdzil Los Alamos National Laboratory Los Alamos, NM 87545 hone: 1-55-667-1367, fax: 1-55-667-6372
More informationA Note on Guaranteed Sparse Recovery via l 1 -Minimization
A Note on Guaranteed Sarse Recovery via l -Minimization Simon Foucart, Université Pierre et Marie Curie Abstract It is roved that every s-sarse vector x C N can be recovered from the measurement vector
More informationTopic 7: Using identity types
Toic 7: Using identity tyes June 10, 2014 Now we would like to learn how to use identity tyes and how to do some actual mathematics with them. By now we have essentially introduced all inference rules
More informationA generalization of Amdahl's law and relative conditions of parallelism
A generalization of Amdahl's law and relative conditions of arallelism Author: Gianluca Argentini, New Technologies and Models, Riello Grou, Legnago (VR), Italy. E-mail: gianluca.argentini@riellogrou.com
More informationANALYTIC NUMBER THEORY AND DIRICHLET S THEOREM
ANALYTIC NUMBER THEORY AND DIRICHLET S THEOREM JOHN BINDER Abstract. In this aer, we rove Dirichlet s theorem that, given any air h, k with h, k) =, there are infinitely many rime numbers congruent to
More informationSpectral Clustering based on the graph p-laplacian
Sectral Clustering based on the grah -Lalacian Thomas Bühler tb@cs.uni-sb.de Matthias Hein hein@cs.uni-sb.de Saarland University Comuter Science Deartment Camus E 663 Saarbrücken Germany Abstract We resent
More informationOn the Fluid Dependence of Rock Compressibility: Biot-Gassmann Refined
Downloaded 0/9/3 to 99.86.4.8. Redistribution subject to SEG license or coyright; see Terms of Use at htt://library.seg.org/ On the luid Deendence of Rock Comressibility: Biot-Gassmann Refined Leon Thomsen,
More informationInformation collection on a graph
Information collection on a grah Ilya O. Ryzhov Warren Powell February 10, 2010 Abstract We derive a knowledge gradient olicy for an otimal learning roblem on a grah, in which we use sequential measurements
More informationBEST CONSTANT IN POINCARÉ INEQUALITIES WITH TRACES: A FREE DISCONTINUITY APPROACH
BEST CONSTANT IN POINCARÉ INEQUALITIES WITH TRACES: A FREE DISCONTINUITY APPROACH DORIN BUCUR, ALESSANDRO GIACOMINI, AND PAOLA TREBESCHI Abstract For Ω R N oen bounded and with a Lischitz boundary, and
More informationInformation collection on a graph
Information collection on a grah Ilya O. Ryzhov Warren Powell October 25, 2009 Abstract We derive a knowledge gradient olicy for an otimal learning roblem on a grah, in which we use sequential measurements
More informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
More informationChurilova Maria Saint-Petersburg State Polytechnical University Department of Applied Mathematics
Churilova Maria Saint-Petersburg State Polytechnical University Deartment of Alied Mathematics Technology of EHIS (staming) alied to roduction of automotive arts The roblem described in this reort originated
More informationMATH 2710: NOTES FOR ANALYSIS
MATH 270: NOTES FOR ANALYSIS The main ideas we will learn from analysis center around the idea of a limit. Limits occurs in several settings. We will start with finite limits of sequences, then cover infinite
More information