High Performance Computing Based Methods for Simulation and Optimisation of Flow Problems

The thesis is concerned with the study of methods in high-performance computing for simulation and optimisation of flow problems that occur in the framework of microflows. We consider the adequate use of techniques in parallel computing by means of finite element based solvers for partial differential equations and by means of sensitivity- and adjoint-based optimisation methods. The main focus is on three-dimensional, low Reynolds number flows described by the instationary Navier-Stokes equations

A newton-type method for fluid computation

Master'sMASTER OF ENGINEERIN

Final report on the Copper Mountain conference on multigrid methods

The Copper Mountain Conference on Multigrid Methods was held on April 6-11, 1997. It took the same format used in the previous Copper Mountain Conferences on Multigrid Method conferences. Over 87 mathematicians from all over the world attended the meeting. 56 half-hour talks on current research topics were presented. Talks with similar content were organized into sessions. Session topics included: fluids; domain decomposition; iterative methods; basics; adaptive methods; non-linear filtering; CFD; applications; transport; algebraic solvers; supercomputing; and student paper winners

Scalability of preconditioners as a strategy for parallel computation of compressible fluid flow

Parallel implementations of a Newton-Krylov-Schwarz algorithm are used to solve a model problem representing low Mach number compressible fluid flow over a backward-facing step. The Mach number is specifically selected to result in a numerically {open_quote}stiff{close_quotes} matrix problem, based on an implicit finite volume discretization of the compressible 2D Navier-Stokes/energy equations using primitive variables. Newton`s method is used to linearize the discrete system, and a preconditioned Krylov projection technique is used to solve the resulting linear system. Domain decomposition enables the development of a global preconditioner via the parallel construction of contributions derived from subdomains. Formation of the global preconditioner is based upon additive and multiplicative Schwarz algorithms, with and without subdomain overlap. The degree of parallelism of this technique is further enhanced with the use of a matrix-free approximation for the Jacobian used in the Krylov technique (in this case, GMRES(k)). Of paramount interest to this study is the implementation and optimization of these techniques on parallel shared-memory hardware, namely the Cray C90 and SGI Challenge architectures. These architectures were chosen as representative and commonly available to researchers interested in the solution of problems of this type. The Newton-Krylov-Schwarz solution technique is increasingly being investigated for computational fluid dynamics (CFD) applications due to the advantages of full coupling of all variables and equations, rapid non-linear convergence, and moderate memory requirements. A parallel version of this method that scales effectively on the above architectures would be extremely attractive to practitioners, resulting in efficient, cost-effective, parallel solutions exhibiting the benefits of the solution technique

Seventh Copper Mountain Conference on Multigrid Methods

The Seventh Copper Mountain Conference on Multigrid Methods was held on April 2-7, 1995 at Copper Mountain, Colorado. This book is a collection of many of the papers presented at the conference and so represents the conference proceedings. NASA Langley graciously provided printing of this document so that all of the papers could be presented in a single forum. Each paper was reviewed by a member of the conference organizing committee under the coordination of the editors. The vibrancy and diversity in this field are amply expressed in these important papers, and the collection clearly shows the continuing rapid growth of the use of multigrid acceleration techniques

The workshop on iterative methods for large scale nonlinear problems

The aim of the workshop was to bring together researchers working on large scale applications with numerical specialists of various kinds. Applications that were addressed included reactive flows (combustion and other chemically reacting flows, tokamak modeling), porous media flows, cardiac modeling, chemical vapor deposition, image restoration, macromolecular modeling, and population dynamics. Numerical areas included Newton iterative (truncated Newton) methods, Krylov subspace methods, domain decomposition and other preconditioning methods, large scale optimization and optimal control, and parallel implementations and software. This report offers a brief summary of workshop activities and information about the participants. Interested readers are encouraged to look into an online proceedings available at http://www.usi.utah.edu/logan.proceedings. In this, the material offered here is augmented with hypertext abstracts that include links to locations such as speakers` home pages, PostScript copies of talks and papers, cross-references to related talks, and other information about topics addresses at the workshop

Convection in the Melt

A physical problem involving the melting/freezing of a phase-change material (PCM) is the applied setting of this research. The development of models that couple the partial differential equations for energy transport and fluid motion with phases of differing densities is a primary goal of the research. In Chapter 2, a general framework is developed for the formulation of conservation laws that admit interfaces. A notion of weak solution is developed and its relation with classical and other weak formulations is discussed. Conditions that hold across various kinds of interfaces are also developed. The formulation is examined for the conservation of mass, momentum and energy in Chapter 3. In Chapter 4, a numerical method for the solution of conservation law equations is given. The method uses a Crank-Nicolson time discretization and solves the implicit equations with a Newton/Approximate Factorization technique. The method captures interfaces and is consistent with the control volume weak formulations of Chapter 2. The numerical solution converges to the distributional solution of the conservation law. In Chapter 5, three applications of the theory are developed and numerical computations are presented. First, a one dimensional problem is studied involving conservation of mass. momentum and energy in a phase-change material with a liquid density larger than that of the solid. The second application is a suction problem in two dimensions. The bulk movement of a liquid and void are simulated with and without the effects of surface tension. The third application is to a three-dimensional simulation of the heating of a cylindrical canister of PCM in 1-g and 0-g. For this simulation the Marangoni stress is the important driving force on the flow

The Sixth Copper Mountain Conference on Multigrid Methods, part 1

The Sixth Copper Mountain Conference on Multigrid Methods was held on 4-9 Apr. 1993, at Copper Mountain, CO. This book is a collection of many of the papers presented at the conference and as such represents the conference proceedings. NASA LaRC graciously provided printing of this document so that all of the papers could be presented in a single forum. Each paper was reviewed by a member of the conference organizing committee under the coordination of the editors. The multigrid discipline continues to expand and mature, as is evident from these proceedings. The vibrancy in this field is amply expressed in these important papers, and the collection clearly shows its rapid trend to further diversity and depth