==== Front Sci Rep Sci Rep Scientific Reports 2045-2322 Nature Publishing Group UK London 77996 10.1038/s41598-020-77996-x Article Dual solutions of nanomaterial flow comprising titanium alloy (Ti6Al4V) suspended in Williamson fluid through a thin moving needle with nonlinear thermal radiation: stability scrutinization Khan Umair 1 Zaib A. 2 Khan Ilyas ilyaskhan@tdtu.edu.vn 3 Nisar Kottakkaran Sooppy n.sooppy@psau.edu.sa 4 1 grid.442838.10000 0004 0609 4757Department of Mathematics and Social Sciences, Sukkur IBA University, Sukkur, 65200 Sindh Pakistan 2 grid.440529.e0000 0004 0607 3470Department of Mathematical Sciences, Federal Urdu University of Arts, Science and Technology, Gulshan-E-Iqbal, Karachi, 75300 Pakistan 3 grid.444812.f0000 0004 5936 4802Faculty of Mathematics and Statistics, Ton Duc Thang University, Ho Chi Minh City, Vietnam 4 grid.449553.aDepartment of Mathematics, College of Arts and Sciences, Prince Sattam bin Abdulaziz University, Wadi Aldawaser, 11991 Saudi Arabia 1 12 2020 1 12 2020 2020 10 209338 12 2019 7 10 2020 © The Author(s) 2020Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.Titanium alloy nanoparticle has a variety of applications in the manufacturing of soap and plastic, microsensors, aerospace design material, nano-wires, optical filters, implantation of surgical, and many biological treatments. Therefore, this research article discussed the influence of nonlinear radiation on magneto Williamson fluid involving titanium alloy particles through a thin needle. The arising system of partial differential equations is exercised by the similarity transformations to get the dimensional form of ordinary differential equations. The dual nature of solutions is obtained by implementing bvp4c. The study of stability has been carried out to check which of the results are physically applicable and stable. Influences of pertinent constraints on the flow field are discussed with the help of graphical representations and the method validation is shown in Table 1. The results imply that more than one result is established when the moving needle and the free-stream travel in the reverse directions. Moreover, the magnetic parameter accelerates the severance of boundary-layer flow, while the separation delays in the absence of the nanoparticle. The velocity gradient of nanofluid decays owing to the Williamson parameter in both branches of the outcome, while the temperature shrinks in the first or upper branch solution (stable one) and uplifts in the second or lower branch solution (unstable one). The size of the needle decreases the velocity in the upper solution and accelerates in the lower solution. The patterns of streamlines are more complicated due to the reverse direction of the free stream and thin needle. Subject terms EngineeringMathematics and computingissue-copyright-statement© The Author(s) 2020 ==== Body Introduction A thin moving needle structure is explained as a paraboloid revolution regarding its direction of axes in addition to the erratic thickness. Lee1 seemed to be very energetic and young first researcher who discussed the flow via the thin moving needle. He perceived that the thickness displacement and drag per-unit length moderate very gradually while increasingly thin needle, but ultimately turn into zero as the needle disappears. It is significant to note that thin moving needle movement disturbs the path of the free-stream in the flow. This observable fact is the key to determine the temperature and velocity profiles in the flow of experimental type research. The handling of the thin moving needle is a progressively more imperative feature in industries of engineering and medicine. For example, problems with blood flow, anemometer hot wire to calculate the velocity of the wind, coating and lubricating of wires and transportation. Narain and Uberoi2 expanded the work of Lee by taking the combination of free and forced convective flow with heat transfer from a thin moving needle. They presented the similarity solution as well as a series solution. Wang3 scrutinized the mixed convection flow from a heated tip through a thin needle. Their results indicated that the solution is unique for aiding flow; however, the dual solutions exist for the opposing flow. Multiple solutions of axisymmetric flow passing from a slim needle moving in the opposite or the same path to the free stream were obtained by the researcher Ishak et al.4. Ahmed et al.5 discussed the influences of assisting flow as well as opposing flow from a thin needle. The Keller-box is utilized to obtain the numerical solutions. They observed that the flow with the characteristics of heat transfer is considerably influenced by the size of the needle and mixed convective parameter. Entropy analysis of flow passing from a thin needle in a parallel stream with radiation influence was investigated by Afridi et al.6. Their results have shown that the entropy generation shrinks due to the needle size. Recently, many researchers utilized the concept of nanofluid regarding the geometry of the thin needle. For an instant, Hayat et al.7 scrutinized the erratic heat flux impact close to a stagnation-point through the phenomena of heat transfer applying water-based carbon nanoliquid. The results indicated that the velocity appreciably augments due to nanoparticle volume fraction. Soid et al.8 presented the stability analysis of MHD radiative flow containing nanofluid through a thin needle and obtained the multiple results using bvp4c solver and they observed that the presence of nanoparticle, the domains of the solution are smaller. The Tiwari–Das nanoliquid model with the characteristic of heat transfer through a thin needle with mixed convection has been developed to achieve the multiple outcomes by Salleh et al.9. They performed the stability analysis and found that the upper branch solution is stable, while the lower one is unstable. Salleh et al.10 scrutinized the heat source on flow comprising nanoliquid through a thin needle with chemical reaction. They explored that dual outcomes exist only when the thin needle travels against the flow direction. The exploration of a magnetic field containing menthol based Al-Cu hybrid nanoliquid through a moving needle submerged horizontally was investigated by Sulochana et al.11. They have shown that due to the size of the needle, the rate of heat transport augments. Tlili et al.12 examined the slip effect on 3D magneto flow comprising menthol based hybrid (Cuo–MgO) nanofluid across an erratic surface. They have examined that the combination of Cuo–MgO works as an excellent insulator. Kumar et al.13 examined the influence of thermal radiation on the time-dependent flow of Casson fluid over a curved exponential heated surface with erratic heat source/sink and Joule heating. The thermal transfer of the time-dependent thin-film flow containing Oldroyd-B ferroliquid suspended in water-based CoFe2O4 (cobalt ferrite) with a magnetic field was scrutinized by Tlili et al.14. They observed that the time-dependent parameter and Deborah number has the command to control the rate of heat transfer. Recently, Tlili et al.15 inspected the 3D magneto flow containing aluminum alloys amalgamated hybrid nanofluid through an irregular thickness sheet with a slip impact. It has been observed that the stimulus of the magnetic field is smaller on hybrid nanoliquid compared to nanoliquid. Properties of thermo-physical particularly the thermal conductivity of regular liquids like glycol, water, oil of the engine can be enhanced if the tiny sized particles (1–100 nm) are merged into the regular liquids. These merged liquids are named as nanoliquid pioneered by Choi et al.16 in line to present the nanoliquids in the regular liquids in the form of engineered colloidal suspensions. Nanoliquids are composed of distinct materials like metals, ceramics, alloys, nanotube, semi-conductors, and composite elements. In the last two decades, nanoliquid has been utilized as a superior liquid in heat transfer, particularly in the production of chemical, transportation, solar collector, cooling of electronics, power generation, and industries in biomedical17–19. The small size of the nanomaterials guides to improved stability of a suspension, the capability to flow efficiently without the obstruction of the system, and provided the improved physical as well as thermal properties. Also, it has been shown that nanoparticles existence escorts to an augmentation up to 15% to 40% in the regular fluid thermal conductivity. Haq et al.20 argued the combination of the thermal as well as velocity slip on MHD flow containing the nanoparticles from a stretched sheet through radiation effect. They observed that the Lorentz force declines the velocity and uplift the temperature, while the slip effect declines the heat transfer. Hayat et al.21 utilized the condition of mass flux to discuss the heat absorption/generation impact on flow comprising nanofluid from a stretching sheet with nonlinear radiation. They found that that heat absorption/generation uplifts the temperature as well as the heat transfer rate. Rehman et al.22 utilized the new condition of mass flux on flow with the heat transfer involving revised second-grade nanoliquid through a stretched nonlinear sheet. They have explored that the concentration of nanofluid became weaker due to the Lewis number. The impact of bio-convection with nanoparticles in second-grade liquid through a thin film with gyrotatic microorganisms in the presence of the passive control condition was discussed by Khan23. He observed that the velocity depreciated due to the Rayleigh number. Ahmed et al.24 scrutinized the nonlinear radiation effect on the time-dependent thin-film flow comprising the Maxwell nanoliquid through a rotating disk with activation energy. The outcomes exposed that the temperature of nanofluid significantly augments owing to the impact of thermophoresis. Wakif et al.25 reported thermo migration of nano/tiny-sized particles through the various fluids motion. Their outcomes have shown that the drag force decreases due to thermophoresis. Riaz et al.26 deliberated the impact of slip on the peristaltic flow containing nano-sized particles through a curved channel in a porous medium. They also added the creeping flow owing to the small Reynolds number. Babazadeh et al.27 reported the modeling of migration of magneto nanoparticle inside a porous space with radiation effect and shape factor. They observed that the greater convection mode can be gained due to the buoyancy force. Majeed et al.28 explored the characteristics of heat transfer comprising the magneto ferrofluid with dipole impact by controlling the thermal and momentum boundary-layer region. They concluded that the temperature of nanofluid decelerates due to nanoparticle. Tarakaramu and Satya Narayana29 studied the influence of heat source on the magneto time-dependent flow of Eyring-Powell nanofluid through a moving surface with chemical reaction and radiation effects. Recently, Tlili et al.30 premeditated the time-dependent free convection flow through the mixture of engine oil based TC4/NiCr nanofluid from a revolving cone in a porous medium. The transmission of the heat transfer rate of hybrid nanofluid is moderately lesser than the regular fluid. The scrutiny of non-Newtonian liquids in the current era has received immense interest due to their industrial and engineering intention. These liquids are indispensable, especially for the manufacturing of materials, production of a chemical and lubricants, mixtures of the polymer, etc. It is prominent that the investigation of non-Newtonian liquids and their properties contain numerous complexities due to the addition of rheological impacts in the constitutive leading equations. These supplementary physical and rheological quantities provide to augment the further complex and non-linear equations. Among different models, the Williamson liquid model is one of the signifying models. This model has an unambiguous benefit over other non-Newtonian models since it holds a minimum as well as maximum viscosities which give better outcomes for pseudoplastic liquids. Khan and Khan31 employed the HAM technique to get the series solution of Williamson liquid. The analysis of two-dimensional flow involving Williamson liquid from an exponential and linear stretching surface was deliberated by Nadeem et al.32,33. Hayat et al.34 explored MHD impact on time-dependent flow from a porous plate with fluid such as Williamson. They found that the velocity gradient decreases due to Williamson parameter. The effect of chemical reaction on the electrical conducting flow of Williamson nanoliquid from a stretched sheet immersed in the porous media with melting heat transfer was scrutinized by Krishnamurthy et al.35. They examined that the drag force and heat transfer are higher due to the presence of nanoparticle. The chemical reaction and radiation on Newtonian as well as non-Newtonian liquids near a stagnation-point from an incline cylindrical sheet in a dual stratified medium were evaluated by Rahman et al.36. They concluded that the heat transfer decreases due to thermophoresis and Brownian parameters. The time-dependent flow comprising Williamson fluid through a porous stretched sheet with MHD suspended nanoliquid was reviewed by Bibi et al.37. The radiation impact on the time-dependent MHD flow of Williamson fluid in the existence of nanoliquid through the radially stretched surface was scrutinized by Hashim et al.38. They scrutinized that the heat transfer declines with augmenting the magnetic field. They concluded that the velocity of nanofluid and temperature uplift due to Brownian motion. Sarojamma et al.39 inspected the homogenous-heterogeneous reaction along with erratic thermal conductivity involving non-Newtonian fluid through an uneven surface. They revealed that the velocity decreases due to the wall parameter. Tarakaramu and Satya Narayana40 investigated the Lorentz motion of a non-Newtonian fluid through a linear stretched surface with heat generation embedded in a porous medium. They found that heat generation and magnetic fields are responsible for greater heat transport in the fluid flow. Recently, Sandhya et al.41 explored the impacts of activation energy and second-order slip on magneto buoyancy flow through an exponential stretched surface with radiation effect. The usage of convection heat transfer along with non-Newtinian nanofluids is crucial to design several kinds of thermal equipments. Thus, the main aim of the current problem is to scrutinize the effect of MHD involving titanium alloy nanoparticle with Williamson fluid past a thin needle. Also, the impact of nonlinear radiation is invoked. Tiwari–Das model is taken to simulate the flow problem. The dual complex nature of solutions are obtained which most of the researchers missed out. Also, the stability analysis is also performed. Moreover, the accumulation of Williamson liquid containing magnetite titanium alloy nanoparticle composes a more complex mixture compared to the host nanoliquid. This explores an incorporate a novel epoch for researchers to determine the heat transfer containing the titanium alloy nanoliquid characteristics. In addition, Impacts of the significant constraints elucidate in detail with the assistance of plots and tabular form. Mathematical formulation As mentioned in Fig. 1, the two-dimensional incompressible flow of titanium alloy nanoparticle in the presence of Williamson liquid through a thin needle is scrutinized. The coordinates are utilized in cylindrical form x1,r1, where the coordinates r1- and x1- signifies radial and axial directions, respectively. Further, we incorporate the impacts of nonlinear radiation and magneto-hydrodynamics (MHD) on flow with characteristics of heat transfer. It is prominent that the needle is deemed as slim whilst the thickness of the needle doesn’t surpass that of the boundary layer flow (BLF) over it. Since the size of the needle is small, the pressure gradient is disregarded. Though, the impact of transverse-curvature is needed. The erratic magnetic field B(x1)=B0/x10.5 is inflicted in the path of the moving fluid while induced magnetic is unnoticed owing to the negligibly small Reynolds number. Also, it is perceived that T∞ and Tw signify the ambient and constant temperature of a thin needle with Tw>T∞. Moreover, the needle travels through steady velocity by uw on the contrary or identical direction to the velocity of free-stream which is u∞. The Cauchy stress-tensor for the Williamson fluid is suggeted by Dapra et al.42 as 1 S1=-pI+τ1, 2 τ1=μ∞+μ0-μ∞1-Γ1γ˙1A1, Figure 1 Geometry of the problem. Here S1 is the extra stress-tensor, μ0 and μ∞, respectively are limiting viscosities at zero and infinite shear rates, p the pressure, I the identity vector, Γ1,A1 represent the time relaxation and the first tensor of Rivlin-Erickson, respectively and γ˙1 is presented in the following form 3 γ˙1=12π, where π=traceA12. Here we discussed the case in which Γ1γ˙1<1 and μ∞. Therefor, τ1 can be simplified as 4 τ1=μ01-Γ1γ˙1A1, After employing the bionomial expansion, one gets 5 τ1=μ01+Γ1γ˙1A1. The expressions of governing flow with these postulations are4,6,38 6 ∂∂r1r1v1+∂∂x1r1u1=0 7 μnanofρnanof1r1∂∂r1r1∂u1∂r1+2Γ1∂u1∂r1∂2u1∂r12+Γ12r1∂u1∂r12-σnanofB2(x1)ρnanofu1=v1∂u1∂r1+u1∂u1∂x1 8 αnanof1r1∂∂r1r1∂T1∂r1+1ρcpnanof16σ∗T133k∗∂2T1∂r12+16σ∗3T123k∗∂T1∂r12=u1∂T1∂x1+v1∂T1∂r1, The substantial boundary-conditions are 9 T1=Tw,v1=0,u1=uwatr1=R(x1),T1→T∞,u1→u∞asr1→∞. Here, v1,u1 signify respectively the radial and axial velocity components, x1,r1 the coordinates in cylindrical form, ρnanof density of nanoliquid, μnanof viscosity of nanoliquid, ρcpnanof heat capacity of nanoliquid, σnanof electric conductivity of nanofluids, σ∗ Stefan Boltzmann constant, k∗ coefficient of mean absorption. The nanoliquid quantities are described as7–9 10 ρcpnanof=ϕρcps+1-ϕρcpf,μnanof=1-ϕ-2.5μf,αnanof=knanofρcpnanof-1,knanof/kf=ks+2kf+ϕkf-ks-2.5ks+2kf-2ϕkf-ks,ρnanof=ϕρs+1-ϕρf,σnanof/σf=1+σs/σf+2-σs/σf-1ϕ-2.53σs/σf-1ϕ. where ϕ called the volume fraction of nanoliquid, ρs,ρf signify the densities of nanoliquid and regular liquid, μf the viscosity of a regular liquid, ks,kf portrays the thermal conductivities of nanoliquid and regular liquid, respectively. We want to ease the computation of the current model owing to imply the following similarity variables as: 11 ψ1=νfx1Fη,θη=T1-T∞Tw-T∞,η=Ur12νfx1. Thus, in the aforementioned Eq. (11), the description of symbols is followed as νf the kinetic viscosity, U the merged velocity classified as U=u∞+uw, and ψ1 the stream function is described as v1=-1r1∂ψ1∂x1 and u1=1r1∂ψ1∂r1. Substituting η=c in the Eq. (11), which portrays the size and as well as the shape of the thin moving needle r1=R(x1) through its surface, which is specified via R(x1)=νfcx1U-0.50.5. Utilizing (11), the governing Eqs. (7)–(9) are converted as 12 F″′2η+8ηWeF″+2F″+6WeF″2+1-ϕ2.51-ϕ+ϕρsρfFF″-1-ϕ2.5σnanofσfMF′=0 13 knanofkfηθ″+θ′+Pr21-ϕ+ϕρcpsρcpfFθ′+43Rdθθr-1+12θθr-1+1ηθ″+θ′2+3θ′2θr-1η=0 14 θ(c)=1,F(c)=0.5λc,F′(c)=0.5λ,θ(η)→0,F′η→0.51-λasη→∞. Here, the dimensionless variables are Prandtl number, magnetic, Weissenberg, local Reynolds numbers, radiation, velocity ratio and temperature ratio parameters, respectively. These parameters are described as 15 Pr=νfαf,M=σfB022Uρf,We=2Γ1U2r1Rex1νfx1Rer1,Rex1=Ux1νf,Rer1=Ur1νf,Rd=kfk∗4T∞3σ∗,λ=uwU,θr=TwT∞. It is important to mention that the domain of the velocity ratio parameter is set to 0<λ<1 which specifies the fluid and needle progress in an identical direction. If λ>1, then the needle is going near a positive x1- direction, whereas the free-stream shifts in the way of a negative x1- direction. Here, we consider only the case λ≤1. Mathematically, the temperature gradient and velocity gradient (heat transfer rate and the surface drag force) are identified as 16 Nux1=x1qwkf(Tw-T∞),Cf=τwU2ρf, where qw and τw are recognized as the heat-flux and the shear stress respectively, which are specified as 17 qw=-knanof∂T1∂r1-16σ∗T133k∗∂T1∂r1r1=c,τw=μnanof∂u1∂r1r1=c, Utilizing (11), we have 18 Rex1-1/2Nux1=-2cθ′(c)knanofkf+43Rdθr-1θc+13,Rex11/2Cf=μnanofμf4cF″(c)(1-ϕ)2.5. Stability analysis of the solution The results numerically confirm that there occurs more than one solution called multiple (the lower branch (LB) and the upper branch (UB)) solutions and this idea is scrutinized in the Weidman et al.43, where they were exercising the stability analysis and established that which result is unstable and which one is stable. The key point to be noticed of this analysis is presented to obtain that which result provided the best path to the fluid flow mean physically realizable or stable result. For this purpose, more than one solution for various models has been investigated by Sharma et al.44, Yasin et al.45, Mansur et al.46, Roşca et al.47. We check the analysis through taking all the governing equations (from(7, 8) in the unsteady form except the continuity equation it will remain unchanged as follow: 19 ρnanof∂u1∂t+u1∂u1∂x1+v1∂u1∂r1=μnanof1r1∂∂r1r1∂u1∂r1+2Γ1∂u1∂r1∂2u1∂r12+Γ12r1∂u1∂r12-σnanofB2(x1)u1, 20 ρcpnanof∂T1∂t+u1∂T1∂x1+v1∂T1∂r1-αnanof1r1∂∂r1r1∂T1∂r1=16σ∗T133k∗∂2T1∂r12+16σ∗3T123k∗∂T1∂r12, The new dimensionless time variable is defined in the form τ=2Ut/x1 , where t stands for time. The new similarity transformations were followed as: 21 ψ1=νfx1Fη,τ,θη,τ=T1-T∞Tw-T∞,η=Ur12νfx1,τ=2Utx1. By exercising the similarity transformations (21) into Eqs. (19) and (20), we get 22 2∂2F∂η2+2η∂3F∂η3+We6∂2F∂η22+8η∂2F∂η2∂3F∂η3-1-ϕ2.5σnanofσfM∂F∂η+1-ϕ+ϕρsρf1-ϕ2.5τ∂2F∂η∂τ∂F∂η+1-ϕ2.5F∂2F∂η2-1-ϕ2.5∂2F∂η∂τ=0, 23 knanofkfη∂2θ∂η2+∂θ∂η+43Rdθθr-1+12θθr-1+1η∂2θ∂η2+12∂θ∂η+3∂θ∂η2θr-1η+Pr21-ϕ+ϕρcpsρcpfF∂θ∂η-∂θ∂τ+τ∂F∂η∂θ∂τ=0,s, along the corresponding condition are 24 θc,τ-1=0,Fc,τ-0.5λc=0,∂Fc,τ∂η-0.5λ=0,θ∞,τ→0,∂F∞,τ∂η=0.51-λ. Then, we consider that27,32 25 Fη,τ-e-ξτHη=f0(η)θ(η,τ)-e-ξτGη=θ0(η). To investigate the temporal stability to specify the answers F=f0(η) and θ=θ0(η) which joined the two-point problem (22) and (23). However, ξ stand for the unknown eigenvalue parameter and as well as the functions Gη and Hη are fairly small as compared to θ0(η) and f0(η). Plugging Eq. (25) into Eqs. (22) and (23) yields the following problem of linearized eigenvalue. 26 2H0″+2H0″′+4We2ηf0″H0″+3f0″H0″′+2ηf0″′H0″-σnanofσf1-ϕ2.5MH0′+1-ϕ+ϕρsρf1-ϕ2.5f0″H0+1-ϕ2.5f0H0″-1-ϕ2.5ξH0′=0, 27 knanofkfG0′+ηG0″+Pr21-ϕ+ϕρcpsρcpff0G0′+H0θ0′+ςG0+2θr-1θ0+1G0θ0′2+2θr-1θ0+12θ0′G0′+123Rdθr-1G0θr-1θ0+12ηθ0″+θ0′2+43Rdθr-1θ0+13ηG0″+G0′2=0, and the subjected conditions (18) are follow as 28 H0′c=0,H0c=0,G0c=0,H0′∞→0,G0∞→0. Working out the Eqs. (26)-(28) numerically, we acquire the infinite eigenvalues ξ1<ξ2<ξ3<ξ4…. The flow is said to be physically realizable (stable) if ξ>0 while it is not physically realizable (unstable) for the choice of ξ<0. The range for the eigenvalue is specified by Harris et al.48 where they have relaxing conditions on G0η and H0η. In this recent work, we take H0′η→0 as η→∞ and then tackled the Eqs. (26) and (27) together with (27) and along with once more new condition which is H0′0=1. Physical explanation The transmuted nonlinear ODEs (12) and (13) along with boundary restriction (14) are tackled numerically with the help of bvp4c solver via the formula of 3-stage Lobatto IIIa. The characteristics and features of different pertinent variables appeared in the problem on the profiles of liquid velocity, temperature, Nusselt number, and the skin friction are portrayed in Figs. 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14 and 15. The values appeared in the problem during the computation are considered fixed as θr=Rd=c=0.1, We=0.05, M=0.1, λ=-3.5, ϕ=0.1. Table 1 depicts the significance of the physical and thermal properties of titanium alloy nanoparticle. A comparison of our results with available outcomes is made in Table 2 to validate our problem and discovered in a tremendous closeness. Also, Table 3 is constructed for comparision between two techniques and they are in an excellent agreement. Besides, in all figures, solid lines explain the first (UBS) solution, whereas the dotted lines display the second (LBS) solution.Figure 2 Influence of ϕ on Rex11/2Cf. Figure 3 Influence of ϕ on Rex1-1/2Nux1. Figure 4 Influence of Rd on Rex1-1/2Nux1. Figure 5 Influence of θr on Rex1-1/2Nux1. Figure 6 Influence of We on Rex11/2Cf. Figure 7 Influence of We on Rex1-1/2Nux1. Figure 8 Influence of M on Rex11/2Cf. Figure 9 Influence of M on Rex1-1/2Nux1. Figure 10 Influence of We on F′η. Figure 11 Influence of We on θη. Figure 12 Influence of c on F′η. Figure 13 Influence of M on F′η. Figure 14 Influence of M on θη. Figure 15 Streamlines pattern for nanofluid. Table 1 Thermo physical properties of base fluid and Ti6Al4V (Makinde et al.49). Material Water Ti6Al4V CpJ/kgK 4179 0.56 ρkg/m3 997.1 4420 kW/mK 0.613 7.2 σS/m 0.005 5.8×106 Pr 6.2 – Table 2 Assessment of the values of F″c when ϕ=M=We=0. c λ=0 λ=-1 λ=-1 Soid et al.8 Present Soid et al.8 Present Present First solution Second solution First solution Second solution 0.01 8.491454 8.4915 26.599394 2.805533 26.6021 2.8031 0.1 1.288778 1.2888 3.703713 0.389103 3.7162 0.3884 0.2 0.751665 0.7515 2.005424 0.227837 2.0055 0.2278 Table 3 Assessment of the values of F″c for different values of ϕ when ϕ=0.1,We=0.05, λ=-3.5,c=0.1. ϕ bvp4c Keller-box First solution Second solution First solution Second solution 0 1.2173 0.6990 1.2196 0.6999 0.03 1.3028 0.7688 1.3052 0.7694 0.06 1.4028 0.8399 1.4066 0.8401 0.1 1.5635 1.9361 1.5651 1.9381 Deviation of the skin coefficient and the local Nusselt number Figures 2 and 3 present the deviation of the skin friction Rex11/2Cf and the Nusselt Rex1-1/2Nux1 for varying nanoparticle volume fraction ϕ versus ratio velocity parameter λ. It is scrutinized from Fig. 2 that the increment in the value of ϕ as a result the UBS, as well as LBS, augment. This is because that the nanoparticles thermal conductivity becomes larger which causes more suspension of particles in the host fluid. Figure 3 reveals that the values of Rex1-1/2Nux1 confirming and to show the escalating behavior in the LBS, as well as the UBS owing to the large value of ϕ. However, the multiple results are existing in the regions of λc<λ<0 while a single result of the phenomenon of λ=λc<0 and in the case λ>0 (the free-stream as well as the needle travel in the identical direction). It is observed from our calculation that the critical points λc are -3.8448, -3.7700, and -3.7670 for ϕ=0,0.05,1, respectively. It is concluded from critical values that the separation of boundary delays in the absence of the nanoparticle ϕ=0. Figures 4 and 5 discussed the behavior of Weissenberg number We on Rex11/2Cf and Rex1-1/2Nux1, respectively. These figures suggest the decreasing tendency in the LBS as well as UBS for augmenting the value of We. Here, the critical values in the case of the Weissenberg effect are − 3.8760, − 3.8650 and − 3.8600 for We=0,0.03,0.05, respectively. Impacts of radiation, and heating parameters on Rex11/2Cf and Rex1-1/2Nux1 are portrayed in Figs. 6 and 7, respectively. Figure 6 depicts that the outputs of the Nusselt number enhance in the LBS, and UBS with the augmentation in the radiation parameter. Due to the influence of radiation, the critical values are obtained as − 3.8610, − 3.8605, and − 3.8605 for the three selected values of Rd=0.1,0.3,0.5, respectively. On the other hand, the impact of the heating parameter in Fig. 7 is quite different (like Fig. 6). More precisely, the values of the Nusselt parameter show a decelerated behavior in the UBS and escalating in the LBS owing to upsurges the value of the heating parameter. In this case, the established critical values are the following − 3.8610, − 3.8600, and − 3.8600 for θr=0.1,0.3,0.5, respectively. The trajectories of Rex11/2Cf and Rex1-1/2Nux1 for selected values of M are demonstrated in Figs. 8 and 9. The graph in Fig. 8 confirms that the values of Rex11/2Cf behaving increasingly in the UBS and contrary in the LBS. However, the contrasting behavior is scrutinized for Rex1-1/2Nux1 as depicted in Fig. 9. The critical values are − 4.2760, − 4.0240, and − 3.8610 for M=0,0.03,0.05, respectively. Thus, the critical values λc decrease, suggesting that the magnetic parameter rises the separation of the boundary. Deviation of the velocity and temperature fields The impact of the Weissenberg parameter We on fluid velocity and temperature distribution is shown in Figs. 10 and 11, respectively. Figure 10 displays that the fluid velocity shrinks with augmenting We in the stable and unstable branch solutions which in turn enhances the thickness of the velocity BLF. Since the Williamson parameter is the ratio of relaxation to retardation time. Therefore, uplifting values of We upsurges the relaxation time. Due to this, the liquid particles need extra time to reinstate their original path. However, the temperature distribution decline with escalating values of We in the upper branch solution, whereas the opposite behavior is seen in the lower branch solution as depicted in Fig. 11. The influence of needle size c on the fluid velocity is depicted in Fig. 12. It is visible that growing the size of the needle disturbs the free-stream, which consecutively causes to moderate the velocity of the fluid in the UBS, while the liquid velocity increases in the LBS. Thus, the magnitude of the momentum boundary-layer increases in the upper solution and shrinks in the lower solution. Thus, the boundary thickness can be organized via the needle size. In Figs. 13 and 14 show the behavior of the magnetic field M on the velocity and temperature profiles, respectively. Figure 13 depicts that the liquid velocity shrinks with increasing M in the UBS and consequently enhances the field of velocity gradient and the magnitude of the boundary layer. On the other hand, the contrary behavior is noticed in the LBS. In physical point of view, the accumulation of a magnetic field will convey a kind of resistive force known as the Lorenz or drag force, which has a tremendous property to reduce or stop the liquid velocity. Figure 14 establishes that the temperature profile upsurges in the UBS and diminutions in the LBS for the larger values of M. Physically, a Lorentz force is created by the applied magnetic field which has an important role to delay the motion of liquid, and consequently, the thermal boundary layer and temperature enhances. Thus, needle temperature can be managed by managing the magnetic field strength. Figure 15 shows the patterns of streamlines for nanofluid. The patterns illustrate that the streamlines are more obscured and divided into two regions owing to due to reverse the direction of free stream and thin needle. Main remarks This paper scrutinized the impact of nonlinear radiation on MHD flow with the characteristic of heat transfer suspended in Williamson fluid comprising titanium alloy nanoparticle through a thin needle. The system of arising PDEs is altered into the highly nonlinear ODEs through similarity variables and then functioned out numerically via the built-in function in Matlab called bvp4c. The dual nature of solutions is found for the current problem where the results are shown in both tabular and graphical form (sees physical explanation). The influences of sundry constraints on the flow field are elaborated with the assistance of plots. Also, the stability analysis is highlighted in the presented work to check which solutions are physically stable or not. The vital outcomes of the paper are summarized as the following:The rate of heat transfer as well as the skin friction augments due to nanoparticle volume fraction in both UB and LB solutions. The liquid velocity gradient is a decaying function of the Weissenberg number in both solutions, while the liquid temperature declines in UBS and upsurges in LBS. The rate of heat transfer increases due to the radiation parameter in both branches of solutions. Multiple solutions are achieved to survive when free stream and needle shows their progress in conflicting directions whilst the result is unique when they move in the same direction. The heat transfer rate declines in the UBS and augments in the LBS as we upsurge the impact of heating parameter. The skin friction is decelerating in the UBS and as well as in the LBS owing to the implementation of the larger values of the Weissenberg number while the same behavior is perceived for the rate of heat transfer. The skin friction augments due to the magnetic field in UBS and diminishes in the LBS, where the reverse impact is seen in the heat transfer rate. List of symbols A1First tensor of Rivlin–Erickson BMagnetic field (kgs−2 A−1) CfSkin friction coefficient Nux1Nusselt number B0Strength of the magnetic field (A/m) cSize of the needle MMagnetic parameter PrPrandtl number IIdentity vector k∗Coefficient of mean absorption kThermal conductivity (Wm−1 K−1) qwHeat-flux (Wm−2) FSimilarity fuction for velocity Rex1,Rer1Local Reynolds number RdRadiation parameter S1Extra stress-tensor T1Temperature (K) TwTemperature of the wall (K) T∞Ambient temperature (K) uwVelocity (m/s) u∞Free-stream velocity (m/s) u1,v1Velocity components (m/s) WeWeissenberg number r1,x1Radial and axial directions (m) Greek symbols λVelocity ratio parameter νfKinematic viscosity of the base fluid (m2/s) αfThermal diffusivity (m2/s) μ0,μ∞Limiting viscosities at zero and infinite shear stresses (Nsm−2) γ˙1Shear stress θDimensionless temperature θrTemperature ratio parameter ϕVolume fraction of nanoliquid σElectric conductivity of nanofluids (Sm−1) σ∗Stefan Boltzmann constant (Wm−2 K−4) ρNanofluid density (kg m−3) μViscosity of nanoliquid (Kg m−1 s−1) ρcpHeat capacity of nanoliquid (J/K) Γ1Time relaxation (s) ψ1Stream function ηSimilarity variable τwWall shear stress Subscript nanofNanofluid fBase fluid Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Acknowledgements This project was supported by the Deanship of Scientific Research at Prince Sattam Bin Abdulaziz University under the research project No. 2020/01/16436. Author contributions U.K. formulated the problem. A.Z solved the problem I.K discussed the results K.S.N revised the manuscript. Funding This project was supported by the Deanship of Scientific Research at Prince Sattam Bin Abdulaziz University under the research Project No. 2020/01/16436. Competing interests The authors declare no competing interest. ==== Refs References 1. Lee LL Boundary layer over a thin needle Phys. Fluids 1967 10 1820 1822 2. Narain JP Uberoi SM Combined forced and free-convection heat transfer from vertical thin needles in a uniform stream Phys. Fluids 1973 15 1879 1882 10.1063/1.1693798 3. Wang CY Mixed convection on a vertical needle with heated tip Phys. Fluids 1990 2 622 625 10.1063/1.857709 4. Ishak A Nazar R Pop I Boundary layer flow over a continuously moving thin needle in a parallel free stream Chin. Phys. Lett 2007 24 2895 2897 10.1088/0256-307X/24/10/051 5. Ahmad S Arifin NM Nazar R Pop I Mixed convection boundary layer flow along vertical thin needles: Assisting and opposing flows Int. Commun. Heat Mass Transfer 2008 35 157 162 10.1016/j.icheatmasstransfer.2007.07.005 6. Afridi MI Qasim M Entropy generation and heat transfer in boundary layer flow over a thin needle moving in a parallel stream in the presence of nonlinear Rosseland radiation Int. J. Therm. Sci 2018 123 117 128 10.1016/j.ijthermalsci.2017.09.014 7. Hayat T Khan MI Farooq M Yasmeen T Alsaedi A Water-carbon nanofluid flow with variable heat flux by a thin needle J. Mol. Liq 2016 224 786 791 10.1016/j.molliq.2016.10.069 8. Soid SK Ishak A Pop I Boundary layer flow past a continuously moving thin needle in a nanofluid Appl. Therm. Eng 2017 114 58 64 10.1016/j.applthermaleng.2016.11.165 9. Salleh SNA Bachok N Arifin NM Ali FM Pop I Stability analysis of mixed convection flow towards a moving thin needle in nanofluid Appl. Sci. 2018 8 842 10.3390/app8060842 10. Afridi MI Qasim M Wakif A Hussanain A Numerical analysis of boundary layer flow adjacent to a thin needle in nanofluid with the presence of heat source and chemical reaction Symmetry 2019 11 543 10.3390/sym11050663 11. Sulochana C Aparna SR Sandeep N Impact of linear/nonlinear radiation on incessantly moving thin needle in MHD quiescent Al-Cu/methanol hybrid nanofluid Int. J. Ambient Energy 2020 10.1080/01430750.2020.1768895 12. Tlili I Nabwey HA Samrat SP Sandeep N 3D MHD nonlinear radiative flow of Cuo–MgO/methanol hybrid nanofluid beyond an irregular dimension surface with slip effect Sci. Rep. 2020 10 9181 10.1038/s41598-020-66102-w 32514160 13. Kumar KA Sugunamma V Sandeep N Effect of thermal radiation on MHD Casson fluid flow over an exponentially stretching curved sheet J. Therm. Anal. Calorim 2020 140 2377 2385 10.1007/s10973-019-08977-0 14. Tlili I Samrat SP Sandeep N Nabwey HA Effect of nanoparticle shape on unsteady liquid film flow of MHD Oldroyd-B ferrofluid Ain Shams Eng. J. 2020 10.1016/j.asej.2020.06.007 15. Tlili I Nabwey HA Ashwinkumar GP Sandeep N 3-D magnetohydrodynamic AA7072-AA7075/methanol hybrid nanofluid flow above an uneven thickness surface with slip effect Sci. Rep. 2020 10 4265 10.1038/s41598-020-61215-8 32144369 16. Choi S Eastman JA Enhancing thermal conductivity of fluids with nanoparticles ASME FED-231/MD 1995 66 99 105 17. Wong KV Leon OD Applications of nanofluids: current and future Adv. Mech. Eng 2010 2 519659 10.1155/2010/519659 18. Saidur R Leong KY Mohammad HA A review on applications and challenges of nanofluids Renew. Sustain. Energy Rev 2011 15 1646 1668 10.1016/j.rser.2010.11.035 19. Huminic G Huminic A Application of nanofluids in heat exchangers: a review Renew. Sustain. Energy Rev 2012 16 5625 5638 10.1016/j.rser.2012.05.023 20. Haq RU Nadeem S Khan ZH Akbar NS Thermal radiation and slip effects on MHD stagnation point flow of nanofluid over a stretching sheet Physica E 2015 65 17 23 10.1016/j.physe.2014.07.013 21. Hayat T Qayyum S Alsaedi A Shafiq A Inclined magnetic field and heat source/sink aspects in flow of nanofluid with nonlinear thermal radiation Int. J. Heat Mass Transfer 2016 103 99 107 10.1016/j.ijheatmasstransfer.2016.06.055 22. Rehman, M. Ur., Khan, M. & Manzur, M. Boundary layer flow and heat transfer of a modified second grade nanofluid with new mass flux condition. Res. Phys 10, 594–600 (2018). 23. Khan NS Bioconvection in second grade nanofluid flow containing nanoparticles and gyrotactic microorganisms Braz. J. Phys 2018 48 3 227 241 10.1007/s13538-018-0567-7 24. Ahmed J Khan M Ahmad L Transient thin-film spin-coating flow of chemically reactive and radiative Maxwell nanofluid over a rotating disk App. Phys. A 2019 125 161 10.1007/s00339-019-2424-0 25. Wakif A Animasaun IL Satya Narayana PV Sarojamma G Meta-analysis on thermo-migration of tiny/nano-sized particles in the motion of various fluids Chin. J. Phys. 2019 10.1016/j.cjph.2019.12.002 26. Riaz A Thermal analysis of peristaltic flow of nanosized particles within a curved channel with second-order partial slip and porous medium J. Therm. Anal. Calorim. 2020 10.1007/s10973-020-09454-9 27. Babazadeh H Numerical modelling for nanoparticle thermal migration with effects of shape of particles and magnetic field inside a porous enclosure Iran. J. Sci. Technol. Trans. Mech. Eng. 2020 10.1007/s40997-020-00354-9 28. Majeed A Zeeshan A Bhatti MM Ellahi R Heat transfer in magnetite (Fe3 O4 ) nanoparticles suspended in convectional fluids: Refrigrant-134A (C2 H2 F4 ), kerosene (C10 H22 ), and water (H2 O) under the impact of dipole Heat Transf. Res. 2020 51 217 232 10.1615/HeatTransRes.2019029919 29. Tarakaramu, N. & Satya Narayana, P.V. Radiation and chemical reaction effects on unsteady Eyring–Powell nanofluid flow over a moving surface, in Advances Fluid Dynamics, 393–405 (2020). 30. Tlili I Sandeep NG Reddy M Nabwey HA Effect of radiation on engine oil-TC4/NiCr mixture nanofluid flow over a revolving cone in mutable permeable medium Ain. Shams Eng. J. 2020 10.1016/j.asej.2020.03.011 31. Alam Khan N Khan H A boundary layer flows of non-Newtonian Williamson fluid Nonlinear Eng 2014 3 2 107 115 32. Nadeem S Hussain ST Lee C Flow of a Williamson fluid over a stretching sheet Braz. J. Chem. Eng 2013 30 3 619 625 10.1590/S0104-66322013000300019 33. Nadeem S Hussain ST Heat transfer analysis of Williamson fluid over exponentially stretching surface Appl. Math. Mech. Engl. Ed 2014 35 4 489 502 10.1007/s10483-014-1807-6 34. Hayat T Khalid U Qasim M Steady flow of a Williamson fluid past a porous plate Asia Pac. J. Chem. Eng 2012 7 302 306 10.1002/apj.496 35. Krishnamurthy MR Prasannakumara BC Gireesha BJ Gorla RSR Effect of chemical reaction on MHD boundary layer flow and melting heat transfer of Williamson nanofluid in porous medium Eng. Sci. Tech. Int. J 2016 19 1 53 61 36. Rehman KU Malik MY Makinde OD Malik AA A comparative study of nanofluids flow yields by an inclined cylindrical surface in a double stratified medium Eur. Phy. J. Plus 2017 132 427 10.1140/epjp/i2017-11679-1 37. Bibi M Rehman KU Malik MY Tahir M Numerical study of unsteady Williamson fluid flow and heat transfer in the presence of MHD through a permeable stretching surface Eur. Phys. J. Plus 2018 133 154 10.1140/epjp/i2018-11991-2 38. Hashim KM Hamid A Convective heat transfer during the flow of Williamson nanofluid with thermal radiation and magnetic effects Eur. Phys. J. Plus 2019 134 50 10.1140/epjp/i2019-12473-9 39. Sarojamma G Sreelakshmi K Jyothi PK Satya Narayana PV Influence of homogeneous and heterogeneous chemical reactions and variable thermal conductivity on the MHD Maxwell fluid flow due to a surface of variable thickness Defect Diffus. 2020 401 148 10.4028/www.scientific.net/DDF.401.148 40. Tarakaramu, N. & Satya Narayana, P.V. Influence of heat generation/absorption on 3D magnetohydrodynamic Casson fluid flow over a porous stretching surface, in Advances Fluid Dynamics, 381–392 (2020). 41. Sandhya G Sarojamma G Satya Narayana PV Venkateswarlu B Buoyancy forces and activation energy on the MHD radiative flow over an exponential stretching sheet with second order slip Heat Transf. 2020 10.1002/htj.21904 42. Dapra I Scarpi G Perturbation solution for pulsatile flow of a non-Newtonian Williamson fluid in a rock fracture Int. J. Rock Mech. Min. Sci. 2007 44 271 278 10.1016/j.ijrmms.2006.07.003 43. Weidman PD Kubitschek DG Davis AMJ The effect of transpiration on self-similar boundary layer flow over moving surfaces Int. J. Eng. Sci. 2006 44 11–12 730 737 10.1016/j.ijengsci.2006.04.005 44. Sharma R Ishak A Pop I Stability analysis of magnetohydrodynamic stagnation-point flow toward a stretching/shrinking sheet Comp. Fluids 2014 102 94 98 10.1016/j.compfluid.2014.06.022 45. Yasin MHM Ishak A Pop I Boundary layer flow and heat transfer past a permeable shrinking surface embedded in a porous medium with a second-order slip: a stability analysis Appl. Thermal Eng. 2017 115 1407 1411 10.1016/j.applthermaleng.2016.08.080 46. Mansur S Ishak A Pop I The magnetohydrodynamic stagnation point flow of a nanofluid over a stretching/shrinking sheet with suction PLoS ONE 2015 10 3 e0117733 10.1371/journal.pone.0117733 25760733 47. Roşca AV Pop I Flow and heat transfer over a vertical permeable stretching/shrinking sheet with a second order slip Int. J. Heat Mass Transf. 2013 60 355 364 10.1016/j.ijheatmasstransfer.2012.12.028 48. Harris SD Ingham DB Pop I Mixed convection boundary-layer flow near the stagnation point on a vertical surface in a porous medium: Brinkman model with slip Tran. Porous Media 2009 77 2 267 285 10.1007/s11242-008-9309-6 49. Makinde OD Mahanthesh B Gireesha BJ Shashikumar NS Monaledi RL Tshehla MS MHD nanofluid flow past a rotating disk with thermal radiation in the presence of aluminum and titanium alloy nanoparticles Defect Diffus. Forum 2018 384 69 79 10.4028/www.scientific.net/DDF.384.69