How to get the cellcell center to center distance dx dy dz

dbxmcf · August 7, 2008, 22:26

Hi, I am trying to implement vonNeumann and Ritchmyer's artificial viscosity for shock capture, which required a term delta^2, which is the center to center distance (or the vector(dx, dy, dz), the grid distance), how to obtain that with OpenFOAM's field operations?

Thanks a lot!

ngj · August 8, 2008, 04:23

Hi Roy

On patches the inverse distance is found as this->patch().deltaCoeffs(). You might find something similar for the interior.
Otherwise, you need to locate owner/neighbour properties for each face and then calculate the distance by hand.

The shock capturing sounds very interesting.

Good luck,

Niels

P.S. Have seen a surfaceScalarField definition of deltaCoeffs at some point, but cannot recall where or when.

hjasak · August 8, 2008, 04:37

mesh.deltaCoeffs() will work. It is in surfaceInterpolation.H.

Enjoy,

Hrv

August 7, 2008, 22:26	Hi, I am trying to implement v	#1
dbxmcf Member roy fokker Join Date: Mar 2009 Posts: 44 Rep Power: 17	Hi, I am trying to implement vonNeumann and Ritchmyer's artificial viscosity for shock capture, which required a term delta^2, which is the center to center distance (or the vector(dx, dy, dz), the grid distance), how to obtain that with OpenFOAM's field operations? Thanks a lot!

August 8, 2008, 04:23	Hi Roy On patches the inver	#2
ngj Senior Member Niels Gjoel Jacobsen Join Date: Mar 2009 Location: Copenhagen, Denmark Posts: 1,900 Rep Power: 37	Hi Roy On patches the inverse distance is found as this->patch().deltaCoeffs(). You might find something similar for the interior. Otherwise, you need to locate owner/neighbour properties for each face and then calculate the distance by hand. The shock capturing sounds very interesting. Good luck, Niels P.S. Have seen a surfaceScalarField definition of deltaCoeffs at some point, but cannot recall where or when. __________________ Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.

August 8, 2008, 04:37	mesh.deltaCoeffs() will work.	#3
hjasak Senior Member Hrvoje Jasak Join Date: Mar 2009 Location: London, England Posts: 1,905 Rep Power: 33	mesh.deltaCoeffs() will work. It is in surfaceInterpolation.H. Enjoy, Hrv mbookin likes this. __________________ Hrvoje Jasak Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
center of pressure	kartik	FLUENT	0	July 2, 2008 10:10
center of pressure	Filippo	FLUENT	2	March 29, 2007 02:48
Center of Pressure	Balder	Siemens	1	October 13, 2006 13:39
Center of pressure	edi	FLUENT	0	March 10, 2006 11:45
About Cell Center	Harry	FLUENT	0	December 3, 2000 04:33