CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

Forces from FVM to FEM conversion (fluid/structure interaction)

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes
  • 1 Post By FMDenaro

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 24, 2017, 07:15
Post Forces from FVM to FEM conversion (fluid/structure interaction)
  #1
New Member
 
A.MAROUF
Join Date: Dec 2016
Posts: 22
Rep Power: 9
Abderahmane is on a distinguished road
Dear all,

i am using an opensource code based on Finite Volume Method, well i want to make an interaction with a FEM software ( Abaqus , hyperworks, Ansys ...) the idea is to have the forces from FVM and go to the FEM do the Modal Analysis, get the mode shapes and come back the the FVM code to run the simulations.

The question is that i want to :
- convert FVM forces are in cells to forces in nodes for FEM

is there any idea how to make it ?


Thanks
Abderahmane is offline   Reply With Quote

Old   February 24, 2017, 08:40
Default
  #2
Senior Member
 
Michael Prinkey
Join Date: Mar 2009
Location: Pittsburgh PA
Posts: 363
Rep Power: 25
mprinkey will become famous soon enough
Your FVM code will (likely) have cell pressures and cell pressure gradients available. The most direct way to approach this would be to use those two piece of information to interpolate from the cell data to the node data. Since your FEM code is (likely) continuous, you will have multiple cells providing pressure data to the same node. You can compute them all and average them. Or you can assemble all of the cell data touching each node and build a multidimensional interpolant (say, least-squares).

The easy way to do it is to just average the pressures in the cells touching the node. Most wall boundary conditions assume dp/dn = 0, so the cell pressure is likely the same as the pressure at the wall. You just really need to average out the wall tangent pressure variations. I'd do that first. Then come back and see if including the wall normal gradient terms would make a difference.
mprinkey is offline   Reply With Quote

Old   February 24, 2017, 09:11
Default
  #3
New Member
 
A.MAROUF
Join Date: Dec 2016
Posts: 22
Rep Power: 9
Abderahmane is on a distinguished road
First of all thank you for your aswer, i already had this idea to make an interpolation between nodes and cells. But we may lose precision. What i know is that for example Workbench ansys does it from Mechanical APDL to Fluent ) , maybe there are some tools for it. Don't you think ?

Thanks
Abderahmane is offline   Reply With Quote

Old   February 24, 2017, 09:33
Default
  #4
Senior Member
 
Michael Prinkey
Join Date: Mar 2009
Location: Pittsburgh PA
Posts: 363
Rep Power: 25
mprinkey will become famous soon enough
There is no magic here. Pressure (as it is applied in FVM) is done via interpolation to the faces and and single point (piecewise constant) quadrature. That is how the hydrodynamics are computed. Perhaps the "most" accurate form is to interpolate the pressure to face centroid and consider that as the face load applied to the adjointing element. That would be the most consistent with the FVM formulation from a momentum conservation point of view. Then, I'd refer you to the finite element literature to find the best way to distributed forces to an element. For example, the node pressures should (in an intergral sense) preserve the total force applied in both formulations. I think you can use area-weighted averaging to accomplish this. Something like:

P_node = sum(A_face P_face) / sum(A_face)

But, you need to verify that the integral over the same surface using your finite element and P_node data matches sum(A_face P_face). My guess is that if you are using linear elements, it should work out. But I haven't done this before, so you really need to investigate it carefully.
mprinkey is offline   Reply With Quote

Old   February 24, 2017, 09:45
Default
  #5
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,793
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
I assume that for FV method you mean that you are working with the integral form of the NS equations. Therefore, you need to evaluate the integral of the total stress along each face of the FV. The total stress is defined by the pressure (-pn) and the tangential stress n.(2mu S).
Therefore you should match such quantities from the FE solution. I don't know if you have a sort of immersed boundary method to link the stresses.
I had some experience with my FV code interfaced with the elastic stresses of a membrane computed in a FE manner.
You could see some details in https://www.researchgate.net/publica...ids?ev=prf_pub
mprinkey likes this.
FMDenaro is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
CFD by anderson, chp 10.... supersonic flow over flat plate varunjain89 Main CFD Forum 18 May 11, 2018 07:31
3D interpolation: From FVM to FEM nodes and back vasilis Main CFD Forum 2 October 15, 2008 03:33
FVM vs FEM elwolv Main CFD Forum 8 June 5, 2007 13:23
comments on FDM, FEM, FVM, SM, SEM, DSEM, BEM kenn Main CFD Forum 2 July 18, 2004 18:28
FEM or FVM for CFD Astrid Main CFD Forum 18 December 15, 2000 00:02


All times are GMT -4. The time now is 14:17.