Controlling magnetic domains in soft underlayer (SUL) of perpendicular magnetic recording (PMR) is an important issue for the application of PMR in HDD. We studied the magnetic domain structures in SUL using the finite element based micromagnetic simulation (FEMM) for the SUL models with different thicknesses. The purpose is to simulate the magnetic domain wall noise when the SUL thickness and saturation magnetization are changed. The simulation results show that a 15 nm SUL forms simpler Neel wall domain wall pattern and 40 nm SUL forms complex Bloch wall. To visualize the effect of these domain walls stray field at a read sensor position, the magnetic stray field of the domain walls at air bearing surface (ABS) which is 50 nm above the SUL was simulated and the results imply that Bloch walls have stronger stray field with more complicated field patterns than Neel walls and this becomes a significant noise source. Therefore, the thickness of the SUL should be controlled to avoid the formation of Bloch walls.