Turbulence and Bubble Break up in Slug Flow with Wall Injection

The present work investigates experimentally the changes on the properties of horizontal slug flows subject to fluid injection at the wall. Measurements include data on global flow rates, pressure drop and local mean and fluctuating velocity profiles for nine different conditions. The properties of the two-phase flow are measured through a Shadow Sizer system and laser-based sensors. Two distinct flow transpiration rates are studied, \(v_{wi}^{++}\) = v _w / U _m = 0.0005 and 0.001. The effects of flow transpiration were observed to induce bubble break-up and large changes in the passage frequency and characteristic lengths of the unit cells. In addition to the two-phase flow results, single-phase flow measurements are presented with a view to compare the different turbulent effects introduced by the second phase. The work also proposes modifications in the models of Dukler and Hubbard (Ind. Eng. Chem. Fund. 14 337–347 (1975)) and Orell (Chem. Eng. Sci. 60 1371–1381 (2005)) so that fluid injection at the wall can be accounted for. All theoretical predictions are compared with the experimental data.     

Near-Wall Stress Balance in Front of a Wall-Mounted Cylinder

The stress balance in the near-wall flow in front of a cylinder mounted on a flat plate at moderate Reynolds number is investigated by applying highly resolved Large-Eddy Simulation (LES). The flow around wall-mounted bluff bodies is subject of research due to its wide relevance for engineering applications. However, the structure of the vortex system in front of such a bluff body is complex, bears strong velocity and pressure gradients in each spatial direction and has rich dynamics. Furthermore, the vortex system is located close to the investigated flat bottom wall (Dargahi, Exp. Fluids 8(1-2):1–12, 1989; Devenport and Simpson, J. Fluid Mech. 210:23–55, 1990). Thus, classical models for the treatment of the near-wall flow based on the logarithmic law of the wall or a power law cannot be expected to suffice in such kind of flow (Pope 2011). This paper assesses which contributors to the stress balance have significant influence on the balances residual and thus have to be considered by an approach to model the investigated near-wall flow. To do so, the momentum equation in streamwise direction is integrated in wall-normal direction and applied to the results gained from the LES. The evaluation of the stress balance along four selected wall-normal profiles indicates that the significance of each single term depends on where the profile is located. Outside the viscous layer, no term except the viscous stresses can be neglected in general. The amplitude of the pressure gradient as well as horizontal gradients of mean and fluctuating velocity are multiples of the estimated wall shear stress. Wall models not including a spatial approach are therefore most likely to fail in such kind of flow.     

Two-subsystem thermodynamics for the mechanics of aging amorphous solids

The effect of physical aging on the mechanics of amorphous solids as well as mechanical rejuvenation is modeled with nonequilibrium thermodynamics, using the concept of two thermal subsystems, namely a kinetic one and a configurational one. Earlier work (Semkiv and Hütter in J Non-Equilib Thermodyn 41(2):79–88, 2016) is extended to account for a fully general coupling of the two thermal subsystems. This coupling gives rise to hypoelastic-type contributions in the expression for the Cauchy stress tensor, that reduces to the more common hyperelastic case for sufficiently long aging. The general model, particularly the reversible and irreversible couplings between the thermal subsystems, is compared in detail with models in the literature (Boyce et al. in Mech Mater 7:15–33, 1988; Buckley et al. in J Mech Phys Solids 52:2355–2377, 2004; Klompen et al. in Macromolecules 38:6997–7008, 2005; Kamrin and Bouchbinder in J Mech Phys Solids 73:269–288 2014; Xiao and Nguyen in J Mech Phys Solids 82:62–81, 2015). It is found that only for the case of Kamrin and Bouchbinder (J Mech Phys Solids 73:269–288, 2014) there is a nontrivial coupling between the thermal subsystems in the reversible dynamics, for which the Jacobi identity is automatically satisfied. Moreover, in their work as well as in Boyce et al. (Mech Mater 7:15–33, 1988), viscoplastic deformation is driven by the deviatoric part of the Cauchy stress tensor, while for Buckley et al. (J Mech Phys Solids 52:2355–2377, 2004) and Xiao and Nguyen (J Mech Phys Solids 82:62–81, 2015) this is not the case.     

Passive Scalar Transport Modeling for Hybrid RANS/LES Simulation

A transport model for hybrid RANS/LES simulation of passive scalars is proposed. It invokes a dynamically computed subgrid Prandtl number. The method is based on computing test-filter fluxes. The formulation proves to be especially effective on coarse grids, as occur in DES. After testing it in a wall resolved LES, the present formulation is applied to the Adaptive DDES model of Yin et al. (Phys. Fluids 27, 025105 2015). It is validated by turbulent channel flow and turbulent boundary layer computations.     

Strichartz Estimates and Moment Bounds for the Relativistic Vlasov–Maxwell System

We consider the relativistic Vlasov–Maxwell system with data of unrestricted size and without compact support in momentum space. In the two-dimensional and the two-and-a-half-dimensional cases, Glassey–Schaeffer proved (Commun Math Phys 185:257–284, 1997; Arch Ration Mech Anal 141:331–354, 1998; Arch Ration Mech Anal. 141:355–374, 1998) that for regular initial data with compact momentum support this system has unique global in time classical solutions. In this work we do not assume compact momentum support for the initial data and instead require only that the data have polynomial decay in momentum space. In the two-dimensional and the two-and-a-half-dimensional cases, we prove the global existence, uniqueness and regularity for solutions arising from this class of initial data. To this end we use Strichartz estimates and prove that suitable moments of the solution remain bounded. Moreover, we obtain a slight improvement of the temporal growth of the \({L^\infty_x}\) norms of the electromagnetic fields compared to Glassey and Schaeffer (Commun Math Phys 185:257–284, 1997; Arch Ration Mech Anal 141:355–374, 1998). In the three-dimensional case, we apply Strichartz estimates and moment bounds to show that a regular solution can be extended as long as \({{\|p_0^{\theta} f \|_{L^{q}_{x}L^1_{p}}}}\) remains bounded for \({\theta > \frac{2}{q}}\), \({2 < q \leqq \infty}\). This improves previous results of Pallard (Indiana Univ Math J 54(5):1395–1409, 2005; Commun Math Sci 13(2):347–354, 2015).     

Thermodynamic model formulation for viscoplastic solids as general equations for non-equilibrium reversible–irreversible coupling

Thermodynamic models for viscoplastic solids are often formulated in the context of continuum thermodynamics and the dissipation principle. The purpose of the current work is to show that models for such material behavior can also be formulated in the form of a General Equation for Non-Equilibrium Reversible–Irreversible Coupling (GENERIC), see, e.g., Grmela and Öttinger (Phys Rev E, 56:6620–6632, 1997), Öttinger and Grmela (Phys Rev E, 56:6633–6655, 1997), Grmela (J Non-Newtonian Fluid Mech, 165:980–986, 2010). A GENERIC combines Hamiltonian-dynamics-based modeling of time-reversible processes with Onsager–Casimir-based modeling of time-irreversible processes. The result is a model for the approach of non-equilibrium systems to thermodynamic equilibrium. Originally developed to model complex fluids, it has recently been applied to anisotropic inelastic solids in Eulerian (Hütter and Tervoort, in J Non-Newtonian Fluid Mech, 152:45–52, 2008; Hütter and Tervoort, in J Non-Newtonian Fluid Mech, 152:53–65, 2008; Hütter and Tervoort, in Adv Appl Mech, 42:254–317, 2008) and Lagrangian (Hütter and Svendsen, in J Elast 104:357–368, 2011) settings, as well as to damage mechanics. For simplicity, attention is focused in the current work on the case of thermoelastic viscoplasticity. Central to this formulation is a GENERIC-based form for the viscoplastic flow rule. A detailed comparison with the formulation based on continuum thermodynamics and the dissipation principle is given.     

Topological and Measure-Theoretical Entropies of Nonautonomous Dynamical Systems

We consider the dynamics of a nonautonomous dynamical system determined by a sequence of continuous self-maps \(f_n:X \rightarrow X,\) where \( n \in {\mathbb {N}},\) defined on a compact metric space X. Applying the theory of the Carathéodory structures, elaborated by Pesin (Dimension Theory in Dynamical Systems. Chicago Lectures in Mathematics. The University of Chicago Press, Chicago, 1997), we construct a Carathéodory structure whose capacity coincides with the topological entropy of the considered system. Generalizing the notion of local measure entropy, introduced by Brin and Katok (in: Palis (ed) Geometric Dynamics, Lecture Notes in Mathematics. Springer, Berlin 1983) for a single map, to a nonautonomous dynamical system we provide a lower and upper estimations of the topological entropy by local measure entropies. The theorems of the paper generalize results of Kawan (Nonautonomous Stoch Dyn Syst 1:26–52, 2013) and of Feng and Huang (J Funct Anal 263:2228–2254, 2012). Also, we construct a new entropy-like invariant such the entropy of a sequence \(\{f_n:X \rightarrow X\}_{n=1}^{\infty }\) of Lipschitz continuous maps with the same Lipschitz constant \(L >1,\) restricted to a subset \(Y\subset X,\) is upper bounded by Hausdorff dimension of Y multiplied by the logarithm of the Lipschitz constant L. This gives a generalizations of results of Dai et al. (Sci China Ser A 41:1068–1075, 1998) and Misiurewicz (Discret Contin Dyn Syst 10:827–833, 2004).     

On a Nonlinear Model for Tumor Growth in a Cellular Medium

We investigate the dynamics of a nonlinear model for tumor growth within a cellular medium. In this setting the “tumor” is viewed as a multiphase flow consisting of cancerous cells in either proliferating phase or quiescent phase and a collection of cells accounting for the “waste” and/or dead cells in the presence of a nutrient. Here, the tumor is thought of as a growing continuum \(\Omega \) with boundary \(\partial \Omega \) both of which evolve in time. In particular, the evolution of the boundary \(\partial \Omega \) is prescibed by a given velocity \({{{\varvec{V}}}.}\) The key characteristic of the present model is that the total density of cancerous cells is allowed to vary, which is often the case within cellular media. We refer the reader to the articles (Enault in Mathematical study of models of tumor growth, 2010; Li and Lowengrub in J Theor Biol, 343:79–91, 2014) where compressible type tumor growth models are investigated. Global-in-time weak solutions are obtained using an approach based on penalization of the boundary behavior, diffusion, viscosity and pressure in the weak formulation, as well as convergence and compactness arguments in the spirit of Lions (Mathematical topics in fluid dynamics. Compressible models, 1998) [see also Donatelli and Trivisa (J Math Fluid Mech 16: 787–803, 2004), Feireisl (Dynamics of viscous compressible fluids, 2014)].     

Bail-Down Test Simulation at Laboratory Scale

This paper presents a comparison of hydraulic oil conductivity obtained from interpreting bail-down test data to values calculated from theory. The bail-down tests were performed at laboratory scale, on a radial portion of a circular domain filled with calibrated sand allowing hydraulic oil conductivity to be calculated using Parker’s theoretical model (Parker et al. in Water Resour Res 23(4):618–624, 1987). The bail-down tests were interpreted using the modified Bouwer and Rice (Huntley in Ground Water 38(1):46–52, 2000) and the modified Cooper methods (Beckett and Lyverse in API Interact LNAPL Guide 2:1–27, 2002). The results show that (1) both interpretation methods from bail-down test data give similar hydraulic oil conductivities, and (2) the hydraulic oil conductivities estimated from bail-down test data agree well with the hydraulic oil conductivity predicted when using the Parker theoretical model. Overall, this paper confirms that the modified Bouwer and Rice (Huntley 2000) and the modified Cooper methods (Beckett and Lyverse 2002) are valid to estimate hydraulic oil conductivity, giving realistic values despite test conditions not meeting all the assumptions and boundary conditions of each analytical solution.     

Shock wave emission during the collapse of cavitation bubbles

Shock wave emission induced by intense laser pulses is investigated experimentally. The present work focuses on the conditions of shock wave emission in glycerine and distilled water during the first bubble collapse. Experimental investigations are carried out in liquids as a function of temperature and viscosity. Comparison is made with the theoretical work of Poritsky (Proc 1st US Natl Congress Appl Mech 813–821, 1952) and Brennen (Cavitation and bubble dynamics, Oxford University Press 1995). To the best knowledge of the authors, this is the first experimental verification of those theories.     

On the Generalization of Reissner Plate Theory to Laminated Plates,Part II: Comparison with the Bending-Gradient Theory

In the first part of this two-part paper (Lebée and Sab in On the generalization of Reissner plate theory to laminated plates, Part I: theory, doi: 10.1007/s10659-016-9581-6, 2015), the original thick plate theory derived by Reissner (J. Math. Phys. 23:184–191, 1944) was rigorously extended to the case of laminated plates. This led to a new plate theory called Generalized-Reissner theory which involves the bending moment, its first and second gradients as static variables. In this second paper, the Bending-Gradient theory (Lebée and Sab in Int. J. Solids Struct. 48(20):2878–2888, 2011 and 2889–2901, 2011) is obtained from the Generalized-Reissner theory and several projections as a Reissner–Mindlin theory are introduced. A comparison with an exact solution for the cylindrical bending of laminated plates is presented. It is observed that the Generalized-Reissner theory converges faster than the Kirchhoff theory for thin plates in terms of deflection. The Bending-Gradient theory does not converge faster but improves considerably the error estimate.     

On the Generalization of Reissner Plate Theory to Laminated Plates,Part I: Theory

This is the first part of a two-part paper presenting the generalization of Reissner thick plate theory (Reissner in J. Math. Phys. 23:184–191, 1944) to laminated plates and its relation with the Bending-Gradient theory (Lebée and Sab in Int. J. Solids Struct. 48(20):2878–2888, 2011 and in Int. J. Solids Struct. 48(20):2889–2901, 2011). The original thick and homogeneous plate theory derived by Reissner (J. Math. Phys. 23:184–191, 1944) is based on the derivation of a statically compatible stress field and the application of the principle of minimum of complementary energy. The static variables of this model are the bending moment and the shear force. In the present paper, the rigorous extension of this theory to laminated plates is presented and leads to a new plate theory called Generalized-Reissner theory which involves the bending moment, its first and second gradients as static variables. When the plate is homogeneous or functionally graded, the original theory from Reissner is retrieved. In the second paper (Lebée and Sab, 2015), the Bending-Gradient theory is obtained from the Generalized-Reissner theory and comparison with an exact solution for the cylindrical bending of laminated plates is presented.     

An Analytical Solution for Capillary Gravity Drainage with Dominant Viscous Forces

The development of the capillary fringe during gravity drainage has a significant influence on saturation and pressure distributions in porous formations (Sarkarfarshi et al. in Int J Greenh Gas Control 23:61–71, 2014). This paper introduces an analytical solution for gravity drainage in an axisymmetric geometry with significant capillary pressure. The drainage process results from the injection of a lighter and less viscous injectant into a porous medium saturated with a heavier and more viscous pore fluid. If the viscous force dominates the capillary and the buoyancy forces, then the flow regime is approximated by differential equations and the admissible solution comprises a front shock wave and a trailing simple wave. In contrast to existing analytical solutions for capillary gravity drainage problems (e.g., Nordbotten and Dahle in 47(2) 2011; Golding et al. in J Fluid Mech 678:248–270 2011), this solution targets the saturation distribution during injection at an earlier point in time. Another contribution of this analytical solution is the incorporation of a completely drained flow regime close to the injection well. The analytical solution demonstrates the strong dependency of the saturation distribution upon relative permeability functions, gas entry capillary pressure, and residual saturation. The analytical results are compared to results from a commercial reservoir engineering software package (\(\hbox {CMG } \hbox {STARS}^{\mathrm{TM}}\)).     

The peculiar elongational viscosity of concentrated solutions of monodisperse PMMA in oligomeric MMA

Concentrated solutions of nearly monodisperse poly(methyl methacrylate), PMMA-270k and PMMA-86k, in oligo(methyl methacrylate), MMA o-4k and MMA o-2k, investigated by Wingstrand et al. (Phys Rev Lett 115:078302, 2015) and Wingstrand (2015) do not follow the linear-viscoelastic scaling relations of monodisperse polystyrenes (PS) dissolved in oligomeric styrene (Wagner in Rheol Acta 53:765–777, 2014a, in Non-Newtonian Fluid Mech 222:121–131, 2014b; Wagner et al. in J Rheol 59:1113–1130, 2015). Rather, PMMA-270k shows an attractive interaction with MMA, in contrast to the interaction of PMMA-86k and MMA. This different behavior can be traced back to different tacticities of the two polymers. The attractive interaction of PMMA-270k with o-4k creates pseudo entanglements, which increase the interchain tube pressure, and therefore, the solution PMMA-270k/o-4k shows, as reported by Wingstrand et al. (Phys Rev Lett 115:078302, 2015), qualitatively a similar scaling of the elongational viscosity with \( {\left(\dot{\varepsilon}{\tau}_R\right)}^{-1/2} \) as observed for polystyrene melts. For the solution PMMA-270/o-2k, this effect is only seen at the highest elongation rates investigated. The elongational viscosity of PMMA-86k dissolved in oligomeric MMA is determined by the Rouse time of the melt, as in the case of polystyrene solutions.

Open image in new window

Graphical abstract ?

     

Linear Pantographic Sheets: Existence and Uniqueness of Weak Solutions

The well-posedness of the boundary value problems for second gradient elasticity has been studied under the assumption of strong ellipticity of the dependence on the second placement gradients (see, e.g., Chambon and Moullet in Comput. Methods Appl. Mech. Eng. 193:2771–2796, 2004 and Mareno and Healey in SIAM J. Math. Anal. 38:103–115, 2006).The study of the equilibrium of planar pantographic lattices has been approached in two different ways: in dell’Isola et al. (Proc. R. Soc. Lond. Ser. A 472:20150, 2016) a discrete model was introduced involving extensional and rotational springs which is also valid in large deformations regimes while in Boutin et al. (Math. Mech. Complex Syst. 5:127–162, 2017) the lattice has been modelled as a set of beam elements interconnected by internal pivots, but the analysis was restricted to the linear case. In both papers a homogenized second gradient deformation energy, quadratic in the neighbourhood of non deformed configuration, is obtained via perturbative methods and the predictions obtained with the obtained continuum model are successfully compared with experiments. This energy is not strongly elliptic in its dependence on second gradients. We consider in this paper also the important particular case of pantographic lattices whose first gradient energy does not depend on shear deformation: this could be considered either a pathological case or an important exceptional case (see Stillwell et al. in Am. Math. Mon. 105:850–858, 1998 and Turro in Angew. Chem., Int. Ed. Engl. 39:2255–2259, 2000). In both cases we believe that such a particular case deserves some attention because of what we can understand by studying it (see Dyson in Science 200:677–678, 1978). This circumstance motivates the present paper, where we address the well-posedness of the planar linearized equilibrium problem for homogenized pantographic lattices. To do so: (i) we introduce a class of subsets of anisotropic Sobolev’s space as the most suitable energy space \(E\) relative to assigned boundary conditions; (ii) we prove that the considered strain energy density is coercive and positive definite in \(E\); (iii) we prove that the set of placements for which the strain energy is vanishing (the so-called floppy modes) must strictly include rigid motions; (iv) we determine the restrictions on displacement boundary conditions which assure existence and uniqueness of linear static problems. The presented results represent one of the first mechanical applications of the concept of Anisotropic Sobolev space, initially introduced only on the basis of purely abstract mathematical considerations.     

Local Fixed Point Indices of Iterations of Planar Maps

Let \({f: U\rightarrow {\mathbb R}^2}\) be a continuous map, where U is an open subset of \({{\mathbb R}^2}\). We consider a fixed point p of f which is neither a sink nor a source and such that {p} is an isolated invariant set. Under these assumption we prove, using Conley index methods and Nielsen theory, that the sequence of fixed point indices of iterations \({\{{\rm ind}(f^n,p)\}_{n=1}^\infty}\) is periodic, bounded from above by 1, and has infinitely many non-positive terms, which is a generalization of Le Calvez and Yoccoz theorem (Annals of Math., 146, 241–293 (1997)) onto the class of non-injective maps. We apply our result to study the dynamics of continuous maps on 2-dimensional sphere.     

Statistical and Numerical Considerations of Backus-Average Product Approximation

In this paper, we examine the applicability of the approximation, \(\overline{f g}\approx \overline{f}\,\overline{g}\), within Backus (J. Geophys. Res. 67(11):4427–4440, 1962) averaging. This approximation is a crucial step in the method proposed by Backus (J. Geophys. Res. 67(11):4427–4440, 1962), which is widely used in studying wave propagation in layered Hookean solids. According to this approximation, the average of the product of a rapidly varying function and a slowly varying function is approximately equal to the product of the averages of those two functions.Considering that the rapidly varying function represents the mechanical properties of layers, we express it as a step function. The slowly varying function is continuous, since it represents the components of the stress or strain tensors. In this paper, beyond the upper bound of the error for that approximation, which is formulated by Bos et al. (J. Elast. 127:179–196, 2017), we provide a statistical analysis of the approximation by allowing the function values to be sampled from general distributions.Even though, according to the upper bound, Backus (J. Geophys. Res. 67(11):4427–4440, 1962) averaging might not appear as a viable approach, we show that—for cases representative of physical scenarios modelled by such an averaging—the approximation is typically quite good. We identify the cases for which there can be a deterioration in its efficacy.In particular, we examine a special case for which the approximation results in spurious values. However, such a case—though physically realizable—is not likely to appear in seismology, where Backus (J. Geophys. Res. 67(11):4427–4440, 1962) averaging is commonly used. Yet, such values might occur in material sciences, in general, for which Backus (J. Geophys. Res. 67(11):4427–4440, 1962) averaging is also considered.     

Second Order Corrections to Mean Field Evolution for Weakly Interacting Bosons in The Case of Three-body Interactions

In this paper, we consider the Hamiltonian evolution of N weakly interacting bosons. Assuming triple collisions, its mean field approximation is given by a quintic Hartree equation. We construct a second order correction to the mean field approximation using a kernel k(t, x, y) and derive an evolution equation for k. We show global existence for the resulting evolution equation for the correction and establish an a priori estimate comparing the approximation to the exact Hamiltonian evolution. Our error estimate is global and uniform in time. Comparing with the work of Rodnianski and Schlein (Commun Math Phys 291:31–61, 2009), and Grillakis, Machedon and Margetis (Commun Math Phys 294:273–301, 2010; Adv Math 288:1788–1815, 2011), where the error estimate grows in time, our approximation tracks the exact dynamics for all time with an error of the order \({O(1/\sqrt{N}).}\)     

Complete Global and Bifurcation Analysis of a Stoichiometric Predator–Prey Model

Loladze et al. (Bull Math Biol 62:1137–1162, 2000) proposed a highly cited stoichiometric predator–prey system, which is nonsmooth, and thus it is extremely difficult to analyze its global dynamics. The main challenge comes from the phase plane fragmentation and parameter space partitioning in order to perform a detailed and complete global stability and bifurcation analysis. Li et al. (J Math Biol 63:901–932, 2011) firstly discussed its global dynamical behavior with Holling type I functional response and found that the system has no limit cycles, and the internal equilibrium is globally asymptotically stable if it exists. Secondly, for the system with Holling type II functional response, Li et al. (2011) fixed all parameters (with realistic values) except K to perform the bifurcation analysis and obtained some interesting phenomena, for instance, the appearance of bistability and many bifurcation types. The aim of this paper is to provide a complete global analysis for the system with Holling type II functional response without fixing any parameter. Our analysis shows that the model has far richer dynamics than those found in the previous paper (Li et al. 2011), for example, four types of bistability appear: besides the bistability between an internal equilibrium and a limit cycle as shown in Li et al. (2011), the other three bistabilities occur between an internal equilibrium and a boundary equilibrium, between two internal equilibria, or between a boundary equilibrium and a limit cycle. In addition, this paper rigorously provides all possible bifurcation passways of this stoichiometric model with Holling type II functional response.     

On Lie symmetries and soliton solutions of $$$$-dimensional Bogoyavlenskii equations

The present article is devoted to find some invariant solutions of the \((2+1)\)-dimensional Bogoyavlenskii equations using similarity transformations method. The system describes \((2+1)\)-dimensional interaction of a Riemann wave propagating along y-axis with a long wave along x-axis. All possible vector fields, commutative relations and symmetry reductions are obtained by using invariance property of Lie group. Meanwhile, the method reduces the number of independent variables by one, which leads to the reduction of Bogoyavlenskii equations into a system of ordinary differential equations. The system so obtained is solved under some parametric restrictions and provides invariant solutions. The derived solutions are much efficient to explain the several physical properties depending upon various existing arbitrary constants and functions. Moreover, some of them are more general than previously established results (Peng and Shen in Pramana 67:449–456, 2006; Malik et al. in Comput Math Appl 64:2850–2859, 2012; Zahran and Khater in Appl Math Model 40:1769–1775, 2016; Zayed and Al-Nowehy in Opt Quant Electron 49(359):1–23, 2017). In order to provide rich physical structures, the solutions are supplemented by numerical simulation, which yield some positons, negatons, kinks, wavefront, multisoliton and asymptotic nature.