Search

A special numerical technique has been developed for identification of solitary wave solutions of Boussinesq and Korteweg–de Vries equations. Stationary localized waves are considered in the frame moving to the right. The original ill-posed problem is transferred into a problem of the unknown coefficient from over-posed boundary data in which the trivial solution is excluded. The Method of Variational Imbedding is used for solving the inverse problem. The generalized sixth-order Boussinesq equation is considered for illustrations.

Generalized Wave Equations containing dispersion, dissipation and energy-production (GDWE) are considered in lieu of dissipative NEE as more suitable models for two-way interaction of localized waves. The quasi-particle behavior and the long-time evolution of localized solutions upon take-over and head-on collisions are investigated numerically by means of an adequate difference scheme which represents faithfully the balance/conservation laws. It is shown that in most cases the balance between energy production/dissipation and nonlinearity plays a similar role to the classical Boussinesq balance between dispersion and nonlinearity, namely it can create and support localized solutions which behave as quasi-particles upon collisions and for a reasonably long time after that.

The inverse problem of identification of boundary-layer thickness is replaced by the higher-order boundary value problem for the Euler–Lagrange equations for minimization of the quadratic functional of the original system (Method of Variational Imbedding – MVI). The imbedding problem is correct in the sense of Hadamard and consists of an explicit differential equation for the boundary-layer thickness. The existence and uniqueness of solution of the linearized imbedding problem is demonstrated and a difference scheme of splitting type is proposed for its numerical solution. The performance of the technique is demonstrated for three different boundary-layer problems: the Blasius problem, flow in the vicinity of plane stagnation point and the flow in the leading stagnation point on a circular cylinder. Comparisons with the self-similar solutions where available are quantitatively very good.

Pattern formation in a finite layer of fluid induced either by buoyancy or by a surface-tension gradient is considered. The fluid is confined between poor conducting horizontal boundaries, leading to patterns with a characteristic horizontal scale much larger than the fluid depth. The evolution of the system is studied by numerical integration of the (1+2)D equation introduced by Knobloch [1990]: [Formula: see text] Here µ is the scaled bifurcation parameter, κ=1, and a represents the effects of a heat transfer finite Biot number. The coefficients β, δ and γ do not vanish when the boundary conditions at top and bottom are not identical (β≠0, δ≠0) or when non-Boussinesq effects are taken into account (γ≠0). When the conductive state becomes unstable due to surface-tension inhomogeneities, it is shown that the system evolves towards a stationary pattern of hexagons with up or down flow depending on the relative value of the coefficients β and δ. In the case of buoyancy-driven convection (β=δ≠0), the system displays a tesselation of steady squares. Knobloch’s equation also describes time-dependent patterns at high thermal gradients, including spatio-temporal chaos, due to the non-variational character of the equation.

Two improved versions of Boussinesq equation (Boussinesq paradigm) have been considered which are well-posed (correct in the sense of Hadamard) as an initial value problem: the Proper Boussinesq Equation (PBE) and the Regularized Long Wave Equation (RLWE). Fully implicit difference schemes have been developed strictly representing, on difference level, the conservation or balance laws for the mass, pseudoenergy or pseudomomentum of the wave system. Thresholds of possible nonlinear blow-up are identified for both PBE and RLWE. The head-on collisions of solitary waves of the sech type (Boussinesq solitons) have been investigated. They are subsonic and negative (surface depressions) for PBE and supersonic and positive (surface elevations) for RLWE. The numerically recovered sign and sizes of the phase shifts are in very good quantitative agreement with analytical results for the two-soliton solution of PBE. The subsonic surface elevations are found to be not permanent but to gradually transform into oscillatory pulses whose support increases and amplitude decreases with time although the total pseudoenergy is conserved within 10−10. The latter allows us to claim that the pulses are solitons despite their “aging” (which is felt on times several times the time-scale of collision). For supersonic phase speeds, the collision of Boussinesq solitons has inelastic character exhibiting not only a significant phase shift but also a residual signal of sizable amplitude but negligible pseudoenergy. The evolution of the residual signal is investigated numerically for very long times.

We consider a random two-phase medium which represents a matrix containing an array of allowed to overlap spherical inclusions with random radii. A statistical theory of transport phenomena in the medium, on the example of heat propagation, is constructed by means of the functional (Volterra-Wiener) series approach. The functional series for the temperature is rendered virial in the sense that its truncation after the p-tuple term yields results for all multipoint correlation functions of the temperature field that are asymptotically correct to the order np, where n is the mean number of spheres per unit volume. The case p=2 is considered in detail and the needed kernels of the factorial series are found to the order n2. In this way not only the effective conductivity, but also the full statistical solution, i.e., all needed correlation functions, can be expressed in a closed form, correct to the said order.