Explicit series solution of travelling waves with a front of Fisher equation

In this paper, an analytic technique, namely the homotopy analysis method, is employed to solve the Fisher equation, which describes a family of travelling waves with a front. The explicit series solution for all possible wave speeds 0 < c < +∞ is given. Such kind of explicit series solution has never been reported, to the best of author’s knowledge. Our series solution indicates that the solution contains an oscillation part when 0 < c < 2. The proposed analytic approach is general, and can be applied to solve other similar nonlinear travelling wave problems.     

A Bijective Approach to the Area of Generalized Motzkin Paths

For fixed positive integer k, let E_n denote the set of lattice paths using the steps (1, 1), (1, − 1), and (k, 0) and running from (0, 0) to (n, 0) while remaining strictly above the x-axis elsewhere. We first prove bijectively that the total area of the regions bounded by the paths of E_n and the x-axis satisfies a four-term recurrence depending only on k. We then give both a bijective and a generating function argument proving that the total area under the paths of E_n equals the total number of lattice points on the x-axis hit by the unrestricted paths running from (0, 0) to (n − 2, 0) and using the same step set as above.     

Coupled nonlinear waves in two-dimensional lattice

For one-dimensional nonlinear lattices, such as Toda lattice, it has been extensively studied. By considering the nonlinear effects of two-dimensional lattice, we set up the equation of motion for each particles (atoms, molecules or ions). For small amplitude and long wavelength nonlinear waves in this system, both the linear dispersion relation and the coupled Korteweg de Vries (KdV) equation are obtained. The simple soliton solution is obtained. If the nonlinear lattice is symmetric in the x and y directions, It is noted that there are two kinds of solitons. one is that propagates in either x or y directions, (1, 0) or (0, 1), the other is that propagates in the direction of (1, 1). It is in agreements with that of one-dimensional lattice. The different properties are investigated for different nonlinear interacting potentials, such as Toda potential, Morse potential and LJ potential.     

Computational modelling of fluid structure interaction at nano-scale boundaries with modified Maxwellian velocity distribution

The soft collisions among fluid–fluid and fluid-wall molecules are modeled from first principles. In particular, the assumption of Maxwellian distribution of velocities for thermalized molecules, in both parallel and perpendicular directions to the wall, has been re-evaluated with supporting experimental and/or numerical evidence.It is proposed that the normal component of molecular velocity post collision is conserved for all fluid molecules. The slip effect at the wall boundary, introduced by the surface roughness, is accounted by an accommodation coefficient f. A moving least square method is used to calculate macroscopic velocity values. The influence of molecular interaction on the macroscopic velocity distribution is investigated at 40 MPa and 300 K for slit pore, inclined and stepped wall configurations. The accommodation coefficient values f = 0, 0.07, 0.257, 0.45, 0.681 and 1; and acceleration values ranging from zero to 1 × 10¹¹ m/s² and 250 × 10¹¹ m/s² are used for comparison.The distribution of macroscopic velocity parallel to the wall is studied to observe the effect of the slip behaviour. The detailed study of average of velocity values at various magnitudes of acceleration has shown an evidence of characteristic low and high speed of molecular flows that is considered as significant and a comparison is sought with an equivalent laminar and turbulent flow style behaviour. The two dimensional vector and contour plots of macroscopic velocity provide further insights in understanding Continuum velocity distributions resulting from molecular fluid-wall interaction at nanoscale. The research has highlighted the need to develop molecular dynamics simulation techniques for non-periodic boundary conditions.     

Approximate solutions of a nonlinear oscillator typified as a mass attached to a stretched elastic wire by the homotopy perturbation method

The homotopy perturbation method is used to solve the nonlinear differential equation that governs the nonlinear oscillations of a system typified as a mass attached to a stretched elastic wire. The restoring force for this oscillator has an irrational term with a parameter λ that characterizes the system (0 ? λ ? 1). For λ = 1 and small values of x, the restoring force does not have a dominant term proportional to x. We find this perturbation method works very well for the whole range of parameters involved, and excellent agreement of the approximate frequencies and periodic solutions with the exact ones has been demonstrated and discussed. Only one iteration leads to high accuracy of the solutions and the maximal relative error for the approximate frequency is less than 2.2% for small and large values of oscillation amplitude. This error corresponds to λ = 1, while for λ < 1 the relative error is much lower. For example, its value is as low as 0.062% for λ = 0.5.     

Global chaos synchronization of new chaotic systems via nonlinear control

Nonlinear control is an effective method for making two identical chaotic systems or two different chaotic systems be synchronized. However, this method assumes that the Lyapunov function of error dynamic (e) of synchronization is always formed as V (e) = 1/2e^Te. In this paper, modification based on Lyapunov stability theory to design a controller is proposed in order to overcome this limitation. The method has been applied successfully to make two identical new systems and two different chaotic systems (new system and Lorenz system) globally asymptotically synchronized. Since the Lyapunov exponents are not required for the calculation, this method is effective and convenient to synchronize two identical systems and two different chaotic systems. Numerical simulations are also given to validate the proposed synchronization approach.     

Construction of analytic solution to chaotic dynamical systems using the Homotopy analysis method

Based on a new kind of analytic method, namely the Homotopy analysis method, an analytic approach to solve non-linear, chaotic system of ordinary differential equations is presented. The method is applied to Lorenz system; this system depends on the three parameters: σ, b and the so-called bifurcation parameter R are real constants. Two cases are considered. The first case is when R = 20.5 which corresponds to the transition region and the second case corresponds to R = 23.5 which corresponds to the chaotic region.The validity of the method is verified by comparing the approximation series solution with the results obtained using the standard numerical techniques such as Runge-Kutta method.     

Inverse problems for parabolic equations 2

Let u_t − u_xx = h(t) in 0 ⩽ x ⩽ π, t ⩾ 0. Assume that u(0, t) = v(t), u(π, t) = 0, and u(x, 0) = g(t). The problem is: what extra data determine the three unknown functions {h, v, g} uniquely? This question is answered and an analytical method for recovery of the above three functions is proposed.     

The applications of a higher-dimensional Lie algebra and its decomposed subalgebras

With the help of invertible linear transformations and the known Lie algebras, a higher-dimensional 6 × 6 matrix Lie algebra sμ(6) is constructed. It follows a type of new loop algebra is presented. By using a (2 + 1)-dimensional partial-differential equation hierarchy we obtain the integrable coupling of the (2 + 1)-dimensional KN integrable hierarchy, then its corresponding Hamiltonian structure is worked out by employing the quadratic-form identity. Furthermore, a higher-dimensional Lie algebra denoted by E, is given by decomposing the Lie algebra sμ(6), then a discrete lattice integrable coupling system is produced. A remarkable feature of the Lie algebras sμ(6) and E is used to directly construct integrable couplings.     

Heat transfer analysis of unsteady boundary layer flow by homotopy analysis method

This paper aims to present complete analytic solution to the unsteady heat transfer flow of an incompressible viscous fluid over a permeable plane wall. The flow is started due to an impulsively stretching porous plate. Homotopy analysis method (HAM) has been used to get accurate and complete analytic solution. The solution is uniformly valid for all time τ ∈ [0, ∞) throughout the spatial domain η ∈ [0, ∞). The accuracy of the present results is shown by giving a comparison between the present results and the results already present in the literature. This comparison proves the validity and accuracy of our present results. Finally, the effects of different parameters on temperature distribution are discussed through graphs.     

An extended rational interpolation method

To interpolate function, f(x), a ? x ? b, when we have some information about the values of f(x) and their derivatives in separate points on {x₀, x₁, … , x_n} ? [a, b], the Hermit interpolation method is usually used. Here, to solve this kind of problems, extended rational interpolation method is presented and it is shown that the suggested method is more efficient and suitable than the Hermit interpolation method, especially when the function f(x) has singular points in interval [a, b]. Also for implementing the extended rational interpolation method, the direct method and the inverse differences method are presented, and with some examples these arguments are examined numerically.     

Vacuum polarization in light two-electron atoms and ions

A general theory of the vacuum polarization in light atomic and muon-atomic systems is considered. We derive the closed analytical expression for the Uehling potential and evaluate corrections on vacuum polarization for the 1¹S-state of the two-electron ³He and ⁴He atoms and for some two-electron ions, including the Li⁺, Be²⁺, B³⁺ and C⁴⁺ ions. The correction for vacuum polarization in two-electron He atoms has been evaluated as ΔE_ueh ≈ 7.253 ± 0.0025 × 10⁻⁷ a.u. The analogous corrections in the two-electron He-like ions rapidly increase with the nuclear charge Q (ΔE_ueh ≈ 2.7061 × 10⁻⁶ a.u. for the Li⁺ ion and ΔE_ueh ≈ 2.3495 × 10⁻⁵ a.u. for the C⁴⁺ ion). The corresponding corrections have also been evaluated for the electron–nucleus and electron–electron interactions.     

A semi-Markov process with an inverse Gaussian distribution as sojourn time

Let X_n denote the state of a device after n repairs. We assume that the time between two repairs is the time τ taken by a Wiener process {W(t), t ? 0}, starting from w₀ and with drift μ < 0, to reach c ∈ [0, w₀). After the nth repair, the process takes on either the value X_n?1 + 1 or X_n?1 + 2. The probability that X_n = X_n?1 + j, for j = 1, 2, depends on whether τ ? t₀ (a fixed constant) or τ > t₀. The device is considered to be worn out when X_n ? k, where k ∈ {1, 2, …}. This model is based on the ones proposed by Rishel (1991) [1] and Tseng and Peng (2007) [2]. We obtain an explicit expression for the mean lifetime of the device. Numerical methods are used to illustrate the analytical findings.     

Multifractal scaling of crack images from pyroligneous acid dried sludge

Tannery effluent (sludge, wastewater) is treated by natural way. The waste sludge has been taken for two treatment process. The alkali chemicals are neutralized by pyroligneous acid which is obtained by pyrolysis process of wood. This sludge is sent out for drying. The dried sludge contains some crack pattern formation. Photographs were used to record two sludge cracking surfaces. Experiment has been performed to study the fracture pattern formation in thin film sludge. We studied changes of crack surface of a sludge by image analysis and also assessed applicability of fractal scaling to crack surfaces. The calculated crack surface dimension shows that the fracture surface exhibit fractal structure. Image size was 256 × 256 pixels. MFA (multifractal analysis) was carried out to the method of moments, i.e., the probability distribution was estimated for moments ranging from ?150 < q < 150 and the generalized dimension were calculated from the log/log slope of the probability distribution for the respective moments over box sizes. Generalized dimension D(q) were attained for this box size range, which are capable of characterizing heterogeneous spatial crack structure. Multifractal spectra analyzed two fracture surface of the image and results were indicated that the width of spectra increases due to pyroligneous acid. Multifractal method is sensitive enough to measure the fracture distribution and can be seen as a different approach for changing research of crack images of manure sludge drying.     

Characterizations of PG(n − 1, q)\PG(k − 1, q) by Numerical and Polynomial Invariants

We show that the simple matroid PG(n − 1, q)\PG(k − 1, q), for n ≥ 4 and 1 ≤ k ≤ n − 2, is characterized by a variety of numerical and polynomial invariants. In particular, any matroid that has the same Tutte polynomial as PG(n − 1, q)\PG(k − 1, q) is isomorphic to PG(n − 1, q)\PG(k − 1, q).     

Dynamics of landslide model with time delay and periodic parameter perturbations

In present paper, we analyze the dynamics of a single-block model on an inclined slope with Dieterich–Ruina friction law under the variation of two new introduced parameters: time delay T_d and initial shear stress μ. It is assumed that this phenomenological model qualitatively simulates the motion along the infinite creeping slope. The introduction of time delay is proposed to mimic the memory effect of the sliding surface and it is generally considered as a function of history of sliding. On the other hand, periodic perturbation of initial shear stress emulates external triggering effect of long-distant earthquakes or some non-natural vibration source. The effects of variation of a single observed parameter, T_d or μ, as well as their co-action, are estimated for three different sliding regimes: β < 1, β = 1 and β > 1, where β stands for the ratio of long-term to short-term stress changes. The results of standard local bifurcation analysis indicate the onset of complex dynamics for very low values of time delay. On the other side, numerical approach confirms an additional complexity that was not observed by local analysis, due to the possible effect of global bifurcations. The most complex dynamics is detected for β < 1, with a complete Ruelle–Takens–Newhouse route to chaos under the variation of T_d, or the co-action of both parameters T_d and μ. These results correspond well with the previous experimental observations on clay and siltstone with low clay fraction. In the same regime, the perturbation of only a single parameter, μ, renders the oscillatory motion of the block. Within the velocity-independent regime, β = 1, the inclusion and variation of T_d generates a transition to equilibrium state, whereas the small oscillations of μ induce oscillatory motion with decreasing amplitude. The co-action of both parameters, in the same regime, causes the decrease of block’s velocity. As for β > 1, highly-frequent, limit-amplitude oscillations of initial stress give rise to oscillatory motion. Also for β > 1, in case of perturbing only the initial shear stress, with smaller amplitude, velocity of the block changes exponentially fast. If the time delay is introduced, besides the stress perturbation, within the same regime, the co-action of T_d (T_d < 0.1) and small oscillations of μ induce the onset of deterministic chaos.     

Simulación numérica del flujo turbulento de aire con gotas dispersas de agua a través de separadores de torres de refrigeración

This work presents a numerical study on the turbulent flow of air with dispersed water droplets in separators of mechanical cooling towers. The averaged Navier-Stokes equations are discretised through a finite volume method, using the Fluent and Phoenics codes, and alternatively employing the turbulence models k ? ?, k ? ω and the Reynolds stress model, with low-Re version and wall enhanced treatment refinements. The results obtained are compared with numerical and experimental results taken from the literature. The degree of accuracy obtained with each of the considered models of turbulence is stated. The influence of considering whether or not the simulation of the turbulent dispersion of droplets is analyzed, as well as the effects of other relevant parameters on the collection efficiency and the coefficient of pressure drop. Focusing on four specific eliminators (‘Belgian wave’, ‘H1-V’, ‘L-shaped’ and ‘Zig-zag’), the following ranges of parameters are outlined: 1 ≤ U_e ≤ 5 m/s for the entrance velocity, 2 ≤ D_p ≤ 50 μm for the droplet diameter, 650 ≤ Re ≤ 8.500 for Reynolds number, and 0.05 ≤ P_i ≤ 5 for the inertial parameter. Results reached alternately with Fluent and Phoenics codes are compared. The best results correspond to the simulations performed with Fluent, using the SST k ? ω turbulence model, with values of the dimensionless scaled distance to wall y⁺ in the range 0.2 to 0.5. Finally, correlations are presented to predict the conditions for maximum collection efficiency (100 %), depending on the geometric parameter of removal efficiency of each of the separators, which is introduced in this work.     

An interval nonlinear program for the planning of waste management systems with economies-of-scale effects—A case study for the region of Hamilton,Ontario, Canada

In practical environmental systems with the effects of economies-of-scale (EOS), most relationships among different system components are nonlinear in nature, which can be described precisely only if a nonlinear model is employed. In this study, an interval nonlinear programming (INLP) model is developed and applied to the planning of a municipal solid waste (MSW) management system with EOS effects on system costs. The INLP has a nonlinear objective function and linear constraints. It handles nonlinearity presented as exponential functions. When exponential term p = 1 (in the INLP’s objective function), the model becomes an interval linear program; when p = 2, it becomes an interval quadratic program. Therefore, the INLP is flexible in reflecting a variety of system complexities. A solution algorithm with satisfactory performance is proposed. Application of the proposed method to the planning of waste management activities in the Hamilton-Wentworth Region, Ontario, Canada, indicated that reasonable solutions have been generated. In general, the INLP model could reflect uncertain and nonlinear characteristics of MSW management systems with EOS effects. The modeling results provided useful decision support for the Region’s waste management activities.     

Numerical simulation of magnetohydrodynamic buoyancy-induced flow in a non-isothermally heated square enclosure

Numerical simulation of magnetohydrodynamic (MHD) buoyancy-induced heat transfer and fluid flow has been analyzed in a non-isothermally heated square enclosure using finite volume method. The bottom wall of enclosure were heated and cooled with a sinusoidal function and top wall was cooled isothermally. Vertical walls of the enclosure were adiabatic. Effects of Rayleigh number (Ra = 10⁴, 10⁵ and 10⁶), Hartman number (Ha = 0, 50 and 100) and amplitude of sinusoidal function (n = 0.25, 0.5 and 1) on temperature and flow fields were analyzed. It was observed that heat transfer was decreased with increasing Hartmann number and decreasing value of amplitude of sinusoidal function.     

Free convection heat and mass transfer in a power law fluid past an inclined surface with thermophoresis

The problem of heat and mass transfer in a power law, two-dimensional, laminar, boundary layer flow of a viscous incompressible fluid over an inclined plate with heat generation and thermophoresis is investigated by the characteristic function method. The governing non-linear partial differential equations describing the flow and heat transfer problem are transformed into a set of coupled non-linear ordinary differential equation which was solved using Runge–Kutta shooting method. Exact solutions for the dimensionless temperature and concentration profiles, are presented graphically for different physical parameters and for the different power law exponents 0 < n < 0.5 and for n > 0.5.