Numerical simulations of fluid-structure interactions in single-reed mouthpieces

Most single-reed woodwind instrument models rely on a quasistationary approximation to describe the relationship between the volume flow and. the pressure difference across the reed channel. Semiempirical models based on the quasistationary approximation are very useful in explaining the fundamental characteristics of this family of instruments such as self-sustained oscillations and threshold of blowing pressure. However, they fail at explaining more complex phenomena associated with the fluid-structure interaction during dynamic flow regimes, such as the transient and steady-state behavior of the system as a function. of the mouthpiece geometry. Previous studies have discussed the accuracy of the quasistationary approximation but the amount of literature on the subject is sparse, mainly due to the difficulties involved in the measurement of dynamic flows in channels with an oscillating reed. In this paper, a numerical technique based on the lattice Boltzmann method and a finite difference scheme is proposed in order to investigate the characteristics of fully coupled fluid-structure interaction in single-reed mouthpieces with different channel configurations. Results obtained for a stationary simulation with a static reed agree very well with those predicted by the literature based on the quasistationary approximation. However, simulations carried out for a dynamic regime with dn oscillating reed show that the phenomenon associated with flow detachment and reattachment diverges considerably frorn the theoretical assumptions. Furthermore, in the case of long reed channels, the results obtained for the vena contracta factor are in significant disagreement with those predicted by theory. For short channels, the assumption of constant vena contracta was found to be valid for only 40% of the duty cycle. (c) 2007 Acoustical Society of America.

Application of a mechanobiological simulation technique to stents used clinically

Many cardiovascular diseases are characterised by the restriction of blood flow through arteries. Stents can be expanded within arteries to remove such restrictions; however, tissue in-growth into the stent can lead to restenosis. In order to predict the long-term efficacy of stenting, a mechanobiological model of the arterial tissue reaction to stress is required. In this study, a computational model of arterial tissue response to stenting is applied to three clinically relevant stent designs. We ask the question whether such a mechanobiological model can differentiate between stents used clinically, and we compare these predictions to a purely mechanical analysis. In doing so, we are testing the hypothesis that a mechanobiological model of arterial tissue response to injury could predict the long-term outcomes of stent design. Finite element analysis of the expansion of three different stent types was performed in an idealised, 3D artery. Injury was calculated in the arterial tissue using a remaining-life damage mechanics approach. The inflammatory response to this initial injury was modelled using equations governing variables which represented tissue-degrading species and growth factors. Three levels of inflammation response were modelled to account for inter-patient variability. A lattice-based model of smooth muscle cell behaviour was implemented, treating cells as discrete agents governed by local rules. The simulations predicted differences between stent designs similar to those found in vivo. It showed that the volume of neointima produced could be quantified, providing a quantitative comparison of stents. In contrast, the differences between stents based on stress alone were highly dependent on the choice of comparison criteria. These results show that the choice of stress criteria for stent comparisons is critical. This study shows that mechanobiological modelling may provide a valuable tool in stent design, allowing predictions of their long-term efficacy. The level of inflammation was shown to affect the sensitivity of the model to stent design. If this finding was verified in patients, this could suggest that high-inflammation patients may require alternative treatments to stenting.

Case study of the uniaxial anisotropic spin-1 bilinear-biquadratic Heisenberg model on a triangular lattice

We study the spin-1 model on a triangular lattice in the presence of a uniaxial anisotropy field using a cluster mean-field (CMF) approach. The interplay among antiferromagnetic exchange, lattice geometry, and anisotropy forces Gutzwiller mean-field approaches to fail in a certain region of the phase diagram. There, the CMF method yields two supersolid phases compatible with those present in the spin-1/2 XXZ model onto which the spin-1 system maps. Between these two supersolid phases, the three-sublattice order is broken and the results of the CMF approach depend heavily on the geometry and size of the cluster. We discuss the possible presence of a spin liquid in this region.

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.

Model calculations for copper clusters on Au(111) surfaces

Using the semi-empirical embedded-atom method, the structure of small copper clusters on Au(111) surfaces has been investigated both by static and dynamic calculations. By varying the size of roughly circular clusters, the edge energy per atom is obtained; it agrees quite well with estimates based on experimental results. Small three-dimensional clusters tend to have the shape of a pyramid, whose sides are oriented in the directions of small surface energy. The presence of a cluster is found to distort the underlying lattice of adsorbed copper atoms. (C) 2002 Published by Elsevier Science B.V.

Complete nonlinear theory of longitudinal-to-transverse dust lattice mode coupling in a single-layer dusty plasma crystal

A comprehensive nonlinear model is put forward for coupled longitudinal to transverse displacements in a horizontal dust mono-layer, levitated under the combined influence of gravity and an electric and/or magnetic sheath field. A set of coupled nonlinear evolution equations are obtained in a discrete description, and a pair of coupled (Boussinesq-like) PDEs are obtained in the continuum approximation. Finally, the amplitude modulation of the coupled modes is discussed, pointing out the importance of the coupling. All these results are generic, i.e. valid for any assumed form of the inter-grain interaction potential U and the sheath potential Phi.

Mass media and polarisation processes in the bounded confidence model of opinion dynamics

This paper presents a social simulation in which we add an additional layer of mass media communication to the social network 'bounded confidence' model of Deffuant et al (2000). A population of agents on a lattice with continuous opinions and bounded confidence adjust their opinions on the basis of binary social network interactions between neighbours or communication with a fixed opinion. There are two mechanisms for interaction. 'Social interaction' occurs between neighbours on a lattice and 'mass communication' adjusts opinions based on an agent interacting with a fixed opinion. Two new variables are added, polarisation: the degree to which two mass media opinions differ, and broadcast ratio: the number of social interactions for each mass media communication. Four dynamical regimes are observed, fragmented, double extreme convergence, a state of persistent opinion exchange leading to single extreme convergence and a disordered state. Double extreme convergence is found where agents are less willing to change opinion and mass media communications are common or where there is moderate willingness to change opinion and a high frequency of mass media communications. Single extreme convergence is found where there is moderate willingness to change opinion and a lower frequency of mass media communication. A period of persistent opinion exchange precedes single extreme convergence, it is characterized by the formation of two opposing groups of opinion separated by a gradient of opinion exchange. With even very low frequencies of mass media communications this results in a move to central opinions followed by a global drift to one extreme as one of the opposing groups of opinion dominates. A similar pattern of findings is observed for Neumann and Moore neighbourhoods.

Detection of Entanglement in Ultracold Lattice Gases

Probing non trivial magnetic ordering in quantum magnets realized with ultracold lattice gases demands detection methods with some spatial resolution built on it. Here we demonstrate that the Faraday matter-light interface provides an experimentally feasible tool to distinguish indubitably different quantum phases of a given many-body system in a non-demolishing way. We illustrate our approach by focussing on the Heisenberg chain for spin-1 bosons in the presence of a SU(2) symmetry breaking field. We explain how using the light signal obtained via homodyne detection one can reconstruct the phase diagram of the model. Further we show that the very same technique that provides a direct experimentally measurable signal of different order parameters can be extended to detect also the presence of multipartite entanglement in such systems.

Ion-acoustic supersolitons in plasmas with two-temperature electrons: Boltzmann and kappa distributions

Acoustic supersolitons arise when a plasma model is able to support three consecutive local extrema of the Sagdeev pseudopotential between the undisturbed conditions and an accessible root. This leads to a characteristic electric field signature, where a simple bipolar shape is enriched by subsidiary maxima. Large-amplitude nonlinear acoustic modes are investigated, using a pseudopotential approach, for plasmas containing two-temperature electrons having Boltzmann or kappa distributions, in the presence of cold fluid ions. The existence domains for positive supersolitons are derived in a methodological way, both for structure velocities and amplitudes, in terms of plasma compositional parameters. In addition, typical pseudopotentials, soliton, and electric field profiles have been given to illustrate that positive supersolitons can be found in the whole range of electron distributions from Maxwellian to a very hard nonthermal spectrum in kappa. However, it is found that the parameter ranges that support supersolitons vary significantly over the wide range of kappa considered. VC 2013 AIP Publishing LLC. [http://dx.doi.org/10.1063/1.4818888]

Low-energy behavior of strongly interacting bosons on a flat-band lattice above the critical filling factor

Bosons interacting repulsively on a lattice with a flat lowest band energy dispersion may, at sufficiently small filling factors, enter into a Wigner-crystal-like phase. This phase is a consequence of the dispersionless nature of the system, which in turn implies the occurrence of single-particle localized eigenstates. We investigate one of these systems-the sawtooth lattice-filled with strongly repulsive bosons at filling factors infinitesimally above the critical point where the crystal phase is no longer the ground state. We find, in the hard-core limit, that the crystal retains its structure in all but one of its cells, where it is broken. The broken cell corresponds to an exotic kind of repulsively bound state, which becomes delocalized. We investigate the excitation spectrum of the system analytically and find that the bound state behaves as a single particle hopping on an effective lattice with reduced periodicity, and is therefore gapless. Thus, the addition of a single particle to a flat-band system at critical filling is found to be enough to make kinetic behavior manifest.

Thermometry precision in strongly correlated ultracold lattice gases

The precise knowledge of the temperature of an ultracold lattice gas simulating a strongly correlated
system is a question of both fundamental and technological importance. Here, we address such
question by combining tools from quantum metrology together with the study of the quantum
correlations embedded in the system at finite temperatures. Within this frame we examine the spin-
1 2 XY chain, first estimating, by means of the quantum Fisher information, the lowest attainable
bound on the temperature precision. We then address the estimation of the temperature of the sample
from the analysis of correlations using a quantum non demolishing Faraday spectroscopy method.
Remarkably, our results show that the collective quantum correlations can become optimal
observables to accurately estimate the temperature of our model in a given range of temperatures.