Search results

The numerical methods are of great importance for approximating the solutions of a system of nonlinear singular ordinary differential equations. In this paper, the authors present the biorthogonal flatlet multiwavelet collocation method (BFMCM) as a numerical scheme for a class of system of Lane–Emden equations with initial or boundary or four-point boundary conditions.

The approach is involved in combining the biorthogonal flatlet multiwavelet (BFM) with the collocation method. The authors investigate the properties and procedure of the BFMCM for first time on this class of equations. By using the BFM and the collocation points, the method is constructed and it transforms the nonlinear differential equations problem into a system of nonlinear algebraic equations. The unknown coefficients of the assuming solution are determined by solving the obtained system. Additionally, convergence analysis and numerical stability of the suggested method are provided.

According to the attained results, the proposed BFMCM has more accurate results in comparison with results of other methods. The maximum absolute errors are calculated by using the BFMCM for comparison purposes provided.

The key desirable properties of BFMCM are its efficiency, simple applicability and minimizes errors. Therefore, the proposed method can be used to solve nonlinear problems or problems with singular points.

The purpose of this study is to solve two unique but difficult partial differential equations: the foam drainage equation and the nonlinear time-fractional fisher’s equation. Through our methods, we aim to provide accurate solutions and gain a deeper understanding of the intricate behaviors exhibited by these systems.

In this study, we use a dual technique that combines the Aboodh residual power series method and the Aboodh transform iteration method, both of which are combined with the Caputo operator.

We develop exact and efficient solutions by merging these unique methodologies. Our results, presented through illustrative figures and data, demonstrate the efficacy and versatility of the Aboodh methods in tackling such complex mathematical models.

Owing to their fractional derivatives and nonlinear behavior, these equations are crucial in modeling complex processes and confront analytical complications in various scientific and engineering contexts.

The purpose of this paper is to analyze numerically the blowup in finite time of the solutions to a one-dimensional, bidirectional, nonlinear wave model equation for the propagation of small-amplitude waves in shallow water, as a function of the relaxation time, linear and nonlinear drift, power of the nonlinear advection flux, viscosity coefficient, viscous attenuation, and amplitude, smoothness and width of three types of initial conditions.

An implicit, first-order accurate in time, finite difference method valid for semipositive relaxation times has been used to solve the equation in a truncated domain for three different initial conditions, a first-order time derivative initially equal to zero and several constant wave speeds.

The numerical experiments show a very rapid transient from the initial conditions to the formation of a leading propagating wave, whose duration depends strongly on the shape, amplitude and width of the initial data as well as on the coefficients of the bidirectional equation. The blowup times for the triangular conditions have been found to be larger than those for the Gaussian ones, and the latter are larger than those for rectangular conditions, thus indicating that the blowup time decreases as the smoothness of the initial conditions decreases. The blowup time has also been found to decrease as the relaxation time, degree of nonlinearity, linear drift coefficient and amplitude of the initial conditions are increased, and as the width of the initial condition is decreased, but it increases as the viscosity coefficient is increased. No blowup has been observed for relaxation times smaller than one-hundredth, viscosity coefficients larger than ten-thousandths, quadratic and cubic nonlinearities, and initial Gaussian, triangular and rectangular conditions of unity amplitude.

The blowup of a one-dimensional, bidirectional equation that is a model for the propagation of waves in shallow water, longitudinal displacement in homogeneous viscoelastic bars, nerve conduction, nonlinear acoustics and heat transfer in very small devices and/or at very high transfer rates has been determined numerically as a function of the linear and nonlinear drift coefficients, power of the nonlinear drift, viscosity coefficient, viscous attenuation, and amplitude, smoothness and width of the initial conditions for nonzero relaxation times.

The study aims to optimize truck routes by minimizing social and economic costs. It introduces a strategy involving diverse drones and their potential for reusing at DNs based on flight range. In HTDRP-DC, trucks can select and transport various drones to LDs to reduce deprivation time. This study estimates the nonlinear deprivation cost function using a linear two-piece-wise function, leading to MILP formulations. A heuristic-based Benders Decomposition approach is implemented to address medium and large instances. Valid inequalities and a heuristic method enhance convergence boundaries, ensuring an efficient solution methodology.

Research has yet to address critical factors in disaster logistics: minimizing the social and economic costs simultaneously and using drones in relief distribution; deprivation as a social cost measures the human suffering from a shortage of relief supplies. The proposed hybrid truck-drone routing problem minimizing deprivation cost (HTDRP-DC) involves distributing relief supplies to dispersed demand nodes with undamaged (LDs) or damaged (DNs) access roads, utilizing multiple trucks and diverse drones. A Benders Decomposition approach is enhanced by accelerating techniques.

Incorporating deprivation and economic costs results in selecting optimal routes, effectively reducing the time required to assist affected areas. Additionally, employing various drone types and their reuse in damaged nodes reduces deprivation time and associated deprivation costs. The study employs valid inequalities and the heuristic method to solve the master problem, substantially reducing computational time and iterations compared to GAMS and classical Benders Decomposition Algorithm. The proposed heuristic-based Benders Decomposition approach is applied to a disaster in Tehran, demonstrating efficient solutions for the HTDRP-DC regarding computational time and convergence rate.

Current research introduces an HTDRP-DC problem that addresses minimizing deprivation costs considering the vehicle’s arrival time as the deprivation time, offering a unique solution to optimize route selection in relief distribution. Furthermore, integrating heuristic methods and valid inequalities into the Benders Decomposition approach enhances its effectiveness in solving complex routing challenges in disaster scenarios.

As to different angles of attack and nonlinear problems caused by high temperatures in coexisting hypersonic aircraft, people mainly rely on fluid software for research but lack analysis of flow mechanisms. Owing to computational difficulties, few people use numerical algorithms to combine them for discussion. Hence, this study aims to make a deep inquiry into the laminar flow and heat transfer of compressible Newtonian fluid in hypersonic aircraft with small attack angles.

In this paper, on the basis of mass, momentum and energy conservation laws, the governing equations of the hypersonic boundary layer are established. Viscosity, specific heat capacity and thermal conductivity are considered nonlinear functions concerning temperature. In virtue of the MacCormack finite difference method, the stationary numerical solutions are solved directly, and the validity of the algorithm is verified.

The results demonstrate that at Mach number 5, compared to the 0° attack angle, the maximum temperature near-wall at the 3° attack angle increases by about 25%. An enjoyable phenomenon is discovered, where the position corresponding to the maximum wall shear force shifts back as the attack angle and Mach number increase. The relationship between the near-wall maximum temperature versus attack angle and Mach number is fitted through numerical calculation results.

Empirical formulas can be used to estimate heat transfer characteristics at small attack angles, which will guide the design of aircraft thermal protection systems.

Ferrofluids are aqueous or non-aqueous solutions with colloidal particles of iron oxide nanoparticles with high magnetic characteristics. Their magnetic characteristics enable them to be controlled and manipulated when ferrofluids are exposed to magnetic fields. This study aims to inspect the features of unsteady stagnation point flow (SPF) and heat flux from the surface by incorporating ferromagnetic particles through a special kind of second-grade fluid (SGF) across a movable sheet with a nonlinear heat source/sink and magnetic field effect. The mass suction/injection and stretching/shrinking boundary conditions are also inspected to calculate the fine points of the features of multiple solutions.

The leading equations that govern the ferrofluid flow are reduced to a group of ordinary differential equations by applying similarity variables. The converted equations are numerically solved through the bvp4c solver. Afterward, study and discussion are carried out to examine the different physical parameters of the characteristics of nanofluid flow and thermal properties.

Multiple solutions are revealed to happen for situations of unsteadiness, shrinking as well as stretching sheets. Greater suction slows the separation of the boundary layers and causes the critical values to expand. The region where the multiple solutions appear is observed to expand with increasing values of the magnetic, non-Newtonian and suction parameters. Moreover, the fluid velocity significantly uplifts while the temperature declines due to the suction parameter.

The novelty of the work is to deliberate the impact of mass suction/injection on the unsteady SPF through the special second-grade ferrofluids across a movable sheet with an erratic heat source/sink. The confirmed results provide a very good consistency with the accepted papers. Previous studies have not yet fully explored the entire analysis of the proposed model.

This study aims to derive a novel spatial numerical method based on multidimensional local Taylor series representations for solving high-order advection-diffusion (AD) equations.

The parabolic AD equations are reduced to the nonhomogeneous elliptic system of partial differential equations by utilizing the Chebyshev spectral collocation method (ChSCM) in the temporal variable. The implicit-explicit local differential transform method (IELDTM) is constructed over two- and three-dimensional meshes using continuity equations of the neighbor representations with either explicit or implicit forms in related directions. The IELDTM yields an overdetermined or underdetermined system of algebraic equations solved in the least square sense.

The IELDTM has proven to have excellent convergence properties by experimentally illustrating both h-refinement and p-refinement outcomes. A distinctive feature of the IELDTM over the existing numerical techniques is optimizing the local spatial degrees of freedom. It has been proven that the IELDTM provides more accurate results with far fewer degrees of freedom than the finite difference, finite element and spectral methods.

This study shows the derivation, applicability and performance of the IELDTM for solving 2D and 3D advection-diffusion equations. It has been demonstrated that the IELDTM can be a competitive numerical method for addressing high-space dimensional-parabolic partial differential equations (PDEs) arising in various fields of science and engineering. The novel ChSCM-IELDTM hybridization has been proven to have distinct advantages, such as continuous utilization of time integration and optimized formulation of spatial approximations. Furthermore, the novel ChSCM-IELDTM hybridization can be adapted to address various other types of PDEs by modifying the theoretical derivation accordingly.

This study aims to analyze the two-dimensional ferrofluid flow in porous media. The effects of changes in parameters such as permeability parameter, buoyancy parameter, Reynolds and Prandtl numbers, radiation parameter, velocity slip parameter, energy dissipation parameter and viscosity parameter on the velocity and temperature profile are displayed numerically and graphically.

By using simplification, nonlinear differential equations are converted into ordinary nonlinear equations. Modeling is done in the Cartesian coordinate system. The finite element method (FEM) and the Akbari-Ganji method (AGM) are used to solve the present problem. The finite element model determines each parameter’s effect on the fluid’s velocity and temperature.

The results show that if the viscosity parameter increases, the temperature of the fluid increases, but the velocity of the fluid decreases. As can be seen in the figures, by increasing the permeability parameter, a reduction in velocity and an enhancement in fluid temperature are observed. When the Reynolds number increases, an increase in fluid velocity and temperature is observed. If the speed slip parameter increases, the speed decreases, and as the energy dissipation parameter increases, the temperature also increases.

When considering factors like thermal conductivity and variable viscosity in this context, they can significantly impact velocity slippage conditions. The primary objective of the present study is to assess the influence of thermal conductivity parameters and variable viscosity within a porous medium on ferrofluid behavior. This particular flow configuration is chosen due to the essential role of ferrofluids and their extensive use in engineering, industry and medicine.

The need for precise synthesis of customized designs has resulted in the development of advanced coating processes for modern nanomaterials. Achieving accuracy in these processes requires a deep understanding of thermophysical behavior, rheology and complex chemical reactions. The manufacturing flow processes for these coatings are intricate and involve heat and mass transfer phenomena. Magnetic nanoparticles are being used to create intelligent coatings that can be externally manipulated, making them highly desirable. In this study, a Keller box calculation is used to investigate the flow of a coating nanofluid containing a viscoelastic polymer over a circular cylinder.

The rheology of the coating polymer nanofluid is described using the viscoelastic model, while the effects of nanoscale are accounted for by using Buongiorno’s two-component model. The nonlinear PDEs are transformed into dimensionless PDEs via a nonsimilar transformation. The dimensionless PDEs are then solved using the Keller box method.

The transport phenomena are analyzed through a comprehensive parametric study that investigates the effects of various emerging parameters, including thermal radiation, Biot number, Eckert number, Brownian motion, magnetic field and thermophoresis. The results of the numerical analysis, such as the physical variables and flow field, are presented graphically. The momentum boundary layer thickness of the viscoelastic polymer nanofluid decreases as fluid parameter increases. An increase in mixed convection parameter leads to a rise in the Nusselt number. The enhancement of the Brinkman number and Biot number results in an increase in the total entropy generation of the viscoelastic polymer nanofluid.

Intelligent materials rely heavily on the critical characteristic of viscoelasticity, which displays both viscous and elastic effects. Viscoelastic models provide a comprehensive framework for capturing a range of polymeric characteristics, such as stress relaxation, retardation, stretching and molecular reorientation. Consequently, they are a valuable tool in smart coating technologies, as well as in various applications like supercapacitor electrodes, solar collector receivers and power generation. This study has practical applications in the field of coating engineering components that use smart magnetic nanofluids. The results of this research can be used to analyze the dimensions of velocity profiles, heat and mass transfer, which are important factors in coating engineering. The study is a valuable contribution to the literature because it takes into account Joule heating, nonlinear convection and viscous dissipation effects, which have a significant impact on the thermofluid transport characteristics of the coating.

The momentum boundary layer thickness of the viscoelastic polymer nanofluid decreases as the fluid parameter increases. An increase in the mixed convection parameter leads to a rise in the Nusselt number. The enhancement of the Brinkman number and Biot number results in an increase in the total entropy generation of the viscoelastic polymer nanofluid. Increasing the strength of the magnetic field promotes an increase in the density of the streamlines. An increase in the mixed convection parameter results in a decrease in the isotherms and isoconcentration.

This work compares the performance of the three boundary element techniques for solving Helmholtz problems: dual reciprocity, multiple reciprocity and direct interpolation. All techniques transform domain integrals into boundary integrals, despite using different principles to reach this purpose.

Comparisons here performed include the solution of eigenvalue and response by frequency scanning, analyzing many features that are not comprehensively discussed in the literature, as follows: the type of boundary conditions, suitable number of degrees of freedom, modal content, number of primitives in the multiple reciprocity method (MRM) and the requirement of internal interpolation points in techniques that use radial basis functions as dual reciprocity and direct interpolation.

Among the other aspects, this work can conclude that the solution of the eigenvalue and response problems confirmed the reasonable accuracy of the dual reciprocity boundary element method (DRBEM) only for the calculation of the first natural frequencies. Concerning the direct interpolation boundary element method (DIBEM), its interpolation characteristic allows more accessibility for solving more elaborate problems. Despite requiring a greater number of interpolating internal points, the DIBEM has presented higher-quality results for the eigenvalue and response problems. The MRM results were satisfactory in terms of accuracy just for the low range of frequencies; however, the neglected higher-order primitives impact the accuracy of the dynamic response as a whole.

There are safe alternatives for solving engineering stationary dynamic problems using the boundary element method (BEM), but there are no suitable comparisons between these different techniques. This paper presents the particularities and detailed comparisons approaching the accuracy of the three important BEM techniques, aiming at response and frequency evaluation, which are not found in the specialized literature.