<div dir="ltr"><div><div><div>Hi All,<br><br></div>In my understanding, reflective boundary condition means that Vn=0 at the boundary layer.<br></div>Here Vn means velocity normal to the boundary.<br><br></div>I noticed that in subroutine velocity_refl_vert_bc in Boundary_condition.c in CitcomCU, the reflective boundary condition is set by making Vx=0 and Szx=0 (taking x direction as an example):<br>
<br>if(E-&gt;parallel.me_loc[1] == 0)<br>                {<br>                    E-&gt;VB[1][node1] = 0.0;<br>                    if((ii != 1) &amp;&amp; (ii != E-&gt;mesh.noz))<br>                        E-&gt;VB[3][node1] = 0.0;<br>
                }<br><br>I wonder why it sets the reflective boundary condition like this.<br><br>How about commenting the lines about E-&gt;VB[3] or aading other lines like setting Syx=0, just like below? <br clear="all">
<div><div><div><div><br>if(E-&gt;parallel.me_loc[1] == 0)<br>
                {<br>
                    E-&gt;VB[1][node1] = 0.0;<br>
                }<br>
<br></div><div>OR<br><br>if(E-&gt;parallel.me_loc[1] == 0)<br>
                {<br>
                    E-&gt;VB[1][node1] = 0.0;<br>
                    if((ii != 1) &amp;&amp; (ii != E-&gt;mesh.noz))<br>
                        E-&gt;VB[3][node1] = 0.0;<br>                    if((XXXXXXXX) &amp;&amp; (XXXXXXXX))<br>

                        E-&gt;VB[2][node1] = 0.0;<br>
                }<br>
<br></div><div>Thank you very much!<br></div><div>-- <br>wish  you  happiness !<br>yours  sincerely,<br>Yang Ting<br>School of Earth and Space Science <br>University of Science and Technology of China<br>
</div></div></div></div></div>