Getting started with CyclicAMI
Hello all,
I'm trying to model the flow interaction between two connected channels, and I've set up the following geometry for the left channel with blockMesh: https://s2.postimg.org/m0wwik5pl/Single_Duct_1.png Where the light blue area represents a coarse mesh, and the darker blue a fine mesh. I'd like to use cyclicAMI to connect the rough and fine areas of the mesh, as it seems impossible to merge the mesh and cyclicAMI is the next best thing. Unfortunatly, I'm quite a novice at using OpenFOAM and I've been unable to find a good example of cyclicAMI being used with blockMesh. This is my blockMeshDict file: Code:
/*--------------------------------*- C++ -*----------------------------------*\ Quote:
Am I even correct in including the cyclicAMI boundary conditions in the blockMeshDict file? |
It seems you have spelt 'neighbour' wrong .. although I'm not sure if there are other more appropriate boundary conditions. Have you had a look at snappyHexMesh ? You will be able to refine the dark blue region without splitting up your blockMesh.
|
Quote:
I've boon considering snappyHexMesh, but my supervisor would like to see if it's possible to do in just blockMesh first. I've got it working below, but there are anomalies around the cyclicAMI patches: https://s1.postimg.org/4k8hruji7/cavity18_velocity.png At first I assumed they might just be the due to some oddity in how it dealt with the initial conditions, but running it for increasing time-steps has shown that while the sudden jump in pressure and velocity sometimes goes away, there are always some discontinuities present along the boundary. Irrespective of whether I should just use snappyHexMesh (and I would welcome any advice on the pros and cons of it compared to the current method), I'm curious as to what causes the discontinuities seen above. Does anyone have an idea? |
The advantage of SHM in this case would be that you don't have to worry about specifying any interior boundary conditions - such as the cyclic part. Also, the mesh that SHM will generate may lead to a more accurate solution since the cell nodes will be shared rather than having to interpolate between them.
I think the image you have printed shows a feature of the paraview interpolation rather than the cell values. If you change the variable to its cell values rather than the interpolated point fields the issues may disappear. |
Edit boundary file so its similar to mine. Obviously without the exact numbers and use whatever patch names you have.
AMI2 { type cyclicAMI; nFaces 9878; startFace 4084895; neighbourPatch AMI1; } AMI1 { type cyclicAMI; nFaces 10174; startFace 4094773; neighbourPatch AMI2; *I know this thread is old but for anyone else who finds this thread. |
All times are GMT -4. The time now is 19:51. |