[Getdp] ... Formulation question ...
mkoch at gvtc.com
mkoch at gvtc.com
Wed Feb 28 23:24:58 CET 2007
Hello Olivier,
and thanks for the quick response. I was afraid you were going to say
something like you did. Is there at least a way to define Norm(Grad
Phi) as a quantity somewhere, say call it NrmGrdPhi, (perhaps in
Quantity under Formulation), and then use it in the Galerkin term. So
something akin to this:
Formulation{
Quantity{NrmGrdPhi; [Norm[Grad Phi]];}
Galerkin{[Spd[]*Dof{NrmGrdPhi},{Phi}];In Vlm;Integration Int;}
}
I suppose this would be cheating, too, but I got to give it a try. Any
thoughts?
Thanks,
Matt Koch
----- Message from castany at quatramaran.ens.fr ---------
Date: Wed, 28 Feb 2007 21:47:53 +0100
From: Olivier Castany <castany at quatramaran.ens.fr>
Reply-To: Olivier Castany <castany at quatramaran.ens.fr>
Subject: Re: [Getdp] ... Formulation question ...
To: getdp at geuz.org
>> I am trying to formulate "s*Norm(Grad Phi)" vs. "v*Grad(Phi)", whereby
>
>> However, the formulation with the speed crashes with a parse error on "[":
>>
>> Galerkin{[Spd[]*Dof{Norm[{d Phi}]},{Phi}];In Vlm;Integration
>> Int;Jacobian Jac;}
>
> My 2 cents :
>
> "Norm(Grad Phi)" is not linear with respect to the field Phi. It is thus
> not possible to build a matrix equation with that term.
>
> You can use something like : "Spd[]*{Norm[{d Phi}]}" in a "Galerkin"
> term. But the "Dof{...}" can not be used here (it would ask to build a
> matrix equation, which is not possible here).
>
> _______________________________________________
> getdp mailing list
> getdp at geuz.org
> http://www.geuz.org/mailman/listinfo/getdp
>
----- End message from castany at quatramaran.ens.fr -----