Stack Exchange Network

Stack Exchange network consists of 175 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

Questions tagged [electromagnetism]

The tag has no usage guidance.

0
votes
1answer
29 views

TMZ TME modes, clarification

I'm refering here to Taflove's "computational electrodynamcis, 3rd ed." He says Let us assume that the structure being modeled extends to infinity in the z-direction with no change in the shape ...
2
votes
0answers
37 views

Open source multiphysics software that can model the interaction of two insulators with a user-specified electrical force law?

Given two solid bodies of masses $m_1$ and $m_2$ with invariable charge distributions $\rho_1(\mathbf{r_1})$ and $\rho_2(\mathbf{r_2})$ in them,* respectively, I would like to model the electrical ...
0
votes
0answers
41 views

How does one do electrodynamics simulations in SU2?

The description of SU2 says The primary applications are computational fluid dynamics and aerodynamic shape optimization, but has been extended to treat more general equations such as ...
1
vote
1answer
44 views

Average value divergence in spectral method for Poisson equation

I'd like to know how to deal with a divergence when trying to solve the Poisson equation for electrostatics with a simple spectral method. I'm not sure how to best state my problem, so I'll explain ...
2
votes
0answers
54 views

Solving electrostatics Poisson equation with Intel MKL routines

I am trying to solve the 3D Poisson equation $$\nabla\cdot(\epsilon(\mathbf{r})\nabla) u(\mathbf{r}) = f(\mathbf{r})$$ I notice intel advertises routines that appear to solve $$\nabla^2u(\mathbf{r}...
0
votes
1answer
85 views

Poincare map for Arnold-Beltrami-Childress Magnetic Field in Python

I want to plot the Poincare map for Arnold-Beltrami-Childress magnetic field for parameters $A=1, B=0.816, C=0.5773$ in Python for the Poincare section $z=0$. Also, I am not able to understand what ...
2
votes
1answer
59 views

Magnetization Vector from XY Model for an AntiFerromagnetic System

I am working on an XY model and I'm trying to calculate the magnetization and direction for an anti-ferromagnetic (AF) system. So I have a collection of spins in the $XY$ plane represented as vectors ...
2
votes
1answer
84 views

Calculation of the EFIE integral

I need help computing the following integral: $$ \int_{}\frac{(1+jk|\vec{r}-\vec{r}^\prime|)e^{-jk|\vec{r}-\vec{r}^\prime|}}{|\vec{r}-\vec{r}^\prime|}d\vec{r}^\prime $$ in this integral $\vec{r}$ ...
1
vote
0answers
40 views

Magnetostatic modelling Radia: Increasing distance between magnets a produces positive force until a certain point, beyond which force goes haywire

I have two sets of magnets. One set consists of two electromagnets (Shown below: Blue) and the other set consists of two NdFeB N40 permanent magnets. Both sets of magnets lie on the same plane. I want ...
0
votes
1answer
40 views

Finite difference time domain and dynamic permittivity

Since the permittivity of any material is usually complex function of temperature, frequency, density, etc. I was wondering if it is possible to use a dynamic permittivity which changes as a function ...
0
votes
1answer
67 views

Help understanding this numerical surface integration technique?

I'm attempting to write a FORTRAN program that calculates the magnetic field, B, at any point outside of a bar magnet. I'm going to use a first order euler scheme, where each side of the bar magnet ...
4
votes
0answers
126 views

How can I solve the wave equation for a circular rod in cylindrical coordinates using finite differences?

I have a problem with the stability of finite difference method for the wave equation in cylindrical coordinates. the equation is: $$ \frac{\partial^2 \omega_n}{\partial r^2}+\frac{1}{r}\frac{\...
1
vote
1answer
91 views

Combine Hydrodynamics and Electromagnetics

Is it possible, in general, to combine hydrodynamical motion and expansion of material with, say, a finite difference time domain method to simulate light-matter interaction? If so, how is this done ...
0
votes
1answer
120 views

Introducing EigenModes from 2D FEM into 3D FEM

This particular FEM question concerns waveguides and FEM 3D simulation. To excite a waveguide with waveport (TE10 and so on), we typically have to solve for eigenvalues ($k$) of helmholtz equation ...
6
votes
0answers
240 views

Understanding Boundary Condition in FEM

I am trying to understand Dirichlet and Neumann boundary conditions in FEM and I wanted to know if my inference is correct. To articulate my understanding, lets consider a simple case of TE and TM ...
1
vote
0answers
89 views

Band structure of nonlinear Schrodinger equation with one dimensional potential

I have a nonlinear Schrodinger equation which reads: $$ \frac{1}{2} \frac{d^2u}{dx^2}+ |u|^2u + V(x)u = -i \frac{du}{dz},$$ where $V(x)=\cos(wx)+ i a \sin(wx)$ and $w$, $a$ are numbers. How to ...
1
vote
0answers
72 views

Numerical scheme to solve Maxwell equations with fixed potential boundaries?

We have a 2D electromagnetic field (in the sense that: $E=(E_x,E_y,0)$, $B=(0,0,B_z)$, and all derivatives with respect to $z$ are $0$), and we are considering a system made up of two walls at $x=-b$ ...
4
votes
1answer
108 views

How to force potential boundary conditions in the Yee scheme for solving Maxwell's equations?

Assume that we have a 2D electromagnetic field (in the sense that: $E=(E_x,E_y,0)$, $B=(0,0,B_z)$, and all derivatives with respect to $z$ are $0$), and that we are considering a system made up of two ...
1
vote
2answers
106 views

Rule of thumb for time-step for solving Maxwell's Equation using 3D-FDTD?

Is there something like a rule of thumb for an adequate time-step size when solving Maxwell's equation for the interaction of light with matter? I guess a single wave oscillation has to be resolved ...
2
votes
1answer
119 views

How to deal with numerical errors in electrostatic field calculations

I want to trace electrostatic field lines emerging from 2D surfaces in 3D space. Eventually I want to find their intersection with an (uncharged) mesh. The charge distribution $\sigma(x), x \in \...
0
votes
1answer
37 views

What software packages are designed towards modelling the radiation from accelerated charges?

I'm interested in modeling the electromagnetic fields radiated from an accelerated charge, but do not want to reinvent stuff if possible. I suspect there are software packages already out there which ...
0
votes
0answers
64 views

How to calculate the CFL number for electrostatic charge transfer?

I am currently implementing the Robin Hood Algorithm for Electrostatics. It roughly works by searching the two triangles in a mesh, which have the biggest positive/negative deviation from an ...
1
vote
0answers
74 views

Can x-ray back-projection be converted to hard-field magnetic induction tomography?

This is a question about hard-field back-projection as used in x-ray tomography, applied magnetic induction tomography. Al-Zeibak and Saunders have shown that x-ray filtered backprojection can be ...
0
votes
1answer
77 views

Questions about implementing an electromagnetism/photonics solver package

I am hoping to start (very slowly) on implementing some form of a computational photonics/electromagnetism package. I know things like Meep, S4, FDTD++, EMPy, and a host of other proprietary/free/...
2
votes
1answer
62 views

Is this a proper implementation of point charge dynamics with ODEs

Since learning about point charges in my physics II class this semester, I want to be able to investigate not only the static force and field distributions but the actual trajectories of movement of ...
1
vote
1answer
629 views

Edge and Nodal finite element methods in MATLAB for Magnetic induction tomography

What is the difference between edge finite elements and nodal finite elements? This for use in modeling the eddy current problem in classical electromagnetism. I am attempting to convert MATLAB code ...
2
votes
1answer
105 views

Why FEM electric analysis gives only access to current density?

A Comsol study using frequency sweep on electric current physics yields only current density as accessible variables. I understand the underlying equation used is Ohm's law, i.e. $$\mathbf{J} = \...
1
vote
0answers
61 views

Dielectric in an external electric field

I want to numerically solve for electric field of an object in an external electric field. The object has linear displacement field such that its permitivity ($\epsilon$) is constant. The object's ...
1
vote
1answer
156 views

FEM or FD eigenvalue equation to get wave number instead of cutoff frequency

To get cutoff frequencies and eigenmode field distributions for a waveguide, one uses following equation: $$1/\epsilon ∇ \times 1/\mu ∇ \times E = \omega^2 E$$ With $ \omega^2 $ as eigenvalues. This ...
1
vote
1answer
54 views

Suitable method for simulation of in-fiber interferometer

I am trying to simulate an optic-fiber sensor (in-fiber interferometer) to study its respond to temperature. The method I am using is finite-difference time-domain (FDTD), and I come out with a large ...
1
vote
0answers
80 views

Representing a 3D system in 2D (Electromagnetic modelling)

Ok so I'm a complete beginner in computational modelling (I use analytical methods of physics typically) but I would like to model an anisotropic, aperiodic (but not random) finite array of metallic ...
0
votes
1answer
45 views

How to introduce a distance vector when decompose a force?

I read reference (1),but I am confused about how to introduce a distance vector in Lorentz force like the author does it in equation (11): $$\rm J\times B = -J\cdot \nabla \left( B\times r\right)+\...
7
votes
1answer
2k views

Simulate electric fields due to surface charges in simple circuits using python

I want to simulate the electric fields in simple circuits using Python and only free software. My first goal is to reproduce the images given in (1) which are made by the commercial ANSYS Maxwell ...
3
votes
1answer
281 views

Data analysis of a magnetic hysteresis loop

I am a physics undergrad and I have MOKE data for a magnetic material (thin permalloy film on a silicon substrate). Here is one of the hysteresis loops I obtained, plotted using Python: The form of ...
0
votes
2answers
3k views

Plot vector field in matlab

I have the function of an electric dipole expressed in cartesian coordinates and I want to create the vector field using Matlab . The function is $$E_z= \frac{p}{4\pi\epsilon_0} \cdot \left(\...
1
vote
0answers
65 views

Calculating tangential electric field intensity on the boundary

I have an object with surface charge $\sigma$ for which I want to caclulate tangential electric field which would correspond to mathematical formulation: $$ \vec{E}(\vec{x}) = \nabla_x \int \frac{\...
2
votes
1answer
268 views

I need to scale variables to solve a 2D PDE. What are the physical considerations of scaling?

I am solving a boundary value problem in 2D via an implicit finite difference scheme. Unfortunately, although the problem is well-posed and should have a unique solution, the condition number of the ...
3
votes
1answer
251 views

Finite Difference Beam Propagation Method problem

I am trying to implement the finite difference beam propagation method to study the propagation of a TE light signal through a waveguide. However, my solutions are exponentially growing, and display ...
5
votes
1answer
342 views

H(curl) conforming Nédélec-Elements to satisfy div(B)=0

Most authors are very clear that it's very dangerous to just use $\mathrm{H}(curl)$ conforming edge elements, which are divergence free, to satisfy $\mathrm{div}(\mathbf{B})=0$ and implement this ...
2
votes
2answers
1k views

How to solve a Poisson equation using the finite difference method when there is an object inside a domain?

I'm interested in solving an electrostatics problem in 2d case in some domain with a conductor placed inside the domain. From a mathematical point of view, I have to solve a Poisson equation with ...
1
vote
0answers
274 views

How to solve a nonlinear diffusion equation?

Consider a thin film with a perpendicular applied magnetic field $H_a$ in $z$-axis. The nonlocal relation between $H_a$, the self-field $H_\text{self}$ (generated by the eddy current $J$) and the ...
2
votes
2answers
114 views

Simulating interaction of Smart Fluids with Solids

Is it possible to simulate interaction of Smart Fluids with Solids? Is there a software capable of doing so? What I Know: The two software that I know of are Autodesk Simulation Mechanical and ...
4
votes
2answers
6k views

Visualising Maxwell's equations using MATLAB

I asked this question to help me understand what is going on in one of Maxwell's equations. I am happy with following through the maths on paper now, but would like to use MATLAB to take it one step ...
1
vote
2answers
321 views

Which type of meshing is more suited for simulation of electromagnetic metamaterial unit cells?

Electromagnetic metamaterials are resonant, periodic structures (repetition of a unit cell) involving both dielectric and conducting elements, and are used and simulated in frequency ranges from ...
1
vote
1answer
40 views

Is it necessary to use Poisson - Boltzmann equation if I only need to build electrostatic potential from a PQR file?

A PQR file of a protein describes the positions in x, y and z and the charges of every atoms in it. If you have them, you can construct the charge density $\rho(x)$ of the protein, then use it to find ...
1
vote
1answer
383 views

PQR files make from pdb2pqr are different to GROMACS

I have a pdb file and I would like to make a pqr one. In GROMACS, I use editconf -mead to produce the pqr file, with topol.top taken from pdb2gmx and use it to ...
3
votes
2answers
898 views

Tool for computing electric field on a 3d surface

I'm looking for software recommendations for an electrostatic simulation: I have a 3d object of uniform, high conductivity, separated from its environment by a high-resistance skin. Inside that ...
1
vote
0answers
1k views

Free HFSS alternative to create geometry using python

I am using a python script to create wirebond interfaces in Ansys HFSS. Variables define the board, wirebond, and package geometry. To simulate a new wirebond, I just edit the geometry variables and ...
1
vote
0answers
171 views

How can simulation be fast when Electric Field is being measured with probe in CST Microwave Studio? [closed]

I am measuring the electric field in the closure box with aperture and using probe to achieve that. But I have a problem. When I use the probe for measurement, results are not accurate in transient ...
8
votes
2answers
4k views

simple MHD simulation code for (self) education and play with

I would like some super simple computational code for solving magnetohydrodynamics problems. High accuracy nor performance is not my concern. I wan't it just to visually explore qualitative behavior ...