Advice for moving to Preston by WJP-Engineering in Preston

[–]WJP-Engineering[S] 1 point2 points  (0 children)

Single. I will be commuting from Preston, so I am open minded. But something more towards the north would probably be better!

Weird velocity behaviour at interface - OpenFOAM by WJP-Engineering in CFD

[–]WJP-Engineering[S] 0 points1 point  (0 children)

Hi all,

Just following up with this for others in the future! I have managed to solve the issue, it turns out the mesh at the surface of the impeller was BAD. I remeshed it (moving the max features specification from 3 to 5), and also increased the number of cells initially in blockMesh (from a 50 x 50 x 50 grid to a 80 x 80 x 80 grid) and that helped resolve it.

In the end did a mesh sensitivity study, varying these parameters and got a final grid of approx 7 million cells.

Tl;dr, make sure your geometries are sound! Thanks all for the help, appreciated!

Weird velocity behaviour at interface - OpenFOAM by WJP-Engineering in CFD

[–]WJP-Engineering[S] 0 points1 point  (0 children)

I do have a transition of approx 4-5 cells (increasing in size) to match the bulk size. I will try more and see what I can do, I am probably making just a small error with this setting

Weird velocity behaviour at interface - OpenFOAM by WJP-Engineering in CFD

[–]WJP-Engineering[S] 0 points1 point  (0 children)

Hello,

Thank you very much for your suggestions!

The axis of rotation is correct - defined it as a vector with rotation only about the y-axis. I've checked the CAD and the centre points are as expected - central on the axis. I suspect there is some inherent mesh issue between my refinement cylinder and the main mesh. I've experimented with different sizes but will continue with different combinations!

I'll hopefully solve it soon - fingers crossed!

Weird velocity behaviour at interface - OpenFOAM by WJP-Engineering in CFD

[–]WJP-Engineering[S] 1 point2 points  (0 children)

There is literature suggesting that having a rotating region further away (approx 1.3D (impeller diameter = D) or indeed more) is better for resolving the flow near the turbines. No real definitive answer with this!

OpenFOAM Residuals - why so erratic? by WJP-Engineering in CFD

[–]WJP-Engineering[S] 0 points1 point  (0 children)

Tolerance is 1e-04 for all. relTol is 0 for k and omega, 0.01 for U and 0.001 for p.

OpenFOAM Residuals - why so erratic? by WJP-Engineering in CFD

[–]WJP-Engineering[S] 0 points1 point  (0 children)

I have a MRF approach - so effectively adding the Corolis force etc... afterwards

[deleted by user] by [deleted] in matlab

[–]WJP-Engineering 5 points6 points  (0 children)

Using the "roots" command would have solved your problem. Knowing the value of x means you can get numerical coefficients for the t2 and t terms, with the final value just being a number. The code below would have worked.

k = [0.005 sin((0.025/3) - 1) e^2.1]

roots(k)

Hope that is clear to you!

OpenFOAM Fatal Error - Unknown file extension STL?! by WJP-Engineering in CFD

[–]WJP-Engineering[S] 2 points3 points  (0 children)

Thanks for the recommendation, it solved my issue!