Potential-based Formulations of the Navier-Stokes Equations and their Application

Based on a Clebsch-like velocity representation and a combination of classical variational principles for the special cases of ideal and Stokes flow a novel discontinuous Lagrangian is constructed; it bypasses the known problems associated with non-physical solutions and recovers the classical Navier-Stokes equations together with the balance of inner energy in the limit when an emerging characteristic frequency parameter tends to infinity. Additionally, a generalized Clebsch transformation for viscous flow is established for the first time. Next, an exact first integral of the unsteady, three-dimensional, incompressible Navier-Stokes equations is derived; following which gauge freedoms are explored leading to favourable reductions in the complexity of the equation set and number of unknowns, enabling a self-adjoint variational principle for steady viscous flow to be constructed. Concurrently, appropriate commonly occurring physical and auxiliary boundary conditions are prescribed, including establishment of a first integral for the dynamic boundary condition at a free surface. Starting from this new formulation, three classical flow problems are considered, the results obtained being in total agreement with solutions in the open literature. A new least-squares finite element method based on the first integral of the steady two-dimensional, incompressible, Navier-Stokes equations is developed, with optimal convergence rates established theoretically. The method is analysed comprehensively, thoroughly validated and shown to be competitive when compared to a corresponding, standard, primitive-variable, finite element formulation. Implementation details are provided, and the well-known problem of mass conservation addressed and resolved via selective weighting. The attractive positive definiteness of the resulting linear systems enables employment of a customized scalable algebraic multigrid method for efficient error reduction. The solution of several engineering related problems from the fields of lubrication and film flow demonstrate the flexibility and efficiency of the proposed method, including the case of unsteady flow, while revealing new physical insights of interest in their own right

Reactive Flow and Transport Through Complex Systems

The meeting focused on mathematical aspects of reactive flow, diffusion and transport through complex systems. The research interest of the participants varied from physical modeling using PDEs, mathematical modeling using upscaling and homogenization, numerical analysis of PDEs describing reactive transport, PDEs from fluid mechanics, computational methods for random media and computational multiscale methods

Institute for Computational Mechanics in Propulsion (ICOMP)

The Institute for Computational Mechanics in Propulsion (ICOMP) is a combined activity of Case Western Reserve University, Ohio Aerospace Institute (OAI) and NASA Lewis. The purpose of ICOMP is to develop techniques to improve problem solving capabilities in all aspects of computational mechanics related to propulsion. The activities at ICOMP during 1991 are described

Statistical Properties of Turbulence: An Overview

We present an introductory overview of several challenging problems in the statistical characterisation of turbulence. We provide examples from fluid turbulence in three and two dimensions, from the turbulent advection of passive scalars, turbulence in the one-dimensional Burgers equation, and fluid turbulence in the presence of polymer additives.Comment: 34 pages, 31 figure

ICASE semiannual report, April 1 - September 30, 1989

The Institute conducts unclassified basic research in applied mathematics, numerical analysis, and computer science in order to extend and improve problem-solving capabilities in science and engineering, particularly in aeronautics and space. The major categories of the current Institute for Computer Applications in Science and Engineering (ICASE) research program are: (1) numerical methods, with particular emphasis on the development and analysis of basic numerical algorithms; (2) control and parameter identification problems, with emphasis on effective numerical methods; (3) computational problems in engineering and the physical sciences, particularly fluid dynamics, acoustics, and structural analysis; and (4) computer systems and software, especially vector and parallel computers. ICASE reports are considered to be primarily preprints of manuscripts that have been submitted to appropriate research journals or that are to appear in conference proceedings

Analysis and optimization of film cooling effectiveness

In the first part, an optimization strategy is described that combines high-fidelity simu- lations with response surface construction, and is applied to pulsed film cooling for turbine blades. The response surface is constructed for the film cooling effectiveness as a function of duty cycle, in the range of DC between 0.05 and 1, and pulsation frequency St in the range of 0.2-2, using a pseudo-spectral projection method. The jet is fully modulated and the blowing ratio, when the jet is on, is 1.5 in all cases. Overall 73 direct numerical sim- ulations (DNS) using spectral element method were performed to sample the film cooling effectiveness on a Clenshaw-Curtis grid in the design space. It is observed that in the parameter space explored a global optimum exists, and in the present study, the best film cooling effectiveness is found at DC = 0.14 and St = 1.03. In the same range of DC and St, four other local optimums were found. The gradient-based optimization algorithms are argued to be unsuitable for the current problem due to the non-convexity of the objective function. In the second part, the effect of randomness of blowing ratio on film cooling performance is investigated by combining direct numerical simulations with a stochastic collocation ap- proach. The blowing ratio variations are assumed to have a truncated Gaussian distribution with mean of 0.3 and the standard variation of approximately 0.1. The parametric space is discretized using Multi-Element general Polynomial Chaos (ME-gPC) with five elements where general polynomial chaos of order 3 is used in each element. Direct numerical simula- tions were carried out using spectral/hp element method to sample the governing equations in space and time. The probability density function of the film cooling effectiveness was obtained and the standard deviation of the adiabatic film cooling effectiveness on the blade surface was calculated. A maximum standard deviation of 15% was observed in the region within a four-jet-diameter distance downstream of the exit hole. The spatially-averaged adiabatic film cooling effectiveness was 0.23 0.02. The calculation of all the statistical properties were carried out as off-line post-processing. Overall the computational strategy is shown to be very effective with the total computational cost being equivalent to solving twenty independent direct numerical simulations that are performed concurrently. In the third part, an accurate and efficient finite difference method for solving the incompressible Navier-Stokes equations on curvilinear grids is developed. This method combines the favorable features of the staggered grid and semi-staggered grid approaches. A novel symmetric finite difference discretization of the Poisson-Neumann problem on curvilinear grids is also presented. The validity of the method is demonstrated on four benchmark problems. The Taylor-Green vortex problem is solved on a curvilinear grid with highly skewed cells and a second-order convergence in .-norm is observed. The mixed convection in a lid-driven cavity is solved on a highly curvilinear grid and excellent agreement with literature is obtained. The results for flow past a cylinder are compared with the existing experimental data in the literature. As the fourth case, three dimensional time-dependent incompressible flow in a curved tube is solved. The predictions agree well with the measured data, and validate the approach used

Second All-Union Seminar on Hydromechanics and Heat and Mass Exchange in Weightlessness, summaries of reports

Abstracts of reports are given which were presented at the Second All Union Seminar on Hydromechanics and Heat-Mass Transfer in Weightlessness. Topics include: (1) features of crystallization of semiconductor materials under conditions of microacceleration; (2) experimental results of crystallization of solid solutions of CDTE-HGTE under conditions of weightlessness; (3) impurities in crystals cultivated under conditions of weightlessness; and (4) a numerical investigation of the distribution of impurities during guided crystallization of a melt