diff options
Diffstat (limited to 'Physics_for_BSc(Paper-3)/Chapter5.ipynb')
| -rwxr-xr-x | Physics_for_BSc(Paper-3)/Chapter5.ipynb | 384 |
1 files changed, 384 insertions, 0 deletions
diff --git a/Physics_for_BSc(Paper-3)/Chapter5.ipynb b/Physics_for_BSc(Paper-3)/Chapter5.ipynb new file mode 100755 index 00000000..3ad752c1 --- /dev/null +++ b/Physics_for_BSc(Paper-3)/Chapter5.ipynb @@ -0,0 +1,384 @@ +{
+ "metadata": {
+ "name": "",
+ "signature": "sha256:4564535166c2472d72e2dcbe6fdb282dd0ebdfe6b75d04aed6f0238329c6aca1"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "5: Magnetostatics"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example number 5.1, Page number 118"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#importing modules\n",
+ "import math\n",
+ "from __future__ import division\n",
+ "\n",
+ "#Variable declaration\n",
+ "m=1.7*10**-27; #mass of proton(kg)\n",
+ "q=1.6*10**-19; #charge of proton(C)\n",
+ "KE=6*10**6*q; #kinetic energy of proton(J) \n",
+ "B=2; #intensity of magnetic field(Wb/m**2)\n",
+ "theta=90*math.pi/180; #angle(radian)\n",
+ "\n",
+ "#Calculation\n",
+ "v=math.sqrt(2*KE/m); #velocity(m/s)\n",
+ "Fb=q*v*B*math.sin(theta); #force acting on proton(N)\n",
+ "\n",
+ "#Result\n",
+ "print \"force acting on proton is\",round(Fb*10**11,3),\"*10**-11 N\"\n",
+ "print \"answer given in the book is wrong\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "force acting on proton is 1.075 *10**-11 N\n",
+ "answer given in the book is wrong\n"
+ ]
+ }
+ ],
+ "prompt_number": 12
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example number 5.2, Page number 119"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#importing modules\n",
+ "import math\n",
+ "from __future__ import division\n",
+ "\n",
+ "#Variable declaration\n",
+ "r=6.4*10**6; #radius of earth(m)\n",
+ "P=6.4*10**21; #magnetic moment(Amp m**2)\n",
+ "\n",
+ "#Calculation\n",
+ "A=math.pi*(r**2); #area(m**2)\n",
+ "i=P/A; #current(amp)\n",
+ "\n",
+ "#Result\n",
+ "print \"current is\",round(i/10**7),\"*10**7 Amp\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "current is 5.0 *10**7 Amp\n"
+ ]
+ }
+ ],
+ "prompt_number": 13
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example number 5.3, Page number 119"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#importing modules\n",
+ "import math\n",
+ "from __future__ import division\n",
+ "\n",
+ "#Variable declaration\n",
+ "B0=1.7*10**-5; #magnetic flux density(Wb/m**2)\n",
+ "mew0=4*math.pi*10**-7; \n",
+ "\n",
+ "#Calculation\n",
+ "H=B0/mew0; #magnetic intensity(Amp/m)\n",
+ "\n",
+ "#Result\n",
+ "print \"magnetic intensity is\",round(H,1),\"Amp/m\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "magnetic intensity is 13.5 Amp/m\n"
+ ]
+ }
+ ],
+ "prompt_number": 14
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example number 5.4, Page number 119"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#importing modules\n",
+ "import math\n",
+ "from __future__ import division\n",
+ "\n",
+ "#Variable declaration\n",
+ "l=30*10**-2; #length(m)\n",
+ "A=10**-4; #area(m**2)\n",
+ "N=300; #number of turns\n",
+ "i=0.032; #current(ampere)\n",
+ "phi=2*10**-6; #magnetic flux(Wb)\n",
+ "mew0=4*math.pi*10**-7; \n",
+ "\n",
+ "#Calculation\n",
+ "B=phi/A; #magnetic flux density(wb/m**2)\n",
+ "H=N*i/l; #magnetising force(amp turns/m)\n",
+ "mew=B/H; #permeability of ion(Wb/Amp m)\n",
+ "mewr=mew/mew0; #relative permeability\n",
+ "\n",
+ "#Result\n",
+ "print \"magnetic flux density is\",B*10**2,\"*10**-2 wb/m**2\"\n",
+ "print \"magnetising force is\",H,\"amp turns/m\"\n",
+ "print \"permeability of ion is\",mew*10**4,\"*10**-4 Wb/Amp m\"\n",
+ "print \"relative permeability is\",round(mewr)\n",
+ "print \"answer for relative permeability differs due to round off errors\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "magnetic flux density is 2.0 *10**-2 wb/m**2\n",
+ "magnetising force is 32.0 amp turns/m\n",
+ "permeability of ion is 6.25 *10**-4 Wb/Amp m\n",
+ "relative permeability is 497.0\n",
+ "answer for relative permeability differs due to round off errors\n"
+ ]
+ }
+ ],
+ "prompt_number": 20
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example number 5.5, Page number 120"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#importing modules\n",
+ "import math\n",
+ "from __future__ import division\n",
+ "\n",
+ "#Variable declaration\n",
+ "M=600; #magnetic moment(cgs units)\n",
+ "V=20; #volume(cm**3)\n",
+ "H=50; #magnetic field(oersteds)\n",
+ "\n",
+ "#Calculation\n",
+ "I=M/V; #intensity of magnetisation(cgs units)\n",
+ "B=H+(4*math.pi*I); #flux density(guass)\n",
+ "\n",
+ "#Result\n",
+ "print \"flux density is\",round(B),\"guass\"\n",
+ "print \"answer given in the book is wrong\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "flux density is 427.0 guass\n",
+ "answer given in the book is wrong\n"
+ ]
+ }
+ ],
+ "prompt_number": 23
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example number 5.6, Page number 120"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#importing modules\n",
+ "import math\n",
+ "from __future__ import division\n",
+ "\n",
+ "#Variable declaration\n",
+ "m=12000; #mass(gm)\n",
+ "d=7.5; #density(g/cm**3)\n",
+ "new=50; #frequency(Hz)\n",
+ "a=3000; #area of hysterisis curve(erg/cm**3)\n",
+ "\n",
+ "#Calculation\n",
+ "V=m/d; #volume of core(cm**3)\n",
+ "El=a*V; #energy loss in 1 cycle(erg)\n",
+ "n=new*60*60; #number of cycles per hour\n",
+ "EL=El*n/10**7; #energy loss per hour(J)\n",
+ "\n",
+ "#Result\n",
+ "print \"energy loss per hour is\",EL/10**4,\"*10**4 J\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "energy loss per hour is 8.64 *10**4 J\n"
+ ]
+ }
+ ],
+ "prompt_number": 24
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example number 5.7, Page number 121"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#importing modules\n",
+ "import math\n",
+ "from __future__ import division\n",
+ "\n",
+ "#Variable declaration\n",
+ "B=8*math.pi; #magnetic flux density(wb/m**2)\n",
+ "H=2000; #magnetic field(amp turn/m)\n",
+ "mew0=4*math.pi*10**-7; \n",
+ "\n",
+ "#Calculation\n",
+ "mewr=B/(mew0*H); #relative permeability\n",
+ "chi=mewr-1; #susceptibility\n",
+ "\n",
+ "#Result\n",
+ "print \"susceptibility is\",chi"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "susceptibility is 9999.0\n"
+ ]
+ }
+ ],
+ "prompt_number": 25
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example number 5.8, Page number 121"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#importing modules\n",
+ "import math\n",
+ "from __future__ import division\n",
+ "\n",
+ "#Variable declaration\n",
+ "q=1.6*10**-19; #charge on electron(C)\n",
+ "i1=3;\n",
+ "j1=2;\n",
+ "k1=3; #components of vector\n",
+ "i2=3;\n",
+ "j2=4; \n",
+ "k2=0; #components of vector\n",
+ "i3=0;\n",
+ "j3=4;\n",
+ "k3=1; #components of vector\n",
+ "\n",
+ "#Calculation\n",
+ "i4=(j2*k3)-(j3*k2);\n",
+ "j4=(i2*k3)-(i3*k2);\n",
+ "k4=(i2*j3)-(i3*j2); #components of derived vector\n",
+ "i5=i1+i4;\n",
+ "j5=j1+j4;\n",
+ "k5=k1+k4;\n",
+ "F=q*math.sqrt((i5**2)+(j5**2)+(k5**2)); #Lorentz force(N)\n",
+ "\n",
+ "#Result\n",
+ "print \"Lorentz force is\",round(F*10**19,2),\"*10**-19 N\"\n",
+ "print \"answer given in the book is wrong\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Lorentz force is 27.67 *10**-19 N\n",
+ "answer given in the book is wrong\n"
+ ]
+ }
+ ],
+ "prompt_number": 29
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
\ No newline at end of file |