Studying the contact point and interface moving in a sinusoidal tube with lattice Boltzmann method

The multicomponent nonideal gas lattice Boltzmann model by Shan and Chen (S-C) is used to study the immiscible displacement in a sinusoidal tube. The movement of interface and the contact point (contact line in three-dimension) is studied. Due to the roughness of the boundary, the contact point shows "stick-slip" mechanics. The "stick-slip" effect decreases as the speed of the interface increases. For fluids that are nonwetting, the interface is almost perpendicular to the boundaries at most time, although its shapes at different position of the tube are rather different. When the tube becomes narrow, the interface turns a complex curves rather than remains simple menisci. The velocity is found to vary considerably between the neighbor nodes close to the contact point, consistent with the experimental observation that the velocity is multi-values on the contact line. Finally, the effect of three boundary conditions is discussed. The average speed is found different for different boundary conditions. The simple bounce-back rule makes the contact point move fastest. Both the simple bounce-back and the no-slip bounce-back rules are more sensitive to the roughness of the boundary in comparison with the half-way bounce-back rule. The simulation results suggest that the S-C model may be a promising tool in simulating the displacement behaviour of two immiscible fluids in complex geometry.

Lattice dynamics study of microstructures of electrorheological fluids

The lattice dynamics method is used to study the stability of the chain structures formed in electrorheological (ER) fluids. The appearance of the soft modes in the phonon dispersion of the structures indicates that the chains tend to distort and aggregate into thicker columns due to the electrostatic attractive forces and thermal generated forces between them. The results show that the stability of the chains relies on their width and the separation between them. The complete chain structures are more stable than the chains with defects. The results can be used to elucidate the densification phenomenon of the chains in the structuring process of ER fluids in the quiescent state.

Lattice Boltzmann simulation of viscous fluid systems with elastic boundaries

A lattice Boltzmann model able to simulate viscous fluid systems with elastic and movable boundaries is proposed. By introducing the virtual distribution function at the boundary, the Galilean invariance is recovered for the full system. As examples of application, the how in elastic vessels is simulated with the pressure-radius relationship similar to that of the pulmonary blood vessels. The numerical results for steady how are in good agreement with the analytical prediction, while the simulation results for pulsative how agree with the experimental observation of the aortic flows qualitatively. The approach has potential application in the study of the complex fluid systems such as the suspension system as well as the arterial blood flow.

Reply to “Comment on ‘Temperature-dependent orientational ordering on a spherical surface modeled with a lattice spin model’ ”

A reply to the comment of S. Romano, Phys. Rev. E 2015 on our previous paper is provided.

Chaotic synchronization in 2D lattice for scene segmentation

Synchronization and chaos play important roles in neural activities and have been applied in oscillatory correlation modeling for scene and data analysis. Although it is an extensively studied topic, there are still few results regarding synchrony in locally coupled systems. In this paper we give a rigorous proof to show that large numbers of coupled chaotic oscillators with parameter mismatch in a 2D lattice can be synchronized by providing a sufficiently large coupling strength. We demonstrate how the obtained result can be applied to construct an oscillatory network for scene segmentation. (C) 2007 Elsevier B.V. All rights reserved.

Clustering and diffusion in a symplectic map lattice with non-local coupling

We study a symplectic chain with a non-local form of coupling by means of a standard map lattice where the interaction strength decreases with the lattice distance as a power-law, in Such a way that one can pass continuously from a local (nearest-neighbor) to a global (mean-field) type of coupling. We investigate the formation of map clusters, or spatially coherent structures generated by the system dynamics. Such clusters are found to be related to stickiness of chaotic phase-space trajectories near periodic island remnants, and also to the behavior of the diffusion coefficient. An approximate two-dimensional map is derived to explain some of the features of this connection. (C) 2008 Elsevier Ltd. All rights reserved.

D*D rho vertex from QCD sum rules

We calculate the form factors and the coupling constant in the D*D rho vertex in the framework of QCD sum rules. We evaluate the three-point correlation functions of the vertex considering D, rho and D* mesons off-shell. The form factors obtained are very different but give the same coupling constant: g(D*D rho) = 4.3 +/- 0.9 GeV(-1). (C) 2011 Elsevier B.V. All rights reserved.

A strongly attractive Fermi gas in an optical lattice

We study strongly attractive fermions in an optical lattice superimposed by a trapping potential. We calculate the densities of fermions and condensed bound molecules at zero temperature. There is a competition between dissociated fermions and molecules leading to a reduction of the density of fermions at the trap center. (C) 2010 Elsevier B.V. All rights reserved.

QCD sum rules for the production of the X(3872) as a mixed molecule-charmonium state in B meson decay

We use QCD sum rules to calculate the branching ratio for the production of the meson X(3872) in the decay B -> X(3872)K, assumed to be a mixture between charmonium and exotic molecular vertical bar c (q) over bar vertical bar vertical bar q (c) over bar vertical bar states with J(PC) = 1(++). We find that in a small range for the values of the mixing angle, 5 degrees <= theta <= 13 degrees, we get the branching ratio B(B -> XK) = (1.00 +/- 0.68) x 10(-5), which is in agreement with the experimental upper limit. This result is compatible with the analysis of the mass and decay width of the mode J/psi(n pi) and the radiative decay mode J/psi gamma performed in the same approach. (C) 2011 Elsevier B.V. All rights reserved.

SU(3) mass-splittings of heavy-baryons in QCD

We extract directly (for the first time) the charmed (C = 1) and bottom (B = -1) heavy-baryons (spin 1/2 and 3/2) mass-splittings due to SU(3) breaking using double ratios of QCD spectral sum rules (QSSR) in full QCD, which are less sensitive to the exact value and definition of the heavy quark mass, to the perturbative radiative corrections and to the QCD continuum contributions than the simple ratios commonly used for determining the heavy baryon masses. Noticing that most of the mass-splittings are mainly controlled by the ratio kappa <(S) over bars >/<(d) over bard > of the condensate, we extract this ratio, by allowing 1 sigma deviation from the observed masses of the Xi(c.b) and of the Omega(c). We obtain: kappa = 0.74(3), which improves the existing estimates: kappa = 0.70(10) from light hadrons. Using this value, we deduce M(Omega b) = 6078.5(27.4) MeV which agrees with the recent CDF data but disagrees by 2.4 sigma with the one from D0. Predictions of the Xi(Q)` and of the spectra of spin 3/2 baryons containing one or two strange quark are given in Table 2. Predictions of the hyperfine splittings Omega(Q)* - Omega(Q) and Xi(Q)* - Xi(Q) are also given in Table 3. Starting for a general choice of the interpolating currents for the spin 1/2 baryons, our analysis favours the optimal value of the mixing angle b similar or equal to (-1/5-0) found from light and non-strange heavy baryons. (C) 2010 Elsevier B.V. All rights reserved.

A QCD sum rule calculation for the Y(4140) narrow structure

We use the QCD sum rules to evaluate the mass of a possible scalar mesonic state that couples to a molecular D(s)*(D) over bar (s)* current. We find a mass m(Ds)*(Ds)* = (4.14 +/- 0.09) GeV, which is in an excellent agreement with the recently observed Y(4140) charmonium state. We consider the contributions of condensates up to dimension-eight, we work at leading order in alpha(s) and we keep terms which are linear in the strange quark mass m(s). We also consider a molecular D*(D) over bar* current and we obtain m m(D)*(D)* = (4.13 +/- 0.10), around 200 MeV above the mass of the Y(3930) charmonium state. We conclude that it is possible to describe the Y(4140) structure as a D(s)*(D) over bar (s)* molecular state or even as a mixture of D(s)*(D) over bar (s)* and D*(D) over bar* molecular states. (C) 2009 Elsevier B.V. All rights reserved.

QCD sum rule calculation for the charmonium-like structures in the J/psi phi and J/psi omega invariant mass spectra

Using the QCD sum rules we test if the charmonium-like structure Y(4274), observed in the J/psi phi invariant mass spectrum, can be described with a D(s)(D) over bar (s0)(2317)+ h.c. molecular current with J(PC) = 0(-+). We consider the contributions of condensates up to dimension ten and we work at leading order in alpha(s). We keep terms which are linear in the strange quark mass m(s). The mass obtained for such state is mD(s)D(s0) = (4.78 +/- 0.54) GeV. We also consider a molecular 0(-+) D (D) over bar (0)(2400)+ h.c. current and we obtain m(DD0) = (4.55 +/- 0.49) GeV. Our study shows that the newly observed Y(4274) in the J/psi phi invariant mass spectrum can be, considering the uncertainties, described using a molecular charmonium current. (C) 2011 Elsevier B.V. All rights reserved.

QCD sum rules study of the J(PC)=1(- -) charmonium Y mesons

We use QCD sum rules to test the nature of the recently observed mesons Y(4260), Y(4350) and Y(4660), assumed to be exotic four-quark (c (c) over barq (q) over bar) or (c (c) over bars (s) over bar) states with J(PC)= 1(--). We work at leading order in alpha(s), consider the contributions of higher dimension condensates and keep terms which are linear in the strange quark mass m(s). We find for the (c (c) over bars (s) over bar) state a mass in m(Y) = (4.65 +/- 0.10) GeV which is compatible with the experimental candidate Y (4660), while for the (c (c) over barq (q) over bar) state we find a mass in m(Y) = (4.49 +/- 0.11) GeV, which is still consistent with the mass of the experimental candidate Y(4350). With the tetraquark structure we are working we cannot explain the Y(4260) as a tetraquark state. We also consider molecular D(s0)(D) over bar (s)* and D(0)(D) over bar* states. For the D(s0)(D) over bar (s)* molecular state we get m(Ds0 (D) over bars*) = (4.42 +/- 0.10) GeV which is consistent, considering the errors, with the mass of the meson Y(4350) and for the D(0)(D) over bar* molecular state we get m(D0 (D) over bar*) = (4.27 +/- 0.10) GeV in excellent agreement with the mass of the meson Y(4260). (C) 2008 Elsevier B.V. All rights reserved.

QCD sum rules study of the meson Z(+)(4430)

We use QCD sum rules to study the recently observed meson Z(+)(4430), considered as a D*D-1 molecule with J(P) = 0(-). We consider the contributions of condensates up to dimension eight and work at leading order in alpha(s). We get m(Z) = (4.40 +/- 0.10) GeV in a very good agreement with the experimental value. We also make predictions for the analogous mesons Z(s) and Z(bb) considered as D-s*D-1 and B*B-1 molecules, respectively. For Z(s) we predict mZ(s) = (4.70 +/- 0.06) GeV, which is above the D-s* D-1 threshold, indicating that it is probably a very broad state and, therefore, difficult to observe experimentally. For Z(bb) we predict m(Zbb) = (10.74 +/- 0.12) GeV, in agreement with quark model predictions. (c) 2008 Elsevier B.V. All rights reserved.

rho D*D* vertex from QCD sum rules

We calculate the form factors and the coupling constant in the rho D*D* vertex in the framework of QCD sum rules. We evaluate the three point correlation functions of the vertex considering both rho and D* mesons off-shell. The form factors obtained are very different but give the same coupling constant: g rho D*D* = 6.60 +/- 0.31. This number is 50% larger than what we would expect from SU(4) estimates. (c) 2007 Elsevier B.V. All rights reserved.