Nonlinear Boundary Value Problems Research Articles

Construction of high order numerical methods for solving fourth order nonlinear boundary value problems

Construction of high order numerical methods for solving fourth order nonlinear boundary value problems

The study of nonlinear fractional boundary value problems involving the p-Laplacian operator

Abstract The $p$-Laplacian has attracted considerable attention in numerous fields such as mechanics, image processing and game theory. It is a nonlinear operator which has been used in the modelling and qualitative aspects in numerous problems. In this research work, we propose a new nonlinear fractional differential equation involving the $p$-Laplacian, which include the generalized Caputo fractional derivatives. We investigate the existence and uniqueness of solutions to our proposed problem through the application using the Banach and Schauder's fixed-point theorems. Additionally, we illustrate the practical applicability of our findings by applying them to a specific example, thereby validating their efficacy.

Mathematical analysis of non-linear boundary-value problems in reaction-diffusion model of chitosan-alginate microsphere using asymptotic methods

This article describes the non-linear reaction-diffusion model, which depicts the behaviour of hydrogen peroxide production and glucose oxidation in the chitosan-alginate microsphere. Analytical solutions of glucose, oxygen, gluconic acid and hydrogen peroxide concentrations in planar coordinates under steady-state circumstances are obtained for all reaction-diffusion parameters. The non-linear boundary value problem's approximate analytical expressions are obtained using asymptotic techniques based on the homotopy perturbation and Akbari-Ganji method. The exact and commonly used MATLAB program provided a numerical simulation. It is demonstrated that the resulting analytical expressions strongly agree with the numerical outcomes reported in the literature. Furthermore, determining hydrogen peroxide flux and pH profiles within the microspheres are also discussed. The theoretical findings make it possible to forecast and enhance the efficiency of enzyme kinetics.

Steady-state triad resonance between a surface gravity wave and two hydroacoustic waves based on the homotopy analysis method

Under water compressibility, resonant triads can occur within the family of acoustic-gravity waves. This work investigates steady-state triad resonance between a surface gravity wave and two hydroacoustic waves. The water-wave equations are solved as a nonlinear boundary value problem using the homotopy analysis method (HAM). Within the HAM framework, a potential singularity resulting from exact triad resonance is avoided by appropriately choosing the auxiliary linear operator. The resonant hydroacoustic wave component, along with the two primary waves (one hydroacoustic wave and one gravity wave), is considered to determine an initial guess for the velocity potential. Additionally, by selecting an optimal “convergence-control parameter,” the steady-state resonance between a surface gravity wave and two hydroacoustic waves is successfully obtained. It is found that steady-state resonant acoustic-gravity waves are ubiquitous under certain circumstances. The two primary wave components and the resonant hydroacoustic wave component take up most of the energy in the steady-state resonant acoustic-gravity wave system. The amplitude of the resonant hydroacoustic wave component is mainly determined by the primary hydroacoustic wave component, and the amplitudes of both hydroacoustic waves are approximately equal in all cases considered. In addition, when the overall amplitude of the wave system increases, both dimensionless angular frequencies decrease, indicating that the nonlinearity of the entire wave system becomes stronger with an increase in the wave system amplitude. The amplitude of the primary hydroacoustic wave has a relatively large effect on the system's nonlinearity. This work will enrich and deepen our understanding of acoustic-gravity waves.

Seismic wave amplification in a weakly nonlinear two-layered soil deposit

Based on numerous observations, local soil site conditions can significantly affect the formation of seismic responses. In this paper, we examine seismic responses caused by shear waves propagating through a two-layered soil deposit. The aim of the research is to study the phenomenon of amplification of shear waves due to the resonant properties of the soil deposit and the nonlinear characteristics of partial layers, which are assumed to be the cubically nonlinear continuous media. The target problem relates to a nonlinear boundary value problem, the solution of which is carried out using a semi-analytical approach. According to the approach, the problem’s solution is sought using the method of separation of variables, while the solution to the corresponding spectral problem is derived in the form of a power series. Additionally, good agreement is shown between theoretical studies and numerical simulations. The developed approach allows one to evaluate the amplification factor for a nonlinear two-layered soil deposit and to detect changes in the shape of amplification factor profiles when the frequency of an incident wave on a rigid bedrock is varied.

On the existence of solutions for a class of nonlinear fully fourth-order differential systems

This paper discusses the existence of solutions to fourth order nonlinear boundary value problem involving systems of differential equations and two point boundary conditions. The nonlinearity f ∈ C ([0,1] × Rn+ × Rn × Rn × Rn, Rn)considered in this paper includes derivatives up to order three. Using recent fixed point results for the sum of two operators, we impose a growth condition on f to establish a new existence criteria that ensure the existence of at least one and the existence of at least two nonnegative solutions.

Training RBF neural networks for solving nonlinear and inverse boundary value problems

Radial basis function neural networks (RBFNN) have been increasingly employed to solve boundary value problems (BVPs). In the current study, we propose such a technique for nonlinear (apparently for the first time) elliptic BVPs of orders two and four in 2-D and 3-D. The method is also extended, in a natural way, to solving 2-D and 3-D inverse BVPs. The RBFNN is trained via the least squares minimization of a nonlinear functional using the MATLAB®routine lsqnonlin. In this way, as well as the solution, appropriate values of the RBF approximation parameters are automatically delivered. The efficacy of the proposed RBFNN is demonstrated through several numerical experiments.

Numerical study of binary mixture and thermophoretic analysis near a solar radiative heat transfer

The main focus of the current article is to investigate the two-dimensional boundary layer flow of non-Newtonian fluid induced by a linear stretching surface in the streamwise direction. The flow features of non-Newtonian fluid (i.e. shear-thinning and thickening nature of the fluids) are explored with Carreau fluid model. In addition, thermophoresis, Brownian motion, thermal radiation (with Rosseland approximation), binary mixture and activation enthalpy are also incorporated in fluid model for better analysis of thermal and solutal properties. The mathematical modelling of under consider physical problem yields nonlinear partial differential system. The governing mathematical system is first simplified with scaling transformations for the high Reynolds number and then transformed into non-dimensional ordinary differential system by using the similarity variables. Since the governing mathematical system comprised on nonlinear boundary value problem, thus shooting method is utilized for the numerical solution of the fluid flow governing equations. The obtained results followed the qualitative manners of the previously reported data. The computed results for important physical quantities (i.e. velocity, temperature, and concentration) are presented through graphs and then interpreted physically. It is observed that the boundary layer thickness increased for shear thickening fluids (n>1) while it is reduced for shear thinning fluids (n<1). Also, the temperature of shear thinning fluid is higher than shear thickening fluid. In addition, heat generation and thermal radiations enhance the thermal boundary layer. Thermophoretic viscosity effect reduces the concentration profile significantly in the flow domain.

A Nonlinear Mixed Finite Element Method for the Analysis of Flexoelectric Semiconductors

Abstract In this paper, we develop a nonlinear mixed finite element method for flexoelectric semiconductors and analyze the mechanically tuned redistributions of free carriers and electric currents through flexoelectric polarization in typical structures. We first present a macroscopic theory for flexoelectric semiconductors by combining flexoelectricity and nonlinear drift-diffusion theory. To use C0 continuous elements, we derive an incremental constrained weak form by introducing Langrage multipliers, in which the kinematic constraints between the displacement and its gradient are guaranteed. Based on the weak form, we established a mixed C0 continuous nine-node quadrilateral finite element as well as an iterative process for solving nonlinear boundary-value problems. The accuracy and convergence of the proposed element are validated by comparing linear finite element method results against analytical solutions for the bending of a beam. Finally, the nonlinear element method is applied to more complex problems, such as a circular ring, a plate with a hole, and an isosceles trapezoid. Results indicate that mechanical loads and doping levels have distinct influences on electric properties.

The improved residual power series method for the solutions of higher-order linear and nonlinear boundary value problems

AbstractIRPSM is the extended form of RPSM that gives an approximate solution to boundary value problems without requiring an exact solution. This method creates a truncated series for the determination of the missing initial conditions. The present study is conducted to evaluate semi-numerical solutions to differential equations using the Improved Residual Power Series Method (IRPSM). Our main objective is to check whether the proposed scheme is efficient for solving such equations or not. The solution of differential equations has been approximated using truncated residual power series. The method is used to solve differential equations for higher-order linear and non-linear boundary value problems. Absolute errors for the solved problems are calculated to ensure accuracy. It is also compared with some known methods, like the Differential Transform Method (DTM) and the Homotopy Perturbation Method (HPM). For the calculations, Mathematica 12.0 software is used. The computed results are compared to the exact solutions as well as the available results in the literature. The results showed that the results obtained by IRPSM are closely related to the exact solution when compared to the DTM and HPM results.

A note on Kirchhoff type boundary value problem involving Riemann-Liouville fractional derivative

In this paper, we study some nonlinear Kirchhoff boundary value problems of fractional differential equations involving Riemann Liouville operator. Under appropriate assumptions on the functions in the given problem, we establish the existence of solutions using variational methods combined with the mountain pass theorem. Moreover, an illustrative example is presented to prove the validity of the main result.

A Novel Meshfree Method for Nonlinear Equations in Flow through Porous Media and Electrohydrodynamic Flows

In this study, an efficient meshfree numerical method is introduced for solving the nonlinear boundary value problems. The method of fundamental solutions (MFS) is one of the most popular among meshfree methods. While traditionally limited to linear and homogeneous problems, this study extends the applicability of the MFS to include nonhomogeneous and nonlinear equations. To achieve this, an extended MFS is combined with a fixed-point iteration scheme. This developed framework is benchmarked to address two different flow problems. The first involves fluid flow through porous media in a channel governed by the nonlinear Brinkman-Forchheimer equation. The second problem pertains to electrohydrodynamic (EHD) flow in a circular conduit. The obtained solutions are compared with the finite element method and the solutions available in the existing literature.

Nonlinear transport and kinetics in electrocatalytic thin film of arbitrary shape for non-Michaelis-Menten reaction kinetics

The mathematical model for non-Michael-Menten kinetics, which describes a substrate that produces a complex with the immobilised catalyst, is discussed. This paper analytically solves the nonlinear reaction–diffusion equation in an electrocatalytic thin film of arbitrary shape. We provide a mathematical process that enables a comprehensive analytical solution to the nonlinear boundary value problem. Closed and simple forms of the approximate expression for substrate concentration profiles for general geometry (planar, cylindrical, and spherical) and corresponding steady-state amperometric current response are presented. The results obtained by three analytical methods are then compared with numerical solutions. The comparison has been represented graphically, and the Tabler form shows the effectiveness and advantages of the featured techniques. The effect of the parameters on concentration is also discussed.

Numerical solutions for magneto-convective boundary layer slip flow from a nonlinear stretching sheet with wall transpiration and thermal radiation effects

A mathematical model is presented for magnetohydrodynamic viscous incompressible flow of an electrically conducting Newtonian polymeric fluid from a moving permeable horizontal stretching sheet with variable magnetic field, momentum and thermal slip boundary conditions. The emerging nonlinear ordinary differential boundary value problem is solved numerically by the Runge-Kutta-Fehlberg fourth-fifth order numerical method in Maple symbolic software. The effects of the governing parameters on the dimensionless stream function, dimensionless flow velocity and the dimensionless temperature have been investigated, displayed graphically and discussed. It is found that greater momentum slip and magnetic field parameters reduce the dimensionless velocity. It is further observed that higher values of Prandtl number, thermal slip, and conductionradiation parameter (weaker radiation contribution) reduce the dimensionless temperature whilst stronger magnetic field increases the temperature due to the supplementary work expended in dragging the polymer against the action of the magnetic field which is dissipated as heat. Comparisons of numerical results with previous published results show an excellent agreement. The study finds applications in electro-conductive polymer (ECP) processing systems.

An Iterative Two-Grid Method for Strongly Nonlinear Elliptic Boundary Value Problems

An Iterative Two-Grid Method for Strongly Nonlinear Elliptic Boundary Value Problems

The stress state of a thick-walled shell with allowance for contact with a hydrogen-containing medium

Numerical and experimental methods are used to solve a multidisciplinary problem on determining the stress state of a steel shell of revolution under mechanical loading and thermal effect with allowance for its contact with a hydrogen-containing medium. The study uses a well-developed mathematical tool for solving heat conduction problems in order to solve the problem of hydrogen diffusion into metal. The effective stresses and their invariants are determined by solving the nonlinear boundary value problem of thermoplasticity of a thick-walled shell of revolution in a three-dimensional formulation. The study takes into account the experimentally found effect of changes in the mechanical properties of steel affected by hydrogen. The correctness of the proposed method and the performed calculations is quantitatively estimated by comparison with a well-known problem having an analytical solution. The paper shows that it is possible and necessary to take into account the change in mechanical properties when determining the stress state of steel structures operating in contact with a hydrogen-containing medium.

Solvability of a nonlinear second order m-point boundary value problem with p-Laplacian at resonance

We study the existence of solutions of the nonlinear second order m-point boundary value problem with p-Laplacian at resonance {(ϕp(x′))′=f(t,x,x′),t∈[0,1],x′(0)=0,x(1)=∑i=1m−2aix(ξi),\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \ extstyle\\begin{cases} (\\phi _{p}(x'))'=f(t,x,x'),\\quad t\\in [0,1],\\\\ x'(0)=0, \\qquad x(1)=\\sum_{i=1}^{m-2}a_{i}x(\\xi _{i}), \\end{cases} $$\\end{document} where ϕp(s)=|s|p−2s\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$\\phi _{p}(s)=|s|^{p-2}s$\\end{document}, p>1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$p>1$\\end{document}, f:[0,1]×R2→R\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$f:[0,1]\ imes \\mathbb{R}^{2}\ o \\mathbb{R}$\\end{document} is a continuous function, ai>0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$a_{i}>0$\\end{document} (i=1,2,…,m−2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$i=1,2,\\ldots ,m-2$\\end{document}) with ∑i=1m−2ai=1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$\\sum_{i=1}^{m-2}a_{i}=1$\\end{document}, 0<ξ1<ξ2<⋯<ξm−2<1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$0<\\xi _{1}<\\xi _{2}<\\cdots <\\xi _{m-2}<1$\\end{document}. Based on the topological transversality method together with the barrier strip technique and the cut-off technique, we obtain new existence results of solutions of the above problem. Meanwhile some examples are also given to illustrate our main results.

A Signed Maximum Principle for Boundary Value Problems for Riemann–Liouville Fractional Differential Equations with Analogues of Neumann or Periodic Boundary Conditions

Sufficient conditions are obtained for a signed maximum principle for boundary value problems for Riemann–Liouville fractional differential equations with analogues of Neumann or periodic boundary conditions in neighborhoods of simple eigenvalues. The primary objective is to exhibit four specific boundary value problems for which the sufficient conditions can be verified. To show an application of the signed maximum principle, a method of upper and lower solutions coupled with monotone methods is developed to obtain sufficient conditions for the existence of a maximal solution and a minimal solution of a nonlinear boundary value problem. A specific example is provided to show that sufficient conditions for the nonlinear problem can be realized.

A high-order B-spline collocation method for solving a class of nonlinear singular boundary value problems

A high-order B-spline collocation method for solving a class of nonlinear singular boundary value problems

Employing CNPS and CPS approaches to calculate numerical roots of ninth-order linear and nonlinear boundary value problems

The extremely accurate findings of ninth-order linear and nonlinear BVPs are described in this paper. CNPS and CPS are used to find out the numerical roots of linear and nonlinear, ninth-order BVPs. The proposed methods transform the ninth-order BVPs into a system of linear equations. The algorithms developed in this study not only provide approximate solutions to the BVPs using CNPS and CPS, but they also estimate the derivatives from the first to the ninth-order of the analytical solution simultaneously. These approaches are implemented across five diverse problems, showcasing the effectiveness of these techniques by utilizing step sizes of [Formula: see text] and [Formula: see text]. In order to gauge the accuracy of the methods, a comparison is drawn between the outcomes and the exact solutions, which are then presented in tabular and graphical format. The precision of these techniques is demonstrated through a detailed investigation and is found to be superior to the CBS, the Petrov–Galerkin Method (PGM) using Splines functions as basis functions and Septic B Splines functions as weight functions, the collocation method for ninth-order BVPs by Quintic B Splines and Sextic B Splines as evidenced by the comparison of AEs of CPS and CNPS with these alternative methods.