Output force distribution for each face on a patch

meh · October 19, 2011, 11:21

Hello.

I'm really sorry if this question has already been asked, I couldn't locate the solution in the forums.

I can output the integrated pressure and viscous forces on a patch using libforces.so. I've also managed to output pressures and stresses per face using the sample function.

What I need is the forces per face on a patch, together with the coordinates and area of each face. I'll be using this to input to a fluid structure interaction model.

Can anyone advise how to output the distribution of forces on a patch, rather than just the integrated value?

Thanks very much for your time

Matt

boger · October 20, 2011, 09:25

I'm not sure I understand. If you're able to output the pressure and viscous stress per face, then don't you just need to multiply by the face area?

meh · October 20, 2011, 09:27

Hi David

Thanks very much for your reply. Using sample i've only managed to output a vtk file, which doesn't have face areas in it. Can I output face area using sample?

Thanks very much

Matt

boger · October 20, 2011, 09:37

Ah - now I understand, but I don't know the answer. I assumed you were modifying a utility like patchIntegrate or wallShearStress, in which case you should be able to modify the utility slightly to write out what you need. Someone else probably has a more elegant solution.

meh · October 20, 2011, 11:07

I'll have a look at these utilities, thanks. If there any other ideas out there I'd be glad to hear them.

October 19, 2011, 11:21	Output force distribution for each face on a patch	#1
meh New Member Matt Join Date: Mar 2009 Posts: 4 Rep Power: 17	Hello. I'm really sorry if this question has already been asked, I couldn't locate the solution in the forums. I can output the integrated pressure and viscous forces on a patch using libforces.so. I've also managed to output pressures and stresses per face using the sample function. What I need is the forces per face on a patch, together with the coordinates and area of each face. I'll be using this to input to a fluid structure interaction model. Can anyone advise how to output the distribution of forces on a patch, rather than just the integrated value? Thanks very much for your time Matt

October 20, 2011, 09:25		#2
boger Senior Member David Boger Join Date: Mar 2009 Location: Penn State Applied Research Laboratory Posts: 146 Rep Power: 17	I'm not sure I understand. If you're able to output the pressure and viscous stress per face, then don't you just need to multiply by the face area? __________________ David A. Boger

October 20, 2011, 09:37		#4
boger Senior Member David Boger Join Date: Mar 2009 Location: Penn State Applied Research Laboratory Posts: 146 Rep Power: 17	Ah - now I understand, but I don't know the answer. I assumed you were modifying a utility like patchIntegrate or wallShearStress, in which case you should be able to modify the utility slightly to write out what you need. Someone else probably has a more elegant solution. __________________ David A. Boger

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[Commercial meshers] Fluent msh and cyclic boundary	cfdengineering	OpenFOAM Meshing & Mesh Conversion	48	January 25, 2013 03:28
[Commercial meshers] Using starToFoam	clo	OpenFOAM Meshing & Mesh Conversion	33	September 26, 2012 04:04
channelFoam for a 3D pipe	AlmostSurelyRob	OpenFOAM	3	June 24, 2011 13:06
gmsh2ToFoam	sarajags_89	OpenFOAM	0	November 24, 2009 22:50
[mesh manipulation] Using createPatch in place of couplePatches	sripplinger	OpenFOAM Meshing & Mesh Conversion	8	November 13, 2009 07:14

October 20, 2011, 09:27		#3
meh New Member Matt Join Date: Mar 2009 Posts: 4 Rep Power: 17	Hi David Thanks very much for your reply. Using sample i've only managed to output a vtk file, which doesn't have face areas in it. Can I output face area using sample? Thanks very much Matt

October 20, 2011, 11:07		#5
meh New Member Matt Join Date: Mar 2009 Posts: 4 Rep Power: 17	I'll have a look at these utilities, thanks. If there any other ideas out there I'd be glad to hear them.