
The development of China's manufacturing industry is still facing the challenge of regional imbalance. To solve the problem of development imbalance, it is necessary to realize regional development. First, we must analyze the development characteristics of different regions. To this end, we consider the requirements of the new development era and design an evaluation index system for the high-quality development level of the manufacturing industry from the dimensions of innovation, green, and efficiency. Then construct a novel hybrid model which combines the grey incidence clustering model and AP algorithm for panel data in this paper. According to the statistical data from 2014 to 2018, we find out the high-quality development of China's manufacturing industry is characterized by obvious regional differences, different development stages and different constraints.
Citation: Zhi-ying Han, Yong Liu, Xue-ge Guo, Jun-qian Xu. Regional differences of high-quality development level for manufacturing industry in China[J]. Mathematical Biosciences and Engineering, 2022, 19(5): 4368-4395. doi: 10.3934/mbe.2022202
[1] | Raimund Bürger, Christophe Chalons, Rafael Ordoñez, Luis Miguel Villada . A multiclass Lighthill-Whitham-Richards traffic model with a discontinuous velocity function. Networks and Heterogeneous Media, 2021, 16(2): 187-219. doi: 10.3934/nhm.2021004 |
[2] | Maya Briani, Emiliano Cristiani . An easy-to-use algorithm for simulating traffic flow on networks: Theoretical study. Networks and Heterogeneous Media, 2014, 9(3): 519-552. doi: 10.3934/nhm.2014.9.519 |
[3] | Mauro Garavello, Roberto Natalini, Benedetto Piccoli, Andrea Terracina . Conservation laws with discontinuous flux. Networks and Heterogeneous Media, 2007, 2(1): 159-179. doi: 10.3934/nhm.2007.2.159 |
[4] | Raimund Bürger, Kenneth H. Karlsen, John D. Towers . On some difference schemes and entropy conditions for a class of multi-species kinematic flow models with discontinuous flux. Networks and Heterogeneous Media, 2010, 5(3): 461-485. doi: 10.3934/nhm.2010.5.461 |
[5] | Helge Holden, Nils Henrik Risebro . Follow-the-Leader models can be viewed as a numerical approximation to the Lighthill-Whitham-Richards model for traffic flow. Networks and Heterogeneous Media, 2018, 13(3): 409-421. doi: 10.3934/nhm.2018018 |
[6] | Adriano Festa, Simone Göttlich, Marion Pfirsching . A model for a network of conveyor belts with discontinuous speed and capacity. Networks and Heterogeneous Media, 2019, 14(2): 389-410. doi: 10.3934/nhm.2019016 |
[7] | Christophe Chalons, Paola Goatin, Nicolas Seguin . General constrained conservation laws. Application to pedestrian flow modeling. Networks and Heterogeneous Media, 2013, 8(2): 433-463. doi: 10.3934/nhm.2013.8.433 |
[8] | Raimund Bürger, Stefan Diehl, M. Carmen Martí, Yolanda Vásquez . A difference scheme for a triangular system of conservation laws with discontinuous flux modeling three-phase flows. Networks and Heterogeneous Media, 2023, 18(1): 140-190. doi: 10.3934/nhm.2023006 |
[9] | Giuseppe Maria Coclite, Lorenzo di Ruvo, Jan Ernest, Siddhartha Mishra . Convergence of vanishing capillarity approximations for scalar conservation laws with discontinuous fluxes. Networks and Heterogeneous Media, 2013, 8(4): 969-984. doi: 10.3934/nhm.2013.8.969 |
[10] | Wen Shen . Traveling wave profiles for a Follow-the-Leader model for traffic flow with rough road condition. Networks and Heterogeneous Media, 2018, 13(3): 449-478. doi: 10.3934/nhm.2018020 |
The development of China's manufacturing industry is still facing the challenge of regional imbalance. To solve the problem of development imbalance, it is necessary to realize regional development. First, we must analyze the development characteristics of different regions. To this end, we consider the requirements of the new development era and design an evaluation index system for the high-quality development level of the manufacturing industry from the dimensions of innovation, green, and efficiency. Then construct a novel hybrid model which combines the grey incidence clustering model and AP algorithm for panel data in this paper. According to the statistical data from 2014 to 2018, we find out the high-quality development of China's manufacturing industry is characterized by obvious regional differences, different development stages and different constraints.
Traffic flow models based on scalar conservation laws with continuous flux functions are widely used in the literature. For a general presentation of the models, we refer to the books [11,12,23] and the references therein. Extensions to road traffic networks have been also established. We mention in particular the contributions [6,15], where the authors introduce the coupled network problem and show the existence of solutions. Within this article, we are concerned with the special case of scalar conservation laws with discontinuous flux in the unknown that are motivated in the traffic flow theory by the observation of a gap between the free flow and the congested flow regime [4,5,8]. This phenomenon generates an interesting dynamical behavior called zero waves, i.e., waves with infinite (negative) speed but zero wave strength, and has been investigated in recent years either from a theoretical or numerical point of view, see for instance [2,19,20,22,24] or more generally [1,3,7,13].
To the best of our knowledge, the study of scalar conservation laws with discontinuous flux functions on networks is still missing in the traffic flow literature. However, in the context of supply chains with discontinuous flux such considerations have been already done [10,14]. We remark that supply chain models differ essentially from traffic flow models due to simpler dynamics and different coupling conditions.
In this work, we aim to derive a traffic network model, where the dynamics on each road are governed by a scalar conservation law with discontinuous flux function in the unknown. For simplicity, we restrict to piecewise linear flux functions. Special emphasis is put on the coupling at junction points to ensure a unique admissible weak solution. In particular, we focus on dispersing junctions where the number of incoming roads does not exceed the number of outgoing roads and merging with two incoming and one outgoing road. The latter type of junction can be extended to the case of multiple incoming roads and a single outgoing one. In order to construct a suitable numerical scheme that is not based on regularization techniques we adapt the splitting algorithm originally introduced in [22]. Therein, the discontinuous flux is decomposed into a Lipschitz continuous flux and a Heaviside flux such that a two-point monotone flux scheme, e.g., Godunov, can be employed in an appropriate manner. This algorithm has been studied in [22] for the case of a single road only. However, in the network case, multiple roads with possibly disjunctive flux functions need to be considered at a junction point to ensure mass conservation. Hence, the key challenge is to determine the correct flux through the junction in an appropriate manner. Therefore, a detailed case distinction in accordance with the theoretical investigations is provided for the different types of junctions. The numerical results validate the proposed algorithm for some relevant network problems.
The paper is organized as follows: in Section 2 we discuss the basic model and Riemann problems which permit to derive an exact solution. We extend the modeling framework to networks in Section 3 and focus on the coupling conditions. In Section 4, we introduce how the splitting algorithm [22] can be extended to also deal with the different types of junctions. Finally, we present a suitable discretization and numerical simulations in Section 5.
In this section, we briefly recall the case of the Lighthill-Whitham-Richards (LWR) model [18,21] on a single edge with a flux function having a single decreasing jump at
Following [22], we consider the scalar conservation law Eq (2.1),
{ut+f(u)x=0,(x,t)∈(a,b)×(0,T)=:ΠT,u(x,0)=u0(x)∈[0,umax],x∈(a,b),u(a,t)=r(t)∈[0,umax],t∈(0,T),u(b,t)=s(t)∈[0,umax],F(t)∈˜f(s(t)),t∈(0,T). | (2.1) |
More precisely the flux function is defined as follows Eq (2.2),
f(u)={f1(u) ˆ=Flux1,ifu∈[0,u∗],f2(u) ˆ=Flux2,ifu∈(u∗,umax]. | (2.2) |
We denote
α:=f(u∗−)−f(u∗+). | (2.3) |
As usual we require for the flux function
As in [22], the multivalued version of
˜f(u)={f(u),u∈[0,u∗),[f(u∗+),f(u∗−)],u=u∗,f(u),u∈(u∗,umax]. | (2.4) |
Finally, we have to discuss the imposed boundary conditions at
F(t)={f(u∗−),ifthetrafficaheadofx=bisfree−flowing,f(u∗+),ifthetrafficaheadofx=biscongested. | (2.5) |
The state of traffic ahead of
Remark 2.1. [22,Remark 1.3] We note that for the boundary condition at the left end the state of traffic ahead of
The following assumptions are important for the proof of existence and uniqueness of solutions.
Assumption 2.2. [22,Assumption 1.1] The initial data satisfies
A weak solution is intended in the following sense:
Definition 2.3 [22,Definition 1.1] A function
v(x,t)∈˜f(u(x,t))a.e. |
such that for each
∫T0∫ba(uψt+vψx)dxdt+∫bau0(x)ψ(x,0)dx=0. |
As usual, weak solutions do not lead to a unique solution and additional criteria are necessary to rule out physically incorrect solutions. In particular, the discontinuity of the flux prohibits from directly using the classical approaches. Note that in [22] an adapted version of Oleinik's entropy condition [9] is used to single out the correct solution, while in [24] the convex hull construction [17] is used to construct solutions to Riemann problems.
Here, we will concentrate on the convex hull construction. For completeness we will shortly recall the solutions to Riemann problems considered in [24] as they are essential in order to construct a Riemann solver at a junction.
We consider a Riemann problem with initial data
We consider the following flux function by Eq (2, 6),
f(u)={d1u+d0,x≤u∗(Flux1),e1u+e0,x>u∗(Flux2) | (2.6) |
with the regularized flux function given by Eq (2.7),
fϵ(u)={f1(u)=d1u+d0,0≤u≤u∗,fϵmid(u)=−1ϵ(f1(u∗)−fϵ2(u∗+ϵ))(u−u∗)+f1(u∗),u∗<u<u∗+ϵ,fϵ2(u)=e1u+e0,u∗+ϵ<u≤umax. | (2.7) |
We define
Case 1: Either
This case corresponds to the classical case of solving Riemann problems, where the solution consists of a single rarefaction wave or shock, see [17].
Case 2:
By using the smallest convex hull approach the solution consists of a contact line following
s=f(uL)−f(u∗)uL−u∗<0, |
where we recall that
u(x,t)={uL,ifx<st,u∗,ifst≤x≤d1t,uR,ifx>d1t. | (2.8) |
Case 3:
Here, the solution is given by a shock connecting
s=f(u∗+)−f(uL)u∗−uL<0. |
Note that due to
u(x,t)={uL,ifx<st,u∗,ifst≤x≤e1t,uR,ifx>e1t. | (2.9) |
Case 4:
In this case, we get only one shock connecting
u(x,t)={uL,ifx<st,uR,ifx≥st. | (2.10) |
Remark 2.4. As aforementioned in [19] Riemann solutions for piecewise quadratic discontinuous flux functions are derived. They also cover the case of a piecewise linear flux function if the quadratic terms are zero. For a general quadratic discontinuous flux, the solutions are more involved since no contact discontinuities occur.
Up to now, we have not addressed the case, where one of the boundary conditions equals the critical density
Next, we focus on networks where we allow for discontinuous flux functions. The key idea is to consider the regularized flux function
We start with a short introduction to the network setting. For more details on traffic flow network models we refer the reader to [11] and the references therein.
Let
We call the couple
In order to derive the network solution, we restrict to the description of a single junction
A=(β1,1⋯β1,m⋮⋮⋮βn,1⋯βn,m). |
To conserve the mass we assume
Definition 3.1 (Weak solution at a junction). Let
n∑i=1(∫T0∫Iini(uini(x,t)∂tϕini(x,t)+wini(x,t)∂xϕini(x,t))dxdt)+m∑j=1(∫T0∫Ioutj(uoutj(x,t)∂tϕoutj(x,t)+woutj(x,t)∂xϕoutj(x,t))dxdt)=0, |
for every collection of test functions
ϕini(bini,⋅)=ϕoutj(aoutj,⋅),∂xϕini(bini,⋅)=∂xϕoutj(aoutj,⋅), |
for
Additionally, in order to get unique solutions, we will consider the following concept of admissible solutions, which adapts ref. [11] (rule (A) and (B), p. 81) to the discontinuous setting:
Definition 3.2 (Admissbile Weak Solution). We call
1.
2.
3.
4.
In particular, the maximization of the inflow with respect to the distributions parameters and the technical assumption [11] guarantee the uniqueness of solutions for a continuous flux.
If
For solving the maximization problems imposed by the definition 3.2 so-called supply and demand functions can be used, see [16]. The demand describes the maximal flux the incoming road wants to send. In contrast, the supply describes the maximal flux the outgoing road is able to absorb. The definition of the supply and demand functions of the regularized function is straightforward. As
Definition 3.3. For a network with flux function
D(u)={f(u),u∈[0,u∗),f(u∗−),u∈[u∗,umax]. | (3.1) |
On the contrary, the supply reads as
S(u)={f(u∗−),u∈[0,u∗),f(u),u∈(u∗,umax] | (3.2) |
and
S(u∗)={f(u∗−),freeflowing,f(u∗+),congested. | (3.3) |
Remark 3.4. If we consider the regularized flux function
In order to show existence and uniqueness in the discontinuous case we need to define an additional function. For a regularized flux function we notice that for every flux value, we get two different density values, see left picture in Figure 4. As different density values lead to different solutions, we need to be able to distinguish them and choose the correct solution. In the continuous or regularized case a mapping usually called
Definition 3.5. Let the function
1.
2. For
Note that if
Remark 3.6. We note that the mapping
Now, we present a Riemann solver for two types of junctions. First, we consider a junction with
Theorem 3.7. Let
uini∈{{uini,0}∪(η(uini,0),umax],ifuini,0∈[0,f−1+],{uini,0}∪[u∗,umax],ifuini,0∈(f−1+,u∗),[u∗,umax],ifuini,0∈[u∗,umax], | (3.4) |
and
uoutj∈{[0,u∗],ifuoutj,0∈[0,u∗),[0,u∗],ifuoutj,0=u∗andfreeflowing,{u∗}∪[0,f−1+),ifuoutj,0=u∗andcongested,{uoutj,0}∪[0,η(uoutj,0)),ifuoutj,0∈(u∗,umax]. | (3.5) |
Proof. Using the definition of the supply and demand functions in definition 3.3 and the results from [12,section 5.2.3] we can follow the proof of [12,theorem 5.1.2] and uniquely determine the inflows which maximize the flux through junctions subject to the distribution parameters. It remains to show that by the choice of the density values the correct waves are induced.
We start with considering the outgoing roads. If
For
On the contrary, considering the incoming roads and
Now, let
s=fini−f(uini,0)u∗−uini,0<0. |
Further, if
Now, let us turn to the remaining case of
s=fini−f(uini,0)u∗−uini,0≤0, |
as
Hence, the choices of the densities induce the correct waves.
Now, we consider the case of more incoming than outgoing roads. Exemplary, we study the 2–to–1 situation, even though the results can be easily extended to the
Theorem 3.8. Let
uini∈{{uini,0}∪(η(uini,0),umax],ifuini,0∈[0,f−1+],{uini,0}∪[u∗,umax],ifuini,0∈(f−1+,u∗),[u∗,umax],ifuini,0∈[u∗,umax], | (3.6) |
and
uout1∈{[0,u∗],ifuout1,0∈[0,u∗),[0,u∗],ifuout1,0=u∗andfreeflowing,{u∗}∪[0,f−1+),ifuout1,0=u∗andcongested,{uout1,0}∪[0,η(uout1,0)),ifuout1,0∈(u∗,umax]. | (3.7) |
Proof. Following [11,Section 3.2.2], the flux values at the junction can be calculated with the following steps:
1. Calculate the maximal possible flux
2. Consider the right of way parameter and the flux maximization and calculate the intersection
3. If
4. The flux values are given by
Completely analogous to theorem 3.7 we can show that the choice of the densities admits the correct wave speeds.
The Riemann solutions proposed in theorem 3.7 and theorem 3.8 are the key ingredients for the splitting algorithm on networks in the next section.
Different problems might occur when designing a numerical scheme for a conservation law with discontinuous flux. However, the main difficulties are induced by the zero waves. Since these waves have infinite speed, the regular CFL condition is scaled by the regularization parameter
We consider a flux function
g(u)=−αH(u−u∗), |
where
This case has been already treated in [22] and will be the basis for the splitting algorithm on networks. When solving the scalar conservation law (2.1) on a single road, the boundary value in the case
˜g(u)={0,u∈[0,u∗),[−α,0],u=u∗,−α,u∈(u∗,umax]. | (4.1) |
Furthermore, we define
F(t)=f(u∗−)⇔G(t)=0,F(t)=f(u∗+)⇔G(t)=−α. | (4.2) |
Additionally to the assumptions 2.2, we assume:
Assumption 4.1. [22,Assumption 1.1] The initial data satisfies
Remark 4.2. We emphasize that the original splitting algorithm for a single road [22] is not limited to piecewise linear discontinuous flux functions. Another prominent example might be concave piecewise quadratic flux functions with discontinuity again at
Then, we are able to handle the flux function
λ=ΔtΔx. |
For an integer
As the algorithm splits the function
We denote the backward spatial difference by
rn=rn+12=r(tn),Un+120=Un0=rnsn=sn+12=s(tn)Un+12K+1=UnK+1=sn. |
The function
gn+12K+1=gnK+1=G(tn)={0,ifs(tn)<u∗,0,ifs(tn)=u∗,trafficaheadofx=bisfree−flowing,−α,ifs(tn)=u∗,trafficaheadofx=biscongested,−α,ifs(tn)>u∗. | (4.3) |
That means, we can describe the boundary value
Definition 4.3 [22,Eq (3.7)] Let
˜G(u)={u,u∈[0,u∗),[u∗,u∗+λα],u=u∗,u+λα,u∈(u∗,umax], ˜G−1(u)={u,u∈[0,u∗),u∗,u∈[u∗,u∗+λα),u−λα,u∈[u∗+λα,umax+λα]. |
The splitting algorithm [22] can be then expressed as
{{Un+1/2k=˜G−1(Unk−λgn+1/2k+1),k=K,K−1,…,1,gn+1/2k=(Un+1/2k−Unk+λgn+1/2k+1)/λ,k=K,K−1,…,1,Un+1k=Un+1/2k−λΔ−pg(Un+1/2k+1,Un+1/2k), k∈K. | (4.4) |
Note that the first half step, which includes the first two equations, is implicit. Nevertheless, instead of solving a nonlinear system of equations, the equation can be solved backwards in space starting with
We note that for the implicit equation a CFL condition is not needed, but it is required for the third step. As
As shown in [22,Theorem 5.1] the splitting algorithm (4.4) converges to a weak solution of Eq (2.1). However, obtaining a similar statement about weak entropy solutions is still an open problem.
The key idea to numerically solve such discontinuous conservation laws on networks is to use the splitting algorithm only on the roads and determine the correct in- and outflows at the boundaries by the help of the Riemann solver established in, e.g., theorem 3.7. As the splitting algorithm works with flux values, there is no need to compute the exact densities at the junction. Instead we need to know how the solution at the junction influences the flux values. The algorithm for a single junction is depicted in algorithm 1. The general description of the algorithm allows for either junctions with given distribution or right of way parameters. For simplicity, we assume that each road is represented by the same interval
Remark 4.4. Note that this simplification enables the use of the same grid points on each road which spares further sub- or superindices. However, the algorithm can be easily adapted to different road lengths.
We assume in the following that the space and time grid is the same as in the previous subsection. The approximate solutions are denoted by
The overall strategy of the splitting algorithm on a network consists of three important steps:
1. Solve the optimization problem induced by definition 3.2 (in particular item 4) at the junction to calculate the flux values
Here, it is crucial to use the discontinuous flux function
These flux values bring us now to:
Require: number of incoming roads Ensure: approximate solutions 1: Initilization: 2: 3: 4: 5: for 6: Solve the by definition 3.2 induced optimization problem at the junction based on the flux 7: Compute the densities at the junction with an appropriate Riemann solver 8: Compute the adjusted flux values for the incoming roads 9: for 10: 11: for 12: 13: 14: end for 15: 16: for 17: 18: end for 19: end for 20: for 21: Compute 22: for 23: 24: 25: end for 26: 27: for 28: 29: end for 30: end for 31: end for |
Using the calculated (unadjusted) flux values from step one, we can determine the densities at the junction with the help of the appropriate Riemann solver (theorem 3.7 and Eqs (3.4)–(3.5) or theorem 3.8 and Eqs (3.6)–(3.7)) at the junction (line 7). Then, these density values can be used to calculate the corresponding flux value of
In addition, the first steps give us all the ingredients for the final step:
The adjusted flux values from the previous step are important for the second half step (line 17 or 28) of the splitting algorithm which uses a Godunov type scheme based on
Further boundary data is needed in the first half step of the algorithm, lines 12–13 and 23–24. Here, we start with
gn+1/2,ini,K+1=gn,ini,K+1=fini−fini,adj. | (4.5) |
Note that the definition of
Furthermore, we can decrease the computational costs of the algorithm: In the second step (or in line 7 of algorithm 1) the density at the junction is computed. This can be very expensive and hence we aim to avoid this. In the third step we have seen that for the missing boundary data only the adjusted flux values are necessary and not the densities at the junction themselves. Therefore, studying first each junction type in detail allows to determine the corresponding flux values based on the density values and the supply and demand functions and the intermediate expensive step for the computation of the exact densities can be skipped. Hence, we combine the first and second step of the strategy in one single step. In the following, we will study a 1–to–1 junction in detail and present the tailored algorithm. As the strategy is completely analogous for the 1–to–2 junction and 2–to–1 junction, we will only present the algorithms and discuss important properties. The algorithms can then be used to replace the lines 6 to 8 in algorithm 1.
Remark 4.5. The extension to
Remark 4.6. Further note that the presented strategy and also algorithm 1 can be used for arbitrary junctions and nonlinear discontinuous flux functions once an appropriate Riemann solver similar to theorem 3.7 and 3.8 is established.
First, we consider a junction with one incoming and one outgoing road in detail. Let
Case A: demand and supply are equal
There are two different situations depending on the density value of the incoming road where demand and supply can be equal.
1.
2.
Case B: supply is restrictive
If the supply is restrictive, i.e.
1.
2.
Case C: demand is restrictive
If the demand is restrictive, the outgoing road is able to take the whole amount of vehicles the demand wants to send. Here, we only have one possible situation for the initial condition on the incoming road:
1.
Note that we have the same flux function on each road. Hence, in the case
The whole procedure is summarized in algorithm 2.
Require: Demand Ensure: Flux values 1: 2: 3: if 4: 5: else if 6: 7: 8: end if |
Remark 4.7. Theoretically, the Riemann solver in theorem 3.7 coincides for a 1–to–1 junction with the Riemann solver on a single road. Hence, the procedure described in algorithm 2 leads to the same solution. In contrast to that on the numerical level, the splitting algorithm for a 1–to–1 junction does not exactly coincide with the splitting algorithm used for a single road [22]. The reason for the computational difference is that the splitting algorithm for the 1–to–1 junction considers the exact solution of the Riemann problem at the junction point and hence uses exact values for
The difference to the 1–to–1 junction is now that we have to consider two supplies. So the case distinctions to determine the flux values are a bit more complex. However, the procedure itself does not change. The details can be seen in algorithm 3.
We remark that if the inflow on the first road equals the demand and the restriction given by at least one of the supplies, the latter road needs to be congested. The Riemann solver states congestion such that the flux needs to be adjusted. Then, the incoming road either stays free flowing or the solution is given by
Require: Demand Ensure: Flux values 1: 2: 3: if 4: if 5: 6: else if 7: 8: end if 9: else if 10: if 11: 12: else 13: 14: end if 15: if 16: 17: 18: else if 19: 20: else if 21: 22: end if 23: end if |
If at least one of the supply restrictions is active, we need to adjust the corresponding flux values as in the 1–to–1 case and the incoming road. Nevertheless, here an interesting case (which is not possible in the 1–to–1 situation) can occur. The solution of the Riemann problem at the incoming road can be given by
Recall that for two incoming and one outgoing roads, the maximal possible flux on the outgoing road is given by
Require: Demand Ensure: Flux values 1: 2: 3: if 4: 6: 7: end if 8: 9: if 10: 11:else if 12: if 13: 14: end if 15: if 16: if 17: 18: else 19: 20: end if 21: end if 22: if 23: if 24: 25: else 26: 27: end if 28: end if 29: end if |
As before, no adjustment is needed when the demand on both roads is smaller than the supply. If the supply is active we might need to adjust the outgoing road and in most cases at least one incoming road. As in the 1–to–2 case,
In this section, we present some numerical examples to compare the splitting algorithm on networks with the Riemann solver. Further, we compute the solution using a regularized flux and the Godunov scheme. We consider the following flux function:
f(u)={u,u∈[0,0.5),0.5(1−u),u∈[0.5,1]. | (5.1) |
The corresponding regularization is given by Eq (2.7).
We consider in particular the 1–to–2 and 2–to–1 situations. For our comparison, we choose constant initial data on each road. The junction is always located at
In both scenarios the supply of the first outgoing road is restrictive. The parameter settings are as follows:
The exact solution is given by
uin1(x,t)={0.4,x≥−32t,0.5,−32t<x<−12t,1315,−12t<x<0,uout1(x,t)=0.9,uout2(x,t)={115β20≤x≤841t,0.7,841t<x. |
Apparently, the solution induces two waves on the incoming road.
The exact solution is given by
uin1(x,t)={0.4,x≥−t,0.5,−t<x<0,uout1(x,t)=0.7,1≤x≤3,uout2(x,t)={0.3β21≤x≤t,0.2,t<x. |
This example generates
In Table 1, the
Example 1 | |||||
Splitt. | Reg. | ||||
33.44e-03 | 46.77e-03 | 82.26e-03 | 51.82e-03 | ||
24.17e-03 | 29.05e-03 | 65.59e-03 | 30.69e-03 | ||
14.16e-03 | 20.12e-03 | 60.86e-03 | 24.45e-03 | ||
8.97e-03 | 12.49e-03 | 58.37e-03 | 20.44e-03 | ||
CR | 0.64695 | 0.62453 | 0.1593 | 0.4353 | |
Example 2 | |||||
Splitt. | Reg. | ||||
4.58e-03 | 7.41e-03 | 44.70e-03 | 14.57e-03 | ||
2.97e-03 | 4.24e-03 | 43.47e-03 | 11.28e-03 | ||
2.03e-03 | 2.89e-03 | 42.50e-03 | 10.06e-03 | ||
1.24e-03 | 1.99e-03 | 41.80e-03 | 9.21e-03 | ||
CR | 0.61911 | 0.62327 | 0.0322 | 0.2150 |
We can see that the error terms obtained by the splitting algorithm are the lowest and so the computational costs. Obviously, the error terms increase with a lower CFL due to numerical diffusion. For a direct comparison with the regularized approach, the CFL condition should be the same. Meaning that the second column in Table 1 for the splitting algorithm should be compared with the first one of the regularized approach. In this case, we can see that the splitting algorithm also performs better in both examples. By choosing a smaller regularization parameter the performance of the regularized approach increases, but also the computational costs. To obtain similar error terms as for the splitting algorithm the regularization parameter needs to be further reduced at very high computational costs.
For
Here, we consider two scenarios, where in the first scenario the demand is restrictive while in the second one the supply. The parameter settings are as follows:
The exact solution is given by
uin1(x,t)=0.2,uout1(x,t)={0.450≤x≤t,0.3,t<xuin2(x,t)=0.25. |
The exact solution is given by
uin1(x,t)={0.5x≤−2t,0.6,−2t<x<0,uout1(x,t)={0.50≤x≤t,0.4,t<x,uin2(x,t)={0.8x≤−0.5t,0.7,−0.5t<x<0. |
In particular, the flux value for the first incoming road needs to be adjusted from
Considering the
Example 1 | |||||
Splitt. | Reg. | ||||
9.25e-03 | 16.22e-03 | 16.22e-03 | 17.01e-03 | ||
5.90e-03 | 11.63e-03 | 11.63e-03 | 12.19e-03 | ||
2.98e-03 | 8.13e-03 | 8.13e-03 | 8.52e-03 | ||
8.97e-03 | 5.71e-03 | 5.71e-03 | 5.99e-03 | ||
CR | 0.53838 | 0.50353 | 0.50353 | 0.50347 | |
Example 2 | |||||
Splitt. | Reg. | ||||
14.12e-03 | 20.10e-03 | 85.69e-03 | 27.55e-03 | ||
9.65e-03 | 13.86e-03 | 79.98e-03 | 21.65e-03 | ||
6.41e-03 | 9.57e-03 | 75.96e-03 | 17.49e-03 | ||
4.51e-03 | 6.69e-03 | 73.22e-03 | 14.65e-03 | ||
CR | 0.55295 | 0.52959 | 0.07551 | 0.30432 |
For
We have presented a Riemann solver at a junction for conservation laws with discontinuous flux. We have adapted the splitting algorithm of [22] to networks and demonstrated its validity in comparison with the exact solution. We have also pointed out that the splitting algorithm on networks is faster and more accurate than the approach with a regularized flux. Future research includes the investigation of other network models, where the flux is discontinuous.
J.F. was supported by the German Research Foundation (DFG) under grant HE 5386/18-1 and S.G. under grant GO 1920/10-1.
The authors declare there is no conflict of interest.
[1] |
F. Yang, Y. Sun, Y. Zhang, T. Wang, Factors affecting the manufacturing industry transformation and upgrading: a case study of Guangdong-Hong Kong-Macao greater bay area, Int. J. Environ. Res. Public Health, 18 (2021), 7157. https://doi.org/10.3390/ijerph18137157 doi: 10.3390/ijerph18137157
![]() |
[2] |
Y. Fu, A. Supriyadi, T. Wang, L. Wang, G. T. Cirella, Effects of regional innovation capability on the green technology efficiency of China's manufacturing industry: evidence from listed companies, Energies, 13 (2020), 5467. https://doi.org/10.3390/en13205467 doi: 10.3390/en13205467
![]() |
[3] |
T. Li, L. Liang, D. Han, Research on the efficiency of green technology innovation in China's provincial high-end manufacturing industry based on the RAGA-PP-SFA model, Math. Probl. Eng., 2018 (2018). https://doi.org/10.1155/2018/9463707 doi: 10.1155/2018/9463707
![]() |
[4] | Y. Guo, L. G. Zhang, An empirical study on the impact of environmental regulation on R&D input of industrial enterprises, China Popul. Resour. Environ., S3 (2014), 104–107. |
[5] |
M. Perkmann, A. Neely, K. Walsh, How should firm evaluate success in university-industry alliances? A performance measurement system, R&D Manage., 41 (2011), 202–216. https://doi.org/10.1111/j.1467-9310.2011.00637.x doi: 10.1111/j.1467-9310.2011.00637.x
![]() |
[6] |
K. Frenken, G. J. Heimeriks, J. Hoekman, What drives university research performance? An analysis using the CWTS Leiden ranking data, J. Informetrics, 11 (2017), 859–872. https://doi.org/10.1016/j.joi.2017.06.006 doi: 10.1016/j.joi.2017.06.006
![]() |
[7] |
B. V. Looy, K. Debackere, P. Andries, Policies to stimulate regional innovation capabilities via university-industry collaboration analysis and an assessment, R&D Manage., 33 (2003), 209–229. https://doi.org/10.1111/1467-9310.00293 doi: 10.1111/1467-9310.00293
![]() |
[8] |
V. Dewangan, M. Godse, Towards a holistic enterprise innovation performance measurement system, Technovation, 34 (2014), 536–545. https://doi.org/10.1016/j.technovation.2014.04.002 doi: 10.1016/j.technovation.2014.04.002
![]() |
[9] |
J. Youtie, P. Shapira, Building an innovation hub: a case study of the transformation of university roles in regional technological and economic development, Res. Policy, 37 (2008), 1188–1204. https://doi.org/10.1016/j.respol.2008.04.012 doi: 10.1016/j.respol.2008.04.012
![]() |
[10] | H. N. Tian, Q. Q. Sun, Research on the evaluation of green technology innovation capability of automobile manufacturing enterprises based on cloud model, Manage. Rev., 32 (2020), 102–114. |
[11] |
Y. Ma, G. Hou, Q. Yin, B. Xin, Y. Pan, The sources of green management innovation: Does internal efficiency demand pull or external knowledge supply push?, J. Cleaner Prod., 202 (2018), 582–590. https://doi.org/10.1016/j.jclepro.2018.08.173 doi: 10.1016/j.jclepro.2018.08.173
![]() |
[12] |
P. Li, Y. Ouyang, Technical change and green productivity, Environ. Resour. Econ., 76 (2020), 271–298. https://doi.org/10.1007/s10640-020-00424-1 doi: 10.1007/s10640-020-00424-1
![]() |
[13] | J. X. Li, Research on the evaluation of the degree of intelligentization of China's manufacturing industry and its influencing factors, China Soft Sci., 1 (2020), 154–163. |
[14] |
Y. Deng, D. You, J. Wang, Optimal strategy for enterprises' green technology innovation from the perspective of political competition, J. Cleaner Prod., 235 (2019), 930–942. https://doi.org/10.1016/j.jclepro.2019.06.248 doi: 10.1016/j.jclepro.2019.06.248
![]() |
[15] |
B. Lu, Expedited patent examination for green inventions: developing countries' policy choices, Energy Policy, 61 (2013), 1529–1538. https://doi.org/10.1016/j.enpol.2013.06.028 doi: 10.1016/j.enpol.2013.06.028
![]() |
[16] |
X. Pan, B. Ai, C. Li, X. Pan, Y. Yan, Dynamic relationship among environmental regulation, technological innovation and energy efficiency based on large scale provincial panel data in China, Technol. Forecast. Soc. Change, 144 (2019), 428–435. https://doi.org/10.1016/j.techfore.2017.12.012 doi: 10.1016/j.techfore.2017.12.012
![]() |
[17] | X. Fan, D. P. Zhao, Y. He, Enterprise innovation efficiencies of university-industry cooperation and their influential factors, Sci. Res. Manage., 2 (2012), 33–39. |
[18] | Y. Y. Ma, X. L. Zhang, Y. T. Sun, Environmental regulations motivate enterprises to make R&D efforts? Empirical evidence from thermal power enterprises, Sci. Res. Manage., 39 (2018), 66–74. |
[19] | J. Y. Liu, J. Xie, Environmental regulation and manufacturing export quality upgrade—based on the perspective of factor input structure heterogeneity, China Popul. Resour. Environ., 28 (2018), 158–167. |
[20] |
S. Roper, P. Vahter, J. H. Love, Externalities of openness in innovation, Res. Policy, 42 (2013), 1544–1554. https://doi.org/10.1016/j.respol.2013.05.006 doi: 10.1016/j.respol.2013.05.006
![]() |
[21] | L. Qian, R. Q. Xiao, Z. W. Chen, Research on green technology innovation efficiency and regional differences of industrial enterprises in my country—based on common frontier theory and DEA model, Econ. Theory Econo. Manage., 1 (2015), 26–43. |
[22] |
D. Li, Y. Zhao, L. Zhang, X. Chen, C. Cao, Impact of quality management on green innovation, J. Cleaner Prod., 170 (2018), 462–470. https://doi.org/10.1016/j.jclepro.2017.09.158 doi: 10.1016/j.jclepro.2017.09.158
![]() |
[23] |
G. F. Zhu, D. Wang, Analysis and prospects of the current development of China's manufacturing industry: based on the evaluation index system of manufacturing powers, J. Manage. Eng., 31 (2017), 1–7. https://doi.org/10.13587/j.cnki.jieem.2017.04.001 doi: 10.13587/j.cnki.jieem.2017.04.001
![]() |
[24] | S. R. Liang, L. W. Luo, The dynamic effect of technology spillovers of international R&D capital on the efficiency of green technology innovation, Sci. Res. Manage., 40 (2019), 21–29. |
[25] | Y. M. Wu, A spatial panel econometrics analysis on industrial R&D, university-industry cooperation and innovation performance, Sci. Res. Manage., 36 (2015), 118–127. |
[26] |
J. Zhang, Y. Chang, L. Zhang, D. Li, Do technological innovations promote urban green development?—A spatial econometric analysis of 105 cities in China, J. Cleaner Prod., 182 (2018), 395–403. https://doi.org/10.1016/j.jclepro.2018.02.067 doi: 10.1016/j.jclepro.2018.02.067
![]() |
[27] |
J. Schleich, R. Walz, M. Ragwitz, Effects of policies on patenting in wind-power technologies, Energy Policy, 108 (2017), 684–695. https://doi.org/10.1016/j.enpol.2017.06.043 doi: 10.1016/j.enpol.2017.06.043
![]() |
[28] | W. H. Li, The spatiotemporal evolution and influencing factors of China's provincial industrial green technological innovation output: an empirical study based on 30 provincial data, Chin. J. Manage. Eng., 31 (2017), 9–19. |
[29] |
R. Leoncini, A. Marzucchi, S. Montresor, F. Rentocchini, U. Rizzo, 'Better late than never': the interplay between green technology and age for firm growth, Small Bus. Econo., 52 (2019), 891–904. https://doi.org/10.1007/s11187-017-9939-6 doi: 10.1007/s11187-017-9939-6
![]() |
[30] |
S. Lin, J. Sun, D. Marinova, D. Zhao, Evaluation of the green technology innovation efficiency of China's manufacturing industries: DEA window analysis with ideal window width, Technol. Anal. Strategic Manage., 30 (2018), 1166–1181. https://doi.org/10.1080/09537325.2018.1457784 doi: 10.1080/09537325.2018.1457784
![]() |
[31] | H. X. Tang, X. Z. Zhang, Y. Wu, Z. C. He, Research on the development quality of China's manufacturing industry and the improvement of international competitiveness, China Soft Sci., 2 (2019), 128–142. |
[32] |
H. Z. Zheng, Research on enterprise green technology innovation and diffusion dynamics under environmental regulation, Sci. Manage. Res., 34 (2016), 77–88. https://doi.org/10.19445/j.cnki.15-1103/g3.2016.05.020 doi: 10.19445/j.cnki.15-1103/g3.2016.05.020
![]() |
[33] |
J. Chen, J. Cheng, S. Dai, Regional eco-innovation in China: an analysis of eco-innovation levels and influencing factors, J. Cleaner Prod., 153 (2017), 1–14. https://doi.org/10.1016/j.jclepro.2017.03.141 doi: 10.1016/j.jclepro.2017.03.141
![]() |
[34] | M. Wei, S. H. Li, Research on the measurement of China's economic high-quality development level in the new era, J. Quant. Econo. Tech. Econo., 35 (2018), 3–20. |
[35] | L. Sun, C. Miao, L. Yang, Ecological-economic efficiency evaluation of green technology innovation in strategic emerging industries based on entropy weighted TOPSIS method, Ecol. Indic., 73 (2017), 554–558. |
[36] |
Q. Luo, C. Miao, L. Sun, X. Meng, M. Duan, Efficiency evaluation of green technology innovation of China's strategic emerging industries: an empirical analysis based on Malmquist-data envelopment analysis index, J. Cleaner Prod., 238 (2019), 117782. https://doi.org/10.1016/j.jclepro.2019.117782 doi: 10.1016/j.jclepro.2019.117782
![]() |
[37] |
Y. Yu, G. Xin, Y. Chen, Deconstruction of the "black box" in industry-university-research institute's technology transfer and evaluation of its efficiency, Sci. Res. Manage., 4 (2017), 28–37. https://doi.org/10.19571/j.cnki.1000-2995.2017.04.004 doi: 10.19571/j.cnki.1000-2995.2017.04.004
![]() |
[38] |
S. Jain, K. P. Triantis, S. Liu, Manufacturing performance measurement and target setting: A data envelopment analysis approach, Eur. J. Oper. Res., 214 (2011), 616–626. https://doi.org/10.1016/j.ejor.2011.05.028 doi: 10.1016/j.ejor.2011.05.028
![]() |
[39] | Y. L. Zhao, J. J. Gu, Measurement and comparative research on the development quality of manufacturing industries in China and America, Quant. Econo. Tech. Econo. Res., 35 (2018), 116–133. |
[40] | W. Zt, Research on the cultivation of practical innovation ability of local college students based on CDIO concept, Res. Explor. Lab., 36 (2017), 220–224. |
[41] | L. X. Wang, X. G. Chen, X. L. Yao, X. Y. Li, C. T. Zhang, Research on the responsiveness of Chinese industrial enterprises to environmental regulation policies, China Soft Sci., 10 (2017), 143–152. |
[42] | J. L. Liu, M. S. Ran, Research on the impact of environmental regulations on industrial production technology progress, Sci. Res. Manage., 36 (2015), 107–114. |
[43] |
M. Buesa, J. Heijs, M. M. Pellitero, T. Baumert, Regional systems of innovation and the knowledge production function: the spanish case, Technovation, 26 (2006), 463–472. https://doi.org/10.1016/j.technovation.2004.11.007 doi: 10.1016/j.technovation.2004.11.007
![]() |
[44] |
L. Luo, S. Liang, Study on the efficiency and regional disparity of green technology innovation in China's industrial companies, Chin. J. Popul. Resour. Environ., 14 (2016), 262–270. https://doi.org/10.1080/10042857.2016.1258799 doi: 10.1080/10042857.2016.1258799
![]() |
[45] |
J. Quan, B. Zeng, D. Liu, Green supplier selection for process industries using weighted grey incidence decision model, Complexity, 2018 (2018). https://doi.org/10.1155/2018/4631670 doi: 10.1155/2018/4631670
![]() |
[46] | X. Tang, N. Xie, Research on the evaluation of tourism development potential of tea intangible cultural heritage based on grey clustering, Grey Syst. Theory Appl., 9 (2019), 295–304. |
[47] |
B. Zeng, X. Ma, M. Zhou, A new-structure grey Verhulst model for China's tight gas production forecasting, Appl. Soft Comput., 96 (2020), 106600. https://doi.org/10.1016/j.asoc.2020.106600 doi: 10.1016/j.asoc.2020.106600
![]() |
[48] |
B. Zeng, H. Li, Prediction of coalbed methane production in China based on an optimized grey system model, Energy Fuels, 35 (2021), 4333-4344. https://doi.org/10.1021/acs.energyfuels.0c04195 doi: 10.1021/acs.energyfuels.0c04195
![]() |
[49] | X. F. Pan, Y. Yang, Evaluation of industrial enterprise's innovation efficiency in China excluding the influence of environment variables, Oper. Res. Manage. Sci., 23 (2014), 244-251. |
[50] |
Y. Deng, Y. Wu, H. Xu, Political connections and firm pollution behaviour: an empirical study, Environ. Resour. Econ., 75 (2020), 867-898. https://doi.org/10.1007/s10640-020-00410-7 doi: 10.1007/s10640-020-00410-7
![]() |
[51] |
W. Gu, T. L. Saaty, L. Wei, Evaluating and optimizing technological innovation efficiency of industrial enterprises based on both data and judgments, Int. J. Inf. Technol. Decis. Making, 17 (2018), 9-43. https://doi.org/10.1142/S0219622017500390 doi: 10.1142/S0219622017500390
![]() |
[52] |
L. Fraccascia, V. Albino, C. A. Garavelli, Technical efficiency measures of industrial symbiosis networks using enterprise input-output analysis, Int. J. Prod. Econ., 183 (2017), 273-286. https://doi.org/10.1016/j.ijpe.2016.11.003 doi: 10.1016/j.ijpe.2016.11.003
![]() |
[53] | L. W. Luo, S. R. Liang, Green technology innovation efficiency and factor decomposition of China's industrial enterprises, China Popul. Resour. Environ., 26 (2016), 149-157. |
[54] |
C. Rusinko, Green manufacturing: An evaluation of environmentally sustainable manufacturing practices and their impact on competitive outcomes, IEEE Trans. Eng. Manage., 54 (2007), 445-454. https://doi.org/10.1109/TEM.2007.900806 doi: 10.1109/TEM.2007.900806
![]() |
[55] |
H. Forsman, Innovation capacity and innovation development in small enterprises. A comparison between the manufacturing and service sectors, Res. Policy, 40 (2011), 739-750. https://doi.org/10.1016/j.respol.2011.02.003 doi: 10.1016/j.respol.2011.02.003
![]() |
[56] | Y. H. Zhou, X. D. He, Y. Shen, An evaluation of the efficiency of Chinese industry enterprises' innovation performance, Econ. Res. J., 5 (2012), 107–119. |
[57] | Z. J. Feng, W. Chen, Research on the efficiency of R&D and innovation in China's high-tech industry: a new perspective based on the resource-constrained two-stage DEA model, Syst. Eng. Theory Prac., 34 (2014), 1202–1212. |
[58] | C. Wu, S. W. Yang, P. C. Tang, T. Wu, S. K. Fu, Construction of a green innovation efficiency improvement model for China's heavy pollution industries, China Popul. Resour. Environ., 28 (2018), 40–48. |
[59] | C. Zhang, B. N. Guo, T. S. Yu, Pollution heterogeneity, optimal environmental regulation intensity and production technology progress, Sci. Res. Manage., 36 (2015), 138–144. |
[60] |
N. Yalcin, A. Bayrakdaroglu, C. Kahraman, Application of fuzzy multi-criteria decision making methods for financial performance evaluation of Turkish manufacturing industries, Expert Syst. Appl., 39 (2011), 350–364. https://doi.org/10.1016/j.eswa.2011.07.024 doi: 10.1016/j.eswa.2011.07.024
![]() |
[61] | Z. Liu, Y. G. Dang, W. Y. Qian, W. J. Zhou, Grey grid correlation model based on panel data, Syst. Eng. Theory Prac., 34 (2014), 991–996. |
[62] |
B. J. Frey, D. Dueck, Clustering by passing messages between data points, Science, 315 (2007), 972–976. https://doi.org/10.1126/science.1136800 doi: 10.1126/science.1136800
![]() |
Require: number of incoming roads Ensure: approximate solutions 1: Initilization: 2: 3: 4: 5: for 6: Solve the by definition 3.2 induced optimization problem at the junction based on the flux 7: Compute the densities at the junction with an appropriate Riemann solver 8: Compute the adjusted flux values for the incoming roads 9: for 10: 11: for 12: 13: 14: end for 15: 16: for 17: 18: end for 19: end for 20: for 21: Compute 22: for 23: 24: 25: end for 26: 27: for 28: 29: end for 30: end for 31: end for |
Require: Demand Ensure: Flux values 1: 2: 3: if 4: 5: else if 6: 7: 8: end if |
Require: Demand Ensure: Flux values 1: 2: 3: if 4: if 5: 6: else if 7: 8: end if 9: else if 10: if 11: 12: else 13: 14: end if 15: if 16: 17: 18: else if 19: 20: else if 21: 22: end if 23: end if |
Require: Demand Ensure: Flux values 1: 2: 3: if 4: 6: 7: end if 8: 9: if 10: 11:else if 12: if 13: 14: end if 15: if 16: if 17: 18: else 19: 20: end if 21: end if 22: if 23: if 24: 25: else 26: 27: end if 28: end if 29: end if |
Example 1 | |||||
Splitt. | Reg. | ||||
33.44e-03 | 46.77e-03 | 82.26e-03 | 51.82e-03 | ||
24.17e-03 | 29.05e-03 | 65.59e-03 | 30.69e-03 | ||
14.16e-03 | 20.12e-03 | 60.86e-03 | 24.45e-03 | ||
8.97e-03 | 12.49e-03 | 58.37e-03 | 20.44e-03 | ||
CR | 0.64695 | 0.62453 | 0.1593 | 0.4353 | |
Example 2 | |||||
Splitt. | Reg. | ||||
4.58e-03 | 7.41e-03 | 44.70e-03 | 14.57e-03 | ||
2.97e-03 | 4.24e-03 | 43.47e-03 | 11.28e-03 | ||
2.03e-03 | 2.89e-03 | 42.50e-03 | 10.06e-03 | ||
1.24e-03 | 1.99e-03 | 41.80e-03 | 9.21e-03 | ||
CR | 0.61911 | 0.62327 | 0.0322 | 0.2150 |
Example 1 | |||||
Splitt. | Reg. | ||||
9.25e-03 | 16.22e-03 | 16.22e-03 | 17.01e-03 | ||
5.90e-03 | 11.63e-03 | 11.63e-03 | 12.19e-03 | ||
2.98e-03 | 8.13e-03 | 8.13e-03 | 8.52e-03 | ||
8.97e-03 | 5.71e-03 | 5.71e-03 | 5.99e-03 | ||
CR | 0.53838 | 0.50353 | 0.50353 | 0.50347 | |
Example 2 | |||||
Splitt. | Reg. | ||||
14.12e-03 | 20.10e-03 | 85.69e-03 | 27.55e-03 | ||
9.65e-03 | 13.86e-03 | 79.98e-03 | 21.65e-03 | ||
6.41e-03 | 9.57e-03 | 75.96e-03 | 17.49e-03 | ||
4.51e-03 | 6.69e-03 | 73.22e-03 | 14.65e-03 | ||
CR | 0.55295 | 0.52959 | 0.07551 | 0.30432 |
Require: number of incoming roads Ensure: approximate solutions 1: Initilization: 2: 3: 4: 5: for 6: Solve the by definition 3.2 induced optimization problem at the junction based on the flux 7: Compute the densities at the junction with an appropriate Riemann solver 8: Compute the adjusted flux values for the incoming roads 9: for 10: 11: for 12: 13: 14: end for 15: 16: for 17: 18: end for 19: end for 20: for 21: Compute 22: for 23: 24: 25: end for 26: 27: for 28: 29: end for 30: end for 31: end for |
Require: Demand Ensure: Flux values 1: 2: 3: if 4: 5: else if 6: 7: 8: end if |
Require: Demand Ensure: Flux values 1: 2: 3: if 4: if 5: 6: else if 7: 8: end if 9: else if 10: if 11: 12: else 13: 14: end if 15: if 16: 17: 18: else if 19: 20: else if 21: 22: end if 23: end if |
Require: Demand Ensure: Flux values 1: 2: 3: if 4: 6: 7: end if 8: 9: if 10: 11:else if 12: if 13: 14: end if 15: if 16: if 17: 18: else 19: 20: end if 21: end if 22: if 23: if 24: 25: else 26: 27: end if 28: end if 29: end if |
Example 1 | |||||
Splitt. | Reg. | ||||
33.44e-03 | 46.77e-03 | 82.26e-03 | 51.82e-03 | ||
24.17e-03 | 29.05e-03 | 65.59e-03 | 30.69e-03 | ||
14.16e-03 | 20.12e-03 | 60.86e-03 | 24.45e-03 | ||
8.97e-03 | 12.49e-03 | 58.37e-03 | 20.44e-03 | ||
CR | 0.64695 | 0.62453 | 0.1593 | 0.4353 | |
Example 2 | |||||
Splitt. | Reg. | ||||
4.58e-03 | 7.41e-03 | 44.70e-03 | 14.57e-03 | ||
2.97e-03 | 4.24e-03 | 43.47e-03 | 11.28e-03 | ||
2.03e-03 | 2.89e-03 | 42.50e-03 | 10.06e-03 | ||
1.24e-03 | 1.99e-03 | 41.80e-03 | 9.21e-03 | ||
CR | 0.61911 | 0.62327 | 0.0322 | 0.2150 |
Example 1 | |||||
Splitt. | Reg. | ||||
9.25e-03 | 16.22e-03 | 16.22e-03 | 17.01e-03 | ||
5.90e-03 | 11.63e-03 | 11.63e-03 | 12.19e-03 | ||
2.98e-03 | 8.13e-03 | 8.13e-03 | 8.52e-03 | ||
8.97e-03 | 5.71e-03 | 5.71e-03 | 5.99e-03 | ||
CR | 0.53838 | 0.50353 | 0.50353 | 0.50347 | |
Example 2 | |||||
Splitt. | Reg. | ||||
14.12e-03 | 20.10e-03 | 85.69e-03 | 27.55e-03 | ||
9.65e-03 | 13.86e-03 | 79.98e-03 | 21.65e-03 | ||
6.41e-03 | 9.57e-03 | 75.96e-03 | 17.49e-03 | ||
4.51e-03 | 6.69e-03 | 73.22e-03 | 14.65e-03 | ||
CR | 0.55295 | 0.52959 | 0.07551 | 0.30432 |