A manufacturing systems modelling using system dynamics:forming a dedicated modelling tool

Discrete event simulation is a popular aid for manufacturing system design; however in application this technique can sometimes be unnecessarily complex. This paper is concerned with applying an alternative technique to manufacturing system design which may well provide an efficient form of rough-cut analysis. This technique is System Dynamics, and the work described in this paper has set about incorporating the principles of this technique into a computer based modelling tool that is tailored to manufacturing system design. This paper is structured to first explore the principles of System Dynamics and how they differ from Discrete Event Simulation. The opportunity for System Dynamics is then explored, and this leads to defining the capabilities that a suitable tool would need. This specification is then transformed into a computer modelling tool, which is then assessed by applying this tool to model an engine production facility. Read More: http://www.worldscientific.com/doi/abs/10.1142/S0219686703000228

Intracavity dynamics in high-power mode-locked fiber lasers

A theoretical model is developed which characterizes the intracavity pulse evolutions in high-power fiber lasers. It is shown that experimentally observed dynamics of the key pulse parameters can be described by a reduced model of ordinary differential equations. Critical in driving the intracavity dynamics is the amplitude and phase modulations generated by the discrete elements in the laser. The theory gives a simple geometrical description of the intracavity dynamics and possible operation modes of the laser cavity. Furthermore, it provides a simple and efficient method for optimizing the performance of complex multiparametric laser systems.

Lattice approach to the dynamics of phase-coded soliton trains

We address the collective dynamics of a soliton train propagating in a medium described by the nonlinear Schrödinger equation. Our approach uses the reduction of train dynamics to the discrete complex Toda chain (CTC) model for the evolution of parameters for each train constituent: such a simplification allows one to carry out an approximate analysis of the dynamics of positions and phases of individual interacting pulses. Here, we employ the CTC model to the problem which has relevance to the field of fibre optics communications where each binary digit of transmitted information is encoded via the phase difference between the two adjacent solitons. Our goal is to elucidate different scenarios of the train distortions and the subsequent information garbling caused solely by the intersoliton interactions. First, we examine how the structure of a given phase pattern affects the initial stage of the train dynamics and explain the general mechanisms for the appearance of unstable collective soliton modes. Then we further discuss the nonlinear regime concentrating on the dependence of the Lax scattering matrix on the input phase distribution; this allows one to classify typical features of the train evolution and determine the distance where the soliton escapes from its slot. In both cases, we demonstrate deep mathematical analogies with the classical theory of crystal lattice dynamics.

Complexity of classical dynamics of molecular systems. II. Finite statistical complexity of a water-Na+ system

The computational mechanics approach has been applied to the orientational behavior of water molecules in a molecular dynamics simulated water–Na + system. The distinctively different statistical complexity of water molecules in the bulk and in the first solvation shell of the ion is demonstrated. It is shown that the molecules undergo more complex orientational motion when surrounded by other water molecules compared to those constrained by the electric field of the ion. However the spatial coordinates of the oxygen atom shows the opposite complexity behavior in that complexity is higher for the solvation shell molecules. New information about the dynamics of water molecules in the solvation shell is provided that is additional to that given by traditional methods of analysis.

Large-scale molecular dynamics simulations of HLA-A*0201 complexed with a tumor-specific antigenic peptide:can the α3 and β2m domains be neglected?

Large-scale massively parallel molecular dynamics (MD) simulations of the human class I major histo-compatibility complex (MHC) protein HLA-A*0201 bound to a decameric tumor-specific antigenic peptide GVY-DGREHTV were performed using a scalable MD code on high-performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available Large-scale massively parallel molecular dynamics (MD) simulations of the human class I major histocompatibility complex (MHC) protein HLA-A*0201 bound to a decameric tumor-specific antigenic peptide GVYDGREHTV were performed using a scalable MD code on high-performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available for this study allow us to compare directly differences in the behavior of very large molecular models; in this case, the entire extracellular portion of the peptide–MHC complex vs. the isolated peptide binding domain. Comparison of the results from the partial and the whole system simulations indicates that the peptide is less tightly bound in the partial system than in the whole system. From a detailed study of conformations, solvent-accessible surface area, the nature of the water network structure, and the binding energies, we conclude that, when considering the conformation of the α1–α2 domain, the α3 and β2m domains cannot be neglected. © 2004 Wiley Periodicals, Inc. J Comput Chem 25: 1803–1813, 2004

Travelling waves in a complex Ginzburg-Landau equation with time-delay feedback

One of the simplest ways to create nonlinear oscillations is the Hopf bifurcation. The spatiotemporal dynamics observed in an extended medium with diffusion (e.g., a chemical reaction) undergoing this bifurcation is governed by the complex Ginzburg-Landau equation, one of the best-studied generic models for pattern formation, where besides uniform oscillations, spiral waves, coherent structures and turbulence are found. The presence of time delay terms in this equation changes the pattern formation scenario, and different kind of travelling waves have been reported. In particular, we study the complex Ginzburg-Landau equation that contains local and global time-delay feedback terms. We focus our attention on plane wave solutions in this model. The first novel result is the derivation of the plane wave solution in the presence of time-delay feedback with global and local contributions. The second and more important result of this study consists of a linear stability analysis of plane waves in that model. Evaluation of the eigenvalue equation does not show stabilisation of plane waves for the parameters studied. We discuss these results and compare to results of other models.

Applied Aspects of Mathematical Modeling and Optimization of Dynamics of Charged Beams

The complex of questions connected with the analysis, estimation and structural-parametrical optimization of dynamic system is considered in this article. Connection of such problems with tasks of control by beams of trajectories is emphasized. The special attention is concentrated on the review and analysis of spent scientific researches, the attention is stressed to their constructability and applied directedness. Efficiency of the developed algorithmic and software is demonstrated on the tasks of modeling and optimization of output beam characteristics in linear resonance accelerators.

Labeled EF-Tus for rapid kinetic studies of pretranslocation complex formation

The universally conserved translation elongation factor EF-Tu delivers aminoacyl(aa)-tRNA in the form of an aa-tRNA·EF-Tu·GTP ternary complex (TC) to the ribosome where it binds to the cognate mRNA codon within the ribosomal A-site, leading to formation of a pretranslocation (PRE) complex. Here we describe preparation of QSY9 and Cy5 derivatives of the variant E348C-EF-Tu that are functional in translation elongation. Together with fluorophore derivatives of aa-tRNA and of ribosomal protein L11, located within the GTPase associated center (GAC), these labeled EF-Tus allow development of two new FRET assays that permit the dynamics of distance changes between EF-Tu and both L11 (Tu-L11 assay) and aa-tRNA (Tu-tRNA assay) to be determined during the decoding process. We use these assays to examine: (i) the relative rates of EF-Tu movement away from the GAC and from aa-tRNA during decoding, (ii) the effects of the misreading-inducing antibiotics streptomycin and paromomycin on tRNA selection at the A-site, and (iii) how strengthening the binding of aa-tRNA to EF-Tu affects the rate of EF-Tu movement away from L11 on the ribosome. These FRET assays have the potential to be adapted for high throughput screening of ribosomal antibiotics.

Histidine hydrogen bonding in MHC at pH 5 and pH 7 modeled by molecular docking and molecular dynamics simulations

Hydrogen bonds play important roles in maintaining the structure of proteins and in the formation of most biomolecular protein-ligand complexes. All amino acids can act as hydrogen bond donors and acceptors. Among amino acids, Histidine is unique, as it can exist in neutral or positively charged forms within the physiological pH range of 5.0 to 7.0. Histidine can thus interact with other aromatic residues as well as forming hydrogen bonds with polar and charged residues. The ability of His to exchange a proton lies at the heart of many important functional biomolecular interactions, including immunological ones. By using molecular docking and molecular dynamics simulation, we examine the influence of His protonation/deprotonation on peptide binding affinity to MHC class II proteins from locus HLA-DP. Peptide-MHC interaction underlies the adaptive cellular immune response, upon which the next generation of commercially-important vaccines will depend. Consistent with experiment, we find that peptides containing protonated His residues bind better to HLA-DP proteins than those with unprotonated His. Enhanced binding at pH 5.0 is due, in part, to additional hydrogen bonds formed between peptide His+ and DP proteins. In acidic endosomes, protein His79β is predominantly protonated. As a result, the peptide binding cleft narrows in the vicinity of His79β, which stabilizes the peptide - HLA-DP protein complex. © 2014 Bentham Science Publishers.

Generalised Riccati solution and pinning control of complex stochastic networks

This paper considers the global synchronisation of a stochastic version of coupled map lattices networks through an innovative stochastic adaptive linear quadratic pinning control methodology. In a stochastic network, each state receives only noisy measurement of its neighbours' states. For such networks we derive a generalised Riccati solution that quantifies and incorporates uncertainty of the forward dynamics and inverse controller in the derivation of the stochastic optimal control law. The generalised Riccati solution is derived using the Lyapunov approach. A probabilistic approximation type algorithm is employed to estimate the conditional distributions of the state and inverse controller from historical data and quantifying model uncertainties. The theoretical derivation is complemented by its validation on a set of representative examples.

Temporal dynamics in an immunological synapse:role of thermal fluctuations in signaling

The article analyzes the contribution of stochastic thermal fluctuations in the attachment times of the immature T-cell receptor TCR: peptide-major-histocompatibility-complex pMHC immunological synapse bond. The key question addressed here is the following: how does a synapse bond remain stabilized in the presence of high-frequency thermal noise that potentially equates to a strong detaching force? Focusing on the average time persistence of an immature synapse, we show that the high-frequency nodes accompanying large fluctuations are counterbalanced by low-frequency nodes that evolve over longer time periods, eventually leading to signaling of the immunological synapse bond primarily decided by nodes of the latter type. Our analysis shows that such a counterintuitive behavior could be easily explained from the fact that the survival probability distribution is governed by two distinct phases, corresponding to two separate time exponents, for the two different time regimes. The relatively shorter timescales correspond to the cohesion:adhesion induced immature bond formation whereas the larger time reciprocates the association:dissociation regime leading to TCR:pMHC signaling. From an estimate of the bond survival probability, we show that, at shorter timescales, this probability PΔ(τ) scales with time τ as a universal function of a rescaled noise amplitude DΔ2, such that PΔ(τ)∼τ-(ΔD+12),Δ being the distance from the mean intermembrane (T cell:Antigen Presenting Cell) separation distance. The crossover from this shorter to a longer time regime leads to a universality in the dynamics, at which point the survival probability shows a different power-law scaling compared to the one at shorter timescales. In biological terms, such a crossover indicates that the TCR:pMHC bond has a survival probability with a slower decay rate than the longer LFA-1:ICAM-1 bond justifying its stability.

The end of stigma? Understanding the dynamics of legitimisation in the context of TV series consumption

This research contributes to prior work on stigmatisation by looking at stigmatisation and legitimisation as social processes in the context of TV series consumption. Using in-depth interviews, we show that the dynamics of legitimisation are complex and accompanied by the reproduction of existing stigmas and creation of new stigmas.

Influence of reaction atmosphere (H2O, N2, H2, CO2, CO) on fluidized-bed fast pyrolysis of biomass using detailed tar vapor chemistry in computational fluid dynamics

Secondary pyrolysis in fluidized bed fast pyrolysis of biomass is the focus of this work. A novel computational fluid dynamics (CFD) model coupled with a comprehensive chemistry scheme (134 species and 4169 reactions, in CHEMKIN format) has been developed to investigate this complex phenomenon. Previous results from a transient three-dimensional model of primary pyrolysis were used for the source terms of primary products in this model. A parametric study of reaction atmospheres (H2O, N2, H2, CO2, CO) has been performed. For the N2 and H2O atmosphere, results of the model compared favorably to experimentally obtained yields after the temperature was adjusted to a value higher than that used in experiments. One notable deviation versus experiments is pyrolytic water yield and yield of higher hydrocarbons. The model suggests a not overly strong impact of the reaction atmosphere. However, both chemical and physical effects were observed. Most notably, effects could be seen on the yield of various compounds, temperature profile throughout the reactor system, residence time, radical concentration, and turbulent intensity. At the investigated temperature (873 K), turbulent intensity appeared to have the strongest influence on liquid yield. With the aid of acceleration techniques, most importantly dimension reduction, chemistry agglomeration, and in-situ tabulation, a converged solution could be obtained within a reasonable time (∼30 h). As such, a new potentially useful method has been suggested for numerical analysis of fast pyrolysis.

Network evolution and the spatiotemporal dynamics of knowledge sourcing

Knowledge accessing from external organisations is important to firms, especially entrepreneurial ones which often cannot generate internally all the knowledge necessary for innovation. There is, however, a lack of evidence concerning the association between the evolution of firms and the evolution of their networks. The aim of this paper is to begin to fill this gap by undertaking an exploratory analysis of the relationship between the vintage of firms and their knowledge sourcing networks. Drawing on an analysis of firms in the UK, the paper finds some evidence of a U-shaped relationship existing between firm age and the frequency of accessing knowledge from certain sources. Emerging entrepreneurial firms tend to be highly active with regard to accessing knowledge for a range of sources and geographic locations, with the rate of networking dropping somewhat during the period of peak firm growth. For instance, it is found that firms tend to less frequently access knowledge sources such as universities and research institutes in their own region during a stage of peak turnover growth. Overall, the results suggest a complex relationship between the lifecycle of the firm and its networking patterns. It is concluded that policymakers need to become more aware that network formation and utilisation by firms is likely to vary dependent upon their lifecycle position.

Complex Oscillations and Limit Cycles in Autonomous Two-Component Incommensurate Fractional Dynamical Systems

MSC 2010: 26A33, 34D05, 37C25