1 package org.opentrafficsim.road.gtu.generator;
2
3 import org.djunits.value.vdouble.scalar.Length;
4 import org.djunits.value.vdouble.scalar.Speed;
5
6 /**
7 * Extends car-following placement with a first-order bounded acceleration (BA) principle. This principle comes down to
8 * reduction of efficiency by increasing headways, as the generated GTU speed is lower than the desired speed. The increased
9 * headways allow for acceleration to occur, allowing faster flow recovery over time.
10 * <p>
11 * Copyright (c) 2013-2024 Delft University of Technology, PO Box 5, 2600 AA, Delft, the Netherlands. All rights reserved. <br>
12 * BSD-style license. See <a href="https://opentrafficsim.org/docs/license.html">OpenTrafficSim License</a>.
13 * </p>
14 * @author <a href="https://github.com/averbraeck">Alexander Verbraeck</a>
15 * @author <a href="https://github.com/peter-knoppers">Peter Knoppers</a>
16 * @author <a href="https://github.com/wjschakel">Wouter Schakel</a>
17 */
18 public class CfBaRoomChecker extends CfRoomChecker
19 {
20
21 /**
22 * Constructor.
23 */
24 public CfBaRoomChecker()
25 {
26 //
27 }
28
29 @Override
30 protected double headwayFactor(final Speed desiredSpeed, final Length desiredHeadway, final Speed generationSpeed,
31 final Length generationHeadway, final Length leaderLength)
32 {
33 if (desiredSpeed.eq(generationSpeed))
34 {
35 return 1.0;
36 }
37 // following state at desired speed (capacity)
38 double k0 = 1.0 / (desiredHeadway.si + leaderLength.si);
39 double q0 = k0 * desiredSpeed.si;
40 // actual state
41 double k = 1.0 / (generationHeadway.si + leaderLength.si);
42 double q = k * generationSpeed.si;
43 // recovery flow qr
44 double qr = q; // if already on free flow branch
45 if (k > k0) // on congestion branch
46 {
47 // recovery wave speed assuming theta = 1 / q0 (theoretically sound and tested to be robust)
48 double rho = generationSpeed.si - q0 / k;
49 // recovery free flow state
50 double kr = (q - rho * k) / (desiredSpeed.si - rho);
51 qr = kr * desiredSpeed.si;
52 }
53 // efficiency factor, > 1 is larger headway, so less efficient
54 return q0 / qr;
55 }
56
57 }