Results 1  10
of
74
Functional Phonology  Formalizing the interactions between articulatory and perceptual drives
, 1998
"... ..."
ImplicitExplicit Methods For TimeDependent PDEs
 SIAM J. Numer. Anal
, 1997
"... . Implicitexplicit (IMEX) schemes have been widely used, especially in conjunction with spectral methods, for the time integration of spatially discretized PDEs of diffusionconvection type. Typically, an implicit scheme is used for the diffusion term and an explicit scheme is used for the convecti ..."
Abstract

Cited by 105 (6 self)
 Add to MetaCart
. Implicitexplicit (IMEX) schemes have been widely used, especially in conjunction with spectral methods, for the time integration of spatially discretized PDEs of diffusionconvection type. Typically, an implicit scheme is used for the diffusion term and an explicit scheme is used for the convection term. Reactiondiffusion problems can also be approximated in this manner. In this work we systematically analyze the performance of such schemes, propose improved new schemes and pay particular attention to their relative performance in the context of fast multigrid algorithms and of aliasing reduction for spectral methods. For the prototype linear advectiondiffusion equation, a stability analysis for first, second, third and fourth order multistep IMEX schemes is performed. Stable schemes permitting large time steps for a wide variety of problems and yielding appropriate decay of high frequency error modes are identified. Numerical experiments demonstrate that weak decay of high freque...
Parallel NewtonKrylovSchwarz Algorithms For The Transonic Full Potential Equation
, 1998
"... We study parallel twolevel overlapping Schwarz algorithms for solving nonlinear finite element problems, in particular, for the full potential equation of aerodynamics discretized in two dimensions with bilinear elements. The overall algorithm, NewtonKrylovSchwarz (NKS), employs an inexact finite ..."
Abstract

Cited by 42 (27 self)
 Add to MetaCart
We study parallel twolevel overlapping Schwarz algorithms for solving nonlinear finite element problems, in particular, for the full potential equation of aerodynamics discretized in two dimensions with bilinear elements. The overall algorithm, NewtonKrylovSchwarz (NKS), employs an inexact finitedifference Newton method and a Krylov space iterative method, with a twolevel overlapping Schwarz method as a preconditioner. We demonstrate that NKS, combined with a density upwinding continuation strategy for problems with weak shocks, is robust and economical for this class of mixed elliptichyperbolic nonlinear partial differential equations, with proper specification of several parameters. We study upwinding parameters, inner convergence tolerance, coarse grid density, subdomain overlap, and the level of fillin in the incomplete factorization, and report their effect on numerical convergence rate, overall execution time, and parallel efficiency on a distributedmemory parallel computer.
Nonlinearly preconditioned inexact Newton algorithms
 SIAM J. Sci. Comput
, 2000
"... Abstract. Inexact Newton algorithms are commonlyused for solving large sparse nonlinear system of equations F (u ∗ ) = 0 arising, for example, from the discretization of partial differential equations. Even with global strategies such as linesearch or trust region, the methods often stagnate at loc ..."
Abstract

Cited by 35 (14 self)
 Add to MetaCart
Abstract. Inexact Newton algorithms are commonlyused for solving large sparse nonlinear system of equations F (u ∗ ) = 0 arising, for example, from the discretization of partial differential equations. Even with global strategies such as linesearch or trust region, the methods often stagnate at local minima of �F �, especiallyfor problems with unbalanced nonlinearities, because the methods do not have builtin machineryto deal with the unbalanced nonlinearities. To find the same solution u ∗ , one maywant to solve instead an equivalent nonlinearlypreconditioned system F(u ∗ ) = 0 whose nonlinearities are more balanced. In this paper, we propose and studya nonlinear additive Schwarzbased parallel nonlinear preconditioner and show numericallythat the new method converges well even for some difficult problems, such as high Reynolds number flows, where a traditional inexact Newton method fails. Key words. nonlinear preconditioning, inexact Newton methods, Krylov subspace methods, nonlinear additive Schwarz, domain decomposition, nonlinear equations, parallel computing, incompressible
Approximate Solutions of Nonlinear Conservation Laws and Related Equations
, 1997
"... During the recent decades there was an enormous amount of activity related to the construction and analysis of modern algorithms for the approximate solution of nonlinear hyperbolic conservation laws and related problems. To present some aspects of this successful activity, we discuss the analytical ..."
Abstract

Cited by 34 (11 self)
 Add to MetaCart
During the recent decades there was an enormous amount of activity related to the construction and analysis of modern algorithms for the approximate solution of nonlinear hyperbolic conservation laws and related problems. To present some aspects of this successful activity, we discuss the analytical tools which are used in the development of convergence theories for these algorithms. These include classical compactness arguments (based on BV a priori estimates), the use of compensated compactness arguments (based on H^1compact entropy production), measure valued solutions (measured by their negative entropy production), and finally, we highlight the most recent addition to this bag of analytical tools  the use of averaging lemmas which yield new compactness and regularity results for nonlinear conservation laws and related equations. We demonstrate how these analytical tools are used in the convergence analysis of approximate solutions for hyperbolic conservation laws and related equations. Our discussion includes examples of Total Variation Diminishing (TVD) finitedifference schemes; error estimates derived from the onesided stability of Godunovtype methods for convex conservation laws (and their multidimensional analogue  viscosity solutions of demiconcave HamiltonJacobi equations); we outline, in the onedimensional case, the convergence proof of finiteelement streamlinediffusion and spectral viscosity schemes based on the divcurl lemma; we also address the questions of convergence and error estimates for multidimensional finitevolume schemes on nonrectangular grids; and finally, we indicate the convergence of approximate solutions with underlying kinetic formulation, e.g., finitevolume and relaxation schemes, once their regularizing effect is quantified in terms of the averaging lemma.
Finite element methods of leastsquares type
, 1998
"... We consider the application of leastsquares variational principles to the numerical solution of partial differential equations. Our main focus is on the development of leastsquares finite element methods for elliptic boundary value problems arising in fields such as fluid flows, linear elasticit ..."
Abstract

Cited by 31 (4 self)
 Add to MetaCart
We consider the application of leastsquares variational principles to the numerical solution of partial differential equations. Our main focus is on the development of leastsquares finite element methods for elliptic boundary value problems arising in fields such as fluid flows, linear elasticity, and convectiondiffusion. For many of these problems, leastsquares principles offer numerous theoretical and computational advantages in the algorithmic design and implementation of corresponding finite element methods, that are not present in standard Galerkin discretizations. Most notably, the use of leastsquares principles leads to symmetric and positive definite algebraic problems and allows one to circumvent stability conditions such as the infsup condition arising in mixed methods for the Stokes and NavierStokes equations. As a result, application of leastsquares principles has led to the development of robust and efficient finite element methods for a large class of problems of practical importance.
Comparison of some Flux Corrected Transport and Total Variation Diminishing Numerical Schemes for Hydrodynamic and Magnetohydrodynamic Problems
, 1996
"... Two versions of flux corrected transport and two versions of total variation diminishing schemes are tested for several one and twodimensional hydrodynamic and magnetohydrodynamic problems. Two of the schemes, YDFCT and TVDLF are tested extensively for the first time. The results give an insight i ..."
Abstract

Cited by 25 (7 self)
 Add to MetaCart
Two versions of flux corrected transport and two versions of total variation diminishing schemes are tested for several one and twodimensional hydrodynamic and magnetohydrodynamic problems. Two of the schemes, YDFCT and TVDLF are tested extensively for the first time. The results give an insight into the limitations of the methods, their relative strengths and weaknesses. Some subtle points of the algorithms and the effects of selecting different options for certain methods are emphasised. 1 INTRODUCTION Many interesting and important problems arise in astrophysical, solar, magnetospheric and thermonuclear research which can be described by the system of magnetohydrodynamic (MHD) equations. The complexity of these problems often prohibits an analytical investigation and/or only some of the variables can be observed or measured experimentally, thus the researcher has to rely on numerical simulations. In many situations, MHD flows develop steep gradients, shock waves, contact disconti...
Krylov subspace acceleration of nonlinear multigrid schemes
 Electronic Transactions in Numerical Analysis, 6:271290
, 1997
"... Abstract. In this paper we present a Krylov acceleration technique for nonlinear PDEs. As a ‘preconditioner’ we use nonlinear multigrid schemes such as the Full Approximation Scheme (FAS) [1]. The benefits of nonlinear multigrid used in combination with the new accelerator are illustrated by difficu ..."
Abstract

Cited by 22 (2 self)
 Add to MetaCart
Abstract. In this paper we present a Krylov acceleration technique for nonlinear PDEs. As a ‘preconditioner’ we use nonlinear multigrid schemes such as the Full Approximation Scheme (FAS) [1]. The benefits of nonlinear multigrid used in combination with the new accelerator are illustrated by difficult nonlinear elliptic scalar problems, such as the Bratu problem, and for systems of nonlinear equations, such as the NavierStokes equations. Key words. nonlinear Krylov acceleration, nonlinear multigrid, robustness, restarting conditions. AMS subject classifications. 65N55, 65H10, 65Bxx. 1. Introduction. It is well known that multigrid solution methods are optimal O(N) solvers, when all components in a method are chosen correctly. For difficult problems, such as some systems of nonlinear equations, it is far from trivial to choose these optimal components. The influence on the multigrid convergence of combinations of complicated factors, like convectiondominance, anisotropies, nonlinearities or non Mmatrix properties (the
Application and Accuracy of the Parallel Diagonal Dominant Algorithm
 Parallel Comput
, 1995
"... The Parallel Diagonal Dominant (PDD) algorithm is an efficient tridiagonal solver. In this paper, a detailed study of the PDD algorithm is given. First the PDD algorithm is extended to solve periodic tridiagonal systems and its scalability is studied. Then the reduced PDD algorithm, which has a smal ..."
Abstract

Cited by 10 (9 self)
 Add to MetaCart
The Parallel Diagonal Dominant (PDD) algorithm is an efficient tridiagonal solver. In this paper, a detailed study of the PDD algorithm is given. First the PDD algorithm is extended to solve periodic tridiagonal systems and its scalability is studied. Then the reduced PDD algorithm, which has a smaller operation count than that of the conventional sequential algorithm for many applications, is proposed. Accuracy analysis is provided for a class of tridiagonal systems, the symmetric and skewsymmetric Toeplitz tridiagonal systems. Implementation results show that the analysis gives a good bound on the relative error, and the PDD and reduced PDD algorithms are good candidates for emerging massively parallel machines. Index Terms: Parallel processing, Parallel numerical algorithms, Scalable computing, Tridiagonal systems, Toeplitz systems Manuscript received April 7, 1993; revised April 7, 1994 and January 27, 1995. This research was supported in part by the National Aeronautics and S...
Numerical solution of the timedomain Maxwell equations using highaccuracy finitedifference methods
 SIAM J. Sci. Comput
"... Abstract. Highaccuracy finitedifference schemes are used to solve the twodimensional timedomain Maxwell equations for electromagnetic wave propagation and scattering. The highaccuracy schemes consist of a sevenpoint spatial operator coupled with a sixstage Runge–Kutta timemarching method. Two ..."
Abstract

Cited by 8 (5 self)
 Add to MetaCart
Abstract. Highaccuracy finitedifference schemes are used to solve the twodimensional timedomain Maxwell equations for electromagnetic wave propagation and scattering. The highaccuracy schemes consist of a sevenpoint spatial operator coupled with a sixstage Runge–Kutta timemarching method. Two methods are studied, one of which produces the maximum order of accuracy and one of which is optimized for propagation distances smaller than roughly 300 wavelengths. Boundary conditions are presented which preserve the accuracy of these schemes when modeling interfaces between different materials. Numerical experiments are performed which demonstrate the utility of the highaccuracy schemes in modeling waves incident on dielectric and perfectconducting scatterers using Cartesian and curvilinear grids. The highaccuracy schemes are shown to be substantially more efficient, in both computing time and memory, than a secondorder and a fourthorder method. The optimized scheme can lead to a reduction in error relative to the maximumorder scheme, with no additional expense, especially when the number of wavelengths of travel is large.