A Phase-field Model for Computer Simulation of Lamellar Structure Formation in y-TiAI

In the present paper, a phase-field model is developed to simulate the formation and evolution of lamellar microstructure in γ-TiAl alloys. The mechanism of formation of TiAl lamellae proposed by Denquin and Naka is incorporated into the model. The model describes the formation and evolution of the face-centered cubic (fcc) stacking lamellar zone followed by the subsequent appearance and growth of the γ-phase, involving both the chemical composition change by atom transfer and the ordering of the fcc lattice. The thermodynamics of the model system and the interaction between the displacive and diffusional transformations are described by a non-equilibrium free energy formulated as a function of concentration and structural order parameter fields. The long-range elastic interactions, arising from the lattice misfit between the α, fcc (A1) and the various orientation variants of the γ-phase are taken into account by incorporating of the elastic strain energy into the total free energy. Simulation studies based on the model successfully predicted some essential features of the lamellar structure. It is found that the formation and evolution of the lamellar structure are predominantly controlled by the minimization of the elastic energy of the interfaces between the different fcc stacking groups, low-symmetry product phase γ and the high-symmetry α-phase, as well as between the various orientation variants of the product phase.

Population-based local search for protein folding simulation in the MJ energy model and cubic lattices

We present experimental results on benchmark problems in 3D cubic lattice structures with the Miyazawa-Jernigan energy function for two local search procedures that utilise the pull-move set: (i) population-based local search (PLS) that traverses the energy landscape with greedy steps towards (potential) local minima followed by upward steps up to a certain level of the objective function; (ii) simulated annealing with a logarithmic cooling schedule (LSA). The parameter settings for PLS are derived from short LSA-runs executed in pre-processing and the procedure utilises tabu lists generated for each member of the population. In terms of the total number of energy function evaluations both methods perform equally well, however. PLS has the potential of being parallelised with an expected speed-up in the region of the population size. Furthermore, both methods require a significant smaller number of function evaluations when compared to Monte Carlo simulations with kink-jump moves. (C) 2009 Elsevier Ltd. All rights reserved.

An Energy Conserving Finite Difference Scheme for Simulation of Collisions

Nonlinear phenomena play an essential role in the sound production process of many musical instruments. A common source of these effects is object collision, the numerical simulation of which is known to give rise to stability
issues. This paper presents a method to construct numerical schemes that conserve the total energy in simulations of one-mass systems involving collisions, with no conditions imposed on any of the physical or numerical parameters.
This facilitates the adaptation of numerical models to experimental data, and allows a more free parameter adjustment in sound synthesis explorations. The energy preservedness of the proposed method is tested and demonstrated though several examples, including a bouncing ball and a non-linear oscillator, and implications regarding the wider applicability are discussed.

Energy conserving schemes for the simulation of musical instrument contact dynamics

Collisions are an innate part of the function of many musical instruments. Due to the nonlinear nature of contact forces, special care has to be taken in the construction of numerical schemes for simulation and sound synthesis. Finite difference schemes and other time-stepping algorithms used for musical instrument modelling purposes are normally arrived at by discretising a Newtonian description of the system. However because impact forces are non-analytic functions of the phase space variables, algorithm stability can rarely be established this way. This paper presents a systematic approach to deriving energy conserving schemes for frictionless impact modelling. The proposed numerical formulations follow from discretising Hamilton׳s equations of motion, generally leading to an implicit system of nonlinear equations that can be solved with Newton׳s method. The approach is first outlined for point mass collisions and then extended to distributed settings, such as vibrating strings and beams colliding with rigid obstacles. Stability and other relevant properties of the proposed approach are discussed and further demonstrated with simulation examples. The methodology is exemplified through a case study on tanpura string vibration, with the results confirming the main findings of previous studies on the role of the bridge in sound generation with this type of string instrument.

A reactive force field simulation of liquid-liquid phase transitions in phosphorus

A force field model of phosphorus has been developed based on density functional (DF) computations and experimental results, covering low energy forms of local tetrahedral symmetry and more compact (simple cubic) structures that arise with increasing pressure. Rules tailored to DF data for the addition, deletion, and exchange of covalent bonds allow the system to adapt the bonding configuration to the thermodynamic state. Monte Carlo simulations in the N-P-T ensemble show that the molecular (P-4) liquid phase, stable at low pressure P and relatively low temperature T, transforms to a polymeric (gel) state on increasing either P or T. These phase changes are observed in recent experiments at similar thermodynamic conditions, as shown by the close agreement of computed and measured structure factors in the molecular and polymer phases. The polymeric phase obtained by increasing pressure has a dominant simple cubic character, while the polymer obtained by raising T at moderate pressure is tetrahedral. Comparison with DF results suggests that the latter is a semiconductor, while the cubic form is metallic. The simulations show that the T-induced polymerization is due to the entropy of the configuration of covalent bonds, as in the polymerization transition in sulfur. The transition observed with increasing P is the continuation at high T of the black P to arsenic (A17) structure observed in the solid state, and also corresponds to a semiconductor to metal transition. (C) 2004 American Institute of Physics.

Numerical simulation of the charge balance in an EBIT

The electron beam ions traps (EBITs) are widely used to study highly charged ions (HCIs). In an EBIT, a high energy electron beam collides with atoms and ions to generate HCIs in the trap region. It is important to study the physics in the trap. The atomic processes, such as electron impact ionisation (EI), radiative recombination (RR), dielectronic recombination (DR) and charge exchange (CX), occur in the trap and numerical simulation can give some parameters for design, predict the composition and describe charge state evolution in an EBIT [Phys. Rev. A 43 (199 1) 4861]. We are presently developing a new code, which additionally includes a description of the overlaps between the ion clouds of the various charge-states. It has been written so that it can simulate experiments where various machine parameters (e.g. beam energy and current) can vary throughout the simulation and will be able to use cross- sections either based on scaling laws or derived from atomic structure calculations. An object-oriented method is used in developing the new software, which is an efficient way to organize and write code. (C) 2003 Elsevier Science B.V. All rights reserved.

Ab initio simulation of charged slabs at constant chemical potential

We present a practical scheme for performing ab initio supercell calculations of charged slabs at constant electron chemical potential mu, rather than at constant number of electrons N-e. To this end, we define the chemical potential relative to a plane (or "reference electrode") at a finite distance from the slab (the distance should reflect the particular geometry of the situation being modeled). To avoid a net charge in the supercell, and thus make possible a standard supercell calculation, we restore the electroneutrality of the periodically repeated unit by means of a compensating charge, whose contribution to the total energy and potential is subtracted afterwards. The "constant mu" mode enables one to perform supercell calculation on slabs, where the slab is kept at a fixed potential relative to the reference electrode. We expect this to be useful in modeling many experimental situations, especially in electro-chemistry. (C) 2001 American Institute of Physics.

Solvation Structure and Transport of Acidic Protons in Ionic Liquids: A First-principles Simulation Study

Ab initio simulations of a single molecule of HCl in liquid dimethyl imidazolium chloride [dmim][Cl] show that the acidic proton exists as a symmetric, linear ClHCl- species. Details of the solvation structure around this molecule are given. The proton-transfer process was investigated by applying a force along the antisymmetric stretch coordinate until the molecule broke. Changes in the free energy and local solvation structure during this process were investigated. In the reaction mechanism identified, a free chloride approaches the proton from the side. As the original ClHCl- distorts and the incoming chloride forms a new bond to the proton, one of the original chlorine atoms is expelled and a new linear molecule is formed.

Robustness of Damping Control Implemented by Energy Storage Systems Installed in Power Systems

An energy storage system (ESS) installed in a power system can effectively damp power system oscillations through controlling exchange of either active or reactive power between the ESS and power system. This paper investigates the robustness of damping control implemented by the ESS to the variations of power system operating conditions. It proposes a new analytical method based on the well-known equal-area criterion and small-signal stability analysis. By using the proposed method, it is concluded in the paper that damping control implemented by the ESS through controlling its active power exchange with the power system is robust to the changes of power system operating conditions. While if the ESS damping control is realized by controlling its reactive power exchange with the power system, effectiveness of damping control changes with variations of power system operating condition. In the paper, an example power system installed with a battery ESS (BESS) is presented. Simulation results confirm the analytical conclusions made in the paper about the robustness of ESS damping control. Laboratory experiment of a physical power system installed with a 35kJ/7kW SMES (Superconducting Magnetic Energy Storage) was carried out to evaluate theoretical study. Results are given in the paper, which demonstrate that effectiveness of SMES damping control realized through regulating active power is robust to changes of load conditions of the physical power system.

Dynamical simulation of inelastic quantum transport

A method for correlated quantum electron-ion dynamics is combined with a method for electronic open boundaries to simulate in real time the heating, and eventual equilibration at an elevated vibrational energy, of a quantum ion under current flow in an atomic wire, together with the response of the current to the ionic heating. The method can also be used to extract inelastic current voltage corrections under steady-state conditions. However, in its present form the open-boundary method contains an approximation that limits the resolution of current-voltage features. The results of the simulations are tested against analytical results from scattering theory. Directions for the improvement of the method are summarized at the end.

Stochastic protein folding simulation in the three-dimensional HP-model

We present results from three-dimensional protein folding simulations in the HP-model on ten benchmark problems. The simulations are executed by a simulated annealing-based algorithm with a time-dependent cooling schedule. The neighbourhood relation is determined by the pull-move set. The results provide experimental evidence that the maximum depth D of local minima of the underlying energy landscape can be upper bounded by D < n(2/3). The local search procedure employs the stopping criterion (In/delta)(D/gamma) where m is an estimation of the average number of neighbouring conformations, gamma relates to the mean of non-zero differences of the objective function for neighbouring conformations, and 1-delta is the confidence that a minimum conformation has been found. The bound complies with the results obtained for the ten benchmark problems. (c) 2008 Elsevier Ltd. All rights reserved.

Glucose solvation by the ionic liquid 1,3-dimethylimidazolium chloride: A simulation study

Simulations of beta-glucose in the ionic liquid 1,3-dimethylimidazoliurn chloride have been performed in order to examine the solvation environment of the carbohydrate. Both single molecule and 1:5 glucose:ionic liquid (16.7 wt %) solutions are Studied, and the hydrogen bonding between sugar and solvent is examined. The primary solvation shell around the perimeter of the glucose ring consists predominantly of chloride anions which hydrogen bond to the hydroxyl groups. A small presence of the cation is also found, with the association Occurring through the weakly acidic hydrogen at the 2-position of the imidazolium ring interacting with the oxygen atoms of the sugar secondary hydroxyls. An average chloride coordination number of 4 is found around the glucose for both the single molecule and high concentration Simulations, despite the reduced chloride:glucose ratio in the latter case. In relation to the cation, the glucose molecules occupy positions above and below the plane of the imidazolium ring. Importantly, even at high glucose concentrations, no significant change in the anion-cation interactions and overall liquid structure of the ionic liquid is found, indicating that the glucose is readily accommodated by the solvent at this concentration. Dominant contributions to the sugar-ionic liquid interaction energy come from favorable hydrogen bonding (electrostatic) interactions between hydroxyls and chlorides, although a small favorable van der Waals energy contribution is also seen between the sugar and cations suggesting that the cation could be tailored in order to further improve the dissolution of glucose/cellulose in ionic liquid systems.

Simulation of coupled heat and moisture transfer in air-conditioned buildings

The simultaneous heat and moisture transfer in the building envelope has an important influence on the indoor environment and the overall performance of buildings. In this paper, a model for predicting whole building heat and moisture transfer was presented. Both heat and moisture transfer in the building envelope and indoor air were simultaneously considered; their interactions were modeled. The coupled model takes into account most of the main hygrothermal effects in buildings. The coupled system model was implemented in MATLAB-Simulink, and validated by using a series of published testing tools. The new program was applied to investigate the moisture transfer effect on indoor air humidity and building energy consumption under different climates. The results show that the use of more detailed simulation routines can result in improvements to the building's design for energy optimisation through the choice of proper hygroscopic materials, which would not be indicated by simpler calculation techniques.

Numerical simulation of the charge balance and temperature evolution in an electron beam ion trap

A computer code has been developed to simulate and study the evolution of ion charge states inside the trap region of an electron beam ion trap. In addition to atomic physics phenomena previously included in similar codes such as electron impact ionization, radiative recombination, and charge exchange, several aspects of the relevant physics such as dielectronic recombination, ionization heating, and ion cloud expansion have been included for the first time in the model. The code was developed using object oriented concepts with database support, making it readable, accurate, and well organized. The simulation results show a good agreement with various experiments, and give useful information for selection of operating conditions and experiment design.