Cylinder Py - 5 x 10 x 20 (nm)

Cylinder Py - 5 x 20 x 40 (nm)

Cylinder Py - 5 x 100 x 200 (nm)

F-AF-F layer coupling

1.0e-8, vert / elem = 1662 / 8038

Amorphous wire plus stress

FOLLOW-UP of Part I and Part II


1.0e-5, vert / elem = 34214 / 184912, rad / len = 0.1 / 1.0

1. NO STRESS



2. R>0.06, STRESS 5.0e7 Pa (50%) 5.0e7 OK for 100% case, nu enough for 50% case


3. R>0.06, STRESS 8.0e7 Pa (50%) 8.0e7 OK for 50% case, no central domain


4. R>0.03, STRESS 5.0e7 Pa (80%) 5.0e7 OK for 80% case, with central domain

Magnetostatic interactions

Material:
K1 = 1.00E05 J/m^3
Js = 2.00 T
A = 1.00E-11


1. same length
L = 1.0e-6 m
r = 1.0e-8 m

One magnetic wire (non-magnetic neighbors)

Nine magnetic wires

Twenty five magnetic wires

2. random length
L around 1.0e-6 m
r = 1.0e-8 m

One magnetic wire (non-magnetic neighbors)

Nine magnetic wires

Twenty five magnetic wires

Interactions

Material:
K1 = 1.00E03 J/m^3
Js = 5.00 T
A = 1.00E-11

2. 9,600 elems.
L = 1.0e-7 m
r = 0.1e-7 m
H = 500 ... -500 kA/m

3. 15,000 elems.
L = 1.0e-7 m
r = 0.1e-7 m
H = 500 ... -500 kA/m

Long wire plus stress

1.0e-6 150638 / 764830 - stress 1e8

Wire plus Stress

1.0e-5 132511 / 723688 - stress 4e7




1.0e-5 132511 / 723688 - stress 4.5e7 - Z




1.0e-5 132511 / 723688 - stress 4.5e7 - Y

Wire

0.5e-6 18248 / 90461



1.0e-6 30650 / 156657



2.0e-6 30650 / 156657





3.0e-6 30650 / 156657


1.0e-5 132511 / 723688

Ring - radiuses 5.0e-7...3.0e-7 - LLG (refine=0)

Ring - radiuses 5.0e-7...3.0e-7 - TAO (refine=2)

 

GiD Procedure

1. start GiD
2. set problem type:
   "Data/Problem type/inp"
3. create object(s)
4. assign "Conditions" (=property ids=grain numbers) - the corresponding material parameters are later defined in project.krn: material properties in the appropriate line (line number=property id):
   "Data/Conditions"
5. set grain number
   "grain number 1"
6. "Assign"
7. select all objects, which should have property id ”1”
8. press the "Finish" button or the Esc key when finished
9. Repeat 5-8 if necessary
10. Just in case:
    "Geometry/Delete/Point"
    "Geometry/Delete/Line"
    "Geometry/Delete/Surface"
    only volumes should be left over!
11. generate a tetrahedral finite element mesh:
    "Meshing/Generate"
12. Enter the size of elements, 
13. save the project:
    "Files/Save" "sphere"
14. generate the inp file:
    Calculate/Calculate
15. rename the inp file:
    mv sphere.dat sphere.inp

New mkinp_vect

#include "stdlib.h"
#include "stdio.h"
#include "string.h"

int main(int argc, char *argv[])
{
    FILE *fp_IN, *fp_OUT;
    char nume_fis[200], old_nume_fis[200], new_nume_fis[200], delete_command[200];
    char *dummy_char_name, *dummy_nr;
    int file_nr;
    int nr_elem, nr_nodes, nr_data, dummy_int;
    int i;
    int dummy;
    double M_x, M_y, M_z, Hdemag_x, Hdemag_y, Hdemag_z, Hani_x, Hani_y, Hani_z, Hexch_x, Hexch_y, Hexch_z, Hext_x, Hext_y, Hext_z;
    char etapa1_line[250];

    strcpy(nume_fis,argv[1]);
    dummy_char_name = strtok (nume_fis,"."); //nume
    dummy_nr        = strtok (NULL,"."); //numar
    file_nr         = atoi (dummy_nr);

    sprintf(old_nume_fis, "%s.%04d.inp", dummy_char_name, file_nr);
    sprintf(new_nume_fis, "%s_v.%04d.inp", dummy_char_name, file_nr);
    while( (fp_IN  = fopen( old_nume_fis, "r" )) != NULL )
    {
        printf("Working on file %04d\n", file_nr);
        fp_OUT  = fopen( new_nume_fis, "w" );
        fscanf(fp_IN, "%d %d %d %d %d\n", &nr_elem, &nr_nodes, &nr_data, &dummy_int, &dummy_int);
        fprintf(fp_OUT, "%d %d %d %d %d\n", nr_elem, nr_nodes, nr_data+3, dummy_int, dummy_int);

        for (i=0;i < nr_elem;i++)
        {
            fscanf(fp_IN, " %[^\n]", etapa1_line);
            fprintf(fp_OUT, "%s\n", etapa1_line);
        }
        for (i=0;i < nr_nodes;i++)
        {
            fscanf(fp_IN, " %[^\n]", etapa1_line);
            fprintf(fp_OUT, "%s\n", etapa1_line);
        }

        fscanf(fp_IN, " %[^\n]", etapa1_line);
        fprintf(fp_OUT, "16 3 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1\n", etapa1_line);

        fprintf(fp_OUT, "M, none\n");

        for (i=0;i < nr_data;i++)
        {
            fscanf(fp_IN, " %[^\n]", etapa1_line);
            fprintf(fp_OUT, "%s\n", etapa1_line);
        }

        for (i=0;i < nr_elem;i++)
        {
            fscanf(fp_IN, "%d %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf\n", &dummy, &M_x, &M_y, &M_z, &Hdemag_x, &Hdemag_y, &Hdemag_z, &Hani_x, &Hani_y, &Hani_z, &Hexch_x, &Hexch_y, &Hexch_z, &Hext_x, &Hext_y, &Hext_z);
            fprintf(fp_OUT, "%d %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf %lf\n",  dummy, M_x, M_y, M_z, M_x, M_y, M_z, Hdemag_x, Hdemag_y, Hdemag_z, Hani_x, Hani_y, Hani_z, Hexch_x, Hexch_y, Hexch_z, Hext_x, Hext_y, Hext_z);
        }

        fclose(fp_IN);
        fclose(fp_OUT);

        sprintf(delete_command, "rm %s", old_nume_fis);
        system(delete_command);

        file_nr++;
        sprintf(old_nume_fis, "%s.%04d.inp", dummy_char_name, file_nr);
        sprintf(new_nume_fis, "%s_v.%04d.inp", dummy_char_name, file_nr);
    }
    return(0);
}

Disk

disk.inp -> disk de raza 5, inaltime 0.5
-> min element size 0.6 -> 3532 elems / 1254 nodes
(-> min element size 0.5 -> 15485 elems / 4200 nodes)

disk.krn -> 0.0 0.0 5.0e2 0.0 1.00 1.30E-11 0.1 uni #(Permalloy)

allopt.txt 
#########################################################################
#
# magpar configuration file: allopt.txt
#
# default settings
#
#########################################################################

-simName disk

-meshtype 1          # read finite element mesh from AVS
#-inp 0232
-refine 1
-size 1.0e-7         # set length scaling factor to 10 nm
-init_mag 1         # Mx
-mode 1              # select energy minimization using TAO
-demag 1             # switch on demagnetizing field

-hextini 100.0        # apply a homogeneous external field of -40 kA/m
-htheta 1.57         # apply homogeneous external field at 0.1 rad from z-axis
-hphi   0.017453293
-hstep -0.5          # change field amplitude in steps of -0.1 kA/m
-hfinal -100.0        # final field value -40 kA/m
-hext_cu 0

-ts_logsteps 1000
-ts_nsteps 1
-condinp_equil 1
-condinp_j 0.1
-condinp_equil_j 0.0
-condinp_t 1e300

#-slice_n 0,0,1
#-slice_p 0,0,2.5e-1
#-slice_g 0        # 0:  plot magnetization of grains with any property id
#-res 200
#########################################################################
#########################################################################
# tolerances
#########################################################################
#########################################################################
-tol 1e-4

-jfinal -1e99
-mfinal -1e99,-1e99,-1e99
#########################################################################
#########################################################################
# options for solvers
#########################################################################
#########################################################################

-magdist 0.0        # magdist==0: do not distort

#########################################################################
# TAO: Energy minimization
-tao_method tao_lmvm
-tao_ls_ftol 1e-20
-tao_ls_rtol 1e+60
-tao_ls_gtol 1.0
-tao_fatol 0.0
-tao_frtol 1e-5
-tao_max_its 2000
-tao_monitor

#########################################################################
# PVODE: LLG time integration
-ts_init_time 0         #Initial time (TSSetInitialTime) (unit: ns)
-ts_max_time 5.0e10  #Time to run to (TSSetDuration) (unit: ns)

-psolve_ksp_atol 1e-7
-psolve_ksp_rtol 0.01
-psolve_ksp_divtol 100

mpdboot --totalnum=4 --file=mpd.hosts.4 --ncpus=8

mpirun -n 10 /storage/fizica/lstoleriu/test/disk/magpar.exe > X.dat &

mpirun -np 12 /storage/fizica/lstoleriu/test/disk/magpar.exe > disk.stdout.log < /dev/null &


nohup mpirun -l -np 12 /storage/fizica/lstoleriu/test/disk/magpar.exe > disk.stdout.log < /dev/null &

mkinp disk.0001.femsh *.gz

mkinp_vect disk.0001.inp

********************************************

primul
-size 1.0e-8

al doilea
-size 1.0e-7