In earlier work, two of the authors developed a high order accurate numerical boundary condition procedure for hyperbolic conservation laws, which allows the computation using high order finite difference schemes on Cartesian meshes to solve problems in arbitrary physical domains whose boundaries do not coincide with grid lines. This procedure is based on the so-called inverse Lax-Wendroff (ILW) procedure for inflow boundary conditions and high order extrapolation for outflow boundary conditions. However, the algebra of the ILW procedure is quite heavy for two dimensional (2D) hyperbolic systems, which makes it difficult to implement the procedure for order of accuracy higher than three. In this paper, we first discuss a simplified and improved implementation for this procedure, which uses the relatively complicated ILW procedure only for the evaluation of the first order normal derivatives. Fifth order WENO type extrapolation is used for all other derivatives, regardless of the direction of the local characteristics and the smoothness of the solution. This makes the implementation of a fifth order boundary treatment practical for 2D systems with source terms. For no-penetration boundary condition of compressible inviscid flows, a further simplification is discussed, in which the evaluation of the tangential derivatives involved in the ILW procedure is avoided. We test our simplified and improved boundary treatment for Euler equations with or without source terms representing chemical reactions in detonations. The results demonstrate the designed fifth order accuracy, stability, and good performance for problems involving complicated interactions between detonation/shock waves and solid boundaries.

We develop a method searching for quasi-conformal maps from point-cloud surfaces to Euclidean plane. The maps can be got by directly solving Beltrami equation with proper boundary condition or solving optimization problem of the variational formu- lation. The point integral methd (PIM) is used for discretization on point clouds. Numerical experiments suggest that the method converges as the density of points in- creases.

This investigation studies nonemptiness problems of plane edge coloring with three colors. In the edge coloring (or Wang tiles) of a plane, unit squares with colored edges that have one of p colors are arranged side by side such that the touching edges of the adjacent tiles have the same colors. Given a basic set B of Wang tiles, the nonemptiness problem is to determine whether or not Σ(B) 6= ∅, where Σ(B) is the set of all global patterns on Z2 that can be constructed from the Wang tiles in B. Wang’s conjecture is that for any B of Wang tiles, Σ(B) 6= ∅ if and only if P(B) 6= ∅, where P(B) is the set of all periodic patterns on Z2 that can be generated by the tiles in B. When p ≥ 5, Wang’s conjecture is known to be wrong. When p = 2, theconjecture is true. This study proves that when p = 3, the conjecture is also true. If P(B) 6= ∅, then B has a subset B′ of minimal cycle generators such that P(B′) 6= ∅ and P(B′′) = ∅ for B′′ $ B′. This study demonstrates that the set C(3) of all minimal cycle generators contains 787, 605 members that can be classified into 2, 906 equivalence classes. N(3) is the set of all maximal non-cycle generators : if B ∈ N(3), then P(B) = ∅ and P(B ˜) 6= ∅ for B ˜ % B. Wang’s conjecture is shown to be true by proving that B ∈ N(3) implies Σ(B) = ∅.

No abstract uploaded!

No abstract uploaded!