I am solving a hybrid nanofluid flow problem in a bifurcated artery using Maple. The governing equations for velocity and temperature are solved using dsolve(..., numeric, method=bvp[midrich]).

My Maple code successfully produces for both the artery  parentartery_and_daughter_artery_error.mw.

The velocity profiles are obtained correctly using odeplot.

However, I want to compute additional physical quantities and generate plots similar to the velocity profiles.

Specifically I want to plot:

1. Flow rate Q versus axial distance z

2. Impedance (flow resistance) λ versus z

3. Wall shear stress τ versus z

for different values of Hartmann number Ha.

The formulas I am using are

Flow rate:

Q=2π(R∫01ηw(η) dη+R2∫01w(η) dη)Q = 2\pi \left( R \int_0^1 \eta w(\eta)\,d\eta + R_2 \int_0^1 w(\eta)\,d\eta \right)Q=2π(R∫01​ηw(η)dη+R2​∫01​w(η)dη)

Wall shear stress:

τ=μ∣dwdr∣\tau = \mu \left|\frac{dw}{dr}\right|τ=μ​drdw​​

Impedance:

λ=∣dp/dz∣Q\lambda = \frac{|dp/dz|}{Q}λ=Q∣dp/dz∣​
Please help me to solve this question.

Hello everyone,

I am trying to reproduce some results from the paper:

Peristaltic flow of a magnetohydrodynamic nanofluid through a bifurcated channel

I want to generate:

• Fig 1 (channel geometry)

• Fig 8 (pressure gradient vs axial coordinate ξ)

• Fig 11 (streamlines)

• Table 1 (resistance values)

I am using Maple 2018.

Below is my Maple code:
p1_ravi_sir.mw

Problem:

The plots do not match the paper

p1_ravi_sir.pdf

Why does my resistance value not change with M

 Could someone please help me:

• Correct the Maple code

• Solve equations (18–22) properly in Maple

• Generate the correct plots 

Thank you very much.

Variation of (a) Skin friction ∂W/∂Z​, (b) Heat Transfer ∂θ/∂Z​, and (c) Mass Transfer ∂ϕ/∂Z​ for γ=10.0, Pr=7.0, ε=1.0, Nt=0.4, and Nb=0.2.

how to get this values by solving the PDE by using the pdsolve method 

Variation of (a)W(b) θ and (c) ϕ for different value of γ when ε = 10.0, Nb = 0.4, ε = 10.0, Sc = 0.5 and Pr = 7.0.
consider X as 0.1

plume_work.mw

How to get a 2D plot for the following command: 
error_in_2d_plot_cf_and_nu.mw

I am trying to plot a contour graph for my problem for (psi) function in the particular boundary, and even though it's working, but the contour  plot is not appearing at the end. Could anyone help me with the code to get proper graph in the specified boundary. 

i have ploted the graph in python i got a plot similar to that i am trying maple but i am not able to plot it. could any one help me to solve.

contour_plots_error_in_wavey_flow.mw