Kinetic Theory of Jet Dynamics in the Stochastic Barotropic and 2D Navier-Stokes Equations

We discuss the dynamics of zonal (or unidirectional) jets for barotropic flows forced by Gaussian stochastic fields with white in time correlation functions. This problem contains the stochastic dynamics of 2D Navier-Stokes equation as a special case. We consider the limit of weak forces and dissipation, when there is a time scale separation between the inertial time scale (fast) and the spin-up or spin-down time (large) needed to reach an average energy balance. In this limit, we show that an adiabatic reduction (or stochastic averaging) of the dynamics can be performed. We then obtain a kinetic equation that describes the slow evolution of zonal jets over a very long time scale, where the effect of non-zonal turbulence has been integrated out. The main theoretical difficulty, achieved in this work, is to analyze the stationary distribution of a Lyapunov equation that describes quasi-Gaussian fluctuations around each zonal jet, in the inertial limit. This is necessary to prove that there is no ultraviolet divergence at leading order, in such a way that the asymptotic expansion is self-consistent. We obtain at leading order a Fokker–Planck equation, associated to a stochastic kinetic equation, that describes the slow jet dynamics. Its deterministic part is related to well known phenomenological theories (for instance Stochastic Structural Stability Theory) and to quasi-linear approximations, whereas the stochastic part allows to go beyond the computation of the most probable zonal jet. We argue that the effect of the stochastic part may be of huge importance when, as for instance in the proximity of phase transitions, more than one attractor of the dynamics is present.     

Particle Creation in Anisotropically Expanding Universe

Using squeezed vacuum states formalism of quantum optics, a homogeneous and massive scalar field minimally coupled to gravity in Bianchi type-I model of the universe is examined in the frame work of semiclassical theory of gravity. Hence an approximate leading solution to the semiclassical Einstein equation is found. The next order solution for each scale factor in their respective direction show power law of expansion. It is further noted that evolution of scale factors are mutually correlated. The phenomena of nonclassical particle creation is also examined in the anisotropic background cosmology.     

Unstrained and strained flamelets for LES of premixed combustion

The unstrained and strained flamelet closures for filtered reaction rate in large eddy simulation (LES) of premixed flames are studied. The required sub-grid scale (SGS) PDF in these closures is presumed using the Beta function. The relative performances of these closures are assessed by comparing numerical results from large eddy simulations of piloted Bunsen flames of stoichiometric methane–air mixture with experimental measurements. The strained flamelets closure is observed to underestimate the burn rate and thus the reactive scalars mass fractions are under-predicted with an over-prediction of fuel mass fraction compared with the unstrained flamelet closure. The physical reasons for this relative behaviour are discussed. The results of unstrained flamelet closure compare well with experimental data. The SGS variance of the progress variable required for the presumed PDF is obtained by solving its transport equation. An order of magnitude analysis of this equation suggests that the commonly used algebraic model obtained by balancing source and sink in this transport equation does not hold. This algebraic model is shown to underestimate the SGS variance substantially and the implications of this variance model for the filtered reaction rate closures are highlighted.     

A priori analysis of sub-grid variance of a reactive scalar using DNS data of high Ka flames

Direct numerical simulations (DNS) of low and high Karlovitz number (Ka) flames are analysed to investigate the behaviour of the reactive scalar sub-grid scale (SGS) variance in premixed combustion under a wide range of combustion conditions (regimes). An order of magnitude analysis is performed to assess the importance of various terms in the variance evolution equation and the analysis is validated using the DNS results. This analysis sheds light on the relative behaviour among turbulent transport and production, scalar dissipation and chemical processes involved in the evolution of the SGS variance at different Ka. The common expectation is that the variance equation shifts from a reaction-dissipation balance at low Ka to a production–dissipation balance at high Ka with diminishing reaction contribution. However, in large eddy simulation (LES), a high Ka alone does not make the reaction term negligible, as the relative importance of the reaction term has a concurrent increase with filter size. The filter size can be relatively large compared with the Kolmogorov length scale in practical LES of high Ka flames, and as a consequence a reaction–production–dissipation balance may prevail in the variance equation even in a high Ka configuration, and this possibility is quantified using the DNS analysis in this work. This has implications from modelling perspectives, and therefore two commonly used closures in LES for the SGS scalar dissipation rate are investigated a priori to estimate the importance of the above balance in LES modelling. The results are explained to highlight the interplay among turbulence, chemistry and dissipation processes as a function of Ka.     

Higher-order quadrature-based moment methods for kinetic equations

Kinetic equations containing terms for spatial transport, body forces, and particle–particle collisions occur in many applications (e.g., rarefied gases, dilute granular gases, fluid-particle flows). The direct numerical solution of the kinetic equation is usually intractable due to the large number of independent variables. A useful alternative is to reformulate the problem in terms of the moments of the velocity distribution function. Closure of the moment equations is challenging for flows sufficiently far away from the Maxwellian limit. In previous work, a quadrature-based third-order moment closure was derived for approximating solutions to the kinetic equation for arbitrary Knudsen number. A key component of quadrature-based closures is the moment-inversion algorithm used to find the non-negative weights and velocity abscissas. Here, a robust inversion procedure is proposed for three-component velocity moments up to ninth order. By reconstructing the velocity distribution function, the spatial fluxes in the moment equations are treated using a kinetic-based finite-volume solver. Because the quadrature-based moment method employs the moment transport equations directly instead of a discretized form of the kinetic equation, the mass, momentum and energy are conserved for arbitrary Knudsen and Mach numbers. The computational algorithm is tested for the Riemann shock problem and, for increasing Knudsen numbers (i.e. larger deviations from the Maxwellian limit), the accuracy of the moment closure is shown to be determined by the discrete representation of the spatial fluxes.     

$\mathcal{PT}$ Invariant Complex E8 Root Spaces

We provide a construction procedure for complex root spaces invariant under antilinear transformations, which may be applied to any Coxeter group. The procedure is based on the factorisation of a chosen element of the Coxeter group into two factors. Each of the factors constitutes an involution and may therefore be deformed in an antilinear fashion. Having the importance of the E ₈-Coxeter group in mind, such as underlying a particular perturbation of the Ising model and the fact that for it no solution could be found previously, we exemplify the procedure for this particular case. As a concrete application of this construction we propose new generalisations of Calogero-Moser-Sutherland models and affine Toda field theories based on the invariant complex root spaces and deformed complex simple roots, respectively.     

Renormalization Group Invariance in Pionless Effective Field Theory for the NN System

We consider the NN interaction in pionless effective field theory (EFT) up to next-to-next-to-leading order (NNLO) and use a recursive subtractive renormalization scheme to describe NN scattering in the ¹ S ₀ channel. We fix the strengths of the contact interactions at a reference scale, chosen to be the one that provides the best fit for the phase-shifts, and then slide the renormalization scale by evolving the driving terms of the subtracted Lippmann?CSchwinger equation through a non-relativistic Callan?CSymanzik equation. The results show that such a systematic renormalization scheme with multiple subtractions is fully renormalization group invariant.     

Convergence to SPDEs in Stratonovich Form

We consider the perturbation of parabolic operators of the form ∂ _t  + P(x, D) by large-amplitude highly oscillatory spatially dependent potentials modeled as Gaussian random fields. The amplitude of the potential is chosen so that the solution to the random equation is affected by the randomness at the leading order. We show that, when the dimension is smaller than the order of the elliptic pseudo-differential operator P(x, D), the perturbed parabolic equation admits a solution given by a Duhamel expansion. Moreover, as the correlation length of the potential vanishes, we show that the latter solution converges in distribution to the solution of a stochastic parabolic equation with multiplicative noise that should be interpreted in the Stratonovich sense. The theory of mild solutions for such stochastic partial differential equations is developed. The behavior described above should be contrasted to the case of dimensions larger than or equal to the order of the elliptic pseudo-differential operator P(x, D). In the latter case, the solution to the random equation converges strongly to the solution of a homogenized (deterministic) parabolic equation as is shown in [2]. A stochastic limit is obtained only for sufficiently small space dimensions in this class of parabolic problems.     

On primitive variable behaviour near shocks in ensemble-averaged methods

Models for averaged shock corrugation effects and the impact of turbulent entropy or acoustic modes on the energy equation are presented, for application in Reynolds-Averaged Navier Stokes(RANS) simulations of shock-turbulence interactions. Unlike previous work that has focused the modification of turbulent statistics by the shock, the proposed models are introduced to capture the effects of the turbulence on the profiles of primitive variables - mean density, velocity, and pressure. By producing accurate profiles for the primitive variables, it is shown that the proposed models improve numerical convergence behaviour with mesh refinement about a shock, and introduce the physical effects of shock asphericity in a converging shock geometry. These effects are achieved by local closures to turbulent statistics in the averaged Navier-Stokes equations, and can be applied in conjunction with existing Reynolds stress closures that have been constructed for broader applications beyond shock-turbulence interactions.     

Master equation — The information gain minimum approximation

Approximate time-dependent solutions of a master equation having unique stationary solution can be obtained by minimizing the information gain functional subject to constraints for mean values of a number of chosen observables. We study mathematical properties of such an approximation. We find the region of applicability, prove that the approximate solutions are globally asymptotically stable, and show how the approximation is related to some exact integrodifferential equation governing the time evolution of the mean values of the chosen observables.     

An integrable 3D model for Heisenberg ferromagnetic spin system with biquadratic interactions

We formulate a model Hamiltonian to a 3D Ferromagnetic spin system incorporating biquadratic interactions. The dynamics is represented by a higher order (3+1) dimensional integrable nonlinear Schrödinger equation. We construct the Lax pair associate with the system and find multisoliton solutions using Darboux transformation(DT). We bilinearize the equation using Hirota’s bilinearization procedure and find one soliton solution.     

Deforming the Travelling Wave Solutions of the KdV Equation to Those of the Generalized Fifth Order KdV Equation

Using some limiting procedures, the solutions of the fifth order KdV equation u_t + (μu²+ υu_xx + αuu_xx + βu_x² + γu³ + δu_xxxx)_x = 0 would degenerate into the solutions of a simple equation, say KdV equation. In this letter, we analyze the possibility of the inverse procedure of the limiting process mentioned above for the travelling wave solutions. The results show that the procedure for deforming a travelling wave solution of the KdV equation to that of the generalized fifth order KdV equation can be accomplished by some pure algebraic tricks. Moreover, this inverse procedure is not unique in general.     

A new method for deriving analytical solutions of partial differential equations — Algebraically explicit analytical solutions of two-buoyancy natural convection in porous media

Analytical solutions of governing equations of various phenomena have their irreplaceable theoretical meanings. In addition, they can also be the benchmark solutions to verify the outcomes and codes of numerical solutions, and even to develop various numerical methods such as their differencing schemes and grid generation skills as well. A hybrid method of separating variables for simultaneous partial differential equation sets is presented. It is proposed that different methods of separating variables for different independent variables in the simultaneous equation set may be used to improve the solution derivation procedure, for example, using the ordinary separating method for some variables and using extraordinary methods of separating variables, such as the separating variables with addition promoted by the first author, for some other variables. In order to prove the ability of the above-mentioned hybrid method, a lot of analytical exact solutions of two-buoyancy convection in porous media are successfully derived with such a method. The physical features of these solutions are given. Supported by the National Natural Science Foundation of China (Grant No. 50576097) and the National Basic Research Development Program of China (Grant No. 2007CB206902)     

Rogue wave dynamics in barotropic relaxing media

In this work, we deal with a nonlinear wave equation, namely the Vakhnenko equation, which models the propagation of nonlinear wave in the barotropic relaxing media. Based on the homoclinic breather limit method, we seek rogue wave solution to the above equation. The results show that rogue wave or giant wave can exist in such a medium.     

Sixth order compact scheme combined with multigrid method and extrapolation technique for 2D poisson equation

We develop a sixth order finite difference discretization strategy to solve the two dimensional Poisson equation, which is based on the fourth order compact discretization, multigrid method, Richardson extrapolation technique, and an operator based interpolation scheme. We use multigrid V-Cycle procedure to build our multiscale multigrid algorithm, which is similar to the full multigrid method (FMG). The multigrid computation yields fourth order accurate solution on both the fine grid and the coarse grid. A sixth order accurate coarse grid solution is computed by using the Richardson extrapolation technique. Then we apply our operator based interpolation scheme to compute sixth order accurate solution on the fine grid. Numerical experiments are conducted to show the solution accuracy and the computational efficiency of our new method, compared to Sun–Zhang’s sixth order Richardson extrapolation compact (REC) discretization strategy using Alternating Direction Implicit (ADI) method and the standard fourth order compact difference (FOC) scheme using a multigrid method.     

Standard and goal-oriented adaptive mesh refinement applied to radiation transport on 2D unstructured triangular meshes

Standard and goal-oriented adaptive mesh refinement (AMR) techniques are presented for the linear Boltzmann transport equation. A posteriori error estimates are employed to drive the AMR process and are based on angular-moment information rather than on directional information, leading to direction-independent adapted meshes. An error estimate based on a two-mesh approach and a jump-based error indicator are compared for various test problems. In addition to the standard AMR approach, where the global error in the solution is diminished, a goal-oriented AMR procedure is devised and aims at reducing the error in user-specified quantities of interest. The quantities of interest are functionals of the solution and may include, for instance, point-wise flux values or average reaction rates in a subdomain. A high-order (up to order 4) Discontinuous Galerkin technique with standard upwinding is employed for the spatial discretization; the discrete ordinates method is used to treat the angular variable.     

Bilinearization and new multisoliton solutions for the (4+1)-dimensional Fokas equation

The (4 + 1)-dimensional Fokas equation is derived in the process of extending the integrable Kadomtsev–Petviashvili and Davey–Stewartson equations to higher-dimensional nonlinear wave equations. This equation is under investigation in this paper. Hirota’s bilinear method is, for the first time, used to solve such a higher-dimensional equation. In order to bilinearize the Fokas equation, some appropriate transformations are adopted. As a result, single-soliton solution, double-soliton solution and three-soliton solution are obtained. A new uniform formula of n-soliton solution is derived from this. It is shown that the transformations adopted in this work play a key role in converting the Fokas equation into Hirota’s bilinear form.     

Separation of variables in Kolmogorov's equation. I

Calculations are made of the coefficients of the nonstationary Kolmogorov equation for which a solution can be found by separation of variables using a complete set of differential symmetry operators of first order; the variables are separated in accordance with this procedure. A solution is found to the equation describing the motion of a dynamical system containing a random parameter.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 8, pp. 16–20, August, 1977.     

Ion-acoustic waves in plasma of warm ions and isothermal electrons using time-fractional KdV equation

The ion-acoustic solitary wave in collisionless unmagnetized plasma consisting of warm ions-fluid and isothermal electrons is studied using the time fractional KdV equation. The reductive perturbation method has been employed to derive the Korteweg-de Vries equation for small but finite amplitude ion-acoustic wave in warm plasma. The Lagrangian of the time fractional KdV equation is used in a similar form to the Lagrangian of the regular KdV equation with fractional derivative for the time differentiation. The variation of the functional of this Lagrangian leads to the Euler-Lagrange equation that gives the time fractional KdV equation. The variational-iteration method is used to solve the derived time fractional KdV equation. The calculations of the solution are carried out for different values of the time fractional order. These calculations show that the time fractional can be used to modulate the electrostatic potential wave instead of adding a higher order dissipation term to the KdV equation. The results of the present investigation may be applicable to some plasma environments,such as the ionosphere plasma.     

Ion-acoustic dressed soliton in electron-ion quantum plasma

Ion acoustic dressed soliton in an unmagnetized two-species electron-ion quantum plasma is studied. Using reductive perturbation technique, a higher order inhomogeneous (KdV-type) differential equation is derived for the second order correction. The nonsecular solution is obtained by using renormalization procedure. A new technique is used to obtain the particular solution of the higher order inhomogeneous equation which is found to be simpler compared to the technique used by previous investigators.