[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