Acoustic wave propagation simulation in a poroelastic medium saturated by two immiscible fluids using a staggered finite-difference with a time partition method

Based on the three-phase theory proposed by Santos, acoustic wave propagation in a poroelastic medium saturated by two immiscible fluids was simulated using a staggered high-order finite-difference algorithm with a time partition method, which is firstly applied to such a three-phase medium. The partition method was used to solve the stiffness problem of the differential equations in the three-phase theory. Considering the effects of capillary pressure, reference pressure and coupling drag of two fluids in pores, three compressional waves and one shear wave predicted by Santos have been correctly simulated. Influences of the parameters, porosity, permeability and gas saturation on the velocities and amplitude of three compressional waves were discussed in detail. Also, a perfectly matched layer (PML) absorbing boundary condition was firstly implemented in the three-phase equations with a staggered-grid high-order finite-difference. Comparisons between the proposed PML method and a commonly used damping method were made to validate the efficiency of the proposed boundary absorption scheme. It was shown that the PML works more efficiently than the damping method in this complex medium. Additionally, the three-phase theory is reduced to the Biot’s theory when there is only one fluid left in the pores, which is shown in Appendix. This reduction makes clear that three-phase equation systems are identical to the typical Biot’s equations if the fluid saturation for either of the two fluids in the pores approaches to zero. Supported by the Key Program of the National Natural Science Foundation of China (Grant No. 10534040) and the National Natural Science Foundation of China (Grant No. 10674148)     

Patch test function for axisymmetric element of conventional and couple stress theory

The enhanced patch test proposed by Chen W J (2006) can be used to assess the convergence of the problem with non-homogeneous differential equations. Based on this theory, we establish the patch test function for axisymmetric elements of conventional and couple stress theories, and reach an important conclusion that the patch test function for axisymmetric elements cannot contain non-zero constant shear. Supported by the National Natural Science Foundation of China (Grant No. 10672032)     

The long time and Brownian limits for tagged particle motion in liquids

Formally exact theories of tagged particle motion in liquids are developed, based upon kinetic theory for hard spheres and mode coupling for smooth potentials. It is shown that the resulting equations are tractable in the long time and Brownian limits. The coefficient of the long time tail of the velocity correlation function is seen to differ from its low-density form by only the replacement of the low-density viscosity and diffusion constant by the true viscosity and diffusion constant. In the Brownian limit, the slip Stokes-Einstein law is obtained, with the true viscosity. The relation to other work is discussed.Supported by NSF Grant No. CHE81-11422 and by a Dreyfus Teacher-Scholar grant to TK.     

Approximate homotopy symmetry method: Homotopy series solutions to the sixth-order Boussinesq equation

An approximate homotopy symmetry method for nonlinear problems is proposed and applied to the sixth-order Boussinesq equation,which arises from fluid dynamics.We summarize the general formulas for similarity reduction solutions and similarity reduction equations of different orders,educing the related homotopy series solutions.Zero-order similarity reduction equations are equivalent to the Painlevé IV type equation or Weierstrass elliptic equation.Higher order similarity solutions can be obtained by solving...     

Compositional Dependence of Static Shear Viscosity of Immiscible PP/PS Blends

Phase structures of immiscible polypropylene (PP)/polystyrene (PS) blends with different volume proportions, PP90/PS10, PP80/PS20, PP70/PS30, PP60/PS40, PP50/PS50, PP40/PS60, PP30/PS70, PP20/PS80, PP10/PS90, were observed by means of scanning electronic microscopy (SEM). The zero shear viscosities of the blends were determined according to a modified Carreau model by fitting the curves of static shear rate sweeps of blends tested at 190°C in a Stress Tech Fluids Rheometer. The results showed that the compositional dependence of zero shear viscosity of PP/PS deviated greatly from linear or log‐linear additivity. When PS was dispersed in a PP continuous phase, the blends showed negative deviation, while for blends with PP dispersed in a PS matrix, positive deviation was generated. When different theoretical equations of Nielsen, Utracki, Taylor, Frankel‐Acrivos (FA), Choi‐Schowalter (CS), and Han‐King (HK) were used to fit the experimental data of zero shear viscosities of blends, none of them was suitable for PP/PS blends. These experimental phenomena may result from the complex phase structures of the blends and their response to shear conditions, which are discussed in detail and compared with the experimental analysis.     

General Zakharov-Shabat equations,multi-time Hamiltonian formalism,and constants of motion

We construct a Hamiltonian formalism for general Zakharov-Shabat equations (zero curvature equations with rational dependence on a parameter) as well as their constants of motion, and prove that the latter are in involution. The field-theoretical (multi-time) Hamiltonian formalism is used.Supported by National Science Foundation, Research Grant DMS 8703407.     

Modified Single-Mode Leonov Rheological Equations for Polymer Melts and Solutions

The shear viscosity and the normal stress coefficients are important parameters in the flow of polymer melts and polymer solutions. Based on the Leonov model, modified single-mode rheological equations are presented by introducing relaxation time and temperature functions, and the shear viscosity and the normal stress coefficients are predicted. Without a complex statistical simulation, the experimental data of a low-density polyethylene melt, a poly(ethylene oxide) solution and a mixed decalin/polybutene oil solution were compared to verify the modified equations in very wide range of deformation rates. Furthermore, based on the equations, the relationship between the stress overshoot and the temperature is discussed. In addition, the predicted shear thinning behavior for the modified equations is also compared with other single-mode models.     

Two-dimensional coupled mathematical modeling of fluvial processes with intense sediment transport and rapid bed evolution

Alluvial rivers may experience intense sediment transport and rapid bed evolution under a high flow regime, for which traditional decoupled mathematical river models based on simplified conservation equations are not applicable. A two-dimensional coupled mathematical model is presented, which is generally applicable to the fluvial processes with either intense or weak sediment transport. The governing equations of the model comprise the complete shallow water hydrodynamic equations closed with Manning roughness for boundary resistance and empirical relationships for sediment exchange with the erodible bed. The second-order Total-Variation-Diminishing version of the Weighted-Average-Flux method, along with the HLLC approximate Riemann Solver, is adapted to solve the governing equations, which can properly resolve shock waves and contact discontinuities. The model is applied to the pilot study of the flooding due to a sudden outburst of a real glacial-lake. Supported by the National Basic Research and Development Program of China (973 Program) (Grant No. 2007CB14106), the National Natural Science Foundation of China (Grant No. 50459001), and the Key Project of Chinese Academy of Sciences (Grant No. KZCX3-SW-357-02)     

Shear viscosity for inelastic hard spheres

The hydrodynamic equations of the Enskog theory for inelastic hard spheres is considered as a model for rapid flow granular fluids at finite densities. A detailed analysis of the shear viscosity of the granular fluid has been done using homogenous cooling state (HCS) and uniform shear flow (USF) models. It is found that shear viscosity is sensitive to the coefficient of restitution α and pair correlation function at contact. The collisional part of the Newtonian shear viscosity is found to be dominant than its kinetic part.     

The vanishing viscosity limit for a dyadic model

A dyadic shell model for the Navier-Stokes equations is studied in the context of turbulence. The model is an infinite nonlinearly coupled system of ODEs. It is proved that the unique fixed point is a global attractor, which converges to the global attractor of the inviscid system as viscosity goes to zero. This implies that the average dissipation rate for the viscous system converges to the anomalous dissipation rate for the inviscid system (which is positive) as viscosity goes to zero. This phenomenon is called the dissipation anomaly predicted by Kolmogorov’s theory for the actual Navier-Stokes equations.     

Steady-state responses of axially accelerating viscoelastic beams: Approximate analysis and numerical confirmation

Nonlinear parametric vibration of axially accelerating viscoelastic beams is investigated via an approximate analytical method with numerical confirmations. Based on nonlinear models of a finite-small-stretching slender beam moving at a speed with a periodic fluctuation, a solvability condition is established via the method of multiple scales for subharmonic resonance. Therefore, the amplitudes of steady-state periodic responses and their existence conditions are derived. The amplitudes of stable steady-state responses increase with the amplitude of the axial speed fluctuation, and decrease with the viscosity coefficient and the nonlinear coefficient. The minimum of the detuning parameter which causes the existence of a stable steady-state periodic response decreases with the amplitude of the axial speed fluctuation, and increases with the viscosity coefficient. Numerical solutions are sought via the finite difference scheme for a nonlinear partial-differential equation and a nonlinear integro-partial-differential equation. The calculation results qualitatively confirm the effects of the related parameters predicted by the approximate analysis on the amplitude and the existence condition of the stable steady-state periodic responses. Quantitative comparisons demonstrate that the approximate analysis results have rather high precision. Supported by the National Outstanding Young Scientists Foundation of China (Grant No. 10725209), the National Natural Science Foundation of China (Grant No. 10672092), Scientific Research Project of Shanghai Municipal Education Commission (Grant No. 07ZZ07), and Shanghai Leading Academic Discipline Project (Grant No. Y0103)     

Vorticity alignment and negative normal stresses in sheared attractive emulsions

Attractive emulsions near the colloidal glass transition are investigated by rheometry and optical microscopy under shear. We find that (i) the apparent viscosity eta drops with increasing shear rate, then remains approximately constant in a range of shear rates, then continues to decay; (ii) the first normal stress difference N1 transitions sharply from nearly zero to negative in the region of constant shear viscosity; and (iii) correspondingly, cylindrical flocs form, align along the vorticity, and undergo a log-rolling movement. An analysis of the interplay between steric constraints, attractive forces, and composition explains this behavior, which seems universal to several other complex systems.     

Transport coefficients in Chiral Perturbation Theory

We present recent results on the calculation of transport coefficients for a pion gas at zero chemical potential in Chiral Perturbation Theory (ChPT) using the Linear Response Theory (LRT). More precisely, we show the behavior of DC conductivity and shear viscosity at low temperatures. To compute transport coefficients, the standard power counting of ChPT has to be modified. The effects derived from imposing unitarity are also analyzed. As physical applications in relativistic heavy-ion collisions, we show the relation of the DC conductivity to soft-photon production and phenomenological effects related to a non-zero shear viscosity. In addition, our values for the shear viscosity to entropy ratio satisfy the KSS bound.     

On the generation of new metrics containing arbitrary functions

We present simple methods by which certain solutions of the Einstein vacuum equation (and some other equations) can be used to generate further solutions of these equations. In some cases the new solutions admit a smaller number of metric automorphisms than the original ones and are, in this sense, more general.Supported by N.R.C. Grant No. A-5205.Supported by N.R.C. Grant No. A-3993.     

Comparison of shear banding in BMGs due to thermal-softening and free volume creation

This paper reports a comparative study of shear banding in BMGs resulting from thermal softening and free volume creation. Firstly, the effects of thermal softening and free volume creation on shear instability are discussed. It is known that thermal softening governs thermal shear banding, hence it is essentially energy related. However, compound free volume creation is the key factor to the other instability, though void-induced softening seems to be the counterpart of thermal softening. So, the driving force for shear instability owing to free volume creation is very different from the thermally assisted one. In particular, long wave perturbations are always unstable owing to compound free volume creation. Therefore, the shear instability resulting from coupled compound free volume creation and thermal softening may start more like that due to free volume creation. Also, the compound free volume creation implies a specific and intrinsic characteristic growth time of shear instability. Finally, the mature shear band width is governed by the corresponding diffusions (thermal or void diffusion) within the band. As a rough guide, the dimensionless numbers: Thermal softening related number B, Deborah number (denoting the relation of instability growth rate owing to compound free volume and loading time) and Lewis number (denoting the competition of different diffusions) show us their relative importance of thermal softening and free volume creation in shear banding. All these results are of particular significance in understanding the mechanism of shear banding in bulk metallic glasses (BMGs). Supported by the Chinese Academy of Sciences under the project “Multi-Scale Complex System” (Grant No. KJCX-SW-L08), the National Natural Science Foundation of China (Grant Nos. 10725211 and 10721202), and the Doctorial Start-up Fund of Hunan University of Science and Technology (Grant No. E50840)     

Negative electrorheological responses of micro-polar fluids in the finite spin viscosity small spin velocity limit. I. Couette flow geometries

Negative electrorheological responses induced by micro-particle electrorotation in two-dimensional Couette flow geometries are analyzed by a set of continuum modeling field equations originating from anti-symmetric/couple stress theories in the finite spin viscosity small spin velocity (FSV) limit. Analytical solutions are obtained for the first time to express the spin velocity, linear velocity, and effective viscosity in terms of the electric field strength, driving shear rate, boundary condition selection parameter, and spin viscosity. Good agreement is achieved between the FSV theoretical predictions presented herein and the experimental measurements reported in recent literature for the effective viscosity for low driving shear rates.     

Novel Rotating Hairy Black Hole in (2+1)-Dimensions and Shear Viscosity to Entropy Ratio

The novel rotating hairy black hole metric in (2 + 1) dimensions, which is an exact solution to the field equations of the Einstein-scalar AdS theory with a non-minimal coupling, considered in this paper and some hydrodynamics quantities such as diffusion constant and shear viscosity investigated. By using thermodynamics quantities such as temperature and entropy we can use diffusion constant to obtain shear viscosity and then calculate shear viscosity to entropy ratio.     

The refined theory of deep rectangular beams for symmetrical deformation

Based on elasticity theory, various one-dimensional equations for symmetrical deformation have been deduced systematically and directly from the two-dimensional theory of deep rectangular beams by using the Papkovich-Neuber solution and the Lur’e method without ad hoc assumptions, and they construct the refined theory of beams for symmetrical deformation. It is shown that the displacements and stresses of the beam can be represented by the transverse normal strain and displacement of the mid-plane. In the case of homogeneous boundary conditions, the exact solutions for the beam are derived, and the exact equations consist of two governing differential equations: the second-order equation and the transcendental equation. In the case of non-homogeneous boundary conditions, the approximate governing differential equations and solutions for the beam under normal loadings only and shear loadings only are derived directly from the refined beam theory, respectively, and the correctness of the stress assumptions in classic extension or compression problems is revised. Meanwhile, as an example, explicit expressions of analytical solutions are obtained for beams subjected to an exponentially distributed load along the length of beams. Supported by the National Natural Science Foundation of China (Grant Nos. 10702077, 10672001, and 10602001), the Beijing Natural Science Foundation (Grant No. 1083012), and the Alexander von Humboldt Foundation in Germany     

The theoretical research of basic function method in incompressible viscous flow and its simulations in three-dimensional aneurysms

Basic function method is developed to treat the incompressible viscous flow. Artificial compressibility coefficient, the technique of flux splitting method and the combination of central and upwind schemes are applied to construct the basic function scheme of trigonometric function type for solving three-dimensional incompressible Navier-Stokes equations numerically. To prove the method, flows in finite-length-pipe are calculated, the velocity and pressure distribution of which solved by our method quite coincide with the exact solutions of Poiseuille flow except in the areas of entrance and exit. After the method is proved elementary, the hemodynamics in two- and three-dimensional aneurysms is researched numerically by using the basic function method of trigonometric function type and unstructured grids generation technique. The distributions of velocity, pressure and shear force in steady flow of aneurysms are calculated, and the influence of the shape of the aneurysms on the hemodynamics is studied. Supported by the National Natural Foundation of China (Grant Nos. 40874077, 40504020, and 40536029) and the National Basic Research Program of China (Grant No. 2006CB806304)     

Pair function for the rectangular Ising ferromagnet

A representation of the pair correlation function for the rectangular Ising model in zero magnetic field is derived using a new spinor technique; this enables the scaling limit to be established, as well as several analytical properties of the scaling functions.Supported by the Fonds National Suisse de la Recherche Scientifique, by the National Science Foundation Grant No. PHY76-17191, and by the National Research Council of Canada Grant No. NRC A9344