Physical and mathematical problems of 1D modelling of transonic two-phase flow in a convergent-divergent nozzle
DOI:
https://doi.org/10.24425/ather.2024.152017Abstract
This work addresses the challenges associated with one-dimensional modelling of steady transonic two-phase flows illustrated through simulations of selected flow cases of the famous Moby Dick experiment. It primarily concentrates on the method for a fast determination of the transonic trajectory. The second proposed approach allows for determining trajectories describing the transonic flow with a normal shock wave. The first method is successfully verified by comparing its results with simulation results obtained from the (widely known and thoroughly verified) Wavefront Algorithm for High-speed Aerodynamics code, utilizing the Delayed Equilibrium Model. The first method mentioned is the author's proposition that is competitive to conventional (time-expensive) approaches such as the Newton Critical Point or achieving a steady flow description by asymptotical convergence of the time-dependent model's solutions, and it is a completely new consistent solution method. The second proposed method is an adaptation of the Rankine–Hugoniot jump conditions to a two-phase flow described by the Delayed Equilibrium Model. In the case of this method, the presented here analysis and results serve only as a proof of concept. Similar methods have been described before, but the results presented in this article, obtained with the Delayed Equilibrium Model, are unique. Also, a limited but coherent model of thermodynamic properties of a superheated liquid is presented and physically justified. This model was formulated earlier and is often used, but its comprehensive derivation has not been presented before.References
[1] Bilicki, Z., & Kestin, J. (1990). Physical aspects of the relaxation model in two-phase flow. Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences, 428(1875), 379–397. https://doi.org/10.1098/rspa.1990.0040
[2] Bilicki, Z., Dafermos, C., Kestin, J., Majda, G., & Zeng, D. L. (1987). Trajectories and singular points in steady-state models of two-phase flows. International Journal of Multiphase Flow, 13(4), 511–533. https://doi.org/10.1016/0301-9322(87)90019-X
[3] De Sterck, H. (2007). Critical point analysis of transonic flow profiles with heat conduction. SIAM Journal on Applied Dynamical Systems, 6(3), 645–662. https://doi.org/10.1137/060677458
[4] Yee, H. C. (1994). Basins of attraction and the time-dependent approach to obtaining steady-state numerical solutions. In Numerical methods for fluid dynamics (pp. 135–162). Oxford Academic. https://doi.org/10.1093/oso/9780198536963.003.0009
[5] Angielczyk, W., Śmierciew, K., & Butrymowicz, D. (2019). Application of a fast transonic trajectory determination approach in 1-D modelling of steady-state two-phase carbon dioxide flow. E3S Web of Conferences, 128, 06005. https://doi.org/10.1051/e3sconf/201912806005
[6] Seynhaeve, J. M., & Giot, M. (2004). The WAHA code: A numerical tool for water hammer in two-phase flow—Some simulations of experiments. Proceedings of the 117th Session of the Scientific and Technical Committee of the French Hydrotechnical Society, Lyon, France.
[7] Banaszkiewicz, M., & Kardaś, D. (1997). Numerical calculations of the Moby Dick experiment by means of unsteady relaxation models. Journal of Theoretical and Applied Mechanics, 35(2), 211–232.
[8] Bartosiewicz, Y., & Seynhaeve, J. M. (2013). Delayed equilibrium model (DEM) of flashing choked flows relevant to LOCA. Multiphase Science and Technology, 25(2–4), 117–131. https://doi.org/10.1615/MultScienTechn.v25.i2-4.50
[9] Bartosiewicz, Y., Giot, M., & Seynhaeve, J. M. (2010). Revisiting modeling techniques and validation experiments for two-phase choked flows relevant to LOCA. Proceedings of NUTHOS-8, Shanghai, China.
[10] Emanuel, G. (1964). A general method for numerical integration through a saddle-point singularity with application to one-dimensional nonequilibrium nozzle flow (Technical Report No. AEDC-TDR-64-29). Stanford University.
[11] De Sterck, H., Rostrup, S., & Tian, F. (2009). A fast and accurate algorithm for computing radial transonic flows. Journal of Computational and Applied Mathematics, 223(2), 916–928. https://doi.org/10.1016/j.cam.2008.03.019
[12] Cheng, K., Meng, T., Zhao, F., & Tan, S. (2019). Development and validation of a thermal hydraulic transient analysis code for offshore floating nuclear reactor based on RELAP5/SCDAPSIM/MOD3.4. Annals of Nuclear Energy, 127, 215–226. https://doi.org/10.1016/j.anucene.2018.12.004
[13] Ohnuki, A., Akimoto, H., & Murao, Y. (1993). Assessment of TRAC-PF1/MOD1 code for thermal-hydraulic behavior in pressure vessel during reflood in SCTF test (JAERI-M-93-139). Japan Atomic Energy Research Institute.
[14] Dor, I., Geffraye, G., Lavialle, G., & Mieusset, T. (2004). Recent improvements of physical models in the CATHARE code and their validation. Proceedings of the 12th International Conference on Nuclear Engineering (pp. 681–691). https://doi.org/10.1115/ICONE12-49408
[15] Sánchez, V. H., Thieme, M., & Tietsch, W. (2012). Validation and application of the thermal hydraulic system code TRACE for analysis of BWR transients. Science and Technology of Nuclear Installations, 2012, 247482. https://doi.org/10.1155/2012/247482
[16] Gale, J., Tiselj, I., & Horvat, A. (2008). Two-fluid model of the WAHA code for simulations of water hammer transients. Multiphase Science and Technology, 20(3–4), 291–322. https://doi.org/10.1615/MultScienTechn.v20.i3-4.40
[17] Va, B., Lee, W.-T., & Roe, P. (1991). Characteristic time-stepping or local preconditioning of the Euler equations. Proceedings of the 10th Computational Fluid Dynamics Conference, Honolulu, USA. https://doi.org/10.2514/6.1991-1552
[18] Angielczyk, W., Seynhaeve, J. M., Gagan, J., Bartosiewicz, Y., & Butrymowicz, D. (2019). Prediction of critical mass rate of flashing carbon dioxide flow in convergent-divergent nozzle. Chemical Engineering and Processing: Process Intensification, 143, 107599. https://doi.org/10.1016/j.cep.2019.107599
[19] Angielczyk, W., Bartosiewicz, Y., & Butrymowicz, D. (2020). Development of delayed equilibrium model for CO₂ convergent-divergent nozzle transonic flashing flow. International Journal of Multiphase Flow, 131, 103351. https://doi.org/10.1016/j.ijmultiphaseflow.2020.103351
[20] Angielczyk, W. (2024). A review of the relaxation models for phase transition flows centered on the topological aspects of the nonequilibrium mass transfer modelling. Acta Mechanica et Automatica, 18(3), 526–535. https://doi.org/10.2478/ama-2024-0056
[21] Nakagawa, M., Berana, M. S., & Kishine, A. (2009). Supersonic two-phase flow of CO₂ through converging–diverging nozzles for the ejector refrigeration cycle. International Journal of Refrigeration, 32(6), 1195–1202. https://doi.org/10.1016/j.ijrefrig.2009.01.015
[22] Gale, J., & Tiselj, I. (2004). Simulation of the critical flashing flow with the transient 1D two-fluid model. ASME Heat Transfer/Fluids Engineering Summer Conference. https://doi.org/10.1115/ht-fed2004-56315
[23] Ebrahimzadeh Azghadi, F. E., Mahpeykar, M. R., Pasandideh Fard, M., & Lakzian, E. (2022). Numerical simulation of nucleating flow and shock capturing in steam turbines by a simple low-dissipation upwind scheme using an Eulerian–Lagrangian model. Computers & Fluids, 249, 105699. https://doi.org/10.1016/j.compfluid.2022.105699
[24] Gaballa, H., Jafari, S., Di-Lella, A., Habchi, C., & de Hemptinne, J. C. (2022). Thermodynamics analysis of CO₂ condensation in supersonic flows for the potential of clean offshore natural gas processing. Applied Energy, 320, 118523. https://doi.org/10.1016/j.apenergy.2022.118523
[25] Chen, Z., & Liu, G. (2024). Analysis of the internal flow features of a CO₂ transonic nozzle and optimization of the nozzle shape profile. Applied Thermal Engineering, 243, 121945. https://doi.org/10.1016/j.applthermaleng.2023.121945
[26] Seynhaeve, J.-M., De Crécy, A., & Bartosiewicz, Y. (2015). Uncertainty analysis of delayed equilibrium model (DEM) using the CIRCE methodology. Proceedings of NURETH-16, Chicago, USA.
[27] Attou, A., Bolle, L., & Seynhaeve, J. M. (2000). Experimental study of the critical flashing flow through a relief line: Evidence of the double-choked flow phenomenon. International Journal of Multiphase Flow, 26(6), 921–947. https://doi.org/10.1016/S0301-9322(99)00077-4
[28] Lorenzo, M. D., Lafon, P., Seynhaeve, J.-M., & Bartosiewicz, Y. (2017). Benchmark of delayed equilibrium model (DEM) and classic two-phase critical flow models against experimental data. International Journal of Multiphase Flow, 92, 112–130. https://doi.org/10.1016/j.ijmultiphaseflow.2017.03.004
[29] Skripov, V. P. (1992). Metastable states. Journal of Non-Equilibrium Thermodynamics, 17(3), 193–236. https://doi.org/10.1515/jnet.1992.17.3.193
[30] Angielczyk, W. (2021). Modelling of CO₂ transonic flashing flow through a convergent-divergent nozzle using the delayed equilibrium model (PhD thesis). Bialystok University of Technology.
[31] Bouré, J. A., Fritte, A. A., Giot, M. M., & Réocreux, M. L. (1976). Highlights of two-phase critical flow. International Journal of Multiphase Flow, 3(1), 1–22. https://doi.org/10.1016/0301-9322(76)90030-6
[32] Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. (2002). Numerical recipes in C++: The art of scientific computing (3rd ed.). Cambridge University Press.
[33] Bell, I. H., Wronski, J., Quoilin, S., & Lemort, V. (2014). Pure and pseudo-pure fluid thermophysical property evaluation and the open-source thermophysical property library CoolProp. Industrial & Engineering Chemistry Research, 53(6), 2498–2508. https://doi.org/10.1021/ie4033999
