If you change your geometry every iteration, you will get negative cell volume very fast. The first few iterations every time step can have very high/low values of pressure and/or velocities which would result in unexpectedly large motion ... nBut anyway, every function in your UDF could be recalled by its name ... so, you can call the grid_motion function from inside define_adjust n