|
[Sponsors] |
2 zones mesh morphing (to keep the boundary layer) |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
May 10, 2013, 06:24 |
2 zones mesh morphing (to keep the boundary layer)
|
#1 |
Senior Member
HECKMANN Frédéric
Join Date: Jul 2010
Posts: 249
Rep Power: 16 |
Dear everybody,
I'm working on a simulation where I have a deforming wall. It is actually an airfoil with a deforming leading edge (the law depends of the flow condition). I would like to know if there is a function / tool in OpenFoam (2.2 if possible) able to keep the boundary layer next to the wall intact and send the deformation in the domain. I know that Fluent uses a 2 zone definition to do so. The first zone next to the wall follow the deformation while the second zone takes the deep compression / extension deformation. In my case the boundary layer is quite important and I need to use a y+ < 1. I know use this options for the dynamicMeshDict: Code:
dynamicFvMesh dynamicMotionSolverFvMesh; motionSolverLibs ("libfvMotionSolvers.so"); solver displacementLaplacian; displacementLaplacianCoeffs { diffusivity inverseDistance (wall); } Does anybody have experience in this kind of simulation ? |
|
May 2, 2016, 11:03 |
|
#2 |
Member
Join Date: Mar 2013
Posts: 98
Rep Power: 13 |
Hi,
I have the same problems that you described and I wonder if you found some reasonable solution for that. I move a boundary based on a calculated field and then I apply a mesh movement solver. It works nicely without the boundary layer mesh but negative cells appears after few iterations if I build a mesh with y+ around 1. I tried all the available solvers in OpenFOAM-2.3.0 but none of them seems to give nice results. Thanks Giacomo |
|
May 11, 2016, 03:12 |
|
#3 |
Senior Member
HECKMANN Frédéric
Join Date: Jul 2010
Posts: 249
Rep Power: 16 |
Hi,
Maybe your boundary motion is too large compared to the first cell height. Try to decrease your time step until the maximum move is less than 20% of the first cell height. You can also try to use other algorithm that sends the deformation further in the domain. For example you can use a strong spring factor. |
|
May 11, 2016, 17:31 |
|
#4 |
Member
Join Date: Mar 2013
Posts: 98
Rep Power: 13 |
Hi,
first of all thanks for your reply. You are right, my boundary motion is certainly larger than the first cell height. The solvers work properly reducing the boundary motion but at a certain moment it crash in any case. Moreover, the boundary motion that I impose to make it works is too small for my application, increasing enormously the computational time. For this reason, I'm looking for a more stable mesh movement solver. The idea to keep a region around the body that follow the boundary deformation (leaving the mesh morphing to the remaining domain) is nice but I didn't manage to make it work (neither with the 6DOF solver). The other option is the one that you pointed out ("to use other algorithm that sends the deformation further in the domain"). When you speak about the strong spring factor, do you refer to some other solver that has not been implemented in OpenFOAM yet?Or that condition can be forced somehow in the case directories? I'll also try all the option available in the extended version. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Wind turbine simulation | Saturn | CFX | 58 | July 3, 2020 01:13 |
[snappyHexMesh] How to Do External Mesh for Airfoil sHM | msuaeronautics | OpenFOAM Meshing & Mesh Conversion | 1 | September 23, 2012 04:00 |
Grid Independence of Boundary layer | Luigi_ | Main CFD Forum | 0 | December 14, 2011 13:42 |
[ANSYS Meshing] Creating Hex Boundary Layer Mesh and Unstructured everywhere else | NLao | ANSYS Meshing & Geometry | 1 | March 15, 2011 06:35 |
errors | Fahad | Main CFD Forum | 0 | March 23, 2004 13:20 |