Study of Nonlinear MHD Tribological Squeeze Film at Generalized Magnetic Reynolds Numbers Using DTM
Identifieur interne : 003B52 ( Ncbi/Merge ); précédent : 003B51; suivant : 003B53Study of Nonlinear MHD Tribological Squeeze Film at Generalized Magnetic Reynolds Numbers Using DTM
Auteurs : Mohammad Mehdi Rashidi [Iran, République populaire de Chine] ; Navid Freidoonimehr [Iran] ; Ebrahim Momoniat [Afrique du Sud] ; Behnam Rostami [Iran]Source :
- PLoS ONE [ 1932-6203 ] ; 2015.
Abstract
In the current article, a combination of the differential transform method (DTM) and Padé approximation method are implemented to solve a system of nonlinear differential equations modelling the flow of a Newtonian magnetic lubricant squeeze film with magnetic induction effects incorporated. Solutions for the transformed radial and tangential momentum as well as solutions for the radial and tangential induced magnetic field conservation equations are determined. The DTM-Padé combined method is observed to demonstrate excellent convergence, stability and versatility in simulating the magnetic squeeze film problem. The effects of involved parameters, i.e. squeeze Reynolds number (
Url:
DOI: 10.1371/journal.pone.0135004
PubMed: 26267247
PubMed Central: 4534046
Links toward previous steps (curation, corpus...)
- to stream Pmc, to step Corpus: 000317
- to stream Pmc, to step Curation: 000317
- to stream Pmc, to step Checkpoint: 000314
Links to Exploration step
PMC:4534046Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Study of Nonlinear MHD Tribological Squeeze Film at Generalized Magnetic Reynolds Numbers Using DTM</title><author><name sortKey="Rashidi, Mohammad Mehdi" sort="Rashidi, Mohammad Mehdi" uniqKey="Rashidi M" first="Mohammad Mehdi" last="Rashidi">Mohammad Mehdi Rashidi</name><affiliation wicri:level="1"><nlm:aff id="aff001"><addr-line>Mechanical Engineering Department, Engineering Faculty of Bu-Ali Sina University, Hamedan, Iran</addr-line></nlm:aff><country xml:lang="fr">Iran</country><wicri:regionArea>Mechanical Engineering Department, Engineering Faculty of Bu-Ali Sina University, Hamedan</wicri:regionArea><wicri:noRegion>Hamedan</wicri:noRegion></affiliation><affiliation wicri:level="1"><nlm:aff id="aff002"><addr-line>University of Michigan-Shanghai Jiao Tong University Joint Institute, Shanghai Jiao Tong University, Shanghai, Peoples Republic of China</addr-line></nlm:aff><country xml:lang="fr">République populaire de Chine</country><wicri:regionArea>University of Michigan-Shanghai Jiao Tong University Joint Institute, Shanghai Jiao Tong University, Shanghai</wicri:regionArea><wicri:noRegion>Shanghai</wicri:noRegion></affiliation></author><author><name sortKey="Freidoonimehr, Navid" sort="Freidoonimehr, Navid" uniqKey="Freidoonimehr N" first="Navid" last="Freidoonimehr">Navid Freidoonimehr</name><affiliation wicri:level="1"><nlm:aff id="aff003"><addr-line>Young Researchers & Elite Club, Hamedan Branch, Islamic Azad University, Hamedan, Iran</addr-line></nlm:aff><country xml:lang="fr">Iran</country><wicri:regionArea>Young Researchers & Elite Club, Hamedan Branch, Islamic Azad University, Hamedan</wicri:regionArea><wicri:noRegion>Hamedan</wicri:noRegion></affiliation></author><author><name sortKey="Momoniat, Ebrahim" sort="Momoniat, Ebrahim" uniqKey="Momoniat E" first="Ebrahim" last="Momoniat">Ebrahim Momoniat</name><affiliation wicri:level="4"><nlm:aff id="aff004"><addr-line>DST/NRF Centre of Excellence in the Mathematical and Statistical Sciences, School of Computer Science and Applied Mathematics, University of the Witwatersrand, Johannesburg, Private Bag 3, Wits 2050, South Africa</addr-line></nlm:aff><country xml:lang="fr">Afrique du Sud</country><wicri:regionArea>DST/NRF Centre of Excellence in the Mathematical and Statistical Sciences, School of Computer Science and Applied Mathematics, University of the Witwatersrand, Johannesburg, Private Bag 3, Wits 2050</wicri:regionArea><orgName type="university">Université du Witwatersrand</orgName><placeName><settlement type="city">Johannesburg</settlement><region type="region" nuts="2">Gauteng</region></placeName></affiliation></author><author><name sortKey="Rostami, Behnam" sort="Rostami, Behnam" uniqKey="Rostami B" first="Behnam" last="Rostami">Behnam Rostami</name><affiliation wicri:level="1"><nlm:aff id="aff003"><addr-line>Young Researchers & Elite Club, Hamedan Branch, Islamic Azad University, Hamedan, Iran</addr-line></nlm:aff><country xml:lang="fr">Iran</country><wicri:regionArea>Young Researchers & Elite Club, Hamedan Branch, Islamic Azad University, Hamedan</wicri:regionArea><wicri:noRegion>Hamedan</wicri:noRegion></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">26267247</idno><idno type="pmc">4534046</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4534046</idno><idno type="RBID">PMC:4534046</idno><idno type="doi">10.1371/journal.pone.0135004</idno><date when="2015">2015</date><idno type="wicri:Area/Pmc/Corpus">000317</idno><idno type="wicri:Area/Pmc/Curation">000317</idno><idno type="wicri:Area/Pmc/Checkpoint">000314</idno><idno type="wicri:Area/Ncbi/Merge">003B52</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Study of Nonlinear MHD Tribological Squeeze Film at Generalized Magnetic Reynolds Numbers Using DTM</title><author><name sortKey="Rashidi, Mohammad Mehdi" sort="Rashidi, Mohammad Mehdi" uniqKey="Rashidi M" first="Mohammad Mehdi" last="Rashidi">Mohammad Mehdi Rashidi</name><affiliation wicri:level="1"><nlm:aff id="aff001"><addr-line>Mechanical Engineering Department, Engineering Faculty of Bu-Ali Sina University, Hamedan, Iran</addr-line></nlm:aff><country xml:lang="fr">Iran</country><wicri:regionArea>Mechanical Engineering Department, Engineering Faculty of Bu-Ali Sina University, Hamedan</wicri:regionArea><wicri:noRegion>Hamedan</wicri:noRegion></affiliation><affiliation wicri:level="1"><nlm:aff id="aff002"><addr-line>University of Michigan-Shanghai Jiao Tong University Joint Institute, Shanghai Jiao Tong University, Shanghai, Peoples Republic of China</addr-line></nlm:aff><country xml:lang="fr">République populaire de Chine</country><wicri:regionArea>University of Michigan-Shanghai Jiao Tong University Joint Institute, Shanghai Jiao Tong University, Shanghai</wicri:regionArea><wicri:noRegion>Shanghai</wicri:noRegion></affiliation></author><author><name sortKey="Freidoonimehr, Navid" sort="Freidoonimehr, Navid" uniqKey="Freidoonimehr N" first="Navid" last="Freidoonimehr">Navid Freidoonimehr</name><affiliation wicri:level="1"><nlm:aff id="aff003"><addr-line>Young Researchers & Elite Club, Hamedan Branch, Islamic Azad University, Hamedan, Iran</addr-line></nlm:aff><country xml:lang="fr">Iran</country><wicri:regionArea>Young Researchers & Elite Club, Hamedan Branch, Islamic Azad University, Hamedan</wicri:regionArea><wicri:noRegion>Hamedan</wicri:noRegion></affiliation></author><author><name sortKey="Momoniat, Ebrahim" sort="Momoniat, Ebrahim" uniqKey="Momoniat E" first="Ebrahim" last="Momoniat">Ebrahim Momoniat</name><affiliation wicri:level="4"><nlm:aff id="aff004"><addr-line>DST/NRF Centre of Excellence in the Mathematical and Statistical Sciences, School of Computer Science and Applied Mathematics, University of the Witwatersrand, Johannesburg, Private Bag 3, Wits 2050, South Africa</addr-line></nlm:aff><country xml:lang="fr">Afrique du Sud</country><wicri:regionArea>DST/NRF Centre of Excellence in the Mathematical and Statistical Sciences, School of Computer Science and Applied Mathematics, University of the Witwatersrand, Johannesburg, Private Bag 3, Wits 2050</wicri:regionArea><orgName type="university">Université du Witwatersrand</orgName><placeName><settlement type="city">Johannesburg</settlement><region type="region" nuts="2">Gauteng</region></placeName></affiliation></author><author><name sortKey="Rostami, Behnam" sort="Rostami, Behnam" uniqKey="Rostami B" first="Behnam" last="Rostami">Behnam Rostami</name><affiliation wicri:level="1"><nlm:aff id="aff003"><addr-line>Young Researchers & Elite Club, Hamedan Branch, Islamic Azad University, Hamedan, Iran</addr-line></nlm:aff><country xml:lang="fr">Iran</country><wicri:regionArea>Young Researchers & Elite Club, Hamedan Branch, Islamic Azad University, Hamedan</wicri:regionArea><wicri:noRegion>Hamedan</wicri:noRegion></affiliation></author></analytic><series><title level="j">PLoS ONE</title><idno type="eISSN">1932-6203</idno><imprint><date when="2015">2015</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><p>In the current article, a combination of the differential transform method (DTM) and Padé approximation method are implemented to solve a system of nonlinear differential equations modelling the flow of a Newtonian magnetic lubricant squeeze film with magnetic induction effects incorporated. Solutions for the transformed radial and tangential momentum as well as solutions for the radial and tangential induced magnetic field conservation equations are determined. The DTM-Padé combined method is observed to demonstrate excellent convergence, stability and versatility in simulating the magnetic squeeze film problem. The effects of involved parameters, i.e. squeeze Reynolds number (<italic>N</italic><sub>1</sub>), dimensionless axial magnetic force strength parameter (<italic>N</italic><sub>2</sub>), dimensionless tangential magnetic force strength parameter (<italic>N</italic><sub>3</sub>), and magnetic Reynolds number (<italic>Re</italic><sub><italic>m</italic></sub>) are illustrated graphically and discussed in detail. Applications of the study include automotive magneto-rheological shock absorbers, novel aircraft landing gear systems and biological prosthetics.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Davidson, P A" uniqKey="Davidson P">P A Davidson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rashidi, M M" uniqKey="Rashidi M">M M Rashidi</name></author><author><name sortKey="Rostami, B" uniqKey="Rostami B">B Rostami</name></author><author><name sortKey="Freidoonimehr, N" uniqKey="Freidoonimehr N">N Freidoonimehr</name></author><author><name sortKey="Abbasbandy, S" uniqKey="Abbasbandy S">S Abbasbandy</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Khan, A" uniqKey="Khan A">A Khan</name></author><author><name sortKey="Khan, I" uniqKey="Khan I">I Khan</name></author><author><name sortKey="Ali, F" uniqKey="Ali F">F Ali</name></author><author><name sortKey="Ulhaq, S" uniqKey="Ulhaq S">S ulhaq</name></author><author><name sortKey="Shafie, S" uniqKey="Shafie S">S Shafie</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hussanan, A" uniqKey="Hussanan A">A Hussanan</name></author><author><name sortKey="Ismail, Z" uniqKey="Ismail Z">Z Ismail</name></author><author><name sortKey="Khan, I" uniqKey="Khan I">I Khan</name></author><author><name sortKey="Hussein, A" uniqKey="Hussein A">A Hussein</name></author><author><name sortKey="Shafie, S" uniqKey="Shafie S">S Shafie</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Samiulhaq" uniqKey="Samiulhaq">Samiulhaq</name></author><author><name sortKey="Ahmad, S" uniqKey="Ahmad S">S Ahmad</name></author><author><name sortKey="Vieru, D" uniqKey="Vieru D">D Vieru</name></author><author><name sortKey="Khan, I" uniqKey="Khan I">I Khan</name></author><author><name sortKey="Shafie, S" uniqKey="Shafie S">S Shafie</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Khan, A" uniqKey="Khan A">A Khan</name></author><author><name sortKey="Khan, I" uniqKey="Khan I">I Khan</name></author><author><name sortKey="Ali, F" uniqKey="Ali F">F Ali</name></author><author><name sortKey="Shafie, S" uniqKey="Shafie S">S Shafie</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Khalid, A" uniqKey="Khalid A">A Khalid</name></author><author><name sortKey="Khan, I" uniqKey="Khan I">I Khan</name></author><author><name sortKey="Khan, A" uniqKey="Khan A">A Khan</name></author><author><name sortKey="Shafie, S" uniqKey="Shafie S">S Shafie</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Park, E" uniqKey="Park E">E Park</name></author><author><name sortKey="Min, K W" uniqKey="Min K">K-W Min</name></author><author><name sortKey="Lee, S K" uniqKey="Lee S">S-K Lee</name></author><author><name sortKey="Lee, S H" uniqKey="Lee S">S-H Lee</name></author><author><name sortKey="Lee, H J" uniqKey="Lee H">H-J Lee</name></author><author><name sortKey="Moon, S J" uniqKey="Moon S">S-J Moon</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Anderson, C G" uniqKey="Anderson C">C G Anderson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chernyakova, Y M" uniqKey="Chernyakova Y">Y M Chernyakova</name></author><author><name sortKey="Pinchuk, L S" uniqKey="Pinchuk L">L S Pinchuk</name></author><author><name sortKey="Lobanovskii, L S" uniqKey="Lobanovskii L">L. S Lobanovskii</name></author><author><name sortKey="Drozd, E S" uniqKey="Drozd E">E S Drozd</name></author><author><name sortKey="Chizhik, S A" uniqKey="Chizhik S">S A Chizhik</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Liu, B" uniqKey="Liu B">B Liu</name></author><author><name sortKey="Li, W H" uniqKey="Li W">W H Li</name></author><author><name sortKey="Kosasih, P B" uniqKey="Kosasih P">P B Kosasih</name></author><author><name sortKey="Zhang, X Z" uniqKey="Zhang X">X Z Zhang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chandra, P" uniqKey="Chandra P">P Chandra</name></author><author><name sortKey="Sinha, P" uniqKey="Sinha P">P Sinha</name></author><author><name sortKey="Kumar, D" uniqKey="Kumar D">D Kumar</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Song, W L" uniqKey="Song W">W L Song</name></author><author><name sortKey="Choi, S B" uniqKey="Choi S">S-B Choi</name></author><author><name sortKey="Choi, J Y" uniqKey="Choi J">J Y Choi</name></author><author><name sortKey="Lee, C H" uniqKey="Lee C">C H Lee</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Duran, J D G" uniqKey="Duran J">J D G Durán</name></author><author><name sortKey="Iglesias, G R" uniqKey="Iglesias G">G R Iglesias</name></author><author><name sortKey="Delgado, A V" uniqKey="Delgado A">A V Delgado</name></author><author><name sortKey="Ruiz Mor N, L F" uniqKey="Ruiz Mor N L">L F Ruiz-Morón</name></author><author><name sortKey="Insa, J" uniqKey="Insa J">J Insa</name></author><author><name sortKey="Gonzalez Caballero, F" uniqKey="Gonzalez Caballero F">F González-Caballero</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Huang, W" uniqKey="Huang W">W Huang</name></author><author><name sortKey="Shen, C" uniqKey="Shen C">C Shen</name></author><author><name sortKey="Liao, S" uniqKey="Liao S">S Liao</name></author><author><name sortKey="Wang, X" uniqKey="Wang X">X Wang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Stolarski, T A" uniqKey="Stolarski T">T A Stolarski</name></author><author><name sortKey="Makida, Y" uniqKey="Makida Y">Y Makida</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Anwar, M I" uniqKey="Anwar M">M I Anwar</name></author><author><name sortKey="Rodkiewicz, C M" uniqKey="Rodkiewicz C">C M Rodkiewicz</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kargulewicz, M" uniqKey="Kargulewicz M">M Kargulewicz</name></author><author><name sortKey="Iordanoff, I" uniqKey="Iordanoff I">I Iordanoff</name></author><author><name sortKey="Marrero, V" uniqKey="Marrero V">V Marrero</name></author><author><name sortKey="Tichy, J" uniqKey="Tichy J">J Tichy</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zueco, J" uniqKey="Zueco J">J Zueco</name></author><author><name sortKey="Beg, O A" uniqKey="Beg O">O A Bég</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhu, G" uniqKey="Zhu G">G Zhu</name></author><author><name sortKey="Ingber, M S" uniqKey="Ingber M">M S Ingber</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gertzos, K P" uniqKey="Gertzos K">K P Gertzos</name></author><author><name sortKey="Nikolakopoulos, P G" uniqKey="Nikolakopoulos P">P G Nikolakopoulos</name></author><author><name sortKey="Papadopoulos, C A" uniqKey="Papadopoulos C">C A Papadopoulos</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Asadi, M" uniqKey="Asadi M">M Asadi</name></author><author><name sortKey="Rashidi, M M" uniqKey="Rashidi M">M M Rashidi</name></author><author><name sortKey="Beg, T A" uniqKey="Beg T">T A Bég</name></author><author><name sortKey="Beg, O A" uniqKey="Beg O">O A Bég</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Moghani, T" uniqKey="Moghani T">T Moghani</name></author><author><name sortKey="Butler, J P" uniqKey="Butler J">J P Butler</name></author><author><name sortKey="Lin, J L W" uniqKey="Lin J">J L-W Lin</name></author><author><name sortKey="Loring, S H" uniqKey="Loring S">S H Loring</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Elshekh, S S" uniqKey="Elshekh S">S S Elshekh</name></author><author><name sortKey="Abd Elhady, M K" uniqKey="Abd Elhady M">M K Abd Elhady</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gul, T" uniqKey="Gul T">T Gul</name></author><author><name sortKey="Islam, S" uniqKey="Islam S">S Islam</name></author><author><name sortKey="Shah, R A" uniqKey="Shah R">R A Shah</name></author><author><name sortKey="Khan, I" uniqKey="Khan I">I Khan</name></author><author><name sortKey="Khalid, A" uniqKey="Khalid A">A Khalid</name></author><author><name sortKey="Shafie, S" uniqKey="Shafie S">S Shafie</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gul, T" uniqKey="Gul T">T Gul</name></author><author><name sortKey="Islam, S" uniqKey="Islam S">S Islam</name></author><author><name sortKey="Shah, R A" uniqKey="Shah R">R A Shah</name></author><author><name sortKey="Khan, I" uniqKey="Khan I">I Khan</name></author><author><name sortKey="Shafie, S" uniqKey="Shafie S">S Shafie</name></author><author><name sortKey="Khan, M A" uniqKey="Khan M">M A Khan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Fatma, A" uniqKey="Fatma A">A Fatma</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Fatma, A" uniqKey="Fatma A">A Fatma</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Liao, S" uniqKey="Liao S">S Liao</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Liao, S" uniqKey="Liao S">S Liao</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Dehghan, M" uniqKey="Dehghan M">M Dehghan</name></author><author><name sortKey="Rahmani, Y" uniqKey="Rahmani Y">Y Rahmani</name></author><author><name sortKey="Ganji, D D" uniqKey="Ganji D">D D Ganji</name></author><author><name sortKey="Saedodin, S" uniqKey="Saedodin S">S Saedodin</name></author><author><name sortKey="Valipour, M S" uniqKey="Valipour M">M S Valipour</name></author><author><name sortKey="Rashidi, S" uniqKey="Rashidi S">S Rashidi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sheikholeslami, M" uniqKey="Sheikholeslami M">M Sheikholeslami</name></author><author><name sortKey="Ganji, D D" uniqKey="Ganji D">D D Ganji</name></author><author><name sortKey="Ashorynejad, H R" uniqKey="Ashorynejad H">H R Ashorynejad</name></author><author><name sortKey="Rokni, H B" uniqKey="Rokni H">H B Rokni</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rashidi, M M" uniqKey="Rashidi M">M M Rashidi</name></author><author><name sortKey="Erfani, E" uniqKey="Erfani E">E Erfani</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hamza, E A" uniqKey="Hamza E">E A Hamza</name></author><author><name sortKey="Macdonald, D A" uniqKey="Macdonald D">D A Macdonald</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ibrahim, F N" uniqKey="Ibrahim F">F N Ibrahim</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ibrahim, F N" uniqKey="Ibrahim F">F N Ibrahim</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hughes, W F" uniqKey="Hughes W">W F Hughes</name></author><author><name sortKey="Elco, R A" uniqKey="Elco R">R A Elco</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Baker, G A" uniqKey="Baker G">G A Baker</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Baker, G A" uniqKey="Baker G">G A Baker</name></author><author><name sortKey="Graves Morris, P R" uniqKey="Graves Morris P">P R Graves-Morris</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rashidi, M M" uniqKey="Rashidi M">M M Rashidi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rashidi, M M" uniqKey="Rashidi M">M M Rashidi</name></author><author><name sortKey="Mohimanian, Pour S A" uniqKey="Mohimanian P">Pour S A Mohimanian</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rashidi, M M" uniqKey="Rashidi M">M M Rashidi</name></author><author><name sortKey="Freidoonimehr, N" uniqKey="Freidoonimehr N">N Freidoonimehr</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hamza, E A" uniqKey="Hamza E">E A Hamza</name></author></analytic></biblStruct></listBibl></div1></back></TEI><pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">PLoS One</journal-id><journal-id journal-id-type="iso-abbrev">PLoS ONE</journal-id><journal-id journal-id-type="publisher-id">plos</journal-id><journal-id journal-id-type="pmc">plosone</journal-id><journal-title-group><journal-title>PLoS ONE</journal-title></journal-title-group><issn pub-type="epub">1932-6203</issn><publisher><publisher-name>Public Library of Science</publisher-name><publisher-loc>San Francisco, CA USA</publisher-loc></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">26267247</article-id><article-id pub-id-type="pmc">4534046</article-id><article-id pub-id-type="doi">10.1371/journal.pone.0135004</article-id><article-id pub-id-type="publisher-id">PONE-D-15-09032</article-id><article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group></article-categories><title-group><article-title>Study of Nonlinear MHD Tribological Squeeze Film at Generalized Magnetic Reynolds Numbers Using DTM</article-title><alt-title alt-title-type="running-head">Nonlinear MHD Tribological Squeeze Film</alt-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Rashidi</surname><given-names>Mohammad Mehdi</given-names></name><xref ref-type="aff" rid="aff001"><sup>1</sup></xref><xref ref-type="aff" rid="aff002"><sup>2</sup></xref><xref rid="cor001" ref-type="corresp">*</xref></contrib><contrib contrib-type="author"><name><surname>Freidoonimehr</surname><given-names>Navid</given-names></name><xref ref-type="aff" rid="aff003"><sup>3</sup></xref></contrib><contrib contrib-type="author"><name><surname>Momoniat</surname><given-names>Ebrahim</given-names></name><xref ref-type="aff" rid="aff004"><sup>4</sup></xref></contrib><contrib contrib-type="author"><name><surname>Rostami</surname><given-names>Behnam</given-names></name><xref ref-type="aff" rid="aff003"><sup>3</sup></xref></contrib></contrib-group><aff id="aff001"><label>1</label><addr-line>Mechanical Engineering Department, Engineering Faculty of Bu-Ali Sina University, Hamedan, Iran</addr-line></aff><aff id="aff002"><label>2</label><addr-line>University of Michigan-Shanghai Jiao Tong University Joint Institute, Shanghai Jiao Tong University, Shanghai, Peoples Republic of China</addr-line></aff><aff id="aff003"><label>3</label><addr-line>Young Researchers & Elite Club, Hamedan Branch, Islamic Azad University, Hamedan, Iran</addr-line></aff><aff id="aff004"><label>4</label><addr-line>DST/NRF Centre of Excellence in the Mathematical and Statistical Sciences, School of Computer Science and Applied Mathematics, University of the Witwatersrand, Johannesburg, Private Bag 3, Wits 2050, South Africa</addr-line></aff><contrib-group><contrib contrib-type="editor"><name><surname>Islam</surname><given-names>Saeed</given-names></name><role>Editor</role><xref ref-type="aff" rid="edit1"></xref></contrib></contrib-group><aff id="edit1"><addr-line>Abdul Wali Khan university Mardan Pakistan, PAKISTAN</addr-line></aff><author-notes><fn fn-type="conflict" id="coi001"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="con" id="contrib001"><p>Conceived and designed the experiments: MMR NF BR. Performed the experiments: NF BR. Analyzed the data: MMR EM. Contributed reagents/materials/analysis tools: MMR NF EM BR. Wrote the paper: MMR NF EM BR.</p></fn><corresp id="cor001">* E-mail: <email>mm_rashidi@sjtu.edu.cn</email></corresp></author-notes><pub-date pub-type="epub"><day>12</day><month>8</month><year>2015</year></pub-date><pub-date pub-type="collection"><year>2015</year></pub-date><volume>10</volume><issue>8</issue><elocation-id>e0135004</elocation-id><history><date date-type="received"><day>2</day><month>3</month><year>2015</year></date><date date-type="accepted"><day>16</day><month>7</month><year>2015</year></date></history><permissions><copyright-year>2015</copyright-year><copyright-holder>Rashidi et al</copyright-holder><license xlink:href="http://creativecommons.org/licenses/by/4.0/"><license-p>This is an open access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited</license-p></license></permissions><self-uri content-type="pdf" xlink:type="simple" xlink:href="pone.0135004.pdf"></self-uri><abstract><p>In the current article, a combination of the differential transform method (DTM) and Padé approximation method are implemented to solve a system of nonlinear differential equations modelling the flow of a Newtonian magnetic lubricant squeeze film with magnetic induction effects incorporated. Solutions for the transformed radial and tangential momentum as well as solutions for the radial and tangential induced magnetic field conservation equations are determined. The DTM-Padé combined method is observed to demonstrate excellent convergence, stability and versatility in simulating the magnetic squeeze film problem. The effects of involved parameters, i.e. squeeze Reynolds number (<italic>N</italic><sub>1</sub>), dimensionless axial magnetic force strength parameter (<italic>N</italic><sub>2</sub>), dimensionless tangential magnetic force strength parameter (<italic>N</italic><sub>3</sub>), and magnetic Reynolds number (<italic>Re</italic><sub><italic>m</italic></sub>) are illustrated graphically and discussed in detail. Applications of the study include automotive magneto-rheological shock absorbers, novel aircraft landing gear systems and biological prosthetics.</p></abstract><funding-group><funding-statement>EM acknowledges support from the National Research foundation of South Africa. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</funding-statement></funding-group><counts><fig-count count="13"></fig-count><table-count count="2"></table-count><page-count count="21"></page-count></counts><custom-meta-group><custom-meta id="data-availability"><meta-name>Data Availability</meta-name><meta-value>All relevant data are within the paper and its Supporting Information files.</meta-value></custom-meta></custom-meta-group></article-meta><notes><title>Data Availability</title><p>All relevant data are within the paper and its Supporting Information files.</p></notes></front><body><sec sec-type="intro" id="sec001"><title>Introduction</title><p>Understanding magneto-hydrodynamics (MHD) is strongly related to the comprehension of physical effects which take place in MHD. When a conductor moves into a magnetic field, electric current is induced in the conductor and creates its own magnetic field (Lenz’s law). Since the induced magnetic field tends to eliminate the original and external supported field, the magnetic field lines will be excluded from the conductor. Conversely, when the magnetic field influences the conductor to move it out of the field, the induced field amplifies the applied field. The net result of this process is that the lines of force appear to be dragged accompanied by the conductor. In this paper the conductor is the fluid with complex motions. To understand the second key effect which is dynamical we should know that when currents are induced by a motion of a conducting fluid through a magnetic field, a Lorentz force acts on the fluid and modifies its motion. In MHD, the motion modifies the field and vice versa. This makes the theory highly non-linear [<xref rid="pone.0135004.ref001" ref-type="bibr">1</xref>, <xref rid="pone.0135004.ref002" ref-type="bibr">2</xref>].</p><p>In recent decades, researchers have performed several studies in the fields of the MHD applications’. Khan <italic>et al</italic>. [<xref rid="pone.0135004.ref003" ref-type="bibr">3</xref>] investigated the effects of an arbitrary wall shear stress on unsteady MHD flow of a Newtonian fluid with conjugate effects of heat and mass transfer using the Laplace transform technique. Hussanan <italic>et al</italic>. [<xref rid="pone.0135004.ref004" ref-type="bibr">4</xref>] analysed the unsteady boundary layer MHD free convection flow past an oscillating vertical plate embedded in a porous medium with constant mass diffusion and Newtonian heating condition using the Laplace transform technique. Samiulhaq <italic>et al</italic>. [<xref rid="pone.0135004.ref005" ref-type="bibr">5</xref>] studied the magnetic field influence on unsteady free convection flow of a second grade fluid near an infinite vertical flat plate with ramped wall temperature embedded in a porous medium. In another study, Khan <italic>et al</italic>. [<xref rid="pone.0135004.ref006" ref-type="bibr">6</xref>] displayed the effects of an arbitrary wall shear stress on unsteady MHD flow of a Newtonian fluid with conjugate effects of heat and mass transfer. Further, Khalid <italic>et al</italic>. [<xref rid="pone.0135004.ref007" ref-type="bibr">7</xref>] illustrated the unsteady MHD free flow of a Casson fluid past an oscillating vertical plate with constant wall temperature.</p><p>Magneto-hydrodynamic lubrication is a type of “smart tribology” which has found increasing applications in diverse branches of engineering in recent years. These include seismic magneto-rheological (MR) shock dampers [<xref rid="pone.0135004.ref008" ref-type="bibr">8</xref>], magnetic-repulsion enhanced hydrostatic bearings for offshore wave energy conversion devices [<xref rid="pone.0135004.ref009" ref-type="bibr">9</xref>], and biomedical systems [<xref rid="pone.0135004.ref010" ref-type="bibr">10</xref>, <xref rid="pone.0135004.ref011" ref-type="bibr">11</xref>]. MHD lubricants respond to the application of magnetic fields and have been presented to enhance load-carrying capacities, reduce wear and achieve more uniform pressure distributions. A tremendous variety of such lubricants has been developed including ferrofluids, magnetic particle-based suspensions, electrically-conducting biopolymers and yield stress magnetic fluids. In parallel with practical and manufacturing developments, there has been a rich contribution from engineering scientists engaged in mathematical and experimental simulations of the behaviour of such fluids in many complex tribological configurations.</p><p>Chandra <italic>et al</italic>. [<xref rid="pone.0135004.ref012" ref-type="bibr">12</xref>] studied electromagnetic lubrication in various journal bearings with cavitation boundary conditions and for regimes where the magnetization vector is oblique to the magnetic field vector, showing that better contact performance is achieved compared with non-magnetic lubricants. Song <italic>et al</italic>. [<xref rid="pone.0135004.ref013" ref-type="bibr">13</xref>] analysed wear and friction characteristics of a magneto-rheological fluid under different magnetic fields, by employing a pin-on-disc tribometer, and showed that MR fluid exhibits improved lubrication characteristics. They also observed that the key wear mechanism for steel and brass specimens was abrasive wear by asperities and MR particles on the worn surfaces, whereas a mixed wear mechanism that included adhesive wear and abrasive wear was observed for the aluminium specimen studied. Durán <italic>et al</italic>. [<xref rid="pone.0135004.ref014" ref-type="bibr">14</xref>] presented a novel formulation of a stable magnetic fluid to show that the yield stress is elevated by several orders of magnitude when the magnetic field strength reaches several hundred microTesla, and furthermore observed that excellent damping of forced oscillations is achievable in automotive magneto-fluid bearings. Huang <italic>et al</italic>. [<xref rid="pone.0135004.ref015" ref-type="bibr">15</xref>] investigated ferrofluid magnetic tribology for lubricants comprising stable colloidal systems consisting of single-domain magnetic particles with a diameter of approximately 10 nm coated with surfactants and dispersed in a carrier liquid. By applying an external magnetic field, they demonstrated that ferrofluid lubricants may be orientated and positioned at optimized locations. They additionally showed that the load capacity of a Fe<sub>3</sub>O<sub>4</sub>-based ferrofluid lubricant film may be significantly boosted with appropriate magnetic field and that these liquids achieve a good friction-reduction performance in the presence of an external magnetic field compared with the carrier liquid with markedly enhanced lifetimes. Stolarsky and Makida [<xref rid="pone.0135004.ref016" ref-type="bibr">16</xref>] conducted experiments on the effect of permanent magnetic fields on the wear of lubricated sliding contact operating at short stroke and high frequency, observing that horizontal magnetic field strongly influences contact performance. They also noted that magnetic field increases the abrasive action by wear particles and allows reduction in wear of the plate specimen. An especially significant regime in magneto-hydrodynamic tribology is the squeeze film. This has attracted considerable attention as it can be simulated using classical methods of elasto-hydrodynamics (EHD). Anwar and Rodkiewicz [<xref rid="pone.0135004.ref017" ref-type="bibr">17</xref>] examined computationally the MHD squeeze lubrication of a slider bearing system, including inertia effects and considering low Hartmann numbers. They found that inertial terms have a reduced influence with greater Hartmann number and that a non-uniform magnetic field achieves noticeably greater load capacity than uniform magnetic fields.</p><p>In recent years many sophisticated numerical and so-called semi-numerical/analytical procedures have been implemented to solve boundary value problems arising in magnetic tribology and also nonlinear squeeze film flows. Kargulewicz <italic>et al</italic>. [<xref rid="pone.0135004.ref018" ref-type="bibr">18</xref>] developed a discrete element algorithm to optimize aircraft ejector seat applications. Zueco and Bég [<xref rid="pone.0135004.ref019" ref-type="bibr">19</xref>] applied the electro-thermal network simulation code to study magneto-elastic hydrodynamic lubrication between rotating disks at generalized magnetic Reynolds numbers, as a model of conceptual spacecraft landing gear systems for Mars NASA missions. Zhu and Ingber [<xref rid="pone.0135004.ref020" ref-type="bibr">20</xref>] utilized a traction-modified boundary element method (BEM) to study Newtonian squeeze films between spherical bodies in locomotive gear systems. Gertzos <italic>et al</italic>. [<xref rid="pone.0135004.ref021" ref-type="bibr">21</xref>] studied the performance characteristics of a hydrodynamic journal bearing lubricated with either magneto-rheological Bingham or electro-rheological Bingham fluids using the commercial CFD software, FLUENT software with a “dynamic meshing” technique. Bég <italic>et al</italic>. [<xref rid="pone.0135004.ref022" ref-type="bibr">22</xref>] analyzed the magneto-hydrodynamic squeezing flow of a microstructural fluid in a porous media biological bearing with the Liao homotopy analysis method (HAM), observing that micro-rotation of lubricant micro-elements is strongly influenced by Hartmann number and medium permeability, and that response time is also enhanced with magnetic field. Moghani <italic>et al</italic>. [<xref rid="pone.0135004.ref023" ref-type="bibr">23</xref>] used a hybrid fluid-solid meshing procedure in the ADINA commercial finite element code to study squeezing lubrication of soft biomaterials. From the above squeezing hydrodynamics studies which have considered magnetic fields, with the exception of Zueco and Bég [<xref rid="pone.0135004.ref019" ref-type="bibr">19</xref>] have generally neglected magnetic induction effects. When magnetic Reynolds number is sizeable, an induced magnetic field is also generated in the flow and a separate magnetic field conservation equation is required. Several researchers have studied magnetic induction effects. Elshekh and Elhady [<xref rid="pone.0135004.ref024" ref-type="bibr">24</xref>] investigated magnetic squeeze film flow between co-rotating disks with induced magnetic field effects, computing the response of radial and azimuthal magnetic fields to squeezing rates and relative disk rotation, although only for a single value of Batchelor number. Gul <italic>et al</italic>. [<xref rid="pone.0135004.ref025" ref-type="bibr">25</xref>] demonstrated the problem of thin film layer flowing on a vertical oscillating belt via two analytical techniques namely Adomian Decomposition Method (ADM) and Optimal Homotopy Asymptotic Method (OHAM). In another study, Gul <italic>et al</italic>. [<xref rid="pone.0135004.ref026" ref-type="bibr">26</xref>] performed an analysis to study the unsteady thin film flow of a second grade fluid over a vertical oscillating belt.</p><p>Nonlinear differential equations are employed to describe some of physical systems. Concurrent with the development of computers, rising use of analytical methods can be observed in comparison with numerical methods. Despite all the benefits, there are a lot of cons for the numerical methods such as the inability to apply infinite boundary condition, etc. There are a lot of analytical methods such as DTM [<xref rid="pone.0135004.ref027" ref-type="bibr">27</xref>, <xref rid="pone.0135004.ref028" ref-type="bibr">28</xref>], HAM [<xref rid="pone.0135004.ref029" ref-type="bibr">29</xref>, <xref rid="pone.0135004.ref030" ref-type="bibr">30</xref>], HPM [<xref rid="pone.0135004.ref031" ref-type="bibr">31</xref>], and ADM [<xref rid="pone.0135004.ref032" ref-type="bibr">32</xref>] applied to solve nonlinear equations. The main advantage of these methods applied to nonlinear differential equations is that no linearization or discretization needs to take place [<xref rid="pone.0135004.ref033" ref-type="bibr">33</xref>]. In the present article we employ DTM-Padé method to analyse two nonlinear magneto-hydrodynamic squeeze film boundary value problems. The present DTM-Padé code is also benchmarked with the numerical method based on shooting technique, illustrating excellent correlation. Excellent convergence and stability characteristics are also observed for the DTM-Padé code. The present simulations find applications in novel aircraft landing gear systems exploiting smart magnetic fluids.</p><p>The paper is divided up as follows: in section 2 we derive the mathematical model we will be investigating in this paper. In section 3 we implement the DTM-Padé method to solve the resulting system of nonlinear differential equations. Results are discussed in section 4. Concluding remarks are presented in section 5.</p></sec><sec id="sec002"><title>Mathematical Model</title><p>In this problem, we assume the axisymmetric flow in a polar coordinate system (<italic>r</italic>, <italic>θ</italic>, <italic>z</italic>) of a thin Newtonian, hydro-magnetic lubricant fluid squeeze film between two disks placed parallel to each other and each rotating at constant velocity in its own plane. The components of the flow velocity (<italic>u</italic>, <italic>v</italic>, <italic>w</italic>) are in the directions of increasing (<italic>r</italic>, <italic>θ</italic>, <italic>z</italic>), respectively. The disks are separated by a distance <italic>D</italic>(1 − <italic>αt</italic>)<sup>1/2</sup> at time <italic>t</italic>, where <italic>D</italic> is a representative length equivalent to the disk separation at <italic>t</italic> = 0 and <italic>t</italic> denotes time. The coordinate system and the physical regime of the problem are shown in <xref rid="pone.0135004.g001" ref-type="fig">Fig 1</xref>.</p><fig id="pone.0135004.g001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0135004.g001</object-id><label>Fig 1</label><caption><title>Hydromagnetic rotating squeeze film lubrication with magnetic induction.</title></caption><graphic xlink:href="pone.0135004.g001"></graphic></fig><p>The upper disk can moves downwards at velocity <italic>d</italic>(<italic>d</italic>(<italic>t</italic>))/<italic>dt</italic> (where <italic>d</italic>(<italic>t</italic>) is separation of the disks at time <italic>t</italic>), towards the constrained lower disk. In the other words, the lower disk is prohibited from moving in the axial direction (along the z-axis). The applied magnetic field (<italic>H</italic>) has two components; include an azimuthal (tangential) component (<italic>H</italic><sub><italic>θ</italic></sub>) and an axial component (<italic>H</italic><sub><italic>z</italic></sub>) which affect the upper disk. These parameters are defined as <italic>H</italic><sub><italic>θ</italic></sub> = <italic>r N</italic><sub>0</sub>/(<italic>μ</italic><sub>2</sub>(1 − <italic>αt</italic>)) and <italic>H</italic><sub><italic>z</italic></sub> = −<italic>α M</italic><sub>0</sub>/(<italic>μ</italic><sub>1</sub>(1 − <italic>αt</italic>)<sup>1/2</sup>), where <italic>M</italic><sub>0</sub> and <italic>N</italic><sub>0</sub> are the magnetic field quantities introduced to render <italic>H</italic><sub><italic>θ</italic></sub> and <italic>H</italic><sub><italic>z</italic></sub>, dimensionless, <italic>μ</italic><sub>2</sub> and <italic>μ</italic><sub>1</sub> are the magnetic permeability’s of the squeeze film and the medium external to the disks, respectively. For liquid metals <italic>μ</italic><sub>2</sub> = <italic>μ</italic><sub>0</sub>, where <italic>μ</italic><sub>0</sub> is the permeability of free space. Following the experimental study, <italic>H</italic><sub><italic>θ</italic></sub> and <italic>H</italic><sub><italic>z</italic></sub> are assumed to be zero on the lower disk [<xref rid="pone.0135004.ref024" ref-type="bibr">24</xref>]. This applied magnetic field (<italic>H</italic>) generates an induced magnetic field <italic>B</italic> with components (<italic>B</italic><sub><italic>r</italic></sub>, <italic>B</italic><sub><italic>θ</italic></sub>, <italic>B</italic><sub><italic>z</italic></sub>) in the squeeze film, between the disks. By neglecting the convective acceleration components in the Navier–Stokes equations and considering the above assumptions, the governing equations of the hydro-magnetic squeeze film regime for the momentum and magnetic field equations in a (<italic>r</italic>, <italic>θ</italic>, <italic>z</italic>) coordinate system are [<xref rid="pone.0135004.ref024" ref-type="bibr">24</xref>, <xref rid="pone.0135004.ref034" ref-type="bibr">34</xref>, <xref rid="pone.0135004.ref035" ref-type="bibr">35</xref>]:
<disp-formula id="pone.0135004.e001"><alternatives><graphic xlink:href="pone.0135004.e001.jpg" id="pone.0135004.e001g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M1"><mml:mrow><mml:mfrac><mml:mrow><mml:mo>∂</mml:mo><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mi>ν</mml:mi><mml:mfrac><mml:mrow><mml:msup><mml:mo>∂</mml:mo><mml:mn>2</mml:mn></mml:msup><mml:mi>u</mml:mi></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:msup><mml:mi>z</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mfrac><mml:mo>−</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mi>ρ</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>μ</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:mfrac><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mn>2</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:mi>r</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>M</mml:mi><mml:mfrac><mml:mrow><mml:msup><mml:mo>∂</mml:mo><mml:mn>2</mml:mn></mml:msup><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:msup><mml:mi>z</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(1)</label></disp-formula><disp-formula id="pone.0135004.e002"><alternatives><graphic xlink:href="pone.0135004.e002.jpg" id="pone.0135004.e002g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M2"><mml:mrow><mml:mfrac><mml:mrow><mml:mo>∂</mml:mo><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mi>ν</mml:mi><mml:mfrac><mml:mrow><mml:msup><mml:mo>∂</mml:mo><mml:mn>2</mml:mn></mml:msup><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:msup><mml:mi>z</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mfrac><mml:mo>−</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mi>ρ</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>μ</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:mfrac><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mn>2</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:mi>r</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>M</mml:mi><mml:mfrac><mml:mrow><mml:msup><mml:mo>∂</mml:mo><mml:mn>2</mml:mn></mml:msup><mml:mi>N</mml:mi></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:msup><mml:mi>z</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(2)</label></disp-formula><disp-formula id="pone.0135004.e003"><alternatives><graphic xlink:href="pone.0135004.e003.jpg" id="pone.0135004.e003g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M3"><mml:mrow><mml:mn>0</mml:mn><mml:mo>=</mml:mo><mml:mi>ν</mml:mi><mml:mfrac><mml:mrow><mml:msup><mml:mo>∂</mml:mo><mml:mn>2</mml:mn></mml:msup><mml:mi>w</mml:mi></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:msup><mml:mi>z</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mfrac><mml:mo>−</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mi>ρ</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>μ</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:mfrac><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msup><mml:mi>r</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mi>N</mml:mi><mml:mfrac><mml:mrow><mml:mo>∂</mml:mo><mml:mi>N</mml:mi></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:mfrac><mml:mo>−</mml:mo><mml:msup><mml:mi>r</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mfrac><mml:mrow><mml:msup><mml:mo>∂</mml:mo><mml:mn>3</mml:mn></mml:msup><mml:mi>N</mml:mi></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:msup><mml:mi>z</mml:mi><mml:mn>3</mml:mn></mml:msup></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(3)</label></disp-formula><disp-formula id="pone.0135004.e004"><alternatives><graphic xlink:href="pone.0135004.e004.jpg" id="pone.0135004.e004g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M4"><mml:mrow><mml:mfrac><mml:mrow><mml:mo>∂</mml:mo><mml:msub><mml:mi>B</mml:mi><mml:mi>r</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mfrac><mml:mrow><mml:msup><mml:mo>∂</mml:mo><mml:mn>3</mml:mn></mml:msup><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:msup><mml:mi>z</mml:mi><mml:mn>3</mml:mn></mml:msup></mml:mrow></mml:mfrac><mml:mo>+</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mi>σ</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>μ</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:mfrac><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:msup><mml:mo>∂</mml:mo><mml:mn>2</mml:mn></mml:msup><mml:msub><mml:mi>B</mml:mi><mml:mi>r</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:msup><mml:mi>z</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(4)</label></disp-formula><disp-formula id="pone.0135004.e005"><alternatives><graphic xlink:href="pone.0135004.e005.jpg" id="pone.0135004.e005g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M5"><mml:mrow><mml:mfrac><mml:mrow><mml:mo>∂</mml:mo><mml:msub><mml:mi>B</mml:mi><mml:mi>θ</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mfrac><mml:mrow><mml:msup><mml:mo>∂</mml:mo><mml:mn>2</mml:mn></mml:msup><mml:mi>N</mml:mi></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:msup><mml:mi>z</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mfrac><mml:mo>+</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mi>σ</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>μ</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:mfrac><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:msup><mml:mo>∂</mml:mo><mml:mn>2</mml:mn></mml:msup><mml:msub><mml:mi>B</mml:mi><mml:mi>θ</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:msup><mml:mi>z</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(5)</label></disp-formula><disp-formula id="pone.0135004.e006"><alternatives><graphic xlink:href="pone.0135004.e006.jpg" id="pone.0135004.e006g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M6"><mml:mrow><mml:mfrac><mml:mrow><mml:mo>∂</mml:mo><mml:msub><mml:mi>B</mml:mi><mml:mi>z</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mi>σ</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>μ</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:mfrac><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:msup><mml:mo>∂</mml:mo><mml:mn>2</mml:mn></mml:msup><mml:msub><mml:mi>B</mml:mi><mml:mi>z</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:msup><mml:mi>z</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(6)</label></disp-formula>where <italic>ν</italic> is the kinematic fluid viscosity, <italic>ρ</italic> is fluid density, <italic>μ</italic><sub>2</sub> is magnetic permeability of the medium between the two disks (squeeze film regime) and <italic>σ</italic> is the electrical conductivity of fluid (squeeze film).</p><p>It should be noticed that the Eqs (<xref rid="pone.0135004.e004" ref-type="disp-formula">4</xref>)–(<xref rid="pone.0135004.e006" ref-type="disp-formula">6</xref>) are the components of magnetic field equation <inline-formula id="pone.0135004.e007"><alternatives><graphic xlink:href="pone.0135004.e007.jpg" id="pone.0135004.e007g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M7"><mml:mrow><mml:mfrac><mml:mrow><mml:mo>∂</mml:mo><mml:mspace width="0.15em"></mml:mspace><mml:mi mathvariant="bold">B</mml:mi></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mi mathvariant="normal">∇</mml:mi><mml:mo>×</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi mathvariant="normal">∇</mml:mi><mml:mo>×</mml:mo><mml:mi mathvariant="bold">B</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mi>σ</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>μ</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:mfrac><mml:msup><mml:mi mathvariant="normal">∇</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mi mathvariant="bold">B</mml:mi></mml:mrow></mml:math></alternatives></inline-formula> in the directions of polar coordinate system, where <bold>B</bold> = <italic>μ</italic><sub>2</sub><bold>H</bold> [<xref rid="pone.0135004.ref024" ref-type="bibr">24</xref>, <xref rid="pone.0135004.ref036" ref-type="bibr">36</xref>]. The governing equations can be reduced from a system of partial differential equations (PDEs) to dimensionless, coupled and nonlinear ordinary differential equations (ODEs) by introducing the following transformations:
<disp-formula id="pone.0135004.e008"><alternatives><graphic xlink:href="pone.0135004.e008.jpg" id="pone.0135004.e008g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M8"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mi>η</mml:mi><mml:mo>=</mml:mo><mml:mfrac><mml:mi>z</mml:mi><mml:mrow><mml:mi>d</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mfrac><mml:mspace width="12em"></mml:mspace><mml:mi>d</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>D</mml:mi><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>α</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>ξ</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mrow><mml:mn>1</mml:mn><mml:mo>/</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:mrow></mml:msup></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mi>u</mml:mi><mml:mo>=</mml:mo><mml:mi>r</mml:mi><mml:mfrac><mml:mrow><mml:mo>∂</mml:mo><mml:mi>F</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.15em"></mml:mspace><mml:mi>t</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>α</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>r</mml:mi><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>η</mml:mi></mml:mrow></mml:mfrac></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>α</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>ξ</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mfrac><mml:mspace width="6em"></mml:mspace><mml:mi>v</mml:mi><mml:mo>=</mml:mo><mml:mi>r</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>G</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.15em"></mml:mspace><mml:mi>t</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>r</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi mathvariant="normal">Ω</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:mi>g</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>η</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>α</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>ξ</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mfrac></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mi>w</mml:mi><mml:mo>=</mml:mo><mml:mo>−</mml:mo><mml:mn>2</mml:mn><mml:mi>F</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.15em"></mml:mspace><mml:mi>t</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mo>−</mml:mo><mml:mfrac><mml:mrow><mml:mi>α</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>D</mml:mi><mml:mspace width="0.1em"></mml:mspace><mml:mi>f</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mrow><mml:mo>(</mml:mo><mml:mi>η</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>α</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>ξ</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mn>1</mml:mn><mml:mo>/</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:mrow></mml:msup></mml:mrow></mml:mfrac><mml:mspace width="2.8em"></mml:mspace><mml:msub><mml:mi>B</mml:mi><mml:mi>r</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mi>r</mml:mi><mml:mfrac><mml:mrow><mml:mo>∂</mml:mo><mml:mi>M</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.15em"></mml:mspace><mml:mi>t</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>α</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>r</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>M</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>m</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>η</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>η</mml:mi></mml:mrow></mml:mfrac></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:mi>D</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>α</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>ξ</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mfrac></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>B</mml:mi><mml:mi>θ</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mi>r</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>N</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.15em"></mml:mspace><mml:mi>t</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>r</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>N</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:mi>n</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>η</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>α</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>ξ</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mfrac><mml:mspace width="4.8em"></mml:mspace><mml:msub><mml:mi>B</mml:mi><mml:mi>z</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mo>−</mml:mo><mml:mn>2</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:mi>M</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.15em"></mml:mspace><mml:mi>t</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mo>−</mml:mo><mml:mfrac><mml:mrow><mml:mi>α</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>M</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:mi>m</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>η</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>α</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>ξ</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mn>1</mml:mn><mml:mo>/</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:mrow></mml:msup></mml:mrow></mml:mfrac></mml:mtd></mml:mtr></mml:mtable></mml:math></alternatives><label>(7)</label></disp-formula></p><p>Substituting <xref rid="pone.0135004.e008" ref-type="disp-formula">Eq 7</xref> in Eqs (<xref rid="pone.0135004.e001" ref-type="disp-formula">1</xref>)–(<xref rid="pone.0135004.e006" ref-type="disp-formula">6</xref>), we obtain the following dimensionless system of highly nonlinear and coupled ordinary differential equations in terms of a single independent space variable (<italic>η</italic>):
<disp-formula id="pone.0135004.e009"><alternatives><graphic xlink:href="pone.0135004.e009.jpg" id="pone.0135004.e009g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M9"><mml:mrow><mml:mfrac><mml:mrow><mml:msup><mml:mi>d</mml:mi><mml:mn>4</mml:mn></mml:msup><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:msup><mml:mi>η</mml:mi><mml:mn>4</mml:mn></mml:msup></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mn>3</mml:mn><mml:mfrac><mml:mrow><mml:msup><mml:mi>d</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:msup><mml:mi>η</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mfrac><mml:mo>+</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>η</mml:mi><mml:mo>−</mml:mo><mml:mn>2</mml:mn><mml:mi>f</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mfrac><mml:mrow><mml:msup><mml:mi>d</mml:mi><mml:mn>3</mml:mn></mml:msup><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:msup><mml:mi>η</mml:mi><mml:mn>3</mml:mn></mml:msup></mml:mrow></mml:mfrac><mml:mo>−</mml:mo><mml:mn>2</mml:mn><mml:msubsup><mml:mi>N</mml:mi><mml:mn>2</mml:mn><mml:mrow><mml:mspace width="0.25em"></mml:mspace><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mi>m</mml:mi><mml:mfrac><mml:mrow><mml:msup><mml:mi>d</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:msup><mml:mi>η</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(8)</label></disp-formula><disp-formula id="pone.0135004.e010"><alternatives><graphic xlink:href="pone.0135004.e010.jpg" id="pone.0135004.e010g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M10"><mml:mrow><mml:mfrac><mml:mrow><mml:msup><mml:mi>d</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mi>g</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:msup><mml:mi>η</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mn>2</mml:mn><mml:mi>g</mml:mi><mml:mo>+</mml:mo><mml:mi>η</mml:mi><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>g</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>η</mml:mi></mml:mrow></mml:mfrac><mml:mo>+</mml:mo><mml:mn>2</mml:mn><mml:mi>g</mml:mi><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>η</mml:mi></mml:mrow></mml:mfrac><mml:mo>−</mml:mo><mml:mn>2</mml:mn><mml:mi>f</mml:mi><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>g</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>η</mml:mi></mml:mrow></mml:mfrac><mml:mo>+</mml:mo><mml:mn>2</mml:mn><mml:msub><mml:mi>N</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:msub><mml:mi>N</mml:mi><mml:mn>3</mml:mn></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>m</mml:mi><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>η</mml:mi></mml:mrow></mml:mfrac><mml:mo>−</mml:mo><mml:mi>n</mml:mi><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>η</mml:mi></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(9)</label></disp-formula><disp-formula id="pone.0135004.e011"><alternatives><graphic xlink:href="pone.0135004.e011.jpg" id="pone.0135004.e011g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M11"><mml:mrow><mml:mfrac><mml:mrow><mml:msup><mml:mi>d</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:msup><mml:mi>η</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mi>R</mml:mi><mml:msub><mml:mi>e</mml:mi><mml:mi>m</mml:mi></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>m</mml:mi><mml:mo>+</mml:mo><mml:mi>η</mml:mi><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>η</mml:mi></mml:mrow></mml:mfrac><mml:mo>−</mml:mo><mml:mn>2</mml:mn><mml:mi>f</mml:mi><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>η</mml:mi></mml:mrow></mml:mfrac><mml:mo>+</mml:mo><mml:mn>2</mml:mn><mml:mi>m</mml:mi><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>η</mml:mi></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(10)</label></disp-formula><disp-formula id="pone.0135004.e012"><alternatives><graphic xlink:href="pone.0135004.e012.jpg" id="pone.0135004.e012g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M12"><mml:mrow><mml:mfrac><mml:mrow><mml:msup><mml:mi>d</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:msup><mml:mi>η</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mi>R</mml:mi><mml:msub><mml:mi>e</mml:mi><mml:mi>m</mml:mi></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mn>2</mml:mn><mml:mi>n</mml:mi><mml:mo>+</mml:mo><mml:mi>η</mml:mi><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>η</mml:mi></mml:mrow></mml:mfrac><mml:mo>−</mml:mo><mml:mn>2</mml:mn><mml:mi>f</mml:mi><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>η</mml:mi></mml:mrow></mml:mfrac><mml:mo>+</mml:mo><mml:mn>2</mml:mn><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>m</mml:mi><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>g</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>η</mml:mi></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(11)</label></disp-formula>where <italic>η</italic> is dimensionless z-coordinate, <italic>f</italic> is dimensionless axial velocity, <italic>g</italic> is dimensionless tangential velocity, <italic>m</italic> is dimensionless axial induced magnetic field component, <italic>n</italic> is dimensionless tangential induced magnetic field component, <italic>N</italic><sub>1</sub>(=<italic>α D</italic><sup>2</sup>/2<italic>ν</italic>) is squeeze Reynolds number (based on the speed of approach of the two disks), represents the ratio between the normal (axial) velocity of the upper disk and kinematic viscosity, <inline-formula id="pone.0135004.e013"><alternatives><graphic xlink:href="pone.0135004.e013.jpg" id="pone.0135004.e013g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M13"><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow><mml:mo>/</mml:mo><mml:mrow><mml:mi>D</mml:mi><mml:msqrt><mml:mrow><mml:msub><mml:mi>μ</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:mi>ρ</mml:mi></mml:mrow></mml:msqrt></mml:mrow></mml:mrow></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:math></alternatives></inline-formula> is dimensionless parameter based on the magnetic force in the axial direction, <inline-formula id="pone.0135004.e014"><alternatives><graphic xlink:href="pone.0135004.e014.jpg" id="pone.0135004.e014g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M14"><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mn>3</mml:mn></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow><mml:mo>/</mml:mo><mml:mrow><mml:msub><mml:mi mathvariant="normal">Ω</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:msqrt><mml:mrow><mml:msub><mml:mi>μ</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:mi>ρ</mml:mi></mml:mrow></mml:msqrt></mml:mrow></mml:mrow></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:math></alternatives></inline-formula> is dimensionless parameter based on magnetic force strength in the azimuthal (tangential) direction, <italic>Re</italic><sub><italic>m</italic></sub>(=<italic>N</italic><sub>1</sub><italic>Bt</italic>) is magnetic Reynolds number and <italic>Bt</italic>(=<italic>σ μ</italic><sub>2</sub><italic>ν</italic>) is the Batchelor number. As it is obvious from Eqs (<xref rid="pone.0135004.e009" ref-type="disp-formula">8</xref>)–(<xref rid="pone.0135004.e012" ref-type="disp-formula">11</xref>) the variation in the magnetic Reynolds number can be varied by keeping <italic>N</italic><sub>1</sub> invariant and altering Batchelor number (<italic>Bt</italic>) or vice versa. According to the Hughes and Elco [<xref rid="pone.0135004.ref037" ref-type="bibr">37</xref>], it is considered that both disks are ideal (perfect) conductors. Electrical forces are much smaller than the magnetic forces and consequently are ignored in the present problem. The boundary conditions for the magnetic induction <bold>B</bold> follow from the fact that the normal component of <bold>B</bold> and the tangential component of <bold>H</bold> are continuous through the two disks. The transformed boundary conditions become:
<disp-formula id="pone.0135004.e015"><alternatives><graphic xlink:href="pone.0135004.e015.jpg" id="pone.0135004.e015g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M15"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mn>0</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mn>0</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>η</mml:mi></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mspace width="0.5em"></mml:mspace><mml:mi>g</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mn>0</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mspace width="0.5em"></mml:mspace><mml:mi>m</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mn>0</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>n</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mn>0</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mspace width="4.3em"></mml:mspace><mml:mtext>at</mml:mtext><mml:mspace width="1.5em"></mml:mspace><mml:mi>η</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>L</mml:mi><mml:mi>o</mml:mi><mml:mi>w</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi><mml:mspace width="0.25em"></mml:mspace><mml:mi>d</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi><mml:mi>k</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mn>0.5</mml:mn><mml:mo>,</mml:mo><mml:mspace width="0.5em"></mml:mspace><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>η</mml:mi></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mspace width="0.5em"></mml:mspace><mml:mi>g</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mspace width="0.5em"></mml:mspace><mml:mi>m</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>n</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mspace width="2em"></mml:mspace><mml:mtext>at</mml:mtext><mml:mspace width="2em"></mml:mspace><mml:mi>η</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>U</mml:mi><mml:mi>p</mml:mi><mml:mi>p</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi><mml:mspace width="0.25em"></mml:mspace><mml:mi>d</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi><mml:mi>k</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></alternatives><label>(12)</label></disp-formula></p><p>In tribological applications, we can further define the dimensionless frictional moment’s i.e. non-dimensional torques exerted on the upper and lower disks. The torque at the upper disk is
<disp-formula id="pone.0135004.e016"><alternatives><graphic xlink:href="pone.0135004.e016.jpg" id="pone.0135004.e016g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M16"><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mi>p</mml:mi><mml:mi>p</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi><mml:mspace width="0.25em"></mml:mspace><mml:mi>d</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>2</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:mi>π</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>μ</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mstyle displaystyle="true"><mml:mrow><mml:msubsup><mml:mi>∫</mml:mi><mml:mn>0</mml:mn><mml:mi>a</mml:mi></mml:msubsup><mml:mrow><mml:msub><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mo>∂</mml:mo><mml:mi>v</mml:mi></mml:mrow><mml:mrow><mml:mo>∂</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>z</mml:mi><mml:mo>=</mml:mo><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:mi>d</mml:mi><mml:mi>r</mml:mi></mml:mrow></mml:mrow></mml:mstyle><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(13)</label></disp-formula>where <italic>μ</italic> m denotes the Newtonian dynamic viscosity. Using (7) and (13), we obtain:
<disp-formula id="pone.0135004.e017"><alternatives><graphic xlink:href="pone.0135004.e017.jpg" id="pone.0135004.e017g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M17"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>T</mml:mi><mml:mo>¯</mml:mo></mml:mover><mml:mrow><mml:mi>U</mml:mi><mml:mi>p</mml:mi><mml:mi>p</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi><mml:mspace width="0.25em"></mml:mspace><mml:mi>d</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mn>2</mml:mn><mml:mi>D</mml:mi><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>α</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>ξ</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mn>3</mml:mn><mml:mo>/</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:mrow></mml:msup></mml:mrow><mml:mrow><mml:mi>π</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>μ</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi mathvariant="normal">Ω</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:msup><mml:mi>a</mml:mi><mml:mn>4</mml:mn></mml:msup></mml:mrow></mml:mfrac><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mi>p</mml:mi><mml:mi>p</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi><mml:mspace width="0.25em"></mml:mspace><mml:mi>d</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>g</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>η</mml:mi></mml:mrow></mml:mfrac><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(14)</label></disp-formula>where <inline-formula id="pone.0135004.e018"><alternatives><graphic xlink:href="pone.0135004.e018.jpg" id="pone.0135004.e018g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M18"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>T</mml:mi><mml:mo>¯</mml:mo></mml:mover><mml:mrow><mml:mi>U</mml:mi><mml:mi>p</mml:mi><mml:mi>p</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi><mml:mspace width="0.25em"></mml:mspace><mml:mi>d</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi><mml:mi>k</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></alternatives></inline-formula> is the dimensionless torque exerted by the fluid on the upper disk and <italic>dg</italic>(1)/<italic>dη</italic> is the azimuthal (tangential) velocity gradient at the upper disk (<italic>η</italic> = 1). Similarly for the lower disk, the dimensionless torque is simply given by the same calculation evaluated at <italic>η</italic> = 0, as
<disp-formula id="pone.0135004.e019"><alternatives><graphic xlink:href="pone.0135004.e019.jpg" id="pone.0135004.e019g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M19"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>T</mml:mi><mml:mo>¯</mml:mo></mml:mover><mml:mrow><mml:mi>L</mml:mi><mml:mi>o</mml:mi><mml:mi>w</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi><mml:mspace width="0.25em"></mml:mspace><mml:mi>d</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>g</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mn>0</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>η</mml:mi></mml:mrow></mml:mfrac><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(15)</label></disp-formula>where <italic>dg</italic>(0)/<italic>dη</italic> is the tangential (azimuthal) velocity gradient at the lower disk (<italic>η</italic> = 0).</p></sec><sec id="sec003"><title>The Differential Transform Method</title><p>DTM is employed to obtain semi- analytical/numerical solutions to the well-posed two-point boundary value problem defined by Eqs (<xref rid="pone.0135004.e009" ref-type="disp-formula">8</xref>)–(<xref rid="pone.0135004.e012" ref-type="disp-formula">11</xref>) and conditions (12). DTM is an extremely strong technique in finding solutions to magneto-hydrodynamic and complex material flow problems. It has also been used very effectively in conjunction with Padé approximants. To provide a summary of the method, the transformation of the <italic>k</italic><sup>th</sup> derivative of a function in one variable is considered which is defined as:
<disp-formula id="pone.0135004.e020"><alternatives><graphic xlink:href="pone.0135004.e020.jpg" id="pone.0135004.e020g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M20"><mml:mrow><mml:mi>F</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>k</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mi>k</mml:mi><mml:mo>!</mml:mo></mml:mrow></mml:mfrac><mml:msub><mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:msup><mml:mi>d</mml:mi><mml:mi>k</mml:mi></mml:msup><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>η</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:msup><mml:mi>η</mml:mi><mml:mi>k</mml:mi></mml:msup></mml:mrow></mml:mfrac></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>η</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mi>η</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow></mml:msub><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(16)</label></disp-formula>where <italic>f</italic> (<italic>η</italic>) is the original function and <italic>F</italic>(<italic>k</italic>) is transformed function. The differential inverse transformation of <italic>F</italic>(<italic>k</italic>) is:
<disp-formula id="pone.0135004.e021"><alternatives><graphic xlink:href="pone.0135004.e021.jpg" id="pone.0135004.e021g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M21"><mml:mrow><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>η</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>∞</mml:mi></mml:munderover><mml:mrow><mml:mi>F</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>k</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>η</mml:mi><mml:mo>−</mml:mo><mml:msub><mml:mi>η</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mi>k</mml:mi></mml:msup></mml:mrow></mml:mstyle><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(17)</label></disp-formula></p><p>The concept of the differential transform is derived from a Taylor series expansion and in actual applications the function <italic>f</italic> (<italic>η</italic>) is expressed by a finite series as follows:
<disp-formula id="pone.0135004.e022"><alternatives><graphic xlink:href="pone.0135004.e022.jpg" id="pone.0135004.e022g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M22"><mml:mrow><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>η</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>≅</mml:mo><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>m</mml:mi></mml:munderover><mml:mrow><mml:mi>F</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>k</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>η</mml:mi><mml:mo>−</mml:mo><mml:msub><mml:mi>η</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mi>k</mml:mi></mml:msup></mml:mrow></mml:mstyle><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(18)</label></disp-formula></p><p>The value of <italic>m</italic> is decided by convergence of the series coefficients.</p><sec id="sec004"><title>3.1. Padé Approximant</title><p>Páde approximants are applied to the problem to increase the convergence of a given series. Suppose that a power series <inline-formula id="pone.0135004.e023"><alternatives><graphic xlink:href="pone.0135004.e023.jpg" id="pone.0135004.e023g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M23"><mml:mrow><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>∞</mml:mi></mml:munderover><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:msup></mml:mrow></mml:mstyle></mml:mrow></mml:math></alternatives></inline-formula> is given, which represents a function <italic>f</italic> (<italic>x</italic>), such that:
<disp-formula id="pone.0135004.e024"><alternatives><graphic xlink:href="pone.0135004.e024.jpg" id="pone.0135004.e024g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M24"><mml:mrow><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>∞</mml:mi></mml:munderover><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:msup></mml:mrow></mml:mstyle><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(19)</label></disp-formula></p><p>The Páde approximant is a rational fraction and the notation for such a Padé approximant is:
<disp-formula id="pone.0135004.e025"><alternatives><graphic xlink:href="pone.0135004.e025.jpg" id="pone.0135004.e025g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M25"><mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mrow><mml:mi>L</mml:mi><mml:mo>/</mml:mo><mml:mi>M</mml:mi></mml:mrow></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi>L</mml:mi></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:msub><mml:mi>Q</mml:mi><mml:mi>M</mml:mi></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mfrac><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(20)</label></disp-formula>where <italic>P</italic><sub><italic>L</italic></sub>(<italic>x</italic>) is a polynomial of degree at most <italic>L</italic> and <italic>Q</italic><sub><italic>M</italic></sub> (<italic>x</italic>) is a polynomial of degree at most <italic>M</italic>. Therefore:
<disp-formula id="pone.0135004.e026"><alternatives><graphic xlink:href="pone.0135004.e026.jpg" id="pone.0135004.e026g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M26"><mml:mrow><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mi>x</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mn>3</mml:mn></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mn>3</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mn>4</mml:mn></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mn>4</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mo>⋯</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(21)</label></disp-formula><disp-formula id="pone.0135004.e027"><alternatives><graphic xlink:href="pone.0135004.e027.jpg" id="pone.0135004.e027g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M27"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi>L</mml:mi></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mi>x</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mn>3</mml:mn></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mn>3</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mo>…</mml:mo><mml:mo>+</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mi>L</mml:mi></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mi>L</mml:mi></mml:msup><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(22)</label></disp-formula><disp-formula id="pone.0135004.e028"><alternatives><graphic xlink:href="pone.0135004.e028.jpg" id="pone.0135004.e028g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M28"><mml:mrow><mml:msub><mml:mi>Q</mml:mi><mml:mi>M</mml:mi></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mi>x</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mn>3</mml:mn></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mn>3</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mo>…</mml:mo><mml:mo>+</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mi>M</mml:mi></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mi>M</mml:mi></mml:msup><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(23)</label></disp-formula>where in Eq (<xref rid="pone.0135004.e025" ref-type="disp-formula">20</xref>) there are <italic>L</italic> + 1 numerator coefficients and <italic>M</italic> + 1 denominator coefficients. Since the numerator and denominator can be multiplied by a constant and [<italic>L</italic>/<italic>M</italic>] left unchanged, the following normalization condition is imposed
<disp-formula id="pone.0135004.e029"><alternatives><graphic xlink:href="pone.0135004.e029.jpg" id="pone.0135004.e029g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M29"><mml:mrow><mml:msub><mml:mi>Q</mml:mi><mml:mi>M</mml:mi></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mn>0</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(24)</label></disp-formula></p><p>So there are <italic>L</italic> + 1 independent numerator coefficients and <italic>M</italic> independent denominator coefficients, which make <italic>L</italic> + <italic>M</italic> + 1 unknown coefficients in all. This number suggests that normally the [<italic>L</italic>/<italic>M</italic>] ought to fit the power series Eq (<xref rid="pone.0135004.e024" ref-type="disp-formula">19</xref>) through the orders 1, <italic>x</italic>, <italic>x</italic><sup>2</sup>,…, <italic>x</italic><sup><italic>L</italic>+<italic>M</italic></sup>. Based on conditions given in [<xref rid="pone.0135004.ref038" ref-type="bibr">38</xref>, <xref rid="pone.0135004.ref039" ref-type="bibr">39</xref>], [<italic>L</italic>/<italic>M</italic>] approximation is uniquely determined. In the notation of formal power series:
<disp-formula id="pone.0135004.e030"><alternatives><graphic xlink:href="pone.0135004.e030.jpg" id="pone.0135004.e030g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M30"><mml:mrow><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>∞</mml:mi></mml:munderover><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:msup><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mi>x</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mn>3</mml:mn></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mn>3</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mo>…</mml:mo><mml:mo>+</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mi>L</mml:mi></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mi>L</mml:mi></mml:msup></mml:mrow><mml:mrow><mml:msub><mml:mi>q</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mi>x</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mn>3</mml:mn></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mn>3</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mo>…</mml:mo><mml:mo>+</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mi>M</mml:mi></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mi>M</mml:mi></mml:msup></mml:mrow></mml:mfrac></mml:mrow></mml:mstyle><mml:mo>+</mml:mo><mml:mi>O</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msup><mml:mi>x</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>+</mml:mo><mml:mi>M</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(25)</label></disp-formula></p><p>By cross-multiplying Eq (<xref rid="pone.0135004.e030" ref-type="disp-formula">25</xref>), one obtains:
<disp-formula id="pone.0135004.e031"><alternatives><graphic xlink:href="pone.0135004.e031.jpg" id="pone.0135004.e031g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M31"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mi>x</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mn>3</mml:mn></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mn>3</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mo>…</mml:mo><mml:mo>+</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mi>L</mml:mi></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mi>L</mml:mi></mml:msup></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>×</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mi>x</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mn>3</mml:mn></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mn>3</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mn>4</mml:mn></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mn>4</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mo>…</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mspace width="9.3em"></mml:mspace><mml:mo>=</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mi>x</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mn>3</mml:mn></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mn>3</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mo>…</mml:mo><mml:mo>+</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mi>M</mml:mi></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mi>M</mml:mi></mml:msup><mml:mo>+</mml:mo><mml:mi>O</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msup><mml:mi>x</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>+</mml:mo><mml:mi>M</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></alternatives><label>(26)</label></disp-formula></p><p>From Eq (<xref rid="pone.0135004.e031" ref-type="disp-formula">26</xref>) the following set of linear equations are obtained
<disp-formula id="pone.0135004.e032"><alternatives><graphic xlink:href="pone.0135004.e032.jpg" id="pone.0135004.e032g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M32"><mml:mrow><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mi>a</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>,</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>a</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>q</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mo>,</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>a</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>q</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>q</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mo>,</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mo>⋮</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>a</mml:mi><mml:mi>L</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>−</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>q</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:mo>⋯</mml:mo><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>q</mml:mi><mml:mi>L</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mi>L</mml:mi></mml:msub><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mrow></mml:math></alternatives><label>(27)</label></disp-formula>and
<disp-formula id="pone.0135004.e033"><alternatives><graphic xlink:href="pone.0135004.e033.jpg" id="pone.0135004.e033g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M33"><mml:mrow><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mi>L</mml:mi></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>q</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:mo>⋯</mml:mo><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>−</mml:mo><mml:mi>M</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>q</mml:mi><mml:mi>M</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>+</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>q</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:mo>⋯</mml:mo><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>−</mml:mo><mml:mi>M</mml:mi><mml:mo>+</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>q</mml:mi><mml:mi>M</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mo>⋮</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>+</mml:mo><mml:mi>M</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>+</mml:mo><mml:mi>M</mml:mi><mml:mo>−</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>q</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:mo>⋯</mml:mo><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mi>L</mml:mi></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>q</mml:mi><mml:mi>M</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mrow></mml:math></alternatives><label>(28)</label></disp-formula>where <italic>a</italic><sub><italic>n</italic></sub> = 0 for <italic>n</italic> < 0 and <italic>q</italic><sub><italic>j</italic></sub> = 0 for <italic>j</italic> > <italic>M</italic>. Eqs (<xref rid="pone.0135004.e032" ref-type="disp-formula">27</xref>) and (<xref rid="pone.0135004.e033" ref-type="disp-formula">28</xref>) can be solved directly if they are non-singular
<disp-formula id="pone.0135004.e034"><alternatives><graphic xlink:href="pone.0135004.e034.jpg" id="pone.0135004.e034g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M34"><mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mi>L</mml:mi><mml:mo>/</mml:mo><mml:mi>M</mml:mi></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mrow><mml:mo>|</mml:mo><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>−</mml:mo><mml:mi>m</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>−</mml:mo><mml:mi>M</mml:mi><mml:mo>+</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mo>…</mml:mo></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mo>⋮</mml:mo></mml:mtd><mml:mtd><mml:mo>⋮</mml:mo></mml:mtd><mml:mtd><mml:mo>⋱</mml:mo></mml:mtd><mml:mtd><mml:mo>⋮</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mi>L</mml:mi></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mo>…</mml:mo></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>+</mml:mo><mml:mi>M</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>j</mml:mi><mml:mo>=</mml:mo><mml:mi>M</mml:mi></mml:mrow><mml:mi>L</mml:mi></mml:munderover><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>j</mml:mi><mml:mo>−</mml:mo><mml:mi>M</mml:mi></mml:mrow></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mi>j</mml:mi></mml:msup></mml:mrow></mml:mstyle></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>j</mml:mi><mml:mo>=</mml:mo><mml:mi>M</mml:mi><mml:mo>−</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>L</mml:mi></mml:munderover><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>j</mml:mi><mml:mo>−</mml:mo><mml:mi>M</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mi>j</mml:mi></mml:msup></mml:mrow></mml:mstyle></mml:mrow></mml:mtd><mml:mtd><mml:mo>…</mml:mo></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>j</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>L</mml:mi></mml:munderover><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mi>j</mml:mi></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mi>j</mml:mi></mml:msup></mml:mrow></mml:mstyle></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>|</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mo>|</mml:mo><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>−</mml:mo><mml:mi>m</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>−</mml:mo><mml:mi>m</mml:mi><mml:mo>+</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mo>…</mml:mo></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mo>⋮</mml:mo></mml:mtd><mml:mtd><mml:mo>⋮</mml:mo></mml:mtd><mml:mtd><mml:mo>⋱</mml:mo></mml:mtd><mml:mtd><mml:mo>⋮</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mi>L</mml:mi></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mo>…</mml:mo></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>+</mml:mo><mml:mi>M</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msup><mml:mi>x</mml:mi><mml:mi>M</mml:mi></mml:msup></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msup><mml:mi>x</mml:mi><mml:mrow><mml:mi>M</mml:mi><mml:mo>−</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:mtd><mml:mtd><mml:mo>…</mml:mo></mml:mtd><mml:mtd><mml:mn>1</mml:mn></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>|</mml:mo></mml:mrow></mml:mrow></mml:mfrac><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(29)</label></disp-formula></p><p>If the lower index on a sum exceeds the upper, the sum is replaced by zero. Alternate forms are:
<disp-formula id="pone.0135004.e035"><alternatives><graphic xlink:href="pone.0135004.e035.jpg" id="pone.0135004.e035g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M35"><mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mrow><mml:mi>L</mml:mi><mml:mo>/</mml:mo><mml:mi>M</mml:mi></mml:mrow></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>j</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mi>L</mml:mi><mml:mo>−</mml:mo><mml:mi>M</mml:mi></mml:mrow></mml:munderover><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mi>j</mml:mi></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mi>j</mml:mi></mml:msup><mml:mo>+</mml:mo><mml:msup><mml:mi>x</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>−</mml:mo><mml:mi>M</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:msubsup><mml:mi>w</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>/</mml:mo><mml:mi>M</mml:mi></mml:mrow><mml:mi>T</mml:mi></mml:msubsup><mml:msubsup><mml:mi>W</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>/</mml:mo><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mo>−</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msubsup></mml:mrow></mml:mstyle><mml:msub><mml:mi>w</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>/</mml:mo><mml:mi>M</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>j</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mi>L</mml:mi><mml:mo>+</mml:mo><mml:mi>n</mml:mi></mml:mrow></mml:munderover><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mi>j</mml:mi></mml:msub><mml:msup><mml:mi>x</mml:mi><mml:mi>j</mml:mi></mml:msup><mml:mo>+</mml:mo><mml:msup><mml:mi>x</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>+</mml:mo><mml:mi>n</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:msubsup><mml:mi>w</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>L</mml:mi><mml:mo>+</mml:mo><mml:mi>M</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>/</mml:mo><mml:mi>M</mml:mi></mml:mrow><mml:mi>T</mml:mi></mml:msubsup><mml:msubsup><mml:mi>W</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>/</mml:mo><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mo>−</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msubsup></mml:mrow></mml:mstyle><mml:msub><mml:mi>w</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>L</mml:mi><mml:mo>+</mml:mo><mml:mi>n</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>/</mml:mo><mml:mi>M</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(30)</label></disp-formula></p><p>For
<disp-formula id="pone.0135004.e036"><alternatives><graphic xlink:href="pone.0135004.e036.jpg" id="pone.0135004.e036g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M36"><mml:mrow><mml:msub><mml:mi>W</mml:mi><mml:mrow><mml:mrow><mml:mi>L</mml:mi><mml:mo>/</mml:mo><mml:mi>M</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>−</mml:mo><mml:mi>M</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>−</mml:mo><mml:mi>x</mml:mi><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>−</mml:mo><mml:mi>M</mml:mi><mml:mo>+</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mo>…</mml:mo></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mi>L</mml:mi></mml:msub><mml:mo>−</mml:mo><mml:mi>x</mml:mi><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mo>⋮</mml:mo></mml:mtd><mml:mtd><mml:mo>⋱</mml:mo></mml:mtd><mml:mtd><mml:mo>⋮</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mi>L</mml:mi></mml:msub><mml:mo>−</mml:mo><mml:mi>x</mml:mi><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mo>…</mml:mo></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>+</mml:mo><mml:mi>M</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>−</mml:mo><mml:mi>x</mml:mi><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>+</mml:mo><mml:mi>M</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(31)</label></disp-formula><disp-formula id="pone.0135004.e037"><alternatives><graphic xlink:href="pone.0135004.e037.jpg" id="pone.0135004.e037g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M37"><mml:mrow><mml:msub><mml:mi>w</mml:mi><mml:mrow><mml:mrow><mml:mi>L</mml:mi><mml:mo>/</mml:mo><mml:mi>M</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>−</mml:mo><mml:mi>M</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>−</mml:mo><mml:mi>M</mml:mi><mml:mo>+</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mo>⋮</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mi>L</mml:mi></mml:msub></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(32)</label></disp-formula></p><p>The construction of [<italic>L</italic>/<italic>M</italic>] approximants involves only algebraic operations [<xref rid="pone.0135004.ref038" ref-type="bibr">38</xref>, <xref rid="pone.0135004.ref039" ref-type="bibr">39</xref>]. Each choice of <italic>L</italic>, degree of the numerator and <italic>M</italic>, degree of the denominator, leads to an approximant. How to direct the choice in order to obtain the best approximant is the major difficulty in applying the technique, which necessitates the need for a criterion for the choice depending on the s<italic>hape</italic> of the solution. A criterion which has worked well here is the choice of [<italic>L</italic>/<italic>M</italic>] approximants such that <italic>L</italic> = <italic>M</italic>.</p></sec><sec id="sec005"><title>3.2. Analytical approximation by means of DTM-Padé</title><p>Taking differential transform of Eqs (<xref rid="pone.0135004.e009" ref-type="disp-formula">8</xref>)–(<xref rid="pone.0135004.e012" ref-type="disp-formula">11</xref>), one can obtain (for more details, see [<xref rid="pone.0135004.ref040" ref-type="bibr">40</xref>–<xref rid="pone.0135004.ref042" ref-type="bibr">42</xref>])
<disp-formula id="pone.0135004.e038"><alternatives><graphic xlink:href="pone.0135004.e038.jpg" id="pone.0135004.e038g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M38"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>2</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>3</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>4</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>4</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>−</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mspace width="3em"></mml:mspace><mml:msub><mml:mi>N</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mspace width="7em"></mml:mspace><mml:mn>3</mml:mn><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>2</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>2</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>r</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>k</mml:mi></mml:munderover><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>2</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>3</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>δ</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>r</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>3</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>−</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mn>2</mml:mn><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>2</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>3</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>r</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>3</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mstyle></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mspace width="1.8em"></mml:mspace><mml:mo>−</mml:mo><mml:mn>2</mml:mn><mml:msubsup><mml:mi>N</mml:mi><mml:mn>2</mml:mn><mml:mn>2</mml:mn></mml:msubsup><mml:mspace width="0.15em"></mml:mspace><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>r</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>k</mml:mi></mml:munderover><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>2</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>m</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>r</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mi>m</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>2</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mstyle></mml:mtd></mml:mtr></mml:mtable><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></alternatives><label>(33)</label></disp-formula><disp-formula id="pone.0135004.e039"><alternatives><graphic xlink:href="pone.0135004.e039.jpg" id="pone.0135004.e039g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M39"><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>2</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>g</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>2</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>−</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mn>2</mml:mn><mml:mspace width="0.1em"></mml:mspace><mml:mi>g</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>k</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>r</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>k</mml:mi></mml:munderover><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>δ</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>r</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mi>g</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mn>2</mml:mn><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>g</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>r</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>−</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mn>2</mml:mn><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>r</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mi>g</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mstyle></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mo>+</mml:mo><mml:mn>2</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>N</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>N</mml:mi><mml:mn>3</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>r</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>k</mml:mi></mml:munderover><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>m</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>r</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mi>n</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>−</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>n</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>r</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mi>m</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mstyle></mml:mtd></mml:mtr></mml:mtable><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(34)</label></disp-formula><disp-formula id="pone.0135004.e040"><alternatives><graphic xlink:href="pone.0135004.e040.jpg" id="pone.0135004.e040g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M40"><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>2</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>m</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>2</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>−</mml:mo><mml:mi>R</mml:mi><mml:msub><mml:mi>e</mml:mi><mml:mi>m</mml:mi></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>m</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>k</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>r</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>k</mml:mi></mml:munderover><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>δ</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>r</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mi>m</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>−</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mn>2</mml:mn><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>r</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mi>m</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mn>2</mml:mn><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>m</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>r</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mstyle></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(35)</label></disp-formula><disp-formula id="pone.0135004.e041"><alternatives><graphic xlink:href="pone.0135004.e041.jpg" id="pone.0135004.e041g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M41"><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>2</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>n</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>2</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>−</mml:mo><mml:mi>R</mml:mi><mml:msub><mml:mi>e</mml:mi><mml:mi>m</mml:mi></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mn>2</mml:mn><mml:mspace width="0.1em"></mml:mspace><mml:mi>n</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>k</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>r</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>k</mml:mi></mml:munderover><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>δ</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>r</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mi>n</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>−</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mn>2</mml:mn><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>r</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mi>n</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mstyle></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mo>+</mml:mo><mml:mn>2</mml:mn><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mo>/</mml:mo><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow></mml:mrow></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>r</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>k</mml:mi></mml:munderover><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>m</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>r</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mi>g</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>−</mml:mo><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mstyle></mml:mtd></mml:mtr></mml:mtable><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(36)</label></disp-formula>where <italic>f</italic> (<italic>k</italic>), <italic>g</italic> (<italic>k</italic>), <italic>m</italic> (<italic>k</italic>) and <italic>n</italic> (<italic>k</italic>) are the differential transforms of <italic>f</italic> (<italic>η</italic>), <italic>g</italic> (<italic>η</italic>), <italic>m</italic> (<italic>η</italic>) and <italic>n</italic> (<italic>η</italic>) are displayed by:
<disp-formula id="pone.0135004.e042"><alternatives><graphic xlink:href="pone.0135004.e042.jpg" id="pone.0135004.e042g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M42"><mml:mrow><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>η</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>∞</mml:mi></mml:munderover><mml:mrow><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>k</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:msup><mml:mi>η</mml:mi><mml:mi>k</mml:mi></mml:msup></mml:mrow></mml:mstyle><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(37)</label></disp-formula><disp-formula id="pone.0135004.e043"><alternatives><graphic xlink:href="pone.0135004.e043.jpg" id="pone.0135004.e043g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M43"><mml:mrow><mml:mi>g</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>η</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>∞</mml:mi></mml:munderover><mml:mrow><mml:mi>g</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>k</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:msup><mml:mi>η</mml:mi><mml:mi>k</mml:mi></mml:msup></mml:mrow></mml:mstyle><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(38)</label></disp-formula><disp-formula id="pone.0135004.e044"><alternatives><graphic xlink:href="pone.0135004.e044.jpg" id="pone.0135004.e044g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M44"><mml:mrow><mml:mi>m</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>η</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>∞</mml:mi></mml:munderover><mml:mrow><mml:mi>m</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>k</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:msup><mml:mi>η</mml:mi><mml:mi>k</mml:mi></mml:msup></mml:mrow></mml:mstyle><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(39)</label></disp-formula><disp-formula id="pone.0135004.e045"><alternatives><graphic xlink:href="pone.0135004.e045.jpg" id="pone.0135004.e045g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M45"><mml:mrow><mml:mi>n</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>η</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mi>∞</mml:mi></mml:munderover><mml:mrow><mml:mi>n</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>k</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:msup><mml:mi>η</mml:mi><mml:mi>k</mml:mi></mml:msup></mml:mrow></mml:mstyle><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(40)</label></disp-formula><disp-formula id="pone.0135004.e046"><alternatives><graphic xlink:href="pone.0135004.e046.jpg" id="pone.0135004.e046g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M46"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mn>0</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mspace width="0.25em"></mml:mspace><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mspace width="0.25em"></mml:mspace><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mn>2</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>α</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.25em"></mml:mspace><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mn>3</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>β</mml:mi><mml:mo>,</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mi>g</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mn>0</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mspace width="0.25em"></mml:mspace><mml:mi>g</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>γ</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.25em"></mml:mspace><mml:mi>m</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mn>0</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mspace width="0.25em"></mml:mspace><mml:mi>m</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>κ</mml:mi><mml:mo>,</mml:mo><mml:mspace width="0.25em"></mml:mspace><mml:mi>n</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mn>0</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mspace width="0.25em"></mml:mspace><mml:mi>n</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>ω</mml:mi><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></alternatives><label>(41)</label></disp-formula>where <italic>α</italic>, <italic>β</italic>, <italic>γ</italic>, <italic>κ</italic> and <italic>ω</italic> are constants. By substituting Eq (<xref rid="pone.0135004.e046" ref-type="disp-formula">41</xref>) into Eqs (<xref rid="pone.0135004.e038" ref-type="disp-formula">33</xref>)–(<xref rid="pone.0135004.e041" ref-type="disp-formula">36</xref>), we obtain the values of <italic>f</italic> (<italic>η</italic>), <italic>g</italic> (<italic>η</italic>), <italic>m</italic> (<italic>η</italic>) and <italic>n</italic>(<italic>η</italic>).
<disp-formula id="pone.0135004.e047"><alternatives><graphic xlink:href="pone.0135004.e047.jpg" id="pone.0135004.e047g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M47"><mml:mrow><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>η</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>α</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:msup><mml:mi>η</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mi>β</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:msup><mml:mi>η</mml:mi><mml:mn>3</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mn>4</mml:mn></mml:mfrac><mml:msub><mml:mi>N</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:mi>α</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:msup><mml:mi>η</mml:mi><mml:mn>4</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mn>5</mml:mn></mml:mfrac><mml:msub><mml:mi>N</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:mi>β</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:msup><mml:mi>η</mml:mi><mml:mn>5</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mo>⋯</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(42)</label></disp-formula><disp-formula id="pone.0135004.e048"><alternatives><graphic xlink:href="pone.0135004.e048.jpg" id="pone.0135004.e048g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M48"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mi>g</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>η</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:mi>γ</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>η</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:msup><mml:mi>η</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mn>6</mml:mn></mml:mfrac><mml:msub><mml:mi>N</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mn>4</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:mi>α</mml:mi><mml:mo>+</mml:mo><mml:mn>3</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:mi>γ</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:msup><mml:mi>η</mml:mi><mml:mn>3</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mn>12</mml:mn></mml:mrow></mml:mfrac><mml:msub><mml:mi>N</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mn>4</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>N</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:mn>6</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:mi>β</mml:mi><mml:mo>+</mml:mo><mml:mn>2</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:mi>α</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>γ</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:msup><mml:mi>η</mml:mi><mml:mn>4</mml:mn></mml:msup></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mspace width="2.5em"></mml:mspace><mml:mo>+</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mn>20</mml:mn></mml:mrow></mml:mfrac><mml:msub><mml:mi>N</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mspace width="3em"></mml:mspace><mml:mn>2</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>N</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:mi>α</mml:mi><mml:mo>+</mml:mo><mml:mn>4</mml:mn><mml:mi>β</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>γ</mml:mi><mml:mo>+</mml:mo><mml:mfrac><mml:mn>5</mml:mn><mml:mn>6</mml:mn></mml:mfrac><mml:msub><mml:mi>N</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mn>4</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:mi>α</mml:mi><mml:mo>+</mml:mo><mml:mn>3</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:mi>γ</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mn>2</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>N</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>N</mml:mi><mml:mn>3</mml:mn></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mo>−</mml:mo><mml:mfrac><mml:mn>2</mml:mn><mml:mn>3</mml:mn></mml:mfrac><mml:mi>R</mml:mi><mml:msub><mml:mi>e</mml:mi><mml:mi>m</mml:mi></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:mi>κ</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>ω</mml:mi><mml:mo>+</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mn>3</mml:mn></mml:mfrac><mml:mi>R</mml:mi><mml:msub><mml:mi>e</mml:mi><mml:mi>m</mml:mi></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:mi>κ</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mn>2</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>N</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:mi>γ</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>κ</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow></mml:mfrac><mml:mo>+</mml:mo><mml:mn>3</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:mi>ω</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable><mml:mo>)</mml:mo></mml:mrow><mml:msup><mml:mi>η</mml:mi><mml:mn>5</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mo>⋯</mml:mo><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></alternatives><label>(43)</label></disp-formula><disp-formula id="pone.0135004.e049"><alternatives><graphic xlink:href="pone.0135004.e049.jpg" id="pone.0135004.e049g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M49"><mml:mrow><mml:mi>m</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>η</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>κ</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>η</mml:mi><mml:mo>+</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mn>3</mml:mn></mml:mfrac><mml:mi>R</mml:mi><mml:msub><mml:mi>e</mml:mi><mml:mi>m</mml:mi></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:mi>κ</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:msup><mml:mi>η</mml:mi><mml:mn>3</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mn>6</mml:mn></mml:mfrac><mml:mi>R</mml:mi><mml:msub><mml:mi>e</mml:mi><mml:mi>m</mml:mi></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:mi>α</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>κ</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:msup><mml:mi>η</mml:mi><mml:mn>4</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mn>20</mml:mn></mml:mrow></mml:mfrac><mml:mi>R</mml:mi><mml:msub><mml:mi>e</mml:mi><mml:mi>m</mml:mi></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mn>4</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:mi>R</mml:mi><mml:msub><mml:mi>e</mml:mi><mml:mi>m</mml:mi></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:mi>κ</mml:mi></mml:mrow><mml:mn>3</mml:mn></mml:mfrac><mml:mo>+</mml:mo><mml:mn>4</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:mi>β</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>κ</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:msup><mml:mi>η</mml:mi><mml:mn>5</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mo>⋯</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:math></alternatives><label>(44)</label></disp-formula><disp-formula id="pone.0135004.e050"><alternatives><graphic xlink:href="pone.0135004.e050.jpg" id="pone.0135004.e050g" position="anchor" mimetype="image" orientation="portrait"></graphic><mml:math id="M50"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mi>n</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>η</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>ω</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>η</mml:mi><mml:mo>+</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mn>6</mml:mn></mml:mfrac><mml:mi>R</mml:mi><mml:msub><mml:mi>e</mml:mi><mml:mi>m</mml:mi></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mn>2</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>N</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:mi>γ</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>κ</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow></mml:mfrac><mml:mo>+</mml:mo><mml:mn>3</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:mi>ω</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:msup><mml:mi>η</mml:mi><mml:mn>3</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mn>12</mml:mn></mml:mrow></mml:mfrac><mml:mi>R</mml:mi><mml:msub><mml:mi>e</mml:mi><mml:mi>m</mml:mi></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mn>4</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>N</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>N</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:mi>κ</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow></mml:mfrac><mml:mo>−</mml:mo><mml:mn>2</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:mi>α</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>ω</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:msup><mml:mi>η</mml:mi><mml:mn>4</mml:mn></mml:msup><mml:mo>+</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mspace width="3em"></mml:mspace><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mn>20</mml:mn></mml:mrow></mml:mfrac><mml:mi>R</mml:mi><mml:msub><mml:mi>e</mml:mi><mml:mi>m</mml:mi></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mn>2</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>N</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:mi>R</mml:mi><mml:msub><mml:mi>e</mml:mi><mml:mi>m</mml:mi></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:mi>γ</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>κ</mml:mi></mml:mrow><mml:mrow><mml:mn>3</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>N</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow></mml:mfrac><mml:mo>+</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>N</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mn>4</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:mi>α</mml:mi><mml:mo>+</mml:mo><mml:mn>3</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:mi>γ</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>κ</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow></mml:mfrac><mml:mo>−</mml:mo><mml:mn>2</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:mi>β</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>ω</mml:mi><mml:mo>+</mml:mo><mml:mfrac><mml:mn>5</mml:mn><mml:mn>6</mml:mn></mml:mfrac><mml:mi>R</mml:mi><mml:msub><mml:mi>e</mml:mi><mml:mi>m</mml:mi></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mn>2</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:msub><mml:mi>N</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mspace width="0.15em"></mml:mspace><mml:mi>γ</mml:mi><mml:mspace width="0.15em"></mml:mspace><mml:mi>κ</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow></mml:mfrac><mml:mo>+</mml:mo><mml:mn>3</mml:mn><mml:mspace width="0.15em"></mml:mspace><mml:mi>ω</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:msup><mml:mi>η</mml:mi><mml:mn>5</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mo>⋯</mml:mo><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></alternatives><label>(45)</label></disp-formula>where the number of required terms is determined by the convergence of the numerical values to one’s desired accuracy. We obtain the approximants using <italic>MATHEMATICA</italic> software. The Padé approximant is employed to extend the convergence radius of the truncated series solution. As it is illustrated in Figs <xref rid="pone.0135004.g002" ref-type="fig">2</xref>–<xref rid="pone.0135004.g003" ref-type="fig">3</xref>, without using the Padé approximant, the different orders of DTM solution cannot satisfy boundary conditions at infinity. Therefore, it is necessary to apply DTM-Padé to provide an effective way to handle infinite boundary value problems. The Padé approximant is applied to Eqs (<xref rid="pone.0135004.e047" ref-type="disp-formula">42</xref>)–(<xref rid="pone.0135004.e050" ref-type="disp-formula">45</xref>) and one can obtain <italic>α</italic>, <italic>β</italic>, <italic>γ</italic>, <italic>κ</italic> and <italic>ω</italic>.</p><fig id="pone.0135004.g002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0135004.g002</object-id><label>Fig 2</label><caption><title>The obtained results of tangential velocity distribution (<italic>g</italic>) for different orders of DTM and DTM-Padé solutions in comparison with the numerical solution when <italic>N</italic><sub>1</sub> = <italic>N</italic><sub>2</sub> = <italic>N</italic><sub>3</sub> = 1 and <italic>Bt</italic> = 6.</title></caption><graphic xlink:href="pone.0135004.g002"></graphic></fig><fig id="pone.0135004.g003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0135004.g003</object-id><label>Fig 3</label><caption><title>The obtained results of tangential induced magnetic field distribution (<italic>n</italic>) for different orders of DTM and DTM-Padé solutions in comparison with the numerical solution when <italic>N</italic><sub>1</sub> = <italic>N</italic><sub>2</sub> = <italic>N</italic><sub>3</sub> = 1 and <italic>Bt</italic> = 6.</title></caption><graphic xlink:href="pone.0135004.g003"></graphic></fig></sec></sec><sec id="sec006"><title>Results and Discussion</title><p>The nonlinear ordinary differential equations subject to the boundary conditions are solved via DTM-Padé method for some values of the four key parameters i.e. squeeze Reynolds number (<italic>N</italic><sub>1</sub>), magnetic Reynolds number (<italic>Re</italic><sub><italic>m</italic></sub> = <italic>N</italic><sub>1</sub><italic>Bt</italic>), dimensionless axial magnetic force parameter (<italic>N</italic><sub>2</sub>), dimensionless tangential magnetic force parameter (<italic>N</italic><sub>3</sub>) on the velocity and induced magnetic field components in the gap between the disks and also on the torques at the lower and upper disk. Computations are performed for the evolution of velocity components (<italic>f</italic>, <italic>g</italic>) and induced magnetic field components (<italic>m</italic>, <italic>n</italic>) with dimensionless axial coordinate (<italic>η</italic>). Representative values are used to simulate physically realistic flows. <xref rid="pone.0135004.t001" ref-type="table">Table 1</xref> and <xref rid="pone.0135004.t002" ref-type="table">Table 2</xref> present the comparison between the DTM-Padé and numerical solution, based on shooting technique, results for torque values at lower (<italic>dg</italic>(0)/<italic>dη</italic>) and upper disk (<italic>dg</italic>(1)/<italic>dη</italic>) for various values of squeeze Reynolds number (<italic>N</italic><sub>1</sub>) and axial magnetic force number (<italic>N</italic><sub>2</sub>).</p><table-wrap id="pone.0135004.t001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0135004.t001</object-id><label>Table 1</label><caption><title>Torque values at lower (<italic>dg</italic>(0)/<italic>dη</italic>) and upper disk (<italic>dg</italic>(1)/<italic>dη</italic>) when <italic>N</italic><sub>2</sub> = 1, <italic>N</italic><sub>3</sub> = 0.5 and <italic>Bt</italic> = 6 for various <italic>N</italic><sub>1</sub>.</title></caption><alternatives><graphic id="pone.0135004.t001g" xlink:href="pone.0135004.t001"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col></colgroup><thead><tr><th rowspan="2" align="center" colspan="1"><italic>N</italic><sub>1</sub></th><th align="center" rowspan="1" colspan="1"></th><th align="center" rowspan="1" colspan="1"><inline-formula id="pone.0135004.e051"><alternatives><graphic id="pone.0135004.e051g" xlink:href="pone.0135004.e051"></graphic><mml:math id="M51"><mml:mrow><mml:mfrac><mml:mrow><mml:mi mathvariant="bold-italic">d</mml:mi><mml:mi mathvariant="bold-italic">g</mml:mi><mml:mrow><mml:mo mathvariant="bold">(</mml:mo><mml:mn mathvariant="bold">0</mml:mn><mml:mo mathvariant="bold">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi mathvariant="bold-italic">d</mml:mi><mml:mi mathvariant="bold-italic">η</mml:mi></mml:mrow></mml:mfrac></mml:mrow></mml:math></alternatives></inline-formula></th><th align="center" rowspan="1" colspan="1"></th><th align="center" rowspan="1" colspan="1"><inline-formula id="pone.0135004.e052"><alternatives><graphic id="pone.0135004.e052g" xlink:href="pone.0135004.e052"></graphic><mml:math id="M52"><mml:mrow><mml:mfrac><mml:mrow><mml:mi mathvariant="bold-italic">d</mml:mi><mml:mi mathvariant="bold-italic">g</mml:mi><mml:mrow><mml:mo mathvariant="bold">(</mml:mo><mml:mn mathvariant="bold">1</mml:mn><mml:mo mathvariant="bold">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi mathvariant="bold-italic">d</mml:mi><mml:mi mathvariant="bold-italic">η</mml:mi></mml:mrow></mml:mfrac></mml:mrow></mml:math></alternatives></inline-formula></th></tr><tr><th align="center" rowspan="1" colspan="1">DTM-Padé Result</th><th align="center" rowspan="1" colspan="1">Numerical Result</th><th align="center" rowspan="1" colspan="1">DTM-Padé Result</th><th align="center" rowspan="1" colspan="1">Numerical Result</th></tr></thead><tbody><tr><td align="center" rowspan="1" colspan="1">0.1</td><td align="char" char="." rowspan="1" colspan="1">− 1.08963495</td><td align="char" char="." rowspan="1" colspan="1">− 1.08963506</td><td align="char" char="." rowspan="1" colspan="1">− 0.95987351</td><td align="char" char="." rowspan="1" colspan="1">− 0.95987349</td></tr><tr><td align="center" rowspan="1" colspan="1">0.2</td><td align="char" char="." rowspan="1" colspan="1">− 1.17203735</td><td align="char" char="." rowspan="1" colspan="1">− 1.17203765</td><td align="char" char="." rowspan="1" colspan="1">− 0.93844853</td><td align="char" char="." rowspan="1" colspan="1">− 0.93844830</td></tr><tr><td align="center" rowspan="1" colspan="1">0.3</td><td align="char" char="." rowspan="1" colspan="1">− 1.25013614</td><td align="char" char="." rowspan="1" colspan="1">− 1.25013649</td><td align="char" char="." rowspan="1" colspan="1">− 0.92615659</td><td align="char" char="." rowspan="1" colspan="1">− 0.92615609</td></tr><tr><td align="center" rowspan="1" colspan="1">0.5</td><td align="char" char="." rowspan="1" colspan="1">− 1.39797802</td><td align="char" char="." rowspan="1" colspan="1">− 1.39797797</td><td align="char" char="." rowspan="1" colspan="1">− 0.91295668</td><td align="char" char="." rowspan="1" colspan="1">− 0.91295593</td></tr><tr><td align="center" rowspan="1" colspan="1">1</td><td align="char" char="." rowspan="1" colspan="1">− 1.73306809</td><td align="char" char="." rowspan="1" colspan="1">− 1.73306821</td><td align="char" char="." rowspan="1" colspan="1">− 0.89280560</td><td align="char" char="." rowspan="1" colspan="1">− 0.89280536</td></tr><tr><td align="center" rowspan="1" colspan="1">2</td><td align="char" char="." rowspan="1" colspan="1">− 2.28925264</td><td align="char" char="." rowspan="1" colspan="1">− 2.28925762</td><td align="char" char="." rowspan="1" colspan="1">− 0.83902287</td><td align="char" char="." rowspan="1" colspan="1">− 0.83902117</td></tr></tbody></table></alternatives></table-wrap><table-wrap id="pone.0135004.t002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0135004.t002</object-id><label>Table 2</label><caption><title>Torque values at lower (<italic>dg</italic>(0)/<italic>dη</italic>) and upper disk (<italic>dg</italic>(1)/<italic>dη</italic>) <italic>N</italic><sub>1</sub> = 1, <italic>N</italic><sub>3</sub> = 0.5 and <italic>Bt</italic> = 6 for various <italic>N</italic><sub>2</sub>.</title></caption><alternatives><graphic id="pone.0135004.t002g" xlink:href="pone.0135004.t002"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col></colgroup><thead><tr><th rowspan="2" align="center" colspan="1"><italic>N</italic><sub>2</sub></th><th colspan="2" align="center" rowspan="1"><inline-formula id="pone.0135004.e053"><alternatives><graphic id="pone.0135004.e053g" xlink:href="pone.0135004.e053"></graphic><mml:math id="M53"><mml:mrow><mml:mfrac><mml:mrow><mml:mi mathvariant="bold-italic">d</mml:mi><mml:mi mathvariant="bold-italic">g</mml:mi><mml:mrow><mml:mo mathvariant="bold">(</mml:mo><mml:mn mathvariant="bold">0</mml:mn><mml:mo mathvariant="bold">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi mathvariant="bold-italic">d</mml:mi><mml:mi mathvariant="bold-italic">η</mml:mi></mml:mrow></mml:mfrac></mml:mrow></mml:math></alternatives></inline-formula></th><th colspan="2" align="center" rowspan="1"><inline-formula id="pone.0135004.e054"><alternatives><graphic id="pone.0135004.e054g" xlink:href="pone.0135004.e054"></graphic><mml:math id="M54"><mml:mrow><mml:mfrac><mml:mrow><mml:mi mathvariant="bold-italic">d</mml:mi><mml:mi mathvariant="bold-italic">g</mml:mi><mml:mrow><mml:mo mathvariant="bold">(</mml:mo><mml:mn mathvariant="bold">1</mml:mn><mml:mo mathvariant="bold">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi mathvariant="bold-italic">d</mml:mi><mml:mi mathvariant="bold-italic">η</mml:mi></mml:mrow></mml:mfrac></mml:mrow></mml:math></alternatives></inline-formula></th></tr><tr><th align="center" rowspan="1" colspan="1">DTM-Padé Result</th><th align="center" rowspan="1" colspan="1">Numerical Result</th><th align="center" rowspan="1" colspan="1">DTM-Padé Result</th><th align="center" rowspan="1" colspan="1">Numerical Result</th></tr></thead><tbody><tr><td align="center" rowspan="1" colspan="1">0</td><td align="char" char="." rowspan="1" colspan="1">− 1.77239023</td><td align="char" char="." rowspan="1" colspan="1">− 1.77238909</td><td align="char" char="." rowspan="1" colspan="1">− 0.60923829</td><td align="char" char="." rowspan="1" colspan="1">− 0.60923369</td></tr><tr><td align="center" rowspan="1" colspan="1">1</td><td align="char" char="." rowspan="1" colspan="1">− 1.73306809</td><td align="char" char="." rowspan="1" colspan="1">− 1.73306821</td><td align="char" char="." rowspan="1" colspan="1">− 0.89280560</td><td align="char" char="." rowspan="1" colspan="1">− 0.89280536</td></tr><tr><td align="center" rowspan="1" colspan="1">2</td><td align="char" char="." rowspan="1" colspan="1">− 1.62765012</td><td align="char" char="." rowspan="1" colspan="1">− 1.62765091</td><td align="char" char="." rowspan="1" colspan="1">− 1.91333401</td><td align="char" char="." rowspan="1" colspan="1">− 1.91333213</td></tr><tr><td align="center" rowspan="1" colspan="1">3</td><td align="char" char="." rowspan="1" colspan="1">− 1.48539942</td><td align="char" char="." rowspan="1" colspan="1">− 1.48539918</td><td align="char" char="." rowspan="1" colspan="1">− 3.92694971</td><td align="char" char="." rowspan="1" colspan="1">− 3.92694747</td></tr><tr><td align="center" rowspan="1" colspan="1">4</td><td align="char" char="." rowspan="1" colspan="1">− 1.34565346</td><td align="char" char="." rowspan="1" colspan="1">− 1.34565395</td><td align="char" char="." rowspan="1" colspan="1">− 6.89037665</td><td align="char" char="." rowspan="1" colspan="1">− 6.89038054</td></tr><tr><td align="center" rowspan="1" colspan="1">5</td><td align="char" char="." rowspan="1" colspan="1">− 1.23080507</td><td align="char" char="." rowspan="1" colspan="1">− 1.23080496</td><td align="char" char="." rowspan="1" colspan="1">− 10.44130841</td><td align="char" char="." rowspan="1" colspan="1">− 10.44131517</td></tr></tbody></table></alternatives></table-wrap><p>Figs <xref rid="pone.0135004.g004" ref-type="fig">4</xref>, <xref rid="pone.0135004.g005" ref-type="fig">5</xref>, <xref rid="pone.0135004.g006" ref-type="fig">6</xref> and <xref rid="pone.0135004.g007" ref-type="fig">7</xref> display the effects of magnetic Reynolds number (<italic>Re</italic><sub><italic>m</italic></sub>) on the axial and tangential velocity distributions and induced magnetic field distributions (<italic>f</italic>, <italic>g</italic>, <italic>m</italic>, <italic>n</italic>). The magnetic Reynolds number defines the ratio of the fluid flux to the magnetic diffusivity. This parameter therefore is instrumental in determining the diffusion of magnetic field along streamlines and is analogous to the classical Reynolds number in viscous hydrodynamics, the latter controlling the vorticity diffusion along the streamlines. When <italic>N</italic><sub>1</sub> is large this implies fast vertical velocity of the upper disk and vice versa for small values of this squeezing parameter. The variation in <italic>Re</italic><sub><italic>m</italic></sub> has almost little effect on the axial velocity distribution. As it is obvious from <xref rid="pone.0135004.g005" ref-type="fig">Fig 5</xref>, the maximum values of <italic>g</italic> arises at the lower disk (<italic>η</italic> = 0) for all cases. In other word, the azimuthal velocity of the fluid decreases as we move from the lower disk towards the upper disk. In addition, the tangential velocity distribution is a decreasing function of magnetic Reynolds number. With an increase in <italic>Re</italic><sub><italic>m</italic></sub> from 1, 2 through to the maximum value of 10, there is a strong decrease in the axial induced magnetic field component (<italic>m</italic>) and tangential induced magnetic field component (<italic>n</italic>). For all cases the maximum values of <italic>m</italic> and <italic>n</italic> arise at the upper disk (<italic>η</italic> = 1). It is to certify that for higher values of <italic>Re</italic><sub><italic>m</italic></sub> the squeeze film must possess much higher electrical conductivities.</p><fig id="pone.0135004.g004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0135004.g004</object-id><label>Fig 4</label><caption><title>Effect of magnetic Reynolds number (<italic>Re</italic><sub><italic>m</italic></sub>) on the axial velocity distribution (<italic>f</italic>) when <italic>N</italic><sub>2</sub> = 1, <italic>N</italic><sub>3</sub> = 0.5 and <italic>Bt</italic> = 6.</title></caption><graphic xlink:href="pone.0135004.g004"></graphic></fig><fig id="pone.0135004.g005" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0135004.g005</object-id><label>Fig 5</label><caption><title>Effect of magnetic Reynolds number (<italic>Re</italic><sub><italic>m</italic></sub>) on the tangential velocity distribution (<italic>g</italic>) when <italic>N</italic><sub>2</sub> = 1, <italic>N</italic><sub>3</sub> = 0.5 and <italic>Bt</italic> = 6.</title></caption><graphic xlink:href="pone.0135004.g005"></graphic></fig><fig id="pone.0135004.g006" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0135004.g006</object-id><label>Fig 6</label><caption><title>Effect of magnetic Reynolds number (<italic>Re</italic><sub><italic>m</italic></sub>) on the axial induced magnetic field distribution (<italic>m</italic>) when <italic>N</italic><sub>2</sub> = 1, <italic>N</italic><sub>3</sub> = 0.5 and <italic>Bt</italic> = 6.</title></caption><graphic xlink:href="pone.0135004.g006"></graphic></fig><fig id="pone.0135004.g007" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0135004.g007</object-id><label>Fig 7</label><caption><title>Effect of magnetic Reynolds number (<italic>Re</italic><sub><italic>m</italic></sub>) on the tangential induced magnetic field distribution (<italic>n</italic>) when <italic>N</italic><sub>2</sub> = 1, <italic>N</italic><sub>3</sub> = 0.5 and <italic>Bt</italic> = 6.</title></caption><graphic xlink:href="pone.0135004.g007"></graphic></fig><p>The effects of dimensionless axial magnetic force parameter (<italic>N</italic><sub>2</sub>) on the axial and tangential velocity distributions and induced magnetic field distributions (<italic>f</italic>, <italic>g</italic>, <italic>m</italic>, <italic>n</italic>) are illustrated in Figs <xref rid="pone.0135004.g008" ref-type="fig">8</xref>, <xref rid="pone.0135004.g009" ref-type="fig">9</xref>, <xref rid="pone.0135004.g010" ref-type="fig">10</xref> and <xref rid="pone.0135004.g011" ref-type="fig">11</xref>. Axial velocity is generated in the two-disk system from the vertical motion of the upper disk and the radial flux far from the axis of rotation in the vicinity of the lower disk. When <italic>N</italic><sub>2</sub> is large this means fast vertical velocity of the upper disk and vice versa for small <italic>N</italic><sub>2</sub>. It can be seen that an increase in <italic>N</italic><sub>2</sub> from 0 to 5 induces a significant decrease in the axial velocity between the disks. The tangential velocity distribution (<italic>g</italic>) is also directly affected by an increase in squeeze Reynolds number. As the axial component of the magnetic force <italic>N</italic><sub>2</sub> increases, the azimuthal velocity <italic>g</italic> increases. Hence, the axial magnetic force <italic>N</italic><sub>2</sub> can be used to increase the azimuthal velocity of the fluid. The increase of <italic>g</italic> with <italic>N</italic><sub>2</sub> agrees with the results obtained by Elshekh and Elhady [<xref rid="pone.0135004.ref024" ref-type="bibr">24</xref>] and Hamza [<xref rid="pone.0135004.ref043" ref-type="bibr">43</xref>]. The effect of the magnetic force <italic>N</italic><sub>2</sub> will be dominant Ω<sub>1</sub>. As the results present, an increase in the dimensionless axial magnetic force parameter causes to decrease the axial induced magnetic field distribution (<italic>m</italic>) and increase in the tangential induced magnetic field distribution (<italic>n</italic>). The tangential velocity profile and tangential induced magnetic field also become significantly more curved with higher <italic>N</italic><sub>2</sub> values. The current results are an expected result since the normal component of the induced magnetic field between the two disks must increase with the increase of the normal component of the external applied magnetic field.</p><fig id="pone.0135004.g008" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0135004.g008</object-id><label>Fig 8</label><caption><title>Effect of axial magnetic force number (<italic>N</italic><sub>2</sub>) on the axial velocity distribution (<italic>f</italic>) when <italic>N</italic><sub>1</sub> = 1, <italic>N</italic><sub>3</sub> = 0.5 and <italic>Bt</italic> = 6.</title></caption><graphic xlink:href="pone.0135004.g008"></graphic></fig><fig id="pone.0135004.g009" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0135004.g009</object-id><label>Fig 9</label><caption><title>Effect of axial magnetic force number (<italic>N</italic><sub>2</sub>) on the tangential velocity distribution (<italic>g</italic>) when <italic>N</italic><sub>1</sub> = 1, <italic>N</italic><sub>3</sub> = 0.5 and <italic>Bt</italic> = 6.</title></caption><graphic xlink:href="pone.0135004.g009"></graphic></fig><fig id="pone.0135004.g010" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0135004.g010</object-id><label>Fig 10</label><caption><title>Effect of axial magnetic force number (<italic>N</italic><sub>2</sub>) on the axial induced magnetic field distribution (<italic>m</italic>) when <italic>N</italic><sub>1</sub> = 1, <italic>N</italic><sub>3</sub> = 0.5 and <italic>Bt</italic> = 6.</title></caption><graphic xlink:href="pone.0135004.g010"></graphic></fig><fig id="pone.0135004.g011" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0135004.g011</object-id><label>Fig 11</label><caption><title>Effect of axial magnetic force number (<italic>N</italic><sub>2</sub>) on the tangential induced magnetic field distribution (<italic>n</italic>) when <italic>N</italic><sub>1</sub> = 1, <italic>N</italic><sub>3</sub> = 0.5 and <italic>Bt</italic> = 6.</title></caption><graphic xlink:href="pone.0135004.g011"></graphic></fig><p>Figs <xref rid="pone.0135004.g012" ref-type="fig">12</xref> and <xref rid="pone.0135004.g013" ref-type="fig">13</xref> depict the effect of dimensionless tangential magnetic force parameter (<italic>N</italic><sub>3</sub>) on the tangential velocity distribution (<italic>g</italic>) and tangential induced magnetic field distribution (<italic>n</italic>), respectively. The tangential velocity distribution slightly increases with the increases in the dimensionless tangential magnetic force parameter from 0.5 to 5. Increasing the dimensionless tangential magnetic force parameter causes to decreases the tangential induced magnetic field distribution.</p><fig id="pone.0135004.g012" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0135004.g012</object-id><label>Fig 12</label><caption><title>Effect of tangential magnetic force number (<italic>N</italic><sub>3</sub>) on the tangential velocity distribution (<italic>g</italic>) when <italic>N</italic><sub>1</sub> = <italic>N</italic><sub>2</sub> = 1, and <italic>Bt</italic> = 6.</title></caption><graphic xlink:href="pone.0135004.g012"></graphic></fig><fig id="pone.0135004.g013" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0135004.g013</object-id><label>Fig 13</label><caption><title>Effect of tangential magnetic force number (<italic>N</italic><sub>3</sub>) on the tangential induced magnetic field distribution (<italic>n</italic>) when <italic>N</italic><sub>1</sub> = <italic>N</italic><sub>2</sub> = 1, and <italic>Bt</italic> = 6.</title></caption><graphic xlink:href="pone.0135004.g013"></graphic></fig></sec><sec sec-type="conclusions" id="sec007"><title>Conclusions</title><p>The present study has displayed novel DTM-Padé solution, the combination of differential transform method and Padé approximation, for two nonlinear magneto-hydrodynamic squeeze film problems. Applications of the study include in automotive magneto-rheological shock absorbers, novel aircraft landing gear systems and biological prosthetics. The transformed dimensionless equations have been formulated and solved with robust boundary conditions. Exceptional stability and convergence characteristics have been demonstrated with the DTM. The physical key parameters emerging have been investigated graphically in detail including dimensionless axial magnetic force strength parameter, dimensionless tangential magnetic force strength parameter and magnetic Reynolds number. The results illustrated that the tangential velocity distribution is a decreasing function of magnetic Reynolds number. Further, the tangential velocity distribution is directly affected by an increase in squeeze Reynolds number. In addition, an increase in the dimensionless axial magnetic force parameter causes to decrease the axial induced magnetic field distribution and increase in the tangential induced magnetic field distribution. Moreover, increasing the dimensionless tangential magnetic force parameter causes to decreases the tangential induced magnetic field distribution.</p></sec></body><back><ack><p>EM acknowledges support from the National Research foundation of South Africa.</p></ack><ref-list><title>References</title><ref id="pone.0135004.ref001"><label>1</label><mixed-citation publication-type="book"><name><surname>Davidson</surname><given-names>P A</given-names></name>. <chapter-title>An introduction to magnetohydrodynamics</chapter-title><publisher-name>Cambridge University Press</publisher-name>, <year>2001</year>.</mixed-citation></ref><ref id="pone.0135004.ref002"><label>2</label><mixed-citation publication-type="journal"><name><surname>Rashidi</surname><given-names>M M</given-names></name>, <name><surname>Rostami</surname><given-names>B</given-names></name>, <name><surname>Freidoonimehr</surname><given-names>N</given-names></name>, <name><surname>Abbasbandy</surname><given-names>S</given-names></name>.<article-title>Free convective heat and mass transfer for MHD fluid flow over a permeable vertical stretching sheet in the presence of the radiation and buoyancy effects</article-title>. <source>Ain Shams Engineering Journal</source><volume>5</volume> (<year>2014</year>) <fpage>901</fpage>–<lpage>912</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref003"><label>3</label><mixed-citation publication-type="journal"><name><surname>Khan</surname><given-names>A</given-names></name>, <name><surname>Khan</surname><given-names>I</given-names></name>, <name><surname>Ali</surname><given-names>F</given-names></name>, <name><surname>ulhaq</surname><given-names>S</given-names></name>, <name><surname>Shafie</surname><given-names>S</given-names></name>, <article-title>Effects of Wall Shear Stress on Unsteady MHD Conjugate Flow in a Porous Medium with Ramped Wall Temperature</article-title>, <source>PLOS ONE</source><volume>9</volume> (<year>2014</year>) <fpage>e90280</fpage><comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1371/journal.pone.0090280">10.1371/journal.pone.0090280</ext-link></comment><pub-id pub-id-type="pmid">24621775</pub-id></mixed-citation></ref><ref id="pone.0135004.ref004"><label>4</label><mixed-citation publication-type="journal"><name><surname>Hussanan</surname><given-names>A</given-names></name>, <name><surname>Ismail</surname><given-names>Z</given-names></name>, <name><surname>Khan</surname><given-names>I</given-names></name>, <name><surname>Hussein</surname><given-names>A</given-names></name>, <name><surname>Shafie</surname><given-names>S</given-names></name>, <article-title>Unsteady boundary layer MHD free convection flow in a porous medium with constant mass diffusion and Newtonian heating</article-title>, <source>The European Physical Journal Plus</source><volume>129</volume> (<year>2014</year>) <fpage>1</fpage>–<lpage>16</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref005"><label>5</label><mixed-citation publication-type="journal"><name><surname>Samiulhaq</surname></name>, <name><surname>Ahmad</surname><given-names>S</given-names></name>, <name><surname>Vieru</surname><given-names>D</given-names></name>, <name><surname>Khan</surname><given-names>I</given-names></name>, <name><surname>Shafie</surname><given-names>S</given-names></name>, <article-title>Unsteady Magnetohydrodynamic Free Convection Flow of a Second Grade Fluid in a Porous Medium with Ramped Wall Temperature</article-title>, <source>PLOS ONE</source><volume>9</volume> (<year>2014</year>) <fpage>e88766</fpage><comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1371/journal.pone.0088766">10.1371/journal.pone.0088766</ext-link></comment><pub-id pub-id-type="pmid">24785147</pub-id></mixed-citation></ref><ref id="pone.0135004.ref006"><label>6</label><mixed-citation publication-type="journal"><name><surname>Khan</surname><given-names>A</given-names></name>, <name><surname>Khan</surname><given-names>I</given-names></name>, <name><surname>Ali</surname><given-names>F</given-names></name>, <name><surname>Shafie</surname><given-names>S</given-names></name>, <article-title>Effects of Wall Shear Stress on MHD Conjugate Flow over an Inclined Plate in a Porous Medium with Ramped Wall Temperature</article-title>, <source>Mathematical Problems in Engineering</source><volume>2014</volume> (<year>2014</year>) <fpage>15</fpage> pages.</mixed-citation></ref><ref id="pone.0135004.ref007"><label>7</label><mixed-citation publication-type="journal"><name><surname>Khalid</surname><given-names>A</given-names></name>, <name><surname>Khan</surname><given-names>I</given-names></name>, <name><surname>Khan</surname><given-names>A</given-names></name>, <name><surname>Shafie</surname><given-names>S</given-names></name>, <article-title>Unsteady MHD free convection flow of Casson fluid past over an oscillating vertical plate embedded in a porous medium</article-title>, <source>Engineering Science and Technology, an International Journal</source>.</mixed-citation></ref><ref id="pone.0135004.ref008"><label>8</label><mixed-citation publication-type="journal"><name><surname>Park</surname><given-names>E</given-names></name>, <name><surname>Min</surname><given-names>K-W</given-names></name>, <name><surname>Lee</surname><given-names>S-K</given-names></name>, <name><surname>Lee</surname><given-names>S-H</given-names></name>, <name><surname>Lee</surname><given-names>H-J</given-names></name>, <name><surname>Moon</surname><given-names>S-J</given-names></name>, <etal>et al</etal>, <article-title>Real-time hybrid test on a semi-actively controlled building structure equipped with full-scale MR dampers</article-title>, <source>Journal of Intelligent Material Systems and Structures</source><volume>21</volume> (<year>2010</year>) <fpage>1831</fpage>–<lpage>1850</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref009"><label>9</label><mixed-citation publication-type="book"><name><surname>Anderson</surname><given-names>C G</given-names></name>, <chapter-title>Magnetic squeeze-films: design proposals for a magnetic-repulsion enhanced hydrostatic bearing for the type N2 wave-power buoy</chapter-title>, <source>Technical Report for Norwegian Institute of Technology-Trondheim</source>, <publisher-name>Dept. Mechanical Engineering, Edinburgh University</publisher-name>, <publisher-loc>Scotland</publisher-loc> (<year>1983</year>).</mixed-citation></ref><ref id="pone.0135004.ref010"><label>10</label><mixed-citation publication-type="journal"><name><surname>Chernyakova</surname><given-names>Y M</given-names></name>, <name><surname>Pinchuk</surname><given-names>L S</given-names></name>, <name><surname>Lobanovskii</surname><given-names>L. S</given-names></name>, <name><surname>Drozd</surname><given-names>E S</given-names></name>, <name><surname>Chizhik</surname><given-names>S A</given-names></name>, <article-title>Structure and magnetic susceptibility of lubricating fluids in synovial joints</article-title>, <source>Journal of Friction and Wear</source><volume>32</volume> (<year>2011</year>) <fpage>54</fpage>–<lpage>60</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref011"><label>11</label><mixed-citation publication-type="journal"><name><surname>Liu</surname><given-names>B</given-names></name>, <name><surname>Li</surname><given-names>W H</given-names></name>, <name><surname>Kosasih</surname><given-names>P B</given-names></name>, <name><surname>Zhang</surname><given-names>X Z</given-names></name>, <article-title>Development of an MR-brake-based haptic device</article-title>, <source>Smart Materials and Structures</source><volume>15</volume> (<year>2006</year>) <fpage>1960</fpage>.</mixed-citation></ref><ref id="pone.0135004.ref012"><label>12</label><mixed-citation publication-type="journal"><name><surname>Chandra</surname><given-names>P</given-names></name>, <name><surname>Sinha</surname><given-names>P</given-names></name>, <name><surname>Kumar</surname><given-names>D</given-names></name>, <article-title>Ferrofluid lubrication of a journal bearing considering cavitation</article-title>, <source>Tribology Transactions</source><volume>35</volume> (<year>1992</year>) <fpage>163</fpage>–<lpage>169</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref013"><label>13</label><mixed-citation publication-type="journal"><name><surname>Song</surname><given-names>W L</given-names></name>, <name><surname>Choi</surname><given-names>S-B</given-names></name>, <name><surname>Choi</surname><given-names>J Y</given-names></name>, <name><surname>Lee</surname><given-names>C H</given-names></name>, <article-title>Wear and friction characteristics of magnetorheological fluid under magnetic field activation</article-title>, <source>Tribology Transactions</source><volume>54</volume> (<year>2011</year>) <fpage>616</fpage>–<lpage>624</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref014"><label>14</label><mixed-citation publication-type="journal"><name><surname>Durán</surname><given-names>J D G</given-names></name>, <name><surname>Iglesias</surname><given-names>G R</given-names></name>, <name><surname>Delgado</surname><given-names>A V</given-names></name>, <name><surname>Ruiz-Morón</surname><given-names>L F</given-names></name>, <name><surname>Insa</surname><given-names>J</given-names></name>, <name><surname>González-Caballero</surname><given-names>F</given-names></name>, <article-title>Stability and Flow Behavior of a Magnetorheological Lubricant in a Magnetic Shock Absorber</article-title>, <source>Tribology Transactions</source><volume>51</volume> (<year>2008</year>) <fpage>271</fpage>–<lpage>277</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref015"><label>15</label><mixed-citation publication-type="journal"><name><surname>Huang</surname><given-names>W</given-names></name>, <name><surname>Shen</surname><given-names>C</given-names></name>, <name><surname>Liao</surname><given-names>S</given-names></name>, <name><surname>Wang</surname><given-names>X</given-names></name>, <article-title>Study on the Ferrofluid Lubrication with an External Magnetic Field</article-title>, <source>Tribology Letters</source><volume>41</volume> (<year>2011</year>) <fpage>145</fpage>–<lpage>151</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref016"><label>16</label><mixed-citation publication-type="journal"><name><surname>Stolarski</surname><given-names>T A</given-names></name>, <name><surname>Makida</surname><given-names>Y</given-names></name>, <article-title>Performance of Lubricated Sliding Contact in Magnetic Field</article-title>, <source>Tribology Letters</source><volume>46</volume> (<year>2012</year>) <fpage>113</fpage>–<lpage>121</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref017"><label>17</label><mixed-citation publication-type="journal"><name><surname>Anwar</surname><given-names>M I</given-names></name>, <name><surname>Rodkiewicz</surname><given-names>C M</given-names></name>, <article-title>Nonuniform Magnetic Field Effects in MHD Slider Bearing</article-title>, <source>Journal of Tribology</source><volume>94</volume> (<year>1972</year>) <fpage>101</fpage>–<lpage>105</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref018"><label>18</label><mixed-citation publication-type="journal"><name><surname>Kargulewicz</surname><given-names>M</given-names></name>, <name><surname>Iordanoff</surname><given-names>I</given-names></name>, <name><surname>Marrero</surname><given-names>V</given-names></name>, <name><surname>Tichy</surname><given-names>J</given-names></name>, <article-title>Modeling of Magnetorheological Fluids by the Discrete Element Method</article-title>, <source>Journal of Tribology</source><volume>134</volume> (<year>2012</year>) <fpage>031706</fpage>–<lpage>031706</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref019"><label>19</label><mixed-citation publication-type="journal"><name><surname>Zueco</surname><given-names>J</given-names></name>, <name><surname>Bég</surname><given-names>O A</given-names></name>, <article-title>Network numerical analysis of hydromagnetic squeeze film flow dynamics between two parallel rotating disks with induced magnetic field effects</article-title>, <source>Tribology International</source><volume>43</volume> (<year>2010</year>) <fpage>532</fpage>–<lpage>543</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref020"><label>20</label><mixed-citation publication-type="journal"><name><surname>Zhu</surname><given-names>G</given-names></name>, <name><surname>Ingber</surname><given-names>M S</given-names></name>, <article-title>Accurate treatment of lubrication forces between rigid spheres in viscous fluids using a traction-corrected boundary element method</article-title>, <source>Engineering Analysis with Boundary Elements</source><volume>33</volume> (<year>2009</year>) <fpage>467</fpage>–<lpage>473</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref021"><label>21</label><mixed-citation publication-type="journal"><name><surname>Gertzos</surname><given-names>K P</given-names></name>, <name><surname>Nikolakopoulos</surname><given-names>P G</given-names></name>, <name><surname>Papadopoulos</surname><given-names>C A</given-names></name>, <article-title>CFD analysis of journal bearing hydrodynamic lubrication by Bingham lubricant</article-title>, <source>Tribology International</source><volume>41</volume> (<year>2008</year>) <fpage>1190</fpage>–<lpage>1204</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref022"><label>22</label><mixed-citation publication-type="journal"><name><surname>Asadi</surname><given-names>M</given-names></name>, <name><surname>Rashidi</surname><given-names>M M</given-names></name>, <name><surname>Bég</surname><given-names>T A</given-names></name>, <name><surname>Bég</surname><given-names>O A</given-names></name>, <article-title>Homotopy analysis of transient magneto-bio-fluid dynamics of micropolar squeeze film in a porous medium: A model for magneto-bio-rheological lubrication</article-title>. <source>Journal of Mechanics in Medicine and Biology</source><volume>12</volume> (<year>2012</year>) <fpage>1250051</fpage>.</mixed-citation></ref><ref id="pone.0135004.ref023"><label>23</label><mixed-citation publication-type="journal"><name><surname>Moghani</surname><given-names>T</given-names></name>, <name><surname>Butler</surname><given-names>J P</given-names></name>, <name><surname>Lin</surname><given-names>J L-W</given-names></name>, <name><surname>Loring</surname><given-names>S H</given-names></name>.<article-title>Finite element simulation of elastohydrodynamic lubrication of soft biological tissues</article-title>. <source>Computers & Structures</source><volume>85</volume> (<year>2007</year>) <fpage>1114</fpage>–<lpage>1120</lpage>.<pub-id pub-id-type="pmid">18037975</pub-id></mixed-citation></ref><ref id="pone.0135004.ref024"><label>24</label><mixed-citation publication-type="journal"><name><surname>Elshekh</surname><given-names>S S</given-names></name>, <name><surname>Abd Elhady</surname><given-names>M K</given-names></name>.<article-title>Fluid film squeezed between two rotating disks in the presence of a magnetic field</article-title>. <source>International Journal of Engineering Science</source><volume>34</volume> (<year>1996</year>) <fpage>1183</fpage>–<lpage>1195</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref025"><label>25</label><mixed-citation publication-type="journal"><name><surname>Gul</surname><given-names>T</given-names></name>, <name><surname>Islam</surname><given-names>S</given-names></name>, <name><surname>Shah</surname><given-names>R A</given-names></name>, <name><surname>Khan</surname><given-names>I</given-names></name>, <name><surname>Khalid</surname><given-names>A</given-names></name>, <name><surname>Shafie</surname><given-names>S</given-names></name>. <article-title>Heat Transfer Analysis of MHD Thin Film Flow of an Unsteady Second Grade Fluid Past a Vertical Oscillating Belt</article-title>. <source>PLOS ONE</source><volume>9</volume> (<year>2014</year>) <fpage>e103843</fpage><comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1371/journal.pone.0103843">10.1371/journal.pone.0103843</ext-link></comment><pub-id pub-id-type="pmid">25383797</pub-id></mixed-citation></ref><ref id="pone.0135004.ref026"><label>26</label><mixed-citation publication-type="journal"><name><surname>Gul</surname><given-names>T</given-names></name>, <name><surname>Islam</surname><given-names>S</given-names></name>, <name><surname>Shah</surname><given-names>R A</given-names></name>, <name><surname>Khan</surname><given-names>I</given-names></name>, <name><surname>Shafie</surname><given-names>S</given-names></name>, <name><surname>Khan</surname><given-names>M A</given-names></name>. <article-title>Analysis of thin film flow over a vertical oscillating belt with a second grade fluid</article-title>. <source>Engineering Science and Technology, an International Journal</source>.</mixed-citation></ref><ref id="pone.0135004.ref027"><label>27</label><mixed-citation publication-type="journal"><name><surname>Fatma</surname><given-names>A</given-names></name>.<article-title>Applications of differential transform method to differential-algebraic equations</article-title>. <source>Applied Mathematics and Computation</source><volume>152</volume> (<year>2004</year>) <fpage>649</fpage>–<lpage>657</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref028"><label>28</label><mixed-citation publication-type="journal"><name><surname>Fatma</surname><given-names>A</given-names></name>. <article-title>Solutions of the system of differential equations by differential transform method</article-title>. <source>Applied Mathematics and Computation</source><volume>147</volume> (<year>2004</year>) <fpage>547</fpage>–<lpage>567</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref029"><label>29</label><mixed-citation publication-type="book"><name><surname>Liao</surname><given-names>S</given-names></name>. <chapter-title>Beyond Perturbation: Introduction to the Homotopy Analysis Method</chapter-title><publisher-name>Taylor & Francis</publisher-name>, <year>2003</year>.</mixed-citation></ref><ref id="pone.0135004.ref030"><label>30</label><mixed-citation publication-type="book"><name><surname>Liao</surname><given-names>S</given-names></name>. <chapter-title>Advances in the Homotopy Analysis Method</chapter-title>, <publisher-name>World Scientific Publishing Company. Shanghai Jiao Tong University</publisher-name>, <publisher-loc>China</publisher-loc>, <year>2014</year>.</mixed-citation></ref><ref id="pone.0135004.ref031"><label>31</label><mixed-citation publication-type="journal"><name><surname>Dehghan</surname><given-names>M</given-names></name>, <name><surname>Rahmani</surname><given-names>Y</given-names></name>, <name><surname>Ganji</surname><given-names>D D</given-names></name>, <name><surname>Saedodin</surname><given-names>S</given-names></name>, <name><surname>Valipour</surname><given-names>M S</given-names></name>, <name><surname>Rashidi</surname><given-names>S</given-names></name>,.<article-title>Convection–radiation heat transfer in solar heat exchangers filled with a porous medium: Homotopy perturbation method versus numerical analysis</article-title>. <source>Renewable Energy</source><volume>74</volume> (<year>2015</year>) <fpage>448</fpage>–<lpage>455</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref032"><label>32</label><mixed-citation publication-type="journal"><name><surname>Sheikholeslami</surname><given-names>M</given-names></name>, <name><surname>Ganji</surname><given-names>D D</given-names></name>, <name><surname>Ashorynejad</surname><given-names>H R</given-names></name>, <name><surname>Rokni</surname><given-names>H B</given-names></name>. <article-title>Analytical investigation of Jeffery-Hamel flow with high magnetic field and nanoparticle by Adomian decomposition method</article-title>. <source>Applied Mathematics and Mechanics</source><volume>33</volume> (<year>2012</year>) <fpage>25</fpage>–<lpage>36</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref033"><label>33</label><mixed-citation publication-type="journal"><name><surname>Rashidi</surname><given-names>M M</given-names></name>, <name><surname>Erfani</surname><given-names>E</given-names></name>. <article-title>A new analytical study of MHD stagnation-point flow in porous media with heat transfer</article-title>. <source>Computers and Fluids</source><volume>40</volume> (<year>2011</year>) <fpage>172</fpage>–<lpage>178</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref034"><label>34</label><mixed-citation publication-type="journal"><name><surname>Hamza</surname><given-names>E A</given-names></name>, <name><surname>Macdonald</surname><given-names>D A</given-names></name>. <article-title>A fluid film squeezed between two parallel plane surfaces</article-title>. <source>Journal of Fluid Mechanics</source><volume>109</volume> (<year>1981</year>) <fpage>147</fpage>–<lpage>160</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref035"><label>35</label><mixed-citation publication-type="journal"><name><surname>Ibrahim</surname><given-names>F N</given-names></name>. <article-title>Unsteady flow between two rotating discs with heat transfer</article-title>. <source>Journal of Physics D: Applied Physics</source><volume>24</volume> (<year>1991</year>) <fpage>1293</fpage>–<lpage>1299</lpage></mixed-citation></ref><ref id="pone.0135004.ref036"><label>36</label><mixed-citation publication-type="journal"><name><surname>Ibrahim</surname><given-names>F N</given-names></name>. <article-title>Unsteady flow between two rotating discs with heat transfer</article-title>. <source>Journal of Physics D: Applied Physics</source><volume>24</volume> (<year>1991</year>) <fpage>1293</fpage>.</mixed-citation></ref><ref id="pone.0135004.ref037"><label>37</label><mixed-citation publication-type="journal"><name><surname>Hughes</surname><given-names>W F</given-names></name>, <name><surname>Elco</surname><given-names>R A</given-names></name>. <article-title>Magnetohydrodynamic lubrication flow between parallel rotating disks</article-title>. <source>Journal of Fluid Mechanics</source><volume>13</volume> (<year>1962</year>) <fpage>21</fpage>–<lpage>32</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref038"><label>38</label><mixed-citation publication-type="book"><name><surname>Baker</surname><given-names>G A</given-names></name>.<chapter-title>Essential of padé approximants</chapter-title><publisher-name>Academic Press</publisher-name>, <publisher-loc>London</publisher-loc>, <year>1975</year>.</mixed-citation></ref><ref id="pone.0135004.ref039"><label>39</label><mixed-citation publication-type="book"><name><surname>Baker</surname><given-names>G A</given-names></name>, <name><surname>Graves-Morris</surname><given-names>P R</given-names></name>, <chapter-title>Padé approximants: Basic theory</chapter-title><publisher-name>Addison-Wesley</publisher-name>, <year>1981</year>.</mixed-citation></ref><ref id="pone.0135004.ref040"><label>40</label><mixed-citation publication-type="journal"><name><surname>Rashidi</surname><given-names>M M</given-names></name>. <article-title>The modified differential transform method for solving MHD boundary-layer equations</article-title>. <source>Computer Physics Communications</source><volume>180</volume> (<year>2009</year>) <fpage>2210</fpage>–<lpage>2217</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref041"><label>41</label><mixed-citation publication-type="journal"><name><surname>Rashidi</surname><given-names>M M</given-names></name>, <name><surname>Mohimanian</surname><given-names>Pour S A</given-names></name>. <article-title>A novel analytical solution of heat transfer of a micropolar fluid through a porous medium with radiation by DTM-Padé</article-title>. <source>Heat Transfer—Asian Research</source><volume>39</volume> (<year>2010</year>) <fpage>575</fpage>–<lpage>589</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref042"><label>42</label><mixed-citation publication-type="journal"><name><surname>Rashidi</surname><given-names>M M</given-names></name>, <name><surname>Freidoonimehr</surname><given-names>N</given-names></name>. <article-title>Series solutions for the flow in the vicinity of the equator of an MHD boundary-layer over a porous rotating sphere with heat transfer</article-title>. <source>Thermal Science</source><volume>18</volume> (<year>2012</year>) <fpage>S527</fpage>—<lpage>S537</lpage>.</mixed-citation></ref><ref id="pone.0135004.ref043"><label>43</label><mixed-citation publication-type="journal"><name><surname>Hamza</surname><given-names>E A</given-names></name>.<article-title>The magnetohydrodynamic effects on a fluid film squeezed between two rotating surfaces</article-title>. <source>Journal of Physics D: Applied Physics</source><volume>24</volume> (<year>1991</year>) <fpage>547</fpage>.</mixed-citation></ref></ref-list></back></pmc><affiliations><list><country><li>Afrique du Sud</li><li>Iran</li><li>République populaire de Chine</li></country><region><li>Gauteng</li></region><settlement><li>Johannesburg</li></settlement><orgName><li>Université du Witwatersrand</li></orgName></list><tree><country name="Iran"><noRegion><name sortKey="Rashidi, Mohammad Mehdi" sort="Rashidi, Mohammad Mehdi" uniqKey="Rashidi M" first="Mohammad Mehdi" last="Rashidi">Mohammad Mehdi Rashidi</name></noRegion><name sortKey="Freidoonimehr, Navid" sort="Freidoonimehr, Navid" uniqKey="Freidoonimehr N" first="Navid" last="Freidoonimehr">Navid Freidoonimehr</name><name sortKey="Rostami, Behnam" sort="Rostami, Behnam" uniqKey="Rostami B" first="Behnam" last="Rostami">Behnam Rostami</name></country><country name="République populaire de Chine"><noRegion><name sortKey="Rashidi, Mohammad Mehdi" sort="Rashidi, Mohammad Mehdi" uniqKey="Rashidi M" first="Mohammad Mehdi" last="Rashidi">Mohammad Mehdi Rashidi</name></noRegion></country><country name="Afrique du Sud"><region name="Gauteng"><name sortKey="Momoniat, Ebrahim" sort="Momoniat, Ebrahim" uniqKey="Momoniat E" first="Ebrahim" last="Momoniat">Ebrahim Momoniat</name></region></country></tree></affiliations></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Ticri/CIDE/explor/HapticV1/Data/Ncbi/Merge
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 003B52 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Ncbi/Merge/biblio.hfd -nk 003B52 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Ticri/CIDE
|area= HapticV1
|flux= Ncbi
|étape= Merge
|type= RBID
|clé= PMC:4534046
|texte= Study of Nonlinear MHD Tribological Squeeze Film at Generalized Magnetic Reynolds Numbers Using DTM
}}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Ncbi/Merge/RBID.i -Sk "pubmed:26267247" \
| HfdSelect -Kh $EXPLOR_AREA/Data/Ncbi/Merge/biblio.hfd \
| NlmPubMed2Wicri -a HapticV1
| This area was generated with Dilib version V0.6.23. |  |