CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Generating grids with pointwise and reproducing propeller rotation calculations with

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 15, 2021, 05:37
Post Generating grids with pointwise and reproducing propeller rotation calculations with
  #1
New Member
 
hirota tomohiro
Join Date: Jun 2021
Posts: 16
Rep Power: 5
hirota is on a distinguished road
Hi Foamers,

Currently, I am trying to reproduce the rotation of a propeller using pimpleDyMFoam.
Specifically, I want to simulate a rotating region (propeller rotation region) inside a stationary region.

First, I assume that I will be using pointwise to generate the grid.
Next, I am considering using the AMI function in consideration of the case where propeller rotation is derived in the future.

As shown in the Fig1, the computational domain is a hexahedral region with a cylinder rotation region inside.

This analysis was output from pointwise as polyMesh to openFoam specification, and after that I had a hard time setting the AMI boundary, and tried to find a solution, but still could not reach a solution.

The procedure is
1, Create a rotated region and a stationary region by pointwise, and output them by polyMesh.
 The boundry files at this time are shown in Figure 2 and Figure 3.
2, Perform renumberMesh.
3, Create a cellzone as a rotation area by topoSet (topoSetDict is shown in Figure 4)
4, MergeMesh to combine two meshes.
5, Use changeDictionary to change the type of AMI1 and AMI2 of the boundary from patch to cyclicAMI (changeDictionaryDict is attached in Figure 5)
6, Describe the name of the cellZone that will be the rotation region in the dynamicMeshDict, and execute the simpleDyMFoam.

If you follow these steps, you will see that the calculation turns, but the velocity and pressure distributions do not change even though you are just rotating an empty mesh inside a rectangular body.
However, the velocity distribution fluctuates and the flow is disturbed. (Figure 6; the circle near the center is the cross section of the cylinder in the rotating region)
After investigating the cause, we found that when step 5 was executed in changeDictionary, an impossible velocity distribution was given to the cellZone, the rotating region, even though it was the initial condition (Fig. 6) Since it is the boundary of AMI, it should be zero.
The pressure and other distributions are still zero, but for some reason only the velocity distribution is like this, and this seems to be the cause of the flow disturbance during the calculation.
I am at a loss and am having a very difficult time. I would be grateful if someone could help me.
Attached Images
File Type: png Fig1.png (47.4 KB, 5 views)
File Type: png Fig2_rotor_Arew_boundry.png (15.9 KB, 2 views)
File Type: png Fig3_steady_Area_boundry.png (46.9 KB, 2 views)
File Type: png Fig_4_topoSetDict.png (53.3 KB, 2 views)
File Type: png Fig5_changeDictionaryDict.png (65.9 KB, 2 views)
hirota is offline   Reply With Quote

Old   July 15, 2021, 05:40
Default Here are some additional images. Fig.6 &Fig.7
  #2
New Member
 
hirota tomohiro
Join Date: Jun 2021
Posts: 16
Rep Power: 5
hirota is on a distinguished road
Here are some additional images.
Fig.6 &Fig.7
hirota is offline   Reply With Quote

Reply


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



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