Physics
Scientific paper
Dec 2004
adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=2004agufmsm53b0420a&link_type=abstract
American Geophysical Union, Fall Meeting 2004, abstract #SM53B-0420
Physics
7827 Kinetic And Mhd Theory, 7835 Magnetic Reconnection, 7843 Numerical Simulation Studies, 2753 Numerical Modeling
Scientific paper
We've developed an algorithm to iteratively solve the linearized Vlasov equation for growth rates and first-order potentials of instabilities such as collisionless ion-electron tearing in a neutral sheet. It requires no assumptions about the nature of the equilibrium trajectories and can give exact (numerical) results. In this work we describe the algorithm, and its application to study the behavior of the tearing-mode starting from a Harris equilibrium over a range of the control parameters: the current sheet half-width w, (Te)/(T_i), (Me)/(M_i). The method starts with computing and saving several thousand equilibrium particle paths to sample the phase space for a chosen equilibrium state. Then for a given k, the coupled equations for ěc{A1} and γ are solved iteratively with ěc{J1} computed from time-integrals over the saved particle paths. A new ěc{A1} is found using a Green's function to ensure the boundary- conditions are satisfied, while the correction to γ is found from a relation involving integrals of the currents and ěc{A1}. φ 1 is found by assuming quasi-neutrality. The method is stable, converges for any starting values we've tried, requires less than a MB of memory, and 10MB--10GB of disk space for the paths. CPU time on a desktop PC is from ˜ 10 minutes to an hour for each k. The algorithm is similar in spirit to previous work [1], but was developed independently so the implementation is different. Typical dispersion curves, plots of potentials and response functions will be shown. We've found the growth-rates can be fit quite well (correlations of 0.999) by the expression γ (k)= ax kb(1/w-k) where parameters ax, b are determined from the peak growth γ x and wavenumber kx. Preliminary results are that for parameter values (Me)/(M_i)=1, (Te)/(T_i)=1 the dependence of γ x, kx on w, over the range 3--20, is given by kx=0.527 w-1.07 and γ x=0.219 w-2.32 where γ x is scaled by Ω i, k and w are scaled by ρ i. For (Me)/(M_i)= 1/1836, (Te)/(T_i)=1/2 the dependence is kx=0.451 w-1.11 and γ x=0.0567 w-2.379. [1] M. Brittnacher, K. B. Quest, and H. Karimabadi, J. Geophys. Res., 100, 3551, (1995); W. Daughton, J. Geophys. Res., 103, 29429, (1998).
No associations
LandOfFree
Using Exact Particle Paths in the Calculation of Linear Ion-Electron Tearing and Dependence on the Equilibrium State Parameters does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with Using Exact Particle Paths in the Calculation of Linear Ion-Electron Tearing and Dependence on the Equilibrium State Parameters, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Using Exact Particle Paths in the Calculation of Linear Ion-Electron Tearing and Dependence on the Equilibrium State Parameters will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-1470727