[lesmodel]

Dear Max,
in aorta.cpp line174, it’s ‘SinusStartScale<T,int> nSinusStartScale( iTperiod,converter.getLatticeU() )’, the velocity at the inlet is a circular poiseuille profile with LatticeU beein the maximum velocity, it’s LatticeU not CharU. Looking forward to your reply.

[lesmodel]

Dear mathias,
Thanks very much for your quick response. I have another question, in the case of aorta3d. the Re is decided by CharU, CharL and CharNu, if I change the CharU to obtain the different Re, what’s the new velocity for inflow?

[lesmodel]

Dear Robin,
Thanks very much for your help, one more question, in the example cylinder3d, a planeReduction is difined to use “gifWriter.write” to generate the ppm files of velocity distribution, but when defining the planeReduction, only a
normal vector (1,0,0) is specified, how can we know which plane is referred to, that is to say, it refers to plan x=0 or someone esle? Thanks very much!!

[lesmodel]

Dear robin,
Thanks very much for your reply, but I have another question, when I run a program in parallel in the Cygwin64,
I use the command
mpirun -np 8 ./cylinder2d.exe
it reports:
[user-PC:98680] [[12126,1],0] tcp_peer_send_blocking: send() to socket 12failed:Transport endpoint is not connected (128)
[user-PC:98680] [[12126,1],0] tcp_peer_send_blocking: send() to socket 12failed:Transport endpoint is not connected (128)

What’s wrong with this, thans very much!!

[lesmodel]

Dear Mathias,
Thanks for you reply. How do we determine this threshold value, that is 0.3 here.

[lesmodel]

Dear Mathias,
Thanks very much for your kind reply. Do you mean that I should only change the charNU in the LBconverter to 0.000001m^2/s? and make sure that the Reynolds number is fine for the calculating, but what about the density, should I write the charRho = 1000 to the LBconverter? if not, can i get the correct force acting on the cylinder? Thank you!

[Author]

Posts