Gravity modeled by customization does not conform with that is added in GUI

boesendorfer · June 15, 2017, 21:28

Gravity model has been included in physics panel of GUI.
However, I disabled this model in GUI and tried to account for gravity force by customization shown below because of some other reasons.

In the forcal.f, I loop over all fluid cell and set fsigz(ijk)=-9.8*rho(ijk)*delx(i)*dely(j)*delz(k)*f(ijk)
Also I've added iffor=1 in prepin file.

Then I run a simulation of stationary liquid with density of 7800 and depth of 0.01m. The mesh size is 0.001m.

The result shows that it does have effort. But the pressure at the bottom is about 15MPa higher than that is simulated by individual gravity model in GUI.

Why?

boy · September 12, 2017, 03:00

Hi,boesendorfer, Can you post your whole forcal subroutine here or send it to 1789078096@qq.com?

June 15, 2017, 21:28	Gravity modeled by customization does not conform with that is added in GUI	#1
boesendorfer New Member Join Date: Apr 2014 Posts: 17 Rep Power: 12	Gravity model has been included in physics panel of GUI. However, I disabled this model in GUI and tried to account for gravity force by customization shown below because of some other reasons. In the forcal.f, I loop over all fluid cell and set fsigz(ijk)=-9.8rho(ijk)delx(i)dely(j)delz(k)*f(ijk) Also I've added iffor=1 in prepin file. Then I run a simulation of stationary liquid with density of 7800 and depth of 0.01m. The mesh size is 0.001m. The result shows that it does have effort. But the pressure at the bottom is about 15MPa higher than that is simulated by individual gravity model in GUI. Why?

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
ANSYS Workbench on "Uncertified" Linux Distros	hsr	CFX	289	April 20, 2023 09:23
Customization using parqsr.f for mass momentum source	boesendorfer	FLOW-3D	0	May 28, 2017 09:30
Segmentation fault in layerAdditionRemoval Foam-extend 3.1	jbettencourt	OpenFOAM Running, Solving & CFD	3	March 10, 2016 05:53

September 12, 2017, 03:00		#2
boy New Member Guangxi Join Date: Sep 2014 Location: China Posts: 4 Rep Power: 11	Hi,boesendorfer, Can you post your whole forcal subroutine here or send it to 1789078096@qq.com?