CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

rotational and inviscid

Register Blogs Community New Posts Updated Threads Search

Like Tree7Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 3, 2023, 16:33
Default
  #21
Senior Member
 
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,761
Rep Power: 66
LuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura about
So now you have the initial value problem.

If the initial value is 0. The change when the initial value is 0 is 0. Hence, when it is 0, it stays 0.


There are numerous analogues. The simple harmonic oscillator for example is y''=y. If you put y(t=0)=0, then y''=y'=y=0. A simple harmonic oscillator doesn't oscillate if there is no initial oscillation even though the system is capable of oscillating.
LuckyTran is offline   Reply With Quote

Old   November 3, 2023, 17:15
Default
  #22
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by ichamail View Post
Let me use some different notation just in case we implie the same thing but we "lost in translation"

\frac{d \vec{\zeta} }{dt}(t) = \left( \vec{\zeta}(t) \cdot \nabla \right) \vec{V}(t)

If we know that \vec{\zeta}(t_0) = 0 then, all I can say is

\frac{d \vec{\zeta} }{dt}(t_0) = \left( \vec{\zeta}(t_0) \cdot \nabla \right) \vec{V}(t_0) = 0

This is all I can understand. If you are telling me that there is something more than that, I'm sorry but maybe I need to go back to the basics, and the come back again to ask for help



Ok, let use the equation and integrate exactly in time




z(t)=z(t0) + Int[t0,t] z(tau).Grad v(tau) dtau




note that we integrate along the path-line dx/dt=v[x(t),t].


Thus, z(t) is a non-vanishing function if along the path-line there is a non-zero contribution of the integral.
ichamail likes this.
FMDenaro is offline   Reply With Quote

Old   November 3, 2023, 17:15
Default
  #23
New Member
 
Join Date: Jul 2022
Posts: 20
Rep Power: 4
ichamail is on a distinguished road
Quote:
Originally Posted by LuckyTran View Post
So now you have the initial value problem.

If the initial value is 0. The change when the initial value is 0 is 0. Hence, when it is 0, it stays 0.


There are numerous analogues. The simple harmonic oscillator for example is y''=y. If you put y(t=0)=0, then y''=y'=y=0. A simple harmonic oscillator doesn't oscillate if there is no initial oscillation even though the system is capable of oscillating.
I have to disagree on that. For the case of simple harmonic oscillator you need two initial conditions. You can have y(0)=0 and y'(0) = 1, where 1 would be the maximum speed. On the other hand y=0 is definately a solution of this ode and it means that when you are at rest you stay at rest. But a solution must satisfy both the equation and and the initial values.

But for the case of vorticity that you need only one initial vvalue I agree
ichamail is offline   Reply With Quote

Old   November 3, 2023, 18:36
Default
  #24
Senior Member
 
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,761
Rep Power: 66
LuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura about
So now it is suddenly obvious how y(t=0)=0 does not necessarily imply y'(t=0)=0?

Interesting.....

So do you get it now?
LuckyTran is offline   Reply With Quote

Old   November 4, 2023, 13:53
Default
  #25
New Member
 
Join Date: Jul 2022
Posts: 20
Rep Power: 4
ichamail is on a distinguished road
Quote:
Originally Posted by LuckyTran View Post
So now it is suddenly obvious how y(t=0)=0 does not necessarily imply y'(t=0)=0?

Interesting.....

So do you get it now?
I never said y(t=0)= 0 \implies y'(t=0) = 0. I said that y = 0 \implies  y' = 0. These things are different. y=0 is the same thing as y(t)=0 that implies y' = y'(t) = 0. Just like y, which is not an independent variable here but a function of time (y = y(t)), vorticity \vec{\zeta} is a function of time and space (\vec{\zeta} =\vec{\zeta}(x,y,z,t)), and not an independent variable. So setting \vec{\zeta} = 0, is like saying that \vec{\zeta}(x,y,z,t) = 0 which definately implies that \frac{d \vec{\zeta}(t)}{dt} = 0. One the other hand setting \vec{\zeta}(t=0,x,y,z) = 0 does not implies that \frac{d\zeta(t)}{dt} = 0, or \frac{d\zeta}{dt}(t=0) = 0 (assuming of course that you don't know the pde)

The point is that now I understand what you say, about knowing the initial value \vec{\zeta}(0, x,y,z) = 0 , and then calculating from the equation of vorticity the time derivative at t=0 \frac{d \vec{\zeta} }{dt}(t_0) = \left( \vec{\zeta}(t_0) \cdot \nabla \right) \vec{V}(t_0) = 0 can lead as to the colclusion that \frac{d \vec{\zeta} }{dt}(t) = 0


I just think, the way I would explain it now is this:

we have this pde
\frac{d \vec{\zeta} }{dt}(t) = \left( \vec{\zeta}(t) \cdot \nabla \right) \vec{V}(t)


\vec{\zeta}(t) = \vec{0} is a solution to this pde, because it satisfies the equation. \frac{d \vec{\zeta} }{dt}(t) - \left( \vec{\zeta}(t) \cdot \nabla \right) \vec{V}(t) =  \frac{d \vec{0} }{dt} - \left( \vec{0} \cdot \nabla \right) \vec{V}(t) = \vec{0} - \vec{0} = \vec{0}
THOUGH, It doesn't satisfies EVERY initial condition.

In case that the initial condition is \vec{\zeta}(t_0) = 0, then this solution satisfies both the pde and the initial condition.
So we can say \vec{\zeta}(t) = \vec{0} \implies  \frac{d \vec{\zeta} }{dt}(t) = 0
ichamail is offline   Reply With Quote

Old   November 4, 2023, 14:11
Default
  #26
Senior Member
 
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,761
Rep Power: 66
LuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura about
That point is that you don't know y'(t=0) is anything until you apply it to the PDE regardless of what y(t=0). You are saying it... I really don't get why you don't understand your own words and why you're confusing yourself.

Note that these two statements are not equivalent:
\vec{\zeta}(t) = \vec{0} \implies  \frac{d \vec{\zeta} }{dt}(t) =  0 is trivial because as you write \vec{\zeta}(t)  = \vec{0} \implies  \frac{d \vec{0} }{dt}(t) = 0. It's trivial because this is true , everywhere, all the time and has nothing to do with PDE's. The derivative of a constant function is 0.

What you want to write is:
\vec{\zeta}(t=t_0) = \vec{0} \implies  \frac{d \vec{\zeta} }{dt}(t) =  0 for any/all real t. This statement is only true due to the vorticity transport equation for this particular case of inviscid irrotational flow and under the condition that the vorticity was zero at least once.

I know I didn't fully write out in detailed explicit notation what was meant by my words, but you cannot fancy it up and improperly insert notation because you can end up writing something inconsistent.
LuckyTran is offline   Reply With Quote

Old   November 4, 2023, 14:14
Default
  #27
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by ichamail View Post
I never said y(t=0)= 0 \implies y'(t=0) = 0. I said that y = 0 \implies  y' = 0. These things are different. y=0 is the same thing as y(t)=0 that implies y' = y'(t) = 0. Just like y, which is not an independent variable here but a function of time (y = y(t)), vorticity \vec{\zeta} is a function of time and space (\vec{\zeta} =\vec{\zeta}(x,y,z,t)), and not an independent variable. So setting \vec{\zeta} = 0, is like saying that \vec{\zeta}(x,y,z,t) = 0 which definately implies that \frac{d \vec{\zeta}(t)}{dt} = 0. One the other hand setting \vec{\zeta}(t=0,x,y,z) = 0 does not implies that \frac{d\zeta(t)}{dt} = 0, or \frac{d\zeta}{dt}(t=0) = 0 (assuming of course that you don't know the pde)

The point is that now I understand what you say, about knowing the initial value \vec{\zeta}(0, x,y,z) = 0 , and then calculating from the equation of vorticity the time derivative at t=0 \frac{d \vec{\zeta} }{dt}(t_0) = \left( \vec{\zeta}(t_0) \cdot \nabla \right) \vec{V}(t_0) = 0 can lead as to the colclusion that \frac{d \vec{\zeta} }{dt}(t) = 0


I just think, the way I would explain it now is this:

we have this pde
\frac{d \vec{\zeta} }{dt}(t) = \left( \vec{\zeta}(t) \cdot \nabla \right) \vec{V}(t)


\vec{\zeta}(t) = \vec{0} is a solution to this pde, because it satisfies the equation. \frac{d \vec{\zeta} }{dt}(t) - \left( \vec{\zeta}(t) \cdot \nabla \right) \vec{V}(t) =  \frac{d \vec{0} }{dt} - \left( \vec{0} \cdot \nabla \right) \vec{V}(t) = \vec{0} - \vec{0} = \vec{0}
THOUGH, It doesn't satisfies EVERY initial condition.

In case that the initial condition is \vec{\zeta}(t_0) = 0, then this solution satisfies both the pde and the initial condition.
So we can say \vec{\zeta}(t) = \vec{0} \implies  \frac{d \vec{\zeta} }{dt}(t) = 0





Why? The derivative is still a function of time, you cannot assess in general this conclusion.
FMDenaro is offline   Reply With Quote

Old   November 4, 2023, 14:53
Default
  #28
New Member
 
Join Date: Jul 2022
Posts: 20
Rep Power: 4
ichamail is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Why? The derivative is still a function of time, you cannot assess in general this conclusion.

It is not a general conclusion. I prove it down below
ichamail is offline   Reply With Quote

Old   November 4, 2023, 14:59
Default
  #29
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by ichamail View Post
It is not a general conclusion. I prove it down below



I don't agree with the part highlighted in red. You can easily see that


z(t)=z(t0)+(t-t0)*dz/dt|t0 + (t-t0)^2/2 *d^2z/dt^2|t0 + ...




Therefore you have to assess that all derivatives are zero, not only the first derivative.
In your case, dz/dt|t0=0 does not imply that dz/dt(t) is zero for all time.
FMDenaro is offline   Reply With Quote

Old   November 4, 2023, 15:10
Default
  #30
New Member
 
Join Date: Jul 2022
Posts: 20
Rep Power: 4
ichamail is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
I don't agree with the part highlighted in red. You can easily see that


z(t)=z(t0)+(t-t0)*dz/dt|t0 + (t-t0)^2/2 *d^2z/dt^2|t0 + ...




Therefore you have to assess that all derivatives are zero, not only the first derivative.
In your case, dz/dt|t0=0 does not imply that dz/dt(t) is zero for all time.
Do you disagree with the proof I wrote under the highlighted part?

In the red highlighted part l didn't mean tha l came to this conclusion without the part bellow.
ichamail is offline   Reply With Quote

Old   November 4, 2023, 15:20
Default
  #31
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
You wrote:


Quote:
Originally Posted by ichamail View Post

The point is that now I understand what you say, about knowing the initial value \vec{\zeta}(0, x,y,z) = 0 , and then calculating from the equation of vorticity the time derivative at t=0 \frac{d \vec{\zeta} }{dt}(t_0) = \left( \vec{\zeta}(t_0) \cdot \nabla  \right) \vec{V}(t_0) = 0 can lead as to the colclusion that \frac{d \vec{\zeta} }{dt}(t) = 0


I just think, the way I would explain it now is this:

we have this pde
\frac{d \vec{\zeta} }{dt}(t) = \left( \vec{\zeta}(t) \cdot \nabla \right) \vec{V}(t)


\vec{\zeta}(t) = \vec{0} is a solution to this pde, because it satisfies the equation. \frac{d \vec{\zeta} }{dt}(t) - \left( \vec{\zeta}(t) \cdot \nabla  \right) \vec{V}(t) =  \frac{d \vec{0} }{dt} - \left( \vec{0} \cdot  \nabla \right) \vec{V}(t) = \vec{0} - \vec{0} = \vec{0}
THOUGH, It doesn't satisfies EVERY initial condition.

In case that the initial condition is \vec{\zeta}(t_0) = 0, then this solution satisfies both the pde and the initial condition.
So we can say \vec{\zeta}(t) = \vec{0} \implies  \frac{d \vec{\zeta} }{dt}(t) = 0


I wrote the Taylor expansion in time that shows you can have a non-vanishing functions z(t) and z'(t) even if z(t0)=0 and z'(t0)=0.
FMDenaro is offline   Reply With Quote

Old   November 4, 2023, 15:25
Default
  #32
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
The key is that you want to start assuming z(t)=0. Actually, we have only the equation and the initial condition at t0, no more conditions are assumed.
FMDenaro is offline   Reply With Quote

Old   November 4, 2023, 15:36
Default
  #33
New Member
 
Join Date: Jul 2022
Posts: 20
Rep Power: 4
ichamail is on a distinguished road
Quote:
Originally Posted by LuckyTran View Post
That point is that you don't know y'(t=0) is anything until you apply it to the PDE regardless of what y(t=0). You are saying it... I really don't get why you don't understand your own words and why you're confusing yourself.

Note that these two statements are not equivalent:
\vec{\zeta}(t) = \vec{0} \implies  \frac{d \vec{\zeta} }{dt}(t) =  0 is trivial because as you write \vec{\zeta}(t)  = \vec{0} \implies  \frac{d \vec{0} }{dt}(t) = 0. It's trivial because this is true , everywhere, all the time and has nothing to do with PDE's. The derivative of a constant function is 0.

What you want to write is:
\vec{\zeta}(t=t_0) = \vec{0} \implies  \frac{d \vec{\zeta} }{dt}(t) =  0 for any/all real t
. This statement is only true due to the vorticity transport equation for this particular case of inviscid irrotational flow and under the condition that the vorticity was zero at least once.

I know I didn't fully write out in detailed explicit notation what was meant by my words, but you cannot fancy it up and improperly insert notation because you can end up writing something inconsistent.
I never said that.

I know saying \vec{\zeta}(t) = \vec{0} is a trivial solution. But if it satisfy my initial condition, it is the solution
ichamail is offline   Reply With Quote

Old   November 4, 2023, 15:43
Default
  #34
New Member
 
Join Date: Jul 2022
Posts: 20
Rep Power: 4
ichamail is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
The key is that you want to start assuming z(t)=0. Actually, we have only the equation and the initial condition at t0, no more conditions are assumed.
I never assumed tha z(t) = 0.
I said z(t)=0 is a solution, because it satisfies the equation.
If z(t)=0 satisfies the initial condition too, it is the solution
ichamail is offline   Reply With Quote

Old   November 4, 2023, 15:44
Default
  #35
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by ichamail View Post
I never said that.

I know saying \vec{\zeta}(t) = \vec{0} is a trivial solution. But if it satisfy my initial condition, it is the solution

No, the mathematical problem is: PDE+intial conditions+boundary conditions. You cannot get a solution only if you satisfy the intial conditions.

The vorticity equation for ideal fluid is hyperbolic and requires BCs on an open space-time domain. Your trivial solution must satisfy also the BC on the frontiers.
FMDenaro is offline   Reply With Quote

Old   November 4, 2023, 15:48
Default
  #36
New Member
 
Join Date: Jul 2022
Posts: 20
Rep Power: 4
ichamail is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
No, the mathematical problem is: PDE+intial conditions+boundary conditions. You cannot get a solution only if you satisfy the intial conditions.

The vorticity equation for ideal fluid is hyperbolic and requires BCs on an open space-time domain. Your trivial solution must satisfy also the BC on the frontiers.
Ok, I see...what type of bc do we need. We know that vorticity of undisturbed flow at boundaries is zero for all time. Can we use that as BC?
ichamail is offline   Reply With Quote

Old   November 4, 2023, 16:14
Default
  #37
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Of course, if you set a zero initial condition and zero BC on the frontier.
FMDenaro is offline   Reply With Quote

Old   November 4, 2023, 16:23
Default
  #38
New Member
 
Join Date: Jul 2022
Posts: 20
Rep Power: 4
ichamail is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Of course, if you set a zero initial condition and zero BC on the frontier.
I had in my mind flows around bodies, so such a BC would always be true I guess.

Thinking of a body starting to move at t=0, then z is everywhere zero at t=0 and after that is zero far away from the body at all directions for t>0.
ichamail is offline   Reply With Quote

Old   November 4, 2023, 16:30
Default
  #39
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by ichamail View Post
I had in my mind flows around bodies, so such a BC would always be true I guess.

Thinking of a body starting to move at t=0, then z is everywhere zero at t=0 and after that is zero far away from the body at all directions for t>0.

You are in the field of the classic aerodynamics, now. Since you have slip condition on a body, the vorticity on a wall is not regular.

You have to think to inviscid flow problems that are solved also with vorticity distribution.

Lift is generated this way.
ichamail likes this.
FMDenaro is offline   Reply With Quote

Old   November 4, 2023, 17:41
Default
  #40
Senior Member
 
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,761
Rep Power: 66
LuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura about
Quote:
Originally Posted by ichamail View Post
I never assumed tha z(t) = 0.
I said z(t)=0 is a solution, because it satisfies the equation.
If z(t)=0 satisfies the initial condition too, it is the solution

The thing is...this is not a proof. You are assuming the conclusion as a premise which is circular reasoning. It is a common mistake in inductive arguments so don't feel too bad about it.

Btw the proof is not actually complicated nor profound. Filippo has basically demonstrated it for you, it falls right out after you apply the fundamental theorem of calculus. The crux is that you need to prove that the lagrangian derivative is zero for all time.
ichamail likes this.
LuckyTran is offline   Reply With Quote

Reply


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
interFoam inviscid and rotational pablodecastillo OpenFOAM Running, Solving & CFD 7 June 16, 2012 09:15
Convergence for inviscid flow over sylinder Rjakk Main CFD Forum 2 March 21, 2007 11:53
Inviscid flow Atit Koonsrisuk CFX 12 January 2, 2003 13:40
Separation by Euler equations Jitendra Main CFD Forum 9 July 5, 2000 15:58


All times are GMT -4. The time now is 17:54.