Award Abstract # 1016332
Fast Interior Penalty Methods

NSF Org: DMS
Division Of Mathematical Sciences
Recipient: LOUISIANA STATE UNIVERSITY
Initial Amendment Date: August 26, 2010
Latest Amendment Date: August 26, 2010
Award Number: 1016332
Award Instrument: Standard Grant
Program Manager: Junping Wang
jwang@nsf.gov
 (703)292-4488
DMS
 Division Of Mathematical Sciences
MPS
 Directorate for Mathematical and Physical Sciences
Start Date: September 1, 2010
End Date: August 31, 2014 (Estimated)
Total Intended Award Amount: $300,951.00
Total Awarded Amount to Date: $300,951.00
Funds Obligated to Date: FY 2010 = $300,951.00
History of Investigator:
  • Susanne Brenner (Principal Investigator)
    brenner@math.lsu.edu
  • Li-yeng Sung (Co-Principal Investigator)
Recipient Sponsored Research Office: Louisiana State University
202 HIMES HALL
BATON ROUGE
LA  US  70803-0001
(225)578-2760
Sponsor Congressional District: 06
Primary Place of Performance: Louisiana State University
202 HIMES HALL
BATON ROUGE
LA  US  70803-0001
Primary Place of Performance
Congressional District:
06
Unique Entity Identifier (UEI): ECQEYCHRNKJ4
Parent UEI:
NSF Program(s): OFFICE OF MULTIDISCIPLINARY AC,
COMPUTATIONAL MATHEMATICS,
Information Technology Researc,
Special Projects - CCF
Primary Program Source: 01001011DB NSF RESEARCH & RELATED ACTIVIT
Program Reference Code(s): 6863, 9150, 9263
Program Element Code(s): 125300, 127100, 164000, 287800
Award Agency Code: 4900
Fund Agency Code: 4900
Assistance Listing Number(s): 47.049

ABSTRACT

This project will develop fast numerical methods for fourth and higher order partial differential equations using the interior penalty approach. The interior penalty approach has advantages over the classical approaches that use conforming, nonconforming or mixed finite elements in terms of the computational complexity, the existence of natuaral hierarchies of elements, the preservation of the symetric positive definiteness of the continuous problem, and the ease of deriving convergent schemes for complicated problems. Another significant advantage of interior penalty methods for higher order problems is due to the fact that discontinuous finite elements for higher order problems are also suitable for lower order problems. Therefore multigrid algorithms for interior penalty methods can be developed recursively through the hierarchy of elliptic problems. Namely, multigrid algorithms for second order problems can be embedded naturally in multigrid algorithms for fourth order problems, which can then be embedded naturally in multigrid algorithms for sixth order problems, and so on. The performance of these multigrid methods for higher order problems is comparable to the performance of multigrid methods for second order problems. This project will initiate a comprehensive study of interior penalty methods for higher order problems together with multigrid, domain decomposition and adaptive algorithms that will provide fast solvers for the resulting discrete problems. The results of this project will make it feasible to solve problems of order six and higher on general domains. Applications of these methods to strain gradient elasticity, plate buckling, the Monge-Ampere equations and the Cahn-Hilliard equations will also be investigated.

The fast algorithms developed in this project will make it practical for scientists and engineers to model complex phenomena by higher order partial differential equations. These algorithms will enhance the performance of numerical simulations in diverse areas such as structural mechanics, fluid mechanics, image processing, nanoscience, geometric optics, meteorology, optimal transport, differential geometry, and crystal growth, among many others.

PUBLICATIONS PRODUCED AS A RESULT OF THIS RESEARCH

Note:  When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

(Showing: 1 - 10 of 40)
A.T. Barker and S.C. Brenner "A mixed finite element method for the Stokes equations based on a weakly over-penalized symmetric interior penalty approach" Journal of Scientific Computing , 2013 DOI:10.1007/s10915-013-9733-9
A.T. Barker and S.C. Brenner "A mixed finite element method for the Stokes equations based on a weakly over-penalized symmetric interior penalty approach" Journal of Scientific Computing , v.58 , 2014 , p.290
A.T. Barker, S.C. Brenner and L.-Y. Sung "Overlapping Schwarz domain decomposition preconditioners for the local discontinuous Galerkin method for elliptic problems" Journal of Numerical Mathematics , v.19 , 2011 , p.165
A.T. Barker, S.C. Brenner and L.-Y. Sung "Overlapping Schwarz domain decomposition preconditionersfor the local discontinuous Galerkin method for ellipticproblems" Journal of Numerical Mathematics , v.19 , 2011 , p.165
A.T. Barker, S.C. Brenner, E.-H. Park and L.-Y. Sung "A nonoverlapping DD preconditioner for a weakly over-penalized symmetric interior penalty method" Domain Decomposition Methods in Science and Engineering XX Lecture Notes in Computational Science and Engineering , v.91 , 2013 , p.257
A.T. Barker, S.C. Brenner, E.-H. Park and L.-Y. Sung "Two-level additive Schwarz preconditioners for a weakly over-penalized symmetric interior penalty method" Journal of Scientific Computing , v.47 , 2011 , p.27 10.1007/s10915-010-9419-5
P.F. Antonietti, B. Ayuso de Dios, S.C. Brenner and L.-Y. Sung "Schwarz methods for a preconditioned WOPSIP method for elliptic problems" Computational Methods in Applied Mathematics , v.12 , 2012 , p.241
P.F. Antonietti, B. Ayuso de Dios, S.C. Brenner and L.-Y. Sung "Schwarz methods for a preconditioned WOPSIP method for elliptic problems" Computational Methods in Applied Mathematics , v.12 , 2012 , p.241
S.C. Brenner "An additive analysis of multiplicative Schwarz methods" Numerische Mathematik , 2012 10.1007/s00211-012-01479-3
S.C. Brenner "An additive analysis of multiplicative Schwarz methods" Numerische Mathematik , v.123 , 2013 , p.1
S.C. Brenner and K. Wang "An iterative substructuring algorithm for a $C^0$ interior penalty method" Electronic Transactions on Numerical Analysis , v.39 , 2012 , p.313
(Showing: 1 - 10 of 40)

PROJECT OUTCOMES REPORT

Disclaimer

This Project Outcomes Report for the General Public is displayed verbatim as submitted by the Principal Investigator (PI) for this award. Any opinions, findings, and conclusions or recommendations expressed in this Report are those of the PI and do not necessarily reflect the views of the National Science Foundation; NSF has not approved or endorsed its content.

New classes of numerical schemes for partial differential equations and related  optimization problems have resulted from this project. One class of schemes can be applied to higher order problems such as the bending of thin plates in mechanics  and the phase separation phenomenon in material science.  A second class of schemes can be applied to fully nonlinear partial differential equations arising from differential geometry and optimal transport.  Another class of schemes can be applied to two dimensional problems in electromagnetics.                                          

New fast solution methods involving multilevel techniques and parallel computing have also been developed for numerical schemes that can handle  complicated meshes and multi-physics couplings, and for numerical schemes that can handle saddle point  problems in structural and fluid mechanics.                                              

Several breakthroughs have been achieved in this project:

(1) A general framework for the numerical analysis of fourth order variational inequalities has been developed,  which fills a void in the literature that has existed for more than three decades.   

(2) A new approach to the numerical solution of two dimensional electromagnetic problems has been developed that would allow many existing methods for scalar equations  to be used in computational electromagnetics.

(3) A new class of multigrid methods  for saddle point problems appearing in mechanics has been developed. These methods      
 converge uniformly in the energy norm that underlies the physics and they also converge for general domains.

(4) A new general approach to the numerical solution of fully nonlinear partial differential equations has been developed.                                                                                                               
 
 The outcomes of this project will provide new tools for numerical simulations in
 science and engineering.


Last Modified: 11/30/2014
Modified by: Susanne C Brenner

Please report errors in award information by writing to: awardsearch@nsf.gov.

Print this page

Back to Top of page