@KIRAN SAJJAN   @dharr

How to get the plots for this govering equation. as suggested by kiran sajjan

Mn := .5; We := 1.5; Omega := 1.5; Grt := 1.5; Grc := 1.5; Grf := 1.5; Pr := 6; Nb := 1.5; Nt := 2; Ntc := .5; Nct := .8; beta := 1.5; d := .8; x := 0; varphi := (1/2)*Pi; b := .6; a := .5; t := 1; lambda := 1; Rd := 3; alpha := 0.5e-1; f := 3; yL := -.8; yR := 1.5; OdeSys := {diff(Gamma(y), y$2)+Nt*(diff(theta(y), y$2))/Nb = 0, diff(phi(y), y$2)+Nct*(diff(theta(y), y$2)) = 0, (Pr*Rd+1)*(diff(theta(y), y$2))+Nb*Pr*(diff(theta(y), y))*(diff(Gamma(y), y))+Nt*Pr*(diff(theta(y), y))^2+Ntc*Pr*(diff(phi(y), y$2))+beta = 0, diff(psi(y), y$4)-Mn^2*(diff(psi(y), y$2))/(Omega^2+1)+Grt*theta(y)+Grc*phi(y)-Grf*Gamma(y)-2*We^2*(Mn^2*(diff(psi(y), y$2))/(Omega^2+1)-Grt*theta(y)-Grc*phi(y)+Grf*Gamma(y))^3 = 0}; Cond := {(D(psi))(1+a*cos(2*Pi*(x-t)/lambda))-alpha*((D@@2)(psi))(1+a*cos(2*Pi*(x-t)/lambda)) = -1, (D(psi))(-d-b*cos(2*Pi*(x-t)/lambda+varphi))+alpha*((D@@2)(psi))(-d-b*cos(2*Pi*(x-t)/lambda+varphi)) = -1, Gamma(1+a*cos(2*Pi*(x-t)/lambda)) = 0, Gamma(-d-b*cos(2*Pi*(x-t)/lambda+varphi)) = 1, phi(1+a*cos(2*Pi*(x-t)/lambda)) = 0, phi(-d-b*cos(2*Pi*(x-t)/lambda+varphi)) = 1, psi(1+a*cos(2*Pi*(x-t)/lambda)) = (1/2)*f, psi(-d-b*cos(2*Pi*(x-t)/lambda+varphi)) = -(1/2)*f, theta(1+a*cos(2*Pi*(x-t)/lambda)) = 0, theta(-d-b*cos(2*Pi*(x-t)/lambda+varphi)) = 1};
sol := dsolve({Cond[], OdeSys[]}, numeric, output = listprocedure, range = yL .. yR);
i  want same plots given in this manner

expected results

@KIRAN SAJJAN 

You are given correct post .

If I am implement to my plots means I am getting circular path only. 

I have one dought psi versus y  2d plots getting correct wave / ossilation form but why I m getting circular path for contour of same psi function

What is the mistake in plot command.getting circular pattern. Where I need to change in this worksheet to get same pattern in stream lines / contour plot

pulatile_flow_error.mw

@acer 

Thank you for your reply sir, 

Dear sir here i have given a work  will all equations they have ploted streamlines like a flow pattern of fluid.

i also need to find similar plots which is given in plots 9, 10 or 11 given in pdf peristatic_flow_paper_demo.pdf

but i have similar equations but in my equations i am not able to get same streamline flow i am getting as circular path. 

it should be osulatary motion.

why it is circular motion in the above plots . in the base work wich is given as oscilatory. if iam using sine or cos function also it is giving the same circular motion only. 

if we mention plots in the circualr motion there will be wrong representation of stream function/ stram lines.

@dharr 

Dear sir howto draw stream lines for the above equation 

Using the stream function u=∂ψ/∂y,v=−∂ψ/∂x.

with(plots):

psi_fun := (x,y) -> eval(psi(y), Ans[1]) + 0.1*sin(2*Pi*x/200)*y;

contourplot(psi_fun(x,y),

            x=0..200, y=0..400,

            contours=20,

            filled=true,

            coloring=[blue, cyan, lightblue],

            axes=boxed,

            label

s=["x","y"]);

Similar plot I want

@sand15 paper2_new_efficiency_plots_2025.mw

for these plots

 i need the animation of that plots with the variation of  X=0..1 in first plot and tau=0..2 in second plot.

please give me solution for it

@Rouben Rostamian  

Download first_order_pde.mw

@Carl Love

Thank you

But what are the changes need to do in this plot sir.

I m not that much familier in this. 

by using this post Pde i have given my equation please check the mistakes and reply

I m waiting for the reply.

Daily i m checking the same post. Please do help for this.

Thank you

please check my question and rectify my errors.

i have posted a new question that is automatically deleted

Download circular_cylider_entropy_plot.mw

Similar problem asked few years back there answer is not there question is incomplete check here Related pot.

Modify my work sheet for solution