A numerical method for multipoint boundary value problems with application to a restricted three body problem

A type of parallel shooting method is proposed for the solution of nonlinear multipoint boundary value problems. It extends the usual quasilinearization method and a previous shooting method developed for such problems, and reduces to usual multiple shooting techniques for two point boundary value problems. The effectiveness of the method for stiff problems is illustrated by an application to the problem of finding periodic solutions of a restricted three body problem with given Jacobian constant and unknown period.     

Analysis of pressure and heat distribution in a dilating or contracting filter chamber with two outlets using multivariate spectral quasilinearization method

The development of efficient filters is an essential part of industrial machinery design, specifically to increase the lifespan of a machine. In the filter chamber design considered in this study, the magnetic material is placed along the horizontal surface of the filter chamber. The inside of the filter chamber is layered with a porous material to restrict the outflow of unwanted particles. This study aims to investigate the flow, pressure, and heat distribution in a dilating or contracting filter chamber with two outlets driven by injection through a permeable surface. The proposed model of the fluid dynamics within the filter chamber follows the conservation equations in the form of partial differential equations. The model equations are further reduced to a steady case through Lie's symmetry group of transformation. They are then solved using a multivariate spectral-based quasilinearization method on the Chebyshev–Gauss–Lobatto nodes. Insights and analyses of the thermophysical parameters that drive optimal outflow during the filtration process are provided through the graphs of the numerical solutions of the differential equations. We find, among other results, that expansion of the filter chamber leads to an overall decrease in internal pressure and an increase in heat distribution inside the filter chamber. The results also show that shrinking the filter chamber increases the internal momentum inside the filter, which leads to more outflow of filtrates.     

An application of paired quasilinearization on double‐diffusive convection flow over a cone embedded in a porous medium in the presence of nanoparticles

Mixed convection flow of a nanofluid near a vertical cone embedded in a a porous medium with Soret and Dufour effects is exercised. The bearing of a porous medium is recounted by the Darcy model. The partial differential equations, modeling the concerned problem, is nondimensionalised by implementing compatible transformations, which results in a similar form. A new paired spectral quasilinearization method is adopted to get the accurate numerical solution. Convergence and accuracy of the solution is elaborated by analyzing the norm of residual and solution errors. Alteration of velocity, temperature, nanoparticle and solute concentration profiles due to flow controlling parameters, namely, Brownian motion, thermophoresis, Soret, Dufour, Lewis number, and buoyancy ratio is outlined by reproducing the obtained numerical solution in graphs and tables. Analysis reveals that the flow profiles are greatly influenced by the physical parameters under investigation.     

MHD mixed convective radiative flow of Eyring‐Powell fluid over an oscillatory stretching sheet using bivariate spectral method on overlapping grids

The bivariate spectral quasilinearization method (BSQLM) on overlapping grids is presented and applied in the analysis of unsteady magnetohydrodynamic mixed convection flow of Eyring‐Powell fluid over an oscillatory stretching sheet embedded in a non‐Darcy porous medium with nonlinear radiative heat flux and variable thermophysical properties. The fluid properties, namely the fluid viscosity, thermal conductivity, and mass diffusivity, are assumed to vary with temperature. It is assumed that the first‐order chemical reaction with heat generation/absorption takes place in the flow. The flow domain is subject to uniform transverse magnetic field perpendicular to the stretching surface. The transformed flow equations are solved numerically using BSQLM on overlapping grids. The convergence properties and accuracy of the method are assessed. The proposed method is computationally efficient, and it gives stable and highly accurate results after few iterations and using few grid points in each subinterval. The improved accuracy rests upon the use of the overlapping grid, which produces sparse coefficient matrices that are easy to invert and have small condition numbers. The effects of physical parameters on the flow fields, local skin friction, the Nusselt number, and the Sherwood number are exhibited through graphs and tables. Amongst other findings, we found that the amplitude of the fluid flow along with flow characteristics may efficiently improve through the utilization of variable fluid viscosity. Heat and mass transportation processes enhance with the inclusion of nonlinear radiative heat flux, temperature‐dependent thermal conductivity, and mass diffusion coefficient, whereas they diminish with the increase in the local inertia coefficient. The current flow analysis can be useful in various engineering applications including paper production, polymer solution, glass blowing, extrusion of thermal system manufacturing process, and heat transportation enhancement.     

Sensitivity analysis of shock wave Burgers’ equation via a novel algorithm based on scale-3 Haar wavelets

In this paper, a novel technique is being formulated for the numerical solutions of Shock wave Burgers' equations for planar and non-planar geometry. It is well known that Burgers' equation is sensitive to the perturbations in the diffusion term. Thus we use robustness of wavelets generated by dilation and translation of Haar wavelets on third scale to capture the sensitivity information. The present approach is an improved form of the scale-2 Haar wavelet method. The scheme is based on the forward finite difference scheme for time integration, scale-3 Haar wavelets for space integration and the nonlinearity has been tackled via quasilinearzation technique. Through scale-3 Haar wavelet analysis once the wavelet coefficient is calculated then we can compute the solutions at near the perturbation point. The computation cost of the present scheme is negligible. The proposed method is tested on six test problems to check its computational efficiency where the convergence analysis of scale-3 Haar wavelet method is the proof of our computational arguments.     

时标上非线性三点边值问题的拟线性方法

利用拟线性方法,研究了时标上非线性动力方程两项和的三点边值问题,得到了两个一致收敛于所讨论问题的唯一解的单调序列,而且是二阶收敛的。     

Convective nanofluid flow over a vertical cone with a rough surface

In recent literature, the analysis of a combined convective flow over a cone has received a lot of attention. To explore the convection effects of flow over a cone in greater detail, in this investigation, we have considered a cone with a rough surface, which is entirely a new flow problem. Recent studies have shown the influence of roughness on fluid flow over several geometries, but flow over a rough conical surface has not been studied so far. In addition, we have analyzed the effects of nanoparticles, magnetohydrodynamic (MHD), and suction/blowing, which could have significant impacts on characteristics of fluid flow over the cone with a rough surface. Initially, the governing equations, which are partial differential equations with a high degree of nonlinearity, are nondimensionalized through Mangler's transformations. Later, linear equations are obtained via the method of quasilinearization, which is then solved numerically through finite difference approximations. The roughness of the cone's surface has notable effects on fluid flow, that too away from the origin. In fact, the roughness increases the friction at the cone's surface. Furthermore, the magnetic field applied at the wall increases the surface friction. Thus, the combination of roughness and MHD helps delay the boundary layer separation. On the other hand, the suction reduces the temperature of the fluid and increases the energy transport strength, while the thermophoresis parameter exhibits the contrary nature. Therefore, the combined consideration of these two could enhance energy transport strength in several industrial applications.     

Thermally Induced Flow in Gas Centrifuge, (I)

Centrifugally driven thermal convection near the cooled rotor conver of a centrifuge is treated.

Boundary layer equations are reduced to ordinary differential equations by the Lees-Levy transformation, and solved numerically by applying the quasilinearization method.

Sample solutions for radial and azimuthal velocity profiles are presented. The values of the axial velocity at the rotor center are calculated as a function of the temperature difference between the rotor cover and the contained gas.     

基于保辛算法的绳系卫星系统闭环反馈控制问题求解

应用非线性系统滚动时域控制的保辛算法求解绳系卫星系统子星释放和回收过程的闭环反馈控制问题.通过第二类Lagrange方程推导出二体绳系卫星系统的动力学方程;通过拟线性化方法将绳系卫星系统闭环反馈控制问题转化为线性非齐次Hamilton系统两端边值问题的迭代求解;通过保辛算法将线性非齐次Hamilton两端边值问题转化为线性方程组的求解;通过递进更新时间步的状态变量和控制变量,完成绳系卫星系统的闭环反馈控制.数值仿真表明:相对于Legendre伪谱方法,用保辛算法求解绳系卫星系统的闭环反馈控制问题的计算速度和收敛速度较快.绳系卫星系统的开环控制和闭环反馈控制问题数值仿真结果表明:在绳系卫星的初始状态存在偏差的情况下,使用开环控制会导致系统在终端无法达到稳定状态,而使用闭环反馈控制则能在一段时间内抵消初始状态向量偏差对系统产生的影响,最终达到稳定状态.