Hi,

I am trying to assign different temperatures for different regions in a 3D simulation box. The temperature calculations should not include any drifting velocity. I am having some troubles to make it work. Would you please help?

Here is what I did,

#assign region:

region 25o cylinder z 0 0 25 -40 40 side out units box

region 26i cylinder z 0 0 26 -40 40 side in units box

region 26o cylinder z 0 0 26 -40 40 side out units box

region 26 intersect 2 26i 25o units box

region 27i cylinder z 0 0 27 -40 40 side in units box

region 27o cylinder z 0 0 27 -40 40 side out units box

region 27 intersect 2 27i 26o units box

region 28i cylinder z 0 0 28 -40 40 side in units box

region 28o cylinder z 0 0 28 -40 40 side out units box

region 28 intersect 2 28i 27o units box

region 29i cylinder z 0 0 29 -40 40 side in units box

region 29 intersect 2 29i 28o units box

#temperature calculations:

fix d26 drop1 temp/rescale 1 1.3 1.3 0.01 1 region 26 # assign temperatures

compute_modify d26_temp extra 3 #take off drifting velocities

compute cr26 drop1 temp/region 26 #checking temperature with compute command

compute_modify cr26 extra 3 #take off drifting velocities

fix d27 drop1 temp/rescale 1 1.4 1.4 0.01 1 region 27

compute_modify d27_temp extra 3

compute cr27 drop1 temp/region 27

compute_modify cr27 extra 3

fix d28 drop1 temp/rescale 1 1.5 1.5 0.01 1 region 28

compute_modify d28_temp extra 3

compute cr28 drop1 temp/region 28

compute_modify cr28 extra 3

fix d29 drop1 temp/rescale 1 1.6 1.6 0.01 1 region 29

compute_modify d29_temp extra 3

compute cr29 drop1 temp/region 29

compute_modify cr29 extra 3

#output results

thermo_style custom step c_d26_temp c_d27_temp c_d28_temp c_d29_temp c_cr26 c_cr27 c_cr28 c_cr29

AND here are the results,

1495000 1.6 1.6 1.6 1.6 1.3736374 1.3389688 1.3441566 1.3349029

1500001 1.6 1.6 1.6 1.6 1.3743383 1.3437465 1.3394622 1.3441895

I thought,

1) _d26_temp c_d27_temp c_d28_temp c_d29 should be about 1.3, 1.4, 1.5 and 1.6

2)c_d2*_temp should have the same value as c_cr2* if they are in the same region.

What did I do wrong? Any help is highly appreciated.

Thanks,

Shihai