BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Shijun Liao

Shijun Liao contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
nlin.CDphysics.comp-phphysics.flu-dynmath-phmath.MPOther Computer Sciencephysics.class-phphysics.gen-phquant-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

11 published item(s)

preprint2026arXiv

Scaling transition in horizontal convection near the density maximum

Horizontal convection (HC) serves as a canonical model for geophysical and industrial flows driven by differential heating along a surface. While the classical Oberbeck-Boussinesq (OB) approximation is well-established, the impact of a nonlinear equation of state, specifically the density maximum of water near $4^\circ\mathrm{C}$, remains underexplored. This study investigates Non-Oberbeck-Boussinesq (NOB) effects on HC via direct numerical simulations (DNS) over a Rayleigh number range of $10^6 \le Ra \le 5\times 10^{10}$. We examine two configurations: Classical HC (CHC) and Symmetric HC (SHC). Our results reveal that the NOB-SHC case undergoes a structural transition, evolving from a bicellular structure to a full-depth, single-roll circulation driven by central `mixing plumes'. This reorganization manifests as transitional anomalies in Reynolds number ($Re$) scaling, whereas the emergence of full-depth plumes fundamentally alters the heat transport mechanism. Consequently, unlike the classical Rossby scaling ($Nu \sim Ra^{1/5}$) observed in reference cases, the NOB-SHC regime exhibits an enhanced heat transport scaling ranging from $Nu \sim Ra^{1/4}$ to $Ra^{1/3}$. To rationalize this behavior, we extend the Shishkina-Grossmann-Lohse (SGL) theory by incorporating a generalized potential energy transfer term ($Φ_{i2}$). The theoretical framework demonstrates that the global scaling law is dictated by the characteristic plume height ($\hat{z}$). Specifically, when plumes penetrate the entire cavity depth ($\hat{z} \sim H$), as observed in the NOB-SHC case, the flow transcends classical bounds for OB HC, accessing a regime analogous to Rayleigh Bénard convection. The proposed theory successfully unifies the scaling laws for both OB and NOB fluids, showing excellent agreement with numerical data.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Avoiding small denominator problems by means of the homotopy analysis method

The so-called ``small denominator problem'' was a fundamental problem of dynamics, as pointed out by Poincaré. Small denominators appear most commonly in perturbative theory. The Duffing equation is the simplest example of a non-integrable system exhibiting all problems due to small denominators. In this paper, using the forced Duffing equation as an example, we illustrate that the famous ``small denominator problems'' never appear if a non-perturbative approach based on the homotopy analysis method (HAM), namely ``the method of directly defining inverse mapping'' (MDDiM), is used. The HAM-based MDDiM provides us great freedom to directly define the inverse operator of an undetermined linear operator so that all small denominators can be completely avoided and besides the convergent series of multiple limit-cycles of the forced Duffing equation with high nonlinearity are successfully obtained. So, from the viewpoint of the HAM, the famous ``small denominator problems'' are only artifacts of perturbation methods. Therefore, completely abandoning perturbation methods but using the HAM-based MDDiM, one would be never troubled by ``small denominators''. The HAM-based MDDiM has general meanings in mathematics and thus can be used to attack many open problems related to the so-called ``small denominators''.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Large-scale influence of numerical noises as artificial stochastic disturbances on a sustained turbulence

We investigate the large-scale influence of numerical noises as tiny artificial stochastic disturbances on a sustained turbulence. Using the two-dimensional (2D) turbulent Rayleigh-Bénard (RB) convection as an example, we numerically solve the NS equations, separately, by means of a traditional algorithm with double precision (marked by RKwD) and the so-called clean numerical simulation (CNS). The numerical simulation given by the RKwD is a mixture of the "true" physical solution and the "false" numerical noises that is random and can be regarded as a kind of artificial stochastic disturbances: unfortunately, the "true" physical solution is mostly at the same level as the "false" numerical noises. By contrast, the CNS can greatly reduce the background numerical noise to any a required level so that the "false" numerical noises are negligible compared with the "true" physical solution and thus the CNS solution can be used as a "clean" benchmark solution for comparison. It is found that the numerical noises as tiny artificial stochastic disturbances could indeed lead to large-scale deviations of simulations not only in spatio-temporal trajectories but also even in statistics. Especially, these numerical noises (as artificial stochastic disturbances) even lead to different types of flows: the shearing convection occurs for the RKwD simulations, and its corresponding flow field turns to a kind of zonal flow thereafter, however the CNS benchmark solution always sustains the non-shearing vortical/roll-like convection during the whole process of simulation. Thus, we provide a rigorous evidence that numerical noises as a kind of small-scale artificial stochastic disturbances have quantitatively and qualitatively large-scale influences on a sustained turbulence, i.e. the 2D turbulent RB convection considered in this paper.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Three-body problem -- from Newton to supercomputer plus machine learning

The famous three-body problem can be traced back to Newton in 1687, but quite few families of periodic orbits were found in 300 years thereafter. In this paper, we propose an effective approach and roadmap to numerically gain planar periodic orbits of three-body systems with arbitrary masses by means of machine learning based on an artificial neural network (ANN) model. Given any a known periodic orbit as a starting point, this approach can provide more and more periodic orbits (of the same family name) with variable masses, while the mass domain having periodic orbits becomes larger and larger, and the ANN model becomes wiser and wiser. Finally we have an ANN model trained by means of all obtained periodic orbits of the same family, which provides a convenient way to give accurate enough predictions of periodic orbits with arbitrary masses for physicists and astronomers. It suggests that the high-performance computer and artificial intelligence (including machine learning) should be the key to gain periodic orbits of the famous three-body problem.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Accurate predictions of chaotic motion of a free fall disk

It is important to know the accurate trajectory of a free fall object in fluid (such as a spacecraft), whose motion might be chaotic in many cases. However, it is impossible to accurately predict its chaotic trajectory in a long enough duration by traditional numerical algorithms in double precision. In this paper, we give the accurate predictions of the same problem by a new strategy, namely the Clean Numerical Simulation (CNS). Without loss of generality, a free fall disk in water is considered, whose motion is governed by the Andersen-Pesavento-Wang model. We illustrate that convergent and reliable trajectories of a chaotic free fall disk in a long enough interval of time can be obtained by means of the CNS, but different traditional algorithms in double precision give disparate trajectories. Besides, unlike the traditional algorithms in double precision, the CNS can predict the accurate posture of the free fall disk near the vicinity of the bifurcation point of some physical parameters in a long duration. Therefore, the CNS can provide reliable prediction of chaotic systems in a long enough interval of time.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Ultra-chaos: an insurmountable objective obstacle of reproducibility and replicability

In this paper, a new concept, i.e. ultra-chaos, is proposed for the first time. Unlike a normal-chaos, statistical properties such as the probability density functions (PDF) of an ultra-chaos are sensitive to tiny disturbances. We illustrate that ultra-chaos is widely existed and thus has general scientific meanings. It is found that statistical non-reproducibility is an inherent property of an ultra-chaos so that an ultra-chaos is at a higher-level of disorder than a normal-chaos. Thus, it is impossible in practice to replicate experimental/numerical results of an ultra-chaos even in statistical meanings, since random environmental noises always exist and are out of control. Thus, the ultra-chaos should be an insurmountable obstacle of reproducibility and replicability. Similar to Gödel's incompleteness theorem, such kind of "incompleteness of reproducibility" reveals a limitation of our traditional scientific paradigm based on reproducible experiments, which can be traced back to Galileo. The ultra-chaos opens a new door and possibility to study chaos theory, turbulence theory, computational fluid dynamics (CFD), the statistical significance, reproducibility crisis, and so on.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the risks of using double precision in numerical simulations of spatio-temporal chaos

Due to the butterfly-effect, computer-generated chaotic simulations often deviate exponentially from the true solution, so that it is very hard to obtain a reliable simulation of chaos in a long-duration time. In this paper, a new strategy of the so-called Clean Numerical Simulation (CNS) in physical space is proposed for spatio-temporal chaos, which is computationally much more efficient than its predecessor (in spectral space). The strategy of the CNS is to reduce both of the truncation and round-off errors to a specified level by implementing high-order algorithms in multiple-precision arithmetic (with sufficient significant digits for all variables and parameters) so as to guarantee that numerical noise is below such a critical level in a temporal interval $t\in[0,T_c]$ that corresponding numerical simulation remains reliable over the whole interval. Without loss of generality, the complex Ginzburg-Landau equation (CGLE) is used to illustrate its validity. As a result, a reliable long-duration numerical simulation of the CGLE is achieved in the whole spatial domain over a long interval of time $t\in[0,3000]$, which is used as a reliable benchmark solution to investigate the influence of numerical noise by comparing it with the corresponding ones given by the 4th-order Runge-Kutta method in double precision (RKwD). Our results demonstrate that the use of double precision in simulations of chaos might lead to huge errors in the prediction of spatio-temporal trajectories and in statistics, not only quantitatively but also qualitatively, particularly in a long interval of time.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

One family of 13315 stable periodic orbits of the non-hierarchical unequal-mass triple system

The three-body problem has been studied for more than three centuries [1,2], and has received much attention in recent years [3-5]. It shows complex dynamical phenomena due to the mutual gravitational interaction of the three bodies. Triple systems are common in astronomy, but all observed periodic triple systems are hierarchical up till now [6-8]. It is traditionally believed that bound non-hierarchical triple systems are almost unstable and disintegrate into a stable binary system and a single star [5], and thus stable periodic orbits of non-hierarchical triple systems are rather scarce. Here we report one family of 13315 stable periodic orbits of the non-hierarchical triple system with unequal mass. Compared with the narrow mass region (only 10E-5) of the stable figure-eight solution [9], our newly-found stable periodic orbits can have fairly large mass region. It is found that many of these newly-found stable periodic orbits have the mass ratios close to those of the hierarchical triple systems that have been measured by the astronomical observation. It implies that these stable periodic orbits of the non-hierarchical triple system with distinctly unequal masses can be quite possibly observed in practice. Our investigation also suggests that there should exist an infinite number of stable periodic orbits of non-hierarchical triple systems with distinctly unequal masses. Obviously, these stable periodic orbits of the non-hierarchical unequal-mass triple system have broad impact for the astrophysical scenario: they could inspire the theoretical and observational study of the non-hierarchical triple system, the formation of triple stars [6], the gravitational waves pattern [10] and the gravitational waves observation [11] of the non-hierarchical triple system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Satellites of the Broucke-Hadjidemetriou-Hénon family of periodic unequal-mass three-body orbits

The Broucke-Hadjidemetriou-Hénon's (BHH) orbits are a family of periodic orbits of the three-body system with the simplest topological free group word $a$, while the BHH satellites have free group words $a^k$ ($k>1$), where $k$ is the topological exponent. Janković and Dmitrašinović [Phy. Rev. Lett. 116, 064301 (2016)] reported 57 new BHH satellites with equal mass and found that at a fixed energy the relationship between the angular momentum ($L$) and the topologically rescaled period ($T/k$) is the same for both of the BHH orbits ($k=1$) and the BHH satellites ($k>1$). In this letter, we report 419,743 new BHH orbits ($k=1$) and 179,253 new BHH satellites ($k>1$) of the three-body system with unequal mass, which have never been reported, to the best of our knowledge. Among these newly-found 598,996 BBH orbits and satellites, about 33.5 % (i.e., 200,686) are linearly stable and thus many among them might be observed in practice. Besides, we discover that, for the three-body system with unequal mass at a fixed energy, relationship between the angular momentum ($L$) and topologically rescaled period ($T/k$) of the BHH satellites ($k>1$) is different from that of the BHH orbits ($k=1$).

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

A new non-perturbative approach in quantum mechanics for time-independent Schrödinger equations

A new non-perturbative approach is proposed to solve time-independent Schrödinger equations in quantum mechanics and chromodynamics (QCD). It is based on the homotopy analysis method (HAM), which was developed by the author for highly nonlinear equations since 1992 and has been widely applied in many fields. Unlike perturbative methods, this HAM-based approach has nothing to do with small/large physical parameters. Besides, convergent series solution can be obtained even if the disturbance is far from the known status. A nonlinear harmonic oscillator is used as an example to illustrate the validity of this approach for disturbances that might be more than hundreds larger than the possible superior limit of the perturbative approach. This HAM-based approach could provide us rigorous theoretical results in quantum mechanics and chromodynamics (QCD), which can be directly compared with experimental data. Obviously, this is of great benefit not only for improving the accuracy of experimental measurements but also for validating physical theories.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Collisionless periodic orbits in the free-fall three-body problem

Although the free-fall three-body problem have been investigated for more than one century, however, only four collisionless periodic orbits have been found. In this paper, we report 234 collisionless periodic orbits of the free-fall three-body system with some mass ratios, including three known collisionless periodic orbits. Thus, 231 collisionless free-fall periodic orbits among them are entirely new. In theory, we can gain periodic orbits of the free-fall three-body system in arbitrary ratio of mass. Besides, it is found that, for a given ratio of masses of two bodies, there exists a generalized Kepler's third law for the periodic three-body system. All of these would enrich our knowledge and deepen our understanding about the famous three-body problem as a whole.

0 citations0 reviewsNo code linkedOpen access