Fluids最新文献

{"title":"Effect of Schmidt Number on Forced Isotropic Turbulence with Passive Scalars","authors":"Paolo Orlandi, Sergio Pirozzoli","doi":"10.3390/fluids8090248","DOIUrl":"https://doi.org/10.3390/fluids8090248","url":null,"abstract":"Traditionally, Fourier spectra have been employed to gain a deeper understanding of turbulence flow structures. The investigation of isotropic forced turbulence with passive scalars offers a straightforward means to examine the disparities between velocity and passive scalar spectra. This flow configuration has been extensively studied in the past, encompassing a range of Reynolds and Schmidt numbers. In this present study, direct numerical simulations (DNS) of this flow are conducted at sufficiently high Reynolds numbers, enabling the formation of a wide inertial range. The primary focus of this investigation is to quantitatively assess the variations in scalar spectra with the Schmidt number (Sc). The spectra exhibit a transition from a k−5/3 scaling for low Sc to a k−4/3 scaling for high Sc. The emergence of the latter power law becomes evident at Sc = 2, with its width expanding as Sc increases. To gain further insights into the underlying flow structures, a statistical analysis is performed by evaluating quantities aligned with the principal axes of the strain field. The study reveals that enstrophy is primarily influenced by the vorticity aligned with the intermediate principal strain axis, while the scalar gradient variance is predominantly controlled by the compressive strain. To provide a clearer understanding of the differences between enstrophy and scalar gradient variance, joint probability density functions (PDFs) and visualizations of the budget terms for both quantities are presented. These visualizations serve to elucidate the distinctions between the two and offer insights into their respective behaviors.","PeriodicalId":12397,"journal":{"name":"Fluids","volume":"363 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-09-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135885896","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Hydrodynamics of Supersonic Steam Jets Injected into Cross-Flowing Water","authors":"Hassan Ali Ghazwani, Khairuddin Sanaullah, Afrasyab Khan","doi":"10.3390/fluids8090250","DOIUrl":"https://doi.org/10.3390/fluids8090250","url":null,"abstract":"High-speed gas/vapour jets injected into a cross-moving sonic liquid signifies a vital phenomenon which bears useful applications in environmental and energy processes. In the present experimental study, a pulsating jet of supersonic steam was injected into cross-flowing water. Circulation zones of opposite vorticity owing to the interaction between the steam jet and cross-water flow were found. However, a large circulation appeared in front of the nozzle exit. Also, most small circulation regions were observed at higher water-flow rates (>2 m3/s). Among the prime mixing variables (i.e., turbulence kinetic energy (TKE) and Reynolds shear stress (RSS)), the RSS estimations backed a small diffusive phenomenon within a region far from the nozzle exit. Further information extracted from the PIV images indicated the existence of Kelvin–Helmholtz (KH) instabilities. The counter-rotating vortex pairs (CVPs) appeared to be significant in the region close to the nozzle exit, and they exhibited leeward side folds. Moreover, the effects of the operating conditions on the pressure recovery and mixing efficiency as well as the penetration and the separation height were evaluated to determine the optimisation of the phenomenon. By applying extreme difference analysis, the mixing efficiency was found as the most influential parameter.","PeriodicalId":12397,"journal":{"name":"Fluids","volume":"11 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-09-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135885297","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A Vaporization Model for Continuous Surface Force Approaches and Subcooled Configurations","authors":"Charles Brissot, Léa Cailly-Brandstäter, Elie Hachem, Rudy Valette","doi":"10.3390/fluids8080233","DOIUrl":"https://doi.org/10.3390/fluids8080233","url":null,"abstract":"The integration of phase change phenomena through an interface is a numerical challenge that requires proper attention. Solutions to properly ensure mass and energy conservation were developed for finite difference and finite volume methods, but not for Finite Element methods. We propose a Finite Element phase change model based on an Eulerian framework with a Continuous Surface Force (CSF) approach. It handles both momentum and energy conservation at the interface for anisotropic meshes in a light an efficient way. To do so, a model based on the Level Set method is developed. A thick interface is considered to fit with the CSF approach. To properly compute the energy conservation, heat fluxes are extended through this interface thanks to the resolution of a transport equation. A dedicated pseudo compressible Navier–Stokes solver is added to compute velocity jumps with a source term at the interface in the velocity divergence equation. Several 1D and 2D benchmarks are considered with increasing complexity to highlight the performances of each feature of the framework. This stresses the capacity of the model to properly tackle phase change problems.","PeriodicalId":12397,"journal":{"name":"Fluids","volume":"59 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-08-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135937369","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Calibration and Verification of Operation Parameters for an Array of Vectrino Profilers Configured for Turbulent Flow Field Measurement around Bridge Piers—Part II","authors":"Gordon Gilja, Robert Fliszar, Antonija Harasti, Manousos Valyrakis","doi":"10.3390/fluids8070199","DOIUrl":"https://doi.org/10.3390/fluids8070199","url":null,"abstract":"High-frequency velocimeters used for flow measurements during laboratory experiments allow the user to select the range for several operation parameters to set up the instrument for optimal velocity measurement. The discrepancies between velocity measurements collected with different instrument configurations can be significant, depending on the flume bed configuration and boundary conditions. The aim of this paper is to quantify the differences in flow velocity profiles measured with Acoustic Doppler Velocimeter Profilers (ADVPs) configured using a combination of profiling parameters: Ping Algorithm (PA), Transmit Pulse Size (TPS), and Cell Size (CS). Whereas in Part I of this research, the goal was to identify the optimal probe configuration for downstream measurement of the complex hydraulic structure (pier protected with riprap) based on a match of the flow rate with measurements from other instruments, in this paper, effect of distinct probe configuration on velocity profile and turbulent kinetic energy (TKE) is demonstrated. Differences between ADVPs’ configurations were analyzed through sensitivity analysis with the intention to calculate and compare any discrepancies in the velocity measurements for all the three measured velocity components: streamwise u, spanwise v and vertical w collected on two characteristic flume cross-sections. The results show that each parameter change has a significant effect on the measured values of each velocity component when compared to the Target Configuration (TC). The largest root-mean-square-error (RMSE) is observed when TPS is changed, followed by CS and PA. Absolute RMSE calculated for TPS change from 4 mm to 1 mm is, on average, 6.30 cm/s, 0.90 cm/s, and 0.82 cm/s for velocity components u, v and w, respectively. Absolute RMSE calculated for CS change from 1 mm to 4 mm is, on average, 4.49 cm/s, 0.88 cm/s, and 0.71 cm/s for velocity components u, v and w, respectively. Absolute RMSE calculated for PA change from Adaptive to Max interval is, on average, 4.04 cm/s, 0.63 cm/s, and 0.68 cm/s for velocity components u, v and w, respectively. For a change in all parameters, RMSE is greater for the cross-section downstream of the pier than for the approach cross-section: on average, 90%, 57% and 54% for a change in the PA, TPS, and CS, respectively.","PeriodicalId":12397,"journal":{"name":"Fluids","volume":"11 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-06-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135099974","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Transition to Equilibrium and Coherent Structure in Ideal MHD Turbulence, Part 2","authors":"John V. Shebalin","doi":"10.3390/fluids8060181","DOIUrl":"https://doi.org/10.3390/fluids8060181","url":null,"abstract":"We continue our study of the transition of ideal, homogeneous, incompressible, magnetohydrodynamic (MHD) turbulence from non-equilibrium initial conditions to equilibrium using long-time numerical simulations on a 1283 periodic grid. A Fourier spectral transform method is used to numerically integrate the dynamical equations forward in time. The six runs that previously went to near equilibrium are here extended into equilibrium. As before, we neglect dissipation as we are primarily concerned with behavior at the largest scale where this behavior has been shown to be essentially the same for ideal and real (forced and dissipative) MHD turbulence. These six runs have various combinations of imposed rotation and mean magnetic field and represent the five cases of ideal, homogeneous, incompressible, and MHD turbulence: Case I (Run 1), with no rotation or mean field; Case II (Runs 2a and 2b), where only rotation is imposed; Case III (Run 3), which has only a mean magnetic field; Case IV (Run 4), where rotation vector and mean magnetic field direction are aligned; and Case V (Run 5), which has non-aligned rotation vector and mean field directions. Statistical mechanics predicts that dynamic Fourier coefficients are zero-mean random variables, but largest-scale coherent magnetic structures emerge and manifest themselves as Fourier coefficients with very large, quasi-steady, mean values compared to their standard deviations, i.e., there is ‘broken ergodicity.’ These magnetic coherent structures appeared in all cases during transition to near equilibrium. Here, we report that, as the runs were continued, these coherent structures remained quasi-steady and energetic only in Cases I and II, while Case IV maintained its coherent structure but at comparatively low energy. The coherent structures that appeared in transition in Cases III and V were seen to collapse as their associated runs extended into equilibrium. The creation of largest-scale, coherent magnetic structure appears to be a dynamo process inherent in ideal MHD turbulence, particularly in Cases I and II, i.e., those cases most pertinent to planets and stars. Furthermore, the statistical theory of ideal MHD turbulence has proven to apply at the largest scale, even when dissipation and forcing are included. This, along with the discovery and explanation of dynamically broken ergodicity, is essentially a solution to the ‘dynamo problem’.","PeriodicalId":12397,"journal":{"name":"Fluids","volume":"77 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-06-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135961472","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Numerical Modelling of Water Flashing at Sub-Atmopsheric Pressure with a Multi-Regime Approach","authors":"Clément Loiseau, Stéphane Mimouni, Didier Colmont, Stéphane Vincent","doi":"10.3390/fluids8050143","DOIUrl":"https://doi.org/10.3390/fluids8050143","url":null,"abstract":"The CFD numerical study of the flash boiling phenomenon of a water film was conducted using an Euler–Euler method, and compared to the experiments on the flashing of a water film. The water film is initially heated at temperatures ranging from 34 to 74 ∘C (frim 1 to 41 ∘C superheat), and the pressure is decreased from 1 bar to 50 mbar during the experiments. This paper shows that the experiments could not be correctly modelled by a simple liquid/bubble model because of the overestimation of the drag force above the water film (in the gas/droplet region). The generalised large interface model (GLIM), however, a multi-regime approach implemented in the version 7.0 of the neptune_cfd software, is able to differentiate the water film, where liquid/bubble interactions are predominant from the gas region where gas/droplet interactions are predominant, and gives nice qualitative results. Finally, this paper shows that the interfacial heat transfer model of Berne for superheated liquids could accurately predict the evolution of the water temperature over time.","PeriodicalId":12397,"journal":{"name":"Fluids","volume":"41 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-04-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"136000472","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Two-Fluid Large-Eddy Simulation of Two-Phase Flow in Air-Sparged Hydrocyclone","authors":"Mustafa Bukhari, Hassan Fayed, Saad Ragab","doi":"10.3390/fluids8050139","DOIUrl":"https://doi.org/10.3390/fluids8050139","url":null,"abstract":"The two-fluid (Euler–Euler) model and large-eddy simulation are used to compute the turbulent two-phase flow of air and water in a cyclonic flotation device known as an Air-Sparged Hydrocyclone (ASH). In the operation of ASH, air is injected through a porous cylindrical wall. The study considers a 48 mm diameter hydrocyclone and uses a block-structured fine mesh of 10.5 million hexagonal elements. The air-to-water injection ratio is 4, and a uniform air bubble diameter of 0.5 mm is specified. The flow field in ASH was investigated for the inlet flow rate of water of 30.6 L/min at different values of underflow exit pressure. The current simulations quantify the effects of the underflow exit pressure on the split ratio and the overall flow physics in ASH, including the distribution of the air volume fraction, water axial velocity, tangential velocity, and swirling-layer thickness. The loci of zero-axial velocity surfaces were determined for different exit pressures. The water split ratio through the overflow opening varies with underflow exit pressure as 6%, 8%, 16%, and 26% for 3, 4, 5, and 6 kPa, respectively. These results indicate that regulating the pressure at the underflow exit can be used to optimize the ASH’s performance. Turbulent energy spectra in different regions of the hydrocyclone were analyzed. Small-scale turbulence spectra at near-wall points exhibit f−4 law, where f is frequency. Whereas for points at the air-column interface, the energy spectra show an inertial subrange f−5/3 followed by a dissipative range of f−7 law.","PeriodicalId":12397,"journal":{"name":"Fluids","volume":"95 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-04-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135068187","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"An Approximate Analytical Model of a Jet Flow with Mach Reflection and Pulsed Energy Supply at the Main Shock","authors":"Mikhail V. Chernyshov, Karina E. Savelova","doi":"10.3390/fluids8040132","DOIUrl":"https://doi.org/10.3390/fluids8040132","url":null,"abstract":"The supersonic flow of a reactive gas mixture with Mach reflection of oblique shocks and pulsed energy supply at the Mach stem is considered within the framework of the Chapman–Jouguet theory. An approximate analytical model is proposed that quickly determines the shape and size of the shock-wave structure as well as the flow parameters in various flow regions. As an example of the application of the proposed analytical model, the “first barrel” of a highly overexpanded jet flow of an air-methane mixture with a high supersonic velocity, is studied. Flows of hydrogen–air and hydrogen–oxygen mixtures were also considered for comparison with preceding numerical results. The height of the triple point of the Mach reflection is determined in the presence of a change in the chemical composition of the mixture and an isobaric pulsed energy supply at the main shock.","PeriodicalId":12397,"journal":{"name":"Fluids","volume":"140 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-04-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135034176","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Vortex Formation Times in the Glottal Jet, Measured in a Scaled-Up Model.","authors":"Michael Krane","doi":"10.3390/fluids6110412","DOIUrl":"https://doi.org/10.3390/fluids6110412","url":null,"abstract":"<p><p>In this paper, the timing of vortex formation on the glottal jet is studied using previously published velocity measurements of flow through a scaled-up model of the human vocal folds. The relative timing of the pulsatile glottal jet and the instability vortices are acoustically important since they determine the harmonic and broadband content of the voice signal. Glottis exit jet velocity time series were extracted from time-resolved planar DPIV measurements. These measurements were acquired at four glottal flow speeds (<i>u</i> _SS = 16.1-38 cm/s) and four glottis open times (<i>T</i> _o = 5.67-23.7 s), providing a Reynolds number range <i>Re</i> = 4100-9700 and reduced vibration frequency <i>f</i>* = 0.01-0.06. Exit velocity waveforms showed temporal behavior on two time scales, one that correlates to the period of vibration and another characterized by short, sharp velocity peaks (which correlate to the passage of instability vortices through the glottis exit plane). The vortex formation time, estimated by computing the time difference between subsequent peaks, was shown to be not well-correlated from one vibration cycle to the next. The principal finding is that vortex formation time depends not only on cycle phase, but varies strongly with reduced frequency of vibration. In all cases, a strong high-frequency burst of vortex motion occurs near the end of the cycle, consistent with perceptual studies using synthesized speech.</p>","PeriodicalId":12397,"journal":{"name":"Fluids","volume":"6 11","pages":""},"PeriodicalIF":1.9,"publicationDate":"2021-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC8627194/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"39941544","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Nanoparticle Delivery in Prostate Tumors Implanted in Mice Facilitated by Either Local or Whole-Body Heating.","authors":"Qimei Gu,&nbsp;Lance Dockery,&nbsp;Marie-Christine Daniel,&nbsp;Charles J Bieberich,&nbsp;Ronghui Ma,&nbsp;Liang Zhu","doi":"10.3390/fluids6080272","DOIUrl":"https://doi.org/10.3390/fluids6080272","url":null,"abstract":"<p><p>This work discusses in vivo experiments that were performed to evaluate whether local or whole-body heating to 40 °C reduced interstitial fluid pressures (IFPs) and enhanced nanoparticle delivery to subcutaneous PC3 human prostate cancer xenograft tumors in mice. After heating, 0.2 mL of a previously developed nanofluid containing gold nanoparticles (10 mg Au/mL) was injected via the tail vein. The induced whole-body hyperthermia led to increases in tumor and mouse body blood perfusion rates of more than 50% and 25%, respectively, while the increases were much smaller in the local heating group. In the whole-body hyperthermia groups, the IFP reduction from the baseline at the tumor center immediately after heating was found to be statistically significant when compared to the control group. The 1 h of local heating group showed IFP reductions at the tumor center, while the IFPs increased in the periphery of the tumor. The intratumoral gold nanoparticle accumulation was quantified using inductively coupled plasma mass spectrometry (ICP-MS). Compared to the control group, 1 h or 4 h of experiencing whole-body hyperthermia resulted in an average increase of 51% or 67% in the gold deposition in tumors, respectively. In the 1 h of local heating group, the increase in the gold deposition was 34%. Our results suggest that 1 h of mild whole-body hyperthermia may be a cost-effective and readily implementable strategy for facilitating nanoparticle delivery to PC3 tumors in mice.</p>","PeriodicalId":12397,"journal":{"name":"Fluids","volume":"6 8","pages":""},"PeriodicalIF":1.9,"publicationDate":"2021-08-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.3390/fluids6080272","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"39519132","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}