Advertise on FreeFEM application

[Denis-CRETE]

Is your feature request related to a problem? Please describe.
After submission of a paper, the publisher rejected my publication on the basis that FreeFEM was not applicable to superconducting circuits only when "dimensions >> coherence length and London penetration depth". I strongly disagree with this statement and I don't know where it comes from.

Describe the solution you'd like
I suggest to add in the examples of the documentation (and manual) a small example (e.g. provided at this end of this message).

Describe alternatives you've considered
Emphasize the application to superconducting circuits by giving references to published papers (e.g. Mikhail Khapaev, M. Yu. Kupriyanov, Sergey Bakurskiy, N. V. Klenov and I.I. Soloviev, "Modeling Superconductor SFN-Structures Using the Finite Element Method" in Differential Equations 56(7):959-967 (July 2020) DOI:10.1134/S0012266120070149).

`// Far field vector potential identical if created by infinite length (along z) filament:

// A=Az0*log(r/r0);
// The superconductor is a single strip with rounded edges (line 26)
real W=25,large=50,loin=500,haut=300, e=2; // linear dimensions (in µm)
real l2= 0.09; // Lambda2 =0.3² is the square of the London penetration depth
bool debug=true;

border bottom(t=-loin,loin) { x=t; y=-haut; label=1; }
border coted(t=-haut,haut) { x=loin; y=t; label=1; }
border coteg(t=haut,-haut) { x=-loin; y=t; label=1; }
border dessus(t=loin,-loin) { x=t; y=haut; label=1; }

border bo(t=-5*W,5*W) { x=t; y=-5*W; label=2; }
border cd(t=-5*W,5*W) { x=5*W; y=t; label=2; }
border cg(t=5*W,-5*W) { x=-5*W; y=t; label=2; }
border de(t=5*W,-5*W) { x=t; y=5*W; label=2; }

mesh Th1=buildmesh(bottom(26)+coted(16)+dessus(26)+coteg(16)+bo(-201)+cd(-121)+de(-201)+cg(-121));
mesh Th2=buildmesh(bo(201)+cd(121)+de(201)+cg(121));
mesh Th=Th1+Th2;
plot(Th,ps="flux_focus_fullstrip.eps",bw=1, wait=true);

fespace Vh(Th,P2);
Vh Aphi, v;
func f=-1/l2*((abs(x)<W-e)*(abs(y)<e)+(x>=W-e)*((x-W+e)^2+y^2<e)+(x<=W-e)*((x+W-e)^2+y^2<e)); // in the superconductor only

problem poissonMagneto(Aphi,v) =
int2d(Th)( dx(Aphi)*dx(v) + dy(Aphi)*dy(v) )
+ int2d(Th) ( -v*Aphi*f ) // the source term: Mu0*J=-Aphi/l2 in the superconductor only
+ on(1, Aphi=log(x*x +y*y)); // far away
poissonMagneto;
plot(Aphi, ps="vector_potential_fullstrip.eps", value=true, wait=debug);
`