SciELO - Scientific Electronic Library Online

 
vol.17 número4The design of a liquefied natural gas (LNG) distribution network of a company operating in MexicoGrid connected DFIG based wind energy conversion system using nine switch converter índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Journal of applied research and technology

versión On-line ISSN 2448-6736versión impresa ISSN 1665-6423

J. appl. res. technol vol.17 no.4 Ciudad de México jul./ago. 2019  Epub 28-Feb-2021

https://doi.org/10.22201/icat.16656423.2019.17.4.794 

Articles

Economic production quantity model with backorders and items with imperfect/perfect quality options

Behrouz Afshar-Nadjafia  * 

Hamed Pourbakhshb 

Mohammad Mirhabibib 

Hasan Khodaeib 

Babak Ghodamib 

Fahimeh Sadighib 

Soheil Azizib 

a Department of Industrial Engineering, Faculty of Industrial and Mechanical Engineering, Qazvin Branch, Islamic Azad University, Qazvin, Iran.

b Industrial Management Department, Faculty of Management and Accounting, Allameh Tabatabaei University, Tehran, Iran.


Abstract

In this paper, an economic production quantity (EPQ) model with backorders considering two options for replenishing of items is proposed; with partially imperfect and perfect quality items. First option assumes a produced shipment contains a fraction of imperfect quality items and supplier does not conduct a full inspection. Therefore, these items are detected by a fully perfect screening process by buyer and are sold as a single batch at a discounted price. While the second one assumes that all items that are produced are fully inspected by supplier and all delivered items are perfect; of course with higher unit production price. Ordering size and backordering level are used as decision variables to derive the closed-form optimal solution. The proposed model is illustrated and discussed by a numerical example. Finally, a sensitivity analysis is done for identifying the impact of crucial parameters on the optimal solution.

Keywords: Inventory Control; Economic Production Quantity; Backorders; Imperfect Items

1. Introduction

The economic production quantity (EPQ) model has been widely used in practice because of its simplicity. However, there are some shortcomings in the assumptions of the basic EPQ model. Two important assumptions of the basic EPQ model are that all the items are of perfect quality and the no shortage constraint. Recently, the classical EPQ model has been generalized in many directions and many researchers have tried to improve it with different viewpoints. Some authors extended the EPQ model by relaxing these restrictions and extend the basic EPQ model to solving inventory models under various realistic situations such as backordering, inflation, varying setup cost, rework process and etc.

Some researchers have studied EPQ model with backordering, including Cárdenas-Barrón (2011), Omar, Zubir, and Moin (2010), Pentico, Drake, and Toews (2009), Sarkar and Moon (2013), Sphicas (2006) and Wee, Huang, Wang, & Cheng (2014). Yang, Teng, and Chern (2010), studied an inventory model under inflation for deteriorating items with stock-dependent consumption rate and partially backordering. Darwish (2008) developed the EPQ model by considering a relationship between the setup cost and the production run length. The effect of learning and forgetting in setups and in product quality is investigated by Jaber and Bonney (2003). Jaber (2006) investigated the lot sizing problem for reduction in setups with reworks and interruptions to restore the process to an “in-control’’ state. Afshar and Abbasi (2009) considered an EPQ model with depreciation cost and process quality cost as continuous functions of time. Tettech, Xu, and Liu (2014) studied the effect of speculation in a dual supply chain in controlling cost of inventory. Wee, Wang, and Cardenas-Barron (2013) generalized an EPQ model with rework process at a single-stage manufacturing system with two backordering costs on imperfect manufacturing system. They proposed other decision variables, the production time and the time to eliminate backorders, instead of the classical decision variables of lot size and backorder quantities. Chiu, Tseng, Wu, and Sung (2014) considered a multi-item EPQ with scrap, rework and multiple deliveries. Also, Chiu, Wu,Chiu, and Chang (2015) applied mathematical modelling for determining the optimal production length of the problem.

Many other researchers have studied EPQ model with imperfect production process including Afshar-Nadjafi et al. (2018), Hemmati and Afshar-Nadjafi (2017), Jaber, Zanoni, and Zavanella (2014), Liao and Sheu (2011), Vörös (2013) and Widyadana and Wee (2012). Yoo, Kim, and Park (2009) developed an EPQ model with imperfect production quality, imperfect inspection and rework. Yassine, Maddah, and Salameh (2012) analyzed shipment of imperfect quality items during a single production runs and over multiple production runs. Rezaei and Salimi (2012) studied the relationship between buyer and supplier with regard to conducting the inspection and resulting in a change the buyer’s economic order quantity and purchasing price. This work will be the starting point of the model developed in this paper.

To the author’s knowledge, none of the above studies considered the economic production quantity (EPQ) model with backorders considering quality-oriented options for replenishing of items. In this paper, the basic EPQ model is extended by considering backorders and two options for replenishing of items; with partially imperfect and perfect quality. A simple method is used to solve the extended model of maximizing the total annual profit. Also, numerical examples are used to show the utility of the proposed models. The organization of the paper is as follows. The model assumptions and notation are presented in Section 2. Then, the model is developed in Section 3. In section 4 a summary of the results given in this paper is explained based on a numerical example. Finally, section 5 contains the conclusions.

2. Assumptions and notations

The basic EPQ model is that of determining a production quantity of an item, subject to the following assumptions:

  • Demand rate is continuous, known and constant.

  • Production rate is greater than or equal to demand rate.

  • All demands must be met.

  • Holding costs are determined by the value of the item.

  • There is no discount.

  • There are no quantity constraints.

  • All items in a batch conform to quality characteristics.

  • No shortages are allowed.

In this section, we derive a mathematical statement for the EPQ model with backordering and quality-oriented options for replenishing of items. Indeed, we consider two scenarios for replenishing of items. First scenario assumes a produced shipment contains a fraction of imperfect quality items and these items are detected by 100% screening by buyer and are sold as a single batch at a discounted price. The second scenario assumes that all produced items are fully inspected by supplier and all delivered items are perfect; of course with higher unit production price. This scenario is the traditional EPQ model in which studied in the literature by relaxing the last assumption (Rezaei & Salimi, 2012).

In the first scenario, we remove the two last assumptions of the basic model considering the following assumptions:

  • A shipment contains a fraction of imperfect items.

  • Imperfect items are identified via a perfect inspection process

  • Shortages are allowed as backlogged.

  • Inspection rate is higher than the production rate.

In the second scenario, the same assumptions of the classical economic production quantity as mentioned above are considered. More specifically, it is assumed that all the items are fully inspected by supplier, therefore all the items are perfect.

In this paper, we examine which scenario is preferred by the buyer?

It is clear that, if the unit price of items is the same, the answer is simply to select the second scenario. However, in reality the unit price of items offered by the supplier in second scenario is higher than the unit price of items in offered by the first scenario, which results in a trade-off between quality and price. To answer the above mentioned question, we need to identify the maximum production price (Mc) the buyer is willing to pay to a supplier belonging to the second scenario. If the unit price is less than Mc, the second scenario is preferred, if not the first scenario is preferred (Chung & Cárdenas-Barrón, 2012).

In order to state the problem mathematically, let:

y

production quantity (real positive decision variable)

b

the maximal backorder level (real positive decision variable)

D

annual demand rate of product

U

annual production rate of product

c

regular unit production cost of product

c'

unit production cost of perfect product

Mc

maximum unit production price of perfect product

d

inspection cost per unit

s

the perfect unit selling price

v

the imperfect unit selling price

π

the unit fixed backorder cost per unit

π^

annual unit backorder cost

h

annual unit holding cost

K

fixed setup cost of production system

p

imperfect rate of production system (random variable)

E(p)

expected value of imperfect rate of production,system

I max

maximum on -hand inventory level

T

the length of the inventory cycle

t 1

time to eliminate the backorder position

t 2

time to build an inventory of I max from a zero position

t 3

time to consume the maximum inventory

t 4

time to build a backorder position

NS(t)

net stock level at time t

TP(y,b)

total profit per ordering cycle

ETPU(y,b)

expected annual total profit (objective function of first scenario)

TPU(y,b)

annual total profit (objective function of second scenario)

3. Model development

We depict the proposed inventory system (first scenario) in Fig. 1. As shown in Fig. 1, an inventory cycle, i.e., T, can be split into four time-intervals.

Fig. 1. The relation between net stock level and time. 

From graphical representation of model in Fig. 1, we have:

t1=bU-D (1)

t2=y 1- DU -bU-D (2)

t3=y 1- DU - p-bD (3)

t4=bD (4)

Imax=y 1- DU - b (5)

The buyer’s expected profit per ordering cycle is as follows:

TP(y,b) = sales of perfect quality items + sales of imperfect quality items - variable production cost - fixed setup cost - inspection cost - holding costs - shortage costs

Equivalently, we have:

TPy,b=sy1-p+vyp-cy-K-dy-πb-π^b22D1-DU-hy1-DU-b22U-D-hy1-DU-b-py22D (6)

Then the buyer’s expected total profit per time unit (annual) can be stated by replacing p by E(p) and dividing it by the expected length of the production cycle ET=(1-E(p))yD as follows:

ETPUy,b=sD1-E(p)+vDE(p)-cD-KDy-dD- πDby-π^b22y1-DU-hDy1-DU-b22yU-D-hy1-DU-b-E(p)y22y/(1-E(p)) (7)

Set ETPU(y,b)b equals zero, it yields:

b*=-πD+hy*1-E(p)1-DUπ^+h (8)

Substituting b * in ETPU(y,b)y=0, we have:

y*=2DKπ^+h1-DU-π2D22hπ^+hEp+0.5+E2p21-DU-h21-E(p)2 (9)

The elements of Hessian matrix are:

2ETPU(y,b)b2=-(π^+h)y1-DU(1-E(p)) (10)

2ETPU(y,b)by=2ETPU(y,b)yb=1y2(1-E(p))πD+(π^+h)b1-DU (11)

And

2ETPU(y,b)y2=-1y3(1-E(p))2D(K+πb)+(π^+h)b21-DU (12)

The values α 1 and α 2 can be computed as follows:

α1=2ETPU(y,b)b20     ;      b , y (13)

α2=det2ETPUy,b=det2ETPU(y,b)b22ETPU(y,b)by2ETPU(y,b)yb2ETPU(y,b)y2=(π^+h)y41-DU(1-E(p))22DK-π2D21-DUπ^+h  ;  free in sign (14)

Evaluation of α 1 and α 2 shows that the ETPU (y, b) is a non-concave function. Therefore, taking the partial derivatives of ETPU (y, b) with respect to the b and y, and solving equations ETPU(y,b)b=0 and ETPU(y,b)y=0 don’t guarantee the necessary conditions for b * and y * to be optimal. However, we can attain the optimal policy for the first scenario as follows:

If π2D2min{2DKπ^+h1-DU,2DKh1-Ep22Ep+11-DU+E2p}, the optimal solution to maximize ETPU (y, b) is (b *, y *) obtained from Eqs. (8) and (9), else, the optimal policy is to fill all demand without backorders.

According to Eq. (14) if π2D22DK(π^+h)1-DU, we have α 2 ≥ 0, which guarantees the concavity of ETPU (y, b). Also, by substituting Eq. (9) in Eq. (8) and setting b * ≥ 0 , we have π2D22DKh(1-E(p))22Ep+11-DU+E2(p) which guarantees feasibility (non-negativity) of b * obtained from Eq (8). If we have π2D2>2DKh(1-E(p))22Ep+11-DU+E2(p), we get b * < 0, then (y *, b *) is never the optimal solution of ETPU (y, b) on y > 0 and b ≥ 0. So, if the optimal solution of ETPU (y, b) on y > 0 and b ≥ 0 exist, then, b * = 0.

Notice that, in special case of π=0, concavity of ETPU (y, b) and feasibility (non-negativity) of b * is guaranteed. In this case, that b * and y * are computed as follows:

b*=hy*1-E(p)1-DUπ^+h (15)

y*=2DKπ^+h2hπ^+hEp+0.51-DU+0.5E2p-h21-E(p)21-DU (16)

To verify our model by reducing it to the model with perfect items and backordering not considering fixed backordering cost, we set π = 0 and E(p) = 0. One can show that b * and y * are computed as follows which are identical to the results in Chung and Cárdenas-Barrón (2012):

b*=2DKh1-DUπ^(π^+h) (17)

y*=2DKh1-DUπ^+hπ^ (18)

For second scenario, we assume the unit price (c´) paid for fully perfect items to be higher than that of a batch with some imperfect items; c'>c. Given fully perfect items (p = 0), which result to zero inspection for the buyer, the total profit per time unit is:

TPUy,b=s-c'D-DyK+πb-π^b2+hy1-DU-b22y1-DU (19)

In this case, the optimal solution to maximize TPU (y, b) is identical to the EPQ inventory model with linear and fixed backorder costs reported in Chung and Cárdenas-Barrón (2012):

b**=hy*-πD1-DUπ^+h (20)

y**=2DKh1-DUπ^+hπ^-π2D2hπ^ (21)

Now the maximum purchasing price for full perfect items should be determined. To this end, we first determine the difference between the total profit per time unit when there are no imperfect items and the expected profit per time unit when there are averagely p% imperfect items in each batch. We consider c' as a variable here.

TPUy**,b**,c'-ETPUy*,b*=s-c'D-Dy**K+πb**-π^b**2+hy**1-DU-b**22y**1-DU-ETPUy*,b* (22)

The buyer accepts to pay more if and only if

TPUy**,b**,c'-ETPUy*,b*0 (23)

Or equivalently:

c's-1y**K+πb**-π^b**2+hy**1-DU-b**22Dy**1-DU-ETPUy*,b*D (24)

The right-hand side of this equation determines the maximum unit purchasing price (Mc) the buyer is willing to pay for batches without imperfect items.

Now, we provide the following rule to answer the research question of this paper: if the item price for the second scenario is less than or equals to the maximum purchasing price Mc computed from Eq. (24), it is wise to select the second scenario (who carries out a full inspection and there are no imperfect items), otherwise the first scenario (here there are averagely p% imperfect items in each batch and the inspection is conducted by the buyer) is preferred.

4. Numerical Example and Discussion

To illustrate the model developed in section 3, the inventory situation with following data is considered: U = 12,000 unit/year, D = 10,000 units/year, K = 450 $/cycle, π^=0.5$/unit/year, π = 1.2$/ unit, h = 75$/unit/year, d=5 $/ unit, $/ unit, s = 22, v=30 $/ unit, c'=125 $/ unit, c = 100 $/ unit and E (p) = 6%. Using these data for the first scenario gives y * = 1691.17, b * = 236.71 and ETPU (y * , b * ) = 1,010,030.28 whereas with the second scenario, we have y * = 10241.09, b ** = 1669.06, TPUy**,b**=947,165.47 and Mc = 110. It is clear that the first scenario is preferred in which there are averagely 15% imperfect items in each batch and the inspection is conducted by the buyer. However, the optimal backorder level, b*, and the optimal order quantity, y*, are decreased in comparison with the second scenario. In this example Mc = 110 means that if the item price for the second scenario is less than or equals to Mc = 110, this scenario will be preferred.

In order to investigate the effect of imperfect rate of production system, p, the sensitivity analysis is performed by changing E (p). The remaining parameters are kept unchanged.

Table 1 gives the optimal solutions for selected values of E (p) ranging from 0 to 0.30. The results show that the economic production quantity, economic backordering level and the optimal buyer’s expected total profit per time unit (annual) decrease as E (p) increases, while Mc is an increasing function of E (p). These results are demonstrated in Figs. 2, 3 and 4. Fig. 2 shows that as E (p) increases, optimal solution goes to infeasibility (negativity of b *). In Fig. 2, the initial point, E (p) = 0, is the optimal solution related to the supplier with full perfect items (second scenario). It can be observed that the supplier with full perfect items is always characterized by a higher lot size and higher backordering level compared to the supplier with partially imperfect items (first scenario).

Table 1. Optimal solutions for different values of E(p). 

E (p) y * b * ETPU * Mc

  • 0

  • 0.02

  • 0.04

  • 0.06

  • 0.08

  • 0.1

  • 0.12

  • 0.14

  • 0.16

  • 0.18

  • 0.2

  • 0.22

  • 0.24

  • 0.26

  • 0.28

  • 0.3

  • 10241.09

  • 2841.45

  • 2053.14

  • 1691.17

  • 1472.65

  • 1322.65

  • 1211.63

  • 1125.24

  • 1055.59

  • 997.93

  • 949.20

  • 907.33

  • 870.87

  • 838.75

  • 810.20

  • 784.61

  • 1669.06

  • 434.54

  • 299.84

  • 236.71

  • 197.82

  • 170.59

  • 150.04

  • 133.73

  • 120.31

  • 108.99

  • 99.23

  • 90.68

  • 83.09

  • 76.27

  • 70.09

  • 64.44

  • 1147165.47

  • 1130996.68

  • 1114372.28

  • 1097141.96

  • 1079215.51

  • 1060524.62

  • 1041005.39

  • 1020592.83

  • 999218.31

  • 976808.08

  • 953281.95

  • 928552.19

  • 902522.27

  • 875085.57

  • 846123.85

  • 815505.57

  • 105

  • 106.6

  • 108.3

  • 110

  • 111.8

  • 113.7

  • 115.6

  • 117.7

  • 119.8

  • 122

  • 124.4

  • 126.9

  • 129.5

  • 132.2

  • 135.1

  • 138.2

Fig. 2. Optimal values of y * 

Fig. 3. Optimal values of Mc against E ( p ). 

Fig. 4. Comparison of two scenarios with different values of E ( p ). 

Fig. 3 shows that as E (p) increases, the buyer is willing to pay more to a supplier belonging to the second scenario to avoid receiving imperfect items. Paying this additional amount implies that the supplier conducts the inspection process. Fig. 4 shows that for low values of E (p) the supplier belonging to the first scenario is preferred. However, as E (p) increases, the profitability of the first scenario decreases rapidly. In our example for p>20%, second scenario is preferred.

An interesting analysis can be carried out by considering changing the purchasing cost of an perfect item, i.e., C', while keeping all other parameters unchanged. The result is reported in Fig. 5 , where it can be seen that for values of C' smaller than Mc = 110 it is more profitable to select a supplier belonging to the second scenario, otherwise it is more profitable to select first one.

Fig. 5. Comparison of two scenarios with different values of C ´. 

5. Conclusions

In this paper, a variation of the EPQ model have been developed considering backorders and two options for replenishing of items; with partially imperfect and perfect quality. The option with perfect quality items is simply the traditional EPQ with backorder in which studied in the literature. In this paper, we have formulated and solved the option with imperfect quality items to determine the economic production quantity, economic backordering level and maximum purchasing price a buyer is willing to pay to a supplier to avoid receiving imperfect items. The problem described with a mathematical model, and the existence of the optimal policy is investigated. Also, taking the limiting parameter values for the optimal solutions, our model can be generalized to the classic EPQ model with backorders and perfect items. The numerical example is provided to illustrate the developed model. From the numerical results, we could clearly see that loss due to using the classical EPQ model with backorders is significant. The sensitivity analysis showed that the profitability of a supplier with partially imperfect quality items decreases with increase in the value of the purchasing cost of a perfect item. Also, the results showed that the imperfect rate of production system, p, has significant effect on optimal policies. The results of this study can help managers make optimal decisions on supplier selection. Further research can investigate the case in which there is rework option for imperfect items and shortages are partially backordered.

References

Afshar-Nadjafi, B., & Abbasi, B. (2009). EPQ model with depreciation cost and process quality cost as continuous functions of time. Journal of Industrial Engineering International, 5(8), 77-89. [ Links ]

Afshar-Nadjafi, B., Hadinejad, F., Khajeh, M., Eftekhar, M., Fakhim-Hashemi, N., & Sadeghian, A. (2018). A Mathematical Approach to Derive the Optimal Ordering Policy with Backorders and Inspection. Chiang Mai Journal of Science, 45(1), 607-616. [ Links ]

Cárdenas-Barrón, L. E. (2011). The derivation of EOQ/EPQ inventory models with two backorders costs using analytic geometry and algebra.Applied Mathematical Modelling,35(5), 2394-2407. [ Links ]

Chiu, Y.S.P., Tseng, C.T., Wu, M.F., & Sung, P.C. (2014). Multi-item EPQ model with scrap, rework, and multi-delivery using common cycle policy. Journal of Applied Research and Technology, 12(3), 615-622. [ Links ]

Chiu, Y.S.P., Wu, M.F., Chiu, S.W., & Chang, H.H. (2015). A simplified approach to the multi-item economic production quantity model with scrap, rework, and multi-delivery. Journal of Applied Research and Technology, 13(4), 472-476. [ Links ]

Chung, K. J., & Cárdenas-Barrón, L. E. (2012). The complete solution procedure for the EOQ and EPQ inventory models with linear and fixed backorder costs.Mathematical and Computer Modelling,55(11-12), 2151-2156. [ Links ]

Darwish, M. A. (2008). EPQ models with varying setup cost.International Journal of Production Economics,113(1), 297-306. [ Links ]

Hemmati, A., & Afshar-Nadjafi, B. (2017). The influence of processing on the optimal inventory policy for deteriorating products.Latin American Applied Research,47(3-4), 137-143. [ Links ]

Jaber, M.Y., & Bonney, M. (2003). Lot sizing with learning and forgetting in set-ups and in product quality. International Journal of Production Economics, 83(1), 95-111. [ Links ]

Jaber, M.Y. (2006). Lot sizing for an imperfect production process with quality corrective interruptions and improvements, and reduction in setups. Computers and Industrial Engineering, 51(4), 781-790. [ Links ]

Jaber, M.Y., Zanoni, S., & Zavanella, L.E. (2014). Economic order quantity models for imperfect items with buy and repair options. International Journal of Production Economics, 155,126-131. [ Links ]

Liao, G.L., & Sheu, S.H. (2011). Economic production quantity model for randomly failing production process with minimal repair and imperfect maintenance. International Journal of Production Economics, 130(1), 118-124. [ Links ]

Omar, M., Zubir, M.B., & Moin, N.H. (2010). An alternative approach to analyze economic ordering quantity and economic production quantity inventory problems using the completing the square method. Computers and Industrial. Engineering, 59(2), 362-364. [ Links ]

Pentico, D.W., Drake, M.J., & Toews, C. (2009). The deterministic EPQ with partial backordering: a new approach. Omega, 37(3), 624-636. [ Links ]

Rezaei, J., & Salimi, N. (2012). Economic order quantity and purchasing price for items with imperfect quality when inspection shifts from buyer to supplier. International Journal of Production Economics, 137(1), 11-18. [ Links ]

Sarkar, B., & Moon, I. (2014). Improved quality, setup cost reduction, and variable backorder costs in an imperfect production process. International Journal of Production Economics, 155, 204-213. [ Links ]

Sphicas, G.P. (2006). EOQ and EPQ with linear and fixed backorder costs: two cases identified and models analyzed without calculus. International Journal of Production Economics, 100(1), 59-64. [ Links ]

Tettech, A., Xu, Q., & Liu, Z. (2014). Inventory control by using speculative strategies in dual channel supply chain. Journal of Applied Research and Technology, 12(2), 296-314. [ Links ]

Vörös, J. (2013). Economic order and production quantity models without constraint on the percentage of defective items. Central European Journal of Operations Research, 21(4), 867-885. [ Links ]

Wee, H.M., Wang, W.T., & Cardenas-Barron, L.E. (2013). An alternative analysis and solution procedure for the EPQ model with rework process at a single-stage manufacturing system with planned backorders. Computers and Industrial Engineering, 64(2), 748-755. [ Links ]

Wee, H.M., Huang, Y.D., Wang, W.T., & Cheng, Y.L. (2014). An EPQ model with partial backorders considering two backordering costs. Appllied Mathematics and Computation, 232, 898-907. [ Links ]

Widyadana, G. A., & Wee, H. M. (2012). An economic production quantity model for deteriorating items with multiple production setups and rework.International Journal of Production Economics,138(1), 62-67. [ Links ]

Yang, H. L., Teng, J. T., & Chern, M. S. (2010). An inventory model under inflation for deteriorating items with stock-dependent consumption rate and partial backlogging shortages.International Journal of Production Economics,123(1), 8-19. [ Links ]

Yassine, A., Maddah, B., & Salameh, M. (2012). Disaggregation and consolidation of imperfect quality shipments in an extended EPQ model. International Journal of Production Economics, 135(1), 345-352. [ Links ]

Yoo, S. H., Kim, D., & Park, M. S. (2009). Economic production quantity model with imperfect-quality items, two-way imperfect inspection and sales return.International Journal of Production Economics,121(1), 255-265. [ Links ]

Peer Review under the responsability on Universidad Nacional Autónoma de México.

* Corresponding author. E-mail address: afsharnb@alum.sharif.edu (Behrouz Afshar-Nadjafi).

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License