Hydrodynamics of Plasma Interaction with an Intense Laser Field

A set of moment relations, which can describe the charged fluids response to an intense pump laser, and a linearization substitution relation, which is more appropriate as compared with the past treatment, are given by theoretical analyses. The relevant equations of state (adiabatic and isothermal), momentum and energy equations are derived self-consistently.The dispersion relations of the electron plasma wave and the ion acoustic wave driven by an intense pump laser field are-obtained. The results show that the frequencies of both the excited electron plasma wave and the excited ion acoustic wave have a great modification in the case of strong pump. The former bears out the theoretical result obtained from Vlasov equation and the later is consistent with experimental observations. It is proved that the zero-frequency component of the laser light wave contribution to the plasma pressure tensor is un-neglected,which implies a greatly change to the wave excitation properties, particularly in the direction of parallel or approximately parallel to the laser field vector.     

Kinetic boundary layers in gas mixtures: Systems described by nonlinearly coupled kinetic and hydrodynamic equations and applications to droplet condensation and evaporation

We consider a mixture of heavy vapor molecules and a light carrier gas surrounding a liquid droplet. The vapor is described by a variant of the Klein-Kramers equation, a kinetic equation for Brownian particles moving in a spatially inhomogeneous background; the gas is described by the Navier-Stokes equations; the droplet acts as a heat source due to the released heat of condensation. The exchange of momentum and energy between the constituents of the mixture is taken into account by force terms in the kinetic equation and source terms in the Navier-Stokes equations. These are chosen to obtain maximal agreement with the irreversible thermodynamics of a gas mixture. The structure of the kinetic boundary layer around the sphere is then determined from the self-consistent solution of this set of coupled equations with appropriate boundary conditions at the surface of the sphere. For this purpose the kinetic equation is rewritten as a set of coupled moment equations. A complete set of solutions of these moment equations is constructed by numerical integration inward from the region far away from the droplet, where the background inhomogeneities are small. A technique developed in an earlier paper is used to deal with the severe numerical instability of the moment equations. The solutions so obtained for given temperature and pressure profiles in the gas are then combined linearly in such a way that they obey the boundary conditions at the droplet surface; from this solution source terms for the Navier-Stokes equation of the gas are constructed and used to determine improved temperature and pressure profiles for the background gas. For not too large temperature differences between the droplet and the gas at infinity, self-consistency is reached after a few iterations. The method is applied to the condensation of droplets from a supersaturated vapor, where small but significant corrections to an earlier, not fully consistent version of the theory are found, as well as to strong evaporation of droplets under the influence of an external heat source, where corrections of up to 40 % are obtained.     

The interaction of transverse domain walls

The interaction between transverse domain walls is calculated analytically using a multipole expansion up to third order. Starting from an analytical expression for the magnetization in the wall, the monopole, dipole, and quadrupole moments are derived and their impact on the interaction is investigated using the surface and volume charges. The surface charges are important for the dipole moment while the volume charges constitute the monopole and quadrupole moments. For domain walls that are situated in different wires it is found that there is a strong deviation from the interaction of two monopoles. This deviation is caused by the interaction of the monopole of the wall in the first wire with the dipole of the wall in the second wire and vice versa. The dipole-dipole and the quadrupole-monopole interactions are found to be also of considerable size and non-negligible. A comparison with micromagnetic simulations shows a good agreement.     

等离子体中非线性朗谬尔波的哈密顿描述

研究了在双离子（H^ ,O^ )成份等离子体中的非线性朗谬尔波的特性,从流体方程出发,考虑低频离子运动的完全非线性和双极势的色散,得到了描述高频电场缓变振幅与低频势扰动的耦合方程组。利用哈密顿方法,在小振幅情况下,对方程组解耦合,利用Sagdeev势方法,对孤立波的性状进行了讨论,结果表明,双离子成份等离子体中双极势的孤立子的幅度相对电子,离子等离子体的双极势孤立子的幅度要大,而相应的高频电场强度的幅度相对要小,相应的孤立子的速度传播区也较小。     

Numerical Simulation of MHD Peristaltic Flow with Variable Electrical Conductivity and Joule Dissipation Using Generalized Differential Quadrature Method

In this paper, the MHD peristaltic flow inside wavy walls of an asymmetric channel is investigated, where the walls of the channel are moving with peristaltic wave velocity along the channel length. During this investigation,the electrical conductivity both in Lorentz force and Joule heating is taken to be temperature dependent. Also, the long wavelength and low Reynolds number assumptions are utilized to reduce the governing partial differential equations into a set of coupled nonlinear ordinary differential equations. The new set of obtained equations is then numerically solved using the generalized differential quadrature method(GDQM). This is the first attempt to solve the nonlinear equations arising in the peristaltic flows using this method in combination with the Newton-Raphson technique. Moreover, in order to check the accuracy of the proposed numerical method, our results are compared with the results of built-in Mathematica command NDSolve. Taking Joule heating and viscous dissipation into account, the effects of various parameters appearing in the problem are used to discuss the fluid flow characteristics and heat transfer in the electrically conducting fluids graphically. In presence of variable electrical conductivity, velocity and temperature profiles are highly decreasing in nature when the intensity of the electrical conductivity parameter is strengthened.     

Classical and semirelativistic magnetohydrodynamics with anisotropic ion pressure

We study the magnetohydrodynamics (MHD) equations with anisotropic ion pressure and isotropic electron pressure under both the classical and semirelativistic approximations in order to develop a numerical model. The dispersion relation as well as the characteristic wave speeds are derived. In addition to the exact wave speed solutions, we also provide efficient approximate formulas for the semirelativistic magnetosonic speeds. The equations are discretized with the Rusanov and Harten-Lax-van Leer numerical schemes and implemented into the BATS-R-US MHD code. We perform a set of verification tests.     

Flow oscillations in radial expansion of an inhomogeneous plasma layer

The cylindrically symmetric radial evolution of an inhomogeneous plasma layer expanding into vacuum is investigated nonperturbatively by first determining the spatial structure of the plasma flow structure. The evolution is then governed by a set of ordinary differential equations. The effect of the plasma inhomogeneity on the nonlinear coupling among the electron and ion flow components and oscillations is investigated.     

Moment Equations for a Granular Material in a Thermal Bath

We compute the moment equations for a granular material under the simplifying assumption of pseudo-Maxwellian particles approximating dissipative hard spheres. We obtain the general moment equations of second and third order and the isotropic moment equations of any order. Our equations describe, in the space homogeneous case, the granular system described by a Boltzmann-like collision term and subject to a Brownian motion due to the interaction with a bath, described by a Fokker–Planck term. The trend to equilibrium is studied in detail.     

The use of polarized coordinates for solving the macroscopic equations of collisionless plasma

A method is developed for the iterative solution of an infinite system of moment equations derived from the collisionless Boltzmann (Vlasov) equation. Since in a strong magnetic field the dominant terms of these equations have a complicated form, the iterative solution after a few steps comes up against difficulties connected with the inversion of matrices of high order. It is found that the transformation of moment equations into the natural polarized coordinate system (Buneman O.: Phys. Fluids4 (1961), 669) diagonalizes the dominant matrices in the magnetic fields with straight lines of force exactly and in adiabatic magnetic fields approximately. In addition, this transformation reveals the symmetries of the moment equations which in some cases permit an exact solution to be found. Apart from the first three moment equations (the equation of continuity, the dynamic (Euler) equation, the equation for pressure) the general moment equation of then-th order is also derived in the Cartesian and general curvilinear coordinate system and in the corresponding polarized natural coordinate systems. The paper deals with some special cases of natural polarized coordinate systems belonging to a magnetic field with straight lines of force (plane, cylindrical and spherical geometry) and with curved lines of force (magnetic trap with mirrors).     

Conditional quadrature method of moments for kinetic equations

Kinetic equations arise in a wide variety of physical systems and efficient numerical methods are needed for their solution. Moment methods are an important class of approximate models derived from kinetic equations, but require closure to truncate the moment set. In quadrature-based moment methods (QBMM), closure is achieved by inverting a finite set of moments to reconstruct a point distribution from which all unclosed moments (e.g. spatial fluxes) can be related to the finite moment set. In this work, a novel moment-inversion algorithm, based on 1-D adaptive quadrature of conditional velocity moments, is introduced and shown to always yield realizable distribution functions (i.e. non-negative quadrature weights). This conditional quadrature method of moments (CQMOM) can be used to compute exact N-point quadratures for multi-valued solutions (also known as the multi-variate truncated moment problem), and provides optimal approximations of continuous distributions. In order to control numerical errors arising in volume averaging and spatial transport, an adaptive 1-D quadrature algorithm is formulated for use with CQMOM. The use of adaptive CQMOM in the context of QBMM for the solution of kinetic equations is illustrated by applying it to problems involving particle trajectory crossing (i.e. collision-less systems), elastic and inelastic particle–particle collisions, and external forces (i.e. fluid drag).     

Equations for finite-difference, time-domain simulation of sound propagation in moving inhomogeneous media and numerical implementation

Finite-difference, time-domain (FDTD) calculations are typically performed with partial differential equations that are first order in time. Equation sets appropriate for FDTD calculations in a moving inhomogeneous medium (with an emphasis on the atmosphere) are derived and discussed in this paper. Two candidate equation sets, both derived from linearized equations of fluid dynamics, are proposed. The first, which contains three coupled equations for the sound pressure, vector acoustic velocity, and acoustic density, is obtained without any approximations. The second, which contains two coupled equations for the sound pressure and vector acoustic velocity, is derived by ignoring terms proportional to the divergence of the medium velocity and the gradient of the ambient pressure. It is shown that the second set has the same or a wider range of applicability than equations for the sound pressure that have been previously used for analytical and numerical studies of sound propagation in a moving atmosphere. Practical FDTD implementation of the second set of equations is discussed. Results show good agreement with theoretical predictions of the sound pressure due to a point monochromatic source in a uniform, high Mach number flow and with Fast Field Program calculations of sound propagation in a stratified moving atmosphere.     

Dust Alfven ordinary and cusp solitons and modulational instability in a self-gravitating magneto-radiative plasma

The effective one-fluid adiabatic magnetohydrodynamic (MHD) equations for a multicomponent plasma comprising of electrons, ions, and dust are used to investigate the nonlinear coupling of dust Alfven and dust acoustic waves in the presence of radiation pressure as well as the Jean’s term that arises in a self-gravitating plasma. In this context, the set of Zakharov equations are derived. The soliton solutions in the presence of radiation pressure and Jeans term are separately discussed. It is found that ordinary solitons are obtained in the absence of Jeans term whereas cusp solitons are found in the absence of radiation pressure. To the best of our knowledge, cusp solitons are obtained for the first time for a self-gravitating plasma with Jeans term for an electromagnetic wave in a dusty plasma. The modulational instability is also investigated in the presence of radiation pressure and Jeans term. It is found that the Jeans term drives the system modulationally unstable provided it dominates the dust acoustic and radiation pressure terms whereas the radiation pressure enhances the stability of the system. The relevance of the present investigation in the photodissociation region that separates the HII region from the dense molecular clouds is also pointed out.     

Acoustic symmetry and the third order elastic constants of tetragonal TII (4/m) crystals

The mode equations which relate the measured hydrostatic and uniaxial pressure derivatives of the “natural” ultrasonic velocity to the second and third order elastic constants for the tetragonal TII (4/m) Laue group have been transformed from the crystallographic (XYZ) axial frame to a reference frame (k, k + π/2, Z) comprised of two acoustic symmetry axes and the fourfold axis. It is found that, if the acoustic symmetry axes do not shift appreciably with pressure, the mode equations for TII crystals in this reference frame are much simpler than those in the crystallographie axial frame (XYZ), having the same form as those for the higher symmetry TI Laue group. Hence 12 transformed third order elastic constants can be obtained experimentally for TII crystals by making measurements of the hydrostatic and uniaxial pressure derivatives of ultrasonic wave velocity of modes propagated in the (k, k + π/2, Z) axial set with the uniaxial stresses also applied in this reference frame. The approach has been used to determine the 12 transformed third order elastic constants of scheelite (CaWO₄).     

Positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations with source terms

In [16], [17], we constructed uniformly high order accurate discontinuous Galerkin (DG) schemes which preserve positivity of density and pressure for the Euler equations of compressible gas dynamics with the ideal gas equation of state. The technique also applies to high order accurate finite volume schemes. For the Euler equations with various source terms (e.g., gravity and chemical reactions), it is more difficult to design high order schemes which do not produce negative density or pressure. In this paper, we first show that our framework to construct positivity-preserving high order schemes in [16], [17] can also be applied to Euler equations with a general equation of state. Then we discuss an extension to Euler equations with source terms. Numerical tests of the third order Runge–Kutta DG (RKDG) method for Euler equations with different types of source terms are reported.     

The steady state out-of-plane response of a Timoshenko curved beam with internal damping

The steady state out-of-plane response of a Timoshenko curved beam with internal damping to a sinusoidally varying point force or moment is determined by use of the transfer matrix approach. For this purpose, the equations of out-of-plane vibration of a curved beam are written as a coupled set of the first order differential equations by using the transfer matrix of the beam. Once the matrix has been determined by numerical integration of the equations, the steady state response of the beam is obtained. The method is applied to free-clamped non-uniform beams with circular, elliptical, catenary and parabolical neutral axes driven at the free end; the driving point impedance and force or moment transmissibility are calculated numerically and the effects of the slenderness ratio, varying cross-section and the function expressing the neutral axis on them are studied.     

New Applications of Langmuir Probes

In this work, two new possibilities for standard probe diagnostics are described. The first one can be used to study isotropic, collisionless low-pressure plasma in which the electron energy distribution function is close to a Maxwellian one. In such plasmas, the Boltzmann law, Bohm effect, and 3/2 power law are valid. Use of corresponding system of equations for cylindrical Langmuir probes allowed for measurements of probe sheath thicknesses and the mean ion mass. The solution of this task was provided by accurate probe diagnostics of inductive xenon plasma at pressure p = 2 mTorr that resulted in the determination of the Bohm coefficient C_BCyl = 1.22. The second possibility of probe diagnostics includes a method and device for evaluation of ion current density to a wall under a floating potential using a radially movable plane wall Langmuir probe simulator. This measurement in the same xenon plasma served as the basis for development of an ion source in which the given wall was represented by an ion extracting electrode of the ion extraction grid system.