Next Issue
Volume 7, March
Previous Issue
Volume 7, January
 
 

Fluids, Volume 7, Issue 2 (February 2022) – 47 articles

Cover Story (view full-size image): Numerical simulations of contaminated spherical drops falling through a stagnant liquid at low Reynolds numbers are carried out. The results are used to describe the behavior of the surfactant concentrations and to understand the surfactant effects on the fluid motions in detail, e.g., peaks appear in the predicted Marangoni stresses near the stagnant-cap angle, which causes similar peaks in the pressure distribution, and high-pressure spots prevent fluid motion along the interface, which results in the formation of the stagnant-cap region and the attenuation of the tangential velocity in the continuous phase. View this paper.
  • Issues are regarded as officially published after their release is announced to the table of contents alert mailing list.
  • You may sign up for e-mail alerts to receive table of contents of newly released issues.
  • PDF is the official format for papers published in both, html and pdf forms. To view the papers in pdf format, click on the "PDF Full-text" link, and use the free Adobe Reader to open them.
Order results
Result details
Section
Select all
Export citation of selected articles as:
15 pages, 16065 KiB  
Article
A Modular Wave Energy Converter for Observational and Navigational Buoys
by Nicholas Vella, Jamie Foley, James Sloat, Alexander Sandoval, Leonardo D’Attile and Masoud Masoumi
Fluids 2022, 7(2), 88; https://doi.org/10.3390/fluids7020088 - 21 Feb 2022
Cited by 9 | Viewed by 5067
Abstract
More than 80% of the ocean is not fully mapped or even observed, even though it covers over 70% of our planet’s surface. One of the primary challenges for ocean observation and monitoring is the required power for exploration and monitoring systems, which [...] Read more.
More than 80% of the ocean is not fully mapped or even observed, even though it covers over 70% of our planet’s surface. One of the primary challenges for ocean observation and monitoring is the required power for exploration and monitoring systems, which often operate in remote areas of the ocean. This work addresses the design and development of an ocean wave energy converter that can be installed on observational buoys to provide enough power for sensors, cameras, data acquisition and recording, as well as data transfer units. The initial simulations of the prototype indicate that this system can produce up to 3.7–3.85 watts of power on average, with greater than 12 watts of maximum power in two selected sites in California and Hawaii. The proposed system is simple and low-cost. Further, multiple energy converters can be installed on one buoy to address higher power needs. Full article
(This article belongs to the Special Issue Wind and Wave Renewable Energy Systems, Volume II)
Show Figures

Graphical abstract

15 pages, 1375 KiB  
Article
General Relativistic Magnetohydrodynamics Mean-Field Dynamos
by Luca Del Zanna , Niccolò Tomei, Kevin Franceschetti, Matteo Bugli and Niccolò Bucciantini
Fluids 2022, 7(2), 87; https://doi.org/10.3390/fluids7020087 - 21 Feb 2022
Cited by 7 | Viewed by 2630
Abstract
Large-scale, ordered magnetic fields in several astrophysical sources are supposed to be originated, and maintained against dissipation, by the combined amplifying action of rotation and small-scale turbulence. For instance, in the solar interior, the so-called αΩ mean-field dynamo is known to [...] Read more.
Large-scale, ordered magnetic fields in several astrophysical sources are supposed to be originated, and maintained against dissipation, by the combined amplifying action of rotation and small-scale turbulence. For instance, in the solar interior, the so-called αΩ mean-field dynamo is known to be responsible for the observed 22-years magnetic cycle. Similar mechanisms could operate in more extreme environments, like proto neutron stars and accretion disks around black holes, for which the physical modelling needs to be translated from the regime of magnetohydrodynamics (MHD) and Newtonian gravity to that of a plasma in a general relativistic curved spacetime (GRMHD). Here we review the theory behind the mean field dynamo in GRMHD, the strategies for the implementation of the relevant equations in numerical conservative schemes, and we show the most important applications to the mentioned astrophysical compact objects obtained by our group in Florence. We also present novel results, such as three-dimensional GRMHD simulations of accretion disks with dynamo and the application of our dynamo model to a super massive neutron star, remnant of a binary neutron star merger as obtained from full numerical relativity simulations. Full article
Show Figures

Figure 1

14 pages, 3237 KiB  
Article
On the Computational Modeling of Inclined Brine Discharges
by Ilias G. Papakonstantis and Panos N. Papanicolaou
Fluids 2022, 7(2), 86; https://doi.org/10.3390/fluids7020086 - 20 Feb 2022
Cited by 3 | Viewed by 2437
Abstract
In this paper, five computational approaches are used to model bulk flow parameters of inclined round negatively buoyant jets. More specifically, an integral model employing Gaussian distributions for velocity and apparent acceleration of gravity, proposed in earlier study, is implemented with two different [...] Read more.
In this paper, five computational approaches are used to model bulk flow parameters of inclined round negatively buoyant jets. More specifically, an integral model employing Gaussian distributions for velocity and apparent acceleration of gravity, proposed in earlier study, is implemented with two different entrainment formulae. The remaining three computational approaches include an integral model known as EMA, which takes into consideration the fluid detachment occurring in the inner side of the flow near the terminal height, the widely known commercial model Corjet and analytical solutions that were proposed in a previous study. Predictions are provided for the maximum centerline height and its horizontal position, the terminal height of the upper jet boundary, the horizontal distance to the points where the jet centerline and the upper jet boundary return to the source level, the centerline dilution at the maximum height and the centerline dilution at the return point. Detailed comparisons are made in dimensionless form between the estimations provided by the models and a wide range of experimental data for discharge angles between 15° and 90°. Conclusions are drawn regarding the performance of the five computational approaches. Full article
(This article belongs to the Collection Feature Paper for Mathematical and Computational Fluid Mechanics)
Show Figures

Graphical abstract

22 pages, 3972 KiB  
Article
Numerical Investigation and Fluid-Structure Interaction (FSI) Analysis on a Double-Element Simplified Formula One (F1) Composite Wing in the Presence of Ground Effect
by Chris Sungkyun Bang, Zeeshan A. Rana, László Könözsy, Veronica Marchante Rodriguez and Clive Temple
Fluids 2022, 7(2), 85; https://doi.org/10.3390/fluids7020085 - 19 Feb 2022
Cited by 5 | Viewed by 6081
Abstract
This research paper focuses on a novel coupling of the aerodynamic and structural behaviour of a double-element composite front wing of a Formula One (F1) vehicle, which was simulated and studied for the first time here. To achieve this goal, a modified two-way [...] Read more.
This research paper focuses on a novel coupling of the aerodynamic and structural behaviour of a double-element composite front wing of a Formula One (F1) vehicle, which was simulated and studied for the first time here. To achieve this goal, a modified two-way coupling method was employed in the context of high performance computing (HPC) to simulate a steady-state fluid-structure interaction (FSI) configuration using the ANSYS software package. The front wing plays a key role in generating aerodynamic forces and controlling the fresh airflow to maximise the aerodynamic performance of an F1 car. Therefore, the composite front wing becomes deflected under aerodynamic loading conditions due to its elastic behaviour which can lead to changes in the flow field and the aerodynamic performance of the wing. To reduce the uncertainty of the simulations, a grid sensitivity study and the assessment of different engineering turbulence models were carried out. The practical contribution of our investigations is the quantification of the coupled effect of the aerodynamic and structural performance of the wing and an understanding of the influence of ride heights on the ground effect. It was found that the obtained numerical surface pressure distributions, the aerodynamic forces, and the wake profiles show an accurate agreement with experimental data taken from the literature. Full article
(This article belongs to the Special Issue Fluid Structure Interaction: Methods and Applications)
Show Figures

Graphical abstract

14 pages, 2711 KiB  
Article
Interannual Bottom-Intensified Current Thickening Observed on the Continental Slope Off the Southeastern Coast of Hokkaido, Japan
by Akira Nagano, Takuya Hasegawa, Keisuke Ariyoshi and Hiroyuki Matsumoto
Fluids 2022, 7(2), 84; https://doi.org/10.3390/fluids7020084 - 19 Feb 2022
Cited by 4 | Viewed by 2450
Abstract
By rotary empirical orthogonal function and coastal-trapped wave mode analyses, we analyzed current velocity data, collected from 2001 to 2016. The data were obtained by an acoustic Doppler current profiler, deployed upward at a location of 41°39.909′ N, 144°20.695′ E, on a 2630-m [...] Read more.
By rotary empirical orthogonal function and coastal-trapped wave mode analyses, we analyzed current velocity data, collected from 2001 to 2016. The data were obtained by an acoustic Doppler current profiler, deployed upward at a location of 41°39.909′ N, 144°20.695′ E, on a 2630-m deep continental slope seabed off the southeastern coast of Hokkaido, Japan. The results indicate that the current intensifies toward the bottom and is directed nearly toward the shore, reaching an average speed of ~2.5 cm s−1 just above the bottom. The thickness of the along-slope northward component of the bottom-intensified current varied within the range of 50–350 m. We found that the current thickness change was caused by oceanic barotropic disturbances, produced by the intensification of the Aleutian Low, largely related to the El Niño–Southern Oscillation and modified through the excitation of bottom-trapped modes of coastal-trapped waves. This finding improves the prediction accuracy of the the bottom-intensified current change, being beneficial for suspended sediment studies, construction and maintenance of marine structures, planning of deep drilling, and so on. Full article
(This article belongs to the Special Issue Boundary Layer Processes in Geophysical/Environmental Flows)
Show Figures

Figure 1

20 pages, 2349 KiB  
Article
Numerical Simulation of Departure from Nucleate Boiling in Rod Bundles under High-Pressure Conditions
by Sai Raja Gopal Vadlamudi and Arun K. Nayak
Fluids 2022, 7(2), 83; https://doi.org/10.3390/fluids7020083 - 18 Feb 2022
Cited by 4 | Viewed by 3290
Abstract
In subcooled boiling flows beyond a certain heat flux, heat transfer is hampered due to a phenomenon known as Departure from Nucleate Boiling (DNB). Conducting DNB experiments at one-to-one nuclear reactor operating conditions is highly challenging and expensive. Another alternative approach is to [...] Read more.
In subcooled boiling flows beyond a certain heat flux, heat transfer is hampered due to a phenomenon known as Departure from Nucleate Boiling (DNB). Conducting DNB experiments at one-to-one nuclear reactor operating conditions is highly challenging and expensive. Another alternative approach is to use Look-up table data. However, its applicability is limited due to its dependence on rod bundle correction factors. In the present investigation, a state-of-the-art Eulerian-Eulerian two-fluid model coupled with an extended heat flux partitioning model is used to predict DNB in tubes and rod bundles with square and hexagonal lattices (relevant to Pressurized Water Reactors). In this approach, bubble departure characteristics are modeled using semi-mechanistic models based on force balance analysis. The predicted DNB values are compared with experimental and Look-up table data and found out to be within 1.8% to 20%. Full article
(This article belongs to the Special Issue Advances in Multiphase Flow Science and Technology)
Show Figures

Figure 1

12 pages, 1112 KiB  
Article
Note on the Bulk Estimate of the Energy Dissipation Rate in the Oceanic Bottom Boundary Layer
by Xiaozhou Ruan
Fluids 2022, 7(2), 82; https://doi.org/10.3390/fluids7020082 - 18 Feb 2022
Cited by 2 | Viewed by 2537
Abstract
The dissipation of the kinetic energy (KE) associated with oceanic flows is believed to occur primarily in the oceanic bottom boundary layer (BBL), where bottom drag converts the KE from mean flows to heat loss through irreversible mixing at molecular scales. Due to [...] Read more.
The dissipation of the kinetic energy (KE) associated with oceanic flows is believed to occur primarily in the oceanic bottom boundary layer (BBL), where bottom drag converts the KE from mean flows to heat loss through irreversible mixing at molecular scales. Due to the practical difficulties associated with direct observations on small-scale turbulence close to the seafloor, most up-to-date estimates on bottom drag rely on a simple bulk formula (CdU3) proposed by G.I. Taylor that relates the integrated BBL dissipation rate to a drag coefficient (Cd) as well as a flow magnitude outside of the BBL (U). Using output from several turbulence-resolving direct numerical simulations, it is shown that the true BBL-integrated dissipation rate is approximately 90% of that estimated using the classic bulk formula, applied here to the simplest scenario where a mean flow is present over a flat and hydrodynamically smooth bottom. It is further argued that Taylor’s formula only provides an upper bound estimate and should be applied with caution in the future quantification of BBL dissipation; the performance of the bulk formula depends on the distribution of velocity and shear stress near the bottom, which, in the real ocean, could be disrupted by bottom roughness. Full article
(This article belongs to the Special Issue Boundary Layer Processes in Geophysical/Environmental Flows)
Show Figures

Figure 1

30 pages, 4327 KiB  
Article
A Hybrid Non-Linear Unsteady Vortex Lattice-Vortex Particle Method for Rotor Blades Aerodynamic Simulations
by Vincent Proulx-Cabana, Minh Tuan Nguyen, Sebastien Prothin, Guilhem Michon and Eric Laurendeau
Fluids 2022, 7(2), 81; https://doi.org/10.3390/fluids7020081 - 16 Feb 2022
Cited by 11 | Viewed by 4092
Abstract
This study presents a hybrid non-linear unsteady vortex lattice method-vortex particle method (NL UVLM-VPM) to investigate the aerodynamics of rotor blades hovering in and out of ground effect. The method is of interest for the fast aerodynamic prediction of helicopter and smaller rotor [...] Read more.
This study presents a hybrid non-linear unsteady vortex lattice method-vortex particle method (NL UVLM-VPM) to investigate the aerodynamics of rotor blades hovering in and out of ground effect. The method is of interest for the fast aerodynamic prediction of helicopter and smaller rotor blades. UVLM models the vorticity along the rotor blades and near field wakes with panels that are then converted into their equivalent vortex particle representations. The standard Vreman subgrid scale model is incorporated in the context of a large eddy simulation for mesh-free VPM to stabilize the wake development via particle strength exchange (PSE). The computation of the pairwise interactions in the VPM are accelerated using the fast-multipole method. Non-linear UVLM is achieved with a low computational cost viscous-inviscid alpha coupling algorithm through a stripwise 2D Reynolds-averaged Navier–Stokes (RANS) or empirical database. The aerodynamics of the scaled S76 rotor blades in and out of ground effect from the hover prediction workshop is investigated with the proposed algorithm. The results are validated with experimental data and various high-fidelity codes. Full article
Show Figures

Figure 1

29 pages, 20806 KiB  
Article
Modeling Transport of SARS-CoV-2 Inside a Charlotte Area Transit System (CATS) Bus
by Matthew Goodson, Jeffrey Feaster, Andy Jones, Gregory McGowan, Lucas Agricola, William Timms and Mesbah Uddin
Fluids 2022, 7(2), 80; https://doi.org/10.3390/fluids7020080 - 16 Feb 2022
Cited by 4 | Viewed by 2710
Abstract
We present in this paper a model of the transport of human respiratory particles on a Charlotte Area Transit System (CATS) bus to examine the efficacy of interventions to limit exposure to SARS-CoV-2, the virus that causes COVID-19. The methods discussed here utilize [...] Read more.
We present in this paper a model of the transport of human respiratory particles on a Charlotte Area Transit System (CATS) bus to examine the efficacy of interventions to limit exposure to SARS-CoV-2, the virus that causes COVID-19. The methods discussed here utilize a commercial Navier–Stokes flow solver, RavenCFD, using a massively parallel supercomputer to model the flow of air through the bus under varying conditions, such as windows being open or the HVAC flow settings. Lagrangian particles are injected into the RavenCFD predicted flow fields to simulate the respiratory droplets from speaking, coughing, or sneezing. These particles are then traced over time and space until they interact with a surface or are removed via the HVAC system. Finally, a volumetric Viral Mean Exposure Time (VMET) is computed to quantify the risk of exposure to the SARS-CoV-2 under various environmental and occupancy scenarios. Comparing the VMET under varying conditions should help identify viable methods to reduce the risk of viral exposure of CATS bus passengers during the COVID-19 pandemic. Full article
(This article belongs to the Special Issue Aerodynamics and Aeroacoustics of Vehicles, Volume II)
Show Figures

Figure 1

12 pages, 15461 KiB  
Article
Rim Breakups of Impacting Drops on a Superhydrophobic Surface and a Superheated Surface
by Minori Shirota, Masaki Kato and Ai Ishio
Fluids 2022, 7(2), 79; https://doi.org/10.3390/fluids7020079 - 15 Feb 2022
Cited by 2 | Viewed by 2471
Abstract
The rim breakup of an impacting drop is experimentally investigated by comparing the impacts on superheated and superhydrophobic surfaces. The objective of the present study is to experimentally examine whether the Bo = 1 criteria holds for the rim breakups of drops impacting [...] Read more.
The rim breakup of an impacting drop is experimentally investigated by comparing the impacts on superheated and superhydrophobic surfaces. The objective of the present study is to experimentally examine whether the Bo = 1 criteria holds for the rim breakups of drops impacting on the surfaces. A transparent sapphire plate was heated to achieve the Leidenfrost impact, which enables us to observe with a high-speed camera from below. The characteristics of the rim breakup were evaluated quantitatively using a particle tracking velocimetry method for both the rim and the drops generated. As a result, we clarified that Bo of the rim increases in the spreading phase and marks the highest value of 0.5 on a superheated surface, which is smaller than that on a pillar, where Bo ≈ 1. On a superhydrophobic surface, the highest Bo was 1.2, which is smaller than that on a wettable solid surface, 2.5, but close to the value on a pillar. We also revealed that diameters of generated drops collapse on a master curve when plotted as a function of pinch-off time for both the impacts on superheated and superhydrophobic surfaces. Full article
(This article belongs to the Special Issue Advances in Multiphase Flow Science and Technology)
Show Figures

Figure 1

20 pages, 4427 KiB  
Article
Experimental Techniques against RANS Method in a Fully Developed Turbulent Pipe Flow: Evolution of Experimental and Computational Methods for the Study of Turbulence
by Gabriela Lopez-Santana, Andrew Kennaugh and Amir Keshmiri
Fluids 2022, 7(2), 78; https://doi.org/10.3390/fluids7020078 - 15 Feb 2022
Cited by 9 | Viewed by 3663
Abstract
Fully developed turbulent flow in a pipe was studied by considering experimental and computational methods. The aim of this work was to build on the legacy of the University of Manchester, which is widely regarded as the birthplace of turbulence due to the [...] Read more.
Fully developed turbulent flow in a pipe was studied by considering experimental and computational methods. The aim of this work was to build on the legacy of the University of Manchester, which is widely regarded as the birthplace of turbulence due to the pioneering work of the prominent academic Professor Osborne Reynolds (1842–1912), by capturing the evolution of fluid turbulence analysis tools over the last 100 years. A classical experimental apparatus was used to measure the mean velocity field and wall shear stress through four historical techniques: static pressure drop; mean square signals measured from a hot-wire; Preston tube; and Clauser plot. Computational Fluid Dynamics (CFD) was used to simulate the pipe flow, utilizing the Reynolds-averaged Navier–Stokes (RANS) method with different two-equation turbulence models. The performance of each approach was assessed to compare the experimental and computational methods. This comparison revealed that the numerical results produced a close agreement with the experiments. The finding shows that, in some cases, CFD simulations could be used as alternative or complementary methods to experimental techniques for analyzing fully developed turbulent pipe flow. Full article
Show Figures

Figure 1

16 pages, 3556 KiB  
Article
Application of an Integral Turbulence Model to Close the Model of an Anisotropic Porous Body as Applied to Rod Structures
by Maksim N. Vlasov and Igor G. Merinov
Fluids 2022, 7(2), 77; https://doi.org/10.3390/fluids7020077 - 14 Feb 2022
Cited by 6 | Viewed by 2284
Abstract
In practice, often devices are ordered rod structures consisting of a large number of rods. Heat exchangers, fuel assemblies of nuclear reactors, and their cores in the case of using caseless assemblies are examples of such devices. Simulation of heat and mass transfer [...] Read more.
In practice, often devices are ordered rod structures consisting of a large number of rods. Heat exchangers, fuel assemblies of nuclear reactors, and their cores in the case of using caseless assemblies are examples of such devices. Simulation of heat and mass transfer processes in such devices in porous-body approximation can significantly reduce the required resources compared to computational fluid dynamics (CFD) approaches. The paper describes an integral turbulence model developed for defining anisotropic model parameters of a porous body. The parameters of the integral turbulence model were determined by numerical simulations for assemblies of smooth rods, assemblies with spacer grids, and wire-wrapped fuel assemblies. The results of modeling the flow of a liquid metal coolant in an experimental fuel assembly with local blocking of its flow section in anisotropic porous-body approximation using an integral turbulence model are described. The possibility of using the model of an anisotropic porous body with the integral model of turbulence to describe thermal-hydraulic processes during fluid flow in rod structures is confirmed. Full article
(This article belongs to the Special Issue Stochastic Equations in Fluid Dynamics)
Show Figures

Figure 1

31 pages, 452 KiB  
Concept Paper
Reduced Numerical Modeling of Turbulent Flow with Fully Resolved Time Advancement. Part 1. Theory and Physical Interpretation
by Alan R. Kerstein
Fluids 2022, 7(2), 76; https://doi.org/10.3390/fluids7020076 - 13 Feb 2022
Cited by 6 | Viewed by 2265
Abstract
A multiscale modeling concept for numerical simulation of multiphysics turbulent flow utilizing map-based advection is described. The approach is outlined with emphasis on its theoretical foundations and physical interpretations in order to establish the context for subsequent presentation of the associated numerical algorithms [...] Read more.
A multiscale modeling concept for numerical simulation of multiphysics turbulent flow utilizing map-based advection is described. The approach is outlined with emphasis on its theoretical foundations and physical interpretations in order to establish the context for subsequent presentation of the associated numerical algorithms and the results of validation studies. The model formulation is a synthesis of existing methods, modified and extended in order to obtain a qualitatively new capability. The salient feature of the approach is that time advancement of the flow is fully resolved both spatially and temporally, albeit with modeled advancement processes restricted to one spatial dimension. This one-dimensional advancement is the basis of a bottom-up modeling approach in which three-dimensional space is discretized into under-resolved mesh cells, each of which contains an instantiation of the modeled one-dimensional advancement. Filtering is performed only to provide inputs to a pressure correction that enforces continuity and to obtain mesh-scale-filtered outputs if desired. The one-dimensional advancement, the pressure correction, and coupling of one-dimensional instantiations using a Lagrangian implementation of mesh-resolved volume fluxes is sufficient to advance the three-dimensional flow without time advancing coarse-grained equations, a feature that motivates the designation of the approach as autonomous microscale evolution (AME). In this sense, the one-dimensional treatment is not a closure because there are no unclosed terms to evaluate. However, the approach is additionally suitable for use as a subgrid-scale closure of existing large-eddy-simulation methods. The potential capabilities and limitations of both of these implementations of the approach are assessed conceptually and with reference to demonstrated capabilities of related methods. Full article
(This article belongs to the Special Issue Turbulent Flow)
Show Figures

Figure 1

15 pages, 13222 KiB  
Article
Efficiency and Aerodynamic Performance of Bristled Insect Wings Depending on Reynolds Number in Flapping Flight
by Felicity O’Callaghan, Amir Sarig, Gal Ribak and Fritz-Olaf Lehmann
Fluids 2022, 7(2), 75; https://doi.org/10.3390/fluids7020075 - 10 Feb 2022
Cited by 6 | Viewed by 3985
Abstract
Insect wings are generally constructed from veins and solid membranes. However, in the case of the smallest flying insects, the wing membrane is often replaced by hair-like bristles. In contrast to large insects, it is possible for both bristled and membranous wings to [...] Read more.
Insect wings are generally constructed from veins and solid membranes. However, in the case of the smallest flying insects, the wing membrane is often replaced by hair-like bristles. In contrast to large insects, it is possible for both bristled and membranous wings to be simultaneously present in small insect species. There is therefore a continuing debate about the advantages and disadvantages of bristled wings for flight. In this study, we experimentally tested bristled robotic wing models on their ability to generate vertical forces and scored aerodynamic efficiency at Reynolds numbers that are typical for flight in miniature insects. The tested wings ranged from a solid membrane to a few bristles. A generic lift-based wing kinematic pattern moved the wings around their root. The results show that the lift coefficients, power coefficients and Froude efficiency decreased with increasing bristle spacing. Skin friction significantly attenuates lift production, which may even result in negative coefficients at elevated bristle spacing and low Reynolds numbers. The experimental data confirm previous findings from numerical simulations. These had suggested that for small insects, flying with bristled instead of membranous wings involved less change in energetic costs than for large insects. In sum, our findings highlight the aerodynamic changes associated with bristled wing designs and are thus significant for assessing the biological fitness and dispersal of flying insects. Full article
(This article belongs to the Special Issue Advances in Biological Flows and Biomimetics, Volume II)
Show Figures

Figure 1

15 pages, 3862 KiB  
Article
Simulation of Natural Convection by Multirelaxation Time Lattice Boltzmann Method in a Triangular Enclosure
by Jamal Baliti, Youssef Elguennouni, Mohamed Hssikou and Mohammed Alaoui
Fluids 2022, 7(2), 74; https://doi.org/10.3390/fluids7020074 - 10 Feb 2022
Cited by 12 | Viewed by 2485
Abstract
The natural convection of incompressible flow confined within an enclosed right-angled triangular and isosceles cavity was investigated numerically using the multirelaxation time lattice Boltzmann method (MRT-LBM). According to the left and inclined walls thermal boundary conditions, two cases were considered in this study. [...] Read more.
The natural convection of incompressible flow confined within an enclosed right-angled triangular and isosceles cavity was investigated numerically using the multirelaxation time lattice Boltzmann method (MRT-LBM). According to the left and inclined walls thermal boundary conditions, two cases were considered in this study. In the first case, the inclined side of the enclosure was adiabatic, and the horizontal wall was heated, while the left one was kept at a cold temperature. However, the states of the left and inclined walls were interchanged in the second case. As the flow is only transported under the convection force, this study was carried out for the Rayleigh number ranging from Ra=103 to 106. The effects of the Rayleigh number on velocity and temperature profiles, streamlines, isotherms, and average Nusselt number were investigated. The position of cold and adiabatic walls had a great effect on the results. The results obtained are in good agreement with those of the literature and show the robustness of the MRT-LBM approach. In both cases, the heat-transfer rate increases with the increase in the Rayleigh number. Full article
(This article belongs to the Section Heat and Mass Transfer)
Show Figures

Figure 1

20 pages, 3422 KiB  
Article
The Capillary Waves’ Contribution to Wind-Wave Generation
by Harald Naeser
Fluids 2022, 7(2), 73; https://doi.org/10.3390/fluids7020073 - 10 Feb 2022
Cited by 2 | Viewed by 4460
Abstract
Published theories and observations have shown that dissipation of gravity waves implies frequency downshifting of wave energy. Hence, for wind-waves, the wind energy input to the highest frequencies is of special interest. Here it is shown that this input is vital, because the [...] Read more.
Published theories and observations have shown that dissipation of gravity waves implies frequency downshifting of wave energy. Hence, for wind-waves, the wind energy input to the highest frequencies is of special interest. Here it is shown that this input is vital, because the direct wind energy input obtained by the air-pressure’s work on most gravity waves is slightly less than what the waves need to grow. Further, the wind’s input of the angular momentum that waves need to grow is found to be absent at most gravity wave frequencies. The capillary waves that appear at the surface of the sea when the wind is blowing solve these problems. To demonstrate this, an extension of linear wave theory is established to study possibilities and limitations for transfer of energy and angular momentum from the wind to waves through these frequencies. The theory describes regular, gravity–capillary waves with constant amplitude under laminar conditions. It includes surface tensions, viscosity, gravity and a wind-generated shear current, and shows that these waves—contrary to most gravity waves—receive more energy from the wind than they dissipate and angular momentum they cannot keep. Hence, the problem of the missing input of energy and angular momentum from wind to gravity waves is solved by transfers through the capillary waves. This implies that capillary waves are vital to obtain growing gravity waves. Full article
(This article belongs to the Section Geophysical and Environmental Fluid Mechanics)
Show Figures

Figure 1

13 pages, 2081 KiB  
Article
Reducing Slipstream Velocities Experienced in Proximity to High-Speed Trains
by Jordan Ashley Dunlop and Mark Christopher Thompson
Fluids 2022, 7(2), 72; https://doi.org/10.3390/fluids7020072 - 9 Feb 2022
Cited by 1 | Viewed by 1941
Abstract
Slipstream wake structures generated by the passing of high-speed rail vehicles represent a hazard to passengers and workers in close proximity. In this article, the possibility of reducing peak slipstream velocities through the implementation of angled fins or swirling flow injection is assessed [...] Read more.
Slipstream wake structures generated by the passing of high-speed rail vehicles represent a hazard to passengers and workers in close proximity. In this article, the possibility of reducing peak slipstream velocities through the implementation of angled fins or swirling flow injection is assessed on the basis of improved delayed detached eddy simulations (IDDES). The key to improving slipstream velocities involves redirecting and/or reducing the internal energy, a pair of meandering counter-rotating vortex cores that are associated with large wake slipstream deviations. It is demonstrated that the danger imposed by slipstream wake structures, as measured by the induced velocity measures recorded at a series of test points adjacent to the passing train, could be significantly reduced, with decreases from 10% up to 23%, recorded across a range of sampling locations. The means by which these reductions are generated and the corresponding changes in the flow are also explored through analysis of the modified wakes. As such, these devices show promise at improving the operational safety of high-speed rail vehicles. Full article
(This article belongs to the Special Issue Aerodynamics of Road Vehicles and Trains)
Show Figures

Figure 1

11 pages, 2252 KiB  
Article
Effect of the Coastline Geometry on the Boundary Currents Intruding through the Gap
by Joseph Kuehl and Vitalii A. Sheremet
Fluids 2022, 7(2), 71; https://doi.org/10.3390/fluids7020071 - 8 Feb 2022
Cited by 6 | Viewed by 1920
Abstract
The problem of a geophysical western boundary current negotiating a gap in its supporting boundary is considered. For traditional straight, parallel gaps, such systems are known to exhibit two dominant states, gap penetrating and leaping, with the transitional dynamics between states displaying hysteresis. [...] Read more.
The problem of a geophysical western boundary current negotiating a gap in its supporting boundary is considered. For traditional straight, parallel gaps, such systems are known to exhibit two dominant states, gap penetrating and leaping, with the transitional dynamics between states displaying hysteresis. However, for more complex geometries, such as angled or offset gap configurations, the question of multiple states and hysteresis is unresolved. In such cases, the inertia of the western boundary current is oriented into the gap, hence the assumption that increased inertia promotes gap penetrating loop current states. Here we address the problem numerically in an idealized setting. It is found that despite the inertia of the current being directed into the gap, for large western boundary current transport values, leaping states will be present. That is, we show here that the presence of multiple states with hysteresis for gap-leaping western boundary current systems is robust to both angled and offset gap geometries. Full article
(This article belongs to the Special Issue Boundary Layer Processes in Geophysical/Environmental Flows)
Show Figures

Figure 1

12 pages, 3622 KiB  
Article
Viscoelastic Representation of the Operation of Sucker Rod Pumps
by Sheldon Wang, Lynn Rowlan, Abbey Henderson, Sean T. Aleman, Trent Creacy and Carrie Anne Taylor
Fluids 2022, 7(2), 70; https://doi.org/10.3390/fluids7020070 - 8 Feb 2022
Cited by 3 | Viewed by 2620
Abstract
The issues of leakage with respect to the clearance between the pump plunger outer surface and the pump barrel inner surface and other operational conditions have been revisited in this paper. Both Poiseuille flow rate due to the pressure difference and Couette flow [...] Read more.
The issues of leakage with respect to the clearance between the pump plunger outer surface and the pump barrel inner surface and other operational conditions have been revisited in this paper. Both Poiseuille flow rate due to the pressure difference and Couette flow rate due to the plunger motion have been considered. The purpose of this study is to explore the possibility of representing the entire downhole pump system with a simple viscoelastic model. We have explored both Kelvin and Maxwell viscoelastic models along with the dynamic behaviors of a mass point attached to the viscoelastic model. By using the time-dependent polished rod force measured with a dynamometer as the input to the viscoelastic models, we have obtained the displacement responses, which match closely with the actual measurements in experiments and operations. Further study and experiments have been planned and partially implemented in the McCoy School of Engineering at Midwestern State University, a member of the Texas Tech University System. Full article
(This article belongs to the Collection Complex Fluids)
Show Figures

Figure 1

14 pages, 2932 KiB  
Article
The Influence of Exit Nozzle Geometry on Sweeping Jet Actuator Performance
by Mobashera Alam and Kursat Kara
Fluids 2022, 7(2), 69; https://doi.org/10.3390/fluids7020069 - 8 Feb 2022
Cited by 8 | Viewed by 3366
Abstract
When pressurized with a fluid, the sweeping jet actuator (SWJA) emits a self-induced and self-sustained temporally continuous, but spatially oscillating bi-stable jet at the outlet. The SWJA adds up local momentum using the Coanda extension without any moving parts and, therefore, can be [...] Read more.
When pressurized with a fluid, the sweeping jet actuator (SWJA) emits a self-induced and self-sustained temporally continuous, but spatially oscillating bi-stable jet at the outlet. The SWJA adds up local momentum using the Coanda extension without any moving parts and, therefore, can be a promising tool for suppressing aerodynamic flow separation. However, the SWJA needs to be integrated into curved aerodynamic surfaces with an angle. The present study focuses on investigating the effects of various exit nozzle geometries on the flow field. The geometric parameters considered were the exit nozzle angle, diffuser arm length, and curvature. The working fluid was air, and the mass flow rate was 0.015 lb/s. A set of time-dependent flow fields was computed using a two-dimensional unsteady Reynolds-averaged Navier–Stokes (URANS) simulation. The time history of pressure was recorded inside the upper and lower feedback channels. The jet oscillation frequency was obtained by employing the fast Fourier transform (FFT) for all datasets. The results were compared against the baseline case and data available in the literature. The results showed that external geometric variations at the nozzle exit had a negligible impact on the oscillation frequency. However, there were notable effects on the pressure and velocity distribution in the flow field, indicating that the actuator had sensitivity towards the geometric variation of the exit nozzle—the wider the exit nozzle, the lower the downstream velocity. Notably, we observed that the mean velocity at the exit nozzle downstream for the curvature case was 40.3% higher than the reference SWJA. Full article
Show Figures

Graphical abstract

20 pages, 4500 KiB  
Article
Application of the Moment Method for Numerical Simulation of Homogeneous-Heterogeneous Condensation
by Igor E. Ivanov, Vladislav S. Nazarov and Igor A. Kryukov
Fluids 2022, 7(2), 68; https://doi.org/10.3390/fluids7020068 - 7 Feb 2022
Cited by 3 | Viewed by 2435
Abstract
The paper considers the numerical modeling of the processes of homogeneous and heterogeneous condensation and evaporation in multiphase flows using the method of moments. Nonstationary processes of gas dynamics and phase transitions in the two-dimensional plane and axisymmetric regions are described by a [...] Read more.
The paper considers the numerical modeling of the processes of homogeneous and heterogeneous condensation and evaporation in multiphase flows using the method of moments. Nonstationary processes of gas dynamics and phase transitions in the two-dimensional plane and axisymmetric regions are described by a general system of equations. The system of equations is expanded by adding two equations. One describes the evolution of the total mass fraction of the condensing substance; the other describes the evolution of the mass fraction of solid particles. An instant wetting model is used to model heterogeneous nucleation. The Gyarmathy model is used for the approximation of the average droplet growth rate. Heterogeneous condensation is modeled based on the distribution function of foreign impurities. An approach to calculating evaporation in the heterogeneous case is proposed. A comparison of the proposed models with a numerical experiment is given. Numerical simulation of homogeneous-heterogeneous condensation in a gas-dynamic ejector is carried out. Full article
Show Figures

Figure 1

15 pages, 13845 KiB  
Article
GIS and Wave Modeling for Establishing a Potential Area of Aquaculture—Case Study: Central Atlantic Part of the Moroccan Coast
by Mohamed Amine Taji, Atika Hilali, Hassan Rhinane, Antoine Mangin, Philippe Bryère, Abdelatif Orbi, Hassan Mabchour, Bendahhou Zourarah and Aïssa Benazzouz
Fluids 2022, 7(2), 67; https://doi.org/10.3390/fluids7020067 - 7 Feb 2022
Cited by 1 | Viewed by 3334
Abstract
Marine aquaculture has proliferated over the past decade, expanding into new, untapped open-water cultivation areas, such as lakes, rivers and deeper offshore environments, in response to increasing demand for seafood by consumers. However, to ensure sustainable development, it is necessary to minimize the [...] Read more.
Marine aquaculture has proliferated over the past decade, expanding into new, untapped open-water cultivation areas, such as lakes, rivers and deeper offshore environments, in response to increasing demand for seafood by consumers. However, to ensure sustainable development, it is necessary to minimize the impact of other ocean activities and the environment through science-based spatial planning. The choice of the primary site (physical carrying capacity) depends mainly on the aquaculture system, which varies around the world. However, the site is considered one of the factors (production, ecological and social) keys to any aquaculture operation, especially in the African continent. This choice affects both the success and sustainability of the products cultivated and the resolution of conflicts between different activities as well as the rational use of space. This study aims to identify suitable areas (primary site selection) for aquaculture in the Moroccan Atlantic continental shelf focused on the sub-area located between Cap Ghir 31.25° and Tarfaya 27.47°, based on the assessment of the dominant wave energy by implementing the hydrodynamical SWAN (Simulating Waves Nearshore) model dedicated for this kind of study. We derived the inputs for the SWAN model from WW3 (WAVEWATCH III model), which the AVISO data-products have extensively validated. The results show that, even if the Atlantic area is known for the agitation of its seas, there is the possibility of having adequate areas for aquaculture with an overall capacity that could extinguish the 389 ha in the study area if aquatic cultivation manages to exploit the offshore areas. At the level of the sub-zone belonging to the sous-Massa region (zone 1), the results show a strong coherence between the values of the surfaces estimated by the study and the actual values resulting from the development plan, with a value of 69 Ha for the first and 75 for the second, i.e., equal to 6 Ha, due to the geomorphology of the coast and natural coastal shelters, which play favorably on the environment for aquaculture development. These areas may attract the greed of investors, although they are in the process of being the subject of an aquaculture development plan. Full article
Show Figures

Figure 1

16 pages, 2169 KiB  
Review
Sublethal Damage to Erythrocytes during Blood Flow
by Mesude Avcı, Edgar A. O’Rear, Kylie M. Foster and Dimitrios V. Papavassiliou
Fluids 2022, 7(2), 66; https://doi.org/10.3390/fluids7020066 - 7 Feb 2022
Cited by 1 | Viewed by 4051
Abstract
Mechanical circulatory support (MCS) devices are designed to perform the functional needs of organs and to meet clinical hemocompability criteria. Critical complications have been reported with their long-term use such as thrombosis, anemia and gastrointestinal bleeding. Damage to red blood cells (RBCs), which [...] Read more.
Mechanical circulatory support (MCS) devices are designed to perform the functional needs of organs and to meet clinical hemocompability criteria. Critical complications have been reported with their long-term use such as thrombosis, anemia and gastrointestinal bleeding. Damage to red blood cells (RBCs), which occurs with nonphysiological blood flow conditions such as contact with foreign surfaces, high shear stress, and turbulence, is a major problem for the design and development of these systems. Even in the absence of hemolysis, cardiovascular devices (CAD) still cause cell injury and shortened RBC lifespans. This review summarizes various effects that occur to erythrocytes exposed to supraphysiological but sublethal stresses. Full article
(This article belongs to the Special Issue Advances in Biological Flows and Biomimetics, Volume II)
Show Figures

Figure 1

15 pages, 514 KiB  
Article
Drag Reduction of Turbulent Boundary Layers by Travelling and Non-Travelling Waves of Spanwise Wall Oscillations
by Martin Skote
Fluids 2022, 7(2), 65; https://doi.org/10.3390/fluids7020065 - 5 Feb 2022
Cited by 5 | Viewed by 2678
Abstract
Turbulence control in the form of a streamwise travelling wave of transverse wall motion was studied numerically by employing direct numerical simulations (DNS). Both total and phase averaging were utilised to examine the statistical behaviour of the turbulence affected by the wall forcing, [...] Read more.
Turbulence control in the form of a streamwise travelling wave of transverse wall motion was studied numerically by employing direct numerical simulations (DNS). Both total and phase averaging were utilised to examine the statistical behaviour of the turbulence affected by the wall forcing, with a focus on the skin friction. Comparison with results from pure temporal and spatial wall forcing are conducted, and a compilation of data is used to explore analogies with drag-reduced channel flow. Full article
(This article belongs to the Special Issue Turbulent Flow)
Show Figures

Figure 1

10 pages, 1352 KiB  
Article
The Development of Forecasting Technique for Cyclic Steam Stimulation Technology Effectiveness in Near-Wellbore Area
by Sergey Krivoshchekov, Alexander Kochnev and Kirill Vyatkin
Fluids 2022, 7(2), 64; https://doi.org/10.3390/fluids7020064 - 3 Feb 2022
Cited by 3 | Viewed by 2294
Abstract
The analytical review has shown that the scientific inquiry for effective technologies for high-viscosity oil field development is a critical task of the present-day oil industry. The paper presents a technique for determining the expediency and effectiveness of deploying the near-wellbore cyclic steam [...] Read more.
The analytical review has shown that the scientific inquiry for effective technologies for high-viscosity oil field development is a critical task of the present-day oil industry. The paper presents a technique for determining the expediency and effectiveness of deploying the near-wellbore cyclic steam stimulation technology for oil recovery enhancement. The method involves the calculation of process parameters of the technology cycle and the comparative analysis of cumulative oil production before the treatment (base case) and after its deployment. Separately, the work focuses on studying the impact of dynamic oil viscosity over the entire temperature range on the technology effectiveness and expediency. The laboratory studies showed dynamic viscosity correlation dependencies for six different oils of the Nozhovskaya group of oil fields (Russian Federation) characterized as viscous and highly viscous. As a case study of the proposed method application, a numerical simulation of the technology deployment was carried out for six oil samples. The calculations determined inexpediency of cyclic steam stimulation for one of the samples since oil well downtime for workover operation prevailed over the time of near-wellbore cooling. Full article
(This article belongs to the Special Issue Multiphase Flow in Pipes with and without Porous Media)
Show Figures

Graphical abstract

22 pages, 2190 KiB  
Article
Recent Upgrades in a 2D Turbulent Transport Solver Based on a Hybrid Discontinuous Galerkin Method for the Simulation of Fusion Plasma in Tokamak
by Giacomo Piraccini, Marcello Capasso, Manuel Scotto D’Abusco, Giorgio Giorgiani, Frédéric Schwander, Eric Serre, Hugo Bufferand, Guido Ciraolo and Patrick Tamain
Fluids 2022, 7(2), 63; https://doi.org/10.3390/fluids7020063 - 2 Feb 2022
Viewed by 2246
Abstract
The simulation of fusion plasmas in realistic magnetic configurations and tokamak geometries still requires the development of advanced numerical algorithms owing to the complexity of the problem. In this context, we propose a Hybrid Discontinuous Galerkin (HDG) method to solve 2D transport fluid [...] Read more.
The simulation of fusion plasmas in realistic magnetic configurations and tokamak geometries still requires the development of advanced numerical algorithms owing to the complexity of the problem. In this context, we propose a Hybrid Discontinuous Galerkin (HDG) method to solve 2D transport fluid equations in realistic magnetic and tokamak wall geometries. This high-order solver can handle magnetic equilibrium free structured and unstructured meshes allowing a much more accurate discretization of the plasma facing components than current solvers based on magnetic field aligned methods associated with finite-differences (volumes) discretization. In addition, the method allows for handling realistic magnetic equilibrium, eventually non steady, a critical point in the modeling of full discharges including ramp up and ramp down phases. In this paper, we introduce the HDG algorithm with a special focus on recent developments related to the treatment of the cross-field diffusive terms, and to an adaptive mesh refinement technique improving the numerical efficiency and robustness of the scheme. The updated solver is verified with a manufactured solution method, and numerical tests are provided to illustrate the new capabilities of the code. Full article
(This article belongs to the Special Issue Complex Fluids and Flows: Algorithms and Applications)
Show Figures

Figure 1

15 pages, 6126 KiB  
Review
Flow Control in Wings and Discovery of Novel Approaches via Deep Reinforcement Learning
by Ricardo Vinuesa, Oriol Lehmkuhl, Adrian Lozano-Durán and Jean Rabault
Fluids 2022, 7(2), 62; https://doi.org/10.3390/fluids7020062 - 1 Feb 2022
Cited by 39 | Viewed by 5790
Abstract
In this review, we summarize existing trends of flow control used to improve the aerodynamic efficiency of wings. We first discuss active methods to control turbulence, starting with flat-plate geometries and building towards the more complicated flow around wings. Then, we discuss active [...] Read more.
In this review, we summarize existing trends of flow control used to improve the aerodynamic efficiency of wings. We first discuss active methods to control turbulence, starting with flat-plate geometries and building towards the more complicated flow around wings. Then, we discuss active approaches to control separation, a crucial aspect towards achieving a high aerodynamic efficiency. Furthermore, we highlight methods relying on turbulence simulation, and discuss various levels of modeling. Finally, we thoroughly revise data-driven methods and their application to flow control, and focus on deep reinforcement learning (DRL). We conclude that this methodology has the potential to discover novel control strategies in complex turbulent flows of aerodynamic relevance. Full article
(This article belongs to the Special Issue External Aerodynamics)
Show Figures

Figure 1

27 pages, 6566 KiB  
Article
Convective Velocity Perturbations and Excess Gain in Flame Response as a Result of Flame-Flow Feedback
by Thomas Steinbacher and Wolfgang Polifke
Fluids 2022, 7(2), 61; https://doi.org/10.3390/fluids7020061 - 31 Jan 2022
Cited by 6 | Viewed by 2896
Abstract
Convective velocity perturbations (CVPs) are known to play an important role in the response of flames to acoustic perturbations and in thermoacoustic combustion instabilities. In order to elucidate the flow-physical origin of CVPs, the present study models the response of laminar premixed slit [...] Read more.
Convective velocity perturbations (CVPs) are known to play an important role in the response of flames to acoustic perturbations and in thermoacoustic combustion instabilities. In order to elucidate the flow-physical origin of CVPs, the present study models the response of laminar premixed slit flames to low amplitude perturbations of the upstream flow velocity with a reduced order flow decomposition approach: A linearized G-equation represents the shape and heat release rate of the perturbed flame, while the velocity perturbation field is decomposed into irrotational and solenoidal contributions. The former are determined with a conformal mapping from geometry and boundary conditions, whereas the latter are governed by flame front curvature and flow expansion across the flame, which generates baroclinic vorticity. High-resolution CFD analysis provides values of model parameters and confirms the plausibility of model results. This flow decomposition approach makes it possible to explicitly evaluate and analyze the respective contributions of irrotational and solenoidal flows to the flame response, and conversely the effect of flame perturbations on the flow. The use of the popular ad hoc hypothesis of convected velocity perturbation is avoided. It is found that convected velocity perturbations do not result from immediate acoustic-to-hydrodynamic mode conversion, but are generated by flame-flow feedback. In this sense, models for flame dynamics that rely on ad-hoc models for CVPs do not respect causality. Furthermore, analysis of the flame impulse response reveals that for the configuration investigated, flame-flow feedback is also responsible for “excess gain” of the flame response, that is, the magnitude of the flame frequency response above unity. Full article
(This article belongs to the Special Issue Stability and Dynamics of Gaseous Flames and Detonations)
Show Figures

Figure 1

7 pages, 241 KiB  
Communication
Can a Dynamo Mechanism Act at the Magnetopauses of Magnetic Rapidly Rotating Exoplanets?
by Elena Belenkaya
Fluids 2022, 7(2), 60; https://doi.org/10.3390/fluids7020060 - 31 Jan 2022
Viewed by 1996
Abstract
An astrophysical dynamo converts the kinetic energy of fluids into magnetic energy. Dynamo is a non-local process. Here, we consider whether a dynamo can operate at the magnetopauses of magnetic rapidly rotating planets. We analyze the main necessary condition for the work of [...] Read more.
An astrophysical dynamo converts the kinetic energy of fluids into magnetic energy. Dynamo is a non-local process. Here, we consider whether a dynamo can operate at the magnetopauses of magnetic rapidly rotating planets. We analyze the main necessary condition for the work of this type of dynamo—the rotation transfer from the planet to the magnetopause. We show the role of the current disc around a rapidly rotating magnetic planet in the redistribution of angular momentum depending on the direction of the external magnetic field, using the example of the Jupiter’s magnetodisc. Full article
(This article belongs to the Special Issue Galactic, Astrophysical and Planetary Dynamos)
18 pages, 130509 KiB  
Article
Effect of Initial Water Saturation on Oil Displacement Efficiency by Nanosuspensions
by Dmitriy Guzei, Vladimir Zhigarev, Valery Rudyak, Sofia Ivanova and Andrey Minakov
Fluids 2022, 7(2), 59; https://doi.org/10.3390/fluids7020059 - 31 Jan 2022
Cited by 5 | Viewed by 2628
Abstract
This article deals with the study of the initial water saturation effect of a porous medium on the oil recovery factor using a water-based nanosuspension. The initial water saturation of the porous medium in the computations varied within the range from 0 to [...] Read more.
This article deals with the study of the initial water saturation effect of a porous medium on the oil recovery factor using a water-based nanosuspension. The initial water saturation of the porous medium in the computations varied within the range from 0 to 90%. The nanoparticle SiO2 concentration varied from 0 to 1 wt%. The particle sizes were equal to 5, 18, 22, and 50 nm. Experimentally measured wetting angles and the interfacial tension coefficient depending on the concentration and size of nanoparticles were used in computations. A mathematical model was developed, describing the transfer and diffusion of nanoparticles within the aqueous phase during immiscible displacement of oil by nanosuspension from a porous medium. Using the developed model, a systematic computational study of the effect of the initial water saturation of the core micromodel on the oil recovery factor using nanosuspension was carried out. It was revealed that with an increase in the initial water saturation, the oil recovery factor monotonically decreased in the case of displacement both by water and nanosuspension. It was shown that with an increase in the concentration of nanoparticles and a decrease in their size, the oil recovery factor increased. At that, the relative increase in the recovery factor had a maximum at an initial water saturation equal to 60%. Full article
(This article belongs to the Collection Advances in Flow of Multiphase Fluids and Granular Materials)
Show Figures

Figure 1

Previous Issue
Next Issue
Back to TopTop