How do I solve an electron density continuity equation for Earth's atmosphere?

Click For Summary
SUMMARY

This discussion focuses on solving the continuity equation for electron density in the E-region of Earth's atmosphere, specifically neglecting vertical transport and considering the ion production rate ceasing at dusk. The equation presented is a Riccati equation, expressed as dNe/dt = p(t) - αd Ne², where Ne represents electron density and αd is the recombination coefficient for ions. The user seeks guidance on integrating this equation with an initial density of approximately 1 x 10^9 m³. Corrections to the equation notation were also noted, emphasizing the correct use of αd.

PREREQUISITES
  • Understanding of Riccati equations
  • Familiarity with electron density concepts in atmospheric physics
  • Knowledge of ion production and recombination coefficients
  • Basic calculus for integration techniques
NEXT STEPS
  • Study the methods for solving Riccati equations
  • Research the impact of ion production rates on electron density
  • Explore numerical integration techniques for differential equations
  • Review atmospheric models related to the E-region
USEFUL FOR

Students and researchers in atmospheric physics, particularly those studying ionospheric processes and electron density dynamics in the E-region of Earth's atmosphere.

Imagin_e
Messages
58
Reaction score
0

Homework Statement


Hi!

I need to solve a continuity equation for electron density as a function of time in the E-region of Earth's atmosphere. I shall neglect vertical transport and that the ion production rate completely shuts off dusk.

Homework Equations


See below

The Attempt at a Solution


Here is my attempt:

(1+λ)dNe/dt=p(t)-(1+λ)(αd-aλi)Ne2-Ne dλ/dt
Ne is the electron density, dy, λ=ratio of negative ions to electrons , p(t) is the electron production rate, αd is the recombination coeff. for ions and ai is the recombination coeff. negative ions with positive ions.

The effective recombination coefficient is: αeff=(1+λ)(αd+aλi).
If we assume that the is no negative ions, we get:
dNe/dt=p(t)-adNe2 (1)
and aeff simply becomes ad

Here comes the issue, I need to have a solution for equation (1) with respect to the density and time. . How should I integrate this? The density from the beginning is ≈ 1*109 m3 , which I calculated
 
Physics news on Phys.org
There is a typo in equation ##(1)##. If I understand correctly you meant to write ##\alpha_d## not ##a_d##.

Equation ##(1)## has the form of a Riccati equation:

$$\frac{dN_e(t)}{dt} = p(t) - \alpha_d N_e^2(t)$$
$$y'(t) = p(t) + g(t)y(t) + f(t)y^2(t)$$

Where ##N_e(t) = y(t)##, ##p(t) = p(t)##, ##g(t) = 0##, and ##f(t) = \alpha_d##.

To learn how to deal with equations of this form, see: https://en.wikipedia.org/wiki/Riccati_equation#Reduction_to_a_second_order_linear_equation
 
Zondrina said:
There is a typo in equation ##(1)##. If I understand correctly you meant to write ##\alpha_d## not ##a_d##.

Equation ##(1)## has the form of a Riccati equation:

$$\frac{dN_e(t)}{dt} = p(t) - \alpha_d N_e^2(t)$$
$$y'(t) = p(t) + g(t)y(t) + f(t)y^2(t)$$

Where ##N_e(t) = y(t)##, ##p(t) = p(t)##, ##g(t) = 0##, and ##f(t) = \alpha_d##.

To learn how to deal with equations of this form, see: https://en.wikipedia.org/wiki/Riccati_equation#Reduction_to_a_second_order_linear_equation
THANK YOU THANK YOU THANK YOU!
 

Similar threads

  • · Replies 5 ·
Replies
5
Views
779
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 5 ·
Replies
5
Views
2K
Replies
19
Views
2K
  • · Replies 10 ·
Replies
10
Views
6K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
Replies
5
Views
7K
  • · Replies 1 ·
Replies
1
Views
5K
  • · Replies 3 ·
Replies
3
Views
3K