An accurate implementation of wetting and pressure drop is crucial to correctly reproduce fluid displacement processes in porous media. Although several strategies have been proposed in the literature, a systematic comparison of them is needed to determine the most suitable for practical applications. Here, we carried out numerical simulations to investigate the performance of two widely used wettability schemes in the lattice Boltzmann color gradient model, namely, the geometrical wetting scheme by Leclaire et al. [Phys. Rev. E 95(3), 033306 (2017)](scheme-I) and the modified direction of the color gradient scheme by Akai et al. [Adv. Water Resour. 116, 56–66 (2018)] (scheme-II). We showed that scheme-II was more accurate in simulating static contact angles of a fluid droplet on a solid surface. However, scheme-I was more accurate in simulating a dynamic case of a binary fluid flow in a horizontal capillary tube described by the Washburn equation. Moreover, we investigated the performance of two popular pressure gradient implementation types. Type-I used the so-called Zou–He pressure boundary conditions at the inlet and the outlet of the domain, while type-II used an external body force as a pressure gradient. We showed that the type-I implementation was slightly more accurate in simulating a neutrally wetting fluid in a horizontal capillary tube described by the Washburn equation. We also investigated the differences between the two types of pressure gradient implementation in simulating two fluid displacement processes in a Bentheimer sandstone rock sample: the primary drainage and the imbibition displacement processes.
要正确再现多孔介质中的流体位移过程，准确实现润湿和压降至关重要。虽然文献中已经提出了几种策略，但需要对它们进行系统比较，以确定最适合实际应用的策略。在此，我们进行了数值模拟，研究了晶格玻尔兹曼颜色梯度模型中两种广泛使用的润湿方案的性能，即 Leclaire 等人的几何润湿方案[Phys. Rev. E 95(3), 033306 (2017)]（方案-I）和 Akai 等人的颜色梯度修正方向方案[Adv. Water Resour. 116, 56-66 (2018)]（方案-II）。我们的研究表明，方案-II 在模拟流体液滴在固体表面的静态接触角方面更为精确。然而，在模拟由瓦什伯恩方程描述的二元流体在水平毛细管中流动的动态情况时，方案一更为精确。此外，我们还研究了两种常用压力梯度实施类型的性能。类型 I 在域的入口和出口处使用了所谓的 Zou-He 压力边界条件，而类型 II 则使用了外部体力作为压力梯度。我们发现，在模拟由瓦什伯恩方程描述的水平毛细管中的中性润湿流体时，I 型实施方案的精确度略高。我们还研究了在模拟本特海默砂岩岩石样本中的两种流体位移过程（原生排水和浸润位移过程）时，两种压力梯度实现方式之间的差异。
{"title":"Wetting and pressure gradient performance in a lattice Boltzmann color gradient model","authors":"M. Sedahmed, R. C. V. Coelho","doi":"10.1063/5.0228835","DOIUrl":"https://doi.org/10.1063/5.0228835","url":null,"abstract":"An accurate implementation of wetting and pressure drop is crucial to correctly reproduce fluid displacement processes in porous media. Although several strategies have been proposed in the literature, a systematic comparison of them is needed to determine the most suitable for practical applications. Here, we carried out numerical simulations to investigate the performance of two widely used wettability schemes in the lattice Boltzmann color gradient model, namely, the geometrical wetting scheme by Leclaire et al. [Phys. Rev. E 95(3), 033306 (2017)](scheme-I) and the modified direction of the color gradient scheme by Akai et al. [Adv. Water Resour. 116, 56–66 (2018)] (scheme-II). We showed that scheme-II was more accurate in simulating static contact angles of a fluid droplet on a solid surface. However, scheme-I was more accurate in simulating a dynamic case of a binary fluid flow in a horizontal capillary tube described by the Washburn equation. Moreover, we investigated the performance of two popular pressure gradient implementation types. Type-I used the so-called Zou–He pressure boundary conditions at the inlet and the outlet of the domain, while type-II used an external body force as a pressure gradient. We showed that the type-I implementation was slightly more accurate in simulating a neutrally wetting fluid in a horizontal capillary tube described by the Washburn equation. We also investigated the differences between the two types of pressure gradient implementation in simulating two fluid displacement processes in a Bentheimer sandstone rock sample: the primary drainage and the imbibition displacement processes.","PeriodicalId":20066,"journal":{"name":"Physics of Fluids","volume":"16 1","pages":""},"PeriodicalIF":4.6,"publicationDate":"2024-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142258835","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
A third-order compact multi-resolution weighted essentially non-oscillatory (CMR-WENO) reconstruction method for three-dimensional (3D) hybrid unstructured grids is developed using the Arbitrary Lagrange–Euler framework. The finite volume method is used to discretize the governing equations, and some turbulent and moving boundary problems are simulated. Only one compact center stencil comprising the neighboring cells of each control cell is required to construct the polynomials in the algorithm. As a result, the number of stencils and stencil cells is significantly reduced when compared with the traditional WENO scheme. This simplifies the code and improves the robustness of the algorithm. By ensuring the cell average and first-order derivatives are consistent with that in stencil cells an over-determined system of equations can be used to reconstruct the polynomials. This system can then be solved using the compact least squares method to avoid an ill-conditioned coefficient matrix. Furthermore, a coupled implicit iteration strategy is used to solve for the unknown coefficients, so no extra determination is required for the derivatives of each control cell. The final interpolation function for discontinuities in the flow field is obtained using CMR-WENO to nonlinearly combine polynomials of different orders, which further improves the stability of the algorithm. The CMR-WENO can be implemented on 3D hybrid unstructured grids and can be used to simulate complex problems such as those involving turbulence and moving boundaries. Finally, the algorithm presented here is verified to be third-order accurate and to exhibit good robustness when used on several representative numerical examples.
{"title":"Three-dimensional compact multi-resolution weighted essentially non-oscillatory reconstruction under the Arbitrary Lagrange–Euler framework","authors":"Ningyu Zhan, Rongqian Chen, Yancheng You","doi":"10.1063/5.0226237","DOIUrl":"https://doi.org/10.1063/5.0226237","url":null,"abstract":"A third-order compact multi-resolution weighted essentially non-oscillatory (CMR-WENO) reconstruction method for three-dimensional (3D) hybrid unstructured grids is developed using the Arbitrary Lagrange–Euler framework. The finite volume method is used to discretize the governing equations, and some turbulent and moving boundary problems are simulated. Only one compact center stencil comprising the neighboring cells of each control cell is required to construct the polynomials in the algorithm. As a result, the number of stencils and stencil cells is significantly reduced when compared with the traditional WENO scheme. This simplifies the code and improves the robustness of the algorithm. By ensuring the cell average and first-order derivatives are consistent with that in stencil cells an over-determined system of equations can be used to reconstruct the polynomials. This system can then be solved using the compact least squares method to avoid an ill-conditioned coefficient matrix. Furthermore, a coupled implicit iteration strategy is used to solve for the unknown coefficients, so no extra determination is required for the derivatives of each control cell. The final interpolation function for discontinuities in the flow field is obtained using CMR-WENO to nonlinearly combine polynomials of different orders, which further improves the stability of the algorithm. The CMR-WENO can be implemented on 3D hybrid unstructured grids and can be used to simulate complex problems such as those involving turbulence and moving boundaries. Finally, the algorithm presented here is verified to be third-order accurate and to exhibit good robustness when used on several representative numerical examples.","PeriodicalId":20066,"journal":{"name":"Physics of Fluids","volume":"40 1","pages":""},"PeriodicalIF":4.6,"publicationDate":"2024-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142258723","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The coexistence of density stratification and surface tension alters the dynamic behaviors of two-phase immiscible wake and jet flows in complex interfacial instability modes. Building on the framework established by Schmidt et al. [“Global stability and nonlinear dynamics of wake flows with a two-fluid interface,” J. Fluid Mech. 915, A96 (2021)], we conduct a global stability analysis to investigate the effects of surface tension and density ratio on the interfacial instabilities of two-phase planar wake and jet flows. Surface tension, acting counterintuitively as either a stabilizer or destabilizer, enhances the self-sustainability of varicose and sinuous disturbances at low levels, while high surface tension ultimately leads to the stabilization of interfacial disturbances. Additionally, sinuous disturbances, characterized by higher oscillation amplitudes, predominate in dense wakes or wakes with strong shear strengths, where surface tension serves exclusively as a stabilizing factor. In contrast, the varicose pattern in jets prevails over a wider range of Weber numbers, exhibiting relatively higher linear growth rates compared to their sinuous counterparts, especially in lighter jets. The temporal-spatial analysis further provides a theoretical demonstration of these findings by offering views into the complex interplay of these factors.
{"title":"On the interfacial instabilities of two-phase wake and jet flows with density stratification and surface tension","authors":"Minjiang Gong, A-Man Zhang, Chengwang Xiong","doi":"10.1063/5.0221984","DOIUrl":"https://doi.org/10.1063/5.0221984","url":null,"abstract":"The coexistence of density stratification and surface tension alters the dynamic behaviors of two-phase immiscible wake and jet flows in complex interfacial instability modes. Building on the framework established by Schmidt et al. [“Global stability and nonlinear dynamics of wake flows with a two-fluid interface,” J. Fluid Mech. 915, A96 (2021)], we conduct a global stability analysis to investigate the effects of surface tension and density ratio on the interfacial instabilities of two-phase planar wake and jet flows. Surface tension, acting counterintuitively as either a stabilizer or destabilizer, enhances the self-sustainability of varicose and sinuous disturbances at low levels, while high surface tension ultimately leads to the stabilization of interfacial disturbances. Additionally, sinuous disturbances, characterized by higher oscillation amplitudes, predominate in dense wakes or wakes with strong shear strengths, where surface tension serves exclusively as a stabilizing factor. In contrast, the varicose pattern in jets prevails over a wider range of Weber numbers, exhibiting relatively higher linear growth rates compared to their sinuous counterparts, especially in lighter jets. The temporal-spatial analysis further provides a theoretical demonstration of these findings by offering views into the complex interplay of these factors.","PeriodicalId":20066,"journal":{"name":"Physics of Fluids","volume":"15 1","pages":""},"PeriodicalIF":4.6,"publicationDate":"2024-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142258724","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Confined jets occur in many engineering applications including combustion chambers, jet pumps, and chemical reactors. The effects of axisymmetric confinement on the vortical structures identified in a turbulent jet are investigated using large eddy simulation at a Reynolds number of 30 000 (based on nozzle exit conditions) and expansion ratio (chamber-to-nozzle diameter ratio) of five. The results obtained from the confined jet are compared with those of a free jet under the same nozzle exit flow conditions. A prominent recirculation zone forms between the expanding jet and the confining wall, resulting in early shear layer distortion and a shorter interaction length in the confined jet (0.85 jet diameters) compared to the free jet (1.15 jet diameters). Using the λ2 criterion for vortex identification, two dominant structural modes are identified in the near-exit region of the free jet: ring and helical modes. However, in the confined jet, the helical mode is absent, and the turbulent confined fluid accelerates the breakup of the ring vortices. The interaction of the secondary line vortices with the breaking structures leads to the formation of new hairpin-like vortices, which also contribute to further vortex breakup. These results explain the enhanced mixing performance of confined jets as the mixing is directly tied to the breakup of large vortical structures. Proper orthogonal decomposition modes are also presented to identify the structures/events with the highest contribution to the total turbulent kinetic energy in both flow fields.
{"title":"Effects of axisymmetric confinement on vortical structures emanating from round turbulent jets","authors":"J. Adjetey, K. Fukuda, R. Balachandar","doi":"10.1063/5.0224877","DOIUrl":"https://doi.org/10.1063/5.0224877","url":null,"abstract":"Confined jets occur in many engineering applications including combustion chambers, jet pumps, and chemical reactors. The effects of axisymmetric confinement on the vortical structures identified in a turbulent jet are investigated using large eddy simulation at a Reynolds number of 30 000 (based on nozzle exit conditions) and expansion ratio (chamber-to-nozzle diameter ratio) of five. The results obtained from the confined jet are compared with those of a free jet under the same nozzle exit flow conditions. A prominent recirculation zone forms between the expanding jet and the confining wall, resulting in early shear layer distortion and a shorter interaction length in the confined jet (0.85 jet diameters) compared to the free jet (1.15 jet diameters). Using the λ2 criterion for vortex identification, two dominant structural modes are identified in the near-exit region of the free jet: ring and helical modes. However, in the confined jet, the helical mode is absent, and the turbulent confined fluid accelerates the breakup of the ring vortices. The interaction of the secondary line vortices with the breaking structures leads to the formation of new hairpin-like vortices, which also contribute to further vortex breakup. These results explain the enhanced mixing performance of confined jets as the mixing is directly tied to the breakup of large vortical structures. Proper orthogonal decomposition modes are also presented to identify the structures/events with the highest contribution to the total turbulent kinetic energy in both flow fields.","PeriodicalId":20066,"journal":{"name":"Physics of Fluids","volume":"17 1","pages":""},"PeriodicalIF":4.6,"publicationDate":"2024-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142258779","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The upcoming energy scarcity problem has driven research toward developing energy-efficient two-phase heat exchangers essential for various cooling applications. This research is rooted in the principles of pool boiling, essential for effective heat transfer in various heat exchangers. A well-known reported problem in heat exchangers is the dry-out phenomena of heated surfaces due to bubble coalescence. To tackle this undesirable problem, an innovative technique has been introduced in this study, which involves the shearing of bubbles through liquid jet impingement over the heated surface. The study has been carried out in a two-dimensional domain numerically, in the wall superheat range of 9–16 K. To study the underlying physics involved in this pool boiling phenomenon, the bubble dynamics parameters such as departure frequency, bubble diameter, cold spot (bubble base) temperature, and vapor volume fraction have been analyzed. The results show that with the jet shearing effect, a maximum enhancement of 25% in heat transfer rate is observed at higher wall superheat. The investigation also highlights that the liquid jet enhances vapor volume fraction, indicating enhanced steam generation, particularly an enhancement of 27% observed at elevated wall superheat. An early onset necking effect is also observed with the shearing effect, which leads to the formation of smaller bubbles with higher departure frequencies. This study is a benchmark to the fundamental physics of enhancing two-phase heat transfer through bubble shearing, offering promising insights for energy conservation in two-phase heat exchanger design, particularly within the context of pool boiling.
{"title":"Unveiling the underlying physics of two-phase boiling heat transfer enhancement through shearing of coalescing bubbles","authors":"Niloy Laskar, Mihir K. Das","doi":"10.1063/5.0227551","DOIUrl":"https://doi.org/10.1063/5.0227551","url":null,"abstract":"The upcoming energy scarcity problem has driven research toward developing energy-efficient two-phase heat exchangers essential for various cooling applications. This research is rooted in the principles of pool boiling, essential for effective heat transfer in various heat exchangers. A well-known reported problem in heat exchangers is the dry-out phenomena of heated surfaces due to bubble coalescence. To tackle this undesirable problem, an innovative technique has been introduced in this study, which involves the shearing of bubbles through liquid jet impingement over the heated surface. The study has been carried out in a two-dimensional domain numerically, in the wall superheat range of 9–16 K. To study the underlying physics involved in this pool boiling phenomenon, the bubble dynamics parameters such as departure frequency, bubble diameter, cold spot (bubble base) temperature, and vapor volume fraction have been analyzed. The results show that with the jet shearing effect, a maximum enhancement of 25% in heat transfer rate is observed at higher wall superheat. The investigation also highlights that the liquid jet enhances vapor volume fraction, indicating enhanced steam generation, particularly an enhancement of 27% observed at elevated wall superheat. An early onset necking effect is also observed with the shearing effect, which leads to the formation of smaller bubbles with higher departure frequencies. This study is a benchmark to the fundamental physics of enhancing two-phase heat transfer through bubble shearing, offering promising insights for energy conservation in two-phase heat exchanger design, particularly within the context of pool boiling.","PeriodicalId":20066,"journal":{"name":"Physics of Fluids","volume":"19 1","pages":""},"PeriodicalIF":4.6,"publicationDate":"2024-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142258689","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Aircraft icing seriously threatens flight safety. This paper describes a modified shallow-water icing thermodynamic model that is applicable to unstructured grids and considers the effects of changes in water physical parameters and sublimation on icing. An icing calculation method enables the automatic determination of whether freezing occurs and the type of icing. An autonomous icing calculation program is developed to extract the geometric coordinates and mesh information of the icing surface and combine this information with the multiphase flow field of air-supercooled water droplets. The icing equations in the Godnov format are discretized to produce a large-scale sparse matrix that is solved iteratively using the biconjugate gradient stabilized method. The output includes parameter distributions for the ice thickness, water film thickness, and equilibrium temperature. Taking the National Advisory Committee for Aeronautic 0012 airfoil as the study object, the results of icing simulations are compared with experimental data from an ice wind tunnel and the ice shapes calculated by the FENSAP-ICE software. The results are found to be in good agreement with the experimental data, and the ice height errors at stagnation points are less than 15%. In the case of rime ice, single-horned ice shapes are simulated for both conventional and large supercooled droplets. For mixed ice and glaze ice, double-horned ice shapes are simulated for both droplet conditions. FENSAP-ICE fails to simulate the ice horns and produces large errors in ice thickness and ice range.
{"title":"Research on icing model and calculation methods","authors":"Yaping Hu, Jiheng Pan, Yong Liu, Changxian Zhang, Yuetao Jiang, Jiangnan Zhu","doi":"10.1063/5.0226037","DOIUrl":"https://doi.org/10.1063/5.0226037","url":null,"abstract":"Aircraft icing seriously threatens flight safety. This paper describes a modified shallow-water icing thermodynamic model that is applicable to unstructured grids and considers the effects of changes in water physical parameters and sublimation on icing. An icing calculation method enables the automatic determination of whether freezing occurs and the type of icing. An autonomous icing calculation program is developed to extract the geometric coordinates and mesh information of the icing surface and combine this information with the multiphase flow field of air-supercooled water droplets. The icing equations in the Godnov format are discretized to produce a large-scale sparse matrix that is solved iteratively using the biconjugate gradient stabilized method. The output includes parameter distributions for the ice thickness, water film thickness, and equilibrium temperature. Taking the National Advisory Committee for Aeronautic 0012 airfoil as the study object, the results of icing simulations are compared with experimental data from an ice wind tunnel and the ice shapes calculated by the FENSAP-ICE software. The results are found to be in good agreement with the experimental data, and the ice height errors at stagnation points are less than 15%. In the case of rime ice, single-horned ice shapes are simulated for both conventional and large supercooled droplets. For mixed ice and glaze ice, double-horned ice shapes are simulated for both droplet conditions. FENSAP-ICE fails to simulate the ice horns and produces large errors in ice thickness and ice range.","PeriodicalId":20066,"journal":{"name":"Physics of Fluids","volume":"18 1","pages":""},"PeriodicalIF":4.6,"publicationDate":"2024-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142258725","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Bowen Hu, Yongjie Ren, Rui Sun, Shengcheng Wang, Shanjie Su
Heterogeneity of shale pores at nano-scale and micrometer-scale is of great significance to gas transport properties. In this study, the pore structure of shale samples from lower Silurian Longmaxi Formation in the Sichuan basin is investigated by field emission-scanning electron microscopy (FE-SEM) and x-ray micro-computed tomography (Xμ-CT) technology. Based on fractal theory, the lacunarity is introduced to describe the clustering degree of pores in shale matrix, which can compensate for the limitations of fractal dimension. Combining lacunarity with fractal dimension allows for quantification of subtle differences in pore spatial distribution. For FE-SEM images at nano-scales, the fractal dimension changes in a “U” shape, while lacunarity changes in a “∩” shape. For Xμ-CT images at micrometer-scale, both the fractal dimension and lacunarity change in a logarithmic function. Lacunarity at both nano-scale and micrometer-scale linearly decreases with the increase in fractal dimension. By three-dimensional (3D) pore network modeling analysis, the structure properties of the connected pores, such as the number of pores and throats, pore diameter, pore volume, pore surface, throat length, and coordination number, are quantitatively calculated, and these structure parameters show strong heterogeneity. The average coordination number of the connected pores ranges in 2.92–4.36. This indicates that these pores in shale matrix have poor connectivity. The permeability varies from 0.06 to 0.17 μm2 in two-dimensional (2D) Xμ-CT images but from 3.20 to 34.99 μm2 in a 3D structure. The permeability in the 3D structure is about two order higher in magnitude than that in the 2D Xμ-CT images.
{"title":"Heterogeneity properties and permeability of shale matrix at nano-scale and micron-scale","authors":"Bowen Hu, Yongjie Ren, Rui Sun, Shengcheng Wang, Shanjie Su","doi":"10.1063/5.0223200","DOIUrl":"https://doi.org/10.1063/5.0223200","url":null,"abstract":"Heterogeneity of shale pores at nano-scale and micrometer-scale is of great significance to gas transport properties. In this study, the pore structure of shale samples from lower Silurian Longmaxi Formation in the Sichuan basin is investigated by field emission-scanning electron microscopy (FE-SEM) and x-ray micro-computed tomography (Xμ-CT) technology. Based on fractal theory, the lacunarity is introduced to describe the clustering degree of pores in shale matrix, which can compensate for the limitations of fractal dimension. Combining lacunarity with fractal dimension allows for quantification of subtle differences in pore spatial distribution. For FE-SEM images at nano-scales, the fractal dimension changes in a “U” shape, while lacunarity changes in a “∩” shape. For Xμ-CT images at micrometer-scale, both the fractal dimension and lacunarity change in a logarithmic function. Lacunarity at both nano-scale and micrometer-scale linearly decreases with the increase in fractal dimension. By three-dimensional (3D) pore network modeling analysis, the structure properties of the connected pores, such as the number of pores and throats, pore diameter, pore volume, pore surface, throat length, and coordination number, are quantitatively calculated, and these structure parameters show strong heterogeneity. The average coordination number of the connected pores ranges in 2.92–4.36. This indicates that these pores in shale matrix have poor connectivity. The permeability varies from 0.06 to 0.17 μm2 in two-dimensional (2D) Xμ-CT images but from 3.20 to 34.99 μm2 in a 3D structure. The permeability in the 3D structure is about two order higher in magnitude than that in the 2D Xμ-CT images.","PeriodicalId":20066,"journal":{"name":"Physics of Fluids","volume":"17 1","pages":""},"PeriodicalIF":4.6,"publicationDate":"2024-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142258776","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The internal to external pressure ratio in a large-span cylindrical roof building with a dominant gable opening fluctuates dramatically between 0 and 1, significantly impacted by the attenuation of internal pressure. Current theories usually assume this ratio equal to 1 and overlook the attenuation effect. This study investigates four cylindrical roof models with varying opening areas, scale ratios, and wind speeds by wind tunnel tests. It analyzes ratios of mean (C¯pi/C¯pe), fluctuating (σpi/σpe), and peak (Ĉpi/Ĉpe) internal to external pressure to pinpoint factors affecting the internal pressure attenuation. The results highlight that the most pronounced internal pressure attenuation is at the sideward opening. The vortex shedding around the opening is induced by the wind direction, scale ratio, and wind speed. The attenuation effect decreases with lower frequencies of periodic vortex shedding. This effect generally vanishes when the windward or leeward opening ratio (A1.5/V0) exceeds 0.57%. Empirical design formulas are proposed to predict ratios of internal to external pressure considering the attenuation effect. The inertia (CI) and loss coefficients (CL) affected by the internal pressure attenuation are analyzed to estimate the air slug inertia and damping through the opening. A governing equation, incorporating reduction coefficients (C¯eddy, C̃eddy) from empirical design formulas, is applied to precisely compute the attenuated internal pressure in the large-span cylindrical roof building with a dominant gable opening for engineering risk assessment.
{"title":"Research on the attenuation of internal pressure in the large-span cylindrical roof building with a dominant gable opening","authors":"Yuhang Ge, Ying Sun, Zhenggang Cao, Qiming Zhu","doi":"10.1063/5.0220765","DOIUrl":"https://doi.org/10.1063/5.0220765","url":null,"abstract":"The internal to external pressure ratio in a large-span cylindrical roof building with a dominant gable opening fluctuates dramatically between 0 and 1, significantly impacted by the attenuation of internal pressure. Current theories usually assume this ratio equal to 1 and overlook the attenuation effect. This study investigates four cylindrical roof models with varying opening areas, scale ratios, and wind speeds by wind tunnel tests. It analyzes ratios of mean (C¯pi/C¯pe), fluctuating (σpi/σpe), and peak (Ĉpi/Ĉpe) internal to external pressure to pinpoint factors affecting the internal pressure attenuation. The results highlight that the most pronounced internal pressure attenuation is at the sideward opening. The vortex shedding around the opening is induced by the wind direction, scale ratio, and wind speed. The attenuation effect decreases with lower frequencies of periodic vortex shedding. This effect generally vanishes when the windward or leeward opening ratio (A1.5/V0) exceeds 0.57%. Empirical design formulas are proposed to predict ratios of internal to external pressure considering the attenuation effect. The inertia (CI) and loss coefficients (CL) affected by the internal pressure attenuation are analyzed to estimate the air slug inertia and damping through the opening. A governing equation, incorporating reduction coefficients (C¯eddy, C̃eddy) from empirical design formulas, is applied to precisely compute the attenuated internal pressure in the large-span cylindrical roof building with a dominant gable opening for engineering risk assessment.","PeriodicalId":20066,"journal":{"name":"Physics of Fluids","volume":"10 1","pages":""},"PeriodicalIF":4.6,"publicationDate":"2024-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142258777","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
In this work, a phase-field-based lattice Boltzmann equation (LBE) model for axisymmetric two-phase flow with phase change is proposed. Two sets of discrete particle distribution functions are employed to match the conserved Allen–Cahn equation and the hydrodynamic equations with phase change effect, respectively. Since phase change occurs at the interface, the divergence-free condition of the velocity field is no longer satisfied due to mass transfer, and the conserved Allen–Cahn equation needs to be equipped with a source term dependent on the phase change model. To deal with these, a novel source term in the hydrodynamic LBE is delicately designed to recover the correct target governing equations. Meanwhile, the LBE for the Allen–Cahn equation is modified with a discrete force term to model mass transfer. In particular, an additional correction term is added into the hydrodynamic LBE to reduce the spurious velocity and improve numerical stability. Several axisymmetric benchmark multiphase problems with phase change, including bubble growing in superheated liquid, D2 law, film boiling, bubble rising in superheated liquid under gravity, and droplet impact on a hot surface, have been conducted to test the performance of the proposed model. Numerical results agree well with analytical solutions and available published data in the literature.
{"title":"Axisymmetric phase-field-based lattice Boltzmann model for incompressible two-phase flow with phase change","authors":"Chunhua Zhang, Wenyuan Hou, Qin Lou, Liang Wang, Hantao Liu, Zhaoli Guo","doi":"10.1063/5.0226691","DOIUrl":"https://doi.org/10.1063/5.0226691","url":null,"abstract":"In this work, a phase-field-based lattice Boltzmann equation (LBE) model for axisymmetric two-phase flow with phase change is proposed. Two sets of discrete particle distribution functions are employed to match the conserved Allen–Cahn equation and the hydrodynamic equations with phase change effect, respectively. Since phase change occurs at the interface, the divergence-free condition of the velocity field is no longer satisfied due to mass transfer, and the conserved Allen–Cahn equation needs to be equipped with a source term dependent on the phase change model. To deal with these, a novel source term in the hydrodynamic LBE is delicately designed to recover the correct target governing equations. Meanwhile, the LBE for the Allen–Cahn equation is modified with a discrete force term to model mass transfer. In particular, an additional correction term is added into the hydrodynamic LBE to reduce the spurious velocity and improve numerical stability. Several axisymmetric benchmark multiphase problems with phase change, including bubble growing in superheated liquid, D2 law, film boiling, bubble rising in superheated liquid under gravity, and droplet impact on a hot surface, have been conducted to test the performance of the proposed model. Numerical results agree well with analytical solutions and available published data in the literature.","PeriodicalId":20066,"journal":{"name":"Physics of Fluids","volume":"206 1","pages":""},"PeriodicalIF":4.6,"publicationDate":"2024-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142258783","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
In this paper, we present a parameter-free particle relaxation technique to improve the accuracy and stability of smoothed particle hydrodynamics (SPH). Instead of imposing a background pressure, particles are regularized following the criteria of 0th-order consistency, i.e., the gradient of a constant to be zero. Specifically, the modifications of particles' position are solved by a gradient decent method according to the error between zero value and the gradient of a constant. This modification decreases the integration error and leads a more uniform particles distribution. A set of challenging benchmarks including lid-driven cavity flow, Taylor-Green vortex, FSI (fluid-solid interaction) problem, 2D (two-dimensional) dam-break case, and water exit of a cylinder are investigated to validate the effectiveness of the present technique for addressing the well-known tensile instability and particle clumping problems. Finally, the study of 3D (three-dimensional) dam-break against an obstacle demonstrates the stability and versatility of the present method.
{"title":"A parameter-free particle relaxation technique for smoothed particle hydrodynamics","authors":"Hualin Zheng, Hongfu Qiang, Yujie Zhu, Chi Zhang","doi":"10.1063/5.0223930","DOIUrl":"https://doi.org/10.1063/5.0223930","url":null,"abstract":"In this paper, we present a parameter-free particle relaxation technique to improve the accuracy and stability of smoothed particle hydrodynamics (SPH). Instead of imposing a background pressure, particles are regularized following the criteria of 0th-order consistency, i.e., the gradient of a constant to be zero. Specifically, the modifications of particles' position are solved by a gradient decent method according to the error between zero value and the gradient of a constant. This modification decreases the integration error and leads a more uniform particles distribution. A set of challenging benchmarks including lid-driven cavity flow, Taylor-Green vortex, FSI (fluid-solid interaction) problem, 2D (two-dimensional) dam-break case, and water exit of a cylinder are investigated to validate the effectiveness of the present technique for addressing the well-known tensile instability and particle clumping problems. Finally, the study of 3D (three-dimensional) dam-break against an obstacle demonstrates the stability and versatility of the present method.","PeriodicalId":20066,"journal":{"name":"Physics of Fluids","volume":"4 1","pages":""},"PeriodicalIF":4.6,"publicationDate":"2024-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142258726","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
京公网安备 11010802042870号
京ICP备2023020795号-1
扫码关注我们
求助内容：
应助结果提醒方式：