Large deformation and electrochemistry of polyelectrolyte gels

Immersed in an ionic solution, a network of polyelectrolytes imbibes the solution and swells, resulting in a polyelectrolyte gel. The swelling is reversible, and the amount of swelling is regulated by ionic concentrations, mechanical forces, and electric potentials. This paper develops a field theory to couple large deformation and electrochemistry. A specific material model is described, including the effects of stretching the network, mixing the polymers with the solvent and ions, and polarizing the gel. We show that the notion of osmotic pressure in a gel has no experimental significance in general, but acquires a physical interpretation within the specific material model. The theory is used to analyze several phenomena: a gel swells freely in an ionic solution, a gel swells under a constraint of a substrate, electric double layer at the interface between the gel and the external solution, and swelling of a gel of a small size.     

Magnetohydrodynamic control of supersonic flow past bodies: Possibilities and undesirable effects

Basic problems of super-and hypersonic magnetohydrodynamics (MHD) associated with the determination of the integral characteristics of bodies and vehicles inside which there are systems generating a uniform magnetic field are considered. Three classes of flows, namely, flow in a hypersonic multimode fixed-geometry air-intake; internal and external flow in a model of a hypersonic vehicle containing an air-intake with an MHD generator, a combustion chamber, and a supersonic nozzle; and hypersonic flow past a blunt cone are studied using numerical simulation and theoretical analysis (on the basis of the complete averaged system of Navier-Stokes equations and the electrodynamic equations). Attention is concentrated on the presence of an additionalmagnetic force acting on the system generating themagnetic field and, consequently, on the body and initiating an additional drag (in the case of a vehicle-reducing its thrust). Attractive possibilities for MHD flow control, namely, an increase in the degree of flow compression in the air-intake, a reduction in the ignition length in the combustion chamber, and a decrease in the heat flux to the nose of the body, are noted, as well as negative effects associated with the action of the magnetic force on the bodies considered.     

A geometric study of shocks in equations that change type

In this paper we validate the generalized geometric entropy criterion for admissibility of shocks in systems which change type. This condition states that a shock between a state in a hyperbolic region and one in a nonhyperbolic region is admissible if the Lax geometric entropy criterion, based on the number of characteristics entering the shock, holds, where now the real part of a complex characteristic replaces the characteristic speed itself. We test this criterion by a nonlinear inviscid perturbation. We prove that the perturbed Cauchy problem in the elliptic region has a solution for a uniform time if the data lie in a suitable class of analytic functions and show that under small perturbations of the data a perturbed shock and a perturbed solution in the hyperbolic region exist, also for a uniform time.     

Creeping motion of a slip spherical particle in a circular cylindrical pore

A combined analytical–numerical study for the creeping flow caused by a spherical fluid or solid particle with a slip-flow surface translating in a viscous fluid along the centerline of a circular cylindrical pore is presented. To solve the axisymmetric Stokes equations for the fluid velocity field, a general solution is constructed from the superposition of the fundamental solutions in both cylindrical and spherical coordinate systems. The boundary conditions are enforced first at the pore wall by the Fourier transforms and then on the particle surface by a collocation technique. Numerical results for the hydrodynamic drag force acting on the particle are obtained with good convergence for various values of the relative viscosity or slip coefficient of the particle, the slip parameter of the pore wall, and the ratio of radii of the particle and pore. For the motion of a fluid sphere along the axis of a cylindrical pore, our drag results are in good agreement with the available solutions in the literature. As expected, the boundary-corrected drag force for all cases is a monotonic increasing function of the ratio of particle-to-pore radii, and approaches infinity in the limit. Except for the case that the cylindrical pore is hardly slip and the value of the ratio of particle-to-pore radii is close to unity, the drag force exerted on the particle increases monotonically with an increase in its relative viscosity or with a decrease in its slip coefficient for a constant ratio of radii. In a comparison for the pore shape effect on the axial translation of a slip sphere, it is found that the particle in a circular cylindrical pore in general acquires a lower hydrodynamic drag than in a spherical cavity, but this trend can be reversed for the case of highly slippery particles and pore walls.     

Spherical inflation of a class of compressible elastic bodies

Using a special model that belongs to a new class of elastic bodies wherein the Cauchy-Green stretch is given in terms of the Cauchy stress and its invariants, within the context of the spherical inflation of a spherical annulus, we show that interesting phenomena like the development of “stress boundary layers” manifest themselves. We consider two cases of boundary value problems, one in which there is a cavity in a sphere and the other in which there is a rigid spherical inclusion in a sphere. We show that in the case of a rigid inclusion, it is possible for a pronounced “stress boundary” layer to develop, in that the values of the stresses within this boundary layer that is adjacent to a spherical inclusion are much larger than external to it. We also show that in the case of both the cavity and a rigid inclusion, the stress concentration is an order of magnitude higher than the increase in the deformation gradient, that is, the stress and the stretch do not scale in a similar manner. While the stress adjacent to a rigid inclusion can be 2500 times the applied radial stress, the maximum stretch, which occurs at the rigid inclusion is about 10. While the variation in the stresses are linear in thin walled annular regions, we find that in thick walled annular regions, the variation of the stresses is non-linear.     

Unidimensional SPH simulations of reactive shock tubes in an astrophysical perspective

Smoothed Particle Hydrodynamics (SPH) is a Lagrangian method widely used for the modelling of a large variety of astrophysical fluid flows in more than one dimension. Simulations of thermonuclear explosions in stars require, besides the hydrodynamic equations, a realistic equation of state, an energy source term, and a set of nuclear kinetic equations to follow the composition changes of the gas during the explosion. The implementation of a realistic stellar equation of state, and the coupling of hydrodynamics and nuclear burning are investigated in the framework of the simple shock tube geometry. We present and discuss the results of a series of SPH simulations of a detonation in the presence of (1) a single exothermic nuclear reaction, and (2) a restricted network of nuclear reactions. Our results are compared to those of identical simulations performed by other authors using a different hydrodynamic method.     

Mechanics and chemical thermodynamics of phase transition in temperature-sensitive hydrogels

This paper uses the thermodynamic data of aqueous solutions of uncrosslinked poly(N-isopropylacrylamide) (PNIPAM) to study the phase transition of PNIPAM hydrogels. At a low temperature, uncrosslinked PNIPAM can be dissolved in water and form a homogenous liquid solution. When the temperature is increased, the solution separates into two liquid phases with different concentrations of the polymer. Covalently crosslinked PNIPAM, however, does not dissolve in water, but can imbibe water and form a hydrogel. When the temperature is changed, the hydrogel undergoes a phase transition: the amount of water in the hydrogel in equilibrium changes with temperature discontinuously. While the aqueous solution is a liquid and cannot sustain any nonhydrostatic stress in equilibrium, the hydrogel is a solid and can sustain nonhydrostatic stress in equilibrium. The nonhydrostatic stress can markedly affect various aspects of the phase transition in the hydrogel. We adopt the Flory-Rehner model, and show that the interaction parameter as a function of temperature and concentration obtained from the PNIPAM-water solution can be used to analyze diverse phenomena associated with the phase transition of the PNIPAM hydrogel. We analyze free swelling, uniaxially and biaxially constrained swelling of a hydrogel, swelling of a core-shell structure, and coexistent phases in a rod. The analysis is related to available experimental observations. Also outlined is a general theory of coexistent phases undergoing inhomogeneous deformation.     

STUDY ON THE FUEL AIR MIXING INDUCED BY A SHOCK WAVE PROPAGATING INTO A H_2-AIR INTERFACE

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｆｕｅｌａｉｒｍｉｘｉｎｇｉｎａｓｈｅａｒｌａｙｅｒｆｌｏｗｉｓａｎｉｍｐｏｒｔａｎｔｐｒｏｂｌｅｍｉｎｓｔｕｄｙｉｎｇｃｏｍｐｒｅｓｓｉｂｌｅｆｌｏｗａｎｄｓｕｐｅｒｓｏｎｉｃｃｏｍｂｕｓｔｉｏｎ ,ｓｕｃｈａｓｆｕｅｌｄｉｆｆｕｓｉｏｎａｎｄｍｉｘｉｎｇｉｎａＳｃｒａｍｊｅｔ[1].Ｒｏｓｈｋｏ[2 ]ｓｔｕｄｉｅｄｅｘｐｅｒｉｍｅｎｔａｌｌｙｔｈｅｐｈｅｎｏｍｅｎｏｎｏｆｆｕｅｌａｉｒｍｉｘｉｎｇｉｎａｓｕｂｓｏｎｉｃｓｈｅａｒｆｌｏｗａｎｄｆｏｕｎｄｌａｒｇｅ ,ｃｏｈ…     

Thin film cracking modulated by underlayer creep

In devices that integrate dissimilar materials in small dimensions, crack extension in one material often accompanies inelastic deformation in another. In this paper we analyze a channel crack advancing in an elastic film, while an underlayer creeps. The film is subject to a tensile stress. As the underlayer creeps, the stress field in the film relaxes in the crack wake, and intensifies around the crack tip. In a blanket film, the crack can attain a steady velocity, set by two rate processes: subcritical decohesion at the crack tip, and creep in the underlayer. In a thin-film microbridge over a viscous stripe, the crack cannot grow when the bridge is short, and can grow at a steady velocity when the bridge is long. We use a two-dimensional shear lag model to approximate the three-dimensional fracture process, and an extended finite element method to simulate the moving crack with an invariant, relatively coarse mesh. On the basis of the theoretical findings, we propose new experiments to measure fracture toughness and creep laws in small structures. As a byproduct, an analytical formula is found for the growth rate per temperature cycle of a channel crack in a brittle film, induced by ratcheting plastic deformation in a metal underlayer.     

A Cosserat shell model for interphases in elastic media

This study is concerned with the modeling of interphases in elastic media in general, and in composite materials in particular. The aim is to replace a boundary value problem consisting of a three-phase configuration, say that of fiber-interphase-matrix, by a simpler problem which involves the fiber and matrix only, plus certain matching conditions which simulate the interphase. The simplest of such known representations replaces a thin interphase by a “perfect contact interface” (a single surface) across which the displacements and tractions are assumed to be continuous. Another classical model replaces a thin and soft interphase by a “spring-type interface”, across which the tractions are continuous, but the displacement field undergoes a discontinuity. In the present paper, a Cosserat shell model of the interphase is derived which successfully models the original interphase in a unified manner, for the full range of its material parameters relative to those of the neighboring media. The model is derived in the setting of three-dimensional linear elasticity with small deformations and displacements. Comparisons with an existing exact solution of a coated fiber in an infinite matrix show that it performs extremely well even for moderately thick interphases.     

Free vibration of two coupled nonlinear oscillators

An analytical investigation is carried out on the free vibration of a two degree of freedom weakly nonlinear oscillator. Namely, the method of multiple time scales is first applied in deriving modulation equations for a van der Pol oscillator coupled with a Duffing oscillator. For the case of non-resonant oscillations, these equations are in standard normal form of a codimension two (Hopf-Hopf) bifurcation, which permits a complete analysis to be performed. Three different types of asymptotic states-corresponding to trivial, periodic and quasiperiodic motions of the original system-are obtained and their stability is analyzed. Transitions between these different solutions are also identified and analyzed in terms of two appropriate parameters. Then, effects of a coupling, a detuning, a nonlinear stiffness and a damping parameter are investigated numerically in a systematic manner. The results are interpreted in terms of classical engineering terminology and are related to some relatively new findings in the area of nonlinear dynamical systems.     

A finite strain model of stress, diffusion, plastic flow, and electrochemical reactions in a lithium-ion half-cell

We formulate the continuum field equations and constitutive equations that govern deformation, stress, and electric current flow in a Li-ion half-cell. The model considers mass transport through the system, deformation and stress in the anode and cathode, electrostatic fields, as well as the electrochemical reactions at the electrode/electrolyte interfaces. It extends existing analyses by accounting for the effects of finite strains and plastic flow in the electrodes, and by exploring in detail the role of stress in the electrochemical reactions at the electrode-electrolyte interfaces. In particular, we find that that stress directly influences the rest potential at the interface, so that a term involving stress must be added to the Nernst equation if the stress in the solid is significant. The model is used to predict the variation of stress and electric potential in a model 1-D half-cell, consisting of a thin film of Si on a rigid substrate, a fluid electrolyte layer, and a solid Li cathode. The predicted cycles of stress and potential are shown to be in good agreement with experimental observations.     

Nonlinear analyses of wrinkles in a film bonded to a compliant substrate

Subject to a compressive membrane force, a film bonded to a compliant substrate often forms a pattern of wrinkles. This paper studies such wrinkles in a layered structure used in several recent experiments. The structure comprises a stiff film bonded to a compliant substrate, which in turn is bonded to a rigid support. Two types of analyses are performed. First, for sinusoidal wrinkles, by minimizing energy, we obtain the wavelength and the amplitude of the wrinkles for substrates of various moduli and thicknesses. Second, we develop a method to simultaneously evolve the two-dimensional pattern in the film and the three-dimensional elastic field in the substrate. The simulations show that the wrinkles can evolve into stripes, labyrinths, or herringbones, depending on the anisotropy of the membrane forces. Statistical averages of the amplitude and wavelength of wrinkles of various patterns correlate well with the analytical solution of the sinusoidal wrinkles.     

Finite element computations of viscoelastic two-phase flows using local projection stabilization

A three-field local projection stabilized (LPS) finite element method is developed for computations of a three-dimensional axisymmetric buoyancy driven liquid drop rising in a liquid column where one of the liquid is viscoelastic. The two-phase flow is described by the time-dependent incompressible Navier-Stokes equations, whereas the viscoelasticity is modeled by the Giesekus constitutive equation in a time-dependent domain. The arbitrary Lagrangian-Eulerian (ALE) formulation with finite elements is used to solve the governing equations in the time-dependent domain. Interface-resolved moving meshes in ALE allows to incorporate the interfacial tension force and jumps in the material parameters accurately. A one-level LPS based on an enriched approximation space and a discontinuous projection space is used to stabilize the numerical scheme. A comprehensive numerical investigation is performed for a Newtonian drop rising in a viscoelastic fluid column and a viscoelastic drop rising in a Newtonian fluid column. The influence of the viscosity ratio, Newtonian solvent ratio, Giesekus mobility factor, and the Eötvös number on the drop dynamics are analyzed. The numerical study shows that beyond a critical Capillary number, a Newtonian drop rising in a viscoelastic fluid column experiences an extended trailing edge with a cusp-like shape and also exhibits a negative wake phenomena. However, a viscoelastic drop rising in a Newtonian fluid column develops an indentation around the rear stagnation point with a dimpled shape.     

On the Need for Control of Nonlinear Oscillations in Machine Systems

Oscillations in machines are invariably nonlinear. This is either because of inertial coupling effects between different motions of the moving components, material and constitutive phenomena giving rise to stiffness modifications, nonlinear dissipation mechanisms, large deflections, or, as is most likely, some sort of combination of all of these. The net effect of nonlinear vibrations is that at best the machine may well behave a little differently from the way the designer intended, or at worst, in a manner which renders it completely unsuitable for the job. The extent of such problems depends on the nature and the scale of the nonlinearities that are present but it is safe to say that nonlinear oscillations can rarely be completely overlooked in precision machinery analysis and design. The unifying theme in this paper is pendulum motion, firstly in the case of a mobile gantry crane for container stacking where we wish to minimise such motion and converge on a target, and then secondly in the case of a vibration absorber in which we choose to initiate pendulum motion within a special absorber, for the purposes of vibration minimisation. The third example involves the potential for pendulum motion at a very much larger scale and summarises the main control problem that is likely to be encountered in a fully deployed momentum exchange propulsion tether operating in space. The paper discusses the general mathematical issues that pertain to pendulum motion in each of the three cases. This motion is investigated initially in the context of the mobile gantry crane, in the form of a basic three dimensional dynamical model. A feedback linearised controller is shown to offer some advantages for the control of such a system and then a simulation based on data from a practical implementation of this within a real-time control system on a 1/10 laboratory scale model is discussed. It is recalled that the real-time effectiveness of the controller can be compromised by relatively slow sensing and data logging hardware but that despite this some useful performance gains can still be obtainable using this sort of control strategy. The second example comprises an autoparametric vibration absorber and here it is shown how even a simple hunting controller can exploit the mode-locking and wide detuning region effects inherent in autoparametric systems. Further experimental results are discussed in the case of a hunting controller for detuning a vertically oriented parametrically excited pendulum in order to exploit and enhance the powerful and persistent absorption available during autoparametric interaction. The paper concludes with a summary review of the third problem in which the theoretical attitude dynamics of a motorised momentum exchange space propulsion tether are summarised and it is shown that they need to be controlled for reliable and optimal payload velocity boost from both circular parking orbits and elliptical transfer orbits about the Earth.     

On well-posedness of the dynamic problem for an anisotropic fluid-saturated solid

It is known that a high degree of anisotropy in the constitutive behaviour of a solid may result in the loss of hyperbolicity of the dynamic equations in the form of either complex-conjugate or purely imaginary characteristic wave speeds (flutter ill-posedness and shear band formation, respectively). In the present paper we investigate the characteristic wave speeds in the dynamic problem for a transversely isotropic fluid-saturated porous solid. Three cases are considered: a dry solid and a saturated solid under locally undrained and drained conditions. It is shown that, for given constitutive parameters of the solid skeleton, the dynamic problem for a drained solid may become ill-posed due to the flutter-type loss of hyperbolicity, while the dynamic equations for a dry and an undrained solids remain hyperbolic. For a given solid skeleton, the characteristic wave speeds are strongly influenced by the pore fluid compressibility which, in turn, is extremely sensitive to the presence of a small amount of free gas.     

Mechanism of unsteady aerodynamic heating with sudden change in surface temperature

The characteristics and mechanism of unsteady aerodynamic heating of a transient hypersonic boundary layer caused by a sudden change in surface temperature are studied. The complete time history of wall heat flux is presented with both analytical and numerical approaches. With the analytical method, the unsteady compressible boundary layer equation is solved. In the neighborhood of the initial and final steady states, the transient responses can be expressed with a steady-state solution plus a perturbation series. By combining these two solutions, a complete solution in the entire time domain is achieved. In the region in which the analytical approach is applicable, numerical results are in good agreement with the analytical results, showing reliability of the methods. The result shows two distinct features of the unsteady response. In a short period just after a sudden increase in the wall temperature, the direction of the wall heat flux is reverted, and a new inflexion near the wall occurs in the profile of the thermal boundary layer. This is a typical unsteady characteristic. However, these unsteady responses only exist in a very short period in hypersonic flows, meaning that, in a long-term aerodynamic heating process considering only unsteady surface temperature, the unsteady characteristics of the flow can be ignored, and the traditional quasi-steady aerodynamic heating prediction methods are still valid.     

分层流体中运动源生成的内波研究进展

针对两类密度分布模型------连续分层流体和间断分层流体, 综述了在运动潜体生成的Kelvin型和非Kelvin型内尾迹研究方面的现状, 内容侧重于运动源生成内波的解析理论和分层拖曳水槽中内尾迹实验方面的研究成果. 介绍了在连续分层流体中运动源生成的Kelvin型非线性内波的一般方程和在间断分层流体中Kelvin型内波的势流分析的一般方法; 概述了运动源诱生的先锋内孤立子、代数孤立子和平孤立波3类特殊非线性内波的研究进展, 其中运动潜体生成的平孤立内波被作者实验证实是一类极限孤立波, 并首次建立了共轭流动模型予以描述; 综合分析了在密度线性分布流体中潜体运动生成内波的动力学过程多样性特征, 其中包括内尾迹近场和远场的时空结构、不稳定结构、涡旋与湍流耦合结构以及湍流与内波相互作用结构等.      

Micromechanics of cleavage fracture initiation in ferritic steels by carbide cracking

Cleavage fracture in ferritic steels is often initiated in brittle carbides randomly distributed in the material. The carbides break as a result of a fibre loading mechanism in which the stress levels in the carbides are raised, as the surrounding ferrite undergoes plastic deformation. The conditions in the vicinity of the nucleated micro-crack will then determine whether the crack will penetrate or be arrested by the ferrite. The ferrite is able to arrest nucleated cracks through the presence of mobile dislocations, which blunt and shield the microcrack and thus lowers the stresses at the crack tip. Hence, the macroscopic toughness of the material directly depends on the ability of the ferrite to arrest nucleated micro-cracks and in turn on the plastic rate sensitivity of the ferrite. The initiation of cleavage fracture is here modelled explicitly in the form of a micro-crack, which nucleates in a brittle carbide and propagates into the surrounding ferrite. The carbide is modelled as an elastic cylinder or in a few cases an elastic sphere and the ferrite as an elastic viscoplastic material. The crack growth is modelled using a cohesive surface, where the tractions are governed by a modified exponential cohesive law. It is shown that the critical stress, required to propagate a microcrack from a broken carbide, increases with decreasing plastic rate sensitivity of the ferrite. The results also show that a low stress triaxiality and a high aspect ratio of the carbide promote the initiation of cleavage fracture from a broken carbide.