A differential approach to suspensions with power-law matrices

We apply the differential method of Roscoe (1952) to the problem of finding the viscosity of a suspension of non-colloidal spheres in a power-law matrix. The results are compared with other theories and published experiments, and reasonable agreement is found up to moderate concentrations (? ～ 0.5) when viscoelasticity and other effects are not important. The Roscoe paper depends on using a “crowding” function in the analysis; here two modified crowding functions are discussed, with a view to explaining the success of the Maron–Pierce formula, which does not reduce to the Einstein form for low concentrations.     

Natural frequencies of a hemispherical shell

An experimental study of the natural frequencies of a hemispherical thin shell with a clamped edge is described. For the particular shell model studied, the four lowest symmetric modes were identified conclusively by establishing the node-line patterns. The associated natural frequencies were found to be in excellent agreement with theoretical predictions. A great number of higher modes were also excited but could not be conclusively identified due to difficulties in establishing the node-line patterns. The symmetric frequency spectrum up to and including the 23rd mode was, however, established by means of a qualitative argument and showed good correlation with the theoretical spectrum. A number of natural frequencies were found to be present in a particular frequency band for the test shell resulting in a severe crowding in this region. This was found to be in agreement with the more general results for hemispherical shells which indicate that the crowding takes place around a frequency parameterp=1.0.     

A simple semiempirical model for the effective viscosity of multicomponent suspensions

We propose a methodology to approximate the viscosity of multicomponent suspensions. The procedure consists of successive applications of expressions for the viscosity of binary mixtures, originally written as the product of monomodal stiffening functions. First, the viscosity of a binary mixture made of the two smallest components is calculated. This allows to extract a volume fraction that will be used, together with the volume fraction of the third component, to feed the next iteration of the procedure to calculate the viscosity of a trimodal mixture and so on. The application of this approach to arbitrary mixtures requires the detailed knowledge of the geometry of the system in the form of size ratios and compositions. When this information is unknown, an approximation of the model can still be used as a fitting tool. With that purpose, the final expression for the viscosity is written in terms of an effective volume fraction that is further approximated by the use of a (1,2) Padé approximant. This approximation allows to incorporate the crowding effects due to different species in a volume fraction-dependent crowding factor that can be used as a fitting parameter to match experimental or simulation data. We have applied the model to mixtures of particles with different sizes and tested its accuracy comparing with experimental results obtaining very good agreement.     

Key factors in chaperonin-assisted protein folding

Incorrect folding of proteins in the macromolecular crowding environment in living cells would cause cellular disasters. All cells have developed their structurally and functionally distinct classes of molecular chaperones to help nonnative proteins fold to their native structures, one of which is the most studied GroEL/ES complex. In the present article, large-scale all-atom explicit solvent molecular dynamics (MD) simulations have been carried out on rhodanese folding in a series of chaperonin mutants for 200 ns to understand the mechanism therein. In accordance with experimental results, two factors have been identified to play a significant role, the geometrical confinement effect of the folding cavity and the charge effect of the inner surface of the cavity. Our analysis of the properties during simulation suggests that the GroEL/ES complex directly exerts force on the contacting residues of the substrate, thus assists substrate folding.     

Improved non-dominated sorting genetic algorithm （NSGA）-II in multi-objective optimization studies of wind turbine blades

The non-dominated sorting genetic algorithm (NSGA) is improved with the controlled elitism and dynamic crowding distance. A novel multi-objective optimization algorithm is obtained for wind turbine blades. As an example, a 5 MW wind turbine blade design is presented by taking the maximum power coefficient and the minimum blade mass as the optimization objectives. The optimal results show that this algorithm has good performance in handling the multi-objective optimization of wind turbines, and it gives a Pareto-optimal solution set rather than the optimum solutions to the conventional multiobjective optimization problems. The wind turbine blade optimization method presented in this paper provides a new and general algorithm for the multi-objective optimization of wind turbines.     

Pedestrian flows in bounded domains with obstacles

In this paper, we systematically apply the mathematical structures by time-evolving measures developed in a previous work to the macroscopic modeling of pedestrian flows. We propose a discrete-time Eulerian model, in which the space occupancy by pedestrians is described via a sequence of Radon-positive measures generated by a push-forward recursive relation. We assume that two fundamental aspects of pedestrian behavior rule the dynamics of the system: on the one hand, the will to reach specific targets, which determines the main direction of motion of the walkers; on the other hand, the tendency to avoid crowding, which introduces interactions among the individuals. The resulting model is able to reproduce several experimental evidences of pedestrian flows pointed out in the specialized literature, being at the same time much easier to handle, from both the analytical and the numerical point of view, than other models relying on nonlinear hyperbolic conservation laws. This makes it suitable to address two-dimensional applications of practical interest, chiefly the motion of pedestrians in complex domains scattered with obstacles.      

Reconstruction of road defects and road roughness classification using vehicle responses with artificial neural networks simulation

The road damage assessment methodology in this paper utilizes an artificial neural network that reconstructs road surface profiles from measured vehicle accelerations. The paper numerically demonstrates the capabilities of such a methodology in the presence of noise, changing vehicle mass, changing vehicle speeds and road defects. In order to avoid crowding out understanding of the methodology, a simple linear pitch-plane model is employed. Initially, road profiles from known roughness classes were applied to a physical model to calculate vehicle responses. The calculated responses and road profiles were used to train an artificial neural network. In this way, the network renders corresponding road profiles on the availability of fresh data on model responses. The results show that the road profiles and associated defects can be reconstructed to within a 20% error at a minimum correlation value of 94%.     

Elasticity and dynamics of particle gels in non-Newtonian melts

We investigate the relation between the structure and the viscoelastic behavior of a model polymer nanocomposite system based on a mixture of titanium dioxide (TiO₂) nanoparticles and polypropylene. Above a critical volume fraction, Φ _c, the elasticity of the hybrids dramatically increases, and the frequency dependence of the elastic and viscous moduli reflects the superposition of the independent responses of the suspending polymer melt and of an elastic particle network. In addition, the elasticity of the hybrids shows critical behavior around Φ _c. We interpret these observations by hypothesizing the formation of a transient network, which forms due to crowding of particle clusters. Consistent with this interpretation, we find a long-time, Φ-dependent, structural relaxation, which emphasizes the transient character of the structure formed by the particle clusters. For times below this characteristic relaxation time, the elasticity of the network is Φ-independent and reminiscent of glassy behavior, with the elastic modulus, G^′, scaling with frequency, ω, as G^′∼ω ^0.3. We expect that our analysis will be useful for understanding the behavior of other complex fluids where the elasticity of the components could be superimposed.     

Darcian Dynamics: A New Approach to the Mobilization of a Trapped Phase in Porous Media

We develop a new approach, which we term Darcian Dynamics, to simulate two-phase (liquid-gas) flow in porous media, when the gas phase is disconnected in the form of ganglia. The method is based on the assumption of homogeneous fluid flow for the liquid, although it does allow for heterogeneous capillary thresholds due to the pore microstructure. Using techniques from potential theory, the hydrodynamic interaction between liquid and gas is expressed through an integral representation over the ganglia interfaces. We use a numerical method to solve the resulting integral equation, and explore conditions for the onset of ganglia mobilization as well as for subsequent events, such as break-up, coalescence and stranding. The interaction between the ganglia and the flowing phase is influenced by the capillary and gravity (Bond) numbers, and by geometric factors, such as size, orientation, and ganglia density. The latter effect depends on the hydrodynamic interaction in addition to the intuitively expected crowding effect.     

A generalized equation for rheology of emulsions and suspensions of deformable particles subjected to simple shear at low Reynolds number

We present analyses to provide a generalized rheological equation for suspensions and emulsions of non-Brownian particles. These multiparticle systems are subjected to a steady straining flow at low Reynolds number. We first consider the effect of a single deformable fluid particle on the ambient velocity and stress fields to constrain the rheological behavior of dilute mixtures. In the homogenization process, we introduce a first volume correction by considering a finite domain for the incompressible matrix. We then extend the solution for the rheology of concentrated system using an incremental differential method operating in a fixed and finite volume, where we account for the effective volume of particles through a crowding factor. This approach provides a self-consistent method to approximate hydrodynamic interactions between bubbles, droplets, or solid particles in concentrated systems. The resultant non-linear model predicts the relative viscosity over particle volume fractions ranging from dilute to the the random close packing in the limit of small deformation (capillary or Weissenberg numbers) for any viscosity ratio between the dispersed and continuous phases. The predictions from our model are tested against published datasets and other constitutive equations over different ranges of viscosity ratio, volume fraction, and shear rate. These comparisons show that our model, is in excellent agreement with published datasets. Moreover, comparisons with experimental data show that the model performs very well when extrapolated to high capillary numbers (C a?1). We also predict the existence of two dimensionless numbers; a critical viscosity ratio and critical capillary numbers that characterize transitions in the macroscopic rheological behavior of emulsions. Finally, we present a regime diagram in terms of the viscosity ratio and capillary number that constrains conditions where emulsions behave like Newtonian or Non-Newtonian fluids.     

基于多目标优化NSGA2改进算法的结构动力学模型确认

传统结构动力学模型确认方法通常采用单目标优化,存在精度不足和稳定性差等缺点,难以满足实际工程需求。基于此,提出一种采用神经网络作为代理模型,建立以马氏距离和鲁棒性为不确定性量化指标的多目标优化模型,并将NSGA2多目标进化算法用于求解。针对NSGA2存在无法有效识别伪非支配解、计算效率低和解集质量较差等设计缺陷,提出一种基于支配强度的NSGA2改进算法INSGA2-DS。INSGA2-DS将支配强度引入非支配排序,采用新型拥挤距离公式和自适应精英保留策略,以提高收敛效率和解集质量。GARTEUR飞机算例的仿真结果表明,INSGA2-DS求解复杂工程问题时具有更好的收敛性和分布性,而考虑鲁棒性的结构动力学模型确认方法可以获得同时满足多种目标要求的Pareto解集,提高了模型确认的精度和稳定性。     

Frontogenesis in turbulent flow through porous media using non-linear Forchheimer equation

We present here both one- and two-dimensional models for turbulent flow through heterogeneous unbounded fluid saturated porous media using non-linear Forchheimer extended Darcy (DF) equation in the presence of gravitational field. The fluid is initially at rest and sets in motion due to a uniform horizontal density gradient. It is shown that a purely horizontal motion develops satisfying non-linear DF equation. Analytical solutions of this non-linear Initial Value Problem are obtained and limiting solutions valid for the Darcy regime in the case of laminar flow are derived. A measure of the stability of the flow is discussed briefly using Richardson number. The comparison between the nature of the solutions satisfying the non-linear and linear initial value problems are made. We found that even in the case of turbulent flow the vertical density gradient varies continuously both with space z and time t but the horizontal density gradient remains unchanged. The existence and uniqueness theorem of the Initial Value Problem is proved. The stability of these solutions are discussed and it is shown that the solutions are qualitatively and quantitatively different for and in the upper and lower half of the region. In particular, we have shown that the solution which is stable for infinitesimal perturbations is also stable for arbitrary perturbations both in time and space.In the case of two-dimensional motion, a piecewise initial density gradient with continuous distribution of density, stream function formulation is used and the solutions are obtained using time-series analysis. In this case solution shows crowding of the density profiles in the lower-half of the channel reflecting an increase in density gradient and incipient of frontogenesis there, because of the increase in circulation of the flow due to piecewise initial density gradient.     

多物理场下FCBGA焊点电迁移失效预测的数值模拟研究

随着微电子封装技术的快速发展, 焊点的电迁移失效问题日益受到关注. 基于有限元法并结合子模型技术对倒装芯片球栅阵列封装(flip chip ball grid array, FCBGA)进行电-热-结构多物理场耦合分析, 详细介绍了封装模型的简化处理方法, 重点分析了易失效关键焊点的电流密度分布、温度分布和应力分布, 发现电子流入口处易产生电流拥挤效应, 而整个焊点的温度梯度较小. 基于综合考虑“电子风力”、温度梯度、应力梯度和原子密度梯度四种电迁移驱动机制的原子密度积分法, 并结合空洞形成/扩散准则及失效判据, 分析FCBGA焊点在不同网格密度下的电迁移空洞演化过程, 发现原子密度积分算法稳定, 不依赖网格密度. 采用原子密度积分法模拟真实 工况下FCBGA关键焊点电迁移空洞形成位置和失效寿命, 重点研究了焊点材料和铜金属层结构对电迁移失效的影响. 结果表明, 电迁移失效寿命随激活能的增加呈指数级增加, 因此Sn3.5Ag焊点的电迁移失效寿命约为63Sn37Pb的2.5倍, 有效电荷数对电迁移寿命也有一定的影响;铜金属层结构的调整会改变电流的流向和焊点的应力分布, 进而影响焊点的电迁移失效寿命.      

Viscosity of suspensions modeled with a shear-dependent maximum packing fraction

A large amount of data from the literature on viscosity of concentrated suspensions of rigid spherical particles are analyzed to support the new concept that the maximum packing fraction ( _M ) is shear-dependent. Incorporation of this behavior in a rheological model for viscosity () as a function of particle volume fraction () succeeds in describing virtually all non-Newtonian effects over the entire concentration range and also accounts for a yield stress. The most successful model is one proposed by Krieger and Dougherty for Newtonian viscosities, (, _M ), but with _M varying from a low-shear limit _M0 to a high-shear limit _M. Microstructural interpretations of this behavior are advanced, with arguments suggesting that similar rheological models should apply to suspensions of nonspherical and irregular particles.Symbols a particle size scale (for spheres, the diameter) - A lumped kinetic parameter in eqs. (23) and (24) - BS butadiene-styrene - C coefficient in Arrhenius model, eq. (2) - D coefficient in Mooney model, eq. (3) - e _i parameter representing one of the three electroviscous effects (i = 1, 2, or 3) - f fraction of total particulates that exist in the dispersed phase, eq. (22) - h solution factor, in Arrhenius model, eq. (2) - k crowding factor, in Mooney model, eq. (3) - k _D ,k _F kinetic rate coefficient for producing particles of dispersed or flocculated type, respectively - K Einstein coefficient for particles of any shape, eq. (1); equal to [] - KD Krieger-Dougherty model, eq. (6) - m exponent to characterize shear-dependence in viscosity models of Cross, eq. (10), and eq. (23), and also in yield stress prediction eq. (24) - N number of monodisperse components in a blend of spheres with different diameters - PD polydispersity (in size) parameter - S generalized shape parameter - T temperature - V _c volume of chamber in figure 6, representing the entire volume of the sample - V _P total volume of particles in the sample - V _D ,V _F sample volumes in which dispersed particles or flocculated particles, respectively, prevail; volumes of the dispersed phase or flocculated phase, containing both particles and carrier fluid - V _PD ,V _PF particle volume within the phase volumeV _D orV _F , respectively Greek coefficient in definition of _c in eq. (8); of order unity - coefficient regulating -sensitivity in eq. (10) - shear rate,dv ₁/dx ₂ in simple shear - shear viscosity of the suspension - ₀, low-shear and high-shear limiting values of - _s viscosity of the suspending fluid - [] intrinsic viscosity, - _r reduced viscosity,/ _s - Boltzmann's constant; in _c - shear stress - _c parameter characterizing sensitivity of viscosity to stress, in eq. (8) - _B dynamic yield stress in the floc model - _y yield stress - volume fraction occupied by solids in a suspension - _M maximum value of attainable by a given collection of particles under given conditions of flow - _M0, _M limiting values of _M at the low- and high- conditions, respectively     

Investigation of the stability problems of elastic bodies using the method of mathematical theory of elasticity

In this paper, using the equilibrium equations and boundary conditions of elastic stability problem of and the method of mathematical theory of elasticity, we solve some elastic stability problems, which were studied by Ишлынский^[2] and Войцеховская^[3,4],and obtained more reasonable results than theirs.     

On Asymptotics of Solutions of Nonlinear Differential Equations of the Second Order

We study the problem of asymptotics of unbounded solutions of differential equations of the form y″ = α₀ p(t)ϕ(y), where α₀ ∈ {−1, 1}, p: [a, ω[→]0, +∞[, −∞ < a < ω ≤ +∞, is a continuous function, and ϕ: [y ₀, +∞[→]0, +∞[ is a twice continuously differentiable function close to a power function in a certain sense.__________Translated from Neliniini Kolyvannya, Vol. 8, No. 1, pp. 18–28, January–March, 2005.     

Experimental investigations of the stability limit of the helical flow of pseudoplastic liquids

The stability of the laminar helical flow of pseudoplastic liquids has been investigated with an indirect method consisting in the measurement of the rate of mass transfer at the surface of the inner rotating cylinder. The experiments have been carried out for different values of the geometric parameter = R ₁/R ₂ (the radius ratio) in the range of small values of the Reynolds number,Re < 200. Water solutions of CMC and MC have been used as pseudoplastic liquids obeying the power law model. The results have been correlated with the Taylor and Reynolds numbers defined with the aid of the mean viscosity value. The stability limit of the Couette flow is described by a functional dependence of the modified critical Taylor number (including geometric factor) on the flow indexn. This dependence, general for pseudoplastic liquids obeying the power law model, is close to the previous theoretical predictions and displays destabilizing influence of pseudoplasticity on the rotational motion. Beyond the initial range of the Reynolds numbers values (Re>20), the stability of the helical flow is not affected considerably by the pseudoplastic properties of liquids. In the range of the monotonic stabilization of the helical flow the stability limit is described by a general dependence of the modified Taylor number on the Reynolds number. The dependence is general for pseudoplastic as well as Newtonian liquids.Nomenclature C _i concentration of reaction ions, kmol/m³ - d = R ₂ –R ₁ gap width, m - F _M () Meksyn's geometric factor (Eq. (1)) - F ₀ Faraday constant, C/kmol - i _l density of limit current, A/m³ - k _c mass transfer coefficient, m/s - n flow index - R ₁,R ₂ inner, outer radius of the gap, m - Re = V _m ·2d·/µ _m Reynolds number - Ta _c = _c ·d^3/2·R ₁ ^1/2 ·/µ _m Taylor number - Z _i number of electrons involved in electrochemical reaction - = R ₁/R ₂ radius ratio - µ apparent viscosity (local), Ns/m² - µ _m mean apparent viscosity value (Eq. (3)), Ns/m² - µ _i apparent viscosity value at a surface of the inner cylinder, Ns/m² - density, kg/m³ - _c angular velocity of the inner cylinder (critical value), 1/s     

On the orders of the best approximations of integrals of functions by integrals of rank σ

We study the values e _σ(f) of the best approximation of integrals of functions from the spaces L _p (A, dμ) by integrals of rank σ. We determine the orders of the least upper bounds of these values as σ → ∞ in the case where the function ƒ is the product of two nonnegative functions one of which is fixed and the other varies on the unit ball U _p (A) of the space L _p (A, dμ). We consider applications of the obtained results to approximation problems in the spaces S _p ^ϕ. __________ Translated from Neliniini Kolyvannya, Vol. 10, No. 4, pp. 528–559, October–December, 2007.     

The effect of random step changes of system parameters on the stability of steady vibration of nonlinear systems

Summary The author proposes a method of analysing the effect of random step changes of system parameters on the stability of steady solution of a nonlinear system in case this steady solution is not the only one. The method is based on the assumption that the time between successive step changes is comparatively long against the period of steady vibration so that transients become stabilized after each parameter change. The parameter is free to vary on a finite interval of values.

---

Übersicht Es wird ein Verfahren vorgeschlagen, um die Auswirkung von zufälligen Sprungstörungen der Systemparameter auf die Stabilität von stationären Lösungen nichtlinearer Systeme zu untersuchen. Dabei wird vorausgesetzt, daß es mehrere stationäre Lösungen gibt. Es wird angenommen, daß die Zeit zwischen den nacheinander folgenden Sprungstörungen vergleichsweise lang gegenüber der Periode der stationären Schwingung ist, so daß die Einschwingvorgänge nach jeder Parameteränderung abgeklungen sind. Der Parameter kann sich dabei in einem endlichen Wertcbcreich ändern.

---

---

Herrn Professor Dr. K. Klotter zum 70. Geburtstag gewidmet.     

Local and global behavior of solutions of quasilinear equations of Emden-Fowler type

We consider the equation div for p≦N, 0<p−1<q. We study the isolated singularities and the behavior near infinity of nonradial positive solutions when q <N(p −1)/(N − p), and give a complete classification of local and global radial solutions of any sign, for any q.