No abstract uploaded!

The influences of non-Oberbeck-Boussinesq (NOB) effects on flow instabilities and bifurcation characteristics of Rayleigh-Benard convection are examined. The working fluid is air with reference Prandtl number Pr = 0.71 and contained in two-dimensional rigid cavities of finite aspect ratios. The fluid flow is governed by the low-Mach-number equations, accounting for the NOB effects due to large temperature difference involving flow compressibility and variations of fluid viscosity and thermal conductivity with temperature. The intensity of NOB effects is measured by the dimensionless temperature differential epsilon. Linear stability analysis of the thermal conduction state is performed. An epsilon(2) scaling of the leading-order corrections of critical Rayleigh number Ra-cr and disturbance growth rate sigma due to NOB effects is identified, which is a consequence of an intrinsic symmetry of the system. The influences of weak NOB effects on flow instabilities are further studied by perturbation expansion of linear stability equations with regard to epsilon, and then the influence of aspect ratio A is investigated in detail. NOB effects are found to enhance (weaken) flow stability in large (narrow) cavities. Detailed contributions of compressibility, viscosity and buoyancy actions on disturbance kinetic energy growth are identified quantitatively by energy analysis. Besides, a weakly nonlinear theory is developed based on centre-manifold reduction to investigate the NOB influences on bifurcation characteristics near convection onset, and amplitude equations are constructed for both codimension-one and -two cases. Rich bifurcation regimes are observed based on amplitude equations and also confirmed by direct numerical simulation. Weakly nonlinear analysis is useful for organizing and understanding these simulation results.

In this paper, we develop, analyze and test the Fourier spectral methods for solving the Degasperis-Procesi(DP) equation which contains nonlinear high order derivatives, and possibly discontinuous or sharp transition solutions. The $L^2$ stability is obtained for general numerical solutions of the Fourier Galerkin method and Fourier collocation (pseudospectral) method.By applying the Gegenbauer reconstruction technique as a post-processing method to the Fourier spectral solution, we reduce the oscillations arising from the discontinuity successfully. The numerical simulation results for different types of solutions of the nonlinear Degasperis-Procesi equation are provided to illustrate the accuracy and capability of the methods.

In this work we investigate the decorated domain wall construction in bosonic group-cohomology symmetry-protected topological (SPT) phases and related quantum anomalies in bosonic topological phases. We first show that a general decorated domain wall construction can be described mathematically as an Atiyah-Hirzebruch spectral sequence, where the terms on the E2 page correspond to decorations by lower-dimensional SPT states at domain wall junctions. For bosonic group-cohomology SPT phases, the spectral sequence becomes the Lyndon-Hochschild-Serre (LHS) spectral sequence for ordinary group cohomology. We then discuss the physical interpretations of the differentials in the spectral sequence, particularly in the context of anomalous SPT phases and symmetry-enriched gauge theories. As the main technical result, we obtain a full description of the LHS spectral sequence concretely at the cochain level. The explicit formulae are then applied to explain Lieb-Schultz-Mattis theorems for SPT phases, and also derive a new LSM theorem for easy-plane spin model in a π flux lattice. We also revisit the classifications of symmetry-enriched 2D and 3D Abelian gauge theories using our results.

Robustness of numerical methods has attracted an increasing interest in the community of computational fluid dynamics. One mathematical aspect of robustness for numerical methods is the positivity-preserving property. At high Mach numbers or for flows near vacuum, solving the conservative Euler equations may generate negative density or internal energy numerically, which may lead to nonlinear instability and crash of the code. This difficulty is particularly profound for high order methods, for multi-material flows and for problems with moving meshes, such as the Lagrangian methods. In this paper, we construct both first order and uniformly high order accurate conservative Lagrangian schemes which preserve positivity of physically positive variables such as density and internal energy in the simulation of compressible multi-material flows with general equations of state (EOS). We first develop a positivity-preserving approximate Riemann solver for the Lagrangian scheme solving compressible Euler equations with both ideal and non ideal EOS. Then we design a class of high order positivity-preserving and conservative Lagrangian schemes by using the essentially non-oscillatory (ENO) reconstruction, the strong stability preserving (SSP) high order time discretizations and the positivity-preserving scaling limiter which can be proven to maintain conservation and uniformly high order accuracy and is easy to implement. One-dimensional and two-dimensional numerical tests for the positivity-preserving Lagrangian schemes are provided to demonstrate the effectiveness of these methods.