How can I calculate in MATLAB-the electric field due to a time varying magentic field for a simple geometry?
23 次查看(过去 30 天)
显示 更早的评论
I have a problem of electromagentic modeling. The problem involves a current carrying circular coil where the current will be varied with time. Calculating the magnetic field was relatively simpler - I used Bio-Savart's law & the superposition theory. But I am facing some difficulty calculating the electric field that arises from the time varying magnetic field, specifically, I didn't know how to implement Maxwell's equation numerically. I would appreciate any help in any form, however, after much speculation, I finished writing my first script of calculating the electric field. I amm attaching the vector plot of the electromagnetic field distribution in a 3D volume. I am also attaching the .fig file of the corresponding plot so anyone interested can maneuver plot as they like. I would like to know if the electric field I calculated is correct.
2 个评论
G A
2021-8-11
Will the FDTD 3D method help you? See e.g.
https://www.mathworks.com/matlabcentral/fileexchange/54557-3d-finite-different-time-domain-first-order-mur-boundary-condition
A.Martusciello
2024-6-21
I calculated magnetic induction of multilyer coil with Biot-Savart, at the moment using Ieff current, after with a set of waveform each at a time. So i have three matrix for Hx, Hy, Hz, and, obviously H(module). My question is relative to use antenna toolbox to calculate and plot E. There is some way to di this? Mi final purpose is the quantization of potenzial on 1mm^2 of bone cells to verify if it is big enough to activate hosteogenesys process, in particulare inside vertebral body. The coil is on the back of patient.
Thank a lot
回答(2 个)
Chandra Amma Palanisamy
2024-5-24
Hi Md Chowdhury,
MATLAB Antenna Toolbox can be used to create circular loop and compute E and H fields. The below code can be used to design circular loop at 1 GHz and calculate the EHfields.
% design a loop circular
l = design(loopCircular,1000e6);
% visualize
figure;show(l);
% specify the Points to compute E and H fields
x = [-0.08 -0.07 -0.06 -0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08];
y = [-0.08 -0.07 -0.06 -0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08];
z = [-0.05 -0.04 -0.03 -0.02 -0.01 0.01 0.02 0.03 0.04 0.05];
[X,Y,Z] = meshgrid(x,y,z);
points = [X(:)';Y(:)';Z(:)'];
% visualize the E & H fields
figure; EHfields(l,1000e6,points,"ScaleFields",[5 5]);
To acccurately model the coil, convert the loopCircular to wireStack
w = wireStack(l);
figure; EHfields(w,1000e6,points,"ScaleFields",[5 5]);
0 个评论
David Goodmanson
2021-8-12
编辑:David Goodmanson
2021-8-12
Hello MC,
From the plot it is not that easy to see what is going on since the two field colors are similar. Could you post a separate plot for each field?
If you can successfully produce the B field by using a Biot-Savart integration then you can produce an E field in the same way, especially since it's easier. Here I use the common assumptions [1] the current is changing slowly, so that there is a negligible amount of change to the system in the time it takes for light to propagate across the loop. Therefore retarded time effects can be neglected, [2] the current has the same value all the way around the loop so that there is no resulting charge density that can be the source of an electric field. In that case the B field is
B = (mu0/4pi) Integral J(r',t) x (r-r')/|r-r'|^3 d^3r'
where r,r' and J are vector quantities and (r-r') extends from the source point to the observation point. If the current is confined to a small filament with constant current I, then
B = (mu0/4pi) I Integral ds' x (r-r')/|r-r'|^3
where vector ds' is a line element that is everywhere tangent to the loop. I assume you are doing an integral like this to find B.
The electric field in this circumstance is
E = -(mu0/4pi) Integral 1/|r-r'| dJ(r',t)/dt d^3r'
and
E = -(mu0/4pi) dI/dt Integral ds'/|r-r'|
For an oscillation at a single frequency, described by e^(-iwt), then d/dt can be replaced by -iw in which case
E = (mu0/4pi) iw Int 1/|r-r'| J(r') d^3r'
and
E = (mu0/4pi) iwI Int ds'/|r-r'|
The factor of i just means that B and E are 90 degrees out of phase as functions of time.
The E field lines go in circles and look similar to the current loop.
0 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Solvers 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!