Simple, efficient, and accurate finite difference methods are introduced for 3D unsteady viscous incompressible flows in the vorticity–vector potential formulation on nonstaggered grids. Two different types of methods are discussed. They differ in the implementation of the normal component of the vorticity boundary condition and consequently the enforcement of the divergence free condition for vorticity. Both second-order and fourth-order accurate schemes are developed. A detailed accuracy test is performed, revealing the structure of the error and the effect of how the convective terms are discretized near the boundary. The influence of the divergence free condition for vorticity to the overall accuracy is studied. Results on the cubic driven cavity flow at Reynolds number 500 and 3200 are shown and compared with that of the MAC scheme.

A new fourth-order accurate finite difference scheme for the computation of unsteady viscous incompressible flows is introduced. The scheme is based on the vorticity-stream function formulation. It is essentially compact and has the nice features of a compact scheme with regard to the treatment of boundary conditions. It is also very efficient, at every time step or Runge–Kutta stage, only two Poisson-like equations have to be solved. The Poisson-like equations are amenable to standard fast Poisson solvers usually designed for second order schemes. Detailed comparison with the second-order scheme shows the clear superiority of this new fourth-order scheme in resolving both the boundary layers and the gross features of the flow. This efficient fourth-order scheme also made it possible to compute the driven cavity flow at Reynolds number 106on a 10242grid at a reasonable cost. Fourth-order convergence is proved under mild regularity requirements. This is the first such result to our knowledge

We propose a multiscale model reduction method for partial differential equations. The main purpose of this method is to derive an effective equation for multiscale problems without scale separation. An essential ingredient of our method is to decompose the harmonic coordinates into a smooth part and a highly oscillatory part so that the smooth part is invertible and the highly oscillatory part is small. Such a decomposition plays a key role in our construction of the effective equation. We show that the solution to the effective equation is in H2, and can be approximated by a regular coarse mesh. When the multiscale problem has scale separation and a periodic structure, our method recovers the traditional homogenized equation. Furthermore, we provide error analysis for our method and show that the solution to the effective equation is close to the original multiscale solution in the H1 norm. Numerical results are presented to demonstrate the accuracy and robustness of the proposed method for several multiscale problems without scale separation, including a problem with a high contrast coefficient.

We present a novel area-preservation mapping/flattening method using the optimal mass transport technique, based on the Monge-Brenier theory. Our optimal transport map approach is rigorous and solid in theory, efficient and parallel in computation, yet general for various applications. By comparison with the conventional Monge-Kantorovich approach, our method reduces the number of variables from O(n²) to O(n), and converts the optimal mass transport problem to a convex optimization problem, which can now be efficiently carried out by Newton's method. Furthermore, our framework includes the area weighting strategy that enables users to completely control and adjust the size of areas everywhere in an accurate and quantitative way. Our method significantly reduces the complexity of the problem, and improves the efficiency, flexibility and scalability during visualization. Our framework, by combining

Photonic graphene, a photonic crystal with honeycomb structures, has been intensively studied in both theoretical and applied fields. Similar to graphene which admits Dirac Fermions and topological edge states, photonic graphene supports novel and subtle propagating modes (edge modes) of electromagnetic waves. These modes have wide applications in many optical systems. In this paper, we propose a novel gradient recovery method based on Bloch theory for the computation of topological edge modes in the honeycomb structure. Compared to standard finite element methods, this method provides higher order accuracy with the help of gradient recovery technique. This high order accuracy is highly desired for constructing the propagating electromagnetic modes in applications. We analyze the accuracy and prove the superconvergence of this method. Numerical examples are presented to show the efficiency by computing the edge mode for the P-symmetry and C-symmetry breaking cases in honeycomb structures.