BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Shijun Liao

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

ResearcherUnclaimed source record
nlin.CDphysics.flu-dynnlin.SIphysics.comp-phmath-phmath.MPastro-ph.EPnlin.PSOther Computer Sciencephysics.class-phphysics.gen-phquant-ph

Catalog footprint

What is connected

26works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

26 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
preprint2016arXiv

On the origin of intrinsic randomness of Rayleigh-Benard turbulence

It is of broad interest to understand how the evolution of non-equilibrium systems can be triggered and the role played by external perturbations. A famous example is the origin of randomness in the laminar-turbulence transition, which is raised in the pipe flow experiment by Reynolds as a century old unresolved problem. Although there exist different hypotheses, it is widely believed that the randomness is "intrinsic", which, however, remains as an open question to be verified. Simulating the modeled Rayleigh-Benard convection system by means of the so-called clean numerical simulation (CNS) with negligible numerical noises that are smaller even than thermal fluctuation, we verify that turbulence can be self-excited from the inherent thermal fluctuation, without any external disturbances, i.e. out of nothing. This reveals a relationship between microscopic physical uncertainty and macroscopic randomness. It is found that in physics the system nonlinearity functions as a channel for randomness information, and energy as well, to transport microscopic uncertainty toward large scales. Such scenario can generally be helpful to understand the various relevant phenomena. In methodology, compared with direct numerical simulation (DNS), CNS opens a new direction to investigate turbulent flows with largely improved accuracy and reliability.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The influence of numerical noises on statistics computation of chaotic dynamic systems

It is well known that chaotic dynamic systems (such as three-body system, turbulent flow and so on) have the sensitive dependance on initial conditions (SDIC). Unfortunately, numerical noises (such as truncation error and round-off error) always exist in practice. Thus, due to the SDIC, long-term accurate prediction of chaotic dynamic systems is practically impossible, and therefore numerical simulations of chaos are only mixtures of "true" solution with physical meanings and numerical noises without physical meanings. However, it is traditionally believed that statistic computations based on such kind of "mixtures" of numerical simulations of chaotic dynamic systems are acceptable. In this paper, using the so-called "clean numerical simulation" (CNS) whose numerical noises might be much smaller even than micro-level physical uncertainty and thus are negligible, we gain accurate prediction of a chaotic dynamic system in a long enough interval of time. Then, based on these reliable simulations, the influence of numerical noises on statistic computations is investigated. It is found that the influence of numerical noises is negligible when statistic results are time-independent. Unfortunately, when a chaotic dynamic system is far from equilibrium state so that its statistics are time-dependent, numerical noises have a great influence even on statistic computations. It suggests that even the direct numerical simulations (DNS) might not give reliable statistic computations for non-equilibrium dynamic systems with the SDIC.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the inherent self-excited macroscopic randomness of chaotic three-body system

What is the origin of macroscopic randomness (uncertainty)? This is one of the most fundamental open questions for human being. In this paper, 10000 samples of reliable (convergent), multiple-scale (from 1.0E-60 to 100) numerical simulations of a chaotic three-body system indicate that, without any external disturbance, the microscopic inherent uncertainty (in the level of 1.0E-60) due to physical fluctuation of initial positions of the three-body system enlarges exponentially into macroscopic randomness (at the level O(1)) until t=T*, the so-called physical limit time of prediction, but propagates algebraically thereafter when accurate prediction of orbit is impossible. Note that these 10000 samples use micro-level, inherent physical fluctuations of initial position, which have nothing to do with human being. Especially, the differences of these 10000 fluctuations are mathematically so small (in the level of 1.0E-60) that they are physically the SAME since a distance shorter than a Planck length does not make physical senses according to the spring theory. It indicates that the macroscopic randomness of the chaotic three-body system is self-excited, say, without any external force or disturbances, from the inherent micro-level uncertainty. This provides us the new concept "self-excited macroscopic randomness (uncertainty)". In addition, it is found that, without any external disturbance, the chaotic three-body system might randomly disrupt with the symmetry-breaking at t=1000 in about 25% probability, which provides us the new concepts "self-excited random disruption", "self-excited random escape" and "self-excited symmetry breaking" of the chaotic three-body system. It suggests that a chaotic three-body system might randomly evolve by itself, without any external forces or disturbance.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A comment on the arguments about the reliability and convergence of chaotic simulations

Yao and Hughes commented (Tellus-A, 60: 803 - 805, 2008) that "all chaotic responses are simply numerical noise and have nothing to do with the solutions of differential equations". However, using 1200 CPUs of the National Supercomputer TH-A1 and a parallel integral algorithm of the so-called "Clean Numerical Simulation" (CNS) based on the 3500th-order Taylor expansion and data in 4180-digit multiple precision, one can gain reliable, convergent chaotic solution of Lorenz equation in a rather long interval [0,10000]. This supports Lorenz's optimistic viewpoint (Tellus-A, 60: 806 - 807, 2008): "numerical approximations can converge to a chaotic true solution throughout any finite range of time".

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Observations of highly localized oscillons with multiple crests and troughs

Two types of stable, highly localized Faraday's resonant standing waves with multiple crests and troughs are observed in the alcoholic solution partly filled in a Hele-Shaw cell vertically oscillated with a single frequency. These oscillons have never been reported. Systematical experiments are performed to investigate the properties of these ocsillons. Especially, some experiments are performed, for the first time, to indicate the high localization of the oscillons, which suggest that these oscillons may be regarded as a combination of the two elementary oscillons discovered by Rajchenbach et al. (Physical Review Letters, 107, 2011), for instance, (2,3) = (1,1) + (1,2), where ($m,n$) denotes oscillon with m crests and n troughs. So, our experiments also reveal an interesting, elegant "arithmetic" of these oscillons.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On cusped solitary waves in finite water depth

It is well-known that the Camassa-Holm (CH) equation admits both of the peaked and cusped solitary waves in shallow water. However, it was an open question whether or not the exact wave equations can admit them in finite water depth. Besides, it was traditionally believed that cusped solitary waves, whose 1st-derivative tends to infinity at crest, are essentially different from peaked solitary ones with finite 1st-derivative. Currently, based on the symmetry and the exact water wave equations, Liao [1] proposed a unified wave model (UWM) for progressive gravity waves in finite water depth. The UWM admits not only all traditional smooth progressive waves but also the peaked solitary waves in finite water depth: in other words, the peaked solitary progressive waves are consistent with the traditional smooth ones. In this paper, in the frame of the linearized UWM, we further give, for the first time, the cusped solitary waves in finite water depth, and besides reveal a close relationship between the cusped and peaked solitary waves: a cusped solitary wave is consist of an infinite number of peaked solitary ones with the same phase speed, so that it can be regarded as a special peaked solitary wave. This also well explains why and how a cuspon has an infinite 1st-derivative at crest. It is found that, like peaked solitary waves, the vertical velocity of a cusped solitary wave in finite water depth is also discontinuous at crest (x=0), and especially its phase speed has nothing to do with wave height, too. In addition, it is unnecessary to consider whether the peaked/cusped solitary waves given by the UWM are weak solution or not, since the governing equation is not necessary to be satisfied at crest. All of these would deepen and enrich our understandings about the cusped solitary waves.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the mathematically reliable long-term simulation of chaos of Lorenz equation in the interval [0,10000]

Using 1200 CPUs of the National Supercomputer TH-A1 and a parallel integral algorithm based on the 3500th-order Taylor expansion and the 4180-digit multiple precision data, we have done a reliable simulation of chaotic solution of Lorenz equation in a rather long interval [0,10000] (Lorenz time unit). Such a kind of mathematically reliable chaotic simulation has never been reported. It provides us a numerical benchmark for mathematically reliable long-term prediction of chaos. Besides, it also proposes a safe method for mathematically reliable simulations of chaos in a finite but long enough interval. In addition, our very fine simulations suggest that such a kind of mathematically reliable long-term prediction of chaotic solution might have no physical meanings, because the inherent physical micro-level uncertainty due to thermal fluctuation might quickly transfer into macroscopic uncertainty so that trajectories for a long enough time would be essentially uncertain in physics.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A commend on "Three Classes of Newtonian Three-Body Planar Periodic Orbits" by Šuvakov and Dmitrašinović (PRL, 2013)

Currently, the fifteen new periodic solutions of Newtonian three-body problem with equal mass were reported by Šuvakov and Dmitrašinović (PRL, 2013) [1]. However, using a reliable numerical approach (namely the Clean Numerical Simulation, CNS) that is based on the arbitrary-order Taylor series method and data in arbitrary-digit precision, it is found that at least seven of them greatly depart from the periodic orbits after a long enough interval of time. Therefore, the reported initial conditions of at least seven of the fifteen orbits reported by Šuvakov and Dmitrašinović [1] are not accurate enough to predict periodic orbits. Besides, it is found that these seven orbits are unstable.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Chaotic motion of three-body problem : an origin of macroscopic randomness of the universe

The famous three-body problem is investigated by means of a numerical approach with negligible numerical noises in a long enough time interval, namely the Clean Numerical Simulation (CNS). From physical viewpoints, position of any bodies contains inherent micro-level uncertainty. The evaluations of such kind of inherent micro-level uncertainty are accurately simulated by means of the CNS. Our reliable, very accurate CNS results indicate that the inherent micro-level uncertainty of position of a star/planet might transfer into macroscopic randomness. Thus, the inherent micro-level uncertainty of a body might be an origin of macroscopic randomness of the universe. In addition, from physical viewpoints, orbits of some three-body systems at large time are inherently random, and thus it has no physical meanings to talk about the accurate long-term prediction of the chaotic orbits. Note that such kind of uncertainty and randomness has nothing to do with the ability of human being. All of these might enrich our knowledge and deepen our understandings about not only the three-body problem but also chaos.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Do peaked solitary water waves indeed exist?

Many models of shallow water waves admit peaked solitary waves. However, it is an open question whether or not the widely accepted peaked solitary waves can be derived from the fully nonlinear wave equations. In this paper, a unified wave model (UWM) based on the symmetry and the fully nonlinear wave equations is put forward for progressive waves with permanent form in finite water depth. Different from traditional wave models, the flows described by the UWM are not necessarily irrotational at crest, so that it is more general. The unified wave model admits not only the traditional progressive waves with smooth crest, but also a new kind of solitary waves with peaked crest that include the famous peaked solitary waves given by the Camassa-Holm equation. Besides, it is proved that Kelvin's theorem still holds everywhere for the newly found peaked solitary waves. Thus, the UWM unifies, for the first time, both of the traditional smooth waves and the peaked solitary waves. In other words, the peaked solitary waves are consistent with the traditional smooth ones. So, in the frame of inviscid fluid, the peaked solitary waves are as acceptable and reasonable as the traditional smooth ones. It is found that the peaked solitary waves have some unusual and unique characteristics. First of all, they have a peaked crest with a discontinuous vertical velocity at crest. Especially, the phase speed of the peaked solitary waves has nothing to do with wave height. In addition, the kinetic energy of the peaked solitary waves either increases or almost keeps the same from free surface to bottom. All of these unusual properties show the novelty of the peaked solitary waves, although it is still an open question whether or not they are reasonable in physics if the viscosity of fluid and surface tension are considered.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Two new standing solitary waves in shallow water

In this paper, the closed-form analytic solutions of two new Faraday's standing solitary waves due to the parametric resonance of liquid in a vessel vibrating vertically with a constant frequency are given for the first time. Using a model based on the symmetry of wave elevation and the linearized Boussinesq equation, we gain the closed-form wave elevations of the two kinds of non-monotonically decaying standing solitary waves with smooth crest and the even or odd symmetry. All of them have never been reported, to the best of our knowledge. Besides, they can well explain some experimental phenomena. All of these are helpful to deepen and enrich our understandings about standing solitary waves and Faraday's wave.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On peaked solitary waves of Camassa-Holm equation

Unlike the Boussinesq, KdV and BBM equations, the celebrated Casamma-Holm (CH) equation can model both phenomena of soliton interaction and wave breaking. Especially, it has peaked solitary waves in case of omega=0. Besides, in case of omega > 0, its solitary wave "becomes $C^\infty$ and there is no derivative discontinuity at its peak", as mentioned by Camassa and Holm in 1993 (PRL). However, it is found in this article that the CH equation has peaked solitary waves even in case of omega > 0. Especially, all of these peaked solitary waves have an unusual property: their phase speeds have nothing to do with the height of peakons or anti-peakons. Therefore, in contrast to the traditional view-points, the peaked solitary waves are a common property of the CH equation: in fact, all mainstream models of shallow water waves admit such kind of peaked solitary waves

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On peaked solitary waves of Degasperis - Procesi equation

The Degasperis - Procesi (DP) equation describing the propagation of shallow water waves contains a physical parameter $ω$, and it is well-known that the DP equation admits solitary waves with a peaked crest when $ω= 0$. In this article, we illustrate, for the first time, that the DP equation admits peaked solitary waves even when $ω\neq 0$. This is helpful to enrich our knowledge and deepen our understandings about peaked solitary waves of the DP equation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Two kinds of peaked solitary waves of the KdV, BBM and Boussinesq equations

It is well-known that the celebrated Camassa-Holm equation has the peaked solitary waves, which have been not reported for other mainstream models of shallow water waves. In this letter, the closed-form solutions of peaked solitary waves of the KdV equation, the BBM equation and the Boussinesq equation are given for the first time. All of them have either a peakon or an anti-peakon. Each of them exactly satisfies the corresponding Rankine-Hogoniot jump condition and should be understood as weak solution. Therefore, the peaked solitary waves might be common for most of shallow water wave models, no matter whether or not they are integrable and/or admit breaking-wave solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the homotopy multiple-variable method and its applications in the interactions of nonlinear gravity waves

The basic ideas of a homotopy-based multiple-variable method is proposed and applied to investigate the nonlinear interactions of periodic traveling waves. Mathematically, this method does not depend upon any small physical parameters at all and thus is more general than the traditional multiple-scale perturbation techniques. Physically, it is found that, for a fully developed wave system, the amplitudes of all wave components are finite even if the wave resonance condition given by Phillips (1960) is exactly satisfied. Besides, it is revealed that there exist multiple resonant waves, and that the amplitudes of resonant wave may be much smaller than those of primary waves so that the resonant waves sometimes contain rather small part of wave energy. Furthermore, a wave resonance condition for arbitrary numbers of traveling waves with large wave amplitudes is given, which logically contains Phillips' four-wave resonance condition but opens a way to investigate the strongly nonlinear interaction of more than four traveling waves with large amplitudes. This work also illustrates that the homotopy multiple-variable method is helpful to gain solutions with important physical meanings of nonlinear problems, if the multiple variables are properly defined with clear physical meanings.

0 citations0 reviewsNo code linkedOpen access