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 @Override
22 protected double headwayFactor(final Speed desiredSpeed, final Length desiredHeadway, final Speed generationSpeed,
23 final Length generationHeadway, final Length leaderLength)
24 {
25 if (desiredSpeed.eq(generationSpeed))
26 {
27 return 1.0;
28 }
29 // following state at desired speed (capacity)
30 double k0 = 1.0 / (desiredHeadway.si + leaderLength.si);
31 double q0 = k0 * desiredSpeed.si;
32 // actual state
33 double k = 1.0 / (generationHeadway.si + leaderLength.si);
34 double q = k * generationSpeed.si;
35 // recovery flow qr
36 double qr = q; // if already on free flow branch
37 if (k > k0) // on congestion branch
38 {
39 // recovery wave speed assuming theta = 1 / q0 (theoretically sound and tested to be robust)
40 double rho = generationSpeed.si - q0 / k;
41 // recovery free flow state
42 double kr = (q - rho * k) / (desiredSpeed.si - rho);
43 qr = kr * desiredSpeed.si;
44 }
45 // efficiency factor, > 1 is larger headway, so less efficient
46 return q0 / qr;
47 }
48
49 }