In this work, the kinetically reduced local Navier-Stokes equations are applied to the simulation of two- and three-dimensional unsteady viscous incompressible flow problems. The reduced differential transform method is used to find the new approximate analytical solutions of these flow problems. The new technique has been tested by using four selected multidimensional unsteady flow problems: two- and three-dimensional Taylor decaying vortices flow, Kovasznay flow, and three-dimensional Beltrami flow. The convergence analysis was discussed for this approach. The numerical results obtained by this approach are compared with other results that are available in previous works. Our results show that this method is efficient to provide new approximate analytic solutions. Moreover, we found that it has highly precise solutions with good convergence, less time consuming, being easily implemented for high Reynolds numbers, and low Mach numbers.

In this paper, optimal homotopy asymptotic method (OHAM) and its implementation on subinterval, called multistage optimal homotopy asymptotic method (MOHAM), are presented for solving linear and nonlinear systems of Volterra integral equations of the second kind. To illustrate these approaches two examples are presented. The results confirm the efficiency and ability of these methods for such equations. The results will be compared to find out which method is more accurate. Advantages of applying MOHAM are also illustrated.

In this paper, a hybrid iterative algorithm is proposed for finding a common element of the set of common fixed points of finite family of nonexpansive mappings and the set of common solutions of the variational inequality for an inverse strongly monotone mapping on the real Hilbert space. We establish the strong convergence of the proposed method for approximating a common element of the above defined sets under some suitable conditions. The results presented in this paper extend and improve some well-known corresponding results in the earlier and recent literature.

In this paper we present a mathematical model for the transmission dynamics of Contagious Bovine Pleuropneumonia (CBPP) by considering antibiotic treatment and vaccination. The model is comprised of susceptible, vaccinated, exposed, infectious, persistently infected, and recovered compartments. We analyse the model by deriving a formula for the control reproduction number ${\mathcal{R}}_c$ and prove that, for , the disease free equilibrium is globally asymptotically stable; thus CBPP dies out, whereas for ${\mathcal{R}}_c>\mathrm{1}$, the unique endemic equilibrium is globally asymptotically stable and hence the disease persists. Thus, ${\mathcal{R}}_c=\mathrm{1}$ acts as a sharp threshold between the disease dying out or causing an epidemic. As a result, the threshold of antibiotic treatment is ${\alpha }_t^{*}=\mathrm{0.1049}$. Thus, without using vaccination, more than $\mathrm{85.45}\mathrm{\%}$ of the infectious cattle should receive antibiotic treatment or the period of infection should be reduced to less than 8.15 days to control the disease. Similarly, the threshold of vaccination is ${\rho }^{*}=\mathrm{0.0084}$. Therefore, we have to vaccinate at least $\mathrm{80}\mathrm{\%}$ of susceptible cattle in less than 49.5 days, to control the disease. Using both vaccination and antibiotic treatment, the threshold value of vaccination depends on the rate of antibiotic treatment, ${\alpha }_t,$ and is denoted by ${\rho }_{{\alpha }_t}$. Hence, if $\mathrm{50}\mathrm{\%}$ of infectious cattle receive antibiotic treatment, then at least $\mathrm{50}\mathrm{\%}$ of susceptible cattle should get vaccination in less than 73.8 days in order to control the disease.

In this work, we investigate various types of solutions for the generalised resonant dispersive nonlinear Schrödinger equation (GRD-NLSE) with power law nonlinearity. Based on simple mathematical techniques, the complicated form of the GRD-NLSE is reduced to an ordinary differential equation (ODE) which has a variety of solutions. The analytic solution of the resulting ODE gives rise to bright soliton, singular soliton, peaked soliton, compacton solutions, solitary pattern solutions, rational solution, Weierstrass elliptic periodic type solutions, and some other types of solutions. Constraint conditions for the existence of solitons and other solutions are given.

In this paper, dysentery diarrhea deterministic compartmental model is proposed. The local and global stability of the disease-free equilibrium is obtained using the stability theory of differential equations. Numerical simulation of the system shows that the backward bifurcation of the endemic equilibrium exists for $R_{\mathrm{0}}>\mathrm{1}$. The system is formulated as a standard nonlinear least squares problem to estimate the parameters. The estimated reproduction number, based on the dysentery diarrhea disease data for Ethiopia in 2017, is $R_{\mathrm{0}}=\mathrm{1.1208}$. This suggests that elimination of the dysentery disease from Ethiopia is not practical. A graphical method is used to validate the model. Sensitivity analysis is carried out to determine the importance of model parameters in the disease dynamics. It is found out that the reproduction number is the most sensitive to the effective transmission rate of dysentery diarrhea (${\beta }_h$). It is also demonstrated that control of the effective transmission rate is essential to stop the spreading of the disease.

In this study, an optimal control theory was applied to a nonautonomous model for Newcastle disease transmission in the village chicken population. A notable feature of this model is the inclusion of environment contamination and wild birds, which act as reservoirs of the disease virus. Vaccination, culling, and environmental hygiene and sanitation time dependent control strategies were adopted in the proposed model. This study proved the existence of an optimal control solution, and the necessary conditions for optimality were determined using Pontryagin’s Maximum Principle. The numerical simulations of the optimal control problem were performed using the forward–backward sweep method. The results showed that the use of only the environmental hygiene and sanitation control strategy has no significant effect on the transmission dynamics of the Newcastle disease. Additionally, the combination of vaccination and environmental hygiene and sanitation strategies reduces more number of infected chickens and the concentration of the Newcastle disease virus in the environment than any other combination of control strategies. Furthermore, a cost-effective analysis was performed using the incremental cost-effectiveness ratio method, and the results showed that the use of vaccination alone as the control measure is less costly compared to other control strategies. Hence, the most effective way to minimize the transmission rate of the Newcastle disease and the operational costs is concluded to be the timely vaccination of the entire population of the village chicken, improvement in the sanitation of facilities, and the maintenance of a hygienically clean environment.

A rainbow $t$-coloring of a $t$-connected graph $G$ is an edge coloring such that for any two distinct vertices $u$ and $v$ of $G$ there are at least $t$ internally vertex-disjoint rainbow $(u,v)$-paths. In this work, we apply a Rank Genetic Algorithm to search for rainbow $t$-colorings of the family of Moore cages with girth six $(t;\mathrm{6})$-cages. We found that an upper bound in the number of colors needed to produce a rainbow 4-coloring of a $(\mathrm{4};\mathrm{6})$-cage is 7, improving the one currently known, which is 13. The computation of the minimum number of colors of a rainbow coloring is known to be NP-Hard and the Rank Genetic Algorithm showed good behavior finding rainbow $t$-colorings with a small number of colors.

After the dawn of the August 2007 financial crisis, banks became more aware of financial risk leading to the appearance of nonnegligible spreads between LIBOR and OIS rates and also between LIBOR of different tenors. This consequently led to the birth of multicurve models. This study establishes a new model; the multicurve cross-currency LIBOR market model (MCCCLMM). The model extends the initial LIBOR Market Model (LMM) from the single-curve cross-currency economy into the multicurve cross-currency economy. The model incorporates both the risk-free OIS rates and the risky forward LIBOR rates of two different currencies. The established model is suitable for pricing different quanto interest rate derivatives. A brief illustration is given on the application of the MCCCLMM on pricing quanto caplets and quanto floorlets using a Black-like formula derived from the MCCCLMM.

In this paper, the magneto hydrodynamic (MHD) squeezing flow of a non-Newtonian, namely, Casson, fluid between parallel plates is studied. The suitable one of similarity transformation conversion laws is proposed to obtain the governing MHD flow nonlinear ordinary differential equation. The resulting equation has been solved by a novel algorithm. Comparisons between the results of the novel algorithm technique and other analytical techniques and one numerical Range-Kutta fourth-order algorithm are provided. The results are found to be in excellent agreement. Also, a novel convergence proof of the proposed algorithm based on properties of convergent series is introduced. Flow behavior under the changing involved physical parameters such as squeeze number, Casson fluid parameter, and magnetic number is discussed and explained in detail with help of tables and graphs.

In this study, we tend to propose a replacement hybrid algorithmic rule which mixes the search directions like Steepest Descent (SD) and Quasi-Newton (QN). First, we tend to develop a replacement search direction for combined conjugate gradient (CG) and QN strategies. Second, we tend to depict a replacement positive CG methodology that possesses the adequate descent property with sturdy Wolfe line search. We tend to conjointly prove a replacement theorem to make sure global convergence property is underneath some given conditions. Our numerical results show that the new algorithmic rule is powerful as compared to different standard high scale CG strategies.

We investigate the filtering problem where the borrower’s time varying credit quality process is estimated using continuous time observation process and her (in this paper we refer to the borrower as female and the lender as male) ego-network data. The hidden credit quality is modeled as a hidden Gaussian mean-reverting process whilst the social network is modeled as a continuous time latent space network model. At discrete times, the network data provides unbiased estimates of the current credit state of the borrower and her ego-network. Combining the continuous time observed behavioral data and network information, we provide filter equations for the hidden credit quality and show how the network information reduces information asymmetry between the borrower and the lender. Further, we consider the case when the network information arrival times are random and solve stochastic optimal control problem for a lender having linear quadratic utility function.

We deem the time-fractional Benjamin-Ono (BO) equation out of the Riemann–Liouville (RL) derivative by applying the Lie symmetry analysis (LSA). By first using prolongation theorem to investigate its similarity vectors and then using these generators to transform the time-fractional BO equation to a nonlinear ordinary differential equation (NLODE) of fractional order, we complete the solutions by utilizing the power series method (PSM).

In this paper we elaborate an algorithm to estimate p-order Random Coefficient Autoregressive Model (RCA(p)) parameters. This algorithm combines quasi-maximum likelihood method, the Kalman filter, and the simulated annealing method. In the aim to generalize the results found for RCA(1), we have integrated a subalgorithm which calculate the theoretical autocorrelation. Simulation results demonstrate that the algorithm is viable and promising.

An interest is often present in knowing evolving variables that are not directly observable; this is the case in aerospace, engineering control, medical imaging, or data assimilation. What is at hand, though, are time-varying measured data, a model connecting them to variables of interest, and a model of how to evolve the variables over time. However, both models are only approximation and the observed data are tainted with noise. This is an ill-posed inverse problem. Methods, such as Kalman filter (KF), have been devised to extract the time-varying quantities of interest. These methods applied to this inverse problem, nonetheless, are slow, computation wise, since they require large matrices multiplications and even matrix inversion. Furthermore, these methods are not usually suitable to impose some constraints. This article introduces a new iterative filtering algorithm based on alternating projections. Experiments were run with simulated moving projectiles and were compared with results using KF. The new optimization algorithm proves to be slightly more accurate than KF, but, more to the point, it is much faster in terms of CPU time.

This paper develops a mapping approach to explore the relationship between different layers of a multilayer air transport network composed of airway, route, and flight network. A two-step methodology is adopted to investigate the hierarchical structure and mapping relationship of the integrated network. First, the relationship between airway and route network is characterized by a multisource multisink shortest path method based on a generalized incidence matrix. Second, the relationship between route and flight network is formulated by a two-dimension array. A case study of an en route airspace in Lanzhou air traffic control area in China verifies the feasibility of the proposed two-step methodology.

In this paper, a new class of two-derivative two-step Runge-Kutta (TDTSRK) methods for the numerical solution of non-stiff initial value problems (IVPs) in ordinary differential equation (ODEs) is considered. The TDTSRK methods are a special case of multi-derivative Runge-Kutta methods proposed by Kastlunger and Wanner (1972). The methods considered herein incorporate only the first and second derivatives terms of ODEs. These methods possess large interval of stability when compared with other existing methods in the literature. The experiments have been performed on standard problems, and comparisons were made with some standard explicit Runge-Kutta methods in the literature.

In the study presented in this paper, we built a nonlinear binary integer programming model of a flexible scheduling problem for the Department of Zhejiang Provincial Local Tax Services. One difference between our model and typical ones is that whereas in the latter the number of open windows within each working day is fixed, in our model it is not. We used a variety of integer programming software in an attempt to solve our scheduling model; however, unfortunately we could not find an optimal solution. Thus, we tested all the combinations of different numbers of employees to construct the optimal solution. When we tested our model in the tax office of Lishui City, China, the average waiting time of taxpayers was less than 15 min and the employees working hours were clearly reduced. Thus, a noteworthy improvement in the quality of the service is achieved by the model.

A mathematical model is proposed and analysed to study the dynamics of two-prey one predator system of fishery model with Holling type II function response. The effect of harvesting was incorporated to both populations and thoroughly analysed. We study the ecological dynamics of the Nile perch, cichlid, and tilapia fishes as prey-predator system of lake Victoria fishery in Tanzania. In both cases, by nondimensionalization of the system, the equilibrium points are computed and conditions for local and global stability of the system are obtained. Condition for local stability was obtained by eigenvalue approach and Routh-Hurwitz Criterion. Moreover, the global stability of the coexistence equilibrium point is proved by defining appropriate Lyapunov function. Bioeconomic equilibrium is analysed and numerical simulations are also carried out to verify the analytical results. The numerical results indicate that the three species would coexist if cichlid and tilapia fishes will not be overharvested as these populations contribute to the growth rates of Nile perch population. The fishery control management should be exercised to avoid overharvesting of cichlid and tilapia fishes.

The nonmonotone alternating direction algorithm (NADA) was recently proposed for effectively solving a class of equality-constrained nonsmooth optimization problems and applied to the total variation minimization in image reconstruction, but the reconstructed images suffer from the artifacts. Though by the $l_{\mathrm{0}}$-norm regularization the edge can be effectively retained, the problem is NP hard. The smoothed $l_{\mathrm{0}}$-norm approximates the $l_{\mathrm{0}}$-norm as a limit of smooth convex functions and provides a smooth measure of sparsity in applications. The smoothed $l_{\mathrm{0}}$-norm regularization has been an attractive research topic in sparse image and signal recovery. In this paper, we present a combined smoothed $l_{\mathrm{0}}$-norm and $l_{\mathrm{1}}$-norm regularization algorithm using the NADA for image reconstruction in computed tomography. We resolve the computation challenge resulting from the smoothed $l_{\mathrm{0}}$-norm minimization. The numerical experiments demonstrate that the proposed algorithm improves the quality of the reconstructed images with the same cost of CPU time and reduces the computation time significantly while maintaining the same image quality compared with the $l_{\mathrm{1}}$-norm regularization in absence of the smoothed $l_{\mathrm{0}}$-norm.

Based on the Hirota bilinear form of the generalized (3+1)-dimensional variable-coefficient B-type Kadomtsev-Petviashvili equation, the lump and lump-type solutions are generated through symbolic computation, whose analyticity can be easily achieved by taking special choices of the involved parameters. The property of solutions is investigated and exhibited vividly by three-dimensional plots and contour plots.

Catastrophe is a loss that has a low probability of occurring but can lead to high-cost claims. This paper uses the data of fire accidents from a reinsurance company in Thailand for an experiment. Our study is in two parts. First, we approximate the parameters of a Weibull distribution. We compare the parameter estimation using a direct search method with other frequently used methods, such as the least squares method, the maximum likelihood estimation, and the method of moments. The results show that the direct search method approximates the parameters more precisely than other frequently used methods (to four-digit accuracy). Second, we approximate the minimum initial capital (MIC) a reinsurance company has to hold under a given ruin probability (insolvency probability) by using parameters from the first part. Finally, we show MIC with varying the premium rate.

Recent studies reveal that Allee effect may play important roles in the growth of tumor. We present one of the first mathematical models of avascular tumor that incorporates the weak Allee effect. The model considers the densities of tumor cells in three stages: proliferating cells, quiescent cells, and necrotic cells. We investigate how Allee effect impacts the growth of the avascular tumor. We also investigate the effect of apoptosis of proliferating cells and necrosis of quiescent cells. The system is numerically solved in 2D using different sets of parameters. We show that Allee effect and apoptosis play important roles in the growth of tumor and the formation of necrotic core.

The nonlinear convective flow of Eyring-Powell nanofluid using Catteneo-Christov model with heat generation or absorption term and chemical reaction rate over nonlinear stretching surface is analyzed. The simultaneous nonlinear partial differential equations governing the boundary layer flow are transformed to the corresponding nonlinear ordinary differential equations using similarity solution and then solved using Galerkin finite element method (GFEM). The impacts of pertinent governing parameters like Brownian diffusion, thermophoresis, mixed convection, heat generation or absorption, chemical reaction rate, Deborah numbers, Prandtl number, magnetic field parameter, Lewis number, nonlinear stretching sheet, and Eyring-Powell fluid parameters on velocity field, temperature, and nanoparticle concentration are given in both figures and tabular form. The result shows that the rise in chemical reaction rate will improve mass transfer rate and reduce heat transfer rate and local buoyancy parameter has quit opposite effect. The attributes of local skin friction coefficient, Nusselt number, and Sheer wood number are investigated and validated with existing literatures.

We prove the polynomial analogues of some Liouville identities from elementary number theory. Consequently several sums defined over the finite fields ${\mathbb{F}}_q[t]$ are evaluated by combining the results obtained and some of the results from sums of reciprocals of polynomials over ${\mathbb{F}}_q[t]$.

In this work, we study the statistical properties of molecular systems of admixtures, which are placed in nanoscale volumes like cube or parallelepiped with impermeable surfaces on their walls. We simulate interactions of free motion of molecules of admixtures with molecules on the surfaces of walls: we modeled them in the form of atomic structure grid using SiO₂ as a material. This type of substance allows us to take into account the manifestation of one of the important quantum effects: Casimir force. We used its general expression in terms of interaction energy, with the assumption of dependence on distance and projected area between atoms of the wall and atoms of the admixtures. To model surface roughness we used uniform random distribution function for surface heights (z coordinates). The results of computational experiments can be used to estimate the distribution of chemical bond lengths, valence, and dihedral angles lengths deviations in polyatomic molecules. Our model can be useful to determine the stable configuration properties of the system, to solve practical problems in the conditions of physical limitations of nanoscale devices, filtration of admixtures in highly dispersed systems, in the development of mechanisms in structures that have parallel plates or membranes, including porous structures. Also we compute radial distribution function for multicomponent admixtures including atoms of inert gases, water vapor, and ethanol. And we took into account the influence of Casimir force on admixtures for cube and parallelepiped type volumes. In results, we showed distributions of deviations of lengths from their equilibrium values.

In the article titled “Study of Two-Sided Similarity Methods Using a Radiation “Switch on” Imploding Shock in a Magnetic Field” [1], there were errors in the equation formatting that should be corrected as follows.

The problem of building a smooth surface containing given points or curves is actual due to development of industry and computer technology. Previously used for those purposes, shells of zero Gaussian curvature and minimal surfaces based on isotropic analytic curves are restricted in their consumer properties. To expand the possibilities regarding the shaping of surfaces we propose the method of constructing surfaces based on isotropic fractional-rational curves. The surfaces are built using flat isothermal and orthogonal grids and on the basis of the Weierstrass method. In the latter case, the surfaces are minimal. Examples of surfaces that were built according to the proposed method are given.

The integrals evaluated are the products of multivariable Aleph-functions with algebraic functions, Jacobi polynomials, Legendre functions, Bessel-Maitland functions, and general class of polynomials. The main results of our paper are quite general in nature and competent at yielding a very large number of integrals involving polynomials and various special functions occurring in the problem of mathematical analysis and mathematical physics.

Given a tree network $T$ with $n$ vertices where each edge has an independent operational probability, we are interested in finding the optimal location of a reliable service provider facility in a shape of subtree with exactly $k$ leaves and with a diameter of at most $l$ which maximizes the expected number of nodes that are reachable from the selected subtree by operational paths. Demand requests for service originate at perfectly reliable nodes. So, the major concern of this paper is to find a location of a reliable tree-shaped facility on the network in order to provide a maximum access to network services by ensuring the highest level of network connectivity between the demand nodes and the facility. An efficient algorithm for finding a reliable $(k,l)$ – tree core of $T$ is developed. The time complexity of the proposed algorithm is $O\left(lkn\right).$ Examples are provided to illustrate the performance of the proposed algorithm.

Two numerical techniques, namely, Haar Wavelet and the product integration methods, have been employed to give an approximate solution of the fractional Volterra integral equation of the second kind. To test the applicability and efficiency of the numerical method, two illustrative examples with known exact solution are presented. Numerical results show clearly that the accuracy of these methods are in a good agreement with the exact solution. A comparison between these methods shows that the product integration method provides more accurate results than its counterpart.

In this study, we develop and analyze a deterministic mathematical model for tuberculosis (TB) transmission dynamics. The model includes vaccination for newborns and treatment for both high-risk latent and active TB patients. The stability of disease-free equilibrium point is discussed in detail. In the numerical simulation, the model parameters are estimated using reported TB incidence data in Ethiopia from the years 2003 to 2017, and ${\mathcal{R}}_0$ is calculated as ${\mathcal{R}}_0\approx 2.13$. Finally, the sensitivity indices of ${\mathcal{R}}_0$ with respect to the model parameters are performed, and their corresponding graphical results are presented. Our results quantify the positive influence of vaccination and the treatment for high-risk latent and active TB patients on the control of tuberculosis.

In this work, we tend to deal within the field of the constrained optimization methods of three-term Conjugate Gradient (CG) technique which is primarily based on Dai–Liao (DL) formula. The new proposed technique satisfies the conjugacy property and the descent conditions of Karush–Kuhn–Tucker (K.K.T.). Our planned constrained technique uses the robust Wolfe line search condition with some assumptions. We tend to prove the global convergence property of the new planned technique. Numeral comparisons for (30-thirty) constrained optimization issues make sure the effectiveness of the new planned formula.

In this paper, a system of deterministic model is presented for the dynamical analysis of the interactional consequence of criminals and criminality on victimisation under two distinguishable forms of rehabilitation—the behavioural reformation of criminals and the emotional psychotherapy of victims. A threshold value, ${\mathcal{R}}_0=\mathrm{m}\mathrm{a}\mathrm{x}\left\{{\mathcal{R}}_K,{\mathcal{R}}_V\right\}$, responsible for the persistence of crime/criminality and victimisation, is obtained and, using it, stability analyses on the model performed. The impact of an effective implementation of the two forms of rehabilitation was found to be substantial on crime and criminality, while an ineffective implementation of same was observed to have a detrimental consequence. The prevention of repeat victimisation was seen to present a more viable option for containing crime than the noncriminalisation of victims. Further, the removal of criminals, either through quitting or death, among others, was also found to have a huge positive impact. Numerical simulations were performed for a variety of mixing criminal scenarios to verify the analytical results obtained.

In this paper, a state-constrained optimal control problem governed by p-Laplacian elliptic equations is studied. The feasible control set or the cost functional may be nonconvex, and the purpose is to obtain the convergence of a solution of the discretized control problem to an optimal control of the relaxed continuous problem.