A level set-based solver for two-phase incompressible flows: extension to magnetic fluids. (English) Zbl 07868312
Summary: Development of a two-phase incompressible solver for magnetic flows in the magnetostatic case is presented. The proposed numerical toolkit couples the Navier-Stokes equations of hydrodynamics with Maxwell’s equations of electromagnetism to model the behaviour of magnetic flows in the presence of a magnetic field. To this end, a rigorous implementation of a second-order two-phase solver for incompressible nonmagnetic flows is introduced first. This solver is implemented in the finite-difference framework, where a fifth-order conservative level set method is employed to capture the evolution of the interface, along with an incompressible solver based on the projection scheme to model the fluids. The solver demonstrates excellent performance even with high density ratios across the interface (Atwood number \(\approx 1\)), while effectively preserving the mass conservation property. Subsequently, the numerical discretisation of Maxwell’s equations under the magnetostatic assumption is described in detail, utilising the vector potential formulation. The primary second-order solver for two-phase flows is extended to the case of magnetic flows, by incorporating the Lorentz force into the momentum equation, accounting for high magnetic permeability ratios across the interface. The implemented solver is then utilised for examining the deformation of ferrofluid droplets in both quiescent and shear flow regimes across various susceptibility values of the droplets. The results suggest that increasing the susceptibility value of the ferrofluid droplet can affect its deformation and rotation in low capillary regimes. In higher capillary flows, increasing the magnetic permeability jump across the interface can further lead to droplet breakup as well. The effect of this property is also investigated for the Rayleigh-Taylor instability growth in magnetic fluids.
MSC:
| 76M99 | Basic methods in fluid mechanics |
| 76M20 | Finite difference methods applied to problems in fluid mechanics |
| 76W05 | Magnetohydrodynamics and electrohydrodynamics |
| 76T10 | Liquid-gas two-phase flows, bubbly flows |
| 76E25 | Stability and instability of magnetohydrodynamic and electrohydrodynamic flows |
Keywords:
finite difference-based fifth-order conservative level set method; Navier-Stokes equations; Maxwell equations; ferrofluid droplet; Rayleigh-Taylor instabilityReferences:
| [1] | Afkhami, S., Renardy, Y., Renardy, M., Riffle, J. S., and St Pierre, T.. 2008. “Field-induced Motion of Ferrofluid Droplets Through Immiscible Viscous Media.” Journal of Fluid Mechanics610:363-380. . · Zbl 1147.76063 |
| [2] | Afkhami, S., Tyler, A., Renardy, Y., Renardy, M., Pierre, T. S., Woodward, R., and Riffle, J. S.. 2010. “Deformation of a Hydrophobic Ferrofluid Droplet Suspended in a Viscous Medium Under Uniform Magnetic Fields.” Journal of Fluid Mechanics663:358-384. . · Zbl 1205.76291 |
| [3] | Awasthi, M. K.2014. “Viscous Potential Flow Analysis of Magnetohydrodynamic Rayleigh-Taylor Instability with Heat and Mass Transfer.” International Journal of Dynamics and Control2 (3): 254-261. . |
| [4] | Bacri, J., and Salin, D.. 1982. “Instability of Ferrofluid Magnetic Drops Under Magnetic Field.” Journal de Physique Lettres43 (17): 649-654. . |
| [5] | Bacri, J.-C., Salin, D., and Massart, R.. 1982. “Study of the Deformation of Ferrofluid Droplets in a Magnetic Field.” Journal de Physique Lettres43 (6): 179-184. . |
| [6] | Balay, S., Gropp, W. D., McInnes, L. C., and Smith, B. F.. 1997. “Efficient Management of Parallelism in Object Oriented Numerical Software Libraries.” In: Modern Software Tools in Scientific Computing, edited by E. Arge, A. M. Bruaset and H. P. Langtangen, 163-202. Boston, MA: Birkhäuser Press. . · Zbl 0882.65154 |
| [7] | Bijarchi, M. A., Favakeh, A., Alborzi, S., and Shafii, M. B.. 2021. “Experimental Investigation of on-demand Ferrofluid Droplet Generation in Microfluidics Using a Pulse-width Modulation Magnetic Field with Proposed Correlation.” Sensors and Actuators B: Chemical329:129274. . |
| [8] | Bijarchi, M. A., and Shafii, M. B.. 2020. “Experimental Investigation on the Dynamics of on-demand Ferrofluid Drop Formation Under a Pulse-width-modulated Nonuniform Magnetic Field.” Langmuir36 (26): 7724-7740. . |
| [9] | Boniou, V., Schmitt, T., and Vié, A.. 2022. “Comparison of Interface Capturing Methods for the Simulation of Two-phase Flow in a Unified Low-Mach Framework.” International Journal of Multiphase Flow149:103957. . · Zbl 07599632 |
| [10] | Brackbill, J. U., Kothe, D. B., and Zemach, C.. 1992. “A Continuum Method for Modeling Surface Tension.” Journal of Computational Physics100 (2): 335-354. . · Zbl 0775.76110 |
| [11] | Bussmann, M., Kothe, D. B., and Sicilian, J. M.. 2002. “Modeling High Density Ratio Incompressible Interfacial Flows.” In Proceedings of the ASME 2002 Joint U.S.-European Fluids Engineering Division Conference. Volume 1: Fora, Parts A and B, 707-713. Montreal, QC: American Society of Mechanical Engineers (ASME). |
| [12] | Chen, D., Tong, X., Xie, B., Xiao, F., and Li, Y.. 2023. “An Accurate and Efficient Multiphase Solver Based on THINC Scheme and Adaptive Mesh Refinement.” International Journal of Multiphase Flow162:104409. . |
| [13] | Chorin, A. J.1997. “A Numerical Method for Solving Incompressible Viscous Flow Problems.” Journal of Computational Physics135 (2): 118-125. . · Zbl 0899.76283 |
| [14] | Clift, R., Grace, J. R., and Weber, M. E.. 2005. Bubbles, Drops, and Particles. Mineola, NY: Dover Publications, Inc. |
| [15] | Cui, X., Habashi, W. G., and Casseau, V.. 2021. “Multiphase SPH Modelling of Supercooled Large Droplets Freezing on Aircraft Surfaces.” International Journal of Computational Fluid Dynamics35 (1-2): 79-92. . · Zbl 1498.76100 |
| [16] | Cunha, L. H., Siqueira, I. R., Oliveira, T. F., and Ceniceros, H. D.. 2018. “Field-induced Control of Ferrofluid Emulsion Rheology and Droplet Break-up in Shear Flows.” Physics of Fluids30 (12): 122110. . |
| [17] | Davidson, P. A.2001. An Introduction to Magnetohydrodynamics. Cambridge, UK: Cambridge University Press. . · Zbl 0974.76002 |
| [18] | Desjardins, O., Blanquart, G., Balarac, G., and Pitsch, H.. 2008. “High Order Conservative Finite Difference Scheme for Variable Density Low Mach Number Turbulent Flows.” Journal of Computational Physics227 (15): 7125-7159. . · Zbl 1201.76139 |
| [19] | Desjardins, O., Moureau, V., and Pitsch, H.. 2008. “An Accurate Conservative Level Set/Ghost Fluid Method for Simulating Turbulent Atomization.” Journal of Computational Physics227 (18): 8395-8416. . · Zbl 1256.76051 |
| [20] | Desjardins, O., and Pitsch, H.. 2009. “A Spectrally Refined Interface Approach for Simulating Multiphase Flows.” Journal of Computational Physics228 (5): 1658-1677. . · Zbl 1409.76096 |
| [21] | Ding, H., Spelt, P. D., and Shu, C.. 2007. “Diffuse Interface Model for Incompressible Two-phase Flows with Large Density Ratios.” Journal of Computational Physics226 (2): 2078-2095. . · Zbl 1388.76403 |
| [22] | El-Dib, Y. O.1994. “Nonlinear Hydrodynamic Rayleigh-Taylor Instability of Viscous Magnetic Fluids: Effect of a Tangential Magnetic Field.” Journal of Plasma Physics51 (1): 1-11. . |
| [23] | Enright, D., Fedkiw, R., Ferziger, J., and Mitchell, I.. 2002. “A Hybrid Particle Level Set Method for Improved Interface Capturing.” Journal of Computational Physics183 (1): 83-116. . · Zbl 1021.76044 |
| [24] | Fedkiw, R. P., Aslam, T., Merriman, B., and Osher, S.. 1999. “A Non-oscillatory Eulerian Approach to Interfaces in Multimaterial Flows (the Ghost Fluid Method).” Journal of Computational Physics152 (2): 457-492. . · Zbl 0957.76052 |
| [25] | Finster, M., Chaves, H., and Schwarze, R.. 2023. “Numerical Simulation of the Accumulation Process in Modulated Liquid Jets of Defined Volume Using the Volume-of-fluid Methodology.” International Journal of Multiphase Flow160:104356. . |
| [26] | Fonty, T., Ferrand, M., Leroy, A., Joly, A., and Violeau, D.. 2019. “Mixture Model for Two-phase Flows with High Density Ratios: A Conservative and Realizable SPH Formulation.” International Journal of Multiphase Flow111:158-174. . |
| [27] | Gao, D., Morley, N. B., and Dhir, V.. 2004. “Understanding Magnetic Field Gradient Effect From a Liquid Metal Droplet Movement.” Journal of Fluids Engineering126 (1): 120-124. . |
| [28] | Gibou, F., Fedkiw, R., and Osher, S.. 2018. “A Review of Level-set Methods and Some Recent Applications.” Journal of Computational Physics353:82-109. . · Zbl 1380.65196 |
| [29] | Gottlieb, S., and Shu, C.-W.. 1998. “Total Variation Diminishing Runge-Kutta Schemes.” Mathematics of Computation67 (221): 73-85. . · Zbl 0897.65058 |
| [30] | Guermond, J.-L., and Quartapelle, L.. 2000. “A Projection FEM for Variable Density Incompressible Flows.” Journal of Computational Physics165 (1): 167-188. . · Zbl 0994.76051 |
| [31] | Guido, S., and Villone, M.. 1998. “Three-dimensional Shape of a Drop Under Simple Shear Flow.” Journal of Rheology42 (2): 395-415. . |
| [32] | Harris, E.1962. “Rayleigh-Taylor Instabilities of a Collapsing Cylindrical Shell in a Magnetic Field.” The Physics of Fluids5 (9): 1057-1062. . · Zbl 0116.19904 |
| [33] | Hassan, M. R., and Wang, C.. 2019. “Magnetic Field Induced Ferrofluid Droplet Breakup in a Simple Shear Flow At a Low Reynolds Number.” Physics of Fluids31 (12). . |
| [34] | Hassan, M. R., Zhang, J., and Wang, C.. 2018. “Deformation of a Ferrofluid Droplet in Simple Shear Flows Under Uniform Magnetic Fields.” Physics of Fluids30 (9). . |
| [35] | Herrmann, M.2008. “A Balanced Force Refined Level Set Grid Method for Two-phase Flows on Unstructured Flow Solver Grids.” Journal of Computational Physics227 (4): 2674-2706. . · Zbl 1388.76252 |
| [36] | Hu, Y., Li, D., and Niu, X.. 2018. “Phase-field-based Lattice Boltzmann Model for Multiphase Ferrofluid Flows.” Physical Review E98 (3): 033301. . |
| [37] | Huang, Z., Lin, G., and Ardekani, A. M.. 2019. “A Mixed Upwind/central WENO Scheme for Incompressible Two-phase Flows.” Journal of Computational Physics387:455-480. . · Zbl 1452.76154 |
| [38] | Huang, Z., Lin, G., and Ardekani, A. M.. 2020. “Consistent, Essentially Conservative and Balanced-force Phase-field Method to Model Incompressible Two-phase Flows.” Journal of Computational Physics406:109192. . · Zbl 1453.76131 |
| [39] | Huang, H., Ying, A., and Abdou, M.. 2002. “3D MHD Free Surface Fluid Flow Simulation Based on Magnetic-field Induction Equations.” Fusion Engineering and Design63-64:361-368. . |
| [40] | Jesus, W. C., Roma, A. M., and Ceniceros, H. D.. 2018. “Deformation of a Sheared Magnetic Droplet in a Viscous Fluid.” Communications in Computational Physics24 (2): 332-355. . · Zbl 1488.76146 |
| [41] | Jiang, G.-S., and Shu, C.-W.. 1996. “Efficient Implementation of Weighted ENO Schemes.” Journal of Computational Physics126 (1): 202-228. . · Zbl 0877.65065 |
| [42] | Jiang, G.-S., and Wu, C.-C.. 1999. “A High-order WENO Finite Difference Scheme for the Equations of Ideal Magnetohydrodynamics.” Journal of Computational Physics150 (2): 561-594. . · Zbl 0937.76051 |
| [43] | Kennedy, M. R., Pozrikidis, C., and Skalak, R.. 1994. “Motion and Deformation of Liquid Drops, and the Rheology of Dilute Emulsions in Simple Shear Flow.” Computers & Fluids23 (2): 251-278. . · Zbl 0810.76003 |
| [44] | Kim, J., and Moin, P.. 1985. “Application of a Fractional-step Method to Incompressible Navier-Stokes Equations.” Journal of Computational Physics59 (2): 308-323. . · Zbl 0582.76038 |
| [45] | Lai, B., Ke, Z., Wang, Z., Zhu, R., Gao, R., Mao, Y., and Zhang, Y.. 2023. “A Lattice Boltzmann Front-tracking Interface Capturing Method Based on Neural Network for Gas-liquid Two-phase Flow.” International Journal of Computational Fluid Dynamics37 (1): 49-66. . · Zbl 07769888 |
| [46] | Li, X., Dong, Z.-Q., Li, Y., Wang, L.-P., Niu, X.-D., Yamaguchi, H., Li, D.-C., and Yu, P.. 2022. “A Fractional-Step Lattice Boltzmann Method for Multiphase Flows with Complex Interfacial Behavior and Large Density Contrast.” International Journal of Multiphase Flow149:103982. . |
| [47] | Li, X., Yu, P., Niu, X.-D., Li, D.-C., and Yamaguchi, H.. 2021. “A Magnetic Field Coupling Lattice Boltzmann Model and Its Application on the Merging Process of Multiple-ferrofluid-droplet System.” Applied Mathematics and Computation393:125769. . · Zbl 1462.76209 |
| [48] | Liu, X.-D., Osher, S., and Chan, T.. 1994. “Weighted Essentially Non-oscillatory Schemes.” Journal of Computational Physics115 (1): 200-212. . · Zbl 0811.65076 |
| [49] | Majidi, M., Bijarchi, M. A., Arani, A. G., Rahimian, M. H., and Shafii, M. B.. 2022. “Magnetic Field-induced Control of a Compound Ferrofluid Droplet Deformation and Breakup in Shear Flow Using a Hybrid Lattice Boltzmann-finite Difference Method.” International Journal of Multiphase Flow146:103846. . |
| [50] | Makaremi-Esfarjani, P., and Najafi-Yazdi, A.. 2022. “Characteristic Boundary Conditions for Magnetohydrodynamic Equations.” Computers & Fluids241:105461. . · Zbl 1521.76566 |
| [51] | Mefford, O. T., Woodward, R. C., Goff, J. D., Vadala, T., Pierre, T. G. S., Dailey, J. P., and Riffle, J. S.. 2007. “Field-induced Motion of Ferrofluids Through Immiscible Viscous Media: Testbed for Restorative Treatment of Retinal Detachment.” Journal of Magnetism and Magnetic Materials311 (1): 347-353. . |
| [52] | Meng, W., Liao, L., Chen, M., Yu, C.-h., Li, J., and An, R.. 2022. “An Enhanced CLSVOF Method with an Algebraic Second-reconstruction Step for Simulating Incompressible Two-phase Flows.” International Journal of Multiphase Flow154:104151. . |
| [53] | Mirjalili, S., Jain, S. S., and Dodd, M.. 2017. “Interface-capturing Methods for Two-phase Flows: An Overview and Recent Developments.” Center for Turbulence Research Annual Research Briefs2017 (117-135): 13. |
| [54] | Mizuno, K., Asahara, M., Kamiya, T., and Miyasaka, T.. 2022. “Finite Difference and Reinitialization Methods with Level Set to Interfacial Area Transport Equations for Gas-liquid Two-phase Flows.” International Journal of Computational Fluid Dynamics36 (5): 361-383. . · Zbl 1506.76120 |
| [55] | Morinishi, Y., Lund, T. S., Vasilyev, O. V., and Moin, P.. 1998. “Fully Conservative Higher Order Finite Difference Schemes for Incompressible Flow.” Journal of Computational Physics143 (1): 90-124. . · Zbl 0932.76054 |
| [56] | Nourgaliev, R., and Theofanous, T.. 2007. “High-fidelity Interface Tracking in Compressible Flows: Unlimited Anchored Adaptive Level Set.” Journal of Computational Physics224 (2): 836-866. . · Zbl 1124.76043 |
| [57] | Olsson, E., and Kreiss, G.. 2005. “A Conservative Level Set Method for Two Phase Flow.” Journal of Computational Physics210 (1): 225-246. . · Zbl 1154.76368 |
| [58] | Olsson, E., Kreiss, G., and Zahedi, S.. 2007. “A Conservative Level Set Method for Two Phase Flow II.” Journal of Computational Physics225 (1): 785-807. . · Zbl 1256.76052 |
| [59] | Pierce, C. D.2001. “Progress-Variable Approach for Large-Eddy Simulation of Turbulent Combustion.” Ph.D. thesis, Stanford University. |
| [60] | Popinet, S.2018. “Numerical Models of Surface Tension.” Annual Review of Fluid Mechanics50 (1): 49-75. . · Zbl 1384.76016 |
| [61] | Prosperetti, A.1981. “Motion of Two Superposed Viscous Fluids.” Physics of Fluids24 (7): 1217-1223. . · Zbl 0469.76035 |
| [62] | Prosperetti, A., and Tryggvason, G.. 2009. Computational Methods for Multiphase Flow. Cambridge, UK: Cambridge University Press. . · Zbl 1166.76004 |
| [63] | Radman, S., Fiorina, C., and Pautz, A.. 2021. “Development of a Novel Two-Phase Flow Solver for Nuclear Reactor Analysis: Algorithms, Verification and Implementation in OpenFOAM.” Nuclear Engineering and Design379:111178. . |
| [64] | Rosensweig, R. E.2013. Ferrohydrodynamics. Mineola, NY: Dover Publications, Inc. |
| [65] | Scardovelli, R., and Zaleski, S.. 1999. “Direct Numerical Simulation of Free-surface and Interfacial Flow.” Annual Review of Fluid Mechanics31 (1): 567-603. . |
| [66] | Sethian, J. A.1999. Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science. Vol. 3. Cambridge, UK: Cambridge University Press. · Zbl 0973.76003 |
| [67] | Su, Y., and Kim, S. H.. 2018. “An Improved Consistent, Conservative, Non-oscillatory and High Order Finite Difference Scheme for Variable Density Low Mach Number Turbulent Flow Simulation.” Journal of Computational Physics372:202-219. . · Zbl 1415.76487 |
| [68] | Sussman, M., and Puckett, E. G.. 2000. “A Coupled Level Set and Volume-of-fluid Method for Computing 3D and Axisymmetric Incompressible Two-phase Flows.” Journal of Computational Physics162 (2): 301-337. . · Zbl 0977.76071 |
| [69] | Sussman, M., Smereka, P., and Osher, S.. 1994. “A Level Set Approach for Computing Solutions to Incompressible Two-phase Flow.” Journal of Computational Physics114 (1): 146-159. . · Zbl 0808.76077 |
| [70] | Tagawa, T.2006. “Numerical Simulation of Two-phase Flows in the Presence of a Magnetic Field.” Mathematics and Computers in Simulation72 (2-6): 212-219. . · Zbl 1116.76460 |
| [71] | Taylor, G. I.1932. “The Viscosity of a Fluid Containing Small Drops of Another Fluid.” Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character138 (834): 41-48. . · Zbl 0005.32201 |
| [72] | Taylor, G. I.1934. “The Formation of Emulsions in Definable Fields of Flow.” Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character146 (858): 501-523. . |
| [73] | Tryggvason, G., Scardovelli, R., and Zaleski, S.. 2011. Direct Numerical Simulations of Gas-Liquid Multiphase Flows. Cambridge, UK: Cambridge University Press. . · Zbl 1226.76001 |
| [74] | Unverdi, S. O., and Tryggvason, G.. 1992. “A Front-tracking Method for Viscous, Incompressible, Multi-fluid Flows.” Journal of Computational Physics100 (1): 25-37. . · Zbl 0758.76047 |
| [75] | Velikovich, A. L., Cochran, F., and Davis, J.. 1996. “Suppression of Rayleigh-Taylor Instability in Z-pinch Loads with Tailored Density Profiles.” Physical Review Letters77 (5): 853-856. . |
| [76] | Velikovich, A. L., and Schmit, P.. 2015. “Bell-Plesset Effects in Rayleigh-Taylor Instability of Finite-thickness Spherical and Cylindrical Shells.” Physics of Plasmas22 (12): 122711. . |
| [77] | Voltairas, P., Fotiadis, D., and Massalas, C.. 2001. “Elastic Stability of Silicone Ferrofluid Internal Tamponade (SFIT) in Retinal Detachment Surgery.” Journal of Magnetism and Magnetic Materials225 (1-2): 248-255. . |
| [78] | Wu, C.-C.2007. “A High Order WENO Finite Difference Scheme for Incompressible Fluids and Magnetohydrodynamics.” Geophysical and Astrophysical Fluid Dynamics101 (1): 37-61. . · Zbl 1505.76068 |
| [79] | Yang, J., Mao, S., He, X., Yang, X., and He, Y.. 2019. “A Diffuse Interface Model and Semi-implicit Energy Stable Finite Element Method for Two-phase Magnetohydrodynamic Flows.” Computer Methods in Applied Mechanics and Engineering356:435-464. . · Zbl 1441.76143 |
| [80] | Zarei Saleh Abad, M., Ebrahimi-Dehshali, M., Bijarchi, M. A., Shafii, M. B., and Moosavi, A.. 2019. “Visualization of Pool Boiling Heat Transfer of Magnetic Nanofluid.” Heat Transfer—Asian Research48 (7): 2700-2713. . |
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.
