Asymptotic problems and numerical schemes for traffic flows with unilateral constraints describing the formation of jams

  • Received: 01 December 2016 Revised: 01 August 2017
  • Primary: 65M08; Secondary: 35L60, 90B20, 35Q91

  • We discuss numerical strategies to deal with PDE systems describing traffic flows, taking into account a density threshold, which restricts the vehicle density in the situation of congestion. These models are obtained through asymptotic arguments. Hence, we are interested in the simulation of approached models that contain stiff terms and large speeds of propagation. We design schemes intended to apply with relaxed stability conditions.

    Citation: Florent Berthelin, Thierry Goudon, Bastien Polizzi, Magali Ribot. Asymptotic problems and numerical schemes for traffic flows with unilateral constraints describing the formation of jams[J]. Networks and Heterogeneous Media, 2017, 12(4): 591-617. doi: 10.3934/nhm.2017024

    Related Papers:

    [1] Marimuthu Mohan Raja, Velusamy Vijayakumar, Anurag Shukla, Kottakkaran Sooppy Nisar, Wedad Albalawi, Abdel-Haleem Abdel-Aty . A new discussion concerning to exact controllability for fractional mixed Volterra-Fredholm integrodifferential equations of order r(1,2) with impulses. AIMS Mathematics, 2023, 8(5): 10802-10821. doi: 10.3934/math.2023548
    [2] Sumbal Ahsan, Rashid Nawaz, Muhammad Akbar, Saleem Abdullah, Kottakkaran Sooppy Nisar, Velusamy Vijayakumar . Numerical solution of system of fuzzy fractional order Volterra integro-differential equation using optimal homotopy asymptotic method. AIMS Mathematics, 2022, 7(7): 13169-13191. doi: 10.3934/math.2022726
    [3] Clement Manga, Alioune Coulibaly, Alassane Diedhiou . On some stochastic differential equations with jumps subject to small positives coefficients. AIMS Mathematics, 2019, 4(5): 1369-1385. doi: 10.3934/math.2019.5.1369
    [4] Subramanian Muthaiah, Manigandan Murugesan, Muath Awadalla, Bundit Unyong, Ria H. Egami . Ulam-Hyers stability and existence results for a coupled sequential Hilfer-Hadamard-type integrodifferential system. AIMS Mathematics, 2024, 9(6): 16203-16233. doi: 10.3934/math.2024784
    [5] Manal Elzain Mohamed Abdalla, Hasanen A. Hammad . Solving functional integrodifferential equations with Liouville-Caputo fractional derivatives by fixed point techniques. AIMS Mathematics, 2025, 10(3): 6168-6194. doi: 10.3934/math.2025281
    [6] Chao Wei . Parameter estimation for partially observed stochastic differential equations driven by fractional Brownian motion. AIMS Mathematics, 2022, 7(7): 12952-12961. doi: 10.3934/math.2022717
    [7] Thabet Abdeljawad, Sabri T. M. Thabet, Imed Kedim, Miguel Vivas-Cortez . On a new structure of multi-term Hilfer fractional impulsive neutral Levin-Nohel integrodifferential system with variable time delay. AIMS Mathematics, 2024, 9(3): 7372-7395. doi: 10.3934/math.2024357
    [8] Mouffak Benchohra, Zohra Bouteffal, Johnny Henderson, Sara Litimein . Measure of noncompactness and fractional integro-differential equations with state-dependent nonlocal conditions in Fréchet spaces. AIMS Mathematics, 2020, 5(1): 15-25. doi: 10.3934/math.2020002
    [9] Guowei Zhang . The lower bound on the measure of sets consisting of Julia limiting directions of solutions to some complex equations associated with Petrenko's deviation. AIMS Mathematics, 2023, 8(9): 20169-20186. doi: 10.3934/math.20231028
    [10] Saima Rashid, Fahd Jarad, Khadijah M. Abualnaja . On fuzzy Volterra-Fredholm integrodifferential equation associated with Hilfer-generalized proportional fractional derivative. AIMS Mathematics, 2021, 6(10): 10920-10946. doi: 10.3934/math.2021635
  • We discuss numerical strategies to deal with PDE systems describing traffic flows, taking into account a density threshold, which restricts the vehicle density in the situation of congestion. These models are obtained through asymptotic arguments. Hence, we are interested in the simulation of approached models that contain stiff terms and large speeds of propagation. We design schemes intended to apply with relaxed stability conditions.



    Mathematical modeling helps to develop a scientific understanding, quantify the current knowledge of a system, and test the impact of changes in a system. It is used in natural sciences and engineering disciplines. So, the mathematical modeling has been attracting the attention of many authors [1,2,3,4].

    Worldwide energy consumption is increasing by around 2.1 percent annually on average [5]. 84 percent of the fuels produced are produced with Fossil fuels, Coal and Natural Gas. Increasing electricity consumption increases CO2 emissions and increases global warming.

    According to 2020 Eurostat data in [6], 26 percent of energy consumption in Europe in 2018 was at homes. 14.1 percent of the total electricity consumption was used for white goods and lighting. The rate of use of dishwashers in Europe is over 44 percent. Dishwashers are among the white goods with a high potential to become widespread. Using a dishwasher provides 50 percent water savings and 28 percent energy savings compared to hand washing [7,8]. The annual energy consumption of dishwashers in Europe reached 31TWH/year in 2015. Efficiency increases in dishwashers have an important place in preventing the global water and energy crisis.

    In the studies about the dishwasher, the factors affecting the efficiency are determined. One of them is to wash with an intermittent water jet instead of applying a conventional continuous water jet from the lower spray arm [9]. Also, the author has investigate the lower energy consumption and the increase in washing performance depending on the operation case. In [10], Tsouknidas and Zhang have presented the relationships among water pressure, rotational speed, and washing performance using motor-controlled spray arm. The nozzles on the standard spray arm are fixed angled. Since these fixed angle nozzles provide rotation, the rotation speed of the spray arm is directly related to the angle of the nozzles. When higher pressure is applied to increase the washing performance, the rotation speed increases, but when the speed increases, the washing performance decreases and energy consumption increases. Since nozzle angle, water pressure, rotational speed are determined according to the priority of maintaining the washing performance in the conventional spray arm, then energy consumption and water consumption remain in the background. The new spray arm design provides a great improvement in energy and water consumption as it stretches conventional design criteria [10].

    Arçelik company has developed a spray arm that sweep the inside of the dishwashers in all four corners. When the spray arm is viewed as a projection, since it sweep a square profile, it can send water vertically to every region. In this way, it can exhibit lower energy, water pressure, sound level and higher washing performance. This new system is called as CornerWash spray arm. In this study, we have obtained a mathematical model of the new spray arm motion. The path followed by the new spray arm corresponds to the spherical timelike helix motion in Minkowski 3space. We examine the new spray arm motion as a two-dimensional projection, we can see that it makes a hypocycloid motion.

    Minkowski 3space is defined with the following inner product

    g(u,v)=u0v0+u1v1+u2v2,

    for the vectors u=(u0,u1,u2) and v=(v0,v1,v2) in the scalar space R31. The exterior product in R31 is defined by

    u×v=|e0e1e2u0u1u2v0v1v2|, (2.1)

    where {ei,i=0,1,2} is canonical bases of R31 [11]. In Minkowski 3space, vectors and curves have casual characters and they are called spacelike, timelike, or lightlike respect to their casual characters. If the observer is moving less than the speed of light along the curve, then the curve is called a timelike curve; if it is moving at a speed of light called a lightlike curve; and if it is moving faster than the speed of light called a spacelike curve. Similarly, the norm of the vector in Minkowski 3space is equal to 1, 0, or 1. Then, it is called a timelike, a lightlike, or a spacelike, respectively. Also, the spheres are defined as follows in Minkowski 3space:

    S21={pR31|g(p,p)=1}de Sitter space,H20={pR31|g(p,p)=1}hyperbolic space,Q2={pR31|g(p,p)=0}lightcone.

    The curve x is a non-null curve with the arc length parameter s in the Minkowski 3space. The Frenet frame apparatus of the curve x are {t(s),n(s),b(s),ϰ(s),τ(s)}, where n(s), b(s) are unit normal and binormal vector field of the curve x, respectively. The derivatives of the Frenet vectors are given in terms of themselves as follows:

    t(s)=κ(s)n(s),n(s)=ε0ε1κ(s)t(s)+τ(s)b(s),b(s)=ε1ε2τ(s)n(s), (2.2)

    where εi (i=0,1,2) is shown signum functions of the Frenet vectors t(s), n(s), and b(s), respectively. Also, the cross product of these vectors is given as b=t×n, ε1t=n×b, and ε0n=b×t.

    Helix curve is one of the important curves in science and nature. The applications of helices are encountered in many applied fields. In differential geometry, a helix is defined as a geometric curve with a non-vanishing constant curvature and a non-vanishing constant torsion. Another definition is that the tangent makes a constant angle with a fixed straight line (the axis of the general helix). Also, slant helix is defined that the normal makes a constant angle with a fixed straight line. In [12], the authors have given parametric representations of the slant helices on de Sitter 2space and hyperbolic 2space.

    In this section, we have examined the geometrical movements of the conventional and new dishwasher spray arms. We give the relationships between the motions of dishwasher spray arms and helices.

    Theorem 1. The unit speed curve x lies on de Sitter 2space with radius r if and only if the curvatures of the curve x are saisfied the equation ε1κ2(s)+ε2[1τ(s)(1κ(s))]2=r2.

    Proof. If the curve x(s) lies on the de Sitter 2space with radius r and center a, then g(x(s)a,x(s)a)=r2 is satisfied. From the last equation, we calculate the followings g(x(s)a,t(s))=0, g(x(s)a,n(s))=ε0κ(s), and g(x(s)a,b(s))=ε0τ(s)(1κ(s)). The position vector of the curve x(s) is expressed by x(s)a=a1(s)t(s)+a2(s)n(s)+a3(s)b(s). By using these equations, we obtain a1(s)=0, a2(s)=ε0ε1κ(s), and a3(s)=ε0ε2τ(s)(1κ(s)). If these equations are replaced in g(x(s)a,x(s)a)=r2, the desired result is obtained.

    Assume that the equation ε1κ2(s)+ε2τ(s)(1κ(s))=r2 is provided. The following curve α is defined as α(s)=x(s)+ε0ε1Rn(s)+ε0ε2TRb(s) where T=1/τ(s) and R=1/κ(s). Now, we take the derivative of α(s), we hold α(s)=0 and so α is a constant. The radius of the pseudo sphere is calculated as g(x(s)α,x(s)α)=r2. This demonstrates that x(s) lies on the de Sitter 2space with radius r.

    Theorem 2. Let x(s):IRR31 be spherical timelike helix on the de Sitter 2space. The parametric representation of the spherical timelike helix is

    x(s)={cossm21,mcosssin(m21sm)m21+sinscos(m21sm),mcosscos(m21sm)m21+sinssin(m21sm)},

    for mR[1,1].

    Proof. The first derivative of the curve x(s) is obtained ˙x(s)=sinsm21{1,sin(m21sm)m,cos(m21sm)m}, and g(˙x(s),˙x(s))=sin2sm2. Since the norm of the vector ˙x(s) is not equal to 1, then the curve x(s) is not a unit speed curve. The tangent vector field t(s) of the curve x with an arbitrary parameter s is calculated by

    t(s)={mm21,1m21sin(m21sm),1m21cos(m21sm)},

    so, we get ˙x(s)=sinsmt(s). The second derivative of the curve x(s) is found

    ¨x(s)={cossm21,cos(m21sm)sins+mm21cosssin(m21sm)m2,sin(m21sm)sinsmm21cosscos(m21sm)m2}.

    To find the binormal vector field b(s)=˙x(s)רx(s)˙x(s)רx(s), we calculate

    ˙x(s)רx(s)={sin2sm3m21,sin2ssin(m21sm)m2m21,sin2scos(m21sm)m2m21},

    and the norm ˙x(s)רx(s)=sin2sm3. Thus, we get

    b(s)={1m21,mm21sin(m21sm),mm21cos(m21sm)}.

    The normal vector field of the timelike curve x(s) is

    n(s)={0,cos(m21sm),sin(m21sm)},

    with the curvature κ(s)=1sins and the torsion τ(s)=msins. Since the ratio of the curvatures is m=constant, then the curve x is a helix. Since g(x(s),x(s))=1 is satisfied, then the curve lies on the de Sitter space S21.

    The natural equation of the curve is the equation that characterizes the curve independent of the choice of coordinates. It is said to be intrinsic because it is invariant under transformations that preserve length and angle. The study of natural equations began with the following problem: given the curvature and torsion of one parameter, find the space curve for two functions. We will now consider the natural equations of spherical helix on the de Sitter 2space projects on a plane perpendicular to its axis.

    Theorem 3. The projection of the spherical timelike helix on a plane perpendicular to its axis is the plane curve with the natural equations

    |ε1R2±r2|=cot2θs2 (3.1)

    where R(s)=1κ(s).

    Proof. A space curve x(s) lies on the de Sitter 2space with radius r if and only if the condition ε1κ2(s)+ε2[1τ(s)(1κ(s))]2=r2 is satisfied. This condition with considering the equation κ(s)=ε1τ(s)tanθ give the following ODE after elimination of τ(s)

    ε1R2+ε2(RR)2tan2θ=r2,

    where the parameter θ is showen the angle between the tangent vector of the helix and its axis. Integrated for R, and by suitable choice of the additive constant in s, gives the desired results. The projection of the spherical timelike helix on a plane perpendicular to their axis is |ε1R2±r2|=cot2θs2. This type of curves is discussed for Euclidean 2sphere in [13].

    Corollary 1. The spherical timelike helix x(s) on the de Sitter 2space projects on a plane perpendicular to its axis in an arc of hypocycloid.

    Proof. The natural equation of the projection on yzplane perpendicular to its axis of the spherical timelike helix x(s) is obtained the following equation cot2θs2+R2(s)=1. From the last equation, we obtain cotθ=m. Since m>1, projection of the curve x(s) on a plane perpendicular to its axis is an hypocycloid.

    The conventional dishwasher spray arm, simply, is in the form of a rotating line at its midpoint. When viewed from the top of the machine, dishwashers have a square area, while conventional spray arm sweep a circular area in a square structure. This structure geometrically corresponds to a circular cylinder. When the movement of the water coming out of the nozzles is examined depending on time, it makes a cylindrical helix movement with the parametric representation α(s)=(coss,sins,ks) where k is constant (red curve in Figure 1 right). Therefore, the spray arm creates regions where the water cannot reach vertically. For water to reach these areas, it is necessary to make angled nozzles on the conventional spray arm. When the water comes at an angle to the area where outside the projection of the conventional spray arm, it loses some of the force required to clean it due to the angled exit. To reduce this loss of force, water must be sprayed with greater pressure from the angled nozzle. So, the pump must be more powerful and consume more energy for high pressure. Also, the amount of water must be higher in order to remove the dirt. For these reasons, the conventional spray arm increases energy and water consumption, while reducing washing performance. As the water coming out of the spray arm with the angled nozzle hits the inner walls of the dishwasher, it also increases the amount of noise.

    Figure 1.  Projection of the paths of the nozzles in the conventional dishwasher spray arm (left) and its geometric interpretation is the cylindrical helix (right).

    New designed spray arm called the CornerWash by Arçelik in Figure 2, since it sweeps projectively the square area of the dishwasher, it is not necessary to use an angled nozzle. Therefore, the water coming out of the nozzle with a lower pressure will be able to provide the necessary force in every region to clean the dishes. Moreover, as the force on the dishes increases, the amount of water to be consumed decreases and the washing performance increases. It is possible to use a smaller pump, thus reducing energy consumption. Since the water does not hit the side walls directly, the operating noise level of the dishwasher is reduced.

    Figure 2.  The projection of the paths of the nozzles in the CornerWash spray arm (left) and its geometric interpretation is the spherical timelike helix for m=1.155 (right).

    The movement of CornerWash begins with the rotation of the orange arm around the red dot on the axis of rotation. The axis of rotation of the three-armed spray arm is off-axis on the orange lever. The axis of rotation of the three-armed spray arm moves on the white circle. Thus, the nozzles on the three-armed spray arm make hypocycloid movement in Figure 3.

    Figure 3.  The movement of the nozzles on the CornerWash spray arm.

    When the movement of CornerWash is examined geometrically, the path of the water coming out of the nozzle on the new spray arm projectively draws a hypocycloid. When the time-dependent variation of this motion in 3-dimensional space is examined, we find that it corresponds to the timelike helix with

    x(s)={cossm21,mcosssin(m21sm)m21+sinscos(m21sm),mcosscos(m21sm)m21+sinssin(m21sm)}

    lying on the de Sitter 2space in Figure 4 for m=1.155 (right). The projection on yzplane perpendicular to its axis of the spherical timelike helix x(s) is obtained cot2θs2+sin2s=1. From the last equation, we calculate cotθ=1.155. Since cotθ>1, this equation gives a hypocycloid.

    Figure 4.  The paths of the nozzles on the CornerWash spray arm (left) and its geometric interpretation is the spherical timelike helix (right).

    Note: CornerWash images are provided from the website: https://www.youtube.com/watch?v=8ertBCEdhlU.

    In this article, the authors created a geometric model of the water flow of a CornerWash arm design used in dishwashers. This design improves the washing performance and sound level while providing efficiency increase in the water and energy consumed by the dishwashers. The study is important because it includes the geometric interpretation of a movement that has an application. Since it can be used as an infrastructure for Fluid Dynamics, it will help scientists doing similar work.

    The authors declare no conflicts of interest in this paper.

    [1] A. Aw, A. Klar, T. Materne and M. Rascle, Derivation of continuum traffic flow models from microscopic follow-the-leader models, SIAM J. Appl. Math., 63 (2002), 259-278, doi: 10.1137/S0036139900380955
    [2] Resurrection of "second order" models of traffic flow. SIAM J. Appl. Math. (2000) 60: 916-938.
    [3] On the multiscale modeling of vehicular traffic: From kinetic to hydrodynamics. Disc. Cont. Dyn. Syst.-B (2014) 19: 1869-1888.
    [4] On the modeling of traffic and crowds: A survey of models, speculations, and perspectives. SIAM Rev. (2011) 53: 409-463.
    [5] S. Benzoni-Gavage and D. Serre, Multidimensional Hyperbolic Partial Differential Equations: First-Order Systems and Applications, Oxford Mathematical Monographs, Oxford University Press, 2007.
    [6] A model for the evolution of traffic jams in multilane. Kinetic and Related Models (2012) 5: 697-728.
    [7] Multifluid flows: A kinetic approach. J. Sci. Comput. (2016) 66: 792-824.
    [8] A model for the formation and evolution of traffic jams. Arch. Rational Mech. Anal. (2008) 187: 185-220.
    [9] F. Bouchut, Nonlinear Stability of Finite Volume Methods for Hyperbolic Conservation Laws and Well-Balanced Schemes for Sources, Frontiers in math., Birkhäuser, 2004. doi: 10.1007/b93802
    [10] Numerical approximations of pressureless and isothermal gas dynamics. SIAM J. Numer. Anal. (2003) 41: 135-158 (electronic).
    [11] Sticky particles and scalar conservation laws. SIAM J. Numer. Anal. (1998) 35: 2317-2328.
    [12] Transport-equilibrium schemes for computing contact discontinuities in traffic flow modeling. Commun. Math. Sci. (2007) 5: 533-551.
    [13] Godunov scheme and sampling technique for computing phase transitions in traffic flow modeling. Interfaces Free Bound. (2008) 10: 197-221.
    [14] Glimm's method for gas dynamics. SIAM J. Sci. Statist. Comput. (1982) 3: 76-110.
    [15] Requiem for second-order fluid approximations of traffic flow. Transportation Research Part B: Methodological (1995) 29: 277-286.
    [16] Modelling and simulation of vehicular traffic jam formation. Kinet. Relat. Models (2008) 1: 279-293.
    [17] All speed scheme for the low Mach number limit of the isentropic Euler equations. Commun. Comput. Phys. (2011) 10: 1-31.
    [18] Self-organized hydrodynamics with congestion and path formation in crowds. J. Comput. Phys. (2013) 237: 299-319.
    [19] Numerical simulations of the Euler system with congestion constraint. J. Comput. Phys. (2011) 230: 8057-8088.
    [20] Nonlinear follow-the-leader models of traffic flow. Operations Research (1961) 9: 545-567.
    [21] Solutions in the large for nonlinear hyperbolic systems of equations. Comm. Pure Applied Math. (1965) 18: 697-715.
    [22] Existence globale pour le systéme des gaz sans pression. Comptes Rendus Acad. Sci. (1995) 321: 171-174.
    [23] J. Jung, Schémas numériques adaptés aux accélérateurs multicoeurs pour les écoulements bifluides, PhD thesis, Univ. Strasbourg, 2014.
    [24] Stability for some equations of gas dynamics. Math. Comput. (1981) 37: 307-320.
    [25] On kinematic waves. Ⅱ. A theory of traffic flow on long crowded roads. Proc. Roy. Soc. London. Ser. A. (1955) 229: 317-345.
    [26] P. -L. Lions and N. Masmoudi, On a free boundary barotropic model, Ann. Inst. H. Poincaré Anal. Non Linéaire, 16 (1999), 373-410, doi: 10.1016/S0294-1449(99)80018-3
    [27] The deterministic version of the Glimm scheme. Comm. Math. Phys. (1977) 57: 135-148.
    [28] B. Maury and A. Preux, Pressureless Euler equations with maximal density constraint: A time-splitting scheme, Technical report, Université Paris-Sud, 2017, 333-355, Available on https://hal.archives-ouvertes.fr/hal-01224008. doi: 10.1515/9783110430417-014
    [29] The Prigogine-Herman kinetic model predicts widely scattered traffic flow data at high concentrations. Transportation Research Part B: Methodological (1998) 32: 589-604.
    [30] On Boltzmann-like treatments for traffic flow: A critical review of the basic model and an alternative proposal for dilute traffic analysis. Transportation Research (1975) 9: 225-235.
    [31] H. J. Payne, Freflo: A Macroscopic Simulation Model of Freeway Traffic, Transportation Research Record.
    [32] I. Prigogine and R. Herman, Kinetic Theory of Vehicular Traffic, American Elsevier Publishing, 1971.
    [33] Fundamental diagrams in traffic flow: The case of heterogeneous kinetic models. Communications in Mathematical Sciences (2016) 14: 643-669.
    [34] J. Smoller, Shock Waves and Reaction-Diffusion Equations, 2nd edition, Springer-Verlag, New York, 1994. doi: 10.1007/978-1-4612-0873-0
    [35] E. F. Toro, Riemann Solvers and Numerical Methods for Fluid Dynamics, 3rd edition, Springer-Verlag, Berlin, 2009. doi: 10.1007/b79761
    [36] A kinetic model for vehicular traffic derived from a stochastic microscopic model. Transport Theory and Stat. Phys. (1996) 25: 785-798.
    [37] A non-equilibrium traffic model devoid of gas-like behavior. Transportation Research Part B: Methodological (2002) 36: 275-290.
  • This article has been cited by:

    1. Haichao Sun, Hui Xu, Yanjun Li, Xikun Wang, Yalin Li, Parametric Analysis and Optimization Design of the Twin-Volute for a New Type of Dishwasher Pump, 2023, 11, 2227-9717, 305, 10.3390/pr11020305
  • Reader Comments
  • © 2017 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(4729) PDF downloads(82) Cited by(6)

Figures and Tables

Figures(13)  /  Tables(1)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog