In this paper, we consider local multiscale model reduction for problems with multiple scales in space and time. We developed our approaches within the framework of the Generalized Multiscale Finite Element Method (GMsFEM) using spacetime coarse cells. The main idea of GMsFEM is to construct a local snapshot space and a local spectral decomposition in the snapshot space. Previous research in developing multiscale spaces within GMsFEM focused on constructing multiscale spaces and relevant ingredients in space only. In this paper, our main objective is to develop a multiscale model reduction framework within GMsFEM that uses spacetime coarse cells. We construct spacetime snapshot and offline spaces. We compute these snapshot solutions by solving local problems. A complete snapshot space will use all possible boundary conditions; however, this can be very expensive. We propose using

Simulating complex processes in fractured media requires some type of model reduction. Well-known approaches include multi-continuum techniques, which have been commonly used in approximating subgrid effects for flow and transport in fractured media. Our goal in this paper is to (1) show a relation between multi-continuum approaches and Generalized Multiscale Finite Element Method (GMsFEM) and (2) to discuss coupling these approaches for solving problems in complex multiscale fractured media. The GMsFEM, a systematic approach, constructs multiscale basis functions via local spectral decomposition in pre-computed snapshot spaces. We show that GMsFEM can automatically identify separate fracture networks via local spectral problems. We discuss the relation between these basis functions and continuums in multi-continuum methods. The GMsFEM can automatically detect each

In this paper, we apply a simple finite element numerical scheme, proposed in an earlier work (Liu in Math Comput 70(234):579-593, 2000), to perform a high resolution numerical simulation of incompressible flow over an irregular domain and analyze its boundary layer separation. Compared with many classical finite element fluid solvers, this numerical method avoids a Stokes solver, and only two Poisson-like equations need to be solved at each time step/stage. In addition, its combination with the fully explicit fourth order Runge-Kutta (RK4) time discretization enables us to compute high Reynolds number flow in a very efficient way. As an application of this robust numerical solver, the dynamical mechanism of the boundary layer separation for a triangular cavity flow with Reynolds numbers \(Re=10^4\) and \(Re=10^5\) , including the precise values of bifurcation location and critical time, are reported in this paper. In addition, we provide a super-convergence analysis for the simple finite element numerical scheme, using linear elements over a uniform triangulation with right triangles.

In this paper, we apply high-order finite difference (FD) schemes for multispecies and multireaction detonations (MMD). In MMD, the density and pressure are positive and the mass fraction of the ith species in the chemical reaction, say zi, is between 0 and 1, with ∑z_i=1. Due to the lack of maximum-principle, most of the previous bound-preserving technique cannot be applied directly. To preserve those bounds, we will use the positivity-preserving technique to all the z′is and enforce ∑z_i=1 by constructing conservative schemes, thanks to conservative time integrations and consistent numerical fluxes in the system. Moreover, detonation is an extreme singular mode of flame propagation in premixed gas, and the model contains a significant stiff source. It is well known that for hyperbolic equations with stiff source, the transition points in the numerical approximations near the shocks may trigger spurious shock speed, leading to wrong shock position. Intuitively, the high-order weighted essentially non-oscillatory (WENO) scheme, which can suppress oscillations near the discontinuities, would be a good choice for spatial discretization. However, with the nonlinear weights, the numerical fluxes are no longer “consistent”, leading to nonconservative numerical schemes and the bound-preserving technique does not work. Numerical experiments demonstrate that, without further numerical techniques such as subcell resolutions, the conservative FD method with linear weights can yield better numerical approximations than the nonconservative WENO scheme.

No abstract uploaded!