Let's Omega= [0, x_max][0, y_max][0, z_max] the domain where some particles move in. I have a question about the update of the particles locations if they across the domain boundaries. More precisely, I want the particles to be reflected inside the domain Omega if they across the limits of any axis. For example if the new location is {0.5,0.2,1.3} with x_max=y_max=z_max=1.0 so here z.location is outside the limit so how to update the particle location in that case with reflection? Thank you in advance Best regards Ilhem

Let's consider Omega= [0, x_max]*[0, y_max]*[0, z_max] the domain in which some particles move. I have a question about the update of the particles locations if they cross the domain boundaries. 

More precisely, I want the particles to be reflected inside the domain Omega if they across the limits of any axis. For example if the new location is {0.5, 0.2, 1.3} and $x_{max} = y_{max}=z_{max}=1.0$, then in this case the $z$ location is outside the limit, so how to update the particle location in that case with reflection?