Numerical simulation of a bubble rising in shear-thinning fluids

The motion of a single bubble rising freely in quiescent non-Newtonian viscous fluids was investigated experimentally and computationally. The non-Newtonian effects in the flow of viscous inelastic fluids are modeled by the Carreau theological model. An improved level set approach for computing the incompressible two-phase flow with deformable free interface is used. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The simulation results demonstrate that the algorithm is robust for shear-thinning liquids with large density (rho(1)/rho(g) up to 10(3)) and high viscosity (eta(1)/eta(g) up to 10(4)). The comparison of the experimental measurements of terminal bubble shape and velocity with the computational results is satisfactory. It is shown that the local change in viscosity around a bubble greatly depends on the bubble shape and the zero-shear viscosity of non-Newtonian shear-thinning liquids. The shear-rate distribution and velocity fields are used to elucidate the formation of a region of large viscosity at the rear of a bubble as a result of the rather stagnant flow behind the bubble. The numerical results provide the basis for further investigations, such as the numerical simulation of viscoelastic fluids. (C) 2010 Elsevier B.V. All rights reserved.

Direct Numerical Simulation of Shock/Turbulent Boundary Layer Interaction in a Supersonic Compression Ramp

A direct numerical simulation of the shock/turbulent boundary layer interaction flow in a supersonic 24-degree compression ramp is conducted with the free stream Mach number 2.9. The blow-and-suction disturbance in the upstream wall boundary is used to trigger the transition. Both the mean wall pressure and the velocity profiles agree with those of the experimental data, which validates the simulation. The turbulent kinetic energy budget in the separation region is analyzed. Results show that the turbulent production term increases fast in the separation region, while the turbulent dissipation term reaches its peak in the near-wall region. The turbulent transport term contributes to the balance of the turbulent conduction and turbulent dissipation. Based on the analysis of instantaneous pressure in the downstream region of the mean shock and that in the separation bubble, the authors suggest that the low frequency oscillation of the shock is not caused by the upstream turbulent disturbance, but rather the instability of separation bubble.

A method on theoretical simulation of chromosome breaks in cells exposed to heavy ions

Background. The aim of this study is to assess an easy and quick method on simulating chromosome breaks in cells exposed to heavy charged particles. Methods. The theoretical value of chromosome break was calculated, and the validated comparison with the experimental value by using a premature chromosome condensation technique was done. Results. A good consistence was found to be appeared between the theoretical and experimental value. Conclusions. This suggested that a higher relative biological effectiveness of heavy ions was closely correlated with its physical characteristics and besides, a safe approach on predicting chromosome breaks in cells exposed to heavy ions at off-line environment come to be considered. Furthermore, three key factors influencing the theoretical simulation was investigated and discussed.

Monte Carlo simulation of the neutron detection wall at RIBLL II in HIRFL-CSR

With the construction of the new Radioactive Ion Beam Line in Lanzhou (RIBLL II) which connects the CSRm and the CSRe, an experimental setup for physics research is highly required. A large area neutron detection wall is the main part of the setup. This paper introduced the detection principle of the neutron detection wall and the Monte-Carlo simulation of its design under the environment of the Geant4 toolkit. We presented the final design with the optimized parameters and the performance of the wall.

the role of software process simulation modeling in software risk management: a systematic review

Univ SE Calif, Ctr Syst & Software Engn, ABB, Microsoft Res, IEEE, ACMSIGSOFT, N Carolina State Univ Comp Sci

Enhancements of the simulation method on the edge effect in resin transfer molding processes

In resin transfer molding processes, small clearances exist between the fiber preform and the mold edges, which result in a preferential resin flow in the edge channel and then disrupt the flow patterns during the mold filling stage. A mathematical model including the effect of cavity thickness on resin flow was developed for flow behavior involving the interface between an edge channel and a porous medium. According to the mathematical analysis of momentum equations in a fully developed rectangular duct and formulations of the equivalent edge permeability, comparing with three-dimensional Navier-Stokes equations, the governing equations were modified in the edge channel. The volume of fluid (VOF) method was applied to track the flow front. A simple case is numerically simulated using the modified governing equations. The effects of edge channel width and cavity thickness on flow front and inlet pressure are analyzed, and the evolution characteristics of simulated results are in agreement with the experimental results. (c) 2007 Elsevier B.V. All rights reserved

Simulation of Banded Spherulite Patterns by Coupled Logistic Map Lattice Model in Three Dimensions

Banded spherulite patterns are simulated in three dimensions by means of a Coupled Logistic map lattice model. The patterns obtained by numerical calculation are consistent with those in experiments. The simulation results also indicate that the hand spacing is decreased with the increase of parameter mu in the Logistic map and increased with the increase of the coupling parameter e for cube lattices, and increased with the increase of the thickness of the lattice for polymer film, which is quite similar to the results in some experiments. Spiral pattern in three dimensions is also shown in this paper, which helps us understand the form of banded spherulite in polymers.

Simulation Study of Aggregate Morphologies Formed by ABC Linear Triblock Copolymers in a Selective Solvent Through the Self-Consistent Field Theory

The effect of the hydrophobic properties of blocks B and C on the aggregate morphologies formed by ABC linear triblock copolymers in selective solvent was studied through the self-consistent field theory. Five typical micelles, such as core-shell-corona, hamburger-like, segmented-wormlike, were obtained by changing the hydrophobic properties of blocks B and C. The simulation results indicate that the shape and size of micelle are basically controlled by the hydrophobic degree of the middle block B, whereas the type of micelle is mainly determined by the hydrophobic degree of the end block C.

Monte Carlo simulation on symmetric ABA/AB copolymer blends in confined thin films

The effects of blend composition on morphology, order-disorder transition (ODT), and chain conformation of symmetric ABA/AB copolymer blends confined between two neutral hard walls have been investigated by lattice Monte Carlo simulation. Only lamellar structure is observed in all the simulation morphologies under thermodynamic equilibrium state, which is supported by theoretical prediction. When the composition of AB diblock copolymer (phi) increases, both lamellar spacing and the corresponding ODT temperature increase, which can be attributed to the variation of conformation distribution of the diblock and the triblock copolymer chains. In addition, both diblock and triblock copolymer, chains with bridge conformation extend dramatically in the direction parallel to the surface when the system is in ordered state. Finally, the copolymer chain conformation depends strongly on both the blend composition and the incompatibility parameter chi N.

Computer simulation of Gd(III) speciation in human interstitial fluid

The speciation and distribution of Gd(III) in human interstitial fluid was studied by computer simulation. The results show that at the background concentration, all the Gd(III) species are soluble and no precipitates appear. However as the total concentration of Gd(III) rises above 2.610 x 10(-9) mol/l the insoluble species become predominant. GdPO4 is formed first as a precipitate and then Gd-2(CO3)(3). Among soluble species, free Gd(III), [Gd(HSA)], [Gd(Ox)] and the ternary complexes of Gd(III) with citrate as the primary ligand are main species when the total concentration of Gd(III) is below 2.074 x 10(-2) mol/l. With the total concentration of Gd(III) further rising, [Gd-3(OH)(4)] begins to appear and gradually becomes a predominant species.