mat_Fe = mag.MagMaterial("Fe",
Ms=1.2,
A=4.0,
extra_H="""H_total_Fe(0) -= mag[0]*1.2;
if(X(0)<0.01 || X(0)>2.99) /* pin the M field at boundaries */
{
H_total_Fe(0)=0.0;
H_total_Fe(1)=0.0;
H_total_Fe(2)=0.0;
}
""")
mag.set_local_magnetic_coupling(mat_Dy, mat_Fe, -80.0)
mag.defregion("DyFe2", nm.box([1.0], [2.0]), mag_mat=[mat_Dy, mat_Fe])
mag.defregion(
"YFe2",
nm.union([nm.box([0.0], [1.0]), nm.box([2.0], [3.0])]),
mag_mat=mat_Fe # can pass either a material or a list of materials
)
mag.set_meshing_parameters(
cache_name="exchange-spring-mesh",
bounding_box=([-1.0], [5.0]),
a0=0.25,
max_steps=100,
)
mag.create_mesh()
mat_Fe = mag.MagMaterial("Fe",Ms=1.2,A=4.0,
extra_H="""H_total_Fe(0) -= mag[0]*1.2;
if(X(0)<0.01 || X(0)>2.99) /* pin the M field at boundaries */
{
H_total_Fe(0)=0.0;
H_total_Fe(1)=0.0;
H_total_Fe(2)=0.0;
}
"""
)
mag.set_local_magnetic_coupling(mat_Dy,mat_Fe,-80.0)
mag.defregion("DyFe2",nm.box([1.0],[2.0]),
mag_mat=[mat_Dy,mat_Fe]
)
mag.defregion("YFe2", nm.union([nm.box([0.0],[1.0]),nm.box([2.0],[3.0])]),
mag_mat=mat_Fe # can pass either a material or a list of materials
)
mag.set_meshing_parameters(cache_name="exchange-spring-mesh",
bounding_box=([-1.0],[5.0]),
a0=0.25,
max_steps=100,
)
mag.create_mesh()
#mag.set_intensive_parameters(["H_x","H_y","H_z"], # may also set T,p,E, etc.
# external_magnetization=["H_x","H_y","H_z"])
# ^ XXX not used yet, as we are going to use the experimental cvode integrator...
mag.set_intensive_parameters([])
# mag.set_default_order(1) # default anyway...
mag.set_features({"demag":True,"exchange":True,"timestep":True}) # XXX change usage!
mat_Py = mag.MagMaterial("Py",Ms=1.0,A=13.0, # A = exchange constant
)
mag.defregion("Py", nm.ellipsoid([3.0,3.0,3.0]), mag_mat=mat_Py)
mag.set_meshing_parameters(cache_name="exchange-spring-mesh",
bounding_box=([-4.0,-4.0,-4.0],[4.0,4.0,4.0]),
a0=1.0,
max_steps=600,
)
mag.create_mesh()
mag.set_magnetization([0.0,0.0,1.0]) # providing a vector rather than a function
import nfem.visual
nfem.visual.fields2vtkfile([mag.default_simulation_context.field_M],'sphere-initial.vtk',mesh=mag.default_simulation_context.mesh)
h_total_Funny[2] += 0.0;
""")
mag.set_default_material(PermAlloy)
mag.set_intensive_parameters(["T", "p", "H_x", "H_y", "H_z"])
#mag.defregion("Ball 1",nm.ellipsoid([3.0,3.0,3.0],transform=[("shift",[-3.0,0.0,0.0])]))
print "OK 1"
sys.stdout.flush()
#mag.defregion("Ball 2",nm.ellipsoid([3.0,3.0,3.0],transform=[("shift",[3.0,0.0,0.0])]))
mag.defregion(
"Ball 2",
nm.ellipsoid([2.0, 2.0, 2.0], transform=[("shift", [3.0, 0.0, 0.0])]))
# Note: clearly, we DO need a better way to specify geometries. Ideally, I would like to be
# able to write instead:
#
# mag.defregion("Ball 1",nm.shifted([-3,0,0],nm.sphere(3)))
# mag.defregion("Ball 2",nm.shifted([ 3,0,0],nm.sphere(3)))
#
# or alternatively:
#
# sphere = nm.sphere(3)
# mag.defregion("Ball 1",nm.shifted([-3,0,0],sphere))
# mag.defregion("Ball 2",nm.shifted([ 3,0,0],sphere))
mag.set_meshing_parameters(cache_name="two-balls")
h_total_Funny[2] += 0.0;
""")
mag.set_default_material(PermAlloy)
mag.set_intensive_parameters(["T","p","H_x","H_y","H_z"])
#mag.defregion("Ball 1",nm.ellipsoid([3.0,3.0,3.0],transform=[("shift",[-3.0,0.0,0.0])]))
print "OK 1"
sys.stdout.flush()
#mag.defregion("Ball 2",nm.ellipsoid([3.0,3.0,3.0],transform=[("shift",[3.0,0.0,0.0])]))
mag.defregion("Ball 2",nm.ellipsoid([2.0,2.0,2.0],transform=[("shift",[3.0,0.0,0.0])]))
# Note: clearly, we DO need a better way to specify geometries. Ideally, I would like to be
# able to write instead:
#
# mag.defregion("Ball 1",nm.shifted([-3,0,0],nm.sphere(3)))
# mag.defregion("Ball 2",nm.shifted([ 3,0,0],nm.sphere(3)))
#
# or alternatively:
#
# sphere = nm.sphere(3)
# mag.defregion("Ball 1",nm.shifted([-3,0,0],sphere))
# mag.defregion("Ball 2",nm.shifted([ 3,0,0],sphere))
mag.set_meshing_parameters(cache_name="two-balls")
