<div dir="ltr"><div><div><div><div><div><div>Dear List Members,<br><br></div>I'm just wondering if it's possible to scale the current density from the time iteration loop. <br><br>Here I introduced a parameter <b>current_factor</b>, <br><br><span style="font-family:monospace,monospace">Function{ <br> ......<br> <b>current_factor</b> = 1.;
<br>}
</span><br><br></div>which appears in the formulation as a scaling factor for the current density:<br><br><span style="font-family:monospace,monospace">Formulation{ // OK
<br> {Name VectorPotentialFormulation; Type FemEquation;
<br> Quantity{
<br> {Name Aphi; Type Local; NameOfSpace Afs;}
<br> {Name js; Type Local; NameOfSpace jfs;}
<br> }
<br> Equation{
<br> Galerkin{ [1/mu0*Dof{Curl Aphi},{Curl Aphi}];
<br> In Omega; Jacobian JVol; Integration I1;}
<br> <br> Galerkin{ DtDof[sigma[]*Dof{Aphi},{Aphi}];
<br> In Omega; Jacobian JVol; Integration I1;}
<br> <br> Galerkin{ [-<b>current_factor</b>*Dof{js},{Aphi}];
<br> In Omega_coil; Jacobian JVol; Integration I1;}
<br> }
<br> }
<br>}
</span><br><br></div>and I'd like to update this parameter during the time iteration:<br><br><span style="font-family:monospace,monospace">Resolution{
<br> {Name Solution;
<br> System{
<br> {Name Syst; NameOfFormulation VectorPotentialFormulation;}
<br> }
<br> Operation{
<br> InitSolution[Syst];SaveSolution[Syst];
<br> TimeLoopTheta[t0,tmax,dt,theta] {<br> <b>current_factor</b> = 1/1.7;
// This has no effect!!<br> Generate[Syst]; Solve[Syst]; SaveSolution[Syst];
<br> }
<br> }
<br> }
<br>}</span><br><br></div>but this has no impact on the solution. Apparently the value of current_vector is still considered one as it's defined in Functions. <br><br>Is there any way to update this parameter from the time loop? <br><br></div>Thanks for your help in advance.<br><br><br></div>Peter<br><div><div><br></div></div></div>