对称结构计算结果不对称原因

[好好科研]

计算一个对称的平板绕流，速度场和涡量大小的云图都对称，利用如下公式计算另一个变量的时候，计算出来的等值面不对称（应该是轻微不对称，但是我现在十分不对称），请问各位大佬，可能有什么原因呢。

	surfaceScalarField phiw
	(
		"phiw",
                fvc::flux(omega) * LL
	);
	

	for (int ntau = 1; ntau <= 5; ntau++)
	{
		H.oldTime() = H;
		fvScalarMatrix HEqn
		(
	             	  fvm::ddt(H)
			+ fvm::div(phiw,H)
			- fvm::laplacian(DH,H)
		);
		HEqn.solve();
	}

[好好科研]

算出来的等值面结果如图