A discrete magnetostatic solver in Matlab/Octave
|
|
- Elwin Bishop
- 6 years ago
- Views:
Transcription
1 discrete magnetostatic solver in Matlab/Octave Tobias Schäfer 1st March 01 Contents 1 Introduction 1 Theory 1 Boundary conditions 5 3 Solution 5 31 Iterative solution Hint 5 3 FEM Solution 6 4 Matlab/Octave code 6 1 Introduction In the open-rtms project a complex magnetic field is applied to the human brain The whole setup consists of several coils placed around the skull I assume, that the properties of the brain are LTI (linear, time invariant), so that a simple superposition of the fields from the different coils can be calculated the get the field distribution inside the brain The coils I assume to be small compared to the distance between them and the field strength is small as well, so that the influence of the coils towards each other can also be neglected These assumptions should work out well, if following types of coils are used: a simple core made from a screw with some turns of copper ferrite core is not needed, because of the low frequency flat coils without core In this case also overlapping coils are not a problem 1
2 THEORY Theory To calculate the magnetic fieldlines we first calculate the vector potential fter is computed we can calculate the magnetic field B = rot() = To calculate we start with the Poisson equation in 3D (: permeability of material, j: current density in coil): ( 1 x ) ( 1 x y ) ( 1 y z ) = j (1) z To simplify the solution of equation (1) the problem is simplified to a D rotational system (Figure ) r Mesh cells ir z Core Coil Figure 1: Rotational symmetric setup in this rotational case the Poisson equation (1) becomes: ( ) 1 (r ) ( 1 r r r z ) z = j () This equation is solved for a grid of cells The terms in () are rewritten as finite differences: The inner partial differentials according to figure will be and (r ) r (r r/) (r r/) (r r/) (r r/) r z (z z/) (z z/) z (3) (4)
3 THEORY (r r,z) (r r,z) (r,z z) (r, z) (r,z z) (r,z z) (r,z) (r,z z) (r r,z) (r r,z) Figure : Cells in the grid For the whole expression () we get ( ) 1 (r ) r r r (r r) (r r) r r (r r) (r r) (r r/) (r r/) r (r r/) (r r/) r (5) and ( 1 z ) z (z z) (z z/) z (z z) (z z/) z z (z z) (z z/) z (z z) (z z/) z (6) s it can be seen from figure is used in equations (5) and (6) on a 45 rotated grid between the cells For this grid is hard to define in Matlab/Octave without making errors, we define at the very same points as (as shown in figure ) and use terms like = (z z/) = (z) (z z) (7) 3
4 THEORY Figure 3: inbetween Equations (5) and (6) become and (r r) (r r) r (r r) (r r) r (r r/) ( (r r)) r (r r/) ( (r r)) r (8) (z z) (z z) ( (z z)) z ( (z z)) z (9) Everything put together and all the different taken apart: (r r) (r r,z) (r r/) ( (r r)) r (r r) (r r,z) (r r/) ( (r r)) r (r,z z) ( (z z)) z (r,z z) ( (z z)) z ( r (r,z) (r r/) ( (r r)) r r (r r/) ( (r r)) r ) ( (z z)) z ( (z z)) z = j (10) 4
5 3 SOLUTION 1 Boundary conditions This lengthy equation (10) is solved for (r,z) 1 Boundary conditions In theory the magnetic field lines (or actually field planes) extent into infinity The simulated volume is finite 3 Solution 31 Iterative solution - and a j-matrix are given The -matrix is initialized to zero The matrix is calculated iteratively For every step k = 1,,3, Calculate k (r,z) from equation (10) using the values of the surrounding cells from the previous step ( k 1 (r r,z), k 1 (r r,z), k 1 (r,z z), k 1 (r,z z)) Calculate the error in cell (r,z): err = k (r,z) k 1 (r,z) Correct the error in cell (r,z): k (r,z) = k 1 (r,z) err V Repeat V is a twiddel-factor: V < 0 instable V 0 exact solution V 1 fast solution 1 < V < oscillates but might become stable eventually V > instable 311 Hint If one looks closely at the error calculation and correction one can see a discrete 1st order lowpass filter (PT-1) here Transfer behavior from x k to y k : Transfer-function: y k = y k 1 V (x k y k 1 ) Y (z) = Y (z) z 1 V (X(z) Y (z) z 1 Y (z) X(z) = V 1 (V 1) z 1 G(z) = V z z (V 1) The filter has a zero at 0 ı0 and a pole at (1 V ) ı0 For V < 0 or V > the pole leaves the unit circle and the filter becomes instable 5
6 3 FEM Solution 4 MTLB/OCTVE CODE 3 FEM Solution The equation 1 can be rewritten in the form C = j with a vector of the vectorpotentials in very cell, j a vector of the current density in every cell and C a matrix of constants With a simple inversion of C we get the solution for : = C 1 j 4 Matlab/Octave code 6
EEL3135: Homework #4
EEL335: Homework #4 Problem : For each of the systems below, determine whether or not the system is () linear, () time-invariant, and (3) causal: (a) (b) (c) xn [ ] cos( 04πn) (d) xn [ ] xn [ ] xn [ 5]
More informationStability Analysis Techniques
Stability Analysis Techniques In this section the stability analysis techniques for the Linear Time-Invarient (LTI) discrete system are emphasized. In general the stability techniques applicable to LTI
More informationCH 19-1 Magnetic Field
CH 19-1 Magnetic Field Important Ideas A moving charged particle creates a magnetic field everywhere in space around it. If the particle has a velocity v, then the magnetic field at this instant is tangent
More informationModelling of Specimen Interaction with Ferrite Cored Coils by Coupling Semi-Analytical and Numerical Techniques
Modelling of Specimen Interaction with Ferrite Cored Coils by Coupling Semi-Analytical and Numerical Techniques A. Skarlatos, E. Demaldent, A. Vigneron and C. Reboud 25 th 28 th of June, Bratislava, Slovak
More informationEE102B Signal Processing and Linear Systems II. Solutions to Problem Set Nine Spring Quarter
EE02B Signal Processing and Linear Systems II Solutions to Problem Set Nine 202-203 Spring Quarter Problem 9. (25 points) (a) 0.5( + 4z + 6z 2 + 4z 3 + z 4 ) + 0.2z 0.4z 2 + 0.8z 3 x[n] 0.5 y[n] -0.2 Z
More informationThe Z transform (2) 1
The Z transform (2) 1 Today Properties of the region of convergence (3.2) Read examples 3.7, 3.8 Announcements: ELEC 310 FINAL EXAM: April 14 2010, 14:00 pm ECS 123 Assignment 2 due tomorrow by 4:00 pm
More informationECE 350 Signals and Systems Spring 2011 Final Exam - Solutions. Three 8 ½ x 11 sheets of notes, and a calculator are allowed during the exam.
ECE 35 Spring - Final Exam 9 May ECE 35 Signals and Systems Spring Final Exam - Solutions Three 8 ½ x sheets of notes, and a calculator are allowed during the exam Write all answers neatly and show your
More informationCourse no. 4. The Theory of Electromagnetic Field
Cose no. 4 The Theory of Electromagnetic Field Technical University of Cluj-Napoca http://www.et.utcluj.ro/cs_electromagnetics2006_ac.htm http://www.et.utcluj.ro/~lcret March 19-2009 Chapter 3 Magnetostatics
More informationDigital Signal Processing. Midterm 1 Solution
EE 123 University of California, Berkeley Anant Sahai February 15, 27 Digital Signal Processing Instructions Midterm 1 Solution Total time allowed for the exam is 8 minutes Some useful formulas: Discrete
More informationMagnetic Fields Due to Currents
PHYS102 Previous Exam Problems CHAPTER 29 Magnetic Fields Due to Currents Calculating the magnetic field Forces between currents Ampere s law Solenoids 1. Two long straight wires penetrate the plane of
More informationM. Popnikolova Radevska, Member, IEEE, M. Cundev, L. Pelkovska
MK0500056 Electromagnetic Characteristics and Static Torque of a Solid Salient Poles Synchronous Motor Computed by 3D-Finite Element Method Magnetics (September 2000) M. Popnikolova Radevska, Member, IEEE,
More informationHomework # Physics 2 for Students of Mechanical Engineering. Part A
Homework #9 203-1-1721 Physics 2 for Students of Mechanical Engineering Part A 5. A 25-kV electron gun in a TV tube fires an electron beam having a diameter of 0.22 mm at the screen. The spot on the screen
More informationFEM: Domain Decomposition and Homogenization for Maxwell s Equations of Large Scale Problems
FEM: and Homogenization for Maxwell s Equations of Large Scale Problems Karl Hollaus Vienna University of Technology, Austria Institute for Analysis and Scientific Computing February 13, 2012 Outline 1
More informationRatio of Charge to Mass for the Electron
Ratio of Charge to Mass for the Electron For a positive charge moving in a uniform magnetic field B with velocity v, the force F on the charge is always perpendicular to the magnetic field and the velocity.
More informationZ Transform (Part - II)
Z Transform (Part - II). The Z Transform of the following real exponential sequence x(nt) = a n, nt 0 = 0, nt < 0, a > 0 (a) ; z > (c) for all z z (b) ; z (d) ; z < a > a az az Soln. The given sequence
More informationIntegrated Square-Shaped Spiral Inductor
Page 1 of 9 Electrical Component Models : Integrated Square-Shaped Spiral Inductor Integrated Square-Shaped Spiral Inductor This example presents a model of a micro-scale square inductor, used for LC bandpass
More informationPHYS102 Previous Exam Problems. Induction
PHYS102 Previous Exam Problems CHAPTER 30 Induction Magnetic flux Induced emf (Faraday s law) Lenz law Motional emf 1. A circuit is pulled to the right at constant speed in a uniform magnetic field with
More informationApplications of Ampere s Law
Applications of Ampere s Law In electrostatics, the electric field due to any known charge distribution ρ(x, y, z) may alwaysbeobtainedfromthecoulomblaw it sauniversal tool buttheactualcalculation is often
More informationStepper Motor. Author: ElectroMagneticWorks. Company: ElectroMagneticWorks Inc.
1. Introduction 2. Model View 3. Materials 4. Load & Restraint Information 5. Coils Information 6. Force and Torque Information 7. Study Properties 8. Results Table 9. Magnetic Flux Density Results 10.
More information(i) Represent discrete-time signals using transform. (ii) Understand the relationship between transform and discrete-time Fourier transform
z Transform Chapter Intended Learning Outcomes: (i) Represent discrete-time signals using transform (ii) Understand the relationship between transform and discrete-time Fourier transform (iii) Understand
More informationEC Objective Paper I (Set - D)
EC-Objective Paper-I ESE-5 www.gateforum.com EC Objective Paper I (Set - D). If a system produces frequencies in the output are not present in the input, then the system cannot be Minimum phase system
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi System Stability - 26 March, 2014
Prof. Dr. Eleni Chatzi System Stability - 26 March, 24 Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can
More informationMagnetostatic Analysis of Solenoid
Magnetostatic Analysis of Solenoid 1. Introduction 2. Model View 3. Materials 4. Load & Restraint Information 5. Coils Information 6. Force and Torque Information 7. Study Properties 8. Results Table 9.
More informationThe Influence of Core Shape and Material Nonlinearities to Corner Losses of Inductive Element
The Influence of Core Shape and Material Nonlinearities to Corner Losses of Inductive Element Magdalena Puskarczyk 1, Brice Jamieson 1, Wojciech Jurczak 1 1 ABB Corporate Research Centre, Kraków, Poland
More information(a) What are the probabilities associated with finding the different allowed values of the z-component of the spin after time T?
1. Quantum Mechanics (Fall 2002) A Stern-Gerlach apparatus is adjusted so that the z-component of the spin of an electron (spin-1/2) transmitted through it is /2. A uniform magnetic field in the x-direction
More informationChapter 5. Magnetostatics
Chapter 5. Magnetostatics 5.4 Magnetic Vector Potential 5.1.1 The Vector Potential In electrostatics, E Scalar potential (V) In magnetostatics, B E B V A Vector potential (A) (Note) The name is potential,
More informationLecture 3 Matlab Simulink Minimum Phase, Maximum Phase and Linear Phase Systems
Lecture 3 Matlab Simulink Minimum Phase, Maximum Phase and Linear Phase Systems Lester Liu October 31, 2012 Minimum Phase, Maximum Phase and Linear Phase LTI Systems In this section, we will explore the
More informationSensibility Analysis of Inductance Involving an E-core Magnetic Circuit for Non Homogeneous Material
Sensibility Analysis of Inductance Involving an E-core Magnetic Circuit for Non Homogeneous Material K. Z. Gomes *1, T. A. G. Tolosa 1, E. V. S. Pouzada 1 1 Mauá Institute of Technology, São Caetano do
More informationA cylinder in a magnetic field (Jackson)
Problem 1. A cylinder in a magnetic field (Jackson) A very long hollow cylinder of inner radius a and outer radius b of permeability µ is placed in an initially uniform magnetic field B o at right angles
More informationPre Lab for Ratio of Mass to. Charge of an Electron
Pre Lab for Ratio of Mass to Charge of an Electron The direction of the magnetic force on a charged particle moving in the magnetic field is given by the right hand rule. Students need practice using the
More informationELEC 1908 The Electric Potential (V) March 28, 2013
ELEC 1908 The Electric Potential (V) March 28, 2013 1 Abstract The objective of this lab is to solve numerically Laplace s equation in order to obtain the electric potential distribution in di erent electric
More informationNumerical and experimental identification of the static characteristics of a combined Journal-Magnetic bearing: Smart Integrated Bearing
Numerical and experimental identification of the static characteristics of a combined Journal-Magnetic bearing: Smart Integrated Bearing M El-Hakim*, A S Dimitri**, T Saqr*, J Mahfoud***, A Adly**, A El-
More informationFinite Element Method (FEM)
Finite Element Method (FEM) The finite element method (FEM) is the oldest numerical technique applied to engineering problems. FEM itself is not rigorous, but when combined with integral equation techniques
More informationComputations of Vector Potential and Toroidal Flux and Applications to Stellarator Simulations
Computations of Vector Potential and Toroidal Flux and Applications to Stellarator Simulations NIMROD Team Meeting Torrin Bechtel April 30, 2017 Outline 1 Project Goals and Progress 2 Vector Potential
More informationMeasurement of electric potential fields
Measurement of electric potential fields Matthew Krupcale, Oliver Ernst Department of Physics, Case Western Reserve University, Cleveland Ohio, 44106-7079 18 November 2012 Abstract In electrostatics, Laplace
More informationMinimizing the Detent Force in Permanent Magnet Linear Synchronous Motor for driving of 2D Laser Marking Table
Minimizing the Detent Force in Permanent Magnet Linear Synchronous Motor for driving of 2D Laser Marking Table Lyubomir Lazov 1, Peter Uzunov 2 Latvia Academy of Technologies, Rezekne, Latvia 1 Electricity
More informationA half submerged metal sphere (UIC comprehensive
Problem 1. exam) A half submerged metal sphere (UIC comprehensive A very light neutral hollow metal spherical shell of mass m and radius a is slightly submerged by a distance b a below the surface of a
More informationName: Class: Date: Multiple Choice Identify the letter of the choice that best completes the statement or answers the question.
Name: Class: Date: AP REVIEW 4 Multiple Choice Identify the letter of the choice that best completes the statement or answers the question. 1. If a positively charged glass rod is used to charge a metal
More informationProblems in Magnetostatics
Problems in Magnetostatics 8th February 27 Some of the later problems are quite challenging. This is characteristic of problems in magnetism. There are trivial problems and there are tough problems. Very
More informationELECTROMAGNETISM. Challenging MCQ questions by The Physics Cafe. Compiled and selected by The Physics Cafe
ELECTROMGNETISM Challenging MCQ questions by The Physics Cafe Compiled and selected by The Physics Cafe 1 Fig 31 shows the top view of a current balance. The wire frame CD is supported by the pivots P
More informationPHYSICS 202 Practice Exam Magnetism, Induction, Simple Harmonic Motion. Name. Constants and Conversion Factors
PHYSICS 202 Practice Exam Magnetism, Induction, Simple Harmonic Motion Name Constants and Conversion Factors "* electron charge, / œ 1.6 10 Coulombs $" electron mass, 7/ œ 9.11 10 kg 1 * # # Coulomb's
More informationA Simple Weak-Field Coupling Benchmark Test of the Electromagnetic-Thermal-Structural Solution Capabilities of LS-DYNA Using Parallel Current Wires
13 th International LS-DYNA Users Conference Session: Electromagnetic A Simple Weak-Field Coupling Benchmark Test of the Electromagnetic-Thermal-Structural Solution Capabilities of LS-DYNA Using Parallel
More informationA Magnetohydrodynamic study af a inductive MHD generator
Excerpt from the Proceedings of the COMSOL Conference 2009 Milan A Magnetohydrodynamic study af a inductive MHD generator Augusto Montisci, Roberto Pintus University of Cagliari, Department of Electrical
More informationMagnetic Field Simulation of a Miniature Circuit Breaker s Coil
Magnetic Field Simulation of a Miniature Circuit Breaker s Coil Comparison of Simulation and analitical results on the solenoid axis. Conecici-Lucian Madalin, Munteanu Calin, Purcar Ioan Marius Technical
More informationPhysics 54 Lecture March 1, Micro-quiz problems (magnetic fields and forces) Magnetic dipoles and their interaction with magnetic fields
Physics 54 Lecture March 1, 2012 OUTLINE Micro-quiz problems (magnetic fields and forces) Magnetic dipoles and their interaction with magnetic fields Electromagnetic induction Introduction to electromagnetic
More informationADAM PIŁAT Department of Automatics, AGH University of Science and Technology Al. Mickiewicza 30, Cracow, Poland
Int. J. Appl. Math. Comput. Sci., 2004, Vol. 14, No. 4, 497 501 FEMLAB SOFTWARE APPLIED TO ACTIVE MAGNETIC BEARING ANALYSIS ADAM PIŁAT Department of Automatics, AGH University of Science and Technology
More informationMutual Inductance. The field lines flow from a + charge to a - change
Capacitors Mutual Inductance Since electrical charges do exist, electric field lines have a starting point and an ending point. For example, if you have a + and a - change, the field lines would look something
More informationNitsche-type Mortaring for Maxwell s Equations
Nitsche-type Mortaring for Maxwell s Equations Joachim Schöberl Karl Hollaus, Daniel Feldengut Computational Mathematics in Engineering at the Institute for Analysis and Scientific Computing Vienna University
More informationFEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment
FEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment Alastair J. Radcliffe Andreas Dedner Timo Betcke Warwick University, Coventry University College of London (UCL) U.K. Radcliffe
More informationThe Helmholtz Decomposition and the Coulomb Gauge
The Helmholtz Decomposition and the Coulomb Gauge 1 Problem Kirk T. McDonald Joseph Henry Laboratories, Princeton University, Princeton, NJ 08544 (April 17, 2008; updated February 7, 2015 Helmholtz showed
More informationFrequency Response of Discrete-Time Systems
Frequency Response of Discrete-Time Systems EE 37 Signals and Systems David W. Graham 6 Relationship of Pole-Zero Plot to Frequency Response Zeros Roots of the numerator Pin the system to a value of ero
More informationResponses of Digital Filters Chapter Intended Learning Outcomes:
Responses of Digital Filters Chapter Intended Learning Outcomes: (i) Understanding the relationships between impulse response, frequency response, difference equation and transfer function in characterizing
More informationIntroduction. HFSS 3D EM Analysis S-parameter. Q3D R/L/C/G Extraction Model. magnitude [db] Frequency [GHz] S11 S21 -30
ANSOFT Q3D TRANING Introduction HFSS 3D EM Analysis S-parameter Q3D R/L/C/G Extraction Model 0-5 -10 magnitude [db] -15-20 -25-30 S11 S21-35 0 1 2 3 4 5 6 7 8 9 10 Frequency [GHz] Quasi-static or full-wave
More informationChapter 22, Magnetism. Magnets
Chapter 22, Magnetism Magnets Poles of a magnet (north and south ) are the ends where objects are most strongly attracted. Like poles repel each other and unlike poles attract each other Magnetic poles
More informationDevelopment of a Compressor for a Miniature Pulse Tube Cryocooler of 2.5 W at 65 K
Development of a Compressor for a Miniature Pulse Tube Cryocooler of 2.5 W at 65 K N. Matsumoto, Y. Yasukawa, K. Ohshima, T. Takeuchi, K. Yoshizawa, T. Matsushita, Y. Mizoguchi, and A. Ikura Fuji Electric
More informationPHYSICS GRADUATE SCHOOL QUALIFYING EXAMINATION
PHYSICS GRADUATE SCHOOL QUALIFYING EXAMINATION Partl INSTRUCTIONS: Work all problems. This is a closed book examination. Start each problem on a new pack of yellow paper and use only one side of each sheet.
More informationUniversity of the Philippines College of Science PHYSICS 72. Summer Second Long Problem Set
University of the Philippines College of Science PHYSICS 72 Summer 2012-2013 Second Long Problem Set INSTRUCTIONS: Choose the best answer and shade the corresponding circle on your answer sheet. To change
More informationTHERMAL FIELD ANALYSIS IN DESIGN AND MANUFACTURING OF A PERMANENT MAGNET LINEAR SYNCHRONOUS MOTOR
THERMAL FIELD ANALYSIS IN DESIGN AND MANUFACTURING OF A PERMANENT MAGNET LINEAR SYNCHRONOUS MOTOR Petar UZUNOV 1 ABSTRACT: The modern Permanent Magnet Linear Synchronous Motors (PMLSM) has a wide range
More informationLecture 2 OKAN UNIVERSITY FACULTY OF ENGINEERING AND ARCHITECTURE
OKAN UNIVERSITY FACULTY OF ENGINEERING AND ARCHITECTURE EEE 43 DIGITAL SIGNAL PROCESSING (DSP) 2 DIFFERENCE EQUATIONS AND THE Z- TRANSFORM FALL 22 Yrd. Doç. Dr. Didem Kivanc Tureli didemk@ieee.org didem.kivanc@okan.edu.tr
More informationElectromagnetic Forces on Parallel Current-
Page 1 of 5 Tutorial Models : Electromagnetic Forces on Parallel Current-Carrying Wires Electromagnetic Forces on Parallel Current- Carrying Wires Introduction One ampere is defined as the constant current
More informationEE402 - Discrete Time Systems Spring Lecture 10
EE402 - Discrete Time Systems Spring 208 Lecturer: Asst. Prof. M. Mert Ankarali Lecture 0.. Root Locus For continuous time systems the root locus diagram illustrates the location of roots/poles of a closed
More informationPRACTICE EXAM 1 for Midterm 1
PRACTICE EXAM 1 for Midterm 1 Multiple Choice Questions 1) The figure shows three electric charges labeled Q 1, Q 2, Q 3, and some electric field lines in the region surrounding the charges. What are the
More informationInverse Problems and Optimal Design in Electricity and Magnetism
Inverse Problems and Optimal Design in Electricity and Magnetism P. Neittaanmäki Department of Mathematics, University of Jyväskylä M. Rudnicki Institute of Electrical Engineering, Warsaw and A. Savini
More informationProblem Set 2: Solution Due on Wed. 25th Sept. Fall 2013
EE 561: Digital Control Systems Problem Set 2: Solution Due on Wed 25th Sept Fall 2013 Problem 1 Check the following for (internal) stability [Hint: Analyze the characteristic equation] (a) u k = 05u k
More informationNumerical Computation of the Restoring Force in a Cylindrical Bearing Containing Magnetic Liquid
SERBIAN JOURNAL OF ELECTRICAL ENGINEERING Vol. 5, No., November 008, 83-90 Numerical Computation of the Restoring Force in a Cylindrical Bearing Containing Magnetic Liquid Marian Greconici, Gheorghe Madescu,
More informationEddy Current Losses in the Tank Wall of Power Transformers
Eddy Current Losses in the Tank Wall of Power Transformers Erich Schmidt Institute of Electrical Drives and Machines, Vienna University of Technology A 14 Vienna, Austria, Gusshausstrasse 25 29 Phone:
More informationElectromagnetic inverse problems in biomedical engineering Jens Haueisen
Electromagnetic inverse problems in biomedical engineering Jens Haueisen Institute of Biomedical Engineering and Informatics Technical University of Ilmenau Germany Biomagnetic Center Department of Neurology
More informationCHARGE TO MASS RATIO FOR THE ELECTRON
CHARGE TO MASS RATIO FOR THE ELECTRON OBJECTIVE: To measure the ratio of the charge of an electron to its mass. METHOD: A stream of electrons is accelerated by having them "fall" through a measured potential
More information( 1 jj) = p + j B (1) j = en (u i u e ) (2) ρu = n (m i u i + m e u e ) (3)
Magnetospheric Physics - Homework, 2/14/2014 11. MHD Equations: a Consider a two component electrons and ions charge neutral ρ c = 0 plasma where the total bulk velocity is defined by ρu = n m i u i +
More informationModeling and Design Optimization of Permanent Magnet Linear Synchronous Motor with Halbach Array
Modeling and Design Optimization of Permanent Magnet Linear Synchronous Motor with Halbach Array N. Roshandel Tavana, and A. Shoulaie nroshandel@ee.iust.ir, and shoulaie@ee.iust.ac.ir Department of Electrical
More informationReading Question 24.1
Reading Question 24.1 A compass in a magnetic field will line up A. With the north pole pointing in the direction of the magnetic field. B. With the north pole pointing opposite the direction of the magnetic
More informationfiziks Institute for NET/JRF, GATE, IIT-JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES
Content-ELECTRICITY AND MAGNETISM 1. Electrostatics (1-58) 1.1 Coulomb s Law and Superposition Principle 1.1.1 Electric field 1.2 Gauss s law 1.2.1 Field lines and Electric flux 1.2.2 Applications 1.3
More informationFaraday's Law ds B B G G ΦB B ds Φ ε = d B dt
Faraday's Law ds ds ε= d Φ dt Φ Global Review Electrostatics» motion of q in external E-field» E-field generated by Σq i Magnetostatics» motion of q and i in external -field» -field generated by I Electrodynamics»
More informationPractical Course Scientific Computing and Visualization
May 3, 24 Page 1 of 18 1 gridpoints for pressure p: j i p, p 1, p 2, p 3, p 4, p,1 p 1,1 p 2,1 p 3,1 p 4,1 p,2 p 1,2 p 2,2 p 3,2 p 4,2 p,3 p 1,3 p 2,3 p 3,3 p 4,3 p,4 p 1,4 p 2,4 p 3,4 p 4,4 Page 2 of
More informationOn Permanent Magnet Material Testing. 1. Permanent Magnet magnetization.
July 3 2006 CBN06-2 On Permanent Magnet Material Testing. Alexander Temnykh, LEPP, Ithaca, NY, 14853 Abstract. Note describes two experiments made to explore basic properties of permanent magnet (PM) material
More informationThe Direction of Magnetic Field. Neil Alberding (SFU Physics) Physics 121: Optics, Electricity & Magnetism Spring / 16
The Direction of Magnetic Field Neil Alberding (SFU Physics) Physics 121: Optics, Electricity & Magnetism Spring 2010 1 / 16 The Magnetic Field We introduced electric field to explain-away long-range electric
More informationE & M Qualifier. January 11, To insure that the your work is graded correctly you MUST:
E & M Qualifier 1 January 11, 2017 To insure that the your work is graded correctly you MUST: 1. use only the blank answer paper provided, 2. use only the reference material supplied (Schaum s Guides),
More informationEXAMPLE 115 NON-ISOTHERMAL FLOW IN A SINGLE SCREW EXTRUDER DESCRIPTION KEYWORDS FILENAMES. Example 115
EXAMPLE 115 NON-ISOTHERMAL FLOW IN A SINGLE SCREW EXTRUDER DESCRIPTION In this example, we evaluate the flow in the last part of a single screw extruder, including the last few conveying screw elements,
More informationSimulation study of the sensing field in electromagnetic tomography for two-phase flow measurement
Flow Measurement and Instrumentation 16 (2005) 199 204 www.elsevier.com/locate/flowmeasinst Simulation study of the sensing field in electromagnetic tomography for two-phase flow measurement Ze Liu a,,minhe
More informationEE123 Digital Signal Processing. M. Lustig, EECS UC Berkeley
EE123 Digital Signal Processing Today Last time: DTFT - Ch 2 Today: Continue DTFT Z-Transform Ch. 3 Properties of the DTFT cont. Time-Freq Shifting/modulation: M. Lustig, EE123 UCB M. Lustig, EE123 UCB
More informationComplex Analysis Topic: Singularities
Complex Analysis Topic: Singularities MA201 Mathematics III Department of Mathematics IIT Guwahati August 2015 Complex Analysis Topic: Singularities 1 / 15 Zeroes of Analytic Functions A point z 0 C is
More informationVersion The diagram below represents lines of magnetic flux within a region of space.
1. The diagram below represents lines of magnetic flux within a region of space. 5. The diagram below shows an electromagnet made from a nail, a coil of insulated wire, and a battery. The magnetic field
More informationPoisson Equation in 2D
A Parallel Strategy Department of Mathematics and Statistics McMaster University March 31, 2010 Outline Introduction 1 Introduction Motivation Discretization Iterative Methods 2 Additive Schwarz Method
More informationDigital Signal Processing IIR Filter Design via Bilinear Transform
Digital Signal Processing IIR Filter Design via Bilinear Transform D. Richard Brown III D. Richard Brown III 1 / 12 Basic Procedure We assume here that we ve already decided to use an IIR filter. The basic
More information0 t < 0 1 t 1. u(t) =
A. M. Niknejad University of California, Berkeley EE 100 / 42 Lecture 13 p. 22/33 Step Response A unit step function is described by u(t) = ( 0 t < 0 1 t 1 While the waveform has an artificial jump (difficult
More informationZero speed sensorless drive capability of fractional slot inset PM machine
Zero speed sensorless drive capability of fractional slot inset PM Adriano Faggion Nicola Bianchi Silverio Bolognani Emanuele Fornasiero Electric Drives Laboratory Department of Industrial Engineering
More informationTransform analysis of LTI systems Oppenheim and Schafer, Second edition pp For LTI systems we can write
Transform analysis of LTI systems Oppenheim and Schafer, Second edition pp. 4 9. For LTI systems we can write yœn D xœn hœn D X kd xœkhœn Alternatively, this relationship can be expressed in the z-transform
More informationChapter 1 Magnetic Circuits
Principles of Electric Machines and Power Electronics Third Edition P. C. Sen Chapter 1 Magnetic Circuits Chapter 1: Main contents i-h relation, B-H relation Magnetic circuit and analysis Property of magnetic
More informationAutomatic Control Systems theory overview (discrete time systems)
Automatic Control Systems theory overview (discrete time systems) Prof. Luca Bascetta (luca.bascetta@polimi.it) Politecnico di Milano Dipartimento di Elettronica, Informazione e Bioingegneria Motivations
More informationContinuous Casting - Comparison between Numerical Simulation and Experimental Measurement
International Scientific Colloquium Modelling for Electromagnetic Processing Hannover, March 24-26, 3 Continuous Casting - Comparison between Numerical Simulation and Experimental Measurement Y. Delannoy,
More informationDon t Copy This. Michael Faraday 3/1/13. Chapter 25: EM Induction and EM Waves. Key Terms:
3/1/13 Chapter 25: EM Induction and EM Waves Key Terms: Induced Induced emf current Flux Lenz s Law waves Photons EM Don t Copy This Last chapter we learned that a current can create a magnetic field.
More informationDesign of Orthonormal Wavelet Filter Banks Using the Remez Exchange Algorithm
Electronics and Communications in Japan, Part 3, Vol. 81, No. 6, 1998 Translated from Denshi Joho Tsushin Gakkai Ronbunshi, Vol. J80-A, No. 9, September 1997, pp. 1396 1402 Design of Orthonormal Wavelet
More informationControl Systems Engineering ( Chapter 6. Stability ) Prof. Kwang-Chun Ho Tel: Fax:
Control Systems Engineering ( Chapter 6. Stability ) Prof. Kwang-Chun Ho kwangho@hansung.ac.kr Tel: 02-760-4253 Fax:02-760-4435 Introduction In this lesson, you will learn the following : How to determine
More informationAtomic and nuclear physics
Atomic and nuclear physics Atomic shell Normal Zeeman effect LEYBOLD Physics Leaflets Observing the normal Zeeman effect in transverse and longitudinal Objects of the experiment Observing the line triplet
More informationIII.Sources of Magnetic Fields - Ampere s Law - solenoids
Magnetism I. Magnetic Field - units, poles - effect on charge II. Magnetic Force on Current - parallel currents, motors III.Sources of Magnetic Fields - Ampere s Law - solenoids IV.Magnetic Induction -
More informationSimplified Analysis Technique for Double Layer Non-overlap Multiphase Slip Permanent Magnet Couplings in Wind Energy Applications
Simplified Analysis Technique for Double Layer Non-overlap Multiphase Slip Permanent Magnet Couplings in Wind Energy Applications Petrus JJ van Wyk - Prof. Maarten J. Kamper Renewable Energy Postgraduate
More informationAnalytical Model for Sizing the Magnets of Permanent Magnet Synchronous Machines
Journal of Electrical Engineering 3 (2015) 134-141 doi: 10.17265/2328-2223/2015.03.004 D DAVID PUBLISHING Analytical Model for Sizing Magnets of Permanent Magnet Synchronous Machines George Todorov and
More informationCharacterization of an Electromagnetic Vibration Isolator
Journal of Electromagnetic Analysis and Applications, 211, 3, 519-528 doi:1.4236/jemaa.211.31279 Published Online December 211 (http://www.scirp.org/journal/jemaa) 519 Characterization of an Electromagnetic
More informationTrilogy of Magnetics
Trilogy of Magnetics Design Guide for EMI Filter Design, MP & RF Circuits Basic principles 11 I Basic principles 1 Basic principles of inductive components Magnetism The basis for understanding inductors
More informationHomework 2 solutions
Homework 2 solutions Section 2.1: Ex 1,2,3,4,6,11; AP 1 (18 points; Ex 3, AP 1 graded, 4 pts each; 2 pts to try the others) 1. Determine if each function has a unique fixed point on the specified interval.
More information