Yonghui Shu is currently a Master’s student under the tutelage of Prof. Peng Hu at the University of Science and Technology of China. His research interests include refrigeration and cryogenics
Peng Hu received his Ph.D. degree in Engineering Thermophysics from the University of Science and Technology of China (USTC). He is currently a professor at USTC. His research interests include thermophysical properties, refrigeration, heat and mass transfer, and solar energy
Capillary tubes have been widely used as expansion devices in small-scale refrigeration and heat-pump systems. However, adiabatic flow through a capillary tube is extremely complicated, despite its simple geometry. This work presents a comparative study on the homogenous flow model and separated flow model, which were used to simulate the flow of isobutene (R600a) through adiabatic capillary tubes. The influence of different combinations of friction factor and two-phase viscosity correlations, and the effect of metastable flow on the flow characteristics were investigated. The predicted mass flow rate was lower when the separated flow model was used. The separated flow model performed better in predicting a mass flow rate over 2 kg·h−1. The Colebrook friction factor correlation combined with the Dukler or McAdams viscosity correlation yielded smaller deviations of 5.43%, 5.49% and 5.44%, 5.43% when ignoring and considering the metastable flow, respectively. Additionally, the homogenous flow model adopting the Bittle and Pate friction factor and Dukler viscosity correlations achieved the highest accuracy with a mass flow rate under 2 kg·h−1. The mean error was 4.12% in the case without metastable flow, and 3.37% in the case with metastable flow.
Graphical Abstract
Simulation on the throttling process of refrigerant flowing through adiabatic capillary tubes.
Abstract
Capillary tubes have been widely used as expansion devices in small-scale refrigeration and heat-pump systems. However, adiabatic flow through a capillary tube is extremely complicated, despite its simple geometry. This work presents a comparative study on the homogenous flow model and separated flow model, which were used to simulate the flow of isobutene (R600a) through adiabatic capillary tubes. The influence of different combinations of friction factor and two-phase viscosity correlations, and the effect of metastable flow on the flow characteristics were investigated. The predicted mass flow rate was lower when the separated flow model was used. The separated flow model performed better in predicting a mass flow rate over 2 kg·h−1. The Colebrook friction factor correlation combined with the Dukler or McAdams viscosity correlation yielded smaller deviations of 5.43%, 5.49% and 5.44%, 5.43% when ignoring and considering the metastable flow, respectively. Additionally, the homogenous flow model adopting the Bittle and Pate friction factor and Dukler viscosity correlations achieved the highest accuracy with a mass flow rate under 2 kg·h−1. The mean error was 4.12% in the case without metastable flow, and 3.37% in the case with metastable flow.
Public Summary
Homogeneous flow models and separated flow models with and without metastable flow are compared in detail.
The optimal hypothesis models for different R600a mass flow rate ranges are investigated.
The optimal combinations of friction factor and two-phase viscosity correlations for each model are recommended.
Barium (Ba) is a highly incompatible element during mantle melting, thus, it is much more enriched in the upper continental crust[1] and marine sediments[2] than in the mantle[3]. Ba is also fluid mobile during slab subduction[4]. Ba has seven stable isotopes: 130Ba (0.11%), 132Ba (0.10%), 134Ba (2.42%), 135Ba (6.59%), 136Ba (7.85%), 137Ba (11.23%), and 138Ba (71.70%)[5]. Previous studies have found that Ba isotope compositions (δ138/134Ba = [(138/134Basample) / (138/134BaSRM3104a) − 1] × 1000‰) of subducted crustal materials such as marine sediment, altered oceanic crust (AOC), and continental rocks[6-9] are distinct from the mantle (0.05‰±0.06‰) [10]. For example, diamictites have a substantial variation in δ138/134Ba from −0.23‰ to +0.47‰ [6]. Therefore, Ba isotopes have recently been suggested as sensitive tracers for crustal material recycling in mantle-derived igneous rocks [7, 8, 11].
Magma produced annually in mid-ocean ridge settings accounts for 75% of global magmatic output[12]. Mid-ocean ridge basalts (MORBs) represent mafic melts derived from a simple melting history of the depleted upper mantle, providing a critical way to constrain the composition of the upper mantle[13]. Substantial chemical heterogeneity has been observed in global MORBs[13-15]. The ratios of highly incompatible trace elements (such as La/Sm) in MORB vary by more than an order of magnitude[13]. A previous study has proven that Ba isotope variation in MORBs is not caused by mantle melting and/or fractional crystallization[7]. Thus, the δ138/134Ba values of the MORBs should reflect their mantle source signatures. It is thus reasonable to use Ba isotope data to identify the contribution of recycled materials to the source of the MORB if the enriched composition in the mantle is related to the addition of materials from surface reservoirs.
Better knowledge of the Ba isotope composition of the upper mantle is essential for the application of Ba isotopes in tracing mantle heterogeneity. However, estimates of δ138/134Ba in the mantle are rare. Nielsen et al.[7] reported δ138/134Ba values of MORB glasses from the Mid-Atlantic Ridge (MAR), East Pacific Rise (EPR), Juan de Fuca Ridge (JdF), Central Indian Ridge, and Southwest Indian Ridge ranging from +0.02‰ to +0.15‰. According to the composition of their most depleted MORBs, they defined the average δ138/134Ba of the depleted MORB mantle (DMM) as ~ +0.14‰. Based on the measurement of mantle-derived carbonatites from Canada, East Africa, Germany and Greenland, Li et al.[10] found that the δ138/134Ba of most carbonatites vary from −0.04‰ to +0.12‰, thus suggesting that the average δ138/134Ba of their mantle source is +0.05‰±0.06‰ (2SD). The average δ138/134Ba of the MORBs (+0.07‰±0.08‰, 2SD) reported by Nielsen et al.[7] is generally indistinguishable from that of the carbonatites reported by Li et al.[10]. However, the δ138/134Ba value of the DMM suggested by Nielsen et al. [7] was much higher. This contradiction hinders the application of Ba isotopes in mantle geochemistry.
In this study, we report new high-precision Ba isotope data for several well-known regions of MORB and back-arc basin basalt (BABB) samples and incorporate these data into a critically compiled database that includes data from a previous study[7]. With these data, we attempted to re-estimate the average Ba isotope composition of the DMM and further constrain the enriched components in the mantle source of the “enriched-type” MORB (E-MORBs) with Ba isotopes.
2.
Sample description
We analyzed 25 MORB samples from different mid-ocean ridge segments. They were from the EPR at 9–10°N, the MAR at 35°N, the Gakkel Ridge (GR) at 82–87°N, and the Southeast Indian Ridge at 50°S (SEIR). These samples were collected from ridge segments with spreading rates varying from lower than 1.5 cm/year to approximately 11 cm/year [12]. In addition, five BABB samples from the Lau Basin at the Australian-Pacific plate boundary were also measured, two of which were andesitic lavas with low MgO content. All samples chosen were fresh lavas[16-18].
According to the relative depletions and enrichments of highly incompatible elements, Arevalo and McDonough[19] divided MORBs into two groups: MORBs with (La/Sm)N[3] (chondrite-normalized La/Sm) < 1 are defined as “normal-type” MORB (N-MORB), and those with (La/Sm)N > 1 are defined as E-MORB. The samples span a broad range of chemical and radiogenic isotope compositions, including the typical N-MORB and E-MORB based on their variations in (La/Sm)N (0.55–3.84). The major and trace element abundances of these samples have been previously characterized [16, 17, 20-22]. In particular, their MgO contents range from 2.7 wt% to 9.1 wt%, SiO2 range from 47.6 wt% to 58.5 wt%, and Ba range from 7 to 453 μg/g. More details (such as the location and chemical composition) are documented in Table 1.
Table
1.
Elemental and isotope compositions for the samples and standard materials in this study.
*87Sr/86Sr values of the samples with asterisks were measured in this study. a La/Sm was normalized to CI-chondrite[3]. b 2SD = two times the standard deviation of the population of n repeat measurements of a sample solution. cn represents the number of repeated measurements of the same solution using MC-ICP-MS. d Data were calculated from the measured δ137/134Ba values, as described in the text. e Replicate = repeat sample dissolution, purified column chemistry, and MC-ICP-MS analysis. PetDB= the Petrological Database website (http://www.earthchem.org/petdb).
Whole rock Ba and Sr isotope analyses were performed in the CAS Key Laboratory of Crust-Mantle and Environments at the University of Science and Technology of China (USTC). Sr was chemically purified following the method reported by Chen et al.[23]. Approximately 120–150 mg of the sample powder was digested in a mixture of concentrated HF-HNO3-HCl. Sr was then separated from the matrix using a quartz column with a 2 mL AG50W-X8 cation exchange resin (200–400 mesh). The total procedural blanks were < 0.5 ng. Ba was chemically purified follow the method reported by Nan et al.[6, 24]. Sample powders containing ~2 μg of Ba were digested in a mixture of concentrated HF-HNO3-HCl. Ba was then separated from the matrix using a Teflon column with a 2 mL AG50W-X12 cation exchange resin (200–400 mesh). To avoid matrix effects from light REEs (e.g., La, Ce, and Nd), a second column with 0.5 mL AG50W-X12 resin was used for further purification. The yield of the purification process was > 99%. The total procedural blanks were < 0.5 ng. Purified samples were dried and diluted with 2% (m/m) HNO3 for instrumental measurements.
Ba and Sr isotope measurements were conducted using a Neptune Plus multi-collector inductively coupled plasma mass spectrometer (MC-ICP-MS). Sr isotope ratios were measured in low-resolution mode under “wet” plasma conditions. The measured Sr isotope ratios were corrected for instrumental mass bias by normalizing them to 86Sr/88Sr = 0.1194. The results for the samples and reference materials are listed in Table 1. The 87Sr/86Sr values of NBS-987 and BCR-2 are 0.710247±0.000011 (2SD, n = 25) and 0.705025±0.000009 (2SD, n = 4), respectively, which are consistent with the previous study[25]. Ba isotope ratios were measured using the double-spike (135Ba-136Ba) method in low-resolution mode under “dry” plasma conditions (Aridus II desolvating nebulizer). The background signal for 137Ba (< 5 mV) is negligible relative to the sample signal (~7 V). The results are reported in δ-notation relative to NIST SRM3104a: δ137/134Ba (‰) = [(137/134Basample)/(137/134BaSRM3104a) − 1] × 1000. For comparison, the δ137/134Ba values of all the samples were calculated to be δ138/134Ba following the mass-dependent fractionation laws (δ138/134Ba ≈ 1.33 × δ137/134Ba). The external precision of δ137/134Ba based on the measurement of two in-house standards (USTC-Ba and ICPUS-Ba) is ≤ 0.04‰ (2SD). We estimate the long-term external precision of δ138/134Ba is ≤ 0.05‰ (2SD), as verified by Deng et al.[26]. The results for the samples and standard materials are listed in Table 1. The δ137/134Ba values of BCR-2 and GSP-2 are in good agreement with previously published values[6, 24].
4.
Results
The Ba-Sr isotope data are presented in Table 1 and Fig. 1. The 87Sr/86Sr values of the MORB and BABB samples range from 0.70243 to 0.70404. The δ138/134Ba of the MORB samples range from −0.06‰ to +0.11‰. E-MORBs (−0.06‰ to +0.11‰, n = 7) show more scattered δ138/134Ba than N-MORBs (−0.04‰ to +0.08‰, n = 18), although the difference is within the analytical uncertainty. There are no systematic cross-regional variations. Therefore, δ138/134Ba varies from −0.03‰ to +0.08‰ in EPR MORBs (n = 13), from −0.02‰ to +0.11‰ in MAR MORBs (n = 6), and from −0.06‰ to +0.10‰ in GR MORBs (n = 6). One MORB sample in SEIR has a value of +0.04‰±0.03‰ (2SD). The δ138/134Ba of the five BABB samples range from +0.01‰ to +0.08‰, which is within the range of MORBs. These values are slightly lower than those of previously reported MORB samples, despite substantial overlapping (+0.02‰ to +0.15‰[7]) (Fig. 2). To our knowledge, the five Lau Basin basalts are the first published Ba isotope data for BABBs. On average, δ138/134Ba of our samples (+0.04‰±0.08‰, 2SD, n = 30) is indistinguishable from their average value (+0.07‰±0.08‰, 2SD, n = 21).
Figure
1.
Correlations between δ138/134Ba and (a) (La/Sm)N, (b) 87Sr/86Sr, and (c) Ba/Th for the mid-ocean ridge basalt (MORB) and back-arc basin basalt (BABB) samples in our study and MORB samples analyzed by Nielsen et al.[7]. Data are from Table 1. Error bars represent 2SD uncertainties. The vertical dotted lines present the defined average (La/Sm)N[16], 87Sr/86Sr[28], and Ba/Th[28] of the depleted MORB mantle. Samples in the orange shade are marked as depleted MORB, and those in the blue shade are marked as normal-type MORB or enriched-type MORB.
Figure
2.
The histogram of Ba isotope compositions of MORBs and BABBs investigated in this study. Literature data of MORBs are from Nielsen et al.[7]. Data are from Table 1.
The basaltic oceanic ridge samples reported herein span a broad range of geographical distributions and chemical compositions. The inconsistent Ba isotope data between this study and that in Nielsen et al.[7] could not reflect the different degrees of rock alteration because all these samples are fresh. Interlaboratory bias is also not considered as Ba isotope analyses of the igneous rock standards in our and their studies display consistent results within error[7, 24]. The most probable explanation is that this difference is due to practical Ba isotope variations among global MORBs.
Barium isotope variation in MORBs can potentially be induced by magmatic differentiation during the formation of these basalts or mantle heterogeneity. Fig. 3 shows that no observable correlations exist between δ138/134Ba and the indicators of magmatic differentiation (MgO and SiO2) or degree of partial melting (Na8)[27]. In addition, as has been verified previously[7], Ba is highly incompatible with the bulk distribution coefficient Dsolid/melt of 0.00012 during mantle melting[28], and ~99% of Ba tends to enter the melt after 1% partial melting. MORBs represent degrees of melting between 6% and 20%[29], as a result, Ba isotope fractionation between the starting mantle composition and these basaltic melts is negligible. Therefore, the Ba isotope compositions of the MORBs principally reflect variations in the δ138/134Ba values of their mantle sources.
Figure
3.
δ138/134Ba versus (a) MgO, (b) SiO2 and (c) Na8 for MORBs and BABBs. Na8 = Na2O + 0.373 × MgO − 2.98 (from Ref. [27]). Data are from Table 1. The error bars represent the 2SD uncertainties.
Together, all the MORB data sets in this study and the literature[7] (51 samples, including BABB samples) provide a range of δ138/134Ba of −0.06‰ to +0.15‰, with an average of +0.05‰±0.09‰ (2SD, n = 51).
5.2
Estimation of the Ba isotope composition of DMM
According to the previous definition of N-MORB [19], 70% of the MORB samples in PetDB are categorized as N-MORB, displaying a considerable range of major and trace element compositions[12]. Only the most depleted samples were representative of the DMM. Therefore, we estimated the δ138/134Ba of DMM based on the most depleted MORB (D-MORB) data. The selected 10 D-MORB samples showed typical depleted compositions with (La/Sm)N < 0.8[16]. We also used 87Sr/86Sr (< 0.70263) and Ba/Th (< 71.3), which are lower than the average DMM[28], as the criteria for recognizing the D-MORBs (Figs. 1b and 1c). Along with the D-MORB samples measured in the literature[7], the average δ138/134Ba of DMM was estimated to be +0.05‰±0.05‰ (2SD, n = 16). The JdF sample (+0.15‰) in Nielsen et al.[7] is beyond the average δ138/134Ba of all MORBs (+0.05‰±0.09‰, 2SD, n = 51) and was thus precluded when calculating the average DMM value. To avoid artificial bias in different laboratories, we also calculated a mean δ138/134Ba solely based on our data of +0.05‰±0.04‰ (2SD, n = 10). Our estimation of δ138/134Ba of DMM is significantly lower than that reported by Nielsen et al. [7] (~ +0.14‰).
It is necessary to evaluate different estimates of the average δ138/134Ba values of the DMM. The estimation of the DMM by Nielsen et al.[7] may be problematic for two reasons. First, the DMM and E-MORB mantle (EMM) were defined using only three MORB samples with the lowest and highest 87Sr/86Sr respectively. The limited sample size leads to a lack of representativeness for the samples that estimate the DMM. Second, as pointed out by Wu et al.[11], using only 87Sr/86Sr is not sufficient to choose appropriate samples to define the DMM. In fact, both the DMM and EMM reservoirs span a large range of radiogenic isotope ratios. In this case, the D-MORBs with 87Sr/86Sr < 0.7024 exhibited a large δ138/134Ba range (+0.03‰ to +0.14‰), overlapping with the δ138/134Ba range (+0.02‰ to +0.11‰) of the MORBs with 87Sr/86Sr > 0.7030. Therefore, one should rigorously choose MORB samples to estimate the average DMM. Based on all Ba isotope data for the D-MORBs, our study provides a more reliable estimate for the average δ138/134Ba of the DMM (+0.05‰±0.05‰, 2SD, n = 16), which leads to new constraints on the nature of the E-MORB source.
5.3
The origin of the enriched source of E-MORB
There are long-term debates on how the mantle source of E-MORB was formed, i.e., melting of the enriched lower mantle [30, 31] , and recycling crustal materials into the mantle (e.g., crustal rocks, metamorphic fluids, and sediments) [3, 13, 15, 32]. Because Ba isotope compositions are unlikely to be modified during mantle melting and fractional crystallization from basaltic melt, Ba isotope variation in MORBs should reflect the introduction of recycled, Ba-rich crustal components into the mantle.
Based on the early estimates of average DMM data, Nielsen et al.[7] proposed that Ba isotope variation in MORBs is derived from pervasive mixing between isotopically heavy DMM and light EMM reservoirs. However, different insights into the formation of the EMM may be gained if the new DMM estimation is applied. The E-MORBs and BABBs show a slightly larger range of δ138/134Ba (−0.06‰ to +0.11‰) relative to D-MORBs. The δ138/134Ba does not show a simple trend with increasing (La/Sm)N, 87Sr/86Sr, and Ba/Th (Fig. 1), which is inconsistent with previously observed hyperbolic mixing relationships[7]. Nielsen et al.[7] attributed the isotopically light Ba signature of EMM (≈ +0.03‰) to the pervasive addition of sedimentary materials to DMM (≈ +0.14‰), as they proposed that sediments have a narrow range of δ138/134Ba (0.01‰±0.04‰, 2SD) and are similar to the E-MORB (+0.03‰ ± 0.02‰, 2SD). However, if the new estimation is applied, the δ138/134Ba values of such sediments will be similar to the DMM value (+0.01‰ vs. +0.05‰). Considering the observations of the MORBs and sediments together, it is unreasonable to infer that the Ba isotope signatures of E-MORBs are derived from pervasive recycled sediments in their mantle source. Additionally, the δ138/134Ba observed in this study did not show a simple trend, but became more scattered with increasing 87Sr/86Sr and Ba/Th (Figs. 1b and 1c).
Significant variations in δ138/134Ba for crustal materials, such as marine sediments (−0.11‰ to +0.10‰)[7-9], AOC (−0.09‰ to +0.33‰)[7], and diamictites[6] have been reported, which are much wider than the range used by Nielsen et al.[7]. Although the AOC does not have much higher Ba content and 87Sr/86Sr relative to the sediment, a study of the Tonga arc revealed that the addition of AOC-derived fluid to the mantle wedge is the most probable source of the high δ138/134Ba in the arc lavas (up to +0.16‰)[11]. In addition, Hao et al.[33] identified the contribution of bulk AOC, aside from the sediment components, to the formation of the Fushui mafic rocks using Ba isotopes. Therefore, we suggest that subduction of sediments and AOC can account for the variation of δ138/134Ba in the E-MORBs (−0.06‰ to +0.11‰). However, Ba isotopes can be considerably fractionated during subduction-zone fluid processes, resulting in a metamorphic fluid with much higher δ138/134Ba than the subducted materials[34]. Therefore, evaluating the contributions of different crustal components to the source of mantle-derived igneous rock requires careful consideration.
6.
Conclusions
This study presents the Ba isotope compositions of a wide range of MORBs and BABBs from different localities, with variable δ138/134Ba ranging from −0.06‰ to +0.11‰. Based on the D-MORB samples, we re-estimated the average δ138/134Ba of the DMM as +0.05‰±0.05‰ (2SD, n = 16). Because of the similar δ138/134Ba values of sediments and the DMM, it is unreasonable to conclude that Ba isotope signatures of E-MORBs are attributed to pervasive sediment recycling in the upper mantle. Instead, the Ba isotope data of the E-MORBs indicate that the EMM could reflect the incorporation of subducted AOC and/or sediments, depending on the Ba isotope composition and other geochemical information of the local mantle.
We are grateful to Charles H. Langmuir and Michael Perfit for providing the MORB samples. This work was financially supported by the National Natural Science Foundation of China (42073007, 41803003) and the Fundamental Research Funds for the Central Universities (WK2080000149). We also appreciate the reviewers’ constructive comments.
Conflict of interest
The authors declare that they have no conflict of interest.
Biographies
Xiaoyun Nan is currently working as an associate research fellow at the University of Science and Technology of China (USTC). She received her PhD in Geology from the USTC in 2017. Her research focuses on metal–stable isotope geochemistry.
Conflict of Interest
The authors declare that they have no conflict of interest.
Homogeneous flow models and separated flow models with and without metastable flow are compared in detail.
The optimal hypothesis models for different R600a mass flow rate ranges are investigated.
The optimal combinations of friction factor and two-phase viscosity correlations for each model are recommended.
Ali M S, Anwar Z, Mujtaba M A, et al. Two-phase frictional pressure drop with pure refrigerants in vertical mini/micro-channels. Case Studies in Thermal Engineering,2021, 23: 100824. DOI: https://doi.org/10.1016/j.csite.2020.100824
[2]
Zhang N C, Li B, Feng L H, et al. Research on the thermophysical properties and cycle performances of R1234yf/R290 and R1234yf/R600a. International Journal of Thermophysics,2021, 42: 123. DOI: 10.1007/s10765-021-02875-0
[3]
Zhang N, Hu P, Chen L X, et al. Measurements of critical properties of the binary mixture of 1,1,1-trifluoroethane (HFC-143a)+ trans-1,3,3,3-Tetrafluoropropene (HFO-1234ze (E)). Journal of Chemical & Engineering Data,2021, 66: 2717–2722. DOI: https://doi.org/10.1021/acs.jced.1c00065
[4]
Bagherzadeh S A, D'Orazio A, Karimipour A, et al. A novel sensitivity analysis model of EANN for F-MWCNTs-Fe3O4/EG nanofluid thermal conductivity: Outputs predicted analytically instead of numerically to more accuracy and less costs. Physica A: Statistical Mechanics and its Applications,2019, 521: 406–415. DOI: 10.1016/j.physa.2019.01.048
[5]
Giwa S O, Sharifpur M, Goodarzi M, et al. Influence of base fluid, temperature, and concentration on the thermophysical properties of hybrid nanofluids of alumina–ferrofluid: Experimental data, modeling through enhanced ANN, ANFIS, and curve fitting. Journal of Thermal Analysis and Calorimetry,2021, 143: 4149–4167. DOI: 10.1007/s10973-020-09372-w
[6]
Sánchez D, Cabello R, Llopis R, et al. Energy performance evaluation of R1234yf, R1234ze(E), R600a, R290 and R152a as low-GWP R134a alternatives. International Journal of Refrigeration,2017, 74: 269–282. DOI: https://doi.org/10.1016/j.ijrefrig.2016.09.020
[7]
Hwang S, Jeong J H. The effects of the parameters of a refrigeration system working with R600a on the non-equilibrium subcooled two-phase flow of the refrigerant. International Journal of Refrigeration,2020, 118: 462–469. DOI: 10.1016/j.ijrefrig.2020.06.026
[8]
Bahmani M H, Sheikhzadeh G, Zarringhalam M, et al. Investigation of turbulent heat transfer and nanofluid flow in a double pipe heat exchanger. Advanced Powder Technology,2018, 29 (2): 273–282. DOI: 10.1016/j.apt.2017.11.013
[9]
Sarafraz M M, Tian Z, Tlili I, et al. Thermal evaluation of a heat pipe working with n-pentane-acetone and n-pentane-methanol binary mixtures. Journal of Thermal Analysis and Calorimetry,2020, 139: 2435–2445. DOI: 10.1007/s10973-019-08414-2
[10]
Heimel M, Lang W, Berger E, et al. A homogeneous capillary tube model - comprehensive parameter studies using isobutane as refrigerant. In: International Refrigeration and Air Conditioning Conference. West Lafayette: Purdue University, 2012: Paper 1233.
[11]
de Lara J F, Melo C, Boeng J, et al. Experimental analysis of HFC-134a expansion through small-bore adiabatic capillary tubes. International Journal of Refrigeration,2020, 112: 37–43. DOI: 10.1016/j.ijrefrig.2019.12.015
[12]
Ardhapurkar P M, Sridharan A, Atrey M D. Investigation of pressure drop in capillary tube for mixed refrigerant Joule-Thomson cryocooler. AIP Conference Proceedings,2014, 1573: 155. DOI: https://doi.org/10.1063/1.4860696
[13]
Kruthiventi S S H, Venkatarathnam G. Studies on capillary tube expansion device used in J-T refrigerators operating with nitrogen-hydrocarbon mixtures. Cryogenics,2017, 87: 76–84. DOI: 10.1016/j.cryogenics.2017.09.002
[14]
Parmar D, Atrey M D. Experimental and numerical investigation on the flow of mixed refrigerants through capillary tubes at cryogenic temperatures. Applied Thermal Engineering,2020, 175: 115339. DOI: 10.1016/j.applthermaleng.2020.115339
[15]
Rocha T T M, de Paula C H, Cangussu V M, et al. Effect of surface roughness on the mass flow rate predictions for adiabatic capillary tubes. International Journal of Refrigeration,2020, 118: 269–278. DOI: 10.1016/j.ijrefrig.2020.05.020
[16]
Khan M K, Kumar R, Sahoo P K. Flow characteristics of refrigerants flowing through capillary tubes-A review. Applied Thermal Engineering,2009, 29: 1426–1439. DOI: 10.1016/j.applthermaleng.2008.08.020
[17]
García-Valladares O. Numerical simulation and experimental validation of coiled adiabatic capillary tubes. Applied Thermal Engineering,2007, 27: 1062–1071. DOI: 10.1016/j.applthermaleng.2006.07.034
[18]
Schenk M, Oellrich L R. Experimental investigation of the refrigerant flow of isobutane (R600a) through adiabatic capillary tubes. International Journal of Refrigeration,2014, 38: 275–280. DOI: 10.1016/j.ijrefrig.2013.08.024
[19]
Chingulpitak S, Wongwises S. Two-phase flow model of refrigerants flowing through helically coiled capillary tubes. Applied Thermal Engineering,2010, 30: 1927–1936. DOI: 10.1016/j.applthermaleng.2010.04.026
[20]
Vinš V, Hrubý J, Vacek V. Numerical simulation of gas-contaminated refrigerant two-phase flow through adiabatic capillary tubes. International Journal of Heat and Mass Transfer,2010, 53: 5430–5439. DOI: 10.1016/j.ijheatmasstransfer.2010.07.013
[21]
Rasti M, Jeong J H. A generalized continuous empirical correlation for the refrigerant mass flow rate through adiabatic straight and helically coiled capillary tubes. Applied Thermal Engineering,2018, 143: 450–460. DOI: 10.1016/j.applthermaleng.2018.07.124
[22]
Hermes C J L, Melo C, Knabben F T. Algebraic solution of capillary tube flows Part I: Adiabatic capillary tubes. Applied Thermal Engineering,2010, 30: 449–457. DOI: 10.1016/j.applthermaleng.2009.10.005
[23]
Dubba S K, Kumar R. Flow of refrigerants through capillary tubes: A state-of-the-art. Experimental Thermal and Fluid Science,2017, 81: 370–381. DOI: 10.1016/j.expthermflusci.2016.09.012
[24]
Jadhav P, Agrawal N. A comparative study in the straight and a spiral adiabatic capillary tube. International Journal of Ambient Energy,2019, 40: 693–698. DOI: 10.1080/01430750.2017.1422146
[25]
Alok P, Sahu D. Numerical simulation of capillary tube for selected refrigerants using homogeneous equilibrium model. International Journal of Air-Conditioning and Refrigeration,2019, 27: 1950001. DOI: 10.1142/S2010132519500019
[26]
Jadhav P, Agrawal N. A comparative study of flow characteristics of adiabatic spiral and helical capillary tube in a CO2 transcritical system. International Journal of Ambient Energy,2021: 1–8. DOI: 10.1080/01430750.2021.1913645
[27]
Zareh M, Heidari M G, Javidmand P. Numerical simulation and experimental comparison of the R12, R22 and R134a flow inside straight and coiled helical capillary tubes. Journal of Mechanical Science and Technology,2016, 30: 1421–1430. DOI: 10.1007/s12206-016-0250-2
[28]
Wang J, Cao F, Wang Z Z, et al. Numerical simulation of coiled adiabatic capillary tubes in CO2 transcritical systems with separated flow model including metastable flow. International Journal of Refrigeration,2012, 35 (8): 2188–2198. DOI: 10.1016/j.ijrefrig.2012.07.012
[29]
Agrawal N, Bhattacharyya S. Homogeneous versus separated two phase flow models: Adiabatic capillary tube flow in a transcritical CO2 heat pump. International Journal of Thermal Sciences,2008, 47 (11): 1555–1562. DOI: 10.1016/j.ijthermalsci.2007.12.008
[30]
Furlong T W, Schmidt D P. A comparison of homogenous and separated flow assumptions for adiabatic capillary flow. Applied Thermal Engineering,2012, 48: 186–193. DOI: 10.1016/j.applthermaleng.2012.05.007
[31]
Lorbek L, Kuhelj A, Dular M, et al. Two-phase flow patterns in adiabatic refrigerant flow through capillary tubes. International Journal of Refrigeration,2020, 115: 107–116. DOI: 10.1016/j.ijrefrig.2020.02.030
[32]
Melo C, Ferreira R T S, Neto C B, et al. An experimental analysis of adiabatic capillary tubes. Applied Thermal Engineering,1999, 19: 669–684. DOI: 10.1016/S1359-4311(98)00062-3
[33]
Collier J G, Thome J R. Convective Boiling and Condensation. 3rd ed. New York: Clarendon Press, 1994.
[34]
Chen Z H, Li R Y, Lin S, et al. A correlation for metastable flow of R-12 through capillary tubes. ASHRAE Transactions,1990, 96: 550–554.
[35]
Feburie V, Giot M, Granger S, et al. A model for choked flow through cracks with inlet subcooling. International Journal of Multiphase Flow,1993, 19: 541–562. DOI: 10.1016/0301-9322(93)90087-B
[36]
Premoli A, Francesco D, Prina A. An empirical correlation for evaluating two-phase mixture density under adiabatic conditions. In: European Two-Phase Flow Group Meeting, Milan, Italy, 1970.
[37]
Chisholm D. Pressure gradients due to friction during the flow of evaporating two-phase mixtures in smooth tubes and channels. International Journal of Heat and Mass Transfer,1973, 16: 347–358. DOI: 10.1016/0017-9310(73)90063-X
[38]
Deodhar S D, Kothadia H B, Iyer K N, et al. Experimental and numerical studies of choked flow through adiabatic and diabatic capillary tubes. Applied Thermal Engineering,2015, 90: 879–894. DOI: 10.1016/j.applthermaleng.2015.07.073
Zhang Y F, Zhou G B, Xie H, et al. An assessment of friction factor and viscosity correlations for model prediction of refrigerant flow in capillary tubes. International Journal of Energy Research,2005, 29 (3): 233–248. DOI: 10.1002/er.1050
[41]
Ahmadi M H, Mohseni-Gharyehsafa B, Ghazvini M, et al. Comparing various machine learning approaches in modeling the dynamic viscosity of CuO/water nanofluid. Journal of Thermal Analysis and Calorimetry,2020, 139: 2585–2599. DOI: 10.1007/s10973-019-08762-z
[42]
Zhou G B, Zhang Y F. Numerical and experimental investigations on the performance of coiled adiabatic capillary tubes. Applied Thermal Engineering,2006, 26: 1106–1114. DOI: 10.1016/j.applthermaleng.2005.11.003
[43]
Bansal P K, Wang G. Numerical analysis of choked refrigerant flow in adiabatic capillary tubes. Applied Thermal Engineering,2004, 24: 851–863. DOI: 10.1016/j.applthermaleng.2003.10.010
[44]
Lemmon E, Huber M, McLinden M. NIST Standard Reference Database 23: NIST thermodynamic and transport properties of refrigerants and refrigerant mixtures-REFPROP, version 9.1. Gaithersburg, MD: National Institute of Standards and Technology, 2013.
Kamaruddin, R., Sulaiman, S., Khalid, A. et al. Computational Fluid Dynamics: Flow Analysis on the Effect of Different Jet Orifice Angle Multi Circular Jet for Fuel and Air Mixing. CFD Letters, 2024, 16(1): 121-137.
DOI:10.37934/cfdl.16.1.121137
Figure
1.
Schematic diagram of adiabatic capillary tube. (a) Case with metastable flow. (b) Case without metastable flow.
Figure
2.
Control volume of adiabatic capillary tube.
Figure
3.
Flowcharts of mass flow rate algorithm in case without metastable flow. (a) Algorithm for data of Melo et al.[32]. (b) Algorithm for data of Schenk and Oellrich[18].
4.
Flowcharts of mass flow rate algorithm in case with metastable flow. (a) Algorithm for data of Melo et al.[32]. (b) Algorithm for data of Schenk and Oellrich[18].
Figure
5.
Comparison of different viscosity correlations.
Figure
6.
Comparison of different friction factor correlations.
Figure
7.
Comparison of mass flow rates predicted by homogenous flow model and separated flow model to experimental data. (a) Case without metastable flow. (b) Case with metastable flow.
Figure
8.
Comparison of pressure distributions along tube.
Figure
9.
Comparison of temperature distributions along tube.
Figure
10.
Comparison of vapor phase mass fraction distributions along two-phase region.
Figure
11.
Comparison of velocity distributions along two-phase region.
References
[1]
Ali M S, Anwar Z, Mujtaba M A, et al. Two-phase frictional pressure drop with pure refrigerants in vertical mini/micro-channels. Case Studies in Thermal Engineering,2021, 23: 100824. DOI: https://doi.org/10.1016/j.csite.2020.100824
[2]
Zhang N C, Li B, Feng L H, et al. Research on the thermophysical properties and cycle performances of R1234yf/R290 and R1234yf/R600a. International Journal of Thermophysics,2021, 42: 123. DOI: 10.1007/s10765-021-02875-0
[3]
Zhang N, Hu P, Chen L X, et al. Measurements of critical properties of the binary mixture of 1,1,1-trifluoroethane (HFC-143a)+ trans-1,3,3,3-Tetrafluoropropene (HFO-1234ze (E)). Journal of Chemical & Engineering Data,2021, 66: 2717–2722. DOI: https://doi.org/10.1021/acs.jced.1c00065
[4]
Bagherzadeh S A, D'Orazio A, Karimipour A, et al. A novel sensitivity analysis model of EANN for F-MWCNTs-Fe3O4/EG nanofluid thermal conductivity: Outputs predicted analytically instead of numerically to more accuracy and less costs. Physica A: Statistical Mechanics and its Applications,2019, 521: 406–415. DOI: 10.1016/j.physa.2019.01.048
[5]
Giwa S O, Sharifpur M, Goodarzi M, et al. Influence of base fluid, temperature, and concentration on the thermophysical properties of hybrid nanofluids of alumina–ferrofluid: Experimental data, modeling through enhanced ANN, ANFIS, and curve fitting. Journal of Thermal Analysis and Calorimetry,2021, 143: 4149–4167. DOI: 10.1007/s10973-020-09372-w
[6]
Sánchez D, Cabello R, Llopis R, et al. Energy performance evaluation of R1234yf, R1234ze(E), R600a, R290 and R152a as low-GWP R134a alternatives. International Journal of Refrigeration,2017, 74: 269–282. DOI: https://doi.org/10.1016/j.ijrefrig.2016.09.020
[7]
Hwang S, Jeong J H. The effects of the parameters of a refrigeration system working with R600a on the non-equilibrium subcooled two-phase flow of the refrigerant. International Journal of Refrigeration,2020, 118: 462–469. DOI: 10.1016/j.ijrefrig.2020.06.026
[8]
Bahmani M H, Sheikhzadeh G, Zarringhalam M, et al. Investigation of turbulent heat transfer and nanofluid flow in a double pipe heat exchanger. Advanced Powder Technology,2018, 29 (2): 273–282. DOI: 10.1016/j.apt.2017.11.013
[9]
Sarafraz M M, Tian Z, Tlili I, et al. Thermal evaluation of a heat pipe working with n-pentane-acetone and n-pentane-methanol binary mixtures. Journal of Thermal Analysis and Calorimetry,2020, 139: 2435–2445. DOI: 10.1007/s10973-019-08414-2
[10]
Heimel M, Lang W, Berger E, et al. A homogeneous capillary tube model - comprehensive parameter studies using isobutane as refrigerant. In: International Refrigeration and Air Conditioning Conference. West Lafayette: Purdue University, 2012: Paper 1233.
[11]
de Lara J F, Melo C, Boeng J, et al. Experimental analysis of HFC-134a expansion through small-bore adiabatic capillary tubes. International Journal of Refrigeration,2020, 112: 37–43. DOI: 10.1016/j.ijrefrig.2019.12.015
[12]
Ardhapurkar P M, Sridharan A, Atrey M D. Investigation of pressure drop in capillary tube for mixed refrigerant Joule-Thomson cryocooler. AIP Conference Proceedings,2014, 1573: 155. DOI: https://doi.org/10.1063/1.4860696
[13]
Kruthiventi S S H, Venkatarathnam G. Studies on capillary tube expansion device used in J-T refrigerators operating with nitrogen-hydrocarbon mixtures. Cryogenics,2017, 87: 76–84. DOI: 10.1016/j.cryogenics.2017.09.002
[14]
Parmar D, Atrey M D. Experimental and numerical investigation on the flow of mixed refrigerants through capillary tubes at cryogenic temperatures. Applied Thermal Engineering,2020, 175: 115339. DOI: 10.1016/j.applthermaleng.2020.115339
[15]
Rocha T T M, de Paula C H, Cangussu V M, et al. Effect of surface roughness on the mass flow rate predictions for adiabatic capillary tubes. International Journal of Refrigeration,2020, 118: 269–278. DOI: 10.1016/j.ijrefrig.2020.05.020
[16]
Khan M K, Kumar R, Sahoo P K. Flow characteristics of refrigerants flowing through capillary tubes-A review. Applied Thermal Engineering,2009, 29: 1426–1439. DOI: 10.1016/j.applthermaleng.2008.08.020
[17]
García-Valladares O. Numerical simulation and experimental validation of coiled adiabatic capillary tubes. Applied Thermal Engineering,2007, 27: 1062–1071. DOI: 10.1016/j.applthermaleng.2006.07.034
[18]
Schenk M, Oellrich L R. Experimental investigation of the refrigerant flow of isobutane (R600a) through adiabatic capillary tubes. International Journal of Refrigeration,2014, 38: 275–280. DOI: 10.1016/j.ijrefrig.2013.08.024
[19]
Chingulpitak S, Wongwises S. Two-phase flow model of refrigerants flowing through helically coiled capillary tubes. Applied Thermal Engineering,2010, 30: 1927–1936. DOI: 10.1016/j.applthermaleng.2010.04.026
[20]
Vinš V, Hrubý J, Vacek V. Numerical simulation of gas-contaminated refrigerant two-phase flow through adiabatic capillary tubes. International Journal of Heat and Mass Transfer,2010, 53: 5430–5439. DOI: 10.1016/j.ijheatmasstransfer.2010.07.013
[21]
Rasti M, Jeong J H. A generalized continuous empirical correlation for the refrigerant mass flow rate through adiabatic straight and helically coiled capillary tubes. Applied Thermal Engineering,2018, 143: 450–460. DOI: 10.1016/j.applthermaleng.2018.07.124
[22]
Hermes C J L, Melo C, Knabben F T. Algebraic solution of capillary tube flows Part I: Adiabatic capillary tubes. Applied Thermal Engineering,2010, 30: 449–457. DOI: 10.1016/j.applthermaleng.2009.10.005
[23]
Dubba S K, Kumar R. Flow of refrigerants through capillary tubes: A state-of-the-art. Experimental Thermal and Fluid Science,2017, 81: 370–381. DOI: 10.1016/j.expthermflusci.2016.09.012
[24]
Jadhav P, Agrawal N. A comparative study in the straight and a spiral adiabatic capillary tube. International Journal of Ambient Energy,2019, 40: 693–698. DOI: 10.1080/01430750.2017.1422146
[25]
Alok P, Sahu D. Numerical simulation of capillary tube for selected refrigerants using homogeneous equilibrium model. International Journal of Air-Conditioning and Refrigeration,2019, 27: 1950001. DOI: 10.1142/S2010132519500019
[26]
Jadhav P, Agrawal N. A comparative study of flow characteristics of adiabatic spiral and helical capillary tube in a CO2 transcritical system. International Journal of Ambient Energy,2021: 1–8. DOI: 10.1080/01430750.2021.1913645
[27]
Zareh M, Heidari M G, Javidmand P. Numerical simulation and experimental comparison of the R12, R22 and R134a flow inside straight and coiled helical capillary tubes. Journal of Mechanical Science and Technology,2016, 30: 1421–1430. DOI: 10.1007/s12206-016-0250-2
[28]
Wang J, Cao F, Wang Z Z, et al. Numerical simulation of coiled adiabatic capillary tubes in CO2 transcritical systems with separated flow model including metastable flow. International Journal of Refrigeration,2012, 35 (8): 2188–2198. DOI: 10.1016/j.ijrefrig.2012.07.012
[29]
Agrawal N, Bhattacharyya S. Homogeneous versus separated two phase flow models: Adiabatic capillary tube flow in a transcritical CO2 heat pump. International Journal of Thermal Sciences,2008, 47 (11): 1555–1562. DOI: 10.1016/j.ijthermalsci.2007.12.008
[30]
Furlong T W, Schmidt D P. A comparison of homogenous and separated flow assumptions for adiabatic capillary flow. Applied Thermal Engineering,2012, 48: 186–193. DOI: 10.1016/j.applthermaleng.2012.05.007
[31]
Lorbek L, Kuhelj A, Dular M, et al. Two-phase flow patterns in adiabatic refrigerant flow through capillary tubes. International Journal of Refrigeration,2020, 115: 107–116. DOI: 10.1016/j.ijrefrig.2020.02.030
[32]
Melo C, Ferreira R T S, Neto C B, et al. An experimental analysis of adiabatic capillary tubes. Applied Thermal Engineering,1999, 19: 669–684. DOI: 10.1016/S1359-4311(98)00062-3
[33]
Collier J G, Thome J R. Convective Boiling and Condensation. 3rd ed. New York: Clarendon Press, 1994.
[34]
Chen Z H, Li R Y, Lin S, et al. A correlation for metastable flow of R-12 through capillary tubes. ASHRAE Transactions,1990, 96: 550–554.
[35]
Feburie V, Giot M, Granger S, et al. A model for choked flow through cracks with inlet subcooling. International Journal of Multiphase Flow,1993, 19: 541–562. DOI: 10.1016/0301-9322(93)90087-B
[36]
Premoli A, Francesco D, Prina A. An empirical correlation for evaluating two-phase mixture density under adiabatic conditions. In: European Two-Phase Flow Group Meeting, Milan, Italy, 1970.
[37]
Chisholm D. Pressure gradients due to friction during the flow of evaporating two-phase mixtures in smooth tubes and channels. International Journal of Heat and Mass Transfer,1973, 16: 347–358. DOI: 10.1016/0017-9310(73)90063-X
[38]
Deodhar S D, Kothadia H B, Iyer K N, et al. Experimental and numerical studies of choked flow through adiabatic and diabatic capillary tubes. Applied Thermal Engineering,2015, 90: 879–894. DOI: 10.1016/j.applthermaleng.2015.07.073
Zhang Y F, Zhou G B, Xie H, et al. An assessment of friction factor and viscosity correlations for model prediction of refrigerant flow in capillary tubes. International Journal of Energy Research,2005, 29 (3): 233–248. DOI: 10.1002/er.1050
[41]
Ahmadi M H, Mohseni-Gharyehsafa B, Ghazvini M, et al. Comparing various machine learning approaches in modeling the dynamic viscosity of CuO/water nanofluid. Journal of Thermal Analysis and Calorimetry,2020, 139: 2585–2599. DOI: 10.1007/s10973-019-08762-z
[42]
Zhou G B, Zhang Y F. Numerical and experimental investigations on the performance of coiled adiabatic capillary tubes. Applied Thermal Engineering,2006, 26: 1106–1114. DOI: 10.1016/j.applthermaleng.2005.11.003
[43]
Bansal P K, Wang G. Numerical analysis of choked refrigerant flow in adiabatic capillary tubes. Applied Thermal Engineering,2004, 24: 851–863. DOI: 10.1016/j.applthermaleng.2003.10.010
[44]
Lemmon E, Huber M, McLinden M. NIST Standard Reference Database 23: NIST thermodynamic and transport properties of refrigerants and refrigerant mixtures-REFPROP, version 9.1. Gaithersburg, MD: National Institute of Standards and Technology, 2013.
Kamaruddin, R., Sulaiman, S., Khalid, A. et al. Computational Fluid Dynamics: Flow Analysis on the Effect of Different Jet Orifice Angle Multi Circular Jet for Fuel and Air Mixing. CFD Letters, 2024, 16(1): 121-137.
DOI:10.37934/cfdl.16.1.121137
*87Sr/86Sr values of the samples with asterisks were measured in this study. a La/Sm was normalized to CI-chondrite[3]. b 2SD = two times the standard deviation of the population of n repeat measurements of a sample solution. cn represents the number of repeated measurements of the same solution using MC-ICP-MS. d Data were calculated from the measured δ137/134Ba values, as described in the text. e Replicate = repeat sample dissolution, purified column chemistry, and MC-ICP-MS analysis. PetDB= the Petrological Database website (http://www.earthchem.org/petdb).