The solar atmosphere is dominated by loops of magnetic flux which connect the multi-million-degree corona to the much cooler chromosphere. The temperature and density structure of quasi-static loops is determined by the continuous flow of energy from the hot corona to the lower solar atmosphere. Loop scaling laws provide relationships between global properties of the loop (such as peak temperature, pressure, and length); they follow from the physical variable dependencies of various terms in the energy equation, and hence the form of the loop scaling law provides insight into the key physics that controls the loop structure. Traditionally, scaling laws have been derived under the assumption of collision-dominated thermal conduction. Here we examine the impact of different regimes of thermal conduction -- collision-dominated, turbulence-dominated, and free-streaming -- on the form of the scaling laws relating the loop temperature and heating rate to its pressure and half-length. We show that the scaling laws for turbulence-dominated conduction are fundamentally different than those for collision-dominated and free-streaming conduction, inasmuch as the form of the scaling laws now depend primarily on conditions at the low-temperature, rather than high-temperature, part of the loop. We also establish regimes in temperature and density space in which each of the applicable scaling laws prevail.

Electromagnetic wave scattering off density inhomogeneities in the solar corona is an important process which determines both the apparent source size and the time profile of radio bursts observed at 1 AU. Here we model the scattering process using a Fokker-Planck equation and apply this formalism to several regimes of interest. In the first regime the density fluctuations are considered quasi-static and diffusion in wavevector space is dominated by angular diffusion on the surface of a constant energy sphere. In the small-angle ("pencil beam") approximation, this diffusion further occurs over a small solid angle in wavevector space. The second regime corresponds to a much later time, by which scattering has rendered the photon distribution near-isotropic resulting in a spatial diffusion of the radiation. The third regime involves time-dependent fluctuations and, therefore, Fermi acceleration of photons. Combined, these results provide a comprehensive theoretical framework within which to understand several important features of propagation of radio burst waves in the solar corona: emitted photons are accelerated in a relatively small inner region and then diffuse outwards to larger distances. En route, angular diffusion results both in source sizes which are substantially larger than the intrinsic source, and in observed intensity-versus-time profiles that are asymmetric, with a sharp rise and an exponential decay. Both of these features are consistent with observations of solar radio bursts.

A considerable fraction of the energy in a solar flare is released as suprathermal electrons; such electrons play a major role in energy deposition in the ambient atmosphere and hence the atmospheric response to flare heating. Historically the transport of these particles has been approximated through a deterministic approach in which first-order secular energy loss to electrons in the ambient target is treated as the dominant effect, with second-order diffusive terms (in both energy and angle) being generally either treated as a small correction or neglected. However, it has recently been pointed out that while neglect of diffusion in energy may indeed be negligible, diffusion in angle is of the same order as deterministic scattering and hence must be included. Here we therefore investigate the effect of angular scattering on the energy deposition profile in the flaring atmosphere. A relatively simple compact expression for the spatial distribution of energy deposition into the ambient plasma is presented and compared with the corresponding deterministic result. For unidirectional injection there is a significant shift in heating from the lower corona to the upper corona; this shift is much smaller for isotropic injection. We also compare the heating profiles due to return current Ohmic heating in the diffusional and deterministic models.

Solar radio observations provide a unique diagnostic of the outer solar atmosphere. However, the inhomogeneous turbulent corona strongly affects the propagation of the emitted radio waves, so decoupling the intrinsic properties of the emitting source from the effects of radio-wave propagation has long been a major challenge in solar physics. Here we report quantitative spatial and frequency characterization of solar radio burst fine structures observed with the LOw Frequency Array (LOFAR), an instrument with high time resolution that also permits imaging at scales much shorter than those corresponding to radio-wave propagation in the corona. The observations demonstrate that radio-wave propagation effects, and not the properties of the intrinsic emission source, dominate the observed spatial characteristics of radio burst images. These results permit more accurate estimates of source brightness temperatures, and open opportunities for quantitative study of the mechanisms that create the turbulent coronal medium through which the emitted radiation propagates.

Using the "enthalpy-based thermal evolution of loops" (EBTEL) model, we investigate the hydrodynamics of the plasma in a flaring coronal loop in which heat conduction is limited by turbulent scattering of the electrons that transport the thermal heat flux. The EBTEL equations are solved analytically in each of the two (conduction-dominated and radiation-dominated) cooling phases. Comparison of the results with typical observed cooling times in solar flares shows that the turbulent mean free-path $\lambda_T$ lies in a range corresponding to a regime in which classical (collision-dominated) conduction plays at most a limited role. We also consider the magnitude and duration of the heat input that is necessary to account for the enhanced values of temperature and density at the beginning of the cooling phase and for the observed cooling times. We find through numerical modeling that in order to produce a peak temperature $\simeq 1.5 \times 10^7$~K and a 200~s cooling time consistent with observations, the flare heating profile must extend over a significant period of time; in particular, its lingering role must be taken into consideration in any description of the cooling phase. Comparison with observationally-inferred values of post-flare loop temperatures, densities, and cooling times thus leads to useful constraints on both the magnitude and duration of the magnetic energy release in the loop, as well as on the value of the turbulent mean free-path $\lambda_T$.

We study the influence of kinetic-scale Alfv\'enic turbulence on the generation of plasma radio emission in the solar coronal regions where the plasma/magnetic pressure ratio $\beta $ is smaller than the electron/ion mass ratio $m_{e}/m_{i}$. The present study is motivated by the phenomenon of solar type I radio storms associated with the strong magnetic field of active regions. The measured brightness temperature of the type I storms can be up to $10^{10}$ K for continuum emission, and can exceed $10^{11}$ K for type I bursts. At present, there is no generally accepted theory explaining such high brightness temperatures and some other properties of the type I storms. We propose the model with the imbalanced turbulence of kinetic-scale Alfv\'en waves producing an asymmetric quasilinear plateau on the upward half of the electron velocity distribution. The Landau damping of resonant Langmuir waves is suppressed and their amplitudes grow spontaneously above the thermal level. The estimated saturation level of Langmuir waves is high enough to generate observed type I radio emission at the fundamental plasma frequency. Harmonic emission does not appear in our model because the backward-propagating Langmuir waves undergo a strong Landau damping. Our model predicts $100\%$ polarization in the sense of the ordinary (o-) mode of type I emission.

The energy released in solar flares derives from a reconfiguration of magnetic fields to a lower energy state, and is manifested in several forms, including bulk kinetic energy of the coronal mass ejection, acceleration of electrons and ions, and enhanced thermal energy that is ultimately radiated away across the electromagnetic spectrum from optical to X-rays. Using an unprecedented set of coordinated observations, from a suite of instruments, we here report on a hitherto largely overlooked energy component -- the kinetic energy associated with small-scale turbulent mass motions. We show that the spatial location of, and timing of the peak in, turbulent kinetic energy together provide persuasive evidence that turbulent energy may play a key role in the transfer of energy in solar flares. Although the kinetic energy of turbulent motions accounts, at any given time, for only $\sim (0.5-1)$\% of the energy released, its relatively rapid ($\sim$$1-10$~s) energization and dissipation causes the associated throughput of energy (i.e., power) to rival that of major components of the released energy in solar flares, and thus presumably in other astrophysical acceleration sites.

The transport of the energy contained in suprathermal electrons in solar flares plays a key role in our understanding of many aspects of flare physics, from the spatial distributions of hard X-ray emission and energy deposition in the ambient atmosphere to global energetics. Historically the transport of these particles has been largely treated through a deterministic approach, in which first-order secular energy loss to electrons in the ambient target is treated as the dominant effect, with second-order diffusive terms (in both energy and angle) being generally either treated as a small correction or even neglected. We here critically analyze this approach, and we show that spatial diffusion through pitch-angle scattering necessarily plays a very significant role in the transport of electrons. We further show that a satisfactory treatment of the diffusion process requires consideration of non-local effects, so that the electron flux depends not just on the local gradient of the electron distribution function but on the value of this gradient within an extended region encompassing a significant fraction of a mean free path. Our analysis applies generally to pitch-angle scattering by a variety of mechanisms, from Coulomb collisions to turbulent scattering. We further show that the spatial transport of electrons along the magnetic field of a flaring loop can be modeled rather effectively as a Continuous Time Random Walk with velocity-dependent probability distribution functions of jump sizes and occurrences, both of which can be expressed in terms of the scattering mean free path.

We investigate the impact of turbulent suppression of parallel heat conduction on the cooling of post-flare coronal loops. Depending on the value of the mean free path $\lambda_T$ associated with the turbulent scattering process, we identify four main cooling scenarios. The overall temperature evolution, from an initial temperature in excess of $10^7$~K, is modeled in each case, highlighting the evolution of the dominant cooling mechanism throughout the cooling process. Comparison with observed cooling times allows the value of $\lambda_T$ to be constrained, and interestingly this range corresponds to situations where collision-dominated conduction plays a very limited role, or even no role at all, in the cooling of post-flare coronal loops.

The transport of the energy contained in electrons, both thermal and suprathermal, in solar flares plays a key role in our understanding of many aspects of the flare phenomenon, from the spatial distribution of hard X-ray emission to global energetics. Motivated by recent {\em RHESSI} observations that point to the existence of a mechanism that confines electrons to the coronal parts of flare loops more effectively than Coulomb collisions, we here consider the impact of pitch-angle scattering off turbulent magnetic fluctuations on the parallel transport of electrons in flaring coronal loops. It is shown that the presence of such a scattering mechanism in addition to Coulomb collisional scattering can significantly reduce the parallel thermal and electrical conductivities relative to their collisional values. We provide illustrative expressions for the resulting thermoelectric coefficients that relate the thermal flux and electrical current density to the temperature gradient and the applied electric field. We then evaluate the effect of these modified transport coefficients on the flare coronal temperature that can be attained, on the post-impulsive-phase cooling of heated coronal plasma, and on the importance of the beam-neutralizing return current on both ambient heating and the energy loss rate of accelerated electrons. We also discuss the possible ways in which anomalous transport processes have an impact on the required overall energy associated with accelerated electrons in solar flares.

X-ray observations are a direct diagnostic of fast electrons produced in solar flares, energized during the energy release process and directed towards the Sun. Since the properties of accelerated electrons can be substantially changed during their transport and interaction with the background plasma, a model must ultimately be applied to X-ray observations in order to understand the mechanism responsible for their acceleration. A cold thick target model is ubiquitously used for this task, since it provides a simple analytic relationship between the accelerated electron spectrum and the emitting electron spectrum in the X-ray source, with the latter quantity readily obtained from X-ray observations. However, such a model is inappropriate for the majority of solar flares in which the electrons propagate in a hot megaKelvin plasma, because it does not take into account the physics of thermalization of fast electrons. The use of a more realistic model, properly accounting for the properties of the background plasma, and the collisional diffusion and thermalization of electrons, can alleviate or even remove many of the traditional problems associated with the cold thick target model and the deduction of the accelerated electron spectrum from X-ray spectroscopy, such as the number problem and the need to impose an ad hoc low energy cut-off.

Extending previous studies of nonthermal electron transport in solar flares which include the effects of collisional energy diffusion and thermalization of fast electrons, we present an analytic method to infer more accurate estimates of the accelerated electron spectrum in solar flares from observations of the hard X-ray spectrum. Unlike for the standard cold-target model, the spatial characteristics of the flaring region, especially the necessity to consider a finite volume of hot plasma in the source, need to be taken into account in order to correctly obtain the injected electron spectrum from the source-integrated electron flux spectrum (a quantity straightforwardly obtained from hard X-ray observations). We show that the effect of electron thermalization can be significant enough to nullify the need to introduce an {\it ad hoc} low-energy cutoff to the injected electron spectrum in order to keep the injected power in non-thermal electrons at a reasonable value. Rather the suppression of the inferred low-energy end of the injected spectrum compared to that deduced from a cold-target analysis allows the inference from hard X-ray observations of a more realistic energy in injected non-thermal electrons in solar flares.

Twisted coronal loops should be ubiquitous in the solar corona. Twisted magnetic fields contain excess magnetic energy, which can be released during magnetic reconnection, causing solar flares. The aim of this work is to investigate magnetic reconnection, and particle acceleration and transport in kink-unstable twisted coronal loops, with a focus on the effects of resistivity, loop geometry and atmospheric stratification. Another aim is to perform forward-modelling of bremsstrahlung emission and determine the structure of hard X-ray sources. We use a combination of magnetohydrodynamic (MHD) and test-particle methods. First, the evolution of the kinking coronal loop is considered using resistive MHD model, incorporating atmospheric stratification and loop curvature. Then, the obtained electric and magnetic fields and density distributions are used to calculate electron and proton trajectories using a guiding-centre approximation, taking into account Coulomb collisions. It is shown that electric fields in twisted coronal loops can effectively accelerate protons and electrons to energies up to 10 MeV. High-energy particles have hard, nearly power-law energy spectra. The volume occupied by high-energy particles demonstrates radial expansion, which results in the expansion of the visible hard X-ray loop and a gradual increase in hard X-ray footpoint area. Synthesised hard X-ray emission reveals strong footpoint sources and the extended coronal source, whose intensity strongly depends on the coronal loop density.

We study a model of particle acceleration coupled with an MHD model of magnetic reconnection in unstable twisted coronal loops. The kink instability leads to the formation of helical currents with strong parallel electric fields resulting in electron acceleration. The motion of electrons in the electric and magnetic fields of the reconnecting loop is investigated using a test-particle approach taking into account collisional scattering. We discuss the effects of Coulomb collisions and magnetic convergence near loop footpoints on the spatial distribution and energy spectra of high-energy electron populations and possible implications on the hard X-ray emission in solar flares.

Driven by recent RHESSI observations of confined loop-top hard X-ray sources in solar flares, we consider stochastic acceleration of electrons in the presence of Coulomb collisions. If electron escape from the acceleration region can be neglected, the electron distribution function is determined by a balance between diffusive acceleration and collisions. Such a scenario admits a stationary solution for the electron distribution function that takes the form of a kappa distribution. We show that the evolution toward this kappa distribution involves a "wave front" propagating forwards in velocity space, so that electrons of higher energy are accelerated later; the acceleration time scales with energy according to $\tau_{\rm acc} \sim E^{3/2}$. At sufficiently high energies escape from the finite-length acceleration region will eventually dominate. For such energies, the electron velocity distribution function is obtained by solving a time-dependent Fokker-Planck equation in the "leaky-box" approximation. Solutions are obtained in the limit of a small escape rate from an acceleration region that can effectively be considered a thick target.

Observations with {\em RHESSI} have enabled the detailed study of the structure of dense hard X-ray coronal sources in solar flares. The variation of source extent with electron energy has been discussed in the context of streaming of non-thermal particles in a one-dimensional cold-target model, and the results used to constrain both the physical extent of, and density within, the electron acceleration region. Here we extend this investigation to a more physically realistic model of electron transport that takes into account the finite temperature of the ambient plasma, the initial pitch-angle distribution of the accelerated electrons, and the effects of collisional pitch-angle scattering. The finite temperature results in the thermal diffusion of electrons, that leads to the observationally-inferred value of the acceleration region volume being an overestimate of its true value. The different directions of the electron trajectories, a consequence of both the non-zero injection pitch-angle and scattering within the target, cause the projected propagation distance parallel to the guiding magnetic field to be reduced, so that a one-dimensional interpretation can overestimate the actual density by a factor of up to $\sim 6$. The implications of these results for the determination of acceleration region properties (specific acceleration rate, filling factor, etc.) are discussed.

Recent observations from {\em RHESSI} have revealed that the number of non-thermal electrons in the coronal part of a flaring loop can exceed the number of electrons required to explain the hard X-ray-emitting footpoints of the same flaring loop. Such sources cannot, therefore, be interpreted on the basis of the standard collisional transport model, in which electrons stream along the loop while losing their energy through collisions with the ambient plasma; additional physical processes, to either trap or scatter the energetic electrons, are required. Motivated by this and other observations that suggest that high energy electrons are confined to the coronal region of the source, we consider turbulent pitch angle scattering of fast electrons off low frequency magnetic fluctuations as a confinement mechanism, modeled as a spatial diffusion parallel to the mean magnetic field. In general, turbulent scattering leads to a reduction of the collisional stopping distance of non-thermal electrons along the loop and hence to an enhancement of the coronal HXR source relative to the footpoints. The variation of source size $L$ with electron energy $E$ becomes weaker than the quadratic behavior pertinent to collisional transport, with the slope of $L(E)$ depending directly on the mean free path $\lambda$ again pitch angle scattering. Comparing the predictions of the model with observations, we find that $\lambda \sim$$(10^8-10^9)$ cm for $\sim30$ keV, less than the length of a typical flaring loop and smaller than, or comparable to, the size of the electron acceleration region.

The acceleration of charged particles in magnetized plasmas is considered during turbulent multi-island magnetic reconnection. The particle acceleration model is constructed for an ensemble of islands which produce adiabatic compression of the particles. The model takes into account the statistical fluctuations in the compression rate experienced by the particles during their transport in the acceleration region. The evolution of the particle distribution function is described as a simultaneous first and second-order Fermi acceleration process. While the efficiency of the first-order process is controlled by the average rate of compression, the second order process involves the variance in the compression rate. Moreover, the acceleration efficiency associated with the second-order process involves both the Eulerian properties of the compression field and the Lagrangian properties of the particles. The stochastic contribution to the acceleration is non-resonant and can dominate the systematic part in the case of a large variance in the compression rate. The model addresses the role of the second-order process, how the latter can be related to the large-scale turbulent transport of particles and explains some features of the numerical simulations of particle acceleration by multi-island contraction during magnetic reconnection.

Non-thermal electron populations are observed throughout the heliosphere. The relaxation of an electron beam is known to produce Langmuir waves which, in turn, may substantially modify the electron distribution function. As the Langmuir waves are refracted by background density gradients and as the solar and heliospheric plasma density is naturally perturbed with various levels of inhomogeneity, the interaction of Langmuir waves with non-thermal electrons in inhomogeneous plasmas is an important topic. We investigate the role played by ambient density fluctuations on the beam-plasma relaxation, focusing on the effect of acceleration of beam electrons. The scattering of Langmuir waves off turbulent density fluctuations is modeled as a wavenumber diffusion process which is implemented in numerical simulations of the one-dimensional quasilinear kinetic equations describing the beam relaxation. The results show that a substantial number of beam electrons are accelerated when the diffusive time scale in wavenumber space \tau_D is of the order of the quasilinear time scale \tau_ql, while when \tau_D << \tau_ql, the beam relaxation is suppressed. Plasma inhomogeneities are therefore an important means of energy redistribution for waves and hence electrons, and so must be taken into account when interpreting, for example, hard X-ray or Type III emission from flare-accelerated electrons.

We establish a classification scheme for stochastic acceleration models involving low-frequency plasma turbulence in a strongly magnetized plasma. This classification takes into account both the properties of the accelerating electromagnetic field, and the nature of the transport of charged particles in the acceleration region. We group the acceleration processes as either resonant, non-resonant or resonant-broadened, depending on whether the particle motion is free-streaming along the magnetic field, diffusive or a combination of the two. Stochastic acceleration by moving magnetic mirrors and adiabatic compressions are addressed as illustrative examples. We obtain expressions for the momentum-dependent diffusion coefficient $D(p)$, both for general forms of the accelerating force and for the situation when the electromagnetic force is wave-like, with a specified dispersion relation $\omega=\omega(k)$. Finally, for models considered, we calculate the energy-dependent acceleration time, a quantity that can be directly compared with observations of the time profile of the radiation field produced by the accelerated particles, such as during solar flares.

Context. High energy electrons accelerated during solar flare are abundant in the solar corona and in the interplanetary space. Commonly, the number and the energy of non-thermal electrons at the Sun is estimated using hard X-ray (HXR) spectral observations (e.g. RHESSI) and a single-particle collisional approximation. Aims. To investigate the role of the spectrally evolving Langmuir turbulence on the population of energetic electrons in the solar corona. Methods. We numerically simulate the relaxation of a power-law non-thermal electron population in a collisional inhomogeneous plasma including wave-particle, and wave-wave interactions. Results. The numerical simulations show that the long-time evolution of electron population above 20 keV deviates substantially from the collisional approximation when wave-particle interactions in non-uniform plasma are taken into account. The evolution of Langmuir wave spectrum towards smaller wavenumbers, due to large-scale density fluctuations and wave-wave interactions, leads to an effective acceleration of electrons. Furthermore, the time-integrated spectrum of non-thermal electrons, which is normally observed with HXR above 20 keV, is noticeably increased due to acceleration of non-thermal electrons by Langmuir waves. Conclusions. The results show that the observed HXR spectrum, when interpreted in terms of collisional relaxation, can lead to an overestimated number and energy of energetic electrons accelerated in the corona.

A fundamental problem in astrophysics is the interaction between magnetic turbulence and charged particles. It is now possible to use \emph{Ramaty High Energy Solar Spectroscopic Imager (RHESSI)} observations of hard X-rays (HXR) emitted by electrons to identify the presence of turbulence and to estimate the magnitude of the magnetic field line diffusion coefficient at least in dense coronal flaring loops.} {We discuss the various possible regimes of cross-field transport of non-thermal electrons resulting from broadband magnetic turbulence in coronal loops. The importance of the Kubo number $K$ as a governing parameter is emphasized and results applicable in both the large and small Kubo number limits are collected.} {Generic models, based on concepts and insights developed in the statistical theory of transport, are applied to the coronal loops and to the interpretation of hard X-ray imaging data in solar flares. The role of trapping effects, which become important in the non-linear regime of transport, is taken into account in the interpretation of the data.} For this flaring solar loop, we constrain the ranges of parallel and perpendicular correlation lengths of turbulent magnetic fields and possible Kubo numbers. We show that a substantial amount of magnetic fluctuations with energy $\sim 1%$ (or more) of the background field can be inferred from the measurements of the magnetic diffusion coefficient inside thick-target coronal loops.

Plasma turbulence is thought to be associated with various physical processes involved in solar flares, including magnetic reconnection, particle acceleration and transport. Using Ramaty High Energy Solar Spectroscopic Imager ({\it RHESSI}) observations and the X-ray visibility analysis, we determine the spatial and spectral distributions of energetic electrons for a flare (GOES M3.7 class, April 14, 2002 23$:$55 UT), which was previously found to be consistent with a reconnection scenario. It is demonstrated that because of the high density plasma in the loop, electrons have to be continuously accelerated about the loop apex of length $\sim 2\times 10^9$cm and width $\sim 7\times 10^8$cm. Energy dependent transport of tens of keV electrons is observed to occur both along and across the guiding magnetic field of the loop. We show that the cross-field transport is consistent with the presence of magnetic turbulence in the loop, where electrons are accelerated, and estimate the magnitude of the field line diffusion coefficient for different phases of the flare. The energy density of magnetic fluctuations is calculated for given magnetic field correlation lengths and is larger than the energy density of the non-thermal electrons. The level of magnetic fluctuations peaks when the largest number of electrons is accelerated and is below detectability or absent at the decay phase. These hard X-ray observations provide the first observational evidence that magnetic turbulence governs the evolution of energetic electrons in a dense flaring loop and is suggestive of their turbulent acceleration.

Previous numerical studies have identified "phase mixing" of low-frequency Alfven waves as a mean of parallel electric field amplification and acceleration of electrons in a collisionless plasma. Theoretical explanations are given of how this produces an amplification of the parallel electric field, and as a consequence, also leads to enhanced collisionless damping of the wave by energy transfer to the electrons. Our results are based on the properties of the Alfven waves in a warm plasma which are obtained from drift-kinetic theory, in particular, the rate of their electron Landau damping. Phase mixing in a collisionless low-$\beta$ plasma proceeds in a manner very similar to the visco-resistive case, except for the fact that electron Landau damping is the primary energy dissipation channel. The time and length scales involved are evaluated. We also focus on the evolution of the parallel electric field and calculate its maximum value in the course of its amplification.

Anisotropic Alfv\'{e}nic fluctuations with $k_{\parallel}/k_{\perp}\ll 1$ remain at frequencies much smaller than the ion cyclotron frequency in the presence of a strong background magnetic field. Based on the simplest truncation of the electromagnetic gyrofluid equations in a homogeneous plasma, a model for the energy cascade produced by Alfv\'{e}nic turbulence is constructed, which smoothly connect the large magnetohydrodynamics (MHD) scales and the small "kinetic" scales. Scaling relations are obtained for the electromagnetic fluctuations, as a function of $k_{\perp}$ and $k_{\parallel}$. Moreover, a particular attention is paid to the spectral structure of the parallel electric field which is produced by Alfv\'{e}nic turbulence. The reason is the potential implication of this parallel electric field in turbulent acceleration and transport of particles. For electromagnetic turbulence, this issue was raised some time ago in [A. Hasegawa, K. Mima, J. Geophys. Res. {\bf 83} 1117 (1978)].

{This work aims to investigate the spectral structure of the parallel electric field generated by strong anisotropic and balanced Alfvenic turbulence in relation with the problem of electron acceleration from the thermal population in solar flare plasma conditions.} {We consider anisotropic Alfvenic fluctuations in the presence of a strong background magnetic field. Exploiting this anisotropy, a set of reduced equations governing non-linear, two-fluid plasma dynamics is derived. The low-$\beta$ limit of this model is used to follow the turbulent cascade of the energy resulting from the non-linear interaction between kinetic Alfven waves, from the large magnetohydrodynamics (MHD) scales with $k_{\perp}\rho_{s}\ll 1$ down to the small "kinetic" scales with $k_{\perp}\rho_{s} \gg 1$, $\rho_{s}$ being the ion sound gyroradius.} {Scaling relations are obtained for the magnitude of the turbulent electromagnetic fluctuations, as a function of $k_{\perp}$ and $k_{\parallel}$, showing that the electric field develops a component parallel to the magnetic field at large MHD scales.} {The spectrum we derive for the parallel electric field fluctuations can be effectively used to model stochastic resonant acceleration and heating of electrons by Alfven waves in solar flare plasma conditions}

The propagation of shear Alfven wave packets in inhomogeneous magnetic fields, at the origin of their distortion, regardless of the occurrence of non-linear coupling, is considered. It is shown that the distortion mechanism can be regarded as mixing process and hence, standard "phase mixing" corresponds to the effect of an "Alfvenic" shear flow while enhanced dissipation at a magnetic X-point corresponds to mixing by an "Alfvenic" strain flow. The evolution of the global wave field is supposed to result from the dynamics of a superposition of wave packets and a kinetic equation for the wave energy is obtained following this eikonal (WKB) description. Since shear Alfven wave packets experience continuous shearing/straining while transported by an inhomogeneous Alfvenic flow $\mathbf{V}_{A}$, their mixing process, in physical space, is also a cascade of wave energy in k-space. The wave energy spectrum resulting from this linear mechanism of energy transfer is determined for the special case of waves propagating along chaotic magnetic field lines, the analog of a chaotic mixing process. The latter follows a $k^{-1}$ power-law, in the energy conserving range in $k$ space.