BH curve

More
1 year 8 months ago #4620 by lahoussaine
BH curve was created by lahoussaine
Good morning community;
is there a way to implment bh curve using ngsolve : 
i do somthing like 
def nuf(u):
     return Norm(grad(u))
nu=CoefficientFunction(nuf)
a+=nu*grad(u)*grad(v)*dx
but always i got nan values !
any tricks ? 
More
1 year 8 months ago #4621 by mrambausek
Replied by mrambausek on topic BH curve
Dear Lahoussaine,

what is the actual model you want to implement?
"Norm(grad(u))" does not seem to a suitable model in first place, as it leads to a singular matrix for grad(u) = 0.

Best,
Matthias
More
1 year 8 months ago #4625 by lahoussaine
Replied by lahoussaine on topic BH curve
I want something like
1/(B**3+B**2+B )
or :
3.8*exp(2.17*B*B)+396.2

I ve tried those expression but it gives always NAN;

Thankes for your reply
More
1 year 8 months ago #4626 by mrambausek
Replied by mrambausek on topic BH curve
TBH, I have never seen any magnetic model like that (which might not mean a lot...)

I can't see how the first one could give something reasonable for B=0.

Before using a nonlinear model, did you have a working version for a constant "nu"?

I fear, we need much more information to be able to help you.


Best
More
1 year 8 months ago #4627 by lahoussaine
Replied by lahoussaine on topic BH curve
def nuf(u):
IfPos(2.28-B,(1/(-0.003995*B**3+0.005956*B**2+0.00646*B + 0.001925)), ((1902.211*B-4266.42343)/0.016694669))

Here is my formule, but it doesn't work
The problem is whenever i put norm(B) it returns NAN
even if the formula is well defined( 1/(B**2+1) for example )


Thank you and Sorry to not be claire
More
1 year 8 months ago #4628 by mrambausek
Replied by mrambausek on topic BH curve
Unfortunately, I still miss some information: how is "B" defined?
Where does the "Norm" enter in "1/(B**2 + 1)
Which command did you use for assembly? Where exactly do you get NaN? Already in the matrix or only after solving, i.e. in the solution?

Could you maybe post a complete minimal example that reproduces your problem?
Time to create page: 0.097 seconds