[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 -----