Sample menu:

Articles written by Gottlieb, D.

  1. Chen, Q.-Y., Gottlieb, D., and Hesthaven, J.S.. "Uncertainty analysis for the steady-state flows in a dual throat nozzle" Journal of Computational Physics. 204 (1). 2005. pp. 378--398.

    It is well known that the steady state of an isentropic flow in a dual-throat nozzle with equal throat areas is not unique. In particular there is a possibility that the flow contains a shock wave, whose location is determined solely by the initial condition. In this paper, we consider cases with uncertainty in this initial condition and use generalized polynomial chaos methods to study the steady-state solutions for stochastic initial conditions. Special interest is given to the statistics of the shock location. The polynomial chaos (PC) expansion modes are shown to be smooth functions of the spatial variable x, although each solution realization is discontinuous in the spatial variable x. When the variance of the initial condition is small, the probability density function of the shock location is computed with high accuracy. Otherwise, many terms are needed in the PC expansion to produce reasonable results due to the slow convergence of the PC expansion, caused by non-smoothness in random space.


  2. Gottlieb, D. and Xiu, D.. "Galerkin method for wave equations with uncertain coefficients" Communications in Computational Physics. 3 (2). FEB 2008. pp. 505--518.

    Polynomial chaos methods (and generalized polynomial chaos methods) have been extensively applied to analyze PDEs that contain uncertainties. However this approach is rarely applied to hyperbolic systems. In this paper we analyze the properties of the resulting deterministic system of equations obtained by stochastic Galerkin projection. We consider a simple model of a scalar wave equation with random wave speed. We show that when uncertainty causes the change of characteristic directions, the resulting deterministic system of equations is a symmetric hyperbolic system with both positive and negative eigenvalues. A consistent method of imposing the boundary conditions is proposed and its convergence is established. Numerical examples are presented to support the analysis.