This article studies the structure and its evolution of incompressible flows with the anti-symmetry using a combination of rigorous analysis and numerical simulations, with an application to an example of oceanic flow. In particular, necessary and sufficient conditions for 2D divergence-free vector fields with anti-symmetry are obtained, and a detailed numerical simulation for a simplified model of Marsigli oceanic flow is provided to explore and verify the structure and its transitions of the flow. It is expected that the study will lead to useful insights to the understanding of the flow dynamics from both the mathematical and physical points of view.

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.

Local discontinuous Galerkin (LDG) methods are popular for convection-diffusion equations. In LDG methods, we introduce an auxiliary variable $p$ to represent the derivative of the primary variable $u$, and solve them on the same mesh. In this paper, we will introduce a new LDG method, and solve $u$ and $p$ on different meshes. The stability and error estimates will be investigated. The new algorithm is more flexible and flux-free for pure diffusion equations without introducing additional computational cost compared with the original LDG methods, since it is not necessary to solve each equation twice. Moreover, it is possible to construct third-order maximum-principle-preserving schemes based on the new algorithm. However, one cannot anticipate optimal accuracy in some special cases. In this paper, we will find out the reason for accuracy degeneration which further leads to several alternatives to obtain optimal convergence rates. Finally, several numerical experiments will be given to demonstrate the stability and optimal accuracy of the new algorithm.

No abstract uploaded!

No abstract uploaded!