Solution of Müller's paradox of frame dependence in mixtures of polymeric liquid crystals

Recently, there appeared in this journal a short review article by Müller ((2002), 14: 227–229), in which it was argued that the internal energy of a reacting mixture of liquid crystals should not be an objective quantity (i.e., a quantity independent of the referential frame). Such a paradoxical conclusion has revealed the urging necessity for a better comprehension of the interactions taking place in structured mixtures, specially when referred to non-inertial observers. This work shows that Müller's paradox is avoided when all inertial effects are carefully accounted for. Further, it predicts interesting phenomena without analogue in reacting mixtures of structureless fluids: internal inertial effects, produced by a combination of mass exchanges (e.g. by chemical reactions or phase changes) with the extra degrees of freedom posed by the microstructure. Such effects clarify the reasons why inertial forces and couples acting on a mixture do not always coincide with the sum of the inertial forces and couples exerted upon its individual constituents. The present conclusions establish a reinterpretation of some fundamental concepts of continuum mechanics and thermodynamics, including a deeper understanding of the manner in which total energy and momenta are conserved in complex media. Received: July 25, 2002 / Published online: February 17, 2003 Communicated by Kolumban Hutter, Darmstadt e-mail: faria@mechanik.tu-darmstadt.de     

On the pressure and stress singularities induced by steady flows of incompressible viscous fluids

Design for structural integrity requires an appreciation of where stress singularities can occur in structural configurations. While there is a rich literature devoted to the identification of such singular behavior in solid mechanics, to date there has been relatively little explicit identification of stress singularities caused by fluid flows. In this study, stress and pressure singularities induced by steady flows of viscous incompressible fluids are asymptotically identified. This is done by taking advantage of an earlier result that the Navier-Stokes equations are locally governed by Stokes flow in angular corners. Findings for power singularities are confirmed by developing and using an analogy with solid mechanics. This analogy also facilitates the identification of flow-induced log singularities. Both types of singularity are further confirmed for two global configurations by applying convergence-divergence checks to numerical results. Even though these flow-induced stress singularities are analogous to singularities in solid mechanics, they nonetheless render a number of structural configurations singular that were not previously appreciated as such from identifications within solid mechanics alone.     

The interactions between wave-currents and offshore structures with consideration of fluid viscosity

Study of the flow field around the large scale offshore structures under the action of waves and viscous currents is of primary importance for the scouring estimation and protection in the vicinity of the structures. But very little has been known in its mechanism when the viscous effects is taken into consideration. As a part of the efforts to tackle the problem, a numerical model is presented for the simulation of the flow field around a fixed vertical truncated circular cylinder subjected to waves and viscous currents based on the depth-averaged Reynolds equations and depth-averagedk-ɛ turbulence model. Finite difference method with a suitable iteration defect correct method and an artificial open boundary condition are adopted in the numerical process. Numerical results presented relate to the interactions of a pure incident viscous current with Reynolds numberRe=10⁵, a pure incident regular sinusoidal wave, and the coexisting of viscous current and wave with a circular cylinder, respectively. Flow fields associated with the hydrodynamic coefficients of the fixed cylinder, as well as corresponding free surface profiles and wave amplitudes, are discussed. The present method is found to be relatively straightforward, computationally effective and numerically stable for treating the problem of interactions among waves, viscous currents and bodies. The project supported by the National Natural Science Foundation of China and Foundation of State Key Laboratory of Ocean Engineering at Shanghai Jiao Tong University.     

On the Chaotic Dynamics of a Spherical Pendulum with a Harmonically Vibrating Suspension

The equations of motion for a lightly damped spherical pendulum are considered. The suspension point is harmonically excited in both vertical and horizontal directions. The equations are approximated in the neighborhood of resonance by including the third order terms in the amplitude. The stability of equilibrium points of the modulation equations in a four-dimensional space is studied. The periodic orbits of the spherical pendulum without base excitations are revisited via the Jacobian elliptic integral to highlight the role played by homoclinic orbits. The homoclinic intersections of the stable and unstable manifolds of the perturbed spherical pendulum are investigated. The physical parameters leading to chaotic solutions in terms of the spherical angles are derived from the vanishing Melnikov–Holmes–Marsden (MHM) integral. The existence of real zeros of the MHM integral implies the possible chaotic motion of the harmonically forced spherical pendulum as a result from the transverse intersection between the stable and unstable manifolds of the weakly disturbed spherical pendulum within the regions of investigated parameters. The chaotic motion of the modulation equations is simulated via the 4th-order Runge–Kutta algorithms for certain cases to verify the analysis.     

Modeling and Algorithm on a Class of Mechanical Systems with Unilateral Constraints

In the present paper, the cable-structures are considered as a class of mechanical complementary-slackness systems. Based on the optimization algorithms used for multi-body dynamics with unilateral contacts, an algorithm by means of artificial neural network (NNW) is developed. The following two classes of cable-structures have been considered force-elongation of cable member follows elastic behavior and work-hardening assumption. Due to simplicity the former is used to prove the method reliability, and the latter, as general cable-structure problem is handled. First, the complementarity problems for those structures have been formulated; then using generalized Gaussian‘ least action principle they are summarized as an optimization problem. Based on Hopfield’s work, an artificial NNW has been designed and used to decide combination of possible constraints at each step in simulation. As examples, two cable-structures have been investigated. The calculated results for a simple suspension structure evidence the reliability and time-economization of the proposed method. An example of guyed mast shows the suitability of the proposed method for practical cable-structures.This paper supported by NSFC, No.10072040     

Solutions for a system of nonlinear random integral and differential equations under weak topology

I.IntroductionTheexistenceandcomparisonresultsofsolutionsfornonlinearVolterraintegralequationsinstrongtopologyofBanachspaceshavebeenobtainedbyVaughnI"'"],LakshmikanthamIl3]andLakshmikantham-Leela114l.TheexistenceresultsofweaksolutionsfortheCauchyproblemof…     

A coupling method of difference with high order accuracy and boundary integral equation for evolutionary equation and its error estimates

In the present paper,a new numerical method for solving initial-boundary valueproblems of evolutionary equations is proposed and studied,combining difference methodwith high accuracy with boundary integral equation method.The numerical approximateschemes for both problems on a bounded or unbounded domain in R~3 are proposed and theirprior error estimates are obtained.     

The effect of NaCl addition on the rheological behavior of cetyltrimethylammonium <Emphasis Type="Italic">p</Emphasis>-toluenesulfonate (CTAT) aqueous solutions and their mixtures with hydrophobically modified polyacrylamide aqueous solutions

In this paper, the influence of NaCl addition, up to very large concentrations, on the rheological properties of cetyltrimethylammonium p-toluenesulfonate (CTAT) solutions and their mixtures with two hydrophobically modified polyacrylamides (HMPAM) has been studied under simple shear. The CTAT concentrations employed were above the critical rod concentration. As salt is added to CTAT aqueous solutions, the zero-shear viscosity first increases, goes through a maximum, and at very high ionic strengths increases once more. The overlap concentration of worm-like micelles decreases as the concentration of NaCl increases. The results are explained by the salt addition-induced growth of worm-like micelles and salting out effects at the highest contents of NaCl. The influence of ionic environment on the rheological properties of CTAT with two HMPAM solutions with different contents of hydrophobic moieties was also studied under simple shear. When NaCl is added to HMPAM/CTAT solutions, the same trends observed in CTAT/NaCl solutions were repeated but the viscosity increases were largely magnified. The large viscosity enhancements with salt increments in HMPAM/CTAT solutions were explained by the formation of an interpenetrated network of hydrophobically modified polymer chains and worm-like micelles with hydrophobic sequences embedded within its structure.     

Thermal instability of a porous medium with relaxation and inertia in the presence of Hall effects

Summary  The thermal instability of a Rivlin–Ericksen fluid in a porous medium is considered in the presence of a uniform vertical magnetic field to include the effect of Hall currents. For the case of stationary convection, the magnetic field has a stabilizing effect on the system, whereas the Hall current has a destabilizing effect on the system. The medium permeability has both stabilizing and destabilizing effects, depending on the Hall parameter M. The kinematic viscoelasticity has no effect on stationary convection. Graphs have been plotted by giving numerical values to the parameters, to depict the stability characteristics. The magnetic field (and corresponding Hall currents) introduces oscillatory modes in the system, which would be nonexistent in their absence. The sufficient conditions for the nonexistence of overstability are also obtained. Received 20 May 1999; accepted for publication 8 March 2000     

Homogenization of scalar wave equations with hysteresis

The paper deals with a scalar wave equation of the form where is a Prandtl–Ishlinskii operator and are given functions. This equation describes longitudinal vibrations of an elastoplastic rod. The mass density and the Prandtl–Ishlinskii distribution function are allowed to depend on the space variable x. We prove existence, uniqueness and regularity of solution to a corresponding initial-boundary value problem. The system is then homogenized by considering a sequence of equations of the above type with spatially periodic data and , where the spatial period tends to 0. We identify the homogenized limits and and prove the convergence of solutions to the solution of the homogenized equation. Received June 17, 1999     

Transient analysis of piles and pile groups in non-homogeneous soil

In this paper, the transient analysis of a single pile and a 3 × 3 pile group is presented for Gibson type non-homogeneous soil by using a hybrid type of boundary and finite element formulation for the soil domain and pile domain, respectively. The formula is presented for a transient point force acting in the interior of a non-homogeneous, isotropic half space. A time stepping boundary element algorithm for soil domain is used together with an implicit time integration scheme for finite pile domain. To investigate the validity of this formulation, a single pile and a pile group are analyzed under Heaviside loading and triangular transient loading. In the analyses, it can be concluded that the results agree well for all cases of the inhomogeneity index by comparing the Laplace domain solutions.     

Viscous and inviscid interactions of an oblique shock with a flexible panel

The complex self-sustained oscillations arising from the interaction of an oblique shock with a flexible panel in both the inviscid and viscous regimes have been investigated numerically. The aeroelastic interactions are simulated using either the Euler or the full compressible Navier–Stokes equations coupled to the nonlinear von Karman plate equations. Results demonstrate that for a sufficiently strong shock limit-cycle oscillations emerge from either subcritical or supercritical bifurcations even in the absence of viscous separated flow effects. The critical dynamic pressure diminishes with increasing shock strength and can be much lower than that corresponding to standard panel flutter. Significant changes in panel dynamics were also found as a function of the shock impingement point and cavity pressure. For viscous laminar flow above the panel without a shock, high-frequency periodic oscillations appear due to the coupling of boundary-layer instabilities with high-mode flexural deflections. For a separated shock laminar boundary layer interaction, non-periodic self-excited oscillations arise which can result in a significant reduction in the extent of the time-averaged separation region. This finding suggests the potential use of an aeroelastically tailored flexible panel as a means of passive flow control. Forced panel oscillations, induced by a specified variable cavity pressure underneath the panel, were also found to be effective in reducing separation. For both inviscid and viscous interactions, the significant unsteadiness generated by the fluttering panel propagates along the complex reflected expansion/recompression wave system.     

A note on probabilistic contractor couple and solutions for a system of nonlinear equations in n a menger pn-spaces

ＩｎｔｒｏｄｕｃｔｉｏｎａｎｄＤｅｆｉｎｉｔｉｏｎｓＩｎ [1 ] ,ＣｈａｎｇＳｈｉｓｈｅｎｇａｎｄＰｅｎｇＹｏｎｇｃｈｅｎｇｉｎｔｒｏｄｕｃｅｄｔｈｅｃｏｎｃｅｐｔｏｆｐｒｏｂａｂｉｌｉｓｔｉｃｃｏｎｔｒａｃｔｏｒｃｏｕｐｌｅｉｎＮｏｎ＿ＡｒｃｈｉｍｅｄｅａｎＰＮ＿ｓｐａｃｅｓ (ｆｏｒｓｈｏｒｔ,Ｎ .Ａ .ＭｅｎｇｅｒＰＮ＿ｓｐａｃｅｓ) .Ｕｓｉｎｇｔｈｉｓｃｏｎｃｅｐｔ,ｕｎｄｅｒｔｈｅｃｏｎｄｉｔｉｏｎｔ＿ｎｏｒｍΔ =ｍｉｎ ,ｔｈｅｙｓｔｕｄｉｅｄｔｈｅｅｘｉｓｔｅｎｃｅａｎｄｕｎｉｑｕｅｎｅｓｓｏ…     

On the use of solid- and shell-type finite elements in creep-damage predictions of thinwalled structures

Summary  Solid- and shell-type finite elements available for plasticity and creep analysis are applied to the creep-damage prediction of a thinwalled pipe bend under uniform internal pressure. Conventional creep-damage material model with scalar damage parameter is used. Based on the comparative numerical study, performed using solid and shell elements, the applicability frame of the shell concept is discussed. Particularly, if a dependence on the stress state is included in the material model, the cross-section assumptions of the first-order shear deformation theory should be refined. The possibilities to modify the through-thickness approximations are demonstrated on the beam equations. The first-order shear-deformation beam theory is discussed in detail. It is shown that if the damage evolution significantly differs for tensile and compressive stresses, the classical parabolic transverse shear-stress distribution and the shear-correction coefficient have to be modified within time-step simulations. Received 30 January 2000; accepted for publication 30 May 2000     

Efficient coupling of direct forcing immersed boundary-lattice Boltzmann method and finite element method to simulate fluid-structure interactions

Fluid-structure interactions (FSI) of rigid and flexible bodies are simulated in this article. For the fluid flow, multidirect forcing immersed boundary method (IBM) is adopted to capture the moving boundary, and lattice Boltzmann method (LBM) is used to evolve the flow field. Compared with our previous no-penetration IBM, less iterations are required in this work. In addition, larger velocity in lattice units can be used and the nonphysical force oscillations are suppressed due to the C³ 6-point kernel. Multi-relaxation-time collision operator and local grid refinement are also adopted in LBM to enhance the numerical stability and efficiency. To improve the efficiency of the FSI coupling algorithm, the mesh of the deformable structure can be coarser than the Lagrangian mesh using Newton-Cotes formulas to integrate the traction on the structure surface. A variety of benchmarks, including flow around a circular cylinder with Reynold number ranging from 20 to 200, forced oscillation of a circular cylinder, vortex-induced vibration (VIV) of an elastically mounted circular cylinder in two degrees of freedom, and VIV of an elastic cantilever beam attached to a circular cylinder, are carried out to evaluate the accuracy and stability of the present coupling algorithm. For the benchmark FSI problem considered in this article, a reduction of 54% of the calculation time is achieved using coarser structure mesh. As an application of the FSI coupling algorithm, the mechanism of an elastic beam in the wake of a circular cylinder is discussed.     

Delay-independent exponential stability of stochastic Cohen–Grossberg neural networks with time-varying delays and reaction–diffusion terms

Different from the approaches used in the earlier papers, in this paper, the Halanay inequality technique, in combination with the Lyapunov method, is exploited to establish a delay-independent sufficient condition for the exponential stability of stochastic Cohen–Grossberg neural networks with time-varying delays and reaction–diffusion terms. Moreover, for the deterministic delayed Cohen–Grossberg neural networks, with or without reaction–diffusion terms, sufficient criteria for their global exponential stability are also obtained. The proposed results improve and extend those in the earlier literature and are easier to verify. An example is also given to illustrate the correctness of our results.     

Effect of Reactive Surface Areas Associated with Different Particle Shapes on Chemical-Dissolution Front Instability in Fluid-Saturated Porous Rocks

In this article, the effect of reactive surface areas associated with different particle shapes on the reactive infiltration instability in a fluid-saturated porous medium is investigated through analytically deriving the dimensionless pore-fluid pressure-gradient of a coupled system between porosity, pore-fluid flow and reactive chemical-species transport within two idealized porous media consisting of spherical and cubic grains respectively. Compared with the critical dimensionless pore-fluid pressure-gradient of the coupled system, the derived dimensionless pore-fluid pressure-gradient can be used to assess the instability of a chemical dissolution front within the fluid-saturated porous medium. The related theoretical analysis has demonstrated that (1) since the shape coefficient of spherical grains is greater than that of cubic grains, the chemical system consisting of spherical grains is more unstable than that consisting of cubic grains, and (2) the instability likelihood of a natural porous medium, which is comprised of irregular grains, is smaller than that of an idealized porous medium, which is comprised of regular spherical grains. To simulate the complicated morphological evolution of a chemical dissolution front in the case of the chemical dissolution system becoming supercritical, a numerical procedure is proposed for solving this kind of problem. The related numerical results have demonstrated that the reactive surface area associated with different particle shapes can have a significant influence on the morphological evolution of an unstable chemical-dissolution front within fluid-saturated porous rocks.     

Modal interactions in primary and subharmonic resonant dynamics of imperfect microplates with geometric nonlinearities

This paper analyses the modal interactions in the nonlinear, size-dependent dynamics of geometrically imperfect microplates. Based on the modified couple stress theory,the equations of motion for the in-plane and out-of-plane motions are obtained employing the von Kármán plate theory as well as Kirchhoff's hypotheses by means of the Lagrange equations. The equations of motions are solved using the pseudo-arclength continuation technique and direct timeintegration method. The system parameters are tuned to the values associated with modal interactions, and then nonlinear resonant responses and energy transfer are analysed.Nonlinear motion characteristics are shown in the form of frequency-response and force-response curves, time histories, phase-plane portraits, and fast Fourier transforms.     

Further application of hybrid solution to another form of Boussinesq equations and comparisons

Recently, a new hybrid scheme is introduced for the solution of the Boussinesq equations. In this study, the hybrid scheme is used to solve another form of the Boussinesq equations. The hybrid solution is composed of finite‐volume and finite difference method. The finite‐volume method is applied to conservative part of the governing equations, whereas the higher order Boussinesq terms are discretized using the finite‐difference scheme. Fourth‐order accuracy is provided in both time and space. The solution is then applied to several test cases, which are taken from the previous studies. The results of this study are compared with experimental and theoretical results as well as those of the previous ones. The comparisons indicate that the Boussinesq equations solved here and in the previous study produce quite similar results. Copyright © 2006 John Wiley & Sons, Ltd.