\r\n\tA number of advanced combustion technologies have been introduced to improve performance, fuel economy and emissions levels. Research in combustion technology has highlighted the importance of new fuels in reducing the petroleum dependence and achieving high efficiency with low pollutant formation. \r\n\tThe purpose of this book is to collect interesting and original studies on combustion methods, advanced combustion strategies and new fuels able to achieve efficiency improvements and environment compliance. \r\n\tContributions in which experimental, theoretical and computation approaches are applied to explore how fuel properties and composition affect advanced combustion systems and how advanced combustion technology can maximize engine efficiency and be environment-friendly are invited and appreciated.
",isbn:null,printIsbn:"979-953-307-X-X",pdfIsbn:null,doi:null,price:0,priceEur:0,priceUsd:0,slug:null,numberOfPages:0,isOpenForSubmission:!1,hash:"025ed42cace07d115f7220a5a7b26f56",bookSignature:"Prof. Ornella Chiavola",publishedDate:null,coverURL:"https://cdn.intechopen.com/books/images_new/8231.jpg",keywords:"Combustion process, Performance, Fuel consumption, Emissions levels, Efficiency improvements, Environmental impact, Sustainable fuels, Fossil fuels, New combustion methods, Combustion analysis, Pollutant formation, Efficiency",numberOfDownloads:null,numberOfWosCitations:0,numberOfCrossrefCitations:0,numberOfDimensionsCitations:0,numberOfTotalCitations:0,isAvailableForWebshopOrdering:!0,dateEndFirstStepPublish:"May 14th 2019",dateEndSecondStepPublish:"September 17th 2019",dateEndThirdStepPublish:"November 16th 2019",dateEndFourthStepPublish:"February 4th 2020",dateEndFifthStepPublish:"April 4th 2020",remainingDaysToSecondStep:"a year",secondStepPassed:!0,currentStepOfPublishingProcess:5,editedByType:null,kuFlag:!1,biosketch:null,coeditorOneBiosketch:null,coeditorTwoBiosketch:null,coeditorThreeBiosketch:null,coeditorFourBiosketch:null,coeditorFiveBiosketch:null,editors:[{id:"201229",title:"Prof.",name:"Ornella",middleName:null,surname:"Chiavola",slug:"ornella-chiavola",fullName:"Ornella Chiavola",profilePictureURL:"https://mts.intechopen.com/storage/users/201229/images/system/201229.jfif",biography:"Ornella Chiavola completed a PhD in Mechanical Engineering at Engineering Department of the University of ‘La Sapienza’ – Rome, Italy in 1999. \r\n\r\nShe is currently working as an Associate Professor at ‘ROMA TRE’ University, Engineering Department, where she teaches ‘Interaction between machines and environment’ at Master Degree of Mechanical Engineering and at Master Degree of Aeronautical Engineering. \r\n\r\nFrom 2012, she is a member of the Teaching Board of Mechanical and Industrial Engineering PhD at ‘ROMA TRE’ University. \r\nShe has published numerous articles in international journals and is a member of the editorial board for several journals concerning her field of research.",institutionString:"Roma Tre University",position:null,outsideEditionCount:0,totalCites:0,totalAuthoredChapters:"2",totalChapterViews:"0",totalEditedBooks:"0",institution:{name:"Roma Tre University",institutionURL:null,country:{name:"Italy"}}}],coeditorOne:null,coeditorTwo:null,coeditorThree:null,coeditorFour:null,coeditorFive:null,topics:[{id:"11",title:"Engineering",slug:"engineering"}],chapters:null,productType:{id:"1",title:"Edited Volume",chapterContentType:"chapter",authoredCaption:"Edited by"},personalPublishingAssistant:{id:"301331",firstName:"Mia",lastName:"Vulovic",middleName:null,title:"Mrs.",imageUrl:"https://mts.intechopen.com/storage/users/301331/images/8498_n.jpg",email:"mia.v@intechopen.com",biography:"As an Author Service Manager, my responsibilities include monitoring and facilitating all publishing activities for authors and editors. From chapter submission and review to approval and revision, copyediting and design, until final publication, I work closely with authors and editors to ensure a simple and easy publishing process. I maintain constant and effective communication with authors, editors and reviewers, which allows for a level of personal support that enables contributors to fully commit and concentrate on the chapters they are writing, editing, or reviewing. I assist authors in the preparation of their full chapter submissions and track important deadlines and ensure they are met. I help to coordinate internal processes such as linguistic review, and monitor the technical aspects of the process. As an ASM I am also involved in the acquisition of editors. Whether that be identifying an exceptional author and proposing an editorship collaboration, or contacting researchers who would like the opportunity to work with IntechOpen, I establish and help manage author and editor acquisition and contact."}},relatedBooks:[{type:"book",id:"1591",title:"Infrared Spectroscopy",subtitle:"Materials Science, Engineering and Technology",isOpenForSubmission:!1,hash:"99b4b7b71a8caeb693ed762b40b017f4",slug:"infrared-spectroscopy-materials-science-engineering-and-technology",bookSignature:"Theophile Theophanides",coverURL:"https://cdn.intechopen.com/books/images_new/1591.jpg",editedByType:"Edited by",editors:[{id:"37194",title:"Dr.",name:"Theophanides",surname:"Theophile",slug:"theophanides-theophile",fullName:"Theophanides Theophile"}],productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"3092",title:"Anopheles mosquitoes",subtitle:"New insights into malaria vectors",isOpenForSubmission:!1,hash:"c9e622485316d5e296288bf24d2b0d64",slug:"anopheles-mosquitoes-new-insights-into-malaria-vectors",bookSignature:"Sylvie Manguin",coverURL:"https://cdn.intechopen.com/books/images_new/3092.jpg",editedByType:"Edited by",editors:[{id:"50017",title:"Prof.",name:"Sylvie",surname:"Manguin",slug:"sylvie-manguin",fullName:"Sylvie Manguin"}],productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"3161",title:"Frontiers in Guided Wave Optics and Optoelectronics",subtitle:null,isOpenForSubmission:!1,hash:"deb44e9c99f82bbce1083abea743146c",slug:"frontiers-in-guided-wave-optics-and-optoelectronics",bookSignature:"Bishnu Pal",coverURL:"https://cdn.intechopen.com/books/images_new/3161.jpg",editedByType:"Edited by",editors:[{id:"4782",title:"Prof.",name:"Bishnu",surname:"Pal",slug:"bishnu-pal",fullName:"Bishnu Pal"}],productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"72",title:"Ionic Liquids",subtitle:"Theory, Properties, New Approaches",isOpenForSubmission:!1,hash:"d94ffa3cfa10505e3b1d676d46fcd3f5",slug:"ionic-liquids-theory-properties-new-approaches",bookSignature:"Alexander Kokorin",coverURL:"https://cdn.intechopen.com/books/images_new/72.jpg",editedByType:"Edited by",editors:[{id:"19816",title:"Prof.",name:"Alexander",surname:"Kokorin",slug:"alexander-kokorin",fullName:"Alexander Kokorin"}],productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"1373",title:"Ionic Liquids",subtitle:"Applications and Perspectives",isOpenForSubmission:!1,hash:"5e9ae5ae9167cde4b344e499a792c41c",slug:"ionic-liquids-applications-and-perspectives",bookSignature:"Alexander Kokorin",coverURL:"https://cdn.intechopen.com/books/images_new/1373.jpg",editedByType:"Edited by",editors:[{id:"19816",title:"Prof.",name:"Alexander",surname:"Kokorin",slug:"alexander-kokorin",fullName:"Alexander Kokorin"}],productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"57",title:"Physics and Applications of Graphene",subtitle:"Experiments",isOpenForSubmission:!1,hash:"0e6622a71cf4f02f45bfdd5691e1189a",slug:"physics-and-applications-of-graphene-experiments",bookSignature:"Sergey Mikhailov",coverURL:"https://cdn.intechopen.com/books/images_new/57.jpg",editedByType:"Edited by",editors:[{id:"16042",title:"Dr.",name:"Sergey",surname:"Mikhailov",slug:"sergey-mikhailov",fullName:"Sergey Mikhailov"}],productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"371",title:"Abiotic Stress in Plants",subtitle:"Mechanisms and Adaptations",isOpenForSubmission:!1,hash:"588466f487e307619849d72389178a74",slug:"abiotic-stress-in-plants-mechanisms-and-adaptations",bookSignature:"Arun Shanker and B. Venkateswarlu",coverURL:"https://cdn.intechopen.com/books/images_new/371.jpg",editedByType:"Edited by",editors:[{id:"58592",title:"Dr.",name:"Arun",surname:"Shanker",slug:"arun-shanker",fullName:"Arun Shanker"}],productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"878",title:"Phytochemicals",subtitle:"A Global Perspective of Their Role in Nutrition and Health",isOpenForSubmission:!1,hash:"ec77671f63975ef2d16192897deb6835",slug:"phytochemicals-a-global-perspective-of-their-role-in-nutrition-and-health",bookSignature:"Venketeshwer Rao",coverURL:"https://cdn.intechopen.com/books/images_new/878.jpg",editedByType:"Edited by",editors:[{id:"82663",title:"Dr.",name:"Venketeshwer",surname:"Rao",slug:"venketeshwer-rao",fullName:"Venketeshwer Rao"}],productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"4816",title:"Face Recognition",subtitle:null,isOpenForSubmission:!1,hash:"146063b5359146b7718ea86bad47c8eb",slug:"face_recognition",bookSignature:"Kresimir Delac and Mislav Grgic",coverURL:"https://cdn.intechopen.com/books/images_new/4816.jpg",editedByType:"Edited by",editors:[{id:"528",title:"Dr.",name:"Kresimir",surname:"Delac",slug:"kresimir-delac",fullName:"Kresimir Delac"}],productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"3621",title:"Silver Nanoparticles",subtitle:null,isOpenForSubmission:!1,hash:null,slug:"silver-nanoparticles",bookSignature:"David Pozo Perez",coverURL:"https://cdn.intechopen.com/books/images_new/3621.jpg",editedByType:"Edited by",editors:[{id:"6667",title:"Dr.",name:"David",surname:"Pozo",slug:"david-pozo",fullName:"David Pozo"}],productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}}]},chapter:{item:{type:"chapter",id:"63437",title:"Interaction of Magnetic Fields on Ferrofluidic Taylor-Couette Flow",doi:"10.5772/intechopen.80301",slug:"interaction-of-magnetic-fields-on-ferrofluidic-taylor-couette-flow",body:'
1. Introduction
Ferrofluids [1] represent one example of magnetic particles with a variety of applications, ranging from liquid seals in rotating systems, their use in computer hard drives, as well as in laboratory experiments, to study geophysical flows [2, 3]. However, deeper fundamental study of their magnetohydrodynamics is inevitable, as there are many modeling assumptions that have been implemented to make theoretical descriptions tractable. Ferrofluids are manufactured fluids which consist of dispersions of magnetized nanoparticles in various kinds of liquid carriers. Stabilizing against agglomeration is typically achieved by the addition of a surfactant monolayer onto the magnetic particles. Being no magnetic field present, the magnetic nanoparticles are randomly orientated, which results in a zero-net magnetization of the fluid. Only effect of the present nanoparticles is a typically small alteration to the fluid’s density and viscosity. Applying a sufficiently strong magnetic field, the ferrofluid flows toward certain regions of the magnetic field. The hydrodynamics of the system can be significantly changed, and properties of the fluid, e.g., the viscosity, are altered. Therefore, nowadays, they are widely used by various car manufactures in order to provide adaptable suspensions, i.e., changing the damping characteristics of an antishock/damper due to modified viscosity of the ferrofluid under influence of an external applied magnetic field. Both numerical simulations and experiments show that any applied magnetic field, whether radial, axial, azimuthal, transversal, or any combination of the previous stabilizes the basic state [4, 5, 6, 7, 8, 9]. This stabilization effect depends, among others, on the particle-particle interactions of the fluid [9]. A strong interaction leads to large/crucial changes in flow behavior, which are driven by the formation of internal flow structures. While an axisymmetric field leaves the flow structures untouched (compared to classical Couette flow), a finite transverse, in-plane magnetic-field component breaks axisymmetry, enforcing new nonlinear effects, modify present flow structures, and can produce new flow structures.
Among various others, one important point to look at is the applied magnetic field, in particular its characteristics while penetrating a moving ferrofluid. In order to describe the hydrodynamics of ferrofluids in containers/boxes typically, models are used which assume that the internal magnetic field and the external applied magnetic field are equal [4, 5, 6]. However, such a simplification is only a leading order approximation, as the magnetic field in the container typically is modified, based on the magnetic susceptibility of the ferrofluid. First work trying to account for modifications [15] includes an expression for the angular velocity of the ferrofluid particles, which differs from the ferrofluid itself due to the magnetic torques acting on the particles. The results are significant differences from the external field, providing among others, a much changed body force in the governing equations. More concretely, the interaction between a uniform external field and the susceptibility of the ferrofluid-filled annulus can be shown to result in a field with an 1/r2 radial dependence superimposed on the external field.
Further common assumption when describing the hydrodynamics of magnetic particles as ferrofluids is to consider them as point-wise objects and ignore their finite size and specific properties. However, in real, such particles typically aggregate to form clusters having the form of chains, and thus they hinder the free flow of the fluid and increase the viscosity [9, 10, 11], and change their rheology and transport properties [14]. In this type of structure’s formation, it is also assumed that the interaction parameter is usually greater than unity [1], and thus the strength of the grain-grain interaction can be measured in terms of the total momentum of a particle. Regarding the magneto-dissipative structure of ferrofluid dynamics, one can derive, based on general principles, an equation for macroscopic ferrofluid dynamics. The magnetization’s relaxation equation includes an additional term that is proportional to the product of the magnetization’s magnitude and the symmetric part of the velocity gradient tensor, which can describe an elongational flow scaled by a so-called transport coefficientλ2. Such description also exists in the dynamics of nematic liquid crystals, as the flow alignment modifies the director field in an applied shear flow [16].
It is worth to mention that there are further crucial and important considerations for describing ferrofluids and in general fluids with magnetic particles which are also important as the discussed ones. However, these are not considered in the present chapter. Therefore, we refer to the huge amount of literature, which can be essentially found by the presented references. Two prominent examples of these further considerations are the modulation of the relaxation times of magnetic particles [7, 13], e.g., Neel’s or Braun’s and the consideration of monodisperse vs. polydisperse materials [7, 17, 18].
In this chapter, we will consider a rotating ferrofluid in Taylor-Couette geometry. First, we are going to describe the modifications in flow solutions and its structural properties and alterations in primary supercritical bifurcation thresholds induced/forced by an external applied magnetic field. Following, we will discuss two different approaches to come to a more realistic description of ferrofluidic flows. Therefore, we will modify and extend the basic ferrohydrodynamical equations. First, considering a modification of the internal magnetic field, which differs from the external applied one due to its internal magnetization. Second incorporating elongational flow. The latter aims to consider finite size effects of the magnetic particles, as well as agglomeration and chain formation procedures.
2. System and theoretical description
As a prototypical system to investigate the influences of a magnetic field on a rotating ferrofluid system, we consider Taylor-Couette flow (Taylor Couette system, TCS) [19, 20] driven by the independent differential rotation of two concentric cylinders (Figure 1).
Figure 1.
Schematic of the Taylor-Couette system with an external applied homogeneous magnetic field H→ext=Hxe→x+Hze→z.
The inner cylinder of radius r1 rotates at ω1, and the outer cylinder of radius r2 rotates at ω2. If not other stated, in this chapter, periodic boundary conditions are considered in the axial direction with periodicity of length λ and no-slip boundary conditions on the cylinder walls. The system is described using a cylindrical polar coordinate system rθz with a velocity field u→=uvw. The radius ratio of the cylinders is β=r1/r2, and the aspect ratio is Γ=λ/r2−r1. The gap between the cylinders is filled with a viscous, incompressible, and isothermal ferrofluid. An external homogeneous magnetic field H→ext=Hxe→x+Hze→z is applied, where x=rcosθ is the transverse direction and z is the axial direction, respectively (Figure 1).
The flow dynamics of an incompressible homogeneous monodisperse ferrofluid with kinematic viscosity ν and density ρ is governed by the incompressible Navier–Stokes equations, including magnetic terms and the continuity equation. Following, an approach based on the model of Niklas [4, 5] is used, where the magnetic fluid is assumed to be incompressible, nonconducting, and to have a constant temperature and a homogeneous distribution of magnetic particles. Using the gap d=r2−r1 as the length scale, the diffusion time τd=d2/ν as the time scale, scaling pressure with ρν2/d2, the magnetic field H→, and the magnetization M→ with ρ/μ00.5ν/d (μ0 is the magnetic constant, i.e., magnetic permeability of free space), the nondimensional governing equations are
∂tu→⋅∇u→−∇2u→∇p=M→⋅∇H→+∇×M→×H→/2,∇⋅u→=0E1
The cylinders are considered to be of no-slip type with velocity boundary conditions
u→r1,θz=0Re1,0andu→r2,θz=0Re20,E2
where the inner and outer Reynolds numbers are Re1=ωr1d/ν and Re2=ωr2d/ν, respectively.
2.1 Magnetization equations
In order to solve (1), one needs additional information about the magnetic particles themselves, in particular their properties and interaction with an applied magnetic field. Therefore, (1) is solved together with an equation that describes the magnetization of the ferrofluid. A first approximation is to use the equilibrium magnetization of an unperturbed state with a homogeneously magnetized ferrofluid at rest with the mean magnetic moments orientated in the direction of the magnetic field, M→eq=χH→, where χ is the magnetic susceptibility of the ferrofluid, determined using Langevin’s formula [21]. The ferrofluid considered for all simulations in this chapter is APG933 [22] with χ=0.9. However, a ferrofluid’s magnetization is also influenced by the flow field. For the presented numerical simulations, an approach based on the model of Niklas [4, 5, 6] is used. Furthermore, a stationary magnetization near equilibrium considering a small perturbation at the equilibrium magnetization M→−M→eq and small relaxation times Ωτ<<1, where Ω is the absolute value of half of the vorticity Ω and τ is the magnetic relaxation time. In the near-equilibrium approximation, Niklas [4] determined the relationship between the magnetization M→, the magnetic field H→, and the velocity u→ to be
M→−Meq=cNΩ→×H→E3
with Ω→=∇×u→/2 being the vorticity and the Niklas coefficient
cN2=τ/1/χ+τμ0H2/6μΦE4
where μ is the dynamic viscosity and Φ is the volume fraction of the magnetic material.
2.2 Ferrohydrodynamic equations of motions
The magnetization can be eliminated from (1) by using (3) in order to obtain the ferrohydrodynamic equation of motion [4]:
∂t+u→⋅∇u→−∇2u→+∇pM=−cN2H→∇⋅F→+H→×∇×F→/2E5
where F→=Ω→×H→ and pM are the dynamic pressure incorporating all magnetic terms which can be written as gradients. Note that the relaxation time τ could typically be defined by τ=τBD31/3, where τB is the Brownian relaxation time and D31/3 is the averaged diameter of the magnetic particles. For simplicity, to investigate the different effects, i.e., internal magnetization and elongational flow (see below), we assume that τ is constant, which is independent of the magnetic field H→. The considered relaxation time is based on the experimental result [22] τ=τAPG933/τd=0.0018.
Consider external field to be equal to internal one H→=H→ext, (5) can be simplified to
∂t+u→⋅∇u→−∇2u→+∇pM=−cN2H→∇⋅F→+H→×∇×F→/2E6
In this approach, the magnetic field and all the magnetic properties of the ferrofluid influence the velocity field only via the magnetic field, the Niklas parameters→N=sxe→x+sze→z [4, 5, 6], whereby magnetic field strengths are:
sx=22+χ/2+χ2−χ2η2HxcNandsz=HzcnE7
In this approach, sx and sz are dimensionless numbers, describing the strength of the applied magnetic field.
3. Ferrofluidic flow under finite magnetic field
3.1 Introduction and motivation
We start our discussion with the investigation of the influence of a homogeneous magnetic field on a ferrofluid. In this first part, we will assume that the internal magnetic field is equal to the external applied magnetic field, i.e., homogeneous. Moreover, we ignore in this section any kind of finite size effects, e.g., possible agglomeration and chain formation, later considered as elongational flow effect (Section 4). Therefore, we consider the ferrofluid particles as point-size objects without any 3D expansion. Structure, dynamics, symmetry properties, bifurcation, and stability behavior of different vortex structures are investigated for axial and transversal magnetic fields, as well as combinations of them, i.e., oblique magnetic field. We will see that a transversal magnetic field has significant influence. It modulates and changes the classical flow structures, as well the basic state (Circular Couette flow, CCF) as the classical primary supercritical ones, Taylor vortex flow (TVF) and the spiral vortex flow (spirals, SPI). They are replaced by wavy-modulated flow structures: wavy Taylor vortex flow (wTVF) and wavy spiral vortex flow (wSPI), as well as they forward bifurcate out of the modified 2-fold circular Couette flow (2-CCF).
3.2 Flow structures
3.2.1 Flow structures in absence of a magnetic field
There are two fundamental basic, topologically different flow structures appearing in (classical) Taylor-Couette flow (see also Nomenclature and abbreviations section). These are either toroidally closed Taylor vortices (Taylor vortex flow, TVF) and helical spiral vortices (spiral vortex flow/spirals, SPI) [20]. Figure 2 shows the structure of these fundamental solutions, TVF and SPI (here for a left-handed SPI with azimuthal wavenumber m=1), by using isosurfaces of the azimuthal vorticity. Because of the system symmetry z→−z, we only discuss left-handed SPI. Without other parameters, e.g., axial through flow, left- and right-handed SPI are mirror images of each other. These pure states can also be found modulated as wavy structures, which are topologically identical with the pure states. However, they show a wavy-like deformation that results from additional and other stimulated modes compared to the classical solutions.
Figure 2.
Azimuthal vorticity isosurfaces η over two axial wavelengths (for visual purpose) in absence of a magnetic field (N) of (a) TVF at Re2=0, Re1=100, and (b) SPI at Re2=−150, Re1=150. Red (yellow) isosurfaces correspond to positive (negative) values as indicated. The bottom row shows mode amplitudes um,n of the radial velocity field u corresponding to structures above over the m-n plane. The values are scaled regarding the maximum mode amplitude to be 1 (TVF: u0,1=1).
While classical wavy flow solutions [27, 28, 29, 30, 31, 32, 33, 34, 35, 36] are well known and have already been widely discussed in literature regarding TCS, the ones generated by a magnetic field are different. In particular, the magnetically induced wTVF are nonrotating structures having a pinned phase, contrary to the classic wTVF, which rotate azimuthally. Regarding the flow structure, the imposed magnetic field shrinks and expands the vortices at different cross sections of constant θ, resulting in an azimuthal alternating wide and narrow vortex belly structure.
3.2.2 Flow structures—mode coupling in magnetic fields
Following, we will have a quick look at the impact on the different applied magnetic field with respect to the flow structures and in particular their mode spectra. The modifications in NSE due to the additional terms incorporating magnetization and magnetic field result in various new phenomena. Most important is probably the fact that magnetic fields with a finite transversal component stimulate additional modes. Depending on different external applied magnetic fields, Figure 3 presents a schematic plot for the different stimulated modes, depending on the applied field.
Figure 3.
Schematic sketch of modes in the presence of different externally imposed (homogeneous) magnetic fields. Here, m represents the azimuthal, and n the axial mode index in the Fourier expansion (Appendix (17) and (18)). (a) Modes either in absence of any magnetic field or in a pure axial magnetic field. (b) Pure transversal magnetic field with additional stimulation of m=n±2 modes. (c) Superposition of axial and transversal fields with even further finite m=n±1 modes. Strength of mode amplitudes, the magnitudes are characterized from dark (large amplitudes) to bright (small amplitudes) colors (legend right). Red and blue lines indicate the dominant mode amplitudes of the flow structures.
Having a pure axial field sx=0sz>0 applied the mode spectra, the flow structures remain qualitatively similar to the no field (N) case (Figure 3(a)). A TVF still contains only the axisymmetric modes m=0 (with dominant n = 1, see also Figures 2 and 6). Similarly, L1-SPI, i.e., a left-handed SPI flow with azimuthal wavenumber m=1 contains modes only on the diagonal m=n (with dominant n=1, see also Figures 2 and 6). However, if a field with finite transversal component is applied additional mode become finite and the flow structures change. A pure transversal field sx>0sz=0 excites m=±2 modes lying on the secondary diagonal m=n±2 (Figure 3(b)). Thus, as a result, the pure states of TVF and SPI vanish in the presence of a transversal magnetic field. They become wavy-like modified, resulting in wTVF and wSPI. Finally, having a magnetic field with either axial and transversal component sx>0sz>0, there are further interactions, which result in stimulation of additional m=n±1 modes (Figure 3(c)). Moreover, these additional modes can induce further nonlinear mode couplings. As a result, the mode spectra contain more nonzero components (Figure 3(c) and Figure 6 for numerical results). Thus all flow structures in a magnetic field with transversal component differ qualitatively from the classical flows found in the absence of magnetic fields [6, 26].
3.3 Effect of finite magnetic fields
3.3.1 Basic state and symmetries
If only a pure axial magnetic field, Hextre→z (see Figure 1), the classical circular Couette flow (CCF) is observed, with u→CCF=0Ar+B/r0, where A=Re2−βRe1/1+β, and B=βRe1−βRe2/1+β1−β2, is a solution of (5) as the associated vorticity Ω→ccf=∇×u→ccf is parallel to the applied magnetic field and therefore all magnetic terms vanish. Although for magnetic fields that are orientated purely in the radial or azimuthal direction, the basic state changes but remains axisymmetric with deviations from CCF only having an azimuthal component [4, 23]. As a result, in all of these cases, the basic state is invariant to a number of symmetries (azimuthal rotation, mirror symmetry, and axial translation) whose actions on a general velocity field are
RΦuvwrθzt=uvwrθ+Φzt,Φ∈02π,E8
Kzuvwrθzt=uv−wrθ−zt,E9
Tαuvwrθzt=uvwrθz+αΓt,α∈NE10
On the contrary, for a field with a finite transverse component, the base state is no longer axisymmetric. Although the basic state remains invariant to Tα and a modified mirror symmetry KzH [24], a transverse magnetic field T breaks the continuous axisymmetry RΦ, resulting in a basic state with discrete symmetry Rπ, i.e., with azimuthal wavenumber m=2 [24]. Following, we will refer to this modified 2-fold symmetric CCF as 2-CCF. It is important to mention that this m=2 state is stationary [24], which is in contrast to the generic breaking of a SO2 symmetry which would result in a rotating wave.
3.3.2 Bifurcation behavior
As seen in Section 3.2, an applied magnetic field can stimulate additional modes and thus changes the Fourier mode spectrum of flow structures resulting in new solutions. The results are wavy-like modified flow states. It is to expect that these mode-stimulated modifications also influence and appear in the primary bifurcation thresholds of the solutions, as well as in their growth rate. Following, we will have a look at the bifurcation behavior of the different flow structures with the additionally field-induced modes explained in Section 3.3. In particular, we will focus on four different field configurations: no field (N)sx=0sz=0; pure axial field (A)sx=0sz=0.6; pure transversal field (T)sx=0.6sz=0; and oblique field (O)sx=0.6sz=0.6 ((7) for parameters sx,sz).
Considering these four, non, pure axial, pure transversal, and oblique magnetic fields, Figure 4 illustrates the stable forward bifurcating branches of toroidally closed TVF and of wTVF solutions, respectively. Figure 5 shows the bifurcation diagrams for helical SPI and wSPI in an analogous way. Note that the shown wavy-like modulated flow states only result from symmetry breaking due to finite transversal field component. In the Figures 4 and 5, the relative distance μ=Re1/Re1,c−1 is used as control parameter for easier comparison of the bifurcation branches (Re1,c stands for the critical onset Reynolds number Re1 of the flow structure, which depends on various system parameters. These are the aspect ratio Γ, the radius ratio r1/r2, the outer Reynolds number Re2, axial boundary conditions, and the field strength sx, sz, respectively). The main reason to use μ is the shift in the bifurcation onset, which varies depending on the applied field direction.
Figure 4.
Bifurcation diagrams of TVF and wTVF, respectively, in different magnetic fields: (a) no field (N), (b) axial (A), (c) transverse (T), and (d) oblique (O). Moduli um,n (18) of the radial flow field amplitudes at mid gap are shown versus the relative distance μ=Re1/Re1,c−1 from the onset of the respective vortex flow. Further parameters are: Re2=0,β=0.5,k=3.927 (cf. Figure 2).
Figure 5.
Bifurcation diagrams of SPI and wSPI, respectively, in different magnetic fields: (a) no field (N), (b) axial (A), (c) transverse (T), and (d) oblique (O). Moduli um,n (18) of the radial flow field amplitudes at mid gap are shown versus the relative distance μ=Re1/Re1,onset−1 from the onset of the respective vortex flow. Further parameters are: Re2=−150,β=0.5,k=3.927 (cf. Figure 2).
Displayed in the Figures 4 and 5 are the moduli um,n of the dominant, as well as the first higher harmonic modes of the corresponding structure (see Appendix (18) and (19)), e.g. mn=01 and mn=02 for TVF. Furthermore, the largest field-induced mode amplitudes u2,1 and u1,±1 (only in oblique configuration) of the wavy structures (cf. schematics in Figure 3) are included for comparison.
Axial field (A): As a pure axial magnetic field sx=0sz≠0 does not stimulate any further modes within the mode-spectra (Figure 3(a)), the flow states TVF and SPI and their structural properties remain qualitatively similar, and they are as without a field. However, quantitative, they eventually change. First important effect is the rise of the onsets to higher values of Re1 [6, 7]. Thereby, the upward shift depends on the various system control parameters (e.g., Re2,β,Γ). For parameters in Figure 4, the onset of TVF rises by about 25%, and for the parameters of Figure 5, the SPI bifurcation threshold moves upward in by about 15%. For the latter, consequently also corresponding SPI frequencies are increased in an axial field. Other change regarding Figure 4 is the amplitude of the modes with increasing the parameter μ (weakest for nonfield and strongest for oblique field configuration).
Transversal field (T): A transversal magnetic field sx≠0sz=0 stimulates additional azimuthal modes m=n±2 (Section 3.2) as illustrated in Figure 3(b). Thus, based on this knowledge, one can already expect a more dramatic influence on the classical TVF and SPI structures. As a pure axial field, a transversal field also delays the onset of the primary bifurcating structures, i.e., it stabilizes the basic state. Furthermore and more important is the fact that such a field also changes the flow structural properties in qualitative manner. The growth of the additional modulus u2,1 (18) (for toroidally closed flow states) with μ=Re1/Re1,onset−1 changes the pure TVF to a wavy-like modulated wTVF. Analog pure SPI changes to wSPI (shown for L1-(w)SPI in Figure 5) due to the additional m=n±2 (see Figure 3(b)), indicated by the growth of the modulus u3,1 in Figure 5. It is worth to mention that due to the mode stimulation and interaction, also the (1,-1) mode becomes excited, which is similar to structure R1-(w)SPI. Note, the same also holds for the symmetry-related mirror image L1-(w)SPI.
Thus, in the presence of a magnetic field with a finite transversal component, the pure/classical TVF and SPI structures are modified. They are replaced by their wavy-like counterparts, wTVF and wSPI, which now bifurcate as primary structured solutions forward out of the CCF-2 ground state; note that already the basic state has an additional mode-2 symmetry due to the symmetry-breaking effect of the transversal field and as such already differ from the classical CCF (Section 3.3.1). As for axial field, the magnitude of stabilization varies with given system control parameters. In general, it is slightly smaller compared to the stabilization effect of an axial field. For instance for parameters presented in Figures 4 and 5, the upward shift for the onset of wTVF (wSPI) is about 20% (12%) in an axial (transversal) magnetic field.
Important to mention that a finite transverse component sx≠0 breaks the system symmetry but does not change the degeneration between left- and right-handed SPIs. Both solutions keep bifurcating at the same critical values. The latter is different to the effect of an axial flow, for instance. Aside that such a flow also breaks the system symmetry, it separates both left- and right-handed SPI solutions. Under axial applied flow, they appear at different critical values.
Oblique field (O): As already seen for both pure field configurations, either sx≠0,sz=0 or sx=0,sz≠0, the onset of (w)TVF and (w)SPI is shifted to higher values. In fact in a field, obliquely orientated to the cylinder sx≠0sz≠0, the magnitude of shift is bigger compared to the pure cases. In addition, as discussed before, such a superposition of an axial and a transversal field stimulates also new, additional modes with higher mode index combinations m=n±1 (Figure 3(c) and for numerical simulations Figures 4(c,d), 5(c,d), and 6(e,f)).
Figure 6.
Azimuthal vorticity isosurfaces η over two axial wavelengths (for visual purpose) with pure transversal applied magnetic field (T): (a) wTVF at Re2=0, Re1=125, (b) wSPI at Re2=−150, Re1=160; and oblique applied magnetic field (O): (c) wTVF at Re2=0, Re1=150, and (d) wSPI at Re2=−150, Re1=190. Red (yellow) isosurfaces correspond to positive (negative) values as indicated. Bottom row shows mode amplitudes um,n of the radial velocity field u corresponding to structures above over the m-n-plane. Values are scaled with the maximum mode amplitude to be 1 (TVF: u0,1=1).
A common finding for (w)TVF and (w)SPI is that their dominant mode amplitudes for m=0 (for (w)TVF) and m=1 (for (w)SPI) grow with respect to the bifurcation onsets. Although the absolute value depends on the different applied magnetic field, comparing the moduli um,1 at the same relative distance μ from the respective onset (Figures 4 and 5, compare dominant mode amplitudes), one finds the following relations
for a given constant value c∈01 and m=0 in case of toroidally closed (w)TVF and m=1 for helical (w)SPI, respectively (Compare upper curves for dominant mode amplitudes um,1 in Figures 4 and 5).
For helical flow structures, the larger variety of mode interactions due to the basic spectra (e.g. m=n [m=−n] for pure L1-SPI [R-SPI]) results in significant variation of the higher order and secondary stimulated modes, depending on the field configuration (see in particular Figures 3(c), 4(d), and 5(d)).
Based on the both bifurcation diagrams presented in Figures 4 and 5, it is worth to remind that in the case of a finite transversal field component sx≠0, no pure solutions, TVF and SPI, exist. They become replaced by wavy-like modulated wTVF and wSPI. However, independent of the field, leading mode amplitude of the respective flow state, (w)TVF and (w)SPI, follows a typical square-root forward bifurcation, whereby the slope of its square grows with the applied field strength.
3.3.3 Vortex structures and mode contents
Having discussed the changes that different oriented magnetic fields can cause in the bifurcation behavior and mode stimulation of (w)TVF and (w)SPI, respectively, one can imagine that these also coincide with flow structural modifications. Now, we will address these structural changes of the vortices and the associated changes in the mode contents of the flows in more detail.
Figure 6 shows isosurfaces of the azimuthal vorticity η over two axial wavelengths (for visual purpose) for supercritical flow states (w)TVF and (w)SPI in the three different considered field configurations: (a,b) no field applied (N), (c,d) pure transverse field (T), and (e,f) oblique field (O) at parameters Re1 and Re2 as indicated.
A good quantity in order to visualize the field induced changes in the 3D vortex structures are isosurfaces of the azimuthal vorticity η. These surfaces appropriately convey structural details of the vortex flows in question [28].
Axial field (A): As a pure axial field does not break any continuous symmetry of the system, there is neither qualitative change in the structure of TVF and SPI in real space nor the mode structure of the flow in the m−n Fourier plane (not shown, but similar as Figure 2(a,b); only quantitative the mode amplitudes are slightly larger). TVF still consists of toroidally closed, rotationally symmetric vortex tubes with a mode spectrum containing only m=0 modes (Figure 2(a)). And also, the helically oriented, open vortex tubes of SPI (here for left-winding spiral L-SPI) containing m=n modes as shown in Figures 2(b) and 3 (mirror image of right-winding R-SPI analog) are unaffected by an axial field.
Transversal and oblique field (T, O): Having a magnetic field with finite transversal component crucially changes the scenario as TVF, and SPI structures are replaced by modified flow states due to the presence of such a magnetic field. In general, the isosurfaces become wavy-like deformed, and additional Fourier modes are stimulated as described before (Section 3.3; see also schematics in Figure 3). For pure transversal field, this mode stimulation results in wTVF with additional m=±2 modes aside the m=0 (Figure 6(c)) and in wSPI with m=n±2 modes on the secondary diagonal(s) in addition to the m=n SPI modes on the diagonal (Figure 6(d)). Beyond these, an oblique field further stimulates the modes m=±1 in wTVF and m=n±1 in wSPI due to more complex nonlinear dynamics/mode interaction. In particular, the spatio-temporal properties of the wTVF generated due to symmetry breaking magnetic field, either transversal or oblique (Figure 6(c,e)), differ significant from those of the classic wTVF. The structure of classical wTVF rotates as a whole in azimuthal direction [27, 32, 33, 34, 35].
The wavy-like deformation on the toroidally closed vortices η introduced by a magnetic field with a finite transversal component is different. Regarding Figure 6(c,e), one sees that the tubes of isosurfaces of the azimuthal vorticity become strong modified and bent at certain azimuthal cross sections of constant θ in comparison with the pure TVF tubes. Classic wTVF arises by deforming the TVF such that the axial stack of closed vortex tubes is wave-like modulated, bending axially upward and downward. Then, this deformation structure of the vortices rotates as a whole with a characteristic frequency. In contrast to these classical rotating wTVF, the wTVF initiated due to the symmetry-breaking transversal field component is stationary and does not rotate in azimuthal direction, instead, developing an 2-fold belly structure. It is worth to mention that such 2-fold symmetry also applies to the basic state. Thus, the classical rotational, axisymmetric CCF is replaced by the 2-CCF due to the symmetry breaking in the transversal magnetic field (see Section 3.2).
Consider wSPI, one can identify qualitative similar wavy characteristics, whereby the topological differences are weaker in comparison to wTVF. As the classical SPI structures already rotate without magnetic field, this also remains for wSPI in the presence of magnetic fields. However, also wSPI remains structural almost identical, and their modulation amplitude becomes enforced at certain cross sections of constant θ in presence of a transversal field component.
The basic state (CCF) loses stability either to (w)TVF (Figure 5) or (w)SPI (Figure 6). In case of wTVF, a supercritical pitchfork of revolution bifurcation breaks the axial translational symmetry Tα (8.3) and parameterizes the axial location of the steady axisymmetric family of solutions. On the other hand, (w)SPIs appear in a supercritical Hopf bifurcation which breaks the mirror symmetry Kz (8.2), resulting in a pair of symmetrically related spiral states, one with left-handed and the other with right-handed winding. It is worth to mention that a finite transverse magnetic field T breaks the axisymmetry Rφ, resulting in a basic state with discrete symmetry Rπ, i.e., with azimuthal wave number m=2 (section 3.3.1). Note that this m=2 state is stationary, which is in contrast to a generic breaking of an SO2 symmetry, resulting in a rotating wave.
3.3.4 Stabilization of basic state: point of higher co-dimension
For TCS, a crucial point of special interest in the Re1−Re2 phase diagram is the bicritical, co-dimension two point, γ, at which the primary bifurcation thresholds of TVF and SPI intersect/converge and as thus exchange the stability of the primary bifurcating solution [20]. For instance, in absence of a magnetic field and given wavenumber k=3.927, the γ-point is given by Re1≈95.25 and Re2≈−73.69 [28] (Figures 7(a) and 8(a)); for smaller Re2 (stronger counter-rotation), SPI bifurcate stable, and for larger Re2 (less counter-rotation and co-rotation), TVF bifurcate stable respectively with increasing Re1 (Figure 8(a)). Thus, with increasing Re1, the 2nd solution bifurcates as unstable solution. Note that the γ-point is quite sensitive to the axial wavenumber; for different k, it moves toward higher Re1 and smaller Re2. This bifurcation scenario also holds for wTVF and wSPI in oblique and transversal magnetic fields.
Figure 7.
(a) Influence of an axial magnetic field on the location of the bifurcation thresholds of TVF and SPI out of CCF in the Re1−Re2 diagram. Blue (full) lines refer to TVF and orange (dashed) ones to SPI. Magenta line with points indicates the moving of the γ-point with sz towards smaller Re2. Further shown are the bifurcation thresholds Re1,c with increasing sz2 for fixed (b) Re2=−100 and (c) Re2=0, respectively, obtained by full non-linear simulations.
Figure 8.
Bifurcation curves for the onset of (w)TVF (blue) and (w)SPI (orange) in Re2−Re1 diagram under (a) no magnetic field (N) sx=0sz=0, (b) an applied transverse field (T) sx≠0sz=0, (c) an applied axial field (a) sx=0sz≠0, and (d) an applied oblique field (O) sx≠0sz≠0. Simulation results of full nonlinear Eqs. (11) consider short periodic TCS with k=3.927 (cf. [6, 28]).
It is worth to mention that the curves in Figures 7(a) and 8 present “marginal stability” thresholds. These are obtained for a fixed wavenumber k=3.927, which differ from the ‘critical’ thresholds. The latter is defined by the (smallest) Reynolds number, for which the flow structure first appears, which crucially depends on k. Thus, the intersection points of the curves for TVF and SPI (e.g., Figure 7(a)) are critical values, with respect to the fixed wavenumber.
Figure 7(a) illustrates an example for the influence of different axial magnetic field strength sz on the bifurcation thresholds for TVF and SPI and the variation on the location of this bicritical γ point. These are results of linear calculations and were obtained with a shooting method [7]. Solid blue (dashed orange) lines refer to TVF (SPI). Note that always the lower (in Re1,c) bifurcating vortex solution is stable and the upper one is unstable. As discussed in Section 3.2, increasing sz stabilizes the CCF basic state, both, against TVF and SPI, i.e., the upward shift of the respective threshold curves. In general, the strength of the stabilization effect is larger for helical SPI vortex structures than for toroidally closed TVF (e.g. Figure 8). As a result, the γ-point moves toward more negative Re2 (i.e., toward the left in Figure 7(a)). Thus, it is observed that the Re2 region of primary, stable bifurcating TVF expands with increasing sz.
The same scenario also applies for transversal and oblique magnetic field configuration, whereby the magnitude of the shift is biggest for a pure axial magnetic field and almost equal for a pure transversal and oblique magnetic field (slightly larger for the latter). Figure 7(b,c) displays the bifurcation thresholds Re1,c for fixed outer cylinder rotations, Re2=−100 and Re2=0 with increasing sz2, resulting from fully nonlinear simulations. The slopes in Figure 7(b,c) depend crucially on Re2. However, in general, the slopes for SPI are steeper than for TVF. For small to moderate field parameters, TVF is the primary stable bifurcating structure while SPI bifurcate secondarily and unstable. For stronger fields, the structures interchange stability and bifurcation order.
Summarizing, for fixed Re2, an axial field shifts the onset of helical SPI more strongly than the one of toroidally closed TVF. We furthermore found this behavior to appear for all parameters and magnetic fields: either for axial, transversal, and oblique magnetic fields (Figure 8). In this case, that for transversal and oblique fields, all structures are wavy-like modulated; wTVF and wSPI, respectively (see also Section 3.2).
In pure axial or pure magnetic fields, the upward shifts of the bifurcation values, the stabilization of the basic state (CCF or 2-CCF), i.e., the onset of the primary stable flow state (w)TVF or (w)SPI grows linearly (Figure 7(b,c)) with the squared magnetic field parameters sx2 or sz2, respectively [6]. Thereby, the magnitude of stabilization shift depends crucially on Re2, in general is largest close to outer cylinder at rest.
3.3.5 Marginal stability thresholds: Stabilization of the basic state
Before going into detail about the stabilization of the basic state it is worth, once more to mention, that the basic state can be different in presence of a magnetic field, whether the field contains a symmetry breaking transversal component or not. For pure axial field, one finds the classical CCF which becomes replaced by the 2-fold symmetric 2-CCF for fields with finite transversal component (see also Section 3.3.1).
Figure 8 presents bifurcation thresholds for vortex flows under the influence of different magnetic fields. Shown is the Re2−Re1 diagram, along which the CCF state loses stability. For comparison, (a) shows the thresholds in absence of any field sx=0sz=0, while the same thresholds are shown for (b) a transversal sx≠0sz=0, (c) an axial sx=0sz≠0, and (d) an oblique field sx≠0sz≠0. These results were obtained by nonlinear simulations with the full field ferrohydrodynamical Eqs. (11), considering a short periodic TCS with k=3.927. See also Figure 7, for linear stability analysis, in which the shift of the higher co-dimension γ-point is shown together with stabilization of CCF in an axial field.
For all bifurcation thresholds of TVF, SPI, wTVF, and wSPI upward shifts of the bifurcation values, Re1,c, as shown in Figure 8, are observed. However, the exact onsets of the different primary bifurcating structures differ, either depending on the particular applied field as well as on Re2.
In absence of a magnetic field (Figure 8(a)) and sufficiently strong counter rotation, Re2≤−73 [28], the basic state, CCF, first becomes unstable to a helical solution, SPI, in a supercritical Hopf bifurcation breaking Kz, which results in a pair of symmetrically related, degenerated spiral states, one with left-handed winding and the other with right-handed winding. For considered Re2-parameter range in Figure 8, the bifurcating SPI has an azimuthal wavenumber m=1 (or m=−1 in case of the right-handed SPI). Note that the wavenumber m of the primary bifurcating SPI increases with stronger counter-rotating cylinders [20]. On the other hand, for Re2≥−73, the basic state loses stability to the toroidally closed TVF. The supercritical pitchfork of revolution bifurcation breaks Tα, and the steady axisymmetric family of solutions are parametrized by their axial location.
Considering a pure axial magnetic field, the scenario remains mainly unchanged (Figure 8(c)); still SPI and TVF appear as primary supercritical flow states with an intersection of the bifurcation thresholds at a specific value Re2. As already seen in Figure 7, increasing field strength sz in a pure axial field tends to shift the co-dimension γ-point toward more negative Re2. The same also holds for increasing sx in a pure transversal field, with significant minor effect on Re2 variation. Thus, the magnetic fields tend to render the corresponding basic states more unstable to (w)TVF and more stable to (w)SPI, respectively.
Having a field with finite transversal component, either pure transversal (Figure 8(b)) or oblique orientated (Figure 8(d)), the primary instabilities, with TVF being replaced by wTVF, and SPI being replaced by wSPI, respectively, just as the basic states were modified from CCF to 2-CCF due to the presence of the magnetic field. Their onset Re1,c occurs at higher Re1 values for any given Re2 value compared to the onset of classical TVF and SPI in the absence of a magnetic field. Thus, the magnetic fields alter the basic states and make them more robust to instabilities. The level of stabilization is greater for a pure axial field compared to a pure transversal field, but strongest, when an oblique field is imposed (bifurcation values Re1,c grow with field parameters sx,sz).
An oblique magnetic field results in the strongest upward shift of bifurcation values Re1,c, and also, this shift depends on Re2, which interestingly is largest about Re2=0, outer cylinder at rest. Furthermore, the topological different flow structures of (w)TVF and (w)SPI are shifted about the same amplitude, with a slightly lager effect on helical (w)SPI for the considered Re2−Re1 diagram. This should be kept in mind regarding later consideration involving field modulation via internal dependence, presented in Section 5.
So, while a pure axial fields shift the bifurcation onsets more strongly (i.e. to larger values Re1) than pure transversal ones (Figure 8), the latter one has a stronger influence on the nonlinear spatio-temporal properties and dynamics of the bifurcating vortex structures (Figure 6). Combining both portions to an oblique field enforces both effects. One finds either larger shifts as well as stronger modulated flow structures (Figure 6(e,f)).
3.4 Resume
An external applied magnetic field (independent of its orientation) stabilizes the CCF (2-CCF) basic state: for a given value of the outer Reynolds number Re2, the bifurcation thresholds Re1,c for vortex structures are shifted to higher values of the inner Reynolds number Re1. This holds likewise for toroidally closed (w)TVF, as well for helical (w)SPI. Thereby, these shifts are typically linear in the squared field parameters sx2 or sz2, respectively. Moreover, stabilization is in general stronger for helical solutions (w)SPI than for toroidally closed ones (w)TVF. Additionally, an external applied magnetic field can, depending on its orientation, change the spatio-temporal structure of the classical TVF and of SPI. The reason is the stimulation of additional modes entering the axial and azimuthal Fourier decomposition of the flow. A pure axial magnetic field keeps the symmetry properties unchanged and does not change the flow structures. It only alter the flow amplitudes, i.e., the magnitude of up-shift in the onsets and magnitude of modes after bifurcation. On the other hand, as soon as a finite transversal field component is present (pure transversal or oblique field), all flow structures become always 3-dimensional and therefore qualitatively different to the classical TVF and SPI states. In particular, they are wavy-like modulated and the pure structures do not exist anymore. Such a field component breaks the system symmetry and therefore alters either basic state (CCF changes to 2-CCF) and supercritical instabilities. As a result, one finds new solutions, field generated, nonrotating wTVF, contrary to the classical rotating ones. The vortex tubes of these wTVF are periodically expanded and constricted in azimuthal θ-direction, and this deformation pattern is stationary. In general, the bifurcation thresholds Re1,c for non-axisymmetric, helical m=1 (w)SPI are stronger modified in presence of a magnetic field than those for toroidally closed m=0 (w)TVF. The upward shift towards larger values of Re1 of the former is bigger than for the latter. One direct effect is the shift of the co-dimension point toward stronger counter-rotating frequencies. This effect has been found to show in axial, transversal, and oblique magnetic fields [6, 9].
Finally it is worth to mention that the magnetic field (transversal or oblique field) generated wTVF differ crucially from the classic wTVF [20, 33, 35, 45] that bifurcates secondary out of TVF at relatively high Re1,c in absence of a magnetic field. The present generated wTVF arises via a primary bifurcation directly out of the 2-CCF basic state.
4. Agglomeration and elongational flow effects
4.1 Introduction and motivation
Aside the general difficulties to incorporate the magnetization of any magnetic field into the governing equations of motion, it is also crucial to consider the specific properties of the magnetic particles. However, most mathematical models describing the flow of complex magnetic fluids typically assume noninteracting magnetic particles with a small volume or only point-size objects. However, real ferrofluids consist of a suspension of particles with a finite size in an almost ellipsoid shape. Moreover, they involve particle-particle interactions resulting in agglomeration; in particular, they tend to form chains of various lengths. One possibility to include such effects and to come close to the realistic situation for ferrofluids is to consider the effect of elongational flow. Mathematically, such flow is incorporated by the symmetric part of the velocity gradient field tensor, which could be scaled by a so-called, material dependent, transport coefficientλ2 [10, 37, 38]. As discussed in the previous section, magnetic fields tend to change the bifurcation threshold of the supercritical unstable flows and, depending on their orientation, resulting in structural modifications. Consider elongational flow, these modifications are expected to change the obtained states at least qualitatively.
4.2 Modified internal magnetic field
In order to consider such elongational flow (i.e. agglomeration, chain formation, etc.), one has to determine the relationship between the magnetization M→, the magnetic field H→, and the velocity u→. One possibility is to consider an additional dependence of the magnetization on the symmetric part of the velocity gradient S=∂iuj+∂jui/2 [10, 38]. Thus, the magnetization (4) becomes modified to
M→−M→eq=cNG→E11
where G=Ω→×H→+λ2SH→.
Onsager symmetry relations then require [12, 39] that additional forces involving λ2 appear in the off-equilibrium momentum balance of the ferrofluid.
Now, using the latest modified magnetization (9), the magnetization part in (1) can be eliminated, similar to previous discussion in Section 2.2. We thus have the following ferrohydrodynamic equation of motion (including elongational flow effects) [12]:
∂t+u→⋅∇u→+∇pM−∇2u→=−cN2H→∇⋅G→+H→×∇×F→/2E12
where pM is the dynamic pressure incorporating all magnetic terms which can be written as gradients, i.e., ∇H→⋅H→ and ∇H→eq−H→⋅M→. See also (5) for comparison without considered elongational flow.
Assuming that the internal magnetic field is equal to the external imposed magnetic field (it is known as a leading-order approximation [24], see also following section 5) which is sufficiently good for a ‘first’ numerical investigation for the effect of elongational flow. Therefore, (10) can be simplified:
In this approach, the magnetic field and all the magnetic properties of the ferrofluid could be influenced by the velocity field, the transport coefficient λ2, and the magnetic field Niklas parameters sx and sz as described before (7). However, it is worth to mention that in this approach, the material dependent transport coefficient λ2 is a parameter [12], which can be scaled by the term describing the elongational flow and is independent of H→. Thus, it only appears in combination with the symmetric component of the velocity field tensor. However, real ferrofluids show H→ dependence in the microscopic aspect of the dynamics of ferrofluids with chains or nonspherical particles [25].
4.3 Linear stability behavior and analysis
Before we come to the influence of λ2 onto the fully nonlinear Eqs. (11), it is worth to have a look at the linear marginal stability threshold of the different flow structures. Therefore, in the following, we first investigate the linear growth behavior of either stationary axisymmetric m=0 (TVF) and nonaxisymmetric m=1 (SPI) oscillatory vortex flow of a ferrofluid in the TCS that is placed in homogeneous axial magnetic fields H→ext=Hze→z (Figure 1). In the literature [4, 5, 7, 40, 41, 42], one can find many theoretical works that analyze (via linear stability analysis) the influence of rotational symmetric magnetic fields, i.e., axial, azimuthal, and radial ones on the flow of a ferrofluid in the TCS. However, for a symmetry breaking (e.g., transverse), magnetic field such linear stability analysis is not known for the corresponding basic state at such a field configuration. For a detailed description regarding the linearization around the CCF and solution of the linearized NSE, we refer to the Appendix in [43].
4.3.1 Marginal stability thresholds
In order to discuss how the transport coefficient λ2 influences the growth of m=0 and m=1 vortex perturbations, we determine the stability of the CCF basic state against these perturbations, considering a magnetic field H→ext=Hze→z, applied in pure axial direction.
Figure 9 presents the bifurcation thresholds in the Re2−Re1 diagram for either axisymmetric, m=0, and nonaxisymmetric, m=1, vortex flow in magnetic fields at (a,b) Hext=67.7kA/m and (c) Hext=135.4kA/m for two different axial wavenumbers. Note that the representative axial wavenumber chosen in (a) k=2.8274 is as detected in experiments [10]; (b) k=3.4558. Insets in Figure 9(b,c) show corresponding SPI frequencies at the bifurcation thresholds. As already discussed in Section 3.3.4, the CCF basic state is stabilized against the two types of vortex perturbations in any magnetic field (the lowest thick solid curve in each plot in Figure 9 refers to the case without any magnetic field H = 0 (c.f. Figure 8(a) for non-linear calculation)). Consider elongational flow (λ2≠0), the basic state becomes destabilized in comparison with the case λ2=0. Regarding Figure 9(a,b), both bifurcation thresholds for m=0 and m=1 vortices are shifted to lower values of Re1. Thus, for the axial wavenumber of k=2.8274, the λ2 term reduces the general stabilization effect caused by an applied magnetic field. For parameters as presented in Figure 9(a,b), this holds for the whole range of presented Re2, even when the changes of the stability boundaries are quite weak, in particular for corotating cylinders. It is important to mention that it never over-compensate the field induced stabilization. Figure 9(c) illustrates in general a strong parameter dependence of the elongational flow effect. While at Re2<22 for TVF and Re2<4 for SPI, increasing λ2 also has destabilizing effect, and for corresponding larger Re2, the effect is just opposite. Increasing λ2 even enforce the stabilization due to the applied magnetic field.
Figure 9.
Stability thresholds of the CCF basic state against vortex growth of (a) m=0 (TVF) and (b) m=1 (SPI) in the Re2−Re1 diagram of control parameters. Thin solid and dashed lines with symbols show the results for λ2=0 and λ2=0.8 for a magnetic field Hext=67.7kA/m. The axial wavenumber is k=2.8274. (c) Same curves, with wavenumber k=3.4558 and magnetic field is Hext=135.4kA/m. The lowest thick curve in (a) and (b) refers in each figure to the threshold in absence of a magnetic field Hext=0. Insets show corresponding SPI frequencies at the bifurcation thresholds.
It is worth to mention that, in contrast to the stabilization of the CCF state by a magnetic field, the λ2-induced modifications of the linear growth of vortices depend significantly on their axial wavenumber k as will be discussed below.
4.3.2 Magnetic field dependence of thresholds
Figure 10(I) illustrates the variation of the bifurcation thresholds m=0 (TVF) and for m=1 (SPI) with the magnetic field H→ext=Hze→z, incorporating the elongational flow effect. The top row of figures refers to Re2=0, and the bottom one to Re2=−100 and wavenumbers are k=2.8274 (left) and k=3.1415 (right), respectively. Ignoring the effect at small Hext, all thresholds increase linearly with Hext. The slopes are given by ∂Re1,stabH/∂H and vary slightly either with λ2 and k. For outer cylinder at rest, Re2=0, (Figure 10(I)(a)) the slopes of the m=0 and m=1 curves are comparable, they differ significantly for counter-rotating cylinders, Re2=−100 (Figure 10(I)(b)). For the latter, this results in intersections of the m=0 and m=1 stability boundaries. Such intersections have been observed also in experiments and also in full nonlinear simulations [6] with λ2=0.
Figure 10.
(I) Variation of the (a) m=0 (TVF) and (b) m=1 (SPI) bifurcation thresholds and (II) corresponding SPI frequencies ω with magnetic field Hext for different λ2 as indicated. The top row refers to Re2=0 and the bottom one to Re2=−100. The wavenumbers are (1) k=2.8274 and (2) k=3.1415, respectively.
In general, modifications of the Hext field dependence of the bifurcation thresholds due to λ2 tend to increase, when the cylinders are stronger counter-rotating (Figure 10(I)). Re2 becomes more negative, and thus, the strain increases. On the other side, they also depend on the axial wavenumber k (see also Section 4.3.4). Considering fixed k, the modifications due to λ2 increase with increasing magnetic field Hext while remaining qualitatively the same. However, the general tendency visible in Figure 11, considering elongational effect, is to shift the stability boundaries downward in Re1. Concrete amounts depend on the parameters k and Re2. Note that this relative small destabilizing elongational effect is always overcompensated by the significant stronger general stabilization effect of the magnetic field itself.
Figure 11.
(1) Stability boundaries of the CCF basic state against the growth of m=0 (TVF) (blue) and m=1 (SPI) vortex flow (orange) in the k−Re1 diagram for λ2=0 (solid lines) and λ2=0.8 dashed lines. (2) corresponding marginal spiral frequencies ω at the m=1 thresholds. Further parameters: (a) Re2=0 and (b) Re2=−100. The (axial) magnetic field is constant with Hext=80kA/m. Results from linearized equations [43].
4.3.3 Frequencies of helical SPI states
Instead of looking at the growth rates—the real part of the eigenvalue of these solutions (previous Sections), one can also consider the marginal spiral frequency ω—the imaginary part of the marginal eigenvalue at the bifurcation threshold of m=1 (SPI) vortex flow. Figure 10(II) shows the frequency, corresponding to the m=1 bifurcation thresholds presented in Figure 10(I). Main effect is the increase of the SPI frequencies with the field Hext, accordingly to the increase of onset Reynolds numbers with Hext (Figure 10(I)). Thus, the threshold frequencies are larger for Re2=−100 than for Re2=0. In general, the effect induced by elongational flow is minor, significant smaller, and only secondary. Analog to the slightly down shift in the onsets (Figure 10(I)) also the SPI frequencies are slightly reduced (Figure 10(II)). Strong counter-rotating cylinders reduce in general the frequency for finite λ2 (Figure 10(II)), while for larger wavenumbers, a finite elongational flow λ2≠0 at Re2=0 it can be slightly decreased.
A general observation is that all modifications of the spiral frequencies due to finite λ2≠0 are most pronounced for the smaller wavenumbers (here k=2.8274 in Figure 10) in counter-rotating cylinders. The field dependence is virtually not much modified due to variations in λ2, and the growth rate of the imaginary part of the marginal eigenvalue with increasing λ2 is quite similar to that of its real part.
4.3.4 Axial wavenumber dependence
Following, we will first discuss the situation for a (moderately) larger wavenumber, taking k=3.4558 as one specific, representative example. For larger value k (Figure 9(c)), the situation becomes more complicate. Incorporating the elongational flow scaled by the parameter λ2 increases and decreases the growth rate, with respect to λ2=0, and thus shifts the stability boundaries, up or down, depending on Re2.
Figure 11 illustrates the relation between the wavenumber dependence of the stability boundaries and of the spiral frequencies. This is of particular interest in experiments, since the wavenumbers of some vortex structures show a strong variation when changing the magnetic field [9, 44]. Reindl et al. [9] also observed hysteresis between vortex flows with different k when increasing and decreasing the Reynolds numbers Re1. The bifurcation thresholds (Figure 11, top row Re2=0 and the bottom row Re2=−100; the (axial) magnetic field is fixed at Hext=80kA/m) for m=0 and m=1 vortex flow structures in the k−Re1 diagram change when the elongational term λ2≠0 is considered in the linearized ferrofluid equations [43]. Right column of Figure 11(2) presents the wavenumber dependence of the corresponding marginal spiral frequencies at the onsets for m=1 SPI.
Typically, the thresholds for λ2≠0 (in Figure 11 the case of λ2=0, lines with squares) are shifted downward in Re1 at smaller wavenumbers k and upward at larger k. Resulting in a fixed magnetic field, the elongational effect destabilizes (stabilizes) the CCF state against growth of vortices with small (large) wavenumbers, either for m=0 or m=1 perturbations. Consequently, the stability onsets cross for different values of λ2 (λ2=0, λ2=0.8 in Figure 11) at a specific value of kcross the wavenumber.
It is worth to mention that the crossover values of kcross for m=0 and m=1 perturbations also depend sensitive on other the control parameters, e.g., Re2. Based on the observations in Figure 11, it is fair to assume that, depending on the wavenumber k, the γ-point (Section 3.3.4) will change.
Similar to the stability thresholds, also the spiral frequencies (Figure 11(2)) increase with increasing λ2 when k is large and decrease for smaller k as a result of switching on the elongational effect. Thereby, the crossing wavenumber, kcross, separating the increase and decrease in ω is slightly smaller than the one for crossing the marginal thresholds (Figure 11(a1,b1)). Following the speculation by Odenbach and Müller [10] on the physical origin of the elongational effect, one can give a reasonable explanation for such an observed wavenumber dependence of the λ2 effect on the growth rates of vortex flow. In [10], the authors argue that the microscopic origin of a finite λ2 might be attributed to the finite asphericity of the colloids, i.e., the presence of (short) particle chains in the ferrofluid. Small axial wavenumber k will allow these chains more easily to remain aligned in the direction of the magnetic field in vortex flows instead in a flow with large axial wavenumber and large azimuthal vorticity. Larger wavenumber means smaller wavelength and therefore more higher probability to unalign in the given field. Thus, since the magneto-rotational dissipation in the latter flow is larger than in the former one, it requires larger centrifugal forces, i.e., larger values of Re1, to drive the growth of a large k flow. However, it is important to state that this argumentation also holds some weakness: In the presence of chains, one could expect also for small k vortices an increase of the threshold Reynolds number Re1. On the contrary, incorporating the elongational term λ2≠0 in the magnetization balance changes the momentum balance of the vortex flow such that the growth of vortices is enhanced for k<kcross (left in Figure 10) and reduced for k>kcross (right in Figure 10) relative to the reference case of λ2=0.
4.4 Nonlinear analysis
4.4.1 Bifurcation thresholds
Following, we will consider the full nonlinear ferrohydrodynamical equation of motions incorporating elongational flow effects (11) in order to investigate the influence of a finite transport coefficient λ2 at different applied magnetic fields.
The general modification in the bifurcation threshold of primary instabilities ((w)TVF) and ((w)SPI) due to finite values λ2≠0 remains qualitative as detected and described before in the linear stability analysis in pure axial field. However, the absolute values are slightly different. A direct comparison for axial field (Figures 9-11) results in variation about one to two percent between the calculated onsets using either the linear analysis or full nonlinear equations [6]. While the different field orientations result in different magnitudes of stabilization of the basic states CCF and 2-CCF, respectively, as discussed in Section 3.3.4, a finite value λ2 does not depend on different field orientations. The variation effect with λ2 is almost identical for all considered field configuration. However, in general, the modifications are slightly larger under oblique magnetic fields. Depending on other system parameters, e.g., wavenumber k, cylinder rotation speeds Re1,Re2, the effect is either stabilizing or destabilizing with respect to the situation λ2=0 (Figure 10(c) and Figure 11). Considering the effect of elongational flow, the critical bifurcating point of primary centrifugal instability will be shifted either up or down. As detected in the linear analysis, the potential destabilizing effect of λ2never overcomes the general stabilizing effect of the external applied magnetic field, i.e., comparing λ2=0 and any value λ2≠0.
4.4.2 Flow structure modifications
Aside the modification of the primary stability thresholds, there are further flow structural changes/modifications which can be observed for finite value λ2. Figure 11 shows, for different field configurations T: sx=0.1sz=0, A: sx=0sz=0.1, O: sx=0.1sz=0.1, how the maxima of vorticity η vary for either the full flow structure Δη≔ηλ2≠0−ηλ2=0 (Figure 12(a)), and reduced 2-fold subspace (Figure 12(b)) varies with λ2 for the basic state configuration Re1=60Re2=0. Remember, that only for sx=0, one has the classical CCF, and for any sx≠0, a modified basic state is present [6, 9, 37] (i.e., symmetry breaking effect of transversal field). Although for a given configuration (A, T, O), the absolute maxima of these values are significantly different, and their relative differences Δ scale is linear with λ2 starting at zero. Lines in Figure 12 are calculated by linear regression fitting to the numerical data points. The effect of λ2 is decoupled from the magnetic field configuration and with this the considered flow structure.
Figure 12.
Variation of vorticity maxima with λ2 for (a) maxη=ηλ2−ηλ2=0 at full solution and (b) maxηm=2=ηλ2−ηλ2=0 at only 2-fold symmetric subspace (empty for A) for basic flow (CCF and 2-CCF, respectively) at Re1=60,Re2=0 and different field orientations: Axial (A) sx=0.0,sz=0.6; transverse (T) sx=0.6,sz=0.0 and oblique (O) sx=0.6,sz=0.6. Lines are linear fits.
Either for classical “pure” CCF or field modified 2-CCF, the variation in maxΔη grows linear with λ2 (Figure 12(a)) (note the significant smaller amplitude in the m=2 subspace). Further confirmation for this is the fact that an analog linear growth can be also found when considering other parameters, at which a supercritical flow state exists [37]. This holds for toroidally closed structures (w)TVF, as well as for helical flow states (w)SPI. Only the absolute values in maxΔη are significant larger.
It is important to mention that the slopes in Figure 12 (for Osx=0.1sz=0.1) not only depend on the applied field orientation (O) but also depend on the strength of the applied magnetic field. With increasing the magnetic field strength sx,sz, the slopes become steeper. This scenario is illustrated, exemplary for 2-CCF basic state in transversal field configuration in Figure 13.
Figure 13.
As Figure 11 for basic state 2-CCF at Re1=60,Re2=0. Variation with λ2 of vorticity maxima maxη=ηλ2−ηλ2=0 for transversal field (T) at different field strength sx as indicated.
Moreover, similar characteristics hold for the different primary bifurcating structures. Figure 14 elucidates this effect for the three different field configurations A, T, and O as indicated (see also Figures 4 and 6 in [37] for another example considering a transversal magnetic field in finite system configuration). For a given field configuration, the slopes increase with λ2, when increasing the magnetic field strength sx. In general, slopes are smallest for transversal and largest for oblique fields, respectively. The variation in the slopes is continuous with increasing the field strength, but not linear; for larger field strength, the variation in the slopes becomes smaller.
Figure 14.
Supercritical (w)TVF solution at Re1=100,Re2=0. Variation of vorticity maxima maxη=ηλ2−ηλ2=0 for (a) transversal (T), (b) axial (A) and (c) oblique (O) magnetic field (strength as indicated). Lines are linear fits.
4.4.3 Wavenumber modification and selection
As seen in the former discussion (section 4.3.2), an elongational flow can have strong effect on the axial wavenumber k. Consider the full equations (and periodic boundary conditions), one obtains similar results as presented in Figure 10 for linear analysis of the marginal stability thresholds for TVF and SPI. As described before, the agreement between the results of linear and full equations is very close (they differ less than 2 percent).
However, in order to see such wavenumber modification from another perspective, it makes sense to study a finite system with sufficient large aspect ratio enclosed by fixed end plates (assumed to be stationary) on top and bottom. As in a finite system, the axial wavenumber k is natural selected by the system itself (depending on the various system parameters), it is possible to see the direct influence, considering a finite value λ2 on to the flow structures and their axial wavenumbers k.
Starting with a classical TVF as initial flow state, with n=22 vortices, which corresponds to an axial wavenumber k=3.41 in the bulk (illustrated on the most left border in Figure 15), we consider a situation with an applied pure transversal magnetic field which is continuously increased for two values, either ignoring any elongational flow effect, λ2=0 (Figure 14(a)), or considering such effect λ2=0.2 (Figure 15(b)) (note 0.2 is quite small). In any case, either for λ2=0 and for λ2=0.2, by increasing the magnetic-field strength, sx, the number of vortices is reduced, which leads to an increase in the axial wavelength (i.e. decrease in axial wavenumber k) due to elimination of vortices in the bulk. Therefore, the corresponding axial wavenumber k in the bulk becomes reduced. Remember that the initial state is a classical, stationary TVF, which becomes modified to stationary wTVF as soon as sx≠0. For the field generated wTVF, the wavy-like modulation is not visible in Figure 15. However qualitative, the same scenario also holds for either axial and oblique magnetic fields, with changes in the corresponding parameters sx,sz at which the wavenumber k changes, i.e., vortices disappear.
Figure 15.
Contours of the radial velocity component u at mid-gap with variation of field strength sx at Re1=88.35,Re2=0 for (a) λ2=0 and (b) λ2=0.2. Left border in each plot gives the initial state, a classical TVF with n = 22 vortices corresponding to an axial wavenumber k = 3.41 in the bulk. Note that for sx=0, only wavy-like modified flow states exist. Red (yellow) [dark gray (light gray)] corresponds to positive (negative) values. The max (min) level is ±9.98.
In general, the exact field strength sx for vortex elimination and/or the number of destroyed vortices depends on λ2. In the example of Figure 15, the initial flow has 22 vortices with an axial wavenumber k=3.41. Neglecting elongational effects λ2=0 (Figure 15(a)) until the parameter sx approximately reaches the value about 0.7, the flow remains stable, and the number of vortices in the flow can be reduced to n=20 with k=3.04. Further increasing sx to about 0.85, the number of vortices can be more reduced to n=16 with k=2.58. Hereafter, no further changes in the number of vortices n appear, before reaching the boundary threshold, sx≈0.9. Consider elongational flow λ2=0.2 (Figure 15(b)), the behavior of the vortex state is similar to the case neglecting such effects. The flow remains stable before the parameter sx reaches about 0.75. At the same time, the number of vortices is also reduced in analog manner to n=20 with k=3.04. Hereafter, as before, the number of vortices n remains unchanged until the parameter sx increases to the boundary threshold, sx≈0.95. Consider values λ2>0.7, the flow structure (for given parameters) is not affected by the magnetic-field strength. A flow with n=22 vortices corresponding to axial wavenumber k=3.41 in the bulk is remaining in the bulk.
An interesting phenomenon, which also has been reported in experimental studies [9] as in numerical simulations, is the fact that one or two vortex pairs can be eliminated. However so far, there is no physical explanation for this observation/phenomenon. Analog to the presented scenario in Figure 15 for pure transversal magnetic field, one also finds similar reduction scenarios for pure axial and oblique magnetic fields. As to expect, the absolute values of field strength, sx,sz, at which the elimination of vortices happen is different, depending on the applied field.
4.5 Resume
Real ferrofluids consist of a suspension of particles with finite size and mainly almost ellipsoid shape, as well as there are particle-particle interactions that tend to form chains of various lengths and tend to agglomerate in general. One possibility in order to make a more realistic approach for ferrofluids and to come close to such a configuration in a real ferrofluids is to consider the effect of elongational flow. Such a flow becomes incorporated by the symmetric part of the velocity gradient field tensor in the nondimensional ferrohydrodynamic equations of motion (11) and scales by a material dependent transport coefficient λ2. By such an assumption, the latter is given by a pure scalar number, describing the strength/influence of all former mentioned properties.
Consider either linear stability analysis or fully nonlinear equations, one finds that incorporating such elongational flow can have either stabilization or destabilization effect (i.e., modification in the supercritical primary bifurcation thresholds) with respect to the reference situation in absence of such flow λ2=0. The modifications due to finite value λ2 are qualitative similar and in the same range for different considered magnetic fields, slightly largest in case of an oblique field compared to pure axial or transversal ones. Typically, the modifications due to λ2≠0 are larger for helical, nonaxisymmetric, m=1, (w)SPIs than for toroidally closed, axisymmetric, m=0, (w)TVF. Thus, the regim of primary stable bifurcating (w)TVF increases. This observation is similar and goes in hand with the previously described general stabilization of any applied magnetic field (Section 3). There we saw that the stability thresholds of helical flow states are shifted further than those of toroidal closed flow structures. However, it is important and fundamental to state that the elongational effect λ2≠0 is never big enough to over compensate the general, field-induced stabilization. Only the magnitude/amplitude in shift can be reduced. Further, it is worth to mention that linear stability analysis is limited to axial magnetic fields; as for fields with finite transversal component, no basic state can be derived, simulations of full nonlinear equations confirm similar effects. The modifications due to λ2 are sensitive to various different system parameters (e.g., Reynolds numbers, radius ratio, etc.), in particular sensitive to the axial wavenumber k. Here, smaller k are strongly effected than larger ones. Thus, considering finite size TCS, a considered elongational flow results in modifications of flow structures in the sense of changing number of vortices in the bulk, which coincides with a different wavenumber selection. Thus, depending on the particular wavenumber, the presence of any elongational effect λ2≠0 can be either stabilizing or destabilizing for the basic state. In finite system, this results in different selection and/or modification of number of vortex pairs within the annulus, depending on the parameter.
Regarding the flow structural properties, a finite transport coefficient, λ2≠0, does not change the qualitative shape of the flow pattern. Only quantitative differences can be observed. Essentially, all structural properties remain as generated/induced due to the external applied magnetic field. Thus for instance, the azimuthal mode-2 symmetry originated in the symmetry-breaking transverse magnetic field is preserved; only the strength of different flow quantities is modified (in general they increase) according to finite λ2. Mainly, the flow states become linearly enforced, e.g., speaking in quantity of the azimuthal vorticity, which holds for both the basic flow (CCF, 2-CCF) and the primary instabilities the supercritical centrifugal unstable flows ((w)TVF, (w)SPI).
5. Modified field dependence - internal magnetization
5.1 Introduction and motivation
Aside the properties and behavior of the magnetic particles itself and their interaction with the applied magnetic field (see last sections), the correct assumption of the applied magnetic field is another crucially point in describing magnetic flows and their dynamics. The frequently used, most common assumption is to consider the internal magnetic field within the magnetic fluid (e.g., a ferrofluid) to be equal to the external applied field. Following, we will see that this is only a leading-order approximation. By accounting for the ferrofluid’s magnetic susceptibility, a uniform externally imposed magnetic field is modified by the presence of the magnetic fluid (ferrofluid) within the annulus. The modification to the magnetic field appears in an radial dependence 1/r2 and their magnitude scales with the susceptibility χ. To be concrete, such modification to the imposed magnetic field can be substantial, also for ferrofluids typically used in laboratory experiments. These have significant consequences on the structural properties and stability of the basic states, as well as on the supercritical (primary) bifurcating solutions.
5.2 Modified internal magnetic field
As already seen in the previous sections, for solving the equation of motion (1), we have to consider the magnetic field within the annular gap between the cylinders. In the last sections, we used the simplest approach with the assumption to take the magnetic field to be identical to the applied external field [4, 6]. However, this simple approach is only a leading-order approximation, and depending on the magnetic susceptibility of the ferrofluid, the magnetic field in the gap is modified.
Assuming infinitely long cylinders, the magnetic boundary conditions are
H→=H→ext−Mre→rE14
at r=r1 and r=r2, where H→ext is the homogeneous external applied magnetic field in the absence of the ferrofluid-filled annulus and Mr is the radial component of the magnetization M→.
The following ansatz satisfies a solenoidal field, ∇⋅H→=0:
and then the boundary conditions defining the external field are.
H→r=r1=H→ext+a1−b1/r12cosθ+a2−b2/r12sinθe→r+a2+b2/r12cosθ−a1−b1/r12sinθe→θ and H→r=r2=H→ext+a1−b1/r22cosθ+a2−b2/r22sinθe→r+a2+b2/r22cosθ−a1−b1/r22sinθe→θ.
To calculate the four unknown constants a1,a2,b1,b2 from the boundary conditions (12), we need the radial component of the magnetization which results from substituting (13) into the magnetization (3):
where HextT stands for the transverse component of H→ext. Using the continuity of the magnetic field condition on the cylinders, the coefficients can be calculated. Substituting them into (13), finally the resulting magnetic field is given by
where K=2+χ2−χ2η2. A special case of (15) when ΩcN=0. This corresponds to an equilibrium magnetization and has been already observed in an earlier work by Odenbach and Müller [38]. They considered the limit of an equilibrium magnetization and used a similar derivation of the modification of the magnetic field for a ferrofluid between two cylinders.
Based on (15), one can use the following magnetic field in the equation of motion:
H→=22+ζKHx[1−ζ/r2cosθe→r−ζ/r2sinθe→θ+Hze→zE18
the dimensionless scalar parameter
ζ=χ2+χr12E19
characterizes the strength of the radial field dependence ζ. The ferrohydrodynamical equation of motion (6) and the Niklas parameters sx and sz (7) remain unchanged/untouched from such considered radial field dependence.
5.3 Flow structural modifications
5.3.1 Primary instabilities
As we learned in previous sections, a magnetic field with a finite transverse component modifies either basic state or primary bifurcating flow structures in a wavy-like manner. More details of these states have been reported in numerical simulations [6], for which a homogeneous internal magnetic field equal to the external applied one has been supposed, as well as in real experimental studies [9]. The question now is how does the former derived radial field dependence ζ affects the stability of the basic states (CCF and 2-CCF) and the characteristics of the bifurcating solutions ((w)TVF and (w)SPI).
We begin our discussion by considering a purely axial applied magnetic field (A), with fixed field strength sz=0.6, which leaves the flows, in particular its structural properties unaffected (TVF and SPI). Figure 16 shows the influence of the radial field parameter, ζ, (17) (exemplary chosen (a) ζ=0, (b) ζ=0.4 and (c) ζ=0.8), on the bifurcation thresholds for TVF and SPI in an axial magnetic field. As in absence of any field dependence, ζ=0 (Figure 16(a)), the bifurcations continue being supercritical. Main effect on the bifurcation curves is that with increasing ζ, the thresholds for the onset of both TVF and SPI are shifted to higher Re1 (Figure 16(b,c)). Thereby, the threshold for the toroidally closed TVF being shifted significantly further (to larger values Re1). Note the different scaling on the abscissa in Figure 16(c). Aside from the enhanced stabilization of the CCF basic state on top of the field introduced stabilization itself, one observes an important consequence to the relative shifts in the two bifurcation curves. Consequently, the bifurcation curve of TVF is shifted more than the one of SPI, the co-dimension-two point, γ, is also shifted to more positive Re2 (c.f. Figure 15(a,b)). Interestingly, this is just the opposite way of the shift with increasing field strength sz in a pure axial magnetic field (the same holds for sx in transversal field), which we saw in Figure 7 with moving the γ-point to the left toward more negative values in Re2. Consider larger values, e.g., ζ=0.8 (Figure 15(c)), one sees that the co-dimension-two γ-point is even shifted to Re2>50 out of the considered parameter space. As a result, the primary instability is always the helical SPI, even with the two cylinders in strong co-rotation. This observation is crucially different from classical TCS, both without magnetic fields and also with assumption of an applied uniform magnetic field. Thus, the radial field dependence, induced by the susceptibility of the ferrofluid-filled annulus, is seen to cause major changes in the quantitative and qualitative characteristics of the flow instabilities.
Figure 16.
Bifurcation curves for the onset of TVF (blue) and SPI (orange) in Re2−Re1 diagram for different ζ-values (a) ζ=0, (b) ζ=0.4 and (c) ζ=0.8 with an applied axial magnetic field (A) with sz=0.6. ((a) is a replot of Figure 7(c) for comparison). Note the different scale on the ordinate in (c). Simulation results of full nonlinear equations (11) consider a short periodic TCS with axial wavenumber k=3.927 (same hold for Figures 15 and 16).
Consider a transversal magnetic field, field strength sx=0.6 (Figure 17) instead of an axial one, the variation by increasing ζ remains virtually and qualitatively the same, aside that now all structures are wavy-like modulated (2-CCF, wTVF, wSPI). One finds the same preference, in the sense of stronger upshifts for toroidally closed wTVF with the same result of a shift in the γ-point to the right to larger values in Re2. Thus also here, for sufficiently strong ζ, e.g., ζ=0.8 in Figure 17(c), the helical wSPI is the primary stable bifurcating solution in all presented parameter space. Comparing both, pure axial and pure transversal field (at same field strength), the up-shift, i.e., the stabilization is slightly larger in axial fields. This holds for toroidally closed (w)TVF and helical (w)SPI solutions.
Figure 17.
Bifurcation curves for the onset of wTVF (blue) and wSPI (orange) in parameter space for different ζ values (a) ζ=0, (b) ζ=0.4, and (c) ζ=0.8 with an applied transversal magnetic field (T) sx=0.6. ((a) is a replot of Figure 7(b) for comparison).
Finally, this effect is even more pronounced when an oblique magnetic field is applied as illustrated in Figure 18 for field strength sx=0.6,sz=0.6. Note also the different scale on the ordinate in Figure 17(c) highlighting the stronger shifts due to such field configuration. One can speculate that this enforcement is due to the axial component of the oblique field leading to a base state with w having axial shear, and the radial field dependence ζ enhances the w component of the 2-CCF base state. The combination of the base state azimuthal shear and this axial shear favors the helical instability. Close to the onset, the variation of the maxima of all flow components u,v,w is almost linear.
Figure 18.
Bifurcation curves for the onset of wTVF (blue) and wSPI (orange) in Re2−Re1 diagram for different ζ-values (a) ζ=0, (b) ζ=0.4, and (c) ζ=0.8 with an applied oblique magnetic field (O) with sx=0.6=s−z. ((a) is a replot of Figure 7(d) for comparison). Note different scaling in (c).
5.3.2 Modified flow structures
Figure 19 illustrates the flow structural properties of wTVF and wSPI in either a transverse magnetic field T (a,b) and an oblique magnetic field O (c,d), both with assumed radial field dependence ζ=0.8. Note, the flow structures of TVF and SPI under a pure axial applied magnetic field are not shown, as they virtually do not differ from the ones in absence of a magnetic field. Only the absolute values (e.g. in η) are different. Comparing these flow structures to the ones for ζ=0 cases as earlier presented in Figure 6 (Section 3.3.3), one sees the modulations in all structures to be enhanced/enforced. In general, the waviness (axial variation in azimuthal direction) of all flow structures is enforced/enhanced by the radial dependent magnetic field. However, the field-induced wTVF always remains a stationary, nonrotating, and phase-pinned structure, irrespective of whether the radial field dependence is considered or not. On the other hand, wSPIs bifurcate as degenerated solutions, either as a left- or right-winding SPI, depending on initial conditions. Their corresponding frequencies remain mainly constant over a wide range and do not differ much (about 20% or less) from those of SPI. Comparing the stimulated modes for flow structures when ζ=0.8 (bottom in Figure 19) to those with ζ=0 (bottom in Figure 6), one sees general larger mode amplitudes in the additional stimulated modes for consideration of an internal field dependence ζ≠0. These larger amplitudes in the field induced modes are also the reason for the increase in waviness.
Figure 19.
(See also Figures 2 and 6 for comparison to flow states without radial field dependence. The radial field dependence is considered to be ζ=0.8. Top row: Azimuthal vorticity isosurfaces η over two axial wavelengths with pure transversal applied magnetic field (T): (a) wTVF at Re2=0, Re1=−130, (b) wSPI at Re2=−150, Re1=190; and with oblique applied magnetic field (O): (c) wTVF at Re2=0, Re1=−150, and (f) wSPI at Re2=−150, Re1=−200 Red (yellow) isosurfaces correspond to positive (negative) values as indicated. The bottom row shows the mode amplitudes um,n of the radial velocity field u corresponding to the structures above over the m-n-plane.
5.4 Resume
Applying a magnetic field across a container of ferrofluid changes the ferrofluid’s susceptibility and results in modification of the magnetic field structure within the container. So far, most theories commonly neglect such modification, either for convenience or simplicity, and at the end, there is also possibility of numerical simulations in appropriate time or feasibility at all.
However, it is a matter of fact that an external applied magnetic field is modified, depending on the magnetic susceptibility. Experimental studies [10, 38] confirmed the modifications also for susceptibilities of commonly used ferrofluids. Consider such modifications to the imposed magnetic field effects the basic state, its stability, as well as the primary (supercritical) bifurcating solutions. This holds for any field configuration, with only differences in the strength of modifications, i.e., the stabilization (its magnitude; up-shift, move of the bifurcation thesholds to larger control parameters). On the other hand, the flow structural properties remain qualitative unaffected by such an radial field dependence, only change in quantitative manner. They are determined only by the general kind and orientation of the applied magnetic field. Thus, a pure axial field does not have any effect on the flow structural properties, leaving them as in absence of any field, i.e., pure TVF and pure SPI. As described before, the presence of a finite transversal field component results in a wavy-like modulated flow state. This is independent/untouched if the magnetic field is considered to be homogeneous or inhomogeneous with a radial field dependence. However, consider stronger modifications to the imposed magnetic field, i.e., larger values of the radial field dependence parameter ζ have the effect of further strengthening the already present waviness of the structures. This holds for (w)TVF and (w)SPI alike. An important observation is the different influence of the radial field dependence ζ≠0 onto topological closed and helical flow states. In general, the effect of modification in the primary bifurcation thresholds is larger for toroidally closed flow structures compared to helical ones. As a result, their intersection point, i.e., co-dimension γ-point is moved toward positive Re2, which is just the opposite direction, as it moves with increasing the magnetic field strength (for any given field configuration, e.g. Figure 7). The parameter space, the area for which (w)SPI bifurcate primary stable out of the basic state, is growing, even into region of co-rotating cylinders, i.e. Re2>0, making wSPI the favorable flow state in the system.
6. Conclusions
Although there has been significant progress in the understanding of ferrofluidic flows in the past decades, their fully theoretical understanding is still limited. A key reason for this is computational high-expensive simulations, which even increase with continuation in refining the models to make them more realistic/appropriate, e.g., finite size effects, internal field modifications, etc. However, these model improvements are essential and necessary to uncover the full potential of ferrofluids, which already are used in a wide field of applications. The use in the car industry within antishocks to provide adjustable damper settings is probably one of the most prominent examples. However, either experiments and numerics are essential to mutually prove each other and combined to provide new inside of potential applications.
This chapter deals with a rotating ferrofluid between two concentric cylinders, i.e., ferrofluid in Taylor-Couette geometry. Therefor we consider different improvements in the model describing such magnetic particles in a flow under external applied field. Concluding the discussions, the main points can be summarized as follows:
Any applied magnetic field changes the bifurcation thresholds of primary solutions—the thresholds of both primary solutions, TVF and SPI, are shifted to larger parameter values—the basic state becomes stabilized against perturbations.
Magnetic fields with a finite transversal component break the system symmetry and renders all flow states to be inherently 3D. As a result, the pure flow states TVF and SPI cease to exist and are substituted by their wave-like cousins, wTVF and wSPI. Thereby, the generated wTVF differs crucially from the classical ones, which rotate in azimuthal direction. Instead, the magnetic field induced wTVFs are nonrotational, phase-pinned flow solutions and develop belly shape structure at specific azimuthal positions.
Assuming elongational flow, in order to consider finite-size effect of the magnetic particles and their interaction (agglomeration, particle-particle interaction and chain formation procedures), one detects further changes in the bifurcation threshold of the solutions. Also, this can result in further stabilization or destabilization. However, the effect of such elongational flow is never big enough to over-compensate the general, field-induced stabilization. Only the magnitude/amplitude varies. The flow states remain structural as in absence of elongational flow.
Accounting for the changes in ferroluid’s susceptibility in the container due to an external magnetic field, the real internal field can be approximated by 1/r2. As before the flow, structural properties remain qualitative similar. Bifurcation thresholds are shifted upward, toward larger Re1, whereby toroidally closed (w)TVFs are in general preferred and shifted more than helical (w)SPI. The result is an increase in the regime of primary bifurcating (w)SPI.
Conflict of interest
The author declares no conflict of interest.
A. Appendix
A.1 Numerical method
There are various different ways to solve the ferrohydrodynamic equations of motions (Eqs. (5) and (10)), the latter including elongational flow). All here presented results have been obtained by a code (G1D3 [6, 7]), which combines a finite difference method of second order in (r,z) and time (explicit) with spectral decomposition in θ [6]:
frθzt=∑m=−mmaxmmaxfmrzteinkθ,E20
where f denotes one of uvwp. For the parameter regimes investigated here, mmax=8 provides adequate accuracy. For discretization, a forward time, centered space (FTCS) algorithm is used on staggered grids in the r−z plane following the procedure of Hirt et al. [45]. In particular, a homogeneous grid with discretization length δr=δz=0.05 and time steps δt<3800 has been used. Here, δt defines the time step between two iterations via FTCS algorithm of the system of coupled equations for the amplitudes fmrzt of the azimuthal normal modes −mmax≤m≤mmax. For diagnostic purposes, the complex mode amplitudes fm,nrt were calculated obtained from a second axial Fourier decomposition:
fmrzt=∑nfm,nrteinkz,E21
where k=2π/λ is the axial wavenumber.
Acknowledgments/Funding
S. Altmeyer is a Serra Húnter Fellow. This work was supported by the Spanish Government grant FIS2017–85794-P.
Nomenclature and abbreviations
TCS
Taylor-Coutte System
NSE
Navier–Stokes equation
CCF
Circular Couette flow [basic flow]
2-CCF
2-fold CCF [basic flow in transversal field]
TVF
Taylor vortex flow
SPI
Spiral vortex flow (spiral)
wTVF
wavy TVF
wSPI
wavy SPI
\n',keywords:"magnetic fluid, complex system, ferrofluid, instabilities, bifurcation, symmetry-breaking, vortex, growth rate, finite-size effect, internal magnetic field, particle interaction, agglomeration, elongational flow",chapterPDFUrl:"https://cdn.intechopen.com/pdfs/63437.pdf",chapterXML:"https://mts.intechopen.com/source/xml/63437.xml",downloadPdfUrl:"/chapter/pdf-download/63437",previewPdfUrl:"/chapter/pdf-preview/63437",totalDownloads:396,totalViews:0,totalCrossrefCites:0,totalDimensionsCites:0,hasAltmetrics:0,dateSubmitted:"April 3rd 2018",dateReviewed:"July 13th 2018",datePrePublished:"November 5th 2018",datePublished:null,dateFinished:null,readingETA:"0",abstract:"When studying ferrofluidic flows, as one example of magnetic flow dynamics, in terms of instability, bifurcation, and properties, one quickly finds out the additional challenges magnetic fluids introduce compared to the investigation of “classical”, “ordinary” shear flows without any kind of particles. Approximation of ferrofluids as fluids including point-size particles or, more realistic fine size particles, the relaxation times of the magnetic particle, their interaction between each other, i.e., the agglomeration and chain forming effects, and the interaction/response between any external applied field and the internal magnetization are just few examples of challenges to overcome. Further dependence on the considered model system, the direction of the external applied magnetic field (homogeneous or inhomogeneous) is crucial, as it can break the system symmetry and thus generate new solutions. As a result, the classical Navier–Stokes equations become modified to the more complex ferrohydrodynamical equation of motion, incorporating magnetic field and magnetization of the fluid itself, which typically makes numerical simulations expensive and challenging. This chapter provides an overview of the tasks/difficulties from describing and simulating magnetic particles, their interaction, and thus finally resulting modification in rotating flow structures and in particular instabilities and bifurcation behavior.",reviewType:"peer-reviewed",bibtexUrl:"/chapter/bibtex/63437",risUrl:"/chapter/ris/63437",book:{slug:"pattern-formation-and-stability-in-magnetic-colloids"},signatures:"Sebastian Altmeyer",authors:null,sections:[{id:"sec_1",title:"1. Introduction",level:"1"},{id:"sec_2",title:"2. System and theoretical description",level:"1"},{id:"sec_2_2",title:"2.1 Magnetization equations",level:"2"},{id:"sec_3_2",title:"2.2 Ferrohydrodynamic equations of motions",level:"2"},{id:"sec_5",title:"3. Ferrofluidic flow under finite magnetic field",level:"1"},{id:"sec_5_2",title:"3.1 Introduction and motivation",level:"2"},{id:"sec_6_2",title:"3.2 Flow structures",level:"2"},{id:"sec_6_3",title:"3.2.1 Flow structures in absence of a magnetic field",level:"3"},{id:"sec_7_3",title:"3.2.2 Flow structures—mode coupling in magnetic fields",level:"3"},{id:"sec_9_2",title:"3.3 Effect of finite magnetic fields",level:"2"},{id:"sec_9_3",title:"3.3.1 Basic state and symmetries",level:"3"},{id:"sec_10_3",title:"3.3.2 Bifurcation behavior",level:"3"},{id:"sec_11_3",title:"3.3.3 Vortex structures and mode contents",level:"3"},{id:"sec_12_3",title:"3.3.4 Stabilization of basic state: point of higher co-dimension",level:"3"},{id:"sec_13_3",title:"3.3.5 Marginal stability thresholds: Stabilization of the basic state",level:"3"},{id:"sec_15_2",title:"3.4 Resume",level:"2"},{id:"sec_17",title:"4. Agglomeration and elongational flow effects",level:"1"},{id:"sec_17_2",title:"4.1 Introduction and motivation",level:"2"},{id:"sec_18_2",title:"4.2 Modified internal magnetic field",level:"2"},{id:"sec_19_2",title:"4.3 Linear stability behavior and analysis",level:"2"},{id:"sec_19_3",title:"4.3.1 Marginal stability thresholds",level:"3"},{id:"sec_20_3",title:"4.3.2 Magnetic field dependence of thresholds",level:"3"},{id:"sec_21_3",title:"4.3.3 Frequencies of helical SPI states",level:"3"},{id:"sec_22_3",title:"4.3.4 Axial wavenumber dependence",level:"3"},{id:"sec_24_2",title:"4.4 Nonlinear analysis",level:"2"},{id:"sec_24_3",title:"4.4.1 Bifurcation thresholds",level:"3"},{id:"sec_25_3",title:"4.4.2 Flow structure modifications",level:"3"},{id:"sec_26_3",title:"4.4.3 Wavenumber modification and selection",level:"3"},{id:"sec_28_2",title:"4.5 Resume",level:"2"},{id:"sec_30",title:"5. Modified field dependence - internal magnetization",level:"1"},{id:"sec_30_2",title:"5.1 Introduction and motivation",level:"2"},{id:"sec_31_2",title:"5.2 Modified internal magnetic field",level:"2"},{id:"sec_32_2",title:"5.3 Flow structural modifications",level:"2"},{id:"sec_32_3",title:"5.3.1 Primary instabilities",level:"3"},{id:"sec_33_3",title:"5.3.2 Modified flow structures",level:"3"},{id:"sec_35_2",title:"5.4 Resume",level:"2"},{id:"sec_37",title:"6. Conclusions",level:"1"},{id:"sec_43",title:"Conflict of interest",level:"1"},{id:"sec_39",title:"A. Appendix",level:"1"},{id:"sec_38",title:"A.1 Numerical method",level:"1"},{id:"sec_39",title:"Acknowledgments/Funding",level:"1"},{id:"sec_42",title:"Nomenclature and abbreviations",level:"1"}],chapterReferences:[{id:"B1",body:'Rosensweig RE. Ferrohydrodynamics. Cambridge University Press, Cambridge; 1985. 344 pp. Journal of Fluid Mechanics, 200, 597-599 DOI: 10.1017/S0022112089220773'},{id:"B2",body:'Hart JE. Dynamics of Atmospheres and Oceans. 2006;41. DOI: 10.1016/j.dynatmoce.2006.03.001'},{id:"B3",body:'Hart JE, Kittelman S. Dynamics of Atmospheres and Oceans. 2006;41:139. DOI: 10.1017/S0022112075001309'},{id:"B4",body:'Niklas M. Zeitschrift fuer Physik B: Condensed Matter. 1987;68:493. DOI: 10.1007/BF01471080'},{id:"B5",body:'Niklas M, Müller-Krumbhaar H, Lücke M. Journal of Magnetism and Magnetic Materials. 1989;81:29. DOI: 10.1016/0304-8853(89)90225-4'},{id:"B6",body:'Altmeyer S, Hoffmann C, Leschhorn A, Lücke M. Physical Review E. 2010;82:016321. DOI: 10.1103/PhysRevE.82.016321'},{id:"B7",body:'Leschhorn A, Lücke M, Hoffmann C, Altmeyer S. Physical Review E. 2009;79:036308. DOI: 10.1103/PhysRevE.79.036308'},{id:"B8",body:'Reindl M, Leschhorn A, Lücke M, Odenbach S. Journal of Physics: Conference Series. 2009;149:012109. DOI: 10.1088/1742-6596/149/1/012109'},{id:"B9",body:'Reindl M, Odenbach S. Physics of Fluids. 2011;23:093102. DOI: 10.1063/1.3633341'},{id:"B10",body:'Odenbach S, Müller HW. Journal of Magnetism and Magnetic Materials. 2005;289:242. DOI: 10.1016/j.jmmm.2004.11.069'},{id:"B11",body:'Mahle S, Ilg P, Liu M. Physical Review E. 2008;77:016305. DOI: 10.1103/PhysRevE.77.016305'},{id:"B12",body:'Müller HW, Liu M. Physical Review E. 2001;64:061405. DOI: 10.1103/PhysRevE.64.061405'},{id:"B13",body:'Debye PJW. Polar Molecules. New York: Dover; 1929 https://hdl.handle.net/2027/uc1.$b45989'},{id:"B14",body:'McTague JP. The Journal of Chemical Physics. 1969;51:133. DOI: 10.1063/1.1671697'},{id:"B15",body:'Martsenyuk MA, Raikher YL, Shliomis MI. Soviet Physics - JETP. 1974;38:413 http://www.jetp.ac.ru/cgi-bin/dn/e_038_02_0413.pdf'},{id:"B16",body:'Gilles de GP, Prost J, editors. The Physics of Liquid Crystals. Oxford: Clarendon; 1983. ISBN: 9780198517856'},{id:"B17",body:'Leschhorn A, Embs JP, Lücke M. Journal of Physics: Condensed Matter. 2006;18:S2633. DOI: 10.1088/0953-8984/18/38/S07'},{id:"B18",body:'Embs JP, May S, Wagner C, Kityk AV, Leschhorn A, Lücke M. Physical Review E. 2006;73:036302. DOI: 10.1103/PhysRevE.73.036302'},{id:"B19",body:'Taylor GI. Philosophical Transactions of the Royal Society of London A. 1923;223:289. DOI: 10.1098/rsta.1923.0008'},{id:"B20",body:'Chossat P, Iooss G. The Couette-Taylor Problem. Berlin: Springer; 1994. DOI: 10.1007/978-1-4612-4300-7'},{id:"B21",body:'Langevin P. Annales de Chimie Physique. 1905;5(70)'},{id:"B22",body:'Embs J, Müller HW, Wagner C, Knorr K, Lücke M. Physical Review E. 2000;61:R2196. DOI: 10.1103/PhysRevE.61.R2196'},{id:"B23",body:'Odenbach S, Gilly H. Journal of Magnetism and Magnetic Materials. 1995;152(123). DOI: 10.1016/0304-8853(95)00442-4'},{id:"B24",body:'Altmeyer S, Lopez JM, Do Y. Physical Review E. 2012;85:066314. DOI: 10.1103/PhysRevE.85.066314'},{id:"B25",body:'Zubarev AY, Skakova Y. Physical Review E. 2000;61:5415. DOI: 10.1103/PhysRevE.61.5415'},{id:"B26",body:'Golubitsky M, Stewart I, Schaeffer D. Singularities and Groups in Bifurcation Theory II. New York: Springer; 1988. ISBN: 978-1-4612-4574-2'},{id:"B27",body:'Golubitsky M, Stewart I. SIAM Journal on Mathematical Analysis. 1986;17:249. DOI: 10.1137/0517023'},{id:"B28",body:'Ch H, Altmeyer S, Pinter A, Lücke M. New Journal of Physics. 2009;11:053002. DOI: 10.1088/1367-2630/11/5/053002'},{id:"B29",body:'Davey A, Di Prima RC, Stuart JT. Journal of Fluid Mechanics. 1968;31:17. DOI: 10.1017/S0022112068000029'},{id:"B30",body:'Iooss G. Journal of Fluid Mechanics. 1986;173:273. DOI: 10.1017/S0022112086001179'},{id:"B31",body:'Altmeyer S, Do Y, Marques F, Lopez JM. Physical Review E. 2012;86:046316. DOI: 10.1103/PhysRevE.86.046316'},{id:"B32",body:'Coles D. Journal of Fluid Mechanics. 1965;21:385. DOI: 10.1017/S0022112065000241'},{id:"B33",body:'Wereley ST, Lueptow RM. Journal of Fluid Mechanics. 1998;364:59. DOI: 10.1017/S0022112098008969'},{id:"B34",body:'Jones CA. Journal of Fluid Mechanics. 1981;102:249. DOI: 10.1017/S0022112081002620'},{id:"B35",body:'Akonur A, Lueptow RM. Physics of Fluids. 2003;15:947. DOI: 10.1063/1.1556615'},{id:"B36",body:'King GP, Li Y, Lee W, Swinney HL, Marcus PS. Journal of Fluid Mechanics. 1984;141:365. DOI: 10.1017/S0022112084000896'},{id:"B37",body:'Altmeyer S, Lopez JM, Do Y. Physical Review E. 2013;88:013003. DOI: 10.1103/PhysRevE.88.013003'},{id:"B38",body:'Odenbach S, Müller HW. Physical Review Letters. 2002;89:037202. DOI: 10.1103/PhysRevLett.89.037202'},{id:"B39",body:'Müller HW, Hahn D, Liu M. Journal of Physics: Condensed Matter. 2006;18:S2623. DOI: 10.1088/0953-8984/18/38/S06'},{id:"B40",body:'Chang MH, Chen CK, Weng HC. International Journal of Engineering Science. 2003;41:103. DOI: 10.1016/S0020-7225(02)00183-0'},{id:"B41",body:'Singh J, Bajaj R. Journal of Magnetism and Magnetic Materials. 2005;294:53. DOI: 10.1016/j.jmmm.2004.10.123'},{id:"B42",body:'Singh J, Bajaj R. International Journal of Mathematics and Mathematical Sciences. 2005;3727. DOI: 10.1155/IJMMS.2005.3727'},{id:"B43",body:'Altmeyer S, Leschhorn A, Ch H, Lücke M. Physical Review E. 2013;87:053010. DOI: 10.1103/PhysRevE.87.053010'},{id:"B44",body:'Reindl M, Odenbach S. Experiments in Fluids. 2010;50:375. DOI: 10.1007/s00348-010-0940-y'},{id:"B45",body:'Hirt CW, Nichols BD, Romero NC. NASA STI/Recon Technical Report N; 1975. http://adsabs.harvard.edu/abs/1975STIN...7532418H'}],footnotes:[],contributors:[{corresp:"yes",contributorFullName:"Sebastian Altmeyer",address:"sebastian.andreas.altmeyer@upc.edu",affiliation:'
Castelldefels School of Telecom and Aerospace Engineering, Universitat Politècnica de Catalunya, Barcelona, Spain
'}],corrections:null},book:{id:"7516",title:"Pattern Formation and Stability in Magnetic Colloids",subtitle:null,fullTitle:"Pattern Formation and Stability in Magnetic Colloids",slug:"pattern-formation-and-stability-in-magnetic-colloids",publishedDate:"January 8th 2020",bookSignature:"Nicola?s O. Rojas",coverURL:"https://cdn.intechopen.com/books/images_new/7516.jpg",licenceType:"CC BY 3.0",editedByType:"Edited by",editors:[{id:"208723",title:"Dr.",name:"Nicolás O.",middleName:null,surname:"Rojas",slug:"nicolas-o.-rojas",fullName:"Nicolás O. Rojas"}],productType:{id:"1",title:"Edited Volume",chapterContentType:"chapter",authoredCaption:"Edited by"},chapters:[{id:"65746",title:"Introductory Chapter: An Example in Superparamagnetic Colloids",slug:"introductory-chapter-an-example-in-superparamagnetic-colloids",totalDownloads:374,totalCrossrefCites:0,signatures:"Nicolás O. Rojas",authors:[{id:"208723",title:"Dr.",name:"Nicolás O.",middleName:null,surname:"Rojas",fullName:"Nicolás O. Rojas",slug:"nicolas-o.-rojas"}]},{id:"63437",title:"Interaction of Magnetic Fields on Ferrofluidic Taylor-Couette Flow",slug:"interaction-of-magnetic-fields-on-ferrofluidic-taylor-couette-flow",totalDownloads:396,totalCrossrefCites:0,signatures:"Sebastian Altmeyer",authors:[null]},{id:"63763",title:"Equilibrium Structures and Stationary Patterns on Magnetic Colloidal Fluids",slug:"equilibrium-structures-and-stationary-patterns-on-magnetic-colloidal-fluids",totalDownloads:394,totalCrossrefCites:0,signatures:"Ricardo Peredo Ortíz, Martin Hernández Contreras and Raquel Hernández Gómez",authors:[null]},{id:"70082",title:"The Study of Magneto-Convection Heat Transfer in a Partially Open Cavity Based on LBM",slug:"the-study-of-magneto-convection-heat-transfer-in-a-partially-open-cavity-based-on-lbm",totalDownloads:217,totalCrossrefCites:0,signatures:"Raoudha Chaabane",authors:[null]},{id:"63782",title:"Convection Flow of MHD Couple Stress Fluid in Vertical Microchannel with Entropy Generation",slug:"convection-flow-of-mhd-couple-stress-fluid-in-vertical-microchannel-with-entropy-generation",totalDownloads:477,totalCrossrefCites:3,signatures:"Abiodun A. Opanuga, Olasunmbo O. Agboola, Hilary I. Okagbue and Sheila A. Bishop",authors:[null]},{id:"63644",title:"Free Convection in a MHD Open Cavity with a Linearly Heated Wall Using LBM",slug:"free-convection-in-a-mhd-open-cavity-with-a-linearly-heated-wall-using-lbm",totalDownloads:285,totalCrossrefCites:0,signatures:"Raoudha Chaabane",authors:[null]}]},relatedBooks:[{type:"book",id:"6514",title:"Microfluidics and Nanofluidics",subtitle:null,isOpenForSubmission:!1,hash:"4ec06fd827f4dc0d3d7653eda88662de",slug:"microfluidics-and-nanofluidics",bookSignature:"Mohsen Sheikholeslami Kandelousi",coverURL:"https://cdn.intechopen.com/books/images_new/6514.jpg",editedByType:"Edited by",editors:[{id:"185811",title:"Dr.",name:"Mohsen",surname:"Sheikholeslami Kandelousi",slug:"mohsen-sheikholeslami-kandelousi",fullName:"Mohsen Sheikholeslami Kandelousi"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"},chapters:[{id:"61556",title:"Microfluidics and Nanofluidics: Science, Fabrication Technology (From Cleanrooms to 3D Printing) and Their Application to Chemical Analysis by Battery-Operated Microplasmas-On-Chips",slug:"microfluidics-and-nanofluidics-science-fabrication-technology-from-cleanrooms-to-3d-printing-and-the",signatures:"Vassili Karanassios",authors:[{id:"60925",title:"Prof.",name:"Vassili",middleName:null,surname:"Karanassios",fullName:"Vassili Karanassios",slug:"vassili-karanassios"}]},{id:"61632",title:"Application of Nanofluids for Thermal Management of Photovoltaic Modules: A Review",slug:"application-of-nanofluids-for-thermal-management-of-photovoltaic-modules-a-review",signatures:"Hafiz Muhammad Ali, Tayyab Raza Shah, Hamza Babar and\nZargham Ahmad Khan",authors:[{id:"187624",title:"Dr.",name:"Hafiz Muhammad",middleName:null,surname:"Ali",fullName:"Hafiz Muhammad Ali",slug:"hafiz-muhammad-ali"},{id:"229676",title:"Mr.",name:"Hamza",middleName:null,surname:"Babar",fullName:"Hamza Babar",slug:"hamza-babar"},{id:"241251",title:"Mr.",name:"Tayyab",middleName:"Raza",surname:"Shah",fullName:"Tayyab Shah",slug:"tayyab-shah"},{id:"241252",title:"Mr.",name:"Zargham Ahmad",middleName:null,surname:"Khan",fullName:"Zargham Ahmad Khan",slug:"zargham-ahmad-khan"}]},{id:"61387",title:"Thermodynamic Mechanism of Nanofluid Minimum Quantity Lubrication Cooling Grinding and Temperature Field Models",slug:"thermodynamic-mechanism-of-nanofluid-minimum-quantity-lubrication-cooling-grinding-and-temperature-f",signatures:"Min Yang, Changhe Li, Yanbin Zhang, Dongzhou Jia, Runze Li and\nWenfeng Ding",authors:[{id:"13294",title:"Prof.",name:"Changhe",middleName:null,surname:"Li",fullName:"Changhe Li",slug:"changhe-li"}]},{id:"61051",title:"Biological Particle Control and Separation using Active Forces in Microfluidic Environments",slug:"biological-particle-control-and-separation-using-active-forces-in-microfluidic-environments",signatures:"Mohd Anuar Md Ali, Aminuddin Bin Ahmad Kayani and\nBurhanuddin Yeop Majlis",authors:[{id:"194385",title:"Prof.",name:"Burhanuddin Yeop",middleName:null,surname:"Majlis",fullName:"Burhanuddin Yeop Majlis",slug:"burhanuddin-yeop-majlis"},{id:"223436",title:"Dr.",name:"Aminuddin Bin",middleName:null,surname:"Ahmad Kayani",fullName:"Aminuddin Bin Ahmad Kayani",slug:"aminuddin-bin-ahmad-kayani"},{id:"223544",title:"Mr.",name:"Mohd Anuar",middleName:null,surname:"Md Ali",fullName:"Mohd Anuar Md Ali",slug:"mohd-anuar-md-ali"}]},{id:"62513",title:"Particle Deposition in Microfluidic Devices at Elevated Temperatures",slug:"particle-deposition-in-microfluidic-devices-at-elevated-temperatures",signatures:"Zhibin Yan, Xiaoyang Huang, Chun Yang, Mingliang Jin and\nLingling Shui",authors:[{id:"228610",title:"Dr.",name:"Zhibin",middleName:null,surname:"Yan",fullName:"Zhibin Yan",slug:"zhibin-yan"},{id:"248748",title:"Prof.",name:"Xiaoyang",middleName:null,surname:"Huang",fullName:"Xiaoyang Huang",slug:"xiaoyang-huang"},{id:"248749",title:"Prof.",name:"Chun",middleName:null,surname:"Yang",fullName:"Chun Yang",slug:"chun-yang"},{id:"248761",title:"Prof.",name:"Mingliang",middleName:null,surname:"Jin",fullName:"Mingliang Jin",slug:"mingliang-jin"},{id:"248764",title:"Prof.",name:"Lingling",middleName:null,surname:"Shui",fullName:"Lingling Shui",slug:"lingling-shui"}]},{id:"59995",title:"Microdroplets Advancement in Newtonian and Non- Newtonian Microfluidic Multiphase System",slug:"microdroplets-advancement-in-newtonian-and-non-newtonian-microfluidic-multiphase-system",signatures:"Kai Seng Koh, Voon Loong Wong and Yong Ren",authors:[{id:"177059",title:"Dr.",name:"Yong",middleName:null,surname:"Ren",fullName:"Yong Ren",slug:"yong-ren"},{id:"222878",title:"Dr.",name:"Kai Seng",middleName:null,surname:"Koh",fullName:"Kai Seng Koh",slug:"kai-seng-koh"},{id:"224066",title:"Dr.",name:"Voon Loong",middleName:null,surname:"Wong",fullName:"Voon Loong Wong",slug:"voon-loong-wong"}]},{id:"59924",title:"Micro/Nanofluids in Sustainable Machining",slug:"micro-nanofluids-in-sustainable-machining",signatures:"Tran The Long and Tran Minh Duc",authors:[{id:"222893",title:"Dr.",name:"TranThe",middleName:null,surname:"Long",fullName:"TranThe Long",slug:"tranthe-long"},{id:"227561",title:"Prof.",name:"Tran Minh",middleName:null,surname:"Duc",fullName:"Tran Minh Duc",slug:"tran-minh-duc"}]},{id:"59808",title:"Precisely Addressed (DNA Gene) Spray Microfluidic Chip Technology",slug:"precisely-addressed-dna-gene-spray-microfluidic-chip-technology",signatures:"Jian-Chiun Liou",authors:[{id:"222901",title:"Prof.",name:"Jian-Chiun",middleName:null,surname:"Liou",fullName:"Jian-Chiun Liou",slug:"jian-chiun-liou"}]},{id:"59009",title:"Thermal Transport and Challenges on Nanofluids Performance",slug:"thermal-transport-and-challenges-on-nanofluids-performance",signatures:"José Jaime Taha-Tijerina",authors:[{id:"182402",title:"Dr.",name:"Jose",middleName:"Jaime",surname:"Taha-Tijerina",fullName:"Jose Taha-Tijerina",slug:"jose-taha-tijerina"}]},{id:"59126",title:"Magnetite Molybdenum Disulphide Nanofluid of Grade Two: A Generalized Model with Caputo-Fabrizio Derivative",slug:"magnetite-molybdenum-disulphide-nanofluid-of-grade-two-a-generalized-model-with-caputo-fabrizio-deri",signatures:"Farhad Ali, Madeha Gohar, Ilyas Khan, Nadeem Ahmad Sheikh,\nSyed Aftab Alam Jan and Muhammad Saqib",authors:[{id:"190576",title:"Dr.",name:"Ilyas",middleName:null,surname:"Khan",fullName:"Ilyas Khan",slug:"ilyas-khan"},{id:"225942",title:"Mr.",name:"Nadeem Ahmad",middleName:null,surname:"Sheikh",fullName:"Nadeem Ahmad Sheikh",slug:"nadeem-ahmad-sheikh"},{id:"225943",title:"Mr.",name:"Muhammad",middleName:null,surname:"Saqib",fullName:"Muhammad Saqib",slug:"muhammad-saqib"},{id:"225983",title:"Dr.",name:"Farhad",middleName:null,surname:"Ali",fullName:"Farhad Ali",slug:"farhad-ali"},{id:"225985",title:"Ms.",name:"Madeha",middleName:null,surname:"Gohar",fullName:"Madeha Gohar",slug:"madeha-gohar"},{id:"225987",title:"Mr.",name:"Syed Aftab Alam",middleName:null,surname:"Jan",fullName:"Syed Aftab Alam Jan",slug:"syed-aftab-alam-jan"}]},{id:"60601",title:"Performance Evaluation Criterion of Nanofluid",slug:"performance-evaluation-criterion-of-nanofluid",signatures:"Sudarmadji Sudarmadji, Bambang SAP and Santoso",authors:[{id:"226665",title:"Dr.",name:"Sudarmadji",middleName:null,surname:"Sudarmadji",fullName:"Sudarmadji Sudarmadji",slug:"sudarmadji-sudarmadji"},{id:"239852",title:"Dr.",name:"Bambang",middleName:null,surname:"Sap",fullName:"Bambang Sap",slug:"bambang-sap"},{id:"239853",title:"MSc.",name:"Santoso",middleName:null,surname:"Santoso",fullName:"Santoso Santoso",slug:"santoso-santoso"}]},{id:"59675",title:"Microfluidic Membrane Filtration Systems to Study Biofouling",slug:"microfluidic-membrane-filtration-systems-to-study-biofouling",signatures:"Ishita Biswas, Aloke Kumar and Mohtada Sadrzadeh",authors:[{id:"215556",title:"Prof.",name:"Mohtada",middleName:null,surname:"Sadrzadeh",fullName:"Mohtada Sadrzadeh",slug:"mohtada-sadrzadeh"},{id:"227070",title:"Ms.",name:"Ishita",middleName:null,surname:"Biswas",fullName:"Ishita Biswas",slug:"ishita-biswas"},{id:"243533",title:"Prof.",name:"Aloke",middleName:null,surname:"Kumar",fullName:"Aloke Kumar",slug:"aloke-kumar"}]}]}]},onlineFirst:{chapter:{type:"chapter",id:"63437",title:"Interaction of Magnetic Fields on Ferrofluidic Taylor-Couette Flow",doi:"10.5772/intechopen.80301",slug:"interaction-of-magnetic-fields-on-ferrofluidic-taylor-couette-flow",body:'
1. Introduction
Ferrofluids [1] represent one example of magnetic particles with a variety of applications, ranging from liquid seals in rotating systems, their use in computer hard drives, as well as in laboratory experiments, to study geophysical flows [2, 3]. However, deeper fundamental study of their magnetohydrodynamics is inevitable, as there are many modeling assumptions that have been implemented to make theoretical descriptions tractable. Ferrofluids are manufactured fluids which consist of dispersions of magnetized nanoparticles in various kinds of liquid carriers. Stabilizing against agglomeration is typically achieved by the addition of a surfactant monolayer onto the magnetic particles. Being no magnetic field present, the magnetic nanoparticles are randomly orientated, which results in a zero-net magnetization of the fluid. Only effect of the present nanoparticles is a typically small alteration to the fluid’s density and viscosity. Applying a sufficiently strong magnetic field, the ferrofluid flows toward certain regions of the magnetic field. The hydrodynamics of the system can be significantly changed, and properties of the fluid, e.g., the viscosity, are altered. Therefore, nowadays, they are widely used by various car manufactures in order to provide adaptable suspensions, i.e., changing the damping characteristics of an antishock/damper due to modified viscosity of the ferrofluid under influence of an external applied magnetic field. Both numerical simulations and experiments show that any applied magnetic field, whether radial, axial, azimuthal, transversal, or any combination of the previous stabilizes the basic state [4, 5, 6, 7, 8, 9]. This stabilization effect depends, among others, on the particle-particle interactions of the fluid [9]. A strong interaction leads to large/crucial changes in flow behavior, which are driven by the formation of internal flow structures. While an axisymmetric field leaves the flow structures untouched (compared to classical Couette flow), a finite transverse, in-plane magnetic-field component breaks axisymmetry, enforcing new nonlinear effects, modify present flow structures, and can produce new flow structures.
Among various others, one important point to look at is the applied magnetic field, in particular its characteristics while penetrating a moving ferrofluid. In order to describe the hydrodynamics of ferrofluids in containers/boxes typically, models are used which assume that the internal magnetic field and the external applied magnetic field are equal [4, 5, 6]. However, such a simplification is only a leading order approximation, as the magnetic field in the container typically is modified, based on the magnetic susceptibility of the ferrofluid. First work trying to account for modifications [15] includes an expression for the angular velocity of the ferrofluid particles, which differs from the ferrofluid itself due to the magnetic torques acting on the particles. The results are significant differences from the external field, providing among others, a much changed body force in the governing equations. More concretely, the interaction between a uniform external field and the susceptibility of the ferrofluid-filled annulus can be shown to result in a field with an 1/r2 radial dependence superimposed on the external field.
Further common assumption when describing the hydrodynamics of magnetic particles as ferrofluids is to consider them as point-wise objects and ignore their finite size and specific properties. However, in real, such particles typically aggregate to form clusters having the form of chains, and thus they hinder the free flow of the fluid and increase the viscosity [9, 10, 11], and change their rheology and transport properties [14]. In this type of structure’s formation, it is also assumed that the interaction parameter is usually greater than unity [1], and thus the strength of the grain-grain interaction can be measured in terms of the total momentum of a particle. Regarding the magneto-dissipative structure of ferrofluid dynamics, one can derive, based on general principles, an equation for macroscopic ferrofluid dynamics. The magnetization’s relaxation equation includes an additional term that is proportional to the product of the magnetization’s magnitude and the symmetric part of the velocity gradient tensor, which can describe an elongational flow scaled by a so-called transport coefficientλ2. Such description also exists in the dynamics of nematic liquid crystals, as the flow alignment modifies the director field in an applied shear flow [16].
It is worth to mention that there are further crucial and important considerations for describing ferrofluids and in general fluids with magnetic particles which are also important as the discussed ones. However, these are not considered in the present chapter. Therefore, we refer to the huge amount of literature, which can be essentially found by the presented references. Two prominent examples of these further considerations are the modulation of the relaxation times of magnetic particles [7, 13], e.g., Neel’s or Braun’s and the consideration of monodisperse vs. polydisperse materials [7, 17, 18].
In this chapter, we will consider a rotating ferrofluid in Taylor-Couette geometry. First, we are going to describe the modifications in flow solutions and its structural properties and alterations in primary supercritical bifurcation thresholds induced/forced by an external applied magnetic field. Following, we will discuss two different approaches to come to a more realistic description of ferrofluidic flows. Therefore, we will modify and extend the basic ferrohydrodynamical equations. First, considering a modification of the internal magnetic field, which differs from the external applied one due to its internal magnetization. Second incorporating elongational flow. The latter aims to consider finite size effects of the magnetic particles, as well as agglomeration and chain formation procedures.
2. System and theoretical description
As a prototypical system to investigate the influences of a magnetic field on a rotating ferrofluid system, we consider Taylor-Couette flow (Taylor Couette system, TCS) [19, 20] driven by the independent differential rotation of two concentric cylinders (Figure 1).
Figure 1.
Schematic of the Taylor-Couette system with an external applied homogeneous magnetic field H→ext=Hxe→x+Hze→z.
The inner cylinder of radius r1 rotates at ω1, and the outer cylinder of radius r2 rotates at ω2. If not other stated, in this chapter, periodic boundary conditions are considered in the axial direction with periodicity of length λ and no-slip boundary conditions on the cylinder walls. The system is described using a cylindrical polar coordinate system rθz with a velocity field u→=uvw. The radius ratio of the cylinders is β=r1/r2, and the aspect ratio is Γ=λ/r2−r1. The gap between the cylinders is filled with a viscous, incompressible, and isothermal ferrofluid. An external homogeneous magnetic field H→ext=Hxe→x+Hze→z is applied, where x=rcosθ is the transverse direction and z is the axial direction, respectively (Figure 1).
The flow dynamics of an incompressible homogeneous monodisperse ferrofluid with kinematic viscosity ν and density ρ is governed by the incompressible Navier–Stokes equations, including magnetic terms and the continuity equation. Following, an approach based on the model of Niklas [4, 5] is used, where the magnetic fluid is assumed to be incompressible, nonconducting, and to have a constant temperature and a homogeneous distribution of magnetic particles. Using the gap d=r2−r1 as the length scale, the diffusion time τd=d2/ν as the time scale, scaling pressure with ρν2/d2, the magnetic field H→, and the magnetization M→ with ρ/μ00.5ν/d (μ0 is the magnetic constant, i.e., magnetic permeability of free space), the nondimensional governing equations are
∂tu→⋅∇u→−∇2u→∇p=M→⋅∇H→+∇×M→×H→/2,∇⋅u→=0E1
The cylinders are considered to be of no-slip type with velocity boundary conditions
u→r1,θz=0Re1,0andu→r2,θz=0Re20,E2
where the inner and outer Reynolds numbers are Re1=ωr1d/ν and Re2=ωr2d/ν, respectively.
2.1 Magnetization equations
In order to solve (1), one needs additional information about the magnetic particles themselves, in particular their properties and interaction with an applied magnetic field. Therefore, (1) is solved together with an equation that describes the magnetization of the ferrofluid. A first approximation is to use the equilibrium magnetization of an unperturbed state with a homogeneously magnetized ferrofluid at rest with the mean magnetic moments orientated in the direction of the magnetic field, M→eq=χH→, where χ is the magnetic susceptibility of the ferrofluid, determined using Langevin’s formula [21]. The ferrofluid considered for all simulations in this chapter is APG933 [22] with χ=0.9. However, a ferrofluid’s magnetization is also influenced by the flow field. For the presented numerical simulations, an approach based on the model of Niklas [4, 5, 6] is used. Furthermore, a stationary magnetization near equilibrium considering a small perturbation at the equilibrium magnetization M→−M→eq and small relaxation times Ωτ<<1, where Ω is the absolute value of half of the vorticity Ω and τ is the magnetic relaxation time. In the near-equilibrium approximation, Niklas [4] determined the relationship between the magnetization M→, the magnetic field H→, and the velocity u→ to be
M→−Meq=cNΩ→×H→E3
with Ω→=∇×u→/2 being the vorticity and the Niklas coefficient
cN2=τ/1/χ+τμ0H2/6μΦE4
where μ is the dynamic viscosity and Φ is the volume fraction of the magnetic material.
2.2 Ferrohydrodynamic equations of motions
The magnetization can be eliminated from (1) by using (3) in order to obtain the ferrohydrodynamic equation of motion [4]:
∂t+u→⋅∇u→−∇2u→+∇pM=−cN2H→∇⋅F→+H→×∇×F→/2E5
where F→=Ω→×H→ and pM are the dynamic pressure incorporating all magnetic terms which can be written as gradients. Note that the relaxation time τ could typically be defined by τ=τBD31/3, where τB is the Brownian relaxation time and D31/3 is the averaged diameter of the magnetic particles. For simplicity, to investigate the different effects, i.e., internal magnetization and elongational flow (see below), we assume that τ is constant, which is independent of the magnetic field H→. The considered relaxation time is based on the experimental result [22] τ=τAPG933/τd=0.0018.
Consider external field to be equal to internal one H→=H→ext, (5) can be simplified to
∂t+u→⋅∇u→−∇2u→+∇pM=−cN2H→∇⋅F→+H→×∇×F→/2E6
In this approach, the magnetic field and all the magnetic properties of the ferrofluid influence the velocity field only via the magnetic field, the Niklas parameters→N=sxe→x+sze→z [4, 5, 6], whereby magnetic field strengths are:
sx=22+χ/2+χ2−χ2η2HxcNandsz=HzcnE7
In this approach, sx and sz are dimensionless numbers, describing the strength of the applied magnetic field.
3. Ferrofluidic flow under finite magnetic field
3.1 Introduction and motivation
We start our discussion with the investigation of the influence of a homogeneous magnetic field on a ferrofluid. In this first part, we will assume that the internal magnetic field is equal to the external applied magnetic field, i.e., homogeneous. Moreover, we ignore in this section any kind of finite size effects, e.g., possible agglomeration and chain formation, later considered as elongational flow effect (Section 4). Therefore, we consider the ferrofluid particles as point-size objects without any 3D expansion. Structure, dynamics, symmetry properties, bifurcation, and stability behavior of different vortex structures are investigated for axial and transversal magnetic fields, as well as combinations of them, i.e., oblique magnetic field. We will see that a transversal magnetic field has significant influence. It modulates and changes the classical flow structures, as well the basic state (Circular Couette flow, CCF) as the classical primary supercritical ones, Taylor vortex flow (TVF) and the spiral vortex flow (spirals, SPI). They are replaced by wavy-modulated flow structures: wavy Taylor vortex flow (wTVF) and wavy spiral vortex flow (wSPI), as well as they forward bifurcate out of the modified 2-fold circular Couette flow (2-CCF).
3.2 Flow structures
3.2.1 Flow structures in absence of a magnetic field
There are two fundamental basic, topologically different flow structures appearing in (classical) Taylor-Couette flow (see also Nomenclature and abbreviations section). These are either toroidally closed Taylor vortices (Taylor vortex flow, TVF) and helical spiral vortices (spiral vortex flow/spirals, SPI) [20]. Figure 2 shows the structure of these fundamental solutions, TVF and SPI (here for a left-handed SPI with azimuthal wavenumber m=1), by using isosurfaces of the azimuthal vorticity. Because of the system symmetry z→−z, we only discuss left-handed SPI. Without other parameters, e.g., axial through flow, left- and right-handed SPI are mirror images of each other. These pure states can also be found modulated as wavy structures, which are topologically identical with the pure states. However, they show a wavy-like deformation that results from additional and other stimulated modes compared to the classical solutions.
Figure 2.
Azimuthal vorticity isosurfaces η over two axial wavelengths (for visual purpose) in absence of a magnetic field (N) of (a) TVF at Re2=0, Re1=100, and (b) SPI at Re2=−150, Re1=150. Red (yellow) isosurfaces correspond to positive (negative) values as indicated. The bottom row shows mode amplitudes um,n of the radial velocity field u corresponding to structures above over the m-n plane. The values are scaled regarding the maximum mode amplitude to be 1 (TVF: u0,1=1).
While classical wavy flow solutions [27, 28, 29, 30, 31, 32, 33, 34, 35, 36] are well known and have already been widely discussed in literature regarding TCS, the ones generated by a magnetic field are different. In particular, the magnetically induced wTVF are nonrotating structures having a pinned phase, contrary to the classic wTVF, which rotate azimuthally. Regarding the flow structure, the imposed magnetic field shrinks and expands the vortices at different cross sections of constant θ, resulting in an azimuthal alternating wide and narrow vortex belly structure.
3.2.2 Flow structures—mode coupling in magnetic fields
Following, we will have a quick look at the impact on the different applied magnetic field with respect to the flow structures and in particular their mode spectra. The modifications in NSE due to the additional terms incorporating magnetization and magnetic field result in various new phenomena. Most important is probably the fact that magnetic fields with a finite transversal component stimulate additional modes. Depending on different external applied magnetic fields, Figure 3 presents a schematic plot for the different stimulated modes, depending on the applied field.
Figure 3.
Schematic sketch of modes in the presence of different externally imposed (homogeneous) magnetic fields. Here, m represents the azimuthal, and n the axial mode index in the Fourier expansion (Appendix (17) and (18)). (a) Modes either in absence of any magnetic field or in a pure axial magnetic field. (b) Pure transversal magnetic field with additional stimulation of m=n±2 modes. (c) Superposition of axial and transversal fields with even further finite m=n±1 modes. Strength of mode amplitudes, the magnitudes are characterized from dark (large amplitudes) to bright (small amplitudes) colors (legend right). Red and blue lines indicate the dominant mode amplitudes of the flow structures.
Having a pure axial field sx=0sz>0 applied the mode spectra, the flow structures remain qualitatively similar to the no field (N) case (Figure 3(a)). A TVF still contains only the axisymmetric modes m=0 (with dominant n = 1, see also Figures 2 and 6). Similarly, L1-SPI, i.e., a left-handed SPI flow with azimuthal wavenumber m=1 contains modes only on the diagonal m=n (with dominant n=1, see also Figures 2 and 6). However, if a field with finite transversal component is applied additional mode become finite and the flow structures change. A pure transversal field sx>0sz=0 excites m=±2 modes lying on the secondary diagonal m=n±2 (Figure 3(b)). Thus, as a result, the pure states of TVF and SPI vanish in the presence of a transversal magnetic field. They become wavy-like modified, resulting in wTVF and wSPI. Finally, having a magnetic field with either axial and transversal component sx>0sz>0, there are further interactions, which result in stimulation of additional m=n±1 modes (Figure 3(c)). Moreover, these additional modes can induce further nonlinear mode couplings. As a result, the mode spectra contain more nonzero components (Figure 3(c) and Figure 6 for numerical results). Thus all flow structures in a magnetic field with transversal component differ qualitatively from the classical flows found in the absence of magnetic fields [6, 26].
3.3 Effect of finite magnetic fields
3.3.1 Basic state and symmetries
If only a pure axial magnetic field, Hextre→z (see Figure 1), the classical circular Couette flow (CCF) is observed, with u→CCF=0Ar+B/r0, where A=Re2−βRe1/1+β, and B=βRe1−βRe2/1+β1−β2, is a solution of (5) as the associated vorticity Ω→ccf=∇×u→ccf is parallel to the applied magnetic field and therefore all magnetic terms vanish. Although for magnetic fields that are orientated purely in the radial or azimuthal direction, the basic state changes but remains axisymmetric with deviations from CCF only having an azimuthal component [4, 23]. As a result, in all of these cases, the basic state is invariant to a number of symmetries (azimuthal rotation, mirror symmetry, and axial translation) whose actions on a general velocity field are
RΦuvwrθzt=uvwrθ+Φzt,Φ∈02π,E8
Kzuvwrθzt=uv−wrθ−zt,E9
Tαuvwrθzt=uvwrθz+αΓt,α∈NE10
On the contrary, for a field with a finite transverse component, the base state is no longer axisymmetric. Although the basic state remains invariant to Tα and a modified mirror symmetry KzH [24], a transverse magnetic field T breaks the continuous axisymmetry RΦ, resulting in a basic state with discrete symmetry Rπ, i.e., with azimuthal wavenumber m=2 [24]. Following, we will refer to this modified 2-fold symmetric CCF as 2-CCF. It is important to mention that this m=2 state is stationary [24], which is in contrast to the generic breaking of a SO2 symmetry which would result in a rotating wave.
3.3.2 Bifurcation behavior
As seen in Section 3.2, an applied magnetic field can stimulate additional modes and thus changes the Fourier mode spectrum of flow structures resulting in new solutions. The results are wavy-like modified flow states. It is to expect that these mode-stimulated modifications also influence and appear in the primary bifurcation thresholds of the solutions, as well as in their growth rate. Following, we will have a look at the bifurcation behavior of the different flow structures with the additionally field-induced modes explained in Section 3.3. In particular, we will focus on four different field configurations: no field (N)sx=0sz=0; pure axial field (A)sx=0sz=0.6; pure transversal field (T)sx=0.6sz=0; and oblique field (O)sx=0.6sz=0.6 ((7) for parameters sx,sz).
Considering these four, non, pure axial, pure transversal, and oblique magnetic fields, Figure 4 illustrates the stable forward bifurcating branches of toroidally closed TVF and of wTVF solutions, respectively. Figure 5 shows the bifurcation diagrams for helical SPI and wSPI in an analogous way. Note that the shown wavy-like modulated flow states only result from symmetry breaking due to finite transversal field component. In the Figures 4 and 5, the relative distance μ=Re1/Re1,c−1 is used as control parameter for easier comparison of the bifurcation branches (Re1,c stands for the critical onset Reynolds number Re1 of the flow structure, which depends on various system parameters. These are the aspect ratio Γ, the radius ratio r1/r2, the outer Reynolds number Re2, axial boundary conditions, and the field strength sx, sz, respectively). The main reason to use μ is the shift in the bifurcation onset, which varies depending on the applied field direction.
Figure 4.
Bifurcation diagrams of TVF and wTVF, respectively, in different magnetic fields: (a) no field (N), (b) axial (A), (c) transverse (T), and (d) oblique (O). Moduli um,n (18) of the radial flow field amplitudes at mid gap are shown versus the relative distance μ=Re1/Re1,c−1 from the onset of the respective vortex flow. Further parameters are: Re2=0,β=0.5,k=3.927 (cf. Figure 2).
Figure 5.
Bifurcation diagrams of SPI and wSPI, respectively, in different magnetic fields: (a) no field (N), (b) axial (A), (c) transverse (T), and (d) oblique (O). Moduli um,n (18) of the radial flow field amplitudes at mid gap are shown versus the relative distance μ=Re1/Re1,onset−1 from the onset of the respective vortex flow. Further parameters are: Re2=−150,β=0.5,k=3.927 (cf. Figure 2).
Displayed in the Figures 4 and 5 are the moduli um,n of the dominant, as well as the first higher harmonic modes of the corresponding structure (see Appendix (18) and (19)), e.g. mn=01 and mn=02 for TVF. Furthermore, the largest field-induced mode amplitudes u2,1 and u1,±1 (only in oblique configuration) of the wavy structures (cf. schematics in Figure 3) are included for comparison.
Axial field (A): As a pure axial magnetic field sx=0sz≠0 does not stimulate any further modes within the mode-spectra (Figure 3(a)), the flow states TVF and SPI and their structural properties remain qualitatively similar, and they are as without a field. However, quantitative, they eventually change. First important effect is the rise of the onsets to higher values of Re1 [6, 7]. Thereby, the upward shift depends on the various system control parameters (e.g., Re2,β,Γ). For parameters in Figure 4, the onset of TVF rises by about 25%, and for the parameters of Figure 5, the SPI bifurcation threshold moves upward in by about 15%. For the latter, consequently also corresponding SPI frequencies are increased in an axial field. Other change regarding Figure 4 is the amplitude of the modes with increasing the parameter μ (weakest for nonfield and strongest for oblique field configuration).
Transversal field (T): A transversal magnetic field sx≠0sz=0 stimulates additional azimuthal modes m=n±2 (Section 3.2) as illustrated in Figure 3(b). Thus, based on this knowledge, one can already expect a more dramatic influence on the classical TVF and SPI structures. As a pure axial field, a transversal field also delays the onset of the primary bifurcating structures, i.e., it stabilizes the basic state. Furthermore and more important is the fact that such a field also changes the flow structural properties in qualitative manner. The growth of the additional modulus u2,1 (18) (for toroidally closed flow states) with μ=Re1/Re1,onset−1 changes the pure TVF to a wavy-like modulated wTVF. Analog pure SPI changes to wSPI (shown for L1-(w)SPI in Figure 5) due to the additional m=n±2 (see Figure 3(b)), indicated by the growth of the modulus u3,1 in Figure 5. It is worth to mention that due to the mode stimulation and interaction, also the (1,-1) mode becomes excited, which is similar to structure R1-(w)SPI. Note, the same also holds for the symmetry-related mirror image L1-(w)SPI.
Thus, in the presence of a magnetic field with a finite transversal component, the pure/classical TVF and SPI structures are modified. They are replaced by their wavy-like counterparts, wTVF and wSPI, which now bifurcate as primary structured solutions forward out of the CCF-2 ground state; note that already the basic state has an additional mode-2 symmetry due to the symmetry-breaking effect of the transversal field and as such already differ from the classical CCF (Section 3.3.1). As for axial field, the magnitude of stabilization varies with given system control parameters. In general, it is slightly smaller compared to the stabilization effect of an axial field. For instance for parameters presented in Figures 4 and 5, the upward shift for the onset of wTVF (wSPI) is about 20% (12%) in an axial (transversal) magnetic field.
Important to mention that a finite transverse component sx≠0 breaks the system symmetry but does not change the degeneration between left- and right-handed SPIs. Both solutions keep bifurcating at the same critical values. The latter is different to the effect of an axial flow, for instance. Aside that such a flow also breaks the system symmetry, it separates both left- and right-handed SPI solutions. Under axial applied flow, they appear at different critical values.
Oblique field (O): As already seen for both pure field configurations, either sx≠0,sz=0 or sx=0,sz≠0, the onset of (w)TVF and (w)SPI is shifted to higher values. In fact in a field, obliquely orientated to the cylinder sx≠0sz≠0, the magnitude of shift is bigger compared to the pure cases. In addition, as discussed before, such a superposition of an axial and a transversal field stimulates also new, additional modes with higher mode index combinations m=n±1 (Figure 3(c) and for numerical simulations Figures 4(c,d), 5(c,d), and 6(e,f)).
Figure 6.
Azimuthal vorticity isosurfaces η over two axial wavelengths (for visual purpose) with pure transversal applied magnetic field (T): (a) wTVF at Re2=0, Re1=125, (b) wSPI at Re2=−150, Re1=160; and oblique applied magnetic field (O): (c) wTVF at Re2=0, Re1=150, and (d) wSPI at Re2=−150, Re1=190. Red (yellow) isosurfaces correspond to positive (negative) values as indicated. Bottom row shows mode amplitudes um,n of the radial velocity field u corresponding to structures above over the m-n-plane. Values are scaled with the maximum mode amplitude to be 1 (TVF: u0,1=1).
A common finding for (w)TVF and (w)SPI is that their dominant mode amplitudes for m=0 (for (w)TVF) and m=1 (for (w)SPI) grow with respect to the bifurcation onsets. Although the absolute value depends on the different applied magnetic field, comparing the moduli um,1 at the same relative distance μ from the respective onset (Figures 4 and 5, compare dominant mode amplitudes), one finds the following relations
for a given constant value c∈01 and m=0 in case of toroidally closed (w)TVF and m=1 for helical (w)SPI, respectively (Compare upper curves for dominant mode amplitudes um,1 in Figures 4 and 5).
For helical flow structures, the larger variety of mode interactions due to the basic spectra (e.g. m=n [m=−n] for pure L1-SPI [R-SPI]) results in significant variation of the higher order and secondary stimulated modes, depending on the field configuration (see in particular Figures 3(c), 4(d), and 5(d)).
Based on the both bifurcation diagrams presented in Figures 4 and 5, it is worth to remind that in the case of a finite transversal field component sx≠0, no pure solutions, TVF and SPI, exist. They become replaced by wavy-like modulated wTVF and wSPI. However, independent of the field, leading mode amplitude of the respective flow state, (w)TVF and (w)SPI, follows a typical square-root forward bifurcation, whereby the slope of its square grows with the applied field strength.
3.3.3 Vortex structures and mode contents
Having discussed the changes that different oriented magnetic fields can cause in the bifurcation behavior and mode stimulation of (w)TVF and (w)SPI, respectively, one can imagine that these also coincide with flow structural modifications. Now, we will address these structural changes of the vortices and the associated changes in the mode contents of the flows in more detail.
Figure 6 shows isosurfaces of the azimuthal vorticity η over two axial wavelengths (for visual purpose) for supercritical flow states (w)TVF and (w)SPI in the three different considered field configurations: (a,b) no field applied (N), (c,d) pure transverse field (T), and (e,f) oblique field (O) at parameters Re1 and Re2 as indicated.
A good quantity in order to visualize the field induced changes in the 3D vortex structures are isosurfaces of the azimuthal vorticity η. These surfaces appropriately convey structural details of the vortex flows in question [28].
Axial field (A): As a pure axial field does not break any continuous symmetry of the system, there is neither qualitative change in the structure of TVF and SPI in real space nor the mode structure of the flow in the m−n Fourier plane (not shown, but similar as Figure 2(a,b); only quantitative the mode amplitudes are slightly larger). TVF still consists of toroidally closed, rotationally symmetric vortex tubes with a mode spectrum containing only m=0 modes (Figure 2(a)). And also, the helically oriented, open vortex tubes of SPI (here for left-winding spiral L-SPI) containing m=n modes as shown in Figures 2(b) and 3 (mirror image of right-winding R-SPI analog) are unaffected by an axial field.
Transversal and oblique field (T, O): Having a magnetic field with finite transversal component crucially changes the scenario as TVF, and SPI structures are replaced by modified flow states due to the presence of such a magnetic field. In general, the isosurfaces become wavy-like deformed, and additional Fourier modes are stimulated as described before (Section 3.3; see also schematics in Figure 3). For pure transversal field, this mode stimulation results in wTVF with additional m=±2 modes aside the m=0 (Figure 6(c)) and in wSPI with m=n±2 modes on the secondary diagonal(s) in addition to the m=n SPI modes on the diagonal (Figure 6(d)). Beyond these, an oblique field further stimulates the modes m=±1 in wTVF and m=n±1 in wSPI due to more complex nonlinear dynamics/mode interaction. In particular, the spatio-temporal properties of the wTVF generated due to symmetry breaking magnetic field, either transversal or oblique (Figure 6(c,e)), differ significant from those of the classic wTVF. The structure of classical wTVF rotates as a whole in azimuthal direction [27, 32, 33, 34, 35].
The wavy-like deformation on the toroidally closed vortices η introduced by a magnetic field with a finite transversal component is different. Regarding Figure 6(c,e), one sees that the tubes of isosurfaces of the azimuthal vorticity become strong modified and bent at certain azimuthal cross sections of constant θ in comparison with the pure TVF tubes. Classic wTVF arises by deforming the TVF such that the axial stack of closed vortex tubes is wave-like modulated, bending axially upward and downward. Then, this deformation structure of the vortices rotates as a whole with a characteristic frequency. In contrast to these classical rotating wTVF, the wTVF initiated due to the symmetry-breaking transversal field component is stationary and does not rotate in azimuthal direction, instead, developing an 2-fold belly structure. It is worth to mention that such 2-fold symmetry also applies to the basic state. Thus, the classical rotational, axisymmetric CCF is replaced by the 2-CCF due to the symmetry breaking in the transversal magnetic field (see Section 3.2).
Consider wSPI, one can identify qualitative similar wavy characteristics, whereby the topological differences are weaker in comparison to wTVF. As the classical SPI structures already rotate without magnetic field, this also remains for wSPI in the presence of magnetic fields. However, also wSPI remains structural almost identical, and their modulation amplitude becomes enforced at certain cross sections of constant θ in presence of a transversal field component.
The basic state (CCF) loses stability either to (w)TVF (Figure 5) or (w)SPI (Figure 6). In case of wTVF, a supercritical pitchfork of revolution bifurcation breaks the axial translational symmetry Tα (8.3) and parameterizes the axial location of the steady axisymmetric family of solutions. On the other hand, (w)SPIs appear in a supercritical Hopf bifurcation which breaks the mirror symmetry Kz (8.2), resulting in a pair of symmetrically related spiral states, one with left-handed and the other with right-handed winding. It is worth to mention that a finite transverse magnetic field T breaks the axisymmetry Rφ, resulting in a basic state with discrete symmetry Rπ, i.e., with azimuthal wave number m=2 (section 3.3.1). Note that this m=2 state is stationary, which is in contrast to a generic breaking of an SO2 symmetry, resulting in a rotating wave.
3.3.4 Stabilization of basic state: point of higher co-dimension
For TCS, a crucial point of special interest in the Re1−Re2 phase diagram is the bicritical, co-dimension two point, γ, at which the primary bifurcation thresholds of TVF and SPI intersect/converge and as thus exchange the stability of the primary bifurcating solution [20]. For instance, in absence of a magnetic field and given wavenumber k=3.927, the γ-point is given by Re1≈95.25 and Re2≈−73.69 [28] (Figures 7(a) and 8(a)); for smaller Re2 (stronger counter-rotation), SPI bifurcate stable, and for larger Re2 (less counter-rotation and co-rotation), TVF bifurcate stable respectively with increasing Re1 (Figure 8(a)). Thus, with increasing Re1, the 2nd solution bifurcates as unstable solution. Note that the γ-point is quite sensitive to the axial wavenumber; for different k, it moves toward higher Re1 and smaller Re2. This bifurcation scenario also holds for wTVF and wSPI in oblique and transversal magnetic fields.
Figure 7.
(a) Influence of an axial magnetic field on the location of the bifurcation thresholds of TVF and SPI out of CCF in the Re1−Re2 diagram. Blue (full) lines refer to TVF and orange (dashed) ones to SPI. Magenta line with points indicates the moving of the γ-point with sz towards smaller Re2. Further shown are the bifurcation thresholds Re1,c with increasing sz2 for fixed (b) Re2=−100 and (c) Re2=0, respectively, obtained by full non-linear simulations.
Figure 8.
Bifurcation curves for the onset of (w)TVF (blue) and (w)SPI (orange) in Re2−Re1 diagram under (a) no magnetic field (N) sx=0sz=0, (b) an applied transverse field (T) sx≠0sz=0, (c) an applied axial field (a) sx=0sz≠0, and (d) an applied oblique field (O) sx≠0sz≠0. Simulation results of full nonlinear Eqs. (11) consider short periodic TCS with k=3.927 (cf. [6, 28]).
It is worth to mention that the curves in Figures 7(a) and 8 present “marginal stability” thresholds. These are obtained for a fixed wavenumber k=3.927, which differ from the ‘critical’ thresholds. The latter is defined by the (smallest) Reynolds number, for which the flow structure first appears, which crucially depends on k. Thus, the intersection points of the curves for TVF and SPI (e.g., Figure 7(a)) are critical values, with respect to the fixed wavenumber.
Figure 7(a) illustrates an example for the influence of different axial magnetic field strength sz on the bifurcation thresholds for TVF and SPI and the variation on the location of this bicritical γ point. These are results of linear calculations and were obtained with a shooting method [7]. Solid blue (dashed orange) lines refer to TVF (SPI). Note that always the lower (in Re1,c) bifurcating vortex solution is stable and the upper one is unstable. As discussed in Section 3.2, increasing sz stabilizes the CCF basic state, both, against TVF and SPI, i.e., the upward shift of the respective threshold curves. In general, the strength of the stabilization effect is larger for helical SPI vortex structures than for toroidally closed TVF (e.g. Figure 8). As a result, the γ-point moves toward more negative Re2 (i.e., toward the left in Figure 7(a)). Thus, it is observed that the Re2 region of primary, stable bifurcating TVF expands with increasing sz.
The same scenario also applies for transversal and oblique magnetic field configuration, whereby the magnitude of the shift is biggest for a pure axial magnetic field and almost equal for a pure transversal and oblique magnetic field (slightly larger for the latter). Figure 7(b,c) displays the bifurcation thresholds Re1,c for fixed outer cylinder rotations, Re2=−100 and Re2=0 with increasing sz2, resulting from fully nonlinear simulations. The slopes in Figure 7(b,c) depend crucially on Re2. However, in general, the slopes for SPI are steeper than for TVF. For small to moderate field parameters, TVF is the primary stable bifurcating structure while SPI bifurcate secondarily and unstable. For stronger fields, the structures interchange stability and bifurcation order.
Summarizing, for fixed Re2, an axial field shifts the onset of helical SPI more strongly than the one of toroidally closed TVF. We furthermore found this behavior to appear for all parameters and magnetic fields: either for axial, transversal, and oblique magnetic fields (Figure 8). In this case, that for transversal and oblique fields, all structures are wavy-like modulated; wTVF and wSPI, respectively (see also Section 3.2).
In pure axial or pure magnetic fields, the upward shifts of the bifurcation values, the stabilization of the basic state (CCF or 2-CCF), i.e., the onset of the primary stable flow state (w)TVF or (w)SPI grows linearly (Figure 7(b,c)) with the squared magnetic field parameters sx2 or sz2, respectively [6]. Thereby, the magnitude of stabilization shift depends crucially on Re2, in general is largest close to outer cylinder at rest.
3.3.5 Marginal stability thresholds: Stabilization of the basic state
Before going into detail about the stabilization of the basic state it is worth, once more to mention, that the basic state can be different in presence of a magnetic field, whether the field contains a symmetry breaking transversal component or not. For pure axial field, one finds the classical CCF which becomes replaced by the 2-fold symmetric 2-CCF for fields with finite transversal component (see also Section 3.3.1).
Figure 8 presents bifurcation thresholds for vortex flows under the influence of different magnetic fields. Shown is the Re2−Re1 diagram, along which the CCF state loses stability. For comparison, (a) shows the thresholds in absence of any field sx=0sz=0, while the same thresholds are shown for (b) a transversal sx≠0sz=0, (c) an axial sx=0sz≠0, and (d) an oblique field sx≠0sz≠0. These results were obtained by nonlinear simulations with the full field ferrohydrodynamical Eqs. (11), considering a short periodic TCS with k=3.927. See also Figure 7, for linear stability analysis, in which the shift of the higher co-dimension γ-point is shown together with stabilization of CCF in an axial field.
For all bifurcation thresholds of TVF, SPI, wTVF, and wSPI upward shifts of the bifurcation values, Re1,c, as shown in Figure 8, are observed. However, the exact onsets of the different primary bifurcating structures differ, either depending on the particular applied field as well as on Re2.
In absence of a magnetic field (Figure 8(a)) and sufficiently strong counter rotation, Re2≤−73 [28], the basic state, CCF, first becomes unstable to a helical solution, SPI, in a supercritical Hopf bifurcation breaking Kz, which results in a pair of symmetrically related, degenerated spiral states, one with left-handed winding and the other with right-handed winding. For considered Re2-parameter range in Figure 8, the bifurcating SPI has an azimuthal wavenumber m=1 (or m=−1 in case of the right-handed SPI). Note that the wavenumber m of the primary bifurcating SPI increases with stronger counter-rotating cylinders [20]. On the other hand, for Re2≥−73, the basic state loses stability to the toroidally closed TVF. The supercritical pitchfork of revolution bifurcation breaks Tα, and the steady axisymmetric family of solutions are parametrized by their axial location.
Considering a pure axial magnetic field, the scenario remains mainly unchanged (Figure 8(c)); still SPI and TVF appear as primary supercritical flow states with an intersection of the bifurcation thresholds at a specific value Re2. As already seen in Figure 7, increasing field strength sz in a pure axial field tends to shift the co-dimension γ-point toward more negative Re2. The same also holds for increasing sx in a pure transversal field, with significant minor effect on Re2 variation. Thus, the magnetic fields tend to render the corresponding basic states more unstable to (w)TVF and more stable to (w)SPI, respectively.
Having a field with finite transversal component, either pure transversal (Figure 8(b)) or oblique orientated (Figure 8(d)), the primary instabilities, with TVF being replaced by wTVF, and SPI being replaced by wSPI, respectively, just as the basic states were modified from CCF to 2-CCF due to the presence of the magnetic field. Their onset Re1,c occurs at higher Re1 values for any given Re2 value compared to the onset of classical TVF and SPI in the absence of a magnetic field. Thus, the magnetic fields alter the basic states and make them more robust to instabilities. The level of stabilization is greater for a pure axial field compared to a pure transversal field, but strongest, when an oblique field is imposed (bifurcation values Re1,c grow with field parameters sx,sz).
An oblique magnetic field results in the strongest upward shift of bifurcation values Re1,c, and also, this shift depends on Re2, which interestingly is largest about Re2=0, outer cylinder at rest. Furthermore, the topological different flow structures of (w)TVF and (w)SPI are shifted about the same amplitude, with a slightly lager effect on helical (w)SPI for the considered Re2−Re1 diagram. This should be kept in mind regarding later consideration involving field modulation via internal dependence, presented in Section 5.
So, while a pure axial fields shift the bifurcation onsets more strongly (i.e. to larger values Re1) than pure transversal ones (Figure 8), the latter one has a stronger influence on the nonlinear spatio-temporal properties and dynamics of the bifurcating vortex structures (Figure 6). Combining both portions to an oblique field enforces both effects. One finds either larger shifts as well as stronger modulated flow structures (Figure 6(e,f)).
3.4 Resume
An external applied magnetic field (independent of its orientation) stabilizes the CCF (2-CCF) basic state: for a given value of the outer Reynolds number Re2, the bifurcation thresholds Re1,c for vortex structures are shifted to higher values of the inner Reynolds number Re1. This holds likewise for toroidally closed (w)TVF, as well for helical (w)SPI. Thereby, these shifts are typically linear in the squared field parameters sx2 or sz2, respectively. Moreover, stabilization is in general stronger for helical solutions (w)SPI than for toroidally closed ones (w)TVF. Additionally, an external applied magnetic field can, depending on its orientation, change the spatio-temporal structure of the classical TVF and of SPI. The reason is the stimulation of additional modes entering the axial and azimuthal Fourier decomposition of the flow. A pure axial magnetic field keeps the symmetry properties unchanged and does not change the flow structures. It only alter the flow amplitudes, i.e., the magnitude of up-shift in the onsets and magnitude of modes after bifurcation. On the other hand, as soon as a finite transversal field component is present (pure transversal or oblique field), all flow structures become always 3-dimensional and therefore qualitatively different to the classical TVF and SPI states. In particular, they are wavy-like modulated and the pure structures do not exist anymore. Such a field component breaks the system symmetry and therefore alters either basic state (CCF changes to 2-CCF) and supercritical instabilities. As a result, one finds new solutions, field generated, nonrotating wTVF, contrary to the classical rotating ones. The vortex tubes of these wTVF are periodically expanded and constricted in azimuthal θ-direction, and this deformation pattern is stationary. In general, the bifurcation thresholds Re1,c for non-axisymmetric, helical m=1 (w)SPI are stronger modified in presence of a magnetic field than those for toroidally closed m=0 (w)TVF. The upward shift towards larger values of Re1 of the former is bigger than for the latter. One direct effect is the shift of the co-dimension point toward stronger counter-rotating frequencies. This effect has been found to show in axial, transversal, and oblique magnetic fields [6, 9].
Finally it is worth to mention that the magnetic field (transversal or oblique field) generated wTVF differ crucially from the classic wTVF [20, 33, 35, 45] that bifurcates secondary out of TVF at relatively high Re1,c in absence of a magnetic field. The present generated wTVF arises via a primary bifurcation directly out of the 2-CCF basic state.
4. Agglomeration and elongational flow effects
4.1 Introduction and motivation
Aside the general difficulties to incorporate the magnetization of any magnetic field into the governing equations of motion, it is also crucial to consider the specific properties of the magnetic particles. However, most mathematical models describing the flow of complex magnetic fluids typically assume noninteracting magnetic particles with a small volume or only point-size objects. However, real ferrofluids consist of a suspension of particles with a finite size in an almost ellipsoid shape. Moreover, they involve particle-particle interactions resulting in agglomeration; in particular, they tend to form chains of various lengths. One possibility to include such effects and to come close to the realistic situation for ferrofluids is to consider the effect of elongational flow. Mathematically, such flow is incorporated by the symmetric part of the velocity gradient field tensor, which could be scaled by a so-called, material dependent, transport coefficientλ2 [10, 37, 38]. As discussed in the previous section, magnetic fields tend to change the bifurcation threshold of the supercritical unstable flows and, depending on their orientation, resulting in structural modifications. Consider elongational flow, these modifications are expected to change the obtained states at least qualitatively.
4.2 Modified internal magnetic field
In order to consider such elongational flow (i.e. agglomeration, chain formation, etc.), one has to determine the relationship between the magnetization M→, the magnetic field H→, and the velocity u→. One possibility is to consider an additional dependence of the magnetization on the symmetric part of the velocity gradient S=∂iuj+∂jui/2 [10, 38]. Thus, the magnetization (4) becomes modified to
M→−M→eq=cNG→E11
where G=Ω→×H→+λ2SH→.
Onsager symmetry relations then require [12, 39] that additional forces involving λ2 appear in the off-equilibrium momentum balance of the ferrofluid.
Now, using the latest modified magnetization (9), the magnetization part in (1) can be eliminated, similar to previous discussion in Section 2.2. We thus have the following ferrohydrodynamic equation of motion (including elongational flow effects) [12]:
∂t+u→⋅∇u→+∇pM−∇2u→=−cN2H→∇⋅G→+H→×∇×F→/2E12
where pM is the dynamic pressure incorporating all magnetic terms which can be written as gradients, i.e., ∇H→⋅H→ and ∇H→eq−H→⋅M→. See also (5) for comparison without considered elongational flow.
Assuming that the internal magnetic field is equal to the external imposed magnetic field (it is known as a leading-order approximation [24], see also following section 5) which is sufficiently good for a ‘first’ numerical investigation for the effect of elongational flow. Therefore, (10) can be simplified:
In this approach, the magnetic field and all the magnetic properties of the ferrofluid could be influenced by the velocity field, the transport coefficient λ2, and the magnetic field Niklas parameters sx and sz as described before (7). However, it is worth to mention that in this approach, the material dependent transport coefficient λ2 is a parameter [12], which can be scaled by the term describing the elongational flow and is independent of H→. Thus, it only appears in combination with the symmetric component of the velocity field tensor. However, real ferrofluids show H→ dependence in the microscopic aspect of the dynamics of ferrofluids with chains or nonspherical particles [25].
4.3 Linear stability behavior and analysis
Before we come to the influence of λ2 onto the fully nonlinear Eqs. (11), it is worth to have a look at the linear marginal stability threshold of the different flow structures. Therefore, in the following, we first investigate the linear growth behavior of either stationary axisymmetric m=0 (TVF) and nonaxisymmetric m=1 (SPI) oscillatory vortex flow of a ferrofluid in the TCS that is placed in homogeneous axial magnetic fields H→ext=Hze→z (Figure 1). In the literature [4, 5, 7, 40, 41, 42], one can find many theoretical works that analyze (via linear stability analysis) the influence of rotational symmetric magnetic fields, i.e., axial, azimuthal, and radial ones on the flow of a ferrofluid in the TCS. However, for a symmetry breaking (e.g., transverse), magnetic field such linear stability analysis is not known for the corresponding basic state at such a field configuration. For a detailed description regarding the linearization around the CCF and solution of the linearized NSE, we refer to the Appendix in [43].
4.3.1 Marginal stability thresholds
In order to discuss how the transport coefficient λ2 influences the growth of m=0 and m=1 vortex perturbations, we determine the stability of the CCF basic state against these perturbations, considering a magnetic field H→ext=Hze→z, applied in pure axial direction.
Figure 9 presents the bifurcation thresholds in the Re2−Re1 diagram for either axisymmetric, m=0, and nonaxisymmetric, m=1, vortex flow in magnetic fields at (a,b) Hext=67.7kA/m and (c) Hext=135.4kA/m for two different axial wavenumbers. Note that the representative axial wavenumber chosen in (a) k=2.8274 is as detected in experiments [10]; (b) k=3.4558. Insets in Figure 9(b,c) show corresponding SPI frequencies at the bifurcation thresholds. As already discussed in Section 3.3.4, the CCF basic state is stabilized against the two types of vortex perturbations in any magnetic field (the lowest thick solid curve in each plot in Figure 9 refers to the case without any magnetic field H = 0 (c.f. Figure 8(a) for non-linear calculation)). Consider elongational flow (λ2≠0), the basic state becomes destabilized in comparison with the case λ2=0. Regarding Figure 9(a,b), both bifurcation thresholds for m=0 and m=1 vortices are shifted to lower values of Re1. Thus, for the axial wavenumber of k=2.8274, the λ2 term reduces the general stabilization effect caused by an applied magnetic field. For parameters as presented in Figure 9(a,b), this holds for the whole range of presented Re2, even when the changes of the stability boundaries are quite weak, in particular for corotating cylinders. It is important to mention that it never over-compensate the field induced stabilization. Figure 9(c) illustrates in general a strong parameter dependence of the elongational flow effect. While at Re2<22 for TVF and Re2<4 for SPI, increasing λ2 also has destabilizing effect, and for corresponding larger Re2, the effect is just opposite. Increasing λ2 even enforce the stabilization due to the applied magnetic field.
Figure 9.
Stability thresholds of the CCF basic state against vortex growth of (a) m=0 (TVF) and (b) m=1 (SPI) in the Re2−Re1 diagram of control parameters. Thin solid and dashed lines with symbols show the results for λ2=0 and λ2=0.8 for a magnetic field Hext=67.7kA/m. The axial wavenumber is k=2.8274. (c) Same curves, with wavenumber k=3.4558 and magnetic field is Hext=135.4kA/m. The lowest thick curve in (a) and (b) refers in each figure to the threshold in absence of a magnetic field Hext=0. Insets show corresponding SPI frequencies at the bifurcation thresholds.
It is worth to mention that, in contrast to the stabilization of the CCF state by a magnetic field, the λ2-induced modifications of the linear growth of vortices depend significantly on their axial wavenumber k as will be discussed below.
4.3.2 Magnetic field dependence of thresholds
Figure 10(I) illustrates the variation of the bifurcation thresholds m=0 (TVF) and for m=1 (SPI) with the magnetic field H→ext=Hze→z, incorporating the elongational flow effect. The top row of figures refers to Re2=0, and the bottom one to Re2=−100 and wavenumbers are k=2.8274 (left) and k=3.1415 (right), respectively. Ignoring the effect at small Hext, all thresholds increase linearly with Hext. The slopes are given by ∂Re1,stabH/∂H and vary slightly either with λ2 and k. For outer cylinder at rest, Re2=0, (Figure 10(I)(a)) the slopes of the m=0 and m=1 curves are comparable, they differ significantly for counter-rotating cylinders, Re2=−100 (Figure 10(I)(b)). For the latter, this results in intersections of the m=0 and m=1 stability boundaries. Such intersections have been observed also in experiments and also in full nonlinear simulations [6] with λ2=0.
Figure 10.
(I) Variation of the (a) m=0 (TVF) and (b) m=1 (SPI) bifurcation thresholds and (II) corresponding SPI frequencies ω with magnetic field Hext for different λ2 as indicated. The top row refers to Re2=0 and the bottom one to Re2=−100. The wavenumbers are (1) k=2.8274 and (2) k=3.1415, respectively.
In general, modifications of the Hext field dependence of the bifurcation thresholds due to λ2 tend to increase, when the cylinders are stronger counter-rotating (Figure 10(I)). Re2 becomes more negative, and thus, the strain increases. On the other side, they also depend on the axial wavenumber k (see also Section 4.3.4). Considering fixed k, the modifications due to λ2 increase with increasing magnetic field Hext while remaining qualitatively the same. However, the general tendency visible in Figure 11, considering elongational effect, is to shift the stability boundaries downward in Re1. Concrete amounts depend on the parameters k and Re2. Note that this relative small destabilizing elongational effect is always overcompensated by the significant stronger general stabilization effect of the magnetic field itself.
Figure 11.
(1) Stability boundaries of the CCF basic state against the growth of m=0 (TVF) (blue) and m=1 (SPI) vortex flow (orange) in the k−Re1 diagram for λ2=0 (solid lines) and λ2=0.8 dashed lines. (2) corresponding marginal spiral frequencies ω at the m=1 thresholds. Further parameters: (a) Re2=0 and (b) Re2=−100. The (axial) magnetic field is constant with Hext=80kA/m. Results from linearized equations [43].
4.3.3 Frequencies of helical SPI states
Instead of looking at the growth rates—the real part of the eigenvalue of these solutions (previous Sections), one can also consider the marginal spiral frequency ω—the imaginary part of the marginal eigenvalue at the bifurcation threshold of m=1 (SPI) vortex flow. Figure 10(II) shows the frequency, corresponding to the m=1 bifurcation thresholds presented in Figure 10(I). Main effect is the increase of the SPI frequencies with the field Hext, accordingly to the increase of onset Reynolds numbers with Hext (Figure 10(I)). Thus, the threshold frequencies are larger for Re2=−100 than for Re2=0. In general, the effect induced by elongational flow is minor, significant smaller, and only secondary. Analog to the slightly down shift in the onsets (Figure 10(I)) also the SPI frequencies are slightly reduced (Figure 10(II)). Strong counter-rotating cylinders reduce in general the frequency for finite λ2 (Figure 10(II)), while for larger wavenumbers, a finite elongational flow λ2≠0 at Re2=0 it can be slightly decreased.
A general observation is that all modifications of the spiral frequencies due to finite λ2≠0 are most pronounced for the smaller wavenumbers (here k=2.8274 in Figure 10) in counter-rotating cylinders. The field dependence is virtually not much modified due to variations in λ2, and the growth rate of the imaginary part of the marginal eigenvalue with increasing λ2 is quite similar to that of its real part.
4.3.4 Axial wavenumber dependence
Following, we will first discuss the situation for a (moderately) larger wavenumber, taking k=3.4558 as one specific, representative example. For larger value k (Figure 9(c)), the situation becomes more complicate. Incorporating the elongational flow scaled by the parameter λ2 increases and decreases the growth rate, with respect to λ2=0, and thus shifts the stability boundaries, up or down, depending on Re2.
Figure 11 illustrates the relation between the wavenumber dependence of the stability boundaries and of the spiral frequencies. This is of particular interest in experiments, since the wavenumbers of some vortex structures show a strong variation when changing the magnetic field [9, 44]. Reindl et al. [9] also observed hysteresis between vortex flows with different k when increasing and decreasing the Reynolds numbers Re1. The bifurcation thresholds (Figure 11, top row Re2=0 and the bottom row Re2=−100; the (axial) magnetic field is fixed at Hext=80kA/m) for m=0 and m=1 vortex flow structures in the k−Re1 diagram change when the elongational term λ2≠0 is considered in the linearized ferrofluid equations [43]. Right column of Figure 11(2) presents the wavenumber dependence of the corresponding marginal spiral frequencies at the onsets for m=1 SPI.
Typically, the thresholds for λ2≠0 (in Figure 11 the case of λ2=0, lines with squares) are shifted downward in Re1 at smaller wavenumbers k and upward at larger k. Resulting in a fixed magnetic field, the elongational effect destabilizes (stabilizes) the CCF state against growth of vortices with small (large) wavenumbers, either for m=0 or m=1 perturbations. Consequently, the stability onsets cross for different values of λ2 (λ2=0, λ2=0.8 in Figure 11) at a specific value of kcross the wavenumber.
It is worth to mention that the crossover values of kcross for m=0 and m=1 perturbations also depend sensitive on other the control parameters, e.g., Re2. Based on the observations in Figure 11, it is fair to assume that, depending on the wavenumber k, the γ-point (Section 3.3.4) will change.
Similar to the stability thresholds, also the spiral frequencies (Figure 11(2)) increase with increasing λ2 when k is large and decrease for smaller k as a result of switching on the elongational effect. Thereby, the crossing wavenumber, kcross, separating the increase and decrease in ω is slightly smaller than the one for crossing the marginal thresholds (Figure 11(a1,b1)). Following the speculation by Odenbach and Müller [10] on the physical origin of the elongational effect, one can give a reasonable explanation for such an observed wavenumber dependence of the λ2 effect on the growth rates of vortex flow. In [10], the authors argue that the microscopic origin of a finite λ2 might be attributed to the finite asphericity of the colloids, i.e., the presence of (short) particle chains in the ferrofluid. Small axial wavenumber k will allow these chains more easily to remain aligned in the direction of the magnetic field in vortex flows instead in a flow with large axial wavenumber and large azimuthal vorticity. Larger wavenumber means smaller wavelength and therefore more higher probability to unalign in the given field. Thus, since the magneto-rotational dissipation in the latter flow is larger than in the former one, it requires larger centrifugal forces, i.e., larger values of Re1, to drive the growth of a large k flow. However, it is important to state that this argumentation also holds some weakness: In the presence of chains, one could expect also for small k vortices an increase of the threshold Reynolds number Re1. On the contrary, incorporating the elongational term λ2≠0 in the magnetization balance changes the momentum balance of the vortex flow such that the growth of vortices is enhanced for k<kcross (left in Figure 10) and reduced for k>kcross (right in Figure 10) relative to the reference case of λ2=0.
4.4 Nonlinear analysis
4.4.1 Bifurcation thresholds
Following, we will consider the full nonlinear ferrohydrodynamical equation of motions incorporating elongational flow effects (11) in order to investigate the influence of a finite transport coefficient λ2 at different applied magnetic fields.
The general modification in the bifurcation threshold of primary instabilities ((w)TVF) and ((w)SPI) due to finite values λ2≠0 remains qualitative as detected and described before in the linear stability analysis in pure axial field. However, the absolute values are slightly different. A direct comparison for axial field (Figures 9-11) results in variation about one to two percent between the calculated onsets using either the linear analysis or full nonlinear equations [6]. While the different field orientations result in different magnitudes of stabilization of the basic states CCF and 2-CCF, respectively, as discussed in Section 3.3.4, a finite value λ2 does not depend on different field orientations. The variation effect with λ2 is almost identical for all considered field configuration. However, in general, the modifications are slightly larger under oblique magnetic fields. Depending on other system parameters, e.g., wavenumber k, cylinder rotation speeds Re1,Re2, the effect is either stabilizing or destabilizing with respect to the situation λ2=0 (Figure 10(c) and Figure 11). Considering the effect of elongational flow, the critical bifurcating point of primary centrifugal instability will be shifted either up or down. As detected in the linear analysis, the potential destabilizing effect of λ2never overcomes the general stabilizing effect of the external applied magnetic field, i.e., comparing λ2=0 and any value λ2≠0.
4.4.2 Flow structure modifications
Aside the modification of the primary stability thresholds, there are further flow structural changes/modifications which can be observed for finite value λ2. Figure 11 shows, for different field configurations T: sx=0.1sz=0, A: sx=0sz=0.1, O: sx=0.1sz=0.1, how the maxima of vorticity η vary for either the full flow structure Δη≔ηλ2≠0−ηλ2=0 (Figure 12(a)), and reduced 2-fold subspace (Figure 12(b)) varies with λ2 for the basic state configuration Re1=60Re2=0. Remember, that only for sx=0, one has the classical CCF, and for any sx≠0, a modified basic state is present [6, 9, 37] (i.e., symmetry breaking effect of transversal field). Although for a given configuration (A, T, O), the absolute maxima of these values are significantly different, and their relative differences Δ scale is linear with λ2 starting at zero. Lines in Figure 12 are calculated by linear regression fitting to the numerical data points. The effect of λ2 is decoupled from the magnetic field configuration and with this the considered flow structure.
Figure 12.
Variation of vorticity maxima with λ2 for (a) maxη=ηλ2−ηλ2=0 at full solution and (b) maxηm=2=ηλ2−ηλ2=0 at only 2-fold symmetric subspace (empty for A) for basic flow (CCF and 2-CCF, respectively) at Re1=60,Re2=0 and different field orientations: Axial (A) sx=0.0,sz=0.6; transverse (T) sx=0.6,sz=0.0 and oblique (O) sx=0.6,sz=0.6. Lines are linear fits.
Either for classical “pure” CCF or field modified 2-CCF, the variation in maxΔη grows linear with λ2 (Figure 12(a)) (note the significant smaller amplitude in the m=2 subspace). Further confirmation for this is the fact that an analog linear growth can be also found when considering other parameters, at which a supercritical flow state exists [37]. This holds for toroidally closed structures (w)TVF, as well as for helical flow states (w)SPI. Only the absolute values in maxΔη are significant larger.
It is important to mention that the slopes in Figure 12 (for Osx=0.1sz=0.1) not only depend on the applied field orientation (O) but also depend on the strength of the applied magnetic field. With increasing the magnetic field strength sx,sz, the slopes become steeper. This scenario is illustrated, exemplary for 2-CCF basic state in transversal field configuration in Figure 13.
Figure 13.
As Figure 11 for basic state 2-CCF at Re1=60,Re2=0. Variation with λ2 of vorticity maxima maxη=ηλ2−ηλ2=0 for transversal field (T) at different field strength sx as indicated.
Moreover, similar characteristics hold for the different primary bifurcating structures. Figure 14 elucidates this effect for the three different field configurations A, T, and O as indicated (see also Figures 4 and 6 in [37] for another example considering a transversal magnetic field in finite system configuration). For a given field configuration, the slopes increase with λ2, when increasing the magnetic field strength sx. In general, slopes are smallest for transversal and largest for oblique fields, respectively. The variation in the slopes is continuous with increasing the field strength, but not linear; for larger field strength, the variation in the slopes becomes smaller.
Figure 14.
Supercritical (w)TVF solution at Re1=100,Re2=0. Variation of vorticity maxima maxη=ηλ2−ηλ2=0 for (a) transversal (T), (b) axial (A) and (c) oblique (O) magnetic field (strength as indicated). Lines are linear fits.
4.4.3 Wavenumber modification and selection
As seen in the former discussion (section 4.3.2), an elongational flow can have strong effect on the axial wavenumber k. Consider the full equations (and periodic boundary conditions), one obtains similar results as presented in Figure 10 for linear analysis of the marginal stability thresholds for TVF and SPI. As described before, the agreement between the results of linear and full equations is very close (they differ less than 2 percent).
However, in order to see such wavenumber modification from another perspective, it makes sense to study a finite system with sufficient large aspect ratio enclosed by fixed end plates (assumed to be stationary) on top and bottom. As in a finite system, the axial wavenumber k is natural selected by the system itself (depending on the various system parameters), it is possible to see the direct influence, considering a finite value λ2 on to the flow structures and their axial wavenumbers k.
Starting with a classical TVF as initial flow state, with n=22 vortices, which corresponds to an axial wavenumber k=3.41 in the bulk (illustrated on the most left border in Figure 15), we consider a situation with an applied pure transversal magnetic field which is continuously increased for two values, either ignoring any elongational flow effect, λ2=0 (Figure 14(a)), or considering such effect λ2=0.2 (Figure 15(b)) (note 0.2 is quite small). In any case, either for λ2=0 and for λ2=0.2, by increasing the magnetic-field strength, sx, the number of vortices is reduced, which leads to an increase in the axial wavelength (i.e. decrease in axial wavenumber k) due to elimination of vortices in the bulk. Therefore, the corresponding axial wavenumber k in the bulk becomes reduced. Remember that the initial state is a classical, stationary TVF, which becomes modified to stationary wTVF as soon as sx≠0. For the field generated wTVF, the wavy-like modulation is not visible in Figure 15. However qualitative, the same scenario also holds for either axial and oblique magnetic fields, with changes in the corresponding parameters sx,sz at which the wavenumber k changes, i.e., vortices disappear.
Figure 15.
Contours of the radial velocity component u at mid-gap with variation of field strength sx at Re1=88.35,Re2=0 for (a) λ2=0 and (b) λ2=0.2. Left border in each plot gives the initial state, a classical TVF with n = 22 vortices corresponding to an axial wavenumber k = 3.41 in the bulk. Note that for sx=0, only wavy-like modified flow states exist. Red (yellow) [dark gray (light gray)] corresponds to positive (negative) values. The max (min) level is ±9.98.
In general, the exact field strength sx for vortex elimination and/or the number of destroyed vortices depends on λ2. In the example of Figure 15, the initial flow has 22 vortices with an axial wavenumber k=3.41. Neglecting elongational effects λ2=0 (Figure 15(a)) until the parameter sx approximately reaches the value about 0.7, the flow remains stable, and the number of vortices in the flow can be reduced to n=20 with k=3.04. Further increasing sx to about 0.85, the number of vortices can be more reduced to n=16 with k=2.58. Hereafter, no further changes in the number of vortices n appear, before reaching the boundary threshold, sx≈0.9. Consider elongational flow λ2=0.2 (Figure 15(b)), the behavior of the vortex state is similar to the case neglecting such effects. The flow remains stable before the parameter sx reaches about 0.75. At the same time, the number of vortices is also reduced in analog manner to n=20 with k=3.04. Hereafter, as before, the number of vortices n remains unchanged until the parameter sx increases to the boundary threshold, sx≈0.95. Consider values λ2>0.7, the flow structure (for given parameters) is not affected by the magnetic-field strength. A flow with n=22 vortices corresponding to axial wavenumber k=3.41 in the bulk is remaining in the bulk.
An interesting phenomenon, which also has been reported in experimental studies [9] as in numerical simulations, is the fact that one or two vortex pairs can be eliminated. However so far, there is no physical explanation for this observation/phenomenon. Analog to the presented scenario in Figure 15 for pure transversal magnetic field, one also finds similar reduction scenarios for pure axial and oblique magnetic fields. As to expect, the absolute values of field strength, sx,sz, at which the elimination of vortices happen is different, depending on the applied field.
4.5 Resume
Real ferrofluids consist of a suspension of particles with finite size and mainly almost ellipsoid shape, as well as there are particle-particle interactions that tend to form chains of various lengths and tend to agglomerate in general. One possibility in order to make a more realistic approach for ferrofluids and to come close to such a configuration in a real ferrofluids is to consider the effect of elongational flow. Such a flow becomes incorporated by the symmetric part of the velocity gradient field tensor in the nondimensional ferrohydrodynamic equations of motion (11) and scales by a material dependent transport coefficient λ2. By such an assumption, the latter is given by a pure scalar number, describing the strength/influence of all former mentioned properties.
Consider either linear stability analysis or fully nonlinear equations, one finds that incorporating such elongational flow can have either stabilization or destabilization effect (i.e., modification in the supercritical primary bifurcation thresholds) with respect to the reference situation in absence of such flow λ2=0. The modifications due to finite value λ2 are qualitative similar and in the same range for different considered magnetic fields, slightly largest in case of an oblique field compared to pure axial or transversal ones. Typically, the modifications due to λ2≠0 are larger for helical, nonaxisymmetric, m=1, (w)SPIs than for toroidally closed, axisymmetric, m=0, (w)TVF. Thus, the regim of primary stable bifurcating (w)TVF increases. This observation is similar and goes in hand with the previously described general stabilization of any applied magnetic field (Section 3). There we saw that the stability thresholds of helical flow states are shifted further than those of toroidal closed flow structures. However, it is important and fundamental to state that the elongational effect λ2≠0 is never big enough to over compensate the general, field-induced stabilization. Only the magnitude/amplitude in shift can be reduced. Further, it is worth to mention that linear stability analysis is limited to axial magnetic fields; as for fields with finite transversal component, no basic state can be derived, simulations of full nonlinear equations confirm similar effects. The modifications due to λ2 are sensitive to various different system parameters (e.g., Reynolds numbers, radius ratio, etc.), in particular sensitive to the axial wavenumber k. Here, smaller k are strongly effected than larger ones. Thus, considering finite size TCS, a considered elongational flow results in modifications of flow structures in the sense of changing number of vortices in the bulk, which coincides with a different wavenumber selection. Thus, depending on the particular wavenumber, the presence of any elongational effect λ2≠0 can be either stabilizing or destabilizing for the basic state. In finite system, this results in different selection and/or modification of number of vortex pairs within the annulus, depending on the parameter.
Regarding the flow structural properties, a finite transport coefficient, λ2≠0, does not change the qualitative shape of the flow pattern. Only quantitative differences can be observed. Essentially, all structural properties remain as generated/induced due to the external applied magnetic field. Thus for instance, the azimuthal mode-2 symmetry originated in the symmetry-breaking transverse magnetic field is preserved; only the strength of different flow quantities is modified (in general they increase) according to finite λ2. Mainly, the flow states become linearly enforced, e.g., speaking in quantity of the azimuthal vorticity, which holds for both the basic flow (CCF, 2-CCF) and the primary instabilities the supercritical centrifugal unstable flows ((w)TVF, (w)SPI).
5. Modified field dependence - internal magnetization
5.1 Introduction and motivation
Aside the properties and behavior of the magnetic particles itself and their interaction with the applied magnetic field (see last sections), the correct assumption of the applied magnetic field is another crucially point in describing magnetic flows and their dynamics. The frequently used, most common assumption is to consider the internal magnetic field within the magnetic fluid (e.g., a ferrofluid) to be equal to the external applied field. Following, we will see that this is only a leading-order approximation. By accounting for the ferrofluid’s magnetic susceptibility, a uniform externally imposed magnetic field is modified by the presence of the magnetic fluid (ferrofluid) within the annulus. The modification to the magnetic field appears in an radial dependence 1/r2 and their magnitude scales with the susceptibility χ. To be concrete, such modification to the imposed magnetic field can be substantial, also for ferrofluids typically used in laboratory experiments. These have significant consequences on the structural properties and stability of the basic states, as well as on the supercritical (primary) bifurcating solutions.
5.2 Modified internal magnetic field
As already seen in the previous sections, for solving the equation of motion (1), we have to consider the magnetic field within the annular gap between the cylinders. In the last sections, we used the simplest approach with the assumption to take the magnetic field to be identical to the applied external field [4, 6]. However, this simple approach is only a leading-order approximation, and depending on the magnetic susceptibility of the ferrofluid, the magnetic field in the gap is modified.
Assuming infinitely long cylinders, the magnetic boundary conditions are
H→=H→ext−Mre→rE14
at r=r1 and r=r2, where H→ext is the homogeneous external applied magnetic field in the absence of the ferrofluid-filled annulus and Mr is the radial component of the magnetization M→.
The following ansatz satisfies a solenoidal field, ∇⋅H→=0:
and then the boundary conditions defining the external field are.
H→r=r1=H→ext+a1−b1/r12cosθ+a2−b2/r12sinθe→r+a2+b2/r12cosθ−a1−b1/r12sinθe→θ and H→r=r2=H→ext+a1−b1/r22cosθ+a2−b2/r22sinθe→r+a2+b2/r22cosθ−a1−b1/r22sinθe→θ.
To calculate the four unknown constants a1,a2,b1,b2 from the boundary conditions (12), we need the radial component of the magnetization which results from substituting (13) into the magnetization (3):
where HextT stands for the transverse component of H→ext. Using the continuity of the magnetic field condition on the cylinders, the coefficients can be calculated. Substituting them into (13), finally the resulting magnetic field is given by
where K=2+χ2−χ2η2. A special case of (15) when ΩcN=0. This corresponds to an equilibrium magnetization and has been already observed in an earlier work by Odenbach and Müller [38]. They considered the limit of an equilibrium magnetization and used a similar derivation of the modification of the magnetic field for a ferrofluid between two cylinders.
Based on (15), one can use the following magnetic field in the equation of motion:
H→=22+ζKHx[1−ζ/r2cosθe→r−ζ/r2sinθe→θ+Hze→zE18
the dimensionless scalar parameter
ζ=χ2+χr12E19
characterizes the strength of the radial field dependence ζ. The ferrohydrodynamical equation of motion (6) and the Niklas parameters sx and sz (7) remain unchanged/untouched from such considered radial field dependence.
5.3 Flow structural modifications
5.3.1 Primary instabilities
As we learned in previous sections, a magnetic field with a finite transverse component modifies either basic state or primary bifurcating flow structures in a wavy-like manner. More details of these states have been reported in numerical simulations [6], for which a homogeneous internal magnetic field equal to the external applied one has been supposed, as well as in real experimental studies [9]. The question now is how does the former derived radial field dependence ζ affects the stability of the basic states (CCF and 2-CCF) and the characteristics of the bifurcating solutions ((w)TVF and (w)SPI).
We begin our discussion by considering a purely axial applied magnetic field (A), with fixed field strength sz=0.6, which leaves the flows, in particular its structural properties unaffected (TVF and SPI). Figure 16 shows the influence of the radial field parameter, ζ, (17) (exemplary chosen (a) ζ=0, (b) ζ=0.4 and (c) ζ=0.8), on the bifurcation thresholds for TVF and SPI in an axial magnetic field. As in absence of any field dependence, ζ=0 (Figure 16(a)), the bifurcations continue being supercritical. Main effect on the bifurcation curves is that with increasing ζ, the thresholds for the onset of both TVF and SPI are shifted to higher Re1 (Figure 16(b,c)). Thereby, the threshold for the toroidally closed TVF being shifted significantly further (to larger values Re1). Note the different scaling on the abscissa in Figure 16(c). Aside from the enhanced stabilization of the CCF basic state on top of the field introduced stabilization itself, one observes an important consequence to the relative shifts in the two bifurcation curves. Consequently, the bifurcation curve of TVF is shifted more than the one of SPI, the co-dimension-two point, γ, is also shifted to more positive Re2 (c.f. Figure 15(a,b)). Interestingly, this is just the opposite way of the shift with increasing field strength sz in a pure axial magnetic field (the same holds for sx in transversal field), which we saw in Figure 7 with moving the γ-point to the left toward more negative values in Re2. Consider larger values, e.g., ζ=0.8 (Figure 15(c)), one sees that the co-dimension-two γ-point is even shifted to Re2>50 out of the considered parameter space. As a result, the primary instability is always the helical SPI, even with the two cylinders in strong co-rotation. This observation is crucially different from classical TCS, both without magnetic fields and also with assumption of an applied uniform magnetic field. Thus, the radial field dependence, induced by the susceptibility of the ferrofluid-filled annulus, is seen to cause major changes in the quantitative and qualitative characteristics of the flow instabilities.
Figure 16.
Bifurcation curves for the onset of TVF (blue) and SPI (orange) in Re2−Re1 diagram for different ζ-values (a) ζ=0, (b) ζ=0.4 and (c) ζ=0.8 with an applied axial magnetic field (A) with sz=0.6. ((a) is a replot of Figure 7(c) for comparison). Note the different scale on the ordinate in (c). Simulation results of full nonlinear equations (11) consider a short periodic TCS with axial wavenumber k=3.927 (same hold for Figures 15 and 16).
Consider a transversal magnetic field, field strength sx=0.6 (Figure 17) instead of an axial one, the variation by increasing ζ remains virtually and qualitatively the same, aside that now all structures are wavy-like modulated (2-CCF, wTVF, wSPI). One finds the same preference, in the sense of stronger upshifts for toroidally closed wTVF with the same result of a shift in the γ-point to the right to larger values in Re2. Thus also here, for sufficiently strong ζ, e.g., ζ=0.8 in Figure 17(c), the helical wSPI is the primary stable bifurcating solution in all presented parameter space. Comparing both, pure axial and pure transversal field (at same field strength), the up-shift, i.e., the stabilization is slightly larger in axial fields. This holds for toroidally closed (w)TVF and helical (w)SPI solutions.
Figure 17.
Bifurcation curves for the onset of wTVF (blue) and wSPI (orange) in parameter space for different ζ values (a) ζ=0, (b) ζ=0.4, and (c) ζ=0.8 with an applied transversal magnetic field (T) sx=0.6. ((a) is a replot of Figure 7(b) for comparison).
Finally, this effect is even more pronounced when an oblique magnetic field is applied as illustrated in Figure 18 for field strength sx=0.6,sz=0.6. Note also the different scale on the ordinate in Figure 17(c) highlighting the stronger shifts due to such field configuration. One can speculate that this enforcement is due to the axial component of the oblique field leading to a base state with w having axial shear, and the radial field dependence ζ enhances the w component of the 2-CCF base state. The combination of the base state azimuthal shear and this axial shear favors the helical instability. Close to the onset, the variation of the maxima of all flow components u,v,w is almost linear.
Figure 18.
Bifurcation curves for the onset of wTVF (blue) and wSPI (orange) in Re2−Re1 diagram for different ζ-values (a) ζ=0, (b) ζ=0.4, and (c) ζ=0.8 with an applied oblique magnetic field (O) with sx=0.6=s−z. ((a) is a replot of Figure 7(d) for comparison). Note different scaling in (c).
5.3.2 Modified flow structures
Figure 19 illustrates the flow structural properties of wTVF and wSPI in either a transverse magnetic field T (a,b) and an oblique magnetic field O (c,d), both with assumed radial field dependence ζ=0.8. Note, the flow structures of TVF and SPI under a pure axial applied magnetic field are not shown, as they virtually do not differ from the ones in absence of a magnetic field. Only the absolute values (e.g. in η) are different. Comparing these flow structures to the ones for ζ=0 cases as earlier presented in Figure 6 (Section 3.3.3), one sees the modulations in all structures to be enhanced/enforced. In general, the waviness (axial variation in azimuthal direction) of all flow structures is enforced/enhanced by the radial dependent magnetic field. However, the field-induced wTVF always remains a stationary, nonrotating, and phase-pinned structure, irrespective of whether the radial field dependence is considered or not. On the other hand, wSPIs bifurcate as degenerated solutions, either as a left- or right-winding SPI, depending on initial conditions. Their corresponding frequencies remain mainly constant over a wide range and do not differ much (about 20% or less) from those of SPI. Comparing the stimulated modes for flow structures when ζ=0.8 (bottom in Figure 19) to those with ζ=0 (bottom in Figure 6), one sees general larger mode amplitudes in the additional stimulated modes for consideration of an internal field dependence ζ≠0. These larger amplitudes in the field induced modes are also the reason for the increase in waviness.
Figure 19.
(See also Figures 2 and 6 for comparison to flow states without radial field dependence. The radial field dependence is considered to be ζ=0.8. Top row: Azimuthal vorticity isosurfaces η over two axial wavelengths with pure transversal applied magnetic field (T): (a) wTVF at Re2=0, Re1=−130, (b) wSPI at Re2=−150, Re1=190; and with oblique applied magnetic field (O): (c) wTVF at Re2=0, Re1=−150, and (f) wSPI at Re2=−150, Re1=−200 Red (yellow) isosurfaces correspond to positive (negative) values as indicated. The bottom row shows the mode amplitudes um,n of the radial velocity field u corresponding to the structures above over the m-n-plane.
5.4 Resume
Applying a magnetic field across a container of ferrofluid changes the ferrofluid’s susceptibility and results in modification of the magnetic field structure within the container. So far, most theories commonly neglect such modification, either for convenience or simplicity, and at the end, there is also possibility of numerical simulations in appropriate time or feasibility at all.
However, it is a matter of fact that an external applied magnetic field is modified, depending on the magnetic susceptibility. Experimental studies [10, 38] confirmed the modifications also for susceptibilities of commonly used ferrofluids. Consider such modifications to the imposed magnetic field effects the basic state, its stability, as well as the primary (supercritical) bifurcating solutions. This holds for any field configuration, with only differences in the strength of modifications, i.e., the stabilization (its magnitude; up-shift, move of the bifurcation thesholds to larger control parameters). On the other hand, the flow structural properties remain qualitative unaffected by such an radial field dependence, only change in quantitative manner. They are determined only by the general kind and orientation of the applied magnetic field. Thus, a pure axial field does not have any effect on the flow structural properties, leaving them as in absence of any field, i.e., pure TVF and pure SPI. As described before, the presence of a finite transversal field component results in a wavy-like modulated flow state. This is independent/untouched if the magnetic field is considered to be homogeneous or inhomogeneous with a radial field dependence. However, consider stronger modifications to the imposed magnetic field, i.e., larger values of the radial field dependence parameter ζ have the effect of further strengthening the already present waviness of the structures. This holds for (w)TVF and (w)SPI alike. An important observation is the different influence of the radial field dependence ζ≠0 onto topological closed and helical flow states. In general, the effect of modification in the primary bifurcation thresholds is larger for toroidally closed flow structures compared to helical ones. As a result, their intersection point, i.e., co-dimension γ-point is moved toward positive Re2, which is just the opposite direction, as it moves with increasing the magnetic field strength (for any given field configuration, e.g. Figure 7). The parameter space, the area for which (w)SPI bifurcate primary stable out of the basic state, is growing, even into region of co-rotating cylinders, i.e. Re2>0, making wSPI the favorable flow state in the system.
6. Conclusions
Although there has been significant progress in the understanding of ferrofluidic flows in the past decades, their fully theoretical understanding is still limited. A key reason for this is computational high-expensive simulations, which even increase with continuation in refining the models to make them more realistic/appropriate, e.g., finite size effects, internal field modifications, etc. However, these model improvements are essential and necessary to uncover the full potential of ferrofluids, which already are used in a wide field of applications. The use in the car industry within antishocks to provide adjustable damper settings is probably one of the most prominent examples. However, either experiments and numerics are essential to mutually prove each other and combined to provide new inside of potential applications.
This chapter deals with a rotating ferrofluid between two concentric cylinders, i.e., ferrofluid in Taylor-Couette geometry. Therefor we consider different improvements in the model describing such magnetic particles in a flow under external applied field. Concluding the discussions, the main points can be summarized as follows:
Any applied magnetic field changes the bifurcation thresholds of primary solutions—the thresholds of both primary solutions, TVF and SPI, are shifted to larger parameter values—the basic state becomes stabilized against perturbations.
Magnetic fields with a finite transversal component break the system symmetry and renders all flow states to be inherently 3D. As a result, the pure flow states TVF and SPI cease to exist and are substituted by their wave-like cousins, wTVF and wSPI. Thereby, the generated wTVF differs crucially from the classical ones, which rotate in azimuthal direction. Instead, the magnetic field induced wTVFs are nonrotational, phase-pinned flow solutions and develop belly shape structure at specific azimuthal positions.
Assuming elongational flow, in order to consider finite-size effect of the magnetic particles and their interaction (agglomeration, particle-particle interaction and chain formation procedures), one detects further changes in the bifurcation threshold of the solutions. Also, this can result in further stabilization or destabilization. However, the effect of such elongational flow is never big enough to over-compensate the general, field-induced stabilization. Only the magnitude/amplitude varies. The flow states remain structural as in absence of elongational flow.
Accounting for the changes in ferroluid’s susceptibility in the container due to an external magnetic field, the real internal field can be approximated by 1/r2. As before the flow, structural properties remain qualitative similar. Bifurcation thresholds are shifted upward, toward larger Re1, whereby toroidally closed (w)TVFs are in general preferred and shifted more than helical (w)SPI. The result is an increase in the regime of primary bifurcating (w)SPI.
Conflict of interest
The author declares no conflict of interest.
A. Appendix
A.1 Numerical method
There are various different ways to solve the ferrohydrodynamic equations of motions (Eqs. (5) and (10)), the latter including elongational flow). All here presented results have been obtained by a code (G1D3 [6, 7]), which combines a finite difference method of second order in (r,z) and time (explicit) with spectral decomposition in θ [6]:
frθzt=∑m=−mmaxmmaxfmrzteinkθ,E20
where f denotes one of uvwp. For the parameter regimes investigated here, mmax=8 provides adequate accuracy. For discretization, a forward time, centered space (FTCS) algorithm is used on staggered grids in the r−z plane following the procedure of Hirt et al. [45]. In particular, a homogeneous grid with discretization length δr=δz=0.05 and time steps δt<3800 has been used. Here, δt defines the time step between two iterations via FTCS algorithm of the system of coupled equations for the amplitudes fmrzt of the azimuthal normal modes −mmax≤m≤mmax. For diagnostic purposes, the complex mode amplitudes fm,nrt were calculated obtained from a second axial Fourier decomposition:
fmrzt=∑nfm,nrteinkz,E21
where k=2π/λ is the axial wavenumber.
Acknowledgments/Funding
S. Altmeyer is a Serra Húnter Fellow. This work was supported by the Spanish Government grant FIS2017–85794-P.
Nomenclature and abbreviations
TCS
Taylor-Coutte System
NSE
Navier–Stokes equation
CCF
Circular Couette flow [basic flow]
2-CCF
2-fold CCF [basic flow in transversal field]
TVF
Taylor vortex flow
SPI
Spiral vortex flow (spiral)
wTVF
wavy TVF
wSPI
wavy SPI
\n',keywords:"magnetic fluid, complex system, ferrofluid, instabilities, bifurcation, symmetry-breaking, vortex, growth rate, finite-size effect, internal magnetic field, particle interaction, agglomeration, elongational flow",chapterPDFUrl:"https://cdn.intechopen.com/pdfs/63437.pdf",chapterXML:"https://mts.intechopen.com/source/xml/63437.xml",downloadPdfUrl:"/chapter/pdf-download/63437",previewPdfUrl:"/chapter/pdf-preview/63437",totalDownloads:396,totalViews:0,totalCrossrefCites:0,dateSubmitted:"April 3rd 2018",dateReviewed:"July 13th 2018",datePrePublished:"November 5th 2018",datePublished:null,dateFinished:null,readingETA:"0",abstract:"When studying ferrofluidic flows, as one example of magnetic flow dynamics, in terms of instability, bifurcation, and properties, one quickly finds out the additional challenges magnetic fluids introduce compared to the investigation of “classical”, “ordinary” shear flows without any kind of particles. Approximation of ferrofluids as fluids including point-size particles or, more realistic fine size particles, the relaxation times of the magnetic particle, their interaction between each other, i.e., the agglomeration and chain forming effects, and the interaction/response between any external applied field and the internal magnetization are just few examples of challenges to overcome. Further dependence on the considered model system, the direction of the external applied magnetic field (homogeneous or inhomogeneous) is crucial, as it can break the system symmetry and thus generate new solutions. As a result, the classical Navier–Stokes equations become modified to the more complex ferrohydrodynamical equation of motion, incorporating magnetic field and magnetization of the fluid itself, which typically makes numerical simulations expensive and challenging. This chapter provides an overview of the tasks/difficulties from describing and simulating magnetic particles, their interaction, and thus finally resulting modification in rotating flow structures and in particular instabilities and bifurcation behavior.",reviewType:"peer-reviewed",bibtexUrl:"/chapter/bibtex/63437",risUrl:"/chapter/ris/63437",signatures:"Sebastian Altmeyer",book:{id:"7516",title:"Pattern Formation and Stability in Magnetic Colloids",subtitle:null,fullTitle:"Pattern Formation and Stability in Magnetic Colloids",slug:"pattern-formation-and-stability-in-magnetic-colloids",publishedDate:"January 8th 2020",bookSignature:"Nicola?s O. Rojas",coverURL:"https://cdn.intechopen.com/books/images_new/7516.jpg",licenceType:"CC BY 3.0",editedByType:"Edited by",editors:[{id:"208723",title:"Dr.",name:"Nicolás O.",middleName:null,surname:"Rojas",slug:"nicolas-o.-rojas",fullName:"Nicolás O. Rojas"}],productType:{id:"1",title:"Edited Volume",chapterContentType:"chapter",authoredCaption:"Edited by"}},authors:null,sections:[{id:"sec_1",title:"1. Introduction",level:"1"},{id:"sec_2",title:"2. System and theoretical description",level:"1"},{id:"sec_2_2",title:"2.1 Magnetization equations",level:"2"},{id:"sec_3_2",title:"2.2 Ferrohydrodynamic equations of motions",level:"2"},{id:"sec_5",title:"3. Ferrofluidic flow under finite magnetic field",level:"1"},{id:"sec_5_2",title:"3.1 Introduction and motivation",level:"2"},{id:"sec_6_2",title:"3.2 Flow structures",level:"2"},{id:"sec_6_3",title:"3.2.1 Flow structures in absence of a magnetic field",level:"3"},{id:"sec_7_3",title:"3.2.2 Flow structures—mode coupling in magnetic fields",level:"3"},{id:"sec_9_2",title:"3.3 Effect of finite magnetic fields",level:"2"},{id:"sec_9_3",title:"3.3.1 Basic state and symmetries",level:"3"},{id:"sec_10_3",title:"3.3.2 Bifurcation behavior",level:"3"},{id:"sec_11_3",title:"3.3.3 Vortex structures and mode contents",level:"3"},{id:"sec_12_3",title:"3.3.4 Stabilization of basic state: point of higher co-dimension",level:"3"},{id:"sec_13_3",title:"3.3.5 Marginal stability thresholds: Stabilization of the basic state",level:"3"},{id:"sec_15_2",title:"3.4 Resume",level:"2"},{id:"sec_17",title:"4. Agglomeration and elongational flow effects",level:"1"},{id:"sec_17_2",title:"4.1 Introduction and motivation",level:"2"},{id:"sec_18_2",title:"4.2 Modified internal magnetic field",level:"2"},{id:"sec_19_2",title:"4.3 Linear stability behavior and analysis",level:"2"},{id:"sec_19_3",title:"4.3.1 Marginal stability thresholds",level:"3"},{id:"sec_20_3",title:"4.3.2 Magnetic field dependence of thresholds",level:"3"},{id:"sec_21_3",title:"4.3.3 Frequencies of helical SPI states",level:"3"},{id:"sec_22_3",title:"4.3.4 Axial wavenumber dependence",level:"3"},{id:"sec_24_2",title:"4.4 Nonlinear analysis",level:"2"},{id:"sec_24_3",title:"4.4.1 Bifurcation thresholds",level:"3"},{id:"sec_25_3",title:"4.4.2 Flow structure modifications",level:"3"},{id:"sec_26_3",title:"4.4.3 Wavenumber modification and selection",level:"3"},{id:"sec_28_2",title:"4.5 Resume",level:"2"},{id:"sec_30",title:"5. Modified field dependence - internal magnetization",level:"1"},{id:"sec_30_2",title:"5.1 Introduction and motivation",level:"2"},{id:"sec_31_2",title:"5.2 Modified internal magnetic field",level:"2"},{id:"sec_32_2",title:"5.3 Flow structural modifications",level:"2"},{id:"sec_32_3",title:"5.3.1 Primary instabilities",level:"3"},{id:"sec_33_3",title:"5.3.2 Modified flow structures",level:"3"},{id:"sec_35_2",title:"5.4 Resume",level:"2"},{id:"sec_37",title:"6. Conclusions",level:"1"},{id:"sec_43",title:"Conflict of interest",level:"1"},{id:"sec_39",title:"A. Appendix",level:"1"},{id:"sec_38",title:"A.1 Numerical method",level:"1"},{id:"sec_39",title:"Acknowledgments/Funding",level:"1"},{id:"sec_42",title:"Nomenclature and abbreviations",level:"1"}],chapterReferences:[{id:"B1",body:'Rosensweig RE. Ferrohydrodynamics. Cambridge University Press, Cambridge; 1985. 344 pp. Journal of Fluid Mechanics, 200, 597-599 DOI: 10.1017/S0022112089220773'},{id:"B2",body:'Hart JE. Dynamics of Atmospheres and Oceans. 2006;41. DOI: 10.1016/j.dynatmoce.2006.03.001'},{id:"B3",body:'Hart JE, Kittelman S. Dynamics of Atmospheres and Oceans. 2006;41:139. DOI: 10.1017/S0022112075001309'},{id:"B4",body:'Niklas M. Zeitschrift fuer Physik B: Condensed Matter. 1987;68:493. DOI: 10.1007/BF01471080'},{id:"B5",body:'Niklas M, Müller-Krumbhaar H, Lücke M. Journal of Magnetism and Magnetic Materials. 1989;81:29. DOI: 10.1016/0304-8853(89)90225-4'},{id:"B6",body:'Altmeyer S, Hoffmann C, Leschhorn A, Lücke M. Physical Review E. 2010;82:016321. DOI: 10.1103/PhysRevE.82.016321'},{id:"B7",body:'Leschhorn A, Lücke M, Hoffmann C, Altmeyer S. Physical Review E. 2009;79:036308. DOI: 10.1103/PhysRevE.79.036308'},{id:"B8",body:'Reindl M, Leschhorn A, Lücke M, Odenbach S. Journal of Physics: Conference Series. 2009;149:012109. DOI: 10.1088/1742-6596/149/1/012109'},{id:"B9",body:'Reindl M, Odenbach S. Physics of Fluids. 2011;23:093102. DOI: 10.1063/1.3633341'},{id:"B10",body:'Odenbach S, Müller HW. Journal of Magnetism and Magnetic Materials. 2005;289:242. DOI: 10.1016/j.jmmm.2004.11.069'},{id:"B11",body:'Mahle S, Ilg P, Liu M. Physical Review E. 2008;77:016305. DOI: 10.1103/PhysRevE.77.016305'},{id:"B12",body:'Müller HW, Liu M. Physical Review E. 2001;64:061405. DOI: 10.1103/PhysRevE.64.061405'},{id:"B13",body:'Debye PJW. Polar Molecules. New York: Dover; 1929 https://hdl.handle.net/2027/uc1.$b45989'},{id:"B14",body:'McTague JP. The Journal of Chemical Physics. 1969;51:133. DOI: 10.1063/1.1671697'},{id:"B15",body:'Martsenyuk MA, Raikher YL, Shliomis MI. Soviet Physics - JETP. 1974;38:413 http://www.jetp.ac.ru/cgi-bin/dn/e_038_02_0413.pdf'},{id:"B16",body:'Gilles de GP, Prost J, editors. The Physics of Liquid Crystals. Oxford: Clarendon; 1983. ISBN: 9780198517856'},{id:"B17",body:'Leschhorn A, Embs JP, Lücke M. Journal of Physics: Condensed Matter. 2006;18:S2633. DOI: 10.1088/0953-8984/18/38/S07'},{id:"B18",body:'Embs JP, May S, Wagner C, Kityk AV, Leschhorn A, Lücke M. Physical Review E. 2006;73:036302. DOI: 10.1103/PhysRevE.73.036302'},{id:"B19",body:'Taylor GI. Philosophical Transactions of the Royal Society of London A. 1923;223:289. DOI: 10.1098/rsta.1923.0008'},{id:"B20",body:'Chossat P, Iooss G. The Couette-Taylor Problem. Berlin: Springer; 1994. DOI: 10.1007/978-1-4612-4300-7'},{id:"B21",body:'Langevin P. Annales de Chimie Physique. 1905;5(70)'},{id:"B22",body:'Embs J, Müller HW, Wagner C, Knorr K, Lücke M. Physical Review E. 2000;61:R2196. DOI: 10.1103/PhysRevE.61.R2196'},{id:"B23",body:'Odenbach S, Gilly H. Journal of Magnetism and Magnetic Materials. 1995;152(123). DOI: 10.1016/0304-8853(95)00442-4'},{id:"B24",body:'Altmeyer S, Lopez JM, Do Y. Physical Review E. 2012;85:066314. DOI: 10.1103/PhysRevE.85.066314'},{id:"B25",body:'Zubarev AY, Skakova Y. Physical Review E. 2000;61:5415. DOI: 10.1103/PhysRevE.61.5415'},{id:"B26",body:'Golubitsky M, Stewart I, Schaeffer D. Singularities and Groups in Bifurcation Theory II. New York: Springer; 1988. ISBN: 978-1-4612-4574-2'},{id:"B27",body:'Golubitsky M, Stewart I. SIAM Journal on Mathematical Analysis. 1986;17:249. DOI: 10.1137/0517023'},{id:"B28",body:'Ch H, Altmeyer S, Pinter A, Lücke M. New Journal of Physics. 2009;11:053002. DOI: 10.1088/1367-2630/11/5/053002'},{id:"B29",body:'Davey A, Di Prima RC, Stuart JT. Journal of Fluid Mechanics. 1968;31:17. DOI: 10.1017/S0022112068000029'},{id:"B30",body:'Iooss G. Journal of Fluid Mechanics. 1986;173:273. DOI: 10.1017/S0022112086001179'},{id:"B31",body:'Altmeyer S, Do Y, Marques F, Lopez JM. Physical Review E. 2012;86:046316. DOI: 10.1103/PhysRevE.86.046316'},{id:"B32",body:'Coles D. Journal of Fluid Mechanics. 1965;21:385. DOI: 10.1017/S0022112065000241'},{id:"B33",body:'Wereley ST, Lueptow RM. Journal of Fluid Mechanics. 1998;364:59. DOI: 10.1017/S0022112098008969'},{id:"B34",body:'Jones CA. Journal of Fluid Mechanics. 1981;102:249. DOI: 10.1017/S0022112081002620'},{id:"B35",body:'Akonur A, Lueptow RM. Physics of Fluids. 2003;15:947. DOI: 10.1063/1.1556615'},{id:"B36",body:'King GP, Li Y, Lee W, Swinney HL, Marcus PS. Journal of Fluid Mechanics. 1984;141:365. DOI: 10.1017/S0022112084000896'},{id:"B37",body:'Altmeyer S, Lopez JM, Do Y. Physical Review E. 2013;88:013003. DOI: 10.1103/PhysRevE.88.013003'},{id:"B38",body:'Odenbach S, Müller HW. Physical Review Letters. 2002;89:037202. DOI: 10.1103/PhysRevLett.89.037202'},{id:"B39",body:'Müller HW, Hahn D, Liu M. Journal of Physics: Condensed Matter. 2006;18:S2623. DOI: 10.1088/0953-8984/18/38/S06'},{id:"B40",body:'Chang MH, Chen CK, Weng HC. International Journal of Engineering Science. 2003;41:103. DOI: 10.1016/S0020-7225(02)00183-0'},{id:"B41",body:'Singh J, Bajaj R. Journal of Magnetism and Magnetic Materials. 2005;294:53. DOI: 10.1016/j.jmmm.2004.10.123'},{id:"B42",body:'Singh J, Bajaj R. International Journal of Mathematics and Mathematical Sciences. 2005;3727. DOI: 10.1155/IJMMS.2005.3727'},{id:"B43",body:'Altmeyer S, Leschhorn A, Ch H, Lücke M. Physical Review E. 2013;87:053010. DOI: 10.1103/PhysRevE.87.053010'},{id:"B44",body:'Reindl M, Odenbach S. Experiments in Fluids. 2010;50:375. DOI: 10.1007/s00348-010-0940-y'},{id:"B45",body:'Hirt CW, Nichols BD, Romero NC. NASA STI/Recon Technical Report N; 1975. http://adsabs.harvard.edu/abs/1975STIN...7532418H'}],footnotes:[],contributors:[{corresp:"yes",contributorFullName:"Sebastian Altmeyer",address:"sebastian.andreas.altmeyer@upc.edu",affiliation:'
Castelldefels School of Telecom and Aerospace Engineering, Universitat Politècnica de Catalunya, Barcelona, Spain
'}],corrections:null},book:{id:"7516",title:"Pattern Formation and Stability in Magnetic Colloids",subtitle:null,fullTitle:"Pattern Formation and Stability in Magnetic Colloids",slug:"pattern-formation-and-stability-in-magnetic-colloids",publishedDate:"January 8th 2020",bookSignature:"Nicola?s O. Rojas",coverURL:"https://cdn.intechopen.com/books/images_new/7516.jpg",licenceType:"CC BY 3.0",editedByType:"Edited by",editors:[{id:"208723",title:"Dr.",name:"Nicolás O.",middleName:null,surname:"Rojas",slug:"nicolas-o.-rojas",fullName:"Nicolás O. Rojas"}],productType:{id:"1",title:"Edited Volume",chapterContentType:"chapter",authoredCaption:"Edited by"}}},profile:{item:{id:"190158",title:"Dr.",name:"Anita-Ioana",middleName:null,surname:"Visan",email:"anita.visan@inflpr.ro",fullName:"Anita-Ioana Visan",slug:"anita-ioana-visan",position:null,biography:null,institutionString:null,profilePictureURL:"//cdnintech.com/web/frontend/www/assets/author.svg",totalCites:0,totalChapterViews:"0",outsideEditionCount:0,totalAuthoredChapters:"1",totalEditedBooks:"0",personalWebsiteURL:null,twitterURL:null,linkedinURL:null,institution:null},booksEdited:[],chaptersAuthored:[{title:"Composite Coatings Based on Renewable Resources Synthesized by Advanced Laser Techniques",slug:"composite-coatings-based-on-renewable-resources-synthesized-by-advanced-laser-techniques",abstract:"This chapter reviews the progress and perspectives of composite materials in the form of thin films based on renewable resources for biofabrication of a new generation of medical implants with antibacterial properties. The chapter starts with an overview of the types of renewable materials that were currently studied and of the unique properties which make them perfect candidates for numerous bio‐related applications. A briefing of recent progresses in the field of advanced laser synthesis of composites from renewable and sustainable materials, as well as the relevant results in our researches is made. The discussion spans composite coatings based on renewable resources, [e.g., chitosan (CHT) and lignin (Lig)] as biomaterials intended for metallic implants. A particular attention is given to lignin synthesis in the form of thin films due to its ability to functionalize the medical implant surface while preserving the similar composition and the structural properties of the raw, renewable biomaterial. We focused on recent technological advancements (e.g., matrix‐assisted pulsed laser evaporation (MAPLE) and Combinatorial‐MAPLE) which have brought the spotlight onto renewable biomaterials, by detailing the relevant engineering data of processing. This chapter concludes that the extensions of advanced laser techniques are viable fabrication methods of a new generation of metallic implants.",signatures:"Anita‐Ioana Visan, Carmen‐Georgeta Ristoscu and Ion N. Mihailescu",authors:[{id:"17636",title:"Dr.",name:"Ion",surname:"Mihailescu",fullName:"Ion Mihailescu",slug:"ion-mihailescu",email:"ion.mihailescu@inflpr.ro"},{id:"60531",title:"Dr.",name:"Carmen-Georgeta",surname:"Ristoscu",fullName:"Carmen-Georgeta Ristoscu",slug:"carmen-georgeta-ristoscu",email:"carmen.ristoscu@inflpr.ro"},{id:"190158",title:"Dr.",name:"Anita-Ioana",surname:"Visan",fullName:"Anita-Ioana Visan",slug:"anita-ioana-visan",email:"anita.visan@inflpr.ro"}],book:{title:"Composites from Renewable and Sustainable Materials",slug:"composites-from-renewable-and-sustainable-materials",productType:{id:"1",title:"Edited Volume"}}}],collaborators:[{id:"20844",title:"Dr.",name:"Edith",surname:"Mäder",slug:"edith-mader",fullName:"Edith Mäder",position:null,profilePictureURL:"//cdnintech.com/web/frontend/www/assets/author.svg",biography:null,institutionString:null,institution:null},{id:"156271",title:"Dr.",name:"Evans",surname:"Egwim",slug:"evans-egwim",fullName:"Evans Egwim",position:null,profilePictureURL:"//cdnintech.com/web/frontend/www/assets/author.svg",biography:null,institutionString:null,institution:{name:"Federal University of Technology Minna",institutionURL:null,country:{name:"Nigeria"}}},{id:"189661",title:"Dr.",name:"Ruey Shan",surname:"Chen",slug:"ruey-shan-chen",fullName:"Ruey Shan Chen",position:null,profilePictureURL:"//cdnintech.com/web/frontend/www/assets/author.svg",biography:null,institutionString:null,institution:{name:"National University of Malaysia",institutionURL:null,country:{name:"Malaysia"}}},{id:"189861",title:"Dr.",name:"Adeshina",surname:"Fadeyibi",slug:"adeshina-fadeyibi",fullName:"Adeshina Fadeyibi",position:null,profilePictureURL:"https://mts.intechopen.com/storage/users/189861/images/4684_n.jpg",biography:null,institutionString:null,institution:{name:"Kwara State University",institutionURL:null,country:{name:"Nigeria"}}},{id:"189904",title:"Prof.",name:"Sahrim",surname:"Ahmad",slug:"sahrim-ahmad",fullName:"Sahrim Ahmad",position:null,profilePictureURL:"//cdnintech.com/web/frontend/www/assets/author.svg",biography:null,institutionString:null,institution:null},{id:"190270",title:"Prof.",name:"Denis",surname:"Rodrigue",slug:"denis-rodrigue",fullName:"Denis Rodrigue",position:null,profilePictureURL:"https://mts.intechopen.com/storage/users/no_image.jpg",biography:null,institutionString:null,institution:{name:"Toronto Western Hospital",institutionURL:null,country:{name:"Canada"}}},{id:"190690",title:"Dr.",name:"Hanna",surname:"Brodowsky",slug:"hanna-brodowsky",fullName:"Hanna Brodowsky",position:null,profilePictureURL:"//cdnintech.com/web/frontend/www/assets/author.svg",biography:null,institutionString:null,institution:{name:"Leibniz Institute of Polymer Research",institutionURL:null,country:{name:"Germany"}}},{id:"194440",title:"Prof.",name:"Zinash",surname:"Delebo Osunde",slug:"zinash-delebo-osunde",fullName:"Zinash Delebo Osunde",position:null,profilePictureURL:"//cdnintech.com/web/frontend/www/assets/author.svg",biography:null,institutionString:null,institution:null},{id:"194441",title:"Dr.",name:"Gbabo",surname:"Agidi",slug:"gbabo-agidi",fullName:"Gbabo Agidi",position:null,profilePictureURL:"//cdnintech.com/web/frontend/www/assets/author.svg",biography:null,institutionString:null,institution:null},{id:"194442",title:"Prof.",name:"Peter",surname:"Aba Idah",slug:"peter-aba-idah",fullName:"Peter Aba Idah",position:null,profilePictureURL:"//cdnintech.com/web/frontend/www/assets/author.svg",biography:null,institutionString:null,institution:null}]},generic:{page:{slug:"translation-policy",title:"Translation Policy",intro:"
All Works published by IntechOpen prior to October 2011 are licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported license (CC BY-BC-SA 3.0). Works published after October 2011 are licensed under a Creative Commons Attribution 3.0 Unported license (CC BY 3.0), the latter allowing for the broadest possible reuse of published material.
",metaTitle:"Translation Policy",metaDescription:"Translation of Works - Book Chapters",metaKeywords:null,canonicalURL:"/page/translation-policy",contentRaw:'[{"type":"htmlEditorComponent","content":"
All Works licensed under CC BY-BC-SA 3.0 can be freely translated and used for non-commercial purposes. Works licensed under CC BY 3.0 license can be freely translated and used for both commercial and non-commercial purposes.
\\n\\n
Book Chapters
\\n\\n
All translated Chapters have to be properly attributed in accordance with the requirements included in IntechOpen's Attribution Policy. Besides proper attribution translated sections of Works must include the following sentence: "This is an unofficial translation of a work published by IntechOpen. The publisher has not endorsed this translation".
\\n\\n
Books and all other compilations
\\n\\n
All rights to Books and other compilations are reserved by IntechOpen. The copyright to Books and other compilations is subject to a Copyright separate from any that exists in the included Works.
\\n\\n
A Book in its entirety, or a significant part of a Book, cannot be translated freely without specific written consent by the publisher. Requests for permission can be made at permissions@intechopen.com.
All Works licensed under CC BY-BC-SA 3.0 can be freely translated and used for non-commercial purposes. Works licensed under CC BY 3.0 license can be freely translated and used for both commercial and non-commercial purposes.
\n\n
Book Chapters
\n\n
All translated Chapters have to be properly attributed in accordance with the requirements included in IntechOpen's Attribution Policy. Besides proper attribution translated sections of Works must include the following sentence: "This is an unofficial translation of a work published by IntechOpen. The publisher has not endorsed this translation".
\n\n
Books and all other compilations
\n\n
All rights to Books and other compilations are reserved by IntechOpen. The copyright to Books and other compilations is subject to a Copyright separate from any that exists in the included Works.
\n\n
A Book in its entirety, or a significant part of a Book, cannot be translated freely without specific written consent by the publisher. Requests for permission can be made at permissions@intechopen.com.
\n\n
Policy last updated: 2016-06-09
\n'}]},successStories:{items:[]},authorsAndEditors:{filterParams:{sort:"featured,name"},profiles:[{id:"6700",title:"Dr.",name:"Abbass A.",middleName:null,surname:"Hashim",slug:"abbass-a.-hashim",fullName:"Abbass A. Hashim",position:null,profilePictureURL:"https://mts.intechopen.com/storage/users/6700/images/1864_n.jpg",biography:"Currently I am carrying out research in several areas of interest, mainly covering work on chemical and bio-sensors, semiconductor thin film device fabrication and characterisation.\nAt the moment I have very strong interest in radiation environmental pollution and bacteriology treatment. The teams of researchers are working very hard to bring novel results in this field. I am also a member of the team in charge for the supervision of Ph.D. students in the fields of development of silicon based planar waveguide sensor devices, study of inelastic electron tunnelling in planar tunnelling nanostructures for sensing applications and development of organotellurium(IV) compounds for semiconductor applications. I am a specialist in data analysis techniques and nanosurface structure. I have served as the editor for many books, been a member of the editorial board in science journals, have published many papers and hold many patents.",institutionString:null,institution:{name:"Sheffield Hallam University",country:{name:"United Kingdom"}}},{id:"54525",title:"Prof.",name:"Abdul Latif",middleName:null,surname:"Ahmad",slug:"abdul-latif-ahmad",fullName:"Abdul Latif Ahmad",position:null,profilePictureURL:"//cdnintech.com/web/frontend/www/assets/author.svg",biography:null,institutionString:null,institution:null},{id:"20567",title:"Prof.",name:"Ado",middleName:null,surname:"Jorio",slug:"ado-jorio",fullName:"Ado Jorio",position:null,profilePictureURL:"//cdnintech.com/web/frontend/www/assets/author.svg",biography:null,institutionString:null,institution:{name:"Universidade Federal de Minas Gerais",country:{name:"Brazil"}}},{id:"47940",title:"Dr.",name:"Alberto",middleName:null,surname:"Mantovani",slug:"alberto-mantovani",fullName:"Alberto Mantovani",position:null,profilePictureURL:"//cdnintech.com/web/frontend/www/assets/author.svg",biography:null,institutionString:null,institution:null},{id:"12392",title:"Mr.",name:"Alex",middleName:null,surname:"Lazinica",slug:"alex-lazinica",fullName:"Alex Lazinica",position:null,profilePictureURL:"https://mts.intechopen.com/storage/users/12392/images/7282_n.png",biography:"Alex Lazinica is the founder and CEO of IntechOpen. After obtaining a Master's degree in Mechanical Engineering, he continued his PhD studies in Robotics at the Vienna University of Technology. Here he worked as a robotic researcher with the university's Intelligent Manufacturing Systems Group as well as a guest researcher at various European universities, including the Swiss Federal Institute of Technology Lausanne (EPFL). During this time he published more than 20 scientific papers, gave presentations, served as a reviewer for major robotic journals and conferences and most importantly he co-founded and built the International Journal of Advanced Robotic Systems- world's first Open Access journal in the field of robotics. Starting this journal was a pivotal point in his career, since it was a pathway to founding IntechOpen - Open Access publisher focused on addressing academic researchers needs. Alex is a personification of IntechOpen key values being trusted, open and entrepreneurial. Today his focus is on defining the growth and development strategy for the company.",institutionString:null,institution:{name:"TU Wien",country:{name:"Austria"}}},{id:"19816",title:"Prof.",name:"Alexander",middleName:null,surname:"Kokorin",slug:"alexander-kokorin",fullName:"Alexander Kokorin",position:null,profilePictureURL:"https://mts.intechopen.com/storage/users/19816/images/1607_n.jpg",biography:"Alexander I. Kokorin: born: 1947, Moscow; DSc., PhD; Principal Research Fellow (Research Professor) of Department of Kinetics and Catalysis, N. Semenov Institute of Chemical Physics, Russian Academy of Sciences, Moscow.\r\nArea of research interests: physical chemistry of complex-organized molecular and nanosized systems, including polymer-metal complexes; the surface of doped oxide semiconductors. He is an expert in structural, absorptive, catalytic and photocatalytic properties, in structural organization and dynamic features of ionic liquids, in magnetic interactions between paramagnetic centers. The author or co-author of 3 books, over 200 articles and reviews in scientific journals and books. He is an actual member of the International EPR/ESR Society, European Society on Quantum Solar Energy Conversion, Moscow House of Scientists, of the Board of Moscow Physical Society.",institutionString:null,institution:{name:"Semenov Institute of Chemical Physics",country:{name:"Russia"}}},{id:"62389",title:"PhD.",name:"Ali Demir",middleName:null,surname:"Sezer",slug:"ali-demir-sezer",fullName:"Ali Demir Sezer",position:null,profilePictureURL:"https://mts.intechopen.com/storage/users/62389/images/3413_n.jpg",biography:"Dr. Ali Demir Sezer has a Ph.D. from Pharmaceutical Biotechnology at the Faculty of Pharmacy, University of Marmara (Turkey). He is the member of many Pharmaceutical Associations and acts as a reviewer of scientific journals and European projects under different research areas such as: drug delivery systems, nanotechnology and pharmaceutical biotechnology. Dr. Sezer is the author of many scientific publications in peer-reviewed journals and poster communications. Focus of his research activity is drug delivery, physico-chemical characterization and biological evaluation of biopolymers micro and nanoparticles as modified drug delivery system, and colloidal drug carriers (liposomes, nanoparticles etc.).",institutionString:null,institution:{name:"Marmara University",country:{name:"Turkey"}}},{id:"61051",title:"Prof.",name:"Andrea",middleName:null,surname:"Natale",slug:"andrea-natale",fullName:"Andrea Natale",position:null,profilePictureURL:"//cdnintech.com/web/frontend/www/assets/author.svg",biography:null,institutionString:null,institution:null},{id:"100762",title:"Prof.",name:"Andrea",middleName:null,surname:"Natale",slug:"andrea-natale",fullName:"Andrea Natale",position:null,profilePictureURL:"//cdnintech.com/web/frontend/www/assets/author.svg",biography:null,institutionString:null,institution:{name:"St David's Medical Center",country:{name:"United States of America"}}},{id:"107416",title:"Dr.",name:"Andrea",middleName:null,surname:"Natale",slug:"andrea-natale",fullName:"Andrea Natale",position:null,profilePictureURL:"//cdnintech.com/web/frontend/www/assets/author.svg",biography:null,institutionString:null,institution:{name:"Texas Cardiac Arrhythmia",country:{name:"United States of America"}}},{id:"64434",title:"Dr.",name:"Angkoon",middleName:null,surname:"Phinyomark",slug:"angkoon-phinyomark",fullName:"Angkoon Phinyomark",position:null,profilePictureURL:"https://mts.intechopen.com/storage/users/64434/images/2619_n.jpg",biography:"My name is Angkoon Phinyomark. I received a B.Eng. degree in Computer Engineering with First Class Honors in 2008 from Prince of Songkla University, Songkhla, Thailand, where I received a Ph.D. degree in Electrical Engineering. My research interests are primarily in the area of biomedical signal processing and classification notably EMG (electromyography signal), EOG (electrooculography signal), and EEG (electroencephalography signal), image analysis notably breast cancer analysis and optical coherence tomography, and rehabilitation engineering. I became a student member of IEEE in 2008. During October 2011-March 2012, I had worked at School of Computer Science and Electronic Engineering, University of Essex, Colchester, Essex, United Kingdom. In addition, during a B.Eng. I had been a visiting research student at Faculty of Computer Science, University of Murcia, Murcia, Spain for three months.\n\nI have published over 40 papers during 5 years in refereed journals, books, and conference proceedings in the areas of electro-physiological signals processing and classification, notably EMG and EOG signals, fractal analysis, wavelet analysis, texture analysis, feature extraction and machine learning algorithms, and assistive and rehabilitative devices. I have several computer programming language certificates, i.e. Sun Certified Programmer for the Java 2 Platform 1.4 (SCJP), Microsoft Certified Professional Developer, Web Developer (MCPD), Microsoft Certified Technology Specialist, .NET Framework 2.0 Web (MCTS). I am a Reviewer for several refereed journals and international conferences, such as IEEE Transactions on Biomedical Engineering, IEEE Transactions on Industrial Electronics, Optic Letters, Measurement Science Review, and also a member of the International Advisory Committee for 2012 IEEE Business Engineering and Industrial Applications and 2012 IEEE Symposium on Business, Engineering and Industrial Applications.",institutionString:null,institution:{name:"Joseph Fourier University",country:{name:"France"}}},{id:"55578",title:"Dr.",name:"Antonio",middleName:null,surname:"Jurado-Navas",slug:"antonio-jurado-navas",fullName:"Antonio Jurado-Navas",position:null,profilePictureURL:"https://mts.intechopen.com/storage/users/55578/images/4574_n.png",biography:"Antonio Jurado-Navas received the M.S. degree (2002) and the Ph.D. degree (2009) in Telecommunication Engineering, both from the University of Málaga (Spain). He first worked as a consultant at Vodafone-Spain. From 2004 to 2011, he was a Research Assistant with the Communications Engineering Department at the University of Málaga. In 2011, he became an Assistant Professor in the same department. From 2012 to 2015, he was with Ericsson Spain, where he was working on geo-location\ntools for third generation mobile networks. Since 2015, he is a Marie-Curie fellow at the Denmark Technical University. His current research interests include the areas of mobile communication systems and channel modeling in addition to atmospheric optical communications, adaptive optics and statistics",institutionString:null,institution:{name:"University of Malaga",country:{name:"Spain"}}}],filtersByRegion:[{group:"region",caption:"North America",value:1,count:5766},{group:"region",caption:"Middle and South America",value:2,count:5227},{group:"region",caption:"Africa",value:3,count:1717},{group:"region",caption:"Asia",value:4,count:10367},{group:"region",caption:"Australia and Oceania",value:5,count:897},{group:"region",caption:"Europe",value:6,count:15790}],offset:12,limit:12,total:118189},chapterEmbeded:{data:{}},editorApplication:{success:null,errors:{}},ofsBooks:{filterParams:{topicId:"19"},books:[{type:"book",id:"10736",title:"Neurotoxicity",subtitle:null,isOpenForSubmission:!0,hash:"f3ae592c3bd56dca45f9ce7d02e06714",slug:null,bookSignature:"",coverURL:"https://cdn.intechopen.com/books/images_new/10736.jpg",editedByType:null,editors:null,productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"10881",title:"Drug Repurposing",subtitle:null,isOpenForSubmission:!0,hash:"8ef09a9da770b582c0c64114a19b29c0",slug:null,bookSignature:"",coverURL:"https://cdn.intechopen.com/books/images_new/10881.jpg",editedByType:null,editors:null,productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"10882",title:"Smart Drug Delivery",subtitle:null,isOpenForSubmission:!0,hash:"70c3ce4256324b3c58db970d446ddac4",slug:null,bookSignature:"Dr. Usama Ahmad",coverURL:"https://cdn.intechopen.com/books/images_new/10882.jpg",editedByType:null,editors:[{id:"255360",title:"Dr.",name:"Usama",surname:"Ahmad",slug:"usama-ahmad",fullName:"Usama Ahmad"}],productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"10883",title:"Pain Management",subtitle:null,isOpenForSubmission:!0,hash:"82abad01d1cffb27e341ffd507117824",slug:null,bookSignature:"",coverURL:"https://cdn.intechopen.com/books/images_new/10883.jpg",editedByType:null,editors:null,productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"10884",title:"Bisphenols",subtitle:null,isOpenForSubmission:!0,hash:"d73ec720cb7577731662ac9d02879729",slug:null,bookSignature:"Prof. Pınar Erkekoglu",coverURL:"https://cdn.intechopen.com/books/images_new/10884.jpg",editedByType:null,editors:[{id:"109978",title:"Prof.",name:"Pınar",surname:"Erkekoglu",slug:"pinar-erkekoglu",fullName:"Pınar Erkekoglu"}],productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}}],filtersByTopic:[{group:"topic",caption:"Agricultural and Biological Sciences",value:5,count:14},{group:"topic",caption:"Biochemistry, Genetics and Molecular Biology",value:6,count:3},{group:"topic",caption:"Business, Management and Economics",value:7,count:1},{group:"topic",caption:"Chemistry",value:8,count:7},{group:"topic",caption:"Computer and Information Science",value:9,count:6},{group:"topic",caption:"Earth and Planetary Sciences",value:10,count:7},{group:"topic",caption:"Engineering",value:11,count:15},{group:"topic",caption:"Environmental Sciences",value:12,count:2},{group:"topic",caption:"Immunology and Microbiology",value:13,count:3},{group:"topic",caption:"Materials Science",value:14,count:5},{group:"topic",caption:"Mathematics",value:15,count:1},{group:"topic",caption:"Medicine",value:16,count:24},{group:"topic",caption:"Neuroscience",value:18,count:1},{group:"topic",caption:"Pharmacology, Toxicology and Pharmaceutical Science",value:19,count:2},{group:"topic",caption:"Physics",value:20,count:2},{group:"topic",caption:"Psychology",value:21,count:4},{group:"topic",caption:"Social Sciences",value:23,count:2},{group:"topic",caption:"Technology",value:24,count:1},{group:"topic",caption:"Veterinary Medicine and Science",value:25,count:1}],offset:12,limit:12,total:5},popularBooks:{featuredBooks:[{type:"book",id:"9385",title:"Renewable Energy",subtitle:"Technologies and Applications",isOpenForSubmission:!1,hash:"a6b446d19166f17f313008e6c056f3d8",slug:"renewable-energy-technologies-and-applications",bookSignature:"Tolga Taner, Archana Tiwari and Taha Selim Ustun",coverURL:"https://cdn.intechopen.com/books/images_new/9385.jpg",editors:[{id:"197240",title:"Associate Prof.",name:"Tolga",middleName:null,surname:"Taner",slug:"tolga-taner",fullName:"Tolga Taner"}],equalEditorOne:{id:"186791",title:"Dr.",name:"Archana",middleName:null,surname:"Tiwari",slug:"archana-tiwari",fullName:"Archana Tiwari",profilePictureURL:"https://mts.intechopen.com/storage/users/186791/images/system/186791.jpg",biography:"Dr. Archana Tiwari is Associate Professor at Amity University, India. Her research interests include renewable sources of energy from microalgae and further utilizing the residual biomass for the generation of value-added products, bioremediation through microalgae and microbial consortium, antioxidative enzymes and stress, and nutraceuticals from microalgae. She has been working on algal biotechnology for the last two decades. She has published her research in many international journals and has authored many books and chapters with renowned publishing houses. She has also delivered talks as an invited speaker at many national and international conferences. Dr. Tiwari is the recipient of several awards including Researcher of the Year and Distinguished Scientist.",institutionString:"Amity University",position:null,outsideEditionCount:0,totalCites:0,totalAuthoredChapters:"3",totalChapterViews:"0",totalEditedBooks:"1",institution:{name:"Amity University",institutionURL:null,country:{name:"India"}}},equalEditorTwo:{id:"197609",title:"Prof.",name:"Taha Selim",middleName:null,surname:"Ustun",slug:"taha-selim-ustun",fullName:"Taha Selim Ustun",profilePictureURL:"https://mts.intechopen.com/storage/users/197609/images/system/197609.jpeg",biography:"Dr. Taha Selim Ustun received a Ph.D. in Electrical Engineering from Victoria University, Melbourne, Australia. He is a researcher with the Fukushima Renewable Energy Institute, AIST (FREA), where he leads the Smart Grid Cybersecurity Laboratory. Prior to that, he was a faculty member with the School of Electrical and Computer Engineering, Carnegie Mellon University, Pittsburgh, PA, USA. His current research interests include power systems protection, communication in power networks, distributed generation, microgrids, electric vehicle integration, and cybersecurity in smart grids. He serves on the editorial boards of IEEE Access, IEEE Transactions on Industrial Informatics, Energies, Electronics, Electricity, World Electric Vehicle and Information journals. Dr. Ustun is a member of the IEEE 2004 and 2800, IEC Renewable Energy Management WG 8, and IEC TC 57 WG17. He has been invited to run specialist courses in Africa, India, and China. He has delivered talks for the Qatar Foundation, the World Energy Council, the Waterloo Global Science Initiative, and the European Union Energy Initiative (EUEI). His research has attracted funding from prestigious programs in Japan, Australia, the European Union, and North America.",institutionString:"Fukushima Renewable Energy Institute, AIST (FREA)",position:null,outsideEditionCount:0,totalCites:0,totalAuthoredChapters:"1",totalChapterViews:"0",totalEditedBooks:"0",institution:{name:"National Institute of Advanced Industrial Science and Technology",institutionURL:null,country:{name:"Japan"}}},equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter"}},{type:"book",id:"8985",title:"Natural Resources Management and Biological Sciences",subtitle:null,isOpenForSubmission:!1,hash:"5c2e219a6c021a40b5a20c041dea88c4",slug:"natural-resources-management-and-biological-sciences",bookSignature:"Edward R. Rhodes and Humood Naser",coverURL:"https://cdn.intechopen.com/books/images_new/8985.jpg",editors:[{id:"280886",title:"Prof.",name:"Edward R",middleName:null,surname:"Rhodes",slug:"edward-r-rhodes",fullName:"Edward R Rhodes"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter"}},{type:"book",id:"9027",title:"Human Blood Group Systems and Haemoglobinopathies",subtitle:null,isOpenForSubmission:!1,hash:"d00d8e40b11cfb2547d1122866531c7e",slug:"human-blood-group-systems-and-haemoglobinopathies",bookSignature:"Osaro Erhabor and Anjana Munshi",coverURL:"https://cdn.intechopen.com/books/images_new/9027.jpg",editors:[{id:"35140",title:null,name:"Osaro",middleName:null,surname:"Erhabor",slug:"osaro-erhabor",fullName:"Osaro Erhabor"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter"}},{type:"book",id:"7841",title:"New Insights Into Metabolic Syndrome",subtitle:null,isOpenForSubmission:!1,hash:"ef5accfac9772b9e2c9eff884f085510",slug:"new-insights-into-metabolic-syndrome",bookSignature:"Akikazu Takada",coverURL:"https://cdn.intechopen.com/books/images_new/7841.jpg",editors:[{id:"248459",title:"Dr.",name:"Akikazu",middleName:null,surname:"Takada",slug:"akikazu-takada",fullName:"Akikazu Takada"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter"}},{type:"book",id:"8558",title:"Aerodynamics",subtitle:null,isOpenForSubmission:!1,hash:"db7263fc198dfb539073ba0260a7f1aa",slug:"aerodynamics",bookSignature:"Mofid Gorji-Bandpy and Aly-Mousaad Aly",coverURL:"https://cdn.intechopen.com/books/images_new/8558.jpg",editors:[{id:"35542",title:"Prof.",name:"Mofid",middleName:null,surname:"Gorji-Bandpy",slug:"mofid-gorji-bandpy",fullName:"Mofid Gorji-Bandpy"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter"}},{type:"book",id:"9668",title:"Chemistry and Biochemistry of Winemaking, Wine Stabilization and Aging",subtitle:null,isOpenForSubmission:!1,hash:"c5484276a314628acf21ec1bdc3a86b9",slug:"chemistry-and-biochemistry-of-winemaking-wine-stabilization-and-aging",bookSignature:"Fernanda Cosme, Fernando M. Nunes and Luís Filipe-Ribeiro",coverURL:"https://cdn.intechopen.com/books/images_new/9668.jpg",editors:[{id:"186819",title:"Prof.",name:"Fernanda",middleName:null,surname:"Cosme",slug:"fernanda-cosme",fullName:"Fernanda Cosme"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter"}},{type:"book",id:"7847",title:"Medical Toxicology",subtitle:null,isOpenForSubmission:!1,hash:"db9b65bea093de17a0855a1b27046247",slug:"medical-toxicology",bookSignature:"Pınar Erkekoglu and Tomohisa Ogawa",coverURL:"https://cdn.intechopen.com/books/images_new/7847.jpg",editors:[{id:"109978",title:"Prof.",name:"Pınar",middleName:null,surname:"Erkekoglu",slug:"pinar-erkekoglu",fullName:"Pınar Erkekoglu"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter"}},{type:"book",id:"8620",title:"Mining Techniques",subtitle:"Past, Present and Future",isOpenForSubmission:!1,hash:"b65658f81d14e9e57e49377869d3a575",slug:"mining-techniques-past-present-and-future",bookSignature:"Abhay Soni",coverURL:"https://cdn.intechopen.com/books/images_new/8620.jpg",editors:[{id:"271093",title:"Dr.",name:"Abhay",middleName:null,surname:"Soni",slug:"abhay-soni",fullName:"Abhay Soni"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter"}},{type:"book",id:"9660",title:"Inland Waters",subtitle:"Dynamics and Ecology",isOpenForSubmission:!1,hash:"975c26819ceb11a926793bc2adc62bd6",slug:"inland-waters-dynamics-and-ecology",bookSignature:"Adam Devlin, Jiayi Pan and Mohammad Manjur Shah",coverURL:"https://cdn.intechopen.com/books/images_new/9660.jpg",editors:[{id:"280757",title:"Dr.",name:"Adam",middleName:"Thomas",surname:"Devlin",slug:"adam-devlin",fullName:"Adam Devlin"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter"}},{type:"book",id:"9122",title:"Cosmetic Surgery",subtitle:null,isOpenForSubmission:!1,hash:"207026ca4a4125e17038e770d00ee152",slug:"cosmetic-surgery",bookSignature:"Yueh-Bih Tang",coverURL:"https://cdn.intechopen.com/books/images_new/9122.jpg",editors:[{id:"202122",title:"Prof.",name:"Yueh-Bih",middleName:null,surname:"Tang",slug:"yueh-bih-tang",fullName:"Yueh-Bih Tang"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter"}},{type:"book",id:"9043",title:"Parenting",subtitle:"Studies by an Ecocultural and Transactional Perspective",isOpenForSubmission:!1,hash:"6d21066c7438e459e4c6fb13217a5c8c",slug:"parenting-studies-by-an-ecocultural-and-transactional-perspective",bookSignature:"Loredana Benedetto and Massimo Ingrassia",coverURL:"https://cdn.intechopen.com/books/images_new/9043.jpg",editors:[{id:"193200",title:"Prof.",name:"Loredana",middleName:null,surname:"Benedetto",slug:"loredana-benedetto",fullName:"Loredana Benedetto"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter"}},{type:"book",id:"9731",title:"Oxidoreductase",subtitle:null,isOpenForSubmission:!1,hash:"852e6f862c85fc3adecdbaf822e64e6e",slug:"oxidoreductase",bookSignature:"Mahmoud Ahmed Mansour",coverURL:"https://cdn.intechopen.com/books/images_new/9731.jpg",editors:[{id:"224662",title:"Prof.",name:"Mahmoud Ahmed",middleName:null,surname:"Mansour",slug:"mahmoud-ahmed-mansour",fullName:"Mahmoud Ahmed Mansour"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter"}}],offset:12,limit:12,total:5229},hotBookTopics:{hotBooks:[],offset:0,limit:12,total:null},publish:{},publishingProposal:{success:null,errors:{}},books:{featuredBooks:[{type:"book",id:"9385",title:"Renewable Energy",subtitle:"Technologies and Applications",isOpenForSubmission:!1,hash:"a6b446d19166f17f313008e6c056f3d8",slug:"renewable-energy-technologies-and-applications",bookSignature:"Tolga Taner, Archana Tiwari and Taha Selim Ustun",coverURL:"https://cdn.intechopen.com/books/images_new/9385.jpg",editors:[{id:"197240",title:"Associate Prof.",name:"Tolga",middleName:null,surname:"Taner",slug:"tolga-taner",fullName:"Tolga Taner"}],equalEditorOne:{id:"186791",title:"Dr.",name:"Archana",middleName:null,surname:"Tiwari",slug:"archana-tiwari",fullName:"Archana Tiwari",profilePictureURL:"https://mts.intechopen.com/storage/users/186791/images/system/186791.jpg",biography:"Dr. Archana Tiwari is Associate Professor at Amity University, India. Her research interests include renewable sources of energy from microalgae and further utilizing the residual biomass for the generation of value-added products, bioremediation through microalgae and microbial consortium, antioxidative enzymes and stress, and nutraceuticals from microalgae. She has been working on algal biotechnology for the last two decades. She has published her research in many international journals and has authored many books and chapters with renowned publishing houses. She has also delivered talks as an invited speaker at many national and international conferences. Dr. Tiwari is the recipient of several awards including Researcher of the Year and Distinguished Scientist.",institutionString:"Amity University",position:null,outsideEditionCount:0,totalCites:0,totalAuthoredChapters:"3",totalChapterViews:"0",totalEditedBooks:"1",institution:{name:"Amity University",institutionURL:null,country:{name:"India"}}},equalEditorTwo:{id:"197609",title:"Prof.",name:"Taha Selim",middleName:null,surname:"Ustun",slug:"taha-selim-ustun",fullName:"Taha Selim Ustun",profilePictureURL:"https://mts.intechopen.com/storage/users/197609/images/system/197609.jpeg",biography:"Dr. Taha Selim Ustun received a Ph.D. in Electrical Engineering from Victoria University, Melbourne, Australia. He is a researcher with the Fukushima Renewable Energy Institute, AIST (FREA), where he leads the Smart Grid Cybersecurity Laboratory. Prior to that, he was a faculty member with the School of Electrical and Computer Engineering, Carnegie Mellon University, Pittsburgh, PA, USA. His current research interests include power systems protection, communication in power networks, distributed generation, microgrids, electric vehicle integration, and cybersecurity in smart grids. He serves on the editorial boards of IEEE Access, IEEE Transactions on Industrial Informatics, Energies, Electronics, Electricity, World Electric Vehicle and Information journals. Dr. Ustun is a member of the IEEE 2004 and 2800, IEC Renewable Energy Management WG 8, and IEC TC 57 WG17. He has been invited to run specialist courses in Africa, India, and China. He has delivered talks for the Qatar Foundation, the World Energy Council, the Waterloo Global Science Initiative, and the European Union Energy Initiative (EUEI). His research has attracted funding from prestigious programs in Japan, Australia, the European Union, and North America.",institutionString:"Fukushima Renewable Energy Institute, AIST (FREA)",position:null,outsideEditionCount:0,totalCites:0,totalAuthoredChapters:"1",totalChapterViews:"0",totalEditedBooks:"0",institution:{name:"National Institute of Advanced Industrial Science and Technology",institutionURL:null,country:{name:"Japan"}}},equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter"}},{type:"book",id:"8985",title:"Natural Resources Management and Biological Sciences",subtitle:null,isOpenForSubmission:!1,hash:"5c2e219a6c021a40b5a20c041dea88c4",slug:"natural-resources-management-and-biological-sciences",bookSignature:"Edward R. Rhodes and Humood Naser",coverURL:"https://cdn.intechopen.com/books/images_new/8985.jpg",editors:[{id:"280886",title:"Prof.",name:"Edward R",middleName:null,surname:"Rhodes",slug:"edward-r-rhodes",fullName:"Edward R Rhodes"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter"}},{type:"book",id:"9027",title:"Human Blood Group Systems and Haemoglobinopathies",subtitle:null,isOpenForSubmission:!1,hash:"d00d8e40b11cfb2547d1122866531c7e",slug:"human-blood-group-systems-and-haemoglobinopathies",bookSignature:"Osaro Erhabor and Anjana Munshi",coverURL:"https://cdn.intechopen.com/books/images_new/9027.jpg",editors:[{id:"35140",title:null,name:"Osaro",middleName:null,surname:"Erhabor",slug:"osaro-erhabor",fullName:"Osaro Erhabor"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter"}},{type:"book",id:"7841",title:"New Insights Into Metabolic Syndrome",subtitle:null,isOpenForSubmission:!1,hash:"ef5accfac9772b9e2c9eff884f085510",slug:"new-insights-into-metabolic-syndrome",bookSignature:"Akikazu Takada",coverURL:"https://cdn.intechopen.com/books/images_new/7841.jpg",editors:[{id:"248459",title:"Dr.",name:"Akikazu",middleName:null,surname:"Takada",slug:"akikazu-takada",fullName:"Akikazu Takada"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter"}},{type:"book",id:"8558",title:"Aerodynamics",subtitle:null,isOpenForSubmission:!1,hash:"db7263fc198dfb539073ba0260a7f1aa",slug:"aerodynamics",bookSignature:"Mofid Gorji-Bandpy and Aly-Mousaad Aly",coverURL:"https://cdn.intechopen.com/books/images_new/8558.jpg",editors:[{id:"35542",title:"Prof.",name:"Mofid",middleName:null,surname:"Gorji-Bandpy",slug:"mofid-gorji-bandpy",fullName:"Mofid Gorji-Bandpy"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter"}},{type:"book",id:"9668",title:"Chemistry and Biochemistry of Winemaking, Wine Stabilization and Aging",subtitle:null,isOpenForSubmission:!1,hash:"c5484276a314628acf21ec1bdc3a86b9",slug:"chemistry-and-biochemistry-of-winemaking-wine-stabilization-and-aging",bookSignature:"Fernanda Cosme, Fernando M. Nunes and Luís Filipe-Ribeiro",coverURL:"https://cdn.intechopen.com/books/images_new/9668.jpg",editors:[{id:"186819",title:"Prof.",name:"Fernanda",middleName:null,surname:"Cosme",slug:"fernanda-cosme",fullName:"Fernanda Cosme"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter"}},{type:"book",id:"7847",title:"Medical Toxicology",subtitle:null,isOpenForSubmission:!1,hash:"db9b65bea093de17a0855a1b27046247",slug:"medical-toxicology",bookSignature:"Pınar Erkekoglu and Tomohisa Ogawa",coverURL:"https://cdn.intechopen.com/books/images_new/7847.jpg",editors:[{id:"109978",title:"Prof.",name:"Pınar",middleName:null,surname:"Erkekoglu",slug:"pinar-erkekoglu",fullName:"Pınar Erkekoglu"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter"}},{type:"book",id:"8620",title:"Mining Techniques",subtitle:"Past, Present and Future",isOpenForSubmission:!1,hash:"b65658f81d14e9e57e49377869d3a575",slug:"mining-techniques-past-present-and-future",bookSignature:"Abhay Soni",coverURL:"https://cdn.intechopen.com/books/images_new/8620.jpg",editors:[{id:"271093",title:"Dr.",name:"Abhay",middleName:null,surname:"Soni",slug:"abhay-soni",fullName:"Abhay Soni"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter"}},{type:"book",id:"9660",title:"Inland Waters",subtitle:"Dynamics and Ecology",isOpenForSubmission:!1,hash:"975c26819ceb11a926793bc2adc62bd6",slug:"inland-waters-dynamics-and-ecology",bookSignature:"Adam Devlin, Jiayi Pan and Mohammad Manjur Shah",coverURL:"https://cdn.intechopen.com/books/images_new/9660.jpg",editors:[{id:"280757",title:"Dr.",name:"Adam",middleName:"Thomas",surname:"Devlin",slug:"adam-devlin",fullName:"Adam Devlin"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter"}},{type:"book",id:"9122",title:"Cosmetic Surgery",subtitle:null,isOpenForSubmission:!1,hash:"207026ca4a4125e17038e770d00ee152",slug:"cosmetic-surgery",bookSignature:"Yueh-Bih Tang",coverURL:"https://cdn.intechopen.com/books/images_new/9122.jpg",editors:[{id:"202122",title:"Prof.",name:"Yueh-Bih",middleName:null,surname:"Tang",slug:"yueh-bih-tang",fullName:"Yueh-Bih Tang"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter"}}],latestBooks:[{type:"book",id:"9550",title:"Entrepreneurship",subtitle:"Contemporary Issues",isOpenForSubmission:!1,hash:"9b4ac1ee5b743abf6f88495452b1e5e7",slug:"entrepreneurship-contemporary-issues",bookSignature:"Mladen Turuk",coverURL:"https://cdn.intechopen.com/books/images_new/9550.jpg",editedByType:"Edited by",editors:[{id:"319755",title:"Prof.",name:"Mladen",middleName:null,surname:"Turuk",slug:"mladen-turuk",fullName:"Mladen Turuk"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"10065",title:"Wavelet Theory",subtitle:null,isOpenForSubmission:!1,hash:"d8868e332169597ba2182d9b004d60de",slug:"wavelet-theory",bookSignature:"Somayeh Mohammady",coverURL:"https://cdn.intechopen.com/books/images_new/10065.jpg",editedByType:"Edited by",editors:[{id:"109280",title:"Dr.",name:"Somayeh",middleName:null,surname:"Mohammady",slug:"somayeh-mohammady",fullName:"Somayeh Mohammady"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"9313",title:"Clay Science and Technology",subtitle:null,isOpenForSubmission:!1,hash:"6fa7e70396ff10620e032bb6cfa6fb72",slug:"clay-science-and-technology",bookSignature:"Gustavo Morari Do Nascimento",coverURL:"https://cdn.intechopen.com/books/images_new/9313.jpg",editedByType:"Edited by",editors:[{id:"7153",title:"Prof.",name:"Gustavo",middleName:null,surname:"Morari Do Nascimento",slug:"gustavo-morari-do-nascimento",fullName:"Gustavo Morari Do Nascimento"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"9888",title:"Nuclear Power Plants",subtitle:"The Processes from the Cradle to the Grave",isOpenForSubmission:!1,hash:"c2c8773e586f62155ab8221ebb72a849",slug:"nuclear-power-plants-the-processes-from-the-cradle-to-the-grave",bookSignature:"Nasser Awwad",coverURL:"https://cdn.intechopen.com/books/images_new/9888.jpg",editedByType:"Edited by",editors:[{id:"145209",title:"Prof.",name:"Nasser",middleName:"S",surname:"Awwad",slug:"nasser-awwad",fullName:"Nasser Awwad"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"8098",title:"Resources of Water",subtitle:null,isOpenForSubmission:!1,hash:"d251652996624d932ef7b8ed62cf7cfc",slug:"resources-of-water",bookSignature:"Prathna Thanjavur Chandrasekaran, Muhammad Salik Javaid, Aftab Sadiq",coverURL:"https://cdn.intechopen.com/books/images_new/8098.jpg",editedByType:"Edited by",editors:[{id:"167917",title:"Dr.",name:"Prathna",middleName:null,surname:"Thanjavur Chandrasekaran",slug:"prathna-thanjavur-chandrasekaran",fullName:"Prathna Thanjavur Chandrasekaran"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"9644",title:"Glaciers and the Polar Environment",subtitle:null,isOpenForSubmission:!1,hash:"e8cfdc161794e3753ced54e6ff30873b",slug:"glaciers-and-the-polar-environment",bookSignature:"Masaki Kanao, Danilo Godone and Niccolò Dematteis",coverURL:"https://cdn.intechopen.com/books/images_new/9644.jpg",editedByType:"Edited by",editors:[{id:"51959",title:"Dr.",name:"Masaki",middleName:null,surname:"Kanao",slug:"masaki-kanao",fullName:"Masaki Kanao"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"10432",title:"Casting Processes and Modelling of Metallic Materials",subtitle:null,isOpenForSubmission:!1,hash:"2c5c9df938666bf5d1797727db203a6d",slug:"casting-processes-and-modelling-of-metallic-materials",bookSignature:"Zakaria Abdallah and Nada Aldoumani",coverURL:"https://cdn.intechopen.com/books/images_new/10432.jpg",editedByType:"Edited by",editors:[{id:"201670",title:"Dr.",name:"Zak",middleName:null,surname:"Abdallah",slug:"zak-abdallah",fullName:"Zak Abdallah"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"9671",title:"Macrophages",subtitle:null,isOpenForSubmission:!1,hash:"03b00fdc5f24b71d1ecdfd75076bfde6",slug:"macrophages",bookSignature:"Hridayesh Prakash",coverURL:"https://cdn.intechopen.com/books/images_new/9671.jpg",editedByType:"Edited by",editors:[{id:"287184",title:"Dr.",name:"Hridayesh",middleName:null,surname:"Prakash",slug:"hridayesh-prakash",fullName:"Hridayesh Prakash"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"8415",title:"Extremophilic Microbes and Metabolites",subtitle:"Diversity, Bioprospecting and Biotechnological Applications",isOpenForSubmission:!1,hash:"93e0321bc93b89ff73730157738f8f97",slug:"extremophilic-microbes-and-metabolites-diversity-bioprospecting-and-biotechnological-applications",bookSignature:"Afef Najjari, Ameur Cherif, Haïtham Sghaier and Hadda Imene Ouzari",coverURL:"https://cdn.intechopen.com/books/images_new/8415.jpg",editedByType:"Edited by",editors:[{id:"196823",title:"Dr.",name:"Afef",middleName:null,surname:"Najjari",slug:"afef-najjari",fullName:"Afef Najjari"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"9731",title:"Oxidoreductase",subtitle:null,isOpenForSubmission:!1,hash:"852e6f862c85fc3adecdbaf822e64e6e",slug:"oxidoreductase",bookSignature:"Mahmoud Ahmed Mansour",coverURL:"https://cdn.intechopen.com/books/images_new/9731.jpg",editedByType:"Edited by",editors:[{id:"224662",title:"Prof.",name:"Mahmoud Ahmed",middleName:null,surname:"Mansour",slug:"mahmoud-ahmed-mansour",fullName:"Mahmoud Ahmed Mansour"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}}]},subject:{topic:{id:"981",title:"Molecular Biology",slug:"biomedicine-molecular-biology",parent:{title:"Biomedicine",slug:"biomedicine"},numberOfBooks:4,numberOfAuthorsAndEditors:86,numberOfWosCitations:37,numberOfCrossrefCitations:26,numberOfDimensionsCitations:59,videoUrl:null,fallbackUrl:null,description:null},booksByTopicFilter:{topicSlug:"biomedicine-molecular-biology",sort:"-publishedDate",limit:12,offset:0},booksByTopicCollection:[{type:"book",id:"9569",title:"Methods in Molecular Medicine",subtitle:null,isOpenForSubmission:!1,hash:"691d3f3c4ac25a8093414e9b270d2843",slug:"methods-in-molecular-medicine",bookSignature:"Yusuf Tutar",coverURL:"https://cdn.intechopen.com/books/images_new/9569.jpg",editedByType:"Edited by",editors:[{id:"158492",title:"Prof.",name:"Yusuf",middleName:null,surname:"Tutar",slug:"yusuf-tutar",fullName:"Yusuf Tutar"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"8400",title:"Molecular Medicine",subtitle:null,isOpenForSubmission:!1,hash:"222a58353415f02de25a212213bddc00",slug:"molecular-medicine",bookSignature:"Sinem Nalbantoglu and Hakima Amri",coverURL:"https://cdn.intechopen.com/books/images_new/8400.jpg",editedByType:"Edited by",editors:[{id:"147712",title:"Dr.",name:"Sinem",middleName:null,surname:"Nalbantoglu",slug:"sinem-nalbantoglu",fullName:"Sinem Nalbantoglu"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"6205",title:"Polypeptide",subtitle:"New Insight into Drug Discovery and Development",isOpenForSubmission:!1,hash:"35f89bf9a197198efc0d44d0ff56d800",slug:"polypeptide-new-insight-into-drug-discovery-and-development",bookSignature:"Usman Sumo Friend Tambunan",coverURL:"https://cdn.intechopen.com/books/images_new/6205.jpg",editedByType:"Edited by",editors:[{id:"70235",title:"Prof.",name:"Usman Sumo Friend",middleName:null,surname:"Tambunan",slug:"usman-sumo-friend-tambunan",fullName:"Usman Sumo Friend Tambunan"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}},{type:"book",id:"1899",title:"Biomedicine",subtitle:null,isOpenForSubmission:!1,hash:"4b1bc3c9f0a151f8c9ffda110c2053d6",slug:"biomedicine",bookSignature:"Chao Lin",coverURL:"https://cdn.intechopen.com/books/images_new/1899.jpg",editedByType:"Edited by",editors:[{id:"109875",title:"Dr.",name:"Chao",middleName:null,surname:"Lin",slug:"chao-lin",fullName:"Chao Lin"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter",authoredCaption:"Edited by"}}],booksByTopicTotal:4,mostCitedChapters:[{id:"33119",doi:"10.5772/38349",title:"Additive Manufacturing Solutions for Improved Medical Implants",slug:"additive-manufacturing-solutions-for-improved-implants",totalDownloads:7517,totalCrossrefCites:7,totalDimensionsCites:21,book:{slug:"biomedicine",title:"Biomedicine",fullTitle:"Biomedicine"},signatures:"Vojislav Petrovic, Juan Vicente Haro, Jose Ramón Blasco and Luis Portolés",authors:[{id:"116774",title:"Dr.",name:"Vojislav",middleName:null,surname:"Petrovic",slug:"vojislav-petrovic",fullName:"Vojislav Petrovic"},{id:"116777",title:"MSc.",name:"Juan",middleName:"Vicente",surname:"Haro González",slug:"juan-haro-gonzalez",fullName:"Juan Haro González"},{id:"116778",title:"BSc.",name:"José Ramón",middleName:null,surname:"Blasco Puchades",slug:"jose-ramon-blasco-puchades",fullName:"José Ramón Blasco Puchades"},{id:"116779",title:"BSc.",name:"Luís",middleName:null,surname:"Portolés Griñán",slug:"luis-portoles-grinan",fullName:"Luís Portolés Griñán"}]},{id:"33113",doi:"10.5772/33951",title:"Encapsulation and Surface Engineering of Pancreatic Islets: Advances and Challenges",slug:"encapsulation-and-surface-engineering-of-pancreatic-islets-advances-and-challenges-",totalDownloads:3020,totalCrossrefCites:4,totalDimensionsCites:8,book:{slug:"biomedicine",title:"Biomedicine",fullTitle:"Biomedicine"},signatures:"Veronika Kozlovskaya, Oleksandra Zavgorodnya and Eugenia Kharlampieva",authors:[{id:"97932",title:"Prof.",name:"Eugenia",middleName:null,surname:"Kharlampieva",slug:"eugenia-kharlampieva",fullName:"Eugenia Kharlampieva"},{id:"101333",title:"Dr.",name:"Veronika",middleName:null,surname:"Kozlovskaya",slug:"veronika-kozlovskaya",fullName:"Veronika Kozlovskaya"},{id:"135852",title:"MSc.",name:"Oleksandra",middleName:null,surname:"Zavgorodnya",slug:"oleksandra-zavgorodnya",fullName:"Oleksandra Zavgorodnya"}]},{id:"68486",doi:"10.5772/intechopen.88563",title:"Metabolomics: Basic Principles and Strategies",slug:"metabolomics-basic-principles-and-strategies",totalDownloads:1557,totalCrossrefCites:7,totalDimensionsCites:8,book:{slug:"molecular-medicine",title:"Molecular Medicine",fullTitle:"Molecular Medicine"},signatures:"Sinem Nalbantoglu",authors:[{id:"147712",title:"Dr.",name:"Sinem",middleName:null,surname:"Nalbantoglu",slug:"sinem-nalbantoglu",fullName:"Sinem Nalbantoglu"}]}],mostDownloadedChaptersLast30Days:[{id:"68486",title:"Metabolomics: Basic Principles and Strategies",slug:"metabolomics-basic-principles-and-strategies",totalDownloads:1567,totalCrossrefCites:7,totalDimensionsCites:9,book:{slug:"molecular-medicine",title:"Molecular Medicine",fullTitle:"Molecular Medicine"},signatures:"Sinem Nalbantoglu",authors:[{id:"147712",title:"Dr.",name:"Sinem",middleName:null,surname:"Nalbantoglu",slug:"sinem-nalbantoglu",fullName:"Sinem Nalbantoglu"}]},{id:"74620",title:"Molecular Mechanisms of Distinct Diseases",slug:"molecular-mechanisms-of-distinct-diseases",totalDownloads:180,totalCrossrefCites:0,totalDimensionsCites:0,book:{slug:"methods-in-molecular-medicine",title:"Methods in Molecular Medicine",fullTitle:"Methods in Molecular Medicine"},signatures:"Adnan Batman, İrem Yalim Camci, Elif Kadioglu, Kezban Uçar Çifçi, Berçem Yeman Kıyak, Servet Tunoglu, Ezgi Nurdan Yenilmez Tunoglu and Yusuf Tutar",authors:[{id:"158492",title:"Prof.",name:"Yusuf",middleName:null,surname:"Tutar",slug:"yusuf-tutar",fullName:"Yusuf Tutar"},{id:"342753",title:"Dr.",name:"Servet",middleName:null,surname:"Tunoglu",slug:"servet-tunoglu",fullName:"Servet Tunoglu"},{id:"342756",title:"Dr.",name:"Adnan",middleName:null,surname:"Batman",slug:"adnan-batman",fullName:"Adnan Batman"},{id:"343421",title:null,name:"Berçem",middleName:null,surname:"Yeman Kıyak",slug:"bercem-yeman-kiyak",fullName:"Berçem Yeman Kıyak"},{id:"344587",title:null,name:"İrem Yalım",middleName:null,surname:"Camcı",slug:"irem-yalim-camci",fullName:"İrem Yalım Camcı"},{id:"344588",title:"Ph.D. Student",name:"Elif",middleName:null,surname:"Kadıoglu",slug:"elif-kadioglu",fullName:"Elif Kadıoglu"},{id:"344589",title:null,name:"Kezban",middleName:null,surname:"Uçar Çiftçi",slug:"kezban-ucar-ciftci",fullName:"Kezban Uçar Çiftçi"}]},{id:"67272",title:"Introductory Chapter: Insight into the OMICS Technologies and Molecular Medicine",slug:"introductory-chapter-insight-into-the-omics-technologies-and-molecular-medicine",totalDownloads:481,totalCrossrefCites:0,totalDimensionsCites:0,book:{slug:"molecular-medicine",title:"Molecular Medicine",fullTitle:"Molecular Medicine"},signatures:"Sinem Nalbantoglu and Abdullah Karadag",authors:[{id:"147712",title:"Dr.",name:"Sinem",middleName:null,surname:"Nalbantoglu",slug:"sinem-nalbantoglu",fullName:"Sinem Nalbantoglu"}]},{id:"72595",title:"Integrating Evolutionary Genetics to Medical Genomics: Evolutionary Approaches to Investigate Disease-Causing Variants",slug:"integrating-evolutionary-genetics-to-medical-genomics-evolutionary-approaches-to-investigate-disease",totalDownloads:171,totalCrossrefCites:0,totalDimensionsCites:0,book:{slug:"methods-in-molecular-medicine",title:"Methods in Molecular Medicine",fullTitle:"Methods in Molecular Medicine"},signatures:"Ugur Sezerman, Tugce Bozkurt and Fatma Sadife Isleyen",authors:[{id:"283538",title:"Prof.",name:"Osman Ugur",middleName:null,surname:"Sezerman",slug:"osman-ugur-sezerman",fullName:"Osman Ugur Sezerman"},{id:"318900",title:"Ph.D. Student",name:"Fatma Sadife",middleName:null,surname:"Isleyen",slug:"fatma-sadife-isleyen",fullName:"Fatma Sadife Isleyen"},{id:"319214",title:"M.Sc.",name:"Tugce",middleName:null,surname:"Bozkurt",slug:"tugce-bozkurt",fullName:"Tugce Bozkurt"}]},{id:"33120",title:"Crossings on Public Perception of Biomedicine: Spain and the European Indicators",slug:"crossings-on-public-perception-of-biomedicine-spain-and-the-european-indicators-",totalDownloads:1359,totalCrossrefCites:0,totalDimensionsCites:0,book:{slug:"biomedicine",title:"Biomedicine",fullTitle:"Biomedicine"},signatures:"Eulalia Pérez Sedeño and María José Miranda Suárez",authors:[{id:"112999",title:"Dr.",name:"Eulalia",middleName:null,surname:"Perez Sedeno",slug:"eulalia-perez-sedeno",fullName:"Eulalia Perez Sedeno"},{id:"117289",title:"MSc.",name:"María José",middleName:null,surname:"Miranda Suárez",slug:"maria-jose-miranda-suarez",fullName:"María José Miranda Suárez"}]},{id:"64435",title:"Precision Medicine: Role of Biomarkers in Early Prediction and Diagnosis of Alzheimer’s Disease",slug:"precision-medicine-role-of-biomarkers-in-early-prediction-and-diagnosis-of-alzheimer-s-disease",totalDownloads:1246,totalCrossrefCites:1,totalDimensionsCites:2,book:{slug:"molecular-medicine",title:"Molecular Medicine",fullTitle:"Molecular Medicine"},signatures:"Liming Shen, Sijian Xia, Huajie Zhang, Fang Yao, Xukun Liu, Yuxi Zhao, Ming Ying, Javed Iqbal and Qiong Liu",authors:[{id:"263601",title:"Associate Prof.",name:"Liming",middleName:null,surname:"Shen",slug:"liming-shen",fullName:"Liming Shen"},{id:"273558",title:"Ms.",name:"Sijian",middleName:null,surname:"Xia",slug:"sijian-xia",fullName:"Sijian Xia"},{id:"273562",title:"Ms.",name:"Huajie",middleName:null,surname:"Zhang",slug:"huajie-zhang",fullName:"Huajie Zhang"},{id:"273563",title:"Ms.",name:"Fang",middleName:null,surname:"Yao",slug:"fang-yao",fullName:"Fang Yao"},{id:"273564",title:"Ms.",name:"Xukun",middleName:null,surname:"Liu",slug:"xukun-liu",fullName:"Xukun Liu"},{id:"273565",title:"Ms.",name:"Yuxi",middleName:null,surname:"Zhao",slug:"yuxi-zhao",fullName:"Yuxi Zhao"},{id:"273566",title:"Dr.",name:"Ming",middleName:null,surname:"Ying",slug:"ming-ying",fullName:"Ming Ying"},{id:"273567",title:"Dr.",name:"Javed",middleName:null,surname:"Iqbal",slug:"javed-iqbal",fullName:"Javed Iqbal"},{id:"273569",title:"Prof.",name:"Qiong",middleName:null,surname:"Liu",slug:"qiong-liu",fullName:"Qiong Liu"}]},{id:"71705",title:"Landscape Genetics: From Classic Molecular Markers to Genomics",slug:"landscape-genetics-from-classic-molecular-markers-to-genomics",totalDownloads:226,totalCrossrefCites:0,totalDimensionsCites:1,book:{slug:"methods-in-molecular-medicine",title:"Methods in Molecular Medicine",fullTitle:"Methods in Molecular Medicine"},signatures:"Enéas Ricardo Konzen and Maria Imaculada Zucchi",authors:[{id:"196963",title:"Dr.",name:"Enéas Ricardo",middleName:null,surname:"Konzen",slug:"eneas-ricardo-konzen",fullName:"Enéas Ricardo Konzen"},{id:"305359",title:"Prof.",name:"Maria Imaculada",middleName:null,surname:"Zucchi",slug:"maria-imaculada-zucchi",fullName:"Maria Imaculada Zucchi"}]},{id:"33113",title:"Encapsulation and Surface Engineering of Pancreatic Islets: Advances and Challenges",slug:"encapsulation-and-surface-engineering-of-pancreatic-islets-advances-and-challenges-",totalDownloads:3024,totalCrossrefCites:4,totalDimensionsCites:8,book:{slug:"biomedicine",title:"Biomedicine",fullTitle:"Biomedicine"},signatures:"Veronika Kozlovskaya, Oleksandra Zavgorodnya and Eugenia Kharlampieva",authors:[{id:"97932",title:"Prof.",name:"Eugenia",middleName:null,surname:"Kharlampieva",slug:"eugenia-kharlampieva",fullName:"Eugenia Kharlampieva"},{id:"101333",title:"Dr.",name:"Veronika",middleName:null,surname:"Kozlovskaya",slug:"veronika-kozlovskaya",fullName:"Veronika Kozlovskaya"},{id:"135852",title:"MSc.",name:"Oleksandra",middleName:null,surname:"Zavgorodnya",slug:"oleksandra-zavgorodnya",fullName:"Oleksandra Zavgorodnya"}]},{id:"61938",title:"Calcitonin-Related Polypeptide Alpha Gene Polymorphisms and Related Diseases",slug:"calcitonin-related-polypeptide-alpha-gene-polymorphisms-and-related-diseases",totalDownloads:457,totalCrossrefCites:0,totalDimensionsCites:0,book:{slug:"polypeptide-new-insight-into-drug-discovery-and-development",title:"Polypeptide",fullTitle:"Polypeptide - New Insight into Drug Discovery and Development"},signatures:"Nevra Alkanli, Arzu Ay and Suleyman Serdar Alkanli",authors:[{id:"234410",title:"Dr.",name:"Nevra",middleName:null,surname:"Alkanli",slug:"nevra-alkanli",fullName:"Nevra Alkanli"},{id:"235954",title:"Dr.",name:"Arzu",middleName:null,surname:"Ay",slug:"arzu-ay",fullName:"Arzu Ay"},{id:"235955",title:"Mr.",name:"Suleyman Serdar",middleName:null,surname:"Alkanli",slug:"suleyman-serdar-alkanli",fullName:"Suleyman Serdar Alkanli"}]},{id:"72817",title:"Clinical Validation of a Whole Exome Sequencing Pipeline",slug:"clinical-validation-of-a-whole-exome-sequencing-pipeline",totalDownloads:183,totalCrossrefCites:0,totalDimensionsCites:0,book:{slug:"methods-in-molecular-medicine",title:"Methods in Molecular Medicine",fullTitle:"Methods in Molecular Medicine"},signatures:"Debra O. Prosser, Indu Raja, Kelly Kolkiewicz, Antonio Milano and Donald Roy Love",authors:[{id:"319362",title:"Dr.",name:"Donald Roy",middleName:null,surname:"Love",slug:"donald-roy-love",fullName:"Donald Roy Love"},{id:"319363",title:"Mrs.",name:"Debra O.",middleName:null,surname:"Prosser",slug:"debra-o.-prosser",fullName:"Debra O. Prosser"}]}],onlineFirstChaptersFilter:{topicSlug:"biomedicine-molecular-biology",limit:3,offset:0},onlineFirstChaptersCollection:[],onlineFirstChaptersTotal:0},preDownload:{success:null,errors:{}},aboutIntechopen:{},privacyPolicy:{},peerReviewing:{},howOpenAccessPublishingWithIntechopenWorks:{},sponsorshipBooks:{sponsorshipBooks:[{type:"book",id:"10176",title:"Microgrids and Local Energy Systems",subtitle:null,isOpenForSubmission:!0,hash:"c32b4a5351a88f263074b0d0ca813a9c",slug:null,bookSignature:"Prof. Nick Jenkins",coverURL:"https://cdn.intechopen.com/books/images_new/10176.jpg",editedByType:null,editors:[{id:"55219",title:"Prof.",name:"Nick",middleName:null,surname:"Jenkins",slug:"nick-jenkins",fullName:"Nick Jenkins"}],equalEditorOne:null,equalEditorTwo:null,equalEditorThree:null,productType:{id:"1",chapterContentType:"chapter"}}],offset:8,limit:8,total:1},route:{name:"profile.detail",path:"/profiles/190158/anita-ioana-visan",hash:"",query:{},params:{id:"190158",slug:"anita-ioana-visan"},fullPath:"/profiles/190158/anita-ioana-visan",meta:{},from:{name:null,path:"/",hash:"",query:{},params:{},fullPath:"/",meta:{}}}},function(){var m;(m=document.currentScript||document.scripts[document.scripts.length-1]).parentNode.removeChild(m)}()