Nonlinear Diffusion Equation Research Articles

An α-robust and new two-grid nonuniform L2-1 FEM for nonlinear time-fractional diffusion equation

An α-robust and new two-grid nonuniform L2-1 FEM for nonlinear time-fractional diffusion equation

Implementation of thermal conduction energy transfer models in the Bifrost solar atmosphere MHD code

Thermal conductivity provides important contributions to the energy evolution of the upper solar atmosphere, behaving as a non-linear concentration-dependent diffusion equation. Recently, different methods have been offered as best-fit solutions to these problems in specific situations, but their effectiveness and limitations are rarely discussed. We have rigorously tested the different implementations of solving the conductivity flux, in the massively parallel magnetohydrodynamics code, Bifrost, with the aim of specifying the best scenarios for the use of each method. We compared the differences and limitations of explicit versus implicit methods, and analyse the convergence of a hyperbolic approximation. Among the tests, we used a newly derived first-order self-similar approximation to compare the efficacy of each method analytically in a 1D pure-thermal test scenario. We find that although the hyperbolic approximation proves the most accurate and the fastest to compute in long-running simulations, there is no optimal method to calculate the mid-term conductivity with both accuracy and efficiency. We also find that the solution of this approximation is sensitive to the initial conditions, and can lead to faster convergence if used correctly. Hyper-diffusivity is particularly useful in aiding the methods to perform optimally. We discuss recommendations for the use of each method within more complex simulations, whilst acknowledging the areas where there are opportunities for new methods to be developed.

Maximal $$L_p$$-regularity for x-dependent fractional heat equations with Dirichlet conditions

AbstractWe prove optimal regularity results in $$L_p$$ L p -based function spaces in space and time for a large class of linear parabolic equations with a nonlocal elliptic operator in bounded domains with limited smoothness. Here the nonlocal operator is given by a strongly elliptic and even pseudodifferential operator P of order 2a ($$0<a<1$$ 0 < a < 1 ) with nonsmooth x-dependent coefficients. This includes the prominent case of the fractional Laplacian $$(-\Delta )^a$$ ( - Δ ) a , as well as elliptic operators $$(-\nabla \cdot A(x)\nabla +b(x))^a$$ ( - ∇ · A ( x ) ∇ + b ( x ) ) a . The proofs are based on general results on maximal $$L_p$$ L p -regularity and its relation to $$\mathcal {R}$$ R -boundedness of the resolvent of the associated (elliptic) operator. Finally, we apply these results to show existence of strong solutions locally in time for a class of nonlinear nonlocal parabolic equations, which include a fractional nonlinear diffusion equation and a fractional porous medium equation after a transformation. The nonlinear results are new in the case of domains with boundary; the linear results are so when P is x-dependent nonsymmetric.

On multidimensional exact solutions of the nonlinear diffusion equation of the pantograph type with variable delay

We consider a multidimensional pantograph-type nonlinear diffusion equation with a linearly increasing time delay and scaling with respect to spatial variables in the source (sink). It is proposed to construct exact solutions by the reduction method using two ansatzes with a quadratic dependence on spatial variables. The dependence of the solution on spatial variables is found from a system of algebraic equations, and the dependence on time is found from a system of ordinary differential equations with a linearly increasing delay of the argument. A number of examples of exact solutions are given, both radially symmetric and anisotropic with respect to spatial variables.

The Second Critical Exponent for a Time-Fractional Reaction-Diffusion Equation

In this paper, we consider the Cauchy problem of a time-fractional nonlinear diffusion equation. According to Kaplan’s first eigenvalue method, we first prove the blow-up of the solutions in finite time under some sufficient conditions. We next provide sufficient conditions for the existence of global solutions by using the results of Zhang and Sun. In conclusion, we find the second critical exponent for the existence of global and non-global solutions via the decay rates of the initial data at spatial infinity.

Weak long‐range diffusion

AbstractWe study the spreading solutions of the nonlinear diffusion equation when the far‐field diffusivity is small. The method of strained coordinates is used to construct a uniform asymptotic correction to the similarity solution of the unperturbed problem. The equation provides a possible analogue to similar models of fluid jets and plumes.

Exact solution of the nonlinear boson diffusion equation for gluon scattering

An exact analytical solution of the nonlinear boson diffusion equation is presented. It accounts for the time evolution toward the Bose–Einstein equilibrium distribution through inelastic and elastic collisions in the case of constant transport coefficients. As a currently interesting application, gluon scattering in relativistic heavy-ion collisions is investigated. An estimate of the time-dependent gluon-condensate formation in overoccupied systems through number-conserving elastic scatterings in Pb–Pb collisions at relativistic energies is given.

Mathematical modeling and machine learning-based optimization for enhancing biofiltration efficiency of volatile organic compounds

Biofiltration is a method of pollution management that utilizes a bioreactor containing live material to absorb and destroy pollutants biologically. In this paper, we investigate mathematical models of biofiltration for mixing volatile organic compounds (VOCs) for instance hydrophilic (methanol) and hydrophobic (α\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha$$\\end{document}-pinene). The system of nonlinear diffusion equations describes the Michaelis-Menten kinetics of the enzymic chemical reaction. These models represent the chemical oxidation in the gas phase and mass transmission within the air-biofilm junction. Furthermore, for the numerical study of the saturation of α\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha$$\\end{document}-pinene and methanol in the biofilm and gas state, we have developed an efficient supervised machine learning algorithm based on the architecture of Elman neural networks (ENN). Moreover, the Levenberg-Marquardt (LM) optimization paradigm is used to find the parameters/ neurons involved in the ENN architecture. The approximation to a solutions found by the ENN-LM technique for methanol saturation and α\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha$$\\end{document}-pinene under variations in different physical parameters are allegorized with the numerical results computed by state-of-the-art techniques. The graphical and statistical illustration of indications of performance relative to the terms of absolute errors, mean absolute deviations, computational complexity, and mean square error validates that our results perfectly describe the real-life situation and can further be used for problems arising in chemical engineering.

Mathematical analysis of batch reactor performance for the enzymatic synthesis of cephalexin: Laplace Homotopy perturbation method and Adomian decomposition method

In this article, a mathematical model of diffusion reaction in kinetically controlled cephalexin synthesis in the batch reactor with penicillin acylase immobilized in glyoxyl-agarose is analyzed. The kinetic model is a non-linear non-stead-state reaction–diffusion equation with non-linear terms related to the Fick’s law. We have presented the approximate analytical expression for the non-steady-state non-linear reaction–diffusion equation by utilizing the Laplace homotopy perturbation method (LHPM) and for steady-state by using LHPM and Adomian decomposition method (ADM). The approximate analytical solution obtained from these methods proved that they are fit for every values of the reaction–diffusion and kinetic parameters. We also present the numerical solution of the considered reaction–diffusion equation by using pdepe tool in MATLAB software. When comparing the semi-analytical solution with the numerical solution, a satisfactory result is noted for all the possible values of the parameters. The closed-form analytical expressions for the effectiveness factor in non-steady and steady-state conditions are obtained and analyzed using LHPM and ADM.

High-precision meshfree methods for solving 3D radiation diffusion problem in irregular spaces

Radiation diffusion is a significant phenomenon in inertial confinement fusion, geology and so on. Solving nonlinear radiation diffusion equation in irregular spaces is challenging. In this paper, the direct linearisation algorithm and the successive permutation iteration algorithm are constructed for solving 3D radiation diffusion equation, which are full-implicit discretization on time. The high-precision meshfree methods are provided by the Kansa's non-symmetric collocation method with the compactly supported radial basis function. Ultimately, the accuracy and efficiency of the presented algorithms are verified through the comparison of the fixed point, Newton and SOR methods in 3D numerical experiments. These novel meshfree methods not only avoid the complexity of grid generation, but also provide high-precision numerical solution for nonlinear radiation diffusion equation.

Mathematical Analysis of Nonlinear Differential Equations in Polymer Coated Microelectrodes

A theoretical analysis is conducted on the electrochemical behaviour of micro disk electrodes coated with thin coatings of electroactive polymers. The main focus of efforts to characterize the planar diffusion chemical reaction within polymer-modified ultra-microelectrodes is the creation of a theoretical model that explicitly takes into account the potential for the polymer film to cover the inlaid micro disc support surface. Approximate formulae for the steady-state amperometric response and formulation of the boundary value problem are given. The impact of substrate concentration, mediator concentration and current responsiveness in the solution besides the polymer film is also investigated. Akbari-Ganji method (AGM) and differential transformation method (DTM), two adequate and widely available analytical methods, were utilized to determine the steady-state non-linear diffusion equation. The approximate analytical solution for the substrate concentration and mediator concentration and the current for the small experimental kinetic values and diffusion coefficients are presented. We additionally determine the problem’s numerical solution using the MATLAB tool. A satisfactory agreement can be seen when the numerical outcomes verify with the analytical findings.

Bayesian reduced-order deep learning surrogate model for dynamic systems described by partial differential equations

We propose a reduced-order deep-learning surrogate model for dynamic systems described by time-dependent partial differential equations. This method employs space–time Karhunen–Loève expansions (KLEs) of the state variables and space-dependent KLEs of space-varying parameters to identify the reduced (latent) dimensions. Subsequently, a deep neural network (DNN) is used to map the parameter latent space to the state variable latent space.An approximate Bayesian method is developed for uncertainty quantification (UQ) in the proposed KL-DNN surrogate model. The KL-DNN method is tested for the linear advection–diffusion and nonlinear diffusion equations, and the Bayesian approach for UQ is compared with the deep ensembling (DE) approach, commonly used for quantifying uncertainty in DNN models. It was found that the approximate Bayesian method provides a more informative distribution of the PDE solutions in terms of the coverage of the reference PDE solutions (the percentage of nodes where the reference solution is within the confidence interval predicted by the UQ methods) and log predictive probability. The DE method is found to underestimate uncertainty and introduce bias.For the nonlinear diffusion equation, we compare the KL-DNN method with the Fourier Neural Operator (FNO) method and find that KL-DNN is 10% more accurate and needs less training time than the FNO method.

The role of adhesion on soft lubrication: A new theory

Recent experiments reveal that adhesive interactions can play a key role in causing surface instability in soft lubrication. Instances of instability include fluid entrapment in isolated pockets upon a soft sphere's normal contact with a hard substrate and surface wrinkling of a soft substrate as a hard sphere slides across it. These phenomena underscore a substantial distinction between hard and soft lubrication. They are of paramount importance from a fundamental standpoint, providing an entirely new explanation for the transition mechanism from elasto-hydrodynamic to the mixed lubrication regimes. Here, we introduce a new theory to elucidate these observations. Our theory modifies the Reynolds elasto-hydrodynamic equation by incorporating adhesive interaction across the fluid layer, investigating the interplay between adhesion, fluid flow and elastic instability. Our analysis proposes the addition of a new dimensionless parameter in lubrication theory, that compares the stiffness of the adhesive interaction to that of the substrate. When this parameter exceeds unity, the soft solid surface exhibits instability to small perturbations in its shape. In mathematical terms, the Reynolds equation undergoes a transition from a nonlinear diffusion equation to a nonlinear wave equation at this critical point. Post-transition, the diffusivity of the nonlinear diffusion equation turns negative, rendering the problem ill-posed. We investigate the transition using the method of characteristics and present an exact analytic solution. This solution offers insights into the occurrence of a vanishing liquid film thickness at specific locations, resulting in dry contact—initiating transition to mixed lubrication.

Investigation of Well-Posedness for a Direct Problem for a Nonlinear Fractional Diffusion Equation and an Inverse Problem

Investigation of Well-Posedness for a Direct Problem for a Nonlinear Fractional Diffusion Equation and an Inverse Problem

Time-dependent concrete chloride diffusion coupled with nonlinear chloride binding: An analytical study

Both the temporal dependence of chloride diffusion coefficient and the effect of nonlinear chloride binding substantially affect the behavior of concrete chloride diffusion. This study develops an exact chloride diffusion solution that accommodates the above two factors by introducing a simplified but reasonably realistic nonlinear chloride binding isotherm. The introduction of this isotherm makes it possible to achieve an analytical solution to a highly nonlinear chloride diffusion equation that involves a free chloride concentration dependent dimensionless factor which reflects the effect of chloride binding. First, a time-dependent chloride diffusion model incorporating the introduced isotherm is established with the corresponding exact solution derived. The proposed solution is successfully validated against two existing analytical solutions and a numerical solution. Parametric study is conducted using the proposed solution to investigate how the parameters regarding chloride ingress influence the corrosion initiation time, and results show that the increase of the initial effective chloride diffusivity or the before-ingress age decreases the corrosion initiation time bringing an acceleration of a different level to chloride diffusion while the increase of the age factor prolongs the corrosion initiation time decelerating chloride diffusion to some extent, and as each of the three nonlinear chloride binding associated parameters (i.e., the deflection-point concentration, the former slope and the latter slope) increases, the corrosion initiation time rises and the retardation effect on chloride diffusion is enhanced to a various extent. Lastly, to demonstrate the practical validity of the proposed solution, the solution is applied to reproducing a reported chloride diffusion test, and it turns out that our solution reproduces the diffusion test reasonably well, justifying its practical validity.

Solutions with moving singularities for a one-dimensional nonlinear diffusion equation

AbstractThe aim of this paper is to study singular solutions for a one-dimensional nonlinear diffusion equation. Due to slow diffusion near singular points, there exists a solution with a singularity at a prescribed position depending on time. To study properties of such singular solutions, we define a minimal singular solution as a limit of a sequence of approximate solutions with large Dirichlet data. Applying the comparison principle and the intersection number argument, we discuss the existence and uniqueness of a singular solution for an initial-value problem, the profile near singular points and large-time behavior of solutions. We also give some results concerning the appearance of a burning core, convergence to traveling waves and the existence of an entire solution.

Solution of Fishers equation using homotopy perturbation Aboodh transform method

In this study, a nonlinear one-dimensional reaction Diffusion Equation called Fisher’s equation is solved using Aboodh Transform coupled with HPM. This method is applied to the general one-dimensional Fishers equation. This technique is demonstrated with three unique examples with different initial conditions solved with the values of β=1 and a=0. This simplicity of Aboodh transform makes the rapid rate of convergence from approximate to exact solution answers makes it easy and dependable.

Level set method via positive parts for optimal design problems of two-material thermal conductors

This paper is concerned with optimal design problems for two-material thermal conductors via level set methods based on (doubly) nonlinear diffusion equations. In level set methods with material representations via characteristic functions, gradient descent methods cannot be applied directly in terms of the differentiability of objective functionals with respect to level set functions, and therefore, appropriate sensitivities need to be constructed. This paper proposes a formulation via the positive parts of level set functions to avoid heuristic derivation of sensitivities and to apply (generalized gradient) descent methods. In particular, some perturbation term, such as a perimeter constraint, is involved in the formulation, and then an existence theorem for minimizers will be proved. Furthermore, convergence of objective functionals for minimizers with respect to a parameter of the perturbation term will also be discussed. In this paper, by deriving so-called weighted sensitivities, two-phase domains are numerically constructed as candidates for optimal configurations to approximate minimum values for classical design problems in two-dimensional cases.

Nonlocal approximation of nonlinear diffusion equations

We show that degenerate nonlinear diffusion equations can be asymptotically obtained as a limit from a class of nonlocal partial differential equations. The nonlocal equations are obtained as gradient flows of interaction-like energies approximating the internal energy. We construct weak solutions as the limit of a (sub)sequence of weak measure solutions by using the Jordan-Kinderlehrer-Otto scheme from the context of 2-Wasserstein gradient flows. Our strategy allows to cover the porous medium equation, for the general slow diffusion case, extending previous results in the literature. As a byproduct of our analysis, we provide a qualitative particle approximation.

Scaling and Localization in Multipole-Conserving Diffusion.

We study diffusion in systems of classical particles whose dynamics conserves the total center of mass. This conservation law leads to several interesting consequences. In finite systems, it allows for equilibrium distributions that are exponentially localized near system boundaries. It also yields an unusual approach to equilibrium, which in d dimensions exhibits scaling with dynamical exponent z=4+d. Similar phenomena occur for dynamics that conserves higher moments of the density, which we systematically classify using a family of nonlinear diffusion equations. In the quantum setting, analogous fermionic systems are shown to form real-space Fermi surfaces, while bosonic versions display a real-space analog of Bose-Einstein condensation.