Cover Image

A simple, robust, and accurate HLLC-type Riemann solver for the compressible Euler equations at various Mach numbers is built. To cure shock instability of the HLLC solver at strong shocks, a pressure-control technique, which plays a role in limiting the propagation of erroneous pressure perturbation, is proposed. With an all Mach correction method for the compressible Euler system, the proposed method is further extended to compute flow problems at low Mach numbers. The proposed all Mach HLLC-type scheme has been implemented and used to compute a variety of flow problems ranging from hypersonic compressible to low Mach incompressible flow regimes. Various numerical results demonstrate that the obtained all Mach HLLC-type scheme is both accurate and stable for all speed ranges.     

On para‐Kähler Lie algebroids and contravariant pseudo‐Hessian structures

In this paper, we generalize all the results obtained on para‐Kähler Lie algebras in [3] to para‐Kähler Lie algebroids. In particular, we study exact para‐Kähler Lie algebroids as a generalization of exact para‐Kähler Lie algebras. This study leads to a natural generalization of pseudo‐Hessian manifolds, we call them contravariant pseudo‐Hessian manifolds. Contravariant pseudo‐Hessian manifolds have many similarities with Poisson manifolds. We explore these similarities which, among others, leads to a powerful machinery to build examples of non trivial pseudo‐Hessian structures. Namely, we will show that given a finite dimensional commutative and associative algebra , the orbits of the action Φ of on given by are pseudo‐Hessian manifolds, where . We illustrate this result by considering many examples of associative commutative algebras and show that the resulting pseudo‐Hessian manifolds are very interesting.     

A high‐order finite difference method for numerical simulations of supersonic turbulent flows

This paper describes a new class of three‐dimensional finite difference schemes for high‐speed turbulent flows in complex geometries based on the high‐order monotonicity‐preserving (MP) method. Simulations conducted for various 1D, 2D, and 3D problems indicate that the new high‐order MP schemes can preserve sharp changes in the flow variables without spurious oscillations and are able to capture the turbulence at the smallest computed scales. Our results also indicate that the MP method has less numerical dissipation and faster grid convergence than the weighted essentially non‐oscillatory method. However, both of these methods are computationally more demanding than the COMP method and are only used for the inviscid fluxes. To reduce the computational cost for reacting flows, the scalar equations are solved by the COMP method, which is shown to yield similar results to those obtained by the MP in supersonic turbulent flows with strong shock waves. Copyright © 2011 John Wiley & Sons, Ltd.     

Development of circular function‐based gas‐kinetic scheme (CGKS) on moving grids for unsteady flows through oscillating cascades

In this paper, the circular function‐based gas‐kinetic scheme (CGKS), which was originally developed for simulation of flows on stationary grids, is extended to solve moving boundary problems on moving grids. Particularly, the unsteady flows through oscillating cascades are our major interests. The main idea of the CGKS is to discretize the macroscopic equations by the finite volume method while the fluxes at the cell interface are evaluated by locally reconstructing the solution of the continuous Boltzmann Bhatnagar–Gross–Krook equation. The present solver is based on the fact that the modified Boltzmann equation, which is expressed in a moving frame of reference, can recover the corresponding macroscopic equations with Chapman–Enskog expansion analysis. Different from the original Maxwellian function‐based gas‐kinetic scheme, in improving the computational efficiency, a simple circular function is used to describe the equilibrium state of distribution function. Considering that the concerned cascade oscillating problems belong to cases that the motion of surface boundary is known a priori, the dynamic mesh method is suitable and is adopted in the present work. In achieving the mesh deformation with high quality and efficiency, a hybrid dynamic mesh method named radial basic functions‐transfinite interpolation is presented and applied for cascade geometries. For validation, several numerical test cases involving a wide range are investigated. Numerical results show that the developed CGKS on moving grids is well applied for cascade oscillating flows. And for some cases where nonlinear effects are strong, the solution accuracy could be effectively improved by using the present method. Copyright © 2017 John Wiley & Sons, Ltd.     

A Few Remarks About the Hilbert Scheme of Smooth Projective Curves

We discuss two conjectures by Francesco Severi and Joe Harris about the irreducibility and the dimension of the Hilbert scheme parameterizing smooth projective curves of given degree and genus.     

Construction of Biorthogonal Multiwavelets Using the Lifting Scheme

The lifting scheme has been found to be a flexible method for constructing scalar wavelets with desirable properties. Here it is extended to the construction of multiwavelets. It is shown that any set of compactly supported biorthogonal multiwavelets can be obtained from the Lazy matrix filters with a finite number of lifting steps. As an illustration of the general theory, compactly supported biorthogonal multiwavelets with optimum time–frequency resolution are constructed. In addition, experimental results of applying these multiwavelets to image compression are presented.     

Developing a unified FVE‐ALE approach to solve unsteady fluid flow with moving boundaries

In this study, an arbitrary Lagrangian–Eulerian (ALE) approach is incorporated with a mixed finite‐volume–element (FVE) method to establish a novel moving boundary method for simulating unsteady incompressible flow on non‐stationary meshes. The method collects the advantages of both finite‐volume and finite‐element (FE) methods as well as the ALE approach in a unified algorithm. In this regard, the convection terms are treated at the cell faces using a physical‐influence upwinding scheme, while the diffusion terms are treated using bilinear FE shape functions. On the other hand, the performance of ALE approach is improved by using the Laplace method to improve the hybrid grids, involving triangular and quadrilateral elements, either partially or entirely. The use of hybrid FE grids facilitates this achievement. To show the robustness of the unified algorithm, we examine both the first‐ and the second‐order temporal stencils. The accuracy and performance of the extended method are evaluated via simulating the unsteady flow fields around a fixed cylinder, a transversely oscillating cylinder, and in a channel with an indented wall. The numerical results presented demonstrate significant accuracy benefits for the new hybrid method on coarse meshes and where large time steps are taken. Of importance, the current method yields the second‐order temporal accuracy when the second‐order stencil is used to discretize the unsteady terms. Copyright © 2009 John Wiley & Sons, Ltd.     

On the stable difference scheme for source identification nonlocal elliptic problem

Approximation of source identification problem for elliptic equation with integral-type nonlocal condition is discussed. The first order of accuracy difference scheme for elliptic nonlocal identification problem is studied. By using spectral resolution of a self-adjoint operator, we establish stability inequalities for solution of constructed scheme. Subsequently, the difference scheme for approximate solution of multidimensional boundary value problem with integral-type nonlocal and first kind boundary conditions is investigated on stability. Numerical test examples are presented.