Source author record

Chaoqun Liu

Chaoqun Liu 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

14works
4topics
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

14 published item(s)

preprint2025arXiv

A Relative Liutex Method for Vortex Identification

A relative Liutex vortex identification method is proposed in this study, together with its explicit mathematical formulation. The method is designed to identify vortical structures based solely on local flow-field information and is inherently Galilean invariant, ensuring robustness under different reference frames. To validate the proposed approach, a three-dimensional flat-plate boundary-layer transition case is investigated, in which the relative Liutex is systematically compared with conventional vortex identification methods, including the Q-criterion and the original Liutex method. The results show that the relative Liutex is capable of simultaneously capturing both strong and weak vortical structures. Importantly, its behavior cannot be interpreted as a simple superposition of Liutex iso-surfaces obtained using different threshold values. Instead, the relative Liutex provides a more selective and physically coherent identification of weak vortices, particularly in regions above the $Λ$-vortex and in the downstream hairpin-vortex structures, while effectively suppressing spurious and noise-induced features. These advantages arise from its formulation based on local velocity-gradient strength rather than a global vortical-intensity measure. Owing to its ability to consistently identify vortices across a wide range of intensities, the relative Liutex demonstrates strong potential for revealing complex vortex structures and underlying flow mechanisms in vortex-dominated flows.

preprint2022arXiv

An objective vortex identification method

Generally, the vortex structures should be independent of the observers who are moving, especially when their coordinates are non-inertial, which may result in confusions in communications between researchers. The property that not being influenced by the choice of coordinate is called objective. Mostly, researchers would like to gather data in the inertial coordinate because most physical laws are valid in the inertial coordinate. As a result, how to transfer the observations in a non-inertial coordinate to those in an inertial coordinate raises interest. In this paper, methods are provided to obtain velocity and velocity gradient tensor in an inertial coordinate from the non-inertial observations based on one known zero-vorticity point which can be chosen from physical fact such as the point is located in the inviscid region.Further, objective vortex structure can be calculated from the objective velocity or velocity gradient tensor. With Liutex chosen as the indicator of the vortex, two numerical examples are used to test the proposed methods. The results show that the methods perform very well and the relative error is extremely small.

preprint2020arXiv

Correlation Analysis among Vorticity, Q method and Liutex

Influenced by the fact that vorticity represents rotation for rigid body, people believe it also works for fluid flow. However, the theoretical predictions by vorticity do not match experiment results, which drove scientists to look for better methods to describe vortex. According to Dr. Liu classification, all methods applied to detect vortex can be categorized into three generations. The vorticity-based method is classified as the first generation. Methods relying on eigenvalues of velocity gradient tensor are considered as the second generation. Although so many methods appeared, people still believe vorticity is vortex since vorticity theory looks perfect in math, and all other methods are only scalars and unable to indicate swirl direction. Recently, Dr. Liu innovated a new vortex identification method called Liutex. Liutex, a vector quantity, which is regarded as the third-generation method, not only overcomes all previous methods drawbacks, but also has a clear physical meaning. The direction of Liutex represents the swirl axis of rotation, and its strength is equal to twice of angular speed. In this paper, we did a correlation analysis between vorticity, Q, Lambda ci, Lambda 2 methods and Liutex based on a DNS case of boundary layer transition. The results show that the correlation between vorticity and Liutex is minimal in strong shear region, which demonstrates the idea that using vorticity to detect vortex lacks a scientific foundation; in other words, vorticity is not vortex.

preprint2020arXiv

Liutex-based Modified Navier-Stokes Equation

The Navier-Stokes (NS) partial differential equations, as the governing equation of fluid dynamics, are based on particles with zero volume. In such a condition, conservation law of moment of momentum is automatically satisfied, and thus NS equations only contain conservation of mass, momentum and energy. Because of the difficulty to get an analytical solution of NS equations, scientists develop finite element method (FEM) and finite volume method (FVM) to obtain approximate solutions. However, these methods are dependent on the size of finite volumes, which is not zero. Considering the size of the finite volume, conservation law of moment of momentum is no longer automatically satisfied and it is necessary to develop a new control equation to take effect of momentum into account. In this paper, new relations of reciprocal shear stresses are derived from conservation law of moment of momentum, and so are new constitutive relations and modified NS equations. The Liutex based NS equation may be the universal governing equations for both laminar and turbulent flow which is dominated by vortices and conservation of moment of momentum must be satisfied.

preprint2020arXiv

Stretching and Shearing Contamination Analysis for Liutex/Rortex and Other Vortex Identification Methods

Although traditional vortex identification methods such as Q, Delta, Lambda2, Lambdaci remain popular in the identification and visualization of vortices, these methods count on shearing and stretching as a part of vortex strength. However, shearing and stretching do not contribute to fluid rotation. In this paper, the contamination effects of stretching and shearing of these methods are investigated and compared with Liutex method. From our investigation, the Liutex is an exact definition of fluid rotation or vortex, while other vortex identification methods are contaminated by stretching and shearing at different levels. The decomposition of the velocity gradient tensor can only be conducted in a so-called Principal Coordinate for uniqueness. The mathematical relation between Liutex and other vortex identification function are derived in this paper and then the effects of shearing and stretching on different vortex identification methods are studied. The mathematical formula and computation of the stretching and shearing effects on different schemes clearly show that the Liutex method has superiority over the other vortex identification methods as it only counts on local fluid rigid rotation while other methods count on stretching/compression and shearing as a part of the fluid rotation or vortex.

preprint2014arXiv

DNS Study for the origin of the flow Randomization in Late Boundary Layer Transition

This paper is devoted to the investigation of the origin and mechanism of randomization in late boundary layer transition over a flat plate without pressure gradient. The flow randomization is a crucial phase before flow transition to the turbulent state. According to existing literatures, the randomization was caused by the big background noises and non-periodic spanwise boundary conditions. It was assumed that the large ring structure is affected by background noises first, and then the change of large ring structure affects the small length scales quickly, which directly leads to randomization and formation of turbulence. However, by careful analysis of our high order DNS results, we believe that the internal instability of multiple ring cycles structure is the main reason. What we observed is that randomization begins when the third cycle overlaps the first and second cycles. A significant asymmetric phenomenon is originated from the second cycle in the middle of both streamwise and spanwise directions. More technically, a visible asymmetric phenomenon in the middle vortex ring cycle starts at time step t=16.25T and x=838.9δin where the top and bottom level rings are still completely symmetric. The non-symmetric structure of middle level ring affects the small length scale in boundary layer bottom quickly. The randomization phenomenon spreads to top level through ejections. Finally, the whole flow domain becomes randomized. A hypothesis of C- and K-types shift is given as a possible mechanism of flow randomization.

preprint2014arXiv

DNS Study on Physics of Late Boundary Layer Transition

This paper serves as a review of our recent new DNS study on physics of late boundary layer transition. This includes mechanism of the large coherent vortex structure formation, small length scale generation and flow randomization. The widely spread concept vortex breakdown to turbulence,which was considered as the last stage of flow transition, is not observed and is found theoretically incorrect. The classical theory on boundary layer transition is challenged and we proposed a new theory with five steps, i.e. receptivity, linear instability, large vortex formation, small length scale generation, loss of symmetry and randomization to turbulence. We have also proposed a new theory about turbulence generation. The new theory shows that all small length scales (turbulence) are generated by shear layer instability which is produced by large vortex structure with multiple level vortex rings, multiple level sweeps and ejections, and multiple level negative and positive spikes near the laminar sub-layers.Therefore, turbulence is not generated by vortex breakdown but rather positive and negative spikes and consequent high shear layers. Shear layer instability is considered as the mother of turbulence. This new theory may give a universal mechanism for turbulence generation and sustenance which shows that the energy is brought by large vortex structure through multiple level sweeps.

preprint2014arXiv

Improvement of Mixing Function for Modified Upwinding Compact Scheme

The compact scheme has high order accuracy and high resolution, but cannot be used to capture the shock. WENO is a great scheme for shock capturing, but is too dissipative for turbulence and small length scales. We developed a modified upwinding compact scheme which uses an effective shock detector to block compact scheme to cross the shock and a control function to mix the flux with WENO scheme near the shock. The new scheme makes the original compact scheme able to capture the shock sharply and, more important, keep high order accuracy and high resolution in the smooth area, which is particularly important for shock boundary layer and shock acoustic interactions. This work is a continuation to modify the control function for the modified up-winding compact scheme (MUCS). Numerical results show the scheme is successful for 2-D Euler.

preprint2014arXiv

LES Study on the Mechanism of Vortex Rings behind Supersonic MVG with Turbulent Inflow

In this study, we investigate the interaction between vortex rings behind MVG and the oblique shocks in the MVG controlled ramp flow at M=2.5 and Req=5760. Implicit large eddy simulation (ILES) method is used by solving the unfiltered form of the Navier-Stokes equations with the 5th order Bandwidth-optimized WENO scheme. The fully developed inflow is given by a series of turbulent profiles obtained from previous DNS simulation. It shows that the ring structure does not break down and keeps its topology after penetrating the strong shock wave and the oblique shocks is influenced a lot by the induced flow field from rings. The bump of the 3D shock wave surface is discovered and its mechanism is explained.

preprint2014arXiv

Numerical and Experimental Studies on the Separation Topology of the MVG Controlled Flow at M=2.5

In this paper, the implicitly implemented LES method and fifth order bandwidth optimized WENO scheme are used to make comprehensive studies on the separation topology of the MVG controlled flow at M=2.5 and Reθ=5760. Experiments are also made to verify the prediction of the computation. Analyses are conducted on three categories of the topology: the surface separation, cross-section separation and the three dimensional structure of the vortices. A complete description about the separation topology and a series of new findings are obtained. Among them, a pair of spiral point is first predicted by the computation and verified by the experiment. A corresponding new vortex model with 7 vortex tubes is presented also.

preprint2014arXiv

Numerical Study of Mechanism of U-shaped Vortex Formation

This paper illustrates the mechanism of U-shaped vortex formation which is found both by experiment and DNS. The main goal of this paper is to explain how the U-shaped vortex is formed and further develops. According to the results obtained by our direct numerical simulation with high order accuracy, the U-shaped vortex is part of the coherent vortex structure and is actually the tertiary streamwise vortices induced by the secondary vortices. The new finding is quite different from existing theories which describe that the U-shaped vortex is newly formed as the head of young turbulence spot and finally break down to small pieces. In addition, we find that the U-shaped vortex has the same vorticity sign as the original λ-shaped vortex tube legs and serves as a second neck to supply vorticity to the ringlike vortex when the original vortex tube is stretched and multiple rings are generated.

preprint2014arXiv

Numerical Study on Mechanism of Multiple Rings Formation

In this paper, the flow around each ring-like vortex is investigated by high order DNS including first sweep, first ejection, second sweep, second ejection, positive spike, momentum deficit, vortex shape, vortex location, strength of sweeps, etc. Meanwhile, the mechanism about formation of momentum deficit is deeply studied. A new mechanism on how the multiple rings are formed one by one found both by experiment and by DNS in late boundary layer transition is presented. It also reveals that the relation between streamwise vortex and spanvise vertex rings, and how the vorticity is transferred between them.

preprint2014arXiv

Numerical Study on Mechanism of Small Vortex Generation in Boundary Layer Transition

The small vortex generation is a key issue of the mechanism for late flow transition and turbulence generation. It was widely accepted that small length vortices were generated by large vortex breakdown. According to our recent DNS, we find that the hairpin vortex structure is very stable and never breaks down to small pieces. On the other hand, we recognize that there are strong positive spikes besides the ring neck in the spanwise direction. The strongly positive spikes are caused by second sweeps which are generated by perfectly circular and perpendicularly standing vortex rings. The second sweep brings energy from the invisid region downdraft to the bottom of the boundary layers, which generates high shear layers around the positive spikes.Since the high shear layer is not stable, all small length scales (turbulence) are generated around high shear layers especially near the wall surface (bottom of boundary layers). This happens near the ring neck in the streamwise direction and besides the original vortex legs in the spanwise direction.The small length scales then rise up from the wall surface and travel to the upper boundary layer.

preprint2014arXiv

Numerical Study on Randomization in Late Boundary Layer Transition

The mechanism of randomization in late boundary layer transition is a key issue of late boundary layer transition and turbulence theory. We studied the mechanism carefully by high order DNS. The randomization was originally considered as a result of large background noise and non-periodic spanwise boundary conditions. It was addressed that the large ring structure is affected by background noises first and then the change of large ring structure affects the small length scale quickly, which directly leads to randomization and formation of turbulence. However, what we observed is that the loss of symmetry starts from the middle level rings while the top and bottom rings are still symmetric. The nonsymmetric structure of second level rings will influence the small length scales at the boundary layer bottom quickly. The symmetry loss at the bottom of the boundary layer is quickly spread to up level through ejections. This will lead to randomization of the whole flow field. Therefore, the internal instability of multiple level ring structure, especially the middle ring cycles, is the main reason for flow randomization, but not the background noise. A hypothesis is given that the loss of symmetry may be caused by the shift from C-type transition to K-type transition or reverses.