A domain decomposition method for time fractional reaction-diffusion equation.
Ontology highlight
ABSTRACT:
SUBMITTER: Gong C
PROVIDER: S-EPMC3977434 | biostudies-other | 2014
REPOSITORIES: biostudies-other
Similar Datasets
Project description:This paper aims to obtain the approximate solution of time-fractional advection-dispersion equation (FADE) involving Jumarie's modification of Riemann-Liouville derivative by the fractional variational iteration method (FVIM). FVIM provides an analytical approximate solution in the form of a convergent series. Some examples are given and the results indicate that the FVIM is of high accuracy, more efficient, and more convenient for solving time FADEs.
| S-EPMC3919056 | biostudies-other
Project description:Spatial stochastic simulation is a valuable technique for studying reactions in biological systems. With the availability of high-performance computing (HPC), the method is poised to allow integration of data from structural, single-molecule and biochemical studies into coherent computational models of cells. Here, we introduce the Lattice Microbes software package for simulating such cell models on HPC systems. The software performs either well-stirred or spatially resolved stochastic simulations with approximated cytoplasmic crowding in a fast and efficient manner. Our new algorithm efficiently samples the reaction-diffusion master equation using NVIDIA graphics processing units and is shown to be two orders of magnitude faster than exact sampling for large systems while maintaining an accuracy of !0.1%. Display of cell models and animation of reaction trajectories involving millions of molecules is facilitated using a plug-in to the popular VMD visualization platform. The Lattice Microbes software is open source and available for download at http://www.scs.illinois.edu/schulten/lm
| S-EPMC3762454 | biostudies-other
Project description:Urticaria is a common skin disorder characterized by the rapid appearance and disappearance of local skin edema and flares with itching. It is characterized by various macroscopic skin eruptions unique to patients and/or subtypes of urticaria with respect to shape, size, color, and/or duration of eruptions. Nevertheless, the mechanism underlying multifarious eruptions in urticaria is largely unknown. The eruptions are believed to be evoked by histamine release from mast cells in the skin. However, the majority of visible characteristics of urticaria cannot be explained by a simple injection of histamine to the skin. To explain the multifarious eruptions of urticaria, we developed a single reaction-diffusion model suggesting the self-activation and self-inhibition regulation of histamine release from mast cells. Using the model, we found that various geometrical shapes of eruptions typically observed in patients can be explained by the model parameters and randomness or strength of the initial stimuli to mast cells. Furthermore, we verified that the wheal-expanding speed of urticaria, which is shown to be much smaller than that of the intradermal injection experimental system may be explained by our model and a simple diffusion equation. Our study suggests that the simple reaction-diffusion dynamics, including the independent self-activating and -inhibitory regulation of histamine release, may account for the essential mechanism underlying the formation of multifarious eruptions in urticaria.
| S-EPMC6961880 | biostudies-literature
Project description:Accurate and efficient forward models of photon migration in heterogeneous geometries are important for many applications of light in medicine because many biological tissues exhibit a layered structure of independent optical properties and thickness. However, closed form analytical solutions are not readily available for layered tissue-models, and often are modeled using computationally expensive numerical techniques or theoretical approximations that limit accuracy and real-time analysis. Here, we develop an open-source accurate, efficient, and stable numerical routine to solve the diffusion equation in the steady-state and time-domain for a layered cylinder tissue model with an arbitrary number of layers and specified thickness and optical coefficients. We show that the steady-state ([Formula: see text] ms) and time-domain ([Formula: see text] ms) fluence (for an 8-layer medium) can be calculated with absolute numerical errors approaching machine precision. The numerical implementation increased computation speed by 3 to 4 orders of magnitude compared to previously reported theoretical solutions in layered media. We verify our solutions asymptotically to homogeneous tissue geometries using closed form analytical solutions to assess convergence and numerical accuracy. Approximate solutions to compute the reflected intensity are presented which can decrease the computation time by an additional 2-3 orders of magnitude. We also compare our solutions for 2, 3, and 5 layered media to gold-standard Monte Carlo simulations in layered tissue models of high interest in biomedical optics (e.g. skin/fat/muscle and brain). The presented routine could enable more robust real-time data analysis tools in heterogeneous tissues that are important in many clinical applications such as functional brain imaging and diffuse optical spectroscopy.
| S-EPMC9643457 | biostudies-literature
Project description:The efficiency of exact simulation methods for the reaction-diffusion master equation (RDME) is severely limited by the large number of diffusion events if the mesh is fine or if diffusion constants are large. Furthermore, inherent properties of exact kinetic-Monte Carlo simulation methods limit the efficiency of parallel implementations. Several approximate and hybrid methods have appeared that enable more efficient simulation of the RDME. A common feature to most of them is that they rely on splitting the system into its reaction and diffusion parts and updating them sequentially over a discrete timestep. This use of operator splitting enables more efficient simulation but it comes at the price of a temporal discretization error that depends on the size of the timestep. So far, existing methods have not attempted to estimate or control this error in a systematic manner. This makes the solvers hard to use for practitioners since they must guess an appropriate timestep. It also makes the solvers potentially less efficient than if the timesteps are adapted to control the error. Here, we derive estimates of the local error and propose a strategy to adaptively select the timestep when the RDME is simulated via a first order operator splitting. While the strategy is general and applicable to a wide range of approximate and hybrid methods, we exemplify it here by extending a previously published approximate method, the Diffusive Finite-State Projection (DFSP) method, to incorporate temporal adaptivity.
| S-EPMC4746020 | biostudies-literature
Project description:Many biochemical processes at the sub-cellular level involve a small number of molecules. The local numbers of these molecules vary in space and time, and exhibit random fluctuations that can only be captured with stochastic simulations. We present a novel stochastic operator-splitting algorithm to model such reaction-diffusion phenomena. The reaction and diffusion steps employ stochastic simulation algorithms and Brownian dynamics, respectively. Through theoretical analysis, we have developed an algorithm to identify if the system is reaction-controlled, diffusion-controlled or is in an intermediate regime. The time-step size is chosen accordingly at each step of the simulation. We have used three examples to demonstrate the accuracy and robustness of the proposed algorithm. The first example deals with diffusion of two chemical species undergoing an irreversible bimolecular reaction. It is used to validate our algorithm by comparing its results with the solution obtained from a corresponding deterministic partial differential equation at low and high number of molecules. In this example, we also compare the results from our method to those obtained using a Gillespie multi-particle (GMP) method. The second example, which models simplified RNA synthesis, is used to study the performance of our algorithm in reaction- and diffusion-controlled regimes and to investigate the effects of local inhomogeneity. The third example models reaction-diffusion of CheY molecules through the cytoplasm of Escherichia coli during chemotaxis. It is used to compare the algorithm's performance against the GMP method. Our analysis demonstrates that the proposed algorithm enables accurate simulation of the kinetics of complex and spatially heterogeneous systems. It is also computationally more efficient than commonly used alternatives, such as the GMP method.
| S-EPMC3505198 | biostudies-other
Project description:The epidemiological aspects of the viral dynamic of the SARS-CoV-2 have become increasingly crucial due to major questions and uncertainties around the unaddressed issues of how corpse burial or the disposal of contaminated waste impacts nearby soil and groundwater. Here, a theoretical framework base on a meshless algorithm using the moving least squares (MLS) shape functions is adopted for solving the time-fractional model of the viral diffusion in and across three different environments including water, tissue, and soil. Our computations predict that by considering the α (order of fractional derivative) best fit to experimental data, the virus has a traveling distance of 1mm in water after 22, regardless of the source of contamination (e.g., from tissue or soil). The outcomes and extrapolations of our study are fundamental for providing valuable benchmarks for future experimentation on this topic and ultimately for the accurate description of viral spread across different environments. In addition to COVID-19 relief efforts, our methodology can be adapted for a wide range of applications such as studying virus ecology and genomic reservoirs in freshwater and marine environments.
| S-EPMC8825624 | biostudies-literature
Project description:Spatial reaction-diffusion models have been employed to describe many emergent phenomena in biological systems. The modelling technique most commonly adopted in the literature implements systems of partial differential equations (PDEs), which assumes there are sufficient densities of particles that a continuum approximation is valid. However, owing to recent advances in computational power, the simulation and therefore postulation, of computationally intensive individual-based models has become a popular way to investigate the effects of noise in reaction-diffusion systems in which regions of low copy numbers exist. The specific stochastic models with which we shall be concerned in this manuscript are referred to as 'compartment-based' or 'on-lattice'. These models are characterized by a discretization of the computational domain into a grid/lattice of 'compartments'. Within each compartment, particles are assumed to be well mixed and are permitted to react with other particles within their compartment or to transfer between neighbouring compartments. Stochastic models provide accuracy, but at the cost of significant computational resources. For models that have regions of both low and high concentrations, it is often desirable, for reasons of efficiency, to employ coupled multi-scale modelling paradigms. In this work, we develop two hybrid algorithms in which a PDE in one region of the domain is coupled to a compartment-based model in the other. Rather than attempting to balance average fluxes, our algorithms answer a more fundamental question: 'how are individual particles transported between the vastly different model descriptions?' First, we present an algorithm derived by carefully redefining the continuous PDE concentration as a probability distribution. While this first algorithm shows very strong convergence to analytical solutions of test problems, it can be cumbersome to simulate. Our second algorithm is a simplified and more efficient implementation of the first, it is derived in the continuum limit over the PDE region alone. We test our hybrid methods for functionality and accuracy in a variety of different scenarios by comparing the averaged simulations with analytical solutions of PDEs for mean concentrations.
| S-EPMC4424691 | biostudies-literature
Project description:We investigate the propagation of one-dimensional and two-dimensional (1D, 2D) Gaussian beams in the fractional Schrödinger equation (FSE) without a potential, analytically and numerically. Without chirp, a 1D Gaussian beam splits into two nondiffracting Gaussian beams during propagation, while a 2D Gaussian beam undergoes conical diffraction. When a Gaussian beam carries linear chirp, the 1D beam deflects along the trajectories z = ±2(x - x0), which are independent of the chirp. In the case of 2D Gaussian beam, the propagation is also deflected, but the trajectories align along the diffraction cone z = 2√(x(2) + y(2)) and the direction is determined by the chirp. Both 1D and 2D Gaussian beams are diffractionless and display uniform propagation. The nondiffracting property discovered in this model applies to other beams as well. Based on the nondiffracting and splitting properties, we introduce the Talbot effect of diffractionless beams in FSE.
| S-EPMC4838869 | biostudies-other
Project description:Reaction-diffusion models have been widely used to model glioma growth. However, it has not been shown how accurately this model can predict future tumor status using model parameters (i.e., tumor cell diffusion and proliferation) estimated from quantitative in vivo imaging data. To this end, we used in silico studies to develop the methods needed to accurately estimate tumor specific reaction-diffusion model parameters, and then tested the accuracy with which these parameters can predict future growth. The analogous study was then performed in a murine model of glioma growth. The parameter estimation approach was tested using an in silico tumor 'grown' for ten days as dictated by the reaction-diffusion equation. Parameters were estimated from early time points and used to predict subsequent growth. Prediction accuracy was assessed at global (total volume and Dice value) and local (concordance correlation coefficient, CCC) levels. Guided by the in silico study, rats (n = 9) with C6 gliomas, imaged with diffusion weighted magnetic resonance imaging, were used to evaluate the model's accuracy for predicting in vivo tumor growth. The in silico study resulted in low global (tumor volume error <8.8%, Dice >0.92) and local (CCC values >0.80) level errors for predictions up to six days into the future. The in vivo study showed higher global (tumor volume error >11.7%, Dice <0.81) and higher local (CCC <0.33) level errors over the same time period. The in silico study shows that model parameters can be accurately estimated and used to accurately predict future tumor growth at both the global and local scale. However, the poor predictive accuracy in the experimental study suggests the reaction-diffusion equation is an incomplete description of in vivo C6 glioma biology and may require further modeling of intra-tumor interactions including segmentation of (for example) proliferative and necrotic regions.
| S-EPMC4486062 | biostudies-literature