The nonlinear and non-stationary nature of Navier-Stokes equations produces fluid flows that can be noticeably different in appearance with subtle changes. In this paper we introduce a method that can analyze the intrinsic multiscale features of flow fields from a decomposition point of view, by using the Hilbert-Huang transform method on 3D fluid simulation. We show how this method can provide insights to flow styles and help modulate the fluid simulation with its internal physical information. We provide easy-toimplement algorithms that can be integrated with standard grid-based fluid simulation methods, and demonstrate how this approach can modulate the flow field and guide the simulation with different flow styles. The modulation is straightforward and relates directly to the flow¡¯s visual effect, with moderate computational overhead.

No abstract uploaded!

In 1960s, Almgren initiated a program to find minimal hypersurfaces in compact manifolds using min-max method. This program was largely advanced by Pitts and Schoen-Simon in 1980s when the manifold has no boundary. In this paper, we finish this program for general compact manifold with nonempty boundary. As a result, we prove the existence of a smooth embedded minimal hypersurface with free boundary in any compact smooth Euclidean domain. An application of our general existence result combined with the work of Marques and Neves shows that for any compact Riemannian manifolds with nonnegative Ricci curvature and convex boundary, there exist infinitely many embedded minimal hypersurfaces with free boundary which are properly embedded.

We quantized the full Einstein equations in a globally hyperbolic spacetime $N=N^{n+1}$, $n\ge 3$, and found solutions of the resulting hyperbolic equation in a fiber bundle $E$ which can be expressed as a product of spatial eigenfunctions (eigendistributions) and temporal eigenfunctions. The spatial eigenfunctions form a basis in an appropriate Hilbert space while the temporal eigenfunctions are solutions to a second order ordinary differential equation in $\R[]_+$. In case $n\ge 17$ and provided the cosmological constant $\Lam$ is negative the temporal eigenfunctions are eigenfunctions of a self-adjoint operator $\hat H_0$ such that the eigenvalues are countable and the eigenfunctions form an orthonormal basis of a Hilbert space.

In this paper, we study the superconvergence of the error for the discontinuous Galerkin (DG) finite element method for linear conservation laws when upwind fluxes are used. We prove that if we apply piecewise $k$-th degree polynomials, the error between the DG solution and the exact solution is ($k+2$)-th order superconvergent at the downwind-biased Radau points with suitable initial discretization. Moreover, we also prove the DG solution is ($k+2$)-th order superconvergent both for the cell averages and for the error to a particular projection of the exact solution. The superconvergence result in this paper leads to a new {\em a posteriori} error estimate. Our analysis is valid for arbitrary regular meshes and for $\mathcal{P}^k$ polynomials with arbitrary $k\geq1$, and for both periodic boundary conditions and for initial-boundary value problems. We perform numerical experiments to demonstrate that the superconvergence rate proved in this paper is optimal.