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}
Large eddy simulations (LES) provide a methodology for both analyzing and simulating multi-scale flows when the smallest scales of motion cannot be resolved. Within environmental flows there exist numerous biogeochemical processes involving tracers undergoing reactions. In this study, we perform an a posteriori LES analysis on a direct numerical simulation of an idealized model reactive tracer subjected to three-dimensional turbulence induced by a Rayleigh–Taylor instability. The governing equations, including an advection–diffusion–reaction equation for the reactive tracer, are filtered, and the resulting sub-filter-scale terms are expressed in terms of interactions between scales. The procedure is demonstrated for a generalized degree N polynomial reaction function. Various spectral filters are applied to the data and compared. The preferential choice is to use the widest filter possible with a smoothed cutoff. The sub-filter-scale reaction term that results from filtering the reaction function is considered for each of the filter choices. When using a particularly harsh filter, local balances are found for the resolved scale and cross-scale components of the sub-filter-scale reaction term. The same result is shown for the vertical sub-filter-scale flux for both a reactive and a passive tracer. The components of the sub-filter-scale reaction and vertical flux terms involving interactions at the sub-filter-scale do not show any evidence of local balances and are distributed around the fine turbulent structures in the flow. This suggests that parameterizations for the sub-filter-scale terms would benefit from considering event specific dynamics.
大型涡流模拟（LES）提供了一种在无法解析最小运动尺度的情况下分析和模拟多尺度流动的方法。在环境流动中，存在着许多涉及示踪剂反应的生物地球化学过程。在本研究中，我们对受到雷利-泰勒不稳定性诱导的三维湍流影响的理想化反应示踪剂模型进行了直接数值模拟，并对 LES 进行了后验分析。对包括反应示踪剂的平流-扩散-反应方程式在内的控制方程进行过滤，由此产生的子过滤尺度项用尺度间的相互作用来表示。该程序针对广义的 N 级多项式反应函数进行了演示。对数据应用了各种光谱滤波器并进行了比较。优先选择使用尽可能宽的滤波器和平滑截止。在每种滤波器选择中，都要考虑对反应函数进行滤波后产生的子滤波器尺度反应项。当使用特别苛刻的滤波器时，会发现子滤波器尺度反应项的分辨尺度和跨尺度分量的局部平衡。对于反应式和被动式示踪剂的垂直亚滤波尺度通量，也显示了同样的结果。亚滤波器尺度反应项和垂直通量项中涉及亚滤波器尺度相互作用的部分没有显示出任何局部平衡的迹象，而是分布在流体中精细湍流结构的周围。这表明，亚滤波尺度项的参数化将受益于考虑特定事件的动力学。
{"title":"Large eddy simulation analysis of a model reactive tracer through spatial filtering","authors":"S. Legare, M. Stastna","doi":"10.1063/5.0226039","DOIUrl":"https://doi.org/10.1063/5.0226039","url":null,"abstract":"Large eddy simulations (LES) provide a methodology for both analyzing and simulating multi-scale flows when the smallest scales of motion cannot be resolved. Within environmental flows there exist numerous biogeochemical processes involving tracers undergoing reactions. In this study, we perform an a posteriori LES analysis on a direct numerical simulation of an idealized model reactive tracer subjected to three-dimensional turbulence induced by a Rayleigh–Taylor instability. The governing equations, including an advection–diffusion–reaction equation for the reactive tracer, are filtered, and the resulting sub-filter-scale terms are expressed in terms of interactions between scales. The procedure is demonstrated for a generalized degree N polynomial reaction function. Various spectral filters are applied to the data and compared. The preferential choice is to use the widest filter possible with a smoothed cutoff. The sub-filter-scale reaction term that results from filtering the reaction function is considered for each of the filter choices. When using a particularly harsh filter, local balances are found for the resolved scale and cross-scale components of the sub-filter-scale reaction term. The same result is shown for the vertical sub-filter-scale flux for both a reactive and a passive tracer. The components of the sub-filter-scale reaction and vertical flux terms involving interactions at the sub-filter-scale do not show any evidence of local balances and are distributed around the fine turbulent structures in the flow. This suggests that parameterizations for the sub-filter-scale terms would benefit from considering event specific dynamics.","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":"142258778","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 general solution of Stokes equations for the problem of an axisymmetric oscillatory flow of an incompressible, viscous fluid past a sphere satisfying general boundary conditions is obtained. The behavior of the magnitude of drag is observed with the variation of the slip parameter. Some more interesting behaviors are detailed, and several existing results pertaining to steady flows and flows with rigid and shear free boundary conditions are recovered as special cases. The corresponding results are discussed for four different axisymmetric oscillatory Stokes flows, namely, uniform flow, flows generated due to a dipole, a source, and a Stokeslet. A few interesting streamline patterns like formation, elongation, and disappearance of viscous eddies in the vicinity of the sphere with a periodic reversal of the flow are observed at different frequencies for different values of the slip parameter.
{"title":"The effect of slip parameter in an axisymmetric oscillatory Stokes flow","authors":"Dadi Dimple S. S., B. Sri Padmavati","doi":"10.1063/5.0226315","DOIUrl":"https://doi.org/10.1063/5.0226315","url":null,"abstract":"A general solution of Stokes equations for the problem of an axisymmetric oscillatory flow of an incompressible, viscous fluid past a sphere satisfying general boundary conditions is obtained. The behavior of the magnitude of drag is observed with the variation of the slip parameter. Some more interesting behaviors are detailed, and several existing results pertaining to steady flows and flows with rigid and shear free boundary conditions are recovered as special cases. The corresponding results are discussed for four different axisymmetric oscillatory Stokes flows, namely, uniform flow, flows generated due to a dipole, a source, and a Stokeslet. A few interesting streamline patterns like formation, elongation, and disappearance of viscous eddies in the vicinity of the sphere with a periodic reversal of the flow are observed at different frequencies for different values of the slip parameter.","PeriodicalId":20066,"journal":{"name":"Physics of Fluids","volume":"21 1","pages":""},"PeriodicalIF":4.6,"publicationDate":"2024-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142269345","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
扫码关注我们
求助内容：
应助结果提醒方式：