Narendra Kumar,Lok Nath Kannaujiya,Amit K. Verma
Abstract:This paper introduces nonlinear fractional Lane-Emden equations of the form, $$ D^{\alpha} y(x) + \frac{\lambda}{x^\beta}~ D^{\beta} y(x) + f(y) =0, ~ ~1 < \alpha \leq 2, ~~ 0< \beta \leq 1, ~~ 0 < x < 1,$$ subject to boundary conditions, $$ y'(0) =\mathbf{a} , ~~~ \mathbf{c}~ y'(1) + \mathbf{d}~ y(1) = \mathbf{b},$$ where, $D^\alpha, D^\beta$ represent Caputo fractional derivative, $\mathbf{a, b,c,d} \in \mathbb{R}$, $ \lambda = 1, 2$, and $f(y)$ is non linear function of $y.$ We have developed collocation method namely, uniform fractional Haar wavelet collocation method and used it to compute solutions. The proposed method combines the quasilinearization method with the Haar wavelet collocation method. In this approach, fractional Haar integrations is used to determine the linear system, which, upon solving, produces the required solution. Our findings suggest that as the values of $(\alpha, \beta)$ approach $(2,1),$ the solutions of the fractional and classical Lane-Emden become identical.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a new numerical method for solving fractional - order regular β - singular boundary value problems (BVPs), which are used to simulate the heat conduction phenomenon in the human head under the action of antipyretic drugs. Specifically, the author introduced a class of nonlinear fractional - order Lane - Emden equations:
\[ D^\alpha y(x) + \lambda x^{-\beta} D^\beta y(x) + f(y) = 0, \quad 1 < \alpha \leq 2, \quad 0 < \beta \leq 1, \quad 0 < x < 1 \]
where \( D^\alpha \) and \( D^\beta \) represent Caputo fractional - order derivatives, \( a, b, c, d \in \mathbb{R} \), and the boundary conditions are:
\[ y'(0) = a, \quad cy'(1) + dy(1) = b \]
### Main research contents
1. **Propose a new method**: The author proposed the uniform fractional - order Haar wavelet collocation method (UFHWCM), combining the quasilinearization technique and the Haar wavelet collocation method to solve the above fractional - order Lane - Emden equations.
2. **Theoretical analysis**: The convergence and stability of the proposed method were analyzed, and it was proved that when \( \alpha \) is close to 2 and \( \beta \) is close to 1, the solution of the fractional - order Lane - Emden equation will approach the solution of the classical Lane - Emden equation.
3. **Numerical experiments**: The effectiveness and accuracy of the proposed method were verified through three test cases, and the residual errors and convergence rates under different parameter settings were shown.
### Research background
Fractional - order differential equations (FrDE) have been widely used in recent years to model phenomena in various applied research, especially since fractional - order derivatives can well capture memory and hereditary characteristics. The classical Lane - Emden equation mentioned in this paper has wide applications in astrophysics, engineering and other fields. However, relatively few studies have been conducted on fractional - order Lane - Emden equations, so this paper fills this gap in this field.
### Method innovation points
- **Introduce fractional - order Haar wavelets**: Utilize the orthogonality and compact support characteristics of Haar wavelets and apply them to the solution of fractional - order differential equations.
- **Combine with quasilinearization techniques**: Transform nonlinear problems into a series of linear problems through quasilinearization, thereby simplifying the calculation process.
### Numerical results
Through numerical experiments on different parameter combinations, the results show that as the resolution level \( J \) increases, the residual error gradually decreases, indicating that the proposed method has high precision and stability.
In summary, the main purpose of this paper is to effectively solve the fractional - order Lane - Emden equation by introducing a new numerical method - the uniform fractional - order Haar wavelet collocation method (UFHWCM), and verify its effectiveness and accuracy in simulating the heat conduction problem in the human head.