|
[Sponsors] |
Heat transfer in a FSI simulation by OpenFoam |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
July 15, 2019, 05:30 |
Heat transfer in a FSI simulation by OpenFoam
|
#1 |
New Member
peyman havaej
Join Date: Jan 2016
Posts: 16
Rep Power: 10 |
Dear Formers
I want to simulate fluid flow and heat transfer in a FSI case. I use Foam-extend 4.0 and FSI - package provided in openfoamwiki. An appropriate form of Energy equation was added to the fluid solver, and it works correctly. Also, for the solid region, Energy equation for solving Temperature was implemented in the main solver, named fsiFoam. Now, my problem is related to boundary condition between the fluid/solid interface for Temperature. After I searched a lot, I found that I can use ggi method for interpolating data between solid and fluid interface, especially in the non-conformal mesh. However, it did not work for me, since when I set ggi condition in /constant/polyMesh/boundary, the shadowPatch for the fluid interface, located in the solid folders, and the code can not access the boundary and faceZone which was in the solid part. Could someone give me a suggestion, please? I appreciate that. Thank you. Here is the boundaries for fluid: 6 ( inlet { type patch; nFaces 40; startFace 23570; } outlet { type patch; nFaces 40; startFace 23610; } bottumWalls { type wall; nFaces 360; startFace 23650; } topWalls { type wall; nFaces 280; startFace 24010; } interface_fluid { type ggi; nFaces 140; startFace 24290; shadowPatch interface_solid; zone interface-zone; bridgeOverlap false; } frontAndBack { type empty; nFaces 24000; startFace 24430; } ) And here for solid 3 ( topWalls { type wall; nFaces 75; startFace 4925; } interface_solid { type ggi; nFaces 75; startFace 5000; shadowPatch interface_fluid; zone interface-zone; bridgeOverlap false; } frontAndBack { type empty; nFaces 5000; startFace 5075; } After running program, this error was appeared: Create time Create dynamic mesh for time = 0 Selecting dynamicFvMesh dynamicMotionSolverFvMesh Selecting motion solver: velocityLaplacian Selecting motion diffusion: quadratic Selecting motion diffusion: inverseDistance --> FOAM FATAL ERROR: Shadow patch name interface_solid not found. Please check your GGI interface definition. From function label ggiPolyPatch::shadowIndex() const in file meshes/polyMesh/polyPatches/constraint/ggi/ggiPolyPatch.C at line 774. FOAM aborting Aborted (core dumped) constant.zip Last edited by peyman.havaej; July 15, 2019 at 06:52. |
|
July 16, 2019, 10:09 |
|
#2 |
Senior Member
Daniel
Join Date: Mar 2013
Location: Noshahr, Iran
Posts: 348
Rep Power: 21 |
Dear Peyman,
I have done exactly what you are trying to do before! You are right about using GGI interpolation, but it's not that easy. I have implemented everything needed in solids4Foam toolkit with a simple tutorial. But it's not yet validated nor completed. For now it is limited to the one of the fluid models but extending it to the other models should be easy and straightforward: 1- fluidModel: buoyantBoussinesqPimpleFluid 2- solidModel: thermalLinGeomSolid (fsi) and thermalSolid (no-fsi) Meanwhile I would suggest you contact Dr. Philip Cardiff (@bigphil) and request for access to the toolkit. There you can check "feature-coupledTemperatureField" branch for more information. Regards, D. Khazaei |
|
October 25, 2019, 02:01 |
|
#3 | |
Senior Member
Hojatollah Gholami
Join Date: Jan 2019
Posts: 171
Rep Power: 7 |
Dear All
I wanna to add a fixed solid region into fluid region of fsi problem to study heat transfer. For coupling heat transfer from fixed solid to fluid and invert, the chtMultiRegionFoam solver algorithm uses a loop and set fluidField, solve fluid region, set solidField and solve solid region for nOuterCorr times without calculating interface residual. As Quote:
In fsiFoam solver, we have elastic solid region and fluid region. the fluid region mesh in runTime will changes and the fixed solid region will have fixed mesh. So the mesh of two conjugate regions will be different. my opinion is to modify fluid solver of fsiFoam, as after solving UEqn and PEqn, start a loop with 3 times and calculate energy equation of fixed solid region, save T boundary, calculate energy equation in fluid region and save T boundary. But as the interface of fluid and solid regions will be different, we should use ggi mapping. My questions are, is it better suggestion? or if I want to use this algorithm, how can I implement ggi class? |
||
November 27, 2019, 02:26 |
|
#4 | |
Senior Member
Hojatollah Gholami
Join Date: Jan 2019
Posts: 171
Rep Power: 7 |
Dear Khazaei
for your recommend solvers, as there aren't any tutorial. I wanna to use coupled condition for solid and fluid regions. For example, for HronTurek benchmark which condition is suitable for plate of solid region and plate of fluid region? I think fixedValue, zeroGradient is not suitable. chtRcTemperature is not work. Thanks. Quote:
|
||
December 8, 2019, 03:56 |
|
#5 |
Senior Member
Hojatollah Gholami
Join Date: Jan 2019
Posts: 171
Rep Power: 7 |
Dear peyman
Do you solve this problem? |
|
Tags |
foam-extend 4.0, fsi, heat transfer, multi zones, openfoam |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Evaporation-Condensation in a porous zone (Heat transfer problem) | maximilian-1 | Fluent Multiphase | 1 | August 22, 2018 10:42 |
conjugate heat transfer in OpenFOAM | skuznet | OpenFOAM Running, Solving & CFD | 99 | March 16, 2017 06:07 |
OpenFOAM v3.0+ ?? | SBusch | OpenFOAM | 22 | December 26, 2016 15:24 |
Heat transfer in an unsteady-state simulation | Raed141 | FLUENT | 11 | August 7, 2009 18:17 |
Natural convection - Inlet boundary condition | max91 | CFX | 1 | July 29, 2008 21:28 |