Simulation of Multi-Phase Flow in Porous Media Through Integrated Upscaling, MPFA Discretization, and Adaptivity

Speaker: 

James Lambers

Institution: 

Stanford University

Time: 

Thursday, February 26, 2009 - 3:00pm

Location: 

RH 340P

In processes involving multi-phase flow in highly heterogeneous media, such as oil recovery by gas injection, mobile phases will seek high-permeability flow paths. Therefore, it is essential that models for such processes effectively account for these paths. For this purpose, we have developed a computational framework for flow solvers based on adapted Cartesian grids that are equipped with multi-point flux approximations obtained with specialized transmissibility upscaling methods.

For gridding, we propose using Cartesian Cell-based Anisotropically Refined (CCAR) grids, which inherit the ease of Cartesian grids while providing rapid transition between coarse and fine scales to resolve fine-scale features accurately and efficiently. We present an iterative algorithm for automatically generating such grids based on geological data and information from global coarse-scale flow simulations.

For upscaling, we discuss a local transmissibility upscaling method, called Variable Compact Multi-Point (VCMP), that uses spatially varying and compact multi-point flux stencils. The stencil weights are chosen so as to reproduce generic local flow problems accurately, while remaining as close as possible to a two-point flux for the sake of robustness. The inherent flexibility of VCMP can also be exploited to ensure that the solution of the resulting system satisfies a discrete maximum principle.

We conclude with application of these gridding, upscaling and discretization methods, originally designed for single-phase flow, to two-phase flow, which requires enhancing our adaptive mesh refinement scheme in order to accurately resolve rapidly advanacing saturation fronts. We show that adaptivity allows such accurate resolution by upscaling only single-phase parameters, thus avoiding the significant computational expense of multi-phase upscaling.

What is Different About the Ergodic Theory of Stochastic PDEs (vs ODEs)?

Speaker: 

Professor Jonathan Mattingly

Institution: 

Duke University

Time: 

Friday, November 14, 2008 - 4:00pm

Location: 

RH 306

I will discuss the difficulties which arise when one considers the long time behavior of a stochastically forced PDE. I will try to highlight that there are different cases which require very different ideas. Some cases can be seen as extensions of what is done in finite
dimensions, others require new tools and ideas. I will concentrate on the case of degenerately forced SPDEs. I will describe an extension of
Hormander's "sum of squares theorem" to hypo-elliptic operators in infinite dimensions. I will discuss the concert examples of the 2D
Navier Stokes equations on the torus and sphere as well as a class of reaction diffusion equations. In these contexts the discussion will center on the transfer of randomness between scales.

Nonlinear Stability of Periodic Traveling-Wave Solutions for the Benjamin-Ono Equation.

Speaker: 

Professor Jaime Angulo Pava

Institution: 

University of Sao Paulo,Brazil

Time: 

Thursday, November 13, 2008 - 3:00pm

Location: 

RH 340P

In this lecture, we present a method which has broad applicability to studies of nonlinear stability of periodic traveling-wave solutions for equations of KdV-type. In particular we obtain the existence and stability of a family of periodic traveling-wave solutions for the Benjamin-Ono equation via the classical Poisson summation theorem and positivity properties of the Fourier transform.

Method for the Linear Schroedinger Equation of N-interacting Particles

Speaker: 

Professor Claude Bardos

Institution: 

University of Paris 7

Time: 

Monday, October 27, 2008 - 4:00pm

Location: 

RH 306

This is a report on a joint work with Isabelle Catto, Norbert Mauser and Saber Trabelsi. The Multiconfiguration time dependent Hartree Fock Method (MCTDHF) is a nonlinear approximation of a linear system of /N/ quantum particles with binary interaction. It combines the principle of the Hartree Fock and the Galerkin approximation. The main difficulty is the introduction of a global (in space) density matrix $\Gamma(t) $ which may degenerate. By construction this approximation formally preserves the mass and the energy of the system. The conservation of energy can be used to balance the singularities Coulomb potential and to provide sufficient conditions for the global in time invertibility of $\Gamma(t)$.

In numerical computations this matrix is very often regularized (changed into $\Gamma(t) +\epsilon(t)$). In this situation the energy is no more conserved
and the mathematical analysis done in $L^2$ relies on Strichartz type estimates.

Nonlinear Diffusions and Image Processing

Speaker: 

Professor Patrick Guidotti

Institution: 

University of California, Irvine

Time: 

Thursday, October 30, 2008 - 3:00pm

Location: 

RH 340P

Since the seminal paper by Perona and Malik nonlinear diffusions have successfully been used for various image processing tasks. They also have attracted steadfast interest in the mathematical community. In this talk we will give an historical overview of the developments on the Perona-Malik equation and describe two new nonlinear diffusions which resolve the main mathematical shortcomings of Perona-Malik without sacrificing but rather enhancing the cherished practical qualities of the Perona-Malik model. The new equations, while
well-posed and purely diffusive, exhibit a non trivial dynamical behavior, which makes them mathematically interesting and practically
effective.

Coherent Vortex Extraction and Coherent Vortex Simulation of Turbulence

Speaker: 

Professor Kai Schneider

Institution: 

Universite de Provence (Aix-Marseille I), France

Time: 

Thursday, December 4, 2008 - 3:00pm

Location: 

RH 340P

Turbulence is characterized by its nonlinear and multiscale behaviour, self-organization into coherent structures and generic randomness. The number of active spatial and temporal scales involved increases with the Reynolds number, therefore it soon becomes prohibitive for direct numerical simulation. However, observations show that for a given flow realization these scales are not homogeneously distributed, neither in space nor in time, which corresponds to the flow intermittency. To be able to benefit from this property, a suitable representation of the flow should reflect the lacunarity of the fine scale activity, in both space and time.

A prominent tool for multiscale decompositions are wavelets. A wavelet is a well localized oscillating smooth function, i.e. a wave packet, which is dilated and translated. The thus obtained wavelet family allows to decompose a flow field into orthogonal scale-space contributions. The flow intermittency is reflected in the sparsity of the wavelet representation, i.e. only few coefficients, the strongest ones, are necessary to represent the dynamically active part of the flow. We will illustrate this by considering different 2D and 3D turbulent flows, either computed by direct numerical simulation (DNS) or measured by particle image velocimetry (PIV).

To compute the evolution of turbulent flows we have proposed the Coherent Vortex Simulation (CVS), which is based on the wavelet filtered Navier-Stokes equations. At each time step the turbulent fluctuations are split into two parts: the first corresponding to the coherent vortices which are kept, and the second to an incoherent background flow corresponding to turbulent dissipation which is discarded. We will present several simulations of 2D and 3D turbulent flows and show that CVS preserves their nonlinear dynamics.

Related publications can be downloaded from the following web pages:

http://wavelets.ens.fr
http://cmi.univ-mrs.fr/~kschneid

Regularity criteria for the 3D Navier-Stokes equations

Speaker: 

Chongshen Cao

Institution: 

Florida International University

Time: 

Thursday, October 2, 2008 - 3:00pm

Location: 

RH 340P

The question of global regularity for the 3D Navier-Stokes equations is a major open problem in applied analysis. It is well-known that the existence and uniqueness of strong solutions could be obtain under suitable additional assumptions. In this talk I will review some old and new results about the sufficient conditions for the global regularity to the 3D Navier-Stokes equations.

Gradient estimates and monotonicity formulas for linear Heat equations on manifolds with negative Ricci curvature

Speaker: 

RH 340P Xiangjin Xu

Institution: 

Binghamton university

Time: 

Thursday, October 9, 2008 - 3:00pm

Location: 

RH 340P

In the first part of the talk, using ideas from Ricci flow, we get a Li-Yau type gradient estimate for positive solutions of heat equation on Riemmannian manifolds with $Ricci(M)\ge -k$, $k\in \mathbb R$.
In the second part, we establish a Perelman type Li-Yau-Hamilton differential Harnack inequality for heat kernels on manifolds with $Ricci(M)\ge -k$. And we obtain various monotonicity formulas of entropy.

Pages

Subscribe to RSS - Nonlinear PDEs