<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="brief-report">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Mater.</journal-id>
<journal-title>Frontiers in Materials</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Mater.</abbrev-journal-title>
<issn pub-type="epub">2296-8016</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmats.2019.00004</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Materials</subject>
<subj-group>
<subject>Perspective</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>The Role of Volume Fraction and Additives on the Rheology of Suspensions of Micron Sized Iron Particles</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name><surname>Bossis</surname> <given-names>Georges</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/593835/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Volkova</surname> <given-names>Olga</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Grasselli</surname> <given-names>Yan</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Ciffreo</surname> <given-names>Alain</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Laboratory Inphyni, Institute of Physics of Nice, CNRS, University of Nice Sophia-Antipolis</institution>, <addr-line>Nice</addr-line>, <country>France</country></aff>
<aff id="aff2"><sup>2</sup><institution>SKEMA Business School, Universit&#x000E9; C&#x000F4;te d&#x00027;Azur</institution>, <addr-line>Sophia Antipolis</addr-line>, <country>France</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Seung-Bok Choi, Inha University, South Korea</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Xufeng Dong, Dalian University of Technology (DUT), China; Evguenia Korobko, A.V. Luikov Heat and Mass Transfer Institute (NASB), Belarus</p></fn>
<corresp id="c001">&#x0002A;Correspondence: Georges Bossis <email>bossis&#x00040;unice.fr</email></corresp>
<fn fn-type="other" id="fn001"><p>This article was submitted to Smart Materials, a section of the journal Frontiers in Materials</p></fn></author-notes>
<pub-date pub-type="epub">
<day>08</day>
<month>02</month>
<year>2019</year>
</pub-date>
<pub-date pub-type="collection">
<year>2019</year>
</pub-date>
<volume>6</volume>
<elocation-id>4</elocation-id>
<history>
<date date-type="received">
<day>29</day>
<month>10</month>
<year>2018</year>
</date>
<date date-type="accepted">
<day>10</day>
<month>01</month>
<year>2019</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2019 Bossis, Volkova, Grasselli and Ciffreo.</copyright-statement>
<copyright-year>2019</copyright-year>
<copyright-holder>Bossis, Volkova, Grasselli and Ciffreo</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract><p>The increase of the yield stress vs. the magnetic field is the most important quantity characterizing the efficiency of a magnetorheological suspension. The theory based on the formation of columnar aggregates predicts a linear variation with the volume fraction of magnetic particles. In this paper we review previous models used to calculate forces and yield stress and will introduce a new model based on rupture at zero strain. Predictions of these models are compared with the experimental data obtained for carbonyl iron particles, by different authors. Whereas, previous analytical prediction strongly overestimates experimental yield stress, those calculated using the Finite Element Method (FEM), together with affine trajectories, reproduce the experiments well and show a linear dependence with the volume fraction and a H<sup>3/2</sup> behavior between 50 and 200 kA/m. Nevertheless, at very high-volume fractions (&#x0003E;55%), where the suspension can only flow in the presence of specific additives, the dependence of the yield stress vs. the volume fraction and the magnetic field is dramatically changed. We observed a jamming transition, which is triggered by the application of a low magnetic field and which depends strongly on the volume of the fraction. Here, we will discuss new perspectives arising from the use of these very high-volume fractions.</p></abstract>
<kwd-group>
<kwd>carbonyl iron</kwd>
<kwd>yield stress</kwd>
<kwd>discontinuous shear thickening</kwd>
<kwd>jamming</kwd>
<kwd>magnetorheology</kwd>
</kwd-group>
<contract-sponsor id="cn001">Centre National de la Recherche Scientifique<named-content content-type="fundref-id">10.13039/501100004794</named-content></contract-sponsor>
<counts>
<fig-count count="2"/>
<table-count count="0"/>
<equation-count count="7"/>
<ref-count count="38"/>
<page-count count="6"/>
<word-count count="5151"/>
</counts>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>Introduction</title>
<p>Magnetorheological (MR) suspensions are smart fluids made of magnetizable particles, with a typical size varying between 100 nm and a few microns dispersed in a carrier liquid, which can be mineral oil, silicone oil, ethylene glycol, etc. The interest of these fluids is their rapid and reversible transformation in a solid when they are submitted to a magnetic field. This transformation occurs due to the attractive interaction of a dipolar nature, generated by the magnetization of the particles. The solid phase is characterized by its resistance to a shear strain, namely the yield stress, over which the solid begins to creep. Most applications look for fluids that present the highest possible yield. Of course, the higher the magnetic field, the higher the yield stress, with the limit of the full magnetization of the particles as well as the inconvenience of weight and place that are necessary for the coils to produce a high magnetic field. Besides the importance of a high yield stress, other requisites are also important for industrial uses such as easy resuspension after sedimentation and a low off state viscosity. In order to obtain a high yield stress and a reversible magnetization, very few materials are available, and in practice industrial MR fluids are based on carbonyl iron particles as they present a high saturation magnetization and are synthesized in large quantities to produce sintered magnetic circuits. Iron particles have the inconvenience of high density and therefore of rapid sedimentation; furthermore, due to the pressure generated by their weight and the attractive Van der Waals forces, they easily aggregate irreversibly. In order to prevent this irreversible aggregation, coating of the particles by some molecules (surfactant, polymer, even ions) is necessary. Many papers and patents describe these formulations and their effects in reducing sedimentation and aggregation [for a recent review see reference (Ashtiani et al., <xref ref-type="bibr" rid="B1">2015</xref>)]. This coating, which prevents the particles to come into close contact can also present the inconvenience of lowering the magnetic force, which strongly depends on the distance between the particles. Ultimately however, despite decades of research, the maximum yield stress of commercial fluids remains in the 40&#x02013;60 kPa range for a 0.3&#x02013;0.5 Tesla magnetic field induction. We first take a look at the yield stresses predicted by the available models and compare these predictions with the experimental data, then, in light of recent experimental results involving frictional stresses, we will focus on the performance of these fluids and how they can be improved.</p>
</sec>
<sec id="s2">
<title>Prediction of the Yield Stress</title>
<p>The cohesion of the solid phase, and therefore its yield stress, depends first on the strength of the magnetic forces developed between the particles, and second on the structure and its evolution under combined magnetic and hydrodynamic forces. We briefly review the main models and then compare their prediction, both for the forces between the particles and for the yield stress. To do this analysis we take a look at the magnetic properties of the carbonyl iron particles, those currently used in practical applications. The magnetization curve is described by the Frolish Kennelly equation: M &#x0003D; M<sub>s</sub> &#x003C7;<sub>i</sub> H/(M<sub>s</sub> &#x0002B; &#x003C7;<sub>i</sub> H) Where M<sub>s</sub> is the saturation magnetization and &#x003C7;<sub>i</sub> the initial permeability. For the saturation magnetization we have taken the pure iron: M<sub>s</sub> &#x0003D; 1,680 kA/m and the initial permeability &#x003C7;<sub>i</sub> &#x0003D; 70 (De Vicente et al., <xref ref-type="bibr" rid="B12">2002</xref>).</p>
<sec>
<title>Magnetic Forces Between Two Particles</title>
<p>Contrary to well-known dipolar approximation, a multipolar approach for two particles, first presented by Klingenberg (Klingenberg and Zukoski, <xref ref-type="bibr" rid="B25">1990</xref>) and then for several particles described in subsequent works (Clercx and Bossis, <xref ref-type="bibr" rid="B10">1993</xref>), gives an exact prediction of the force for a linear magnetization. The result for two particles reads:</p>
<disp-formula id="E1"><label>(1)</label><mml:math id="M1"><mml:mtable columnalign='left'><mml:mtr><mml:mtd><mml:msub><mml:mtext>F</mml:mtext><mml:mtext>m</mml:mtext></mml:msub><mml:mo stretchy='false'>(</mml:mo><mml:mover><mml:mtext>r</mml:mtext><mml:mo>&#x02192;</mml:mo></mml:mover><mml:mo stretchy='false'>)</mml:mo><mml:mtext>&#x000A0;</mml:mtext><mml:mo>=</mml:mo><mml:mn>12</mml:mn><mml:mo>&#x003C0;</mml:mo><mml:msub><mml:mo>&#x003BC;</mml:mo><mml:mn>0</mml:mn></mml:msub><mml:msup><mml:mtext>a</mml:mtext><mml:mn>2</mml:mn></mml:msup><mml:msup><mml:mi>&#x003B2;</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:msup><mml:mtext>H</mml:mtext><mml:mn>2</mml:mn></mml:msup><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mtext>a</mml:mtext><mml:mtext>r</mml:mtext></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mn>4</mml:mn></mml:msup><mml:mtext>&#x000A0;</mml:mtext><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mn>2</mml:mn><mml:msub><mml:mtext>f</mml:mtext><mml:mrow><mml:mo>/</mml:mo><mml:mo>/</mml:mo></mml:mrow></mml:msub><mml:msup><mml:mrow><mml:mtext>cos</mml:mtext></mml:mrow><mml:mn>2</mml:mn></mml:msup><mml:mo stretchy='false'>(</mml:mo><mml:mo>&#x003B8;</mml:mo><mml:mo stretchy='false'>)</mml:mo><mml:mo>-</mml:mo><mml:msub><mml:mtext>f</mml:mtext><mml:mo>&#x022A5;</mml:mo></mml:msub><mml:msup><mml:mrow><mml:mtext>sin</mml:mtext></mml:mrow><mml:mn>2</mml:mn></mml:msup><mml:mo stretchy='false'>(</mml:mo><mml:mo>&#x003B8;</mml:mo><mml:mo stretchy='false'>)</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mover><mml:mrow><mml:msub><mml:mtext>e</mml:mtext><mml:mtext>r</mml:mtext></mml:msub></mml:mrow><mml:mo>&#x02192;</mml:mo></mml:mover></mml:mrow></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtext>&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mrow><mml:mrow><mml:mo>+</mml:mo><mml:msub><mml:mtext>f</mml:mtext><mml:mo>&#x00393;</mml:mo></mml:msub><mml:mtext>sin</mml:mtext><mml:mo stretchy='false'>(</mml:mo><mml:mn>2</mml:mn><mml:mo>&#x003B8;</mml:mo><mml:mo stretchy='false'>)</mml:mo><mml:msub><mml:mrow><mml:mover><mml:mtext>e</mml:mtext><mml:mo>&#x02192;</mml:mo></mml:mover></mml:mrow><mml:mi>&#x003B8;</mml:mi></mml:msub></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mtext>&#x02003;with&#x000A0;</mml:mtext><mml:mi>&#x003B2;</mml:mi><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mo>&#x003BC;</mml:mo><mml:mtext>p</mml:mtext></mml:msub><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:msub><mml:mo>&#x003BC;</mml:mo><mml:mtext>p</mml:mtext></mml:msub><mml:mo>+</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:mfrac></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>In Equation (1) a is the radius of the particles, &#x003BC;<sub>p</sub> their relative permeability given by &#x003BC;<sub>p</sub> &#x0003D; 1&#x0002B;M/H and f<sub>//</sub>, f<sub>&#x022A5;</sub>, f<sub><italic>&#x00393;</italic></sub> are functions of the initial permeability: &#x003BC;<sub>i</sub> &#x0003D; &#x003C7;<sub>i</sub> &#x0002B;1 and of the normalized separation between the surfaces of the two particles &#x003B5; &#x0003D; r/a-2 (cf <xref ref-type="supplementary-material" rid="SM1">Figure S1</xref>). A <xref ref-type="supplementary-material" rid="SM5">Table S2</xref> of these functions for different separations and different values of &#x003C7;<sub>i</sub> is provided in the <xref ref-type="supplementary-material" rid="SM1">Supplementary Material</xref>.</p>
<p>A popular model derived in Ginder et al. (<xref ref-type="bibr" rid="B18">1996</xref>) uses, on the contrary, the approximation that the pole of the particles is completely saturated (cf <xref ref-type="supplementary-material" rid="SM1">Figure S1</xref>) in a zone defined by &#x003C1; &#x0003D; &#x003B4; where &#x003C1; is the polar coordinate. The interparticle force coming from the integral of the magnetic stress tensor is then given by:</p>
<disp-formula id="E2"><label>(2)</label><mml:math id="M2"><mml:mtable columnalign='left'><mml:mtr><mml:mtd><mml:msub><mml:mtext>F</mml:mtext><mml:mtext>r</mml:mtext></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>&#x003BC;</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow><mml:mn>2</mml:mn></mml:mfrac><mml:mstyle displaystyle='true'><mml:mrow><mml:munderover><mml:mo>&#x0222B;</mml:mo><mml:mn>0</mml:mn><mml:mtext>W</mml:mtext></mml:munderover><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mtext>H</mml:mtext><mml:mo stretchy='false'>(</mml:mo><mml:mi>&#x003C1;</mml:mi><mml:mo stretchy='false'>)</mml:mo><mml:mo>-</mml:mo><mml:mtext>H</mml:mtext></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mn>2</mml:mn></mml:msup><mml:mtext>&#x000A0;</mml:mtext><mml:mn>2</mml:mn><mml:mi>&#x003C0;</mml:mi><mml:mi>&#x003C1;</mml:mi><mml:mtext>d</mml:mtext><mml:mi>&#x003C1;</mml:mi></mml:mrow></mml:mrow></mml:mstyle><mml:mtext>&#x02003;with&#x000A0;H</mml:mtext><mml:mo stretchy='false'>(</mml:mo><mml:mi>&#x003C1;</mml:mi><mml:mo stretchy='false'>)</mml:mo><mml:mo>=</mml:mo><mml:mtext>H</mml:mtext><mml:mfrac><mml:mrow><mml:mo stretchy='false'>(</mml:mo><mml:mn>2</mml:mn><mml:mo>+</mml:mo><mml:mi>&#x003B5;</mml:mi><mml:mo stretchy='false'>)</mml:mo></mml:mrow><mml:mrow><mml:mi>&#x003B5;</mml:mi><mml:mo>+</mml:mo><mml:msup><mml:mi>&#x003C1;</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mfrac></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtext>&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;for&#x000A0;&#x000A0;</mml:mtext><mml:mi>&#x003C1;</mml:mi><mml:mo>&#x0003E;</mml:mo><mml:mi>&#x003B4;</mml:mi><mml:mtext>&#x000A0;&#x000A0;and&#x000A0;&#x000A0;H</mml:mtext><mml:mo stretchy='false'>(</mml:mo><mml:mi>&#x003C1;</mml:mi><mml:mo stretchy='false'>)</mml:mo><mml:mo>=</mml:mo><mml:msub><mml:mtext>M</mml:mtext><mml:mtext>s</mml:mtext></mml:msub><mml:mtext>&#x000A0;&#x000A0;for&#x000A0;&#x000A0;</mml:mtext><mml:mi>&#x003C1;</mml:mi><mml:mo>&#x0003C;</mml:mo><mml:mi>&#x003B4;</mml:mi></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>The value of &#x003B4; is obtained from the equality between the two expressions for &#x003C1; &#x0003D; &#x003B4; which gives: <inline-formula><mml:math id="M3"><mml:mi>&#x003B4;</mml:mi><mml:mo>=</mml:mo><mml:msqrt><mml:mrow><mml:mn>2</mml:mn><mml:mtext>H</mml:mtext><mml:mo>/</mml:mo><mml:msub><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mrow><mml:mtext>s</mml:mtext></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:mi>&#x003B5;</mml:mi></mml:mrow></mml:msqrt></mml:math></inline-formula></p>
<p>In these expressions the distances are normalized by the radius of the sphere and W is of order unity. In the limit where &#x003B5;/2 &#x0003C; <italic>H</italic>/<italic>M</italic><sub><italic>s</italic></sub> &#x0003C; &#x0003C; 1, the upper bound of the integral is negligible and we obtain:</p>
<disp-formula id="E3"><label>(3)</label><mml:math id="M4"><mml:msub><mml:mtext>F</mml:mtext><mml:mrow><mml:mtext>ra</mml:mtext></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mo>&#x003C0;</mml:mo><mml:msub><mml:mo>&#x003BC;</mml:mo><mml:mn>0</mml:mn></mml:msub><mml:msup><mml:mtext>a</mml:mtext><mml:mn>2</mml:mn></mml:msup><mml:msubsup><mml:mtext>M</mml:mtext><mml:mtext>s</mml:mtext><mml:mn>2</mml:mn></mml:msubsup><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mn>2</mml:mn><mml:mtext>H</mml:mtext></mml:mrow><mml:mrow><mml:msub><mml:mtext>M</mml:mtext><mml:mtext>s</mml:mtext></mml:msub></mml:mrow></mml:mfrac><mml:mo>-</mml:mo><mml:mfrac><mml:mi>&#x003B5;</mml:mi><mml:mn>2</mml:mn></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>
<p>It is possible to check the validity of this approximation, using a finite element method (FEM) for the calculation of the magnetic field. We have used the FEMM software to calculate the force between two spheres placed in a constant magnetic field H:</p>
<disp-formula id="E4"><label>(4)</label><mml:math id="M5"><mml:mtable columnalign='left'><mml:mtr><mml:mtd><mml:msub><mml:mtext>F</mml:mtext><mml:mrow><mml:mtext>rFEM</mml:mtext></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>2</mml:mn><mml:mi>&#x003C0;</mml:mi><mml:msup><mml:mtext>a</mml:mtext><mml:mn>2</mml:mn></mml:msup><mml:mstyle displaystyle='true'><mml:mrow><mml:munderover><mml:mo>&#x0222B;</mml:mo><mml:mn>0</mml:mn><mml:mi>&#x003C0;</mml:mi></mml:munderover><mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:msub><mml:mi>&#x003C3;</mml:mi><mml:mrow><mml:mtext>rr</mml:mtext></mml:mrow></mml:msub><mml:mtext>&#x000A0;cos</mml:mtext><mml:mo stretchy='false'>(</mml:mo><mml:mi>&#x003B8;</mml:mi><mml:mo stretchy='false'>)</mml:mo><mml:mo>-</mml:mo><mml:msub><mml:mi>&#x003C3;</mml:mi><mml:mrow><mml:mtext>r</mml:mtext><mml:mi>&#x003B8;</mml:mi></mml:mrow></mml:msub><mml:mtext>&#x000A0;sin</mml:mtext><mml:mo stretchy='false'>(</mml:mo><mml:mi>&#x003B8;</mml:mi><mml:mo stretchy='false'>)</mml:mo></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mrow></mml:mrow></mml:mstyle><mml:mtext>sin</mml:mtext><mml:mo stretchy='false'>(</mml:mo><mml:mi>&#x003B8;</mml:mi><mml:mo stretchy='false'>)</mml:mo><mml:mtext>d</mml:mtext><mml:mi>&#x003B8;</mml:mi></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtext>&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;with&#x000A0;</mml:mtext><mml:msub><mml:mi>&#x003C3;</mml:mi><mml:mrow><mml:mtext>rr</mml:mtext></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>&#x003BC;</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo stretchy='false'>(</mml:mo><mml:msub><mml:mi>&#x003BC;</mml:mi><mml:mtext>f</mml:mtext></mml:msub><mml:msubsup><mml:mtext>H</mml:mtext><mml:mtext>r</mml:mtext><mml:mn>2</mml:mn></mml:msubsup><mml:mo>-</mml:mo><mml:mn>0.5</mml:mn><mml:msup><mml:mtext>H</mml:mtext><mml:mn>2</mml:mn></mml:msup><mml:mo stretchy='false'>)</mml:mo><mml:mtext>&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:msub><mml:mi>&#x003C3;</mml:mi><mml:mrow><mml:mtext>r</mml:mtext><mml:mi>&#x003B8;</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>&#x003BC;</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:msub><mml:mi>&#x003BC;</mml:mi><mml:mtext>f</mml:mtext></mml:msub><mml:msub><mml:mtext>H</mml:mtext><mml:mtext>r</mml:mtext></mml:msub><mml:msub><mml:mtext>H</mml:mtext><mml:mi>&#x003B8;</mml:mi></mml:msub></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>In our case the carrier liquid was non-magnetic and so &#x003BC;<sub>f</sub> &#x0003D; 1. The fields H<sub>r</sub> and H<sub>&#x003B8;</sub> are the field normal and tangential to the internal surface of the particle, respectively. In <xref ref-type="supplementary-material" rid="SM2">Figure S2</xref>, we have compared the dependence of these three models on the amplitude of the magnetic field for a separation distance of 0.01a and a radius of 1 mm. As shown in this figure, the multipolar force (Equation 1) and the FEM calculation (Equation 4) gave the same result at low field as expected but above H &#x0003D; 25 kA/m, the multipolar model overestimates the force. The analytical model (Equation 3) largely overestimates the force. This force calculation is necessary to predict the yield stress, but we also need to know the structure and its deformation under strain.</p>
</sec>
<sec>
<title>Models for the Yield Stress</title>
<sec>
<title>Affine Trajectories</title>
<p>In the standard model the particles are supposed to follow trajectories which are affine with the imposed shear flow. The magnetic field is perpendicular to the imposed flow (cf <xref ref-type="supplementary-material" rid="SM1">Figure S1</xref>). The equilibrium position of the particles is obtained by balancing the shear stress on the top particle with the magnetic restoring force. When the two representative particles are aligned with the magnetic field (&#x003B8; &#x0003D; 0 on <xref ref-type="supplementary-material" rid="SM1">Figure S1</xref>) the gap between the surfaces is given by the thickness of the coating layer, typically of the order of a few nanometers. When particle 2 moves along the stream line, the gap increases with the strain: &#x003B3; &#x0003D; tg(&#x003B8;) as: <inline-formula><mml:math id="M6"><mml:mrow><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>&#x003B5;</mml:mi></mml:mstyle><mml:mo>&#x02212;</mml:mo><mml:msub><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>&#x003B5;</mml:mi></mml:mstyle><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mn>0</mml:mn></mml:mstyle></mml:msub><mml:mo>=</mml:mo><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mn>2</mml:mn></mml:mstyle><mml:mo stretchy='false'>(</mml:mo><mml:msqrt><mml:mrow><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mn>1</mml:mn></mml:mstyle><mml:mo>+</mml:mo><mml:msup><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>&#x003B3;</mml:mi></mml:mstyle><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mn>2</mml:mn></mml:mstyle></mml:msup></mml:mrow></mml:msqrt><mml:mo>&#x02212;</mml:mo><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mn>1</mml:mn></mml:mstyle><mml:mo stretchy='false'>)</mml:mo></mml:mrow></mml:math></inline-formula>. On the other hand, the magnetic force along the x axis is F<sub>r</sub> sin&#x003B8; cos<sup>2</sup>&#x003B8; where F<sub>r</sub> is the radial force calculated in the preceding section for particles aligned in the direction of the field. In a model where all particles form chains spanning the gap with a surface density n<sub>c</sub>/S &#x0003D; 2a N/V &#x0003D; (3/2)&#x003A6;/(&#x003C0;a<sup>2</sup>), with &#x003A6; the volume fraction, the shear stress is given by:</p>
<disp-formula id="E5"><label>(5)</label><mml:math id="M7"><mml:mi>&#x003C4;</mml:mi><mml:mo stretchy='false'>(</mml:mo><mml:mi>&#x003B3;</mml:mi><mml:mo stretchy='false'>)</mml:mo><mml:mo>=</mml:mo><mml:mfrac><mml:mtext>3</mml:mtext><mml:mn>2</mml:mn></mml:mfrac><mml:mo>&#x003A6;</mml:mo><mml:mtext>&#x000A0;</mml:mtext><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi>x</mml:mi><mml:mo stretchy='false'>(</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>F</mml:mi><mml:mi>r</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:mi>&#x003C0;</mml:mi><mml:msup><mml:mi>a</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mfrac><mml:mtext>sin</mml:mtext><mml:mo>&#x003B8;</mml:mo><mml:msup><mml:mtext>cos</mml:mtext><mml:mn>2</mml:mn></mml:msup><mml:mo>&#x003B8;</mml:mo><mml:mo stretchy='false'>)</mml:mo><mml:mo>=</mml:mo><mml:mfrac><mml:mn>3</mml:mn><mml:mn>2</mml:mn></mml:mfrac><mml:mo>&#x003A6;</mml:mo><mml:mtext>&#x000A0;</mml:mtext><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi>x</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:msub><mml:mi>F</mml:mi><mml:mi>r</mml:mi></mml:msub><mml:mo stretchy='false'>(</mml:mo><mml:mi>&#x003B3;</mml:mi><mml:mo stretchy='false'>)</mml:mo></mml:mrow><mml:mrow><mml:mi>&#x003C0;</mml:mi><mml:msup><mml:mi>a</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mfrac><mml:mfrac><mml:mi>&#x003B3;</mml:mi><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:msup><mml:mi>&#x003B3;</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mn>3</mml:mn><mml:mo>/</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>
<p>For the analytical expression of the force Equation (3), we get, (always for &#x003B5;/2 &#x0003C;H/Ms&#x0003C;&#x0003C;1):</p>
<disp-formula id="E6"><label>(6)</label><mml:math id="M8"><mml:msub><mml:mi>&#x003C4;</mml:mi><mml:mrow><mml:mtext>ya</mml:mtext></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msqrt><mml:mrow><mml:mn>16</mml:mn><mml:mo>/</mml:mo><mml:mn>3</mml:mn></mml:mrow></mml:msqrt><mml:mtext>&#x000A0;&#x000A0;</mml:mtext><mml:mo>&#x003A6;</mml:mo><mml:msub><mml:mi>&#x003BC;</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:msub><mml:mtext>M</mml:mtext><mml:mtext>s</mml:mtext></mml:msub><mml:msup><mml:mtext>H</mml:mtext><mml:mrow><mml:mn>3</mml:mn><mml:mo>/</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:math></disp-formula>
<p>Equation (6) differs slightly from that given in Ginder et al. (<xref ref-type="bibr" rid="B18">1996</xref>) where the prefactor is <inline-formula><mml:math id="M9"><mml:msqrt><mml:mrow><mml:mn>6</mml:mn></mml:mrow></mml:msqrt></mml:math></inline-formula>. It is worth noting that when the magnetization of the particles is saturated, the force between the particles is obtained from interactions between rigid dipoles m &#x0003D; M<sub>s</sub>.v, where v is the volume of the particle. The results for the yield stress is readily derived, for non-interacting particle chains (Ginder et al., <xref ref-type="bibr" rid="B18">1996</xref>) as: &#x003C4;<sub>ymax</sub> &#x0003D; 0.086 &#x003A6;.&#x003BC;<sub>0</sub>.<inline-formula><mml:math id="M10"><mml:msubsup><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mrow><mml:mtext>s</mml:mtext></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup></mml:math></inline-formula>.</p>
<p>For the multipolar force we need to interpolate the tables f<sub>//</sub>, f<sub>&#x022A5;</sub>, f<sub><italic>&#x00393;</italic></sub> as a function of the separation and, in the non-linear case derived by FEM, it is also needed to do this interpolation for each value of the field; the data used are given in <xref ref-type="supplementary-material" rid="SM6">Table S3</xref>. The results are presented in <xref ref-type="fig" rid="F1">Figure 1A</xref> but, before discussing them, it is worth introducing another model, which to our knowledge, has never been used for MR (or ER) fluids.</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p><bold>(A)</bold> Theoretical yield stress vs. the field acting on the particles. Red squares: calculated from FEMM; green triangles calculated from multipolar forces; black solid line: Equation (6); blue crosses: Equation (7) with F<sub>r</sub> calculated with FEM. The red dot corresponds to the saturation of the magnetization. <bold>(B)</bold> Comparison theory-experiment for &#x003C4;<sub>y</sub>/&#x003A6;. Blue diamonds: experimental points taken from <xref ref-type="supplementary-material" rid="SM4">Table S1</xref>. Red line, FEM result; black line, Equation (6); Blue line, Equation (7); Purple star, jamming transition.</p></caption>
<graphic xlink:href="fmats-06-00004-g0001.tif"/>
</fig>
</sec>
<sec>
<title>Rupture at Zero Strain</title>
<p>Here we suppose that the particles are packed in some disordered way, and that the solid phase will yield suddenly when the shear stress is high enough to separate the particles on a given &#x0201C;cutting&#x0201D; plane which, in a simple shear, is perpendicular to the velocity gradient. This kind of model is used to predict the yield stress of suspensions of particles aggregated in the presence of Van der Waals forces (Scales et al., <xref ref-type="bibr" rid="B35">1998</xref>). The only difference in our case is that the interparticle force is proportional to (Hcos&#x003B8;)<sup>2</sup> instead of being independent of &#x003B8;. Otherwise the interparticle distances are supposed to be distributed isotropically, which is mainly true at a high-volume fraction. Following the same steps as in Scales et al. (<xref ref-type="bibr" rid="B35">1998</xref>) we obtain:</p>
<disp-formula id="E7"><label>(7)</label><mml:math id="M11"><mml:msub><mml:mi>&#x003C4;</mml:mi><mml:mtext>y</mml:mtext></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mtext>3</mml:mtext><mml:mo>&#x003A6;</mml:mo><mml:mtext>&#x000A0;K</mml:mtext><mml:mo stretchy='false'>(</mml:mo><mml:mo>&#x003A6;</mml:mo><mml:mo stretchy='false'>)</mml:mo></mml:mrow><mml:mrow><mml:mtext>20</mml:mtext><mml:mi>&#x003C0;</mml:mi><mml:msup><mml:mtext>R</mml:mtext><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mfrac><mml:msub><mml:mtext>F</mml:mtext><mml:mtext>r</mml:mtext></mml:msub><mml:mo stretchy='false'>(</mml:mo><mml:msub><mml:mi>&#x003B5;</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo stretchy='false'>)</mml:mo><mml:mtext>&#x02003;with&#x02003;K</mml:mtext><mml:mo stretchy='false'>(</mml:mo><mml:mo>&#x003A6;</mml:mo><mml:mo stretchy='false'>)</mml:mo><mml:mo>=</mml:mo><mml:mn>36</mml:mn><mml:mo>&#x003A6;</mml:mo><mml:mo>/</mml:mo><mml:mo>&#x003C0;</mml:mo></mml:math></disp-formula>
<p>The factor 3/20 comes from the cos<sup>2</sup>&#x003B8; instead of 1/4 for an isotropic force. The function K(&#x003A6;) is the coordination number taken from Suzuki et al. (<xref ref-type="bibr" rid="B36">1981</xref>) and &#x003B5;<sub>0</sub> is the minimum gap between two particles.</p>
</sec>
<sec>
<title>Comparison of the Different Models</title>
<p>The predictions of the different models are reported in <xref ref-type="fig" rid="F1">Figure 1A</xref>. Firstly if we compare, in the affine model, the predictions related to the 3 different approximations of forces: Equations (1, 3, 4),we recover the fact that the multipolar model gives a rather good prediction for low fields but, due to the linear magnetization hypothesis, highly overestimates the yield stress at high fields. The analytic expression overestimates the yield stress given by FEM, by about one order of a magnitude. Lastly, the model of rupture at zero strain (Equation 7), is calculated with the force F<sub>r</sub> obtained by FEM for &#x003B5; &#x0003D; 0.002, which is a reasonable value considering a coating layer of 2 nm and a radius of 1 &#x003BC;m. The corresponding prediction is well-above the FEM prediction and is, by chance, quite close to that of the analytical model. It should also be noted that the fit of the yield stress by a power law:&#x003C4;<sub>y</sub> &#x0003D; K.H<sup>n</sup> provides a good fit with <italic>n</italic> &#x0003D; 1.94 (<italic>R</italic><sup>2</sup> &#x0003D; 0.9997) only for H &#x02264; 50 kA/m (<xref ref-type="supplementary-material" rid="SM3">Figure S3</xref>) which is three times smaller than the field used for the transition from H<sup>2</sup> to H<sup>3/2</sup> (H<sub>c</sub> &#x0003D; 151 kA/m) in Fang et al. (<xref ref-type="bibr" rid="B17">2009</xref>). It is also worth noting that the radial force, at contact between two particles, increases steadily in the range of &#x003B5; &#x0003D; 0&#x02013;0.01(cf <xref ref-type="supplementary-material" rid="SM5">Table S2</xref>), which may explain, at least in the rupture model, the decrease of the yield stress in the presence of a coating layer (Cvek et al., <xref ref-type="bibr" rid="B11">2018</xref>). This comparison of the models was done for two spheres in a given field, which we consider as the average field, in order to compare with the experiments in the next section. Some numerical results for the yield stress of a set of particles can be obtained with a good precision in a multipolar approach, even for high permeability, but for particles distributed on a periodic network (Clercx and Bossis, <xref ref-type="bibr" rid="B10">1993</xref>) with an affine deformation. In this case all the particles are separated by the same gap, which is an unphysical and unstable structure. It is also possible to use a double chain structure so that it can deform and still keep contacts and stability, but it underestimates the yield stress (L&#x000F3;pez-L&#x000F3;pez et al., <xref ref-type="bibr" rid="B31">2012</xref>) and also can&#x00027;t be used in numerical simulations. In fact, almost all numerical simulations of MR suspensions use the dipolar approximation, which largely underestimates the magnetic force (Clercx and Bossis, <xref ref-type="bibr" rid="B10">1993</xref>). We have found one paper (Lacis and Gosko, <xref ref-type="bibr" rid="B27">2009</xref>) using an expression similar to Equation (1), with some tables derived from FEM simulations for two particles and another based on FEM simulations between several particles in a single chain (Kang et al., <xref ref-type="bibr" rid="B23">2012</xref>). A simulation method, called smooth particle dynamics (Hashemi et al., <xref ref-type="bibr" rid="B19">2016</xref>, <xref ref-type="bibr" rid="B20">2018</xref>), can also be used to solve the magnetostatics equation between several particles, but is restricted to a few particles with a rather small (&#x0003C;10) susceptibility. As discussed in the next section, the use of a magnetic force derived from FEM, calculated with a very fine mesh for two particles, is quite realistic and certainly much more efficient in terms of computing time.</p>
</sec>
</sec>
</sec>
<sec id="s3">
<title>Comparison With Experiments</title>
<p>Comparing these predictions with experimental data is not so simple as, very often, the knowledge of the experimental conditions is incomplete. We have selected data concerning carbonyl iron particles (CI) whose magnetization curve of the powder is rather well-established with a saturation magnetization of 200 emu/g. Usually the experiments are made in plate-plate geometry, with the field perpendicular to the plate. In order to compare with the theory, we must first consider that the field inside the suspension is H &#x0003D; H<sub>0</sub>/&#x003BC;(H), where H<sub>0</sub> is the external field measured in the absence of suspension and &#x003BC;(H) the permeability of the suspension. We have deduced the magnetization curve of CI particles from a measurement of M(H) at a volume fraction: &#x003A6; &#x0003D; 64%. This magnetization curve of CI particles is well-fitted in the range 0&#x02013;400 kA/m by a Langevin curve:M<sub>1</sub>(H) &#x0003D; (coth(c.H)-1/c.H)M<sub>sr</sub> with c &#x0003D; 0.014 and M<sub>sr</sub> &#x0003D; 1750 kA/m. This last value is slightly higher than the saturation magnetization of CI (M<sub>s</sub> &#x02248; 1600 kA/m) due to a shift of the fit for H&#x0003E;400 kA/m. The magnetization at any volume fraction of the suspension is obtained by replacing M<sub>sr</sub> by M<sub>sr</sub>.&#x003A6; and we have checked that it works well-even at low &#x003A6;, for instance for &#x003A6; &#x0003D; 0.2 and &#x003A6; &#x0003D; 0.32 BASF. The permeability for any given volume fraction is then: &#x003BC;(H) &#x0003D; 1&#x0002B;M<sub>1</sub>(H).&#x003A6;/H The second point we need to take into account is that, in the plate-plate configuration, the yield stress provided by the software is overestimated by a factor of 4/3 (Bossis et al., <xref ref-type="bibr" rid="B7">2002</xref>; Laun et al., <xref ref-type="bibr" rid="B28">2010</xref>). The data collected in different publications (Kordonski and Golini, <xref ref-type="bibr" rid="B26">1999</xref>; Chin et al., <xref ref-type="bibr" rid="B9">2001</xref>; Lim et al., <xref ref-type="bibr" rid="B29">2004</xref>; Jang et al., <xref ref-type="bibr" rid="B22">2005</xref>; Cheng et al., <xref ref-type="bibr" rid="B8">2009</xref>; Ierardi and Bombard, <xref ref-type="bibr" rid="B21">2009</xref>; Laun et al., <xref ref-type="bibr" rid="B28">2010</xref>; Kim et al., <xref ref-type="bibr" rid="B24">2012</xref>; Bombard et al., <xref ref-type="bibr" rid="B4">2015</xref>; Esmaeilnezhad et al., <xref ref-type="bibr" rid="B14">2017</xref>; Dong et al., <xref ref-type="bibr" rid="B13">2018</xref>; Plachy et al., <xref ref-type="bibr" rid="B34">2018</xref>) are reported in <xref ref-type="supplementary-material" rid="SM4">Table S1</xref> and plotted in <xref ref-type="fig" rid="F1">Figure 1B</xref> (blue diamonds). There is a considerable scatter of the data, but it appears clearly that the affine model with the calculation of forces through FEM (red curve), provides the best prediction compared to both Equation (6) and Equation (7). Now, we draw the attention on the single point (purple cross) above all the curves in <xref ref-type="fig" rid="F1">Figure 1B</xref> which is the result of a jamming transition.</p>
</sec>
<sec id="s4">
<title>Jamming Transition at High Volume Fraction</title>
<p>Jamming transition, also known as discontinuous shear thickening (DST), was observed in highly concentrated suspensions of solid non-Brownian particles like ceramics (Bergstr&#x000F6;m, <xref ref-type="bibr" rid="B3">1998</xref>), cornstarch (Fall et al., <xref ref-type="bibr" rid="B16">2008</xref>, <xref ref-type="bibr" rid="B15">2015</xref>) and also in model suspensions made of PMMA (Pan et al., <xref ref-type="bibr" rid="B33">2015</xref>) or silica particles (Lootens et al., <xref ref-type="bibr" rid="B30">2003</xref>) and calcium carbonate (Bossis et al., <xref ref-type="bibr" rid="B5">2017</xref>). In the two last systems, superplasticizer molecules, common in the cement industry, were used as surfactants. We have recently found that by using, for instance, a polyoxyethylene polyphosphonate (Bossis et al., <xref ref-type="bibr" rid="B5">2017</xref>), we can reach volume fractions of CI particles in water larger than 60%, still keeping a low yield stress and a viscosity of 1-2Pa.s. Now the rheogram at the zero field presents a singularity at a critical stress, where the shear rate abruptly decreases during a ramp of stress (<xref ref-type="fig" rid="F2">Figure 2A</xref>). This DST phenomenon is due to the transition of a lubricated regime to a frictional one (Mari et al., <xref ref-type="bibr" rid="B32">2014</xref>; Wyart and Cates, <xref ref-type="bibr" rid="B38">2014</xref>). Above this point the shear rate fluctuates around some low values. We also see on this graph that, by increasing the volume fraction, the critical stress increases and the critical shear rate decreases. These critical values are obtained at zero magnetic field, but when the field is increased for a given volume fraction, the critical shear rate will continuously decrease to zero in a field range typically of 0&#x02013;20 kA/m (Bossis et al., <xref ref-type="bibr" rid="B6">2016</xref>). If instead of applying a ramp of stress we apply a ramp of shear rate, at the critical point we then obtain a jump of stress to very high values which overcomes the range of conventional rheometers. This is the purple point in <xref ref-type="fig" rid="F1">Figure 1B</xref>. which was obtained in the presence of a field of 19.2 kA/m. The corresponding evolution of stress is plotted in <xref ref-type="fig" rid="F2">Figure 2B</xref> vs. the strain. It was obtained in a cylindrical Couette flow with a vane tool and a homemade rheometer to be able to measure high torques (Bossis et al., <xref ref-type="bibr" rid="B6">2016</xref>). The shear rate was increased from 0 to 1.6 s<sup>&#x02212;1</sup> in 2 s then kept constant for 25 s. The stress does not grow until the strain reaches a value around one and then increases up to 150 kPa where it remains constant afterwards. The interesting point is that this transition from a normal regime to this frictional regime can be controlled by the application of a low magnetic field (a few kA/m), thus opening the door for new applications of MR fluids. It is worth noting that a very high value of yield stress in a MR fluid was previously attained in another way, which also mobilizes the friction between particles, but with the help of a lateral pressure obtained by the compression of the suspension of the initial volume fraction &#x003A6; &#x0003D; 45% by sliding wedges (Tao, <xref ref-type="bibr" rid="B37">2001</xref>). As reported in numerous previous works, [for instance (Ashtiani et al., <xref ref-type="bibr" rid="B1">2015</xref>; Cvek et al., <xref ref-type="bibr" rid="B11">2018</xref>)] the additives are used to stabilize the particles against corrosion and sedimentation; they can either increase the yield stress through a better compaction of the local structure [term K[&#x003A6;] in Equation (7)], or decrease it if the coating layer is too thick [term F<sub>r</sub>(&#x003B5;<sub>0</sub>) in Equation (7)]. In our case it plays a key role due its expulsion from the surface, which provokes the jamming transition.</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p><bold>(A)</bold> Rheogram of carbonyl iron suspensions with 0.2% of superplasticizer at high-volume fraction (blue:60%, orange 61%, black 62%,green 63%). <bold>(B)</bold> Stress vs. strain in cylindrical Couette geometry with vane tool, &#x003A6; &#x0003D; 0.61;H<sub>0</sub> &#x0003D; 19.2 kA/m, The vertical red line corresponds to a time of 2 s where the shear rate has reached its final value of 1.6 s<sup>&#x02212;1</sup>.</p></caption>
<graphic xlink:href="fmats-06-00004-g0002.tif"/>
</fig>
</sec>
<sec sec-type="conclusions" id="s5">
<title>Conclusion</title>
<p>In this paper we have briefly reviewed the models of yield stress under the hypothesis of affine displacement of the particles and have also proposed a new approach based on a rupture model. It appears that the affine motion, with a calculation of the forces through FEM simulation, provides the best agreement with the experimental data of &#x003C4;<sub>y</sub>/&#x003A6;, collected for carbonyl iron suspensions. The commonly used Equation (6) strongly overestimates the experimental results but its prediction of a regime &#x003C4;<sub>y</sub>&#x0221D; H<sup>3/2</sup> is verified between &#x0007E;50 and 200 kPa. No model is able to predict the high value of stress at low field obtained in suspensions of very high-volume fractions: &#x003A6; &#x0003E;0.6. In this case, it is the abrupt transition toward a frictional regime after expulsion of a superplasticizer layer, which leads to very high values of stress. This transition is obtained either above a critical magnetic field or above a critical shear rate, which opens up new perspectives for applications that require a strong response with a small field.</p>
</sec>
<sec id="s6">
<title>Author Contributions</title>
<p>GB, OV, YG, and AC contributed conception and design of the study. GB wrote the first draft of the paper. OV found and organized the references and the figures. YG filled the database of the previous works on the subject and AC contributed to the theoretical part. All authors contributed to manuscript revision, read and approved the submitted version.</p>
<sec>
<title>Conflict of Interest Statement</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
</sec>
</body>
<back>
<ack><p>This work was supported by the Center National d&#x00027;Etudes Spatiales CNES. We also thank the mechanical workshop of our institute for fabricating part of our home made rheometer.</p>
</ack>
<sec sec-type="supplementary-material" id="s8">
<title>Supplementary Material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fmats.2019.00004/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fmats.2019.00004/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Image_1.tif" id="SM1" mimetype="image/tif" xmlns:xlink="http://www.w3.org/1999/xlink">
<label>Figure S1</label>
<caption><p>Schematic of the notations for two particles in a magnetic field.</p></caption> </supplementary-material>
<supplementary-material xlink:href="Image_2.tif" id="SM2" mimetype="image/tif" xmlns:xlink="http://www.w3.org/1999/xlink">
<label>Figure S2</label>
<caption><p>Force vs. magnetic field for a separation &#x003B5; &#x0003D; 0.01; black line: analytic expression: Equation (3); green triangles: multipolar approach: Equation (1); red triangles: FEM, Equation (4).</p></caption> </supplementary-material>
<supplementary-material xlink:href="Image_3.tif" id="SM3" mimetype="image/tif" xmlns:xlink="http://www.w3.org/1999/xlink">
<label>Figure S3</label>
<caption><p>Power law fit of the FEM result for the normalized yield stress &#x003C4;<sub>y</sub>/&#x003A6;. Red solid line:power law with <italic>n</italic> &#x0003D; 1.94; blue solid line: power law with <italic>n</italic> &#x0003D; 3/2;red crosses: FEM result; black solid line: Equation (6).</p></caption> </supplementary-material>
<supplementary-material xlink:href="Table_1.XLS" id="SM4" mimetype="application/vnd.ms-excel" xmlns:xlink="http://www.w3.org/1999/xlink">
<label>Table S1</label>
<caption><p>Corrected experimental yield stress normalized by the volume fraction (column E) in kPa vs. the Maxwell field in the suspension (column H) in kA/m.</p></caption> </supplementary-material>
<supplementary-material xlink:href="Table_2.XLS" id="SM5" mimetype="application/vnd.ms-excel" xmlns:xlink="http://www.w3.org/1999/xlink">
<label>Table S2</label>
<caption><p>Coefficients of the force in Equation (1) for differents values of &#x003B1; &#x0003D; &#x003BC;<sub>p</sub>/&#x003BC;<sub>f</sub> vs. the gap &#x003B5;.</p></caption> </supplementary-material>
<supplementary-material xlink:href="Table_3.XLS" id="SM6" mimetype="application/vnd.ms-excel" xmlns:xlink="http://www.w3.org/1999/xlink">
<label>Table S3</label>
<caption><p>Values of the force obtained with the free software FEMM for different fields and separations between two spheres. The parameters for the fit vs. the normalized separation, &#x003BE;, (&#x0003C;0.2) are given for each field.</p></caption> </supplementary-material>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ashtiani</surname> <given-names>M.</given-names></name> <name><surname>Hashemabadi</surname> <given-names>S. H.</given-names></name> <name><surname>Ghaffari</surname> <given-names>A.</given-names></name></person-group> (<year>2015</year>). <article-title>A review on the magnetorheological fluid preparation and stabilization</article-title>. <source>J. Magn. Magn. Mater.</source> <volume>374</volume>, <fpage>716</fpage>&#x02013;<lpage>730</lpage>. <pub-id pub-id-type="doi">10.1016/j.jmmm.2014.09.020</pub-id></citation></ref>
<ref id="B2">
<citation citation-type="other"><person-group person-group-type="author"><collab>BASF</collab></person-group> (<year>2010</year>). <source>Technical Data Sheet MR2899</source>. <pub-id pub-id-type="pmid">15296025</pub-id></citation></ref>
<ref id="B3">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bergstr&#x000F6;m</surname> <given-names>L.</given-names></name></person-group> (<year>1998</year>). <article-title>Shear thinning and shear thickening of concentrated ceramic suspensions</article-title>. <source>Colloids Surf. A</source> <volume>133</volume>, <fpage>151</fpage>&#x02013;<lpage>155</lpage>. <pub-id pub-id-type="doi">10.1016/S0927-7757(97)00133-7</pub-id></citation></ref>
<ref id="B4">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bombard</surname> <given-names>A. J.</given-names></name> <name><surname>Gon&#x000E7;alves</surname> <given-names>F. R.</given-names></name> <name><surname>de Vicente</surname> <given-names>J.</given-names></name></person-group> (<year>2015</year>). <article-title>Magnetorheology of carbonyl iron dispersions in 1-Alkyl-3-methylimidazolium ionic liquids</article-title>. <source>Ind. Eng. Chem. Res.</source> <volume>54</volume>, <fpage>9956</fpage>&#x02013;<lpage>9963</lpage>. <pub-id pub-id-type="doi">10.1021/acs.iecr.5b02824</pub-id></citation></ref>
<ref id="B5">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bossis</surname> <given-names>G.</given-names></name> <name><surname>Boustingorry</surname> <given-names>P.</given-names></name> <name><surname>Grasselli</surname> <given-names>Y.</given-names></name> <name><surname>Meunier</surname> <given-names>A.</given-names></name> <name><surname>Morini</surname> <given-names>R.</given-names></name> <name><surname>Zubarev</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Discontinuous shear thickening in the presence of polymers adsorbed on the surface of calcium carbonate particles</article-title>. <source>Rheol. Acta</source> <volume>56</volume>, <fpage>415</fpage>&#x02013;<lpage>430</lpage>. <pub-id pub-id-type="doi">10.1007/s00397-017-1005-4</pub-id></citation></ref>
<ref id="B6">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bossis</surname> <given-names>G.</given-names></name> <name><surname>Grasselli</surname> <given-names>Y.</given-names></name> <name><surname>Meunier</surname> <given-names>A.</given-names></name> <name><surname>Volkova</surname> <given-names>O.</given-names></name></person-group> (<year>2016</year>). <article-title>Outstanding magnetorheological effect based on discontinuous shear thickening in the presence of a superplastifier molecule</article-title>. <source>Appl. Phys. Lett.</source> <volume>109</volume>:<fpage>111902</fpage>. <pub-id pub-id-type="doi">10.1063/1.4962467</pub-id></citation></ref>
<ref id="B7">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Bossis</surname> <given-names>G.</given-names></name> <name><surname>Volkova</surname> <given-names>O.</given-names></name> <name><surname>Lacis</surname> <given-names>S.</given-names></name> <name><surname>Meunier</surname> <given-names>A.</given-names></name></person-group> (<year>2002</year>). <article-title>Magnetorheology: fluids, structures and rheology</article-title>, in <source>Ferrofluids, Part of Lecture Notes in Physics, Vol. 594</source> (<publisher-loc>Springer</publisher-loc>), <fpage>202</fpage>&#x02013;<lpage>230</lpage>.</citation></ref>
<ref id="B8">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cheng</surname> <given-names>H. B.</given-names></name> <name><surname>Wang</surname> <given-names>J. M.</given-names></name> <name><surname>Zhang</surname> <given-names>Q. J.</given-names></name> <name><surname>Wereley</surname> <given-names>N. M.</given-names></name></person-group> (<year>2009</year>). <article-title>Preparation of composite magnetic particles and aqueous magnetorheological fluids</article-title>. <source>Smart Mater Struct.</source> <volume>18</volume>:<fpage>085009</fpage>. <pub-id pub-id-type="doi">10.1088/0964-1726/18/8/085009</pub-id></citation></ref>
<ref id="B9">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chin</surname> <given-names>B. D.</given-names></name> <name><surname>Park</surname> <given-names>J. H.</given-names></name> <name><surname>Kwon</surname> <given-names>M. H.</given-names></name> <name><surname>Park</surname> <given-names>O. O.</given-names></name></person-group> (<year>2001</year>). <article-title>Rheological properties and dispersion stability of magnetorheological (MR) suspensions</article-title>. <source>Rheologica Acta</source> <volume>40</volume>, <fpage>211</fpage>&#x02013;<lpage>219</lpage>. <pub-id pub-id-type="doi">10.1007/s003970000150</pub-id></citation></ref>
<ref id="B10">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Clercx</surname> <given-names>H. J. H.</given-names></name> <name><surname>Bossis</surname> <given-names>G.</given-names></name></person-group> (<year>1993</year>). <article-title>Electrostatic interactions in slabs of polarizable particles</article-title>. <source>J. Chem. Phys.</source> <volume>98</volume>, <fpage>8284</fpage>&#x02013;<lpage>8293</lpage>. <pub-id pub-id-type="doi">10.1063/1.464534</pub-id></citation></ref>
<ref id="B11">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cvek</surname> <given-names>M.</given-names></name> <name><surname>Mrl&#x000ED;k</surname> <given-names>M.</given-names></name> <name><surname>Moucka</surname> <given-names>R.</given-names></name> <name><surname>Sedlacik</surname> <given-names>M.</given-names></name></person-group> (<year>2018</year>). <article-title>A systematical study of the overall influence of carbon allotrope additives on performance, stability and redispersibility of magnetorheological fluids</article-title>. <source>Colloids Surf. A</source> <volume>543</volume>, <fpage>83</fpage>&#x02013;<lpage>92</lpage>. <pub-id pub-id-type="doi">10.1016/j.colsurfa.2018.01.046</pub-id></citation></ref>
<ref id="B12">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>De Vicente</surname> <given-names>J.</given-names></name> <name><surname>Bossis</surname> <given-names>G.</given-names></name> <name><surname>Lacis</surname> <given-names>S.</given-names></name> <name><surname>Guyot</surname> <given-names>M.</given-names></name></person-group> (<year>2002</year>). <article-title>Permeability measurements in cobalt ferrite and carbonyl iron powders and suspensions</article-title>. <source>J. Magn. Magn. Mater.</source> <volume>251</volume>, <fpage>100</fpage>&#x02013;<lpage>108</lpage>. <pub-id pub-id-type="doi">10.1016/S0304-8853(02)00484-5</pub-id></citation></ref>
<ref id="B13">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dong</surname> <given-names>Y. Z.</given-names></name> <name><surname>Piao</surname> <given-names>S. H.</given-names></name> <name><surname>Zhang</surname> <given-names>K.</given-names></name> <name><surname>Choi</surname> <given-names>H. J.</given-names></name></person-group> (<year>2018</year>). <article-title>Effect of CoFe2O4 nanoparticles on a carbonyl iron based magnetorheological suspension</article-title>. <source>Colloids Surf. A</source> <volume>537</volume>, <fpage>102</fpage>&#x02013;<lpage>108</lpage>. <pub-id pub-id-type="doi">10.1016/j.colsurfa.2017.10.017</pub-id></citation></ref>
<ref id="B14">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Esmaeilnezhad</surname> <given-names>E.</given-names></name> <name><surname>Jin Choi</surname> <given-names>H.</given-names></name> <name><surname>Schaffie</surname> <given-names>M.</given-names></name> <name><surname>Gholizadeh</surname> <given-names>M.</given-names></name> <name><surname>Ranjbar</surname> <given-names>M.</given-names></name> <name><surname>Hyuk Kwon</surname> <given-names>S.</given-names></name></person-group> (<year>2017</year>). <article-title>Rheological analysis of magnetite added carbonyl iron based magnetorheological fluid</article-title>. <source>J. Magnetism Magnetic Mater.</source> <volume>444</volume>, <fpage>161</fpage>&#x02013;<lpage>167</lpage>. <pub-id pub-id-type="doi">10.1016/j.jmmm.2017.08.023</pub-id></citation></ref>
<ref id="B15">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fall</surname> <given-names>A.</given-names></name> <name><surname>Bertrand</surname> <given-names>F.</given-names></name> <name><surname>Hautemayou</surname> <given-names>D.</given-names></name> <name><surname>M&#x000E9;zi&#x000E8;re</surname> <given-names>C.</given-names></name> <name><surname>Moucheront</surname> <given-names>P.</given-names></name> <name><surname>Lemaitre</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2015</year>). <article-title>Macroscopic discontinuous shear thickening versus local shear jamming in cornstarch</article-title>. <source>Phys. Rev. Lett.</source> <volume>114</volume>:<fpage>098301</fpage>. <pub-id pub-id-type="doi">10.1103/PhysRevLett.114.098301</pub-id><pub-id pub-id-type="pmid">25793857</pub-id></citation></ref>
<ref id="B16">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fall</surname> <given-names>A.</given-names></name> <name><surname>Huang</surname> <given-names>N.</given-names></name> <name><surname>Bertrand</surname> <given-names>F.</given-names></name> <name><surname>Ovarlez</surname> <given-names>G.</given-names></name> <name><surname>Bonn</surname> <given-names>D.</given-names></name></person-group> (<year>2008</year>). <article-title>Shear thickening of cornstarch suspensions as a reentrant jamming transition</article-title>. <source>Phys. Rev. Lett.</source> <volume>100</volume>:<fpage>018301</fpage>. <pub-id pub-id-type="doi">10.1103/PhysRevLett.100.018301</pub-id><pub-id pub-id-type="pmid">18232829</pub-id></citation></ref>
<ref id="B17">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fang</surname> <given-names>F. F.</given-names></name> <name><surname>Choi</surname> <given-names>H. J.</given-names></name> <name><surname>Jhon</surname> <given-names>M. S.</given-names></name></person-group> (<year>2009</year>). <article-title>Magnetorheology of soft magnetic carbonyl iron suspension with single-walled carbon nanotube additive and its yield stress scaling function</article-title>. <source>Colloids Surf. A</source> <volume>351</volume>, <fpage>46</fpage>&#x02013;<lpage>51</lpage>. <pub-id pub-id-type="doi">10.1016/j.colsurfa.2009.09.032</pub-id></citation></ref>
<ref id="B18">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ginder</surname> <given-names>J. M.</given-names></name> <name><surname>Davis</surname> <given-names>L. C.</given-names></name> <name><surname>Elie</surname> <given-names>L. D.</given-names></name></person-group> (<year>1996</year>). <article-title>Rheology of magnetorheological fluids: models and measurements</article-title>. <source>Int. J. Mod. Phys. B</source> <volume>10</volume>, <fpage>3293</fpage>&#x02013;<lpage>3303</lpage>. <pub-id pub-id-type="doi">10.1142/S0217979296001744</pub-id></citation></ref>
<ref id="B19">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hashemi</surname> <given-names>M. R.</given-names></name> <name><surname>Manzari</surname> <given-names>M. T.</given-names></name> <name><surname>Fatehi</surname> <given-names>R.</given-names></name></person-group> (<year>2016</year>). <article-title>A SPH solver for simulating paramagnetic solid fluid interaction in the presence of an external magnetic field</article-title>. <source>Appl. Math. Model.</source> <volume>40</volume>, <fpage>4341</fpage>&#x02013;<lpage>4369</lpage>. <pub-id pub-id-type="doi">10.1016/j.apm.2015.11.020</pub-id></citation></ref>
<ref id="B20">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hashemi</surname> <given-names>M. R.</given-names></name> <name><surname>Manzari</surname> <given-names>M. T.</given-names></name> <name><surname>Fatehi</surname> <given-names>R.</given-names></name></person-group> (<year>2018</year>). <article-title>Direct numerical simulation of magnetic particles suspended in a Newtonian fluid exhibiting finite inertia under SAOS</article-title>. <source>J. Non-Newtonian Fluid Mech.</source> <volume>256</volume>, <fpage>8</fpage>&#x02013;<lpage>22</lpage>. <pub-id pub-id-type="doi">10.1016/j.jnnfm.2018.03.004</pub-id></citation></ref>
<ref id="B21">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ierardi</surname> <given-names>R. F.</given-names></name> <name><surname>Bombard</surname> <given-names>A. J.</given-names></name></person-group> (<year>2009</year>). <article-title>Off-state viscosity and yield stress optimization of magneto-rheological fluids: a mixture design of experiments approach</article-title>. <source>J. Phys. Conf. Series</source> <volume>149</volume>:<fpage>012037</fpage>. <pub-id pub-id-type="doi">10.1088/1742-6596/149/1/012037</pub-id></citation></ref>
<ref id="B22">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jang</surname> <given-names>I. B.</given-names></name> <name><surname>Kim</surname> <given-names>H. B.</given-names></name> <name><surname>Lee</surname> <given-names>J. Y.</given-names></name> <name><surname>You</surname> <given-names>J. L.</given-names></name> <name><surname>Choi</surname> <given-names>H. J.</given-names></name> <name><surname>Jhon</surname> <given-names>M. S.</given-names></name></person-group> (<year>2005</year>). <article-title>Role of organic coating on carbonyl iron suspended particles in magnetorheological fluids</article-title>. <source>J. Appl. Phys.</source> <volume>97</volume>:<fpage>10Q</fpage>912. <pub-id pub-id-type="doi">10.1063/1.1853835</pub-id></citation></ref>
<ref id="B23">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kang</surname> <given-names>T. G.</given-names></name> <name><surname>Hulsen</surname> <given-names>M. A.</given-names></name> <name><surname>den Toonder</surname> <given-names>J. M.</given-names></name></person-group> (<year>2012</year>). <article-title>Dynamics of magnetic chains in a shear flow under the influence of a uniform magnetic field</article-title>. <source>Phys. Fluids</source> <volume>24</volume>:<fpage>042001</fpage>. <pub-id pub-id-type="doi">10.1063/1.4704822</pub-id></citation></ref>
<ref id="B24">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kim</surname> <given-names>M. S.</given-names></name> <name><surname>Liu</surname> <given-names>Y. D.</given-names></name> <name><surname>Park</surname> <given-names>B. J.</given-names></name> <name><surname>You</surname> <given-names>C.-Y.</given-names></name> <name><surname>Choi</surname> <given-names>H. J.</given-names></name></person-group> (<year>2012</year>). <article-title>Carbonyl iron particles dispersed in a polymer solution and their rheological characteristics under applied magnetic field</article-title>. <source>J. Indus. Eng. Chem.</source> <volume>18</volume>, <fpage>664</fpage>&#x02013;<lpage>667</lpage>. <pub-id pub-id-type="doi">10.1016/j.jiec.2011.11.062</pub-id></citation></ref>
<ref id="B25">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Klingenberg</surname> <given-names>D. J.</given-names></name> <name><surname>Zukoski</surname> <given-names>C. F.</given-names></name></person-group> (<year>1990</year>). <article-title>Studies on the steady-shear behavior of electrorheological suspensions</article-title>. <source>Langmuir</source> <volume>6</volume>, <fpage>15</fpage>&#x02013;<lpage>24</lpage>. <pub-id pub-id-type="doi">10.1021/la00091a003</pub-id></citation></ref>
<ref id="B26">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kordonski</surname> <given-names>W. I.</given-names></name> <name><surname>Golini</surname> <given-names>D.</given-names></name></person-group> (<year>1999</year>). <article-title>Fundamentals of Magnetorheological fluid utilization in high precision finishing</article-title>. <source>J. Intell. Mater. Syst. Struct.</source> <volume>10</volume>, <fpage>683</fpage>&#x02013;<lpage>689</lpage>. <pub-id pub-id-type="doi">10.1106/011M-CJ25-64QC-F3A6</pub-id></citation></ref>
<ref id="B27">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lacis</surname> <given-names>S.</given-names></name> <name><surname>Gosko</surname> <given-names>D.</given-names></name></person-group> (<year>2009</year>). <article-title>Direct numerical simulation of MR suspension: the role of viscous and magnetic interactions between particles</article-title>. <source>J. Phys.</source> <volume>149</volume>:<fpage>012066</fpage>. <pub-id pub-id-type="doi">10.1088/1742-6596/149/1/012066</pub-id></citation></ref>
<ref id="B28">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Laun</surname> <given-names>H. M.</given-names></name> <name><surname>Gabriel</surname> <given-names>C.</given-names></name> <name><surname>Kieburg</surname> <given-names>C.</given-names></name></person-group> (<year>2010</year>). <article-title>Twin gap magnetorheometer using ferromagnetic steel plates&#x02014;performance and validation</article-title>. <source>J. Rheol.</source> <volume>54</volume>, <fpage>327</fpage>&#x02013;<lpage>354</lpage>. <pub-id pub-id-type="doi">10.1122/1.3302804</pub-id></citation></ref>
<ref id="B29">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lim</surname> <given-names>S. T.</given-names></name> <name><surname>Cho</surname> <given-names>M. S.</given-names></name> <name><surname>Jang</surname> <given-names>I. B.</given-names></name> <name><surname>Choi</surname> <given-names>H. J.</given-names></name></person-group> (<year>2004</year>). <article-title>Magnetorheological characterization of carbonyl iron based suspension stabilized by fumed silica</article-title>. <source>J. Magnetism Magnetic Mater.</source> <volume>282</volume>, <fpage>170</fpage>&#x02013;<lpage>173</lpage>. <pub-id pub-id-type="doi">10.1016/j.jmmm.2004.04.040</pub-id></citation></ref>
<ref id="B30">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lootens</surname> <given-names>D.</given-names></name> <name><surname>Van Damme</surname> <given-names>H.</given-names></name> <name><surname>H&#x000E9;braud</surname> <given-names>P.</given-names></name></person-group> (<year>2003</year>). <article-title>Giant stress fluctuations at the jamming transition</article-title>. <source>Phys. Rev. Lett.</source> <volume>90</volume>:<fpage>178301</fpage>. <pub-id pub-id-type="doi">10.1103/PhysRevLett.90.178301</pub-id><pub-id pub-id-type="pmid">12786109</pub-id></citation></ref>
<ref id="B31">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>L&#x000F3;pez-L&#x000F3;pez</surname> <given-names>M. T.</given-names></name> <name><surname>Kuzhir</surname> <given-names>P.</given-names></name> <name><surname>Caballero-Hernandez</surname> <given-names>J.</given-names></name> <name><surname>Rodr&#x000ED;guez-Arco</surname> <given-names>L.</given-names></name> <name><surname>Dur&#x000E1;n</surname> <given-names>J. D.</given-names></name> <name><surname>Bossis</surname> <given-names>G.</given-names></name></person-group> (<year>2012</year>). <article-title>Yield stress in magnetorheological suspensions near the limit of maximum-packing fraction</article-title>. <source>J. Rheol.</source> <volume>56</volume>:<fpage>1209</fpage>. <pub-id pub-id-type="doi">10.1122/1.4731659</pub-id></citation></ref>
<ref id="B32">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mari</surname> <given-names>R.</given-names></name> <name><surname>Seto</surname> <given-names>R.</given-names></name> <name><surname>Morris</surname> <given-names>J. F.</given-names></name> <name><surname>Denn</surname> <given-names>M. M.</given-names></name></person-group> (<year>2014</year>). <article-title>Shear thickening, frictionless and frictional rheologies in non-Brownian suspensions</article-title>. <source>J. Rheol.</source> <volume>58</volume>, <fpage>1693</fpage>&#x02013;<lpage>1724</lpage>. <pub-id pub-id-type="doi">10.1122/1.4890747</pub-id></citation></ref>
<ref id="B33">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pan</surname> <given-names>Z.</given-names></name> <name><surname>de Cagny</surname> <given-names>H.</given-names></name> <name><surname>Weber</surname> <given-names>B.</given-names></name> <name><surname>Bonn</surname> <given-names>D.</given-names></name></person-group> (<year>2015</year>). <article-title>S-shaped flow curves of shear thickening suspensions: Direct observation of frictional rheology</article-title>. <source>Phys. Rev. E</source> <volume>92</volume>:<fpage>032202</fpage>. <pub-id pub-id-type="doi">10.1103/PhysRevE.92.032202</pub-id><pub-id pub-id-type="pmid">26465464</pub-id></citation></ref>
<ref id="B34">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Plachy</surname> <given-names>T.</given-names></name> <name><surname>Kutalkova</surname> <given-names>E.</given-names></name> <name><surname>Sedlacik</surname> <given-names>M.</given-names></name> <name><surname>Vesel</surname> <given-names>A.</given-names></name> <name><surname>Masar</surname> <given-names>M.</given-names></name> <name><surname>Kuritka</surname> <given-names>I.</given-names></name></person-group> (<year>2018</year>). <article-title>Impact of corrosion process of carbonyl iron particles on magnetorheological behavior of their suspensions</article-title>. <source>J. Indus. Eng. Chem.</source> <volume>66</volume>, <fpage>362</fpage>&#x02013;<lpage>369</lpage>. <pub-id pub-id-type="doi">10.1016/j.jiec.2018.06.002</pub-id></citation></ref>
<ref id="B35">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Scales</surname> <given-names>P. J.</given-names></name> <name><surname>Johnson</surname> <given-names>S. B.</given-names></name> <name><surname>Healy</surname> <given-names>T. W.</given-names></name> <name><surname>Kapur</surname> <given-names>P. C.</given-names></name></person-group> (<year>1998</year>). <article-title>Shear yield stress of partially flocculated colloidal suspensions</article-title>. <source>AIChE J.</source> <volume>44</volume>, <fpage>538</fpage>&#x02013;<lpage>544</lpage>. <pub-id pub-id-type="doi">10.1002/aic.690440305</pub-id></citation></ref>
<ref id="B36">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Suzuki</surname> <given-names>M.</given-names></name> <name><surname>Makino</surname> <given-names>K.</given-names></name> <name><surname>Yamada</surname> <given-names>M.</given-names></name> <name><surname>Linoya</surname> <given-names>K.</given-names></name></person-group> (<year>1981</year>). <article-title>A study on the coordination number in a system of randomly packed uniform sized spheres</article-title>. <source>Int. Chem. Eng.</source> <volume>21</volume>, <fpage>482</fpage>&#x02013;<lpage>488</lpage>.</citation></ref>
<ref id="B37">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tao</surname> <given-names>R.</given-names></name></person-group> (<year>2001</year>). <article-title>Super-strong magnetorheological fluids</article-title>. <source>J. Phys.</source> <volume>13</volume>, <fpage>R979</fpage>&#x02013;<lpage>R99</lpage>. <pub-id pub-id-type="doi">10.1088/0953-8984/13/50/202</pub-id></citation></ref>
<ref id="B38">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wyart</surname> <given-names>M.</given-names></name> <name><surname>Cates</surname> <given-names>M. E.</given-names></name></person-group> (<year>2014</year>). <article-title>Discontinuous shear thickening without inertia in dense non-Brownian suspensions</article-title>. <source>Phys. Rev. Lett.</source> <volume>112</volume>:<fpage>098302</fpage>. <pub-id pub-id-type="doi">10.1103/PhysRevLett.112.098302</pub-id><pub-id pub-id-type="pmid">24655284</pub-id></citation></ref>
</ref-list>
<fn-group>
<fn fn-type="financial-disclosure"><p><bold>Funding.</bold> Center National d&#x00027;Etudes Spatiales (CNES, the French Space Agency).</p></fn>
</fn-group>
</back>
</article>