An Advanced Galerkin Approach to Solve the Nonlinear \\[6pt]Reaction-Diffusion Equations With Different Boundary Conditions

This study proposed a scheme originated from the Galerkin finite element method (GFEM) for solving nonlinear parabolic partial differential equations (PDEs) numerically with initial and different types of boundary conditions. The scheme is applied generally handling the nonlinear terms in a simple way and throwing over restrictive assumptions. The convergence and stability analysis of the method are derived. The error of the method is estimated. In the series, eminent problems are solved, such as  Fisher's equation, Newell-Whitehead-Segel equation, Burger's equation, and  Burgers-Huxley equation to demonstrate the validity, efficiency, accuracy, simplicity and applicability of this scheme. In each example, the comparison results are presented both numerically and graphically

The effect of landscape fragmentation on Turing-pattern formation

<abstract><p>Diffusion-driven instability and Turing pattern formation are a well-known mechanism by which the local interaction of species, combined with random spatial movement, can generate stable patterns of population densities in the absence of spatial heterogeneity of the underlying medium. Some examples of such patterns exist in ecological interactions between predator and prey, but the conditions required for these patterns are not easily satisfied in ecological systems. At the same time, most ecological systems exist in heterogeneous landscapes, and landscape heterogeneity can affect species interactions and individual movement behavior. In this work, we explore whether and how landscape heterogeneity might facilitate Turing pattern formation in predator–prey interactions. We formulate reaction-diffusion equations for two interacting species on an infinite patchy landscape, consisting of two types of periodically alternating patches. Population dynamics and movement behavior differ between patch types, and individuals may have a preference for one of the two habitat types. We apply homogenization theory to derive an appropriately averaged model, to which we apply stability analysis for Turing patterns. We then study three scenarios in detail and find mechanisms by which diffusion-driven instabilities may arise even if the local interaction and movement rates do not indicate it.</p></abstract>

Bounds on the critical times for the Fisher-KPP equation

The Fisher–Kolmogorov–Petrovsky–Piskunov (Fisher–KPP) equation is one of the prototypical reaction–diffusion equations and is encountered in many areas, primarily in population dynamics. An important consideration for the phenomena modelled by diffusion equations is the length of the diffusive process. In this paper, three definitions of the critical time are given, and bounds are obtained by a careful construction of the upper and lower solutions. The comparison functions satisfy the nonlinear, but linearizable, partial differential equations of Fisher–KPP type. Results of the numerical simulations are displayed. Extensions to some classes of reaction–diffusion systems and an application to a spatially heterogeneous harvesting model are also presented. doi:10.1017/S1446181121000365

THREE-DIMENSIONAL COMPUTATIONAL MODEL OF EARLY UPPER LIMB DEVELOPMENT

Limb development begins during embryogenesis when a series of biochemical interactions are triggered between a particular region of the mesoderm and the ectoderm. These processes affect the morphogenesis and growth of bones, joints, and all the other constituent elements of limbs; nevertheless, how the biochemical regulation affects mesenchymal condensation is not entirely clear. In this study, a three-dimensional computational model is designed to predict the appearance and location of the mesenchymal condensation in the stylopod and zeugopod; the biochemical events were described with reaction–diffusion equations that were solved using the finite elements method. The result of the gene expression in our model was consistent with the one reported in literature; the obtained patterns of Fgf8, Fgf10, and Wnt3a can predict the shape of the mesenchymal condensation of early upper limb development; the simple diffusive patterns of molecules were suitable to explain the areas where sox9 is expressed. Furthermore, our results suggest that the expression of Tgf-[Formula: see text] in the upper limb could be due to the inhibition of retinoic acid. These results suggest the importance of building computational scenarios where pathologies may be comprehensively examined.

A Mathematical Model for Transport and Growth of Microbes in Unsaturated Porous Soil

In this work, we develop a mathematical model for transport and growth of microbes by natural (rain) water infiltration and flow through unsaturated porous soil along the vertical direction under gravity and capillarity by coupling a system of advection diffusion equations (for concentration of microbes and their growth-limiting substrate) with the Richards equation. The model takes into consideration several major physical, chemical, and biological mechanisms. The resulting coupled system of PDEs together with their boundary conditions is highly nonlinear and complicated to solve analytically. We present both a partial analytic approach towards solving the nonlinear system and finding the main type of dynamics of microbes, and a full-scale numerical simulation. Following the auxiliary equation method for nonlinear reaction-diffusion equations, we obtain a closed form traveling wave solution for the Richards equation. Using the propagating front solution for the pressure head, we reduce the transport equation to an ODE along the moving frame and obtain an analytic solution for the history of bacteria concentration for a specific test case. To solve the system numerically, we employ upwind finite volume method for the transport equations and stabilized explicit Runge–Kutta–Legendre super-time-stepping scheme for the Richards equation. Finally, some numerical simulation results of an infiltration experiment are presented, providing a validation and backup to the analytic partial solutions for the transport and growth of bacteria in the soil, stressing the occurrence of front moving solitons in the nonlinear dynamics.