The minimal speeds (c) of the KolmogorovPetrovskyPiskunov (KPP) fronts at small diffusion ( 1) in a class of time-periodic cellular flows with chaotic streamlines is investigated in this paper. The variational principle of c reduces the computation to that of a principle eigenvalue problem on a periodic domain of a linear advectiondiffusion operator with spacetime periodic coefficients and small diffusion. To solve the advection dominated time-dependent eigenvalue problem efficiently over large time, a combination of spectral methods and finite element, as well as the associated fast solvers, are utilized to accelerate computation. In contrast to the scaling c= O ( 1/4) in steady cellular flows, a new relation c= O (1) as 1 is revealed in the time-periodic cellular flows due to the presence of chaotic streamlines. Residual propagation speed emerges from the Lagrangian chaos which is quantified as a sub

In this paper, we propose a rigorous and accurate non-local (in the oversampled region) upscaling framework based on some recently developed multiscale methods [10]. Our proposed method consists of identifying multi-continua parameters via local basis functions and constructing non-local (in the oversampled region) transfer and effective properties. To achieve this, we significantly modify our recent work proposed within Generalized Multiscale Finite Element Method (GMsFEM) in [10] and derive appropriate local problems in oversampled regions once we identify important modes representing each continuum. We use piecewise constant functions in each fracture network and in the matrix to write an upscaled equation. Thus, the resulting upscaled equation is of minimal size and the unknowns are average pressures in the fractures and the matrix. Note that the use of non-local upscaled model for porous

In the long-time scale, we consider the fluid dynamical limits for the kinetic equations when the fluctuation is decomposed into even and odd parts with respect to the microscopic velocity with different scalings. It is shown that when the background state is an absolute Maxwellian, the limit fluid dynamical equations are the incompressible Navier-Stokes equations with viscous heating. This is different from the case when the even and odd parts of the fluctuation have the same scaling where the standard incompressible Navier-Stokes equations without viscous heating are obtained. On the other hand, when the background is a local Maxwellian, it is shown that the above even-odd decomposition leads to a non-classical fluid dynamical system without viscous heating which has been used to describe the ghost effect in the kinetic theory. In addition, the above even-odd decomposition is justified rigorously for the Boltzmann equation for the former case when the background is an absolute Maxwellian.

We introduce a unified particle framework which integrates the phase-field method with multi-material simulation to allow modeling of both liquids and solids, as well as phase transitions between them. A simple elastoplastic model is used to capture the behavior of various kinds of solids, including deformable bodies, granular materials, and cohesive soils. States of matter or phases, particularly liquids and solids, are modeled using the non-conservative Allen-Cahn equation. In contrast, materials—made of different substances—are advected by the conservative Cahn-Hilliard equation. The distributions of phases and materials are represented by a phase variable and a concentration variable, respectively, allowing us to represent commonly observed fluid-solid interactions. Our multi-phase, multi-material system is governed by a unified Helmholtz free energy density. This framework provides the first method in computer graphics capable of modeling a continuous interface between phases. It is versatile and can be readily used in many scenarios that are challenging to simulate. Examples are provided to demonstrate the capabilities and effectiveness of this approach.

We develop an efficient numerical method to compute single slit or double slit diffraction patterns from high frequency wave in inhomogeneous media. We approximate the high frequency asymptotic solution to the Helmholtz equation using the Eulerian Gaussian beam summation proposed in [20, 21]. The emitted rays from a slit are embedded in the phase space using an open segment. The evolution of this open curve is accurately computed using the recently developed Grid Based Particle Method [24] which results in a very efficient computational algorithm. Following the grid based particle method we proposed in [23, 24], we represent the open curve or the open surface by meshless Lagrangian particles sampled according to an underlying fixed Eulerian mesh. The end-points of the open curve are tracked explicitly and consistently with interior particles. To construct the overall wavefield, each of these sampling particles also carry necessary quantities that are obtained by solving advection-reaction equations. Numerical experiments show that the resulting method can model diffraction patterns in inhomogeneous media accurately, even in the occurrence of caustics.