On the singular perturbations for fractional differential equation.
Ontology highlight
ABSTRACT:
SUBMITTER: Atangana A
PROVIDER: S-EPMC3934389 | biostudies-other | 2014
REPOSITORIES: biostudies-other
Similar Datasets
Project description:Fractional derivative or fractional calculus plays a significant role in theoretical modeling of scientific and engineering problems. However, only relatively low order fractional derivatives are used at present. In general, it is not obvious what role a high fractional derivative can play and how to make use of arbitrarily high-order fractional derivatives. This work introduces arbitrarily high-order fractional partial differential equations (PDEs) to describe fractional hyperdiffusions. The fractional PDEs are constructed via fractional variational principle. A fast fractional Fourier transform (FFFT) is proposed to numerically integrate the high-order fractional PDEs so as to avoid stringent stability constraints in solving high-order evolution PDEs. The proposed high-order fractional PDEs are applied to the surface generation of proteins. We first validate the proposed method with a variety of test examples in two and three-dimensional settings. The impact of high-order fractional derivatives to surface analysis is examined. We also construct fractional PDE transform based on arbitrarily high-order fractional PDEs. We demonstrate that the use of arbitrarily high-order derivatives gives rise to time-frequency localization, the control of the spectral distribution, and the regulation of the spatial resolution in the fractional PDE transform. Consequently, the fractional PDE transform enables the mode decomposition of images, signals, and surfaces. The effect of the propagation time on the quality of resulting molecular surfaces is also studied. Computational efficiency of the present surface generation method is compared with the MSMS approach in Cartesian representation. We further validate the present method by examining some benchmark indicators of macromolecular surfaces, i.e., surface area, surface enclosed volume, surface electrostatic potential and solvation free energy. Extensive numerical experiments and comparison with an established surface model indicate that the proposed high-order fractional PDEs are robust, stable and efficient for biomolecular surface generation.
| S-EPMC3869418 | biostudies-literature
Project description:To boost the handful of nutrient-dense individuals in the societal structure, adequate health care documentation and comprehension are permitted. This will strengthen and optimize the well-being of the community, particularly the girls and women of the community that are welcoming the new generation. In this article, we extensively explored a deterministic-stochastic malnutrition model involving nonlinear perturbation via piecewise fractional operators techniques. This novel concept leads us to analyze and predict the process from the beginning to the end of the well-being growth, as it offers the possibility to observe many behaviors from cross over to stochastic processes. Moreover, the piecewise differential operators, which can be constructed with operators such as classical, Caputo, Caputo-Fabrizio, Atangana-Baleanu and stochastic derivative. The threshold parameter is developed and the role of malnutrition in society is examined. Through a rigorous analysis, we first demonstrated that the stochastic model's solution is positive and global. Then, using appropriate stochastic Lyapunov candidates, we examined whether the stochastic system acknowledges a unique ergodic stationary distribution. The objective of this investigation is to design a nutritional deficiency in pregnant women using a piecewise fractional differential equation scheme. We examined multiple options and outlined numerical methods of coping with problems. To exemplify the effectiveness of the suggested concept, graphical conclusions, including chaotic and random perturbation patterns, are supplied. Consequently, fractional calculus' innovative aspects provide more powerful and flexible layouts, enabling us to more effectively adapt to the system dynamics tendencies of real-world representations. This has opened new doors to readers in different disciplines and enabled them to capture different behaviors at different time intervals.
| S-EPMC10491687 | biostudies-literature
Project description:Gene expression is controlled by multiple regulators and their interactions. Data from genome-wide gene expression assays can be used to estimate molecular activities of regulators within a model organism and extrapolate them to biological processes in humans. This approach is valuable in studies to better understand complex human biological systems which may be involved in diseases and hence, have potential clinical relevance. In order to achieve this, it is necessary to infer gene interactions that are not directly observed (i.e. latent or hidden) by way of structural equation modeling (SEM) on the expression levels or activities of the downstream targets of regulator genes. Here we developed an R Shiny application, termed "Structural Equation Modeling of In silico Perturbations (SEMIPs)" to compute a two-sided t-statistic (T-score) from analysis of gene expression data, as a surrogate to gene activity in a given human specimen. SEMIPs can be used in either correlational studies between outcome variables of interest or subsequent model fitting on multiple variables. This application implements a 3-node SEM model that consists of two upstream regulators as input variables and one downstream reporter as an outcome variable to examine the significance of interactions among these variables. SEMIPs enables scientists to investigate gene interactions among three variables through computational and mathematical modeling (i.e. in silico). In a case study using SEMIPs, we have shown that putative direct downstream genes of the GATA Binding Protein 2 (GATA2) transcription factor are sufficient to infer its activities in silico for the conserved progesterone receptor (PGR)-GATA2-SRY-box transcription factor 17 (SOX17) genetic network in the human uterine endometrium.
| S-EPMC8652139 | 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: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:Existing estimation methods for ordinary differential equation (ODE) models are not applicable to discrete data. The generalized ODE (GODE) model is therefore proposed and investigated for the first time. We develop the likelihood-based parameter estimation and inference methods for GODE models. We propose robust computing algorithms and rigorously investigate the asymptotic properties of the proposed estimator by considering both measurement errors and numerical errors in solving ODEs. The simulation study and application of our methods to an influenza viral dynamics study suggest that the proposed methods have a superior performance in terms of accuracy over the existing ODE model estimation approach and the extended smoothing-based (ESB) method.
| S-EPMC4274811 | biostudies-literature
Project description:The computational complexity of one-dimensional time fractional reaction-diffusion equation is O(N²M) compared with O(NM) for classical integer reaction-diffusion equation. Parallel computing is used to overcome this challenge. Domain decomposition method (DDM) embodies large potential for parallelization of the numerical solution for fractional equations and serves as a basis for distributed, parallel computations. A domain decomposition algorithm for time fractional reaction-diffusion equation with implicit finite difference method is proposed. The domain decomposition algorithm keeps the same parallelism but needs much fewer iterations, compared with Jacobi iteration in each time step. Numerical experiments are used to verify the efficiency of the obtained algorithm.
| S-EPMC3977434 | biostudies-other
Project description:Integrating genetic perturbations with gene expression data not only improves accuracy of regulatory network topology inference, but also enables learning of causal regulatory relations between genes. Although a number of methods have been developed to integrate both types of data, the desiderata of efficient and powerful algorithms still remains. In this paper, sparse structural equation models (SEMs) are employed to integrate both gene expression data and cis-expression quantitative trait loci (cis-eQTL), for modeling gene regulatory networks in accordance with biological evidence about genes regulating or being regulated by a small number of genes. A systematic inference method named sparsity-aware maximum likelihood (SML) is developed for SEM estimation. Using simulated directed acyclic or cyclic networks, the SML performance is compared with that of two state-of-the-art algorithms: the adaptive Lasso (AL) based scheme, and the QTL-directed dependency graph (QDG) method. Computer simulations demonstrate that the novel SML algorithm offers significantly better performance than the AL-based and QDG algorithms across all sample sizes from 100 to 1,000, in terms of detection power and false discovery rate, in all the cases tested that include acyclic or cyclic networks of 10, 30 and 300 genes. The SML method is further applied to infer a network of 39 human genes that are related to the immune function and are chosen to have a reliable eQTL per gene. The resulting network consists of 9 genes and 13 edges. Most of the edges represent interactions reasonably expected from experimental evidence, while the remaining may just indicate the emergence of new interactions. The sparse SEM and efficient SML algorithm provide an effective means of exploiting both gene expression and perturbation data to infer gene regulatory networks. An open-source computer program implementing the SML algorithm is freely available upon request.
| S-EPMC3662697 | biostudies-literature
Project description:Partial differential equation (PDE) models are commonly used to model complex dynamic systems in applied sciences such as biology and finance. The forms of these PDE models are usually proposed by experts based on their prior knowledge and understanding of the dynamic system. Parameters in PDE models often have interesting scientific interpretations, but their values are often unknown, and need to be estimated from the measurements of the dynamic system in the present of measurement errors. Most PDEs used in practice have no analytic solutions, and can only be solved with numerical methods. Currently, methods for estimating PDE parameters require repeatedly solving PDEs numerically under thousands of candidate parameter values, and thus the computational load is high. In this article, we propose two methods to estimate parameters in PDE models: a parameter cascading method and a Bayesian approach. In both methods, the underlying dynamic process modeled with the PDE model is represented via basis function expansion. For the parameter cascading method, we develop two nested levels of optimization to estimate the PDE parameters. For the Bayesian method, we develop a joint model for data and the PDE, and develop a novel hierarchical model allowing us to employ Markov chain Monte Carlo (MCMC) techniques to make posterior inference. Simulation studies show that the Bayesian method and parameter cascading method are comparable, and both outperform other available methods in terms of estimation accuracy. The two methods are demonstrated by estimating parameters in a PDE model from LIDAR data.
| S-EPMC3867159 | biostudies-literature
Project description:To better describe its constitutive relation, we need a new constitutive equation for an important nonlinear elastic material, Mn-Cu damping alloy. In this work, we studied the nonlinear and hysteretic characteristics of the stress-strain curve of the M2052 alloy with the uniaxial cyclic tensile test with constant strain rate. The strain rate and amplitude correlations of M2052 resembled those of nonlinear viscoelastic material. Therefore, we created a new constitutive equation for the M2052 damping alloy by modifying the fractional Maxwell model, and we used the genetic algorithm to carry out numerical fitting with MATLAB. By comparing with the experimental data, we confirmed that the new constitutive equation could accurately depict the nonlinear constitutive relation and hysteretic property of the damping alloy. Taken together, this new constitutive equation for Mn-Cu damping alloy based on the fractional Maxwell model can serve as an effective tool for further studies of the constitutive relation of the Mn-Cu damping alloys.
| S-EPMC7254341 | biostudies-literature