I put the excel file in google drive. Here is the link.

https://drive.google.com/file/d/0BycN9tiDv9kmODlZRXl2alJ1SEE/view?usp=sharing

Best regards,

ryan

hi,

i am very interested in the Excel sheet which you used for schroedinger eq.

Kind regards

Wolfgang Schmidt

germany

anyone get that kind of information in such a perfect approach of

writing? I have a presentation next week, and I am at the look for such info. ]]>

Nonetheless I am here now and would just like to say

kudos for a marvelous post and a all round entertaining blog (I also love the theme/design),

I don’t have time to look over it all at the minute but I have book-marked it and also added your RSS feeds,

so when I have time I will be back to read a lot more, Please do keep

up the fantastic b. ]]>

I’m a student. I have an exercise that it looks like your post. Could you help me? ]]>

Below you find the context I included the equation.

4. Tools 14

4.1 Biot Savart and Lorentz 14

4.2 Finite straight wire 15

..1 Field 15

..2 Force 16

..3 Force between round wires 17

4.3 Finite flat plane 18

..1 Field 18

..2 Force 21

..3 Non-cuboidal magnets 23

4.4 Infinite bar 25

..1 Field 25

..2 Force 25

4.5 Finite bar 27

..1 Field 27

..2 Force 28

4.6 Cuboidal magnet 33

..1 Field 33

..2 Force 34

4.7 Round coil or magnet 38

..1 Field 38

..2 Force 42

..3 The mutual inductance, coaxial 43

..4 Self-inductance 44

..5 Mutual inductance between non-coaxial thin wall coils 45

4.8 Self-, mutual inductance and forces of a coils with rectangular shape and cross section 49

..1 Self-inductance 49

..2 Mutual Inductance 55

..3 Force 56

..4 Neumannâs formula 58

4.9 Presence of iron 63

..5 Reluctance force 65

..6 Attraction force on iron 66

4.7 Round coil or magnet

Field

To determine the field of a circular magnet like a disc or cylinder one needs at first the field of a circular loop with

an infinitely thin conductor. One has to apply Biot-Savart in cylinder coordinates as done in [19] and [20].

It is helpful to use the full elliptic integrals of the first and second kind, defined as:

(4.7.1)

The derivation via the vector potential is described in [21]. The vector potential is derived as:

(4.7.2)

with , R as radius and I as current at (r,z) as observation location.

The current has only one component directed tangentially, so in the direction Î¸.

Consequently the same holds for the vector potential. The field components are obtained via:

(4.7.3)

(4.7.4)

As example a picture of a single loop of infinitely thin conductor carrying with 1 Amp:

Fig. 4.7.1, the field lines and flux density of a single loop with its centre as reference

The field lines are drawn in the left-hand figure. The same amount of flux can be found between two field lines.

To get the field lines in an axially symmetric case one should plot r.AÎ¸ .

Prove: from follows with the contour C enclosing surface S,

as the normal unit vector on S and the unit vector along C.

For the flux holds , so .

The vector potential only has a tangential component and this component is not a function

of the tangential position in this case of axial symmetry.

Doing an integration at a fixed value of z over a coaxial circle r leads to

Consequently one should plot lines with a constant value of to get the field lines.

The previous expressions are the base for the field of a coil with an infinitely thin wall.

By integration of the expressions for a single loop over the height of the coil one obtains the field.

The link to a round magnet is that a magnet can be considered as a number of thin coils as e.g. a

round magnet with a hole in its centre.

It is modelled according the following model as two coaxial coils with opposite current direction.

Fig. 4.7.2, a ring magnet as a set of coils

The number of ampere turns is equal to the HcJ times the magnet height.

Representation of a magnet by coils as indicated is fully allowed when:

– the relative permeability Î¼r is near 1 (as holds for NdFeB, SmCo and ferroxdure)

– the point of operation can be found in the linear part of the BH-curve

– the magnetization is uniform and unidirectional.

At first one should define the function J as the complete elliptic integral of the third kind:

(4.7.5)

The full exercise to get the field for a coil is done in [22]. A fast numerical procedure for this integral is given in [23].

(4.7.6)

with: R radius coil,

H half height of the coil,

r, z coordinates,

N number of turns,

I current,

Z=Z(i)

k=k(i) .

The elliptic integral E, K and J are not always available as âin houseâ procedure. Appendix 10 gives a Mathematica implementation.

Fig. 4.7.3 radial field at N.I=-50 amp

Fig. 4.7.4 axial field at N.I=-50 amp

Another recent source for the field of a permanent magnet ring is given in [24].

T.L. Tang published in [25] the following equation for the modified vector potential AÏ(r,z) of a cylindrical coil with N turns, the current I, the radius R and height 2H:

(4.7.7)

with

The combination of the field components and vector potential leads to the following figures.

Fig. 4.7.5 potential lines and field direction of a cylindrical coil

Fig. 4.7.6 potential lines and field direction for two cylindrical coils with opposing current direction

The field of thin wall cylindrical coils and disk and ring magnets can be described well with the previous expressions.

Expressions for the field of thick wall cylindrical coils are still under study in publications,

because a compact formulation is not found yet. Of course, Finite Element Analysis will give correct values,

but a compact analytical expression would be very nice for optimization.

One way to solve this is to cut the thick wall coils in layers and sum the contributions to get the field.

This is demonstrated with the following figures, showing the field of a coil with a circular and rectangular

cross section respectively. The number of ampere-turns equals 1000 and the number of layers equals 40.

Fig. 4.7.7 potential lines, field strength and field direction of a circular loop with a round cross section and uniform current density

Fig. 4.7.8 potential lines, field strength and field direction of a circular loop with a rectangular cross section and uniform current density

Flat coils are applied on printed circuit boards and in wireless energy transfer systems.

Analytical descriptions of the magnetic field of a disk coil have become available in the last 10 years,

but the intended advantages are lost to a high extend because its solution process becomes more

complex and/or time consuming than adding the fields of e.g. 10 concentric rings as given in Fig. 4.7.1.

As example 10 rings with each 100 Amp. with a radius between 0.1 and 0.2 m. leads to the following figure.

Fig. 4.7.9 potential lines, field strength and field direction of a disk coil with 10 wires carrying 100 Amp. each.

]]>