Orientador: Gilberto Medeiros Kremer / Tese (doutorado) - Universidade Federal do Paraná / Resumo: Ao longo deste trabalho tivemos três objetivos: (i) primeiramente determinamos as equações de Burnett em segunda ordem para um gás único relativístico com base na teoria cinética e no método dos momentos de Grad. Fizemos uma análise detalhada para todos os coeficientes da pressão dinâmica, fluxo de calor e do deviante do tensor pressão, expandindo estes coeficientes para os casos de temperaturas altas e baixas respectivamente. A seguir verificamos que os resultados obtidos para todos os coeficientes, no limite não-relativístico, concordam com os resultados obtidos por Wang Chang e Uhlenbeck; (ii) em seguida determinamos novamente as equações de Burnett relativísticas empregando agora o modelo cinético de Anderson e Witting e comparamos os resultados obtidos para os coeficientes de transporte com o método dos momentos de Grad. O modelo cinético de Anderson e Witting e uma boa aproximação para os resultados encontrados pelo método de Grad no caso ultra-relativístico; (iii) e finalmente, como uma aplicação dos resultados encontrados, fizemos uma análise para a propagação de ondas sonoras para os casos de cinco campos e quatorze campos. A seguir fizemos uma análise gráfica para a velocidade de fase e para a absorção das ondas sonoras analisando os limites não-relativístico e ultra-relativístico. / Abstract: In the course of this work we have three objectives: (i) firstly we obtained the relativistic second order Burnett equations for a single-component gas on the basis of kinetic theory and Grad’s moment method. We did a detailed analysis of the coefficients of dynamic pressure, heat flux and pressure deviator, making their expansions for the non-relativistic and ultra-relativistic cases. Then we verified that our results agree with Wang Chang and Uhlenbeck’s results, on the non-relativistic limit; (ii) next, we obtained the relativistic second order Burnett equations for a single-component gas on the basis of the kinetic model of Anderson and Witting and we compared the transport coefficients with results obtained by the Grad’s moment method. The kinetic model of Anderson and Witting is a good approximation for the Grad’s moment method in the ultra-relativistic limit. (iii) finally, as an application of our results, we analysed the propagation of sound in the case of five and fourteen fields. Plots of the phase velocity and absorption of sound are given together with their non-relativistic and ultra-relativistic limits.

O tokamak é um dos mais promissores modelos de reatores de fusão. O mesmo consiste em um sistema de confinamento magnético de forma toroidal, onde coexistem um campo magnético toroidal gerado por bobinas externas e um campo poloidal gerado por uma corrente que percorre o plasma confinado, formando um campo resultante helicoidal. A operação do tokamak é sujeita a várias perturbações e instabilidades. Os modos tearing neoclássicos (NTM) são perturbações perigosas que interferem na estrutura magnética do plasma, pela formação de estruturas chamadas ilhas magnéticas. Acredita-se que os NTM’s são gerados por uma queda na corrente de bootstrap dentro das ilhas magnéticas, onde o perfil de pressão fica achatado. Isso ocorre em superfícies de fluxo onde o fator de segurança é racional, que assim, pode ser escrito como q = m/n, onde o número poloidal m e o número toroidal n são ambos inteiros. O controle desses modos é essencial para a operação dos futuros reatores de fusão. Um dos métodos considerados para esse controle é a geração de calor e corrente localizada na região da ilha, usando ondas do tipo de cíclotron eletrônica, em um processo chamado electron cyclotron current drive (ECCD). Estudamos características importantes sobre a eficiência do processo, via um modelo simplificado do tokamak e da estrutura da ilha, acoplado a um tratamento auto-consistente da interação da ondas com o plasma no tokamak e, consequentemente, geração de corrente, considerando a existência de uma corrente inicial gerada por um transformador (loop voltage), e levando em conta efeitos de transporte radial de partículas e a existência de efeitos indutivos. Todos esses efeitos formam uma equação de Fokker-Planck para evolução temporal da função de distribuição eletrônica, que é resolvida numericamente. O formalismo também inclui uma descrição auto-consistente da largura da ilha magnética, como função da corrente no plasma, atraves de uma equação chamada equação de Rutherford modificada. O formalismo inclui também efeitos clássicos relacionados a um pulo da derivada do logarítmo do campo magnético radial, e dados basicamente pelo gradiente da corrente de equilíbrio e efeitos neoclássicos relacionados a perturbações na corrente de bootstrap. Os resultados são obtidos na análise numérica para diversos parâmetros iniciais e para ondas de cíclotron, e mostram que a largura das ilhas é reduzida pela geração de corrente gerada pelas ondas, o que confirma o ECCD como um bom método para estabilização dos NTM. / The tokamak is one of the most promising models for a fusion reactor. It consists in a system of magnetic confinement of toroidal form, in which coexist a toroidal magnetic field generated by a system of external coils and a poloidal magnetic field generated by the current in the confined plasma, forming a helicoidal magnetic field. The operation of a tokamak is subject to the occurrence of perturbations and instabilities. Neoclassical tearing modes (NTM) are dangerous perturbations which interfere in the magnetic structure of the plasma, by formation of structures known as magnetic islands. It is believed that the NTM are driven by a lack of the bootstrap current inside the magnetic island where the pressure profile is flattened. This occurs in flux surfaces with rational safety factor, which can be writen as q = m/n, where the poloidal number m and the toroidal number n are both integers. The control of these modes is essential to the operation of future fusion reactors. One of the methods considered for control is the heating and generation of current in the island location, using electron cyclotron (EC) waves, in a process known as electron cyclotron current drive (ECCD). We study important features related to the efficiency of such method by the use of a simplified model for a tokamak and for the structure of magnetic islands, coupled to a self-consistent treatment of the interaction of EC waves with the plasma in the tokamak and the consequent current generation, considering the occurrence of a loop voltage in the tokamak and taking into account radial transport of particles and the existence of induced effects. All of these effects are included in a general Fokker-Plank equation for the evolution of the electron distribution function, which is numerically solved. The formalism also includes a self-consistent description of the width of the magnetic islands, as a function of the current in the plasma, through the so-called modified Rutherford equation. The formalism includes classical effects are given by the jump in the logarithmic derivative of the radial magnetic field and driven mainly by the equilibrium current gradient, and neoclassical effects related to the effect of the perturbed bootstrap current. The results obtained in our numerical analyses, for several parameters associated to EC waves, show that the width of the islands can be substantially reduced by the EC generated currents, which seems to confirm ECCD as a promising method for NTM stabilization.

Neste trabalho estudamos a interação feixe-plasma aplicando a teoria de turbulência fraca. Primeiramente fazemos uma introdução `a teoria cinética de plasmas, estudando aspectos fundamentais como: a abordagem estatística das equações do plasma, o tratamento a ser dado `as funções de correlação que aparecem nessa abordagem, e o sistema de equações de Vlasov-Maxwell. Em seguida estudamos a aproximação de Vlasov e a solução do sistema de Vlasov- Maxwell na aproximação linear, enfocando a descrição de ondas no plasma e o amortecimento de Landau. Ao longo desse desenvolvimento fazemos uma apresentação relativamente detalhada dos procedimentos introduzidos por Landau ao tratar o sistema Vlasov-Maxwell como um problema de valor inicial, incluindo uma discussão sobre a resolução de integrais no plano complexo, com polos no denominador, e a solução da relação de dispersão para encontrar os modos normais de oscilação no plasma. Apresentamos também uma breve revisão a respeito da aproximação quase-linear do sistema de Vlasov- Maxwell, na qual abordamos a obtenção da equação quase-linear de difusão no espaço das velocidades para partículas e estudamos suas propriedades de conservação. Nesse contexto da teoria quase-linear, apresentamos uma revisão a respeito da evolução temporal do amortecimento de Landau no caso de distribuição de velocidades Maxwelliana e da instabilidade que pode ocorrer no caso de uma função distribuição com um feixe de partículas. Depois dessa introdução `a aproximação linear e `a teoria quase-linear, abordamos a teoria de turbulência fraca num plasma não magnetizado, incluindo interações não lineares entre ondas e partículas. O formalismo apresentado inclui efeitos como emissão espontânea e induzida, decaimento e espalhamento de ondas. O sistema de equações acopladas da teoria de turbulência fraca é então reduzido a uma aplicação a um sistema considerando duas dimensões, depois reescrito em termos de coordenadas polares e adaptado para solução numérica. Apresentamos então uma descrição do código numérico desenvolvido usando linguagem Fortran, abordando um plasma com elétrons descritos por uma função distribuição Maxwelliana com um feixe tênue, e íons descritos por uma função distribuição Maxwelliana. Finalmente, comparamos nossos resultados com outros obtidos por outros autores em trabalhos anteriores, desenvolvidos usando coordenadas cartesianas, avaliando nosso trabalho. Por último, discutimos algumas perspectivas para o desenvolvimento futuro do trabalho. / In the present work we study the beam-plasma interaction using the weak turbulence theory. We start with an introduction to the kinetic theory of plasmas, studying fundamental features, like the statistical approach to the plasma equations, the procedures to be employed to deal with the correlation functions appearing in the statistical approach, and the system of Vlasov-Maxwell equations. In the sequence we discuss the Vlasov approximation and the solution of the Vlasov- Maxwell system in the linear approximation, emphasizing the description of waves in the plasma and the Landau damping. Along the development we present a relatively detailed description of the procedures introduced by Landau to treat the Vlasov-Maxwell system as an initial value problem, including a discussion about the resolution of integrals in the complex plane, with poles in the denominator, and the solution of the dispersion relation to find the normal mode of oscillations in the plasma. We also present a short review about the quasilinear approximation of the Vlasov-Maxwell system, in which we discuss the derivation of the quasilinear diffusion equation in the space of particle velocities, and study its properties of conservation. In the context of the quasilinear theory, we present a short review about the time evolution of the Landau damping in the case of Maxwellian velocity distribution, and about the instability which can occur in the case of a distribution function with a beam of particles. After the introduction to the linear approximation and to the quasilinear theory, we present the equations of weak turbulence theory for a unmagnetized plasma, including non-linear interactions between waves and particles. The formalism which is presented includes effects and spontaneous and induced emission, decay and scattering of waves. The system of coupled equations of the weak turbulence theory is then reduced to application to a bi-dimensional case, and then re-written in terms of polar coordinates and adapted to numerical solution. We then present a description of the numerical code developed using Fortran language, suitable to describe a plasma with electrons described by a Maxwellian distribution function with a tenuous beam, and ions described by a Maxwellian distribution. Finally, we compare our results with results obtained by other authors in previous works, developed using cartesian coordinates, as a validation of our work. Lastly we discuss some perspectives for future developments.

A situação em que ocorre a interação entre um feixe de elétrons rápidos e um plasma é considerada um modelo fundamental para o estudo de instabilidades cinéticas em plasmas. A dinâmica não linear desta interação vem sendo amplamente estudada, com a intenção de entender melhor os processos por trás da geração de turbulência de Langmuir, observada durante os chamados type II radio bursts e type III radio bursts. Nesses estudos, é comum que seja suposto um feixe de elétrons de baixa densidade, o que leva a um longo intervalo de evolução temporal, muito maior do que o período das ondas que estão sendo consideradas. Para longos períodos de evolução temporal, uma investigação a respeito da atuação dos efeitos colisionais se faz necessária. No presente trabalho, incluímos os efeitos das colisões binárias ao conjunto completo de processos não lineares da instabilidade bump-on-tail, na presença de oscilações eletrostáticas, e obtemos a evolução temporal bidimensional desse sistema. Para isso começamos com uma revisão teórica que inclui teoria cinética de plasmas, teoria de turbulência fraca e um estudo sobre colisões em plasmas. A seguir, adaptamos o operador colisional de Fokker-Planck, na aproximação de Landau, de maneira que pudéssemos agrupá-lo à equação cinética bidimensional para a evolução das partículas. Feito isso, passamos essas equações para a forma de diferenças finitas e as integramos em conjunto com o restante das equações da teoria de turbulência fraca. Uma comparação dos resultados obtidos com a evolução do mesmo sistema na ausência de colisões, nos indica que os primeiros sinais da atuação colisional na função de distribuição, embora muito sutis, surgem logo no início da evolução temporal. No entanto, a evolução das ondas só será afetada pelas colisões bem mais tarde, quando, aparentemente, os processos colisionais passam a dominar a dinâmica do sistema, superando os processos coletivos não lineares. Essa predominância da atuação colisional é percebida como uma lenta, mas progressiva simetrização das distribuições no espaço de velocidades e no espaço de vetor de onda. Os resultados estão de acordo com o esperado e as perspectivas futuras são promissoras. / The situation in which occurs the interaction between a beam of fast electrons and a plasma is considered a fundamental model for the study of kinetic instabilities in plasmas. The nonlinear dynamics of this interaction has been extensively studied with the intention to better understand the processes behind the Langmuir turbulence generation, observed during the so-called type II and type III radio bursts. In these studies, is usual to assume a low density electron beam, which leads to a long time development interval; much longer than the period of the considered waves. For long periods of time evolution, an investigation concerning the action of collisional effects is needed. In this study, we include the effects of binary collisions to the complete set of non-linear processes of the bump-on-tail instability in the presence of electrostatic oscillations, and obtain the solution for the time evolution of the system in two dimensions. For that, we start with a literature review which includes kinetic theory of plasmas, weak turbulence theory and a study of collisions in plasmas. Then we adapt the Fokker-Planck collisional operator, in Landau approach, so we could group it to the two-dimensional kinetic equation for the evolution of the particles. Then, we write these equations in the form of finite difference equations, and integrate them together with the rest of the weak turbulence equations. A comparison of the results obtained with the evolution of the same system in the absence of collisions, indicates that the first signs of collisional activity in the distribution function, although very subtle, appear at the beginning of time evolution. However, the evolution of the waves will only be affected by collisions much later, when the collisional processes apparently dominate the dynamics of the system, overcoming the nonlinear collective processes. This predominance of collisional acting is perceived as a slow but progressive symmetrization of the distributions in the velocity space and the wave vector space. The results are consistent with expectations, and future perspectives are promising.

Neste trabalho fizemos uma revisão da teoria cinética de plasmas e obtivemos as equações cinéticas para um plasma empoeirado magnetizado com carga variável das partículas de poeira. Além das correntes de carregamento por absorção de partículas do plasma, incluímos a corrente de emissão de fotoelétrons a corrente total, e consideramos que o plasma é povoado por poeira com partículas com diferentes tamanhos. / In this work we have made a review of the plasmas kinetic theory and obtained the kinetic equations for a magnetized dusty plasma with variable charges of the dust particles. In addition to currents due to absorption of plasma particles, we have included in the total current the photoelectrons emission current and considered that the dusty plasma is populated by dust with particles with di erent sizes.

Neste trabalho, uma versão analítica do método de ordenadas discretas (ADO) é utilizada no desenvolvimento de soluções para uma ampla classe de problemas de gases rarefeitos em semi-espa»co e canal plano. A modelagem dos problemas baseada em modelos cinéticos derivados da equação linearizada de Boltzmann, tais como os modelos BGK, o S, o Gross-Jackson e o MRS. Em particular, resultados para o modelo MRS são originais. A solução ADO se mostrou eficiente e precisa e uma série de resultados são apresentados no sentido de estabelecer uma comparação entre os modelos propriamente ditos e resultados obtidos a partir da ELB. Além disso, uma análise da influência de alguns parâmetros é apresentada, para todos os problemas. / In this work, an analytical version of the discrete-ordinates method (the ADO method) is used to develop solutions for a wide class of rare¯ed gas problems in half- space and plane channel. The modelling of the problems is based on kinetic models of the linearized Boltzmann equation, such that the BGK, the S, the Gross-Jackson and the MRS models. In particular, results for the MRS model are originals. The ADO solution was founded to be e±cient and accurate and a series of results are presented in order to establish a comparison between the kinetic models themselves and results provided by the ELB. In addition, an analysis of the in°uence of some parameters was presented, for all problems.

A situação em que ocorre a interação entre um feixe de elétrons rápidos e um plasma é considerada um modelo fundamental para o estudo de instabilidades cinéticas em plasmas. A dinâmica não linear desta interação vem sendo amplamente estudada, com a intenção de entender melhor os processos por trás da geração de turbulência de Langmuir, observada durante os chamados type II radio bursts e type III radio bursts. Nesses estudos, é comum que seja suposto um feixe de elétrons de baixa densidade, o que leva a um longo intervalo de evolução temporal, muito maior do que o período das ondas que estão sendo consideradas. Para longos períodos de evolução temporal, uma investigação a respeito da atuação dos efeitos colisionais se faz necessária. No presente trabalho, incluímos os efeitos das colisões binárias ao conjunto completo de processos não lineares da instabilidade bump-on-tail, na presença de oscilações eletrostáticas, e obtemos a evolução temporal bidimensional desse sistema. Para isso começamos com uma revisão teórica que inclui teoria cinética de plasmas, teoria de turbulência fraca e um estudo sobre colisões em plasmas. A seguir, adaptamos o operador colisional de Fokker-Planck, na aproximação de Landau, de maneira que pudéssemos agrupá-lo à equação cinética bidimensional para a evolução das partículas. Feito isso, passamos essas equações para a forma de diferenças finitas e as integramos em conjunto com o restante das equações da teoria de turbulência fraca. Uma comparação dos resultados obtidos com a evolução do mesmo sistema na ausência de colisões, nos indica que os primeiros sinais da atuação colisional na função de distribuição, embora muito sutis, surgem logo no início da evolução temporal. No entanto, a evolução das ondas só será afetada pelas colisões bem mais tarde, quando, aparentemente, os processos colisionais passam a dominar a dinâmica do sistema, superando os processos coletivos não lineares. Essa predominância da atuação colisional é percebida como uma lenta, mas progressiva simetrização das distribuições no espaço de velocidades e no espaço de vetor de onda. Os resultados estão de acordo com o esperado e as perspectivas futuras são promissoras. / The situation in which occurs the interaction between a beam of fast electrons and a plasma is considered a fundamental model for the study of kinetic instabilities in plasmas. The nonlinear dynamics of this interaction has been extensively studied with the intention to better understand the processes behind the Langmuir turbulence generation, observed during the so-called type II and type III radio bursts. In these studies, is usual to assume a low density electron beam, which leads to a long time development interval; much longer than the period of the considered waves. For long periods of time evolution, an investigation concerning the action of collisional effects is needed. In this study, we include the effects of binary collisions to the complete set of non-linear processes of the bump-on-tail instability in the presence of electrostatic oscillations, and obtain the solution for the time evolution of the system in two dimensions. For that, we start with a literature review which includes kinetic theory of plasmas, weak turbulence theory and a study of collisions in plasmas. Then we adapt the Fokker-Planck collisional operator, in Landau approach, so we could group it to the two-dimensional kinetic equation for the evolution of the particles. Then, we write these equations in the form of finite difference equations, and integrate them together with the rest of the weak turbulence equations. A comparison of the results obtained with the evolution of the same system in the absence of collisions, indicates that the first signs of collisional activity in the distribution function, although very subtle, appear at the beginning of time evolution. However, the evolution of the waves will only be affected by collisions much later, when the collisional processes apparently dominate the dynamics of the system, overcoming the nonlinear collective processes. This predominance of collisional acting is perceived as a slow but progressive symmetrization of the distributions in the velocity space and the wave vector space. The results are consistent with expectations, and future perspectives are promising.

Neste trabalho fizemos uma revisão da teoria cinética de plasmas e obtivemos as equações cinéticas para um plasma empoeirado magnetizado com carga variável das partículas de poeira. Além das correntes de carregamento por absorção de partículas do plasma, incluímos a corrente de emissão de fotoelétrons a corrente total, e consideramos que o plasma é povoado por poeira com partículas com diferentes tamanhos. / In this work we have made a review of the plasmas kinetic theory and obtained the kinetic equations for a magnetized dusty plasma with variable charges of the dust particles. In addition to currents due to absorption of plasma particles, we have included in the total current the photoelectrons emission current and considered that the dusty plasma is populated by dust with particles with di erent sizes.

Neste trabalho um modelo é apresentado que permite-nos prever quantitativamente o estado estacionário atingido por um plasma não-neutro durante um processo de relaxação nãocolisional. Como uma aplicação específica, a teoria é usada para estudar a relaxação de feixes de partículas carregadas na qual, mostra-se que, um feixe inicialmente casado relaxa para a distribuição de Lynden-Bell. No entanto, quando existe um descasamento inicial, o mesmo oscila e ressonâncias paramétricas conduzem-no a uma separação de fases: caroço e halo. A abordagem desenvolvida vale tanto para a densidade como para a distribuição de velocidade no estado final estacionário. / In this work a theoretical framework is presented which allows us to quantitatively predict the final stationary state achieved by a non-neutral plasma during a process of collisionless relaxation. As a specific application, the theory is used to study relaxation of charged-particle beams. It is shown that a fully matched beam relaxes to the Lynden-Bell distribution. However, when a mismatch is present and the beam oscillates, parametric resonances lead to a corehalo phase separation. The approach developed accounts for both the density and the velocity distributions in the final stationary state.

Neste trabalho estudamos a interação feixe-plasma aplicando a teoria de turbulência fraca. Primeiramente fazemos uma introdução `a teoria cinética de plasmas, estudando aspectos fundamentais como: a abordagem estatística das equações do plasma, o tratamento a ser dado `as funções de correlação que aparecem nessa abordagem, e o sistema de equações de Vlasov-Maxwell. Em seguida estudamos a aproximação de Vlasov e a solução do sistema de Vlasov- Maxwell na aproximação linear, enfocando a descrição de ondas no plasma e o amortecimento de Landau. Ao longo desse desenvolvimento fazemos uma apresentação relativamente detalhada dos procedimentos introduzidos por Landau ao tratar o sistema Vlasov-Maxwell como um problema de valor inicial, incluindo uma discussão sobre a resolução de integrais no plano complexo, com polos no denominador, e a solução da relação de dispersão para encontrar os modos normais de oscilação no plasma. Apresentamos também uma breve revisão a respeito da aproximação quase-linear do sistema de Vlasov- Maxwell, na qual abordamos a obtenção da equação quase-linear de difusão no espaço das velocidades para partículas e estudamos suas propriedades de conservação. Nesse contexto da teoria quase-linear, apresentamos uma revisão a respeito da evolução temporal do amortecimento de Landau no caso de distribuição de velocidades Maxwelliana e da instabilidade que pode ocorrer no caso de uma função distribuição com um feixe de partículas. Depois dessa introdução `a aproximação linear e `a teoria quase-linear, abordamos a teoria de turbulência fraca num plasma não magnetizado, incluindo interações não lineares entre ondas e partículas. O formalismo apresentado inclui efeitos como emissão espontânea e induzida, decaimento e espalhamento de ondas. O sistema de equações acopladas da teoria de turbulência fraca é então reduzido a uma aplicação a um sistema considerando duas dimensões, depois reescrito em termos de coordenadas polares e adaptado para solução numérica. Apresentamos então uma descrição do código numérico desenvolvido usando linguagem Fortran, abordando um plasma com elétrons descritos por uma função distribuição Maxwelliana com um feixe tênue, e íons descritos por uma função distribuição Maxwelliana. Finalmente, comparamos nossos resultados com outros obtidos por outros autores em trabalhos anteriores, desenvolvidos usando coordenadas cartesianas, avaliando nosso trabalho. Por último, discutimos algumas perspectivas para o desenvolvimento futuro do trabalho. / In the present work we study the beam-plasma interaction using the weak turbulence theory. We start with an introduction to the kinetic theory of plasmas, studying fundamental features, like the statistical approach to the plasma equations, the procedures to be employed to deal with the correlation functions appearing in the statistical approach, and the system of Vlasov-Maxwell equations. In the sequence we discuss the Vlasov approximation and the solution of the Vlasov- Maxwell system in the linear approximation, emphasizing the description of waves in the plasma and the Landau damping. Along the development we present a relatively detailed description of the procedures introduced by Landau to treat the Vlasov-Maxwell system as an initial value problem, including a discussion about the resolution of integrals in the complex plane, with poles in the denominator, and the solution of the dispersion relation to find the normal mode of oscillations in the plasma. We also present a short review about the quasilinear approximation of the Vlasov-Maxwell system, in which we discuss the derivation of the quasilinear diffusion equation in the space of particle velocities, and study its properties of conservation. In the context of the quasilinear theory, we present a short review about the time evolution of the Landau damping in the case of Maxwellian velocity distribution, and about the instability which can occur in the case of a distribution function with a beam of particles. After the introduction to the linear approximation and to the quasilinear theory, we present the equations of weak turbulence theory for a unmagnetized plasma, including non-linear interactions between waves and particles. The formalism which is presented includes effects and spontaneous and induced emission, decay and scattering of waves. The system of coupled equations of the weak turbulence theory is then reduced to application to a bi-dimensional case, and then re-written in terms of polar coordinates and adapted to numerical solution. We then present a description of the numerical code developed using Fortran language, suitable to describe a plasma with electrons described by a Maxwellian distribution function with a tenuous beam, and ions described by a Maxwellian distribution. Finally, we compare our results with results obtained by other authors in previous works, developed using cartesian coordinates, as a validation of our work. Lastly we discuss some perspectives for future developments.