A method searching for optimum fractional order and its application in self-phase modulation induced nonlinear phase noise estimation in coherent optical fiber transmission systems

NASA Astrophysics Data System (ADS)

Huang, Chuan; Guo, Peng; Yang, Aiying; Qiao, Yaojun

2018-07-01

In single channel systems, the nonlinear phase noise only comes from the channel itself through self-phase modulation (SPM). In this paper, a fast-nonlinear effect estimation method is proposed based on fractional Fourier transformation (FrFT). The nonlinear phase noise caused by Self-phase modulation effect is accurately estimated for single model 10Gbaud OOK and RZ-QPSK signals with the fiber length range of 0-200 km and the launch power range of 1-10 mW. The pulse windowing is adopted to search the optimum fractional order for the OOK and RZ-QPSK signals. Since the nonlinear phase shift caused by the SPM effect is very small, the accurate optimum fractional order of the signal cannot be found based on the traditional method. In this paper, a new method magnifying the phase shift is proposed to get the accurate optimum order and thus the nonlinear phase shift is calculated. The simulation results agree with the theoretical analysis and the method is applicable to signals whose pulse type has the similar characteristics with Gaussian pulse.

The Effects of Core Composition on Iron Isotope Fractionation During Planetary Differentiation

NASA Astrophysics Data System (ADS)

Elardo, S. M.; Shahar, A.; Caracas, R.; Mock, T. D.; Sio, C. K. I.

2018-05-01

High pressure and temperature isotope exchange experiments and density functional theory calculations show how the composition of planetary cores affects the fractionation of iron isotopes during planetary differentiation.

Parareal algorithms with local time-integrators for time fractional differential equations

NASA Astrophysics Data System (ADS)

Wu, Shu-Lin; Zhou, Tao

2018-04-01

It is challenge work to design parareal algorithms for time-fractional differential equations due to the historical effect of the fractional operator. A direct extension of the classical parareal method to such equations will lead to unbalance computational time in each process. In this work, we present an efficient parareal iteration scheme to overcome this issue, by adopting two recently developed local time-integrators for time fractional operators. In both approaches, one introduces auxiliary variables to localized the fractional operator. To this end, we propose a new strategy to perform the coarse grid correction so that the auxiliary variables and the solution variable are corrected separately in a mixed pattern. It is shown that the proposed parareal algorithm admits robust rate of convergence. Numerical examples are presented to support our conclusions.

Exploring inductive linearization for pharmacokinetic-pharmacodynamic systems of nonlinear ordinary differential equations.

PubMed

Hasegawa, Chihiro; Duffull, Stephen B

2018-02-01

Pharmacokinetic-pharmacodynamic systems are often expressed with nonlinear ordinary differential equations (ODEs). While there are numerous methods to solve such ODEs these methods generally rely on time-stepping solutions (e.g. Runge-Kutta) which need to be matched to the characteristics of the problem at hand. The primary aim of this study was to explore the performance of an inductive approximation which iteratively converts nonlinear ODEs to linear time-varying systems which can then be solved algebraically or numerically. The inductive approximation is applied to three examples, a simple nonlinear pharmacokinetic model with Michaelis-Menten elimination (E1), an integrated glucose-insulin model and an HIV viral load model with recursive feedback systems (E2 and E3, respectively). The secondary aim of this study was to explore the potential advantages of analytically solving linearized ODEs with two examples, again E3 with stiff differential equations and a turnover model of luteinizing hormone with a surge function (E4). The inductive linearization coupled with a matrix exponential solution provided accurate predictions for all examples with comparable solution time to the matched time-stepping solutions for nonlinear ODEs. The time-stepping solutions however did not perform well for E4, particularly when the surge was approximated by a square wave. In circumstances when either a linear ODE is particularly desirable or the uncertainty in matching the integrator to the ODE system is of potential risk, then the inductive approximation method coupled with an analytical integration method would be an appropriate alternative.

Analytical and Numerical solutions of a nonlinear alcoholism model via variable-order fractional differential equations

NASA Astrophysics Data System (ADS)

Gómez-Aguilar, J. F.

2018-03-01

In this paper, we analyze an alcoholism model which involves the impact of Twitter via Liouville-Caputo and Atangana-Baleanu-Caputo fractional derivatives with constant- and variable-order. Two fractional mathematical models are considered, with and without delay. Special solutions using an iterative scheme via Laplace and Sumudu transform were obtained. We studied the uniqueness and existence of the solutions employing the fixed point postulate. The generalized model with variable-order was solved numerically via the Adams method and the Adams-Bashforth-Moulton scheme. Stability and convergence of the numerical solutions were presented in details. Numerical examples of the approximate solutions are provided to show that the numerical methods are computationally efficient. Therefore, by including both the fractional derivatives and finite time delays in the alcoholism model studied, we believe that we have established a more complete and more realistic indicator of alcoholism model and affect the spread of the drinking.

Optimal Energy Measurement in Nonlinear Systems: An Application of Differential Geometry

NASA Technical Reports Server (NTRS)

Fixsen, Dale J.; Moseley, S. H.; Gerrits, T.; Lita, A.; Nam, S. W.

2014-01-01

Design of TES microcalorimeters requires a tradeoff between resolution and dynamic range. Often, experimenters will require linearity for the highest energy signals, which requires additional heat capacity be added to the detector. This results in a reduction of low energy resolution in the detector. We derive and demonstrate an algorithm that allows operation far into the nonlinear regime with little loss in spectral resolution. We use a least squares optimal filter that varies with photon energy to accommodate the nonlinearity of the detector and the non-stationarity of the noise. The fitting process we use can be seen as an application of differential geometry. This recognition provides a set of well-developed tools to extend our work to more complex situations. The proper calibration of a nonlinear microcalorimeter requires a source with densely spaced narrow lines. A pulsed laser multi-photon source is used here, and is seen to be a powerful tool for allowing us to develop practical systems with significant detector nonlinearity. The combination of our analysis techniques and the multi-photon laser source create a powerful tool for increasing the performance of future TES microcalorimeters.

A homotopy analysis method for the nonlinear partial differential equations arising in engineering

NASA Astrophysics Data System (ADS)

Hariharan, G.

2017-05-01

In this article, we have established the homotopy analysis method (HAM) for solving a few partial differential equations arising in engineering. This technique provides the solutions in rapid convergence series with computable terms for the problems with high degree of nonlinear terms appearing in the governing differential equations. The convergence analysis of the proposed method is also discussed. Finally, we have given some illustrative examples to demonstrate the validity and applicability of the proposed method.

Lie symmetry analysis, exact solutions and conservation laws for the time fractional Caudrey-Dodd-Gibbon-Sawada-Kotera equation

NASA Astrophysics Data System (ADS)

Baleanu, Dumitru; Inc, Mustafa; Yusuf, Abdullahi; Aliyu, Aliyu Isa

2018-06-01

In this work, we investigate the Lie symmetry analysis, exact solutions and conservation laws (Cls) to the time fractional Caudrey-Dodd-Gibbon-Sawada-Kotera (CDGDK) equation with Riemann-Liouville (RL) derivative. The time fractional CDGDK is reduced to nonlinear ordinary differential equation (ODE) of fractional order. New exact traveling wave solutions for the time fractional CDGDK are obtained by fractional sub-equation method. In the reduced equation, the derivative is in Erdelyi-Kober (EK) sense. Ibragimov's nonlocal conservation method is applied to construct Cls for time fractional CDGDK.

Fractional calculus phenomenology in two-dimensional plasma models

NASA Astrophysics Data System (ADS)

Gustafson, Kyle; Del Castillo Negrete, Diego; Dorland, Bill

2006-10-01

Transport processes in confined plasmas for fusion experiments, such as ITER, are not well-understood at the basic level of fully nonlinear, three-dimensional kinetic physics. Turbulent transport is invoked to describe the observed levels in tokamaks, which are orders of magnitude greater than the theoretical predictions. Recent results show the ability of a non-diffusive transport model to describe numerical observations of turbulent transport. For example, resistive MHD modeling of tracer particle transport in pressure-gradient driven turbulence for a three-dimensional plasma reveals that the superdiffusive (2̂˜t^α where α> 1) radial transport in this system is described quantitatively by a fractional diffusion equation Fractional calculus is a generalization involving integro-differential operators, which naturally describe non-local behaviors. Our previous work showed the quantitative agreement of special fractional diffusion equation solutions with numerical tracer particle flows in time-dependent linearized dynamics of the Hasegawa-Mima equation (for poloidal transport in a two-dimensional cold-ion plasma). In pursuit of a fractional diffusion model for transport in a gyrokinetic plasma, we now present numerical results from tracer particle transport in the nonlinear Hasegawa-Mima equation and a planar gyrokinetic model. Finite Larmor radius effects will be discussed. D. del Castillo Negrete, et al, Phys. Rev. Lett. 94, 065003 (2005).

Regarding on the prototype solutions for the nonlinear fractional-order biological population model

DOE Office of Scientific and Technical Information (OSTI.GOV)

Baskonus, Haci Mehmet, E-mail: hmbaskonus@gmail.com; Bulut, Hasan

2016-06-08

In this study, we have submitted to literature a method newly extended which is called as Improved Bernoulli sub-equation function method based on the Bernoulli Sub-ODE method. The proposed analytical scheme has been expressed with steps. We have obtained some new analytical solutions to the nonlinear fractional-order biological population model by using this technique. Two and three dimensional surfaces of analytical solutions have been drawn by wolfram Mathematica 9. Finally, a conclusion has been submitted by mentioning important acquisitions founded in this study.

Nonlinear dynamics under varying temperature conditions of the resonating beams of a differential resonant accelerometer

NASA Astrophysics Data System (ADS)

Zhang, Jing; Wang, Yagang; Zega, Valentina; Su, Yan; Corigliano, Alberto

2018-07-01

In this work the nonlinear dynamic behaviour under varying temperature conditions of the resonating beams of a differential resonant accelerometer is studied from the theoretical, numerical and experimental points of view. A complete analytical model based on the Hamilton’s principle is proposed to describe the nonlinear behaviour of the resonators under varying temperature conditions and numerical solutions are presented in comparison with experimental data. This provides a novel perspective to examine the relationship between temperature and nonlinearity, which helps predicting the dynamic behaviour of resonant devices and can guide their optimal design.

Monitoring inter-channel nonlinearity based on differential pilot

NASA Astrophysics Data System (ADS)

Wang, Wanli; Yang, Aiying; Guo, Peng; Lu, Yueming; Qiao, Yaojun

2018-06-01

We modify and simplify the inter-channel nonlinearity (NL) estimation method by using differential pilot. Compared to previous works, the inter-channel NL estimation method we propose has much lower complexity and does not need modification of the transmitter. The performance of inter-channel NL monitoring with different launch power is tested. For both QPSK and 16QAM systems with 9 channels, the estimation error of inter-channel NL is lower than 1 dB when the total launch power is bigger than 12 dBm after 1000 km optical transmission. At last, we compare our inter-channel NL estimation method with other methods.

Synthesis of robust nonlinear autopilots using differential game theory

NASA Technical Reports Server (NTRS)

Menon, P. K. A.

1991-01-01

A synthesis technique for handling unmodeled disturbances in nonlinear control law synthesis was advanced using differential game theory. Two types of modeling inaccuracies can be included in the formulation. The first is a bias-type error, while the second is the scale-factor-type error in the control variables. The disturbances were assumed to satisfy an integral inequality constraint. Additionally, it was assumed that they act in such a way as to maximize a quadratic performance index. Expressions for optimal control and worst-case disturbance were then obtained using optimal control theory.

On new classes of solutions of nonlinear partial differential equations in the form of convergent special series

NASA Astrophysics Data System (ADS)

Filimonov, M. Yu.

2017-12-01

The method of special series with recursively calculated coefficients is used to solve nonlinear partial differential equations. The recurrence of finding the coefficients of the series is achieved due to a special choice of functions, in powers of which the solution is expanded in a series. We obtain a sequence of linear partial differential equations to find the coefficients of the series constructed. In many cases, one can deal with a sequence of linear ordinary differential equations. We construct classes of solutions in the form of convergent series for a certain class of nonlinear evolution equations. A new class of solutions of generalized Boussinesque equation with an arbitrary function in the form of a convergent series is constructed.

Novel Numerical Methods for Optimal Control Problems Involving Fractional-Order Differential Equations

DTIC Science & Technology

2018-03-14

pricing, Appl. Math . Comp. Vol.305, 174-187 (2017) 5. W. Li, S. Wang, Pricing European options with proportional transaction costs and stochastic...for fractional differential equation. Numer. Math . Theor. Methods Appl. 5, 229–241, 2012. [23] Kilbas A.A. and Marzan, S.A., Cauchy problem for...numerical technique for solving fractional optimal control problems, Comput. Math . Appl., 62, Issue 3, 1055–1067, 2011. [26] Lotfi A., Yousefi SA., Dehghan M

High-order fractional partial differential equation transform for molecular surface construction.

PubMed

Hu, Langhua; Chen, Duan; Wei, Guo-Wei

2013-01-01

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

High-order fractional partial differential equation transform for molecular surface construction

PubMed Central

Hu, Langhua; Chen, Duan; Wei, Guo-Wei

2013-01-01

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

A method for exponential propagation of large systems of stiff nonlinear differential equations

NASA Technical Reports Server (NTRS)

Friesner, Richard A.; Tuckerman, Laurette S.; Dornblaser, Bright C.; Russo, Thomas V.

1989-01-01

A new time integrator for large, stiff systems of linear and nonlinear coupled differential equations is described. For linear systems, the method consists of forming a small (5-15-term) Krylov space using the Jacobian of the system and carrying out exact exponential propagation within this space. Nonlinear corrections are incorporated via a convolution integral formalism; the integral is evaluated via approximate Krylov methods as well. Gains in efficiency ranging from factors of 2 to 30 are demonstrated for several test problems as compared to a forward Euler scheme and to the integration package LSODE.

Maglev Train Signal Processing Architecture Based on Nonlinear Discrete Tracking Differentiator.

PubMed

Wang, Zhiqiang; Li, Xiaolong; Xie, Yunde; Long, Zhiqiang

2018-05-24

In a maglev train levitation system, signal processing plays an important role for the reason that some sensor signals are prone to be corrupted by noise due to the harsh installation and operation environment of sensors and some signals cannot be acquired directly via sensors. Based on these concerns, an architecture based on a new type of nonlinear second-order discrete tracking differentiator is proposed. The function of this signal processing architecture includes filtering signal noise and acquiring needed signals for levitation purposes. The proposed tracking differentiator possesses the advantages of quick convergence, no fluttering, and simple calculation. Tracking differentiator's frequency characteristics at different parameter values are studied in this paper. The performance of this new type of tracking differentiator is tested in a MATLAB simulation and this tracking-differentiator is implemented in Very-High-Speed Integrated Circuit Hardware Description Language (VHDL). In the end, experiments are conducted separately on a test board and a maglev train model. Simulation and experiment results show that the performance of this novel signal processing architecture can fulfill the real system requirement.

From stochastic processes to numerical methods: A new scheme for solving reaction subdiffusion fractional partial differential equations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Angstmann, C.N.; Donnelly, I.C.; Henry, B.I., E-mail: B.Henry@unsw.edu.au

We have introduced a new explicit numerical method, based on a discrete stochastic process, for solving a class of fractional partial differential equations that model reaction subdiffusion. The scheme is derived from the master equations for the evolution of the probability density of a sum of discrete time random walks. We show that the diffusion limit of the master equations recovers the fractional partial differential equation of interest. This limiting procedure guarantees the consistency of the numerical scheme. The positivity of the solution and stability results are simply obtained, provided that the underlying process is well posed. We also showmore » that the method can be applied to standard reaction–diffusion equations. This work highlights the broader applicability of using discrete stochastic processes to provide numerical schemes for partial differential equations, including fractional partial differential equations.« less

Unsteady Solution of Non-Linear Differential Equations Using Walsh Function Series

NASA Technical Reports Server (NTRS)

Gnoffo, Peter A.

2015-01-01

Walsh functions form an orthonormal basis set consisting of square waves. The discontinuous nature of square waves make the system well suited for representing functions with discontinuities. The product of any two Walsh functions is another Walsh function - a feature that can radically change an algorithm for solving non-linear partial differential equations (PDEs). The solution algorithm of non-linear differential equations using Walsh function series is unique in that integrals and derivatives may be computed using simple matrix multiplication of series representations of functions. Solutions to PDEs are derived as functions of wave component amplitude. Three sample problems are presented to illustrate the Walsh function series approach to solving unsteady PDEs. These include an advection equation, a Burgers equation, and a Riemann problem. The sample problems demonstrate the use of the Walsh function solution algorithms, exploiting Fast Walsh Transforms in multi-dimensions (O(Nlog(N))). Details of a Fast Walsh Reciprocal, defined here for the first time, enable inversion of aWalsh Symmetric Matrix in O(Nlog(N)) operations. Walsh functions have been derived using a fractal recursion algorithm and these fractal patterns are observed in the progression of pairs of wave number amplitudes in the solutions. These patterns are most easily observed in a remapping defined as a fractal fingerprint (FFP). A prolongation of existing solutions to the next highest order exploits these patterns. The algorithms presented here are considered a work in progress that provide new alternatives and new insights into the solution of non-linear PDEs.

A nonlinear ordinary differential equation associated with the quantum sojourn time

NASA Astrophysics Data System (ADS)

Benguria, Rafael D.; Duclos, Pierre; Fernández, Claudio; Sing-Long, Carlos

2010-11-01

We study a nonlinear ordinary differential equation on the half-line, with the Dirichlet boundary condition at the origin. This equation arises when studying the local maxima of the sojourn time for a free quantum particle whose states belong to an adequate subspace of the unit sphere of the corresponding Hilbert space. We establish several results concerning the existence and asymptotic behavior of the solutions.

The Non-Linear Relationship between BMI and Health Care Costs and the Resulting Cost Fraction Attributable to Obesity.

PubMed

Laxy, Michael; Stark, Renée; Peters, Annette; Hauner, Hans; Holle, Rolf; Teuner, Christina M

2017-08-30

This study aims to analyse the non-linear relationship between Body Mass Index (BMI) and direct health care costs, and to quantify the resulting cost fraction attributable to obesity in Germany. Five cross-sectional surveys of cohort studies in southern Germany were pooled, resulting in data of 6757 individuals (31-96 years old). Self-reported information on health care utilisation was used to estimate direct health care costs for the year 2011. The relationship between measured BMI and annual costs was analysed using generalised additive models, and the cost fraction attributable to obesity was calculated. We found a non-linear association of BMI and health care costs with a continuously increasing slope for increasing BMI without any clear threshold. Under the consideration of the non-linear BMI-cost relationship, a shift in the BMI distribution so that the BMI of each individual is lowered by one point is associated with a 2.1% reduction of mean direct costs in the population. If obesity was eliminated, and the BMI of all obese individuals were lowered to 29.9 kg/m², this would reduce the mean direct costs by 4.0% in the population. Results show a non-linear relationship between BMI and health care costs, with very high costs for a few individuals with high BMI. This indicates that population-based interventions in combination with selective measures for very obese individuals might be the preferred strategy.

Analyzing Lie symmetry and constructing conservation laws for time-fractional Benny-Lin equation

NASA Astrophysics Data System (ADS)

Rashidi, Saeede; Hejazi, S. Reza

This paper investigates the invariance properties of the time fractional Benny-Lin equation with Riemann-Liouville and Caputo derivatives. This equation can be reduced to the Kawahara equation, fifth-order Kdv equation, the Kuramoto-Sivashinsky equation and Navier-Stokes equation. By using the Lie group analysis method of fractional differential equations (FDEs), we derive Lie symmetries for the Benny-Lin equation. Conservation laws for this equation are obtained with the aid of the concept of nonlinear self-adjointness and the fractional generalization of the Noether’s operators. Furthermore, by means of the invariant subspace method, exact solutions of the equation are also constructed.

A mass-energy preserving Galerkin FEM for the coupled nonlinear fractional Schrödinger equations

NASA Astrophysics Data System (ADS)

Zhang, Guoyu; Huang, Chengming; Li, Meng

2018-04-01

We consider the numerical simulation of the coupled nonlinear space fractional Schrödinger equations. Based on the Galerkin finite element method in space and the Crank-Nicolson (CN) difference method in time, a fully discrete scheme is constructed. Firstly, we focus on a rigorous analysis of conservation laws for the discrete system. The definitions of discrete mass and energy here correspond with the original ones in physics. Then, we prove that the fully discrete system is uniquely solvable. Moreover, we consider the unconditionally convergent properties (that is to say, we complete the error estimates without any mesh ratio restriction). We derive L2-norm error estimates for the nonlinear equations and L^{∞}-norm error estimates for the linear equations. Finally, some numerical experiments are included showing results in agreement with the theoretical predictions.

Fractional Gaussian noise-enhanced information capacity of a nonlinear neuron model with binary signal input

NASA Astrophysics Data System (ADS)

Gao, Feng-Yin; Kang, Yan-Mei; Chen, Xi; Chen, Guanrong

2018-05-01

This paper reveals the effect of fractional Gaussian noise with Hurst exponent H ∈(1 /2 ,1 ) on the information capacity of a general nonlinear neuron model with binary signal input. The fGn and its corresponding fractional Brownian motion exhibit long-range, strong-dependent increments. It extends standard Brownian motion to many types of fractional processes found in nature, such as the synaptic noise. In the paper, for the subthreshold binary signal, sufficient conditions are given based on the "forbidden interval" theorem to guarantee the occurrence of stochastic resonance, while for the suprathreshold binary signal, the simulated results show that additive fGn with Hurst exponent H ∈(1 /2 ,1 ) could increase the mutual information or bits count. The investigation indicated that the synaptic noise with the characters of long-range dependence and self-similarity might be the driving factor for the efficient encoding and decoding of the nervous system.

A Fractional Differential Kinetic Equation and Applications to Modelling Bursts in Turbulent Nonlinear Space Plasmas

NASA Astrophysics Data System (ADS)

Watkins, N. W.; Rosenberg, S.; Sanchez, R.; Chapman, S. C.; Credgington, D.

2008-12-01

Since the 1960s Mandelbrot has advocated the use of fractals for the description of the non-Euclidean geometry of many aspects of nature. In particular he proposed two kinds of model to capture persistence in time (his Joseph effect, common in hydrology and with fractional Brownian motion as the prototype) and/or prone to heavy tailed jumps (the Noah effect, typical of economic indices, for which he proposed Lévy flights as an exemplar). Both effects are now well demonstrated in space plasmas, notably in the turbulent solar wind. Models have, however, typically emphasised one of the Noah and Joseph parameters (the Lévy exponent μ and the temporal exponent β) at the other's expense. I will describe recent work in which we studied a simple self-affine stable model-linear fractional stable motion, LFSM, which unifies both effects and present a recently-derived diffusion equation for LFSM. This replaces the second order spatial derivative in the equation of fBm with a fractional derivative of order μ, but retains a diffusion coefficient with a power law time dependence rather than a fractional derivative in time. I will also show work in progress using an LFSM model and simple analytic scaling arguments to study the problem of the area between an LFSM curve and a threshold. This problem relates to the burst size measure introduced by Takalo and Consolini into solar-terrestrial physics and further studied by Freeman et al [PRE, 2000] on solar wind Poynting flux near L1. We test how expressions derived by other authors generalise to the non-Gaussian, constant threshold problem. Ongoing work on extension of these LFSM results to multifractals will also be discussed.

Experimentally determined sulfur isotope fractionation between metal and silicate and implications for planetary differentiation

NASA Astrophysics Data System (ADS)

Labidi, J.; Shahar, A.; Le Losq, C.; Hillgren, V. J.; Mysen, B. O.; Farquhar, J.

2016-02-01

The Earth's mantle displays a subchondritic 34S/32S ratio. Sulfur is a moderately siderophile element (i.e. iron-loving), and its partitioning into the Earth's core may have left such a distinctive isotope composition on the terrestrial mantle. In order to constrain the sulfur isotope fractionation occurring during core-mantle differentiation, high-pressure and temperature experiments were conducted with synthetic mixtures of metal and silicate melts. With the purpose to identify the mechanism(s) responsible for the S isotope fractionations, we performed our experiments in different capsules - namely, graphite and boron nitride capsules - and thus at different fO2, with varying major element chemistry of the silicate and metal fractions. The S isotope fractionations Δ34Smetal-silicate of equilibrated metal alloys versus silicate melts is +0.2 ± 0.1‰ in a boron-free and aluminum-poor system quenched at 1-1.5 GPa and 1650 °C. The isotope fractionation increases linearly with increasing boron and aluminum content, up to +1.4 ± 0.2‰, and is observed to be independent of the silicon abundance as well as of the fO2 over ∼3.5 log units of variations explored here. The isotope fractionations are also independent of the graphite or nitride saturation of the metal. Only the melt structural changes associated with aluminum and boron concentration in silicate melts have been observed to affect the strength of sulfur bonding. These results establish that the structure of silicate melts has a direct influence on the S2- average bonding strengths. These results can be interpreted in the context of planetary differentiation. Indeed, the structural environments of silicate evolve strongly with pressure. For example, the aluminum, iron or silicon coordination numbers increase under the effect of pressure. Consequently, based on our observations, the sulfur-bonding environment is likely to be affected. In this scheme, we tentatively hypothesize that S isotope fractionations

Extremely low order time-fractional differential equation and application in combustion process

NASA Astrophysics Data System (ADS)

Xu, Qinwu; Xu, Yufeng

2018-11-01

Fractional blow-up model, especially which is of very low order of fractional derivative, plays a significant role in combustion process. The order of time-fractional derivative in diffusion model essentially distinguishes the super-diffusion and sub-diffusion processes when it is relatively high or low accordingly. In this paper, the blow-up phenomenon and condition of its appearance are theoretically proved. The blow-up moment is estimated by using differential inequalities. To numerically study the behavior around blow-up point, a mixed numerical method based on adaptive finite difference on temporal direction and highly effective discontinuous Galerkin method on spatial direction is proposed. The time of blow-up is calculated accurately. In simulation, we analyze the dynamics of fractional blow-up model under different orders of fractional derivative. It is found that the lower the order, the earlier the blow-up comes, by fixing the other parameters in the model. Our results confirm the physical truth that a combustor for explosion cannot be too small.

A nonlinear and fractional derivative viscoelastic model for rail pads in the dynamic analysis of coupled vehicle-slab track systems

NASA Astrophysics Data System (ADS)

Zhu, Shengyang; Cai, Chengbiao; Spanos, Pol D.

2015-01-01

A nonlinear and fractional derivative viscoelastic (FDV) model is used to capture the complex behavior of rail pads. It is implemented into the dynamic analysis of coupled vehicle-slab track (CVST) systems. The vehicle is treated as a multi-body system with 10 degrees of freedom, and the slab track is represented by a three layer Bernoulli-Euler beam model. The model for the rail pads is one dimensional, and the force-displacement relation is based on a superposition of elastic, friction, and FDV forces. This model takes into account the influences of the excitation frequency and of the displacement amplitude through a fractional derivative element, and a nonlinear friction element, respectively. The Grünwald representation of the fractional derivatives is employed to numerically solve the fractional and nonlinear equations of motion of the CVST system by means of an explicit integration algorithm. A dynamic analysis of the CVST system exposed to excitations of rail harmonic irregularities is carried out, pointing out the stiffness and damping dependence on the excitation frequency and the displacement amplitude. The analysis indicates that the dynamic stiffness and damping of the rail pads increase with the excitation frequency while they decrease with the displacement amplitude. Furthermore, comparisons between the proposed model and ordinary Kelvin model adopted for the CVST system, under excitations of welded rail joint irregularities and of random track irregularities, are conducted in the time domain as well as in the frequency domain. The proposed model is shown to possess several modeling advantages over the ordinary Kelvin element which overestimates both the stiffness and damping features at high frequencies.

Reproducibility of Differential Proteomic Technologies in CPTAC Fractionated Xenografts

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tabb, David L.; Wang, Xia; Carr, Steven A.

2016-03-04

The NCI Clinical Proteomic Tumor Analysis Consortium (CPTAC) employed a pair of reference xenograft proteomes for initial platform validation and ongoing quality control of its data collection for The Cancer Genome Atlas (TCGA) tumors. These two xenografts, representing basal and luminal-B human breast cancer, were fractionated and analyzed on six mass spectrometers in a total of 46 replicates divided between iTRAQ and label-free technologies, spanning a total of 1095 LC-MS/MS experiments. These data represent a unique opportunity to evaluate the stability of proteomic differentiation by mass spectrometry over many months of time for individual instruments or across instruments running dissimilarmore » workflows. We evaluated iTRAQ reporter ions, label-free spectral counts, and label-free extracted ion chromatograms as strategies for data interpretation. From these assessments we found that differential genes from a single replicate were confirmed by other replicates on the same instrument from 61-93% of the time. When comparing across different instruments and quantitative technologies, differential genes were reproduced by other data sets from 67-99% of the time. Projecting gene differences to biological pathways and networks increased the similarities. These overlaps send an encouraging message about the maturity of technologies for proteomic differentiation.« less

A lattice Boltzmann model with an amending function for simulating nonlinear partial differential equations

NASA Astrophysics Data System (ADS)

Chen, Lin-Jie; Ma, Chang-Feng

2010-01-01

This paper proposes a lattice Boltzmann model with an amending function for one-dimensional nonlinear partial differential equations (NPDEs) in the form ut + αuux + βunux + γuxx + δuxxx + ζuxxxx = 0. This model is different from existing models because it lets the time step be equivalent to the square of the space step and derives higher accuracy and nonlinear terms in NPDEs. With the Chapman-Enskog expansion, the governing evolution equation is recovered correctly from the continuous Boltzmann equation. The numerical results agree well with the analytical solutions.

ESTIMATION OF CONSTANT AND TIME-VARYING DYNAMIC PARAMETERS OF HIV INFECTION IN A NONLINEAR DIFFERENTIAL EQUATION MODEL.

PubMed

Liang, Hua; Miao, Hongyu; Wu, Hulin

2010-03-01

Modeling viral dynamics in HIV/AIDS studies has resulted in deep understanding of pathogenesis of HIV infection from which novel antiviral treatment guidance and strategies have been derived. Viral dynamics models based on nonlinear differential equations have been proposed and well developed over the past few decades. However, it is quite challenging to use experimental or clinical data to estimate the unknown parameters (both constant and time-varying parameters) in complex nonlinear differential equation models. Therefore, investigators usually fix some parameter values, from the literature or by experience, to obtain only parameter estimates of interest from clinical or experimental data. However, when such prior information is not available, it is desirable to determine all the parameter estimates from data. In this paper, we intend to combine the newly developed approaches, a multi-stage smoothing-based (MSSB) method and the spline-enhanced nonlinear least squares (SNLS) approach, to estimate all HIV viral dynamic parameters in a nonlinear differential equation model. In particular, to the best of our knowledge, this is the first attempt to propose a comparatively thorough procedure, accounting for both efficiency and accuracy, to rigorously estimate all key kinetic parameters in a nonlinear differential equation model of HIV dynamics from clinical data. These parameters include the proliferation rate and death rate of uninfected HIV-targeted cells, the average number of virions produced by an infected cell, and the infection rate which is related to the antiviral treatment effect and is time-varying. To validate the estimation methods, we verified the identifiability of the HIV viral dynamic model and performed simulation studies. We applied the proposed techniques to estimate the key HIV viral dynamic parameters for two individual AIDS patients treated with antiretroviral therapies. We demonstrate that HIV viral dynamics can be well characterized and

Comment on "Applications of homogenous balanced principle on investigating exact solutions to a series of time fractional nonlinear PDEs", [Commun Nonlinear Sci Numer Simulat 47 (2017) 253-266

NASA Astrophysics Data System (ADS)

Li, Xiangzheng

2018-06-01

A counterexample is given to show that the product rule of the Caputo fractional derivatives does not hold except on a special point. The function-expansion method of separation variable proposed by Rui[Commun Nonlinear Sci Numer Simulat 47 (2017) 253-266] based on the product rule must be modified.

Remarks on the Non-Linear Differential Equation the Second Derivative of Theta Plus A Sine Theta Equals 0.

ERIC Educational Resources Information Center

Fay, Temple H.; O'Neal, Elizabeth A.

1985-01-01

The authors draw together a variety of facts concerning a nonlinear differential equation and compare the exact solution with approximate solutions. Then they provide an expository introduction to the elliptic sine function suitable for presentation in undergraduate courses on differential equations. (MNS)

Systems-based decomposition schemes for the approximate solution of multi-term fractional differential equations

NASA Astrophysics Data System (ADS)

Ford, Neville J.; Connolly, Joseph A.

2009-07-01

We give a comparison of the efficiency of three alternative decomposition schemes for the approximate solution of multi-term fractional differential equations using the Caputo form of the fractional derivative. The schemes we compare are based on conversion of the original problem into a system of equations. We review alternative approaches and consider how the most appropriate numerical scheme may be chosen to solve a particular equation.

A new multi-step technique with differential transform method for analytical solution of some nonlinear variable delay differential equations.

PubMed

Benhammouda, Brahim; Vazquez-Leal, Hector

2016-01-01

This work presents an analytical solution of some nonlinear delay differential equations (DDEs) with variable delays. Such DDEs are difficult to treat numerically and cannot be solved by existing general purpose codes. A new method of steps combined with the differential transform method (DTM) is proposed as a powerful tool to solve these DDEs. This method reduces the DDEs to ordinary differential equations that are then solved by the DTM. Furthermore, we show that the solutions can be improved by Laplace-Padé resummation method. Two examples are presented to show the efficiency of the proposed technique. The main advantage of this technique is that it possesses a simple procedure based on a few straight forward steps and can be combined with any analytical method, other than the DTM, like the homotopy perturbation method.

The existence of periodic solutions for nonlinear beam equations on Td by a para-differential method

NASA Astrophysics Data System (ADS)

Chen, Bochao; Li, Yong; Gao, Yixian

2018-05-01

This paper focuses on the construction of periodic solutions of nonlinear beam equations on the $d$-dimensional tori. For a large set of frequencies, we demonstrate that an equivalent form of the nonlinear equations can be obtained by a para-differential conjugation. Given the non-resonant conditions on each finite dimensional subspaces, it is shown that the periodic solutions can be constructed for the block diagonal equation by a classical iteration scheme.

Nonlinear dynamics analysis of a low-temperature-differential kinematic Stirling heat engine

NASA Astrophysics Data System (ADS)

Izumida, Yuki

2018-03-01

The low-temperature-differential (LTD) Stirling heat engine technology constitutes one of the important sustainable energy technologies. The basic question of how the rotational motion of the LTD Stirling heat engine is maintained or lost based on the temperature difference is thus a practically and physically important problem that needs to be clearly understood. Here, we approach this problem by proposing and investigating a minimal nonlinear dynamic model of an LTD kinematic Stirling heat engine. Our model is described as a driven nonlinear pendulum where the motive force is the temperature difference. The rotational state and the stationary state of the engine are described as a stable limit cycle and a stable fixed point of the dynamical equations, respectively. These two states coexist under a sufficient temperature difference, whereas the stable limit cycle does not exist under a temperature difference that is too small. Using a nonlinear bifurcation analysis, we show that the disappearance of the stable limit cycle occurs via a homoclinic bifurcation, with the temperature difference being the bifurcation parameter.

Series expansion solutions for the multi-term time and space fractional partial differential equations in two- and three-dimensions

NASA Astrophysics Data System (ADS)

Ye, H.; Liu, F.; Turner, I.; Anh, V.; Burrage, K.

2013-09-01

Fractional partial differential equations with more than one fractional derivative in time describe some important physical phenomena, such as the telegraph equation, the power law wave equation, or the Szabo wave equation. In this paper, we consider two- and three-dimensional multi-term time and space fractional partial differential equations. The multi-term time-fractional derivative is defined in the Caputo sense, whose order belongs to the interval (1,2],(2,3],(3,4] or (0, m], and the space-fractional derivative is referred to as the fractional Laplacian form. We derive series expansion solutions based on a spectral representation of the Laplacian operator on a bounded region. Some applications are given for the two- and three-dimensional telegraph equation, power law wave equation and Szabo wave equation.

Study of coupled nonlinear partial differential equations for finding exact analytical solutions.

PubMed

Khan, Kamruzzaman; Akbar, M Ali; Koppelaar, H

2015-07-01

Exact solutions of nonlinear partial differential equations (NPDEs) are obtained via the enhanced (G'/G)-expansion method. The method is subsequently applied to find exact solutions of the Drinfel'd-Sokolov-Wilson (DSW) equation and the (2+1)-dimensional Painlevé integrable Burgers (PIB) equation. The efficiency of this method for finding these exact solutions is demonstrated. The method is effective and applicable for many other NPDEs in mathematical physics.

A Differential Evolution Algorithm Based on Nikaido-Isoda Function for Solving Nash Equilibrium in Nonlinear Continuous Games

PubMed Central

He, Feng; Zhang, Wei; Zhang, Guoqiang

2016-01-01

A differential evolution algorithm for solving Nash equilibrium in nonlinear continuous games is presented in this paper, called NIDE (Nikaido-Isoda differential evolution). At each generation, parent and child strategy profiles are compared one by one pairwisely, adapting Nikaido-Isoda function as fitness function. In practice, the NE of nonlinear game model with cubic cost function and quadratic demand function is solved, and this method could also be applied to non-concave payoff functions. Moreover, the NIDE is compared with the existing Nash Domination Evolutionary Multiplayer Optimization (NDEMO), the result showed that NIDE was significantly better than NDEMO with less iterations and shorter running time. These numerical examples suggested that the NIDE method is potentially useful. PMID:27589229

Adaptive neural output-feedback control for nonstrict-feedback time-delay fractional-order systems with output constraints and actuator nonlinearities.

PubMed

Zouari, Farouk; Ibeas, Asier; Boulkroune, Abdesselem; Cao, Jinde; Mehdi Arefi, Mohammad

2018-06-01

This study addresses the issue of the adaptive output tracking control for a category of uncertain nonstrict-feedback delayed incommensurate fractional-order systems in the presence of nonaffine structures, unmeasured pseudo-states, unknown control directions, unknown actuator nonlinearities and output constraints. Firstly, the mean value theorem and the Gaussian error function are introduced to eliminate the difficulties that arise from the nonaffine structures and the unknown actuator nonlinearities, respectively. Secondly, the immeasurable tracking error variables are suitably estimated by constructing a fractional-order linear observer. Thirdly, the neural network, the Razumikhin Lemma, the variable separation approach, and the smooth Nussbaum-type function are used to deal with the uncertain nonlinear dynamics, the unknown time-varying delays, the nonstrict feedback and the unknown control directions, respectively. Fourthly, asymmetric barrier Lyapunov functions are employed to overcome the violation of the output constraints and to tune online the parameters of the adaptive neural controller. Through rigorous analysis, it is proved that the boundedness of all variables in the closed-loop system and the semi global asymptotic tracking are ensured without transgression of the constraints. The principal contributions of this study can be summarized as follows: (1) based on Caputo's definitions and new lemmas, methods concerning the controllability, observability and stability analysis of integer-order systems are extended to fractional-order ones, (2) the output tracking objective for a relatively large class of uncertain systems is achieved with a simple controller and less tuning parameters. Finally, computer-simulation studies from the robotic field are given to demonstrate the effectiveness of the proposed controller. Copyright © 2018 Elsevier Ltd. All rights reserved.

A class of fractional differential hemivariational inequalities with application to contact problem

NASA Astrophysics Data System (ADS)

Zeng, Shengda; Liu, Zhenhai; Migorski, Stanislaw

2018-04-01

In this paper, we study a class of generalized differential hemivariational inequalities of parabolic type involving the time fractional order derivative operator in Banach spaces. We use the Rothe method combined with surjectivity of multivalued pseudomonotone operators and properties of the Clarke generalized gradient to establish existence of solution to the abstract inequality. As an illustrative application, a frictional quasistatic contact problem for viscoelastic materials with adhesion is investigated, in which the friction and contact conditions are described by the Clarke generalized gradient of nonconvex and nonsmooth functionals, and the constitutive relation is modeled by the fractional Kelvin-Voigt law.

Variational Solutions and Random Dynamical Systems to SPDEs Perturbed by Fractional Gaussian Noise

PubMed Central

Zeng, Caibin; Yang, Qigui; Cao, Junfei

2014-01-01

This paper deals with the following type of stochastic partial differential equations (SPDEs) perturbed by an infinite dimensional fractional Brownian motion with a suitable volatility coefficient Φ: dX(t) = A(X(t))dt+Φ(t)dB H(t), where A is a nonlinear operator satisfying some monotonicity conditions. Using the variational approach, we prove the existence and uniqueness of variational solutions to such system. Moreover, we prove that this variational solution generates a random dynamical system. The main results are applied to a general type of nonlinear SPDEs and the stochastic generalized p-Laplacian equation. PMID:24574903

A comparison/validation of a fractional derivative model with an empirical model of non-linear shock waves in swelling shales

NASA Astrophysics Data System (ADS)

Droghei, Riccardo; Salusti, Ettore

2013-04-01

Control of drilling parameters, as fluid pressure, mud weight, salt concentration is essential to avoid instabilities when drilling through shale sections. To investigate shale deformation, fundamental for deep oil drilling and hydraulic fracturing for gas extraction ("fracking"), a non-linear model of mechanic and chemo-poroelastic interactions among fluid, solute and the solid matrix is here discussed. The two equations of this model describe the isothermal evolution of fluid pressure and solute density in a fluid saturated porous rock. Their solutions are quick non-linear Burger's solitary waves, potentially destructive for deep operations. In such analysis the effect of diffusion, that can play a particular role in fracking, is investigated. Then, following Civan (1998), both diffusive and shock waves are applied to fine particles filtration due to such quick transients , their effect on the adjacent rocks and the resulting time-delayed evolution. Notice how time delays in simple porous media dynamics have recently been analyzed using a fractional derivative approach. To make a tentative comparison of these two deeply different methods,in our model we insert fractional time derivatives, i.e. a kind of time-average of the fluid-rocks interactions. Then the delaying effects of fine particles filtration is compared with fractional model time delays. All this can be seen as an empirical check of these fractional models.

Penalized Nonlinear Least Squares Estimation of Time-Varying Parameters in Ordinary Differential Equations

PubMed Central

Cao, Jiguo; Huang, Jianhua Z.; Wu, Hulin

2012-01-01

Ordinary differential equations (ODEs) are widely used in biomedical research and other scientific areas to model complex dynamic systems. It is an important statistical problem to estimate parameters in ODEs from noisy observations. In this article we propose a method for estimating the time-varying coefficients in an ODE. Our method is a variation of the nonlinear least squares where penalized splines are used to model the functional parameters and the ODE solutions are approximated also using splines. We resort to the implicit function theorem to deal with the nonlinear least squares objective function that is only defined implicitly. The proposed penalized nonlinear least squares method is applied to estimate a HIV dynamic model from a real dataset. Monte Carlo simulations show that the new method can provide much more accurate estimates of functional parameters than the existing two-step local polynomial method which relies on estimation of the derivatives of the state function. Supplemental materials for the article are available online. PMID:23155351

Negative Differential Conductance in Polyporphyrin Oligomers with Nonlinear Backbones.

PubMed

Kuang, Guowen; Chen, Shi Zhang; Yan, Linghao; Chen, Ke Qiu; Shang, Xuesong; Liu, Pei Nian; Lin, Nian

2018-01-17

We study negative differential conductance (NDC) effects in polyporphyrin oligomers with nonlinear backbones. Using a low-temperature scanning tunneling microscope, we selectively controlled the charge transport path in single oligomer wires. We observed robust NDC when charge passed through a T-shape junction, bistable NDC when charge passed through a 90° kink and no NDC when charge passed through a 120° kink. Aided by density functional theory with nonequilibrium Green's functions simulations, we attributed this backbone-dependent NDC to bias-modulated hybridization of the electrode states with the resonant transport molecular orbital. We argue this mechanism is generic in molecular systems, which opens a new route of designing molecular NDC devices.

Numerical method for solution of systems of non-stationary spatially one-dimensional nonlinear differential equations

NASA Technical Reports Server (NTRS)

Morozov, S. K.; Krasitskiy, O. P.

1978-01-01

A computational scheme and a standard program is proposed for solving systems of nonstationary spatially one-dimensional nonlinear differential equations using Newton's method. The proposed scheme is universal in its applicability and its reduces to a minimum the work of programming. The program is written in the FORTRAN language and can be used without change on electronic computers of type YeS and BESM-6. The standard program described permits the identification of nonstationary (or stationary) solutions to systems of spatially one-dimensional nonlinear (or linear) partial differential equations. The proposed method may be used to solve a series of geophysical problems which take chemical reactions, diffusion, and heat conductivity into account, to evaluate nonstationary thermal fields in two-dimensional structures when in one of the geometrical directions it can take a small number of discrete levels, and to solve problems in nonstationary gas dynamics.

Nonlinear zero-sum differential game analysis by singular perturbation methods

NASA Technical Reports Server (NTRS)

Sinar, J.; Farber, N.

1982-01-01

A class of nonlinear, zero-sum differential games, exhibiting time-scale separation properties, can be analyzed by singular-perturbation techniques. The merits of such an analysis, leading to an approximate game solution, as well as the 'well-posedness' of the formulation, are discussed. This approach is shown to be attractive for investigating pursuit-evasion problems; the original multidimensional differential game is decomposed to a 'simple pursuit' (free-stream) game and two independent (boundary-layer) optimal-control problems. Using multiple time-scale boundary-layer models results in a pair of uniformly valid zero-order composite feedback strategies. The dependence of suboptimal strategies on relative geometry and own-state measurements is demonstrated by a three dimensional, constant-speed example. For game analysis with realistic vehicle dynamics, the technique of forced singular perturbations and a variable modeling approach is proposed. Accuracy of the analysis is evaluated by comparison with the numerical solution of a time-optimal, variable-speed 'game of two cars' in the horizontal plane.

Cubication of Conservative Nonlinear Oscillators

ERIC Educational Resources Information Center

Belendez, Augusto; Alvarez, Mariela L.; Fernandez, Elena; Pascual, Immaculada

2009-01-01

A cubication procedure of the nonlinear differential equation for conservative nonlinear oscillators is analysed and discussed. This scheme is based on the Chebyshev series expansion of the restoring force, and this allows us to approximate the original nonlinear differential equation by a Duffing equation in which the coefficients for the linear…

Study of coupled nonlinear partial differential equations for finding exact analytical solutions

PubMed Central

Khan, Kamruzzaman; Akbar, M. Ali; Koppelaar, H.

2015-01-01

Exact solutions of nonlinear partial differential equations (NPDEs) are obtained via the enhanced (G′/G)-expansion method. The method is subsequently applied to find exact solutions of the Drinfel'd–Sokolov–Wilson (DSW) equation and the (2+1)-dimensional Painlevé integrable Burgers (PIB) equation. The efficiency of this method for finding these exact solutions is demonstrated. The method is effective and applicable for many other NPDEs in mathematical physics. PMID:26587256

Design of an all-optical fractional-order differentiator with terahertz bandwidth based on a fiber Bragg grating in transmission.

PubMed

Liu, Xin; Shu, Xuewen

2017-08-20

All-optical fractional-order temporal differentiators with bandwidths reaching terahertz (THz) values are demonstrated with transmissive fiber Bragg gratings. Since the designed fractional-order differentiator is a minimum phase function, the reflective phase of the designed function can be chosen arbitrarily. As examples, we first design several 0.5th-order differentiators with bandwidths reaching the THz range for comparison. The reflective phases of the 0.5th-order differentiators are chosen to be linear phase, quadratic phase, cubic phase, and biquadratic phase, respectively. We find that both the maximum coupling coefficient and the spatial resolution of the designed grating increase when the reflective phase varies from quadratic function to cubic function to biquadratic function. Furthermore, when the reflective phase is chosen to be a quadratic function, the obtained grating coupling coefficient and period are more likely to be achieved in practice. Then we design fractional-order differentiators with different orders when the reflective phase is chosen to be a quadratic function. We see that when the designed order of the differentiator increases, the obtained maximum coupling coefficient also increases while the oscillation of the coupling coefficient decreases. Finally, we give the numerical performance of the designed 0.5th-order differentiator by showing its temporal response and calculating its cross-correlation coefficient.

Stability Analysis of Distributed Order Fractional Chen System

PubMed Central

Aminikhah, H.; Refahi Sheikhani, A.; Rezazadeh, H.

2013-01-01

We first investigate sufficient and necessary conditions of stability of nonlinear distributed order fractional system and then we generalize the integer-order Chen system into the distributed order fractional domain. Based on the asymptotic stability theory of nonlinear distributed order fractional systems, the stability of distributed order fractional Chen system is discussed. In addition, we have found that chaos exists in the double fractional order Chen system. Numerical solutions are used to verify the analytical results. PMID:24489508

State-dependent differential Riccati equation to track control of time-varying systems with state and control nonlinearities.

PubMed

Korayem, M H; Nekoo, S R

2015-07-01

This work studies an optimal control problem using the state-dependent Riccati equation (SDRE) in differential form to track for time-varying systems with state and control nonlinearities. The trajectory tracking structure provides two nonlinear differential equations: the state-dependent differential Riccati equation (SDDRE) and the feed-forward differential equation. The independence of the governing equations and stability of the controller are proven along the trajectory using the Lyapunov approach. Backward integration (BI) is capable of solving the equations as a numerical solution; however, the forward solution methods require the closed-form solution to fulfill the task. A closed-form solution is introduced for SDDRE, but the feed-forward differential equation has not yet been obtained. Different ways of solving the problem are expressed and analyzed. These include BI, closed-form solution with corrective assumption, approximate solution, and forward integration. Application of the tracking problem is investigated to control robotic manipulators possessing rigid or flexible joints. The intention is to release a general program for automatic implementation of an SDDRE controller for any manipulator that obeys the Denavit-Hartenberg (D-H) principle when only D-H parameters are received as input data. Copyright © 2015 ISA. Published by Elsevier Ltd. All rights reserved.

Zinc isotope fractionation during magmatic differentiation and the isotopic composition of the bulk Earth

USGS Publications Warehouse

Chen, Heng; Savage, Paul S.; Teng, Fang-Zehn; Helz, Rosalind T.; Moynier, Frédéric

2013-01-01

he zinc stable isotope system has been successfully applied to many and varied fields in geochemistry, but to date it is still not completely clear how this isotope system is affected by igneous processes. In order to evaluate the potential application of Zn isotopes as a proxy for planetary differentiation and volatile history, it is important to constrain the magnitude of Zn isotopic fractionation induced by magmatic differentiation. In this study we present high-precision Zn isotope analyses of two sets of chemically diverse, cogenetic samples from Kilauea Iki lava lake, Hawaii, and Hekla volcano, Iceland, which both show clear evidence of having undergone variable and significant degrees of magmatic differentiation. The Kilauea Iki samples display small but resolvable variations in Zn isotope composition (0.26‰66Zn66Zn defined as the per mille deviation of a sample's 66Zn/64Zn compositional ratio from the JMC-Lyon standard), with the most differentiated lithologies exhibiting more positive δ66Zn values. This fractionation is likely a result of the crystallization of olivine and/or Fe–Ti oxides, which can both host Zn in their crystal structures. Samples from Hekla have a similar range of isotopic variation (0.22‰66Zn66Zn=0.28±0.05‰ (2s.d.).

Unveiling hidden properties of young star clusters: differential reddening, star-formation spread, and binary fraction

NASA Astrophysics Data System (ADS)

Bonatto, C.; Lima, E. F.; Bica, E.

2012-04-01

Context. Usually, important parameters of young, low-mass star clusters are very difficult to obtain by means of photometry, especially when differential reddening and/or binaries occur in large amounts. Aims: We present a semi-analytical approach (ASAmin) that, when applied to the Hess diagram of a young star cluster, is able to retrieve the values of mass, age, star-formation spread, distance modulus, foreground and differential reddening, and binary fraction. Methods: The global optimisation method known as adaptive simulated annealing (ASA) is used to minimise the residuals between the observed and simulated Hess diagrams of a star cluster. The simulations are realistic and take the most relevant parameters of young clusters into account. Important features of the simulations are a normal (Gaussian) differential reddening distribution, a time-decreasing star-formation rate, the unresolved binaries, and the smearing effect produced by photometric uncertainties on Hess diagrams. Free parameters are cluster mass, age, distance modulus, star-formation spread, foreground and differential reddening, and binary fraction. Results: Tests with model clusters built with parameters spanning a broad range of values show that ASAmin retrieves the input values with a high precision for cluster mass, distance modulus, and foreground reddening, but they are somewhat lower for the remaining parameters. Given the statistical nature of the simulations, several runs should be performed to obtain significant convergence patterns. Specifically, we find that the retrieved (absolute minimum) parameters converge to mean values with a low dispersion as the Hess residuals decrease. When applied to actual young clusters, the retrieved parameters follow convergence patterns similar to the models. We show how the stochasticity associated with the early phases may affect the results, especially in low-mass clusters. This effect can be minimised by averaging out several twin clusters in the

Dynamic assessment of nonlinear typical section aeroviscoelastic systems using fractional derivative-based viscoelastic model

NASA Astrophysics Data System (ADS)

Sales, T. P.; Marques, Flávio D.; Pereira, Daniel A.; Rade, Domingos A.

2018-06-01

Nonlinear aeroelastic systems are prone to the appearance of limit cycle oscillations, bifurcations, and chaos. Such problems are of increasing concern in aircraft design since there is the need to control nonlinear instabilities and improve safety margins, at the same time as aircraft are subjected to increasingly critical operational conditions. On the other hand, in spite of the fact that viscoelastic materials have already been successfully used for the attenuation of undesired vibrations in several types of mechanical systems, a small number of research works have addressed the feasibility of exploring the viscoelastic effect to improve the behavior of nonlinear aeroelastic systems. In this context, the objective of this work is to assess the influence of viscoelastic materials on the aeroelastic features of a three-degrees-of-freedom typical section with hardening structural nonlinearities. The equations of motion are derived accounting for the presence of viscoelastic materials introduced in the resilient elements associated to each degree-of-freedom. A constitutive law based on fractional derivatives is adopted, which allows the modeling of temperature-dependent viscoelastic behavior in time and frequency domains. The unsteady aerodynamic loading is calculated based on the classical linear potential theory for arbitrary airfoil motion. The aeroelastic behavior is investigated through time domain simulations, and subsequent frequency transformations, from which bifurcations are identified from diagrams of limit cycle oscillations amplitudes versus airspeed. The influence of the viscoelastic effect on the aeroelastic behavior, for different values of temperature, is also investigated. The numerical simulations show that viscoelastic damping can increase the flutter speed and reduce the amplitudes of limit cycle oscillations. These results prove the potential that viscoelastic materials have to increase aircraft components safety margins regarding aeroelastic

Boundedness of the solutions for certain classes of fractional differential equations with application to adaptive systems.

PubMed

Aguila-Camacho, Norelys; Duarte-Mermoud, Manuel A

2016-01-01

This paper presents the analysis of three classes of fractional differential equations appearing in the field of fractional adaptive systems, for the case when the fractional order is in the interval α ∈(0,1] and the Caputo definition for fractional derivatives is used. The boundedness of the solutions is proved for all three cases, and the convergence to zero of the mean value of one of the variables is also proved. Applications of the obtained results to fractional adaptive schemes in the context of identification and control problems are presented at the end of the paper, including numerical simulations which support the analytical results. Copyright © 2015 ISA. Published by Elsevier Ltd. All rights reserved.

Multilevel Modeling of Two Cyclical Processes: Extending Differential Structural Equation Modeling to Nonlinear Coupled Systems

ERIC Educational Resources Information Center

Butner, Jonathan; Amazeen, Polemnia G.; Mulvey, Genna M.

2005-01-01

The authors present a dynamical multilevel model that captures changes over time in the bidirectional, potentially asymmetric influence of 2 cyclical processes. S. M. Boker and J. Graham's (1998) differential structural equation modeling approach was expanded to the case of a nonlinear coupled oscillator that is common in bimanual coordination…

Accessible solitons of fractional dimension

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zhong, Wei-Ping, E-mail: zhongwp6@126.com; Texas A&M University at Qatar, P.O. Box 23874, Doha; Belić, Milivoj

We demonstrate that accessible solitons described by an extended Schrödinger equation with the Laplacian of fractional dimension can exist in strongly nonlocal nonlinear media. The soliton solutions of the model are constructed by two special functions, the associated Legendre polynomials and the Laguerre polynomials in the fraction-dimensional space. Our results show that these fractional accessible solitons form a soliton family which includes crescent solitons, and asymmetric single-layer and multi-layer necklace solitons. -- Highlights: •Analytic solutions of a fractional Schrödinger equation are obtained. •The solutions are produced by means of self-similar method applied to the fractional Schrödinger equation with parabolic potential.more » •The fractional accessible solitons form crescent, asymmetric single-layer and multilayer necklace profiles. •The model applies to the propagation of optical pulses in strongly nonlocal nonlinear media.« less

Paradox of enrichment: A fractional differential approach with memory

NASA Astrophysics Data System (ADS)

Rana, Sourav; Bhattacharya, Sabyasachi; Pal, Joydeep; N'Guérékata, Gaston M.; Chattopadhyay, Joydev

2013-09-01

The paradox of enrichment (PoE) proposed by Rosenzweig [M. Rosenzweig, The paradox of enrichment, Science 171 (1971) 385-387] is still a fundamental problem in ecology. Most of the solutions have been proposed at an individual species level of organization and solutions at community level are lacking. Knowledge of how learning and memory modify behavioral responses to species is a key factor in making a crucial link between species and community levels. PoE resolution via these two organizational levels can be interpreted as a microscopic- and macroscopic-level solution. Fractional derivatives provide an excellent tool for describing this memory and the hereditary properties of various materials and processes. The derivatives can be physically interpreted via two time scales that are considered simultaneously: the ideal, equably flowing homogeneous local time, and the cosmic (inhomogeneous) non-local time. Several mechanisms and theories have been proposed to resolve the PoE problem, but a universally accepted theory is still lacking because most studies have focused on local effects and ignored non-local effects, which capture memory. Here we formulate the fractional counterpart of the Rosenzweig model and analyze the stability behavior of a system. We conclude that there is a threshold for the memory effect parameter beyond which the Rosenzweig model is stable and may be used as a potential agent to resolve PoE from a new perspective via fractional differential equations.

An approximation theory for nonlinear partial differential equations with applications to identification and control

NASA Technical Reports Server (NTRS)

Banks, H. T.; Kunisch, K.

1982-01-01

Approximation results from linear semigroup theory are used to develop a general framework for convergence of approximation schemes in parameter estimation and optimal control problems for nonlinear partial differential equations. These ideas are used to establish theoretical convergence results for parameter identification using modal (eigenfunction) approximation techniques. Results from numerical investigations of these schemes for both hyperbolic and parabolic systems are given.

A pertinent approach to solve nonlinear fuzzy integro-differential equations.

PubMed

Narayanamoorthy, S; Sathiyapriya, S P

2016-01-01

Fuzzy integro-differential equations is one of the important parts of fuzzy analysis theory that holds theoretical as well as applicable values in analytical dynamics and so an appropriate computational algorithm to solve them is in essence. In this article, we use parametric forms of fuzzy numbers and suggest an applicable approach for solving nonlinear fuzzy integro-differential equations using homotopy perturbation method. A clear and detailed description of the proposed method is provided. Our main objective is to illustrate that the construction of appropriate convex homotopy in a proper way leads to highly accurate solutions with less computational work. The efficiency of the approximation technique is expressed via stability and convergence analysis so as to guarantee the efficiency and performance of the methodology. Numerical examples are demonstrated to verify the convergence and it reveals the validity of the presented numerical technique. Numerical results are tabulated and examined by comparing the obtained approximate solutions with the known exact solutions. Graphical representations of the exact and acquired approximate fuzzy solutions clarify the accuracy of the approach.

Measurements of the S-wave fraction in B 0 → K + π - μ + μ - decays and the B 0 → K ∗(892)0 μ + μ - differential branching fraction

NASA Astrophysics Data System (ADS)

Aaij, R.; Adeva, B.; Adinolfi, M.; Ajaltouni, Z.; Akar, S.; Albrecht, J.; Alessio, F.; Alexander, M.; Ali, S.; Alkhazov, G.; Alvarez Cartelle, P.; Alves, A. A.; Amato, S.; Amerio, S.; Amhis, Y.; An, L.; Anderlini, L.; Andreassi, G.; Andreotti, M.; Andrews, J. E.; Appleby, R. B.; Aquines Gutierrez, O.; Archilli, F.; d'Argent, P.; Artamonov, A.; Artuso, M.; Aslanides, E.; Auriemma, G.; Baalouch, M.; Bachmann, S.; Back, J. J.; Badalov, A.; Baesso, C.; Baldini, W.; Barlow, R. J.; Barschel, C.; Barsuk, S.; Barter, W.; Batozskaya, V.; Battista, V.; Bay, A.; Beaucourt, L.; Beddow, J.; Bedeschi, F.; Bediaga, I.; Bel, L. J.; Bellee, V.; Belloli, N.; Belous, K.; Belyaev, I.; Ben-Haim, E.; Bencivenni, G.; Benson, S.; Benton, J.; Berezhnoy, A.; Bernet, R.; Bertolin, A.; Bettler, M.-O.; van Beuzekom, M.; Bifani, S.; Billoir, P.; Bird, T.; Birnkraut, A.; Bitadze, A.; Bizzeti, A.; Blake, T.; Blanc, F.; Blouw, J.; Blusk, S.; Bocci, V.; Boettcher, T.; Bondar, A.; Bondar, N.; Bonivento, W.; Borghi, S.; Borisyak, M.; Borsato, M.; Bossu, F.; Boubdir, M.; Bowcock, T. J. V.; Bowen, E.; Bozzi, C.; Braun, S.; Britsch, M.; Britton, T.; Brodzicka, J.; Buchanan, E.; Burr, C.; Bursche, A.; Buytaert, J.; Cadeddu, S.; Calabrese, R.; Calvi, M.; Calvo Gomez, M.; Campana, P.; Campora Perez, D.; Capriotti, L.; Carbone, A.; Carboni, G.; Cardinale, R.; Cardini, A.; Carniti, P.; Carson, L.; Carvalho Akiba, K.; Casse, G.; Cassina, L.; Castillo Garcia, L.; Cattaneo, M.; Cauet, Ch.; Cavallero, G.; Cenci, R.; Charles, M.; Charpentier, Ph.; Chatzikonstantinidis, G.; Chefdeville, M.; Chen, S.; Cheung, S.-F.; Chobanova, V.; Chrzaszcz, M.; Cid Vidal, X.; Ciezarek, G.; Clarke, P. E. L.; Clemencic, M.; Cliff, H. V.; Closier, J.; Coco, V.; Cogan, J.; Cogneras, E.; Cogoni, V.; Cojocariu, L.; Collazuol, G.; Collins, P.; Comerma-Montells, A.; Contu, A.; Cook, A.; Coquereau, S.; Corti, G.; Corvo, M.; Couturier, B.; Cowan, G. A.; Craik, D. C.; Crocombe, A.; Cruz Torres, M.; Cunliffe, S.; Currie, R.; D'Ambrosio, C.; Dall'Occo, E.; Dalseno, J.; David, P. N. Y.; Davis, A.; De Aguiar Francisco, O.; De Bruyn, K.; De Capua, S.; De Cian, M.; De Miranda, J. M.; De Paula, L.; De Simone, P.; Dean, C.-T.; Decamp, D.; Deckenhoff, M.; Del Buono, L.; Demmer, M.; Derkach, D.; Deschamps, O.; Dettori, F.; Dey, B.; Di Canto, A.; Dijkstra, H.; Dordei, F.; Dorigo, M.; Dosil Suárez, A.; Dovbnya, A.; Dreimanis, K.; Dufour, L.; Dujany, G.; Dungs, K.; Durante, P.; Dzhelyadin, R.; Dziurda, A.; Dzyuba, A.; Déléage, N.; Easo, S.; Egede, U.; Egorychev, V.; Eidelman, S.; Eisenhardt, S.; Eitschberger, U.; Ekelhof, R.; Eklund, L.; Elsasser, Ch.; Ely, S.; Esen, S.; Evans, H. M.; Evans, T.; Falabella, A.; Farley, N.; Farry, S.; Fay, R.; Ferguson, D.; Fernandez Albor, V.; Ferrari, F.; Ferreira Rodrigues, F.; Ferro-Luzzi, M.; Filippov, S.; Fiore, M.; Fiorini, M.; Firlej, M.; Fitzpatrick, C.; Fiutowski, T.; Fleuret, F.; Fohl, K.; Fontana, M.; Fontanelli, F.; Forshaw, D. C.; Forty, R.; Frank, M.; Frei, C.; Frosini, M.; Fu, J.; Furfaro, E.; Färber, C.; Gallas Torreira, A.; Galli, D.; Gallorini, S.; Gambetta, S.; Gandelman, M.; Gandini, P.; Gao, Y.; García Pardiñas, J.; Garra Tico, J.; Garrido, L.; Garsed, P. J.; Gascon, D.; Gaspar, C.; Gavardi, L.; Gazzoni, G.; Gerick, D.; Gersabeck, E.; Gersabeck, M.; Gershon, T.; Ghez, Ph.; Gianì, S.; Gibson, V.; Girard, O. G.; Giubega, L.; Gizdov, K.; Gligorov, V. V.; Golubkov, D.; Golutvin, A.; Gomes, A.; Gorelov, I. V.; Gotti, C.; Grabalosa Gándara, M.; Graciani Diaz, R.; Granado Cardoso, L. A.; Graugés, E.; Graverini, E.; Graziani, G.; Grecu, A.; Griffith, P.; Grillo, L.; Grünberg, O.; Gushchin, E.; Guz, Yu.; Gys, T.; Göbel, C.; Hadavizadeh, T.; Hadjivasiliou, C.; Haefeli, G.; Haen, C.; Haines, S. C.; Hall, S.; Hamilton, B.; Han, X.; Hansmann-Menzemer, S.; Harnew, N.; Harnew, S. T.; Harrison, J.; He, J.; Head, T.; Heister, A.; Hennessy, K.; Henrard, P.; Henry, L.; Hernando Morata, J. A.; van Herwijnen, E.; Heß, M.; Hicheur, A.; Hill, D.; Hombach, C.; Hulsbergen, W.; Humair, T.; Hushchyn, M.; Hussain, N.; Hutchcroft, D.; Idzik, M.; Ilten, P.; Jacobsson, R.; Jaeger, A.; Jalocha, J.; Jans, E.; Jawahery, A.; John, M.; Johnson, D.; Jones, C. R.; Joram, C.; Jost, B.; Jurik, N.; Kandybei, S.; Kanso, W.; Karacson, M.; Karbach, T. M.; Karodia, S.; Kecke, M.; Kelsey, M.; Kenyon, I. R.; Kenzie, M.; Ketel, T.; Khairullin, E.; Khanji, B.; Khurewathanakul, C.; Kirn, T.; Klaver, S.; Klimaszewski, K.; Kolpin, M.; Komarov, I.; Koopman, R. F.; Koppenburg, P.; Kozachuk, A.; Kozeiha, M.; Kravchuk, L.; Kreplin, K.; Kreps, M.; Krokovny, P.; Kruse, F.; Krzemien, W.; Kucewicz, W.; Kucharczyk, M.; Kudryavtsev, V.; Kuonen, A. K.; Kurek, K.; Kvaratskheliya, T.; Lacarrere, D.; Lafferty, G.; Lai, A.; Lambert, D.; Lanfranchi, G.; Langenbruch, C.; Langhans, B.; Latham, T.; Lazzeroni, C.; Le Gac, R.; van Leerdam, J.; Lees, J.-P.; Leflat, A.; Lefrançois, J.; Lefèvre, R.; Lemaitre, F.; Lemos Cid, E.; Leroy, O.; Lesiak, T.; Leverington, B.; Li, Y.; Likhomanenko, T.; Lindner, R.; Linn, C.; Lionetto, F.; Liu, B.; Liu, X.; Loh, D.; Longstaff, I.; Lopes, J. H.; Lucchesi, D.; Lucio Martinez, M.; Luo, H.; Lupato, A.; Luppi, E.; Lupton, O.; Lusiani, A.; Lyu, X.; Machefert, F.; Maciuc, F.; Maev, O.; Maguire, K.; Malde, S.; Malinin, A.; Maltsev, T.; Manca, G.; Mancinelli, G.; Manning, P.; Maratas, J.; Marchand, J. F.; Marconi, U.; Marin Benito, C.; Marino, P.; Marks, J.; Martellotti, G.; Martin, M.; Martinelli, M.; Martinez Santos, D.; Martinez Vidal, F.; Martins Tostes, D.; Massacrier, L. M.; Massafferri, A.; Matev, R.; Mathad, A.; Mathe, Z.; Matteuzzi, C.; Mauri, A.; Maurin, B.; Mazurov, A.; McCann, M.; McCarthy, J.; McNab, A.; McNulty, R.; Meadows, B.; Meier, F.; Meissner, M.; Melnychuk, D.; Merk, M.; Michielin, E.; Milanes, D. A.; Minard, M.-N.; Mitzel, D. S.; Molina Rodriguez, J.; Monroy, I. A.; Monteil, S.; Morandin, M.; Morawski, P.; Mordà, A.; Morello, M. J.; Moron, J.; Morris, A. B.; Mountain, R.; Muheim, F.; Mulder, M.; Mussini, M.; Müller, D.; Müller, J.; Müller, K.; Müller, V.; Naik, P.; Nakada, T.; Nandakumar, R.; Nandi, A.; Nasteva, I.; Needham, M.; Neri, N.; Neubert, S.; Neufeld, N.; Neuner, M.; Nguyen, A. D.; Nguyen-Mau, C.; Niess, V.; Nieswand, S.; Niet, R.; Nikitin, N.; Nikodem, T.; Novoselov, A.; O'Hanlon, D. P.; Oblakowska-Mucha, A.; Obraztsov, V.; Ogilvy, S.; Oldeman, R.; Onderwater, C. J. G.; Otalora Goicochea, J. M.; Otto, A.; Owen, P.; Oyanguren, A.; Palano, A.; Palombo, F.; Palutan, M.; Panman, J.; Papanestis, A.; Pappagallo, M.; Pappalardo, L. L.; Pappenheimer, C.; Parker, W.; Parkes, C.; Passaleva, G.; Patel, G. D.; Patel, M.; Patrignani, C.; Pearce, A.; Pellegrino, A.; Penso, G.; Pepe Altarelli, M.; Perazzini, S.; Perret, P.; Pescatore, L.; Petridis, K.; Petrolini, A.; Petrov, A.; Petruzzo, M.; Picatoste Olloqui, E.; Pietrzyk, B.; Pikies, M.; Pinci, D.; Pistone, A.; Piucci, A.; Playfer, S.; Plo Casasus, M.; Poikela, T.; Polci, F.; Poluektov, A.; Polyakov, I.; Polycarpo, E.; Pomery, G. J.; Popov, A.; Popov, D.; Popovici, B.; Potterat, C.; Price, E.; Price, J. D.; Prisciandaro, J.; Pritchard, A.; Prouve, C.; Pugatch, V.; Puig Navarro, A.; Punzi, G.; Qian, W.; Quagliani, R.; Rachwal, B.; Rademacker, J. H.; Rama, M.; Ramos Pernas, M.; Rangel, M. S.; Raniuk, I.; Raven, G.; Redi, F.; Reichert, S.; dos Reis, A. C.; Remon Alepuz, C.; Renaudin, V.; Ricciardi, S.; Richards, S.; Rihl, M.; Rinnert, K.; Rives Molina, V.; Robbe, P.; Rodrigues, A. B.; Rodrigues, E.; Rodriguez Lopez, J. A.; Rodriguez Perez, P.; Rogozhnikov, A.; Roiser, S.; Romanovskiy, V.; Romero Vidal, A.; Ronayne, J. W.; Rotondo, M.; Ruf, T.; Ruiz Valls, P.; Saborido Silva, J. J.; Sagidova, N.; Saitta, B.; Salustino Guimaraes, V.; Sanchez Mayordomo, C.; Sanmartin Sedes, B.; Santacesaria, R.; Santamarina Rios, C.; Santimaria, M.; Santovetti, E.; Sarti, A.; Satriano, C.; Satta, A.; Saunders, D. M.; Savrina, D.; Schael, S.; Schellenberg, M.; Schiller, M.; Schindler, H.; Schlupp, M.; Schmelling, M.; Schmelzer, T.; Schmidt, B.; Schneider, O.; Schopper, A.; Schubert, K.; Schubiger, M.; Schune, M.-H.; Schwemmer, R.; Sciascia, B.; Sciubba, A.; Semennikov, A.; Sergi, A.; Serra, N.; Serrano, J.; Sestini, L.; Seyfert, P.; Shapkin, M.; Shapoval, I.; Shcheglov, Y.; Shears, T.; Shekhtman, L.; Shevchenko, V.; Shires, A.; Siddi, B. G.; Silva Coutinho, R.; Silva de Oliveira, L.; Simi, G.; Sirendi, M.; Skidmore, N.; Skwarnicki, T.; Smith, E.; Smith, I. T.; Smith, J.; Smith, M.; Snoek, H.; Sokoloff, M. D.; Soler, F. J. P.; Souza, D.; Souza De Paula, B.; Spaan, B.; Spradlin, P.; Sridharan, S.; Stagni, F.; Stahl, M.; Stahl, S.; Stefko, P.; Stefkova, S.; Steinkamp, O.; Stenyakin, O.; Stevenson, S.; Stoica, S.; Stone, S.; Storaci, B.; Stracka, S.; Straticiuc, M.; Straumann, U.; Sun, L.; Sutcliffe, W.; Swientek, K.; Syropoulos, V.; Szczekowski, M.; Szumlak, T.; T'Jampens, S.; Tayduganov, A.; Tekampe, T.; Tellarini, G.; Teubert, F.; Thomas, C.; Thomas, E.; van Tilburg, J.; Tisserand, V.; Tobin, M.; Tolk, S.; Tomassetti, L.; Tonelli, D.; Topp-Joergensen, S.; Tournefier, E.; Tourneur, S.; Trabelsi, K.; Traill, M.; Tran, M. T.; Tresch, M.; Trisovic, A.; Tsaregorodtsev, A.; Tsopelas, P.; Tuning, N.; Ukleja, A.; Ustyuzhanin, A.; Uwer, U.; Vacca, C.; Vagnoni, V.; Valat, S.; Valenti, G.; Vallier, A.; Vazquez Gomez, R.; Vazquez Regueiro, P.; Vecchi, S.; van Veghel, M.; Velthuis, J. J.; Veltri, M.; Veneziano, G.; Venkateswaran, A.; Vesterinen, M.; Viaud, B.; Vieira, D.; Vieites Diaz, M.; Vilasis-Cardona, X.; Volkov, V.; Vollhardt, A.; Voneki, B.; Voong, D.; Vorobyev, A.; Vorobyev, V.; Voß, C.; de Vries, J. A.; Vázquez Sierra, C.; Waldi, R.; Wallace, C.; Wallace, R.; Walsh, J.; Wang, J.; Ward, D. R.; Watson, N. K.; Websdale, D.; Weiden, A.; Whitehead, M.; Wicht, J.; Wilkinson, G.; Wilkinson, M.; Williams, M.; Williams, M. P.; Williams, M.; Williams, T.; Wilson, F. F.; Wimberley, J.; Wishahi, J.; Wislicki, W.; Witek, M.; Wormser, G.; Wotton, S. A.; Wraight, K.; Wright, S.; Wyllie, K.; Xie, Y.; Xing, Z.; Xu, Z.; Yang, Z.; Yin, H.; Yu, J.; Yuan, X.; Yushchenko, O.; Zangoli, M.; Zarebski, K. A.; Zavertyaev, M.; Zhang, L.; Zhang, Y.; Zhang, Y.; Zhelezov, A.; Zheng, Y.; Zhokhov, A.; Zhukov, V.; Zucchelli, S.

2016-11-01

A measurement of the differential branching fraction of the decay B 0 → K ∗(892)0 μ + μ - is presented together with a determination of the S-wave fraction of the K + π - system in the decay B 0 → K +π- μ + μ -. The analysis is based on pp-collision data corresponding to an integrated luminosity of 3 fb-1 collected with the LHCb experiment. The measurements are made in bins of the invariant mass squared of the dimuon system, q 2. Precise theoretical predictions for the differential branching fraction of B 0 → K ∗(892)0 μ + μ - decays are available for the q 2 region 1 .1 < q 2 < 6 .0 GeV2 /c 4. In this q 2 region, for the K +π- invariant mass range 796 < m Kπ < 996 MeV /c 2, the S-wave fraction of the K +π- system in B 0 → K +π- μ + μ - decays is found to be {F}S=0.101± 0.017(stat)± 0.009(syst), and the differential branching fraction of B 0 → K ∗(892)0 μ + μ - decays is determined to be dB/d{q}^2=(0.{392}_{-0.019}^{+0.020}(stat)± 0.010(syst)± 0.027(norm))× 1{0}^{-7}{c}^4/{GeV}^2.

Pseudo almost periodic solutions to some systems of nonlinear hyperbolic second-order partial differential equations

NASA Astrophysics Data System (ADS)

Al-Islam, Najja Shakir

In this Dissertation, the existence of pseudo almost periodic solutions to some systems of nonlinear hyperbolic second-order partial differential equations is established. For that, (Al-Islam [4]) is first studied and then obtained under some suitable assumptions. That is, the existence of pseudo almost periodic solutions to a hyperbolic second-order partial differential equation with delay. The second-order partial differential equation (1) represents a mathematical model for the dynamics of gas absorption, given by uxt+a x,tux=Cx,t,u x,t , u0,t=4 t, 1 where a : [0, L] x RR , C : [0, L] x R x RR , and ϕ : RR are (jointly) continuous functions ( t being the greatest integer function) and L > 0. The results in this Dissertation generalize those of Poorkarimi and Wiener [22]. Secondly, a generalization of the above-mentioned system consisting of the non-linear hyperbolic second-order partial differential equation uxt+a x,tux+bx,t ut+cx,tu=f x,t,u, x∈ 0,L,t∈ R, 2 equipped with the boundary conditions ux,0 =40x, u0,t=u 0t, uxx,0=y 0x, x∈0,L, t∈R, 3 where a, b, c : [0, L ] x RR and f : [0, L] x R x RR are (jointly) continuous functions is studied. Under some suitable assumptions, the existence and uniqueness of pseudo almost periodic solutions to particular cases, as well as the general case of the second-order hyperbolic partial differential equation (2) are studied. The results of all studies contained within this text extend those obtained by Aziz and Meyers [6] in the periodic setting.

A Nonlinear differential equation model of Asthma effect of environmental pollution using LHAM

NASA Astrophysics Data System (ADS)

Joseph, G. Arul; Balamuralitharan, S.

2018-04-01

In this paper, we investigated a nonlinear differential equation mathematical model to study the spread of asthma in the environmental pollutants from industry and mainly from tobacco smoke from smokers in different type of population. Smoking is the main cause to spread Asthma in the environment. Numerical simulation is also discussed. Finally by using Liao’s Homotopy analysis Method (LHAM), we found that the approximate analytical solution of Asthmatic disease in the environmental.

Fractional differential equations based modeling of microbial survival and growth curves: model development and experimental validation.

PubMed

Kaur, A; Takhar, P S; Smith, D M; Mann, J E; Brashears, M M

2008-10-01

A fractional differential equations (FDEs)-based theory involving 1- and 2-term equations was developed to predict the nonlinear survival and growth curves of foodborne pathogens. It is interesting to note that the solution of 1-term FDE leads to the Weibull model. Nonlinear regression (Gauss-Newton method) was performed to calculate the parameters of the 1-term and 2-term FDEs. The experimental inactivation data of Salmonella cocktail in ground turkey breast, ground turkey thigh, and pork shoulder; and cocktail of Salmonella, E. coli, and Listeria monocytogenes in ground beef exposed at isothermal cooking conditions of 50 to 66 degrees C were used for validation. To evaluate the performance of 2-term FDE in predicting the growth curves-growth of Salmonella typhimurium, Salmonella Enteritidis, and background flora in ground pork and boneless pork chops; and E. coli O157:H7 in ground beef in the temperature range of 22.2 to 4.4 degrees C were chosen. A program was written in Matlab to predict the model parameters and survival and growth curves. Two-term FDE was more successful in describing the complex shapes of microbial survival and growth curves as compared to the linear and Weibull models. Predicted curves of 2-term FDE had higher magnitudes of R(2) (0.89 to 0.99) and lower magnitudes of root mean square error (0.0182 to 0.5461) for all experimental cases in comparison to the linear and Weibull models. This model was capable of predicting the tails in survival curves, which was not possible using Weibull and linear models. The developed model can be used for other foodborne pathogens in a variety of food products to study the destruction and growth behavior.

Multiple positive solutions to nonlinear boundary value problems of a system for fractional differential equations.

PubMed

Zhai, Chengbo; Hao, Mengru

2014-01-01

By using Krasnoselskii's fixed point theorem, we study the existence of at least one or two positive solutions to a system of fractional boundary value problems given by -D(0+)(ν1)y1(t) = λ1a1(t)f(y1(t), y2(t)), - D(0+)(ν2)y2(t) = λ2a2(t)g(y1(t), y2(t)), where D(0+)(ν) is the standard Riemann-Liouville fractional derivative, ν1, ν2 ∈ (n - 1, n] for n > 3 and n ∈ N, subject to the boundary conditions y1((i))(0) = 0 = y ((i))(0), for 0 ≤ i ≤ n - 2, and [D(0+)(α)y1(t)] t=1 = 0 = [D(0+ (α)y2(t)] t=1, for 1 ≤ α ≤ n - 2, or y1((i))(0) = 0 = y ((i))(0), for 0 ≤ i ≤ n - 2, and [D(0+)(α)y1(t)] t=1 = ϕ1(y1), [D(0+)(α)y2(t)] t=1 = ϕ2(y2), for 1 ≤ α ≤ n - 2, ϕ1, ϕ2 ∈ C([0,1], R). Our results are new and complement previously known results. As an application, we also give an example to demonstrate our result.

Generalized Nonlinear Yule Models

NASA Astrophysics Data System (ADS)

Lansky, Petr; Polito, Federico; Sacerdote, Laura

2016-11-01

With the aim of considering models related to random graphs growth exhibiting persistent memory, we propose a fractional nonlinear modification of the classical Yule model often studied in the context of macroevolution. Here the model is analyzed and interpreted in the framework of the development of networks such as the World Wide Web. Nonlinearity is introduced by replacing the linear birth process governing the growth of the in-links of each specific webpage with a fractional nonlinear birth process with completely general birth rates. Among the main results we derive the explicit distribution of the number of in-links of a webpage chosen uniformly at random recognizing the contribution to the asymptotics and the finite time correction. The mean value of the latter distribution is also calculated explicitly in the most general case. Furthermore, in order to show the usefulness of our results, we particularize them in the case of specific birth rates giving rise to a saturating behaviour, a property that is often observed in nature. The further specialization to the non-fractional case allows us to extend the Yule model accounting for a nonlinear growth.

Detection of Differential Item Functioning with Nonlinear Regression: A Non-IRT Approach Accounting for Guessing

ERIC Educational Resources Information Center

Drabinová, Adéla; Martinková, Patrícia

2017-01-01

In this article we present a general approach not relying on item response theory models (non-IRT) to detect differential item functioning (DIF) in dichotomous items with presence of guessing. The proposed nonlinear regression (NLR) procedure for DIF detection is an extension of method based on logistic regression. As a non-IRT approach, NLR can…

Analytical solutions to time-fractional partial differential equations in a two-dimensional multilayer annulus

NASA Astrophysics Data System (ADS)

Chen, Shanzhen; Jiang, Xiaoyun

2012-08-01

In this paper, analytical solutions to time-fractional partial differential equations in a multi-layer annulus are presented. The final solutions are obtained in terms of Mittag-Leffler function by using the finite integral transform technique and Laplace transform technique. In addition, the classical diffusion equation (α=1), the Helmholtz equation (α→0) and the wave equation (α=2) are discussed as special cases. Finally, an illustrative example problem for the three-layer semi-circular annular region is solved and numerical results are presented graphically for various kind of order of fractional derivative.

One Adaptive Synchronization Approach for Fractional-Order Chaotic System with Fractional-Order 1 < q < 2

PubMed Central

Zhou, Ping; Bai, Rongji

2014-01-01

Based on a new stability result of equilibrium point in nonlinear fractional-order systems for fractional-order lying in 1 < q < 2, one adaptive synchronization approach is established. The adaptive synchronization for the fractional-order Lorenz chaotic system with fractional-order 1 < q < 2 is considered. Numerical simulations show the validity and feasibility of the proposed scheme. PMID:25247207

Some operational tools for solving fractional and higher integer order differential equations: A survey on their mutual relations

NASA Astrophysics Data System (ADS)

Kiryakova, Virginia S.

2012-11-01

The Laplace Transform (LT) serves as a basis of the Operational Calculus (OC), widely explored by engineers and applied scientists in solving mathematical models for their practical needs. This transform is closely related to the exponential and trigonometric functions (exp, cos, sin) and to the classical differentiation and integration operators, reducing them to simple algebraic operations. Thus, the classical LT and the OC give useful tool to handle differential equations and systems with constant coefficients. Several generalizations of the LT have been introduced to allow solving, in a similar way, of differential equations with variable coefficients and of higher integer orders, as well as of fractional (arbitrary non-integer) orders. Note that fractional order mathematical models are recently widely used to describe better various systems and phenomena of the real world. This paper surveys briefly some of our results on classes of such integral transforms, that can be obtained from the LT by means of "transmutations" which are operators of the generalized fractional calculus (GFC). On the list of these Laplace-type integral transforms, we consider the Borel-Dzrbashjan, Meijer, Krätzel, Obrechkoff, generalized Obrechkoff (multi-index Borel-Dzrbashjan) transforms, etc. All of them are G- and H-integral transforms of convolutional type, having as kernels Meijer's G- or Fox's H-functions. Besides, some special functions (also being G- and H-functions), among them - the generalized Bessel-type and Mittag-Leffler (M-L) type functions, are generating Gel'fond-Leontiev (G-L) operators of generalized differentiation and integration, which happen to be also operators of GFC. Our integral transforms have operational properties analogous to those of the LT - they do algebrize the G-L generalized integrations and differentiations, and thus can serve for solving wide classes of differential equations with variable coefficients of arbitrary, including non-integer order

Analytical study of exact solutions of the nonlinear Korteweg-de Vries equation with space-time fractional derivatives

NASA Astrophysics Data System (ADS)

Liu, Jiangen; Zhang, Yufeng

2018-01-01

This paper gives an analytical study of dynamic behavior of the exact solutions of nonlinear Korteweg-de Vries equation with space-time local fractional derivatives. By using the improved (G‧ G )-expansion method, the explicit traveling wave solutions including periodic solutions, dark soliton solutions, soliton solutions and soliton-like solutions, are obtained for the first time. They can better help us further understand the physical phenomena and provide a strong basis. Meanwhile, some solutions are presented through 3D-graphs.

On a solution of the nonlinear differential equation for transonic flow past a wave-shaped wall

NASA Technical Reports Server (NTRS)

Kaplan, Carl

1952-01-01

The Prandtl-Busemann small-perturbation method is utilized to obtain the flow of a compressible fluid past an infinitely long wave-shaped wall. When the essential assumption for transonic flow (that all Mach numbers in the region of flow are nearly unity) is introduced, the expression for the velocity potential takes the form of a power series in the transonic similarity parameter. On the basis of this form of the solution, an attempt is made to solve the nonlinear differential equation for transonic flow past the wavy wall. The analysis utilized exhibits clearly the difficulties inherent in nonlinear-flow problems.

Sufficient conditions for asymptotic stability and stabilization of autonomous fractional order systems

NASA Astrophysics Data System (ADS)

Lenka, Bichitra Kumar; Banerjee, Soumitro

2018-03-01

We discuss the asymptotic stability of autonomous linear and nonlinear fractional order systems where the state equations contain same or different fractional orders which lie between 0 and 2. First, we use the Laplace transform method to derive some sufficient conditions which ensure asymptotic stability of linear fractional order systems. Then by using the obtained results and linearization technique, a stability theorem is presented for autonomous nonlinear fractional order system. Finally, we design a control strategy for stabilization of autonomous nonlinear fractional order systems, and apply the results to the chaotic fractional order Lorenz system in order to verify its effectiveness.

Limited knowledge of fraction representations differentiates middle school students with mathematics learning disability (dyscalculia) vs. low mathematics achievement

PubMed Central

Mazzocco, Michèle M. M.; Myers, Gwen F.; Lewis, Katherine E.; Hanich, Laurie B.; Murphy, Melissa M.

2014-01-01

Fractions pose significant challenges for many children, but for some children those challenges persist into high school. Here we administered a fractions magnitude comparison test to 122 children, from Grades 4 to 8, to test whether their knowledge of fractions typically learned early in the sequence of formal math instruction (e.g., fractions equivalent to “one-half,” and fraction pairs with common denominators) differentiates those with mathematical learning disability (MLD) versus low achievement (LA) or typical achievement (TA) in mathematics, and whether long term learning trajectories of this knowledge also differentiate these groups. We confirmed that although 4th graders with LA (n = 18) or TA (n = 93) are more accurate evaluating one-half vs. non-half fractions (until they reach ceiling performance levels on both types of fractions), children with MLD (n=11) do not show a one-half advantage until Grade 7 and do not reach ceiling performance even by Grade 8. Both the MLD and LA groups have early difficulties with fractions, but by Grade 5 the LA group approaches performance levels of the TA group and deviates from the MLD group. All groups showed a visual model advantage over Arabic number representation of fractions, but this advantage was short lived for the TA group (because ceiling level was achieved across formats), slightly more persistent for the LA group, and persisted through Grade 8 for children with MLD. Thus, difficulties with fractions persist through Grade 8 for many students, but the nature and trajectories of those difficulties varies across children with math difficulties (MLD or LA). PMID:23587941

Controllability of fractional higher order stochastic integrodifferential systems with fractional Brownian motion.

PubMed

Sathiyaraj, T; Balasubramaniam, P

2017-11-30

This paper presents a new set of sufficient conditions for controllability of fractional higher order stochastic integrodifferential systems with fractional Brownian motion (fBm) in finite dimensional space using fractional calculus, fixed point technique and stochastic analysis approach. In particular, we discuss the complete controllability for nonlinear fractional stochastic integrodifferential systems under the proved result of the corresponding linear fractional system is controllable. Finally, an example is presented to illustrate the efficiency of the obtained theoretical results. Copyright © 2017 ISA. Published by Elsevier Ltd. All rights reserved.

All-optical temporal fractional order differentiator using an in-fiber ellipsoidal air-microcavity

NASA Astrophysics Data System (ADS)

Zhang, Lihong; Sun, Shuqian; Li, Ming; Zhu, Ninghua

2017-12-01

An all-optical temporal fractional order differentiator with ultrabroad bandwidth (~1.6 THz) and extremely simple fabrication is proposed and experimentally demonstrated based on an in-fiber ellipsoidal air-microcavity. The ellipsoidal air-microcavity is fabricated by splicing a single mode fiber (SMF) and a photonic crystal fiber (PCF) together using a simple arc-discharging technology. By changing the arc-discharging times, the propagation loss can be adjusted and then the differentiation order is tuned. A nearly Gaussian-like optical pulse with 3 dB bandwidth of 8 nm is launched into the differentiator and a 0.65 order differentiation of the input pulse is achieved with a processing error of 2.55%. Project supported by the the National Natural Science Foundation of China (Nos. 61522509, 61377002, 61535012), the National High-Tech Research & Development Program of China (No. SS2015AA011002), and the Beijing Natural Science Foundation (No. 4152052). Ming Li was supported in part by the Thousand Young Talent Program.

Investigation of magnesium isotope fractionation during basalt differentiation: Implications for a chondritic composition of the terrestrial mantle

USGS Publications Warehouse

Teng, F.-Z.; Wadhwa, M.; Helz, R.T.

2007-01-01

To investigate whether magnesium isotopes are fractionated during basalt differentiation, we have performed high-precision Mg isotopic analyses by multi-collector inductively coupled plasma mass spectrometry (MC-ICP-MS) on a set of well-characterized samples from Kilauea Iki lava lake, Hawaii, USA. Samples from the Kilauea Iki lava lake, produced by closed-system crystal-melt fractionation, range from olivine-rich cumulates to highly differentiated basalts with MgO content ranging from 2.37 to 26.87??wt.%. Our results demonstrate that although these basalts have diverse chemical compositions, mineralogies, crystallization temperatures and degrees of differentiation, their Mg isotopic compositions display no measurable variation within the limits of our external precision (average ??26Mg = - 0.36 ?? 0.10 and ??25Mg = - 0.20 ?? 0.07; uncertainties are 2SD). This indicates that Mg isotopic fractionation during crystal-melt fractionation at temperatures of ??? 1055????C is undetectable at the level of precision of the current investigation. Calculations based on our data suggest that at near-magmatic temperatures the maximum fractionation in the 26Mg/24Mg ratio between olivine and melt is 0.07???. Two additional oceanic basalts, two continental basalts (BCR-1 and BCR-2), and two primitive carbonaceous chondrites (Allende and Murchison) analyzed in this study have Mg isotopic compositions similar to the Kilauea Iki lava lake samples. In contrast to a recent report [U. Wiechert, A.N. Halliday, Non-chondritic magnesium and the origins of the inner terrestrial planets, Earth and Planetary Science Letters 256 (2007) 360-371], the results presented here suggest that the Bulk Silicate Earth has a chondritic Mg isotopic composition. ?? 2007.

Limited knowledge of fraction representations differentiates middle school students with mathematics learning disability (dyscalculia) versus low mathematics achievement.

PubMed

Mazzocco, Michèle M M; Myers, Gwen F; Lewis, Katherine E; Hanich, Laurie B; Murphy, Melissa M

2013-06-01

Fractions pose significant challenges for many children, but for some children those challenges persist into high school. Here we administered a fractions magnitude comparison test to 122 children, from Grades 4 to 8, to test whether their knowledge of fractions typically learned early in the sequence of formal math instruction (e.g., fractions equivalent to one-half, fraction pairs with common denominators) differentiates those with mathematics learning disability (MLD) versus low achievement (LA) or typical achievement (TA) in mathematics and whether long-term learning trajectories of this knowledge also differentiate these groups. We confirmed that although fourth graders with TA (n=93) were more accurate in evaluating "one-half" fractions than in evaluating "non-half" fractions (until they reached ceiling performance levels on both types of fractions), children with MLD (n=11) did not show a one-half advantage until Grade 7 and did not reach ceiling performance even by Grade 8. Both the MLD and LA groups had early difficulties with fractions, but by Grade 5 the LA group approached performance levels of the TA group and deviated from the MLD group. All groups showed a visual model advantage over Arabic number representation of fractions, but this advantage was short-lived for the TA group (because ceiling level was achieved across formats), whereas it was slightly more persistent for the LA group and persisted through Grade 8 for children with MLD. Thus, difficulties with fractions persist through Grade 8 for many students, but the nature and trajectories of those difficulties vary across children with math difficulties (MLD or LA). Copyright © 2013 Elsevier Inc. All rights reserved.

Varieties of operator manipulation. [for solving differential equations and calculating finite differences

NASA Technical Reports Server (NTRS)

Doohovskoy, A.

1977-01-01

A change in MACSYMA syntax is proposed to accommodate the operator manipulators necessary to implement direct and indirect methods for the solution of differential equations, calculus of finite differences, and the fractional calculus, as well as their modern counterparts. To illustrate the benefits and convenience of this syntax extension, an example is given to show how MACSYMA's pattern-matching capability can be used to implement a particular set of operator identities which can then be used to obtain exact solutions to nonlinear differential equations.

Fractional Diffusion Equations and Anomalous Diffusion

NASA Astrophysics Data System (ADS)

Evangelista, Luiz Roberto; Kaminski Lenzi, Ervin

2018-01-01

Preface; 1. Mathematical preliminaries; 2. A survey of the fractional calculus; 3. From normal to anomalous diffusion; 4. Fractional diffusion equations: elementary applications; 5. Fractional diffusion equations: surface effects; 6. Fractional nonlinear diffusion equation; 7. Anomalous diffusion: anisotropic case; 8. Fractional Schrödinger equations; 9. Anomalous diffusion and impedance spectroscopy; 10. The Poisson–Nernst–Planck anomalous (PNPA) models; References; Index.

Nonlinear differential equations for the wavefront surface at arbitrary Hartmann-plane distances.

PubMed

Téllez-Quiñones, Alejandro; Malacara-Doblado, Daniel; Flores-Hernández, Ricardo; Gutiérrez-Hernández, David A; León-Rodríguez, Miguel

2016-03-20

In the Hartmann test, a wave aberration function W is estimated from the information of the spot diagram drawn in an observation plane. The distance from a reference plane to the observation plane, the Hartmann-plane distance, is typically chosen as z=f, where f is the radius of a reference sphere. The function W and the transversal aberrations {X,Y} calculated at the plane z=f are related by two well-known linear differential equations. Here, we propose two nonlinear differential equations to denote a more general relation between W and the transversal aberrations {U,V} calculated at any arbitrary Hartmann-plane distance z=r. We also show how to directly estimate the wavefront surface w from the information of {U,V}. The use of arbitrary r values could improve the reliability of the measurements of W, or w, when finding difficulties in adequate ray identification at z=f.

Quantitative label-free multimodality nonlinear optical imaging for in situ differentiation of cancerous lesions

NASA Astrophysics Data System (ADS)

Xu, Xiaoyun; Li, Xiaoyan; Cheng, Jie; Liu, Zhengfan; Thrall, Michael J.; Wang, Xi; Wang, Zhiyong; Wong, Stephen T. C.

2013-03-01

The development of real-time, label-free imaging techniques has recently attracted research interest for in situ differentiation of cancerous lesions from normal tissues. Molecule-specific intrinsic contrast can arise from label-free imaging techniques such as Coherent Anti-Stokes Raman Scattering (CARS), Two-Photon Excited AutoFluorescence (TPEAF), and Second Harmonic Generation (SHG), which, in combination, would hold the promise of a powerful label-free tool for cancer diagnosis. Among cancer-related deaths, lung carcinoma is the leading cause for both sexes. Although early treatment can increase the survival rate dramatically, lesion detection and precise diagnosis at an early stage is unusual due to its asymptomatic nature and limitations of current diagnostic techniques that make screening difficult. We investigated the potential of using multimodality nonlinear optical microscopy that incorporates CARS, TPEAF, and SHG techniques for differentiation of lung cancer from normal tissue. Cancerous and non-cancerous lung tissue samples from patients were imaged using CARS, TPEAF, and SHG techniques for comparison. These images showed good pathology correlation with hematoxylin and eosin (H and E) stained sections from the same tissue samples. Ongoing work includes imaging at various penetration depths to show three-dimensional morphologies of tumor cell nuclei using CARS, elastin using TPEAF, and collagen using SHG and developing classification algorithms for quantitative feature extraction to enable lung cancer diagnosis. Our results indicate that via real-time morphology analyses, a multimodality nonlinear optical imaging platform potentially offers a powerful minimally-invasive way to differentiate cancer lesions from surrounding non-tumor tissues in vivo for clinical applications.

Fractional Poisson Fields and Martingales

NASA Astrophysics Data System (ADS)

Aletti, Giacomo; Leonenko, Nikolai; Merzbach, Ely

2018-02-01

We present new properties for the Fractional Poisson process (FPP) and the Fractional Poisson field on the plane. A martingale characterization for FPPs is given. We extend this result to Fractional Poisson fields, obtaining some other characterizations. The fractional differential equations are studied. We consider a more general Mixed-Fractional Poisson process and show that this process is the stochastic solution of a system of fractional differential-difference equations. Finally, we give some simulations of the Fractional Poisson field on the plane.

3D early embryogenesis image filtering by nonlinear partial differential equations.

PubMed

Krivá, Z; Mikula, K; Peyriéras, N; Rizzi, B; Sarti, A; Stasová, O

2010-08-01

We present nonlinear diffusion equations, numerical schemes to solve them and their application for filtering 3D images obtained from laser scanning microscopy (LSM) of living zebrafish embryos, with a goal to identify the optimal filtering method and its parameters. In the large scale applications dealing with analysis of 3D+time embryogenesis images, an important objective is a correct detection of the number and position of cell nuclei yielding the spatio-temporal cell lineage tree of embryogenesis. The filtering is the first and necessary step of the image analysis chain and must lead to correct results, removing the noise, sharpening the nuclei edges and correcting the acquisition errors related to spuriously connected subregions. In this paper we study such properties for the regularized Perona-Malik model and for the generalized mean curvature flow equations in the level-set formulation. A comparison with other nonlinear diffusion filters, like tensor anisotropic diffusion and Beltrami flow, is also included. All numerical schemes are based on the same discretization principles, i.e. finite volume method in space and semi-implicit scheme in time, for solving nonlinear partial differential equations. These numerical schemes are unconditionally stable, fast and naturally parallelizable. The filtering results are evaluated and compared first using the Mean Hausdorff distance between a gold standard and different isosurfaces of original and filtered data. Then, the number of isosurface connected components in a region of interest (ROI) detected in original and after the filtering is compared with the corresponding correct number of nuclei in the gold standard. Such analysis proves the robustness and reliability of the edge preserving nonlinear diffusion filtering for this type of data and lead to finding the optimal filtering parameters for the studied models and numerical schemes. Further comparisons consist in ability of splitting the very close objects which

Time-Domain Evaluation of Fractional Order Controllers’ Direct Discretization Methods

NASA Astrophysics Data System (ADS)

Ma, Chengbin; Hori, Yoichi

Fractional Order Control (FOC), in which the controlled systems and/or controllers are described by fractional order differential equations, has been applied to various control problems. Though it is not difficult to understand FOC’s theoretical superiority, realization issue keeps being somewhat problematic. Since the fractional order systems have an infinite dimension, proper approximation by finite difference equation is needed to realize the designed fractional order controllers. In this paper, the existing direct discretization methods are evaluated by their convergences and time-domain comparison with the baseline case. Proposed sampling time scaling property is used to calculate the baseline case with full memory length. This novel discretization method is based on the classical trapezoidal rule but with scaled sampling time. Comparative studies show good performance and simple algorithm make the Short Memory Principle method most practically superior. The FOC research is still at its primary stage. But its applications in modeling and robustness against non-linearities reveal the promising aspects. Parallel to the development of FOC theories, applying FOC to various control problems is also crucially important and one of top priority issues.

Modelling groundwater fractal flow with fractional differentiation via Mittag-Leffler law

NASA Astrophysics Data System (ADS)

Ahokposi, D. P.; Atangana, Abdon; Vermeulen, D. P.

2017-04-01

Modelling the flow of groundwater within a network of fractures is perhaps one of the most difficult exercises within the field of geohydrology. This physical problem has attracted the attention of several scientists across the globe. Already two different types of differentiations have been used to attempt modelling this problem including the classical and the fractional differentiation. In this paper, we employed the most recent concept of differentiation based on the non-local and non-singular kernel called the generalized Mittag-Leffler function, to reshape the model of groundwater fractal flow. We presented the existence of positive solution of the new model. Using the fixed-point approach, we established the uniqueness of the positive solution. We solve the new model with three different numerical schemes including implicit, explicit and Crank-Nicholson numerical methods. Experimental data collected from four constant discharge tests conducted in a typical fractured crystalline rock aquifer of the Northern Limb (Bushveld Complex) in the Limpopo Province (South Africa) are compared with the numerical solutions. It is worth noting that the four boreholes (BPAC1, BPAC2, BPAC3, and BPAC4) are located on Faults.

A study of fractional Schrödinger equation composed of Jumarie fractional derivative

NASA Astrophysics Data System (ADS)

Banerjee, Joydip; Ghosh, Uttam; Sarkar, Susmita; Das, Shantanu

2017-04-01

In this paper we have derived the fractional-order Schrödinger equation composed of Jumarie fractional derivative. The solution of this fractional-order Schrödinger equation is obtained in terms of Mittag-Leffler function with complex arguments, and fractional trigonometric functions. A few important properties of the fractional Schrödinger equation are then described for the case of particles in one-dimensional infinite potential well. One of the motivations for using fractional calculus in physical systems is that the space and time variables, which we often deal with, exhibit coarse-grained phenomena. This means infinitesimal quantities cannot be arbitrarily taken to zero - rather they are non-zero with a minimum spread. This type of non-zero spread arises in the microscopic to mesoscopic levels of system dynamics, which means that, if we denote x as the point in space and t as the point in time, then limit of the differentials d x (and d t) cannot be taken as zero. To take the concept of coarse graining into account, use the infinitesimal quantities as (Δ x) α (and (Δ t) α ) with 0 < α < 1; called as `fractional differentials'. For arbitrarily small Δ x and Δ t (tending towards zero), these `fractional' differentials are greater than Δ x (and Δ t), i.e. (Δ x) α > Δ x and (Δ t) α > Δ t. This way of defining the fractional differentials helps us to use fractional derivatives in the study of dynamic systems.

Non-linear patterns in age-related DNA methylation may reflect CD4+ T cell differentiation

PubMed Central

Johnson, Nicholas D.; Wiener, Howard W.; Smith, Alicia K.; Nishitani, Shota; Absher, Devin M.; Arnett, Donna K.; Aslibekyan, Stella; Conneely, Karen N.

2017-01-01

ABSTRACT DNA methylation (DNAm) is an important epigenetic process involved in the regulation of gene expression. While many studies have identified thousands of loci associated with age, few have differentiated between linear and non-linear DNAm trends with age. Non-linear trends could indicate early- or late-life gene regulatory processes. Using data from the Illumina 450K array on 336 human peripheral blood samples, we identified 21 CpG sites that associated with age (P<1.03E-7) and exhibited changing rates of DNAm change with age (P<1.94E-6). For 2 of these CpG sites (cg07955995 and cg22285878), DNAm increased with age at an increasing rate, indicating that differential DNAm was greatest among elderly individuals. We observed significant replication for both CpG sites (P<5.0E-8) in a second set of peripheral blood samples. In 8 of 9 additional data sets comprising samples of monocytes, T cell subtypes, and brain tissue, we observed a pattern directionally consistent with DNAm increasing with age at an increasing rate, which was nominally significant in the 3 largest data sets (4.3E-15differentiation via the repression of FOXP3. These findings may suggest a possible role for cg07955995 and cg22285878 in immunosenescence. PMID:28387568

Continuum Modeling and Control of Large Nonuniform Wireless Networks via Nonlinear Partial Differential Equations

DOE PAGES

Zhang, Yang; Chong, Edwin K. P.; Hannig, Jan; ...

2013-01-01

We inmore » troduce a continuum modeling method to approximate a class of large wireless networks by nonlinear partial differential equations (PDEs). This method is based on the convergence of a sequence of underlying Markov chains of the network indexed by N , the number of nodes in the network. As N goes to infinity, the sequence converges to a continuum limit, which is the solution of a certain nonlinear PDE. We first describe PDE models for networks with uniformly located nodes and then generalize to networks with nonuniformly located, and possibly mobile, nodes. Based on the PDE models, we develop a method to control the transmissions in nonuniform networks so that the continuum limit is invariant under perturbations in node locations. This enables the networks to maintain stable global characteristics in the presence of varying node locations.« less

Parachute dynamics and stability analysis. [using nonlinear differential equations of motion

NASA Technical Reports Server (NTRS)

Ibrahim, S. K.; Engdahl, R. A.

1974-01-01

The nonlinear differential equations of motion for a general parachute-riser-payload system are developed. The resulting math model is then applied for analyzing the descent dynamics and stability characteristics of both the drogue stabilization phase and the main descent phase of the space shuttle solid rocket booster (SRB) recovery system. The formulation of the problem is characterized by a minimum number of simplifying assumptions and full application of state-of-the-art parachute technology. The parachute suspension lines and the parachute risers can be modeled as elastic elements, and the whole system may be subjected to specified wind and gust profiles in order to assess their effects on the stability of the recovery system.

A novel interval type-2 fractional order fuzzy PID controller: Design, performance evaluation, and its optimal time domain tuning.

PubMed

Kumar, Anupam; Kumar, Vijay

2017-05-01

In this paper, a novel concept of an interval type-2 fractional order fuzzy PID (IT2FO-FPID) controller, which requires fractional order integrator and fractional order differentiator, is proposed. The incorporation of Takagi-Sugeno-Kang (TSK) type interval type-2 fuzzy logic controller (IT2FLC) with fractional controller of PID-type is investigated for time response measure due to both unit step response and unit load disturbance. The resulting IT2FO-FPID controller is examined on different delayed linear and nonlinear benchmark plants followed by robustness analysis. In order to design this controller, fractional order integrator-differentiator operators are considered as design variables including input-output scaling factors. A new hybridized algorithm named as artificial bee colony-genetic algorithm (ABC-GA) is used to optimize the parameters of the controller while minimizing weighted sum of integral of time absolute error (ITAE) and integral of square of control output (ISCO). To assess the comparative performance of the IT2FO-FPID, authors compared it against existing controllers, i.e., interval type-2 fuzzy PID (IT2-FPID), type-1 fractional order fuzzy PID (T1FO-FPID), type-1 fuzzy PID (T1-FPID), and conventional PID controllers. Furthermore, to show the effectiveness of the proposed controller, the perturbed processes along with the larger dead time are tested. Moreover, the proposed controllers are also implemented on multi input multi output (MIMO), coupled, and highly complex nonlinear two-link robot manipulator system in presence of un-modeled dynamics. Finally, the simulation results explicitly indicate that the performance of the proposed IT2FO-FPID controller is superior to its conventional counterparts in most of the cases. Copyright © 2017 ISA. Published by Elsevier Ltd. All rights reserved.

Fractional dynamics pharmacokinetics–pharmacodynamic models

PubMed Central

2010-01-01

While an increasing number of fractional order integrals and differential equations applications have been reported in the physics, signal processing, engineering and bioengineering literatures, little attention has been paid to this class of models in the pharmacokinetics–pharmacodynamic (PKPD) literature. One of the reasons is computational: while the analytical solution of fractional differential equations is available in special cases, it this turns out that even the simplest PKPD models that can be constructed using fractional calculus do not allow an analytical solution. In this paper, we first introduce new families of PKPD models incorporating fractional order integrals and differential equations, and, second, exemplify and investigate their qualitative behavior. The families represent extensions of frequently used PK link and PD direct and indirect action models, using the tools of fractional calculus. In addition the PD models can be a function of a variable, the active drug, which can smoothly transition from concentration to exposure, to hyper-exposure, according to a fractional integral transformation. To investigate the behavior of the models we propose, we implement numerical algorithms for fractional integration and for the numerical solution of a system of fractional differential equations. For simplicity, in our investigation we concentrate on the pharmacodynamic side of the models, assuming standard (integer order) pharmacokinetics. PMID:20455076

Transport equations for subdiffusion with nonlinear particle interaction.

PubMed

Straka, P; Fedotov, S

2015-02-07

We show how the nonlinear interaction effects 'volume filling' and 'adhesion' can be incorporated into the fractional subdiffusive transport of cells and individual organisms. To this end, we use microscopic random walk models with anomalous trapping and systematically derive generic non-Markovian and nonlinear governing equations for the mean concentrations of the subdiffusive cells or organisms. We uncover an interesting interaction between the nonlinearities and the non-Markovian nature of the transport. In the subdiffusive case, this interaction manifests itself in a nontrivial combination of nonlinear terms with fractional derivatives. In the long time limit, however, these equations simplify to a form without fractional operators. This provides an easy method for the study of aggregation phenomena. In particular, this enables us to show that volume filling can prevent "anomalous aggregation," which occurs in subdiffusive systems with a spatially varying anomalous exponent. Copyright © 2014 Elsevier Ltd. All rights reserved.

Darboux transformation and solitons for an integrable nonautonomous nonlinear integro-differential Schrödinger equation

NASA Astrophysics Data System (ADS)

Yong, Xuelin; Fan, Yajing; Huang, Yehui; Ma, Wen-Xiu; Tian, Jing

2017-10-01

By modifying the scheme for an isospectral problem, the non-isospectral Ablowitz-Kaup-Newell-Segur (AKNS) hierarchy is constructed via allowing the time varying spectrum. In this paper, we consider an integrable nonautonomous nonlinear integro-differential Schrödinger equation discussed before in “Multi-soliton management by the integrable nonautonomous nonlinear integro-differential Schrödinger equation” [Y. J. Zhang, D. Zhao and H. G. Luo, Ann. Phys. 350 (2014) 112]. We first analyze the integrability conditions and identify the model. Second, we modify the existing Darboux transformation (DT) for such a non-isospectral problem. Third, the nonautonomous soliton solutions are obtained via the resulting DT and basic properties of these solutions in the inhomogeneous media are discussed graphically to illustrate the influences of the variable coefficients. In the process, a technique by selecting appropriate spectral parameters instead of the variable inhomogeneities is employed to realize a different type of one-soliton management. Several novel optical solitons are constructed and their features are shown by some specific figures. In addition, four kinds of the special localized two-soliton solutions are obtained. The solitonic excitations localized both in space and time, which exhibit the feature of the so-called rogue waves but with a zero background, are discussed.

Analysis of blood pressure signal in patients with different ventricular ejection fraction using linear and non-linear methods.

PubMed

Arcentales, Andres; Rivera, Patricio; Caminal, Pere; Voss, Andreas; Bayes-Genis, Antonio; Giraldo, Beatriz F

2016-08-01

Changes in the left ventricle function produce alternans in the hemodynamic and electric behavior of the cardiovascular system. A total of 49 cardiomyopathy patients have been studied based on the blood pressure signal (BP), and were classified according to the left ventricular ejection fraction (LVEF) in low risk (LR: LVEF>35%, 17 patients) and high risk (HR: LVEF≤35, 32 patients) groups. We propose to characterize these patients using a linear and a nonlinear methods, based on the spectral estimation and the recurrence plot, respectively. From BP signal, we extracted each systolic time interval (STI), upward systolic slope (BPsl), and the difference between systolic and diastolic BP, defined as pulse pressure (PP). After, the best subset of parameters were obtained through the sequential feature selection (SFS) method. According to the results, the best classification was obtained using a combination of linear and nonlinear features from STI and PP parameters. For STI, the best combination was obtained considering the frequency peak and the diagonal structures of RP, with an area under the curve (AUC) of 79%. The same results were obtained when comparing PP values. Consequently, the use of combined linear and nonlinear parameters could improve the risk stratification of cardiomyopathy patients.

The numerical dynamic for highly nonlinear partial differential equations

NASA Technical Reports Server (NTRS)

Lafon, A.; Yee, H. C.

1992-01-01

Problems associated with the numerical computation of highly nonlinear equations in computational fluid dynamics are set forth and analyzed in terms of the potential ranges of spurious behaviors. A reaction-convection equation with a nonlinear source term is employed to evaluate the effects related to spatial and temporal discretizations. The discretization of the source term is described according to several methods, and the various techniques are shown to have a significant effect on the stability of the spurious solutions. Traditional linearized stability analyses cannot provide the level of confidence required for accurate fluid dynamics computations, and the incorporation of nonlinear analysis is proposed. Nonlinear analysis based on nonlinear dynamical systems complements the conventional linear approach and is valuable in the analysis of hypersonic aerodynamics and combustion phenomena.

On new fractional Hermite-Hadamard type inequalities for n-time differentiable quasi-convex functions and P-functions

NASA Astrophysics Data System (ADS)

Set, Erhan; Özdemir, M. Emin; Alan, E. Aykan

2017-04-01

In this article, by using the Hölder's inequality and power mean inequality the authors establish several inequalities of Hermite-Hadamard type for n- time differentiable quasi-convex functions and P- functions involving Riemann-Liouville fractional integrals.

Adaptive neural network backstepping control for a class of uncertain fractional-order chaotic systems with unknown backlash-like hysteresis

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wu, Yimin; Lv, Hui, E-mail: lvhui207@gmail.com

In this paper, we consider the control problem of a class of uncertain fractional-order chaotic systems preceded by unknown backlash-like hysteresis nonlinearities based on backstepping control algorithm. We model the hysteresis by using a differential equation. Based on the fractional Lyapunov stability criterion and the backstepping algorithm procedures, an adaptive neural network controller is driven. No knowledge of the upper bound of the disturbance and system uncertainty is required in our controller, and the asymptotical convergence of the tracking error can be guaranteed. Finally, we give two simulation examples to confirm our theoretical results.

Constructing general partial differential equations using polynomial and neural networks.

PubMed

Zjavka, Ladislav; Pedrycz, Witold

2016-01-01

Sum fraction terms can approximate multi-variable functions on the basis of discrete observations, replacing a partial differential equation definition with polynomial elementary data relation descriptions. Artificial neural networks commonly transform the weighted sum of inputs to describe overall similarity relationships of trained and new testing input patterns. Differential polynomial neural networks form a new class of neural networks, which construct and solve an unknown general partial differential equation of a function of interest with selected substitution relative terms using non-linear multi-variable composite polynomials. The layers of the network generate simple and composite relative substitution terms whose convergent series combinations can describe partial dependent derivative changes of the input variables. This regression is based on trained generalized partial derivative data relations, decomposed into a multi-layer polynomial network structure. The sigmoidal function, commonly used as a nonlinear activation of artificial neurons, may transform some polynomial items together with the parameters with the aim to improve the polynomial derivative term series ability to approximate complicated periodic functions, as simple low order polynomials are not able to fully make up for the complete cycles. The similarity analysis facilitates substitutions for differential equations or can form dimensional units from data samples to describe real-world problems. Copyright © 2015 Elsevier Ltd. All rights reserved.

HL-60 differentiating activity and flavonoid content of the readily extractable fraction prepared from citrus juices.

PubMed

Kawaii, S; Tomono, Y; Katase, E; Ogawa, K; Yano, M

1999-01-01

Citrus plants are rich sources of various bioactive flavonoids. To eliminate masking effects caused by hesperidin, naringin, and neoeriocitrin, the abundant flavonoid glycosides which make up 90% of the conventionally prepared sample, the readily extractable fraction from Citrus juice was prepared by adsorbing on HP-20 resin and eluting with EtOH and acetone from the resin and was subjected to HL-60 differentiation assay and quantitative analysis of major flavonoids. Screening of 34 Citrus juices indicated that King (C. nobilis) had a potent activity for inducing differentiation of HL-60, and the active principles were isolated and identified as four polymethoxylated flavonoids, namely, nobiletin, 3,3',4',5,6,7, 8-heptamethoxyflavone, natsudaidain, and tangeretin. HPLC analysis of the readily extractable fraction also indicated that King contained high amounts of these polymethoxylated flavonoids among the Citrus juices examined. Principal component and cluster analyses of the readily extractable flavonoids indicated peculiarities of King and Bergamot.

Geometrically nonlinear resonance of higher-order shear deformable functionally graded carbon-nanotube-reinforced composite annular sector plates excited by harmonic transverse loading

NASA Astrophysics Data System (ADS)

Gholami, Raheb; Ansari, Reza

2018-02-01

This article presents an attempt to study the nonlinear resonance of functionally graded carbon-nanotube-reinforced composite (FG-CNTRC) annular sector plates excited by a uniformly distributed harmonic transverse load. To this purpose, first, the extended rule of mixture including the efficiency parameters is employed to approximately obtain the effective material properties of FG-CNTRC annular sector plates. Then, the focus is on presenting the weak form of discretized mathematical formulation of governing equations based on the variational differential quadrature (VDQ) method and Hamilton's principle. The geometric nonlinearity and shear deformation effects are considered based on the von Kármán assumptions and Reddy's third-order shear deformation plate theory, respectively. The discretization process is performed via the generalized differential quadrature (GDQ) method together with numerical differential and integral operators. Then, an efficient multi-step numerical scheme is used to obtain the nonlinear dynamic behavior of the FG-CNTRC annular sector plates near their primary resonance as the frequency-response curve. The accuracy of the present results is first verified and then a parametric study is presented to show the impacts of CNT volume fraction, CNT distribution pattern, geometry of annular sector plate and sector angle on the nonlinear frequency-response curve of FG-CNTRC annular sector plates with different edge supports.

Analysis of a class of boundary value problems depending on left and right Caputo fractional derivatives

NASA Astrophysics Data System (ADS)

Antunes, Pedro R. S.; Ferreira, Rui A. C.

2017-07-01

In this work we study boundary value problems associated to a nonlinear fractional ordinary differential equation involving left and right Caputo derivatives. We discuss the regularity of the solutions of such problems and, in particular, give precise necessary conditions so that the solutions are C1([0, 1]). Taking into account our analytical results, we address the numerical solution of those problems by the augmented -RBF method. Several examples illustrate the good performance of the numerical method.

On the comparison of perturbation-iteration algorithm and residual power series method to solve fractional Zakharov-Kuznetsov equation

NASA Astrophysics Data System (ADS)

Şenol, Mehmet; Alquran, Marwan; Kasmaei, Hamed Daei

2018-06-01

In this paper, we present analytic-approximate solution of time-fractional Zakharov-Kuznetsov equation. This model demonstrates the behavior of weakly nonlinear ion acoustic waves in a plasma bearing cold ions and hot isothermal electrons in the presence of a uniform magnetic field. Basic definitions of fractional derivatives are described in the Caputo sense. Perturbation-iteration algorithm (PIA) and residual power series method (RPSM) are applied to solve this equation with success. The convergence analysis is also presented for both methods. Numerical results are given and then they are compared with the exact solutions. Comparison of the results reveal that both methods are competitive, powerful, reliable, simple to use and ready to apply to wide range of fractional partial differential equations.

Multiplicative noise removal through fractional order tv-based model and fast numerical schemes for its approximation

NASA Astrophysics Data System (ADS)

Ullah, Asmat; Chen, Wen; Khan, Mushtaq Ahmad

2017-07-01

This paper introduces a fractional order total variation (FOTV) based model with three different weights in the fractional order derivative definition for multiplicative noise removal purpose. The fractional-order Euler Lagrange equation which is a highly non-linear partial differential equation (PDE) is obtained by the minimization of the energy functional for image restoration. Two numerical schemes namely an iterative scheme based on the dual theory and majorization- minimization algorithm (MMA) are used. To improve the restoration results, we opt for an adaptive parameter selection procedure for the proposed model by applying the trial and error method. We report numerical simulations which show the validity and state of the art performance of the fractional-order model in visual improvement as well as an increase in the peak signal to noise ratio comparing to corresponding methods. Numerical experiments also demonstrate that MMAbased methodology is slightly better than that of an iterative scheme.

Mass dependent fractionation of stable chromium isotopes in mare basalts: Implications for the formation and the differentiation of the Moon

NASA Astrophysics Data System (ADS)

Bonnand, Pierre; Parkinson, Ian J.; Anand, Mahesh

2016-02-01

We present the first stable chromium isotopic data from mare basalts in order to investigate the similarity between the Moon and the Earth's mantle. A double spike technique coupled with MC-ICP-MS measurements was used to analyse 19 mare basalts, comprising high-Ti, low-Ti and KREEP-rich varieties. Chromium isotope ratios (δ53Cr) for mare basalts are positively correlated with indices of magmatic differentiation such as Mg# and Cr concentration which suggests that Cr isotopes were fractionated during magmatic differentiation. Modelling of the results provides evidence that spinel and pyroxene are the main phases controlling the Cr isotopic composition during fractional crystallisation. The most evolved samples have the lightest isotopic compositions, complemented by cumulates that are isotopically heavy. Two hypotheses are proposed to explain this fractionation: (i) equilibrium fractionation where heavy isotopes are preferentially incorporated into the spinel lattice and (ii) a difference in isotopic composition between Cr2+ and Cr3+ in the melt. However, both processes require magmatic temperatures below 1200 °C for appreciable Cr3+ to be present at the low oxygen fugacities found in the Moon (IW -1 to -2 log units). There is no isotopic difference between the most primitive high-Ti, low-Ti and KREEP basalts, which suggest that the sources of these basalts were homogeneous in terms of stable Cr isotopes. The least differentiated sample in our sample set is the low-Ti basalt 12016, characterised by a Cr isotopic composition of -0.222 ± 0.025‰, which is within error of the current BSE value (-0.124 ± 0.101‰). The similarity between the mantles of the Moon and Earth is consistent with a terrestrial origin for a major fraction of the lunar Cr. This similarity also suggests that Cr isotopes were not fractionated by core formation on the Moon.

Alexander fractional differential window filter for ECG denoising.

PubMed

Verma, Atul Kumar; Saini, Indu; Saini, Barjinder Singh

2018-06-01

The electrocardiogram (ECG) non-invasively monitors the electrical activities of the heart. During the process of recording and transmission, ECG signals are often corrupted by various types of noises. Minimizations of these noises facilitate accurate detection of various anomalies. In the present paper, Alexander fractional differential window (AFDW) filter is proposed for ECG signal denoising. The designed filter is based on the concept of generalized Alexander polynomial and the R-L differential equation of fractional calculus. This concept is utilized to formulate a window that acts as a forward filter. Thereafter, the backward filter is constructed by reversing the coefficients of the forward filter. The proposed AFDW filter is then obtained by averaging of the forward and backward filter coefficients. The performance of the designed AFDW filter is validated by adding the various type of noise to the original ECG signal obtained from MIT-BIH arrhythmia database. The two non-diagnostic measure, i.e., SNR, MSE, and one diagnostic measure, i.e., wavelet energy based diagnostic distortion (WEDD) have been employed for the quantitative evaluation of the designed filter. Extensive experimentations on all the 48-records of MIT-BIH arrhythmia database resulted in average SNR of 22.014 ± 3.806365, 14.703 ± 3.790275, 13.3183 ± 3.748230; average MSE of 0.001458 ± 0.00028, 0.0078 ± 0.000319, 0.01061 ± 0.000472; and average WEDD value of 0.020169 ± 0.01306, 0.1207 ± 0.061272, 0.1432 ± 0.073588, for ECG signal contaminated by the power line, random, and the white Gaussian noise respectively. A new metric named as morphological power preservation measure (MPPM) is also proposed that account for the power preservance (as indicated by PSD plots) and the QRS morphology. The proposed AFDW filter retained much of the original (clean) signal power without any significant morphological distortion as validated by MPPM measure that were 0

Fast and high-order numerical algorithms for the solution of multidimensional nonlinear fractional Ginzburg-Landau equation

NASA Astrophysics Data System (ADS)

Mohebbi, Akbar

2018-02-01

In this paper we propose two fast and accurate numerical methods for the solution of multidimensional space fractional Ginzburg-Landau equation (FGLE). In the presented methods, to avoid solving a nonlinear system of algebraic equations and to increase the accuracy and efficiency of method, we split the complex problem into simpler sub-problems using the split-step idea. For a homogeneous FGLE, we propose a method which has fourth-order of accuracy in time component and spectral accuracy in space variable and for nonhomogeneous one, we introduce another scheme based on the Crank-Nicolson approach which has second-order of accuracy in time variable. Due to using the Fourier spectral method for fractional Laplacian operator, the resulting schemes are fully diagonal and easy to code. Numerical results are reported in terms of accuracy, computational order and CPU time to demonstrate the accuracy and efficiency of the proposed methods and to compare the results with the analytical solutions. The results show that the present methods are accurate and require low CPU time. It is illustrated that the numerical results are in good agreement with the theoretical ones.

A Comparison of Two-Stage Approaches for Fitting Nonlinear Ordinary Differential Equation Models with Mixed Effects.

PubMed

Chow, Sy-Miin; Bendezú, Jason J; Cole, Pamela M; Ram, Nilam

2016-01-01

Several approaches exist for estimating the derivatives of observed data for model exploration purposes, including functional data analysis (FDA; Ramsay & Silverman, 2005 ), generalized local linear approximation (GLLA; Boker, Deboeck, Edler, & Peel, 2010 ), and generalized orthogonal local derivative approximation (GOLD; Deboeck, 2010 ). These derivative estimation procedures can be used in a two-stage process to fit mixed effects ordinary differential equation (ODE) models. While the performance and utility of these routines for estimating linear ODEs have been established, they have not yet been evaluated in the context of nonlinear ODEs with mixed effects. We compared properties of the GLLA and GOLD to an FDA-based two-stage approach denoted herein as functional ordinary differential equation with mixed effects (FODEmixed) in a Monte Carlo (MC) study using a nonlinear coupled oscillators model with mixed effects. Simulation results showed that overall, the FODEmixed outperformed both the GLLA and GOLD across all the embedding dimensions considered, but a novel use of a fourth-order GLLA approach combined with very high embedding dimensions yielded estimation results that almost paralleled those from the FODEmixed. We discuss the strengths and limitations of each approach and demonstrate how output from each stage of FODEmixed may be used to inform empirical modeling of young children's self-regulation.

Robust variable selection method for nonparametric differential equation models with application to nonlinear dynamic gene regulatory network analysis.

PubMed

Lu, Tao

2016-01-01

The gene regulation network (GRN) evaluates the interactions between genes and look for models to describe the gene expression behavior. These models have many applications; for instance, by characterizing the gene expression mechanisms that cause certain disorders, it would be possible to target those genes to block the progress of the disease. Many biological processes are driven by nonlinear dynamic GRN. In this article, we propose a nonparametric differential equation (ODE) to model the nonlinear dynamic GRN. Specially, we address following questions simultaneously: (i) extract information from noisy time course gene expression data; (ii) model the nonlinear ODE through a nonparametric smoothing function; (iii) identify the important regulatory gene(s) through a group smoothly clipped absolute deviation (SCAD) approach; (iv) test the robustness of the model against possible shortening of experimental duration. We illustrate the usefulness of the model and associated statistical methods through a simulation and a real application examples.

A Numerical Scheme for Ordinary Differential Equations Having Time Varying and Nonlinear Coefficients Based on the State Transition Matrix

NASA Technical Reports Server (NTRS)

Bartels, Robert E.

2002-01-01

A variable order method of integrating initial value ordinary differential equations that is based on the state transition matrix has been developed. The method has been evaluated for linear time variant and nonlinear systems of equations. While it is more complex than most other methods, it produces exact solutions at arbitrary time step size when the time variation of the system can be modeled exactly by a polynomial. Solutions to several nonlinear problems exhibiting chaotic behavior have been computed. Accuracy of the method has been demonstrated by comparison with an exact solution and with solutions obtained by established methods.

On mixed derivatives type high dimensional multi-term fractional partial differential equations approximate solutions

NASA Astrophysics Data System (ADS)

Talib, Imran; Belgacem, Fethi Bin Muhammad; Asif, Naseer Ahmad; Khalil, Hammad

2017-01-01

In this research article, we derive and analyze an efficient spectral method based on the operational matrices of three dimensional orthogonal Jacobi polynomials to solve numerically the mixed partial derivatives type multi-terms high dimensions generalized class of fractional order partial differential equations. We transform the considered fractional order problem to an easily solvable algebraic equations with the aid of the operational matrices. Being easily solvable, the associated algebraic system leads to finding the solution of the problem. Some test problems are considered to confirm the accuracy and validity of the proposed numerical method. The convergence of the method is ensured by comparing our Matlab software simulations based obtained results with the exact solutions in the literature, yielding negligible errors. Moreover, comparative results discussed in the literature are extended and improved in this study.

Numerical solution of nonlinear partial differential equations of mixed type. [finite difference approximation

NASA Technical Reports Server (NTRS)

Jameson, A.

1976-01-01

A review is presented of some recently developed numerical methods for the solution of nonlinear equations of mixed type. The methods considered use finite difference approximations to the differential equation. Central difference formulas are employed in the subsonic zone and upwind difference formulas are used in the supersonic zone. The relaxation method for the small disturbance equation is discussed and a description is given of difference schemes for the potential flow equation in quasi-linear form. Attention is also given to difference schemes for the potential flow equation in conservation form, the analysis of relaxation schemes by the time dependent analogy, the accelerated iterative method, and three-dimensional calculations.

Formal Integrals and Noether Operators of Nonlinear Hyperbolic Partial Differential Systems Admitting a Rich Set of Symmetries

NASA Astrophysics Data System (ADS)

Startsev, Sergey Ya.

2017-05-01

The paper is devoted to hyperbolic (generally speaking, non-Lagrangian and nonlinear) partial differential systems possessing a full set of differential operators that map any function of one independent variable into a symmetry of the corresponding system. We demonstrate that a system has the above property if and only if this system admits a full set of formal integrals (i.e., differential operators which map symmetries into integrals of the system). As a consequence, such systems possess both direct and inverse Noether operators (in the terminology of a work by B. Fuchssteiner and A.S. Fokas who have used these terms for operators that map cosymmetries into symmetries and perform transformations in the opposite direction). Systems admitting Noether operators are not exhausted by Euler-Lagrange systems and the systems with formal integrals. In particular, a hyperbolic system admits an inverse Noether operator if a differential substitution maps this system into a system possessing an inverse Noether operator.

Fractional vector calculus and fluid mechanics

NASA Astrophysics Data System (ADS)

Lazopoulos, Konstantinos A.; Lazopoulos, Anastasios K.

2017-04-01

Basic fluid mechanics equations are studied and revised under the prism of fractional continuum mechanics (FCM), a very promising research field that satisfies both experimental and theoretical demands. The geometry of the fractional differential has been clarified corrected and the geometry of the fractional tangent spaces of a manifold has been studied in Lazopoulos and Lazopoulos (Lazopoulos KA, Lazopoulos AK. Progr. Fract. Differ. Appl. 2016, 2, 85-104), providing the bases of the missing fractional differential geometry. Therefore, a lot can be contributed to fractional hydrodynamics: the basic fractional fluid equations (Navier Stokes, Euler and Bernoulli) are derived and fractional Darcy's flow in porous media is studied.

Numerical simulation of the nonlinear dynamics of harmonically driven Riesz-fractional extensions of the Fermi-Pasta-Ulam chains

NASA Astrophysics Data System (ADS)

Macías-Díaz, J. E.

2018-02-01

In this work, we introduce a spatially discrete model that is a modification of the well-known α-Fermi-Pasta-Ulam chain with damping. The system is perturbed at one end by a harmonic disturbance irradiating at a frequency in the forbidden band-gap of the classical regime, and a nonlocal coupling between the oscillators is considered using discrete Riesz fractional derivatives. We propose fully discrete expressions to approximate an energy functional of the system, and we use them to calculate the total energy of fractional chains over a relatively long period of time [Fract. Diff. Appl. 4 (2004) 153-162]. The approach is thoroughly tested in the case of local couplings against known qualitative results, including simulations of the process of nonlinear recurrence in the traditional chains of anharmonic oscillators. As an application, we provide evidence that the process of supratransmission is present in spatially discrete Fermi-Pasta-Ulam lattices with Riesz fractional derivatives in space. Moreover, we perform numerical experiments for small and large amplitudes of the harmonic disturbance. In either case, we establish the dependency of the critical amplitude at which supratransmission begins as a function of the driving frequency. Our results are in good agreement with the analytic predictions for the classical Fermi-Pasta-Ulam chain.

Jump Resonance in Fractional Order Circuits

NASA Astrophysics Data System (ADS)

Buscarino, Arturo; Caponetto, Riccardo; Famoso, Carlo; Fortuna, Luigi

The occurrence of an hysteretic loop in the frequency response of a driven nonlinear system is a phenomenon deeply investigated in nonlinear control theory. Such a phenomenon, which is linked to the multistable behavior of the system, is called jump resonance, since the magnitude of the frequency response is subjected to an abrupt jump up/down with respect to the increasing/decreasing of the frequency of the driving signal. In this paper, we aim at investigating fractional order nonlinear systems showing jump resonance, that is systems in which the order of the derivative is noninteger and their frequency response has a magnitude that is a multivalued function in a given range of frequencies. Furthermore, a strategy for designing fractional order systems showing jump resonance is presented along with the procedure to design and implement an analog circuit based on the approximation of the fractional order derivative. An extensive numerical analysis allows one to assess that the phenomenon is robust to the difference in the derivative order, enlightening the first example of a system with order lower than two which is able to demonstrate a jump resonance behavior.

Nonlinear differential system applied of a mechanical plan model of the automotives used for the nonlinear stability analysis

NASA Astrophysics Data System (ADS)

Simniceanu, Loreta; Mihaela, Bogdan; Otat, Victor; Trotea, Mario

2017-10-01

This paper proposes a plan mechanical model for the vehicles with two axles, taking into account the lateral deflection of the tire. For this mechanical model are determined two mathematical models under the nonlinear differential equations systems form without taking into account the action of the driver and taking into account. The analysis of driver-vehicle system consists in the mathematical description of vehicle dynamics, coupled with the possibilities and limits of the human factor. Description seeks to emphasize the significant influence of the driver in handling and stability analyzes of vehicles and vehicle-driver system stability until the advent of skidding. These mathematical models are seen as very useful tools to analyzing the vehicles stability. The paper analyzes the influence of some parameters of the vehicle on its behavior in terms of stability of dynamic systems.

Spline approximations for nonlinear hereditary control systems

NASA Technical Reports Server (NTRS)

Daniel, P. L.

1982-01-01

A sline-based approximation scheme is discussed for optimal control problems governed by nonlinear nonautonomous delay differential equations. The approximating framework reduces the original control problem to a sequence of optimization problems governed by ordinary differential equations. Convergence proofs, which appeal directly to dissipative-type estimates for the underlying nonlinear operator, are given and numerical findings are summarized.

Design of minimum multiplier fractional order differentiator based on lattice wave digital filter.

PubMed

Barsainya, Richa; Rawat, Tarun Kumar; Kumar, Manjeet

2017-01-01

In this paper, a novel design of fractional order differentiator (FOD) based on lattice wave digital filter (LWDF) is proposed which requires minimum number of multiplier for its structural realization. Firstly, the FOD design problem is formulated as an optimization problem using the transfer function of lattice wave digital filter. Then, three optimization algorithms, namely, genetic algorithm (GA), particle swarm optimization (PSO) and cuckoo search algorithm (CSA) are applied to determine the optimal LWDF coefficients. The realization of FOD using LWD structure increases the design accuracy, as only N number of coefficients are to be optimized for Nth order FOD. Finally, two design examples of 3rd and 5th order lattice wave digital fractional order differentiator (LWDFOD) are demonstrated to justify the design accuracy. The performance analysis of the proposed design is carried out based on magnitude response, absolute magnitude error (dB), root mean square (RMS) magnitude error, arithmetic complexity, convergence profile and computation time. Simulation results are attained to show the comparison of the proposed LWDFOD with the published works and it is observed that an improvement of 29% is obtained in the proposed design. The proposed LWDFOD approximates the ideal FOD and surpasses the existing ones reasonably well in mid and high frequency range, thereby making the proposed LWDFOD a promising technique for the design of digital FODs. Copyright © 2016 ISA. Published by Elsevier Ltd. All rights reserved.

Parametric Identification of Nonlinear Dynamical Systems

NASA Technical Reports Server (NTRS)

Feeny, Brian

2002-01-01

In this project, we looked at the application of harmonic balancing as a tool for identifying parameters (HBID) in a nonlinear dynamical systems with chaotic responses. The main idea is to balance the harmonics of periodic orbits extracted from measurements of each coordinate during a chaotic response. The periodic orbits are taken to be approximate solutions to the differential equations that model the system, the form of the differential equations being known, but with unknown parameters to be identified. Below we summarize the main points addressed in this work. The details of the work are attached as drafts of papers, and a thesis, in the appendix. Our study involved the following three parts: (1) Application of the harmonic balance to a simulation case in which the differential equation model has known form for its nonlinear terms, in contrast to a differential equation model which has either power series or interpolating functions to represent the nonlinear terms. We chose a pendulum, which has sinusoidal nonlinearities; (2) Application of the harmonic balance to an experimental system with known nonlinear forms. We chose a double pendulum, for which chaotic response were easily generated. Thus we confronted a two-degree-of-freedom system, which brought forth challenging issues; (3) A study of alternative reconstruction methods. The reconstruction of the phase space is necessary for the extraction of periodic orbits from the chaotic responses, which is needed in this work. Also, characterization of a nonlinear system is done in the reconstructed phase space. Such characterizations are needed to compare models with experiments. Finally, some nonlinear prediction methods can be applied in the reconstructed phase space. We developed two reconstruction methods that may be considered if the common method (method of delays) is not applicable.

Solving Fractional Programming Problems based on Swarm Intelligence

NASA Astrophysics Data System (ADS)

Raouf, Osama Abdel; Hezam, Ibrahim M.

2014-04-01

This paper presents a new approach to solve Fractional Programming Problems (FPPs) based on two different Swarm Intelligence (SI) algorithms. The two algorithms are: Particle Swarm Optimization, and Firefly Algorithm. The two algorithms are tested using several FPP benchmark examples and two selected industrial applications. The test aims to prove the capability of the SI algorithms to solve any type of FPPs. The solution results employing the SI algorithms are compared with a number of exact and metaheuristic solution methods used for handling FPPs. Swarm Intelligence can be denoted as an effective technique for solving linear or nonlinear, non-differentiable fractional objective functions. Problems with an optimal solution at a finite point and an unbounded constraint set, can be solved using the proposed approach. Numerical examples are given to show the feasibility, effectiveness, and robustness of the proposed algorithm. The results obtained using the two SI algorithms revealed the superiority of the proposed technique among others in computational time. A better accuracy was remarkably observed in the solution results of the industrial application problems.

Stable Numerical Approach for Fractional Delay Differential Equations

NASA Astrophysics Data System (ADS)

Singh, Harendra; Pandey, Rajesh K.; Baleanu, D.

2017-12-01

In this paper, we present a new stable numerical approach based on the operational matrix of integration of Jacobi polynomials for solving fractional delay differential equations (FDDEs). The operational matrix approach converts the FDDE into a system of linear equations, and hence the numerical solution is obtained by solving the linear system. The error analysis of the proposed method is also established. Further, a comparative study of the approximate solutions is provided for the test examples of the FDDE by varying the values of the parameters in the Jacobi polynomials. As in special case, the Jacobi polynomials reduce to the well-known polynomials such as (1) Legendre polynomial, (2) Chebyshev polynomial of second kind, (3) Chebyshev polynomial of third and (4) Chebyshev polynomial of fourth kind respectively. Maximum absolute error and root mean square error are calculated for the illustrated examples and presented in form of tables for the comparison purpose. Numerical stability of the presented method with respect to all four kind of polynomials are discussed. Further, the obtained numerical results are compared with some known methods from the literature and it is observed that obtained results from the proposed method is better than these methods.

Nonlinear Conservation Laws and Finite Volume Methods

NASA Astrophysics Data System (ADS)

Leveque, Randall J.

Introduction Software Notation Classification of Differential Equations Derivation of Conservation Laws The Euler Equations of Gas Dynamics Dissipative Fluxes Source Terms Radiative Transfer and Isothermal Equations Multi-dimensional Conservation Laws The Shock Tube Problem Mathematical Theory of Hyperbolic Systems Scalar Equations Linear Hyperbolic Systems Nonlinear Systems The Riemann Problem for the Euler Equations Numerical Methods in One Dimension Finite Difference Theory Finite Volume Methods Importance of Conservation Form - Incorrect Shock Speeds Numerical Flux Functions Godunov's Method Approximate Riemann Solvers High-Resolution Methods Other Approaches Boundary Conditions Source Terms and Fractional Steps Unsplit Methods Fractional Step Methods General Formulation of Fractional Step Methods Stiff Source Terms Quasi-stationary Flow and Gravity Multi-dimensional Problems Dimensional Splitting Multi-dimensional Finite Volume Methods Grids and Adaptive Refinement Computational Difficulties Low-Density Flows Discrete Shocks and Viscous Profiles Start-Up Errors Wall Heating Slow-Moving Shocks Grid Orientation Effects Grid-Aligned Shocks Magnetohydrodynamics The MHD Equations One-Dimensional MHD Solving the Riemann Problem Nonstrict Hyperbolicity Stiffness The Divergence of B Riemann Problems in Multi-dimensional MHD Staggered Grids The 8-Wave Riemann Solver Relativistic Hydrodynamics Conservation Laws in Spacetime The Continuity Equation The 4-Momentum of a Particle The Stress-Energy Tensor Finite Volume Methods Multi-dimensional Relativistic Flow Gravitation and General Relativity References

Mathematical nonlinear optics

NASA Astrophysics Data System (ADS)

McLaughlin, David W.

1995-08-01

The principal investigator, together with a post-doctoral fellows Tetsuji Ueda and Xiao Wang, several graduate students, and colleagues, has applied the modern mathematical theory of nonlinear waves to problems in nonlinear optics and to equations directly relevant to nonlinear optics. Projects included the interaction of laser light with nematic liquid crystals and chaotic, homoclinic, small dispersive, and random behavior of solutions of the nonlinear Schroedinger equation. In project 1, the extremely strong nonlinear response of a continuous wave laser beam in a nematic liquid crystal medium has produced striking undulation and filamentation of the laser beam which has been observed experimentally and explained theoretically. In project 2, qualitative properties of the nonlinear Schroedinger equation (which is the fundamental equation for nonlinear optics) have been identified and studied. These properties include optical shocking behavior in the limit of very small dispersion, chaotic and homoclinic behavior in discretizations of the partial differential equation, and random behavior.

Fractional Stochastic Field Theory

NASA Astrophysics Data System (ADS)

Honkonen, Juha

2018-02-01

Models describing evolution of physical, chemical, biological, social and financial processes are often formulated as differential equations with the understanding that they are large-scale equations for averages of quantities describing intrinsically random processes. Explicit account of randomness may lead to significant changes in the asymptotic behaviour (anomalous scaling) in such models especially in low spatial dimensions, which in many cases may be captured with the use of the renormalization group. Anomalous scaling and memory effects may also be introduced with the use of fractional derivatives and fractional noise. Construction of renormalized stochastic field theory with fractional derivatives and fractional noise in the underlying stochastic differential equations and master equations and the interplay between fluctuation-induced and built-in anomalous scaling behaviour is reviewed and discussed.

Fitting Nonlinear Ordinary Differential Equation Models with Random Effects and Unknown Initial Conditions Using the Stochastic Approximation Expectation-Maximization (SAEM) Algorithm.

PubMed

Chow, Sy-Miin; Lu, Zhaohua; Sherwood, Andrew; Zhu, Hongtu

2016-03-01

The past decade has evidenced the increased prevalence of irregularly spaced longitudinal data in social sciences. Clearly lacking, however, are modeling tools that allow researchers to fit dynamic models to irregularly spaced data, particularly data that show nonlinearity and heterogeneity in dynamical structures. We consider the issue of fitting multivariate nonlinear differential equation models with random effects and unknown initial conditions to irregularly spaced data. A stochastic approximation expectation-maximization algorithm is proposed and its performance is evaluated using a benchmark nonlinear dynamical systems model, namely, the Van der Pol oscillator equations. The empirical utility of the proposed technique is illustrated using a set of 24-h ambulatory cardiovascular data from 168 men and women. Pertinent methodological challenges and unresolved issues are discussed.

FITTING NONLINEAR ORDINARY DIFFERENTIAL EQUATION MODELS WITH RANDOM EFFECTS AND UNKNOWN INITIAL CONDITIONS USING THE STOCHASTIC APPROXIMATION EXPECTATION–MAXIMIZATION (SAEM) ALGORITHM

PubMed Central

Chow, Sy- Miin; Lu, Zhaohua; Zhu, Hongtu; Sherwood, Andrew

2014-01-01

The past decade has evidenced the increased prevalence of irregularly spaced longitudinal data in social sciences. Clearly lacking, however, are modeling tools that allow researchers to fit dynamic models to irregularly spaced data, particularly data that show nonlinearity and heterogeneity in dynamical structures. We consider the issue of fitting multivariate nonlinear differential equation models with random effects and unknown initial conditions to irregularly spaced data. A stochastic approximation expectation–maximization algorithm is proposed and its performance is evaluated using a benchmark nonlinear dynamical systems model, namely, the Van der Pol oscillator equations. The empirical utility of the proposed technique is illustrated using a set of 24-h ambulatory cardiovascular data from 168 men and women. Pertinent methodological challenges and unresolved issues are discussed. PMID:25416456

A Comparison of Two-Stage Approaches for Fitting Nonlinear Ordinary Differential Equation (ODE) Models with Mixed Effects

PubMed Central

Chow, Sy-Miin; Bendezú, Jason J.; Cole, Pamela M.; Ram, Nilam

2016-01-01

Several approaches currently exist for estimating the derivatives of observed data for model exploration purposes, including functional data analysis (FDA), generalized local linear approximation (GLLA), and generalized orthogonal local derivative approximation (GOLD). These derivative estimation procedures can be used in a two-stage process to fit mixed effects ordinary differential equation (ODE) models. While the performance and utility of these routines for estimating linear ODEs have been established, they have not yet been evaluated in the context of nonlinear ODEs with mixed effects. We compared properties of the GLLA and GOLD to an FDA-based two-stage approach denoted herein as functional ordinary differential equation with mixed effects (FODEmixed) in a Monte Carlo study using a nonlinear coupled oscillators model with mixed effects. Simulation results showed that overall, the FODEmixed outperformed both the GLLA and GOLD across all the embedding dimensions considered, but a novel use of a fourth-order GLLA approach combined with very high embedding dimensions yielded estimation results that almost paralleled those from the FODEmixed. We discuss the strengths and limitations of each approach and demonstrate how output from each stage of FODEmixed may be used to inform empirical modeling of young children’s self-regulation. PMID:27391255

Clinical Nonlinear Laser Imaging of Human Skin: A Review

PubMed Central

Pavone, Francesco Saverio

2014-01-01

Nonlinear optical microscopy has the potential of being used in vivo as a noninvasive imaging modality for both epidermal and dermal imaging. This paper reviews the capabilities of nonlinear microscopy as a noninvasive high-resolution tool for clinical skin inspection. In particular, we show that two-photon fluorescence microscopy can be used as a diagnostic tool for characterizing epidermal layers by means of a morphological examination. Additional functional information on the metabolic state of cells can be provided by measuring the fluorescence decay of NADH. This approach allows differentiating epidermal layers having different structural and cytological features and has the potential of diagnosing pathologies in a very early stage. Regarding therapy follow-up, we demonstrate that nonlinear microscopy could be successfully used for monitoring the effect of a treatment. In particular, combined two-photon fluorescence and second-harmonic generation microscopy were used in vivo for monitoring collagen remodeling after microablative fractional laser resurfacing and for quantitatively monitoring psoriasis on the basis of the morphology of epidermal cells and dermal papillae. We believe that the described microscopic modalities could find in the near future a stable place in a clinical dermatological setting for quantitative diagnostic purposes and as a monitoring method for various treatments. PMID:25250337

Bayesian parameter estimation for nonlinear modelling of biological pathways.

PubMed

Ghasemi, Omid; Lindsey, Merry L; Yang, Tianyi; Nguyen, Nguyen; Huang, Yufei; Jin, Yu-Fang

2011-01-01

The availability of temporal measurements on biological experiments has significantly promoted research areas in systems biology. To gain insight into the interaction and regulation of biological systems, mathematical frameworks such as ordinary differential equations have been widely applied to model biological pathways and interpret the temporal data. Hill equations are the preferred formats to represent the reaction rate in differential equation frameworks, due to their simple structures and their capabilities for easy fitting to saturated experimental measurements. However, Hill equations are highly nonlinearly parameterized functions, and parameters in these functions cannot be measured easily. Additionally, because of its high nonlinearity, adaptive parameter estimation algorithms developed for linear parameterized differential equations cannot be applied. Therefore, parameter estimation in nonlinearly parameterized differential equation models for biological pathways is both challenging and rewarding. In this study, we propose a Bayesian parameter estimation algorithm to estimate parameters in nonlinear mathematical models for biological pathways using time series data. We used the Runge-Kutta method to transform differential equations to difference equations assuming a known structure of the differential equations. This transformation allowed us to generate predictions dependent on previous states and to apply a Bayesian approach, namely, the Markov chain Monte Carlo (MCMC) method. We applied this approach to the biological pathways involved in the left ventricle (LV) response to myocardial infarction (MI) and verified our algorithm by estimating two parameters in a Hill equation embedded in the nonlinear model. We further evaluated our estimation performance with different parameter settings and signal to noise ratios. Our results demonstrated the effectiveness of the algorithm for both linearly and nonlinearly parameterized dynamic systems. Our proposed

Analysis of differential detergent fractions of an AtT-20 cellular homogenate using one- and two-dimensional capillary electrophoresis.

PubMed

Fazal, Md Abul; Palmer, Vanessa R; Dovichi, Norman J

2006-10-20

Differential detergent fractionation was used to sequentially extract cytosolic, membrane, nuclear, and cytoskeletal fractions from AtT-20 cells. Extracted components were denatured by sodium dodecyl sulfate (SDS) and then labeled with the fluorogenic reagent 3-(2-furoyl) quinoline-1-carboxaldehyde. Both capillary sieving electrophoresis (CSE) and micellar electrokinetic capillary chromatography (MECC) were used to separate labeled components by one-dimensional (1D) electrophoresis. Labeled components were also separated by two-dimensional (2D) capillary electrophoresis; CSE was employed in the first dimension and MECC in the second dimension. Roughly 150 fractions were transferred from the first to the second capillary for this comprehensive analysis in 2.5 h.

Numerical analysis for trajectory controllability of a coupled multi-order fractional delay differential system via the shifted Jacobi method

NASA Astrophysics Data System (ADS)

Priya, B. Ganesh; Muthukumar, P.

2018-02-01

This paper deals with the trajectory controllability for a class of multi-order fractional linear systems subject to a constant delay in state vector. The solution for the coupled fractional delay differential equation is established by the Mittag-Leffler function. The necessary and sufficient condition for the trajectory controllability is formulated and proved by the generalized Gronwall's inequality. The approximate trajectory for the proposed system is obtained through the shifted Jacobi operational matrix method. The numerical simulation of the approximate solution shows the theoretical results. Finally, some remarks and comments on the existing results of constrained controllability for the fractional dynamical system are also presented.

On the origins of generalized fractional calculus

NASA Astrophysics Data System (ADS)

Kiryakova, Virginia

2015-11-01

In Fractional Calculus (FC), as in the (classical) Calculus, the notions of derivatives and integrals (of first, second, etc. or arbitrary, incl. non-integer order) are basic and co-related. One of the most frequent approach in FC is to define first the Riemann-Liouville (R-L) integral of fractional order, and then by means of suitable integer-order differentiation operation applied over it (or under its sign) a fractional derivative is defined - in the R-L sense (or in Caputo sense). The first mentioned (R-L type) is closer to the theoretical studies in analysis, but has some shortages - from the point of view of interpretation of the initial conditions for Cauchy problems for fractional differential equations (stated also by means of fractional order derivatives/ integrals), and also for the analysts' confusion that such a derivative of a constant is not zero in general. The Caputo (C-) derivative, arising first in geophysical studies, helps to overcome these problems and to describe models of applied problems with physically consistent initial conditions. The operators of the Generalized Fractional Calculus - GFC (integrals and derivatives) are based on commuting m-tuple (m = 1, 2, 3, …) compositions of operators of the classical FC with power weights (the so-called Erdélyi-Kober operators), but represented in compact and explicit form by means of integral, integro-differential (R-L type) or differential-integral (C-type) operators, where the kernels are special functions of most general hypergeometric kind. The foundations of this theory are given in Kiryakova 18. In this survey we present the genesis of the definitions of the GFC - the generalized fractional integrals and derivatives (of fractional multi-order) of R-L type and Caputo type, analyze their properties and applications. Their special cases are all the known operators of classical FC, their generalizations introduced by other authors, the hyper-Bessel differential operators of higher integer

A differentiable reformulation for E-optimal design of experiments in nonlinear dynamic biosystems.

PubMed

Telen, Dries; Van Riet, Nick; Logist, Flip; Van Impe, Jan

2015-06-01

Informative experiments are highly valuable for estimating parameters in nonlinear dynamic bioprocesses. Techniques for optimal experiment design ensure the systematic design of such informative experiments. The E-criterion which can be used as objective function in optimal experiment design requires the maximization of the smallest eigenvalue of the Fisher information matrix. However, one problem with the minimal eigenvalue function is that it can be nondifferentiable. In addition, no closed form expression exists for the computation of eigenvalues of a matrix larger than a 4 by 4 one. As eigenvalues are normally computed with iterative methods, state-of-the-art optimal control solvers are not able to exploit automatic differentiation to compute the derivatives with respect to the decision variables. In the current paper a reformulation strategy from the field of convex optimization is suggested to circumvent these difficulties. This reformulation requires the inclusion of a matrix inequality constraint involving positive semidefiniteness. In this paper, this positive semidefiniteness constraint is imposed via Sylverster's criterion. As a result the maximization of the minimum eigenvalue function can be formulated in standard optimal control solvers through the addition of nonlinear constraints. The presented methodology is successfully illustrated with a case study from the field of predictive microbiology. Copyright © 2015. Published by Elsevier Inc.

Usefulness of the troponin-ejection fraction product to differentiate stress cardiomyopathy from ST-segment elevation myocardial infarction.

PubMed

Nascimento, Francisco O; Yang, Solomon; Larrauri-Reyes, Maiteder; Pineda, Andres M; Cornielle, Vertilio; Santana, Orlando; Heimowitz, Todd B; Stone, Gregg W; Beohar, Nirat

2014-02-01

The presentation of stress cardiomyopathy (SC) with nonobstructive coronary artery disease mimics that of ST-segment elevation myocardial infarction (STEMI) due to coronary occlusion. No single parameter has been successful in differentiating the 2 entities. We thus sought to develop a noninvasive clinical tool to discriminate between these 2 conditions. We retrospectively reviewed 59 consecutive cases of SC at our institution from July 2005 through June 2011 and compared those with 60 consecutives cases of angiographically confirmed STEMI treated with primary percutaneous coronary intervention in the same period. All patients underwent acute echocardiography, and the peak troponin I level was determined. The troponin-ejection fraction product (TEFP) was derived by multiplying the peak troponin I level and the echocardiographically derived left ventricular ejection fraction. Comparing the SC and STEMI groups, the mean left ventricular ejection fraction at the time of presentation was 30 ± 9% versus 44 ± 11%, respectively (p <0.001), and the peak troponin I was 7.6 ± 18 versus 102.2 ± 110.3 ng/dl, respectively (p <0.001). The mean TEFP was thus 182 ± 380 and 4,088 ± 4,244 for the SC and STEMI groups, respectively (p <0.001). Receiver operating characteristic curve analysis showed that a TEFP value ≥250 had a sensitivity of 95%, a specificity of 87%, a negative predictive value of 94%, a positive predictive value of 88%, and an overall accuracy of 91% to differentiate a true STEMI from SC (C-statistic 0.91 ± 0.02, p <0.001). In conclusion, for patients not undergoing emergent angiography, the TEFP may be used with high accuracy to differentiate SC with nonobstructive coronary artery disease from true STEMI due to coronary occlusion. Copyright © 2014 Elsevier Inc. All rights reserved.

Dynamical Approach Study of Spurious Steady-State Numerical Solutions of Nonlinear Differential Equations. Part 2; Global Asymptotic Behavior of Time Discretizations

NASA Technical Reports Server (NTRS)

Yee, H. C.; Sweby, P. K.

1995-01-01

The global asymptotic nonlinear behavior of 11 explicit and implicit time discretizations for four 2 x 2 systems of first-order autonomous nonlinear ordinary differential equations (ODEs) is analyzed. The objectives are to gain a basic understanding of the difference in the dynamics of numerics between the scalars and systems of nonlinear autonomous ODEs and to set a baseline global asymptotic solution behavior of these schemes for practical computations in computational fluid dynamics. We show how 'numerical' basins of attraction can complement the bifurcation diagrams in gaining more detailed global asymptotic behavior of time discretizations for nonlinear differential equations (DEs). We show how in the presence of spurious asymptotes the basins of the true stable steady states can be segmented by the basins of the spurious stable and unstable asymptotes. One major consequence of this phenomenon which is not commonly known is that this spurious behavior can result in a dramatic distortion and, in most cases, a dramatic shrinkage and segmentation of the basin of attraction of the true solution for finite time steps. Such distortion, shrinkage and segmentation of the numerical basins of attraction will occur regardless of the stability of the spurious asymptotes, and will occur for unconditionally stable implicit linear multistep methods. In other words, for the same (common) steady-state solution the associated basin of attraction of the DE might be very different from the discretized counterparts and the numerical basin of attraction can be very different from numerical method to numerical method. The results can be used as an explanation for possible causes of error, and slow convergence and nonconvergence of steady-state numerical solutions when using the time-dependent approach for nonlinear hyperbolic or parabolic PDEs.

Evidence for high-temperature fractionation of lithium isotopes during differentiation of the Moon

NASA Astrophysics Data System (ADS)

Day, James M. D.; Qiu, Lin; Ash, Richard D.; McDonough, William F.; Teng, Fang-Zhen; Rudnick, Roberta L.; Taylor, Lawrence A.

2016-06-01

Lithium isotope and abundance data are reported for Apollo 15 and 17 mare basalts and the LaPaz low-Ti mare basalt meteorites, along with lithium isotope data for carbonaceous, ordinary, and enstatite chondrites, and chondrules from the Allende CV3 meteorite. Apollo 15 low-Ti mare basalts have lower Li contents and lower δ7Li (3.8 ± 1.2‰; all uncertainties are 2 standard deviations) than Apollo 17 high-Ti mare basalts (δ7Li = 5.2 ± 1.2‰), with evolved LaPaz mare basalts having high Li contents, but similar low δ7Li (3.7 ± 0.5‰) to Apollo 15 mare basalts. In low-Ti mare basalt 15555, the highest concentrations of Li occur in late-stage tridymite (>20 ppm) and plagioclase (11 ± 3 ppm), with olivine (6.1 ± 3.8 ppm), pyroxene (4.2 ± 1.6 ppm), and ilmenite (0.8 ± 0.7 ppm) having lower Li concentrations. Values of δ7Li in low- and high-Ti mare basalt sources broadly correlate negatively with 18O/16O and positively with 56Fe/54Fe (low-Ti: δ7Li ≤4‰; δ56Fe ≤0.04‰; δ18O ≥5.7‰; high-Ti: δ7Li >6‰ δ56Fe >0.18‰ δ18O <5.4‰). Lithium does not appear to have acted as a volatile element during planetary formation, with subequal Li contents in mare basalts compared with terrestrial, martian, or vestan basaltic rocks. Observed Li isotopic fractionations in mare basalts can potentially be explained through large-degree, high-temperature igneous differentiation of their source regions. Progressive magma ocean crystallization led to enrichment in Li and δ7Li in late-stage liquids, probably as a consequence of preferential retention of 7Li and Li in the melt relative to crystallizing solids. Lithium isotopic fractionation has not been observed during extensive differentiation in terrestrial magmatic systems and may only be recognizable during extensive planetary magmatic differentiation under volatile-poor conditions, as expected for the lunar magma ocean. Our new analyses of chondrites show that they have δ7Li ranging between -2.5‰ and 4�

Proton density fat fraction (PDFF) MRI for differentiation of benign and malignant vertebral lesions.

PubMed

Schmeel, Frederic Carsten; Luetkens, Julian Alexander; Wagenhäuser, Peter Johannes; Meier-Schroers, Michael; Kuetting, Daniel Lloyd; Feißt, Andreas; Gieseke, Jürgen; Schmeel, Leonard Christopher; Träber, Frank; Schild, Hans Heinz; Kukuk, Guido Matthias

2018-06-01

To investigate whether proton density fat fraction (PDFF) measurements using a six-echo modified Dixon sequence can help to differentiate between benign and malignant vertebral bone marrow lesions. Sixty-six patients were prospectively enrolled in our study. In addition to conventional MRI at 3.0-Tesla including at least sagittal T2-weighted/spectral attenuated inversion recovery and T1-weighted sequences, all patients underwent a sagittal six-echo modified Dixon sequence of the spine. The mean PDFF was calculated using regions of interest and compared between vertebral lesions. A cut-off value of 6.40% in PDFF was determined by receiver operating characteristic curves and used to differentiate between malignant (< 6.40%) and benign (≥ 6.40%) vertebral lesions. There were 77 benign and 44 malignant lesions. The PDFF of malignant lesions was statistically significant lower in comparison with benign lesions (p < 0.001) and normal vertebral bone marrow (p < 0.001). The areas under the curves (AUC) were 0.97 for differentiating benign from malignant lesions (p < 0.001) and 0.95 for differentiating acute vertebral fractures from malignant lesions (p < 0.001). This yielded a diagnostic accuracy of 96% in the differentiation of both benign lesions and acute vertebral fractures from malignancy. PDFF derived from six-echo modified Dixon allows for differentiation between benign and malignant vertebral lesions with a high diagnostic accuracy. • Establishing a diagnosis of indeterminate vertebral lesions is a common clinical problem • Benign bone marrow processes may mimic the signal alterations observed in malignancy • PDFF differentiates between benign and malignant lesions with a high diagnostic accuracy • PDFF of non-neoplastic vertebral lesions is significantly higher than that of malignancy • PDFF from six-echo modified Dixon may help avoid potentially harmful bone biopsy.

Linearization of Conservative Nonlinear Oscillators

ERIC Educational Resources Information Center

Belendez, A.; Alvarez, M. L.; Fernandez, E.; Pascual, I.

2009-01-01

A linearization method of the nonlinear differential equation for conservative nonlinear oscillators is analysed and discussed. This scheme is based on the Chebyshev series expansion of the restoring force which allows us to obtain a frequency-amplitude relation which is valid not only for small but also for large amplitudes and, sometimes, for…

FAST TRACK COMMUNICATION Quasi self-adjoint nonlinear wave equations

NASA Astrophysics Data System (ADS)

Ibragimov, N. H.; Torrisi, M.; Tracinà, R.

2010-11-01

In this paper we generalize the classification of self-adjoint second-order linear partial differential equation to a family of nonlinear wave equations with two independent variables. We find a class of quasi self-adjoint nonlinear equations which includes the self-adjoint linear equations as a particular case. The property of a differential equation to be quasi self-adjoint is important, e.g. for constructing conservation laws associated with symmetries of the differential equation.

Second-order numerical methods for multi-term fractional differential equations: Smooth and non-smooth solutions

NASA Astrophysics Data System (ADS)

Zeng, Fanhai; Zhang, Zhongqiang; Karniadakis, George Em

2017-12-01

Starting with the asymptotic expansion of the error equation of the shifted Gr\\"{u}nwald--Letnikov formula, we derive a new modified weighted shifted Gr\\"{u}nwald--Letnikov (WSGL) formula by introducing appropriate correction terms. We then apply one special case of the modified WSGL formula to solve multi-term fractional ordinary and partial differential equations, and we prove the linear stability and second-order convergence for both smooth and non-smooth solutions. We show theoretically and numerically that numerical solutions up to certain accuracy can be obtained with only a few correction terms. Moreover, the correction terms can be tuned according to the fractional derivative orders without explicitly knowing the analytical solutions. Numerical simulations verify the theoretical results and demonstrate that the new formula leads to better performance compared to other known numerical approximations with similar resolution.

Zirconium and hafnium fractionation in differentiation of alkali carbonatite magmatic systems

NASA Astrophysics Data System (ADS)

Kogarko, L. N.

2016-05-01

Zirconium and hafnium are valuable strategic metals which are in high demand in industry. The Zr and Hf contents are elevated in the final products of magmatic differentiation of alkali carbonatite rocks in the Polar Siberia region (Guli Complex) and Ukraine (Chernigov Massif). Early pyroxene fractionation led to an increase in the Zr/Hf ratio in the evolution of the ultramafic-alkali magmatic system due to a higher distribution coefficient of Hf in pyroxene with respect to Zr. The Rayleigh equation was used to calculate a quantitative model of variation in the Zr/Hf ratio in the development of the Guli magmatic system. Alkali carbonatite rocks originated from rare element-rich mantle reservoirs, in particular, the metasomatized mantle. Carbonated mantle xenoliths are characterized by a high Zr/Hf ratio due to clinopyroxene development during metasomatic replacement of orthopyroxene by carbonate fluid melt.

Analysis of tristable energy harvesting system having fractional order viscoelastic material

DOE Office of Scientific and Technical Information (OSTI.GOV)

Oumbé Tékam, G. T.; Woafo, P.; Kitio Kwuimy, C. A.

2015-01-15

A particular attention is devoted to analyze the dynamics of a strongly nonlinear energy harvester having fractional order viscoelastic flexible material. The strong nonlinearity is obtained from the magnetic interaction between the end free of the flexible material and three equally spaced magnets. Periodic responses are computed using the KrylovBogoliubov averaging method, and the effects of fractional order damping on the output electric energy are analyzed. It is obtained that the harvested energy is enhanced for small order of the fractional derivative. Considering the order and strength of the fractional viscoelastic property as control parameter, the complexity of the systemmore » response is investigated through the Melnikov criteria for horseshoes chaos, which allows us to derive the mathematical expression of the boundary between intra-well motion and bifurcations appearance domain. We observe that the order and strength of the fractional viscoelastic property can be effectively used to control chaos in the system. The results are confirmed by the smooth and fractal shape of the basin of attraction as the order of derivative decreases. The bifurcation diagrams and the corresponding Lyapunov exponents are plotted to get insight into the nonlinear response of the system.« less

Fractional diffusion on bounded domains

DOE PAGES

Defterli, Ozlem; D'Elia, Marta; Du, Qiang; ...

2015-03-13

We found that the mathematically correct specification of a fractional differential equation on a bounded domain requires specification of appropriate boundary conditions, or their fractional analogue. In this paper we discuss the application of nonlocal diffusion theory to specify well-posed fractional diffusion equations on bounded domains.

The Ground Flash Fraction Retrieval Algorithm Employing Differential Evolution: Simulations and Applications

NASA Technical Reports Server (NTRS)

Koshak, William; Solakiewicz, Richard

2012-01-01

The ability to estimate the fraction of ground flashes in a set of flashes observed by a satellite lightning imager, such as the future GOES-R Geostationary Lightning Mapper (GLM), would likely improve operational and scientific applications (e.g., severe weather warnings, lightning nitrogen oxides studies, and global electric circuit analyses). A Bayesian inversion method, called the Ground Flash Fraction Retrieval Algorithm (GoFFRA), was recently developed for estimating the ground flash fraction. The method uses a constrained mixed exponential distribution model to describe a particular lightning optical measurement called the Maximum Group Area (MGA). To obtain the optimum model parameters (one of which is the desired ground flash fraction), a scalar function must be minimized. This minimization is difficult because of two problems: (1) Label Switching (LS), and (2) Parameter Identity Theft (PIT). The LS problem is well known in the literature on mixed exponential distributions, and the PIT problem was discovered in this study. Each problem occurs when one allows the numerical minimizer to freely roam through the parameter search space; this allows certain solution parameters to interchange roles which leads to fundamental ambiguities, and solution error. A major accomplishment of this study is that we have employed a state-of-the-art genetic-based global optimization algorithm called Differential Evolution (DE) that constrains the parameter search in such a way as to remove both the LS and PIT problems. To test the performance of the GoFFRA when DE is employed, we applied it to analyze simulated MGA datasets that we generated from known mixed exponential distributions. Moreover, we evaluated the GoFFRA/DE method by applying it to analyze actual MGAs derived from low-Earth orbiting lightning imaging sensor data; the actual MGA data were classified as either ground or cloud flash MGAs using National Lightning Detection Network[TM] (NLDN) data. Solution error

Nonlinear dissipative and dispersive electrostatic structures in unmagnetized relativistic electron-ion plasma with warm ions and trapped electrons

NASA Astrophysics Data System (ADS)

Masood, W.; Hamid, Naira; Ilyas, Iffat; Siddiq, M.

2017-06-01

In this paper, we have investigated electrostatic solitary and shock waves in an unmagnetized relativistic electron-ion (ei) plasma in the presence of warm ions and trapped electrons. In this regard, we have derived the trapped Korteweg-de Vries Burgers (TKdVB) equation using the small amplitude approximation method, which to the best of our knowledge has not been investigated in plasmas. Since the TKdVB equation involves fractional nonlinearity on account of trapped electrons, we have employed a smartly crafted extension of the tangent hyperbolic method and presented the solution of the TKdVB equation in this paper. The limiting cases of the TKdVB equation yield trapped Burgers (TB) and trapped Korteweg-de Vries (TKdV) equations. We have also presented the solutions of TB and TKdV equations. We have also explored how the plasma parameters affect the propagation characteristics of the nonlinear structures obtained for these modified nonlinear partial differential equations. We hope that the present work will open new vistas of research in the nonlinear plasma theory both in classical and quantum plasmas.

Finite elements of nonlinear continua.

NASA Technical Reports Server (NTRS)

Oden, J. T.

1972-01-01

The finite element method is extended to a broad class of practical nonlinear problems, treating both theory and applications from a general and unifying point of view. The thermomechanical principles of continuous media and the properties of the finite element method are outlined, and are brought together to produce discrete physical models of nonlinear continua. The mathematical properties of the models are analyzed, and the numerical solution of the equations governing the discrete models is examined. The application of the models to nonlinear problems in finite elasticity, viscoelasticity, heat conduction, and thermoviscoelasticity is discussed. Other specific topics include the topological properties of finite element models, applications to linear and nonlinear boundary value problems, convergence, continuum thermodynamics, finite elasticity, solutions to nonlinear partial differential equations, and discrete models of the nonlinear thermomechanical behavior of dissipative media.

Chirped femtosecond pulses in the higher-order nonlinear Schrödinger equation with non-Kerr nonlinear terms and cubic-quintic-septic nonlinearities

NASA Astrophysics Data System (ADS)

Triki, Houria; Biswas, Anjan; Milović, Daniela; Belić, Milivoj

2016-05-01

We consider a high-order nonlinear Schrödinger equation with competing cubic-quintic-septic nonlinearities, non-Kerr quintic nonlinearity, self-steepening, and self-frequency shift. The model describes the propagation of ultrashort (femtosecond) optical pulses in highly nonlinear optical fibers. A new ansatz is adopted to obtain nonlinear chirp associated with the propagating femtosecond soliton pulses. It is shown that the resultant elliptic equation of the problem is of high order, contains several new terms and is more general than the earlier reported results, thus providing a systematic way to find exact chirped soliton solutions of the septic model. Novel soliton solutions, including chirped bright, dark, kink and fractional-transform soliton solutions are obtained for special choices of parameters. Furthermore, we present the parameter domains in which these optical solitons exist. The nonlinear chirp associated with each of the solitonic solutions is also determined. It is shown that the chirping is proportional to the intensity of the wave and depends on higher-order nonlinearities. Of special interest is the soliton solution of the bright and dark type, determined for the general case when all coefficients in the equation have nonzero values. These results can be useful for possible chirped-soliton-based applications of highly nonlinear optical fiber systems.

A method based on the Jacobi tau approximation for solving multi-term time-space fractional partial differential equations

NASA Astrophysics Data System (ADS)

Bhrawy, A. H.; Zaky, M. A.

2015-01-01

In this paper, we propose and analyze an efficient operational formulation of spectral tau method for multi-term time-space fractional differential equation with Dirichlet boundary conditions. The shifted Jacobi operational matrices of Riemann-Liouville fractional integral, left-sided and right-sided Caputo fractional derivatives are presented. By using these operational matrices, we propose a shifted Jacobi tau method for both temporal and spatial discretizations, which allows us to present an efficient spectral method for solving such problem. Furthermore, the error is estimated and the proposed method has reasonable convergence rates in spatial and temporal discretizations. In addition, some known spectral tau approximations can be derived as special cases from our algorithm if we suitably choose the corresponding special cases of Jacobi parameters θ and ϑ. Finally, in order to demonstrate its accuracy, we compare our method with those reported in the literature.

Stochastic P-bifurcation and stochastic resonance in a noisy bistable fractional-order system

NASA Astrophysics Data System (ADS)

Yang, J. H.; Sanjuán, Miguel A. F.; Liu, H. G.; Litak, G.; Li, X.

2016-12-01

We investigate the stochastic response of a noisy bistable fractional-order system when the fractional-order lies in the interval (0, 2]. We focus mainly on the stochastic P-bifurcation and the phenomenon of the stochastic resonance. We compare the generalized Euler algorithm and the predictor-corrector approach which are commonly used for numerical calculations of fractional-order nonlinear equations. Based on the predictor-corrector approach, the stochastic P-bifurcation and the stochastic resonance are investigated. Both the fractional-order value and the noise intensity can induce an stochastic P-bifurcation. The fractional-order may lead the stationary probability density function to turn from a single-peak mode to a double-peak mode. However, the noise intensity may transform the stationary probability density function from a double-peak mode to a single-peak mode. The stochastic resonance is investigated thoroughly, according to the linear and the nonlinear response theory. In the linear response theory, the optimal stochastic resonance may occur when the value of the fractional-order is larger than one. In previous works, the fractional-order is usually limited to the interval (0, 1]. Moreover, the stochastic resonance at the subharmonic frequency and the superharmonic frequency are investigated respectively, by using the nonlinear response theory. When it occurs at the subharmonic frequency, the resonance may be strong and cannot be ignored. When it occurs at the superharmonic frequency, the resonance is weak. We believe that the results in this paper might be useful for the signal processing of nonlinear systems.

Differential Response of Human Nasal and Bronchial Epithelial Cells upon Exposure to Size-fractionated Dairy Dust

PubMed Central

Hawley, Brie; Schaeffer, Joshua; Poole, Jill A.; Dooley, Gregory P.; Reynolds, Stephen; Volckens, John

2015-01-01

Exposure to organic dusts is associated with increased respiratory morbidity and mortality in agricultural workers. Organic dusts in dairy farm environments are complex, polydisperse mixtures of toxic and immunogenic compounds. Previous toxicological studies focused primarily on exposures to the respirable size fraction, however, organic dusts in dairy farm environments are known to contain larger particles. Given the size distribution of dusts from dairy farm environments, the nasal and bronchial epithelia represent targets of agricultural dust exposures. In this study, well-differentiated normal human bronchial epithelial cells and human nasal epithelial cells were exposed to two different size fractions (PM10 and PM>10) of dairy parlor dust using a novel aerosol-to-cell exposure system. Levels of pro-inflammatory transcripts (IL-8, IL-6, and TNF-α) were measured two hr after exposure. Lactate dehydrogenase (LDH) release was also measured as an indicator of cytotoxicity. Cell exposure to dust was measured in each size fraction as a function of mass, endotoxin, and muramic acid levels. To our knowledge, this is the first study to evaluate the effects of distinct size fractions of agricultural dust on human airway epithelial cells. Our results suggest that both PM10 and PM>10 size fractions elicit a pro-inflammatory response in airway epithelial cells and that the entire inhalable size fraction needs to be considered when assessing potential risks from exposure to agricultural dusts. Further, data suggest that human bronchial cells respond differently to these dusts than human nasal cells and, therefore, the two cell types need to be considered separately in airway cell models of agricultural dust toxicity. PMID:25965193

Application of the enhanced homotopy perturbation method to solve the fractional-order Bagley-Torvik differential equation

NASA Astrophysics Data System (ADS)

Zolfaghari, M.; Ghaderi, R.; Sheikhol Eslami, A.; Ranjbar, A.; Hosseinnia, S. H.; Momani, S.; Sadati, J.

2009-10-01

The enhanced homotopy perturbation method (EHPM) is applied for finding improved approximate solutions of the well-known Bagley-Torvik equation for three different cases. The main characteristic of the EHPM is using a stabilized linear part, which guarantees the stability and convergence of the overall solution. The results are finally compared with the Adams-Bashforth-Moulton numerical method, the Adomian decomposition method (ADM) and the fractional differential transform method (FDTM) to verify the performance of the EHPM.

Nonlinear fishbone dynamics in spherical tokamaks

DOE Data Explorer

Wang, Feng [Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States); Dalian Univ Technol, Sch Phys & Optoelect Technol, Minist Educ, Key Lab Mat Modificat Laser Ion & Electron Beams, Dalian 116024, Peoples R China.; Fu, G.Y. [Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States); Institute for Fusion Theory and Simulation and Department of Physics Hangzhou, Zhejiang University, Hangzhou, 310027, People's Republic of China; Shen, Wei [Institute of Plasma Physics, Chinese Academy of Science, Hefei 230031, People's Republic of China

2017-01-01

Linear and nonlinear kinetic-MHD hybrid simulations have been carried out to investigate linear stability and nonlinear dynamics of beam-driven fishbone instability in spherical tokamak plasmas. Realistic NSTX parameters with finite toroidal rotation were used. The results show that the fishbone is driven by both trapped and passing particles. The instability drive of passing particles is comparable to that of trapped particles in the linear regime. The effects of rotation are destabilizing and a new region of instability appears at higher q min (>1.5) values, q min being the minimum of safety factor profile. In the nonlinear regime, the mode saturates due to flattening of beam ion distribution, and this persists after initial saturation while mode frequency chirps down in such a way that the resonant trapped particles move out radially and keep in resonance with the mode. Correspondingly, the flattening region of beam ion distribution expands radially outward. A substantial fraction of initially non-resonant trapped particles become resonant around the time of mode saturation and keep in resonance with the mode as frequency chirps down. On the other hand, the fraction of resonant passing particles is significantly smaller than that of trapped particles. Our analysis shows that trapped particles provide the main drive to the mode in the nonlinear regime.

Soliton structure versus singularity analysis: Third-order completely intergrable nonlinear differential equations in 1 + 1-dimensions

NASA Astrophysics Data System (ADS)

Fuchssteiner, Benno; Carillo, Sandra

1989-01-01

Bäcklund transformations between all known completely integrable third-order differential equations in (1 + 1)-dimensions are established and the corresponding transformations formulas for their hereditary operators and Hamiltonian formulations are exhibited. Some of these Bäcklund transformations are not injective; therefore additional non-commutative symmetry groups are found for some equations. These non-commutative symmetry groups are classified as having a semisimple part isomorphic to the affine algebra A(1)1. New completely integrable third-order integro-differential equations, some depending explicitly on x, are given. These new equations give rise to nonin equation. Connections between the singularity equations (from the Painlevé analysis) and the nonlinear equations for interacting solitons are established. A common approach to singularity analysis and soliton structure is introduced. The Painlevé analysis is modified in such a sense that it carries over directly and without difficulty to the time evolution of singularity manifolds of equations like the sine-Gordon and nonlinear Schrödinger equation. A method to recover the Painlevé series from its constant level term is exhibit. The soliton-singularity transform is recognized to be connected to the Möbius group. This gives rise to a Darboux-like result for the spectral properties of the recursion operator. These connections are used in order to explain why poles of soliton equations move like trajectories of interacting solitons. Furthermore it is explicitly computed how solitons of singularity equations behave under the effect of this soliton-singularity transform. This then leads to the result that only for scaling degrees α = -1 and α = -2 the usual Painlevé analysis can be carried out. A new invariance principle, connected to kernels of differential operators is discovered. This new invariance, for example, connects the explicit solutions of the Liouville equation with the Miura transform

Peakompactons: Peaked compact nonlinear waves

DOE PAGES

Christov, Ivan C.; Kress, Tyler; Saxena, Avadh

2017-04-20

This paper is meant as an accessible introduction to/tutorial on the analytical construction and numerical simulation of a class of nonstandard solitary waves termed peakompactons. We present that these peaked compactly supported waves arise as solutions to nonlinear evolution equations from a hierarchy of nonlinearly dispersive Korteweg–de Vries-type models. Peakompactons, like the now-well-known compactons and unlike the soliton solutions of the Korteweg–de Vries equation, have finite support, i.e., they are of finite wavelength. However, unlike compactons, peakompactons are also peaked, i.e., a higher spatial derivative suffers a jump discontinuity at the wave’s crest. Here, we construct such solutions exactly bymore » reducing the governing partial differential equation to a nonlinear ordinary differential equation and employing a phase-plane analysis. Lastly, a simple, but reliable, finite-difference scheme is also designed and tested for the simulation of collisions of peakompactons. In addition to the peakompacton class of solutions, the general physical features of the so-called K #(n,m) hierarchy of nonlinearly dispersive Korteweg–de Vries-type models are discussed as well.« less

Finite element analysis of steady and transiently moving/rolling nonlinear viscoelastic structure. Part 1: Theory

NASA Technical Reports Server (NTRS)

Padovan, Joe

1986-01-01

In a three part series of papers, a generalized finite element analysis scheme is developed to handle the steady and transient response of moving/rolling nonlinear viscoelastic structure. This paper considers the development of the moving/rolling element strategy, including the effects of large deformation kinematics and viscoelasticity modelled by fractional integro-differential operators. To improve the solution strategy, a special hierarchical constraint procedure is developed for the case of steady rolling/translating as well as a transient scheme involving the use of a Grunwaldian representation of the fractional operator. In the second and third parts of the paper, 3-D extensions are developed along with transient contact strategies enabling the handling of impacts with obstructions. Overall, the various developments are benchmarked via comprehensive 2- and 3-D simulations. These are correlated with experimental data to define modelling capabilities.

Investigation of a Nonlinear Control System

NASA Technical Reports Server (NTRS)

Flugge-Lotz, I; Taylor, C F; Lindberg, H E

1958-01-01

A discontinuous variation of coefficients of the differential equation describing the linear control system before nonlinear elements are added is studied in detail. The nonlinear feedback is applied to a second-order system. Simulation techniques are used to study performance of the nonlinear control system and to compare it with the linear system for a wide variety of inputs. A detailed quantitative study of the influence of relay delays and of a transport delay is presented.

Non-linear feedback control of the p53 protein-mdm2 inhibitor system using the derivative-free non-linear Kalman filter.

PubMed

Rigatos, Gerasimos G

2016-06-01

It is proven that the model of the p53-mdm2 protein synthesis loop is a differentially flat one and using a diffeomorphism (change of state variables) that is proposed by differential flatness theory it is shown that the protein synthesis model can be transformed into the canonical (Brunovsky) form. This enables the design of a feedback control law that maintains the concentration of the p53 protein at the desirable levels. To estimate the non-measurable elements of the state vector describing the p53-mdm2 system dynamics, the derivative-free non-linear Kalman filter is used. Moreover, to compensate for modelling uncertainties and external disturbances that affect the p53-mdm2 system, the derivative-free non-linear Kalman filter is re-designed as a disturbance observer. The derivative-free non-linear Kalman filter consists of the Kalman filter recursion applied on the linearised equivalent of the protein synthesis model together with an inverse transformation based on differential flatness theory that enables to retrieve estimates for the state variables of the initial non-linear model. The proposed non-linear feedback control and perturbations compensation method for the p53-mdm2 system can result in more efficient chemotherapy schemes where the infusion of medication will be better administered.

Well-conditioned fractional collocation methods using fractional Birkhoff interpolation basis

NASA Astrophysics Data System (ADS)

Jiao, Yujian; Wang, Li-Lian; Huang, Can

2016-01-01

The purpose of this paper is twofold. Firstly, we provide explicit and compact formulas for computing both Caputo and (modified) Riemann-Liouville (RL) fractional pseudospectral differentiation matrices (F-PSDMs) of any order at general Jacobi-Gauss-Lobatto (JGL) points. We show that in the Caputo case, it suffices to compute F-PSDM of order μ ∈ (0 , 1) to compute that of any order k + μ with integer k ≥ 0, while in the modified RL case, it is only necessary to evaluate a fractional integral matrix of order μ ∈ (0 , 1). Secondly, we introduce suitable fractional JGL Birkhoff interpolation problems leading to new interpolation polynomial basis functions with remarkable properties: (i) the matrix generated from the new basis yields the exact inverse of F-PSDM at "interior" JGL points; (ii) the matrix of the highest fractional derivative in a collocation scheme under the new basis is diagonal; and (iii) the resulted linear system is well-conditioned in the Caputo case, while in the modified RL case, the eigenvalues of the coefficient matrix are highly concentrated. In both cases, the linear systems of the collocation schemes using the new basis can be solved by an iterative solver within a few iterations. Notably, the inverse can be computed in a very stable manner, so this offers optimal preconditioners for usual fractional collocation methods for fractional differential equations (FDEs). It is also noteworthy that the choice of certain special JGL points with parameters related to the order of the equations can ease the implementation. We highlight that the use of the Bateman's fractional integral formulas and fast transforms between Jacobi polynomials with different parameters, is essential for our algorithm development.

On the removal of boundary errors caused by Runge-Kutta integration of non-linear partial differential equations

NASA Technical Reports Server (NTRS)

Abarbanel, Saul; Gottlieb, David; Carpenter, Mark H.

1994-01-01

It has been previously shown that the temporal integration of hyperbolic partial differential equations (PDE's) may, because of boundary conditions, lead to deterioration of accuracy of the solution. A procedure for removal of this error in the linear case has been established previously. In the present paper we consider hyperbolic (PDE's) (linear and non-linear) whose boundary treatment is done via the SAT-procedure. A methodology is present for recovery of the full order of accuracy, and has been applied to the case of a 4th order explicit finite difference scheme.

Nanoparticles and nonlinear thermal radiation properties in the rheology of polymeric material

NASA Astrophysics Data System (ADS)

Awais, M.; Hayat, T.; Muqaddass, N.; Ali, A.; Aqsa; Awan, Saeed Ehsan

2018-03-01

The present analysis is related to the dynamics of polymeric liquids (Oldroyd-B model) with the presence of nanoparticles. The rheological system is considered under the application of nonlinear thermal radiations. Energy and concentration equations are presented when thermophoresis and Brownian motion effects are present. Bidirectional form of stretching is considered to interpret the three-dimensional flow dynamics of polymeric liquid. Making use of the similarity transformations, problem is reduced into ordinary differential system which is approximated by using HAM. Influence of physical parameters including Deborah number, thermophoresis and Brownian motion on velocity, temperature and mass fraction expressions are plotted and analyzed. Numerical values for local Sherwood and Nusselt numbers are presented and discussed.

Spectral decomposition of nonlinear systems with memory

NASA Astrophysics Data System (ADS)

Svenkeson, Adam; Glaz, Bryan; Stanton, Samuel; West, Bruce J.

2016-02-01

We present an alternative approach to the analysis of nonlinear systems with long-term memory that is based on the Koopman operator and a Lévy transformation in time. Memory effects are considered to be the result of interactions between a system and its surrounding environment. The analysis leads to the decomposition of a nonlinear system with memory into modes whose temporal behavior is anomalous and lacks a characteristic scale. On average, the time evolution of a mode follows a Mittag-Leffler function, and the system can be described using the fractional calculus. The general theory is demonstrated on the fractional linear harmonic oscillator and the fractional nonlinear logistic equation. When analyzing data from an ill-defined (black-box) system, the spectral decomposition in terms of Mittag-Leffler functions that we propose may uncover inherent memory effects through identification of a small set of dynamically relevant structures that would otherwise be obscured by conventional spectral methods. Consequently, the theoretical concepts we present may be useful for developing more general methods for numerical modeling that are able to determine whether observables of a dynamical system are better represented by memoryless operators, or operators with long-term memory in time, when model details are unknown.

Correcting the initialization of models with fractional derivatives via history-dependent conditions

NASA Astrophysics Data System (ADS)

Du, Maolin; Wang, Zaihua

2016-04-01

Fractional differential equations are more and more used in modeling memory (history-dependent, non-local, or hereditary) phenomena. Conventional initial values of fractional differential equations are defined at a point, while recent works define initial conditions over histories. We prove that the conventional initialization of fractional differential equations with a Riemann-Liouville derivative is wrong with a simple counter-example. The initial values were assumed to be arbitrarily given for a typical fractional differential equation, but we find one of these values can only be zero. We show that fractional differential equations are of infinite dimensions, and the initial conditions, initial histories, are defined as functions over intervals. We obtain the equivalent integral equation for Caputo case. With a simple fractional model of materials, we illustrate that the recovery behavior is correct with the initial creep history, but is wrong with initial values at the starting point of the recovery. We demonstrate the application of initial history by solving a forced fractional Lorenz system numerically.

Nonlinear Acoustic Propagation into the Seafloor

NASA Astrophysics Data System (ADS)

McDonald, B. Edward

2006-05-01

Explosions near the seafloor result in shock waves entering a much more complicated medium than water or air. Nonlinearities may be increased by two processes inherent to granular media: (1) a poroelastic nonlinearity comparable to the addition of bubbles to water, and (2) the Hertz force resulting from elastic deformation of grains, proportional to the Youngs modulus of the grains times the strain rate to the power 3/2. These two types of nonlinearity for shock propagation into the seafloor are investigated using a variant of the NPE model. The traditional Taylor series expansion of the equation of state (pressure as a function of density) is not appropriate to the Hertz force in the limit of small strain. We present a simple nonlinear wave equation model for compressional waves in marine sediments that retains the Hertz force explicitly with overdensity to the power 3/2. Numerical results for shock propagation are compared with similarity solutions for quadratic nonlinearity and for the fractional nonlinearity of the Hertz force.

Discrete fractional solutions of a Legendre equation

NASA Astrophysics Data System (ADS)

Yılmazer, Resat

2018-01-01

One of the most popular research interests of science and engineering is the fractional calculus theory in recent times. Discrete fractional calculus has also an important position in fractional calculus. In this work, we acquire new discrete fractional solutions of the homogeneous and non homogeneous Legendre differential equation by using discrete fractional nabla operator.

On the fractional Eulerian numbers and equivalence of maps with long term power-law memory (integral Volterra equations of the second kind) to Grünvald-Letnikov fractional difference (differential) equations.

PubMed

Edelman, Mark

2015-07-01

In this paper, we consider a simple general form of a deterministic system with power-law memory whose state can be described by one variable and evolution by a generating function. A new value of the system's variable is a total (a convolution) of the generating functions of all previous values of the variable with weights, which are powers of the time passed. In discrete cases, these systems can be described by difference equations in which a fractional difference on the left hand side is equal to a total (also a convolution) of the generating functions of all previous values of the system's variable with the fractional Eulerian number weights on the right hand side. In the continuous limit, the considered systems can be described by the Grünvald-Letnikov fractional differential equations, which are equivalent to the Volterra integral equations of the second kind. New properties of the fractional Eulerian numbers and possible applications of the results are discussed.

Bayesian analysis of non-linear differential equation models with application to a gut microbial ecosystem.

PubMed

Lawson, Daniel J; Holtrop, Grietje; Flint, Harry

2011-07-01

Process models specified by non-linear dynamic differential equations contain many parameters, which often must be inferred from a limited amount of data. We discuss a hierarchical Bayesian approach combining data from multiple related experiments in a meaningful way, which permits more powerful inference than treating each experiment as independent. The approach is illustrated with a simulation study and example data from experiments replicating the aspects of the human gut microbial ecosystem. A predictive model is obtained that contains prediction uncertainty caused by uncertainty in the parameters, and we extend the model to capture situations of interest that cannot easily be studied experimentally. Copyright © 2011 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

A Nonlinear Modal Aeroelastic Solver for FUN3D

NASA Technical Reports Server (NTRS)

Goldman, Benjamin D.; Bartels, Robert E.; Biedron, Robert T.; Scott, Robert C.

2016-01-01

A nonlinear structural solver has been implemented internally within the NASA FUN3D computational fluid dynamics code, allowing for some new aeroelastic capabilities. Using a modal representation of the structure, a set of differential or differential-algebraic equations are derived for general thin structures with geometric nonlinearities. ODEPACK and LAPACK routines are linked with FUN3D, and the nonlinear equations are solved at each CFD time step. The existing predictor-corrector method is retained, whereby the structural solution is updated after mesh deformation. The nonlinear solver is validated using a test case for a flexible aeroshell at transonic, supersonic, and hypersonic flow conditions. Agreement with linear theory is seen for the static aeroelastic solutions at relatively low dynamic pressures, but structural nonlinearities limit deformation amplitudes at high dynamic pressures. No flutter was found at any of the tested trajectory points, though LCO may be possible in the transonic regime.

Period of vibration of axially vibrating truly nonlinear rod

NASA Astrophysics Data System (ADS)

Cveticanin, L.

2016-07-01

In this paper the axial vibration of a muscle whose fibers are parallel to the direction of muscle compression is investigated. The model is a clamped-free rod with a strongly nonlinear elastic property. Axial vibration is described by a nonlinear partial differential equation. A solution of the equation is constructed for special initial conditions by using the method of separation of variables. The partial differential equation is separated into two uncoupled strongly nonlinear second order differential equations. Both equations, with displacement function and with time function are exactly determined. Exact solutions are given in the form of inverse incomplete and inverse complete Beta function. Using boundary and initial conditions, the frequency of vibration is obtained. It has to be mentioned that the determined frequency represents the exact analytic description for the axially vibrating truly nonlinear clamped-free rod. The procedure suggested in this paper is applied for calculation of the frequency of the longissimus dorsi muscle of a cow. The influence of elasticity order and elasticity coefficient on the frequency property is tested.

Bound-Electron Nonlinearity Beyond the Ionization Threshold.

PubMed

Wahlstrand, J K; Zahedpour, S; Bahl, A; Kolesik, M; Milchberg, H M

2018-05-04

We present absolute space- and time-resolved measurements of the ultrafast laser-driven nonlinear polarizability in argon, krypton, xenon, nitrogen, and oxygen up to ionization fractions of a few percent. These measurements enable determination of the strongly nonperturbative bound-electron nonlinear polarizability well beyond the ionization threshold, where it is found to remain approximately quadratic in the laser field, a result normally expected at much lower intensities where perturbation theory applies.

Bound-Electron Nonlinearity Beyond the Ionization Threshold

NASA Astrophysics Data System (ADS)

Wahlstrand, J. K.; Zahedpour, S.; Bahl, A.; Kolesik, M.; Milchberg, H. M.

2018-05-01

We present absolute space- and time-resolved measurements of the ultrafast laser-driven nonlinear polarizability in argon, krypton, xenon, nitrogen, and oxygen up to ionization fractions of a few percent. These measurements enable determination of the strongly nonperturbative bound-electron nonlinear polarizability well beyond the ionization threshold, where it is found to remain approximately quadratic in the laser field, a result normally expected at much lower intensities where perturbation theory applies.

Rotating Flow of Magnetite-Water Nanofluid over a Stretching Surface Inspired by Non-Linear Thermal Radiation.

PubMed

Mustafa, M; Mushtaq, A; Hayat, T; Alsaedi, A

2016-01-01

Present study explores the MHD three-dimensional rotating flow and heat transfer of ferrofluid induced by a radiative surface. The base fluid is considered as water with magnetite-Fe3O4 nanoparticles. Novel concept of non-linear radiative heat flux is considered which produces a non-linear energy equation in temperature field. Conventional transformations are employed to obtain the self-similar form of the governing differential system. The arising system involves an interesting temperature ratio parameter which is an indicator of small/large temperature differences in the flow. Numerical simulations with high precision are determined by well-known shooting approach. Both uniform stretching and rotation have significant impact on the solutions. The variation in velocity components with the nanoparticle volume fraction is non-monotonic. Local Nusselt number in Fe3O4-water ferrofluid is larger in comparison to the pure fluid even at low particle concentration.

Learning automata-based solutions to the nonlinear fractional knapsack problem with applications to optimal resource allocation.

PubMed

Granmo, Ole-Christoffer; Oommen, B John; Myrer, Svein Arild; Olsen, Morten Goodwin

2007-02-01

This paper considers the nonlinear fractional knapsack problem and demonstrates how its solution can be effectively applied to two resource allocation problems dealing with the World Wide Web. The novel solution involves a "team" of deterministic learning automata (LA). The first real-life problem relates to resource allocation in web monitoring so as to "optimize" information discovery when the polling capacity is constrained. The disadvantages of the currently reported solutions are explained in this paper. The second problem concerns allocating limited sampling resources in a "real-time" manner with the purpose of estimating multiple binomial proportions. This is the scenario encountered when the user has to evaluate multiple web sites by accessing a limited number of web pages, and the proportions of interest are the fraction of each web site that is successfully validated by an HTML validator. Using the general LA paradigm to tackle both of the real-life problems, the proposed scheme improves a current solution in an online manner through a series of informed guesses that move toward the optimal solution. At the heart of the scheme, a team of deterministic LA performs a controlled random walk on a discretized solution space. Comprehensive experimental results demonstrate that the discretization resolution determines the precision of the scheme, and that for a given precision, the current solution (to both problems) is consistently improved until a nearly optimal solution is found--even for switching environments. Thus, the scheme, while being novel to the entire field of LA, also efficiently handles a class of resource allocation problems previously not addressed in the literature.

Nonlinear fishbone dynamics in spherical tokamaks

DOE PAGES

Wang, Feng; Fu, G. Y.; Shen, Wei

2016-11-22

Linear and nonlinear kinetic-MHD hybrid simulations have been carried out to investigate linear stability and nonlinear dynamics of beam-driven fishbone instability in spherical tokamak plasmas. Realistic NSTX parameters with finite toroidal rotation were used. Our results show that the fishbone is driven by both trapped and passing particles. The instability drive of passing particles is comparable to that of trapped particles in the linear regime. The effects of rotation are destabilizing and a new region of instability appears at higher q min (>1.5) values, q min being the minimum of safety factor profile. In the nonlinear regime, the mode saturatesmore » due to flattening of beam ion distribution, and this persists after initial saturation while mode frequency chirps down in such a way that the resonant trapped particles move out radially and keep in resonance with the mode. Correspondingly, the flattening region of beam ion distribution expands radially outward. Furthermore, a substantial fraction of initially non-resonant trapped particles become resonant around the time of mode saturation and keep in resonance with the mode as frequency chirps down. On the other hand, the fraction of resonant passing particles is significantly smaller than that of trapped particles. Finally, our analysis shows that trapped particles provide the main drive to the mode in the nonlinear regime.« less

Numerical analysis for the fractional diffusion and fractional Buckmaster equation by the two-step Laplace Adam-Bashforth method

NASA Astrophysics Data System (ADS)

Jain, Sonal

2018-01-01

In this paper, we aim to use the alternative numerical scheme given by Gnitchogna and Atangana for solving partial differential equations with integer and non-integer differential operators. We applied this method to fractional diffusion model and fractional Buckmaster models with non-local fading memory. The method yields a powerful numerical algorithm for fractional order derivative to implement. Also we present in detail the stability analysis of the numerical method for solving the diffusion equation. This proof shows that this method is very stable and also converges very quickly to exact solution and finally some numerical simulation is presented.

A chaotic secure communication scheme using fractional chaotic systems based on an extended fractional Kalman filter

NASA Astrophysics Data System (ADS)

Kiani-B, Arman; Fallahi, Kia; Pariz, Naser; Leung, Henry

2009-03-01

In recent years chaotic secure communication and chaos synchronization have received ever increasing attention. In this paper, for the first time, a fractional chaotic communication method using an extended fractional Kalman filter is presented. The chaotic synchronization is implemented by the EFKF design in the presence of channel additive noise and processing noise. Encoding chaotic communication achieves a satisfactory, typical secure communication scheme. In the proposed system, security is enhanced based on spreading the signal in frequency and encrypting it in time domain. In this paper, the main advantages of using fractional order systems, increasing nonlinearity and spreading the power spectrum are highlighted. To illustrate the effectiveness of the proposed scheme, a numerical example based on the fractional Lorenz dynamical system is presented and the results are compared to the integer Lorenz system.

Nonlinear Problems in Fluid Dynamics and Inverse Scattering

DTIC Science & Technology

1993-05-31

nonlinear Kadomtsev - Petviashvili (KP) equations , have solutions which will become infinite in finite time. This phenomenon is sometimes referred to as...40 (November 1992). 4 7. Wave Collapse and Instability of Solitary Waves of a Generalized Nonlinear Kaoiomtsev- Petviashvili Equation , X.P. Wang, M.J...words) The inverse scattering of a class of differential-difference equations and multidimensional operators has been constructed. Solutions of nonlinear

The Vertical Linear Fractional Initialization Problem

NASA Technical Reports Server (NTRS)

Lorenzo, Carl F.; Hartley, Tom T.

1999-01-01

This paper presents a solution to the initialization problem for a system of linear fractional-order differential equations. The scalar problem is considered first, and solutions are obtained both generally and for a specific initialization. Next the vector fractional order differential equation is considered. In this case, the solution is obtained in the form of matrix F-functions. Some control implications of the vector case are discussed. The suggested method of problem solution is shown via an example.

Subcellular fractionation by differential and zonal centrifugation of aerobically grown glucose-de-repressed Saccharomyces carlsbergensis

PubMed Central

Cartledge, T. G.; Lloyd, D.

1972-01-01

1. Homogenates were prepared from sphaeroplasts of aerobically grown glucose-de-repressed Saccharomyces carlsbergensis and the distributions of marker enzymes were investigated after differential centrifugation. Cytochrome c oxidase and cytochrome c were sedimented almost completely at 105g-min, and this fraction also contained 37% of the catalase, 27% of the acid p-nitrophenyl phosphatase, 53 and 54% respectively of the NADH– and NADPH–cytochrome c oxidoreductases. 2. Zonal centrifugation indicated complex density distributions of the sedimentable portions of these enzymes and of adenosine triphosphatases and suggested the presence of two mitochondrial populations, as well as a bimodal distribution of peroxisomes and heterogeneity of the acid p-nitrophenyl phosphatase-containing particles. 3. Several different adenosine triphosphatases were distinguished in a post-mitochondrial supernatant that contained no mitochondrial fragments; these enzymes varied in their sensitivities to oligomycin and ouabain and their distributions were different from those of pyrophosphatase, adenosine phosphatase and adenosine pyrophosphatase. 4. The distribution of NADPH–cytochrome c oxidoreductase demonstrated that it cannot be used in S. carlsbergensis as a specific marker enzyme for the microsomal fraction. Glucose 6-phosphatase, inosine pyrophosphatase, cytochrome P-450 and five other enzymes frequently assigned to microsomal fractions of mammalian origin were not detected in yeast under these growth conditions. ImagesPLATE 2PLATE 1 (cont.)PLATE 1PLATE 2 (cont.) PMID:4400904

Isotopomers as a method for differentiating between bacterial and fungal production of nitrous oxide

NASA Astrophysics Data System (ADS)

Sutka, R. L.; Adams, G.; Ostrom, N.; Ostrom, P.

2007-12-01

In order to study the importance of fungi to nitrous oxide (N2O) production in the environment it is critical to have a non-intrusive method for differentiating between fungal and bacterial N2O production. Site preference (SP), the difference in d15N between the central and outer N atoms in N2O, has been used to differentiate between bacterial nitrification and denitrification. In this study we compare the SP, d15N and d18O of N2O produced by the two best-studied fungal denitrifiers, Fusarium oxysporum and Cylindrocarpon tonkinense, to data from our previous bacterial studies. Both d18O and SP values remained fairly constant during the course of nitrite reduction which likely reflects isotopic exchange with water in the case of d18O and conservative behavior in SP that has been observed previously (Sutka et al., 2006). We observed a wide range of fractionation factors for fungal denitrification, -74.7 to -6.6 ‰, and non-linear behavior indicating that fractionation was controlled by more than one step. We interpret the small degree of fractionation as reflecting fractionation during diffusion and the more negative values as being controlled by enzymatic fractionation. Data from this and our previous study of bacterial production (Sutka et al., 2006) reveals that N2O produced via nitrification by fungi can be differentiated from N2O produced by bacterial denitrification primarily on the basis of d18O. The site preference of N2O produced by F. oxysporum and C. tonkinense was 37.1 ± 2.5 ‰ and 36.9 ± 2.8 ‰, respectively. These results indicate that isotopomers can be used as a basis for differentiating bacterial and fungal denitrification. Our work further reveals the role that fungal and bacterial nitric oxide reductases have in determining site preference during N2O production.

The role of nonlinear critical layers in boundary layer transition

NASA Technical Reports Server (NTRS)

Goldstein, M.E.

1995-01-01

Asymptotic methods are used to describe the nonlinear self-interaction between pairs of oblique instability modes that eventually develops when initially linear spatially growing instability waves evolve downstream in nominally two-dimensional laminar boundary layers. The first nonlinear reaction takes place locally within a so-called 'critical layer', with the flow outside this layer consisting of a locally parallel mean flow plus a pair of oblique instability waves - which may or may not be accompanied by an associated plane wave. The amplitudes of these waves, which are completely determined by nonlinear effects within the critical layer, satisfy either a single integro-differential equation or a pair of integro-differential equations with quadratic to quartic-type nonlinearities. The physical implications of these equations are discussed.

The nonlinear evolution of modes on unstable stratified shear layers

NASA Technical Reports Server (NTRS)

Blackaby, Nicholas; Dando, Andrew; Hall, Philip

1993-01-01

The nonlinear development of disturbances in stratified shear flows (having a local Richardson number of value less than one quarter) is considered. Such modes are initially fast growing but, like related studies, we assume that the viscous, non-parallel spreading of the shear layer results in them evolving in a linear fashion until they reach a position where their amplitudes are large enough and their growth rates have diminished sufficiently so that amplitude equations can be derived using weakly nonlinear and non-equilibrium critical-layer theories. Four different basic integro-differential amplitude equations are possible, including one due to a novel mechanism; the relevant choice of amplitude equation, at a particular instance, being dependent on the relative sizes of the disturbance amplitude, the growth rate of the disturbance, its wavenumber, and the viscosity of the fluid. This richness of choice of possible nonlinearities arises mathematically from the indicial Frobenius roots of the governing linear inviscid equation (the Taylor-Goldstein equation) not, in general, differing by an integer. The initial nonlinear evolution of a mode will be governed by an integro-differential amplitude equations with a cubic nonlinearity but the resulting significant increase in the size of the disturbance's amplitude leads on to the next stage of the evolution process where the evolution of the mode is governed by an integro-differential amplitude equations with a quintic nonlinearity. Continued growth of the disturbance amplitude is expected during this stage, resulting in the effects of nonlinearity spreading to outside the critical level, by which time the flow has become fully nonlinear.

Finite-time H∞ filtering for non-linear stochastic systems

NASA Astrophysics Data System (ADS)

Hou, Mingzhe; Deng, Zongquan; Duan, Guangren

2016-09-01

This paper describes the robust H∞ filtering analysis and the synthesis of general non-linear stochastic systems with finite settling time. We assume that the system dynamic is modelled by Itô-type stochastic differential equations of which the state and the measurement are corrupted by state-dependent noises and exogenous disturbances. A sufficient condition for non-linear stochastic systems to have the finite-time H∞ performance with gain less than or equal to a prescribed positive number is established in terms of a certain Hamilton-Jacobi inequality. Based on this result, the existence of a finite-time H∞ filter is given for the general non-linear stochastic system by a second-order non-linear partial differential inequality, and the filter can be obtained by solving this inequality. The effectiveness of the obtained result is illustrated by a numerical example.

On some new properties of fractional derivatives with Mittag-Leffler kernel

NASA Astrophysics Data System (ADS)

Baleanu, Dumitru; Fernandez, Arran

2018-06-01

We establish a new formula for the fractional derivative with Mittag-Leffler kernel, in the form of a series of Riemann-Liouville fractional integrals, which brings out more clearly the non-locality of fractional derivatives and is easier to handle for certain computational purposes. We also prove existence and uniqueness results for certain families of linear and nonlinear fractional ODEs defined using this fractional derivative. We consider the possibility of a semigroup property for these derivatives, and establish extensions of the product rule and chain rule, with an application to fractional mechanics.

Differential morphology and image processing.

PubMed

Maragos, P

1996-01-01

Image processing via mathematical morphology has traditionally used geometry to intuitively understand morphological signal operators and set or lattice algebra to analyze them in the space domain. We provide a unified view and analytic tools for morphological image processing that is based on ideas from differential calculus and dynamical systems. This includes ideas on using partial differential or difference equations (PDEs) to model distance propagation or nonlinear multiscale processes in images. We briefly review some nonlinear difference equations that implement discrete distance transforms and relate them to numerical solutions of the eikonal equation of optics. We also review some nonlinear PDEs that model the evolution of multiscale morphological operators and use morphological derivatives. Among the new ideas presented, we develop some general 2-D max/min-sum difference equations that model the space dynamics of 2-D morphological systems (including the distance computations) and some nonlinear signal transforms, called slope transforms, that can analyze these systems in a transform domain in ways conceptually similar to the application of Fourier transforms to linear systems. Thus, distance transforms are shown to be bandpass slope filters. We view the analysis of the multiscale morphological PDEs and of the eikonal PDE solved via weighted distance transforms as a unified area in nonlinear image processing, which we call differential morphology, and briefly discuss its potential applications to image processing and computer vision.

On the reliability of computed chaotic solutions of non-linear differential equations

NASA Astrophysics Data System (ADS)

Liao, Shijun

2009-08-01

A new concept, namely the critical predictable time Tc, is introduced to give a more precise description of computed chaotic solutions of non-linear differential equations: it is suggested that computed chaotic solutions are unreliable and doubtable when t > Tc. This provides us a strategy to detect reliable solution from a given computed result. In this way, the computational phenomena, such as computational chaos (CC), computational periodicity (CP) and computational prediction uncertainty, which are mainly based on long-term properties of computed time-series, can be completely avoided. Using this concept, the famous conclusion `accurate long-term prediction of chaos is impossible' should be replaced by a more precise conclusion that `accurate prediction of chaos beyond the critical predictable time Tc is impossible'. So, this concept also provides us a timescale to determine whether or not a particular time is long enough for a given non-linear dynamic system. Besides, the influence of data inaccuracy and various numerical schemes on the critical predictable time is investigated in details by using symbolic computation software as a tool. A reliable chaotic solution of Lorenz equation in a rather large interval 0 <= t < 1200 non-dimensional Lorenz time units is obtained for the first time. It is found that the precision of the initial condition and the computed data at each time step, which is mathematically necessary to get such a reliable chaotic solution in such a long time, is so high that it is physically impossible due to the Heisenberg uncertainty principle in quantum physics. This, however, provides us a so-called `precision paradox of chaos', which suggests that the prediction uncertainty of chaos is physically unavoidable, and that even the macroscopical phenomena might be essentially stochastic and thus could be described by probability more economically.

A fast direct method for block triangular Toeplitz-like with tri-diagonal block systems from time-fractional partial differential equations

NASA Astrophysics Data System (ADS)

Ke, Rihuan; Ng, Michael K.; Sun, Hai-Wei

2015-12-01

In this paper, we study the block lower triangular Toeplitz-like with tri-diagonal blocks system which arises from the time-fractional partial differential equation. Existing fast numerical solver (e.g., fast approximate inversion method) cannot handle such linear system as the main diagonal blocks are different. The main contribution of this paper is to propose a fast direct method for solving this linear system, and to illustrate that the proposed method is much faster than the classical block forward substitution method for solving this linear system. Our idea is based on the divide-and-conquer strategy and together with the fast Fourier transforms for calculating Toeplitz matrix-vector multiplication. The complexity needs O (MNlog2 ⁡ M) arithmetic operations, where M is the number of blocks (the number of time steps) in the system and N is the size (number of spatial grid points) of each block. Numerical examples from the finite difference discretization of time-fractional partial differential equations are also given to demonstrate the efficiency of the proposed method.

Nonlinear behaviour of cantilevered carbon nanotube resonators based on a new nonlinear electrostatic load model

NASA Astrophysics Data System (ADS)

Farokhi, Hamed; Païdoussis, Michael P.; Misra, Arun K.

2018-04-01

The present study examines the nonlinear behaviour of a cantilevered carbon nanotube (CNT) resonator and its mass detection sensitivity, employing a new nonlinear electrostatic load model. More specifically, a 3D finite element model is developed in order to obtain the electrostatic load distribution on cantilevered CNT resonators. A new nonlinear electrostatic load model is then proposed accounting for the end effects due to finite length. Additionally, a new nonlinear size-dependent continuum model is developed for the cantilevered CNT resonator, employing the modified couple stress theory (to account for size-effects) together with the Kelvin-Voigt model (to account for nonlinear damping); the size-dependent model takes into account all sources of nonlinearity, i.e. geometrical and inertial nonlinearities as well as nonlinearities associated with damping, small-scale, and electrostatic load. The nonlinear equation of motion of the cantilevered CNT resonator is obtained based on the new models developed for the CNT resonator and the electrostatic load. The Galerkin method is then applied to the nonlinear equation of motion, resulting in a set of nonlinear ordinary differential equations, consisting of geometrical, inertial, electrical, damping, and size-dependent nonlinear terms. This high-dimensional nonlinear discretized model is solved numerically utilizing the pseudo-arclength continuation technique. The nonlinear static and dynamic responses of the system are examined for various cases, investigating the effect of DC and AC voltages, length-scale parameter, nonlinear damping, and electrostatic load. Moreover, the mass detection sensitivity of the system is examined for possible application of the CNT resonator as a nanosensor.

Mango (Mangifera indica L.) peel extract fractions from different cultivars differentially affect lipid accumulation in 3T3-L1 adipocyte cells.

PubMed

Taing, Meng-Wong; Pierson, Jean-Thomas; Shaw, Paul N; Dietzgen, Ralf G; Roberts-Thomson, Sarah J; Gidley, Michael J; Monteith, Gregory R

2013-02-26

Plant phytochemicals are increasingly recognised as sources of bioactive molecules which may have potential benefit in many health conditions. In mangoes, peel extracts from different cultivars exhibit varying effects on adipogenesis in the 3T3-L1 adipocyte cell line. In this study, the effects of preparative HPLC fractions of methanol peel extracts from Irwin, Nam Doc Mai and Kensington Pride mangoes were evaluated. Fraction 1 contained the most hydrophilic components while subsequent fractions contained increasingly more hydrophobic components. High content imaging was used to assess mango peel fraction effects on lipid accumulation, nuclei count and nuclear area in differentiating 3T3-L1 cells. For all three mango cultivars, the more hydrophilic peel fractions 1-3 inhibited lipid accumulation with greater potency than the more hydrophobic peel fractions 4. For all three cultivars, the more lipophilic fraction 4 had concentrations that enhanced lipid accumulation greater than fractions 1-3 as assessed by lipid droplet integrated intensity. The potency of this fraction 4 varied significantly between cultivars. Using mass spectrometry, five long chain free fatty acids were detected in fraction 4; these were not present in any other peel extract fractions. Total levels varied between cultivars, with Irwin fraction 4 containing the highest levels of these free fatty acids. Lipophilic components appear to be responsible for the lipid accumulation promoting effects of some mango extracts and are the likely cause of the diverse effects of peel extracts from different mango cultivars on lipid accumulation.

A new perspective for quintic B-spline based Crank-Nicolson-differential quadrature method algorithm for numerical solutions of the nonlinear Schrödinger equation

NASA Astrophysics Data System (ADS)

Başhan, Ali; Uçar, Yusuf; Murat Yağmurlu, N.; Esen, Alaattin

2018-01-01

In the present paper, a Crank-Nicolson-differential quadrature method (CN-DQM) based on utilizing quintic B-splines as a tool has been carried out to obtain the numerical solutions for the nonlinear Schrödinger (NLS) equation. For this purpose, first of all, the Schrödinger equation has been converted into coupled real value differential equations and then they have been discretized using both the forward difference formula and the Crank-Nicolson method. After that, Rubin and Graves linearization techniques have been utilized and the differential quadrature method has been applied to obtain an algebraic equation system. Next, in order to be able to test the efficiency of the newly applied method, the error norms, L2 and L_{∞}, as well as the two lowest invariants, I1 and I2, have been computed. Besides those, the relative changes in those invariants have been presented. Finally, the newly obtained numerical results have been compared with some of those available in the literature for similar parameters. This comparison clearly indicates that the currently utilized method, namely CN-DQM, is an effective and efficient numerical scheme and allows us to propose to solve a wide range of nonlinear equations.

Single shot, double differential spectral measurements of inverse Compton scattering in the nonlinear regime

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sakai, Y.; Gadjev, I.; Hoang, P.

Inverse Compton scattering (ICS) is a unique mechanism for producing fast pulses$-$picosecond and below$-$of bright photons, ranging from x to γ rays. These nominally narrow spectral bandwidth electromagnetic radiation pulses are efficiently produced in the interaction between intense, well-focused electron and laser beams. The spectral characteristics of such sources are affected by many experimental parameters, with intense laser effects often dominant. A laser field capable of inducing relativistic oscillatory motion may give rise to harmonic generation and, importantly for the present work, nonlinear redshifting, both of which dilute the spectral brightness of the radiation. As the applications enabled by thismore » source often depend sensitively on its spectra, it is critical to resolve the details of the wavelength and angular distribution obtained from ICS collisions. With this motivation, we present an experimental study that greatly improves on previous spectral measurement methods based on x-ray K -edge filters, by implementing a multilayer bent-crystal x-ray spectrometer. In tandem with a collimating slit, this method reveals a projection of the double differential angular-wavelength spectrum of the ICS radiation in a single shot. The measurements enabled by this diagnostic illustrate the combined off-axis and nonlinear-field-induced redshifting in the ICS emission process. The spectra obtained illustrate in detail the strength of the normalized laser vector potential, and provide a nondestructive measure of the temporal and spatial electron-laser beam overlap.« less

Single shot, double differential spectral measurements of inverse Compton scattering in the nonlinear regime

DOE PAGES

Sakai, Y.; Gadjev, I.; Hoang, P.; ...

2017-06-05

Inverse Compton scattering (ICS) is a unique mechanism for producing fast pulses$-$picosecond and below$-$of bright photons, ranging from x to γ rays. These nominally narrow spectral bandwidth electromagnetic radiation pulses are efficiently produced in the interaction between intense, well-focused electron and laser beams. The spectral characteristics of such sources are affected by many experimental parameters, with intense laser effects often dominant. A laser field capable of inducing relativistic oscillatory motion may give rise to harmonic generation and, importantly for the present work, nonlinear redshifting, both of which dilute the spectral brightness of the radiation. As the applications enabled by thismore » source often depend sensitively on its spectra, it is critical to resolve the details of the wavelength and angular distribution obtained from ICS collisions. With this motivation, we present an experimental study that greatly improves on previous spectral measurement methods based on x-ray K -edge filters, by implementing a multilayer bent-crystal x-ray spectrometer. In tandem with a collimating slit, this method reveals a projection of the double differential angular-wavelength spectrum of the ICS radiation in a single shot. The measurements enabled by this diagnostic illustrate the combined off-axis and nonlinear-field-induced redshifting in the ICS emission process. The spectra obtained illustrate in detail the strength of the normalized laser vector potential, and provide a nondestructive measure of the temporal and spatial electron-laser beam overlap.« less

Fractional Bateman—Feshbach Tikochinsky Oscillator

NASA Astrophysics Data System (ADS)

Dumitru, Baleanu; Jihad, H. Asad; Ivo, Petras

2014-02-01

In the last few years the numerical methods for solving the fractional differential equations started to be applied intensively to real world phenomena. Having these things in mind in this manuscript we focus on the fractional Lagrangian and Hamiltonian of the complex Bateman—Feshbach Tikochinsky oscillator. The numerical analysis of the corresponding fractional Euler-Lagrange equations is given within the Grünwald—Letnikov approach, which is power series expansion of the generating function.

Boundary value problems for multi-term fractional differential equations

NASA Astrophysics Data System (ADS)

Daftardar-Gejji, Varsha; Bhalekar, Sachin

2008-09-01

Multi-term fractional diffusion-wave equation along with the homogeneous/non-homogeneous boundary conditions has been solved using the method of separation of variables. It is observed that, unlike in the one term case, solution of multi-term fractional diffusion-wave equation is not necessarily non-negative, and hence does not represent anomalous diffusion of any kind.

Automated reverse engineering of nonlinear dynamical systems.

PubMed

Bongard, Josh; Lipson, Hod

2007-06-12

Complex nonlinear dynamics arise in many fields of science and engineering, but uncovering the underlying differential equations directly from observations poses a challenging task. The ability to symbolically model complex networked systems is key to understanding them, an open problem in many disciplines. Here we introduce for the first time a method that can automatically generate symbolic equations for a nonlinear coupled dynamical system directly from time series data. This method is applicable to any system that can be described using sets of ordinary nonlinear differential equations, and assumes that the (possibly noisy) time series of all variables are observable. Previous automated symbolic modeling approaches of coupled physical systems produced linear models or required a nonlinear model to be provided manually. The advance presented here is made possible by allowing the method to model each (possibly coupled) variable separately, intelligently perturbing and destabilizing the system to extract its less observable characteristics, and automatically simplifying the equations during modeling. We demonstrate this method on four simulated and two real systems spanning mechanics, ecology, and systems biology. Unlike numerical models, symbolic models have explanatory value, suggesting that automated "reverse engineering" approaches for model-free symbolic nonlinear system identification may play an increasing role in our ability to understand progressively more complex systems in the future.

Newton's method: A link between continuous and discrete solutions of nonlinear problems

NASA Technical Reports Server (NTRS)

Thurston, G. A.

1980-01-01

Newton's method for nonlinear mechanics problems replaces the governing nonlinear equations by an iterative sequence of linear equations. When the linear equations are linear differential equations, the equations are usually solved by numerical methods. The iterative sequence in Newton's method can exhibit poor convergence properties when the nonlinear problem has multiple solutions for a fixed set of parameters, unless the iterative sequences are aimed at solving for each solution separately. The theory of the linear differential operators is often a better guide for solution strategies in applying Newton's method than the theory of linear algebra associated with the numerical analogs of the differential operators. In fact, the theory for the differential operators can suggest the choice of numerical linear operators. In this paper the method of variation of parameters from the theory of linear ordinary differential equations is examined in detail in the context of Newton's method to demonstrate how it might be used as a guide for numerical solutions.

A predictive universal fractional-order differential model of wall-turbulence

NASA Astrophysics Data System (ADS)

Song, Fangying; Karniadakis, George

2017-11-01

Fractional calculus has been around for centuries but its use in computational since and engineering has emerged only recently. Here we develop a relatively simple one-dimensional model for fully-developed wall-turbulence that involves a fractional operator with variable fractional order. We use available DNS data bases to ``learn'' the function that describes the fractional order, which has a high value at the wall and decays monotonically to an asymptotic value at the centerline. We show that this function is universal upon re-scaling and hence it can be used to predict the mean velocity profile at all Reynolds numbers. We demonstrate the accuracy of our universal fractional model for channel flow at high Reynolds number as well as for pipe flow and we obtain good agreement with the Princeton super-pipe data up to Reynolds numbers 35,000,000. This work was supported by an ARO MURI Number: W911NF-15-1-0562.

The study of nonlinear almost periodic differential equations without recourse to the H-classes of these equations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Slyusarchuk, V. E., E-mail: V.E.Slyusarchuk@gmail.com, E-mail: V.Ye.Slyusarchuk@NUWM.rv.ua

2014-06-01

The well-known theorems of Favard and Amerio on the existence of almost periodic solutions to linear and nonlinear almost periodic differential equations depend to a large extent on the H-classes and the requirement that the bounded solutions of these equations be separated. The present paper provides different conditions for the existence of almost periodic solutions. These conditions, which do not depend on the H-classes of the equations, are formulated in terms of a special functional on the set of bounded solutions of the equations under consideration. This functional is used, in particular, to test whether solutions are separated. Bibliography: 24more » titles. (paper)« less

General purpose nonlinear system solver based on Newton-Krylov method.

DOE Office of Scientific and Technical Information (OSTI.GOV)

2013-12-01

KINSOL is part of a software family called SUNDIALS: SUite of Nonlinear and Differential/Algebraic equation Solvers [1]. KINSOL is a general-purpose nonlinear system solver based on Newton-Krylov and fixed-point solver technologies [2].

Fractional Order and Dynamic Simulation of a System Involving an Elastic Wide Plate

NASA Astrophysics Data System (ADS)

David, S. A.; Balthazar, J. M.; Julio, B. H. S.; Oliveira, C.

2011-09-01

Numerous researchers have studied about nonlinear dynamics in several areas of science and engineering. However, in most cases, these concepts have been explored mainly from the standpoint of analytical and computational methods involving integer order calculus (IOC). In this paper we have examined the dynamic behavior of an elastic wide plate induced by two electromagnets of a point of view of the fractional order calculus (FOC). The primary focus of this study is on to help gain a better understanding of nonlinear dynamic in fractional order systems.

Optical solitons, nonlinear self-adjointness and conservation laws for the cubic nonlinear Shrödinger's equation with repulsive delta potential

NASA Astrophysics Data System (ADS)

Baleanu, Dumitru; Inc, Mustafa; Aliyu, Aliyu Isa; Yusuf, Abdullahi

2017-11-01

In this paper, the complex envelope function ansatz method is used to acquire the optical solitons to the cubic nonlinear Shrödinger's equation with repulsive delta potential (δ-NLSE). The method reveals dark and bright optical solitons. The necessary constraint conditions which guarantee the existence of the solitons are also presented. We studied the δ-NLSE by analyzing a system of partial differential equations (PDEs) obtained by decomposing the equation into real and imaginary components. We derive the Lie point symmetry generators of the system and prove that the system is nonlinearly self-adjoint with an explicit form of a differential substitution satisfying the nonlinear self-adjoint condition. Then we use these facts to establish a set of conserved vectors for the system using the general Cls theorem presented by Ibragimov. Some interesting figures for the acquired solutions are also presented.

Modified harmonic balance method for the solution of nonlinear jerk equations

NASA Astrophysics Data System (ADS)

Rahman, M. Saifur; Hasan, A. S. M. Z.

2018-03-01

In this paper, a second approximate solution of nonlinear jerk equations (third order differential equation) can be obtained by using modified harmonic balance method. The method is simpler and easier to carry out the solution of nonlinear differential equations due to less number of nonlinear equations are required to solve than the classical harmonic balance method. The results obtained from this method are compared with those obtained from the other existing analytical methods that are available in the literature and the numerical method. The solution shows a good agreement with the numerical solution as well as the analytical methods of the available literature.

Multiplicity of solutions for the noncooperative Schrödinger-Kirchhoff system involving the fractional p-Laplacian in R^N

NASA Astrophysics Data System (ADS)

Liang, Sihua; Zhang, Jihui

2017-06-01

In this paper, we investigate the existence of solutions for the noncooperative Schrödinger-Kirchhoff-type system involving the fractional p-Laplacian and critical nonlinearities in RN. By applying the Limit Index Theory due to Li (Nonlinear Anal 25:1371-1389, 1995) and the fractional version of concentration-compactness principle, we obtain the existence and multiplicity of solutions for the above systems under some suitable assumptions. To our best knowledge, it seems that this is the first time to exploit the existence of solutions for the noncooperative Schrödinger-Kirchhoff-type system involving the fractional p-Laplacian and critical nonlinearity in RN.

Bifurcation and chaos of a new discrete fractional-order logistic map

NASA Astrophysics Data System (ADS)

Ji, YuanDong; Lai, Li; Zhong, SuChuan; Zhang, Lu

2018-04-01

The fractional-order discrete maps with chaotic behaviors based on the theory of ;fractional difference; are proposed in recent years. In this paper, instead of using fractional difference, a new fractionalized logistic map is proposed based on the numerical algorithm of fractional differentiation definition. The bifurcation diagrams of this map with various differential orders are given by numerical simulation. The simulation results show that the fractional-order logistic map derived in this manner holds rich dynamical behaviors because of its memory effect. In addition, new types of behaviors of bifurcation and chaos are found, which are different from those of the integer-order and the previous fractional-order logistic maps.

Frequency-noise cancellation in semiconductor lasers by nonlinear heterodyne detection.

PubMed

Bondurant, R S; Welford, D; Alexander, S B; Chan, V W

1986-12-01

The bit-error-rate (BER) performance of conventional noncoherent, heterodyne frequency-shift-keyed (FSK) optical communications systems can be surpassed by the use of a differential FSK modulation format and nonlinear postdetection processing at the receiver. A BER floor exists for conventional frequency-shift keying because of the frequency noise of the transmitter and local oscillator. The use of differential frequency-shift keying with nonlinear postdetection processing suppresses this BER floor for the semiconductor laser system considered here.

Bounding solutions of geometrically nonlinear viscoelastic problems

NASA Technical Reports Server (NTRS)

Stubstad, J. M.; Simitses, G. J.

1985-01-01

Integral transform techniques, such as the Laplace transform, provide simple and direct methods for solving viscoelastic problems formulated within a context of linear material response and using linear measures for deformation. Application of the transform operator reduces the governing linear integro-differential equations to a set of algebraic relations between the transforms of the unknown functions, the viscoelastic operators, and the initial and boundary conditions. Inversion either directly or through the use of the appropriate convolution theorem, provides the time domain response once the unknown functions have been expressed in terms of sums, products or ratios of known transforms. When exact inversion is not possible approximate techniques may provide accurate results. The overall problem becomes substantially more complex when nonlinear effects must be included. Situations where a linear material constitutive law can still be productively employed but where the magnitude of the resulting time dependent deformations warrants the use of a nonlinear kinematic analysis are considered. The governing equations will be nonlinear integro-differential equations for this class of problems. Thus traditional as well as approximate techniques, such as cited above, cannot be employed since the transform of a nonlinear function is not explicitly expressible.

Bounding solutions of geometrically nonlinear viscoelastic problems

NASA Technical Reports Server (NTRS)

Stubstad, J. M.; Simitses, G. J.

1986-01-01

Integral transform techniques, such as the Laplace transform, provide simple and direct methods for solving viscoelastic problems formulated within a context of linear material response and using linear measures for deformation. Application of the transform operator reduces the governing linear integro-differential equations to a set of algebraic relations between the transforms of the unknown functions, the viscoelastic operators, and the initial and boundary conditions. Inversion either directly or through the use of the appropriate convolution theorem, provides the time domain response once the unknown functions have been expressed in terms of sums, products or ratios of known transforms. When exact inversion is not possible approximate techniques may provide accurate results. The overall problem becomes substantially more complex when nonlinear effects must be included. Situations where a linear material constitutive law can still be productively employed but where the magnitude of the resulting time dependent deformations warrants the use of a nonlinear kinematic analysis are considered. The governing equations will be nonlinear integro-differential equations for this class of problems. Thus traditional as well as approximate techniques, such as cited above, cannot be employed since the transform of a nonlinear function is not explicitly expressible.

Generalized Functions for the Fractional Calculus

NASA Technical Reports Server (NTRS)

Lorenzo, Carl F.; Hartley, Tom T.

1999-01-01

Previous papers have used two important functions for the solution of fractional order differential equations, the Mittag-Leffler functionE(sub q)[at(exp q)](1903a, 1903b, 1905), and the F-function F(sub q)[a,t] of Hartley & Lorenzo (1998). These functions provided direct solution and important understanding for the fundamental linear fractional order differential equation and for the related initial value problem (Hartley and Lorenzo, 1999). This paper examines related functions and their Laplace transforms. Presented for consideration are two generalized functions, the R-function and the G-function, useful in analysis and as a basis for computation in the fractional calculus. The R-function is unique in that it contains all of the derivatives and integrals of the F-function. The R-function also returns itself on qth order differ-integration. An example application of the R-function is provided. A further generalization of the R-function, called the G-function brings in the effects of repeated and partially repeated fractional poles.

An efficient nonlinear finite-difference approach in the computational modeling of the dynamics of a nonlinear diffusion-reaction equation in microbial ecology.

PubMed

Macías-Díaz, J E; Macías, Siegfried; Medina-Ramírez, I E

2013-12-01

In this manuscript, we present a computational model to approximate the solutions of a partial differential equation which describes the growth dynamics of microbial films. The numerical technique reported in this work is an explicit, nonlinear finite-difference methodology which is computationally implemented using Newton's method. Our scheme is compared numerically against an implicit, linear finite-difference discretization of the same partial differential equation, whose computer coding requires an implementation of the stabilized bi-conjugate gradient method. Our numerical results evince that the nonlinear approach results in a more efficient approximation to the solutions of the biofilm model considered, and demands less computer memory. Moreover, the positivity of initial profiles is preserved in the practice by the nonlinear scheme proposed. Copyright © 2013 Elsevier Ltd. All rights reserved.

Automated reverse engineering of nonlinear dynamical systems

PubMed Central

Bongard, Josh; Lipson, Hod

2007-01-01

Complex nonlinear dynamics arise in many fields of science and engineering, but uncovering the underlying differential equations directly from observations poses a challenging task. The ability to symbolically model complex networked systems is key to understanding them, an open problem in many disciplines. Here we introduce for the first time a method that can automatically generate symbolic equations for a nonlinear coupled dynamical system directly from time series data. This method is applicable to any system that can be described using sets of ordinary nonlinear differential equations, and assumes that the (possibly noisy) time series of all variables are observable. Previous automated symbolic modeling approaches of coupled physical systems produced linear models or required a nonlinear model to be provided manually. The advance presented here is made possible by allowing the method to model each (possibly coupled) variable separately, intelligently perturbing and destabilizing the system to extract its less observable characteristics, and automatically simplifying the equations during modeling. We demonstrate this method on four simulated and two real systems spanning mechanics, ecology, and systems biology. Unlike numerical models, symbolic models have explanatory value, suggesting that automated “reverse engineering” approaches for model-free symbolic nonlinear system identification may play an increasing role in our ability to understand progressively more complex systems in the future. PMID:17553966

A nonlinear Kalman filtering approach to embedded control of turbocharged diesel engines

NASA Astrophysics Data System (ADS)

Rigatos, Gerasimos; Siano, Pierluigi; Arsie, Ivan

2014-10-01

The development of efficient embedded control for turbocharged Diesel engines, requires the programming of elaborated nonlinear control and filtering methods. To this end, in this paper nonlinear control for turbocharged Diesel engines is developed with the use of Differential flatness theory and the Derivative-free nonlinear Kalman Filter. It is shown that the dynamic model of the turbocharged Diesel engine is differentially flat and admits dynamic feedback linearization. It is also shown that the dynamic model can be written in the linear Brunovsky canonical form for which a state feedback controller can be easily designed. To compensate for modeling errors and external disturbances the Derivative-free nonlinear Kalman Filter is used and redesigned as a disturbance observer. The filter consists of the Kalman Filter recursion on the linearized equivalent of the Diesel engine model and of an inverse transformation based on differential flatness theory which enables to obtain estimates for the state variables of the initial nonlinear model. Once the disturbances variables are identified it is possible to compensate them by including an additional control term in the feedback loop. The efficiency of the proposed control method is tested through simulation experiments.

Dark Soliton Solutions of Space-Time Fractional Sharma-Tasso-Olver and Potential Kadomtsev-Petviashvili Equations

NASA Astrophysics Data System (ADS)

Guner, Ozkan; Korkmaz, Alper; Bekir, Ahmet

2017-02-01

Dark soliton solutions for space-time fractional Sharma-Tasso-Olver and space-time fractional potential Kadomtsev-Petviashvili equations are determined by using the properties of modified Riemann-Liouville derivative and fractional complex transform. After reducing both equations to nonlinear ODEs with constant coefficients, the \\tanh ansatz is substituted into the resultant nonlinear ODEs. The coefficients of the solutions in the ansatz are calculated by algebraic computer computations. Two different solutions are obtained for the Sharma-Tasso-Olver equation as only one solution for the potential Kadomtsev-Petviashvili equation. The solution profiles are demonstrated in 3D plots in finite domains of time and space.

Trigonometric Integrals via Partial Fractions

ERIC Educational Resources Information Center

Chen, H.; Fulford, M.

2005-01-01

Parametric differentiation is used to derive the partial fractions decompositions of certain rational functions. Those decompositions enable us to integrate some new combinations of trigonometric functions.

Linear approximations of nonlinear systems

NASA Technical Reports Server (NTRS)

Hunt, L. R.; Su, R.

1983-01-01

The development of a method for designing an automatic flight controller for short and vertical take off aircraft is discussed. This technique involves transformations of nonlinear systems to controllable linear systems and takes into account the nonlinearities of the aircraft. In general, the transformations cannot always be given in closed form. Using partial differential equations, an approximate linear system called the modified tangent model was introduced. A linear transformation of this tangent model to Brunovsky canonical form can be constructed, and from this the linear part (about a state space point x sub 0) of an exact transformation for the nonlinear system can be found. It is shown that a canonical expansion in Lie brackets about the point x sub 0 yields the same modified tangent model.

Flatness-Based Tracking Control and Nonlinear Observer for a Micro Aerial Quadcopter

NASA Astrophysics Data System (ADS)

Rivera, G.; Sawodny, O.

2010-09-01

This paper deals with the design of a nonlinear observer and a differential flat based path tracking controller for a mini aerial quadcopter. Taking into account that only the inertial coordinates and the yaw angle are available for measurements, it is shown, that the system is differentially flat, allowing a systematic design of a nonlinear tracking control in open and closed loop. A nonlinear observer is carried out to estimate the roll and pitch angle as well as all the linear and angular velocities. Finally the performance of the feedback controller and observer are illustrated in a computer simulation.

R-Function Relationships for Application in the Fractional Calculus

NASA Technical Reports Server (NTRS)

Lorenzo, Carl F.; Hartley, Tom T.

2000-01-01

The F-function, and its generalization the R-function, are of fundamental importance in the fractional calculus. It has been shown that the solution of the fundamental linear fractional differential equation may be expressed in terms of these functions. These functions serve as generalizations of the exponential function in the solution of fractional differential equations. Because of this central role in the fractional calculus, this paper explores various intrarelationships of the R-function, which will be useful in further analysis. Relationships of the R-function to the common exponential function, e(t), and its fractional derivatives are shown. From the relationships developed, some important approximations are observed. Further, the inverse relationships of the exponential function, el, in terms of the R-function are developed. Also, some approximations for the R-function are developed.

R-function relationships for application in the fractional calculus.

PubMed

Lorenzo, Carl F; Hartley, Tom T

2008-01-01

The F-function, and its generalization the R-function, are of fundamental importance in the fractional calculus. It has been shown that the solution of the fundamental linear fractional differential equation may be expressed in terms of these functions. These functions serve as generalizations of the exponential function in the solution of fractional differential equations. Because of this central role in the fractional calculus, this paper explores various intrarelationships of the R-function, which will be useful in further analysis. Relationships of the R-function to the common exponential function, et, and its fractional derivatives are shown. From the relationships developed, some important approximations are observed. Further, the inverse relationships of the exponential function, et, in terms of the R-function are developed. Also, some approximations for the R-function are developed.

Modern aspects of nonlinear convection and magnetic field in flow of thixotropic nanofluid over a nonlinear stretching sheet with variable thickness

NASA Astrophysics Data System (ADS)

Hayat, Tasawar; Qayyum, Sajid; Alsaedi, Ahmed; Ahmad, Bashir

2018-05-01

Main objective of present analysis is to study the magnetohydrodynamic (MHD) nonlinear convective flow of thixotropic nanofluid. Flow is due to nonlinear stretching surface with variable thickness. Nonlinear thermal radiation and heat generation/absorption are utilized in the energy expression. Convective conditions and zero mass flux at sheet are considered. Intention in present analysis is to develop a model for nanomaterial comprising Brownian motion and thermophoresis phenomena. Appropriate transformations are implemented for the conversion of partial differential systems into a sets of ordinary differential equations. The transformed expressions have been scrutinized through homotopic algorithm. Behavior of various sundry variables on velocity, temperature, nanoparticle concentration, skin friction coefficient and local Nusselt number are displayed through graphs. It is concluded that qualitative behaviors of temperature and thermal layer thickness are similar for radiation and temperature ratio variables. Moreover an enhancement in heat generation/absorption show rise to thermal field.

Characterizing nonlinearity in invasive EEG recordings from temporal lobe epilepsy

NASA Astrophysics Data System (ADS)

Casdagli, M. C.; Iasemidis, L. D.; Sackellares, J. C.; Roper, S. N.; Gilmore, R. L.; Savit, R. S.

Invasive electroencephalographic (EEG) recordings from depth and subdural electrodes, performed in eight patients with temporal lobe epilepsy, are analyzed using a variety of nonlinear techniques. A surrogate data technique is used to find strong evidence for nonlinearities in epileptogenic regions of the brain. Most of these nonlinearities are characterized as “spiking” by a wavelet analysis. A small fraction of the nonlinearities are characterized as “recurrent” by a nonlinear prediction algorithm. Recurrent activity is found to occur in spatio-temporal patterns related to the location of the epileptogenic focus. Residual delay maps, used to characterize “lag-one nonlinearity”, are remarkably stationary for a given electrode, and exhibit striking variations among electrodes. The clinical and theoretical implications of these results are discussed.

Dynamical Approach Study of Spurious Steady-State Numerical Solutions of Nonlinear Differential Equations. 2; Global Asymptotic Behavior of Time Discretizations; 2. Global Asymptotic Behavior of time Discretizations

NASA Technical Reports Server (NTRS)

Yee, H. C.; Sweby, P. K.

1995-01-01

The global asymptotic nonlinear behavior of 1 1 explicit and implicit time discretizations for four 2 x 2 systems of first-order autonomous nonlinear ordinary differential equations (ODES) is analyzed. The objectives are to gain a basic understanding of the difference in the dynamics of numerics between the scalars and systems of nonlinear autonomous ODEs and to set a baseline global asymptotic solution behavior of these schemes for practical computations in computational fluid dynamics. We show how 'numerical' basins of attraction can complement the bifurcation diagrams in gaining more detailed global asymptotic behavior of time discretizations for nonlinear differential equations (DEs). We show how in the presence of spurious asymptotes the basins of the true stable steady states can be segmented by the basins of the spurious stable and unstable asymptotes. One major consequence of this phenomenon which is not commonly known is that this spurious behavior can result in a dramatic distortion and, in most cases, a dramatic shrinkage and segmentation of the basin of attraction of the true solution for finite time steps. Such distortion, shrinkage and segmentation of the numerical basins of attraction will occur regardless of the stability of the spurious asymptotes, and will occur for unconditionally stable implicit linear multistep methods. In other words, for the same (common) steady-state solution the associated basin of attraction of the DE might be very different from the discretized counterparts and the numerical basin of attraction can be very different from numerical method to numerical method. The results can be used as an explanation for possible causes of error, and slow convergence and nonconvergence of steady-state numerical solutions when using the time-dependent approach for nonlinear hyperbolic or parabolic PDES.

Application of the fractional Fourier transform to image reconstruction in MRI.

PubMed

Parot, Vicente; Sing-Long, Carlos; Lizama, Carlos; Tejos, Cristian; Uribe, Sergio; Irarrazaval, Pablo

2012-07-01

The classic paradigm for MRI requires a homogeneous B(0) field in combination with linear encoding gradients. Distortions are produced when the B(0) is not homogeneous, and several postprocessing techniques have been developed to correct them. Field homogeneity is difficult to achieve, particularly for short-bore magnets and higher B(0) fields. Nonlinear magnetic components can also arise from concomitant fields, particularly in low-field imaging, or intentionally used for nonlinear encoding. In any of these situations, the second-order component is key, because it constitutes the first step to approximate higher-order fields. We propose to use the fractional Fourier transform for analyzing and reconstructing the object's magnetization under the presence of quadratic fields. The fractional fourier transform provides a precise theoretical framework for this. We show how it can be used for reconstruction and for gaining a better understanding of the quadratic field-induced distortions, including examples of reconstruction for simulated and in vivo data. The obtained images have improved quality compared with standard Fourier reconstructions. The fractional fourier transform opens a new paradigm for understanding the MR signal generated by an object under a quadratic main field or nonlinear encoding. Copyright © 2011 Wiley Periodicals, Inc.

Nonlinear ring resonator: spatial pattern generation

NASA Astrophysics Data System (ADS)

Ivanov, Vladimir Y.; Lachinova, Svetlana L.; Irochnikov, Nikita G.

2000-03-01

We consider theoretically spatial pattern formation processes in a unidirectional ring cavity with thin layer of Kerr-type nonlinear medium. Our method is based on studying of two coupled equations. The first is a partial differential equation for temporal dynamics of phase modulation of light wave in the medium. It describes nonlinear interaction in the Kerr-type lice. The second is a free propagation equation for the intracavity field complex amplitude. It involves diffraction effects of light wave in the cavity.

Exact solutions for STO and (3+1)-dimensional KdV-ZK equations using (G‧/G2) -expansion method

NASA Astrophysics Data System (ADS)

Bibi, Sadaf; Mohyud-Din, Syed Tauseef; Ullah, Rahmat; Ahmed, Naveed; Khan, Umar

This article deals with finding some exact solutions of nonlinear fractional differential equations (NLFDEs) by applying a relatively new method known as (G‧/G2) -expansion method. Solutions of space-time fractional Sharma-Tasso-Olever (STO) equation of fractional order and (3+1)-dimensional KdV-Zakharov Kuznetsov (KdV-ZK) equation of fractional order are reckoned to demonstrate the validity of this method. The fractional derivative version of modified Riemann-Liouville, linked with Fractional complex transform is employed to transform fractional differential equations into the corresponding ordinary differential equations.

Bell-polynomial approach and Wronskian determinant solutions for three sets of differential-difference nonlinear evolution equations with symbolic computation

NASA Astrophysics Data System (ADS)

Qin, Bo; Tian, Bo; Wang, Yu-Feng; Shen, Yu-Jia; Wang, Ming

2017-10-01

Under investigation in this paper are the Belov-Chaltikian (BC), Leznov and Blaszak-Marciniak (BM) lattice equations, which are associated with the conformal field theory, UToda(m_1,m_2) system and r-matrix, respectively. With symbolic computation, the Bell-polynomial approach is developed to directly bilinearize those three sets of differential-difference nonlinear evolution equations (NLEEs). This Bell-polynomial approach does not rely on any dependent variable transformation, which constitutes the key step and main difficulty of the Hirota bilinear method, and thus has the advantage in the bilinearization of the differential-difference NLEEs. Based on the bilinear forms obtained, the N-soliton solutions are constructed in terms of the N × N Wronskian determinant. Graphic illustrations demonstrate that those solutions, more general than the existing results, permit some new properties, such as the solitonic propagation and interactions for the BC lattice equations, and the nonnegative dark solitons for the BM lattice equations.

Evolutionary algorithm based heuristic scheme for nonlinear heat transfer equations.

PubMed

Ullah, Azmat; Malik, Suheel Abdullah; Alimgeer, Khurram Saleem

2018-01-01

In this paper, a hybrid heuristic scheme based on two different basis functions i.e. Log Sigmoid and Bernstein Polynomial with unknown parameters is used for solving the nonlinear heat transfer equations efficiently. The proposed technique transforms the given nonlinear ordinary differential equation into an equivalent global error minimization problem. Trial solution for the given nonlinear differential equation is formulated using a fitness function with unknown parameters. The proposed hybrid scheme of Genetic Algorithm (GA) with Interior Point Algorithm (IPA) is opted to solve the minimization problem and to achieve the optimal values of unknown parameters. The effectiveness of the proposed scheme is validated by solving nonlinear heat transfer equations. The results obtained by the proposed scheme are compared and found in sharp agreement with both the exact solution and solution obtained by Haar Wavelet-Quasilinearization technique which witnesses the effectiveness and viability of the suggested scheme. Moreover, the statistical analysis is also conducted for investigating the stability and reliability of the presented scheme.

Heterodyne interferometer with subatomic periodic nonlinearity.

PubMed

Wu, C M; Lawall, J; Deslattes, R D

1999-07-01

A new, to our knowledge, heterodyne interferometer for differential displacement measurements is presented. It is, in principle, free of periodic nonlinearity. A pair of spatially separated light beams with different frequencies is produced by two acousto-optic modulators, avoiding the main source of periodic nonlinearity in traditional heterodyne interferometers that are based on a Zeeman split laser. In addition, laser beams of the same frequency are used in the measurement and the reference arms, giving the interferometer theoretically perfect immunity from common-mode displacement. We experimentally demonstrated a residual level of periodic nonlinearity of less than 20 pm in amplitude. The remaining periodic error is attributed to unbalanced ghost reflections that drift slowly with time.

Nonlinear Tollmien-Schlichting/vortex interaction in boundary layers

NASA Technical Reports Server (NTRS)

Hall, P.; Smith, F. T.

1988-01-01

The nonlinear reaction between two oblique 3-D Tollmein-Schlichting (TS) waves and their induced streamwise-vortex flow is considered theoretically for an imcompressible boundary layer. The same theory applies to the destabilization of an incident vortex motion by subharmonic TS waves, followed by interaction. The scales and flow structure involved are addressed for high Reynolds numbers. The nonlionear interaction is powerful, starting at quite low amplitudes with a triple-deck structure for the TS waves but a large-scale structure for the induced vortex, after which strong nonlinear amplification occurs. This includes nonparallel-flow effects. The nonlinear interaction is governed by a partial differential system for the vortex flow coupled with an ordinary-differential one for the TS pressure. The solution properties found sometimes produce a breakup within a finite distance and sometimes further downstream, depending on the input amplitudes upstream and on the wave angles, and that then leads to the second stages of interaction associated with higher amplitudes, the main second stages giving either long-scale phenomena significantly affected by nonparallelism or shorter quasi-parallel ones governed by the full nonlinear triple-deck response.

The coupled nonlinear dynamics of a lift system

DOE Office of Scientific and Technical Information (OSTI.GOV)

Crespo, Rafael Sánchez, E-mail: rafael.sanchezcrespo@northampton.ac.uk, E-mail: stefan.kaczmarczyk@northampton.ac.uk, E-mail: phil.picton@northampton.ac.uk, E-mail: huijuan.su@northampton.ac.uk; Kaczmarczyk, Stefan, E-mail: rafael.sanchezcrespo@northampton.ac.uk, E-mail: stefan.kaczmarczyk@northampton.ac.uk, E-mail: phil.picton@northampton.ac.uk, E-mail: huijuan.su@northampton.ac.uk; Picton, Phil, E-mail: rafael.sanchezcrespo@northampton.ac.uk, E-mail: stefan.kaczmarczyk@northampton.ac.uk, E-mail: phil.picton@northampton.ac.uk, E-mail: huijuan.su@northampton.ac.uk

2014-12-10

Coupled lateral and longitudinal vibrations of suspension and compensating ropes in a high-rise lift system are often induced by the building motions due to wind or seismic excitations. When the frequencies of the building become near the natural frequencies of the ropes, large resonance motions of the system may result. This leads to adverse coupled dynamic phenomena involving nonplanar motions of the ropes, impact loads between the ropes and the shaft walls, as well as vertical vibrations of the car, counterweight and compensating sheave. Such an adverse dynamic behaviour of the system endangers the safety of the installation. This papermore » presents two mathematical models describing the nonlinear responses of a suspension/ compensating rope system coupled with the elevator car / compensating sheave motions. The models accommodate the nonlinear couplings between the lateral and longitudinal modes, with and without longitudinal inertia of the ropes. The partial differential nonlinear equations of motion are derived using Hamilton Principle. Then, the Galerkin method is used to discretise the equations of motion and to develop a nonlinear ordinary differential equation model. Approximate numerical solutions are determined and the behaviour of the system is analysed.« less

Quantification and parametrization of non-linearity effects by higher-order sensitivity terms in scattered light differential optical absorption spectroscopy

NASA Astrophysics Data System (ADS)

Puķīte, Jānis; Wagner, Thomas

2016-05-01

We address the application of differential optical absorption spectroscopy (DOAS) of scattered light observations in the presence of strong absorbers (in particular ozone), for which the absorption optical depth is a non-linear function of the trace gas concentration. This is the case because Beer-Lambert law generally does not hold for scattered light measurements due to many light paths contributing to the measurement. While in many cases linear approximation can be made, for scenarios with strong absorptions non-linear effects cannot always be neglected. This is especially the case for observation geometries, for which the light contributing to the measurement is crossing the atmosphere under spatially well-separated paths differing strongly in length and location, like in limb geometry. In these cases, often full retrieval algorithms are applied to address the non-linearities, requiring iterative forward modelling of absorption spectra involving time-consuming wavelength-by-wavelength radiative transfer modelling. In this study, we propose to describe the non-linear effects by additional sensitivity parameters that can be used e.g. to build up a lookup table. Together with widely used box air mass factors (effective light paths) describing the linear response to the increase in the trace gas amount, the higher-order sensitivity parameters eliminate the need for repeating the radiative transfer modelling when modifying the absorption scenario even in the presence of a strong absorption background. While the higher-order absorption structures can be described as separate fit parameters in the spectral analysis (so-called DOAS fit), in practice their quantitative evaluation requires good measurement quality (typically better than that available from current measurements). Therefore, we introduce an iterative retrieval algorithm correcting for the higher-order absorption structures not yet considered in the DOAS fit as well as the absorption dependence on

Differential branching fraction and angular moments analysis of the decay B 0 → K +π - μ + μ - in the K 0,2 * (1430)⁰ region

DOE Office of Scientific and Technical Information (OSTI.GOV)

Aaij, R.; Adeva, B.; Adinolfi, M.

Here, measurements of the differential branching fraction and angular moments of the decay B 0 → K +π - μ + μ - in the K +π - invariant mass range 1330 < m(K +π -) < 1530 MeV/c 2 are presented. Proton-proton collision data are used, corresponding to an integrated luminosity of 3 fb -1 collected by the LHCb experiment. Differential branching fraction measurements are reported in five bins of the invariant mass squared of the dimuon system, q 2, between 0.1 and 8.0 GeV 2/c 4. For the first time, an angular analysis sensitive to the S-, P-more » and D-wave contributions of this rare decay is performed. The set of 40 normalised angular moments describing the decay is presented for the q 2 range 1.1-6.0 GeV 2/c 4.« less

Differential branching fraction and angular moments analysis of the decay B 0 → K +π - μ + μ - in the K 0,2 * (1430)⁰ region

DOE PAGES

Aaij, R.; Adeva, B.; Adinolfi, M.; ...

2016-12-01

Here, measurements of the differential branching fraction and angular moments of the decay B 0 → K +π - μ + μ - in the K +π - invariant mass range 1330 < m(K +π -) < 1530 MeV/c 2 are presented. Proton-proton collision data are used, corresponding to an integrated luminosity of 3 fb -1 collected by the LHCb experiment. Differential branching fraction measurements are reported in five bins of the invariant mass squared of the dimuon system, q 2, between 0.1 and 8.0 GeV 2/c 4. For the first time, an angular analysis sensitive to the S-, P-more » and D-wave contributions of this rare decay is performed. The set of 40 normalised angular moments describing the decay is presented for the q 2 range 1.1-6.0 GeV 2/c 4.« less

The absence of lithium isotope fractionation during basalt differentiation: New measurements by multicollector sector ICP-MS

USGS Publications Warehouse

Tomascak, P.B.; Tera, F.; Helz, R.T.; Walker, R.J.

1999-01-01

We report measurements of the isotopic composition of lithium in basalts using a multicollector magnetic sector plasma-source mass spectrometer (MC-ICP-MS). This is the first application of this analytical technique to Li isotope determination. External precision of multiple replicate and duplicate measurements for a variety of sample types averages ??1.1??? (2?? population). The method allows for the rapid (???8 min/sample) analysis of small samples (???40 ng Li) relative to commonly used thermal ionization methods. The technique has been applied to a suite of samples from Kilauea Iki lava lake, Hawaii. The samples range from olivine-rich cumulitic lava to SiO2 - and K2O-enriched differentiated liquids, and have ??7Li (per mil deviation of sample 7Li/6Li relative to the L-SVEC standard) of +3.0 to +4.8. The data indicate a lack of per mil-level Li isotope fractionation as a result of crystal-liquid fractionation at temperatures greater than 1050??C. This conclusion has been tacitly assumed but never demonstrated, and is important to the interpretation of Li isotope results from such geochemically complex environments as island arcs. Copyright ?? 1999 Elsevier Science Ltd.

LYSOSOMAL FRACTIONS FROM TRANSITIONAL EPITHELIUM

PubMed Central

Kanczak, Norbert M.; Krall, Joseph I.; Hayes, E. Russell; Elliott, Willard B.

1965-01-01

Histochemical data suggested that the so called lipoid granules of transitional epithelium in some species are equivalent to lysosomes. Scrapings of bovine and canine transitional epithelium were subjected to differential centrifugation to confirm this identification biochemically. Fractions of rat liver, the classic source of lysosomes, were also prepared by the same methods to compare with the fractions obtained from urinary epithelium. In contrast to rat liver, uroepithelial fractions with a high relative specific activity for hydrolases were sedimented before the heavy mitochondria. Microscopically, these fractions contained the highest proportion of lipoid granules. The size and sedimentation characteristics of lysosomes from transitional epithelium more closely resembled those of lysosomes derived from rat kidney than those isolated from liver. PMID:14326111

Model-order reduction of lumped parameter systems via fractional calculus

NASA Astrophysics Data System (ADS)

Hollkamp, John P.; Sen, Mihir; Semperlotti, Fabio

2018-04-01

This study investigates the use of fractional order differential models to simulate the dynamic response of non-homogeneous discrete systems and to achieve efficient and accurate model order reduction. The traditional integer order approach to the simulation of non-homogeneous systems dictates the use of numerical solutions and often imposes stringent compromises between accuracy and computational performance. Fractional calculus provides an alternative approach where complex dynamical systems can be modeled with compact fractional equations that not only can still guarantee analytical solutions, but can also enable high levels of order reduction without compromising on accuracy. Different approaches are explored in order to transform the integer order model into a reduced order fractional model able to match the dynamic response of the initial system. Analytical and numerical results show that, under certain conditions, an exact match is possible and the resulting fractional differential models have both a complex and frequency-dependent order of the differential operator. The implications of this type of approach for both model order reduction and model synthesis are discussed.

SIMULTANEOUS DIFFERENTIAL EQUATION COMPUTER

DOEpatents

Collier, D.M.; Meeks, L.A.; Palmer, J.P.

1960-05-10

A description is given for an electronic simulator for a system of simultaneous differential equations, including nonlinear equations. As a specific example, a homogeneous nuclear reactor system including a reactor fluid, heat exchanger, and a steam boiler may be simulated, with the nonlinearity resulting from a consideration of temperature effects taken into account. The simulator includes three operational amplifiers, a multiplier, appropriate potential sources, and interconnecting R-C networks.

Toward lattice fractional vector calculus

NASA Astrophysics Data System (ADS)

Tarasov, Vasily E.

2014-09-01

An analog of fractional vector calculus for physical lattice models is suggested. We use an approach based on the models of three-dimensional lattices with long-range inter-particle interactions. The lattice analogs of fractional partial derivatives are represented by kernels of lattice long-range interactions, where the Fourier series transformations of these kernels have a power-law form with respect to wave vector components. In the continuum limit, these lattice partial derivatives give derivatives of non-integer order with respect to coordinates. In the three-dimensional description of the non-local continuum, the fractional differential operators have the form of fractional partial derivatives of the Riesz type. As examples of the applications of the suggested lattice fractional vector calculus, we give lattice models with long-range interactions for the fractional Maxwell equations of non-local continuous media and for the fractional generalization of the Mindlin and Aifantis continuum models of gradient elasticity.

Fractional calculus in bioengineering, part 3.

PubMed

Magin, Richard L

2004-01-01

Fractional calculus (integral and differential operations of noninteger order) is not often used to model biological systems. Although the basic mathematical ideas were developed long ago by the mathematicians Leibniz (1695), Liouville (1834), Riemann (1892), and others and brought to the attention of the engineering world by Oliver Heaviside in the 1890s, it was not until 1974 that the first book on the topic was published by Oldham and Spanier. Recent monographs and symposia proceedings have highlighted the application of fractional calculus in physics, continuum mechanics, signal processing, and electromagnetics, but with few examples of applications in bioengineering. This is surprising because the methods of fractional calculus, when defined as a Laplace or Fourier convolution product, are suitable for solving many problems in biomedical research. For example, early studies by Cole (1933) and Hodgkin (1946) of the electrical properties of nerve cell membranes and the propagation of electrical signals are well characterized by differential equations of fractional order. The solution involves a generalization of the exponential function to the Mittag-Leffler function, which provides a better fit to the observed cell membrane data. A parallel application of fractional derivatives to viscoelastic materials establishes, in a natural way, hereditary integrals and the power law (Nutting/Scott Blair) stress-strain relationship for modeling biomaterials. In this review, I will introduce the idea of fractional operations by following the original approach of Heaviside, demonstrate the basic operations of fractional calculus on well-behaved functions (step, ramp, pulse, sinusoid) of engineering interest, and give specific examples from electrochemistry, physics, bioengineering, and biophysics. The fractional derivative accurately describes natural phenomena that occur in such common engineering problems as heat transfer, electrode/electrolyte behavior, and sub

Nonlinearly preconditioned semismooth Newton methods for variational inequality solution of two-phase flow in porous media

NASA Astrophysics Data System (ADS)

Yang, Haijian; Sun, Shuyu; Yang, Chao

2017-03-01

Most existing methods for solving two-phase flow problems in porous media do not take the physically feasible saturation fractions between 0 and 1 into account, which often destroys the numerical accuracy and physical interpretability of the simulation. To calculate the solution without the loss of this basic requirement, we introduce a variational inequality formulation of the saturation equilibrium with a box inequality constraint, and use a conservative finite element method for the spatial discretization and a backward differentiation formula with adaptive time stepping for the temporal integration. The resulting variational inequality system at each time step is solved by using a semismooth Newton algorithm. To accelerate the Newton convergence and improve the robustness, we employ a family of adaptive nonlinear elimination methods as a nonlinear preconditioner. Some numerical results are presented to demonstrate the robustness and efficiency of the proposed algorithm. A comparison is also included to show the superiority of the proposed fully implicit approach over the classical IMplicit Pressure-Explicit Saturation (IMPES) method in terms of the time step size and the total execution time measured on a parallel computer.

Equilibrium Tin Isotope Fractionation during Metal-Sulfide-Silicate Differentiation: A Nuclear Resonant Inelastic X-ray Scattering Approach

NASA Astrophysics Data System (ADS)

Roskosz, M.; Amet, Q.; Fitoussi, C.; Laporte, D.; Hu, M. Y.; Alp, E. E.

2016-12-01

Metal-silicate differentiation was recently addressed through the insight of the isotopic composition of siderophile elements (mainly Fe, Si and Cr isotopes) of planetary and extraterrestrial bodies. A key limitation of this approach is however the knowledge of equilibrium fractionation factors between coexisting phases (metal alloys, silicates and sulfides) used to interpret data on natural samples. These properties are difficult to determine experimentally. In this context, tin is generally classified as a chalcophile element but it is also siderophile and volatile. We applied a synchrotron-based method to circumvent difficulties related to determination of equilibrium isotope fractionation. The nuclear resonant inelastic x-ray scattering (NRIXS) was used to measure the phonon excitation spectrum and then to derive the force constant and finally the fractionation factors of Sn-bearing geomaterials. Spectroscopic measurements were carried out at room pressure at Sector 30-ID (APS, USA). A range of Fe-Ni alloys, rhyolitic and basaltic glasses and iron sulfides containing isotopically enriched 119Sn were synthesized. The tin content and the redox conditions prevailing during the synthesis were varied. The data evaluation was carried out using PHOENIX and SciPhon programs. A strong effect of both the redox state and the tin content was measured. In addition, the composition of the silicate glasses was found to be another important factor determining the tin isotope metal-silicate-sulfide fractionation factors. Our results are consistent with trends previously observed in the case of iron isotopes [1,2]. We will discuss the implications of our experimental results in terms of tin isotope planetary signatures. References: [1] Dauphas et al. (2014), EPSL, 398, 127-140; [2] Roskosz et al. (2015), GCA, 169, 184-199.

Fractional power-law spatial dispersion in electrodynamics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tarasov, Vasily E., E-mail: tarasov@theory.sinp.msu.ru; Departamento de Análisis Matemático, Universidad de La Laguna, 38271 La Laguna, Tenerife; Trujillo, Juan J., E-mail: jtrujill@ullmat.es

2013-07-15

Electric fields in non-local media with power-law spatial dispersion are discussed. Equations involving a fractional Laplacian in the Riesz form that describe the electric fields in such non-local media are studied. The generalizations of Coulomb’s law and Debye’s screening for power-law non-local media are characterized. We consider simple models with anomalous behavior of plasma-like media with power-law spatial dispersions. The suggested fractional differential models for these plasma-like media are discussed to describe non-local properties of power-law type. -- Highlights: •Plasma-like non-local media with power-law spatial dispersion. •Fractional differential equations for electric fields in the media. •The generalizations of Coulomb’s lawmore » and Debye’s screening for the media.« less

Nonlinear optical properties of semiconductor nanocrystals

NASA Astrophysics Data System (ADS)

Ricard, Gianpiero Banfi Vittorio Degiorgio Daniel

1998-05-01

This review is devoted to the description of recent experimental results concerning the nonlinear optical properties of semiconductor-doped glasses SDGs with particular emphasis on the regime in which the energy of the incident photon is smaller than the energy gap. A considerable theoretical and experimental effort has been devoted in the last 10years to the fundamental aspects of quantumconfined structures, which have properties somewhat intermediate between the bulk crystals and atoms or molecules. From this point of view, SDGs represent an easily available test system, and optical techniques have been a major diagnostic tool. Luminescence and absorption spectroscopy were extensively used to characterize the electronic states. The experiments aimed at the measurement of the real and imaginary parts of the third-order optical susceptibility of SDGs below the bandgap are described in some detail, and the results obtained with different techniques are compared. Besides the intrinsic fast nonlinearity due to bound electrons, SDGs may present a larger but much slower nonlinearity due to the free carriers generated by two-photon absorption. This implies that experiments have to be properly designed for separation of the two effects. In this article we stress the importance of a detailed structural characterization of the samples. Knowledge of the volume fraction occupied by the nanocrystals is necessary in order to derive from the experimental data the intrinsic nonlinearity and to compare it with the bulk nonlinearity. We discuss recent experiments in which the dependence of the intrinsic nonlinearity on the crystal size is derived by performing, on the samples, measurements of the real part and imaginary part of the nonlinear optical susceptibility and measurements of crystal size and volume fraction. Structural characterization is of interest also for a better understanding of the physical processes underlying the growth of crystallites in SDGs. The average size of

Modeling of long-range memory processes with inverse cubic distributions by the nonlinear stochastic differential equations

NASA Astrophysics Data System (ADS)

Kaulakys, B.; Alaburda, M.; Ruseckas, J.

2016-05-01

A well-known fact in the financial markets is the so-called ‘inverse cubic law’ of the cumulative distributions of the long-range memory fluctuations of market indicators such as a number of events of trades, trading volume and the logarithmic price change. We propose the nonlinear stochastic differential equation (SDE) giving both the power-law behavior of the power spectral density and the long-range dependent inverse cubic law of the cumulative distribution. This is achieved using the suggestion that when the market evolves from calm to violent behavior there is a decrease of the delay time of multiplicative feedback of the system in comparison to the driving noise correlation time. This results in a transition from the Itô to the Stratonovich sense of the SDE and yields a long-range memory process.

Linear and nonlinear stability of the Blasius boundary layer

NASA Technical Reports Server (NTRS)

Bertolotti, F. P.; Herbert, TH.; Spalart, P. R.

1992-01-01

Two new techniques for the study of the linear and nonlinear instability in growing boundary layers are presented. The first technique employs partial differential equations of parabolic type exploiting the slow change of the mean flow, disturbance velocity profiles, wavelengths, and growth rates in the streamwise direction. The second technique solves the Navier-Stokes equation for spatially evolving disturbances using buffer zones adjacent to the inflow and outflow boundaries. Results of both techniques are in excellent agreement. The linear and nonlinear development of Tollmien-Schlichting (TS) waves in the Blasius boundary layer is investigated with both techniques and with a local procedure based on a system of ordinary differential equations. The results are compared with previous work and the effects of non-parallelism and nonlinearity are clarified. The effect of nonparallelism is confirmed to be weak and, consequently, not responsible for the discrepancies between measurements and theoretical results for parallel flow.

Universal nonlinear small-scale dynamo.

PubMed

Beresnyak, A

2012-01-20

We consider astrophysically relevant nonlinear MHD dynamo at large Reynolds numbers (Re). We argue that it is universal in a sense that magnetic energy grows at a rate which is a constant fraction C(E) of the total turbulent dissipation rate. On the basis of locality bounds we claim that this "efficiency of the small-scale dynamo", C(E), is a true constant for large Re and is determined only by strongly nonlinear dynamics at the equipartition scale. We measured C(E) in numerical simulations and observed a value around 0.05 in the highest resolution simulations. We address the issue of C(E) being small, unlike the Kolmogorov constant which is of order unity. © 2012 American Physical Society

A recurrent neural network for nonlinear optimization with a continuously differentiable objective function and bound constraints.

PubMed

Liang, X B; Wang, J

2000-01-01

This paper presents a continuous-time recurrent neural-network model for nonlinear optimization with any continuously differentiable objective function and bound constraints. Quadratic optimization with bound constraints is a special problem which can be solved by the recurrent neural network. The proposed recurrent neural network has the following characteristics. 1) It is regular in the sense that any optimum of the objective function with bound constraints is also an equilibrium point of the neural network. If the objective function to be minimized is convex, then the recurrent neural network is complete in the sense that the set of optima of the function with bound constraints coincides with the set of equilibria of the neural network. 2) The recurrent neural network is primal and quasiconvergent in the sense that its trajectory cannot escape from the feasible region and will converge to the set of equilibria of the neural network for any initial point in the feasible bound region. 3) The recurrent neural network has an attractivity property in the sense that its trajectory will eventually converge to the feasible region for any initial states even at outside of the bounded feasible region. 4) For minimizing any strictly convex quadratic objective function subject to bound constraints, the recurrent neural network is globally exponentially stable for almost any positive network parameters. Simulation results are given to demonstrate the convergence and performance of the proposed recurrent neural network for nonlinear optimization with bound constraints.

Global Regularity for the Fractional Euler Alignment System

NASA Astrophysics Data System (ADS)

Do, Tam; Kiselev, Alexander; Ryzhik, Lenya; Tan, Changhui

2018-04-01

We study a pressureless Euler system with a non-linear density-dependent alignment term, originating in the Cucker-Smale swarming models. The alignment term is dissipative in the sense that it tends to equilibrate the velocities. Its density dependence is natural: the alignment rate increases in the areas of high density due to species discomfort. The diffusive term has the order of a fractional Laplacian {(-partial _{xx})^{α/2}, α \\in (0, 1)}. The corresponding Burgers equation with a linear dissipation of this type develops shocks in a finite time. We show that the alignment nonlinearity enhances the dissipation, and the solutions are globally regular for all {α \\in (0, 1)}. To the best of our knowledge, this is the first example of such regularization due to the non-local nonlinear modulation of dissipation.

Stochastic Evolution Equations Driven by Fractional Noises

DTIC Science & Technology

2016-11-28

rate of convergence to zero or the error and the limit in distribution of the error fluctuations. We have studied time discrete numerical schemes...error fluctuations. We have studied time discrete numerical schemes based on Taylor expansions for rough differential equations and for stochastic...variations of the time discrete Taylor schemes for rough differential equations and for stochastic differential equations driven by fractional Brownian

On the limits of probabilistic forecasting in nonlinear time series analysis II: Differential entropy.

PubMed

Amigó, José M; Hirata, Yoshito; Aihara, Kazuyuki

2017-08-01

In a previous paper, the authors studied the limits of probabilistic prediction in nonlinear time series analysis in a perfect model scenario, i.e., in the ideal case that the uncertainty of an otherwise deterministic model is due to only the finite precision of the observations. The model consisted of the symbolic dynamics of a measure-preserving transformation with respect to a finite partition of the state space, and the quality of the predictions was measured by the so-called ignorance score, which is a conditional entropy. In practice, though, partitions are dispensed with by considering numerical and experimental data to be continuous, which prompts us to trade off in this paper the Shannon entropy for the differential entropy. Despite technical differences, we show that the core of the previous results also hold in this extended scenario for sufficiently high precision. The corresponding imperfect model scenario will be revisited too because it is relevant for the applications. The theoretical part and its application to probabilistic forecasting are illustrated with numerical simulations and a new prediction algorithm.

Differential geometry techniques for sets of nonlinear partial differential equations

NASA Technical Reports Server (NTRS)

Estabrook, Frank B.

1990-01-01

An attempt is made to show that the Cartan theory of partial differential equations can be a useful technique for applied mathematics. Techniques for finding consistent subfamilies of solutions that are generically rich and well-posed and for introducing potentials or other usefully consistent auxiliary fields are introduced. An extended sample calculation involving the Korteweg-de Vries equation is given.

Two-dimensional linear and nonlinear Talbot effect from rogue waves.

PubMed

Zhang, Yiqi; Belić, Milivoj R; Petrović, Milan S; Zheng, Huaibin; Chen, Haixia; Li, Changbiao; Lu, Keqing; Zhang, Yanpeng

2015-03-01

We introduce two-dimensional (2D) linear and nonlinear Talbot effects. They are produced by propagating periodic 2D diffraction patterns and can be visualized as 3D stacks of Talbot carpets. The nonlinear Talbot effect originates from 2D rogue waves and forms in a bulk 3D nonlinear medium. The recurrences of an input rogue wave are observed at the Talbot length and at the half-Talbot length, with a π phase shift; no other recurrences are observed. Differing from the nonlinear Talbot effect, the linear effect displays the usual fractional Talbot images as well. We also find that the smaller the period of incident rogue waves, the shorter the Talbot length. Increasing the beam intensity increases the Talbot length, but above a threshold this leads to a catastrophic self-focusing phenomenon which destroys the effect. We also find that the Talbot recurrence can be viewed as a self-Fourier transform of the initial periodic beam that is automatically performed during propagation. In particular, linear Talbot effect can be viewed as a fractional self-Fourier transform, whereas the nonlinear Talbot effect can be viewed as the regular self-Fourier transform. Numerical simulations demonstrate that the rogue-wave initial condition is sufficient but not necessary for the observation of the effect. It may also be observed from other periodic inputs, provided they are set on a finite background. The 2D effect may find utility in the production of 3D photonic crystals.

Analysis of Nonlinear Periodic and Aperiodic Media: Application to Optical Logic Gates

NASA Astrophysics Data System (ADS)

Yu, Yisheng

This dissertation is about the analysis of nonlinear periodic and aperiodic media and their application to the design of intensity controlled all optical logic gates: AND, OR, and NOT. A coupled nonlinear differential equation that characterizes the electromagnetic wave propagation in a nonlinear periodic (and aperiodic) medium has been derived from the first principle. The equations are general enough that it reflects the effect of transverse modal fields and can be used to analyze both co-propagating and counter propagating waves. A numerical technique based on the finite differences method and absorbing boundary condition has been developed to solve the coupled differential equations here. The numerical method is simple and accurate. Unlike the method based on characteristics that has been reported in the literature, this method does not involve integration and step sizes of time and space coordinates are decoupled. The decoupling provides independent choice for time and space step sizes. The concept of "gap soliton" has also been re-examined. The dissertation consists of four manuscripts. Manuscript I reports on the design of all optical logic gates: AND, OR, and NOT based on the bistability property of nonlinear periodic and aperiodic waveguiding structures. The functioning of the logic gates has been shown by analysis. The numerical technique that has been developed to solve the nonlinear differential equations are addressed in manuscript II. The effect of transverse modal fields on the bistable property of nonlinear periodic medium is reported in manuscript III. The concept of "gap soliton" that are generated in a nonlinear periodic medium has been re-examined. The details on the finding of the re-examination are discussed in manuscript IV.

Application of the Homotopy Perturbation Method to the Nonlinear Pendulum

ERIC Educational Resources Information Center

Belendez, A.; Hernandez, A.; Belendez, T.; Marquez, A.

2007-01-01

The homotopy perturbation method is used to solve the nonlinear differential equation that governs the nonlinear oscillations of a simple pendulum, and an approximate expression for its period is obtained. Only one iteration leads to high accuracy of the solutions and the relative error for the approximate period is less than 2% for amplitudes as…

Bayesian Estimation of Panel Data Fractional Response Models with Endogeneity: An Application to Standardized Test Rates

ERIC Educational Resources Information Center

Kessler, Lawrence M.

2013-01-01

In this paper I propose Bayesian estimation of a nonlinear panel data model with a fractional dependent variable (bounded between 0 and 1). Specifically, I estimate a panel data fractional probit model which takes into account the bounded nature of the fractional response variable. I outline estimation under the assumption of strict exogeneity as…

Colloquium: Fractional calculus view of complexity: A tutorial

NASA Astrophysics Data System (ADS)

West, Bruce J.

2014-10-01

The fractional calculus has been part of the mathematics and science literature for 310 years. However, it is only in the past decade or so that it has drawn the attention of mainstream science as a way to describe the dynamics of complex phenomena with long-term memory, spatial heterogeneity, along with nonstationary and nonergodic statistics. The most recent application encompasses complex networks, which require new ways of thinking about the world. Part of the new cognition is provided by the fractional calculus description of temporal and topological complexity. Consequently, this Colloquium is not so much a tutorial on the mathematics of the fractional calculus as it is an exploration of how complex phenomena in the physical, social, and life sciences that have eluded traditional mathematical modeling become less mysterious when certain historical assumptions such as differentiability are discarded and the ordinary calculus is replaced with the fractional calculus. Exemplars considered include the fractional differential equations describing the dynamics of viscoelastic materials, turbulence, foraging, and phase transitions in complex social networks.

Systematic Computation of Nonlinear Cellular and Molecular Dynamics with Low-Power CytoMimetic Circuits: A Simulation Study

PubMed Central

Papadimitriou, Konstantinos I.; Stan, Guy-Bart V.; Drakakis, Emmanuel M.

2013-01-01

This paper presents a novel method for the systematic implementation of low-power microelectronic circuits aimed at computing nonlinear cellular and molecular dynamics. The method proposed is based on the Nonlinear Bernoulli Cell Formalism (NBCF), an advanced mathematical framework stemming from the Bernoulli Cell Formalism (BCF) originally exploited for the modular synthesis and analysis of linear, time-invariant, high dynamic range, logarithmic filters. Our approach identifies and exploits the striking similarities existing between the NBCF and coupled nonlinear ordinary differential equations (ODEs) typically appearing in models of naturally encountered biochemical systems. The resulting continuous-time, continuous-value, low-power CytoMimetic electronic circuits succeed in simulating fast and with good accuracy cellular and molecular dynamics. The application of the method is illustrated by synthesising for the first time microelectronic CytoMimetic topologies which simulate successfully: 1) a nonlinear intracellular calcium oscillations model for several Hill coefficient values and 2) a gene-protein regulatory system model. The dynamic behaviours generated by the proposed CytoMimetic circuits are compared and found to be in very good agreement with their biological counterparts. The circuits exploit the exponential law codifying the low-power subthreshold operation regime and have been simulated with realistic parameters from a commercially available CMOS process. They occupy an area of a fraction of a square-millimetre, while consuming between 1 and 12 microwatts of power. Simulations of fabrication-related variability results are also presented. PMID:23393550

Nonlinear Acoustical Assessment of Precipitate Nucleation

NASA Technical Reports Server (NTRS)

Cantrell, John H.; Yost, William T.

2004-01-01

The purpose of the present work is to show that measurements of the acoustic nonlinearity parameter in heat treatable alloys as a function of heat treatment time can provide quantitative information about the kinetics of precipitate nucleation and growth in such alloys. Generally, information on the kinetics of phase transformations is obtained from time-sequenced electron microscopical examination and differential scanning microcalorimetry. The present nonlinear acoustical assessment of precipitation kinetics is based on the development of a multiparameter analytical model of the effects on the nonlinearity parameter of precipitate nucleation and growth in the alloy system. A nonlinear curve fit of the model equation to the experimental data is then used to extract the kinetic parameters related to the nucleation and growth of the targeted precipitate. The analytical model and curve fit is applied to the assessment of S' precipitation in aluminum alloy 2024 during artificial aging from the T4 to the T6 temper.

Double symbolic joint entropy in nonlinear dynamic complexity analysis

NASA Astrophysics Data System (ADS)

Yao, Wenpo; Wang, Jun

2017-07-01

Symbolizations, the base of symbolic dynamic analysis, are classified as global static and local dynamic approaches which are combined by joint entropy in our works for nonlinear dynamic complexity analysis. Two global static methods, symbolic transformations of Wessel N. symbolic entropy and base-scale entropy, and two local ones, namely symbolizations of permutation and differential entropy, constitute four double symbolic joint entropies that have accurate complexity detections in chaotic models, logistic and Henon map series. In nonlinear dynamical analysis of different kinds of heart rate variability, heartbeats of healthy young have higher complexity than those of the healthy elderly, and congestive heart failure (CHF) patients are lowest in heartbeats' joint entropy values. Each individual symbolic entropy is improved by double symbolic joint entropy among which the combination of base-scale and differential symbolizations have best complexity analysis. Test results prove that double symbolic joint entropy is feasible in nonlinear dynamic complexity analysis.

Nonlinear random response prediction using MSC/NASTRAN

NASA Technical Reports Server (NTRS)

Robinson, J. H.; Chiang, C. K.; Rizzi, S. A.

1993-01-01

An equivalent linearization technique was incorporated into MSC/NASTRAN to predict the nonlinear random response of structures by means of Direct Matrix Abstract Programming (DMAP) modifications and inclusion of the nonlinear differential stiffness module inside the iteration loop. An iterative process was used to determine the rms displacements. Numerical results obtained for validation on simple plates and beams are in good agreement with existing solutions in both the linear and linearized regions. The versatility of the implementation will enable the analyst to determine the nonlinear random responses for complex structures under combined loads. The thermo-acoustic response of a hexagonal thermal protection system panel is used to highlight some of the features of the program.

Evans functions and bifurcations of nonlinear waves of some nonlinear reaction diffusion equations

NASA Astrophysics Data System (ADS)

Zhang, Linghai

2017-10-01

The main purposes of this paper are to accomplish the existence, stability, instability and bifurcation of the nonlinear waves of the nonlinear system of reaction diffusion equations ut =uxx + α [ βH (u - θ) - u ] - w, wt = ε (u - γw) and to establish the existence, stability, instability and bifurcation of the nonlinear waves of the nonlinear scalar reaction diffusion equation ut =uxx + α [ βH (u - θ) - u ], under different conditions on the model constants. To establish the bifurcation for the system, we will study the existence and instability of a standing pulse solution if 0 < 2 (1 + αγ) θ < αβγ; the existence and stability of two standing wave fronts if 2 (1 + αγ) θ = αβγ and γ2 ε > 1; the existence and instability of two standing wave fronts if 2 (1 + αγ) θ = αβγ and 0 <γ2 ε < 1; the existence and instability of an upside down standing pulse solution if 0 < (1 + αγ) θ < αβγ < 2 (1 + αγ) θ. To establish the bifurcation for the scalar equation, we will study the existence and stability of a traveling wave front as well as the existence and instability of a standing pulse solution if 0 < 2 θ < β; the existence and stability of two standing wave fronts if 2 θ = β; the existence and stability of a traveling wave front as well as the existence and instability of an upside down standing pulse solution if 0 < θ < β < 2 θ. By the way, we will also study the existence and stability of a traveling wave back of the nonlinear scalar reaction diffusion equation ut =uxx + α [ βH (u - θ) - u ] -w0, where w0 = α (β - 2 θ) > 0 is a positive constant, if 0 < 2 θ < β. To achieve the main goals, we will make complete use of the special structures of the model equations and we will construct Evans functions and apply them to study the eigenvalues and eigenfunctions of several eigenvalue problems associated with several linear differential operators. It turns out that a complex number λ0 is an eigenvalue of the linear

Darcy-Forchheimer flow of Maxwell nanofluid flow with nonlinear thermal radiation and activation energy

NASA Astrophysics Data System (ADS)

Sajid, T.; Sagheer, M.; Hussain, S.; Bilal, M.

2018-03-01

The present article is about the study of Darcy-Forchheimer flow of Maxwell nanofluid over a linear stretching surface. Effects like variable thermal conductivity, activation energy, nonlinear thermal radiation is also incorporated for the analysis of heat and mass transfer. The governing nonlinear partial differential equations (PDEs) with convective boundary conditions are first converted into the nonlinear ordinary differential equations (ODEs) with the help of similarity transformation, and then the resulting nonlinear ODEs are solved with the help of shooting method and MATLAB built-in bvp4c solver. The impact of different physical parameters like Brownian motion, thermophoresis parameter, Reynolds number, magnetic parameter, nonlinear radiative heat flux, Prandtl number, Lewis number, reaction rate constant, activation energy and Biot number on Nusselt number, velocity, temperature and concentration profile has been discussed. It is viewed that both thermophoresis parameter and activation energy parameter has ascending effect on the concentration profile.

Anomalous Transport of Cosmic Rays in a Nonlinear Diffusion Model

DOE Office of Scientific and Technical Information (OSTI.GOV)

Litvinenko, Yuri E.; Fichtner, Horst; Walter, Dominik

2017-05-20

We investigate analytically and numerically the transport of cosmic rays following their escape from a shock or another localized acceleration site. Observed cosmic-ray distributions in the vicinity of heliospheric and astrophysical shocks imply that anomalous, superdiffusive transport plays a role in the evolution of the energetic particles. Several authors have quantitatively described the anomalous diffusion scalings, implied by the data, by solutions of a formal transport equation with fractional derivatives. Yet the physical basis of the fractional diffusion model remains uncertain. We explore an alternative model of the cosmic-ray transport: a nonlinear diffusion equation that follows from a self-consistent treatmentmore » of the resonantly interacting cosmic-ray particles and their self-generated turbulence. The nonlinear model naturally leads to superdiffusive scalings. In the presence of convection, the model yields a power-law dependence of the particle density on the distance upstream of the shock. Although the results do not refute the use of a fractional advection–diffusion equation, they indicate a viable alternative to explain the anomalous diffusion scalings of cosmic-ray particles.« less

Incorporating nonlinearity into mediation analyses.

PubMed

Knafl, George J; Knafl, Kathleen A; Grey, Margaret; Dixon, Jane; Deatrick, Janet A; Gallo, Agatha M

2017-03-21

Mediation is an important issue considered in the behavioral, medical, and social sciences. It addresses situations where the effect of a predictor variable X on an outcome variable Y is explained to some extent by an intervening, mediator variable M. Methods for addressing mediation have been available for some time. While these methods continue to undergo refinement, the relationships underlying mediation are commonly treated as linear in the outcome Y, the predictor X, and the mediator M. These relationships, however, can be nonlinear. Methods are needed for assessing when mediation relationships can be treated as linear and for estimating them when they are nonlinear. Existing adaptive regression methods based on fractional polynomials are extended here to address nonlinearity in mediation relationships, but assuming those relationships are monotonic as would be consistent with theories about directionality of such relationships. Example monotonic mediation analyses are provided assessing linear and monotonic mediation of the effect of family functioning (X) on a child's adaptation (Y) to a chronic condition by the difficulty (M) for the family in managing the child's condition. Example moderated monotonic mediation and simulation analyses are also presented. Adaptive methods provide an effective way to incorporate possibly nonlinear monotonicity into mediation relationships.

Voltage-Induced Nonlinear Conduction Properties of Epoxy Resin/Micron-Silver Particles Composites

NASA Astrophysics Data System (ADS)

Qu, Zhaoming; Lu, Pin; Yuan, Yang; Wang, Qingguo

2018-01-01

The nonlinear conduction properties of epoxy resin (ER)/micron-silver particles (MP) composites were investigated. Under sufficient high intensity applied constant voltage, the obvious nonlinear conduction properties of the samples with volume fraction 25% were found. With increments in the voltage, the conductive switching effect was observed. The nonlinear conduction mechanism of the ER/MP composites under high applied voltages could be attributed to the electrical current conducted via discrete paths of conductive particles induced by the electric field. The test results show that the ER/MP composites with nonlinear conduction properties are of great potential application in electromagnetic protection of electron devices and systems.

Chaotic vibrations of the duffing system with fractional damping

DOE Office of Scientific and Technical Information (OSTI.GOV)

Syta, Arkadiusz; Litak, Grzegorz; Lenci, Stefano

2014-03-15

We examined the Duffing system with a fractional damping term. Calculating the basins of attraction, we demonstrate a broad spectrum of non-linear behaviour connected with sensitivity to the initial conditions and chaos. To quantify dynamical response of the system, we propose the statistical 0-1 test as well as the maximal Lyapunov exponent; the application of the latter encounter a few difficulties because of the memory effect due to the fractional derivative. The results are confirmed by bifurcation diagrams, phase portraits, and Poincaré sections.

A pseudo energy-invariant method for relativistic wave equations with Riesz space-fractional derivatives

NASA Astrophysics Data System (ADS)

Macías-Díaz, J. E.; Hendy, A. S.; De Staelen, R. H.

2018-03-01

In this work, we investigate a general nonlinear wave equation with Riesz space-fractional derivatives that generalizes various classical hyperbolic models, including the sine-Gordon and the Klein-Gordon equations from relativistic quantum mechanics. A finite-difference discretization of the model is provided using fractional centered differences. The method is a technique that is capable of preserving an energy-like quantity at each iteration. Some computational comparisons against solutions available in the literature are performed in order to assess the capability of the method to preserve the invariant. Our experiments confirm that the technique yields good approximations to the solutions considered. As an application of our scheme, we provide simulations that confirm, for the first time in the literature, the presence of the phenomenon of nonlinear supratransmission in Riesz space-fractional Klein-Gordon equations driven by a harmonic perturbation at the boundary.

A review and evaluation of numerical tools for fractional calculus and fractional order controls

NASA Astrophysics Data System (ADS)

Li, Zhuo; Liu, Lu; Dehghan, Sina; Chen, YangQuan; Xue, Dingyü

2017-06-01

In recent years, as fractional calculus becomes more and more broadly used in research across different academic disciplines, there are increasing demands for the numerical tools for the computation of fractional integration/differentiation, and the simulation of fractional order systems. Time to time, being asked about which tool is suitable for a specific application, the authors decide to carry out this survey to present recapitulative information of the available tools in the literature, in hope of benefiting researchers with different academic backgrounds. With this motivation, the present article collects the scattered tools into a dashboard view, briefly introduces their usage and algorithms, evaluates the accuracy, compares the performance, and provides informative comments for selection.

Fuzzy model-based servo and model following control for nonlinear systems.

PubMed

Ohtake, Hiroshi; Tanaka, Kazuo; Wang, Hua O

2009-12-01

This correspondence presents servo and nonlinear model following controls for a class of nonlinear systems using the Takagi-Sugeno fuzzy model-based control approach. First, the construction method of the augmented fuzzy system for continuous-time nonlinear systems is proposed by differentiating the original nonlinear system. Second, the dynamic fuzzy servo controller and the dynamic fuzzy model following controller, which can make outputs of the nonlinear system converge to target points and to outputs of the reference system, respectively, are introduced. Finally, the servo and model following controller design conditions are given in terms of linear matrix inequalities. Design examples illustrate the utility of this approach.

Nonlinear microscopy of collagen fibers

NASA Astrophysics Data System (ADS)

Strupler, M.; Pena, A.-M.; Hernest, M.; Tharaux, P.-L.; Fabre, A.; Marchal-Somme, J.; Crestani, B.; Débarre, D.; Martin, J.-L.; Beaurepaire, E.; Schanne-Klein, M.-C.

2007-02-01

We used intrinsic Second Harmonic Generation (SHG) by fibrillar collagen to visualize the three-dimensional architecture of collagen fibrosis at the micrometer scale using laser scanning nonlinear microscopy. We showed that SHG signals are highly specific to fibrillar collagen and provide a sensitive probe of the micrometer-scale structural organization of collagen in tissues. Moreover, recording simultaneously other nonlinear optical signals in a multimodal setup, we visualized the tissue morphology using Two-Photon Excited Fluorescence (2PEF) signals from endogenous chromophores such as NADH or elastin. We then compared different methods to determine accurate indexes of collagen fibrosis using nonlinear microscopy, given that most collagen fibrils are smaller than the microscope resolution and that second harmonic generation is a coherent process. In order to define a robust method to process our three-dimensional images, we either calculated the fraction of the images occupied by a significant SHG signal, or averaged SHG signal intensities. We showed that these scores provide an estimation of the extension of renal and pulmonary fibrosis in murine models, and that they clearly sort out the fibrotic mice.

Geometric Theory of Reduction of Nonlinear Control Systems

NASA Astrophysics Data System (ADS)

Elkin, V. I.

2018-02-01

The foundations of a differential geometric theory of nonlinear control systems are described on the basis of categorical concepts (isomorphism, factorization, restrictions) by analogy with classical mathematical theories (of linear spaces, groups, etc.).

New solitary wave solutions of the time-fractional Cahn-Allen equation via the improved (G'/G)-expansion method

NASA Astrophysics Data System (ADS)

Batool, Fiza; Akram, Ghazala

2018-05-01

An improved (G'/G)-expansion method is proposed for extracting more general solitary wave solutions of the nonlinear fractional Cahn-Allen equation. The temporal fractional derivative is taken in the sense of Jumarie's fractional derivative. The results of this article are generalized and extended version of previously reported solutions.

Multiple model self-tuning control for a class of nonlinear systems

NASA Astrophysics Data System (ADS)

Huang, Miao; Wang, Xin; Wang, Zhenlei

2015-10-01

This study develops a novel nonlinear multiple model self-tuning control method for a class of nonlinear discrete-time systems. An increment system model and a modified robust adaptive law are proposed to expand the application range, thus eliminating the assumption that either the nonlinear term of the nonlinear system or its differential term is global-bounded. The nonlinear self-tuning control method can address the situation wherein the nonlinear system is not subject to a globally uniformly asymptotically stable zero dynamics by incorporating the pole-placement scheme. A novel, nonlinear control structure based on this scheme is presented to improve control precision. Stability and convergence can be confirmed when the proposed multiple model self-tuning control method is applied. Furthermore, simulation results demonstrate the effectiveness of the proposed method.

On Partial Fraction Decompositions by Repeated Polynomial Divisions

ERIC Educational Resources Information Center

Man, Yiu-Kwong

2017-01-01

We present a method for finding partial fraction decompositions of rational functions with linear or quadratic factors in the denominators by means of repeated polynomial divisions. This method does not involve differentiation or solving linear equations for obtaining the unknown partial fraction coefficients, which is very suitable for either…

Partial Fractions via Calculus

ERIC Educational Resources Information Center

Bauldry, William C.

2018-01-01

The standard technique taught in calculus courses for partial fraction expansions uses undetermined coefficients to generate a system of linear equations; we present a derivative-based technique that calculus and differential equations instructors can use to reinforce connections to calculus. Simple algebra shows that we can use the derivative to…

Solution of time fractional Black-Scholes European option pricing equation arising in financial market

NASA Astrophysics Data System (ADS)

Ravi Kanth, A. S. V.; Aruna, K.

2016-12-01

In this paper, we present fractional differential transform method (FDTM) and modified fractional differential transform method (MFDTM) for the solution of time fractional Black-Scholes European option pricing equation. The method finds the solution without any discretization, transformation, or restrictive assumptions with the use of appropriate initial or boundary conditions. The efficiency and exactitude of the proposed methods are tested by means of three examples.

Dynamic evolution characteristics of a fractional order hydropower station system

NASA Astrophysics Data System (ADS)

Gao, Xiang; Chen, Diyi; Yan, Donglin; Xu, Beibei; Wang, Xiangyu

2018-01-01

This paper investigates the dynamic evolution characteristics of the hydropower station by introducing the fractional order damping forces. A careful analysis of the dynamic characteristics of the generator shaft system is carried out under different values of fractional order. It turns out the vibration state of the axis coordinates has a certain evolution law with the increase of the fractional order. Significantly, the obtained law exists in the horizontal evolution and vertical evolution of the dynamical behaviors. Meanwhile, some interesting dynamical phenomena were found in this process. The outcomes of this study enrich the nonlinear dynamic theory from the engineering practice of hydropower stations.

Investigation of nonlinear motion simulator washout schemes

NASA Technical Reports Server (NTRS)

Riedel, S. A.; Hofmann, L. G.

1978-01-01

An overview is presented of some of the promising washout schemes which have been devised. The four schemes presented fall into two basic configurations; crossfeed and crossproduct. Various nonlinear modifications further differentiate the four schemes. One nonlinear scheme is discussed in detail. This washout scheme takes advantage of subliminal motions to speed up simulator cab centering. It exploits so-called perceptual indifference thresholds to center the simulator cab at a faster rate whenever the input to the simulator is below the perceptual indifference level. The effect is to reduce the angular and translational simulation motion by comparison with that for the linear washout case. Finally, the conclusions and implications for further research in the area of nonlinear washout filters are presented.

Collective effect of personal behavior induced preventive measures and differential rate of transmission on spread of epidemics

NASA Astrophysics Data System (ADS)

Sagar, Vikram; Zhao, Yi

2017-02-01

In the present work, the effect of personal behavior induced preventive measures is studied on the spread of epidemics over scale free networks that are characterized by the differential rate of disease transmission. The role of personal behavior induced preventive measures is parameterized in terms of variable λ, which modulates the number of concurrent contacts a node makes with the fraction of its neighboring nodes. The dynamics of the disease is described by a non-linear Susceptible Infected Susceptible model based upon the discrete time Markov Chain method. The network mean field approach is generalized to account for the effect of non-linear coupling between the aforementioned factors on the collective dynamics of nodes. The upper bound estimates of the disease outbreak threshold obtained from the mean field theory are found to be in good agreement with the corresponding non-linear stochastic model. From the results of parametric study, it is shown that the epidemic size has inverse dependence on the preventive measures (λ). It has also been shown that the increase in the average degree of the nodes lowers the time of spread and enhances the size of epidemics.

Fractional calculus in hydrologic modeling: A numerical perspective

PubMed Central

Benson, David A.; Meerschaert, Mark M.; Revielle, Jordan

2013-01-01

Fractional derivatives can be viewed either as handy extensions of classical calculus or, more deeply, as mathematical operators defined by natural phenomena. This follows the view that the diffusion equation is defined as the governing equation of a Brownian motion. In this paper, we emphasize that fractional derivatives come from the governing equations of stable Lévy motion, and that fractional integration is the corresponding inverse operator. Fractional integration, and its multi-dimensional extensions derived in this way, are intimately tied to fractional Brownian (and Lévy) motions and noises. By following these general principles, we discuss the Eulerian and Lagrangian numerical solutions to fractional partial differential equations, and Eulerian methods for stochastic integrals. These numerical approximations illuminate the essential nature of the fractional calculus. PMID:23524449

Building Blocks for Reliable Complex Nonlinear Numerical Simulations

NASA Technical Reports Server (NTRS)

Yee, H. C.

2005-01-01

This chapter describes some of the building blocks to ensure a higher level of confidence in the predictability and reliability (PAR) of numerical simulation of multiscale complex nonlinear problems. The focus is on relating PAR of numerical simulations with complex nonlinear phenomena of numerics. To isolate sources of numerical uncertainties, the possible discrepancy between the chosen partial differential equation (PDE) model and the real physics and/or experimental data is set aside. The discussion is restricted to how well numerical schemes can mimic the solution behavior of the underlying PDE model for finite time steps and grid spacings. The situation is complicated by the fact that the available theory for the understanding of nonlinear behavior of numerics is not at a stage to fully analyze the nonlinear Euler and Navier-Stokes equations. The discussion is based on the knowledge gained for nonlinear model problems with known analytical solutions to identify and explain the possible sources and remedies of numerical uncertainties in practical computations.

Hipergeometric solutions to some nonhomogeneous equations of fractional order

NASA Astrophysics Data System (ADS)

Olivares, Jorge; Martin, Pablo; Maass, Fernando

2017-12-01

In this paper a study is performed to the solution of the linear non homogeneous fractional order alpha differential equation equal to I 0(x), where I 0(x) is the modified Bessel function of order zero, the initial condition is f(0)=0 and 0 < alpha < 1. Caputo definition for the fractional derivatives is considered. Fractional derivatives have become important in physical and chemical phenomena as visco-elasticity and visco-plasticity, anomalous diffusion and electric circuits. In particular in this work the values of alpha=1/2, 1/4 and 3/4. are explicitly considered . In these cases Laplace transform is applied, and later the inverse Laplace transform leads to the solutions of the differential equation, which become hypergeometric functions.

Liquid fractionation. Part I: Basic principles and experimental simulations

NASA Astrophysics Data System (ADS)

McBirney, Alexander R.; Baker, Brian H.; Nilson, Robert H.

1985-03-01

A possible explanation for the closely associated magmas of contrasting compositions erupted from many mature volcanic centers can be found in the large differences of density produced by relatively small compositional variations in liquids that evolve by crystallization or melting at the walls of shallow magma chambers. A mechanism of liquid fractionation in which differentiated liquids segragate gravitationally to form compositionally graded columns of magma may surmount the long-standing problem of explaining large volumes of highly evolved liquids that reach advanced degrees of differentiation in times that are too short to be consistent with conventional models of crystal fractionation based on crystal settling. In those types of magmas that decrease in density as they differentiate, a fractionated liquid next to a wall may form a buoyant compositional boundary layer that flows up the wall and accumulates as a separate zone in the upper levels of the reservoir. Magmas that increase in density as they differentiate will have the opposite behavior; they descend along the wall and pond on the floor. Both types of systems can be modeled using simple aqueous solutions and techniques similar to those developed by Chen and Turner (1980). The insights gained through experiments of this kind suggest a number of processes that may be responsible for common types of volcanic behavior and patterns of differentiation in shallow plutons.

Nonlinear Bubble Dynamics And The Effects On Propagation Through Near-Surface Bubble Layers

NASA Astrophysics Data System (ADS)

Leighton, Timothy G.

2004-11-01

Nonlinear bubble dynamics are often viewed as the unfortunate consequence of having to use high acoustic pressure amplitudes when the void fraction in the near-surface oceanic bubble layer is great enough to cause severe attenuation (e.g. >50 dB/m). This is seen as unfortunate since existing models for acoustic propagation in bubbly liquids are based on linear bubble dynamics. However, the development of nonlinear models does more than just allow quantification of the errors associated with the use of linear models. It also offers the possibility of propagation modeling and acoustic inversions which appropriately incorporate the bubble nonlinearity. Furthermore, it allows exploration and quantification of possible nonlinear effects which may be exploited. As a result, high acoustic pressure amplitudes may be desirable even in low void fractions, because they offer opportunities to gain information about the bubble cloud from the nonlinearities, and options to exploit the nonlinearities to enhance communication and sonar in bubbly waters. This paper presents a method for calculating the nonlinear acoustic cross-sections, scatter, attenuations and sound speeds from bubble clouds which may be inhomogeneous. The method allows prediction of the time dependency of these quantities, both because the cloud may vary and because the incident acoustic pulse may have finite and arbitrary time history. The method can be readily adapted for bubbles in other environments (e.g. clouds of interacting bubbles, sediments, structures, in vivo, reverberant conditions etc.). The possible exploitation of bubble acoustics by marine mammals, and for sonar enhancement, is explored.

Differential behavior of amino-imino constitutional isomers in nonlinear optical processes.

PubMed

Latorre, Sonia; Moreira, Ibério de P R; Villacampa, Belén; Julià, Lluís; Velasco, Dolores; Bofill, Josep Maria; López-Calahorra, Francisco

2010-03-15

A detailed study of the "blocked" amino-imino tautomers derived from N-acridine-substituted 2-aminobenzothiazole--and their effect on the nonlinear optical response--is presented. The synthesis, characterization, and nonlinear optical properties of these frozen tautomers, namely, N-methyl-N-(2-nitroacridin-6-yl)-2-aminobenzothia-zole and 3-methyl-N-(7-nitroacridin-3-yl)-2-iminobenzothiazole, are reported. A theoretical model based on valence-bond theory is also proposed and used to analyze the effects of the nuclear configuration corresponding to each frozen tautomer structure. In the present case, the aromatic form and the allylic-anion-like system of the -N-C-N- group inherent to each isomer are crucial for understanding and analyzing the different responses of each "blocked" tautomer.

Modified Method of Simplest Equation Applied to the Nonlinear Schrödinger Equation

NASA Astrophysics Data System (ADS)

Vitanov, Nikolay K.; Dimitrova, Zlatinka I.

2018-03-01

We consider an extension of the methodology of the modified method of simplest equation to the case of use of two simplest equations. The extended methodology is applied for obtaining exact solutions of model nonlinear partial differential equations for deep water waves: the nonlinear Schrödinger equation. It is shown that the methodology works also for other equations of the nonlinear Schrödinger kind.

Differential geometric methods in system theory.

NASA Technical Reports Server (NTRS)

Brockett, R. W.

1971-01-01

Discussion of certain problems in system theory which have been or might be solved using some basic concepts from differential geometry. The problems considered involve differential equations, controllability, optimal control, qualitative behavior, stochastic processes, and bilinear systems. The main goal is to extend the essentials of linear theory to some nonlinear classes of problems.

Davidenko’s Method for the Solution of Nonlinear Operator Equations.

DTIC Science & Technology

NONLINEAR DIFFERENTIAL EQUATIONS, NUMERICAL INTEGRATION), OPERATORS(MATHEMATICS), BANACH SPACE , MAPPING (TRANSFORMATIONS), NUMERICAL METHODS AND PROCEDURES, INTEGRALS, SET THEORY, CONVERGENCE, MATRICES(MATHEMATICS)

Prolongation structures of nonlinear evolution equations

NASA Technical Reports Server (NTRS)

Wahlquist, H. D.; Estabrook, F. B.

1975-01-01

A technique is developed for systematically deriving a 'prolongation structure' - a set of interrelated potentials and pseudopotentials - for nonlinear partial differential equations in two independent variables. When this is applied to the Korteweg-de Vries equation, a new infinite set of conserved quantities is obtained. Known solution techniques are shown to result from the discovery of such a structure: related partial differential equations for the potential functions, linear 'inverse scattering' equations for auxiliary functions, Backlund transformations. Generalizations of these techniques will result from the use of irreducible matrix representations of the prolongation structure.

The existence of solutions of q-difference-differential equations.

PubMed

Wang, Xin-Li; Wang, Hua; Xu, Hong-Yan

2016-01-01

By using the Nevanlinna theory of value distribution, we investigate the existence of solutions of some types of non-linear q-difference differential equations. In particular, we generalize the Rellich-Wittich-type theorem and Malmquist-type theorem about differential equations to the case of q-difference differential equations (system).

Description of a computer program and numerical techniques for developing linear perturbation models from nonlinear systems simulations

NASA Technical Reports Server (NTRS)

Dieudonne, J. E.

1978-01-01

A numerical technique was developed which generates linear perturbation models from nonlinear aircraft vehicle simulations. The technique is very general and can be applied to simulations of any system that is described by nonlinear differential equations. The computer program used to generate these models is discussed, with emphasis placed on generation of the Jacobian matrices, calculation of the coefficients needed for solving the perturbation model, and generation of the solution of the linear differential equations. An example application of the technique to a nonlinear model of the NASA terminal configured vehicle is included.

Power-law spatial dispersion from fractional Liouville equation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tarasov, Vasily E.

2013-10-15

A microscopic model in the framework of fractional kinetics to describe spatial dispersion of power-law type is suggested. The Liouville equation with the Caputo fractional derivatives is used to obtain the power-law dependence of the absolute permittivity on the wave vector. The fractional differential equations for electrostatic potential in the media with power-law spatial dispersion are derived. The particular solutions of these equations for the electric potential of point charge in this media are considered.

Discrete time learning control in nonlinear systems

NASA Technical Reports Server (NTRS)

Longman, Richard W.; Chang, Chi-Kuang; Phan, Minh

1992-01-01

In this paper digital learning control methods are developed primarily for use in single-input, single-output nonlinear dynamic systems. Conditions for convergence of the basic form of learning control based on integral control concepts are given, and shown to be satisfied by a large class of nonlinear problems. It is shown that it is not the gross nonlinearities of the differential equations that matter in the convergence, but rather the much smaller nonlinearities that can manifest themselves during the short time interval of one sample time. New algorithms are developed that eliminate restrictions on the size of the learning gain, and on knowledge of the appropriate sign of the learning gain, for convergence to zero error in tracking a feasible desired output trajectory. It is shown that one of the new algorithms can give guaranteed convergence in the presence of actuator saturation constraints, and indicate when the requested trajectory is beyond the actuator capabilities.

Nonlinear bleaching, absorption, and scattering of 532-nm-irradiated plasmonic nanoparticles

DOE Office of Scientific and Technical Information (OSTI.GOV)

Liberman, V.; Sworin, M.; Kingsborough, R. P.

2013-02-07

Single-pulse irradiation of Au and Ag suspensions of nanospheres and nanodisks with 532-nm 4-ns pulses has identified complex optical nonlinearities while minimizing material damage. For all materials tested, we observe competition between saturable absorption (SA) and reverse SA (RSA), with RSA behavior dominating for intensities above {approx}50 MW/cm{sup 2}. Due to reduced laser damage in single-pulse experiments, the observed intrinsic nonlinear absorption coefficients are the highest reported to date for Au nanoparticles. We find size dependence to the nonlinear absorption enhancement for Au nanoparticles, peaking in magnitude for 80-nm nanospheres and falling off at larger sizes. The nonlinear absorption coefficientsmore » for Au and Ag spheres are comparable in magnitude. On the other hand, the nonlinear absorption for Ag disks, when corrected for volume fraction, is several times higher. These trends in nonlinear absorption are correlated to local electric field enhancement through quasi-static mean-field theory. Through variable size aperture measurements, we also separate nonlinear scattering from nonlinear absorption. For all materials tested, we find that nonlinear scattering is highly directional and that its magnitude is comparable to that of nonlinear absorption. These results indicate methods to improve the efficacy of plasmonic nanoparticles as optical limiters in pulsed laser systems.« less

The nonlinear differential equations governing a hierarchy of self-exciting coupled Faraday-disk homopolar dynamos

NASA Astrophysics Data System (ADS)

Hide, Raymond

1997-02-01

This paper discusses the derivation of the autonomous sets of dimensionless nonlinear ordinary differential equations (ODE's) that govern the behaviour of a hierarchy of related electro-mechanical self-exciting Faraday-disk homopolar dynamo systems driven by steady mechanical couples. Each system comprises N interacting units which could be arranged in a ring or lattice. Within each unit and connected in parallel or in series with the coil are electric motors driven into motion by the dynamo, all having linear characteristics, so that nonlinearity arises entirely through the coupling between components. By introducing simple extra terms into the equations it is possible to represent biasing effects arising from impressed electromotive forces due to thermoelectric or chemical processes and from the presence of ambient magnetic fields. Dissipation in the system is due not only to ohmic heating but also to mechanical friction in the disk and the motors, with the latter agency, no matter how weak, playing an unexpectedly crucial rôle in the production of régimes of chaotic behaviour. This has already been demonstrated in recent work on a case of a single unit incorporating just one series motor, which is governed by a novel autonomous set of nonlinear ODE's with three time-dependent variables and four control parameters. It will be of mathematical as well as geophysical and astrophysical interest to investigate systematically phase and amplitude locking and other types of behaviour in the more complicated cases that arise when N > 1, which can typically involve up to 6 N dependent variables and 19 N-5 control parameters. Even the simplest members of the hierarchy, with N as low as 1, 2 or 3, could prove useful as physically-realistic low-dimensional models in theoretical studies of fluctuating stellar and planetary magnetic fields. Geomagnetic polarity reversals could be affected by the presence of the Earth's solid metallic inner core, driven like an electric motor

Sampled-data chain-observer design for a class of delayed nonlinear systems

NASA Astrophysics Data System (ADS)

Kahelras, M.; Ahmed-Ali, T.; Giri, F.; Lamnabhi-Lagarrigue, F.

2018-05-01

The problem of observer design is addressed for a class of triangular nonlinear systems with not-necessarily small delay and sampled output measurements. One more difficulty is that the system state matrix is dependent on the un-delayed output signal which is not accessible to measurement, making existing observers inapplicable. A new chain observer, composed of m elementary observers in series, is designed to compensate for output sampling and arbitrary large delays. The larger the time-delay the larger the number m. Each elementary observer includes an output predictor that is conceived to compensate for the effects of output sampling and a fractional delay. The predictors are defined by first-order ordinary differential equations (ODEs) much simpler than those of existing predictors which involve both output and state predictors. Using a small gain type analysis, sufficient conditions for the observer to be exponentially convergent are established in terms of the minimal number m of elementary observers and the maximum sampling interval.

New non-linear model of groundwater recharge: Inclusion of memory, heterogeneity and visco-elasticity