Chapman–Enskog solutions to arbitrary order in Sonine polynomials V: Summational expressions for the viscosity-related bracket integrals

The Chapman–Enskog solutions of the Boltzmann equations provide a basis for the computation of important transport coefficients for both simple gases and gas mixtures. These coefficients include the viscosity, the thermal conductivity, and the diffusion coefficient. In a preceding paper on simple gases (I), we have shown that the use of higher-order Sonine polynomial expansions enables one to obtain results of arbitrary precision that are free of numerical error. In two subsequent papers (II–III), we extended our original simple gas work to encompass binary gas mixture computations of the viscosity, thermal conductivity, diffusion, and thermal diffusion coefficients to high-order. In a fourth paper (IV) we derived general summational representations for the diffusion- and thermal conductivity-related bracket integrals and provided compact, explicit expressions for all of these bracket integrals needed to compute the diffusion- and thermal conductivity-related transport coefficients up to order 5 in the Sonine polynomial expansions used. In all of this previous work we retained the full dependence of our solutions on the molecular masses, the molecular sizes, the mole fractions, and the intermolecular potential model via the omega integrals up to the final point of solution via matrix inversion. The elements of the matrices to be inverted are, in each case, determined by appropriate combinations of bracket integrals which contain, in general form, all of the various dependencies. Since accurate expressions for the needed bracket integrals have not previously been available in the literature beyond orders 2 or 3, and since such expressions are necessary for any extensive program of computations of the transport coefficients involving Sonine polynomial expansions to higher orders, we have investigated alternative methods of constructing appropriately general bracket integral expressions that do not rely on the term-by-term, expansion and pattern matching techniques that we developed for our previous work. It is our purpose in this paper to report the results of our efforts to obtain useful, alternative, general expressions for the bracket integrals associated with the viscosity-related Chapman–Enskog solutions for gas mixtures. Specifically, we have obtained such expressions in summational form that are conducive to use in high-order viscosity coefficient computations for arbitrary gas mixtures and have computed and reported explicit expressions for all of the orders up to 5.     

Chapman–Enskog solutions to arbitrary order in Sonine polynomials III: Diffusion,thermal diffusion,and thermal conductivity in a binary,rigid-sphere,gas mixture

The Chapman–Enskog solutions of the Boltzmann equation provide a basis for the computation of important transport coefficients for both simple gases and gas mixtures. These coefficients include the viscosity, the thermal conductivity, and the diffusion coefficient. In a preceding paper (I), for simple, rigid-sphere gases (i.e. single-component, monatomic gases) we have shown that the use of higher-order Sonine polynomial expansions enables one to obtain results of arbitrary precision that are free of numerical error and, in a second paper (II), we have extended our initial simple gas work to modeling the viscosity in a binary, rigid-sphere, gas mixture. In this latter paper we reported an extensive set of order 60 results which are believed to constitute the best currently available benchmark viscosity values for binary, rigid-sphere, gas mixtures. It is our purpose in this paper to similarly report the results of our investigation of relatively high-order (order 70), standard, Sonine polynomial expansions for the diffusion- and thermal conductivity-related Chapman–Enskog solutions for binary gas mixtures of rigid-sphere molecules. We note that in this work, as in our previous work, we have retained the full dependence of the solution on the molecular masses, the molecular sizes, the mole fractions, and the intermolecular potential model via the omega integrals. For rigid-sphere gases, all of the relevant omega integrals needed for these solutions are analytically evaluated and, thus, results to any desired precision can be obtained. The values of the transport coefficients obtained using Sonine polynomial expansions for the Chapman–Enskog solutions converge and, therefore, the exact diffusion and thermal conductivity solutions to a given degree of convergence can be determined with certainty by expanding to sufficiently high an order. We have used Mathematica® for its versatility in permitting both symbolic and high-precision computations. Our results also establish confidence in the results reported recently by other authors who used direct numerical techniques to solve the relevant Chapman–Enskog equations. While in all of the direct numerical methods more-or-less full calculations need to be carried out with each variation in molecular parameters, our work has utilized explicit, general expressions for the necessary matrix elements that retain the complete parametric dependence of the problem and, thus, only a matrix inversion at the final step is needed as a parameter is varied. This work also indicates how similar results may be obtained for more realistic intermolecular potential models and how other gas-mixture problems may also be addressed with some additional effort.     

Competitive modes for the Baier–Sahle hyperchaotic flow in arbitrary dimensions

In this article, the stretch-twist-fold (STF) flow is numerically studied using phase portraits, sensitive dependence on initial conditions, Lyapunov exponents, power spectrum, and the Poincaré map. The stretch-twist-fold flow is a two-parameter family of Stokes flows defined in a unit sphere that is associated with the fluid particle motion that naturally arises in the dynamo theory, which proposes a mechanism by which celestial bodies, such as earth and sun can maintain and amplify the magnetic field continuously. For this continuous growth of magnetic field, scientists are interested to invent new tools for the nonfuel consumption magnetism propulsion for the low earth orbit of spacecrafts or satellites. General properties of a chaotic dynamical system reference to the stretch-twist-fold flow model are addressed and numerical solutions are generated to explain some of these properties. Analytically, we studied the local behavior at equilibrium points. The predictability of chaos in the STF flow with the numerical calculation of Lyapunov exponents and Poincaré map is presented in this paper.     

Static equilibrium of a thermoelastic half-plane: Green’s functions and solutions in integrals

This article presents in a closed form new influence functions of a unit point heat source on the displacements for three boundary value problems of thermoelasticity for a half-plane. We also obtain the corresponding new integral formulas of Green’s and Poisson’s types that directly determine the thermoelastic displacements and stresses in the form of integrals of the products of specified internal heat sources or prescribed boundary temperature and constructed already thermoelastic influence functions (kernels). All these results are presented in terms of elementary functions in the form of three theorems. Based on these theorems and on derived early by author the general Green-type integral formula, we obtain in elementary functions new solutions to two particular boundary value problems of thermoelasticity for half-plane. The graphical presentation of the temperature and thermal stresses of one concrete boundary value problems of thermoelasticity for half-plane also is included. The proposed method of constructing thermoelastic Green’s functions and integral formulas is applicable not only for a half-plane, but also for many other two- and three-dimensional canonical domains of different orthogonal coordinate systems.     

Nonlinear wave solutions and their relations for the modified Benjamin–Bona–Mahony equation

Nonlinear Dynamics - In this paper, we use the bifurcation method of dynamical systems to investigate the nonlinear wave solutions of the modified Benjamin–Bona–Mahony equation. These...     

The higher-order lump,breather and hybrid solutions for the generalized Konopelchenko–Dubrovsky–Kaup–Kupershmidt equation in fluid mechanics

Nonlinear Dynamics - Under investigation in this paper is the (2 + 1)-dimensional generalized Konopelchenko–Dubrovsky–Kaup–Kupershmidt equation, which can be utilized to describe...     

Exact harmonic solutions to Guyer–Krumhansl-type equation and application to heat transport in thin films

A system of hyperbolic-type inhomogeneous differential equations (DE) is considered for non-Fourier heat transfer in thin films. Exact harmonic solutions to Guyer–Krumhansl-type heat equation and to the system of inhomogeneous DE are obtained in Cauchy- and Dirichlet-type conditions. The contribution of the ballistic-type heat transport, of the Cattaneo heat waves and of the Fourier heat diffusion is discussed and compared with each other in various conditions. The application of the study to the ballistic heat transport in thin films is performed. Rapid evolution of the ballistic quasi-temperature component in low-dimensional systems is elucidated and compared with slow evolution of its diffusive counterpart. The effect of the ballistic quasi-temperature component on the evolution of the complete quasi-temperature is explored. In this context, the influence of the Knudsen number and of Cauchy- and Dirichlet-type conditions on the evolution of the temperature distribution is explored. The comparative analysis of the obtained solutions is performed.     

On the directional approach in constitutive modelling: A general thermomechanical framework and exact solutions for Mooney–Rivlin type elasticity in each direction

In order to represent process-induced anisotropies in continuum mechanics or to transfer one-dimensional material models to three spatial dimensions the directional approach is a helpful technique. Since the essential equations are defined in the orientation space it is also denoted as microsphere approach. In the current article, the relation for the directional stress tensor of the second Piola–Kirchhoff type is motivated using the volumetric/isochoric split of the deformation gradient and the Clausius–Duhem inequality. Owing to inherent nonlinearities, numerical discretisation techniques are usually applied to calculate the total stress by averaging the directional stress tensors over the unit sphere. In order to investigate the accuracy of such simulations, the availability of exact solutions in closed form is essential. To this end, the tension/compression behaviour which belongs to a certain direction in the orientation space is modelled by an elasticity relation of the Mooney Rivlin type. The exact solutions are calculated, visualized and discussed for uniaxial tension and compression as well as for equibiaxial tension.     

Explicit solutions for the coupled stretching–bending problems of holes in composite laminates

Although the classical lamination theory was developed long time ago, it is still not easy to apply this theory to find the analytical solutions for the curvilinear boundary value problems especially when the stretching and bending are coupled each other. To overcome the difficulties, recently we developed a Stroh-like formalism for the general composite laminates. By using this formalism, most of the relations for the coupled stretching–bending problems can be organized into the forms of Stroh formalism for two-dimensional anisotropic elasticity problems. With this newly developed Stroh-like formalism, it becomes easier to obtain an analytical solution for the coupled stretching–bending problems of holes in composite laminates. Because the Stroh-like formalism is a complex variable formalism, the analytical solutions for the whole field are expressed in complex form. Through the use of some identities derived in this paper, the resultant forces and moments around the hole boundary are obtained explicitly in real form. Due to the lack of analytical solutions for the general cases, the comparison is made with the existing analytical solutions for some special cases. In addition, to show the generality of our analytical solutions, several numerical examples are presented to discuss the coupling effect of the laminates and the shape effect of the holes.     

Soliton solutions for the reduced Maxwell–Bloch system in nonlinear optics via the N-fold Darboux transformation

Under investigation in this paper is the reduced Maxwell?CBloch system, which describes the propagation of the intense ultra-short optical pulses through a two-level dielectric medium. Through symbolic computation, conservation laws are derived and N-fold Darboux transformation (DT) is constructed for that system. By virtue of the DT obtained, multi-soliton solutions are generated. Figures are plotted to reveal the following dynamic features of the solitons: (1) Elastic interactions between two bright one-peak solitons, between two bight two-peak solitons and between two dark two-peak solitons; (2) Parallel propagations between two bright one-peak solitons, between two bright two-peak solitons and between two dark two-peak solitons; (3) Periodic propagations of hump solitons, of a pair of bound hump solitons with the same amplitude and of dark solitons.     

Novel characteristics of lump and lump–soliton interaction solutions to the generalized variable-coefficient Kadomtsev–Petviashvili equation

Nonlinear Dynamics - With the inhomogeneities of media taken into account, a generalized variable-coefficient Kadomtsev–Petviashvili (vcKP) equation is proposed to model nonlinear waves in...     

On the evolutions of the discontinuities of orders ⩾ m in solutions for quasi-linear systems of PDEs of order m,or for the linearizations of these systems

Summary One considers a system L[u]=0 of PDEs, quasi-linear (according to [1]) and of order m, which possesses a bicharacteristic line , as it happens in the hyperbolic case. For v=0, , –m (>0) let u(^v) be a discontinuity wave of order m+v that solves the system above and whose discontinuity hypersurface includes . The corresponding transport equations along are considered. Furthermore some interesting cases are pointed out, in which these equations turn out to be mutually equivalent in a suitable sense. Some theorems are stated to compare the transport equations for the discontinuities of the above kinds, that are connected with the systems d^hL[u]/dt^h=0 (h=0, , –m) and/or the linearization of the system L[u]=0 around any regular solution of it.

---

Sommario Si considera un sistema L[u]=0 di equazioni alle derivate parziali, quasi lineare (secondo [1]) e di ordine m, il quale sia dotato di qualche bicaratteristica , come accade nel caso iperbolico. Per v=0, , –m(>0) sia u(^v) un'onda di discontinuità di ordine m+v risolvente il detto sistema e avente ipersuperficie di discontinuità contenente Si considerano le relative equazioni di trasporto lungo e si determinano casi interessanti in cui queste equazioni sono mutuamente equivalenti in senso opportuno. Si stabiliscono teoremi di confronto per il trasporto delle discontinuità del tipo suddetto, relative ai sistemi d^hL[u]/dt^h=0 (h=0, , –m) e/o alla linearizazione del sistema L[u]=0 attorno a qualche sua soluzione regolare.

     

An empirical model for the thermal conductivity of compacted bentonite and a bentonite–sand mixture

The thermal conductivities of compacted bentonite and a bentonite–sand mixture were measured to investigate the effects of dry density, water content and sand fraction on the thermal conductivity. A single expression has been proposed to describe the thermal conductivity of the compacted bentonite and the bentonite–sand mixture once their primary parameters such as dry density, water content and sand fraction are known.     

Darboux transformation,exact solutions and conservation laws for the reverse space-time Fokas–Lenells equation

Nonlinear Dynamics - In this paper, we firstly deduce a reverse space-time Fokas–Lenells equation which can be derived from a rather simple but extremely important symmetry reduction of...     

Correction to: Darcy–Brinkman–Forchheimer Model for Film Boiling in Porous Media

Transport in Porous Media - A correction to this paper has been published: https://doi.org/10.1007/s11242-021-01631-0     

Interaction solutions between lump and stripe soliton to the (2+1)-dimensional Date–Jimbo–Kashiwara–Miwa equation

Nonlinear Dynamics - In this paper, we apply the ansatz method to the multi-linear form of the (2+1)-dimensional Date–Jimbo–Kashiwara–Miwa equation for constructing interaction...     

Eulerian–Lagrangian bridge for the energy and dissipation spectra in isotropic turbulence

We study, numerically and analytically, the relationship between the Eulerian spectrum of kinetic energy, E _E(k, t), in isotropic turbulence and the corresponding Lagrangian frequency energy spectrum, E _L(ω, t), for which we derive an evolution equation. Our DNS results show that not only E _L(ω, t) but also the Lagrangian frequency spectrum of the dissipation rate ${\varepsilon_{\rm L} (\omega, t)}$ has its maximum at low frequencies (about the turnover frequency of energy-containing eddies) and decays exponentially at large frequencies ω (about a half of the Kolmogorov microscale frequency) for both stationary and decaying isotropic turbulence. Our main analytical result is the derivation of equations that bridge the Eulerian and Lagrangian spectra and allow the determination of the Lagrangian spectrum, E _L (ω) for a given Eulerian spectrum, E _E (k), as well as the Lagrangian dissipation, ${\varepsilon_{\rm L}(\omega)}$ , for a given Eulerian counterpart, ${\varepsilon_{\rm E} (k)=2\nu k^2 E_{\rm E}(k)}$ . These equations were derived from the Navier–Stokes equations in the sweeping-free coordinate system (intermediate between the Eulerian and Lagrangian frameworks) which eliminates the effect of the kinematic sweeping of the small eddies by the larger eddies. We show that both analytical relationships between E _L (ω) and E _E (k) and between ${\varepsilon_{\rm L} (\omega)}$ and ${\varepsilon_{\rm E} (k)}$ are in very good quantitative agreement with our DNS results and explain how ${\varepsilon_{\rm L} (\omega, t)}$ has its maximum at low frequencies and decays exponentially at large frequencies.