Mixing describes the process by which solutes evolve from an initial heterogeneous state to uniformity under the stirring action of a fluid flow. Fluid stretching forms thin scalar lamellae that coalesce due to molecular diffusion. Owing to the linearity of the advection–diffusion equation, coalescence can be envisioned as an aggregation process. Here, we demonstrate that in smooth two-dimensional chaotic flows, mixing obeys a correlated aggregation process, where the spatial distribution of the number of lamellae in aggregates is highly correlated with their elongation, and is set by the fractal properties of the advected material lines. We show that the presence of correlations makes mixing less efficient than a completely random aggregation process because lamellae with similar elongations and scalar levels tend to remain isolated from each other. We show that correlated aggregation is uniquely determined by a single exponent that quantifies the effective number of random aggregation events. These findings expand aggregation theories to a larger class of systems, which have relevance to various fundamental and applied mixing problems.
{"title":"Mixing as a correlated aggregation process","authors":"J. Heyman, E. Villermaux, P. Davy, T. Le Borgne","doi":"10.1017/jfm.2024.537","DOIUrl":"https://doi.org/10.1017/jfm.2024.537","url":null,"abstract":"Mixing describes the process by which solutes evolve from an initial heterogeneous state to uniformity under the stirring action of a fluid flow. Fluid stretching forms thin scalar lamellae that coalesce due to molecular diffusion. Owing to the linearity of the advection–diffusion equation, coalescence can be envisioned as an aggregation process. Here, we demonstrate that in smooth two-dimensional chaotic flows, mixing obeys a correlated aggregation process, where the spatial distribution of the number of lamellae in aggregates is highly correlated with their elongation, and is set by the fractal properties of the advected material lines. We show that the presence of correlations makes mixing less efficient than a completely random aggregation process because lamellae with similar elongations and scalar levels tend to remain isolated from each other. We show that correlated aggregation is uniquely determined by a single exponent that quantifies the effective number of random aggregation events. These findings expand aggregation theories to a larger class of systems, which have relevance to various fundamental and applied mixing problems.","PeriodicalId":15853,"journal":{"name":"Journal of Fluid Mechanics","volume":"19 1","pages":""},"PeriodicalIF":3.7,"publicationDate":"2024-08-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142207465","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}
Dung Viet Duong, Luc Van Nguyen, Duc Van Nguyen, Truong Cong Dinh, Lavi Rizki Zuhal, Long Ich Ngo
Comprehensive coherent structures around a surface-mounted low aspect ratio square cylinder in uniform flow with an oblique angle of $45^{circ }$ were investigated for cylinder-width-based Reynolds numbers of 3000 and 10 000 by direct numerical simulation based on a topology-confined mesh refinement framework. High-resolution simulations and the critical-point concept were scrutinized to reveal for the first time the reasonable and compatible topologies of flow separation and complete near-wall structures, due to their extensive impact on various engineering applications. Large-scale horseshoe vortices are observed at two notable foci in the viscous sublayer. Within this layer, a wall-parallel jet is formed by downflow intruding into the bottom surface at the half-saddle point, then deflecting in the upstream direction and finally penetrating the bottom surface until the half-saddle point. A pair of conical vortices on the cylinder's top surface switch themselves on two sides of the square cylinder, where the switching frequency is identical with that of the sway of the side shear layer. The undulation of the Kelvin–Helmholtz instability is identified in the instantaneous development of a conical vortex and side shear layer, where the scaling of the ratio of the Kelvin–Helmholtz and von Kármán frequencies follows the power-law relation obtained by Lander et al. (J. Fluid Mech., vol. 849, 2018, pp. 1096–1119). Large-scale arch-shaped vortex is often detected in the intermediate wake region of a square cylinder, involving two interconnected portions, such as the leg portion separated from leeward surfaces and head portion rolled up from the top surface. The leg portion of the arch-shaped vortex was rooted by two foci near the bottom-surface plane.
{"title":"Direct numerical simulation of 45° oblique flow past surface-mounted square cylinder","authors":"Dung Viet Duong, Luc Van Nguyen, Duc Van Nguyen, Truong Cong Dinh, Lavi Rizki Zuhal, Long Ich Ngo","doi":"10.1017/jfm.2024.554","DOIUrl":"https://doi.org/10.1017/jfm.2024.554","url":null,"abstract":"Comprehensive coherent structures around a surface-mounted low aspect ratio square cylinder in uniform flow with an oblique angle of <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" mime-subtype=\"png\" xlink:href=\"S0022112024005548_inline2.png\"/> <jats:tex-math>$45^{circ }$</jats:tex-math> </jats:alternatives> </jats:inline-formula> were investigated for cylinder-width-based Reynolds numbers of 3000 and 10 000 by direct numerical simulation based on a topology-confined mesh refinement framework. High-resolution simulations and the critical-point concept were scrutinized to reveal for the first time the reasonable and compatible topologies of flow separation and complete near-wall structures, due to their extensive impact on various engineering applications. Large-scale horseshoe vortices are observed at two notable foci in the viscous sublayer. Within this layer, a wall-parallel jet is formed by downflow intruding into the bottom surface at the half-saddle point, then deflecting in the upstream direction and finally penetrating the bottom surface until the half-saddle point. A pair of conical vortices on the cylinder's top surface switch themselves on two sides of the square cylinder, where the switching frequency is identical with that of the sway of the side shear layer. The undulation of the Kelvin–Helmholtz instability is identified in the instantaneous development of a conical vortex and side shear layer, where the scaling of the ratio of the Kelvin–Helmholtz and von Kármán frequencies follows the power-law relation obtained by Lander <jats:italic>et al.</jats:italic> (<jats:italic>J. Fluid Mech.</jats:italic>, vol. 849, 2018, pp. 1096–1119). Large-scale arch-shaped vortex is often detected in the intermediate wake region of a square cylinder, involving two interconnected portions, such as the leg portion separated from leeward surfaces and head portion rolled up from the top surface. The leg portion of the arch-shaped vortex was rooted by two foci near the bottom-surface plane.","PeriodicalId":15853,"journal":{"name":"Journal of Fluid Mechanics","volume":"25 1","pages":""},"PeriodicalIF":3.7,"publicationDate":"2024-08-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142207463","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}
Stylianos Varchanis, Eliane Younes, Simon J. Haward, Amy Q. Shen
The motion of rigid particles in complex fluids is ubiquitous in natural and industrial processes. The most popular toy model for understanding the physics of such systems is the settling of a solid sphere in a viscoelastic fluid. There is general agreement that an elastic wake develops downstream of the sphere, causing the breakage of fore-and-aft symmetry, while the flow remains axisymmetric, independent of fluid viscoelasticity and flow conditions. Using a continuum mechanics model, we reveal that axisymmetry holds only for weak viscoelastic flows. Beyond a critical value of the settling velocity, steady, non-axisymmetric disturbances develop peripherally of the rear pole of the sphere, giving rise to a four-lobed fingering instability. The transition from axisymmetric to non-axisymmetric flow fields is characterized by a regular bifurcation and depends solely on the interplay between shear and extensional properties of the viscoelastic fluid under different flow regimes. At higher settling velocities, each lobe tip is split into two new lobes, resembling fractal fingering in interfacial flows. For the first time, we capture an elastic fingering instability under steady-state conditions, and provide the missing information for understanding and predicting such instabilities in the response of viscoelastic fluids and soft media.
{"title":"Emergence of lobed wakes during the sedimentation of spheres in viscoelastic fluids","authors":"Stylianos Varchanis, Eliane Younes, Simon J. Haward, Amy Q. Shen","doi":"10.1017/jfm.2024.459","DOIUrl":"https://doi.org/10.1017/jfm.2024.459","url":null,"abstract":"The motion of rigid particles in complex fluids is ubiquitous in natural and industrial processes. The most popular toy model for understanding the physics of such systems is the settling of a solid sphere in a viscoelastic fluid. There is general agreement that an elastic wake develops downstream of the sphere, causing the breakage of fore-and-aft symmetry, while the flow remains axisymmetric, independent of fluid viscoelasticity and flow conditions. Using a continuum mechanics model, we reveal that axisymmetry holds only for weak viscoelastic flows. Beyond a critical value of the settling velocity, steady, non-axisymmetric disturbances develop peripherally of the rear pole of the sphere, giving rise to a four-lobed fingering instability. The transition from axisymmetric to non-axisymmetric flow fields is characterized by a regular bifurcation and depends solely on the interplay between shear and extensional properties of the viscoelastic fluid under different flow regimes. At higher settling velocities, each lobe tip is split into two new lobes, resembling fractal fingering in interfacial flows. For the first time, we capture an elastic fingering instability under steady-state conditions, and provide the missing information for understanding and predicting such instabilities in the response of viscoelastic fluids and soft media.","PeriodicalId":15853,"journal":{"name":"Journal of Fluid Mechanics","volume":"61 1","pages":""},"PeriodicalIF":3.7,"publicationDate":"2024-08-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142226603","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}
Alexander O. Korotkevich, Sergey V. Nazarenko, Yulin Pan, Jalal Shatah
The $k^{-23/6}$ wave action spectrum with an inverse cascade is one of the fundamental Kolmogorov–Zakharov solutions for gravity wave turbulence, which is part of the citation for the Dirac Medal in 2003. Instead of confirming this solution, however, several existing simulations and experiments suggest a spectrum of $k^{-3}$ in set-ups corresponding to the inverse cascade. We provide a theoretical explanation for the latter, considering the condensate that naturally forms in finite domains of experiments/simulations. Our new theory hinges on: (1) derivation of a spectral diffusion equation when non-local interactions with the condensate become dominant, for the first time systematically formulated for quartet-interaction systems; and (2) careful analysis of the asymptotics of interaction coefficient with a remarkable cancellation of all leading-order terms.
{"title":"Non-local gravity wave turbulence in presence of condensate","authors":"Alexander O. Korotkevich, Sergey V. Nazarenko, Yulin Pan, Jalal Shatah","doi":"10.1017/jfm.2024.423","DOIUrl":"https://doi.org/10.1017/jfm.2024.423","url":null,"abstract":"The <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" mime-subtype=\"png\" xlink:href=\"S0022112024004233_inline1.png\"/> <jats:tex-math>$k^{-23/6}$</jats:tex-math> </jats:alternatives> </jats:inline-formula> wave action spectrum with an inverse cascade is one of the fundamental Kolmogorov–Zakharov solutions for gravity wave turbulence, which is part of the citation for the Dirac Medal in 2003. Instead of confirming this solution, however, several existing simulations and experiments suggest a spectrum of <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" mime-subtype=\"png\" xlink:href=\"S0022112024004233_inline2.png\"/> <jats:tex-math>$k^{-3}$</jats:tex-math> </jats:alternatives> </jats:inline-formula> in set-ups corresponding to the inverse cascade. We provide a theoretical explanation for the latter, considering the condensate that naturally forms in finite domains of experiments/simulations. Our new theory hinges on: (1) derivation of a spectral diffusion equation when non-local interactions with the condensate become dominant, for the first time systematically formulated for quartet-interaction systems; and (2) careful analysis of the asymptotics of interaction coefficient with a remarkable cancellation of all leading-order terms.","PeriodicalId":15853,"journal":{"name":"Journal of Fluid Mechanics","volume":"7 1","pages":""},"PeriodicalIF":3.7,"publicationDate":"2024-08-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142226605","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}
David Barreiro-Villaverde, Anne Gosset, Marcos Lema, Miguel A. Mendez
We investigate the dynamics of a gas jet impinging perpendicular to a thin liquid film dragged by a rising vertical substrate. This configuration is relevant to the jet-wiping process in hot-dip galvanization and it is unstable. Previous studies analysed the dynamics of the instability in the case of liquids with low Kapitza numbers (highly viscous liquids), more amenable to experimental and numerical investigations. This work extends the previous investigations by focusing on the wiping at much higher Kapitza numbers, which are more relevant to the galvanizing process. The simulations are carried out by combining volume of fluid and large-eddy simulations, and the dynamics of the gas–liquid interaction is analysed using extended multiscale proper orthogonal decomposition. The simulations allowed for analysing the jet-wiping instability in new flow conditions. Despite the largely different conditions, the results show that the interaction between the gas jet and the liquid film is qualitatively similar, featuring two-dimensional waves in the liquid correlated with oscillations and deflections of the gas jet in all cases. The wave characteristics (e.g. frequency and propagation speed) scale remarkably well using the Shkadov-like scaling based on the liquid, suggesting a dominant role of the liquid film in the coupling, and potentially enabling extrapolation of the results to a broader range of wiping conditions. Finally, we use the numerical results to discuss the limitations of liquid-film models, which constitute currently the only possible approach to study the jet-wiping process in industrial conditions.
{"title":"On the coupling instability of a gas jet impinging on a liquid film","authors":"David Barreiro-Villaverde, Anne Gosset, Marcos Lema, Miguel A. Mendez","doi":"10.1017/jfm.2024.553","DOIUrl":"https://doi.org/10.1017/jfm.2024.553","url":null,"abstract":"We investigate the dynamics of a gas jet impinging perpendicular to a thin liquid film dragged by a rising vertical substrate. This configuration is relevant to the jet-wiping process in hot-dip galvanization and it is unstable. Previous studies analysed the dynamics of the instability in the case of liquids with low Kapitza numbers (highly viscous liquids), more amenable to experimental and numerical investigations. This work extends the previous investigations by focusing on the wiping at much higher Kapitza numbers, which are more relevant to the galvanizing process. The simulations are carried out by combining volume of fluid and large-eddy simulations, and the dynamics of the gas–liquid interaction is analysed using extended multiscale proper orthogonal decomposition. The simulations allowed for analysing the jet-wiping instability in new flow conditions. Despite the largely different conditions, the results show that the interaction between the gas jet and the liquid film is qualitatively similar, featuring two-dimensional waves in the liquid correlated with oscillations and deflections of the gas jet in all cases. The wave characteristics (e.g. frequency and propagation speed) scale remarkably well using the Shkadov-like scaling based on the liquid, suggesting a dominant role of the liquid film in the coupling, and potentially enabling extrapolation of the results to a broader range of wiping conditions. Finally, we use the numerical results to discuss the limitations of liquid-film models, which constitute currently the only possible approach to study the jet-wiping process in industrial conditions.","PeriodicalId":15853,"journal":{"name":"Journal of Fluid Mechanics","volume":"16 1","pages":""},"PeriodicalIF":3.7,"publicationDate":"2024-08-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142207464","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}
Thin-film equations are utilised in many different areas of fluid dynamics when there exists a direction in which the aspect ratio can be considered small. We consider thin free films with Marangoni effects in the extensional flow regime, where velocity gradients occur predominantly along the film. In practice, because of the local deposition of surfactants or input of energy, asymmetric distributions of surfactants or surface tension more generally, are possible. Such examples include the surface of bubbles and the rupture of thin films. In this study, we consider the asymmetric thin-film equations for extensional flow with Marangoni effects. Concentrating on the case of small Reynolds number $ Re $, we study the deposition of insoluble surfactants on one side of a liquid sheet otherwise at rest and the resulting thinning and rupture of the sheet. The analogous problem with a uniformly thinning liquid sheet is also considered. In addition, the centreline deformation is discussed. In particular, we show analytically that if the surface tension isotherm $sigma = sigma (varGamma )$ is nonlinear (surface tension $sigma$ varies with surfactant concentration $varGamma$), then accounting for top–bottom asymmetry leads to slower (faster) thinning and pinching if $sigma = sigma (varGamma )$ is convex (concave). The analytical progress reported in this paper allows us to discuss the production of satellite drops from rupture via Marangoni effects, which, if relevant to surface bubbles, would be an aerosol production mechanism that is distinct from jet drops and film drops.
{"title":"Thin-film flow due to an asymmetric distribution of surface tension and applications to surfactant deposition","authors":"Jun Eshima, Luc Deike, Howard A. Stone","doi":"10.1017/jfm.2024.501","DOIUrl":"https://doi.org/10.1017/jfm.2024.501","url":null,"abstract":"Thin-film equations are utilised in many different areas of fluid dynamics when there exists a direction in which the aspect ratio can be considered small. We consider thin free films with Marangoni effects in the extensional flow regime, where velocity gradients occur predominantly along the film. In practice, because of the local deposition of surfactants or input of energy, asymmetric distributions of surfactants or surface tension more generally, are possible. Such examples include the surface of bubbles and the rupture of thin films. In this study, we consider the asymmetric thin-film equations for extensional flow with Marangoni effects. Concentrating on the case of small Reynolds number <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" mime-subtype=\"png\" xlink:href=\"S0022112024005019_inline1.png\"/> <jats:tex-math>$ Re $</jats:tex-math> </jats:alternatives> </jats:inline-formula>, we study the deposition of insoluble surfactants on one side of a liquid sheet otherwise at rest and the resulting thinning and rupture of the sheet. The analogous problem with a uniformly thinning liquid sheet is also considered. In addition, the centreline deformation is discussed. In particular, we show analytically that if the surface tension isotherm <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" mime-subtype=\"png\" xlink:href=\"S0022112024005019_inline2.png\"/> <jats:tex-math>$sigma = sigma (varGamma )$</jats:tex-math> </jats:alternatives> </jats:inline-formula> is nonlinear (surface tension <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" mime-subtype=\"png\" xlink:href=\"S0022112024005019_inline3.png\"/> <jats:tex-math>$sigma$</jats:tex-math> </jats:alternatives> </jats:inline-formula> varies with surfactant concentration <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" mime-subtype=\"png\" xlink:href=\"S0022112024005019_inline4.png\"/> <jats:tex-math>$varGamma$</jats:tex-math> </jats:alternatives> </jats:inline-formula>), then accounting for top–bottom asymmetry leads to slower (faster) thinning and pinching if <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" mime-subtype=\"png\" xlink:href=\"S0022112024005019_inline5.png\"/> <jats:tex-math>$sigma = sigma (varGamma )$</jats:tex-math> </jats:alternatives> </jats:inline-formula> is convex (concave). The analytical progress reported in this paper allows us to discuss the production of satellite drops from rupture via Marangoni effects, which, if relevant to surface bubbles, would be an aerosol production mechanism that is distinct from jet drops and film drops.","PeriodicalId":15853,"journal":{"name":"Journal of Fluid Mechanics","volume":"63 1","pages":""},"PeriodicalIF":3.7,"publicationDate":"2024-08-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142226604","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}
Si Suo, Seyed Morteza Habibi Khorasani, Shervin Bagheri
In this study, we investigate the stability of a film that is attached to a corner between a cylinder and a substrate, using a combination of theoretical and numerical approaches. Notably, we place our focus on flat and thin films where the contact line is almost perpendicular to the cylinder wall whereas a small angle forms between the contact line and the substrate, and the film size is smaller than the cylinder radius. The film stability, which depends on the film size and the wall wettability, is first predicted by a standard linear stability analysis (LSA) within the long-wave theoretical framework. We find that the film size plays the most important role in controlling the film stability. Specifically, the thicker the film is, the less sensitive it becomes to the large-wavenumber perturbation. The wall wettability mainly impacts the growth rates of perturbations and slightly influences the marginal stability and postinstability patterns of wrapping films. We compare the LSA predictions with numerical results obtained from a disjoining pressure model (DPM) and volume-of-fluid (VOF) simulations, which provide more insights into the film breakup process. At the early stage there is a strong agreement between the LSA predictions and the DPM results. Notably, as the perturbation grows, thin film regions connecting two neighbouring satellite droplets form which may eventually lead to a stable or temporary secondary droplet, an aspect which the LSA is incapable of capturing. In addition, the VOF simulations suggest that beyond a critical film size, merging between two neighbouring drops becomes involved during the breakup stage. Therefore, the LSA predictions are able to provide only an upper limit on the final number of satellite droplets.
{"title":"Dewetting of a corner film wrapping a wall-mounted cylinder","authors":"Si Suo, Seyed Morteza Habibi Khorasani, Shervin Bagheri","doi":"10.1017/jfm.2024.416","DOIUrl":"https://doi.org/10.1017/jfm.2024.416","url":null,"abstract":"In this study, we investigate the stability of a film that is attached to a corner between a cylinder and a substrate, using a combination of theoretical and numerical approaches. Notably, we place our focus on flat and thin films where the contact line is almost perpendicular to the cylinder wall whereas a small angle forms between the contact line and the substrate, and the film size is smaller than the cylinder radius. The film stability, which depends on the film size and the wall wettability, is first predicted by a standard linear stability analysis (LSA) within the long-wave theoretical framework. We find that the film size plays the most important role in controlling the film stability. Specifically, the thicker the film is, the less sensitive it becomes to the large-wavenumber perturbation. The wall wettability mainly impacts the growth rates of perturbations and slightly influences the marginal stability and postinstability patterns of wrapping films. We compare the LSA predictions with numerical results obtained from a disjoining pressure model (DPM) and volume-of-fluid (VOF) simulations, which provide more insights into the film breakup process. At the early stage there is a strong agreement between the LSA predictions and the DPM results. Notably, as the perturbation grows, thin film regions connecting two neighbouring satellite droplets form which may eventually lead to a stable or temporary secondary droplet, an aspect which the LSA is incapable of capturing. In addition, the VOF simulations suggest that beyond a critical film size, merging between two neighbouring drops becomes involved during the breakup stage. Therefore, the LSA predictions are able to provide only an upper limit on the final number of satellite droplets.","PeriodicalId":15853,"journal":{"name":"Journal of Fluid Mechanics","volume":"9 1","pages":""},"PeriodicalIF":3.7,"publicationDate":"2024-08-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142207460","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}
Recently, direct numerical simulations (DNS) of stably stratified turbulence have shown that as the Prandtl number ($Pr$) is increased from 1 to 7, the mean turbulent potential energy dissipation rate (TPE-DR) drops dramatically, while the mean turbulent kinetic energy dissipation rate (TKE-DR) increases significantly. Through an analysis of the equations governing the fluctuating velocity and density gradients we provide a mechanistic explanation for this surprising behaviour and test the predictions using DNS. We show that the mean density gradient gives rise to a mechanism that opposes the production of fluctuating density gradients, and this is connected to the emergence of ramp cliffs. The same term appears in the velocity gradient equation but with the opposite sign, and is the contribution from buoyancy. This term is ultimately the reason why the TPE-DR reduces while the TKE-DR increases with increasing $Pr$. Our analysis also predicts that the effects of buoyancy on the smallest scales of the flow become stronger as $Pr$ is increased, and this is confirmed by our DNS data. A consequence of this is that the standard buoyancy Reynolds number does not correctly estimate the impact of buoyancy at the smallest scales when $Pr$ deviates from 1, and we derive a suitable alternative parameter. Finally, an analysis of the filtered gradient equations reveals that the mean density gradient term changes sign at sufficiently large scales, such that buoyancy acts as a source for velocity gradients at small scales, but as a sink at large scales.
{"title":"Understanding the effect of Prandtl number on momentum and scalar mixing rates in neutral and stably stratified flows using gradient field dynamics","authors":"Andrew D. Bragg, Stephen M. de Bruyn Kops","doi":"10.1017/jfm.2024.548","DOIUrl":"https://doi.org/10.1017/jfm.2024.548","url":null,"abstract":"Recently, direct numerical simulations (DNS) of stably stratified turbulence have shown that as the Prandtl number (<jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" mime-subtype=\"png\" xlink:href=\"S0022112024005482_inline1.png\"/> <jats:tex-math>$Pr$</jats:tex-math> </jats:alternatives> </jats:inline-formula>) is increased from 1 to 7, the mean turbulent potential energy dissipation rate (TPE-DR) drops dramatically, while the mean turbulent kinetic energy dissipation rate (TKE-DR) increases significantly. Through an analysis of the equations governing the fluctuating velocity and density gradients we provide a mechanistic explanation for this surprising behaviour and test the predictions using DNS. We show that the mean density gradient gives rise to a mechanism that opposes the production of fluctuating density gradients, and this is connected to the emergence of ramp cliffs. The same term appears in the velocity gradient equation but with the opposite sign, and is the contribution from buoyancy. This term is ultimately the reason why the TPE-DR reduces while the TKE-DR increases with increasing <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" mime-subtype=\"png\" xlink:href=\"S0022112024005482_inline2.png\"/> <jats:tex-math>$Pr$</jats:tex-math> </jats:alternatives> </jats:inline-formula>. Our analysis also predicts that the effects of buoyancy on the smallest scales of the flow become stronger as <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" mime-subtype=\"png\" xlink:href=\"S0022112024005482_inline3.png\"/> <jats:tex-math>$Pr$</jats:tex-math> </jats:alternatives> </jats:inline-formula> is increased, and this is confirmed by our DNS data. A consequence of this is that the standard buoyancy Reynolds number does not correctly estimate the impact of buoyancy at the smallest scales when <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" mime-subtype=\"png\" xlink:href=\"S0022112024005482_inline4.png\"/> <jats:tex-math>$Pr$</jats:tex-math> </jats:alternatives> </jats:inline-formula> deviates from 1, and we derive a suitable alternative parameter. Finally, an analysis of the filtered gradient equations reveals that the mean density gradient term changes sign at sufficiently large scales, such that buoyancy acts as a source for velocity gradients at small scales, but as a sink at large scales.","PeriodicalId":15853,"journal":{"name":"Journal of Fluid Mechanics","volume":"46 1","pages":""},"PeriodicalIF":3.7,"publicationDate":"2024-08-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142207466","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}
Ashish Mishra, George Mamatsashvili, Martin Seilmayer, Frank Stefani
The effects of thermal convection on turbulence in accretion discs, and particularly its interplay with the magnetorotational instability (MRI), are of significant astrophysical interest. Despite extensive theoretical and numerical studies, such an interplay has not been explored experimentally. We conduct linear analysis of the azimuthal version of MRI (AMRI) in the presence of thermal convection and compare the results with our experimental data published before. We show that the critical Hartmann number ($Ha$) for the onset of AMRI is reduced by convection. Importantly, convection breaks symmetry between $m = pm 1$ instability modes ($m$ is the azimuthal wavenumber). This preference for one mode over the other makes the AMRI wave appear as a ‘one-winged butterfly’.
{"title":"One-winged butterflies: mode selection for azimuthal magnetorotational instability by thermal convection","authors":"Ashish Mishra, George Mamatsashvili, Martin Seilmayer, Frank Stefani","doi":"10.1017/jfm.2024.517","DOIUrl":"https://doi.org/10.1017/jfm.2024.517","url":null,"abstract":"The effects of thermal convection on turbulence in accretion discs, and particularly its interplay with the magnetorotational instability (MRI), are of significant astrophysical interest. Despite extensive theoretical and numerical studies, such an interplay has not been explored experimentally. We conduct linear analysis of the azimuthal version of MRI (AMRI) in the presence of thermal convection and compare the results with our experimental data published before. We show that the critical Hartmann number (<jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" mime-subtype=\"png\" xlink:href=\"S0022112024005172_inline1.png\"/> <jats:tex-math>$Ha$</jats:tex-math> </jats:alternatives> </jats:inline-formula>) for the onset of AMRI is reduced by convection. Importantly, convection breaks symmetry between <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" mime-subtype=\"png\" xlink:href=\"S0022112024005172_inline2.png\"/> <jats:tex-math>$m = pm 1$</jats:tex-math> </jats:alternatives> </jats:inline-formula> instability modes (<jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" mime-subtype=\"png\" xlink:href=\"S0022112024005172_inline3.png\"/> <jats:tex-math>$m$</jats:tex-math> </jats:alternatives> </jats:inline-formula> is the azimuthal wavenumber). This preference for one mode over the other makes the AMRI wave appear as a ‘one-winged butterfly’.","PeriodicalId":15853,"journal":{"name":"Journal of Fluid Mechanics","volume":"11 1","pages":""},"PeriodicalIF":3.7,"publicationDate":"2024-08-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142207469","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 theory of incompressible turbulent plane jets (TPJs) is proposed by advancing an improved boundary layer approximation over the limiting classical – retaining more terms in the momentum balance equations. A pressure deficit inside the jet (with respect to the ambient) must exist due to transverse turbulence (Miller & Comings, <jats:italic>J. Fluid Mech.</jats:italic>, vol. 3, 1957, pp. 1–16; Hussain & Clarke, <jats:italic>Phys. Fluids</jats:italic>, vol. 20, 1977, pp. 1416–1426). Contrary to the universally accepted invariance of the total momentum flux <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" mime-subtype="png" xlink:href="S0022112024005275_inline1.png"/> <jats:tex-math>$J_T(x)$</jats:tex-math> </jats:alternatives> </jats:inline-formula> (non-dimensionalized by its inlet value) as a function of the streamwise distance <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" mime-subtype="png" xlink:href="S0022112024005275_inline2.png"/> <jats:tex-math>$x$</jats:tex-math> </jats:alternatives> </jats:inline-formula>, we prove that <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" mime-subtype="png" xlink:href="S0022112024005275_inline3.png"/> <jats:tex-math>$J_T(x) >1$</jats:tex-math> </jats:alternatives> </jats:inline-formula> – a condition that all TPJs must satisfy; surprisingly, prior theories and most experiments do not satisfy this condition. This motivated us to apply Lie symmetry analysis with translational and dilatational transformations of the modified equations (incorporating <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" mime-subtype="png" xlink:href="S0022112024005275_inline4.png"/> <jats:tex-math>$J_T>1$</jats:tex-math> </jats:alternatives> </jats:inline-formula>), which yields scaling laws for key jet measures: the mean streamwise and transverse velocities <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" mime-subtype="png" xlink:href="S0022112024005275_inline5.png"/> <jats:tex-math>$U(x,y)$</jats:tex-math> </jats:alternatives> </jats:inline-formula> and <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" mime-subtype="png" xlink:href="S0022112024005275_inline6.png"/> <jats:tex-math>$V(x,y)$</jats:tex-math> </jats:alternatives> </jats:inline-formula>, the turbulence intensities, the Reynolds shear stress <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink="http://www.w3.org/1999/xlink" mime-subtype="png" xlink:href="S0022112024005275_inline8.png"/> <jats:tex-math>$-rho ,overline {u'v'}(x,y)$</jats:tex-math> </jats:alternatives> </jats:inline-formula>, the mean pressure <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink="http://www.w3.org
{"title":"Turbulent jet theory via Lie symmetry analysis: the free plane jet","authors":"Nadeem A. Malik, Fazle Hussain","doi":"10.1017/jfm.2024.527","DOIUrl":"https://doi.org/10.1017/jfm.2024.527","url":null,"abstract":"A theory of incompressible turbulent plane jets (TPJs) is proposed by advancing an improved boundary layer approximation over the limiting classical – retaining more terms in the momentum balance equations. A pressure deficit inside the jet (with respect to the ambient) must exist due to transverse turbulence (Miller & Comings, <jats:italic>J. Fluid Mech.</jats:italic>, vol. 3, 1957, pp. 1–16; Hussain & Clarke, <jats:italic>Phys. Fluids</jats:italic>, vol. 20, 1977, pp. 1416–1426). Contrary to the universally accepted invariance of the total momentum flux <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" mime-subtype=\"png\" xlink:href=\"S0022112024005275_inline1.png\"/> <jats:tex-math>$J_T(x)$</jats:tex-math> </jats:alternatives> </jats:inline-formula> (non-dimensionalized by its inlet value) as a function of the streamwise distance <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" mime-subtype=\"png\" xlink:href=\"S0022112024005275_inline2.png\"/> <jats:tex-math>$x$</jats:tex-math> </jats:alternatives> </jats:inline-formula>, we prove that <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" mime-subtype=\"png\" xlink:href=\"S0022112024005275_inline3.png\"/> <jats:tex-math>$J_T(x) >1$</jats:tex-math> </jats:alternatives> </jats:inline-formula> – a condition that all TPJs must satisfy; surprisingly, prior theories and most experiments do not satisfy this condition. This motivated us to apply Lie symmetry analysis with translational and dilatational transformations of the modified equations (incorporating <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" mime-subtype=\"png\" xlink:href=\"S0022112024005275_inline4.png\"/> <jats:tex-math>$J_T>1$</jats:tex-math> </jats:alternatives> </jats:inline-formula>), which yields scaling laws for key jet measures: the mean streamwise and transverse velocities <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" mime-subtype=\"png\" xlink:href=\"S0022112024005275_inline5.png\"/> <jats:tex-math>$U(x,y)$</jats:tex-math> </jats:alternatives> </jats:inline-formula> and <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" mime-subtype=\"png\" xlink:href=\"S0022112024005275_inline6.png\"/> <jats:tex-math>$V(x,y)$</jats:tex-math> </jats:alternatives> </jats:inline-formula>, the turbulence intensities, the Reynolds shear stress <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" mime-subtype=\"png\" xlink:href=\"S0022112024005275_inline8.png\"/> <jats:tex-math>$-rho ,overline {u'v'}(x,y)$</jats:tex-math> </jats:alternatives> </jats:inline-formula>, the mean pressure <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org","PeriodicalId":15853,"journal":{"name":"Journal of Fluid Mechanics","volume":"46 1","pages":""},"PeriodicalIF":3.7,"publicationDate":"2024-08-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142226606","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
扫码关注我们
求助内容：
应助结果提醒方式：