Trying to integrate waves2Foam with overInterDyMFoam
Hello all,
I am currently trying to integrate waves2Foam with the overInterDyMFoam solver in OpenFOAM-v1712. I have previously achieved to do the same with interDyMFoam, so I thought I could give this also a try. My solver (overWaveDyMFoam) compiles successully, but when I run a simple case with a floating box, I get the following error in the very first time step: Code:
Starting time loop Thank you! Pascal. |
Hi Pascal,
Interesting that you try waves2Foam with the overset functionality. It has long been on my to-do list. Based on the error, I guess that you have forgotten to modify createFields.H or include the required header files in the *.C file. An alternative explanation, on second thought, is that the overset mesh functionality modifies some of the member variables in the boundary conditions. If this deletes all pointers, it means that you get the missing autoPtrI-error. One possible work-around is the substitute waveAlpha with zeroGradient and waveVelocity with zeroGradient and rely solely on the relaxationZones. If you make it work, I would very much like you to share some first experiences here in this thread. Kind regards, Niels |
First experiences with waves2Foam and overset grids
2 Attachment(s)
Hi Niels,
I am really happy to have you in this discussion! I can confirm that things start to run correctly when I substitute waveAlpha and waveVelocity with zeroGradient types. I run the simulation for a few seconds and it seems to behave well... I attach my little test case and my overWaveDyM solver as it has evolved so far. I do not know what else should be modified in the createFields.H file, so any hint is appreciated. It would be great if you had the chance to check the solver and the test case (to eliminate a possibly wrong setup). I tried to track the error further down and found, that things run until psi.correctBoundaryConditions(); is called in the MULES correction in file OpenFOAM-v1712/src/finiteVolume/fvMatrices/solvers/MULES/CMULESTemplates.C Beyond this point I am not able to follow the code any further... unfortunately... @InterestedFolks: The attached solver files actually contain my waveDyMFoam solver as well, so everything should go to waves2Foam/applications/solvers/solvers1712_PLUS I am looking forward to reading your comments! Best regards, Pascal. |
Interesting. It seems that something is not working in the boundary conditions.
I will try to get time to look at the bug; however, no promises at present. Kind regards, Niels |
First experiences with waves2Foam and overset grids - 2
Hi Niels,
I am looking forward to your findings! With regard to your request for sharing first experiences - have a short laugh with the attached animation of the case I posted before. It does not behave well at all... :-/ I suppose I have some boundary condition settings wrong... https://drive.google.com/file/d/1w1H...ew?usp=sharing (hope this works) Anyway, next week I will be out of office, so I will not be able to respond to postings... just in case. Best regards, Pascal. |
Hi Pascal,
That is actually pretty funny! My first guess is that your acceleration relaxation is too large or you are using too few pimple-type corrections per time step. Are you initializing with setWaveFields? If yes, does it seem to recognize the multiple domains? Kind regards, Niels |
1 Attachment(s)
Hi Niels,
back again... As often, it turned out that the error sits in front of the screen... ;-) On close inspection of the boundary conditions, I found some mistakes. Now things run smoothly (at least for the first 10 seconds) and the floating body does not seem to be jumping out of the domain again. Anyway, the problem with the boundary condition for alpha.water at the inlet still persists. Things only work when I set the inlet condition to: Code:
type zeroGradient; Here is an animation of the corrected case: https://drive.google.com/open?id=1KG...ihI3K5DEbykeGk I am looking forward to hearing from you about the boundary condition issue! Thank you and best wishes, Pascal. |
Hi Pascal and Niels,
I am facing the same bug now. Do you find it out? I can run correctly when I substitute waveAlpha and waveVelocity with zeroGradient types, as you said. However, I am just wondering that does this subsitute that rely solely on the relaxationZones affect the finally wave generation? Because, I think, zerogradient at Inlet Boundary condition cannot gives the correctly wave elevation as waveAlpha BC. Maybe I am worry. best, Shiyu:) |
Possible (little) bug found in waves2Foam
Hello Niels,
I think I might have stumbled across a weird behaviour of the relaxationZoneLayout tool in waves2Foam. When I use it in a non decomposed case, I get the output: Code:
... snippet ... Now, when I run relaxationZoneLayout on the same case, but decomposed and so in parallel, I get: Code:
... snippet ... At the time being I think this does not affect the actual simulation, but it is confusing. Am I doing something wrong here or might the tool be flawed? I am using a version downloaded on 2018-06-26. Thank you and best regards, Pascal. -- |
Hallo,
You have now posted the same question in multiple threads, but you provide no information. Please refrain from cross-posting and please inform us, what it is you try to accomplish. You state that waves2Foam is compiled, so I can only assume that you try to make a custom solver. Kind regards Niels |
Quote:
I have tried your test case, and it gives out the error message as below: --> FOAM FATAL IO ERROR: Dictionary entry for patch XMIN not found May I know how can I solve it or which directory should I look at, Thank you! __________________________________________________ __________________ A quick update, the problem is solved. The error came from the incomplete installation of overWaveDyMFoam. Best, Venus |
overwavedymfoam
2 Attachment(s)
Hi, did you succeed in coupling the slovers of waves2foam and overinterdymfoam? I am also trying this work and have encountered an unsolvable problem:
I finished compiling overwavedymfoam and no error was reported. But when I run a case, I find the case diverged after 35s which is close to stability. I tried various modifications but couldn't change the divergence phenomenon and the version of OpenFoam I used is OF2206. Can you give me some advice? Thanks. I will attach my case and solver below. |
Quote:
I have a question about wave2foam. Does wave2foam support RBFMeshMotionSolver? From my test, wave2foam can't use RBFMeshMotionSolver. If not, I want to integrate RBFMeshMotionSolver into wave2foam. But I have no idea about how to do this. Could you share your experience and give me any advice? |
Quote:
|
All times are GMT -4. The time now is 13:49. |