[Getdp] Issue with constraint from pre-resolution

Ruth Vazquez Sabariego ruth.sabariego at kuleuven.be
Thu Feb 8 15:48:53 CET 2018


Dear Guillaume,

You need to include the ‘Infinity’ surface in the support of your basis function:

{ Name Hrot_Hphi_3D ; Type Form1 ;
    BasisFunction {
      { Name sn ; NameOfCoef phi ; Function BF_GradNode ;
        Support Region[{Domain, Infinity}] ; Entity NodesOf[ DomainCC ] ; }
      { Name se ; NameOfCoef he ; Function BF_Edge ;
        Support DomainC ; Entity EdgesOf[ All , Not BndDomainC ]; }
    }
    Constraint {
      { NameOfCoef phi ;
        EntityType NodesOf ;
        NameOfConstraint ExtField ; }
    }
  }

Otherwise, it does not find the degrees of freedom as the research is done starting with that Support region.

If I do that I get
Info    : System 1/1: 60 Dofs
Info    : Generate[SystemBoundaryConditions]
Info    : Solve[SystemBoundaryConditions]
Info    : N: 60 - preonly lu mumps
Info    : TransferSolution[SystemBoundaryConditions]
E n d   P r e - R e s o l u t i o n  (1/1)
Info    : (Wall = 1.29631s, CPU = 0.027512s, Mem = 9.85547Mb)
E n d   P r e - P r o c e s s i n g
P r o c e s s i n g . . .
Info    : CreateDir[Results]
Info    : InitSolution[SystemMagnetics]
Info    : SaveSolution[SystemMagnetics]
Info    : TimeLoopTheta ...
Info    : Theta Time = 1e-06 s (TimeStep 1, DTime 1e-06)
Info    : GenerateJac[SystemMagnetics]
Info    : SolveJac[SystemMagnetics]
Info    : N: 1348 - preonly lu mumps
Info    :   0 Nonlinear Residual norm 0.000000000000e+00
Info    : GenerateJac[SystemMagnetics]
Info    : GetResidual[SystemMagnetics]

that does not converge, but at least the BC are there…

HTH,
Ruth




—
Prof. Ruth V. Sabariego
KU Leuven
Dept. Electrical Engineering ESAT/Electa, EnergyVille
http://www.esat.kuleuven.be/electa
http://www.energyville.be

Free software: http://gmsh.info | http://getdp.info | http://onelab.info







On 8 Feb 2018, at 15:23, DILASSER Guillaume <Guillaume.DILASSER at cea.fr<mailto:Guillaume.DILASSER at cea.fr>> wrote:

Dear GetDP user,

I am writing to call for some help concerning the constraints of type AssignFromResolution.

I am using an H-Phi magnetodynamic model and I would like to control the field on the outer boundaries of the modelling domain to match a specific time-dependant profile. Because the external field I want to apply is too complex for the usual Assign type of constraint (it is the sum of two functions with different time dependences), I have been trying to use AssignFromResolution. My current model (not working) is attached with this email. The issue is, during pre-processing, I get the following warning about the system for the pre-resolution, with is apparently empty :
Info    : Generate[SystemBoundaryConditions]
Warning : Generated system is of dimension zero
Info    : Solve[SystemBoundaryConditions]
Warning : Zero-size system: skipping solve!
Unfortunately, I do not understand where those errors come from... I have tried to identify the problem by comparing my code with the benchmark examples antennas/dipole.pro and antennas/mstrip.pro that also use AssignFromResolution constraints but I can’t seem to find it… Thank you in advance for your help !

Yours faithfully,

Guillaume DILASSER
Ingénieur de recherche / Research engineer

CEA Saclay
Bât.123 – P. 213c
Point courrier n°11
91191 Gif-sur-Yvette Cedex
France

guillaume.dilasser at cea.fr<mailto:guillaume.dilasser at cea.fr>

<Vortex_shaking_Fazilleau_1_display.geo><Vortex_shaking_Fazilleau_1_geometry.geo><Vortex_shaking_Fazilleau_1_model.pro><Vortex_shaking_Fazilleau_1_parameters.geo>_______________________________________________
getdp mailing list
getdp at onelab.info<mailto:getdp at onelab.info>
http://onelab.info/mailman/listinfo/getdp

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://onelab.info/pipermail/getdp/attachments/20180208/fa3c1a98/attachment-0001.html>


More information about the getdp mailing list