
For a spatial twisted central configuration of the Newtonian (2N+1)-body problem where 2N masses are at the vertices of two paralleled regular N-polygons with distance h>0, and the twist angle between the two regular N-polygons is 0≤θ<2π, we study the sufficient and necessary conditions for the existence of the spatial twisted central configuration. Additionally, we obtain the uniqueness of the spatial twisted central configuration.
Citation: Liang Ding, Jinrong Wang, Jinlong Wei. Spatial twisted central configuration for Newtonian (2N+1)-body problem[J]. Communications in Analysis and Mechanics, 2024, 16(2): 388-415. doi: 10.3934/cam.2024018
[1] | Muhammad Imran Asjad, Naeem Ullah, Asma Taskeen, Fahd Jarad . Study of power law non-linearity in solitonic solutions using extended hyperbolic function method. AIMS Mathematics, 2022, 7(10): 18603-18615. doi: 10.3934/math.20221023 |
[2] | Mahmoud Soliman, Hamdy M. Ahmed, Niveen Badra, M. Elsaid Ramadan, Islam Samir, Soliman Alkhatib . Influence of the β-fractional derivative on optical soliton solutions of the pure-quartic nonlinear Schrödinger equation with weak nonlocality. AIMS Mathematics, 2025, 10(3): 7489-7508. doi: 10.3934/math.2025344 |
[3] | Chun Huang, Zhao Li . New soliton solutions of the conformal time derivative generalized q-deformed sinh-Gordon equation. AIMS Mathematics, 2024, 9(2): 4194-4204. doi: 10.3934/math.2024206 |
[4] | Tianyong Han, Ying Liang, Wenjie Fan . Dynamics and soliton solutions of the perturbed Schrödinger-Hirota equation with cubic-quintic-septic nonlinearity in dispersive media. AIMS Mathematics, 2025, 10(1): 754-776. doi: 10.3934/math.2025035 |
[5] | Emad H. M. Zahran, Omar Abu Arqub, Ahmet Bekir, Marwan Abukhaled . New diverse types of soliton solutions to the Radhakrishnan-Kundu-Lakshmanan equation. AIMS Mathematics, 2023, 8(4): 8985-9008. doi: 10.3934/math.2023450 |
[6] | Kun Zhang, Tianyong Han, Zhao Li . New single traveling wave solution of the Fokas system via complete discrimination system for polynomial method. AIMS Mathematics, 2023, 8(1): 1925-1936. doi: 10.3934/math.2023099 |
[7] | Da Shi, Zhao Li, Dan Chen . New traveling wave solutions, phase portrait and chaotic patterns for the dispersive concatenation model with spatio-temporal dispersion having multiplicative white noise. AIMS Mathematics, 2024, 9(9): 25732-25751. doi: 10.3934/math.20241257 |
[8] | Mahmoud Soliman, Hamdy M. Ahmed, Niveen Badra, Taher A. Nofal, Islam Samir . Highly dispersive gap solitons for conformable fractional model in optical fibers with dispersive reflectivity solutions using the modified extended direct algebraic method. AIMS Mathematics, 2024, 9(9): 25205-25222. doi: 10.3934/math.20241229 |
[9] | Hamood Ur Rehman, Aziz Ullah Awan, Sayed M. Eldin, Ifrah Iqbal . Study of optical stochastic solitons of Biswas-Arshed equation with multiplicative noise. AIMS Mathematics, 2023, 8(9): 21606-21621. doi: 10.3934/math.20231101 |
[10] | Dumitru Baleanu, Kamyar Hosseini, Soheil Salahshour, Khadijeh Sadri, Mohammad Mirzazadeh, Choonkil Park, Ali Ahmadian . The (2+1)-dimensional hyperbolic nonlinear Schrödinger equation and its optical solitons. AIMS Mathematics, 2021, 6(9): 9568-9581. doi: 10.3934/math.2021556 |
For a spatial twisted central configuration of the Newtonian (2N+1)-body problem where 2N masses are at the vertices of two paralleled regular N-polygons with distance h>0, and the twist angle between the two regular N-polygons is 0≤θ<2π, we study the sufficient and necessary conditions for the existence of the spatial twisted central configuration. Additionally, we obtain the uniqueness of the spatial twisted central configuration.
Hydraulic jumps are frequently observed in lab studies, river flows, and coastal environments. The shallow water equations, commonly known as Saint-Venant equations, have important applications in oceanography and hydraulics. These equations are used to represent fluid flow when the depth of the fluid is minor in comparison to the horizontal scale of the flow field fluctuations [1]. They have been used to simulate flow at the atmospheric and ocean scales, and they have been used to forecast tsunamis, storm surges, and flow around constructions, among other things [2,3]. Shallow water is a thin layer of constant density fluid in hydrostatic equilibrium that is bordered on the bottom by a rigid surface and on the top by a free surface [4].
Nonlinear waves have recently acquired importance due to their capacity to highlight several complicated phenomena in nonlinear research with spirited applications [5,6,7,8]. Teshukov in [9] developed the equation system that describes multi-dimensional shear shallow water (SSW) flows. Some novel solitary wave solutions for the ill-posed Boussinesq dynamic wave model under shallow water beneath gravity were constructed [10]. The effects of vertical shear, which are disregarded in the traditional shallow water model, are included in this system to approximate shallow water flows. It is a non-linear hyperbolic partial differential equations (PDE) system with non-conservative products. Shocks, rarefactions, shear, and contact waves are all allowed in this model. The SSW model can capture the oscillatory character of turbulent hydraulic leaps, which corrects the conventional non-linear shallow water equations' inability to represent such occurrences [11]. Analytical and numerical techniques for solving such problems have recently advanced; for details, see [12,13,14,15,16,17] and the references therein.
The SSW model consists of six non-conservative hyperbolic equations, which makes its numerical solution challenging since the concept of a weak solution necessitates selecting a path that is often unknown. The physical regularization mechanism determines the suitable path, and even if one knows the correct path, it is challenging to construct a numerical scheme that converges to the weak solution since the solution is susceptible to numerical viscosity [18].
In the present work, we constructed the generalized Rusanov (G. Rusanov) scheme to solve the SSW model in 1D of space. This technique is divided into steps for predictors and correctors [19,20,21,22,23]. The first one includes a numerical diffusion control parameter. In the second stage, the balance conservation equation is retrieved. Riemann solutions were employed to determine the numerical flow in the majority of the typical schemes. Unlike previous schemes, the interesting feature of the G. Rusanov scheme is that it can evaluate the numerical flow even when the Riemann solution is not present, which is a very fascinating advantage. In actuality, this approach may be applied as a box solver for a variety of non-conservative law models.
The remainder of the article is structured as follows. Section 2 offers the non-conservative shear shallow water model. Section 3 presents the structure of the G. Rusanov scheme to solve the 1D SSW model. Section 4 shows that the G. Rusanov satisfies the C-property. Section 5 provides several test cases to show the validation of the G. Rusanov scheme versus Rusanov, Lax-Friedrichs, and reference solutions. Section 6 summarizes the work and offers some conclusions.
The one-dimensional SSW model without a source term is
∂W∂t+∂F(W)∂x+K(W)∂h∂x=S(W), | (2.1) |
W=(hhv1hv2E11E12E22),F(W)=(hv1R11+hv21+gh22R12+hv1v2(E11+R11)v1E12v1+12(R11v2+R12v1)E22v1+R12v2), |
K(W)=(000ghv112ghv20),S(W)=(0−gh∂B∂x−Cf|v|v1−Cf|v|v2−ghv1∂B∂x+12D11−Cf|v|v21−12ghv2∂B∂x+12D12−Cf|v|v1v212D22−Cf|v|v22), |
while R11 and R12 are defined by the following tensor Rij=hpij, and also, E11, E12 and E22 are defined by the following tensor Eij=12Rij+12hvivj with i≥1,j≤2. We write the previous system 2.1 in nonconservative form with non-dissipative (Cf=0,D11=D12=D22=0), see [24], which can be written as follows
∂W∂t+(∇F(W)+C(W))∂W∂x=S1(W). | (2.2) |
Also, we can write the previous system as the following
∂W∂t+A(W)∂W∂x=S1(W) | (2.3) |
with
A(W)=[010000gh00200000020−3E11v1h+2v31+ghv23E11h−3v2103v100−2E12v1h−E11v2h+2v21v2+ghv222E12h−2v1v2E11h−v21v22v10−E22v1h−2E12v2h+v22v1E22h+v222E12h−2v2v102v2v1] |
and
S1(W)=(0−gh∂B∂x0−ghv1∂B∂x−12ghv2∂B∂x0). |
System (2.3) is a hyperbolic system and has the following eigenvalues
λ1=λ2=v1,λ3=v1−√2E11h−v21,λ4=v1+√2E11h−v21, |
λ5=v1−√gh+3(v21−2E11h),λ6=v1+√gh+3(v21−2E11h). |
In order deduce the G. Rusanov scheme, we rewrite the system (2.1) as follows
∂W∂t+∂F(W)∂x=−K(W)∂h∂x+S1(W)=Q(W). | (3.1) |
Integrating Eq (3.1) over the domain [tn,tn+1]×[xi−12,xi+12] gives the finite-volume scheme
Wn+1i=Wni−ΔtΔx(F(Wni+1/2)−F(Wni−1/2))+ΔtQin | (3.2) |
in the interval [xi−1/2,xi+1/2] at time tn. Wni represents the average value of the solution W as follows
Wni=1Δx∫xi+1/2xi−1/2W(tn,x)dx, |
and F(Wni±1/2) represent the numerical flux at time tn and space x=xi±1/2. It is necessary to solve the Riemann problem at xi+1/2 interfaces due to the organization of the numerical fluxes. Assume that for the first scenario given below, there is a self-similar Riemann problem solution associated with Eq (3.1)
W(x,0)={WL,ifx<0,WR,ifx>0 | (3.3) |
is supplied by
W(t,x)=Rs(xt,WL,WR). | (3.4) |
The difficulties with discretization of the source term in (3.2) could arise from singular values of the Riemann solution at the interfaces. In order to overcome these challenges and reconstruct a Wni+1/2 approximation, we created a finite-volume scheme in [19,20,21,22,23,25,26,27] for numerical solutions of conservation laws including source terms and without source terms. The principal objective here is building the intermediate states Wni±1/2 to be utilized in the corrector stage (3.2). The following is obtained by integrating Eq (3.1) through a control volume [tn,tn+θni+1/2]×[x−,x+], that includes the point (tn,xi+1/2), with the objective to accomplish this:
∫x+x−W(tn+θni+1/2,x)dx=Δx−Wni+Δx+Wni+1−θni+1/2(F(Wni+1)−F(Wni))+θni+1/2(Δx−−Δx+)Qni+12 | (3.5) |
with Wni±1/2 denoting the approximation of Riemann solution Rs across the control volume [x−,x+] at time tn+θni+1/2. While calculating distances Δx− and Δx+ as
Δx−=|x−−xi+1/2|,Δx+=|x+−xi+1/2|, |
Qni+12 closely resembles the average source term Q and is given by the following
Qni+12=1θni+12(Δx−+Δx+)∫tn+θni+12tn∫x+x−Q(W)dtdx. | (3.6) |
Selecting x−=xi, x+=xi+1 causes the Eq (3.5) to be reduced to the intermediate state, which given as follows
Wni+1/2=12(Wni+Wni+1)−θni+1/2Δx(F(Wni+1)−F(Wni))+θni+12Qni+12, | (3.7) |
where the approximate average value of the solution W in the control domain [tn,tn+θni+1/2]×[xi,xi+1] is expressed by Wni+1/2 as the following
Wni+1/2=1Δx∫xi+1xiW(x,tn+θni+1/2)dx | (3.8) |
by selecting θni+1/2 as follows
θni+1/2=αni+1/2ˉθi+1/2,ˉθi+1/2=Δx2Sni+1/2. | (3.9) |
This choice is contingent upon the results of the stability analysis[19]. The following represents the local Rusanov velocity
Sni+1/2=maxk=1,...,K(max(∣λnk,i∣,∣λnk,i+1∣)), | (3.10) |
where αni+1/2 is a positive parameter, and λnk,i represents the kth eigenvalues in (2.3) evaluated at the solution state Wni. Here, in our case, k=6 for the shear shallow water model. One can use the Lax-Wendroff technique again for αni+1/2=ΔtΔxSni+1/2. Choosing the slope αni+1/2=˜αni+1/2, the proposed scheme became a first-order scheme, where
˜αni+1/2=Sni+1/2sni+1/2, | (3.11) |
and
sni+1/2=mink=1,...,K(max(∣λnk,i∣,∣λnk,i+1∣)). | (3.12) |
In this case, the control parameter can be written as follows:
αni+1/2=˜αni+1/2+σni+1/2Φ(rmi+1/2), | (3.13) |
where ˜αni+1/2 is given by (3.11), and Φi+1/2=Φ(ri+1/2) is a function that limits the slope. While for
ri+1/2=Wi+1−q−Wi−qWi+1−Wi,q=sgn[F′(Xn+1,Wni+1/2)] |
and
σni+1/2=ΔtΔxSni+1/2−Sni+1/2sni+1/2. |
One may use any slope limiter function, including minmod, superbee, and Van leer [28] and [29]. At the end, the G. Rusanov scheme for Eq (2.3) can be written as follows
{Wni+12=12(Wni+Wni+1)−αni+122Snj+12[F(Wni+1)−F(Wni)]+αni+122Snj+12Qni+12,Wn+1i=Wni−rn[F(Wni+12)−F(Wni−12)]+ΔtnQin. | (3.14) |
With a definition of the C-property [30], the source term in (2.1) is discretized in the G. Rusanov method in a way that is well-balanced with the discretization of the flux gradients. According to [30,31], if the following formulas hold, a numerical scheme is said to achieve the C-property for the system (3.1).
(hp11)n= constant,hn+B= constantandvn1=vn2=0. | (4.1) |
When we set v1=v2=0 in the stationary flow at rest, we get system (3.1), which can be expressed as follows.
∂∂t(h0012hp1112hp1212hp22)+∂∂x(0hp11+12gh2hp12000)=(0−gh∂B∂x0000)=Q(x,t). | (4.2) |
We can express the predictor stage (3.14) as follows after applying the G. Rusanov scheme to the previous system
Wni+12=(12(hni+hni+1)−αni+124Sni+12(hni+hni+1)[(hni+1+Bi+1)−(hni+Bi)]−αni+122Sni+12((hp11)ni+1−(hp11)ni−αni+122Sni+12((hp12)ni+1−(hp12)ni)14((hp11)ni+1+(hp11)ni)14((hp12)ni+1+(hp12)ni)14((hp22)ni+1+(hp22)ni)). | (4.3) |
Also, we can write the previous equations as follows
Wni+12=(hni+120−αni+122Sni+12((hp12)ni+1−(hp12)ni)12(hp11)ni+1212(hp12)ni+1212(hp22)ni+12), | (4.4) |
and the stage of the corrector updates the solution to take on the desired form.
((h)n+1i(hv1)n+1i(hv2)n+1i(E11)n+1i(E12)n+1i(E22)n+1i)=((h)ni(hv1)ni(hv2)ni(E11)ni(E12)ni(E22)ni)−ΔtnΔx(0g2((hni+12)2−(hni−12)2)(hp12)ni+12−(hp12)ni−12000) | (4.5) |
+(0ΔtnQni0000). |
The solution is stationary when Wn+1i=Wni and this leads us to rewrite the previous system (4.5) as follows
Δtn2Δxg((hni+12)2−(hni−12)2)−ΔtnQni. |
This lead to the following
Qni=−g8Δx(hni+1+2hni+hni−1)(Bi+1−Bi−1). | (4.6) |
Afterward, when the source term in the corrector stage is discretized in the earlier equations, this leads to the G. Rusanov satisfying the C-property.
In order to simulate the SSW system numerically, we provide five test cases without a source term using the G. Rusanov, Rusanov, and Lax-Friedrichs schemes to illustrate the effectiveness and precision of the suggested G. Rusanov scheme. The computational domain for all cases is L=[0,1] divided into 300 gridpoints, and the final time is t=0.5s, except for test case 2, where the final time is t=10s. We evaluate the three schemes with the reference solution on the extremely fine mesh of 30000 cells that was produced by the traditional Rusanov scheme. Also, we provide one test case with a source term using the G. Rusanov scheme in the domain [0,1] divided into 200 gridpoints and final time t=0.5. We select the stability condition [19] in the following sense
Δt=CFLΔxmaxi(|αni+12Sni+12|), | (5.1) |
where a constant CFL=0.5.
This test case was studied in [24] and [32], and the initial condition is given by
(h,v1,v2,p11,p22,p12)={(0.01,0.1,0.2,0.04,0.04,1×10−8)ifx≤L2,(0.02,0.1,−0.2,0.04,0.04,1×10−8)ifx>L2. | (5.2) |
The solution include five waves, which are represented by 1− shock and 6− rarefaction. We compare the numerical results with the reference solution, which is calculated with the Rusanov scheme with very fine mesh of 30000 cells. Figures 1–3 show the numerical results and variation of parameter of control. We note that all waves are captured by this scheme and the numerical solution agrees with the reference solution. Table 1 shows the CPU time of computation for the three schemes by using a Dell i5 laptop, CPU 2.5 GHZ.
Scheme | G. Rusanov | Rusanov | Lax-Friedrichs |
CPU time | 0.068005 | 0.056805 | 0.0496803 |
This test case (the shear waves problem) was studied in [24] and [32], and the initial condition is given by
(h,v1,v2,p11,p22,p12)={(0.01,0,0.2,1×10−4,1×10−4,0)ifx≤L2,(0.01,0,−0.2,1×10−4,1×10−4,0)ifx>L2. | (5.3) |
The solution is represented by two shear waves with only transverse velocity discontinuities, the p12 and p22 elements of a stress tensor. We compare the numerical results by G. Rusanov, Rusanov, and Lax-Friedrichs with the reference solution, which is calculated with the Rusanov scheme with very fine mesh of 30000 cells and a final time t=10s, see Figures 4 and 5. We observe that Figure 5 displays a spurious rise in the central of the domain of computation, which has also been seen in the literature using different methods, see [24,32].
This test case was studied in [24] with the following initial condition
(h,v1,v2,p11,p22,p12)={(0.02,0,0,4×10−2,4×10−2,0)ifx≤L2,(0.01,0,0,4×10−2,4×10−2,0)ifx>L2. | (5.4) |
The solution for h, v1, and p11 consists of two shear rarefaction waves, one moving to the right and the other moving to the left, and a shock wave to the right. See Figures 6 and 7 which show the behavior of the height of the water, velocity, pressure, and variation of the parameter of control.
We take the same example, but with another variant of the modified test case, where p12=1×10−8 is set to a small non-zero value. We notice that the behavior of h, v1, and p11 is the same, but there is a change in the behavior of p12 such that it displays all of the five waves (four rarefaction and one shock) of the shear shallow water model, see Figures 8 and 9.
This test case was studied in [24] with the following initial condition
(h,v1,v2,p11,p22,p12)={(0.02,0,0,1×10−4,1×10−4,0)ifx≤L2,(0.03,−0.221698,0.0166167,1×10−4,1×10−4,0)ifx>L2. | (5.5) |
The solution for this test case consists of a single shock wave. We compare the numerical results with the reference solution, which is calculated with the Rusanov scheme with very fine mesh of 30000 cells and a final time t=0.5s. Figures 10 and 11 show the numerical results, and we note that all waves are captured by this scheme and the numerical solution agrees with the reference solution, but the solution of p11 does not agree with the exact solution [24]. There is an oscillation with Lax-Friedrichs scheme in the water height h and the pressure p11.
This test case (the shear waves problem) was studied in [24] and [32] with the following initial condition
(h,v1,v2,p11,p22,p12)={(0.02,0,0,1×10−4,1×10−4,0)ifx≤L2,(0.01,0,0,1×10−4,1×10−4,0)ifx>L2. | (5.6) |
The solution of this test case has a single shock wave and a single rarefaction wave, separated by a contact discontinuity. We compare the numerical results with the reference solution, which is calculated with the Rusanov scheme with very fine mesh of 30000 cells and a final time t=0.5s. Figures 12 and 13 show the numerical results, and we note that all waves are captured by this scheme and the numerical solution agrees with the reference solution, but the solution of p11 does not agree with the reference solution [24].
In this test case, we take the previous test case 3 and we added the discontinuous fond topography with the following initial condition
(h,v1,v2,p11,p22,p12,B)={(0.02,0,0,4×10−2,4×10−2,1×10−8,0)ifx≤L2,(0.01,0,0,4×10−2,4×10−2,1×10−8,0.01)ifx>L2. | (5.7) |
We simulate this test case with the G. Rusanov scheme with 200 gridpoints on the domain L=[0,1] at the final time t=0.5s, and we note that the behavior of the solution for h consists of one shear rarefaction wave moving to the left, a contact discontinuity after that two shear rarefaction waves, one moving to the right and the other moving to the left (see the left side of Figure 14. The right side of the Figure 14 shows the behavior of the velocity v1, which it consists of a rarefaction wave, a contact discontinuity after that, and two shear rarefactions. The left side of Figure 14 displays the behavior of the velocity v1. It is composed of two shear rarefactions followed by a contact discontinuity and three rarefaction waves moving to the right. The the right side of Figure 14 displays the behavior of the velocity v2. It is composed of two shear rarefactions moving to the left followed by a contact discontinuity after that, and two rarefaction waves moving to the right.
In summary, we have reported five numerical examples to solve the SSW model. Namely, we implemented the G. Rusanov, Rusanov, and Lax-Friedrichs schemes. We also compared the numerical solutions with the reference solution obtained by the classical Rusanov scheme on the very fine mesh of 30000 gridpoints. The three schemes were capable of capturing shock waves and rarefaction. Also, we have given the last numerical test with a source term. Finally, we discovered that the G. Rusanov scheme was more precise than the Rusanov and Lax-Friedrichs schemes.
The current study is concerned with the SSW model, which is a higher order variant of the traditional shallow water model since it adds vertical shear effects. The model features a non-conservative structure, which makes numerical solutions problematic. The 1D SSW model was solved using the G. Rusanov technique. We clarify that this scheme satisfied the C-property. Several numerical examples were given to solve the SSW model using the G. Rusanov, Rusanov, Lax-Friedrichs, and reference solution techniques. The simulations given verified the G. Rusanov technique's high resolution and validated its capabilities and efficacy in dealing with such models. This approach may be expanded in a two-dimensional space.
H. S. Alayachi: Conceptualization, Data curation, Formal analysis, Writing - original draft; Mahmoud A. E. Abdelrahman: Conceptualization, Data curation, Formal analysis, Writing - original draft; K. Mohamed: Conceptualization, Software, Formal analysis, Writing - original draft.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The authors extend their appreciation to the Deputyship for Research & Innovation, Ministry of Education in Saudi Arabia for funding this research work through the project number 445-9-753.
The authors declare that they have no competing interests.
[1] | A. Wintner, The analytical foundations of celestial mechanics, Princeton University Press, New Jersey, 1941. |
[2] | L. Euler, De moto rectilineo trium corporum se mutuo attrahentium, Novi. Comm. Acad. Sci. Imp. Petrop., 11 (1767), 144–151. |
[3] | J. Lagrange, Essai sur le problème des trois crops, Œuvres, 6 (1772), 229–324. |
[4] |
L. Perko, E. Walter, Regular polygon solutions of the N-body problem, Proc. Amer. Math. Soc., 94 (1985), 301–309. https://doi.org/10.1090/S0002-9939-1985-0784183-1 doi: 10.1090/S0002-9939-1985-0784183-1
![]() |
[5] |
M. Corbera, C. Valls, On centered co-circular central configurations of the n-body problem, J. Dyn. Differ. Equ., 31 (2019), 2053–2060. https://doi.org/10.1007/s10884-018-9699-2 doi: 10.1007/s10884-018-9699-2
![]() |
[6] |
W. Li, Z. Q. Wang, The relationships between regular polygon central configurations and masses for Newtonian N-body problems, Phys. Lett. A, 377 (2013), 1875–1880. https://doi.org/10.1016/j.physleta.2013.05.044 doi: 10.1016/j.physleta.2013.05.044
![]() |
[7] | J. Llibre, R. Moeckel, C. Simó, Central configurations, periodic orbits, and Hamiltonian systems, Birkhäuser, New York, 2015. https://doi.org/10.1007/978-3-0348-0933-7 |
[8] |
D. G. Saari, On the role and the properties of n body central configurations, Celest. Mech., 21 (1980), 9–20. https://doi.org/10.1007/bf01230241 doi: 10.1007/bf01230241
![]() |
[9] |
S. Q. Zhang, Q. Zhou, Periodic solutions for planar 2N-body problems, Proc. Amer. Math. Soc., 131 (2003), 2161–2170. https://doi.org/10.1090/S0002-9939-02-06795-3 doi: 10.1090/S0002-9939-02-06795-3
![]() |
[10] |
R. Moeckel, C. Simó, Bifurcation of spatial central configurations from planar ones, SIAM. J. Math. Anal., 26 (1995), 978–998. https://doi.org/10.1137/S0036141093248414 doi: 10.1137/S0036141093248414
![]() |
[11] | E. Barrabés, J. M. Cors, On central configurations of the κn-body problem, J. Math. Anal. Appl., 476 (2) (2019), 720–736. https://doi.org/10.1016/j.jmaa.2019.04.010 |
[12] |
M. Corbera, J. Delgado, J. Llibre, On the existence of central configurations of p nested n-gons, Qual. Theory Dyn. Syst., 8 (2009), 255–265. https://doi.org/10.1007/s12346-010-0004-y doi: 10.1007/s12346-010-0004-y
![]() |
[13] | F. Zhao, J. Chen, Central configurations for (pN+gN)-body problems, Celest. Mech. Dyn. Astr., 121 (1) (2015), 101–106. https://10.1007/s10569-014-9593-0 |
[14] |
J. Lei, M. Santoprete, Rosette central configurations, degenerate central configurations and bifurcations, Celest. Mech. Dyn. Astr., 94 (2006), 271–287. https://doi.org/10.1007/s10569-005-5534-2 doi: 10.1007/s10569-005-5534-2
![]() |
[15] |
J. Llibre, L. F. Mello, Triple and quadruple nested central configurations for the planar n-body problem, Physica D, 238 (2009), 563–571. https://doi.org/10.1016/j.physd.2008.12.014 doi: 10.1016/j.physd.2008.12.014
![]() |
[16] |
M. Marchesin, A family of three nested regular polygon central configurations, Astrophys. Space Sci., 364 (2019), 1–12. https://doi.org/10.1007/s10509-019-3648-3 doi: 10.1007/s10509-019-3648-3
![]() |
[17] |
A. Siluszyk, On a class of central configurations in the planar 3n-body problem, Math. Comput. Sci., 11 (2017), 457–467. http://dx.doi.org/10.1007/s11786-017-0309-1 doi: 10.1007/s11786-017-0309-1
![]() |
[18] |
Z. Q. Wang, F. Y. Li, A note on the two nested regular polygonal central configurations, Proc. Amer. Math. Soc., 143 (2015), 4817–4822. https://doi.org/10.1090/S0002-9939-2015-12618-4 doi: 10.1090/S0002-9939-2015-12618-4
![]() |
[19] |
Z. F. Xie, G. Bhusal, H. Tahir, Central configurations in the planar 6-body problem forming two equilateral triangles, J. Geom. Phys., 153 (2020), 1–11. https://doi.org/10.1016/j.geomphys.2020.103645 doi: 10.1016/j.geomphys.2020.103645
![]() |
[20] |
S. Q. Zhang, C. R. Zhu, Central configurations consist of two layer twisted regular polygons, Sci. China Ser. A, 45 (2002), 1428–1438. https://doi.org/10.1007/BF02880037 doi: 10.1007/BF02880037
![]() |
[21] |
X. Yu, S. Q. Zhang, Central configurations formed by two twisted regular polygons, J. Math. Anal. Appl., 425 (2015), 372–380. https://doi.org/10.1016/j.jmaa.2014.12.023 doi: 10.1016/j.jmaa.2014.12.023
![]() |
[22] |
J. Chen, J. B. Luo, Solutions of regular polygon with an inner particle for Newtonian N+1-body problem, J. Differ. Equ., 265 (2018), 1248–1258. https://doi.org/10.1016/j.jde.2018.04.004 doi: 10.1016/j.jde.2018.04.004
![]() |
[23] | T. Ouyang, Z. F. Xie, S. Q. Zhang, Pyramidal central configurations and perverse solutions, Electron. J. Differ. Equ., 1 (2004), 1–9. |
[24] |
X. Yu, S. Q. Zhang, Twisted angles for central configurations formed by two twisted regular polygons, J. Differ. Equ., 253 (2012), 2106–2122. https://doi.org/10.1016/j.jde.2012.06.017 doi: 10.1016/j.jde.2012.06.017
![]() |
[25] | M. Marcus, H. Minc, A survey of matrix theory and matrix inequalities, Allyn & Bacon, Boston, 1964. |
[26] |
L. Ding, J. L. Wei, S. Q. Zhang, The central configuration of the planar (N+1)-body problem with a regular N-polygon for homogeneous force laws, Astrophys. Space Sci., 367 (2022), 1–9. https://doi.org/10.1007/s10509-022-04095-w doi: 10.1007/s10509-022-04095-w
![]() |
1. | Georges Chamoun, Finite Volume Analysis of the Two Competing-species Chemotaxis Models with General Diffusive Functions, 2025, 22, 2224-2902, 232, 10.37394/23208.2025.22.24 |
Scheme | G. Rusanov | Rusanov | Lax-Friedrichs |
CPU time | 0.068005 | 0.056805 | 0.0496803 |