CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Community Contributions (https://www.cfd-online.com/Forums/openfoam-community-contributions/)
-   -   [waves2Foam] floating object inside the flume tutorial with interDyMFoam (https://www.cfd-online.com/Forums/openfoam-community-contributions/207971-floating-object-inside-flume-tutorial-interdymfoam.html)

Rising Star April 2, 2015 21:37

floating object inside the flume tutorial with interDyMFoam
 
hi everybody

I'm completely new to openfoam and CFD world, so please be patient if my request seems stupid

for studying purposes I've to simulate a wavetank including a floating plate. During my researches I've found that waves2foam is more or less what I need, but I don't know what I've to do to add the floating object, for example to the waveflume tutorial.

I'm running openfoam 2.3.0, and I've found the tutorial linked here (https://sites.google.com/site/jordim...ic-mesh-motion) but it's not working at all, since it's designed for OF 1.7.1

could you please help me?

tnx in advance

Best Regards

Rising Star April 13, 2015 02:21

floating object inside the flume tutorial with interDyMFoam
 
hi everybody

I think I've coupled waves2foaam with interdymfoam (i'm a newbie), so I've tried to put a floating object inside the flume tutorial, but the simulation crashes immediately:(

this is the error I get

Quote:

Centre of mass: (10.125 -0.0881379 0.05)
Linear velocity: (-0.034629 -43.6092 6.78388e-05)
Angular velocity: (-13839 8.03976 -10.6603)
Execution time for mesh.update() = 0.07 s
GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 4.64817e-06, No Iterations 8
GAMGPCG: Solving for pcorr, Initial residual = 0.0490237, Final residual = 4.62368e-06, No Iterations 4
time step continuity errors : sum local = 3.06048e-06, global = 2.01192e-07, cumulative = 1.18549e-07
smoothSolver: Solving for alpha.water, Initial residual = 0.00160084, Final residual = 6.28579e-09, No Iterations 17
Phase-1 volume fraction = 0.5637 Min(alpha1) = 0 Max(alpha1) = 1.892
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.5637 Min(alpha1) = -6.15106 Max(alpha1) = 76.4838
Relaxing time: 0.03 s
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
#4 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#8
at ??:?
#9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10
at ??:?
those are the files I'm using
https://www.dropbox.com/s/z6h5nhm8sr...me.tar.gz?dl=0


please help me, I'm stuck and I don't know how to solve this situation

thanks in advance
best regards

Olbap April 20, 2015 23:42

@Rising : I ran your case and found the same error, although changing the time scheme seems to change the error time. So I tried to output what happens before the crash, and it seems your floating object is not well defined, the wave field just crashes into it: I don't know what happens exactly but the Courant number explodes and the maximum water fraction is way above 1. It looks like and obstacle.
I have no solid advice, but have you tried either 1) your waveDymFoam on a known tutorial like waveFlume 2) to put a floating object using just waveFoam (or interFoam)? What were the results?
Best,
Pablo


All times are GMT -4. The time now is 01:53.