Source author record

Jiequan Li

Jiequan Li appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2021arXiv

On a supersonic-sonic patch arising from the Frankl problem in transonic flows

We construct a supersonic-sonic smooth patch solution for the two dimensional steady Euler equations in gas dynamics. This patch is extracted from the Frankl problem in the study of transonic flow with local supersonic bubble over an airfoil. Based on the methodology of characteristic decompositions, we establish the global existence and regularity of solutions in a partial hodograph coordinate system in terms of angle variables. The original problem is solved by transforming the solution in the partial hodograph plane back to that in the physical plane. Moreover, the uniform regularity of the solution and the regularity of an associated sonic curve are also verified.

preprint2021arXiv

One-sided GRP Solver and Numerical Boundary Conditions for compressible fluid flows

In the computation of compressible fluid flows, numerical boundary conditions are always necessary for all physical variables at computational boundaries while just partial physical variables are often prescribed as physical boundary conditions. Certain extrapolation technique or ghost cells are often employed traditionally for this issue but spurious wave reflections often arise to cause numerical instability. In this paper, we associate this issue with the one-sided generalized Riemann problem (GRP) solver motivated by the accelerated piston problem in gas dynamics so that the extrapolation technique can be actually avoided. In fact, the compatibility arguments naturally requires to formulate the one-sided generalized Riemann problem and incorporate it into the numerical procedure of boundary conditions. As far as the interaction of nonlinear waves with physical boundaries, such a one-sided GRP solver shows significant effects, as numerical experiments demonstrate, on avoiding spurious wave reflections at the computational boundaries.

preprint2021arXiv

Two-stage Fourth-order Gas Kinetic Solver-based Compact Subcell Finite Volume Method for Compressible Flows over Triangular Meshes

To meet the demand for complex geometries and high resolutions of small-scale flow structures, a two-stage fourth-order subcell finite volume (SCFV) method combining the gas-kinetic solver (GKS) with subcell techniques for compressible flows over (unstructured) triangular meshes was developed to improve the compactness and efficiency. Compared to the fourth-order GKS-based traditional finite volume (FV) method, the proposed method realizes compactness effectively by subdividing each cell into a set of subcells or control volumes (CVs) and selecting only face-neighboring cells for high-order compact reconstruction. Because a set of CVs share a solution polynomial, the reconstruction is more efficient than that for traditional FV-GKS, where each CV needs to be separately reconstructed. Unlike in the single-stage third-order SCFV-GKS, both accuracy and efficiency are improved significantly by two-stage fourth-order temporal discretization, for which only a second-order gas distribution function is needed to simplify the construction of the flux function and reduce computational costs. For viscous flows, it is not necessary to compute the viscous term with GKS. Compared to the fourth-stage Runge--Kutta method, one half of the stage is saved for achieving fourth-order time accuracy, which also helps to improve the efficiency. Therefore, a new high-order method with compactness, efficiency, and robustness is proposed by combining the SCFV method with the two-stage gas-kinetic flux. Several benchmark cases were tested to demonstrate the performance of the method in compressible flow simulations.

preprint2020arXiv

A spacetime outlook on CFD: Spacetime correlated models and spacetime coupled algorithms

A spacetime outlook on Computational Fluid Dynamics is advocated: models in fluid mechanics often have the spacetime correlation property, which should be inherited and preserved in the corresponding numerical algorithms. Starting from the fundamental formulation of fluid mechanics under continuum hypothesis, this paper defines the meaning of spacetime correlation of the models, establishes the fundamental principle of finite volume schemes, expounds the necessity of spacetime coupling of algorithms, as well as realizes the physical and mathematical unification of basic governing equations of fluid mechanics and finite volume schemes. In practice, the design methodology of spacetime coupling high order numerical algorithms is presented, and the difference from spacetime decoupling method is compared. It should be pointed out that most of the contents in this paper are suitable for computational fluid dynamics under the assumption of continuous medium, and some are only suitable for compressible flow.

preprint2020arXiv

Consistency and Convergence of Finite Volume Approximations to Nonlinear Hyperbolic Balance Laws

This paper addresses the three concepts of \textit{ consistency, stability and convergence } in the context of compact finite volume schemes for systems of nonlinear hyperbolic conservation laws. The treatment utilizes the framework of "balance laws". Such laws express the relevant physical conservation laws in the presence of discontinuities. Finite volume approximations employ this viewpoint, and the present paper can be regarded as being in this category. It is first shown that under very mild conditions a weak solution is indeed a solution to the balance law. The schemes considered here allow the computation of several quantities per mesh cell (e.g., slopes) and the notion of consistency must be extended to this framework. Then a suitable convergence theorem is established, generalizing the classical convergence theorem of Lax and Wendroff. Finally, the limit functions are shown to be entropy solutions by using a notion of "Godunov compatibility", which serves as a substitute to the entropy condition.

preprint2020arXiv

Regularity of Fluxes in Nonlinear Hyperbolic Balance Laws

This paper addresses the issue of the formulation of weak solutions to systems of nonlinear hyperbolic conservation laws as integral balance laws. The basic idea is that the "meaningful objects" are the fluxes, evaluated across domain boundaries over time intervals. The fundamental result in this treatment is the regularity of the flux trace in the multi-dimensional setting. It implies that a weak solution indeed satisfies the balance law. In fact, it is shown that the flux is Lipschitz continuous with respect to suitable perturbations of the boundary.

preprint2016arXiv

A Few Benchmark Test Cases for Higher-order Euler Solvers

There have been great efforts on the development of higher-order numerical schemes for compressible Euler equations. The traditional tests mostly targeting on the strong shock interactions alone may not be adequate to test the performance of higher-order schemes. This study will introduce a few test cases with a wide range of wave structures for testing higher-order schemes. As reference solutions, all test cases will be calculated by our recently developed two-stage fourth-order gas-kinetic scheme (GKS). All examples are selected so that the numerical settings are very simple and any high order accurate scheme can be straightly used for these test cases, and compare their performance with the GKS solutions. The examples include highly oscillatory solutions and the large density ratio problem in one dimensional case; hurricane-like solutions, interactions of planar contact discontinuities (the composite of entropy wave and vortex sheets) sheets with large Mach number asymptotic, interaction of planar rarefaction waves with transition from continuous flows to the presence of shocks, and other types of interactions of two-dimensional planar waves. The numerical results from the fourth-order gas-kinetic scheme provide reference solutions only. These benchmark test cases will help CFD developers to validate and further develop their schemes to a higher level of accuracy and robustness.

preprint2016arXiv

An Efficient and Accurate Two-Stage Fourth-order Gas-kinetic Scheme for the Navier-Stokes Equations

For computational fluid dynamics (CFD), the generalized Riemann problem (GRP) solver and the gas-kinetic kinetic scheme (GKS) provide a time-accurate flux function starting from a discontinuous piecewise linear flow distributions around each cell interface. With the use of time derivative of the flux function, a two-stage Lax-Wendroff-type (L-W for short) time stepping method has been recently proposed in the design of a fourth-order time accurate method [18]. In this paper, based on the same time-stepping method and the second-order GKS flux function [34], a fourth-order gas-kinetic scheme is constructed for the Euler and Navier-Stokes equations. In comparison with the formal one-stage time-stepping third-order gas-kinetic solver [21], the current fourth-order method not only reduces the complexity of the flux function, but also improves the accuracy of the scheme, even though the third- and fourth-order schemes have similar computation cost. Most importantly, the robustness of the fourth-order GKS is as good as the second-order one. Perfect numerical solutions can be obtained from the high Reynolds number boundary layer solutions to the hypersonic viscous heat conducting flow computations. Many numerical tests, including many difficult ones for the Navier-Stokes solvers, have been used to validate the current fourth-order method. Following the two-stage time-stepping framework, the one-stage third-order GKS can be easily extended to a fifth-order method with the usage of both first-order and second-order time derivatives of the flux function. The use of time-accurate flux function may have great impact on the development of higher-order CFD methods.

preprint2015arXiv

A Two-Stage Fourth Order Time-Accurate Discretization for Lax-Wendroff Type Flow Solvers. I. Hyperbolic Conservation Laws

In this paper we develop a novel two-stage fourth order time-accurate discretization for time-dependent flow problems, particularly for hyperbolic conservation laws. Different from the classical Runge-Kutta (R-K) temporal discretization for first order Riemann solvers as building blocks, the current approach is solely associated with Lax-Wendroff (L-W) type schemes as the building blocks. As a result, a two-stage procedure can be constructed to achieve a fourth order temporal accuracy, rather than using well-developed four stages for R-K methods. The generalized Riemann problem (GRP) solver is taken as a representative of L-W type schemes for the construction of a two-stage fourth order scheme.

preprint2013arXiv

The generalized Riemann problems for hyperbolic balance laws: A unified formulation towards high order

The Generalized Riemann Problems (GRP) for nonlinear hyperbolic systems of balance laws in one space dimension are now well-known and can be formulated as follows: Given initial-data which are smooth on two sides of a discontinuity, determine the time evolution of the solution near the discontinuity. While the classical Riemann problem serves as a primary building block in the construction of many numerical schemes (most notably the Godunov scheme), the analytic study of GRP will lead to an array of GRP schemes, which extend the Godunov scheme. Currently there are extensive studies on the second-order GRP scheme, which proves to be robust and is capable of resolving complex multidimensional fluid dynamic problems [M. Ben-Artzi and J. Falcovitz, "Generalized Riemann Problems in Computational Fluid Dynamics", Cambridge University Press, 2003]. A more general formulation of the second-order GRP solver is still confined with a class of weakly coupled systems [Numer. Math. (2007) 106:369-425]. This paper provides a unified approach for solving the GRP in the general context of hyperbolic balance laws, without weakly coupled constraint, towards high order accuracy. The derivation of the second-order GRP solver is more concise compared to those in previous works and the third-order quadratic GRP is resolved for the first time. The latter is shown to be necessary through numerical experiments with strong discontinuities. Our method relies heavily on the new treatment of the rarefaction wave by deriving the L(Q)-equations, an ODE system capturing the "evolution" of the characteristic derivatives in x-t space for generalized Riemann invariants. The case of a sonic point is incorporated into a general treatment. The accuracy of the derived GRP solvers are justified and numerical examples are presented for the performance of the resulting scheme.