Numerical modeling for primary atomization of liquid jets