Generalized Langevin equation: An efficient approach to nonequilibrium molecular dynamics of open systems

The generalized Langevin equation (GLE) has been recently suggested to simulate the time evolution of classical solid and molecular systems when considering general nonequilibrium processes. In this approach, a part of the whole system (an open system), which interacts and exchanges energy with its dissipative environment, is studied. Because the GLE is derived by projecting out exactly the harmonic environment, the coupling to it is realistic, while the equations of motion are non-Markovian. Although the GLE formalism has already found promising applications, e. g., in nanotribology and as a powerful thermostat for equilibration in classical molecular dynamics simulations, efficient algorithms to solve the GLE for realistic memory kernels are highly nontrivial, especially if the memory kernels decay nonexponentially. This is due to the fact that one has to generate a colored noise and take account of the memory effects in a consistent manner. In this paper, we present a simple, yet efficient, algorithm for solving the GLE for practical memory kernels and we demonstrate its capability for the exactly solvable case of a harmonic oscillator coupled to a Debye bath.

External validation of the 2008 Framingham cardiovascular risk equation for CHD and stroke events in a European population of middle-aged men. The PRIME study

To test the applicability of the sex-specific 2008 Framingham general cardiovascular risk equation for coronary heart disease (CHD) and stroke in European middle-aged men from Ireland and France.

Material Logistics in Urban Construction Sites: A Structural Equation Model

Aim: The aim of this paper is to identify best practice relating to the effective management of materials in an urban, confined construction site, using structural equation modelling.

Methodology: A literature review, case study analysis and questionnaire survey are employed, with the results scrutinised using confirmatory factor analysis in the form of structural equation modelling.

Results: The following are the leading strategies in the management of materials in a confined urban site environment; (1) Consult and review the project programme, (2) Effective communication and delivery, (3) Implement site safety management plans, and (4) Proactive spatial monitoring and control.

Implication for Practice: With the relentless expansion of urban centres and the increasing high cost of materials, any potential savings made on-site would translate into significant monetary concessions on completion of a development.

Originality/Value: As on-site project management professionals successfully identify and implement the various strategies in the management of plant and materials on a confined urban site, successful resource management in this restrictive environment is attainable.

Innovative Aspect of Paper: An empirical study of three different construction sites in three different countries (Ireland, England and USA) together with a questionnaire survey from the industry, investigating the managerial strategies in the management of plant and material in confined urban site environments

SOLVING THE INHOMOGENEOUS SCHRODINGER-EQUATION ON A BIDIRECTIONAL PIPELINE OF TRANSPUTERS

In this paper we describe the design of a parallel solution of the inhomogeneous Schrodinger equation, which arises in the construction of continuum orbitals in the R-matrix theory of atomic continuum processes. A prototype system is described which has been programmed in occam2 and implemented on a bi-directional pipeline of transputers. Some timing results for the prototype system are presented, and the development of a full production system is discussed.

Difference equation approach to two-thermocouple sensor characterization in constant velocity flow environments

Thermocouples are one of the most popular devices for temperature measurement due to their robustness, ease of manufacture and installation, and low cost. However, when used in certain harsh environments, for example, in combustion systems and engine exhausts, large wire diameters are required, and consequently the measurement bandwidth is reduced. This article discusses a software compensation technique to address the loss of high frequency fluctuations based on measurements from two thermocouples. In particular, a difference equation (DE) approach is proposed and compared with existing methods both in simulation and on experimental test rig data with constant flow velocity. It is found that the DE algorithm, combined with the use of generalized total least squares for parameter identification, provides better performance in terms of time constant estimation without any a priori assumption on the time constant ratios of the thermocouples.

Applications of the generalized Langevin equation: Towards a realistic description of the baths

The generalized Langevin equation (GLE) method, as developed previously [L. Stella et al., Phys. Rev. B 89, 134303 (2014)], is used to calculate the dissipative dynamics of systems described at the atomic level. The GLE scheme goes beyond the commonly used bilinear coupling between the central system and the bath, and permits us to have a realistic description of both the dissipative central system and its surrounding bath. We show how to obtain the vibrational properties of a realistic bath and how to convey such properties into an extended Langevin dynamics by the use of the mapping of the bath vibrational properties onto a set of auxiliary variables. Our calculations for a model of a Lennard-Jones solid show that our GLE scheme provides a stable dynamics, with the dissipative/relaxation processes properly described. The total kinetic energy of the central system always thermalizes toward the expected bath temperature, with appropriate fluctuation around the mean value. More importantly, we obtain a velocity distribution for the individual atoms in the central system which follows the expected canonical distribution at the corresponding temperature. This confirms that both our GLE scheme and our mapping procedure onto an extended Langevin dynamics provide the correct thermostat. We also examined the velocity autocorrelation functions and compare our results with more conventional Langevin dynamics.

Modeling relativistic soliton interactions in overdense plasmas: A perturbed nonlinear Schrodinger equation framework

We investigate the dynamics of localized solutions of the relativistic cold-fluid plasma model in the small but finite amplitude limit, for slightly overcritical plasma density. Adopting a multiple scale analysis, we derive a perturbed nonlinear Schrodinger equation that describes the evolution of the envelope of circularly polarized electromagnetic field. Retaining terms up to fifth order in the small perturbation parameter, we derive a self-consistent framework for the description of the plasma response in the presence of localized electromagnetic field. The formalism is applied to standing electromagnetic soliton interactions and the results are validated by simulations of the full cold-fluid model. To lowest order, a cubic nonlinear Schrodinger equation with a focusing nonlinearity is recovered. Classical quasiparticle theory is used to obtain analytical estimates for the collision time and minimum distance of approach between solitons. For larger soliton amplitudes the inclusion of the fifth-order terms is essential for a qualitatively correct description of soliton interactions. The defocusing quintic nonlinearity leads to inelastic soliton collisions, while bound states of solitons do not persist under perturbations in the initial phase or amplitude

GPU Acceleration of Partial Differential Equation Solvers

Differential equations are often directly solvable by analytical means only in their one dimensional version. Partial differential equations are generally not solvable by analytical means in two and three dimensions, with the exception of few special cases. In all other cases, numerical approximation methods need to be utilized. One of the most popular methods is the finite element method. The main areas of focus, here, are the Poisson heat equation and the plate bending equation. The purpose of this paper is to provide a quick walkthrough of the various approaches that the authors followed in pursuit of creating optimal solvers, accelerated with the use of graphical processing units, and comparing them in terms of accuracy and time efficiency with existing or self-made non-accelerated solvers.

Experimental investigation of flow and segregation behaviour of bulk solids in silos under high gravity conditions

Many researchers have investigated the flow and segregation behaviour in model scale experimental silos at normal gravity conditions. However it is known that the stresses experienced by the bulk solid in industrial silos are high when compared to model silos. Therefore it is important to understand the effect of stress level on flow and segregation behaviour and establish the scaling laws governing this behaviour. The objective of this paper is to understand the effect of gravity on the flow and segregation behaviour of bulk solids in a silo centrifuge model. The materials used were two mixtures composed of Polyamide and glass beads. The discharge of two bi-disperse bulk solids in a silo centrifuge model were recorded under accelerations ranging from 1g to 15g. The velocity distribution during discharge was evaluated using Particle Image Velocimetry (PIV) techniques and the concentration distribution of large and small particles were obtained by imaging processing techniques. The flow and segregation behaviour at high gravities were then quantified and compared with the empirical equations available in the literature.

Solutions of fourth-order parabolic equation modeling thin film growth

In this paper we study the well-posedness for a fourth-order parabolic equation modeling epitaxial thin film growth. Using Kato's Method [1], [2] and [3] we establish existence, uniqueness and regularity of the solution to the model, in suitable spaces, namelyC⁰([0,T];L^p(Ω)) where  with 1<α<2, n∈N and n≥2. We also show the global existence solution to the nonlinear parabolic equations for small initial data. Our main tools are L_p–L_q-estimates, regularization property of the linear part of e^−tΔ2 and successive approximations. Furthermore, we illustrate the qualitative behavior of the approximate solution through some numerical simulations. The approximate solutions exhibit some favorable absorption properties of the model, which highlight the stabilizing effect of our specific formulation of the source term associated with the upward hopping of atoms. Consequently, the solutions describe well some experimentally observed phenomena, which characterize the growth of thin film such as grain coarsening, island formation and thickness growth.