GDSCalc: A Web-Based Application for Evaluating Discrete Graph Dynamical Systems
Elmeligy Abdelhamid, Sherif H.; Kuhlman, Chris J.; Marathe, Madhav V.; Mortveit, Henning S.; Ravi, S. S.
2015-01-01
Discrete dynamical systems are used to model various realistic systems in network science, from social unrest in human populations to regulation in biological networks. A common approach is to model the agents of a system as vertices of a graph, and the pairwise interactions between agents as edges. Agents are in one of a finite set of states at each discrete time step and are assigned functions that describe how their states change based on neighborhood relations. Full characterization of state transitions of one system can give insights into fundamental behaviors of other dynamical systems. In this paper, we describe a discrete graph dynamical systems (GDSs) application called GDSCalc for computing and characterizing system dynamics. It is an open access system that is used through a web interface. We provide an overview of GDS theory. This theory is the basis of the web application; i.e., an understanding of GDS provides an understanding of the software features, while abstracting away implementation details. We present a set of illustrative examples to demonstrate its use in education and research. Finally, we compare GDSCalc with other discrete dynamical system software tools. Our perspective is that no single software tool will perform all computations that may be required by all users; tools typically have particular features that are more suitable for some tasks. We situate GDSCalc within this space of software tools. PMID:26263006
NASA Astrophysics Data System (ADS)
Wuensche, Andrew
DDLab is interactive graphics software for creating, visualizing, and analyzing many aspects of Cellular Automata, Random Boolean Networks, and Discrete Dynamical Networks in general and studying their behavior, both from the time-series perspective — space-time patterns, and from the state-space perspective — attractor basins. DDLab is relevant to research, applications, and education in the fields of complexity, self-organization, emergent phenomena, chaos, collision-based computing, neural networks, content addressable memory, genetic regulatory networks, dynamical encryption, generative art and music, and the study of the abstract mathematical/physical/dynamical phenomena in their own right.
Systoles in discrete dynamical systems
NASA Astrophysics Data System (ADS)
Fernandes, Sara; Grácio, Clara; Ramos, Carlos Correia
2013-01-01
The fruitful relationship between Geometry and Graph Theory has been explored by several authors benefiting also the Theory of discrete dynamical systems seen as Markov chains in graphs. In this work we will further explore the relation between these areas, giving a geometrical interpretation of notions from dynamical systems. In particular, we relate the topological entropy with the systole, here defined in the context of discrete dynamical systems. We show that for continuous interval maps the systole is trivial; however, for the class of interval maps with one discontinuity point the systole acquires relevance from the point of view of the dynamical behavior. Moreover, we define the geodesic length spectrum associated to a Markov interval map and we compute the referred spectrum in several examples.
Discreteness effects in population dynamics
NASA Astrophysics Data System (ADS)
Guevara Hidalgo, Esteban; Lecomte, Vivien
2016-05-01
We analyse numerically the effects of small population size in the initial transient regime of a simple example population dynamics. These effects play an important role for the numerical determination of large deviation functions of additive observables for stochastic processes. A method commonly used in order to determine such functions is the so-called cloning algorithm which in its non-constant population version essentially reduces to the determination of the growth rate of a population, averaged over many realizations of the dynamics. However, the averaging of populations is highly dependent not only on the number of realizations of the population dynamics, and on the initial population size but also on the cut-off time (or population) considered to stop their numerical evolution. This may result in an over-influence of discreteness effects at initial times, caused by small population size. We overcome these effects by introducing a (realization-dependent) time delay in the evolution of populations, additional to the discarding of the initial transient regime of the population growth where these discreteness effects are strong. We show that the improvement in the estimation of the large deviation function comes precisely from these two main contributions.
Terminal Dynamics Approach to Discrete Event Systems
NASA Technical Reports Server (NTRS)
Zak, Michail; Meyers, Ronald
1995-01-01
This paper presents and discusses a mathematical formalism for simulation of discrete event dynamic (DED)-a special type of 'man-made' systems to serve specific purposes of information processing. The main objective of this work is to demonstrate that the mathematical formalism for DED can be based upon a terminal model of Newtonian dynamics which allows one to relax Lipschitz conditions at some discrete points.!.
Dynamic discretization method for solving Kepler's equation
NASA Astrophysics Data System (ADS)
Feinstein, Scott A.; McLaughlin, Craig A.
2006-09-01
Kepler’s equation needs to be solved many times for a variety of problems in Celestial Mechanics. Therefore, computing the solution to Kepler’s equation in an efficient manner is of great importance to that community. There are some historical and many modern methods that address this problem. Of the methods known to the authors, Fukushima’s discretization technique performs the best. By taking more of a system approach and combining the use of discretization with the standard computer science technique known as dynamic programming, we were able to achieve even better performance than Fukushima. We begin by defining Kepler’s equation for the elliptical case and describe existing solution methods. We then present our dynamic discretization method and show the results of a comparative analysis. This analysis will demonstrate that, for the conditions of our tests, dynamic discretization performs the best.
Dynamic Web Pages: Performance Impact on Web Servers.
ERIC Educational Resources Information Center
Kothari, Bhupesh; Claypool, Mark
2001-01-01
Discussion of Web servers and requests for dynamic pages focuses on experimentally measuring and analyzing the performance of the three dynamic Web page generation technologies: CGI, FastCGI, and Servlets. Develops a multivariate linear regression model and predicts Web server performance under some typical dynamic requests. (Author/LRW)
Discrete dislocation dynamics simulations in a cylinder
NASA Astrophysics Data System (ADS)
Li, Maosheng; Gao, Chan; Xu, Jianing
2015-02-01
Mechanical properties of material are closely related to the motion of dislocations, and predicting the interactions and resulting collective motion of dislocations is a major task in understanding and modelling plastically deforming materials. A discrete dislocation dynamics model is used to describe the orientation substructure within the microstructure. Discrete dislocation dynamics simulations in three dimensions have been used to examine the role of dislocation multiplication and mobility on the plasticity in small samples under uniaxial compression. In this paper we describe the application of the dislocation dynamics simulations in a cylindrical geometry. The boundary conditions for the simulation were estimated from the distribution of the geometrically necessary dislocation density which was obtained from the orientation map. Numerical studies benchmark could validate the accuracy of the algorithms and the importance of handling the singularity correctly. The results of the simulation explain the formation of the experimentally observed substructure.
Time Discretization Approach to Dynamic Localization Conditions
NASA Astrophysics Data System (ADS)
Papp, E.
An alternative wavefunction to the description of the dynamic localization of a charged particle moving on a one-dimensional lattice under the influence of a periodic time dependent electric field is written down. For this purpose the method of characteristics such as applied by Dunlap and Kenkre [Phys. Rev. B 34, 3625 (1986)] has been modified by using a different integration variable. Handling this wavefunction one is faced with the selection of admissible time values. This results in a conditionally exactly solvable problem, now by accounting specifically for the implementation of a time discretization working in conjunction with a related dynamic localization condition. In addition, one resorts to the strong field limit, which amounts to replace, to leading order, the large order zeros of the Bessel function J0(z), used before in connection with the cosinusoidal modulation, by integral multiples of π. Here z stands for the ratio between the field amplitude and the frequency. The modulation function of the electric field vanishes on the nodal points of the time grid, which stands for an effective field-free behavior. This opens the way to propose quickly tractable dynamic localization conditions for arbitrary periodic modulations. We have also found that the present time discretization approach produces the minimization of the mean square displacement characterizing the usual exact wavefunction. Other realizations and comparisons have also been presented.
Discrete dynamics and non-Markovianity
NASA Astrophysics Data System (ADS)
Luoma, Kimmo; Piilo, Jyrki
2016-06-01
We study discrete quantum dynamics where a single evolution step consists of unitary system transformation followed by decoherence via coupling to an environment. Often, non-Markovian memory effects are attributed to structured environments, whereas, here, we take a more general approach within a discrete setting. In addition of controlling the structure of the environment, we are interested in how local unitaries on the open system allow the appearance and control of memory effects. Our first simple qubit model where local unitary is followed by dephasing illustrates how memory effects arise, despite having no structure in the environment the system is coupled with. We, then, elaborate on this observation by constructing a model for an open quantum walk where the unitary coin and transfer operation is augmented with the dephasing of the coin. The results demonstrate tha,t in the limit of strong dephasing within each evolution step, the combined coin-position open system always displays memory effects, and their quantities are independent of the structure of the environment. Our construction makes possible an experimentally realizable open quantum walk with photons exhibiting non-Markovian features.
Non-Lipschitz Dynamics Approach to Discrete Event Systems
NASA Technical Reports Server (NTRS)
Zak, M.; Meyers, R.
1995-01-01
This paper presents and discusses a mathematical formalism for simulation of discrete event dynamics (DED) - a special type of 'man- made' system designed to aid specific areas of information processing. A main objective is to demonstrate that the mathematical formalism for DED can be based upon the terminal model of Newtonian dynamics which allows one to relax Lipschitz conditions at some discrete points.
Reducing Neuronal Networks to Discrete Dynamics
Terman, David; Ahn, Sungwoo; Wang, Xueying; Just, Winfried
2008-01-01
We consider a general class of purely inhibitory and excitatory-inhibitory neuronal networks, with a general class of network architectures, and characterize the complex firing patterns that emerge. Our strategy for studying these networks is to first reduce them to a discrete model. In the discrete model, each neuron is represented as a finite number of states and there are rules for how a neuron transitions from one state to another. In this paper, we rigorously demonstrate that the continuous neuronal model can be reduced to the discrete model if the intrinsic and synaptic properties of the cells are chosen appropriately. In a companion paper [1], we analyze the discrete model. PMID:18443649
Dynamics of discrete screw dislocations on glide directions
NASA Astrophysics Data System (ADS)
Alicandro, R.; De Luca, L.; Garroni, A.; Ponsiglione, M.
2016-07-01
We consider a simple discrete model for screw dislocations in crystals. Using a variational discrete scheme we study the motion of a configuration of dislocations toward low energy configurations. We deduce an effective fully overdamped dynamics that follows the maximal dissipation criterion introduced in Cermelli and Gurtin (1999) and predicts motion along the glide directions of the crystal.
Dynamics of a discrete auroral arc
NASA Technical Reports Server (NTRS)
Bruening, K.; Goertz, C. K.
1986-01-01
Porcupine Flight 4 data were used to determine the field-aligned currents associated with a southward moving discrete auroral arc in the postmidnight sector. Three different methods were used for determining the field-aligned current which should give identical results if the arcs are quasi-stationary and no parallel electric field exists between the payload and the dynamo region of the ionosphere. As long as the rocket is above the arc, the three methods agree. The integral of precipitating electron flux, the local magnetic field perturbations, and the divergence of the horizontal Pedersen current all indicate an upward current of 5 + or - 3 microamperes/sq m. Immediately north of the arc a strong downward current of about 10-20 microamperes/sq m is detected. The magnitude, however, is not well known because the rocket's velocity relative to the arc cannot be clearly established. Further north of the southward moving arc, the two methods that can be applied (magnetic field perturbations and divergence of the horizontal Pedersen current) yield contradictory results not only about the magnitude of the current but also about the direction of the current. It is suggested that this discrepancy is due to time-dependent electric field.
Dynamics of a discrete auroral arc
NASA Astrophysics Data System (ADS)
Bruening, K.; Goertz, C. K.
1986-06-01
Porcupine Flight 4 data were used to determine the field-aligned currents associated with a southward moving discrete auroral arc in the postmidnight sector. Three different methods were used for determining the field-aligned current which should give identical results if the arcs are quasi-stationary and no parallel electric field exists between the payload and the dynamo region of the ionosphere. As long as the rocket is above the arc, the three methods agree. The integral of precipitating electron flux, the local magnetic field perturbations, and the divergence of the horizontal Pedersen current all indicate an upward current of 5 + or - 3 microamperes/sq m. Immediately north of the arc a strong downward current of about 10-20 microamperes/sq m is detected. The magnitude, however, is not well known because the rocket's velocity relative to the arc cannot be clearly established. Further north of the southward moving arc, the two methods that can be applied (magnetic field perturbations and divergence of the horizontal Pedersen current) yield contradictory results not only about the magnitude of the current but also about the direction of the current. It is suggested that this discrepancy is due to time-dependent electric field.
Symbolic dynamics and the discrete variational principle
NASA Astrophysics Data System (ADS)
Dullin, H. R.
1998-11-01
We show how to construct symbolic dynamics for the class of 2d-dimensional twist mappings generated by piecewise strictly convex/concave generating functions. The method is constructive and gives an efficient way to find all periodic orbits of these high-dimensional symplectic mappings. It is illustrated with the cardioid and the stadium billiard.
On the Importance of the Dynamics of Discretizations
NASA Technical Reports Server (NTRS)
Sweby, Peter K.; Yee, H. C.; Rai, ManMohan (Technical Monitor)
1995-01-01
It has been realized recently that the discrete maps resulting from numerical discretizations of differential equations can possess asymptotic dynamical behavior quite different from that of the original systems. This is the case not only for systems of Ordinary Differential Equations (ODEs) but in a more complicated manner for Partial Differential Equations (PDEs) used to model complex physics. The impact of the modified dynamics may be mild and even not observed for some numerical methods. For other classes of discretizations the impact may be pronounced, but not always obvious depending on the nonlinear model equations, the time steps, the grid spacings and the initial conditions. Non-convergence or convergence to periodic solutions might be easily recognizable but convergence to incorrect but plausible solutions may not be so obvious - even for discretized parameters within the linearized stability constraint. Based on our past four years of research, we will illustrate some of the pathology of the dynamics of discretizations, its possible impact and the usage of these schemes for model nonlinear ODEs, convection-diffusion equations and grid adaptations.
Global asymptotic coherence in discrete dynamical systems
Earn, David J. D.; Levin, Simon A.
2006-01-01
Spatial synchrony (coherence) in dynamical systems is of both theoretical and applied importance. We address this problem for a generalization of coupled map lattices (CMLs). In the systems we study, which we term “meta-CMLs,” the map at each lattice point may be multidimensional (corresponding, for example, to multispecies ecological systems in which all species have the same dispersal pattern). Most previous work on coherence of CMLs has focused on local stability. Here, we prove a global theorem that provides a useful sufficient condition guaranteeing decay of incoherence in meta-CMLs regardless of initial conditions and regardless of the nature of the attractors of the system. This result facilitates useful analyses of a variety of applied problems, including conservation of endangered species and eradication of pests or infectious diseases. PMID:16537469
Analyzing neuronal networks using discrete-time dynamics
NASA Astrophysics Data System (ADS)
Ahn, Sungwoo; Smith, Brian H.; Borisyuk, Alla; Terman, David
2010-05-01
We develop mathematical techniques for analyzing detailed Hodgkin-Huxley like models for excitatory-inhibitory neuronal networks. Our strategy for studying a given network is to first reduce it to a discrete-time dynamical system. The discrete model is considerably easier to analyze, both mathematically and computationally, and parameters in the discrete model correspond directly to parameters in the original system of differential equations. While these networks arise in many important applications, a primary focus of this paper is to better understand mechanisms that underlie temporally dynamic responses in early processing of olfactory sensory information. The models presented here exhibit several properties that have been described for olfactory codes in an insect’s Antennal Lobe. These include transient patterns of synchronization and decorrelation of sensory inputs. By reducing the model to a discrete system, we are able to systematically study how properties of the dynamics, including the complex structure of the transients and attractors, depend on factors related to connectivity and the intrinsic and synaptic properties of cells within the network.
Analyzing Neuronal Networks Using Discrete-Time Dynamics
Ahn, Sungwoo; Smith, Brian H.; Borisyuk, Alla; Terman, David
2010-01-01
We develop mathematical techniques for analyzing detailed Hodgkin-Huxley like models for excitatory-inhibitory neuronal networks. Our strategy for studying a given network is to first reduce it to a discrete-time dynamical system. The discrete model is considerably easier to analyze, both mathematically and computationally, and parameters in the discrete model correspond directly to parameters in the original system of differential equations. While these networks arise in many important applications, a primary focus of this paper is to better understand mechanisms that underlie temporally dynamic responses in early processing of olfactory sensory information. The models presented here exhibit several properties that have been described for olfactory codes in an insect's Antennal Lobe. These include transient patterns of synchronization and decorrelation of sensory inputs. By reducing the model to a discrete system, we are able to systematically study how properties of the dynamics, including the complex structure of the transients and attractors, depend on factors related to connectivity and the intrinsic and synaptic properties of cells within the network. PMID:20454529
Some remarks on discrete physics as an ultimate dynamical theory
Noyes, H.P.
1995-10-01
The standard model of quarks and leptons currently fails to meet `t Hooft`s criterion for an ``Ultimate Dynamical Theory`` in that it contains 19 parameters which must be taken from experiment. Evaluating discrete physics in the same way we find that it requires 17 parameters and has already succeeded in computing 7 of them. While we are confident that the rest can also be computed, the very concept of an ultimate theory is incompatible with our attitude toward physics.
Pinning synchronization of discrete dynamical networks with delay coupling
NASA Astrophysics Data System (ADS)
Cheng, Ranran; Peng, Mingshu; Zuo, Jun
2016-05-01
The purpose of this paper is to investigate the pinning synchronization analysis for nonlinear coupled delayed discrete dynamical networks with the identical or nonidentical topological structure. Based on the Lyapunov stability theory, pinning control method and linear matrix inequalities, several adaptive synchronization criteria via two kinds of pinning control method are obtained. Two examples based on Rulkov chaotic system are included to illustrate the effectiveness and verification of theoretical analysis.
Dynamic modelling and analysis of space webs
NASA Astrophysics Data System (ADS)
Yu, Yang; Baoyin, HeXi; Li, JunFeng
2011-04-01
Future space missions demand operations on large flexible structures, for example, space webs, the lightweight cable nets deployable in space, which can serve as platforms for very large structures or be used to capture orbital objects. The interest in research on space webs is likely to increase in the future with the development of promising applications such as Furoshiki sat-ellite of JAXA, Robotic Geostationary Orbit Restorer (ROGER) of ESA and Grapple, Retrieve And Secure Payload (GRASP) of NASA. Unlike high-tensioned nets in civil engineering, space webs may be low-tensioned or tensionless, and extremely flexible, owing to the microgravity in the orbit and the lack of support components, which may cause computational difficulties. Mathematical models are necessary in the analysis of space webs, especially in the conceptual design and evaluation for prototypes. A full three-dimensional finite element (FE) model was developed in this work. Trivial truss elements were adopted to reduce the computational complexity. Considering cable is a compression-free material and its tensile stiffness is also variable, we introduced the cable material constitutive relationship to work out an accurate and feasible model for prototype analysis and design. In the static analysis, the stress distribution and global deformation of the webs were discussed to get access to the knowledge of strength of webs with different types of meshes. In the dynamic analysis, special attention was paid to the impact problem. The max stress and global deformation were investigated. The simulation results indicate the interesting phenomenon which may be worth further research.
Applications of Discrete Molecular Dynamics in biology and medicine.
Proctor, Elizabeth A; Dokholyan, Nikolay V
2016-04-01
Discrete Molecular Dynamics (DMD) is a physics-based simulation method using discrete energetic potentials rather than traditional continuous potentials, allowing microsecond time scale simulations of biomolecular systems to be performed on personal computers rather than supercomputers or specialized hardware. With the ongoing explosion in processing power even in personal computers, applications of DMD have similarly multiplied. In the past two years, researchers have used DMD to model structures of disease-implicated protein folding intermediates, study assembly of protein complexes, predict protein-protein binding conformations, engineer rescue mutations in disease-causative protein mutants, design a protein conformational switch to control cell signaling, and describe the behavior of polymeric dispersants for environmental cleanup of oil spills, among other innovative applications. PMID:26638022
Survey of recent work on the analysis of discretely attached corrugated shear webs.
NASA Technical Reports Server (NTRS)
Libove, C.
1972-01-01
Summary and progress report of more recent work by the author and his students on the theoretical analysis of stiffness, stresses, and deformations of corrugated shear webs with discrete, rather than continuous, attachment between the ends of the corrugations and the surrounding members. Various kinds of discrete attachment are considered, and two kinds of corrugation cross section: the trapezoidal and the curvilinear, the latter having crests and valleys made up of identical circular arcs. The more recent analyses employ the method of minimum total potential energy and the calculus of variations to obtain differential equations and boundary conditions governing the longitudinal variation of various component modes used to describe the deformations of a cross section. They are believed to be more accurate than earlier analyses in that they generally permit more degrees of freedom in the assumed deformations. In particular, they abandon the assumption, characteristic of the earlier analyses, that the straight-line generators of the corrugation remain straight during the shearing of the web.
Discrete dynamic modeling of T cell survival signaling networks
NASA Astrophysics Data System (ADS)
Zhang, Ranran
2009-03-01
Biochemistry-based frameworks are often not applicable for the modeling of heterogeneous regulatory systems that are sparsely documented in terms of quantitative information. As an alternative, qualitative models assuming a small set of discrete states are gaining acceptance. This talk will present a discrete dynamic model of the signaling network responsible for the survival and long-term competence of cytotoxic T cells in the blood cancer T-LGL leukemia. We integrated the signaling pathways involved in normal T cell activation and the known deregulations of survival signaling in leukemic T-LGL, and formulated the regulation of each network element as a Boolean (logic) rule. Our model suggests that the persistence of two signals is sufficient to reproduce all known deregulations in leukemic T-LGL. It also indicates the nodes whose inactivity is necessary and sufficient for the reversal of the T-LGL state. We have experimentally validated several model predictions, including: (i) Inhibiting PDGF signaling induces apoptosis in leukemic T-LGL. (ii) Sphingosine kinase 1 and NFκB are essential for the long-term survival of T cells in T-LGL leukemia. (iii) T box expressed in T cells (T-bet) is constitutively activated in the T-LGL state. The model has identified potential therapeutic targets for T-LGL leukemia and can be used for generating long-term competent CTL necessary for tumor and cancer vaccine development. The success of this model, and of other discrete dynamic models, suggests that the organization of signaling networks has an determining role in their dynamics. Reference: R. Zhang, M. V. Shah, J. Yang, S. B. Nyland, X. Liu, J. K. Yun, R. Albert, T. P. Loughran, Jr., Network Model of Survival Signaling in LGL Leukemia, PNAS 105, 16308-16313 (2008).
Convergence Time towards Periodic Orbits in Discrete Dynamical Systems
San Martín, Jesús; Porter, Mason A.
2014-01-01
We investigate the convergence towards periodic orbits in discrete dynamical systems. We examine the probability that a randomly chosen point converges to a particular neighborhood of a periodic orbit in a fixed number of iterations, and we use linearized equations to examine the evolution near that neighborhood. The underlying idea is that points of stable periodic orbit are associated with intervals. We state and prove a theorem that details what regions of phase space are mapped into these intervals (once they are known) and how many iterations are required to get there. We also construct algorithms that allow our theoretical results to be implemented successfully in practice. PMID:24736594
Convergence time towards periodic orbits in discrete dynamical systems.
San Martín, Jesús; Porter, Mason A
2014-01-01
We investigate the convergence towards periodic orbits in discrete dynamical systems. We examine the probability that a randomly chosen point converges to a particular neighborhood of a periodic orbit in a fixed number of iterations, and we use linearized equations to examine the evolution near that neighborhood. The underlying idea is that points of stable periodic orbit are associated with intervals. We state and prove a theorem that details what regions of phase space are mapped into these intervals (once they are known) and how many iterations are required to get there. We also construct algorithms that allow our theoretical results to be implemented successfully in practice. PMID:24736594
Power-law creep from discrete dislocation dynamics.
Keralavarma, Shyam M; Cagin, T; Arsenlis, A; Benzerga, A Amine
2012-12-28
We report two-dimensional discrete dislocation dynamics simulations of combined dislocation glide and climb leading to "power-law" creep in a model aluminum crystal. The approach fully accounts for matter transport due to vacancy diffusion and its coupling with dislocation motion. The existence of quasiequilibrium or jammed states under the applied creep stresses enables observations of diffusion and climb over time scales relevant to power-law creep. The predictions for the creep rates and stress exponents fall within experimental ranges, indicating that the underlying physics is well captured. PMID:23368581
Experiments of reconstructing discrete atmospheric dynamic models from data (I)
NASA Astrophysics Data System (ADS)
Lin, Zhenshan; Zhu, Yanyu; Deng, Ziwang
1995-03-01
In this paper, we give some experimental results of our study in reconstructing discrete atmospheric dynamic models from data. After a great deal of numerical experiments, we found that the logistic map, x n + 1 = 1- μx {2/n}, could be used in monthly mean temperature prediction when it was approaching the chaotic region, and its predictive results were in reverse states to the practical data. This means that the nonlinear developing behavior of the monthly mean temperature system is bifurcating back into the critical chaotic states from the chaotic ones.
On control of continuous dynamical polysystems in discrete times
Khryashchev, S. M.
2015-03-10
This paper considers control systems with a finite number of control parameters, i.e. dynamical polysystems. It is assumed that control switchings can occur only at certain discrete times. Statistical methods in number theory are used to investigate the controllability of these systems. Existence of control switching times is established through analysis of Diophantine equations. The values of control switching times are found through suitable numerical methods. Certain assertions on the controllability of the polysystems of the class under consideration are proved. Some examples are considered.
Novel coupling scheme to control dynamics of coupled discrete systems
NASA Astrophysics Data System (ADS)
Shekatkar, Snehal M.; Ambika, G.
2015-08-01
We present a new coupling scheme to control spatio-temporal patterns and chimeras on 1-d and 2-d lattices and random networks of discrete dynamical systems. The scheme involves coupling with an external lattice or network of damped systems. When the system network and external network are set in a feedback loop, the system network can be controlled to a homogeneous steady state or synchronized periodic state with suppression of the chaotic dynamics of the individual units. The control scheme has the advantage that its design does not require any prior information about the system dynamics or its parameters and works effectively for a range of parameters of the control network. We analyze the stability of the controlled steady state or amplitude death state of lattices using the theory of circulant matrices and Routh-Hurwitz criterion for discrete systems and this helps to isolate regions of effective control in the relevant parameter planes. The conditions thus obtained are found to agree well with those obtained from direct numerical simulations in the specific context of lattices with logistic map and Henon map as on-site system dynamics. We show how chimera states developed in an experimentally realizable 2-d lattice can be controlled using this scheme. We propose this mechanism can provide a phenomenological model for the control of spatio-temporal patterns in coupled neurons due to non-synaptic coupling with the extra cellular medium. We extend the control scheme to regulate dynamics on random networks and adapt the master stability function method to analyze the stability of the controlled state for various topologies and coupling strengths.
The discrete dynamics of symmetric competition in the plane.
Jiang, H; Rogers, T D
1987-01-01
We consider the generalized Lotka-Volterra two-species system xn + 1 = xn exp(r1(1 - xn) - s1yn) yn + 1 = yn exp(r2(1 - yn) - s2xn) originally proposed by R. M. May as a model for competitive interaction. In the symmetric case that r1 = r2 and s1 = s2, a region of ultimate confinement is found and the dynamics therein are described in some detail. The bifurcations of periodic points of low period are studied, and a cascade of period-doubling bifurcations is indicated. Within the confinement region, a parameter region is determined for the stable Hopf bifurcation of a pair of symmetrically placed period-two points, which imposes a second component of oscillation near the stable cycles. It is suggested that the symmetric competitive model contains much of the dynamical complexity to be expected in any discrete two-dimensional competitive model. PMID:3437226
A Jini-based dynamic service WebGIS model
NASA Astrophysics Data System (ADS)
Xuan, Wenling; Chen, Xiuwan; Huang, Zhaoqiang; Zhao, Gang
2007-06-01
The development of current GIS technology has evolved from single platform GIS system into WebGIS. However, The Geographic Information Services (GIServices) provision and application manner cannot meet the requirement of pervasive computing environment. Jini/JAVA technique, a dynamic distributed architecture for providing spontaneous network of services, might be a tool/solution to improve the GIService performance of current WebGIS. This paper studies and analyses Jini infrastructure and its dynamic service mechanism, designs a new WebGIS architecture with Jini-based dynamic service model. The experiment shows that Jini technique can be integrated into WebGIS and to realize the dynamic services organization and management.
Dynamics of Sequence -Discrete Bacterial Populations Inferred Using Metagenomes
Stevens, Sarah; Bendall, Matthew; Kang, Dongwan; Froula, Jeff; Egan, Rob; Chan, Leong-Keat; Tringe, Susannah; McMahon, Katherine; Malmstrom, Rex
2014-03-14
From a multi-year metagenomic time series of two dissimilar Wisconsin lakes we have assembled dozens of genomes using a novel approach that bins contigs into distinct genome based on sequence composition, e.g. kmer frequencies, and contig coverage patterns at various times points. Next, we investigated how these genomes, which represent sequence-discrete bacterial populations, evolved over time and used the time series to discover the population dynamics. For example, we explored changes in single nucleotide polymorphism (SNP) frequencies as well as patterns of gene gain and loss in multiple populations. Interestingly, SNP diversity was purged at nearly every genome position in some populations during the course of this study, suggesting these populations may have experienced genome-wide selective sweeps. This represents the first direct, time-resolved observations of periodic selection in natural populations, a key process predicted by the ecotype model of bacterial diversification.
Bauschinger effect in thin metal films: Discrete dislocation dynamics study
NASA Astrophysics Data System (ADS)
Davoudi, Kamyar M.; Nicola, Lucia; Vlassak, Joost J.
2014-01-01
The effects of dislocation climb on plastic deformation during loading and unloading are studied using a two-dimensional discrete dislocation dynamics model. Simulations are performed for polycrystalline thin films passivated on both surfaces. Dislocation climb lowers the overall level of the stress inside thin films and reduces the work hardening rate. Climb decreases the density of dislocations in pile-ups and reduces back stresses. These factors result in a smaller Bauschinger effect on unloading compared to simulations without climb. As dislocations continue to climb at the onset of unloading and the dislocation density continues to increase, the initial unloading slope increases with decreasing unloading rate. Because climb disperses dislocations, fewer dislocations are annihilated during unloading, leading to a higher dislocation density at the end of the unloading step.
Displaying R spatial statistics on Google dynamic maps with web applications created by Rwui
2012-01-01
Background The R project includes a large variety of packages designed for spatial statistics. Google dynamic maps provide web based access to global maps and satellite imagery. We describe a method for displaying directly the spatial output from an R script on to a Google dynamic map. Methods This is achieved by creating a Java based web application which runs the R script and then displays the results on the dynamic map. In order to make this method easy to implement by those unfamiliar with programming Java based web applications, we have added the method to the options available in the R Web User Interface (Rwui) application. Rwui is an established web application for creating web applications for running R scripts. A feature of Rwui is that all the code for the web application being created is generated automatically so that someone with no knowledge of web programming can make a fully functional web application for running an R script in a matter of minutes. Results Rwui can now be used to create web applications that will display the results from an R script on a Google dynamic map. Results may be displayed as discrete markers and/or as continuous overlays. In addition, users of the web application may select regions of interest on the dynamic map with mouse clicks and the coordinates of the region of interest will automatically be made available for use by the R script. Conclusions This method of displaying R output on dynamic maps is designed to be of use in a number of areas. Firstly it allows statisticians, working in R and developing methods in spatial statistics, to easily visualise the results of applying their methods to real world data. Secondly, it allows researchers who are using R to study health geographics data, to display their results directly onto dynamic maps. Thirdly, by creating a web application for running an R script, a statistician can enable users entirely unfamiliar with R to run R coded statistical analyses of health geographics
Consequences of symbiosis for food web dynamics.
Kooi, B W; Kuijper, L D J; Kooijman, S A L M
2004-09-01
Basic Lotka-Volterra type models in which mutualism (a type of symbiosis where the two populations benefit both) is taken into account, may give unbounded solutions. We exclude such behaviour using explicit mass balances and study the consequences of symbiosis for the long-term dynamic behaviour of a three species system, two prey and one predator species in the chemostat. We compose a theoretical food web where a predator feeds on two prey species that have a symbiotic relationships. In addition to a species-specific resource, the two prey populations consume the products of the partner population as well. In turn, a common predator forages on these prey populations. The temporal change in the biomass and the nutrient densities in the reactor is described by ordinary differential equations (ODE). Since products are recycled, the dynamics of these abiotic materials must be taken into account as well, and they are described by odes in a similar way as the abiotic nutrients. We use numerical bifurcation analysis to assess the long-term dynamic behaviour for varying degrees of symbiosis. Attractors can be equilibria, limit cycles and chaotic attractors depending on the control parameters of the chemostat reactor. These control parameters that can be experimentally manipulated are the nutrient density of the inflow medium and the dilution rate. Bifurcation diagrams for the three species web with a facultative symbiotic association between the two prey populations, are similar to that of a bi-trophic food chain; nutrient enrichment leads to oscillatory behaviour. Predation combined with obligatory symbiotic prey-interactions has a stabilizing effect, that is, there is stable coexistence in a larger part of the parameter space than for a bi-trophic food chain. However, combined with a large growth rate of the predator, the food web can persist only in a relatively small region of the parameter space. Then, two zero-pair bifurcation points are the organizing centers. In
Constant pressure and temperature discrete-time Langevin molecular dynamics
Grønbech-Jensen, Niels; Farago, Oded
2014-11-21
We present a new and improved method for simultaneous control of temperature and pressure in molecular dynamics simulations with periodic boundary conditions. The thermostat-barostat equations are built on our previously developed stochastic thermostat, which has been shown to provide correct statistical configurational sampling for any time step that yields stable trajectories. Here, we extend the method and develop a set of discrete-time equations of motion for both particle dynamics and system volume in order to seek pressure control that is insensitive to the choice of the numerical time step. The resulting method is simple, practical, and efficient. The method is demonstrated through direct numerical simulations of two characteristic model systems—a one-dimensional particle chain for which exact statistical results can be obtained and used as benchmarks, and a three-dimensional system of Lennard-Jones interacting particles simulated in both solid and liquid phases. The results, which are compared against the method of Kolb and Dünweg [J. Chem. Phys. 111, 4453 (1999)], show that the new method behaves according to the objective, namely that acquired statistical averages and fluctuations of configurational measures are accurate and robust against the chosen time step applied to the simulation.
Constant pressure and temperature discrete-time Langevin molecular dynamics.
Grønbech-Jensen, Niels; Farago, Oded
2014-11-21
We present a new and improved method for simultaneous control of temperature and pressure in molecular dynamics simulations with periodic boundary conditions. The thermostat-barostat equations are built on our previously developed stochastic thermostat, which has been shown to provide correct statistical configurational sampling for any time step that yields stable trajectories. Here, we extend the method and develop a set of discrete-time equations of motion for both particle dynamics and system volume in order to seek pressure control that is insensitive to the choice of the numerical time step. The resulting method is simple, practical, and efficient. The method is demonstrated through direct numerical simulations of two characteristic model systems-a one-dimensional particle chain for which exact statistical results can be obtained and used as benchmarks, and a three-dimensional system of Lennard-Jones interacting particles simulated in both solid and liquid phases. The results, which are compared against the method of Kolb and Dünweg [J. Chem. Phys. 111, 4453 (1999)], show that the new method behaves according to the objective, namely that acquired statistical averages and fluctuations of configurational measures are accurate and robust against the chosen time step applied to the simulation. PMID:25416875
Discrete Dynamic Bayesian Network Analysis of fMRI Data
Burge, John; Lane, Terran; Link, Hamilton; Qiu, Shibin; Clark, Vincent P.
2010-01-01
We examine the efficacy of using discrete Dynamic Bayesian Networks (dDBNs), a data-driven modeling technique employed in machine learning, to identify functional correlations among neuroanatomical regions of interest. Unlike many neuroimaging analysis techniques, this method is not limited by linear and/or Gaussian noise assumptions. It achieves this by modeling the time series of neuroanatomical regions as discrete, as opposed to continuous, random variables with multinomial distributions. We demonstrated this method using an fMRI dataset collected from healthy and demented elderly subjects and identify correlates based on a diagnosis of dementia. The results are validated in three ways. First, the elicited correlates are shown to be robust over leave-one-out cross-validation and, via a Fourier bootstrapping method, that they were not likely due to random chance. Second, the dDBNs identified correlates that would be expected given the experimental paradigm. Third, the dDBN's ability to predict dementia is competitive with two commonly employed machine-learning classifiers: the support vector machine and the Gaussian naïve Bayesian network. We also verify that the dDBN selects correlates based on non-linear criteria. Finally, we provide a brief analysis of the correlates elicited from Buckner et al.'s data that suggests that demented elderly subjects have reduced involvement of entorhinal and occipital cortex and greater involvement of the parietal lobe and amygdala in brain activity compared with healthy elderly (as measured via functional correlations among BOLD measurements). Limitations and extensions to the dDBN method are discussed. PMID:17990301
Discretized torsional dynamics and the folding of an RNA chain.
Fernández, A; Salthú, R; Cendra, H
1999-08-01
The aim of this work is to implement a discrete coarse codification of local torsional states of the RNA chain backbone in order to explore the long-time limit dynamics and ultimately obtain a coarse solution to the RNA folding problem. A discrete representation of the soft-mode dynamics is turned into an algorithm for a rough structure prediction. The algorithm itself is inherently parallel, as it evaluates concurrent folding possibilities by pattern recognition, but it may be implemented in a personal computer as a chain of perturbation-translation-renormalization cycles performed on a binary matrix of local topological constraints. This requires suitable representational tools and a periodic quenching of the dynamics for system renormalization. A binary coding of local topological constraints associated with each structural motif is introduced, with each local topological constraint corresponding to a local torsional state. This treatment enables us to adopt a computation time step far larger than hydrodynamic drag time scales. Accordingly, the solvent is no longer treated as a hydrodynamic drag medium. Instead we incorporate its capacity for forming local conformation-dependent dielectric domains. Each translation of the matrix of local topological constraints (LTM's) depends on the conformation-dependent local dielectric created by a confined solvent. Folding pathways are resolved as transitions between patterns of locally encoded structural signals which change within the 1 ns-100 ms time scale range. These coarse folding pathways are generated by a search at regular intervals for structural patterns in the LTM. Each pattern is recorded as a base-pairing pattern (BPP) matrix, a consensus-evaluation operation subject to a renormalization feedback loop. Since several mutually conflicting consensus evaluations might occur at a given time, the need arises for a probabilistic approach appropriate for an ensemble of RNA molecules. Thus, a statistical dynamics of
An implicit finite element method for discrete dynamic fracture
Jobie M. Gerken
1999-12-01
A method for modeling the discrete fracture of two-dimensional linear elastic structures with a distribution of small cracks subject to dynamic conditions has been developed. The foundation for this numerical model is a plane element formulated from the Hu-Washizu energy principle. The distribution of small cracks is incorporated into the numerical model by including a small crack at each element interface. The additional strain field in an element adjacent to this crack is treated as an externally applied strain field in the Hu-Washizu energy principle. The resulting stiffness matrix is that of a standard plane element. The resulting load vector is that of a standard plane element with an additional term that includes the externally applied strain field. Except for the crack strain field equations, all terms of the stiffness matrix and load vector are integrated symbolically in Maple V so that fully integrated plane stress and plane strain elements are constructed. The crack strain field equations are integrated numerically. The modeling of dynamic behavior of simple structures was demonstrated within acceptable engineering accuracy. In the model of axial and transverse vibration of a beam and the breathing mode of vibration of a thin ring, the dynamic characteristics were shown to be within expected limits. The models dominated by tensile forces (the axially loaded beam and the pressurized ring) were within 0.5% of the theoretical values while the shear dominated model (the transversely loaded beam) is within 5% of the calculated theoretical value. The constant strain field of the tensile problems can be modeled exactly by the numerical model. The numerical results should therefore, be exact. The discrepancies can be accounted for by errors in the calculation of frequency from the numerical results. The linear strain field of the transverse model must be modeled by a series of constant strain elements. This is an approximation to the true strain field, so some
Folding a protein by discretizing its backbone torsional dynamics
NASA Astrophysics Data System (ADS)
Fernández, Ariel
1999-05-01
The aim of this work is to provide a coarse codification of local conformational constraints associated with each folding motif of a peptide chain in order to obtain a rough solution to the protein folding problem. This is accomplished by implementing a discretized version of the soft-mode dynamics on a personal computer (PC). Our algorithm mimics a parallel process as it evaluates concurrent folding possibilities by pattern recognition. It may be implemented in a PC as a sequence of perturbation-translation-renormalization (p-t-r) cycles performed on a matrix of local topological constraints (LTM). This requires suitable representational tools and a periodic quenching of the dynamics required for renormalization. We introduce a description of the peptide chain based on a local discrete variable the values of which label the basins of attraction of the Ramachandran map for each residue. Thus, the local variable indicates the basin in which the torsional coordinates of each residue lie at a given time. In addition, a coding of local topological constraints associated with each secondary and tertiary structural motif is introduced. Our treatment enables us to adopt a computation time step of 81 ps, a value far larger than hydrodynamic drag time scales. Folding pathways are resolved as transitions between patterns of locally encoded structural signals that change within the 10 μs-100 ms time scale range. These coarse folding pathways are generated by the periodic search for structural patterns in the time-evolving LTM. Each pattern is recorded as a contact matrix, an operation subject to a renormalization feedback loop. The validity of our approach is tested vis-a-vis experimentally-probed folding pathways eventually generating tertiary interactions in proteins which recover their active structure under in vitro renaturation conditions. As an illustration, we focus on determining significant folding intermediates and late kinetic bottlenecks that occur within the
Modelling machine ensembles with discrete event dynamical system theory
NASA Technical Reports Server (NTRS)
Hunter, Dan
1990-01-01
Discrete Event Dynamical System (DEDS) theory can be utilized as a control strategy for future complex machine ensembles that will be required for in-space construction. The control strategy involves orchestrating a set of interactive submachines to perform a set of tasks for a given set of constraints such as minimum time, minimum energy, or maximum machine utilization. Machine ensembles can be hierarchically modeled as a global model that combines the operations of the individual submachines. These submachines are represented in the global model as local models. Local models, from the perspective of DEDS theory , are described by the following: a set of system and transition states, an event alphabet that portrays actions that takes a submachine from one state to another, an initial system state, a partial function that maps the current state and event alphabet to the next state, and the time required for the event to occur. Each submachine in the machine ensemble is presented by a unique local model. The global model combines the local models such that the local models can operate in parallel under the additional logistic and physical constraints due to submachine interactions. The global model is constructed from the states, events, event functions, and timing requirements of the local models. Supervisory control can be implemented in the global model by various methods such as task scheduling (open-loop control) or implementing a feedback DEDS controller (closed-loop control).
Dynamic hybrid algorithms for MAP inference in discrete MRFs.
Alahari, Karteek; Kohli, Pushmeet; Torr, Philip H S
2010-10-01
In this paper, we present novel techniques that improve the computational and memory efficiency of algorithms for solving multilabel energy functions arising from discrete mrfs or crfs. These methods are motivated by the observations that the performance of minimization algorithms depends on: 1) the initialization used for the primal and dual variables and 2) the number of primal variables involved in the energy function. Our first method (dynamic alpha-expansion) works by "recycling" results from previous problem instances. The second method simplifies the energy function by "reducing" the number of unknown variables present in the problem. Further, we show that it can also be used to generate a good initialization for the dynamic alpha-expansion algorithm by "reusing" dual variables. We test the performance of our methods on energy functions encountered in the problems of stereo matching and color and object-based segmentation. Experimental results show that our methods achieve a substantial improvement in the performance of alpha-expansion, as well as other popular algorithms such as sequential tree-reweighted message passing and max-product belief propagation. We also demonstrate the applicability of our schemes for certain higher order energy functions, such as the one described in [1], for interactive texture-based image and video segmentation. In most cases, we achieve a 10-15 times speed-up in the computation time. Our modified alpha-expansion algorithm provides similar performance to Fast-PD, but is conceptually much simpler. Both alpha-expansion and Fast-PD can be made orders of magnitude faster when used in conjunction with the "reduce" scheme proposed in this paper. PMID:20724761
Signatures of discrete breathers in coherent state quantum dynamics
Igumenshchev, Kirill; Ovchinnikov, Misha; Prezhdo, Oleg; Maniadis, Panagiotis
2013-02-07
In classical mechanics, discrete breathers (DBs) - a spatial time-periodic localization of energy - are predicted in a large variety of nonlinear systems. Motivated by a conceptual bridging of the DB phenomena in classical and quantum mechanical representations, we study their signatures in the dynamics of a quantum equivalent of a classical mechanical point in phase space - a coherent state. In contrast to the classical point that exhibits either delocalized or localized motion, the coherent state shows signatures of both localized and delocalized behavior. The transition from normal to local modes have different characteristics in quantum and classical perspectives. Here, we get an insight into the connection between classical and quantum perspectives by analyzing the decomposition of the coherent state into system's eigenstates, and analyzing the spacial distribution of the wave-function density within these eigenstates. We find that the delocalized and localized eigenvalue components of the coherent state are separated by a mixed region, where both kinds of behavior can be observed. Further analysis leads to the following observations. Considered as a function of coupling, energy eigenstates go through avoided crossings between tunneling and non-tunneling modes. The dominance of tunneling modes in the high nonlinearity region is compromised by the appearance of new types of modes - high order tunneling modes - that are similar to the tunneling modes but have attributes of non-tunneling modes. Certain types of excitations preferentially excite higher order tunneling modes, allowing one to study their properties. Since auto-correlation functions decrease quickly in highly nonlinear systems, short-time dynamics are sufficient for modeling quantum DBs. This work provides a foundation for implementing modern semi-classical methods to model quantum DBs, bridging classical and quantum mechanical signatures of DBs, and understanding spectroscopic experiments that
Web-based dynamic Delphi: a new survey instrument
NASA Astrophysics Data System (ADS)
Yao, JingTao; Liu, Wei-Ning
2006-04-01
We present a mathematical model for a dynamic Delphi survey method which takes advantages of Web technology. A comparative study on the performance of the conventional Delphi method and the dynamic Delphi instrument is conducted. It is suggested that a dynamic Delphi survey may form a consensus quickly. However, the result may not be robust due to the judgement leaking issues.
DIVE in the cosmic web: voids with Delaunay triangulation from discrete matter tracer distributions
NASA Astrophysics Data System (ADS)
Zhao, Cheng; Tao, Charling; Liang, Yu; Kitaura, Francisco-Shu; Chuang, Chia-Hsun
2016-07-01
We present a novel parameter-free cosmological void finder (DIVE, Delaunay TrIangulation Void findEr) based on Delaunay Triangulation (DT), which efficiently computes the empty spheres constrained by a discrete set of tracers. We define the spheres as DT voids, and describe their properties, including a universal density profile together with an intrinsic scatter. We apply this technique on 100 halo catalogues with volumes of 2.5 h-1Gpc side each, with a bias and number density similar to the Baryon Oscillation Spectroscopic Survey CMASS luminous red galaxies, performed with the PATCHY code. Our results show that there are two main species of DT voids, which can be characterized by the radius: they have different responses to halo redshift space distortions, to number density of tracers, and reside in different dark matter environments. Based on dynamical arguments using the tidal field tensor, we demonstrate that large DT voids are hosted in expanding regions, whereas the haloes used to construct them reside in collapsing ones. Our approach is therefore able to efficiently determine the troughs of the density field from galaxy surveys, and can be used to study their clustering. We further study the power spectra of DT voids, and find that the bias of the two populations are different, demonstrating that the small DT voids are essentially tracers of groups of haloes.
Enabling a GeoWeb with Net-Centric Fusion on a Discrete Global Grid (Invited)
NASA Astrophysics Data System (ADS)
Peterson, P. R.
2010-12-01
There is a pressing expectation for general real-time access to multi-source geo-spatial content in support of evidence-based decision-making. Earth location promises to be a decentralized organizational structure for such net-centric decision support systems - the GeoWeb, Digital Earth, GEOINT2, and Planetary Skin are some terms in use. However, these platforms assume a critical provision for access and integration of multi-sources of geo-data on-demand and unassisted by the unanticipated unsophisticated end-use decision-maker. How can this occur when geo-data integration is a complex time-consuming problem? We present a solution. A discrete global grid system (DGGS) incorporates an Earth partitioning that acts as a unifying structure for encoding and integrating/fusing multi-source location-based information necessary for this class of location-based platforms. As a global reference model the DGGS is uniform over the entire planet at any resolution - from continents to bird baths. The DGGS provides fast, seamless assimilation of new, numerous, and disparate geo-data sources regardless of scale, origin, resolution, legacy formats, datum, or projection - allowing any content to reside at its own level of granularity at any location on the globe. The DGGS renders data fused, ubiquitous, searchable, and ready for analysis. Testbed development of a DGGS using the optimized Icosahedral Snyder Equal Area aperture 3 Hexagonal grid (ISEA3H) demonstrated solutions to challenging aspects of multi-source data exploitation and decision support within military geospatial-intelligence. The ISEA3H tessellation is optimized to use the fine increments and the close packed equal area partitioning properties of a square root three (hexagonal) subdivision. The investigations advanced the ISEA3H grid development to include cell indexing, quantization strategy and numeric functions required for a formal digital Earth reference model (DERM). Notably, the global index that was selected
Technology Transfer Automated Retrieval System (TEKTRAN)
Composed of linear difference equations, a discrete dynamic system model was designed to reconstruct transcriptional regulations in gene regulatory networks in response to 5-hydroxymethylfurfural, a bioethanol conversion inhibitor for ethanologenic yeast Saccharomyces cerevisiae. The modeling aims ...
Discrete Dynamical Systems Meet the Classic Monkey-and-the-Bananas Problem.
ERIC Educational Resources Information Center
Gannon, Gerald E.; Martelli, Mario U.
2001-01-01
Presents a solution of the three-sailors-and-the-bananas problem and attempts a generalization. Introduces an interesting way of looking at the mathematics with an idea drawn from discrete dynamical systems. (KHR)
A Model of Longitudinal Tension Dynamics in Paper Webs and Implications for Web Breaks
NASA Astrophysics Data System (ADS)
Hristopulos, Dionissios; Uesaka, Tetsu
2002-03-01
Web breaks are rare but costly events with significant impact on printing presses and papermaking machines. Macroscopic paper defects have been traditionally blamed for causing web breaks, but improvements in paper quality have created interest in the microscopic causes of strength reliability of paper and the role of web tension variations. We present a continuum model of longitudinal tension dynamics in paper webs. The model predicts high straining rates at nip contacts, which may have significant impact on strength reliability. From the general model we obtain the low-frequency web response in terms of an ordinary differential equation. We show that the web is a low-pass filter that reduces the impact of dynamic speed perturbations on strain variations. We prove that the empirical relation between mechanical draw and strain increments is valid only in the steady state. We also estimate the strain variations caused by shape deformations of paper rolls, and we derive an explicit expression for the resulting frequency of web breaks based on a weak-link model of fracture. Finally, we briefly discuss the high-frequency response.
Development and evaluation of a dynamic web-based application.
Hsieh, Yichuan; Brennan, Patricia Flatley
2007-01-01
Traditional consumer health informatics (CHI) applications that were developed for lay public on the Web were commonly written in a Hypertext Markup Language (HTML). As genetics knowledge rapidly advances and requires updating information in a timely fashion, a different content structure is therefore needed to facilitate information delivery. This poster will present the process of developing a dynamic database-driven Web CHI application. PMID:18694081
Dynamical behavior of fractional-order Hastings-Powell food chain model and its discretization
NASA Astrophysics Data System (ADS)
Matouk, A. E.; Elsadany, A. A.; Ahmed, E.; Agiza, H. N.
2015-10-01
In this work, the dynamical behavior of fractional-order Hastings-Powell food chain model is investigated and a new discretization method of the fractional-order system is introduced. A sufficient condition for existence and uniqueness of the solution of the proposed system is obtained. Local stability of the equilibrium points of the fractional-order system is studied. Furthermore, the necessary and sufficient conditions of stability of the discretized system are also studied. It is shown that the system's fractional parameter has effect on the stability of the discretized system which shows rich variety of dynamical behaviors such as Hopf bifurcation, an attractor crisis and chaotic attractors. Numerical simulations show the tea-cup chaotic attractor of the fractional-order system and the richer dynamical behavior of the corresponding discretized system.
Discrete neural dynamic programming in wheeled mobile robot control
NASA Astrophysics Data System (ADS)
Hendzel, Zenon; Szuster, Marcin
2011-05-01
In this paper we propose a discrete algorithm for a tracking control of a two-wheeled mobile robot (WMR), using an advanced Adaptive Critic Design (ACD). We used Dual-Heuristic Programming (DHP) algorithm, that consists of two parametric structures implemented as Neural Networks (NNs): an actor and a critic, both realized in a form of Random Vector Functional Link (RVFL) NNs. In the proposed algorithm the control system consists of the DHP adaptive critic, a PD controller and a supervisory term, derived from the Lyapunov stability theorem. The supervisory term guaranties a stable realization of a tracking movement in a learning phase of the adaptive critic structure and robustness in face of disturbances. The discrete tracking control algorithm works online, uses the WMR model for a state prediction and does not require a preliminary learning. Verification has been conducted to illustrate the performance of the proposed control algorithm, by a series of experiments on the WMR Pioneer 2-DX.
On the Dynamics of Some Discretizations of Convection-Diffusion Equations
NASA Technical Reports Server (NTRS)
Sweby, Peter K.; Yee, H. C.; Rai, Man Mohan (Technical Monitor)
1995-01-01
Numerical discretizations of differential equations which model physical processes can possess dynamics quite different from that of the equations themselves. Recently the emphasis has been on the the dynamics of numerical discretizations for Ordinary Differential Equations (ODEs). For Partial Differential Equations (PDEs) using a method of lines approach the situation is more complex. First, the spatial discretisation may introduce dynamics not present in the original equations; second, the solution of the resulting system of ODEs is open to the modified dynamics of the ODE solver used. These two effects may interact in a complex manner. In this talk we present some results of our recent work on the dynamics of discretizations of convection-diffusion equations, including those produced using Total Variation Diminishing (TVD) schemes and adaptive grid techniques. A more general overview of the area may be found on our accompanying poster presentation.
Astronomy Education via The Dynamic Web
NASA Astrophysics Data System (ADS)
Flurchick, K. M.; Avery, W.; Griego, B. F.; Culver, R.
2011-05-01
The ability of web applications to provide students the ability to explore and investigate astronomical concepts presented in class in a way which can help student understanding. In this presentation we report on the results of students making use of the computational tools in webMathematicaTM to analyze and investigate a variety of astronomical phenomena, including topics such as the Runge-Lenz vector, descriptions of the orbits of the exo-planets and other topics related to celestial mechanics. Using the exercise described herein, students at the North Carolina A & T State University and Colorado State University investigated via computational simulations the creation and characteristics and the effects of various parameters on these systems being studied.
Dislocation dynamics of web type silicon ribbon
NASA Technical Reports Server (NTRS)
Dillon, O. W., Jr.; Tsai, C. T.; De Angelis, R. J.
1987-01-01
Silicon ribbon grown by the dendritic web process passes through a rapidly changing thermal profile in the growth direction. This rapidly changing profile induces stresses which produce changes in the dislocation density in the ribbon. A viscoplastic material response function (Haasen-Sumino model) is used herein to calculate the stresses and the dislocation density at each point in the silicon ribbon. The residual stresses are also calculated.
Simulating food web dynamics along a gradient: quantifying human influence.
Jordán, Ferenc; Gjata, Nerta; Mei, Shu; Yule, Catherine M
2012-01-01
Realistically parameterized and dynamically simulated food-webs are useful tool to explore the importance of the functional diversity of ecosystems, and in particular relations between the dynamics of species and the whole community. We present a stochastic dynamical food web simulation for the Kelian River (Borneo). The food web was constructed for six different locations, arrayed along a gradient of increasing human perturbation (mostly resulting from gold mining activities) along the river. Along the river, the relative importance of grazers, filterers and shredders decreases with increasing disturbance downstream, while predators become more dominant in governing eco-dynamics. Human activity led to increased turbidity and sedimentation which adversely impacts primary productivity. Since the main difference between the study sites was not the composition of the food webs (structure is quite similar) but the strengths of interactions and the abundance of the trophic groups, a dynamical simulation approach seemed to be useful to better explain human influence. In the pristine river (study site 1), when comparing a structural version of our model with the dynamical model we found that structurally central groups such as omnivores and carnivores were not the most important ones dynamically. Instead, primary consumers such as invertebrate grazers and shredders generated a greater dynamical response. Based on the dynamically most important groups, bottom-up control is replaced by the predominant top-down control regime as distance downstream and human disturbance increased. An important finding, potentially explaining the poor structure to dynamics relationship, is that indirect effects are at least as important as direct ones during the simulations. We suggest that our approach and this simulation framework could serve systems-based conservation efforts. Quantitative indicators on the relative importance of trophic groups and the mechanistic modeling of eco-dynamics
Hamiltonian dynamics for complex food webs.
Kozlov, Vladimir; Vakulenko, Sergey; Wennergren, Uno
2016-03-01
We investigate stability and dynamics of large ecological networks by introducing classical methods of dynamical system theory from physics, including Hamiltonian and averaging methods. Our analysis exploits the topological structure of the network, namely the existence of strongly connected nodes (hubs) in the networks. We reveal new relations between topology, interaction structure, and network dynamics. We describe mechanisms of catastrophic phenomena leading to sharp changes of dynamics and hence completely altering the ecosystem. We also show how these phenomena depend on the structure of interaction between species. We can conclude that a Hamiltonian structure of biological interactions leads to stability and large biodiversity. PMID:27078396
Hamiltonian dynamics for complex food webs
NASA Astrophysics Data System (ADS)
Kozlov, Vladimir; Vakulenko, Sergey; Wennergren, Uno
2016-03-01
We investigate stability and dynamics of large ecological networks by introducing classical methods of dynamical system theory from physics, including Hamiltonian and averaging methods. Our analysis exploits the topological structure of the network, namely the existence of strongly connected nodes (hubs) in the networks. We reveal new relations between topology, interaction structure, and network dynamics. We describe mechanisms of catastrophic phenomena leading to sharp changes of dynamics and hence completely altering the ecosystem. We also show how these phenomena depend on the structure of interaction between species. We can conclude that a Hamiltonian structure of biological interactions leads to stability and large biodiversity.
NASA Astrophysics Data System (ADS)
Wang, Tong; Ding, Yongsheng; Zhang, Lei; Hao, Kuangrong
2016-08-01
This paper considered the synchronisation of continuous complex dynamical networks with discrete-time communications and delayed nodes. The nodes in the dynamical networks act in the continuous manner, while the communications between nodes are discrete-time; that is, they communicate with others only at discrete time instants. The communication intervals in communication period can be uncertain and variable. By using a piecewise Lyapunov-Krasovskii function to govern the characteristics of the discrete communication instants, we investigate the adaptive feedback synchronisation and a criterion is derived to guarantee the existence of the desired controllers. The globally exponential synchronisation can be achieved by the controllers under the updating laws. Finally, two numerical examples including globally coupled network and nearest-neighbour coupled networks are presented to demonstrate the validity and effectiveness of the proposed control scheme.
NASA Astrophysics Data System (ADS)
Zañudo, Jorge G. T.; Albert, Réka
2013-06-01
Discrete dynamic models are a powerful tool for the understanding and modeling of large biological networks. Although a lot of progress has been made in developing analysis tools for these models, there is still a need to find approaches that can directly relate the network structure to its dynamics. Of special interest is identifying the stable patterns of activity, i.e., the attractors of the system. This is a problem for large networks, because the state space of the system increases exponentially with network size. In this work, we present a novel network reduction approach that is based on finding network motifs that stabilize in a fixed state. Notably, we use a topological criterion to identify these motifs. Specifically, we find certain types of strongly connected components in a suitably expanded representation of the network. To test our method, we apply it to a dynamic network model for a type of cytotoxic T cell cancer and to an ensemble of random Boolean networks of size up to 200. Our results show that our method goes beyond reducing the network and in most cases can actually predict the dynamical repertoire of the nodes (fixed states or oscillations) in the attractors of the system.
Food-web dynamics in a large river discontinuum
Cross, Wyatt F.; Baxter, Colden V.; Rosi-Marshall, Emma J.; Hall, Robert O., Jr.; Kennedy, Theodore A.; Donner, Kevin C.; Kelly, Holly A. Wellard; Seegert, Sarah E.Z.; Behn, Kathrine E.; Yard, Michael D.
2013-01-01
Nearly all ecosystems have been altered by human activities, and most communities are now composed of interacting species that have not co-evolved. These changes may modify species interactions, energy and material flows, and food-web stability. Although structural changes to ecosystems have been widely reported, few studies have linked such changes to dynamic food-web attributes and patterns of energy flow. Moreover, there have been few tests of food-web stability theory in highly disturbed and intensely managed freshwater ecosystems. Such synthetic approaches are needed for predicting the future trajectory of ecosystems, including how they may respond to natural or anthropogenic perturbations. We constructed flow food webs at six locations along a 386-km segment of the Colorado River in Grand Canyon (Arizona, USA) for three years. We characterized food-web structure and production, trophic basis of production, energy efficiencies, and interaction-strength distributions across a spatial gradient of perturbation (i.e., distance from Glen Canyon Dam), as well as before and after an experimental flood. We found strong longitudinal patterns in food-web characteristics that strongly correlated with the spatial position of large tributaries. Above tributaries, food webs were dominated by nonnative New Zealand mudsnails (62% of production) and nonnative rainbow trout (100% of fish production). The simple structure of these food webs led to few dominant energy pathways (diatoms to few invertebrate taxa to rainbow trout), large energy inefficiencies (i.e., Below large tributaries, invertebrate production declined ∼18-fold, while fish production remained similar to upstream sites and comprised predominately native taxa (80–100% of production). Sites below large tributaries had increasingly reticulate and detritus-based food webs with a higher prevalence of omnivory, as well as interaction strength distributions more typical of theoretically stable food webs (i
Dynamic Reasoning under Probabilistic Uncertainty in the Semantic Web
NASA Astrophysics Data System (ADS)
Chen, Limin; Shi, Zhongzhi
Towards a full realization of the Semantic Web as a source of processable information and services, we extend description logic ALCO@ to admit dynamic reasoning under probabilistic uncertainty. In our formalism, the crisp knowledge is encoded in DL-knowledgebase, and the uncertain one is in conditional constraints; atomic actions are represented in terms of their preconditions (conditional constraints) and effects (possibly negated ALCO@-assertions involving only atomic concepts). We also summarize some reasoning tasks in the new formalism and provide algorithms for them. This work is strictly related to the practical need of DL-based reasoning in the Semantic Web and is much closer than the prior formalisms to a practical formalism in the Web setting full of uncertain information and dynamic processing.
NASA Astrophysics Data System (ADS)
Lehtinen, Arttu; Granberg, Fredric; Laurson, Lasse; Nordlund, Kai; Alava, Mikko J.
2016-01-01
The stress-driven motion of dislocations in crystalline solids, and thus the ensuing plastic deformation process, is greatly influenced by the presence or absence of various pointlike defects such as precipitates or solute atoms. These defects act as obstacles for dislocation motion and hence affect the mechanical properties of the material. Here we combine molecular dynamics studies with three-dimensional discrete dislocation dynamics simulations in order to model the interaction between different kinds of precipitates and a 1/2 <111 > {110 } edge dislocation in BCC iron. We have implemented immobile spherical precipitates into the ParaDis discrete dislocation dynamics code, with the dislocations interacting with the precipitates via a Gaussian potential, generating a normal force acting on the dislocation segments. The parameters used in the discrete dislocation dynamics simulations for the precipitate potential, the dislocation mobility, shear modulus, and dislocation core energy are obtained from molecular dynamics simulations. We compare the critical stresses needed to unpin the dislocation from the precipitate in molecular dynamics and discrete dislocation dynamics simulations in order to fit the two methods together and discuss the variety of the relevant pinning and depinning mechanisms.
Lehtinen, Arttu; Granberg, Fredric; Laurson, Lasse; Nordlund, Kai; Alava, Mikko J
2016-01-01
The stress-driven motion of dislocations in crystalline solids, and thus the ensuing plastic deformation process, is greatly influenced by the presence or absence of various pointlike defects such as precipitates or solute atoms. These defects act as obstacles for dislocation motion and hence affect the mechanical properties of the material. Here we combine molecular dynamics studies with three-dimensional discrete dislocation dynamics simulations in order to model the interaction between different kinds of precipitates and a 1/2〈111〉{110} edge dislocation in BCC iron. We have implemented immobile spherical precipitates into the ParaDis discrete dislocation dynamics code, with the dislocations interacting with the precipitates via a Gaussian potential, generating a normal force acting on the dislocation segments. The parameters used in the discrete dislocation dynamics simulations for the precipitate potential, the dislocation mobility, shear modulus, and dislocation core energy are obtained from molecular dynamics simulations. We compare the critical stresses needed to unpin the dislocation from the precipitate in molecular dynamics and discrete dislocation dynamics simulations in order to fit the two methods together and discuss the variety of the relevant pinning and depinning mechanisms. PMID:26871192
Navier-Stokes Dynamics by a Discrete Boltzmann Model
NASA Technical Reports Server (NTRS)
Rubinstein, Robet
2010-01-01
This work investigates the possibility of particle-based algorithms for the Navier-Stokes equations and higher order continuum approximations of the Boltzmann equation; such algorithms would generalize the well-known Pullin scheme for the Euler equations. One such method is proposed in the context of a discrete velocity model of the Boltzmann equation. Preliminary results on shock structure are consistent with the expectation that the shock should be much broader than the near discontinuity predicted by the Pullin scheme, yet narrower than the prediction of the Boltzmann equation. We discuss the extension of this essentially deterministic method to a stochastic particle method that, like DSMC, samples the distribution function rather than resolving it completely.
Discrete-time pilot model. [human dynamics and digital simulation
NASA Technical Reports Server (NTRS)
Cavalli, D.
1978-01-01
Pilot behavior is considered as a discrete-time process where the decision making has a sequential nature. This model differs from both the quasilinear model which follows from classical control theory and from the optimal control model which considers the human operator as a Kalman estimator-predictor. An additional factor considered is that the pilot's objective may not be adequately formulated as a quadratic cost functional to be minimized, but rather as a more fuzzy measure of the closeness with which the aircraft follows a reference trajectory. All model parameters, in the digital program simulating the pilot's behavior, were successfully compared in terms of standard-deviation and performance with those of professional pilots in IFR configuration. The first practical application of the model was in the study of its performance degradation when the aircraft model static margin decreases.
Discrete Element Method Simulations of Ice Floe Dynamics
NASA Astrophysics Data System (ADS)
Calantoni, J.; Bateman, S. P.; Shi, F.; Orzech, M.; Veeramony, J.
2014-12-01
Ice floes were modeled using LIGGGHTS, an open source discrete element method (DEM) software, where individual elements were bonded together to make floes. The bonds were allowed to break with a critical stress calibrated to existing laboratory measurements for the compressive, tensile, and flexural strength of ice floes. The DEM allows for heterogeneous shape and size distributions of the ice floes to evolve over time. We simulated the interaction between sea ice and ocean waves in the marginal ice zone using a coupled wave-ice system. The waves were modeled with NHWAVE, a non-hydrostatic wave model that predicts instantaneous surface elevation and the three-dimensional flow field. The ice floes and waves were coupled through buoyancy and drag forces. Preliminary comparisons with field and laboratory measurements for coupled simulations will be presented.
Extinction dynamics of a discrete population in an oasis
NASA Astrophysics Data System (ADS)
Berti, Stefano; Cencini, Massimo; Vergni, Davide; Vulpiani, Angelo
2015-07-01
Understanding the conditions ensuring the persistence of a population is an issue of primary importance in population biology. The first theoretical approach to the problem dates back to the 1950s with the Kierstead, Slobodkin, and Skellam (KiSS) model, namely a continuous reaction-diffusion equation for a population growing on a patch of finite size L surrounded by a deadly environment with infinite mortality, i.e., an oasis in a desert. The main outcome of the model is that only patches above a critical size allow for population persistence. Here we introduce an individual-based analog of the KiSS model to investigate the effects of discreteness and demographic stochasticity. In particular, we study the average time to extinction both above and below the critical patch size of the continuous model and investigate the quasistationary distribution of the number of individuals for patch sizes above the critical threshold.
Studying Human Dynamics Through Web Analytics
NASA Astrophysics Data System (ADS)
Ramasco, Jose; Goncalves, Bruno
2008-03-01
When Tim Berners Lee, a physicist at the European Center for Nuclear Research (CERN) first conceived the World Wide Web (WWW) in 1990 as a way to facilitate the sharing of scientific information and results among the centers different researchers and groups, even the most ingenious of science fiction writers could not have imagined the role it would come to play in the following decades. The increasing ubiquitousness of Internet access and the frequency with which people interact with it raise the possibility of using it to better observe, understand, and even monitor several aspects of human social behavior. Websites with large numbers of frequently returning users, such as search engines, company or university websites, are ideal for this task. The properly anonymized logs detailing the access history to Emory University's website is studied. We find that a small number of users is responsible for a finite fraction of the total activity. A saturation phenomenon is observed where, certain connections age, becoming less attractive to new activity over time. Finally, by measuring the average activity as a function of the day of the week, we find that productivity seems to be higher on Tuesdays and Wednesdays, with Sundays being the least active day.
Discrete and continuous dynamics modeling of a mass moving on a flexible structure
NASA Technical Reports Server (NTRS)
Herman, Deborah Ann
1992-01-01
A general discrete methodology for modeling the dynamics of a mass that moves on the surface of a flexible structure is developed. This problem was motivated by the Space Station/Mobile Transporter system. A model reduction approach is developed to make the methodology applicable to large structural systems. To validate the discrete methodology, continuous formulations are also developed. Three different systems are examined: (1) simply-supported beam, (2) free-free beam, and (3) free-free beam with two points of contact between the mass and the flexible beam. In addition to validating the methodology, parametric studies were performed to examine how the system's physical properties affect its dynamics.
Dynamics of Discrete Breathers in a Pt3Al Crystal
NASA Astrophysics Data System (ADS)
Starostenkov, M. D.; Potekaev, A. I.; Dmitriev, S. V.; Zakharov, P. V.; Eremin, A. M.; Kulagina, V. V.
2016-01-01
The discrete breathers in a Pt3Al crystal, which exhibit soft (DB1) and hard (DB2) nonlinearity, are shown to possess a number of principal differences. Unlike an immobile and stable DB1, a DB2 breather is mainly localized on four Al atoms and is stretched along one of the close-packed rows of crystals. On the other hand, DB2 can displace hundreds of nanometers along one of the directions of close packing. Having localized a considerable amount of energy, both DB1 and DB2 breathers slowly emit it during their lifetime. A collision of DB1 and DB2 results in part of their energy being released into the Al sublattice, the larger part lost by DB2 that is destroyed faster than DB1. The DB2 breather can effectively transport the energy throughout the crystal, and a collision of DBs results in its considerable localization in the crystal. A capability of transferring the energy can thus give rise to structural transformations far from the focus of excitation of such localized objects.
Simulating continuous-time Hamiltonian dynamics by way of a discrete-time quantum walk
NASA Astrophysics Data System (ADS)
Schmitz, A. T.; Schwalm, W. A.
2016-03-01
Much effort has been made to connect the continuous-time and discrete-time quantum walks. We present a method for making that connection for a general graph Hamiltonian on a bigraph. Furthermore, such a scheme may be adapted for simulating discretized quantum models on a quantum computer. A coin operator is found for the discrete-time quantum walk which exhibits the same dynamics as the continuous-time evolution. Given the spectral decomposition of the graph Hamiltonian and certain restrictions, the discrete-time evolution is solved for explicitly and understood at or near important values of the parameters. Finally, this scheme is connected to past results for the 1D chain.
Network Analysis of the State Space of Discrete Dynamical Systems
NASA Astrophysics Data System (ADS)
Shreim, Amer; Grassberger, Peter; Nadler, Walter; Samuelsson, Björn; Socolar, Joshua E. S.; Paczuski, Maya
2007-05-01
We study networks representing the dynamics of elementary 1D cellular automata (CA) on finite lattices. We analyze scaling behaviors of both local and global network properties as a function of system size. The scaling of the largest node in-degree is obtained analytically for a variety of CA including rules 22, 54, and 110. We further define the path diversity as a global network measure. The coappearance of nontrivial scaling in both the hub size and the path diversity separates simple dynamics from the more complex behaviors typically found in Wolfram’s class IV and some class III CA.
Network analysis of the state space of discrete dynamical systems.
Shreim, Amer; Grassberger, Peter; Nadler, Walter; Samuelsson, Björn; Socolar, Joshua E S; Paczuski, Maya
2007-05-11
We study networks representing the dynamics of elementary 1D cellular automata (CA) on finite lattices. We analyze scaling behaviors of both local and global network properties as a function of system size. The scaling of the largest node in-degree is obtained analytically for a variety of CA including rules 22, 54, and 110. We further define the path diversity as a global network measure. The coappearance of nontrivial scaling in both the hub size and the path diversity separates simple dynamics from the more complex behaviors typically found in Wolfram's class IV and some class III CA. PMID:17677672
A Voice Web Application Based on Dynamic Navigation of VXML
NASA Astrophysics Data System (ADS)
Bhingarkar, Sukhada P.
2010-11-01
Voice Web, as the name suggests, accesses web resources via voice. VoiceXML is the markup language used to develop speech applications. VoiceXML is interactive and allows voice input to be received and processed by voice browser. Unfortunately, the navigation of VoiceXML document is completely controlled by application developer. Also, the user does not have flexibility to utter random word from currently executing dialog. The aim of the paper is to address the weakness of VoiceXML and develop an application, which dynamically detects recognition candidates in user content, in contrast with recognition candidates of the existing voice web, which depend on the application developer. In this application, a news service is implemented along with dictionary of IT-specific terms and dictionary of words from currently executing news.
Emerging hierarchies in dynamically adapting webs
NASA Astrophysics Data System (ADS)
Katifori, Eleni; Graewer, Johannes; Magnasco, Marcelo; Modes, Carl
Transport networks play a key role across four realms of eukaryotic life: slime molds, fungi, plants, and animals. In addition to the developmental algorithms that build them, many also employ adaptive strategies to respond to stimuli, damage, and other environmental changes. We model these adapting network architectures using a generic dynamical system on weighted graphs and find in simulation that these networks ultimately develop a hierarchical organization of the final weighted architecture accompanied by the formation of a system-spanning backbone. We quantify the hierarchical organization of the networks by developing an algorithm that decomposes the architecture to multiple scales and analyzes how the organization in each scale relates to that of the scale above and below it. The methodologies developed in this work are applicable to a wide range of systems including the slime mold physarum polycephalum, human microvasculature, and force chains in granular media.
Bifurcations and dynamics of a discrete predator–prey system
Asheghi, Rasoul
2014-01-01
In this paper, we study the dynamics behaviour of a stratum of plant–herbivore which is modelled through the following F(x, y)=(f(x, y), g(x, y)) two-dimensional map with four parameters defined by where x≥0, y≥0, and the real parameters a, b, r, k are all positive. We will focus on the case a≠b. We study the stability of fixed points and do the analysis of the period-doubling and the Neimark–Sacker bifurcations in a standard way. PMID:24963984
ERIC Educational Resources Information Center
Winkel, Brian
2012-01-01
We give an example of cross coursing in which a subject or approach in one course in undergraduate mathematics is used in a completely different course. This situation crosses falling body modelling in an upper level differential equations course into a modest discrete dynamical systems unit of a first-year mathematics course. (Contains 1 figure.)
GA-based discrete dynamic programming approach for scheduling in FMS environments.
Yang, J B
2001-01-01
The paper presents a new genetic algorithm (GA)-based discrete dynamic programming (DDP) approach for generating static schedules in a flexible manufacturing system (FMS) environment. This GA-DDP approach adopts a sequence-dependent schedule generation strategy, where a GA is employed to generate feasible job sequences and a series of discrete dynamic programs are constructed to generate legal schedules for a given sequence of jobs. In formulating the GA, different performance criteria could be easily included. The developed DDF algorithm is capable of identifying locally optimized partial schedules and shares the computation efficiency of dynamic programming. The algorithm is designed In such a way that it does not suffer from the state explosion problem inherent in pure dynamic programming approaches in FMS scheduling. Numerical examples are reported to illustrate the approach. PMID:18244848
A discrete dynamical system as a model of a neural network with generalized Hebbian synapses
NASA Astrophysics Data System (ADS)
Pieniazek, Leszek; Wójtowicz, Bernardeta
2005-08-01
We study the dynamical behavior of a discrete time dynamical system which can serve as a model of a learning process. We determine fixed points of this system and basins of attraction of attracting points. This system was studied by Fernanda Botelho and James J. Jamison in [A learning rule with generalized Hebbian synapses, J. Math. Anal. Appl. 273 (2002) 529-547] but authors used its continuous counterpart to describe basins of attraction.
Toccafondo, Iacopo; Taki, Mohammad; Signorini, Alessandro; Zaidi, Farhan; Nannipieri, Tiziano; Faralli, Stefano; Di Pasquale, Fabrizio
2012-11-01
We propose and experimentally demonstrate a hybrid fiber optic sensing technique that effectively combines Raman optical time domain reflectometry and in-line time-division-multiplexing for fiber Bragg grating (FBG) dynamic interrogation. The highly integrated proposed scheme employs broadband apodized low reflectivity FBGs with a single narrowband optical source and a shared receiver block, allowing for simultaneous measurements of distributed static temperature and discrete dynamic strain, over the same sensing fiber. PMID:23114320
Dynamics in a three species food-web system
NASA Astrophysics Data System (ADS)
Gupta, K.; Gakkhar, S.
2016-04-01
In this paper, the dynamics of a three species food-web system is discussed. The food-web comprises of one predator and two logistically growing competing species. The predator species is taking food from one of the competitors with Holling type II functional response. Another competitor is the amensal species for the predator of first species. The system is shown to be positive and bounded. The stability of various axial points, boundary points and interior point has been investigated. The persistence of the system has been studied. Numerical simulation has been performed to show the occurrence of Hopf bifurcation and stable limit cycle about the interior point. The presence of second competitor and its interaction with predator gives more complex dynamics than the simple prey-predator system. The existence of transcritical bifurcation has been established about two axial points. The existence of periodic attractor having period-2 solution has been shown, when amensal coefficient is chosen as bifurcation parameter.
Čech cohomology of attractors of discrete dynamical systems
NASA Astrophysics Data System (ADS)
Ruiz del Portal, Francisco R.; Sánchez-Gabites, J. J.
2014-10-01
Let f:Rn→Rn be a homeomorphism and K an asymptotically stable attractor for f. The aim of this paper is to study when the inclusion of K in its basin of attraction A(K) induces isomorphisms in Čech cohomology. We show that (i) this is true if coefficients are taken in Q or Zp (p prime) and (ii) it is true for integral cohomology if and only if the Čech cohomology of K or A(K) is finitely generated. We compute the Čech cohomology of periodic point free attractors of volume-contracting R3-homeomorphisms and present applications to quite general models in population dynamics.
Dynamics of quantal heating in electron systems with discrete spectra
NASA Astrophysics Data System (ADS)
Dietrich, Scott; Mayer, William; Vitkalov, Sergey; Bykov, A. A.
2015-05-01
The temporal evolution of quantal Joule heating of two-dimensional (2D) electrons in a GaAs quantum well placed in quantizing magnetic fields is studied using a difference-frequency method. The method is based on measurements of the electron conductivity oscillating at the beat frequency f =f1-f2 between two microwaves applied to the 2D system at frequencies f1 and f2. The method provides direct access to the dynamical characteristics of the heating and yields the inelastic-scattering time τi n of 2D electrons. The obtained τi n is strongly temperature dependent, varying from 0.13 ns at 5.5 K to 1 ns at 2.4 K in magnetic field B =0.333 T . When the temperature T exceeds the Landau-level separation, the relaxation rate 1 /τi n is proportional to T2, indicating electron-electron interaction as the dominant mechanism limiting the quantal heating. At lower temperatures, the rate tends to be proportional to T3, indicating considerable contribution from electron-phonon scattering.
Discrete Adjoint-Based Design Optimization of Unsteady Turbulent Flows on Dynamic Unstructured Grids
NASA Technical Reports Server (NTRS)
Nielsen, Eric J.; Diskin, Boris; Yamaleev, Nail K.
2009-01-01
An adjoint-based methodology for design optimization of unsteady turbulent flows on dynamic unstructured grids is described. The implementation relies on an existing unsteady three-dimensional unstructured grid solver capable of dynamic mesh simulations and discrete adjoint capabilities previously developed for steady flows. The discrete equations for the primal and adjoint systems are presented for the backward-difference family of time-integration schemes on both static and dynamic grids. The consistency of sensitivity derivatives is established via comparisons with complex-variable computations. The current work is believed to be the first verified implementation of an adjoint-based optimization methodology for the true time-dependent formulation of the Navier-Stokes equations in a practical computational code. Large-scale shape optimizations are demonstrated for turbulent flows over a tiltrotor geometry and a simulated aeroelastic motion of a fighter jet.
Mathematics of Web science: structure, dynamics and incentives.
Chayes, Jennifer
2013-03-28
Dr Chayes' talk described how, to a discrete mathematician, 'all the world's a graph, and all the people and domains merely vertices'. A graph is represented as a set of vertices V and a set of edges E, so that, for instance, in the World Wide Web, V is the set of pages and E the directed hyperlinks; in a social network, V is the people and E the set of relationships; and in the autonomous system Internet, V is the set of autonomous systems (such as AOL, Yahoo! and MSN) and E the set of connections. This means that mathematics can be used to study the Web (and other large graphs in the online world) in the following way: first, we can model online networks as large finite graphs; second, we can sample pieces of these graphs; third, we can understand and then control processes on these graphs; and fourth, we can develop algorithms for these graphs and apply them to improve the online experience. PMID:23419846
A discrete geometric approach for simulating the dynamics of thin viscous threads
Audoly, B.; Clauvelin, N.; Brun, P.-T.; Bergou, M.; Grinspun, E.; Wardetzky, M.
2013-11-15
We present a numerical model for the dynamics of thin viscous threads based on a discrete, Lagrangian formulation of the smooth equations. The model makes use of a condensed set of coordinates, called the centerline/spin representation: the kinematic constraints linking the centerline's tangent to the orientation of the material frame is used to eliminate two out of three degrees of freedom associated with rotations. Based on a description of twist inspired from discrete differential geometry and from variational principles, we build a full-fledged discrete viscous thread model, which includes in particular a discrete representation of the internal viscous stress. Consistency of the discrete model with the classical, smooth equations for thin threads is established formally. Our numerical method is validated against reference solutions for steady coiling. The method makes it possible to simulate the unsteady behavior of thin viscous threads in a robust and efficient way, including the combined effects of inertia, stretching, bending, twisting, large rotations and surface tension.
NASA Astrophysics Data System (ADS)
Alicandro, Roberto; De Luca, Lucia; Garroni, Adriana; Ponsiglione, Marcello
2014-10-01
This paper aims at building a variational approach to the dynamics of discrete topological singularities in two dimensions, based on Γ-convergence. We consider discrete systems, described by scalar functions defined on a square lattice and governed by periodic interaction potentials. Our main motivation comes from XY spin systems, described by the phase parameter, and screw dislocations, described by the displacement function. For these systems, we introduce a discrete notion of vorticity. As the lattice spacing tends to zero we derive the first order Γ-limit of the free energy which is referred to as renormalized energy and describes the interaction of vortices. As a byproduct of this analysis, we show that such systems exhibit increasingly many metastable configurations of singularities. Therefore, we propose a variational approach to the depinning and dynamics of discrete vortices, based on minimizing movements. We show that, letting first the lattice spacing and then the time step of the minimizing movements tend to zero, the vortices move according with the gradient flow of the renormalized energy, as in the continuous Ginzburg-Landau framework.
Minimax terminal approach problem in two-level hierarchical nonlinear discrete-time dynamical system
NASA Astrophysics Data System (ADS)
Shorikov, A. F.
2015-11-01
We consider a discrete-time dynamical system consisting of three controllable objects. The motions of all objects are given by the corresponding vector nonlinear or linear discrete-time recurrent vector relations, and control system for its has two levels: basic (first or I level) that is dominating and subordinate level (second or II level) and both have different criterions of functioning and united a priori by determined informational and control connections defined in advance. For the dynamical system in question, we propose a mathematical formalization in the form of solving a multistep problem of two-level hierarchical minimax program control over the terminal approach process with incomplete information and give a general scheme for its solving.
Distributing flight dynamics products via the World Wide Web
NASA Technical Reports Server (NTRS)
Woodard, Mark; Matusow, David
1996-01-01
The NASA Flight Dynamics Products Center (FDPC), which make available selected operations products via the World Wide Web, is reported on. The FDPC can be accessed from any host machine connected to the Internet. It is a multi-mission service which provides Internet users with unrestricted access to the following standard products: antenna contact predictions; ground tracks; orbit ephemerides; mean and osculating orbital elements; earth sensor sun and moon interference predictions; space flight tracking data network summaries; and Shuttle transport system predictions. Several scientific data bases are available through the service.
Stochastic aspects of one-dimensional discrete dynamical systems: Benford's law
NASA Astrophysics Data System (ADS)
Snyder, Mark A.; Curry, James H.; Dougherty, Anne M.
2001-08-01
Benford's law owes its discovery to the ``Grubby Pages Hypothesis,'' a 19th century observation made by Simon Newcomb that the beginning pages of logarithm books were grubbier than the last few pages, implying that scientists referenced the values toward the front of the books more frequently. If a data set satisfies Benford's law, then it's significant digits will have a logarithmic distribution, which favors smaller significant digits. In this article we demonstrate two ways of creating discrete one-dimensional dynamical systems that satisfy Benford's law. We also develop a numerical simulation methodology that we use to study dynamical systems when analytical results are not readily available.
Stochastic aspects of one-dimensional discrete dynamical systems: Benford's law.
Snyder, M A; Curry, J H; Dougherty, A M
2001-08-01
Benford's law owes its discovery to the "Grubby Pages Hypothesis," a 19th century observation made by Simon Newcomb that the beginning pages of logarithm books were grubbier than the last few pages, implying that scientists referenced the values toward the front of the books more frequently. If a data set satisfies Benford's law, then it's significant digits will have a logarithmic distribution, which favors smaller significant digits. In this article we demonstrate two ways of creating discrete one-dimensional dynamical systems that satisfy Benford's law. We also develop a numerical simulation methodology that we use to study dynamical systems when analytical results are not readily available. PMID:11497692
Eisenhauer, Philipp; Heckman, James J.; Mosso, Stefano
2015-01-01
We compare the performance of maximum likelihood (ML) and simulated method of moments (SMM) estimation for dynamic discrete choice models. We construct and estimate a simplified dynamic structural model of education that captures some basic features of educational choices in the United States in the 1980s and early 1990s. We use estimates from our model to simulate a synthetic dataset and assess the ability of ML and SMM to recover the model parameters on this sample. We investigate the performance of alternative tuning parameters for SMM. PMID:26494926
A microcosm model with discrete time: an approach and dynamical regimes
NASA Astrophysics Data System (ADS)
Evdokimov, E. V.; Shapovalov, A. V.
A dynamical model of a closed 4-chain microcosm (with non-organic component, producers, consumers, and detritus) with discrete time is proposed. The model is based on a simple 3-chain integrable system of biological growth - eating away (non-organic component, producers, consumers) with continuous time. To construct the 4-chain model we integrate the 3-chain system and find the population numbers as explicit functions of time. By application of the time-shift operation for these functions we construct a map which describes the dynamics of the system with discrete time. This map is modified by introduction of additional terms responsible for producer and consumers death-rate, intra-species competition, and a relation describing detritus dynamics. Computer simulation of the map dynamics is performed. The model is shown to have stable regimes for correspondent values of the model parameters and initial data. The following types of dynamical regimes are also revealed: stable stationary states, decaying oscillations, limit cycles, and quasiharmonic oscillations. The detritus chain plays the role of a "friction" factor of the system oscillations. The work is supported by Russian Foundation for the Humanities.
An Isometric Dynamics for a Causal Set Approach to Discrete Quantum Gravity
NASA Astrophysics Data System (ADS)
Gudder, S.
2015-12-01
We consider a covariant causal set approach to discrete quantum gravity. We first review the microscopic picture of this approach. In this picture a universe grows one element at a time and its geometry is determined by a sequence of integers called the shell sequence. We next present the macroscopic picture which is described by a sequential growth process. We introduce a model in which the dynamics is governed by a quantum transition amplitude. The amplitude satisfies a stochastic and unitary condition and the resulting dynamics becomes isometric. We show that the dynamics preserves stochastic states. By "doubling down" on the dynamics we obtain a unitary group representation and a natural energy operator. These unitary operators are employed to define canonical position and momentum operators.
A network of discrete events for the representation and analysis of diffusion dynamics
NASA Astrophysics Data System (ADS)
Pintus, Alberto M.; Pazzona, Federico G.; Demontis, Pierfranco; Suffritti, Giuseppe B.
2015-11-01
We developed a coarse-grained description of the phenomenology of diffusive processes, in terms of a space of discrete events and its representation as a network. Once a proper classification of the discrete events underlying the diffusive process is carried out, their transition matrix is calculated on the basis of molecular dynamics data. This matrix can be represented as a directed, weighted network where nodes represent discrete events, and the weight of edges is given by the probability that one follows the other. The structure of this network reflects dynamical properties of the process of interest in such features as its modularity and the entropy rate of nodes. As an example of the applicability of this conceptual framework, we discuss here the physics of diffusion of small non-polar molecules in a microporous material, in terms of the structure of the corresponding network of events, and explain on this basis the diffusivity trends observed. A quantitative account of these trends is obtained by considering the contribution of the various events to the displacement autocorrelation function.
Lin, Shih-Wei; Ying, Kuo-Ching; Wan, Shu-Yen
2014-01-01
Berth allocation is the forefront operation performed when ships arrive at a port and is a critical task in container port optimization. Minimizing the time ships spend at berths constitutes an important objective of berth allocation problems. This study focuses on the discrete dynamic berth allocation problem (discrete DBAP), which aims to minimize total service time, and proposes an iterated greedy (IG) algorithm to solve it. The proposed IG algorithm is tested on three benchmark problem sets. Experimental results show that the proposed IG algorithm can obtain optimal solutions for all test instances of the first and second problem sets and outperforms the best-known solutions for 35 out of 90 test instances of the third problem set. PMID:25295295
NASA Astrophysics Data System (ADS)
Chang, Wei-Der; Yan, Jun-Juh
2006-10-01
In this paper, we propose a novel genetic algorithm (GA) with a multi-crossover fashion to estimate the associated coefficients for a class of nonlinear discrete-time multivariable dynamical systems. Unlike the traditional crossover method of using two chromosomes, the proposed method uses three chromosomes to achieve a crossover. According to the adjusting direction by crossing three chromosomes, more excellent offspring can be produced. To solve the identification problem of multivariable nonlinear discrete-time systems, each of estimated system coefficients represents a gene, and a collection of genes is referred to as a chromosome in the view of GA. The chromosomes in the population are then evolved using the proposed multi-crossover method. An illustrative example of multivariable nonlinear systems is given to demonstrate the effectiveness, as compared with the traditional crossover method, of the proposed method.
Discrete Adjoint-Based Design for Unsteady Turbulent Flows On Dynamic Overset Unstructured Grids
NASA Technical Reports Server (NTRS)
Nielsen, Eric J.; Diskin, Boris
2012-01-01
A discrete adjoint-based design methodology for unsteady turbulent flows on three-dimensional dynamic overset unstructured grids is formulated, implemented, and verified. The methodology supports both compressible and incompressible flows and is amenable to massively parallel computing environments. The approach provides a general framework for performing highly efficient and discretely consistent sensitivity analysis for problems involving arbitrary combinations of overset unstructured grids which may be static, undergoing rigid or deforming motions, or any combination thereof. General parent-child motions are also accommodated, and the accuracy of the implementation is established using an independent verification based on a complex-variable approach. The methodology is used to demonstrate aerodynamic optimizations of a wind turbine geometry, a biologically-inspired flapping wing, and a complex helicopter configuration subject to trimming constraints. The objective function for each problem is successfully reduced and all specified constraints are satisfied.
2014-01-01
Berth allocation is the forefront operation performed when ships arrive at a port and is a critical task in container port optimization. Minimizing the time ships spend at berths constitutes an important objective of berth allocation problems. This study focuses on the discrete dynamic berth allocation problem (discrete DBAP), which aims to minimize total service time, and proposes an iterated greedy (IG) algorithm to solve it. The proposed IG algorithm is tested on three benchmark problem sets. Experimental results show that the proposed IG algorithm can obtain optimal solutions for all test instances of the first and second problem sets and outperforms the best-known solutions for 35 out of 90 test instances of the third problem set. PMID:25295295
Modeling and 2-D discrete simulation of dislocation dynamics for plastic deformation of metal
NASA Astrophysics Data System (ADS)
Liu, Juan; Cui, Zhenshan; Ou, Hengan; Ruan, Liqun
2013-05-01
Two methods are employed in this paper to investigate the dislocation evolution during plastic deformation of metal. One method is dislocation dynamic simulation of two-dimensional discrete dislocation dynamics (2D-DDD), and the other is dislocation dynamics modeling by means of nonlinear analysis. As screw dislocation is prone to disappear by cross-slip, only edge dislocation is taken into account in simulation. First, an approach of 2D-DDD is used to graphically simulate and exhibit the collective motion of a large number of discrete dislocations. In the beginning, initial grains are generated in the simulation cells according to the mechanism of grain growth and the initial dislocation is randomly distributed in grains and relaxed under the internal stress. During the simulation process, the externally imposed stress, the long range stress contribution of all dislocations and the short range stress caused by the grain boundaries are calculated. Under the action of these forces, dislocations begin to glide, climb, multiply, annihilate and react with each other. Besides, thermal activation process is included. Through the simulation, the distribution of dislocation and the stress-strain curves can be obtained. On the other hand, based on the classic dislocation theory, the variation of the dislocation density with time is described by nonlinear differential equations. Finite difference method (FDM) is used to solve the built differential equations. The dislocation evolution at a constant strain rate is taken as an example to verify the rationality of the model.
Mesoscale dynamic coupling of finite- and discrete-element methods for fluid-particle interactions.
Srivastava, S; Yazdchi, K; Luding, S
2014-08-01
A new method for two-way fluid-particle coupling on an unstructured mesoscopically coarse mesh is presented. In this approach, we combine a (higher order) finite-element method (FEM) on the moving mesh for the fluid with a soft sphere discrete-element method for the particles. The novel feature of the proposed scheme is that the FEM mesh is a dynamic Delaunay triangulation based on the positions of the moving particles. Thus, the mesh can be multi-purpose: it provides (i) a framework for the discretization of the Navier-Stokes equations, (ii) a simple tool for detecting contacts between moving particles, (iii) a basis for coarse-graining or upscaling, and (iv) coupling with other physical fields (temperature, electromagnetic, etc.). This approach is suitable for a wide range of dilute and dense particulate flows, because the mesh resolution adapts with particle density in a given region. Two-way momentum exchange is implemented using semi-empirical drag laws akin to other popular approaches; for example, the discrete particle method, where a finite-volume solver on a coarser, fixed grid is used. We validate the methodology with several basic test cases, including single- and double-particle settling with analytical and empirical expectations, and flow through ordered and random porous media, when compared against finely resolved FEM simulations of flow through fixed arrays of particles. PMID:24982251
Parallel Discrete Molecular Dynamics Simulation With Speculation and In-Order Commitment*†
Khan, Md. Ashfaquzzaman; Herbordt, Martin C.
2011-01-01
Discrete molecular dynamics simulation (DMD) uses simplified and discretized models enabling simulations to advance by event rather than by timestep. DMD is an instance of discrete event simulation and so is difficult to scale: even in this multi-core era, all reported DMD codes are serial. In this paper we discuss the inherent difficulties of scaling DMD and present our method of parallelizing DMD through event-based decomposition. Our method is microarchitecture inspired: speculative processing of events exposes parallelism, while in-order commitment ensures correctness. We analyze the potential of this parallelization method for shared-memory multiprocessors. Achieving scalability required extensive experimentation with scheduling and synchronization methods to mitigate serialization. The speed-up achieved for a variety of system sizes and complexities is nearly 6× on an 8-core and over 9× on a 12-core processor. We present and verify analytical models that account for the achieved performance as a function of available concurrency and architectural limitations. PMID:21822327
Chaos analysis and delayed-feedback control in a discrete dynamic coupled map traffic model
NASA Astrophysics Data System (ADS)
Fang, Yaling; Shi, Zhongke
2015-03-01
The presence of chaos in traffic flow is studied using a modified discrete dynamic coupled map model which is derived from both the flow-density-speed fundamental diagram and Del Castillo's speed-density model. The modified model employs occupancy as its new variable and introduces a coupling strength with the consideration of effect of the front adjacent vehicle. And we analyze its stability of the control system and provide a procedure to design the decentralized delayed-feedback controllers for the traffic control system. These theoretical results are illustrated by numerical simulations.
Dissipative soliton dynamics in a discrete magnetic nano-dot chain
Lee, Kyeong-Dong; You, Chun-Yeol; Song, Hyon-Seok; Shin, Sung-Chul; Park, Byong-Guk
2014-02-03
Soliton dynamics is studied in a discrete magnetic nano-dot chain by means of micromagnetic simulations together with an analytic model equation. A soliton under a dissipative system is driven by an applied field. The field-driven dissipative soliton enhances its mobility nonlinearly, as the characteristic frequency and the intrinsic Gilbert damping decrease. During the propagation, the soliton emits spin waves which act as an extrinsic damping channel. The characteristic frequency, the maximum velocity, and the localization length of the soliton are found to be proportional to the threshold field, the threshold velocity, and the initial mobility, respectively.
An easy implementation of displacement calculations in 3D discrete dislocation dynamics codes
NASA Astrophysics Data System (ADS)
Fivel, Marc; Depres, Christophe
2014-10-01
Barnett's coordinate-free expression of the displacement field of a triangular loop in an isotropic media is revisited in a view to be implemented in 3D discrete dislocation dynamics codes. A general meshing procedure solving the problems of non-planar loops is presented. The method is user-friendly and can be used in numerical simulations since it gives the contribution of each dislocation segment to the global displacement field without defining the connectivity of closed loops. Easy to implement in parallel calculations, this method is successfully applied to large-scale simulations.
Recent Progress in Discrete Dislocation Dynamics and Its Applications to Micro Plasticity
NASA Astrophysics Data System (ADS)
Po, Giacomo; Mohamed, Mamdouh S.; Crosby, Tamer; Erel, Can; El-Azab, Anter; Ghoniem, Nasr
2014-09-01
We present a self-contained review of the discrete dislocation dynamics (DDD) method for the numerical investigation of plasticity in crystals, focusing on recent development and implementation progress. The review covers the theoretical foundations of DDD within the framework of incompatible elasticity, its numerical implementation via the nodal method, the extension of the method to finite domains and several implementation details. Applications of the method to current topics in micro-plasticity are presented, including the size effects in nano-indentation, the evolution of the dislocation microstructure in persistent slip bands, and the phenomenon of dislocation avalanches in micro-pillar compression.
Recent Progress in Discrete Dislocation Dynamics and Its Applications to Micro Plasticity
NASA Astrophysics Data System (ADS)
Po, Giacomo; Mohamed, Mamdouh S.; Crosby, Tamer; Erel, Can; El-Azab, Anter; Ghoniem, Nasr
2014-10-01
We present a self-contained review of the discrete dislocation dynamics (DDD) method for the numerical investigation of plasticity in crystals, focusing on recent development and implementation progress. The review covers the theoretical foundations of DDD within the framework of incompatible elasticity, its numerical implementation via the nodal method, the extension of the method to finite domains and several implementation details. Applications of the method to current topics in micro-plasticity are presented, including the size effects in nano-indentation, the evolution of the dislocation microstructure in persistent slip bands, and the phenomenon of dislocation avalanches in micro-pillar compression.
Davidchack, Ruslan L.
2010-12-10
We investigate the influence of numerical discretization errors on computed averages in a molecular dynamics simulation of TIP4P liquid water at 300 K coupled to different deterministic (Nose-Hoover and Nose-Poincare) and stochastic (Langevin) thermostats. We propose a couple of simple practical approaches to estimating such errors and taking them into account when computing the averages. We show that it is possible to obtain accurate measurements of various system quantities using step sizes of up to 70% of the stability threshold of the integrator, which for the system of TIP4P liquid water at 300 K corresponds to the step size of about 7 fs.
Neural networks for tracking of unknown SISO discrete-time nonlinear dynamic systems.
Aftab, Muhammad Saleheen; Shafiq, Muhammad
2015-11-01
This article presents a Lyapunov function based neural network tracking (LNT) strategy for single-input, single-output (SISO) discrete-time nonlinear dynamic systems. The proposed LNT architecture is composed of two feedforward neural networks operating as controller and estimator. A Lyapunov function based back propagation learning algorithm is used for online adjustment of the controller and estimator parameters. The controller and estimator error convergence and closed-loop system stability analysis is performed by Lyapunov stability theory. Moreover, two simulation examples and one real-time experiment are investigated as case studies. The achieved results successfully validate the controller performance. PMID:26456201
World-trade web: Topological properties, dynamics, and evolution
NASA Astrophysics Data System (ADS)
Fagiolo, Giorgio; Reyes, Javier; Schiavo, Stefano
2009-03-01
This paper studies the statistical properties of the web of import-export relationships among world countries using a weighted-network approach. We analyze how the distributions of the most important network statistics measuring connectivity, assortativity, clustering, and centrality have coevolved over time. We show that all node-statistic distributions and their correlation structure have remained surprisingly stable in the last 20years —and are likely to do so in the future. Conversely, the distribution of (positive) link weights is slowly moving from a log-normal density towards a power law. We also characterize the autoregressive properties of network-statistics dynamics. We find that network-statistics growth rates are well-proxied by fat-tailed densities like the Laplace or the asymmetric exponential power. Finally, we find that all our results are reasonably robust to a few alternative, economically meaningful, weighting schemes.
Dynamic Space for Rent: Using Commercial Web Hosting to Develop a Web 2.0 Intranet
ERIC Educational Resources Information Center
Hodgins, Dave
2010-01-01
The explosion of Web 2.0 into libraries has left many smaller academic libraries (and other libraries with limited computing resources or support) to work in the cloud using free Web applications. The use of commercial Web hosting is an innovative approach to the problem of inadequate local resources. While the idea of insourcing IT will seem…
Optimal discrete-time dynamic output-feedback design - A w-domain approach
NASA Technical Reports Server (NTRS)
Ha, Cheolkeun; Ly, Uy-Loi
1991-01-01
An alternative method for optimal digital control design is described in this paper. The method is based on the usage of the w-transform and has many attractive design features. One of these is its immediate connection with frequency loop-shaping techniques that are now popular and effective for multivariable control synthesis in continuous-time domain. Furthermore, any design algorithms originally developed for continuous-time systems can now be immediately extended to the discrete-time domain. The main results presented in this paper are the exact problem formulation and solution of an optimal discrete-time dynamic output-feedback design in the w-domain involving a quadratic performance index to random disturbances. In addition, necessary conditions for optimality are obtained for the numerical solution of the optimal output-feedback compensator design. A numerical example is presented illustrating its application to the design of a low-order dynamic compensator in a stability augmentation system of a commercial transport.
NASA Astrophysics Data System (ADS)
Doebrich, Marcus; Markstaller, Klaus; Karmrodt, Jens; Kauczor, Hans-Ulrich; Eberle, Balthasar; Weiler, Norbert; Thelen, Manfred; Schreiber, Wolfgang G.
2005-04-01
In this study, an algorithm was developed to measure the distribution of pulmonary time constants (TCs) from dynamic computed tomography (CT) data sets during a sudden airway pressure step up. Simulations with synthetic data were performed to test the methodology as well as the influence of experimental noise. Furthermore the algorithm was applied to in vivo data. In five pigs sudden changes in airway pressure were imposed during dynamic CT acquisition in healthy lungs and in a saline lavage ARDS model. The fractional gas content in the imaged slice (FGC) was calculated by density measurements for each CT image. Temporal variations of the FGC were analysed assuming a model with a continuous distribution of exponentially decaying time constants. The simulations proved the feasibility of the method. The influence of experimental noise could be well evaluated. Analysis of the in vivo data showed that in healthy lungs ventilation processes can be more likely characterized by discrete TCs whereas in ARDS lungs continuous distributions of TCs are observed. The temporal behaviour of lung inflation and deflation can be characterized objectively using the described new methodology. This study indicates that continuous distributions of TCs reflect lung ventilation mechanics more accurately compared to discrete TCs.
NASA Astrophysics Data System (ADS)
Crone, Joshua C.; Chung, Peter W.; Leiter, Kenneth W.; Knap, Jaroslaw; Aubry, Sylvie; Hommes, Gregg; Arsenlis, Athanasios
2014-04-01
Discrete dislocation dynamics (DD) approaches have proven useful in modeling the dynamics of large ensembles of dislocations. Continuing interest in finite body effects via image stresses has extended DD numerical approaches to improve the handling of surfaces. However, a physically accurate, yet computationally scalable, implementation has been elusive. This paper presents a new framework and implementation of a finite element-based discrete DD code that (1) treats arbitrarily shaped non-convex surfaces through image tractions, (2) allows for systematic refinement of the finite element mesh both in the bulk and on the surface and (3) provides a platform to scale to relatively larger and lengthier simulations. The approach is based on the capabilities of the Parallel Dislocation Simulator coupled through a distributed shared memory implementation for the calculation of large numbers of dislocation segments interacting with an independently large number of surface finite elements. Surface tracking approaches enable topological features at surfaces to be modeled. We verify the computed results via comparisons with analytical solutions for an infinite screw dislocation and prismatic loop near a surface and examine surface effects on a Frank-Read source. Convergence of the image force error with h- and p-refinement is shown to indicate the computational robustness. Additionally, through larger numerical experiments, we demonstrate the new capabilities in a three-dimensional elastic body of finite extent.
NASA Astrophysics Data System (ADS)
Martin, Hugo; Mangeney, Anne; Farin, Maxime; Richard, Patrick
2016-04-01
The mechanical behavior of granular flows is still an open issue. In particular, quantitative agreement between the detailed dynamics of the flow and laboratory experiments is necessary to better constrain the performance and limits of the models. We propose here to compare quantitatively the flow profiles and the force during granular column collapse simulated using Discrete Element Models and laboratory experiments. These small scale experiments are performed with dry granular material released initially from a cylinder on a sloping plane. The flow profiles and the acoustic signal generated by the granular impacts and stresses on the plane are recorded systematically [Farin et al., 2015]. These experiments are simulated using the Discrete Element Method Modys [Richard et al., 2000]. We show that the effect of the removing gate should be taken into account in the model in order to quantatively reproduce the flow dynamics. Furthermore we compare the simulated and observed acoustic signals that are generated by the fluctuating stresses exerted by the grains on the substrate in different frequency bands. [1] P. Richard et Luc Oger. 2000 Etude de la géométrie de milieux granulaires modèles tridimensionnels par simulation numérique. [2] Farin, M., Mangeney, A., Toussaint, R., De Rosny, J., Shapiro, N., Dewez, T., Hibert, C., Mathon, C., Sedan, O., Berger. 2015, Characterization of rockfalls from seismic signal: insights from laboratory experiments
Complex dynamics of a delayed discrete neural network of two nonidentical neurons
Chen, Yuanlong; Huang, Tingwen; Huang, Yu
2014-03-15
In this paper, we discover that a delayed discrete Hopfield neural network of two nonidentical neurons with self-connections and no self-connections can demonstrate chaotic behaviors. To this end, we first transform the model, by a novel way, into an equivalent system which has some interesting properties. Then, we identify the chaotic invariant set for this system and show that the dynamics of this system within this set is topologically conjugate to the dynamics of the full shift map with two symbols. This confirms chaos in the sense of Devaney. Our main results generalize the relevant results of Huang and Zou [J. Nonlinear Sci. 15, 291–303 (2005)], Kaslik and Balint [J. Nonlinear Sci. 18, 415–432 (2008)] and Chen et al. [Sci. China Math. 56(9), 1869–1878 (2013)]. We also give some numeric simulations to verify our theoretical results.
Zhong, Xiangnan; He, Haibo; Zhang, Huaguang; Wang, Zhanshan
2014-12-01
In this paper, we develop and analyze an optimal control method for a class of discrete-time nonlinear Markov jump systems (MJSs) with unknown system dynamics. Specifically, an identifier is established for the unknown systems to approximate system states, and an optimal control approach for nonlinear MJSs is developed to solve the Hamilton-Jacobi-Bellman equation based on the adaptive dynamic programming technique. We also develop detailed stability analysis of the control approach, including the convergence of the performance index function for nonlinear MJSs and the existence of the corresponding admissible control. Neural network techniques are used to approximate the proposed performance index function and the control law. To demonstrate the effectiveness of our approach, three simulation studies, one linear case, one nonlinear case, and one single link robot arm case, are used to validate the performance of the proposed optimal control method. PMID:25420238
Interplay between topology and dynamics in the World Trade Web
NASA Astrophysics Data System (ADS)
Garlaschelli, D.; di Matteo, T.; Aste, T.; Caldarelli, G.; Loffredo, M. I.
2007-05-01
We present an empirical analysis of the network formed by the trade relationships between all world countries, or World Trade Web (WTW). Each (directed) link is weighted by the amount of wealth flowing between two countries, and each country is characterized by the value of its Gross Domestic Product (GDP). By analysing a set of year-by-year data covering the time interval 1950 2000, we show that the dynamics of all GDP values and the evolution of the WTW (trade flow and topology) are tightly coupled. The probability that two countries are connected depends on their GDP values, supporting recent theoretical models relating network topology to the presence of a `hidden' variable (or fitness). On the other hand, the topology is shown to determine the GDP values due to the exchange between countries. This leads us to a new framework where the fitness value is a dynamical variable determining, and at the same time depending on, network topology in a continuous feedback.
Dynamic selection mechanism for quality of service aware web services
NASA Astrophysics Data System (ADS)
D'Mello, Demian Antony; Ananthanarayana, V. S.
2010-02-01
A web service is an interface of the software component that can be accessed by standard Internet protocols. The web service technology enables an application to application communication and interoperability. The increasing number of web service providers throughout the globe have produced numerous web services providing the same or similar functionality. This necessitates the use of tools and techniques to search the suitable services available over the Web. UDDI (universal description, discovery and integration) is the first initiative to find the suitable web services based on the requester's functional demands. However, the requester's requirements may also include non-functional aspects like quality of service (QoS). In this paper, the authors define a QoS model for QoS aware and business driven web service publishing and selection. The authors propose a QoS requirement format for the requesters, to specify their complex demands on QoS for the web service selection. The authors define a tree structure called quality constraint tree (QCT) to represent the requester's variety of requirements on QoS properties having varied preferences. The paper proposes a QoS broker based architecture for web service selection, which facilitates the requesters to specify their QoS requirements to select qualitatively optimal web service. A web service selection algorithm is presented, which ranks the functionally similar web services based on the degree of satisfaction of the requester's QoS requirements and preferences. The paper defines web service provider qualities to distinguish qualitatively competitive web services. The paper also presents the modelling and selection mechanism for the requester's alternative constraints defined on the QoS. The authors implement the QoS broker based system to prove the correctness of the proposed web service selection mechanism.
Sparsified-dynamics modeling of discrete point vortices with graph theory
NASA Astrophysics Data System (ADS)
Taira, Kunihiko; Nair, Aditya
2014-11-01
We utilize graph theory to derive a sparsified interaction-based model that captures unsteady point vortex dynamics. The present model builds upon the Biot-Savart law and keeps the number of vortices (graph nodes) intact and reduces the number of inter-vortex interactions (graph edges). We achieve this reduction in vortex interactions by spectral sparsification of graphs. This approach drastically reduces the computational cost to predict the dynamical behavior, sharing characteristics of reduced-order models. Sparse vortex dynamics are illustrated through an example of point vortex clusters interacting amongst themselves. We track the centroids of the individual vortex clusters to evaluate the error in bulk motion of the point vortices in the sparsified setup. To further improve the accuracy in predicting the nonlinear behavior of the vortices, resparsification strategies are employed for the sparsified interaction-based models. The model retains the nonlinearity of the interaction and also conserves the invariants of discrete vortex dynamics; namely the Hamiltonian, linear impulse, and angular impulse as well as circulation. Work supported by US Army Research Office (W911NF-14-1-0386) and US Air Force Office of Scientific Research (YIP: FA9550-13-1-0183).
A dynamic model of mobile concrete pump boom based on discrete time transfer matrix method
NASA Astrophysics Data System (ADS)
Ren, Wu; Wu, Yunxin; Zhang, Zhaowei
2013-12-01
Mobile concrete pump boom is typical multibody large-scale motion manipulator. Due to posture constantly change in working process, kinematic rule and dynamic characteristic are difficult to solve. A dynamics model of a mobile concrete pump boom is established based on discrete time transfer matrix method (DTTMM). The boom system is divided into sub-structure A and substructure B. Sub-structure A is composed by the 1st boom and hydraulic actuator as well as the support. And substructure B is consists of the other three booms and corresponding hydraulic actuators. In the model, the booms and links are regarded as rigid elements and the hydraulic cylinders are equivalent to spring-damper. The booms are driven by the controllable hydraulic actuators. The overall dynamic equation and transfer matrix of the model can be assembled by sub-structures A and B. To get a precise result, step size and integration parameters are studied then. Next the tip displacement is calculated and compared with the result of ADAMS software. The displacement and rotation angle curves of the proposed method fit well with the ADAMS model. Besides it is convenient in modeling and saves time. So it is suitable for mobile concrete pump boom real-time monitoring and dynamic analysis. All of these provide reference to boom optimize and engineering application of such mechanisms.
Event-driven Monte Carlo: Exact dynamics at all time scales for discrete-variable models
NASA Astrophysics Data System (ADS)
Mendoza-Coto, Alejandro; Díaz-Méndez, Rogelio; Pupillo, Guido
2016-06-01
We present an algorithm for the simulation of the exact real-time dynamics of classical many-body systems with discrete energy levels. In the same spirit of kinetic Monte Carlo methods, a stochastic solution of the master equation is found, with no need to define any other phase-space construction. However, unlike existing methods, the present algorithm does not assume any particular statistical distribution to perform moves or to advance the time, and thus is a unique tool for the numerical exploration of fast and ultra-fast dynamical regimes. By decomposing the problem in a set of two-level subsystems, we find a natural variable step size, that is well defined from the normalization condition of the transition probabilities between the levels. We successfully test the algorithm with known exact solutions for non-equilibrium dynamics and equilibrium thermodynamical properties of Ising-spin models in one and two dimensions, and compare to standard implementations of kinetic Monte Carlo methods. The present algorithm is directly applicable to the study of the real-time dynamics of a large class of classical Markovian chains, and particularly to short-time situations where the exact evolution is relevant.
Hunt, H. B.; Rosenkrantz, D. J.; Barrett, C. L.; Marathe, M. V.; Ravi, S. S.
2001-01-01
We identify several simple but powerful concepts, techniques, and results; and we use them to characterize the complexities of a number of basic problems II, that arise in the analysis and verification of the following models M of communicating automata and discrete dynamical systems: systems of communicating automata including both finite and infinite cellular automata, transition systems, discrete dynamical systems, and succinctly-specified finite automata. These concepts, techniques, and results are centered on the following: (1) reductions Of STATE-REACHABILITY problems, especially for very simple systems of communicating copies of a single simple finite automaton, (2) reductions of generalized CNF satisfiability problems [Sc78], especially to very simple communicating systems of copies of a few basic acyclic finite sequential machines, and (3) reductions of the EMPTINESS and EMPTINESS-OF-INTERSECTION problems, for several kinds of regular set descriptors. For systems of communicating automata and transition systems, the problems studied include: all equivalence relations and simulation preorders in the Linear-time/Branching-time hierarchies of equivalence relations and simulation preorders of [vG90, vG93], both without and with the hiding abstraction. For discrete dynamical systems, the problems studied include the INITIAL and BOUNDARY VALUE PROBLEMS (denoted IVPs and BVPs, respectively), for nonlinear difference equations over many different algebraic structures, e.g. all unitary rings, all finite unitary semirings, and all lattices. For succinctly specified finite automata, the problems studied also include the several problems studied in [AY98], e.g. the EMPTINESS, EMPTINESS-OF-INTERSECTION, EQUIVALENCE and CONTAINMENT problems. The concepts, techniques, and results presented unify and significantly extend many of the known results in the literature, e.g. [Wo86, Gu89, BPT91, GM92, Ra92, HT94, SH+96, AY98, AKY99, RH93, SM73, Hu73, HRS76, HR78], for
H. B. HUNT; D. J. ROSENKRANTS; ET AL
2001-03-01
We identify several simple but powerful concepts, techniques, and results; and we use them to characterize the complexities of a number of basic problems II, that arise in the analysis and verification of the following models M of communicating automata and discrete dynamical systems: systems of communicating automata including both finite and infinite cellular automata, transition systems, discrete dynamical systems, and succinctly-specified finite automata. These concepts, techniques, and results are centered on the following: (i) reductions Of STATE-REACHABILITY problems, especially for very simple systems of communicating copies of a single simple finite automaton, (ii) reductions of generalized CNF satisfiability problems [Sc78], especially to very simple communicating systems of copies of a few basic acyclic finite sequential machines, and (iii) reductions of the EMPTINESS and EMPTINESS-OF-INTERSECTION problems, for several kinds of regular set descriptors. For systems of communicating automata and transition systems, the problems studied include: all equivalence relations and simulation preorders in the Linear-time/Branching-time hierarchies of equivalence relations and simulation preorders of [vG90, vG93], both without and with the hiding abstraction. For discrete dynamical systems, the problems studied include the INITIAL and BOUNDARY VALUE PROBLEMS (denoted IVPs and BVPS, respectively), for nonlinear difference equations over many different algebraic structures, e.g. all unitary rings, all finite unitary semirings, and all lattices. For succinctly-specified finite automata, the problems studied also include the several problems studied in [AY98], e.g. the EMPTINESS, EMPTINESS-OF-INTERSECTION, EQUIVALENCE and CONTAINMENT problems. The concepts, techniques, and results presented unify and significantly extend many of the known results in the literature, e.g. [Wo86, Gu89, BPT91, GM92, Ra92, HT94, SH+96, AY98, AKY99, RH93, SM73, Hu73, HRS76, HR78], for
Kroenlein, Kenneth; Muzny, Chris D; Diky, Vladimir; Kazakov, Andrei F; Chirico, Robert D; Magee, Joseph W; Abdulagatov, Ilmutdin; Frenkel, Michael
2011-06-27
ThermoData Engine (TDE) is the first full-scale software implementation of the dynamic data evaluation concept, as reported recently in this journal. In the present paper, we describe the development of a World Wide Web-based interface to TDE evaluations of pure compound properties, including critical properties, phase boundary equilibria (vapor pressures, sublimation pressures, and crystal-liquid boundary pressures), densities, energetic properties, and transport properties. This includes development of a system for caching evaluation results to maintain high availability and an advanced window-in-window interface that leverages modern Web-browser technologies. Challenges associated with bringing the principal advantages of the TDE technology to the Web are described, as are compromises to maintain general access and speed of interaction while remaining true to the tenets of dynamic data evaluation. Future extensions of the interface and associated Web-services are outlined. PMID:21517125
Discrete states and carrier-phonon scattering in quantum dot population dynamics
Man, Minh Tan; Lee, Hong Seok
2015-01-01
The influence of the growth conditions of multilayer CdTe/ZnTe quantum dots (QDs) on Si substrate upon their carrier dynamics is studied using intensity integration and broadening photoluminescence. The unusual temperature dependence of the line broadening is explained using a model for interband transitions that involves a lowest discrete electronic state (1Se) with different discrete hole states (1S3/2 and 2S3/2) and a 1P transition. These transitions are expected to play a critical role in both the thermally activated energy and the line broadening of the QDs. We also demonstrate that a thermally activated transition between two different states occurs with band low-temperature quenching, with values separated by 5.8–16 meV. The main nonradiative process is thermal escape assisted by carrier scattering via emission of longitudinal phonons through the hole states at high temperature, with an average energy of 19.3–20.2 meV. PMID:25652600
Chan, Kung-Sik; Mysterud, Atle; Øritsland, Nils Are; Severinsen, Torbjørn; Stenseth, Nils Chr
2005-10-01
Climate at northern latitudes are currently changing both with regard to the mean and the temporal variability at any given site, increasing the frequency of extreme events such as cold and warm spells. Here we use a conceptually new modelling approach with two different dynamic terms of the climatic effects on a Svalbard reindeer population (the Brøggerhalvøya population) which underwent an extreme icing event ("locked pastures") with 80% reduction in population size during one winter (1993/94). One term captures the continuous and linear effect depending upon the Arctic Oscillation and another the discrete (rare) "event" process. The introduction of an "event" parameter describing the discrete extreme winter resulted in a more parsimonious model. Such an approach may be useful in strongly age-structured ungulate populations, with young and very old individuals being particularly prone to mortality factors during adverse conditions (resulting in a population structure that differs before and after extreme climatic events). A simulation study demonstrates that our approach is able to properly detect the ecological effects of such extreme climate events. PMID:16010537
NASA Astrophysics Data System (ADS)
Gurrutxaga-Lerma, Benat; Sutton, Adrian; Eakins, Daniel; Balint, Daniel; Dini, Daniele
2013-06-01
This talk intends to offer some insight as to how Discrete Dislocation Plasticity (DDP) can be adapted to simulate plastic relaxation processes under weak shock loading and high strain rates. In those circumstances, dislocations are believed to be the main cause of plastic relaxation in crystalline solids. Direct simulation of dislocations as the dynamic agents of plastic relaxation in those cases remains a challenge. DDP, where dislocations are modelled as discrete discontinuities in elastic continuum media, is often unable to adequately simulate plastic relaxation because it treats dislocation motion quasi-statically, thus neglecting the time-dependent nature of the elastic fields and assuming that they instantaneously acquire the shape and magnitude predicted by elastostatics. Under shock loading, this assumption leads to several artefacts that can only be overcome with a fully time-dependent formulation of the elastic fields. In this talk one of such formulations for the creation, annihilation and arbitrary motion of straight edge dislocations will be presented. These solutions are applied in a two-dimensional model of time-dependent plastic relaxation under shock loading, and some relevant results will be presented. EPSRC CDT in Theory and Simulation of Materials
NASA Astrophysics Data System (ADS)
Derakhshani, S. M.; Schott, D. L.; Lodewijks, G.
2013-06-01
Dust emissions can have significant effects on the human health, environment and industry equipment. Understanding the dust generation process helps to select a suitable dust preventing approach and also is useful to evaluate the environmental impact of dust emission. To describe these processes, numerical methods such as Computational Fluid Dynamics (CFD) are widely used, however nowadays particle based methods like Discrete Element Method (DEM) allow researchers to model interaction between particles and fluid flow. In this study, air flow over a stockpile, dust emission, erosion and surface deformation of granular material in the form of stockpile are studied by using DEM and CFD as a coupled method. Two and three dimensional simulations are respectively developed for CFD and DEM methods to minimize CPU time. The standard κ-ɛ turbulence model is used in a fully developed turbulent flow. The continuous gas phase and the discrete particle phase link to each other through gas-particle void fractions and momentum transfer. In addition to stockpile deformation, dust dispersion is studied and finally the accuracy of stockpile deformation results obtained by CFD-DEM modelling will be validated by the agreement with the existing experimental data.
Policy iteration adaptive dynamic programming algorithm for discrete-time nonlinear systems.
Liu, Derong; Wei, Qinglai
2014-03-01
This paper is concerned with a new discrete-time policy iteration adaptive dynamic programming (ADP) method for solving the infinite horizon optimal control problem of nonlinear systems. The idea is to use an iterative ADP technique to obtain the iterative control law, which optimizes the iterative performance index function. The main contribution of this paper is to analyze the convergence and stability properties of policy iteration method for discrete-time nonlinear systems for the first time. It shows that the iterative performance index function is nonincreasingly convergent to the optimal solution of the Hamilton-Jacobi-Bellman equation. It is also proven that any of the iterative control laws can stabilize the nonlinear systems. Neural networks are used to approximate the performance index function and compute the optimal control law, respectively, for facilitating the implementation of the iterative ADP algorithm, where the convergence of the weight matrices is analyzed. Finally, the numerical results and analysis are presented to illustrate the performance of the developed method. PMID:24807455
Food web dynamics affect Northeast Arctic cod recruitment
Hjermann, Dag Ø; Bogstad, Bjarte; Eikeset, Anne Maria; Ottersen, Geir; Gjøsæter, Harald; Stenseth, Nils Chr
2006-01-01
Proper management of ecosystems requires an understanding of both the species interactions as well as the effect of climate variation. However, a common problem is that the available time-series are of different lengths. Here, we present a general approach for studying the dynamic structure of such interactions. Specifically, we analyse the recruitment of the world's largest cod stock, the Northeast Arctic cod. Studies based on data starting in the 1970–1980s indicate that this stock is affected by temperature through a variety of pathways. However, the value of such studies is somewhat limited by the fact that they are based on a quite specific ecological and climatic situation. Recently, this stock has consisted of fairly young fish and the spawning stock has consisted of relatively few age groups. In this study, we develop a model for the effect of capelin (the cod's main prey) and herring on cod recruitment since 1973. Based on this model, we analyse data on cod, herring and temperature going back to 1921 and find that food-web effects explain a significant part of the cod recruitment variation back to around 1950. PMID:17254990
Web Services for Dynamic Coloring of UAVSAR Images
NASA Astrophysics Data System (ADS)
Wang, Jun; Pierce, Marlon; Donnellan, Andrea; Parker, Jay
2015-08-01
QuakeSim has implemented a service-based Geographic Information System to enable users to access large amounts of Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) data through an online interface. The QuakeSim Interferometric Synthetic Aperture Radar (InSAR) profile tool calculates radar-observed displacement (from an unwrapped interferogram product) along user-specified lines. Pre-rendered thumbnails with InSAR fringe patterns are used to display interferogram and unwrapped phase images on a Google Map in the InSAR profile tool. One challenge with this tool lies in the user visually identifying regions of interest when drawing the profile line. This requires that the user correctly interpret the InSAR imagery, which currently uses fringe patterns. The mapping between pixel color and pixel value is not a one-to-one relationship from the InSAR fringe pattern, and it causes difficulty in understanding general displacement information for QuakeSim users. The goal of this work is to generate color maps that directly reflect the pixel values (displacement) as an addition to the pre-rendered images. Because of an extremely uneven distribution of pixel values on an InSAR image, a histogram-based, nonlinear color template generation algorithm is currently under development. A web service enables on-the-fly coloring of UAVSAR images with dynamically generated color templates.
Penalised logistic regression and dynamic prediction for discrete-time recurrent event data.
Elgmati, Entisar; Fiaccone, Rosemeire L; Henderson, R; Matthews, John N S
2015-10-01
We consider methods for the analysis of discrete-time recurrent event data, when interest is mainly in prediction. The Aalen additive model provides an extremely simple and effective method for the determination of covariate effects for this type of data, especially in the presence of time-varying effects and time varying covariates, including dynamic summaries of prior event history. The method is weakened for predictive purposes by the presence of negative estimates. The obvious alternative of a standard logistic regression analysis at each time point can have problems of stability when event frequency is low and maximum likelihood estimation is used. The Firth penalised likelihood approach is stable but in removing bias in regression coefficients it introduces bias into predicted event probabilities. We propose an alterative modified penalised likelihood, intermediate between Firth and no penalty, as a pragmatic compromise between stability and bias. Illustration on two data sets is provided. PMID:25626559
Dynamical Localization for Discrete and Continuous Random Schrödinger Operators
NASA Astrophysics Data System (ADS)
Germinet, F.; De Bièvre, S.
We show for a large class of random Schrödinger operators Ho on and on that dynamical localization holds, i.e. that, with probability one, for a suitable energy interval I and for q a positive real,
NASA Astrophysics Data System (ADS)
Rocha, M. C.; Saraiva, J. T.
2012-10-01
The basic objective of Transmission Expansion Planning (TEP) is to schedule a number of transmission projects along an extended planning horizon minimizing the network construction and operational costs while satisfying the requirement of delivering power safely and reliably to load centres along the horizon. This principle is quite simple, but the complexity of the problem and the impact on society transforms TEP on a challenging issue. This paper describes a new approach to solve the dynamic TEP problem, based on an improved discrete integer version of the Evolutionary Particle Swarm Optimization (EPSO) meta-heuristic algorithm. The paper includes sections describing in detail the EPSO enhanced approach, the mathematical formulation of the TEP problem, including the objective function and the constraints, and a section devoted to the application of the developed approach to this problem. Finally, the use of the developed approach is illustrated using a case study based on the IEEE 24 bus 38 branch test system.
Discrete Molecular Dynamics Can Predict Helical Prestructured Motifs in Disordered Proteins
Han, Kyou-Hoon; Dokholyan, Nikolay V.; Tompa, Péter; Kalmár, Lajos; Hegedűs, Tamás
2014-01-01
Intrinsically disordered proteins (IDPs) lack a stable tertiary structure, but their short binding regions termed Pre-Structured Motifs (PreSMo) can form transient secondary structure elements in solution. Although disordered proteins are crucial in many biological processes and designing strategies to modulate their function is highly important, both experimental and computational tools to describe their conformational ensembles and the initial steps of folding are sparse. Here we report that discrete molecular dynamics (DMD) simulations combined with replica exchange (RX) method efficiently samples the conformational space and detects regions populating α-helical conformational states in disordered protein regions. While the available computational methods predict secondary structural propensities in IDPs based on the observation of protein-protein interactions, our ab initio method rests on physical principles of protein folding and dynamics. We show that RX-DMD predicts α-PreSMos with high confidence confirmed by comparison to experimental NMR data. Moreover, the method also can dissect α-PreSMos in close vicinity to each other and indicate helix stability. Importantly, simulations with disordered regions forming helices in X-ray structures of complexes indicate that a preformed helix is frequently the binding element itself, while in other cases it may have a role in initiating the binding process. Our results indicate that RX-DMD provides a breakthrough in the structural and dynamical characterization of disordered proteins by generating the structural ensembles of IDPs even when experimental data are not available. PMID:24763499
Dynamics of Investor Attention on the Social Web
ERIC Educational Resources Information Center
Li, Xian
2013-01-01
The World Wide Web has been revolutionizing how investors produce and consume information while participating in financial markets. Both the amount of information and the speed it flows around have achieved unprecedented magnitudes. The preeminent change is the growth of investor communities on the social web, which give rise to multidimensional…
Web-based Toolkit for Dynamic Generation of Data Processors
NASA Astrophysics Data System (ADS)
Patel, J.; Dascalu, S.; Harris, F. C.; Benedict, K. K.; Gollberg, G.; Sheneman, L.
2011-12-01
All computation-intensive scientific research uses structured datasets, including hydrology and all other types of climate-related research. When it comes to testing their hypotheses, researchers might use the same dataset differently, and modify, transform, or convert it to meet their research needs. Currently, many researchers spend a good amount of time performing data processing and building tools to speed up this process. They might routinely repeat the same process activities for new research projects, spending precious time that otherwise could be dedicated to analyzing and interpreting the data. Numerous tools are available to run tests on prepared datasets and many of them work with datasets in different formats. However, there is still a significant need for applications that can comprehensively handle data transformation and conversion activities and help prepare the various processed datasets required by the researchers. We propose a web-based application (a software toolkit) that dynamically generates data processors capable of performing data conversions, transformations, and customizations based on user-defined mappings and selections. As a first step, the proposed solution allows the users to define various data structures and, in the next step, can select various file formats and data conversions for their datasets of interest. In a simple scenario, the core of the proposed web-based toolkit allows the users to define direct mappings between input and output data structures. The toolkit will also support defining complex mappings involving the use of pre-defined sets of mathematical, statistical, date/time, and text manipulation functions. Furthermore, the users will be allowed to define logical cases for input data filtering and sampling. At the end of the process, the toolkit is designed to generate reusable source code and executable binary files for download and use by the scientists. The application is also designed to store all data
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.
A Dynamic Recommender System for Improved Web Usage Mining and CRM Using Swarm Intelligence.
Alphy, Anna; Prabakaran, S
2015-01-01
In modern days, to enrich e-business, the websites are personalized for each user by understanding their interests and behavior. The main challenges of online usage data are information overload and their dynamic nature. In this paper, to address these issues, a WebBluegillRecom-annealing dynamic recommender system that uses web usage mining techniques in tandem with software agents developed for providing dynamic recommendations to users that can be used for customizing a website is proposed. The proposed WebBluegillRecom-annealing dynamic recommender uses swarm intelligence from the foraging behavior of a bluegill fish. It overcomes the information overload by handling dynamic behaviors of users. Our dynamic recommender system was compared against traditional collaborative filtering systems. The results show that the proposed system has higher precision, coverage, F1 measure, and scalability than the traditional collaborative filtering systems. Moreover, the recommendations given by our system overcome the overspecialization problem by including variety in recommendations. PMID:26229978
A Dynamic Recommender System for Improved Web Usage Mining and CRM Using Swarm Intelligence
Alphy, Anna; Prabakaran, S.
2015-01-01
In modern days, to enrich e-business, the websites are personalized for each user by understanding their interests and behavior. The main challenges of online usage data are information overload and their dynamic nature. In this paper, to address these issues, a WebBluegillRecom-annealing dynamic recommender system that uses web usage mining techniques in tandem with software agents developed for providing dynamic recommendations to users that can be used for customizing a website is proposed. The proposed WebBluegillRecom-annealing dynamic recommender uses swarm intelligence from the foraging behavior of a bluegill fish. It overcomes the information overload by handling dynamic behaviors of users. Our dynamic recommender system was compared against traditional collaborative filtering systems. The results show that the proposed system has higher precision, coverage, F1 measure, and scalability than the traditional collaborative filtering systems. Moreover, the recommendations given by our system overcome the overspecialization problem by including variety in recommendations. PMID:26229978
Dorn, Martin; Hekmat, Dariusch
2016-03-01
Preparative packed-bed chromatography using polymer-based, compressible, porous resins is a powerful method for purification of macromolecular bioproducts. During operation, a complex, hysteretic, thus, history-dependent packed bed behavior is often observed but theoretical understanding of the causes is limited. Therefore, a rigorous modeling approach of the chromatography column on the particle scale has been made which takes into account interparticle micromechanics and fluid-particle interactions for the first time. A three-dimensional deterministic model was created by applying Computational Fluid Dynamics (CFD) coupled with the Discrete Element Method (DEM). The column packing behavior during either flow or mechanical compression was investigated in-silico and in laboratory experiments. A pronounced axial compression-relaxation profile was identified that differed for both compression strategies. Void spaces were clearly visible in the packed bed after compression. It was assumed that the observed bed inhomogeneity was because of a force-chain network at the particle scale. The simulation satisfactorily reproduced the measured behavior regarding packing compression as well as pressure-flow dependency. Furthermore, the particle Young's modulus and particle-wall friction as well as interparticle friction were identified as crucial parameters affecting packing dynamics. It was concluded that compaction of the chromatographic bed is rather because of particle rearrangement than particle deformation. © 2015 American Institute of Chemical Engineers Biotechnol. Prog., 32:363-371, 2016. PMID:26588806
Discrete-time nonlinear HJB solution using approximate dynamic programming: convergence proof.
Al-Tamimi, Asma; Lewis, Frank L; Abu-Khalaf, Murad
2008-08-01
Convergence of the value-iteration-based heuristic dynamic programming (HDP) algorithm is proven in the case of general nonlinear systems. That is, it is shown that HDP converges to the optimal control and the optimal value function that solves the Hamilton-Jacobi-Bellman equation appearing in infinite-horizon discrete-time (DT) nonlinear optimal control. It is assumed that, at each iteration, the value and action update equations can be exactly solved. The following two standard neural networks (NN) are used: a critic NN is used to approximate the value function, whereas an action network is used to approximate the optimal control policy. It is stressed that this approach allows the implementation of HDP without knowing the internal dynamics of the system. The exact solution assumption holds for some classes of nonlinear systems and, specifically, in the specific case of the DT linear quadratic regulator (LQR), where the action is linear and the value quadratic in the states and NNs have zero approximation error. It is stressed that, for the LQR, HDP may be implemented without knowing the system A matrix by using two NNs. This fact is not generally appreciated in the folklore of HDP for the DT LQR, where only one critic NN is generally used. PMID:18632382
NASA Astrophysics Data System (ADS)
Liu, Wenyan; Tang, Z. P.; Liu, Yunxin
2000-04-01
In recent years, more attention has been paid to a better understanding of the failure behavior and mechanism of heterogeneous materials at the meso-scale level. In this paper, the crack initiation and development in epoxy composites reinforced with short steel fibers under dynamic loading were simulated and analyzed with the 2D Discrete Meso-Element Dynamic Method. Results show that the damage process depends greatly on the binding property between matrix and fibers.
Discrete two-sex models of population dynamics: On modelling the mating function
NASA Astrophysics Data System (ADS)
Bessa-Gomes, Carmen; Legendre, Stéphane; Clobert, Jean
2010-09-01
Although sexual reproduction has long been a central subject of theoretical ecology, until recently its consequences for population dynamics were largely overlooked. This is now changing, and many studies have addressed this issue, showing that when the mating system is taken into account, the population dynamics depends on the relative abundance of males and females, and is non-linear. Moreover, sexual reproduction increases the extinction risk, namely due to the Allee effect. Nevertheless, different studies have identified diverse potential consequences, depending on the choice of mating function. In this study, we investigate the consequences of three alternative mating functions that are frequently used in discrete population models: the minimum; the harmonic mean; and the modified harmonic mean. We consider their consequences at three levels: on the probability that females will breed; on the presence and intensity of the Allee effect; and on the extinction risk. When we consider the harmonic mean, the number of times the individuals of the least abundant sex mate exceeds their mating potential, which implies that with variable sex-ratios the potential reproductive rate is no longer under the modeller's control. Consequently, the female breeding probability exceeds 1 whenever the sex-ratio is male-biased, which constitutes an obvious problem. The use of the harmonic mean is thus only justified if we think that this parameter should be re-defined in order to represent the females' breeding rate and the fact that females may reproduce more than once per breeding season. This phenomenon buffers the Allee effect, and reduces the extinction risk. However, when we consider birth-pulse populations, such a phenomenon is implausible because the number of times females can reproduce per birth season is limited. In general, the minimum or modified harmonic mean mating functions seem to be more suitable for assessing the impact of mating systems on population dynamics.
Discrete Molecular Dynamics Distinguishes Nativelike Binding Poses from Decoys in Difficult Targets
Proctor, Elizabeth A.; Yin, Shuangye; Tropsha, Alexander; Dokholyan, Nikolay V.
2012-01-01
Virtual screening is one of the major tools used in computer-aided drug discovery. In structure-based virtual screening, the scoring function is critical to identifying the correct docking pose and accurately predicting the binding affinities of compounds. However, the performance of existing scoring functions has been shown to be uneven for different targets, and some important drug targets have proven especially challenging. In these targets, scoring functions cannot accurately identify the native or near-native binding pose of the ligand from among decoy poses, which affects both the accuracy of the binding affinity prediction and the ability of virtual screening to identify true binders in chemical libraries. Here, we present an approach to discriminating native poses from decoys in difficult targets for which several scoring functions failed to correctly identify the native pose. Our approach employs Discrete Molecular Dynamics simulations to incorporate protein-ligand dynamics and the entropic effects of binding. We analyze a collection of poses generated by docking and find that the residence time of the ligand in the native and nativelike binding poses is distinctly longer than that in decoy poses. This finding suggests that molecular simulations offer a unique approach to distinguishing the native (or nativelike) binding pose from decoy poses that cannot be distinguished using scoring functions that evaluate static structures. The success of our method emphasizes the importance of protein-ligand dynamics in the accurate determination of the binding pose, an aspect that is not addressed in typical docking and scoring protocols. PMID:22225808
NASA Technical Reports Server (NTRS)
Saether, Erik; Hochhalter, Jacob D.; Glaessgen, Edward H.; Mishin, Yuri
2014-01-01
A multiscale modeling methodology is developed for structurally-graded material microstructures. Molecular dynamic (MD) simulations are performed at the nanoscale to determine fundamental failure mechanisms and quantify material constitutive parameters. These parameters are used to calibrate material processes at the mesoscale using discrete dislocation dynamics (DD). Different grain boundary interactions with dislocations are analyzed using DD to predict grain-size dependent stress-strain behavior. These relationships are mapped into crystal plasticity (CP) parameters to develop a computationally efficient finite element-based DD/CP model for continuum-level simulations and complete the multiscale analysis by predicting the behavior of macroscopic physical specimens. The present analysis is focused on simulating the behavior of a graded microstructure in which grain sizes are on the order of nanometers in the exterior region and transition to larger, multi-micron size in the interior domain. This microstructural configuration has been shown to offer improved mechanical properties over homogeneous coarse-grained materials by increasing yield stress while maintaining ductility. Various mesoscopic polycrystal models of structurally-graded microstructures are generated, analyzed and used as a benchmark for comparison between multiscale DD/CP model and DD predictions. A final series of simulations utilize the DD/CP analysis method exclusively to study macroscopic models that cannot be analyzed by MD or DD methods alone due to the model size.
Sensor Web Dynamic Measurement Techniques and Adaptive Observing Strategies
NASA Technical Reports Server (NTRS)
Talabac, Stephen J.
2004-01-01
Sensor Web observing systems may have the potential to significantly improve our ability to monitor, understand, and predict the evolution of rapidly evolving, transient, or variable environmental features and events. This improvement will come about by integrating novel data collection techniques, new or improved instruments, emerging communications technologies and protocols, sensor mark-up languages, and interoperable planning and scheduling systems. In contrast to today's observing systems, "event-driven" sensor webs will synthesize real- or near-real time measurements and information from other platforms and then react by reconfiguring the platforms and instruments to invoke new measurement modes and adaptive observation strategies. Similarly, "model-driven" sensor webs will utilize environmental prediction models to initiate targeted sensor measurements or to use a new observing strategy. The sensor web concept contrasts with today's data collection techniques and observing system operations concepts where independent measurements are made by remote sensing and in situ platforms that do not share, and therefore cannot act upon, potentially useful complementary sensor measurement data and platform state information. This presentation describes NASA's view of event-driven and model-driven Sensor Webs and highlights several research and development activities at the Goddard Space Flight Center.
Using Web Services and XML Harvesting to Achieve a Dynamic Web Site. Computers in Small Libraries
ERIC Educational Resources Information Center
Roberts, Gary
2005-01-01
Exploiting and contextualizing free information is a natural part of library culture. In this column, Gary Roberts, the information systems and reference librarian at Herrick Library, Alfred University in Alfred, NY, describes how to use XML content on a Web site to link to hundreds of free and useful resources. He gives a general overview of the…
Principles of Discrete Time Mechanics
NASA Astrophysics Data System (ADS)
Jaroszkiewicz, George
2014-04-01
1. Introduction; 2. The physics of discreteness; 3. The road to calculus; 4. Temporal discretization; 5. Discrete time dynamics architecture; 6. Some models; 7. Classical cellular automata; 8. The action sum; 9. Worked examples; 10. Lee's approach to discrete time mechanics; 11. Elliptic billiards; 12. The construction of system functions; 13. The classical discrete time oscillator; 14. Type 2 temporal discretization; 15. Intermission; 16. Discrete time quantum mechanics; 17. The quantized discrete time oscillator; 18. Path integrals; 19. Quantum encoding; 20. Discrete time classical field equations; 21. The discrete time Schrodinger equation; 22. The discrete time Klein-Gordon equation; 23. The discrete time Dirac equation; 24. Discrete time Maxwell's equations; 25. The discrete time Skyrme model; 26. Discrete time quantum field theory; 27. Interacting discrete time scalar fields; 28. Space, time and gravitation; 29. Causality and observation; 30. Concluding remarks; Appendix A. Coherent states; Appendix B. The time-dependent oscillator; Appendix C. Quaternions; Appendix D. Quantum registers; References; Index.
Investigation on Two-phase Flow Dynamics with Discrete Bubble Model
NASA Astrophysics Data System (ADS)
Ami, Takeyuki; Umekawa, Hisashi; Ozawa, Mamoru; Shoji, Masahiro
Conventional modeling including drift-flux model and two-fluid model is based on “continuous flow hypothesis”, being constructed by time-averaging, and thus both phases are defined in every spatio-temporal space. This makes it possible to apply to a variety of two-phase flow dynamics, while the intrinsic void fraction fluctuations, typically observed in slug and churn flows, are hardly simulated. In order to break through such a problem caused by time-averaging, discrete bubble model based on one-dimensional mass conservation equation, i.e. void propagation equation, has been developed. This model takes into account, as momentum effects, the wake effect induced by preceding bubbles, the local pressure fluctuation and the compressibility of gas phase together with the phase re-distribution due to geometrical constrains. Thus obtained spatio-temporal fluctuation characteristics of void fraction well simulated inherent two-phase behavior not only in a steady flow but also in an oscillatory flow.
Prostate segmentation algorithm using dyadic wavelet transform and discrete dynamic contour
NASA Astrophysics Data System (ADS)
Chiu, Bernard; Freeman, George H.; Salama, M. M. A.; Fenster, Aaron
2004-11-01
Knowing the location and the volume of the prostate is important for ultrasound-guided prostate brachytherapy, a commonly used prostate cancer treatment method. The prostate boundary must be segmented before a dose plan can be obtained. However, manual segmentation is arduous and time consuming. This paper introduces a semi-automatic segmentation algorithm based on the dyadic wavelet transform (DWT) and the discrete dynamic contour (DDC). A spline interpolation method is used to determine the initial contour based on four user-defined initial points. The DDC model then refines the initial contour based on the approximate coefficients and the wavelet coefficients generated using the DWT. The DDC model is executed under two settings. The coefficients used in these two settings are derived using smoothing functions with different sizes. A selection rule is used to choose the best contour based on the contours produced in these two settings. The accuracy of the final contour produced by the proposed algorithm is evaluated by comparing it with the manual contour outlined by an expert observer. A total of 114 2D TRUS images taken for six different patients scheduled for brachytherapy were segmented using the proposed algorithm. The average difference between the contour segmented using the proposed algorithm and the manually outlined contour is less than 3 pixels.
Prostate segmentation algorithm using dyadic wavelet transform and discrete dynamic contour.
Chiu, Bernard; Freeman, George H; Salama, M M A; Fenster, Aaron
2004-11-01
Knowing the location and the volume of the prostate is important for ultrasound-guided prostate brachytherapy, a commonly used prostate cancer treatment method. The prostate boundary must be segmented before a dose plan can be obtained. However, manual segmentation is arduous and time consuming. This paper introduces a semi-automatic segmentation algorithm based on the dyadic wavelet transform (DWT) and the discrete dynamic contour (DDC). A spline interpolation method is used to determine the initial contour based on four user-defined initial points. The DDC model then refines the initial contour based on the approximate coefficients and the wavelet coefficients generated using the DWT. The DDC model is executed under two settings. The coefficients used in these two settings are derived using smoothing functions with different sizes. A selection rule is used to choose the best contour based on the contours produced in these two settings. The accuracy of the final contour produced by the proposed algorithm is evaluated by comparing it with the manual contour outlined by an expert observer. A total of 114 2D TRUS images taken for six different patients scheduled for brachytherapy were segmented using the proposed algorithm. The average difference between the contour segmented using the proposed algorithm and the manually outlined contour is less than 3 pixels. PMID:15584529
Dynamical and statistical behavior of discrete combustion waves: a theoretical and numerical study.
Bharath, Naine Tarun; Rashkovskiy, Sergey A; Tewari, Surya P; Gundawar, Manoj Kumar
2013-04-01
We present a detailed theoretical and numerical study of combustion waves in a discrete one-dimensional disordered system. The distances between neighboring reaction cells were modeled with a gamma distribution. The results show that the random structure of the microheterogeneous system plays a crucial role in the dynamical and statistical behavior of the system. This is a consequence of the nonlinear interaction of the random structure of the system with the thermal wave. An analysis of the experimental data on the combustion of a gasless system (Ti + xSi) and a wide range of thermite systems was performed in view of the developed model. We have shown that the burning rate of the powder system sensitively depends on its internal structure. The present model allows for reproducing theoretically the experimental data for a wide range of pyrotechnic mixtures. We show that Arrhenius' macrokinetics at combustion of disperse systems can take place even in the absence of Arrhenius' microkinetics; it can have a purely thermal nature and be related to their heterogeneity and to the existence of threshold temperature. It is also observed that the combustion of disperse systems always occurs in the microheterogeneous mode according to the relay-race mechanism. PMID:23679470
Dynamical and statistical behavior of discrete combustion waves: A theoretical and numerical study
NASA Astrophysics Data System (ADS)
Bharath, Naine Tarun; Rashkovskiy, Sergey A.; Tewari, Surya P.; Gundawar, Manoj Kumar
2013-04-01
We present a detailed theoretical and numerical study of combustion waves in a discrete one-dimensional disordered system. The distances between neighboring reaction cells were modeled with a gamma distribution. The results show that the random structure of the microheterogeneous system plays a crucial role in the dynamical and statistical behavior of the system. This is a consequence of the nonlinear interaction of the random structure of the system with the thermal wave. An analysis of the experimental data on the combustion of a gasless system (Ti + xSi) and a wide range of thermite systems was performed in view of the developed model. We have shown that the burning rate of the powder system sensitively depends on its internal structure. The present model allows for reproducing theoretically the experimental data for a wide range of pyrotechnic mixtures. We show that Arrhenius’ macrokinetics at combustion of disperse systems can take place even in the absence of Arrhenius’ microkinetics; it can have a purely thermal nature and be related to their heterogeneity and to the existence of threshold temperature. It is also observed that the combustion of disperse systems always occurs in the microheterogeneous mode according to the relay-race mechanism.
Finite-approximation-error-based discrete-time iterative adaptive dynamic programming.
Wei, Qinglai; Wang, Fei-Yue; Liu, Derong; Yang, Xiong
2014-12-01
In this paper, a new iterative adaptive dynamic programming (ADP) algorithm is developed to solve optimal control problems for infinite horizon discrete-time nonlinear systems with finite approximation errors. First, a new generalized value iteration algorithm of ADP is developed to make the iterative performance index function converge to the solution of the Hamilton-Jacobi-Bellman equation. The generalized value iteration algorithm permits an arbitrary positive semi-definite function to initialize it, which overcomes the disadvantage of traditional value iteration algorithms. When the iterative control law and iterative performance index function in each iteration cannot accurately be obtained, for the first time a new "design method of the convergence criteria" for the finite-approximation-error-based generalized value iteration algorithm is established. A suitable approximation error can be designed adaptively to make the iterative performance index function converge to a finite neighborhood of the optimal performance index function. Neural networks are used to implement the iterative ADP algorithm. Finally, two simulation examples are given to illustrate the performance of the developed method. PMID:25265640
Food web dynamics in a seasonally varying wetland
DeAngelis, D.L.; Trexler, J.C.; Donalson, D.D.
2008-01-01
A spatially explicit model is developed to simulate the small fish community and its underlying food web, in the freshwater marshes of the Everglades. The community is simplified to a few small fish species feeding on periphyton and invertebrates. Other compartments are detritus, crayfish, and a piscivorous fish species. This unit food web model is applied to each of the 10,000 spatial cells on a 100 x 100 pixel landscape. Seasonal variation in water level is assumed and rules are assigned for fish movement in response to rising and falling water levels, which can cause many spatial cells to alternate between flooded and dry conditions. It is shown that temporal variations of water level on a spatially heterogeneous landscape can maintain at least three competing fish species. In addition, these environmental factors can strongly affect the temporal variation of the food web caused by top-down control from the piscivorous fish.
NASA Astrophysics Data System (ADS)
Bertin, N.; Upadhyay, M. V.; Pradalier, C.; Capolungo, L.
2015-09-01
In this paper, we propose a novel full-field approach based on the fast Fourier transform (FFT) technique to compute mechanical fields in periodic discrete dislocation dynamics (DDD) simulations for anisotropic materials: the DDD-FFT approach. By coupling the FFT-based approach to the discrete continuous model, the present approach benefits from the high computational efficiency of the FFT algorithm, while allowing for a discrete representation of dislocation lines. It is demonstrated that the computational time associated with the new DDD-FFT approach is significantly lower than that of current DDD approaches when large number of dislocation segments are involved for isotropic and anisotropic elasticity, respectively. Furthermore, for fine Fourier grids, the treatment of anisotropic elasticity comes at a similar computational cost to that of isotropic simulation. Thus, the proposed approach paves the way towards achieving scale transition from DDD to mesoscale plasticity, especially due to the method’s ability to incorporate inhomogeneous elasticity.
Bengoetxea, Ana; Leurs, Françoise; Hoellinger, Thomas; Cebolla, Ana M.; Dan, Bernard; McIntyre, Joseph; Cheron, Guy
2014-01-01
In this study we employed a dynamic recurrent neural network (DRNN) in a novel fashion to reveal characteristics of control modules underlying the generation of muscle activations when drawing figures with the outstretched arm. We asked healthy human subjects to perform four different figure-eight movements in each of two workspaces (frontal plane and sagittal plane). We then trained a DRNN to predict the movement of the wrist from information in the EMG signals from seven different muscles. We trained different instances of the same network on a single movement direction, on all four movement directions in a single movement plane, or on all eight possible movement patterns and looked at the ability of the DRNN to generalize and predict movements for trials that were not included in the training set. Within a single movement plane, a DRNN trained on one movement direction was not able to predict movements of the hand for trials in the other three directions, but a DRNN trained simultaneously on all four movement directions could generalize across movement directions within the same plane. Similarly, the DRNN was able to reproduce the kinematics of the hand for both movement planes, but only if it was trained on examples performed in each one. As we will discuss, these results indicate that there are important dynamical constraints on the mapping of EMG to hand movement that depend on both the time sequence of the movement and on the anatomical constraints of the musculoskeletal system. In a second step, we injected EMG signals constructed from different synergies derived by the PCA in order to identify the mechanical significance of each of these components. From these results, one can surmise that discrete-rhythmic movements may be constructed from three different fundamental modules, one regulating the co-activation of all muscles over the time span of the movement and two others elliciting patterns of reciprocal activation operating in orthogonal directions
The Spider's Web: Creativity and Survival in Dynamic Balance.
ERIC Educational Resources Information Center
Cohen, Bill
2001-01-01
The spider's web is presented as a model for Indigenous education and community transformation, grounded in Okanagan philosophy. Children are at the center and benefit from the influence of extended family and community. The model's relevance for language revitalization, cultural maintenance, and educational planning and assessment is discussed.…
Using a Simulation To Teach Food Web Dynamics.
ERIC Educational Resources Information Center
Rueter, John G.; Perrin, Nancy A.
1999-01-01
Reports on research that tested the effect of using a computer simulation to teach the concept of a food web to nonbiology majors in a large introductory course. Concludes that the use of the simulation resulted in significantly better performance on an open-ended essay question for those students who used the software, particularly for average…
NASA Astrophysics Data System (ADS)
Vivaldi, Franco
2015-12-01
The concept of resonance has been instrumental to the study of Hamiltonian systems with divided phase space. One can also define such systems over discrete spaces, which have a finite or countable number of points, but in this new setting the notion of resonance must be re-considered from scratch. I review some recent developments in the area of arithmetic dynamics which outline some salient features of linear and nonlinear stable (elliptic) orbits over a discrete space, and also underline the difficulties that emerge in their analysis.
NASA Astrophysics Data System (ADS)
Vivaldi, Franco
The concept of resonance has been instrumental to the study of Hamiltonian systems with divided phase space. One can also define such systems over discrete spaces, which have a finite or countable number of points, but in this new setting the notion of resonance must be re-considered from scratch. I review some recent developments in the area of arithmetic dynamics which outline some salient features of linear and nonlinear stable (elliptic) orbits over a discrete space, and also underline the difficulties that emerge in their analysis.
There must be a better way! Managing a corporate web site dynamically from a database
j.z. cohen
1998-10-21
This document is a set of slides available from http://www1.y12.org/lmes_sti/html/ycsdinf-98-8/index.htm that describes limitations of static web pages for conveying information, a plan for overcoming these limitations by generating web pages dynamically from a database, expected advantages and disadvantages of this method, design for a system using the method, and future plans.
Shorikov, A. F.
2014-11-18
We consider a discrete-time dynamical system consisting of three controllable objects. The motions of all objects are given by the corresponding vector linear or convex discrete-time recurrent vector relations, and control system for its has two levels: basic (first or I level) that is dominating and subordinate level (second or II level) and both have different criterions of functioning and united a priori by determined informational and control connections defined in advance. For the dynamical system in question, we propose a mathematical formalization in the form of solving a multistep problem of two-level hierarchical minimax program control over the terminal approach process with incomplete information and give a general scheme for its solution.
Martinez-Pedrero, Fernando; Tierno, Pietro; Johansen, Tom H.; Straube, Arthur V.
2016-01-01
The emergence of wave fronts in dissipative driven systems is a fascinating phenomenon which can be found in a broad range of physical and biological disciplines. Here we report the direct experimental observation of discrete fronts propagating along chains of paramagnetic colloidal particles, the latter propelled above a traveling wave potential generated by a structured magnetic substrate. We develop a rigorously reduced theoretical framework and describe the dynamics of the system in terms of a generalized one-dimensional dissipative Frenkel-Kontorova model. The front dynamics is explored in a wide range of field parameters close to and far from depinning, where the discrete and continuum limits apply. We show how symmetry breaking and finite size of chains are used to control the direction of front propagation, a universal feature relevant to different systems and important for real applications. PMID:26837286
NASA Astrophysics Data System (ADS)
Martinez-Pedrero, Fernando; Tierno, Pietro; Johansen, Tom H.; Straube, Arthur V.
2016-02-01
The emergence of wave fronts in dissipative driven systems is a fascinating phenomenon which can be found in a broad range of physical and biological disciplines. Here we report the direct experimental observation of discrete fronts propagating along chains of paramagnetic colloidal particles, the latter propelled above a traveling wave potential generated by a structured magnetic substrate. We develop a rigorously reduced theoretical framework and describe the dynamics of the system in terms of a generalized one-dimensional dissipative Frenkel-Kontorova model. The front dynamics is explored in a wide range of field parameters close to and far from depinning, where the discrete and continuum limits apply. We show how symmetry breaking and finite size of chains are used to control the direction of front propagation, a universal feature relevant to different systems and important for real applications.
Martinez-Pedrero, Fernando; Tierno, Pietro; Johansen, Tom H; Straube, Arthur V
2016-01-01
The emergence of wave fronts in dissipative driven systems is a fascinating phenomenon which can be found in a broad range of physical and biological disciplines. Here we report the direct experimental observation of discrete fronts propagating along chains of paramagnetic colloidal particles, the latter propelled above a traveling wave potential generated by a structured magnetic substrate. We develop a rigorously reduced theoretical framework and describe the dynamics of the system in terms of a generalized one-dimensional dissipative Frenkel-Kontorova model. The front dynamics is explored in a wide range of field parameters close to and far from depinning, where the discrete and continuum limits apply. We show how symmetry breaking and finite size of chains are used to control the direction of front propagation, a universal feature relevant to different systems and important for real applications. PMID:26837286
A Web Interface for the Quantification of Microtubule Dynamics
Kong, Koon Yin; Marcus, Adam I.; Giaanakakou, Paraskevi; Wang, May D.
2016-01-01
We propose a web interface that allows researchers to quantify and analyze microtubule confocal images online. Most analyses of microtubule confocal images are performed manually using very simple software or tools. Analysis results are stored locally within each collaborator with different styles and formats. This has limited the sharing of data and results when collaborating among different research parties. A web interface provides a simple way for users to process data online. It also allows easy sharing of both data and results among different participating groups. Analysis workflow of the interface is made similar to existing manual protocols. We demonstrate the integration of image processing algorithm in the current workflow to aid the analysis. Our design also allows integration of novel automated analysis algorithms and modules to re-evaluate existing data. This interface can provide a validation platform for new automated algorithm and allow collaboration on microtubule image analysis from different locations.
NASA Astrophysics Data System (ADS)
Zhu, Ling; van de Ven, Glenn; Watkins, Laura L.; Posti, Lorenzo
2016-08-01
We present a new discrete chemo-dynamical axisymmetric modeling technique, which we apply to the dwarf spheroidal galaxy Sculptor. The major improvement over previous Jeans models is that realistic chemical distributions are included directly in the dynamical modelling of the discrete data. This avoids loss of information due to spatial binning and eliminates the need for hard cuts to remove contaminants and to separate stars based on their chemical properties. Using a combined likelihood in position, metallicity and kinematics, we find that our models naturally separate Sculptor stars into a metal-rich and a metal-poor population. Allowing for non-spherical symmetry, our approach provides a central slope of the dark matter density of γ = 0.5 ± 0.3. The metal-rich population is nearly isotropic (with β _r^{red} = 0.0± 0.1) while the metal-poor population is tangentially anisotropic (with β _r^{blue} = -0.2± 0.1) around the half light radius of 0.26 kpc. A weak internal rotation of the metal-rich population is revealed with vmax/σ0 = 0.15 ± 0.15. We run tests using mock data to show that a discrete dataset with ˜6000 stars is required to distinguish between a core (γ = 0) and cusp (γ = 1), and to constrain the possible internal rotation to better than 1 σ confidence with our model. We conclude that our discrete chemo-dynamical modelling technique provides a flexible and powerful tool to robustly constrain the internal dynamics of multiple populations, and the total mass distribution in a stellar system.
A discrete momentum-conserving explicit algorithm for rigid body dynamics analysis
NASA Technical Reports Server (NTRS)
Park, K. C.; Chiou, J. C.
1993-01-01
A discrete momentum-conserving explicit time integration is presented. The accurate feature and simplicity of the present algorithm are realized by a mid-point implicit formula for integrating the Euler parameters and a second-order discrete momentum-conserving form of the central difference algorithm, respectively. The accuracy and robustness of the algorithm is demonstrated by example problems which exhibit large overall rigid motions under holonomic constraints.
2015-01-01
When simulating molecular systems using deterministic equations of motion (e.g., Newtonian dynamics), such equations are generally numerically integrated according to a well-developed set of algorithms that share commonly agreed-upon desirable properties. However, for stochastic equations of motion (e.g., Langevin dynamics), there is still broad disagreement over which integration algorithms are most appropriate. While multiple desiderata have been proposed throughout the literature, consensus on which criteria are important is absent, and no published integration scheme satisfies all desiderata simultaneously. Additional nontrivial complications stem from simulating systems driven out of equilibrium using existing stochastic integration schemes in conjunction with recently developed nonequilibrium fluctuation theorems. Here, we examine a family of discrete time integration schemes for Langevin dynamics, assessing how each member satisfies a variety of desiderata that have been enumerated in prior efforts to construct suitable Langevin integrators. We show that the incorporation of a novel time step rescaling in the deterministic updates of position and velocity can correct a number of dynamical defects in these integrators. Finally, we identify a particular splitting (related to the velocity Verlet discretization) that has essentially universally appropriate properties for the simulation of Langevin dynamics for molecular systems in equilibrium, nonequilibrium, and path sampling contexts. PMID:24555448
Stabilization of chaotic and non-permanent food-web dynamics
NASA Astrophysics Data System (ADS)
Williams, R. J.; Martinez, N. D.
2004-03-01
Several decades of dynamical analyses of food-web networks[CITE] have led to important insights into the effects of complexity, omnivory and interaction strength on food-web stability[CITE]. Several recent insights[CITE] are based on nonlinear bioenergetic consumer-resource models[CITE] that display chaotic behavior in three species food chains[CITE] which can be stabilized by omnivory[CITE] and weak interaction of a fourth species[CITE]. We slightly relax feeding on low-density prey in these models by modifying standard food-web interactions known as “typeII” functional responses[CITE]. This change drastically alters the dynamics of realistic systems containing up to ten species. Our modification stabilizes chaotic dynamics in three species systems and reduces or eliminates extinctions and non-persistent chaos[CITE] in ten species systems. This increased stability allows analysis of systems with greater biodiversity than in earlier work and suggests that dynamic stability is not as severe a constraint on the structure of large food webs as previously thought. The sensitivity of dynamical models to small changes in the predator-prey functional response well within the range of what is empirically observed suggests that functional response is a crucial aspect of species interactions that must be more precisely addressed in empirical studies.
Cheng, Cheng; Zhang, Xiaobing
2013-05-01
In conventional models for two-phase reactive flow of interior ballistic, the dynamic collision phenomenon of particles is neglected or empirically simplified. However, the particle collision between particles may play an important role in dilute two-phase flow because the distribution of particles is extremely nonuniform. The collision force may be one of the key factors to influence the particle movement. This paper presents the CFD-DEM approach for simulation of interior ballistic two-phase flow considering the dynamic collision process. The gas phase is treated as a Eulerian continuum and described by a computational fluid dynamic method (CFD). The solid phase is modeled by discrete element method (DEM) using a soft sphere approach for the particle collision dynamic. The model takes into account grain combustion, particle-particle collisions, particle-wall collisions, interphase drag and heat transfer between gas and solid phases. The continuous gas phase equations are discretized in finite volume form and solved by the AUSM+-up scheme with the higher order accurate reconstruction method. Translational and rotational motions of discrete particles are solved by explicit time integrations. The direct mapping contact detection algorithm is used. The multigrid method is applied in the void fraction calculation, the contact detection procedure, and CFD solving procedure. Several verification tests demonstrate the accuracy and reliability of this approach. The simulation of an experimental igniter device in open air shows good agreement between the model and experimental measurements. This paper has implications for improving the ability to capture the complex physics phenomena of two-phase flow during the interior ballistic cycle and to predict dynamic collision phenomena at the individual particle scale. PMID:24891728
Dynamics of the Lake Michigan food web, 1970-2000
Madenjian, Charles P.; Fahnenstiel, Gary L.; Johengen, Thomas H.; Nalepa, Thomas F.; Vanderploeg, Henry A.; Fleischer, Guy W.; Schneeberger, Philip J.; Benjamin, Darren M.; Smith, Emily B.; Bence, James R.; Rutherford, Edward S.; Lavis, Dennis S.; Robertson, Dale M.; Jude, David J.; Ebener, Mark P.
2002-01-01
Herein, we document changes in the Lake Michigan food web between 1970 and 2000 and identify the factors responsible for these changes. Control of sea lamprey (Petromyzon marinus) and alewife (Alosa pseudoharengus) populations in Lake Michigan, beginning in the 1950s and 1960s, had profound effects on the food web. Recoveries of lake whitefish (Coregonus clupeaformis) and burbot (Lota lota) populations, as well as the buildup of salmonine populations, were attributable, at least in part, to sea lamprey control. Based on our analyses, predation by salmonines was primarily responsible for the reduction in alewife abundance during the 1970s and early 1980s. In turn, the decrease in alewife abundance likely contributed to recoveries of deepwater sculpin (Myoxocephalus thompsoni), yellow perch (Perca flavescens), and burbot populations during the 1970s and 1980s. Decrease in the abundance of all three dominant benthic macroinvertebrate groups, including Diporeia, oligochaetes, and sphaeriids, during the 1980s in nearshore waters (50 m deep) of Lake Michigan, was attributable to a decrease in primary production linked to a decline in phosphorus loadings. Continued decrease in Diporeia abundance during the 1990s was associated with the zebra mussel (Dreissena polymorpha) invasion, but specific mechanisms for zebra mussels affecting Diporeia abundance remain unidentified.
Howard, Ian S; Ingram, James N; Wolpert, Daniel M
2011-04-01
Rhythmic and discrete arm movements occur ubiquitously in everyday life, and there is a debate as to whether these two classes of movements arise from the same or different underlying neural mechanisms. Here we examine interference in a motor-learning paradigm to test whether rhythmic and discrete movements employ at least partially separate neural representations. Subjects were required to make circular movements of their right hand while they were exposed to a velocity-dependent force field that perturbed the circularity of the movement path. The direction of the force-field perturbation reversed at the end of each block of 20 revolutions. When subjects made only rhythmic or only discrete circular movements, interference was observed when switching between the two opposing force fields. However, when subjects alternated between blocks of rhythmic and discrete movements, such that each was uniquely associated with one of the perturbation directions, interference was significantly reduced. Only in this case did subjects learn to corepresent the two opposing perturbations, suggesting that different neural resources were employed for the two movement types. Our results provide further evidence that rhythmic and discrete movements employ at least partially separate control mechanisms in the motor system. PMID:21273324
Systematic generation of nonlinear discretized dynamic equilibrium equations of spinning cantilevers
NASA Technical Reports Server (NTRS)
El-Essawi, M.; Utku, S.; Salama, M.
1982-01-01
General nonlinear discretized governing equations of motion of spinning elastic solids and structures are adjusted for the case of a spinning cantilever with initial geometric imperfections. Consideration is given to second degree nonlinearities in the strain-displacement and velocity-displacement relationships. Parameters of the discretization are developed to include the type and number of the coordinate functions used in the admissible trial solution in order to unify the discretization approaches associated with stationarity principles. The coordinate functions comprise both sets of continuous and piecewise continuous functions employed in the Rayleigh-Ritz and the finite element methods, respectively. Coefficient matrices are provided which contain the energy density expressions and which are adaptable to computer programming.
NASA Technical Reports Server (NTRS)
Housner, J. M.; Mcgowan, P. E.
1986-01-01
A computerized procedure is presented for calculating instantaneous velocity changes due to discrete imposition of constraints, impact or discrete inertial changes. Such discrete changes are impulsive in nature, involve wave propagation and therefore usually require detailed modeling for accurate prediction. It is demonstrated that when modeling is coarse, physically unrealistic results can be produced. The present procedure, which concentrates on finite element analysis, reduces the need for detailed modeling by assuming that only those velocities in the neighborhood of grid points involved in the constraint or those rigidly connected to grid points involved in the constraint need be considered in calculating velocity changes. The procedure utilizes Newton's second law and as such conserves momentum where valid. Sample problems are provided.
From a discrete to a continuum model of cell dynamics in one dimension.
Murray, Philip J; Edwards, Carina M; Tindall, Marcus J; Maini, Philip K
2009-09-01
Multiscale modeling is emerging as one of the key challenges in mathematical biology. However, the recent rapid increase in the number of modeling methodologies being used to describe cell populations has raised a number of interesting questions. For example, at the cellular scale, how can the appropriate discrete cell-level model be identified in a given context? Additionally, how can the many phenomenological assumptions used in the derivation of models at the continuum scale be related to individual cell behavior? In order to begin to address such questions, we consider a discrete one-dimensional cell-based model in which cells are assumed to interact via linear springs. From the discrete equations of motion, the continuous Rouse [P. E. Rouse, J. Chem. Phys. 21, 1272 (1953)] model is obtained. This formalism readily allows the definition of a cell number density for which a nonlinear "fast" diffusion equation is derived. Excellent agreement is demonstrated between the continuum and discrete models. Subsequently, via the incorporation of cell division, we demonstrate that the derived nonlinear diffusion model is robust to the inclusion of more realistic biological detail. In the limit of stiff springs, where cells can be considered to be incompressible, we show that cell velocity can be directly related to cell production. This assumption is frequently made in the literature but our derivation places limits on its validity. Finally, the model is compared with a model of a similar form recently derived for a different discrete cell-based model and it is shown how the different diffusion coefficients can be understood in terms of the underlying assumptions about cell behavior in the respective discrete models. PMID:19905151
Technology Transfer Automated Retrieval System (TEKTRAN)
A linear discrete dynamic system model is constructed to represent the temporal interactions among significantly expressed genes in response to bioethanol conversion inhibitor 5-hydroxymethylfurfural for ethanologenic yeast Saccharomyces cerevisiae. This study identifies the most significant linear...
NASA Astrophysics Data System (ADS)
Hancock, W.; Weatherley, D.; Wruck, B.; Chitombo, G. P.
2012-04-01
The flow dynamics of granular materials is of broad interest in both the geosciences (e.g. landslides, fault zone evolution, and brecchia pipe formation) and many engineering disciplines (e.g chemical engineering, food sciences, pharmaceuticals and materials science). At the interface between natural and human-induced granular media flow, current underground mass-mining methods are trending towards the induced failure and subsequent gravitational flow of large volumes of broken rock, a method known as cave mining. Cave mining relies upon the undercutting of a large ore body, inducement of fragmentation of the rock and subsequent extraction of ore from below, via hopper-like outlets. Design of such mines currently relies upon a simplified kinematic theory of granular flow in hoppers, known as the ellipsoid theory of mass movement. This theory assumes that the zone of moving material grows as an ellipsoid above the outlet of the silo. The boundary of the movement zone is a shear band and internal to the movement zone, the granular material is assumed to have a uniformly high bulk porosity compared with surrounding stagnant regions. There is however, increasing anecdotal evidence and field measurements suggesting this theory fails to capture the full complexity of granular material flow within cave mines. Given the practical challenges obstructing direct measurement of movement both in laboratory experiments and in-situ, the Discrete Element Method (DEM [1]) is a popular alternative to investigate granular media flow. Small-scale DEM studies (c.f. [3] and references therein) have confirmed that movement within DEM silo flow models matches that predicted by ellipsoid theory, at least for mono-disperse granular material freely outflowing at a constant rate. A major draw-back of these small-scale DEM studies is that the initial bulk porosity of the simulated granular material is significantly higher than that of broken, prismatic rock. In this investigation, more
Prototype of a Mobile Social Network for Education Using Dynamic Web Service
ERIC Educational Resources Information Center
Hoentsch, Sandra Costa Pinto; Carvalho, Felipe Oliveira; Santos, Luiz Marcus Monteiro de Almeida; Ribeiro, Admilson de Ribamar Lima
2012-01-01
This article presents the proposal of a social network site SocialNetLab that belongs to the Department of Computing-Federal University of Sergipe and which aims to locate and notify users of a nearby friend independently of the location technology available in the equipment through dynamic Web Service; to serve as a laboratory for research in…
Sociometry: An Approach for Assessing Group Dynamics in Web-Based Courses
ERIC Educational Resources Information Center
Daugherty, Martha; Turner, Jeff
2003-01-01
Student interactivity in web-based educational environments has shown to increase academic learning and motivation (Jiang, 1998; Petraglia, 1998). However, instructors often find it difficult to assess the quality of online group dynamics without visual observations of student behaviors. The purpose of this study was to investigate the use of…
ERIC Educational Resources Information Center
Flaherty, Brian P.
2008-01-01
Developmental research often involves studying change across 2 or more processes or constructs simultaneously. A natural question in this work is whether change in these 2 processes is related or independent. Associative latent transition analysis (ALTA) was designed to test hypotheses about the degree to which change in 2 discrete latent…
Regime shifts in marine communities: a complex systems perspective on food web dynamics
Yletyinen, Johanna; Bodin, Örjan; Weigel, Benjamin; Nordström, Marie C.; Bonsdorff, Erik; Blenckner, Thorsten
2016-01-01
Species composition and habitats are changing at unprecedented rates in the world's oceans, potentially causing entire food webs to shift to structurally and functionally different regimes. Despite the severity of these regime shifts, elucidating the precise nature of their underlying processes has remained difficult. We address this challenge with a new analytic approach to detect and assess the relative strength of different driving processes in food webs. Our study draws on complexity theory, and integrates the network-centric exponential random graph modelling (ERGM) framework developed within the social sciences with community ecology. In contrast to previous research, this approach makes clear assumptions of direction of causality and accommodates a dynamic perspective on the emergence of food webs. We apply our approach to analysing food webs of the Baltic Sea before and after a previously reported regime shift. Our results show that the dominant food web processes have remained largely the same, although we detect changes in their magnitudes. The results indicate that the reported regime shift may not be a system-wide shift, but instead involve a limited number of species. Our study emphasizes the importance of community-wide analysis on marine regime shifts and introduces a novel approach to examine food webs. PMID:26888032
Regime shifts in marine communities: a complex systems perspective on food web dynamics.
Yletyinen, Johanna; Bodin, Örjan; Weigel, Benjamin; Nordström, Marie C; Bonsdorff, Erik; Blenckner, Thorsten
2016-02-24
Species composition and habitats are changing at unprecedented rates in the world's oceans, potentially causing entire food webs to shift to structurally and functionally different regimes. Despite the severity of these regime shifts, elucidating the precise nature of their underlying processes has remained difficult. We address this challenge with a new analytic approach to detect and assess the relative strength of different driving processes in food webs. Our study draws on complexity theory, and integrates the network-centric exponential random graph modelling (ERGM) framework developed within the social sciences with community ecology. In contrast to previous research, this approach makes clear assumptions of direction of causality and accommodates a dynamic perspective on the emergence of food webs. We apply our approach to analysing food webs of the Baltic Sea before and after a previously reported regime shift. Our results show that the dominant food web processes have remained largely the same, although we detect changes in their magnitudes. The results indicate that the reported regime shift may not be a system-wide shift, but instead involve a limited number of species. Our study emphasizes the importance of community-wide analysis on marine regime shifts and introduces a novel approach to examine food webs. PMID:26888032
Scaling dynamics of deep orbits in a periodic stochastic web
NASA Astrophysics Data System (ADS)
Lowenstein, J. H.
1991-04-01
In a fourfold-symmetric stochastic web generated by a two-dimensional area-preserving map, those phase-space orbits that remain forever in the central part of the chaotic channels are selected by a suitable cutoff condition. By means of a horseshoe construction, these deep orbits are classified according to (i) their paths on the square lattice of hyperbolic fixed points, (ii) their transit times between successive lattice sites, and (iii) their horseshoe branch (left or ) at each lattice site. The set of initial conditions for the deep orbits can be placed on a tree that branches with infinite multiplicity at each level: level N of the tree corresponds to paths of N+1 steps on the lattice. Empirically, one finds two scaling relations, one for the level number and the other for the number of map iterations between successive saddle points. Using a transfer-matrix technique, the scaling relations are exploited to calculate the multifractal properties [Hausdorff dimension, Hausdorff measure, the f(α) curve] of the set of deep orbits, to obtain an invariant probability measure and to derive a Green function for the diffusion process on the lattice. Numerical calculations of many of these quantities are carried out for the value 0.5 of the control parameter.
NASA Astrophysics Data System (ADS)
Chang, K.-J.; Chen, R.-F.; Chan, Y.-C.; Kuo, C.-Y.; Weng, C.-H.
2012-04-01
transported debris was about 30 - 90 m thick, covered on the preexisting debris deposit hill and around the river channel. The debris formed a dammed lake, with a maximum volume of 45 Mm3. Based on DTMs data sets, field observations, the discrete element method - PFC3D is adapts to analyze the triggering mechanism and sling dynamic process. The presence and the coupling effect from the strong ground excitation and high pore water pressure is the essential factor to triggering the landslide event. The results shows that the best fit between the deposit topography of the post-event DTM and numerical simulations, the fictional coefficient of the sliding surface is as low as 0.087. The maximum sliding speed is as high as 87.2 m/s, the result coincide with the seismic record from the nearby strong motion seismic record.
The effects of host-feeding on stability of discrete-time host-parasitoid population dynamic models.
Emerick, Brooks; Singh, Abhyudai
2016-02-01
Discrete-time models are the traditional approach for capturing population dynamics of a host-parasitoid system. Recent work has introduced a semi-discrete framework for obtaining model update functions that connect host-parasitoid population levels from year-to-year. In particular, this framework uses differential equations to describe the host-parasitoid interaction during the time of year when they come in contact, allowing specific behaviors to be mechanistically incorporated. We use the semi-discrete approach to study the effects of host-feeding, which occurs when a parasitoid consumes a potential host larva without ovipositing. We find that host-feeding by itself cannot stabilize the system, and both populations exhibit behavior similar to the Nicholson-Bailey model. However, when combined with stabilizing mechanisms such as density-dependent host mortality, host-feeding contracts the region of parameter space that allows for a stable host-parasitoid equilibrium. In contrast, when combined with a density-dependent parasitoid attack rate, host-feeding expands the non-zero equilibrium stability region. Our results show that host-feeding causes inefficiency in the parasitoid population, which yields a higher population of hosts per generation. This suggests that host-feeding may have limited long-term impact in terms of suppressing host levels for biological control applications. PMID:26686008
Discrete simulation of the dynamics of spread of extreme opinions in a society
NASA Astrophysics Data System (ADS)
Stauffer, Dietrich; Sahimi, Muhammad
2006-05-01
We propose a discrete model for how opinions about a given “extreme” subject, about which various groups of a population have different degrees of enthusiasm for or susceptibility to, such as fanaticism, extreme social and political positions, and terrorism, may spread. The model, in a certain limit, is the discrete analogue of a deterministic continuum model suggested by others. We carry out extensive computer simulation of the model by utilizing it on lattices with infinite- or short-range interactions, and on symmetric and hierarchical (or directed) Barabási-Albert scale-free networks. Several interesting features of the model are demonstrated, and comparison is made with the deterministic continuum model.
NASA Technical Reports Server (NTRS)
Jameson, Antony
1994-01-01
The effect of artificial diffusion on discrete shock structures is examined for a family of schemes which includes scalar diffusion, convective upwind and split pressure (CUSP) schemes, and upwind schemes with characteristics splitting. The analysis leads to conditions on the diffusive flux such that stationary discrete shocks can contain a single interior point. The simplest formulation which meets these conditions is a CUSP scheme in which the coefficients of the pressure differences is fully determined by the coefficient of convective diffusion. It is also shown how both the characteristic and CUSP schemes can be modified to preserve constant stagnation enthalpy in steady flow, leading to four variants, the E and H-characteristic schemes, and the E and H-CUSP schemes. Numerical results are presented which confirm the properties of these schemes.
Biological vs. physical mixing effects on benthic food web dynamics.
Braeckman, Ulrike; Provoost, Pieter; Moens, Tom; Soetaert, Karline; Middelburg, Jack J; Vincx, Magda; Vanaverbeke, Jan
2011-01-01
Biological particle mixing (bioturbation) and solute transfer (bio-irrigation) contribute extensively to ecosystem functioning in sediments where physical mixing is low. Macrobenthos transports oxygen and organic matter deeper into the sediment, thereby likely providing favourable niches to lower trophic levels (i.e., smaller benthic animals such as meiofauna and bacteria) and thus stimulating mineralisation. Whether this biological transport facilitates fresh organic matter assimilation by the metazoan lower part of the food web through niche establishment (i.e., ecosystem engineering) or rather deprives them from food sources, is so far unclear. We investigated the effects of the ecosystem engineers Lanice conchilega (bio-irrigator) and Abra alba (bioturbator) compared to abiotic physical mixing events on survival and food uptake of nematodes after a simulated phytoplankton bloom. The (13)C labelled diatom Skeletonema costatum was added to 4 treatments: (1) microcosms containing the bioturbator, (2) microcosms containing the bio-irrigator, (3) control microcosms and (4) microcosms with abiotic manual surface mixing. Nematode survival and subsurface peaks in nematode density profiles were most pronounced in the bio-irrigator treatment. However, nematode specific uptake (Δδ(13)C) of the added diatoms was highest in the physical mixing treatment, where macrobenthos was absent and the diatom (13)C was homogenised. Overall, nematodes fed preferentially on bulk sedimentary organic material rather than the added diatoms. The total C budget (µg C m(-2)), which included TO(13)C remaining in the sediment, respiration, nematode and macrobenthic uptake, highlighted the limited assimilation by the metazoan benthos and the major role of bacterial respiration. In summary, bioturbation and especially bio-irrigation facilitated the lower trophic levels mainly over the long-term through niche establishment. Since the freshly added diatoms represented only a limited food source
Networked Physics Curriculum:. From Static Web to Dynamic Java
NASA Astrophysics Data System (ADS)
Bothun, G. D.; Kevan, S. D.; Micklavzina, S.; Mason, D.
We describe our efforts at the University of Oregon to use Web-based Instructional Technology (IT) supplemented with interactive Java virtual experiments to change the standard pedagogy associated with large, introductory undergraduate classes in physics and astronomy. We begin by examining some of the problems associated with the standard pedagogy in these classes and how these problems motivated our development of networked courseware. Although we identify and describe five empirical positive outcomes associated with IT, we conclude that the use of HTML-based course material and assignments does not substantially alter the standard pedagogy as this medium alone is not conducive to interactive exercises. To build interactivity into our courseware, we have undertaken a vigorous effort of creating Java-based experiments which are grounded in physical reality and duplicate the kinds of experiments that are done in the physical lab. In so doing, we build experimentation into a curriculum for large lecture-based classes in which the standard pedagogy and resource constraints normally preclude lab sections. The main goal is to create a networked environment where the student can easily retrieve the notes and the demonstrations that were done in class as well as to engage in experiments that are designed to illustrate basic principles. In so doing, we hope to move to a more learner-centered environment which is driven by student inquiry. Five specific Java experiments are described here and each is accompanied by a snapshot of the experimental apparatus and controls. An appendix contains the relevant URLs of the experiments, courseware, and animation described herein.
Nutrients, food web and tropho-dynamic processes
NASA Astrophysics Data System (ADS)
Wilson, James G.; Devlin, Michelle
2013-12-01
This volume presents the papers on the theme of nutrients and tropho-dynamic processes from the ECSA50 Venice 2012: Today's Science for Tomorrow's Management meeting. They illustrate how the application of cutting-edge techniques such as stable isotope analysis (SIA) together with improved analytical and synthetic procedures allow a much finer degree of definition in understanding system function. To this may be added increasingly sophisticated models, underpinned by ever more precise and detailed data, through which we can predict with greater and greater certainty the consequences of the changes Drivers, Pressures, States and Impacts on the system to shape the management Responses.
NASA Astrophysics Data System (ADS)
Piccardi, Carlo; Soncini-Sessa, Rodolfo
1991-05-01
The solution via dynamic programming (DP) of a reservoir optimal control problem is often computationally prohibitive when the proper description of the inflow process leads to a system model having several state variables and/or when a sufficiently dense state discretization is required to achieve numerical accuracy. Thus, to simplify, the inflow correlation is usually neglected and/or a coarse state discretization is adopted. However, these simplifications may significantly affect the reliability of the solution of the optimization problem. Nowadays, the availability of very powerful computers based on innovative architectures (vector and parallel machines), even in the domain of personal computers (transputer architectures), stimulates the reformulation of the standard dynamic programming algorithm in a form able to exploit these new machine architectures. The reformulated DP algorithm and new machines enable faster and less costly solution of optimization problems involving a system model having two state variables (storage and previous period inflow, then taking into account the inflow correlation) and a number of states (of the order of 104) such as to guarantee a high numerical accuracy.
NASA Astrophysics Data System (ADS)
Kengne, E.; Lakhssassi, A.
2015-03-01
We consider a lossless one-dimensional nonlinear discrete bi-inductance electrical transmission line made of N identical unit cells. When lattice effects are considered, we use the reductive perturbation method in the semidiscrete limit to show that the dynamics of modulated waves can be modeled by the classical nonlinear Schrödinger (CNLS) equation, which describes the modulational instability and the propagation of bright and dark solitons on a continuous-wave background. Our theoretical analysis based on the CNLS equation predicts either two or four frequency regions with different behavior concerning the modulational instability of a plane wave. With the help of the analytical solutions of the CNLS equation, we investigate analytically the effects of the linear capacitance CS on the dynamics of matter-wave solitons in the network. Our results reveal that the linear parameter CS can be used to manipulate the motion of bright, dark, and kink soliton in the network.
Kengne, E; Lakhssassi, A
2015-03-01
We consider a lossless one-dimensional nonlinear discrete bi-inductance electrical transmission line made of N identical unit cells. When lattice effects are considered, we use the reductive perturbation method in the semidiscrete limit to show that the dynamics of modulated waves can be modeled by the classical nonlinear Schrödinger (CNLS) equation, which describes the modulational instability and the propagation of bright and dark solitons on a continuous-wave background. Our theoretical analysis based on the CNLS equation predicts either two or four frequency regions with different behavior concerning the modulational instability of a plane wave. With the help of the analytical solutions of the CNLS equation, we investigate analytically the effects of the linear capacitance CS on the dynamics of matter-wave solitons in the network. Our results reveal that the linear parameter CS can be used to manipulate the motion of bright, dark, and kink soliton in the network. PMID:25871172
Avalanche dynamics of magnetic flux in a two-dimensional discrete superconductor
Ginzburg, S. L.; Nakin, A. V.; Savitskaya, N. E.
2006-11-15
The critical state of a two-dimensional discrete superconductor in an external magnetic field is studied. This state is found to be self-organized in the generalized sense, i.e., is a set of metastable states that transform to each other by means of avalanches. An avalanche is characterized by the penetration of a magnetic flux to the system. The sizes of the occurring avalanches, i.e., changes in the magnetic flux, exhibit the power-law distribution. It is also shown that the size of the avalanche occurring in the critical state and the external magnetic field causing its change are statistically independent quantities.
Shaitan, K V; Orshanskiy, I A
2015-01-01
In this study we suggested a dynamics simulation for the formation of protofiber of spider web nanofiber. It was shown that a bundle of parallel polyalanine β-strands of sufficient length is arranged through self-assembly into a stable right-handed super helix. By numerical analysis we investigated the rheological properties and provided in nonlinear regime a generalization of the model of Singer for description of the rheological behaviour of super helix. PMID:26394464
Kooistra, Lammert; Bergsma, Aldo; Chuma, Beatus; de Bruin, Sytze
2009-01-01
This paper describes the development of a sensor web based approach which combines earth observation and in situ sensor data to derive typical information offered by a dynamic web mapping service (WMS). A prototype has been developed which provides daily maps of vegetation productivity for the Netherlands with a spatial resolution of 250 m. Daily available MODIS surface reflectance products and meteorological parameters obtained through a Sensor Observation Service (SOS) were used as input for a vegetation productivity model. This paper presents the vegetation productivity model, the sensor data sources and the implementation of the automated processing facility. Finally, an evaluation is made of the opportunities and limitations of sensor web based approaches for the development of web services which combine both satellite and in situ sensor sources. PMID:22574019
Feng, Rui; Xenos, Michalis; Girdhar, Gaurav; Kang, Wei; Davenport, James W; Deng, Yuefan; Bluestein, Danny
2012-01-01
Flow and stresses induced by blood flow acting on the blood cellular constituents can be represented to a certain extent by a continuum mechanics approach down to the order of the μm level. However, the molecular effects of, e.g., adhesion/aggregation bonds of blood clotting can be on the order of nm. The coupling of the disparate length and timescales between such molecular levels and macroscopic transport represents a major computational challenge. To address this challenge, a multiscale numerical approach based on discrete particle dynamics (DPD) methodology derived from molecular dynamics (MD) principles is proposed. The feasibility of the approach was firstly tested for its ability to simulate viscous flow conditions. Simulations were conducted in low Reynolds numbers flows (Re = 25-33) through constricted tubes representing blood vessels with various degrees of stenosis. Multiple discrete particles interacting with each other were simulated, with 1.24-1.36 million particles representing the flow domain and 0.4 million particles representing the vessel wall. The computation was carried out on the massive parallel supercomputer NY BlueGene/L employing NAMD-a parallel MD package for high performance computing (HPC). Typical recirculation zones were formed distal to the stenoses. The velocity profiles and recirculation zones were in excellent agreement with computational fluid dynamics (CFD) 3D Navier-Stokes viscous fluid flow simulations and with classic numerical and experimental results by YC Fung in constricted tubes. This feasibility analysis demonstrates the potential of a methodology that widely departs from a continuum approach to simulate multiscale phenomena such as flow induced blood clotting. PMID:21369918
The multi-stream flows and the dynamics of the cosmic web
Shandarin, Sergei F.
2011-05-01
A new numerical technique to identify the cosmic web is proposed. It is based on locating multi-stream flows, i.e. the places where the velocity field is multi-valued. The method is local in Eulerian space, simple and computationally efficient. This technique uses the velocities of particles and thus takes into account the dynamical information. This is in contrast with the majority of standard methods that use the coordinates of particles only. Two quantities are computed in every mesh cell: the mean and variance of the velocity field. Ideally in the cells where the velocity is single-valued the variance must be equal to zero exactly, therefore the cells with non-zero variance are identified as multi-stream flows. The technique has been tested in the Zel'dovich approximation and in the N-body simulation of the ΛCDM model. The effect of numerical noise is discussed. The web identified by the new method has been compared with the web identified by the standard technique using only the particle coordinates. The comparison has shown overall similarity of two webs as expected, however they by no means are identical. For example, the isocontours of the corresponding fields have significantly different shapes and some density peaks of similar heights exhibit significant differences in the velocity variance and vice versa. This suggests that the density and velocity variance have a significant degree of independence. The shape of the two-dimensional pdf of density and velocity variance confirms this proposition. Thus, we conclude that the dynamical information probed by this technique introduces an additional dimension into analysis of the web.
Dynamic web cache publishing for IaaS clouds using Shoal
NASA Astrophysics Data System (ADS)
Gable, Ian; Chester, Michael; Armstrong, Patrick; Berghaus, Frank; Charbonneau, Andre; Leavett-Brown, Colin; Paterson, Michael; Prior, Robert; Sobie, Randall; Taylor, Ryan
2014-06-01
We have developed a highly scalable application, called Shoal, for tracking and utilizing a distributed set of HTTP web caches. Our application uses the Squid HTTP cache. Squid servers advertise their existence to the Shoal server via AMQP messaging by running Shoal Agent. The Shoal server provides a simple REST interface that allows clients to determine their closest Squid cache. Our goal is to dynamically instantiate Squid caches on IaaS clouds in response to client demand. Shoal provides the VMs on IaaS clouds with the location of the nearest dynamically instantiated Squid Cache.
2015-01-01
By utilizing Graphics Processing Units, we show that constant pH molecular dynamics simulations (CpHMD) run in Generalized Born (GB) implicit solvent for long time scales can yield poor pKa predictions as a result of sampling unrealistic conformations. To address this shortcoming, we present a method for performing constant pH molecular dynamics simulations (CpHMD) in explicit solvent using a discrete protonation state model. The method involves standard molecular dynamics (MD) being propagated in explicit solvent followed by protonation state changes being attempted in GB implicit solvent at fixed intervals. Replica exchange along the pH-dimension (pH-REMD) helps to obtain acceptable titration behavior with the proposed method. We analyzed the effects of various parameters and settings on the titration behavior of CpHMD and pH-REMD in explicit solvent, including the size of the simulation unit cell and the length of the relaxation dynamics following protonation state changes. We tested the method with the amino acid model compounds, a small pentapeptide with two titratable sites, and hen egg white lysozyme (HEWL). The proposed method yields superior predicted pKa values for HEWL over hundreds of nanoseconds of simulation relative to corresponding predicted values from simulations run in implicit solvent. PMID:24803862
Attractors of relaxation discrete-time systems with chaotic dynamics on a fast time scale
NASA Astrophysics Data System (ADS)
Maslennikov, Oleg V.; Nekorkin, Vladimir I.
2016-07-01
In this work, a new type of relaxation systems is considered. Their prominent feature is that they comprise two distinct epochs, one is slow regular motion and another is fast chaotic motion. Unlike traditionally studied slow-fast systems that have smooth manifolds of slow motions in the phase space and fast trajectories between them, in this new type one observes, apart the same geometric objects, areas of transient chaos. Alternating periods of slow regular motions and fast chaotic ones as well as transitions between them result in a specific chaotic attractor with chaos on a fast time scale. We formulate basic properties of such attractors in the framework of discrete-time systems and consider several examples. Finally, we provide an important application of such systems, the neuronal electrical activity in the form of chaotic spike-burst oscillations.
Discrete approach to stochastic parametrization and dimension reduction in nonlinear dynamics.
Chorin, Alexandre J; Lu, Fei
2015-08-11
Many physical systems are described by nonlinear differential equations that are too complicated to solve in full. A natural way to proceed is to divide the variables into those that are of direct interest and those that are not, formulate solvable approximate equations for the variables of greater interest, and use data and statistical methods to account for the impact of the other variables. In the present paper we consider time-dependent problems and introduce a fully discrete solution method, which simplifies both the analysis of the data and the numerical algorithms. The resulting time series are identified by a NARMAX (nonlinear autoregression moving average with exogenous input) representation familiar from engineering practice. The connections with the Mori-Zwanzig formalism of statistical physics are discussed, as well as an application to the Lorenz 96 system. PMID:26216975
The global dynamics of a discrete juvenile-adult model with continuous and seasonal reproduction.
Ackleh, Azmy S; Chiquet, Ross A
2009-03-01
A general discrete juvenile-adult population model with time-dependent birth rate and nonlinear survivorship rates is studied. When breeding is continuous, it is shown that the model has a unique globally asymptotically stable positive equilibrium provided the net reproductive number is larger than one. If it is smaller than one, then the extinction equilibrium is globally asymptotically stable. When breeding is seasonal, it is shown that there exists a unique globally asymptotically stable periodic solution provided the net reproductive number is larger than one. When this value is less than one, the population goes to extinction. Conditions on the birth rate where the population with seasonal breeding survives while the population with continuous breeding becomes extinct are provided. PMID:22880823
Discrete approach to stochastic parametrization and dimension reduction in nonlinear dynamics
Chorin, Alexandre J.; Lu, Fei
2015-01-01
Many physical systems are described by nonlinear differential equations that are too complicated to solve in full. A natural way to proceed is to divide the variables into those that are of direct interest and those that are not, formulate solvable approximate equations for the variables of greater interest, and use data and statistical methods to account for the impact of the other variables. In the present paper we consider time-dependent problems and introduce a fully discrete solution method, which simplifies both the analysis of the data and the numerical algorithms. The resulting time series are identified by a NARMAX (nonlinear autoregression moving average with exogenous input) representation familiar from engineering practice. The connections with the Mori–Zwanzig formalism of statistical physics are discussed, as well as an application to the Lorenz 96 system. PMID:26216975
Invited Article: VEDA: A web-based virtual environment for dynamic atomic force microscopy
NASA Astrophysics Data System (ADS)
Melcher, John; Hu, Shuiqing; Raman, Arvind
2008-06-01
We describe here the theory and applications of virtual environment dynamic atomic force microscopy (VEDA), a suite of state-of-the-art simulation tools deployed on nanoHUB (www.nanohub.org) for the accurate simulation of tip motion in dynamic atomic force microscopy (dAFM) over organic and inorganic samples. VEDA takes advantage of nanoHUB's cyberinfrastructure to run high-fidelity dAFM tip dynamics computations on local clusters and the teragrid. Consequently, these tools are freely accessible and the dAFM simulations are run using standard web-based browsers without requiring additional software. A wide range of issues in dAFM ranging from optimal probe choice, probe stability, and tip-sample interaction forces, power dissipation, to material property extraction and scanning dynamics over hetereogeneous samples can be addressed.
A new dimension: Evolutionary food web dynamics in two dimensional trait space.
Ritterskamp, Daniel; Bearup, Daniel; Blasius, Bernd
2016-09-21
Species within a habitat are not uniformly distributed. However this aspect of community structure, which is fundamental to many conservation activities, is neglected in the majority of models of food web assembly. To address this issue, we introduce a model which incorporates a second dimension, which can be interpreted as space, into the trait space used in evolutionary food web models. Our results show that the additional trait axis allows the emergence of communities with a much greater range of network structures, similar to the diversity observed in real ecological communities. Moreover, the network properties of the food webs obtained are in good agreement with those of empirical food webs. Community emergence follows a consistent pattern with spread along the second trait axis occurring before the assembly of higher trophic levels. Communities can reach either a static final structure, or constantly evolve. We observe that the relative importance of competition and predation is a key determinant of the network structure and the evolutionary dynamics. The latter are driven by the interaction-competition and predation-between small groups of species. The model remains sufficiently simple that we are able to identify the factors, and mechanisms, which determine the final community state. PMID:27060671
Wang, Fei-Yue; Jin, Ning; Liu, Derong; Wei, Qinglai
2011-01-01
In this paper, we study the finite-horizon optimal control problem for discrete-time nonlinear systems using the adaptive dynamic programming (ADP) approach. The idea is to use an iterative ADP algorithm to obtain the optimal control law which makes the performance index function close to the greatest lower bound of all performance indices within an ε-error bound. The optimal number of control steps can also be obtained by the proposed ADP algorithms. A convergence analysis of the proposed ADP algorithms in terms of performance index function and control policy is made. In order to facilitate the implementation of the iterative ADP algorithms, neural networks are used for approximating the performance index function, computing the optimal control policy, and modeling the nonlinear system. Finally, two simulation examples are employed to illustrate the applicability of the proposed method. PMID:20876014
NASA Astrophysics Data System (ADS)
Vattré, A.; Devincre, B.; Feyel, F.; Gatti, R.; Groh, S.; Jamond, O.; Roos, A.
2014-02-01
A unified model coupling 3D dislocation dynamics (DD) simulations with the finite element (FE) method is revisited. The so-called Discrete-Continuous Model (DCM) aims to predict plastic flow at the (sub-)micron length scale of materials with complex boundary conditions. The evolution of the dislocation microstructure and the short-range dislocation-dislocation interactions are calculated with a DD code. The long-range mechanical fields due to the dislocations are calculated by a FE code, taking into account the boundary conditions. The coupling procedure is based on eigenstrain theory, and the precise manner in which the plastic slip, i.e. the dislocation glide as calculated by the DD code, is transferred to the integration points of the FE mesh is described in full detail. Several test cases are presented, and the DCM is applied to plastic flow in a single-crystal Nickel-based superalloy.
Existence, stability and dynamics of discrete solitary waves in a binary waveguide array
NASA Astrophysics Data System (ADS)
Shen, Y.; Kevrekidis, P. G.; Srinivasan, G.; Aceves, A. B.
2016-07-01
Recent work has explored binary waveguide arrays in the long-wavelength, near-continuum limit, here we examine the opposite limit, namely the vicinity of the so-called anti-continuum limit. We provide a systematic discussion of states involving one, two and three excited waveguides, and provide comparisons that illustrate how the stability of these states differ from the monoatomic limit of a single type of waveguide. We do so by developing a general theory which systematically tracks down the key eigenvalues of the linearized system. When we find the states to be unstable, we explore their dynamical evolution through direct numerical simulations. The latter typically illustrate, for the parameter values considered herein, the persistence of localized dynamics and the emergence for the duration of our simulations of robust quasi-periodic states for two excited sites. As the number of excited nodes increases, the unstable dynamics feature less regular oscillations of the solution’s amplitude.
Markov state modeling and dynamical coarse-graining via discrete relaxation path sampling.
Fačkovec, B; Vanden-Eijnden, E; Wales, D J
2015-07-28
A method is derived to coarse-grain the dynamics of complex molecular systems to a Markov jump process (MJP) describing how the system jumps between cells that fully partition its state space. The main inputs are relaxation times for each pair of cells, which are shown to be robust with respect to positioning of the cell boundaries. These relaxation times can be calculated via molecular dynamics simulations performed in each cell separately and are used in an efficient estimator for the rate matrix of the MJP. The method is illustrated through applications to Sinai billiards and a cluster of Lennard-Jones discs. PMID:26233119
Data compression of discrete sequence: A tree based approach using dynamic programming
NASA Technical Reports Server (NTRS)
Shivaram, Gurusrasad; Seetharaman, Guna; Rao, T. R. N.
1994-01-01
A dynamic programming based approach for data compression of a ID sequence is presented. The compression of an input sequence of size N to that of a smaller size k is achieved by dividing the input sequence into k subsequences and replacing the subsequences by their respective average values. The partitioning of the input sequence is carried with the intention of reducing the mean squared error in the reconstructed sequence. The complexity involved in finding the partitions which would result in such an optimal compressed sequence is reduced by using the dynamic programming approach, which is presented.
Polyakov, Pavel D; Duval, Jérôme F L
2014-02-01
We report a comprehensive theory to evaluate the kinetics of complex formation between metal ions and charged spherical nanoparticles. The latter consist of an ion-impermeable core surrounded by a soft shell layer characterized by a discrete axisymmetric 2D distribution of charged sites that bind metal ions. The theory explicitly integrates the conductive diffusion of metal ions from bulk solution toward the respective locations of the reactive sites within the particle shell volume. The kinetic constant k for outer-sphere nanoparticle-metal association is obtained from the sum of the contributions stemming from all reactive sites, each evaluated from the corresponding incoming flux of metal ions derived from steady-state Poisson-Nernst-Planck equations. Illustrations are provided to capture the basic intertwined impacts of particle size, overall particle charge, spatial heterogeneity in site distribution, type of particle (hard, core-shell or porous) and concentration of the background electrolyte on k. As a limit, k converges with predictions from previously reported analytical expressions derived for porous particles with low and high charge density, cases that correspond to coulombic and mean-field (smeared-out) electrostatic treatments, respectively. The conditions underlying the applicability of these latter approaches are rigorously identified in terms of (i) the extent of overlap between electric double layers around charged neighbouring sites, and (ii) the magnitude of the intraparticulate metal concentration gradient. For the first time, the proposed theory integrates the differentiated impact of the local potential around the charged binding sites amidst the overall particle field, together with that of the so-far discarded intraparticulate flux of metal ions. PMID:24336523
Stratification of discharge in noble gases from the viewpoint of the discrete dynamics
Golubovskii, Yu. Pelyukhova, E.; Sigeneger, F.; Nekuchaev, V.
2015-03-15
Based on the analysis of electron phase trajectories in sinusoidal electric fields, a new point of view on discharge stratification is proposed. It is shown that the positive column can be considered as a spatial resonator in which electric fields with a fundamental period length L{sub S} or higher mode length q/p L{sub S} establish, where p and q are integers and p > q. The fundamental mode length L{sub S} is equivalent to the distance on which electrons gain energy equal to the lowest excitation threshold. This distance determines a length of the S-striation. Unlike kinetic theory, in the presented model resonance properties of the discharge column are not connected with elastic collision energy losses. A point map is used to obtain the resonance trajectories of electrons in the phase plane. Stable points for the positions of inelastic collisions in the resonance trajectories have been found at the positions of field maxima in the case of integer ratios p/q . For non-integer ratios p/q , multiple resonance trajectories arise according to a more complex stability criterion. From this point of view, S-, P-, and R-striations in noble gas discharges can be explained. Due to energy losses in elastic collisions, initial electron energy distribution functions converge to the resonance trajectories (the so-called “bunch effect”). The findings of the discrete model agree with results of kinetic theory and experiment. The new approach avoids difficulties of the kinetic theory in the case of exceptionally large relaxation lengths which can even exceed the positive column length.
Minimax terminal approach problem in two-level hierarchical nonlinear discrete-time dynamical system
Shorikov, A. F.
2015-11-30
We consider a discrete–time dynamical system consisting of three controllable objects. The motions of all objects are given by the corresponding vector nonlinear or linear discrete–time recurrent vector relations, and control system for its has two levels: basic (first or I level) that is dominating and subordinate level (second or II level) and both have different criterions of functioning and united a priori by determined informational and control connections defined in advance. For the dynamical system in question, we propose a mathematical formalization in the form of solving a multistep problem of two-level hierarchical minimax program control over the terminal approach process with incomplete information and give a general scheme for its solving.
pH-replica exchange molecular dynamics in proteins using a discrete protonation method.
Sabri Dashti, Danial; Meng, Yilin; Roitberg, Adrian E
2012-08-01
Protonation equilibria in biological molecules modulates structure, dynamics, and function. A pH-replica exchange molecular dynamics (pH-REMD) method is described here to improve the coupling between conformational and protonation sampling. Under a Hamiltonian replica exchange setup, conformations are swapped between two neighboring replicas, which themselves are at different pHs. The method has been validated on a series of biological systems. We applied pH-REMD to a series of model compounds, to an terminally charged ADFDA pentapeptide, and to a heptapeptide derived from the ovomucoid third domain (OMTKY3). In all of those systems, the predicted pK(a) by pH-REMD is very close to the experimental value and almost identical to the ones obtained by constant pH molecular dynamics (CpH MD). The method presented here, pH-REMD, has the advantage of faster convergence properties due to enhanced sampling of both conformation and protonation spaces. PMID:22694266
NASA Astrophysics Data System (ADS)
Alonso-Marroquín, Fernando; Galindo-Torres, Sergio-Andres; Tordesillas, Antoinette; Wang, Yucang
2009-06-01
One of the most challenging problems in the realistic modeling of granular materials is how to capture the real shape of the particles. Here we present a method to simulate systems with complex-shaped particles. This method integrates developments in two traditionally separate research areas: computational geometry and molecular dynamics. The computational geometry involves the implementation of techniques of computer graphics to represent particle shape and collision detection. Traditional techniques from molecular dynamics are used to integrate the equations of motion and to perform an efficient calculation of contact forces. The algorithm to solve the dynamics of the system is much more efficient, accurate and easier to implement than other models. The algorithm is used to simulate quasistatic deformation of granular materials using two different models. The first model consists of non-circular particles interacting via frictional forces. The second model consists of circular particles interacting via rolling and sliding resistance. The comparison of both models help us to understand and quantify the extend to which the effects of particle shape can be captured by the introduction of artificial rolling resistance on circular particles. Biaxial test simulation show that the overall response of the system and the collapse of force chains at the critical state is qualitatively similar in both 2D and 3D simulations.
Dynamically Coupled Food-web and Hydrodynamic Modeling with ADH-CASM
NASA Astrophysics Data System (ADS)
Piercy, C.; Swannack, T. M.
2012-12-01
Oysters and freshwater mussels are "ecological engineers," modifying the local water quality by filtering zooplankton and other suspended particulate matter from the water column and flow hydraulics by impinging on the near-bed flow environment. The success of sessile, benthic invertebrates such as oysters depends on environmental factors including but not limited to temperature, salinity, and flow regime. Typically food-web and other types of ecological models use flow and water quality data as direct input without regard to the feedback between the ecosystem and the physical environment. The USACE-ERDC has developed a coupled hydrodynamic-ecological modeling approach that dynamically couples a 2-D hydrodynamic and constituent transport model, Adaptive Hydraulics (ADH), with a bioenergetics food-web model, the Comprehensive Aquatics Systems Model (CASM), which captures the dynamic feedback between aquatic ecological systems and the environment. We present modeling results from restored oyster reefs in the Great Wicomico River on the western shore of the Chesapeake Bay, which quantify ecosystem services such as the influence of the benthic ecosystem on water quality. Preliminary results indicate that while the influence of oyster reefs on bulk flow dynamics is limited due to the localized influence of oyster reefs, large reefs and the associated benthic ecosystem can create measurable changes in the concentrations of nitrogen, phosphorus, and carbon in the areas around reefs. We also present a sensitivity analysis to quantify the relative sensitivity of the coupled ADH-CASM model to both hydrodynamic and ecological parameter choice.
Web-based experiments for the study of collective social dynamics in cultural markets
Salganik, Matthew J.; Watts, Duncan J.
2013-01-01
Social scientists are often interested in understanding how the dynamics of social systems are driven by the behavior of individuals that make up those systems. However, this process is hindered by the difficulty of experimentally studying how individual behavioral tendencies lead to collective social dynamics in large groups of people interacting over time. In this paper we investigate the role of social influence, a process well studied at the individual level, on the puzzling nature of success for cultural products such as books, movies, and music. Using a “multiple-worlds” experimental design we are able to isolate the causal effect of an individual level mechanism on collective social outcomes. We employ this design in a web-based experiment in which 2,930 participants listened to, rated, and download 48 songs by up-and-coming bands. Surprisingly, despite relatively large differences in the demographics, behavior, and preferences of participants, the experimental results at both the individual and collective level were similar to those found in Salganik, Dodds, and Watts (2006). Further, by comparing results from two distinct pools of participants we are able to gain new insights into the role of individual behavior on collective outcomes. We conclude with a discussion of the strengths and weaknesses of web-based experiments to address questions of collective social dynamics. PMID:25164996
Web-based experiments for the study of collective social dynamics in cultural markets.
Salganik, Matthew J; Watts, Duncan J
2009-07-01
Social scientists are often interested in understanding how the dynamics of social systems are driven by the behavior of individuals that make up those systems. However, this process is hindered by the difficulty of experimentally studying how individual behavioral tendencies lead to collective social dynamics in large groups of people interacting over time. In this study, we investigate the role of social influence, a process well studied at the individual level, on the puzzling nature of success for cultural products such as books, movies, and music. Using a "multiple-worlds" experimental design, we are able to isolate the causal effect of an individual-level mechanism on collective social outcomes. We employ this design in a Web-based experiment in which 2,930 participants listened to, rated, and downloaded 48 songs by up-and-coming bands. Surprisingly, despite relatively large differences in the demographics, behavior, and preferences of participants, the experimental results at both the individual and collective levels were similar to those found in Salganik, Dodds, and Watts (2006). Further, by comparing results from two distinct pools of participants, we are able to gain new insights into the role of individual behavior on collective outcomes. We conclude with a discussion of the strengths and weaknesses of Web-based experiments to address questions of collective social dynamics. PMID:25164996
Development of a laboratory niche Web site.
Dimenstein, Izak B; Dimenstein, Simon I
2013-10-01
This technical note presents the development of a methodological laboratory niche Web site. The "Grossing Technology in Surgical Pathology" (www.grossing-technology.com) Web site is used as an example. Although common steps in creation of most Web sites are followed, there are particular requirements for structuring the template's menu on methodological laboratory Web sites. The "nested doll principle," in which one object is placed inside another, most adequately describes the methodological approach to laboratory Web site design. Fragmentation in presenting the Web site's material highlights the discrete parts of the laboratory procedure. An optimally minimal triad of components can be recommended for the creation of a laboratory niche Web site: a main set of media, a blog, and an ancillary component (host, contact, and links). The inclusion of a blog makes the Web site a dynamic forum for professional communication. By forming links and portals, cloud computing opens opportunities for connecting a niche Web site with other Web sites and professional organizations. As an additional source of information exchange, methodological laboratory niche Web sites are destined to parallel both traditional and new forms, such as books, journals, seminars, webinars, and internal educational materials. PMID:23769601
Zhang, Hao; Douglas, Jack F.
2012-01-01
Recent studies of the dynamics of diverse condensed amorphous materials have indicated significant heterogeneity in the local mobility and a progressive increase in collective particle motion upon cooling that takes the form of string-like particle rearrangements. In a previous paper (Part I), we examined the possibility that fluctuations in potential energy E and particle mobility μ associated with this ‘dynamic heterogeneity’ might offer information about the scale of collective motion in glassy materials based on molecular dynamics simulations of the glassy interfacial region of Ni nanoparticles (NPs) at elevated temperatures. We found that the noise exponent associated with fluctuations in the Debye-Waller factor, a mobility related quantity, was directly proportional to the scale of collective motion L under a broad range of conditions, but the noise exponent associated with E(t) fluctuations was seemingly unrelated to L. In the present work, we focus on this unanticipated difference between potential energy and mobility fluctuations by examining these quantities at an atomic scale. We find that the string atoms exhibit a jump-like motion between two well-separated bands of energy states and the rate at which these jumps occur seems to be consistent with the phenomenology of the ‘slow-beta’ relaxation process of glass-forming liquids. Concurrently with these local E(t) jumps, we also find ‘quake-like’ particle displacements having a power-law distribution in magnitude so that particle displacement fluctuations within the strings are strikingly different from local E(t) fluctuations. An analysis of these E(t) fluctuations suggests that we are dealing with ‘discrete breather’ excitations in which large energy fluctuations develop in arrays of non-linear oscillators by virtue of large anharmonicity in the interparticle interactions and discreteness effects associated with particle packing. We quantify string collective motions on a fast caging
Dynamics of Leslie-Gower type generalist predator in a tri-trophic food web system
NASA Astrophysics Data System (ADS)
Priyadarshi, A.; Gakkhar, S.
2013-11-01
In this paper, the dynamics of a tri-trophic food web system consists of Leslie-Gower type generalist predator has been explored. The system is bounded under certain conditions. The Hopf-bifurcation has been established in the phase planes. The bifurcation diagrams exhibit coexistence of all three species in the form of periodic/chaotic solutions. The "snail-shell" chaotic attractor has very high Lyapunov exponents. The coexistence in the form of stable equilibrium is also possible for lower values of parameters. The two-parameter bifurcation diagrams are drawn for critical parameters.
PCB Food Web Dynamics Quantify Nutrient and Energy Flow in Aquatic Ecosystems.
McLeod, Anne M; Paterson, Gordon; Drouillard, Ken G; Haffner, G Douglas
2015-11-01
Measuring in situ nutrient and energy flows in spatially and temporally complex aquatic ecosystems represents a major ecological challenge. Food web structure, energy and nutrient budgets are difficult to measure, and it is becoming more important to quantify both energy and nutrient flow to determine how food web processes and structure are being modified by multiple stressors. We propose that polychlorinated biphenyl (PCB) congeners represent an ideal tracer to quantify in situ energy and nutrient flow between trophic levels. Here, we demonstrate how an understanding of PCB congener bioaccumulation dynamics provides multiple direct measurements of energy and nutrient flow in aquatic food webs. To demonstrate this novel approach, we quantified nitrogen (N), phosphorus (P) and caloric turnover rates for Lake Huron lake trout, and reveal how these processes are regulated by both growth rate and fish life history. Although minimal nutrient recycling was observed in young growing fish, slow growing, older lake trout (>5 yr) recycled an average of 482 Tonnes·yr(-1) of N, 45 Tonnes·yr(-1) of P and assimilated 22 TJ yr(-1) of energy. Compared to total P loading rates of 590 Tonnes·yr(-1), the recycling of primarily bioavailable nutrients by fish plays an important role regulating the nutrient states of oligotrophic lakes. PMID:26437236
ERIC Educational Resources Information Center
Wang, Tzu-Hua
2010-01-01
This research combines the idea of cake format dynamic assessment defined by Sternberg and Grigorenko (2001) and the "graduated prompt approach" proposed by (Campione and Brown, 1985) and (Campione and Brown, 1987) to develop a multiple-choice Web-based dynamic assessment system. This research adopts a quasi-experimental design to investigate the…
Turing pattern dynamics and adaptive discretization for a super-diffusive Lotka-Volterra model.
Bendahmane, Mostafa; Ruiz-Baier, Ricardo; Tian, Canrong
2016-05-01
In this paper we analyze the effects of introducing the fractional-in-space operator into a Lotka-Volterra competitive model describing population super-diffusion. First, we study how cross super-diffusion influences the formation of spatial patterns: a linear stability analysis is carried out, showing that cross super-diffusion triggers Turing instabilities, whereas classical (self) super-diffusion does not. In addition we perform a weakly nonlinear analysis yielding a system of amplitude equations, whose study shows the stability of Turing steady states. A second goal of this contribution is to propose a fully adaptive multiresolution finite volume method that employs shifted Grünwald gradient approximations, and which is tailored for a larger class of systems involving fractional diffusion operators. The scheme is aimed at efficient dynamic mesh adaptation and substantial savings in computational burden. A numerical simulation of the model was performed near the instability boundaries, confirming the behavior predicted by our analysis. PMID:26219250
NASA Astrophysics Data System (ADS)
Yang, Hui; Li, Zhenhuan; Huang, Minsheng
2014-12-01
Unlike common single crystals, the nickel-based single crystal superalloy shows surprisingly anomalous flow strength (i.e. with the increase of temperature, the yield strength first increases to a peak value and then decreases) and tension-compression (TC) asymmetry. A comprehensive three-dimensional discrete dislocation dynamics (3D-DDD) procedure was developed to model these abnormal mechanical properties. For this purpose, a series of complicated dynamic evolution details of Kear-Wilsdorf (KW) locks, which are closely related to the flow strength anomaly and TC asymmetry, were incorporated into this 3D-DDD framework. Moreover, the activation of the cubic slip system, which is the origin of the decrease in yield strength with increasing temperature at relatively high temperatures, was especially taken into account by introducing a competition criterion between the unlocking of the KW locks and the activation of the cubic slip system. To test our framework, a series of 3D-DDD simulations were performed on a representative volume cell model with a cuboidal Ni3Al precipitate phase embedded in a nickel matrix. Results show that the present 3D-DDD procedure can successfully capture the dynamic evolution of KW locks, the flow strength anomaly and TC asymmetry. Then, the underlying dislocation mechanisms leading to these abnormal mechanical responses were investigated and discussed in detail. Finally, a cyclic deformation of the nickel-based single crystal superalloy was modeled by using the present DDD model, with a special focus on the influence of KW locks on the Bauschinger effect and cyclic softening.
The effects of mixotrophy on the stability and dynamics of a simple planktonic food web
Jost, Christian; Lawrence, Cathryn A.; Campolongo, Francesca; Wouter, van de Bund; Hill, Sheryl; DeAngelis, Donald L.
2004-01-01
Recognition of the microbial loop as an important part of aquatic ecosystems disrupted the notion of simple linear food chains. However, current research suggests that even the microbial loop paradigm is a gross simplification of microbial interactions due to the presence of mixotrophs—organisms that both photosynthesize and graze. We present a simple food web model with four trophic species, three of them arranged in a food chain (nutrients–autotrophs–herbivores) and the fourth as a mixotroph with links to both the nutrients and the autotrophs. This model is used to study the general implications of inclusion of the mixotrophic link in microbial food webs and the specific predictions for a parameterization that describes open ocean mixed layer plankton dynamics. The analysis indicates that the system parameters reside in a region of the parameter space where the dynamics converge to a stable equilibrium rather than displaying periodic or chaotic solutions. However, convergence requires weeks to months, suggesting that the system would never reach equilibrium in the ocean due to alteration of the physical forcing regime. Most importantly, the mixotrophic grazing link seems to stabilize the system in this region of the parameter space, particularly when nutrient recycling feedback loops are included.
Identification of continuous-time models for nonlinear dynamic systems from discrete data
NASA Astrophysics Data System (ADS)
Guo, Yuzhu; Guo, Ling Zhong; Billings, Stephen A.; Wei, Hua-Liang
2016-09-01
A new iOFR-MF (iterative orthogonal forward regression--modulating function) algorithm is proposed to identify continuous-time models from noisy data by combining the MF method and the iOFR algorithm. In the new method, a set of candidate terms, which describe different dynamic relationships among the system states or between the input and output, are first constructed. These terms are then modulated using the MF method to generate the data matrix. The iOFR algorithm is next applied to build the relationships between these modulated terms, which include detecting the model structure and estimating the associated parameters. The relationships between the original variables are finally recovered from the model of the modulated terms. Both nonlinear state-space models and a class of higher order nonlinear input-output models are considered. The new direct method is compared with the traditional finite difference method and results show that the new method performs much better than the finite difference method. The new method works well even when the measurements are severely corrupted by noise. The selection of appropriate MFs is also discussed.
NASA Astrophysics Data System (ADS)
Xie, Zheng; Ye, Zheng; Ma, Yu-Jie
2009-12-01
Numerical simulation of antennae is a topic in computational electromagnetism, which is concerned with the numerical study of Maxwell equations. By discrete exterior calculus and the lattice gauge theory with coefficient R, we obtain the Bianchi identity on prism lattice. By defining an inner product of discrete differential forms, we derive the source equation and continuity equation. Those equations compose the discrete Maxwell equations in vacuum case on discrete manifold, which are implemented on Java development platform to simulate the Gaussian pulse radiation on antennaes.
NASA Astrophysics Data System (ADS)
D'Agostino, Stefania; Della Sala, Fabio; Andreani, Lucio Claudio
2013-05-01
A theoretical control of the electromagnetic coupling between localized surface plasmons and pointlike sources of radiation is a relevant topic in nanoscience and nanophotonics. In this paper a numerical approach based on the discrete dipole approximation is presented as a practical and reliable computational tool to study the decay dynamics of a dipole when it is located in the near proximities of metallic nanoparticles whose shapes do not allow a fully analytical treatment. The method is first applied to Ag nanospheres and nanoshells, which represent two analytically solvable cases, and it is shown to lead to a very good agreement with exact results. The approach is then used to consider the response, in terms of perturbations induced on the radiative and nonradiative decay rates, of elongated nanoparticles, like Ag prolate spheroids and nanocones. Results demonstrate how the optical response of conically shaped nanoparticles can be affected by the distance and the orientation of the emitter of radiation, as well as by other geometrical parameters. The particular symmetry of these plasmonic objects results in peculiar features: the absorption efficiencies of the modes depend on the distance of the source of radiation in a counterintuitive way, and this is explained in terms of the excited charge density distributions. The possibility to simulate arbitrary-shaped nanostructures and several dipole-metal configurations presented here, could thus open new avenues for an aware use of surface plasmons in fluorescence spectroscopy applications or single photon emission studies.
Assmann, Sarah M; Albert, Réka
2009-01-01
A major aim of systems biology is the study of the inter-relationships found within and between large biological data sets. Here we describe one systems biology method, in which the tools of network analysis and discrete dynamic (Boolean) modeling are used to develop predictive models of cellular signaling in cases where detailed temporal and kinetic information regarding the propagation of the signal through the system is lacking. This approach is also applicable to data sets derived from some other types of biological systems, such as transcription factor-mediated regulation of gene expression during the control of developmental fate, or host defense responses following pathogen attack, and is equally applicable to plant and non-plant systems. The method also allows prediction of how elimination of one or more individual signaling components will affect the ultimate outcome, thus allowing the researcher to model the effects of genetic knockout or pharmacological block. The method also serves as a starting point from which more quantitative models can be developed as additional information becomes available. PMID:19588107
Fan, Haidong; Aubry, Sylvie; Arsenlis, Athanasios; El-Awady, Jaafar A.
2015-04-13
The mechanical response of micro-twinned polycrystalline magnesium was studied through three-dimensional discrete dislocation dynamics (DDD). A systematic interaction model between dislocations and (1012) tension twin boundaries (TBs) was proposed and introduced into the DDD framework. In addition, a nominal grain boundary (GB) model agreeing with experimental results was also introduced to mimic the GB’s barrier effect. The current simulation results show that TBs act as a strong obstacle to gliding dislocations, which contributes significantly to the hardening behavior of magnesium. On the other hand, the deformation accommodated by twinning plays a softening role. Therefore, the concave shape of the Mg stress-strain curve results from the competition between dislocation-TB induced hardening and twinning deformation induced softening. At low strain levels, twinning deformation induced softening dominates and a decreasing hardening rate is observed in Stage-I. In Stage-II, both the hardening and softening effects decline, but twinning deformation induced softening declines faster, which leads to an increasing hardening rate.
Fan, Haidong; Aubry, Sylvie; Arsenlis, Athanasios; El-Awady, Jaafar A.
2015-04-13
The mechanical response of micro-twinned polycrystalline magnesium was studied through three-dimensional discrete dislocation dynamics (DDD). A systematic interaction model between dislocations and (1012) tension twin boundaries (TBs) was proposed and introduced into the DDD framework. In addition, a nominal grain boundary (GB) model agreeing with experimental results was also introduced to mimic the GB’s barrier effect. The current simulation results show that TBs act as a strong obstacle to gliding dislocations, which contributes significantly to the hardening behavior of magnesium. On the other hand, the deformation accommodated by twinning plays a softening role. Therefore, the concave shape of the Mgmore » stress-strain curve results from the competition between dislocation-TB induced hardening and twinning deformation induced softening. At low strain levels, twinning deformation induced softening dominates and a decreasing hardening rate is observed in Stage-I. In Stage-II, both the hardening and softening effects decline, but twinning deformation induced softening declines faster, which leads to an increasing hardening rate.« less
NASA Astrophysics Data System (ADS)
Akiyama, Terunobu; de Rooij, Nicolaas F.; Staufer, Urs; Detterbeck, Manfred; Braendlin, Dominik; Waldmeier, Simon; Scheidiger, Martin
2010-06-01
The quartz tuning fork based probe {e.g., Akiyama et al. [Appl. Surf. Sci. 210, 18 (2003)]}, termed "A-Probe," is a self-sensing and self-actuating (exciting) probe for dynamic mode atomic force microscope (AFM) operation. It is an oscillatory force sensor consisting of the two discrete resonators. This paper presents the investigations on an improved A-Probe: its batch fabrication and assembly, mounting on an AFM head, electrical setup, characterization, and AFM imaging. The fundamental features of the A-Probe are electrically and optically characterized in "approach-withdraw" experiments. Further investigations include the frequency response of an A-Probe to small mechanical vibrations externally applied to the tip and the effective loading force yielding between the tip and the sample during the periodic contact. Imaging of an electronic chip, a compact disk stamper, carbon nanotubes, and Si beads is demonstrated with this probe at ambient conditions in the so-called frequency modulation mode. A special probe substrate, which can snap on a receptacle fixed on an AFM head, and a special holder including a preamplifier electronic are introduced. We hope that the implementation and characterization of the A-Probe described in this paper will provide hints for new scanning probe techniques.
NASA Astrophysics Data System (ADS)
Gao, Siwen; Fivel, Marc; Ma, Anxin; Hartmaier, Alexander
2015-03-01
In the characteristic γ / γ ‧ microstructure of single crystal superalloys, misfit stresses occur due to a significant lattice mismatch of those two phases. The magnitude of this lattice mismatch depends on the chemical composition of both phases as well as on temperature. Furthermore, the lattice mismatch of γ and γ ‧ phases can be either positive or negative in sign. The internal stresses caused by such lattice mismatch play a decisive role for the micromechanical processes that lead to the observed macroscopic athermal deformation behavior of these high-temperature alloys. Three-dimensional discrete dislocation dynamics (DDD) simulations are applied to investigate dislocation glide in γ matrix channels and shearing of γ ‧ precipitates by superdislocations under externally applied uniaxial stresses, by fully taking into account internal misfit stresses. Misfit stress fields are calculated by the fast Fourier transformation (FFT) method and hybridized with DDD simulations. For external loading along the crystallographic [001] direction of the single crystal, it was found that the different internal stress states for negative and positive lattice mismatch result in non-uniform dislocation movement and different dislocation patterns in horizontal and vertical γ matrix channels. Furthermore, positive lattice mismatch produces a lower deformation rate than negative lattice mismatch under the same tensile loading, but for an increasing magnitude of lattice mismatch, the deformation resistance always diminishes. Hence, the best deformation performance is expected to result from alloys with either small positive, or even better, vanishing lattice mismatch between γ and γ ‧ phase.
NASA Astrophysics Data System (ADS)
Liang, Xing; Jiang, Jifa
The asymptotic behavior of discrete type-K monotone dynamical systems and reaction-diffusion equations is investigated. The studying content includes the index theory for fixed points, permanence, global stability, convergence everywhere and coexistence. It is shown that the system has a globally asymptotically stable fixed point if every fixed point is locally asymptotically stable with respect to the face it belongs to and at this point the principal eigenvalue of the diagonal partial derivative about any component not belonging to the face is not one. A nice result presented is the sufficient and necessary conditions for the system to have a globally asymptotically stable positive fixed point. It can be used to establish the sufficient conditions for the system to persist uniformly and the convergent result for all orbits. Applications are made to time-periodic Lotka-Volterra systems with diffusion, and sufficient conditions for such systems to have a unique positive periodic solution attracting all positive initial value functions are given. For more general time-periodic type-K monotone reaction-diffusion systems with spatial homogeneity, a simple condition is given to guarantee the convergence of all positive solutions.
Discrete Inverse and State Estimation Problems
NASA Astrophysics Data System (ADS)
Wunsch, Carl
2006-06-01
The problems of making inferences about the natural world from noisy observations and imperfect theories occur in almost all scientific disciplines. This book addresses these problems using examples taken from geophysical fluid dynamics. It focuses on discrete formulations, both static and time-varying, known variously as inverse, state estimation or data assimilation problems. Starting with fundamental algebraic and statistical ideas, the book guides the reader through a range of inference tools including the singular value decomposition, Gauss-Markov and minimum variance estimates, Kalman filters and related smoothers, and adjoint (Lagrange multiplier) methods. The final chapters discuss a variety of practical applications to geophysical flow problems. Discrete Inverse and State Estimation Problems is an ideal introduction to the topic for graduate students and researchers in oceanography, meteorology, climate dynamics, and geophysical fluid dynamics. It is also accessible to a wider scientific audience; the only prerequisite is an understanding of linear algebra. Provides a comprehensive introduction to discrete methods of inference from incomplete information Based upon 25 years of practical experience using real data and models Develops sequential and whole-domain analysis methods from simple least-squares Contains many examples and problems, and web-based support through MIT opencourseware
Synchronous Discrete Harmonic Oscillator
Antippa, Adel F.; Dubois, Daniel M.
2008-10-17
We introduce the synchronous discrete harmonic oscillator, and present an analytical, numerical and graphical study of its characteristics. The oscillator is synchronous when the time T for one revolution covering an angle of 2{pi} in phase space, is an integral multiple N of the discrete time step {delta}t. It is fully synchronous when N is even. It is pseudo-synchronous when T/{delta}t is rational. In the energy conserving hyperincursive representation, the phase space trajectories are perfectly stable at all time scales, and in both synchronous and pseudo-synchronous modes they cycle through a finite number of phase space points. Consequently, both the synchronous and the pseudo-synchronous hyperincursive modes of time-discretization provide a physically realistic and mathematically coherent, procedure for dynamic, background independent, discretization of spacetime. The procedure is applicable to any stable periodic dynamical system, and provokes an intrinsic correlation between space and time, whereby space-discretization is a direct consequence of background-independent time-discretization. Hence, synchronous discretization moves the formalism of classical mechanics towards that of special relativity. The frequency of the hyperincursive discrete harmonic oscillator is ''blue shifted'' relative to its continuum counterpart. The frequency shift has the precise value needed to make the speed of the system point in phase space independent of the discretizing time interval {delta}t. That is the speed of the system point is the same on the polygonal (in the discrete case) and the circular (in the continuum case) phase space trajectories.
The Impact of 850,000 Years of Climate Changes on the Structure and Dynamics of Mammal Food Webs
Nenzén, Hedvig K.; Montoya, Daniel; Varela, Sara
2014-01-01
Most evidence of climate change impacts on food webs comes from modern studies and little is known about how ancient food webs have responded to climate changes in the past. Here, we integrate fossil evidence from 71 fossil sites, body-size relationships and actualism to reconstruct food webs for six large mammal communities that inhabited the Iberian Peninsula at different times during the Quaternary. We quantify the long-term dynamics of these food webs and study how their structure changed across the Quaternary, a period for which fossil data and climate changes are well known. Extinction, immigration and turnover rates were correlated with climate changes in the last 850 kyr. Yet, we find differences in the dynamics and structural properties of Pleistocene versus Holocene mammal communities that are not associated with glacial-interglacial cycles. Although all Quaternary mammal food webs were highly nested and robust to secondary extinctions, general food web properties changed in the Holocene. These results highlight the ability of communities to re-organize with the arrival of phylogenetically similar species without major structural changes, and the impact of climate change and super-generalist species (humans) on Iberian Holocene mammal communities. PMID:25207754
Microbial food web dynamics along a soil chronosequence of a glacier forefield
NASA Astrophysics Data System (ADS)
Esperschütz, J.; Pérez-de-Mora, A.; Schreiner, K.; Welzl, G.; Buegger, F.; Zeyer, J.; Hagedorn, F.; Munch, J. C.; Schloter, M.
2011-11-01
Microbial food webs are critical for efficient nutrient turnover providing the basis for functional and stable ecosystems. However, the successional development of such microbial food webs and their role in "young" ecosystems is unclear. Due to a continuous glacier retreat since the middle of the 19th century, glacier forefields have expanded offering an excellent opportunity to study food web dynamics in soils at different developmental stages. In the present study, litter degradation and the corresponding C fluxes into microbial communities were investigated along the forefield of the Damma glacier (Switzerland). 13C-enriched litter of the pioneering plant Leucanthemopsis alpina (L.) Heywood was incorporated into the soil at sites that have been free from ice for approximately 10, 60, 100 and more than 700 years. The structure and function of microbial communities were identified by 13C analysis of phospholipid fatty acids (PLFA) and phospholipid ether lipids (PLEL). Results showed increasing microbial diversity and biomass, and enhanced proliferation of bacterial groups as ecosystem development progressed. Initially, litter decomposition proceeded faster at the more developed sites, but at the end of the experiment loss of litter mass was similar at all sites, once the more easily-degradable litter fraction was processed. As a result incorporation of 13C into microbial biomass was more evident during the first weeks of litter decomposition. 13C enrichments of both PLEL and PLFA biomarkers following litter incorporation were observed at all sites, suggesting similar microbial foodwebs at all stages of soil development. Nonetheless, the contribution of bacteria, especially actinomycetes to litter turnover became more pronounced as soil age increased in detriment of archaea, fungi and protozoa, more prominent in recently deglaciated terrain.
Bioaccumulation Dynamics of Arsenate at the Base of Aquatic Food Webs.
Lopez, Adeline R; Hesterberg, Dean R; Funk, David H; Buchwalter, David B
2016-06-21
Periphyton is an important food source at the base of freshwater ecosystems that tends to bioconcentrate trace elements making them trophically available. The potential for arsenic-a trace element of particular concern due to its widespread occurrence, toxicity, and carcinogenicity-to bioconcentrate in periphyton and thus be available to benthic grazers is less well characterized. To better understand arsenate bioaccumulation dynamics in lotic food webs, we used a radiotracer approach to characterize accumulation in periphyton and subsequent trophic transfer to benthic grazers. Periphyton bioconcentrated As between 3,200-9,700-fold (dry weight) over 8 days without reaching steady state, suggesting that periphyton is a major sink for arsenate. However, As-enriched periphyton as a food source for the mayfly Neocloeon triangulifer resulted in negligible As accumulation in a full lifecycle exposure. Additional studies estimate dietary assimilation efficiency in several primary consumers ranging from 22% in the mayfly N. triangulifer to 75% in the mayfly Isonychia sp. X-ray fluorescence mapping revealed that As was predominantly associated with iron oxides in periphyton. We speculate that As adsorption to Fe in periphyton may play a role in reducing dietary bioavailability. Together, these results suggest that trophic movement of As in lotic food webs is relatively low, though species differences in bioaccumulation patterns are important. PMID:27223406