skip to main content
Department of Computer Science University of Colorado Boulder
cu: home | engineering | mycuinfo | about | cu a-z | search cu | contact cu cs: about | calendar | directory | catalog | schedules | mobile | contact cs
home · events · colloquia · 2002-2003 · 

Colloquium - Paraschivoiu

Eaton South

Bounds for the Output of the Steady Incompressible Navier-Stokes Equations in Three Space Dimensions
Concordia University, Canada

A design engineer is often interested in optimizing a system based on some specific characteristic values of that system. For fluid mechanics or heat/mass transfer problems, these characteristic values, namely `outputs', may be the temperature at some location, or the flow rate of the underlying system, and can be expressed as linear functionals of the system's field variables, which are temperature, density, and velocities. For simulation based engineering design, the relevant system may be modeled by partial differential equations (PDEs) which yield field variable solutions. The accuracy of numerical solution of the PDEs depends of the discretization size and may required large computational resources.

Marius Paraschivoiu photo

In this work, an implicit a posteriori error estimation method is used to provide bounds to the linear functional output of the numerical solutions to the three dimensional Navier-Stokes equations at a fraction of the cost of calculating the output directly. The finite element bound method is based on a domain decomposition technique for efficiently generating lower and upper bounds. The novelty of this research is to utilize the posteriori error estimation method with Lagrange multipliers and the finite element tearing and interconnecting (FETI) procedure to extend the bound method to three space dimensions for Navier-Stokes (N-S) problems. The approach herein exploits a coarse mesh linearized N-S problem and demonstrates that the omitted nonlinear term is very small.

The computational advantage of the bound procedure is that a single coupled nonsymmetric/symmetric large problem can be decomposed into several uncoupled symmetric small problems. A simple benchmark problem, which is selected here to illustrate the procedure, is to find the lower and upper bounds of the average velocity of a pressure driven, incompressible, steady Newtonian fluid flow moving at low Reynolds numbers through an endless square channel which has an array of rectangular obstacles. Numerical results from the three dimensional Navier--Stokes equations show that the bounds are sharp and that the required computational resources decrease significantly. Parallel implementation on a Beowulf cluster is also reported. The bound method has the potential to be an inexpensive and reliable approach to incorporate numerical simulation of large-scale engineering problems in engineering design.

Hosted by Xiao-Chuan Cai.

The Department holds colloquia throughout the Fall and Spring semesters. These colloquia, open to the public, are typically held on Thursday afternoons, but sometimes occur at other times as well. If you would like to receive email notification of upcoming colloquia, subscribe to our Colloquia Mailing List. If you would like to schedule a colloquium, see Colloquium Scheduling.

Sign language interpreters are available upon request. Please contact Stephanie Morris at least five days prior to the colloquium.

See also:
Department of Computer Science
College of Engineering and Applied Science
University of Colorado Boulder
Boulder, CO 80309-0430 USA
Send email to

Engineering Center Office Tower
ECOT 717
FAX +1-303-492-2844
XHTML 1.0/CSS2 ©2012 Regents of the University of Colorado
Privacy · Legal · Trademarks
May 5, 2012 (13:29)