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

Search Results

Register Blogs Community New Posts Updated Threads Search

Showing results 1 to 6 of 6
Search took 0.00 seconds.
Search: Posts Made By: bjohn
Forum: Main CFD Forum September 10, 2010, 23:08
Replies: 12
Views: 3,201
Posted By bjohn
>u_L=1.0 and u_R=-1.0 i.e. two opposing flows at...

>u_L=1.0 and u_R=-1.0 i.e. two opposing flows at the discontinuity

This is probably not a steady shock satisfying the Rankine-Hugoniot relation.
So, waves will come out. Sod' shock tube...
Forum: Main CFD Forum September 8, 2010, 10:57
Replies: 12
Views: 3,201
Posted By bjohn
The sample code looks good to me. Maybe, errors...

The sample code looks good to me.
Maybe, errors have been introduced when you convert it into C code...

Looks like you're trying to run Sod's shock tube problem.
There is a plot in the Wiki...
Forum: Main CFD Forum September 8, 2010, 05:50
Replies: 12
Views: 3,201
Posted By bjohn
>ρ_R = ρ_L = 1.0, p_L = p_R = 1.0, u = v = 0.0...

>ρ_R = ρ_L = 1.0, p_L = p_R = 1.0, u = v = 0.0 through out.

In this case, you have

f1Left = 0
f2Left = nx
f3Left = ny
f4Left = 0

Right flux
Forum: Main CFD Forum September 8, 2010, 05:02
Replies: 12
Views: 3,201
Posted By bjohn
>qR[i][j] =...

>qR[i][j] = u[i][j]*deltayR[i][j]-v[i][j]*deltaxR[i][j]; //Normal and
>rR[i][j] = u[i][j]*deltaxR[i][j]+v[i][j]*deltayR[i][j]; //tangent.

These quantities must be based on the unit vectors:
...
Forum: Main CFD Forum September 6, 2010, 06:26
Replies: 12
Views: 3,201
Posted By bjohn
1) Your (deltax,deltay) is tangent to the face....

1) Your (deltax,deltay) is tangent to the face. If you define deltax and deltay as you describe, normal will be (deltay, -deltax).

2) The vector multiplied by pressure should be the normal:
...
Forum: Main CFD Forum September 4, 2010, 23:25
Replies: 12
Views: 3,201
Posted By bjohn
I think it's correct that you say. Some...

I think it's correct that you say.

Some things to check:

1) Add up face normals for each cell. The result should be zero.
2) Set up a uniform flow over the domain and run the code.
The...
Showing results 1 to 6 of 6

 
Forum Jump

All times are GMT -4. The time now is 08:58.