What is a Self Consistent Electric Field?

AI Thread Summary
A self-consistent electric field in one-dimensional plasma simulations is determined by the relationship between the electric field and the charge density, which influences each other. To solve for this field numerically, start with an initial guess for the charge distribution, compute the resulting electric field, and then update the charge distribution based on the new field. This iterative process continues until convergence is achieved, meaning the results stabilize to a desired accuracy. The parameter lambda, which balances the old and new charge distributions, must be chosen carefully during iterations. The specific approach may vary depending on the simulation context.
fys iks!
Messages
40
Reaction score
0
Hi,

Ive been doing some reading into 1 dimensional plasma numerical simulations and they keep referring to solving for a "self-consistent" field. If the simulation is in one dimension with periodic boundary conditions, how would I go about solving this electric field?

Example:

dE/dx = n - ρ(x)

where: n = const = 1
ρ(x) is the charge density and I want to solve for E numerically where E is "self consistent"


Thanks for your input.
 
Physics news on Phys.org
The electric field depends on a distribution of charges - but the distribution of charges depends on the electric field. This creates a chicken-and-egg situation.

A "self consistent" field is one which makes the charges distributed so that they generate the field. We can compute them using an iterative procedure.

You start with a guess for a charge distribution ρ0, compute the field that distribution gives rise to. That field will push the charges into a new configuration ρ' - so work out that new distribution as if the field were fixed at what you calculated.

Now repeat the procedure for ρ1=(1-λ)ρ0+λρ' where 0<λ<1.
You have to guess lambda.

Keep going until you keep getting the same result to the desired level of accuracy.

The exact method will depend on the context.
 
This is from Griffiths' Electrodynamics, 3rd edition, page 352. I am trying to calculate the divergence of the Maxwell stress tensor. The tensor is given as ##T_{ij} =\epsilon_0 (E_iE_j-\frac 1 2 \delta_{ij} E^2)+\frac 1 {\mu_0}(B_iB_j-\frac 1 2 \delta_{ij} B^2)##. To make things easier, I just want to focus on the part with the electrical field, i.e. I want to find the divergence of ##E_{ij}=E_iE_j-\frac 1 2 \delta_{ij}E^2##. In matrix form, this tensor should look like this...
Thread 'Applying the Gauss (1835) formula for force between 2 parallel DC currents'
Please can anyone either:- (1) point me to a derivation of the perpendicular force (Fy) between two very long parallel wires carrying steady currents utilising the formula of Gauss for the force F along the line r between 2 charges? Or alternatively (2) point out where I have gone wrong in my method? I am having problems with calculating the direction and magnitude of the force as expected from modern (Biot-Savart-Maxwell-Lorentz) formula. Here is my method and results so far:- This...
Back
Top