Re: [Gfs-users] singular matrix problem !
Brought to you by:
popinet
Menu
▾
▴
Re: [Gfs-users] singular matrix problem !
|
From: sylv s. <syl...@ho...> - 2008-09-27 11:24:25
|
NOBODY ?
>
> hello !
> i have got a problem with a warning "singular matrix".
> I suppose i should have done something wrong in the file !
>
> Could someone chekh my new user error ?
>
>
>
>
>
>
> # gerris3D bubble.gfs 1>out 2>err
> #Wombbling bubble in a pipe
> #moving walls at Vgliss=-0.2
> #mu_l=4.e-3 rho_l=997.48
> #mu_g=18.57e-6 rho_g=1.161
> #surface tension=0.0717
>
>
> 4 3 GfsSimulation GfsBox GfsGEdge {} {
> Time { end = 0.5 dtmax=1e-5}
>
> Global {
> #define var(T,min,max) (CLAMP(T,0,1)*(max - min) + min)
> #define rho(T) var(T, 1.161,997.48)
> #define mu(T) var(T, 18.57e-6,4.e-3)
> #define level 5
> #define radius 0.0045
> #define surface_tension 0.0717
> #define Vgliss -0.2
> #define init_pos 0.0775
> }
>
> Refine 5
>
> ProjectionParams { }
> ApproxProjectionParams { }
>
> # We also need gravity
> Source {} V -9.81
>
> VariableTracerVOF T
> InitFraction T ( (x*x+(y-init_pos)*(y-init_pos)+z*z-radius*radius))
> Init {} { U = 0 V=Vgliss W=0}
>
> PhysicalParams { alpha = 1./rho(T) L=0.025}
> SourceViscosity mu(T)
> VariableCurvature K T
> SourceTension T surface_tension K
>
> #MESH REFINEMENT
> AdaptVorticity { istep = 1 } { cmax = 0.05 maxlevel = level minlevel = 3}
> AdaptGradient { istep = 1 } { cmax = 1e-3 maxlevel = level minlevel = 3} T
>
> #outputs
> GfsOutputBalance { istep = 1 } stdout
> GfsOutputTime { istep = 1 } stdout
>
> OutputScalarSum { istep = 1 } Positionx { v = x*(1-T) }
> OutputScalarSum { istep = 1 } Positiony { v = y*(1-T) }
> OutputScalarSum { istep = 1 } Positionz { v = z* (1-T) }
> OutputScalarSum { istep = 1 } Volume { v = (1-T) }
>
> GfsOutputSimulation { step = 0.01 } case_02-%g.gfs { }
> }
> #MESH
> GfsBox { bottom = GfsBoundaryOutflow
> left = GfsBoundary { GfsBcDirichlet U 0 GfsBcDirichlet V Vgliss GfsBcDirichlet W 0}
> right = GfsBoundary { GfsBcDirichlet U 0 GfsBcDirichlet V Vgliss GfsBcDirichlet W 0 }
> front = GfsBoundary { GfsBcDirichlet U 0 GfsBcDirichlet V Vgliss GfsBcDirichlet W 0}
> back = GfsBoundary { GfsBcDirichlet U 0 GfsBcDirichlet V Vgliss GfsBcDirichlet W 0 } }
>
> GfsBox { left = GfsBoundary { GfsBcDirichlet U 0 GfsBcDirichlet V Vgliss GfsBcDirichlet W 0 }
> right = GfsBoundary { GfsBcDirichlet U 0 GfsBcDirichlet V Vgliss GfsBcDirichlet W 0 }
> front = GfsBoundary { GfsBcDirichlet U 0 GfsBcDirichlet V Vgliss GfsBcDirichlet W 0}
> back = GfsBoundary { GfsBcDirichlet U 0 GfsBcDirichlet V Vgliss GfsBcDirichlet W 0 } }
>
> GfsBox { left = GfsBoundary { GfsBcDirichlet U 0 GfsBcDirichlet V Vgliss GfsBcDirichlet W 0 }
> right = GfsBoundary { GfsBcDirichlet U 0 GfsBcDirichlet V Vgliss GfsBcDirichlet W 0 }
> front = GfsBoundary { GfsBcDirichlet U 0 GfsBcDirichlet V Vgliss GfsBcDirichlet W 0}
> back = GfsBoundary { GfsBcDirichlet U 0 GfsBcDirichlet V Vgliss GfsBcDirichlet W 0 } }
>
> GfsBox { top = GfsBoundary { GfsBcDirichlet U 0 GfsBcDirichlet V Vgliss GfsBcDirichlet W 0 }
> left = GfsBoundary { GfsBcDirichlet U 0 GfsBcDirichlet V Vgliss GfsBcDirichlet W 0 }
> right = GfsBoundary { GfsBcDirichlet U 0 GfsBcDirichlet V Vgliss GfsBcDirichlet W 0 }
> front = GfsBoundary { GfsBcDirichlet U 0 GfsBcDirichlet V Vgliss GfsBcDirichlet W 0}
> back = GfsBoundary { GfsBcDirichlet U 0 GfsBcDirichlet V Vgliss GfsBcDirichlet W 0 } }
> 1 2 top
> 2 3 top
> 3 4 top
>
>
> thank you
> sylv1
> _________________________________________________________________
> Téléphonez gratuitement à tous vos proches avec Windows Live Messenger ! Téléchargez-le maintenant !
> http://www.windowslive.fr/messenger/1.asp
> -------------------------------------------------------------------------
> This SF.Net email is sponsored by the Moblin Your Move Developer's challenge
> Build the coolest Linux based applications with Moblin SDK & win great prizes
> Grand prize is a trip for two to an Open Source event anywhere in the world
> http://moblin-contest.org/redirect.php?banner_id=100&url=/
> _______________________________________________
> Gfs-users mailing list
> Gfs...@li...
> https://lists.sourceforge.net/lists/listinfo/gfs-users
_________________________________________________________________
Installez gratuitement les 20 émôticones Windows Live Messenger les plus fous ! Cliquez ici !
http://www.emoticones-messenger.fr/ |