Advanced Search

pooyan1990

35 Reputation

3 Badges

5 years, 128 days
Contact pooyan1990

MaplePrimes Activity


These are replies submitted by pooyan1990

@Carl Love  Thanks. it worked....

pooyan1990 35
April 05 2023
0 0

@Carl Love 

Thanks. it worked.

@tomleslie  Thank you so much. Thi...

pooyan1990 35
November 09 2019
0 0

@tomleslie 

Thank you so much. This was exactly what I need...

@pooyan1990 at first, it was immedi...

pooyan1990 35
November 03 2019
0 0

@pooyan1990 

at first, it was immediately (2 second) after I press the enter button.

but when I impose assume to the entire code, it took several seconds (about 20).

@Carl Love several seconds about 20...

pooyan1990 35
November 03 2019
0 0

@Carl Love 

several seconds about 20

@Rouben Rostamian   thank you B,H...

pooyan1990 35
November 03 2019
0 0

@Rouben Rostamian  

thank you

B,H,L__1 are independent parameters and they are not relating to any parameter and I will give them different values for my work. they are actually my entrance parameters for my problem and the only way for me is to solve the ODE for symbolic B,H,L__1.

didn't work...

pooyan1990 35
November 03 2019
0 0

@Carl Love 

thank you for your answer

I did exactly what you said and I get the same error again!! "Error, (in dsolve) give the main variable as a second argument"

is it possible that the problem is with my maple?

Mechrect.mw

@Carl Love  thank you so much, wor...

pooyan1990 35
October 30 2019
0 0

@Carl Love 

thank you so much, worked for me

pic from the paper...

pooyan1990 35
October 21 2019
0 0

@tomleslie 

you can see the initial condition from the paper.

 

Initial Condition...

pooyan1990 35
October 21 2019
0 0

@tomleslie 

The initial condition from the paper for this differential equation is :

v__r(R__max(theta,beta),theta,beta)=0

but as I stated in my first post, I can't implement this initial condition to this DE.

The error is like that:
Error, (in dsolve) unexpected occurence of the variables {beta,theta} in the forst operand... 
and you can see that completely in my first post


Is there anyway to implement this initial condition to the equation?

thank you

Help plz...

pooyan1990 35
October 20 2019
0 0

@tomleslie 

I have tried many times and again I have done it in your way and doesn't obtain the correct answer.

when solving the differential equation, an integral comes up and doesn't solve.

I have attached the worksheet.

would you please take a look at this?

eq48a is from the paper and sol2 should be equal to eq48a.

the paper's link in removed because of copyright violaton. if you need the paper let me know to send it to you.

thank you for taking the time


Test2.mw
 

> 

d := 10; c := 10; L__1 := 4; R__i := (1/2)*d+L__1; R__f := c+(1/2)*d+L__1;
eq43:= v__beta(r,theta,beta)= v__m(beta)*(1-(r^2/(R__max(theta,beta)^2)));
eq37:=diff(v__r(r, theta, beta),r)*r*(R__f-r*cos(theta))+v__r(r,theta,beta)*(R__f-2*r*cos(theta))+r*diff(v__beta(r,theta,beta),beta)=0;

d := 10

c := 10

L__1 := 4

R__i := 9

R__f := 19

eq43 := v__beta(r, theta, beta) = v__m(beta)*(1-r^2/R__max(theta, beta)^2)

eq37 := (diff(v__r(r, theta, beta), r))*r*(19-r*cos(theta))+v__r(r, theta, beta)*(19-2*r*cos(theta))+r*(diff(v__beta(r, theta, beta), beta)) = 0

 (1)
> 

eq39:=R(beta)=R__i+2*beta*(R__f-R__i)/Pi;

eq39 := R(beta) = 9+20*beta/Pi

 (2)
> 

eq38:=R__max(theta,beta)=R__max(theta,0)*(R(beta)/R__i);
eq41:=R__max(theta,0)=((2*L__1*cos(theta)) + sqrt(d^2 - 2*(L__1)^2 + 2*(L__1)^2*cos(2*theta)))/2;

eq38 := R__max(theta, beta) = (1/9)*R__max(theta, 0)*R(beta)

eq41 := R__max(theta, 0) = 4*cos(theta)+(17+8*cos(2*theta))^(1/2)

 (3)
> 

eq38a:= subs(eq41, eq38);

eq38a := R__max(theta, beta) = (1/9)*(4*cos(theta)+(17+8*cos(2*theta))^(1/2))*R(beta)

 (4)
> 

eq38b:=simplify(subs(eq39,eq38a));

eq38b := R__max(theta, beta) = (1/9)*(9*Pi+20*beta)*(4*cos(theta)+(17+8*cos(2*theta))^(1/2))/Pi

 (5)
> 

eq42:=R__max(theta,beta)=( (((d+2*L__1)*Pi) - 2*beta*(d+2*L__1-2*R__f))*(((2*L__1*cos(theta)) + sqrt(d^2 - 2*(L__1)^2 + 2*(L__1)^2*cos(2*theta)))))/(2*Pi*(d+2*L__1));

eq42 := R__max(theta, beta) = (1/36)*(18*pi+40*beta)*(8*cos(theta)+2*(17+8*cos(2*theta))^(1/2))/pi

 (6)
> 

eq44 := v__m(beta) = 1/(int(int((1-r^2/R__max(theta, beta)^2)*r, r = 0 .. R__max(theta, beta)), theta = 0 .. 2*Pi))

eq44 := v__m(beta) = 4/(int(R__max(theta, beta)^2, theta = 0 .. 2*Pi))

 (7)
> 

eq44a:=simplify(expand(subs(eq38b,eq44)));

eq44a := v__m(beta) = (162/25)*Pi/(9*Pi+20*beta)^2

 (8)
> 

eq47:=subs([eq44a,eq38b],eq43);

eq47 := v__beta(r, theta, beta) = (162/25)*Pi*(1-81*r^2*Pi^2/((9*Pi+20*beta)^2*(4*cos(theta)+(17+8*cos(2*theta))^(1/2))^2))/(9*Pi+20*beta)^2

 (9)
> 

eq37a:=simplify(expand(subs(eq47, eq37)));

eq37a := (1/5)*(-2361960*(r*cos(theta)-19)*(cos(theta)*(9+16*cos(theta)^2)^(1/2)+4*cos(theta)^2+9/8)*r*(Pi+(20/9)*beta)^5*(diff(v__r(r, theta, beta), r))-4723920*(r*cos(theta)-19/2)*(cos(theta)*(9+16*cos(theta)^2)^(1/2)+4*cos(theta)^2+9/8)*(Pi+(20/9)*beta)^5*v__r(r, theta, beta)-839808*(4*(Pi+(20/9)*beta)^2*cos(theta)^2+(9+16*cos(theta)^2)^(1/2)*(Pi+(20/9)*beta)^2*cos(theta)+(-(1/4)*r^2+9/8)*Pi^2+5*Pi*beta+(50/9)*beta^2)*Pi*r)/((9*Pi+20*beta)^5*(32*cos(theta)^2+8*cos(theta)*(9+16*cos(theta)^2)^(1/2)+9)) = 0

 (10)
> 

sol1:=simplify(dsolve(eq37a, v__r(r, theta,beta)));

sol1 := v__r(r, theta, beta) = (1/5)*(-1296*Pi*(Int(-162*(-2*(cos(2*theta)*r+r-38*cos(theta))*(Pi+(20/9)*beta)^2*(17+8*cos(2*theta))^(1/2)+152*(Pi+(20/9)*beta)^2*cos(2*theta)-4*r*(Pi+(20/9)*beta)^2*cos(3*theta)+(Pi^2*r^2-(33/2)*(Pi+(20/9)*beta)^2)*r*cos(theta)-19*Pi^2*r^2+(475/2)*(Pi+(20/9)*beta)^2)*r/((4*cos(2*theta)*r-152*cos(theta)+4*r)*(17+8*cos(2*theta))^(1/2)+33*r*cos(theta)+8*r*cos(3*theta)-304*cos(2*theta)-475), r))+295245*_F1(theta, beta)*(Pi+(20/9)*beta)^5)/((9*Pi+20*beta)^5*(r*cos(theta)-19)*r)

 (11)
> 

sol2:=eval(sol1, _F1(theta, beta)=0);

sol2 := v__r(r, theta, beta) = -(1296/5)*Pi*(Int(-162*(-2*(cos(2*theta)*r+r-38*cos(theta))*(Pi+(20/9)*beta)^2*(17+8*cos(2*theta))^(1/2)+152*(Pi+(20/9)*beta)^2*cos(2*theta)-4*r*(Pi+(20/9)*beta)^2*cos(3*theta)+(Pi^2*r^2-(33/2)*(Pi+(20/9)*beta)^2)*r*cos(theta)-19*Pi^2*r^2+(475/2)*(Pi+(20/9)*beta)^2)*r/((4*cos(2*theta)*r-152*cos(theta)+4*r)*(17+8*cos(2*theta))^(1/2)+33*r*cos(theta)+8*r*cos(3*theta)-304*cos(2*theta)-475), r))/((9*Pi+20*beta)^5*(r*cos(theta)-19)*r)

 (12)
> 

eq0 := o__1 = (1-2*L__1/d)^2*(1+2*L__1/d+4*c*beta/(d*Pi))^2; eq00 := o__2 = 1-4*L__1*(cos(theta)*sqrt(1-4*L__1^2*sin(theta)^2/d^2)-L__1*cos(theta)/d)/d

eq0 := o__1 = (1/25)*(9/5+4*beta/Pi)^2

eq00 := o__2 = 1-(8/25)*cos(theta)*(25-16*sin(theta)^2)^(1/2)+(16/25)*cos(theta)

 (13)
> 

eq48 := v__r(r, theta, beta) = 8*(1+2*L__1/d)^2*c*r*(o__1-4*(r/d)^2*o__2)/(d^2*Pi*(1-2*L__1/d)^2*(1+2*L__1/d+4*c*beta/(d*Pi))^5*(R__f/d-r*cos(theta)/d))

eq48 := v__r(r, theta, beta) = (324/5)*r*(o__1-(1/25)*r^2*o__2)/(Pi*(9/5+4*beta/Pi)^5*(19/10-(1/10)*r*cos(theta)))

 (14)
> 

eq48a := simplify(expand(subs([eq0, eq00], eq48)))

eq48a := v__r(r, theta, beta) = -25920*Pi^2*(Pi^2*(9+16*cos(theta)^2)^(1/2)*cos(theta)*r^2-2*Pi^2*cos(theta)*r^2+(-(25/8)*r^2+81/8)*Pi^2+45*Pi*beta+50*beta^2)*r/((9*Pi+20*beta)^5*(r*cos(theta)-19))

 (15)
> 

%?

 

@tomleslie  Thank you for taking t...

pooyan1990 35
October 18 2019
0 0

@tomleslie 

Thank you for taking the time

but there some mistakes in this solution.

In this paper there two different mechanism and equations and you use equations for the first one but I want the second one.

V(beta) is not eq1 in the paper, it is eq.43

Our final answer for V(r) is not eq.2, it is eq.48

we also have an initial condition in page 10 of the paper which I highlighted that above eq.48.

I substitute these equations and didn't get the right answer.

would you please do this for me with these reforms?

I will be grateful for your help.

@ecterrab Thanks for your answerI U...

pooyan1990 35
October 17 2019
0 0

@ecterrab 

Thanks for your answer

I Understand this but I can't change my variable names. I extracted these formulations from a valid paper so I think they are correct. I have uploaded the relevant article pdf. you can see all of these equations in pages 9 and 10 from eq. 37 to 48. the boundry condition is also highlighted in page 10. I would be grateful if you take a brief look at this article and if you have any idea to solve this problem.

thanks again

[Edit: .pdf deleted because of copyright violation.]

@Christian Wolinski  Thanks for yo...

pooyan1990 35
September 11 2019
0 0

@Christian Wolinski 

Thanks for your answer

I have solved the problem. I just forgot to Multiply Beta with other terms and there was no sign of multiplying between them. sorry for the inconvenience.

Page 1 of 1