混合层中柱状粒子取向分布函数的研究

数值模拟了平面流固混合层中柱状粒子在平面取向的分布函数，提出了描述流场中某点直粒子取向占优程度的参数Dψ的分析，说明在混合层中心区域粒子的取向排列是有序的，在混合层相邻两个大涡之间的区域，是粒子最可几角占优程度最高的地方，而且最可几取向就是沿流场的拉伸方向，说明流场的拉伸作用及方向对粒子的取向排列具有重要影响。与其它模型相比，粒子取向的分布函数模型能定量反映流场大涡对粒子取向控制的强弱程度。     

Numerical simulation of flow kinematics and fiber orientation for multi-disperse suspension

The development of flow kinematics and fiber orientation distribution from the parabolic velocity profile and isotropic orientation at the channel inlet was computed in multi-disperse suspension flow through a parallel plate channel and their predictions were compared with those of mono- and bi-disperse suspensions. A statistical scheme (orientations of a large number of fibers are evaluated from the solution of the Jeffery equation along the streamlines) was confirmed to be very useful and feasible method to analyze accurately the orientation distribution of fibers in multi-disperse fiber suspension flow as well as mono- and bi-dispersions, instead of direct solutions of the orientation distribution function of fibers or the evolution equation of the orientation tensor which involves a closure equation. It was found that the flow kinematics and the fiber orientation depend completely on both the fiber aspect-ratio and the fiber parameter for multi-disperse suspension when the fiber–fiber and fiber-wall interactions are neglected. Furthermore, the addition of large aspect-ratio fibers as well as an increase in the fiber parameter related to the large aspect-ratio fibers could suppress the complex velocity field and stress distributions which are observed in suspensions containing small aspect-ratio fibers. From a practical point of view, therefore, the mechanical and physical properties of fiber composites should be improved with an increase in the volume fraction of large aspect-ratio fibers.     

Anisotropic work hardening of steel sheets under plane stress

This work is a follow-up of the previous report by Kim and Yin [Kim, K.H., Yin, J.J., 1997. Evolution of anisotropy under plane stress. J. Mech. Phys. Solids 45, 841–851] regarding the anisotropic work hardening of cold rolled steel sheets. Tensile prestrain has been applied at angles to the rolling direction and then tensile uniaxial yield stress and R-value distributions are measured. As reported earlier, the orientations of local maxima and minima in the yield stress are altered when the prestrain axis is not in the rolling direction. This led Kim and Yin [Kim and Yin (1997)] to suggest that the orientations of orthotropy axes are altered by the tensile prestrain at angles to the rolling direction. However, R-value distribution is found to be hardly affected by the prestrain. The unchanging R-value distribution shows that the material remembers the rolling direction even after the prestrain. An attempt is made to approximate the observed yield and flow behavior based upon isotropic-kinematic hardening with the quadratic yield function (Hill, 1948). The degree of approximation raises the issues of yield point definition, flexibility of yield function, non-associated flow rules, distortional hardening and others.     

Fiber suspension flow in a tapered channel: The effect of flow/fiber coupling

A numerical model for predicting the flow and orientation state of semi-dilute, rigid fiber suspensions in a tapered channel is presented. The effect of the two-way flow/fiber coupling is investigated for low Reynolds number flow using the constitutive model of Shaqfeh and Fredrickson. An orientation distribution function is used to describe the local orientation state of the suspension and evolves according to a Fokker–Plank type equation. The planar orientation distribution function is determined along streamlines of the flow and is coupled with the fluid momentum equations through a fourth-order orientation tensor. The coupling term accounts for the two-way interaction and momentum exchange between the fluid and fiber phases. The fibers are free to interact through long range hydrodynamic fiber–fiber interactions which are modeled using a rotary diffusion coefficient, an approach outlined by Folgar and Tucker. Numerical predictions are made for two different orientation states at the inlet to the contraction, namely a fully random and a partially aligned fiber orientation state. Results from these numerical predictions show that the streamlines of the flow are altered and that velocity profiles change from Jeffery–Hamel, to something resembling a plug flow when the fiber phase is considered in the fluid momentum equations. This phenomenon was found when the suspension enters the channel in either a pre-aligned, or in a fully random orientation state. When the suspension enters the channel in an aligned orientation state, fiber orientation is shown to be only marginally changed when the two-way coupling is included. However, significant differences between coupled and uncoupled predictions of fiber orientation were found when the suspension enters the channel in a random orientation state. In this case, the suspension was shown to align much more quickly when the mutual coupling was accounted for and profiles of the orientation anisotropy were considerably different both qualitatively and quantitatively.     

纳米多晶NiTi形状记忆合金超弹性的晶粒取向依赖性相场研究

本文基于Ginzburg-Landau理论,建立了一个反映纳米多晶NiTi形状记忆合金取向依赖性的二维多晶相场模型,研究了晶粒取向对其超弹性性能的影响。结果表明,纳米多晶NiTi形状记忆合金的超弹性行为依赖于晶粒取向分布,即：多晶模型中在所研究的参数变化范围内,晶粒取向分布范围越广、晶粒间取向差越大（无明显织构）,超弹性性能越差;而晶粒取向分布范围越窄、晶粒间取向差越小（织构越强）,超弹性性能越好。该晶粒取向依赖性可解释为：由于晶粒取向的不同,马氏体相变过程中相邻晶粒间的变形不匹配程度不同,因此,多晶模型中在所研究的参数变化范围内,晶粒间取向差异越大,晶界处的变形失配越严重,由此而产生的局部内应力将阻碍其附近马氏体相变的扩展,进而导致纳米多晶NiTi形状记忆合金超弹性性能下降。     

Theory of non-isotropic spatial resolution in PIV

The spatial resolution of the PIV interrogation technique is discussed from an analytical standpoint and assessed with Monte Carlo numerical simulation of particle image motion. The PIV measurement error associated with lack of spatial resolution is modelled associating the cross-correlation operator to a moving average filter. The error associated with the "low-pass filtering" effect is investigated by adopting a second-order polynomial expression for the velocity spatial distribution. According to the present error analysis, the measurement error is proportional to the second-order spatial derivative of the velocity field and increases with the square of the window linear size. The strategy for the selection of the window size and properties (aspect ratio and orientation) so as to minimize the error is discussed. The principle is based on nonisotropic interrogation windows of elliptical shape, with a constant area and elongated in the direction of the largest curvature radius. The nonisotropic parameters are defined as eccentricity and orientation, which are based on the local eigenvalues/vectors of the Hessian tensor of the displacement spatial distribution. The technique is implemented in a recursive PIV interrogation method. The performance of nonisotropic interrogation technique is assessed by means of synthetic PIV images, which simulate three situations: first, a one-dimensional sinusoidal shear displacement, which allows comparison of the cross-correlation spatial response with the transfer function of linear filters. Second, the stream-wise exponential velocity decay is simulated, which simulates the particle tracers decelerating downstream of a shock wave and gives an example of a flow with main velocity differences aligned with the velocity direction. The results show that keeping the image density fixed, the error caused by insufficient spatial resolution can be reduced by a factor two when a preferential direction is found in the flow field. Finally, a Lamb–Oseen vortex flow is presented, which shows the complex pattern formed by the interrogation windows in a two-dimensional case. In this case, the improvement in interrogation performance is limited due to the isotropic nature of the velocity spatial fluctuation.     

Grid Orientation Revisited: Near-well, Early-time Effects and Solution Coupling Methods

The grid orientation effect is a phenomenon which leads to the computation of fundamentally different solutions on grids oriented diagonal and parallel to the principal flow direction. Grid orientation remains an important consideration for many practical simulation studies, and renewed interest in gas injection processes motivates the revisiting of this classical problem. In this article, we show that there are aspects of the grid orientation effect that can be traced back directly to the treatment of early-time, near-well flow and therefore have a major impact on adverse mobility ratio displacements such as miscible or immiscible gas injection. The details of this effect mean that any uncertainty quantification study should account for the interaction of the near-well heterogeneity and the grid orientation effect. We also show how two possible methods—a well-sponge method and a local embedding technique—are able to produce a solution largely independent of grid orientation for single phase two-component miscible flow. Both methods are versatile in that they can be implemented on general grid topologies. They are illustrated on Cartesian grids for both the standard quarter five spot problem with two different grid orientations, and a problem with a single injection well and two producing wells at different angles to the grid lines. Our results show that it is possible to reduce grid-orientation effects for challenging adverse mobility ratio miscible displacements with only local treatments around the injection wells.     

Self-organized domain microstructures in a plate-like particle suspension subjected to rapid simple shear

The evolution of the microstructure and rheological properties of plate-like particle suspensions subjected to rapid simple shear is studied numerically. In response to the shear-induced strain, particles in the suspensions rearrange to form a steady-state microstructure, and the suspension viscosity reaches a steady value. Under this condition, the microstructure is composed of two domains having different particle fractions and particle orientations. In the matrix (particle-poor) and cluster (particle-rich) domains, the particles’ long axes are oriented subparallel to the shear plane and normal to the maximum compressive principal direction, respectively. A higher particle concentration and friction coefficient enhance the development of cluster domains relative to matrix domains leading the intensity of the preferred particle orientation to decrease and the number of contacting particles, the aspect ratio of clusters, the inter-particle force, and the suspension viscosity to increase. The domain microstructure is governed by two factors: (1) geometric relations between the particle orientation and the maximum compressive axes and (2) the magnitude of particle–fluid and particle–particle interactions. The first factor results in the coupling of the particle orientation and the local fraction of particles, which is an important character of the domain microstructure. The second factor controls the relative development of the cluster and matrix domains through the change in the particles’ rotational behavior. Our results suggest that the microstructure of plate-like suspensions subjected to rapid shear is predictable in terms of the cluster stability, which has important implications for the kinematics of flow-related microstructures in nature and manufacturing.     

Structure of supersonic inviscid nonsymmetric conical flow around a V-wing

The shock wave structure of flow around a V-wing and its properties determining the conical flow topology are numerically investigated within the framework of the inviscid gas model on a wide range of the angles of attack and yaw when in the disturbed supersonic flow either nonsymmetric Mach interaction between the shocks attached to the leading edges of the wing or a shockless flow in the compressed layer on the windward cantilever is realized. The subranges of the angles of attack and yaw with the disturbed flow properties characteristic of the wing of the given geometry are determined. It is found that at high angles of attack, when the branching point of the bow shock beneath the leeward cantilever generates an intense contact discontinuity, the structure of the conical flow in the shock layer on the windward cantilever involves a singularity of a new type which can be characterized as a “vortical” Ferri singularity. It is located above the point of convergence of the streamlines proceeding from the leading edges of the wing, at the vertex of the corresponding contact discontinuity. Flow patterns with the point of convergence of the streamlines proceeding from the leading edges located in the elliptical flow region, which is placed at a local maximum of the pressure distribution over the surface are also found. The range of the angles of attack and yaw on which this new property of supersonic conical flows is realized in the presence of a branched shock system is determined.     

Bonded half planes containing an arbitrarily oriented crack

The plane elastostatic problem for two bonded half planes containing an arbitrarily oriented crack in the neighborhood of the interface is considered. Using Mellin Transforms the problem is formulated as a system of singular integral equations. The equations are solved for various crack orientations, material combinations, and external loads. The numerical results given in the paper include the stress intensity factors, the strain energy release rates, and the probable cleavage angles giving the direction of crack propagation.     

Coupled simulation of the spherical angles of rigid fibres by using a fibre orientation probability distribution model

We study numerically orientation of rigid fibres in a turbulent flow. Instead of a Lagrangian approach, i.e. following fibres along trajectories and analyzing orientations after simulations, we are using a fibre orientation probability distribution (FOPD) model. The mathematical model predicting the two spherical orientation angles (in a three-dimensional computational domain) is presented and solved in a coupled manner, which is the novelty of our approach. Furthermore, numerical predictions are validated with experiments. The origin of the problem is in the papermaking industry, that is, distributing and orienting fibres optimally into a paper sheet produced.     

Adaptive reorientation of cardiac myofibers: The long-term effect of initial and boundary conditions

On the basis of results from modeling and experimental studies it has been hypothesized that myocytes adapt their orientation to achieve a preferred mechanical load. In a previous computational model study in which fiber reorientation was considered as a local response to local fiber cross-fiber shear strain, we have shown that predicted left ventricular (LV) myofiber orientations agreed well with experimental data. In this study, we investigated in the latter model the effect of initial and boundary conditions on predicted fiber orientations on the long term. After adaptation, predicted fiber orientation and deformation became more realistic, irrespective of initial and boundary conditions. As adaptation proceeded, the effect of initial conditions was found to disappear, suggesting that one single optimal fiber orientation field exists for the heart. In contrast, the effect of the boundary conditions persisted, indicating that modeling of in particular the interaction between myocardium and valvular annulus is relevant for predicting LV myofiber reorientation.     

Orientation tensors in simple flows of dilute suspensions of non-Brownian rigid ellipsoids,comparison of analytical and approximate solutions

General analytical solutions are obtained for the planar orientation structure of rigid ellipsoid of revolutions subjected to an arbitrary homogeneous flow in a Newtonian fluid. Both finite and infinite aspect ratio particles are considered. The orientation structure is described in terms of two-dimensional, time-dependent tensors that are commonly employed in constitutive equations for anisotropic fluids such as fiber suspensions. The effect of particle aspect ratio on the evolution of orientation structure is studied in simple shear and planar elongational flows. With the availability of analytical solutions, accuracies of quadratic closure approximations used for nonhomogeneous flows are analyzed, avoiding numerical integration of orientation distribution function. In general, fourth-order orientation evolution equations with sixth-order quadratic closure approximations yield more accurate representations compared to the commonly used second-order evolution equations with fourth-order quadratic closure approximations. However, quadratic closure approximations of any order are found to give correct maximum orientation angle (i.e., preferred direction) results for all particle aspect ratios and flow cases.     

Rotational dynamics of bottom-heavy rods in turbulence from experiments and numerical simulations

We successfully perform the three-dimensional tracking in a turbulent fluid flow of small axisymmetrical particles that are neutrally-buoyant and bottom-heavy, i.e., they have a non-homogeneous mass distribution along their symmetry axis. We experimentally show how a tiny mass inhomogeneity can affect the particle orientation along the preferred vertical direction and modify its tumbling rate. The experiment is complemented by a series of simulations based on realistic Navier–Stokes turbulence and on a point-like particle model that is capable to explore the full range of parameter space characterized by the gravitational torque stability number and by the particle aspect ratio. We propose a theoretical perturbative prediction valid in the high bottom-heaviness regime that agrees well with the observed preferential orientation and tumbling rate of the particles. We also show that the heavy-tail shape of the probability distribution function of the tumbling rate is weakly affected by the bottom-heaviness of the particles.     

刚杆状聚合物和混合层中拟序结构的相互影响

用港方法（即Galerkin方法）求解多球刚杆模型位型空间分布函数的扩散方程，用拟谱方法求解物理空间的N－S方程，同时得到了加入刚杆状聚合物分子的二维混合层中拟序结构的演变过程和聚合物分子的最可几取向．计算结果较为合理地描述了聚合物分子和拟序结构的相互影响．||关键词##4多球刚杆分子模型;;混合层;;拟序结构;;谱方法     

Numerical study of flows of complex fluids between eccentric cylinders using transformation functions

In this paper, we investigate fluid flows between eccentric cylinders by means of two stream‐tube analyses. The first method considers a one‐to‐one global transformation function that allows the physical domain to be transformed into a mapped domain, used as computational domain, that involves concentric streamlines. The second approach uses local transformations and domain decomposition techniques to deal with mixed flow regimes. Both formulations are particularly adapted for handling time‐dependent constitutive equations, since particle‐tracking problems are avoided. Mass conservation is verified in both formulations and the relevant numerical procedure can be carried out using simple meshes built on the mapped streamlines. Fluids obeying anelastic and viscoelastic constitutive equations are considered in the calculations. The numerical results are consistent with those in the literature for the flow rates tested. Application of the method to the K‐BKZ memory‐integral constitutive equation highlights significant differences between the model predictions and those provided by more simple rheological models. Copyright © 2002 John Wiley & Sons, Ltd.     

Continuum-mechanical, Anisotropic Flow model for polar ice masses, based on an anisotropic Flow Enhancement factor

A complete theoretical presentation of the Continuum-mechanical, Anisotropic Flow model, based on an anisotropic Flow Enhancement factor (CAFFE model) is given. The CAFFE model is an application of the theory of mixtures with continuous diversity for the case of large polar ice masses in which induced anisotropy occurs. The anisotropic response of the polycrystalline ice is described by a generalization of Glen’s flow law, based on a scalar anisotropic enhancement factor. The enhancement factor depends on the orientation mass density, which is closely related to the orientation distribution function and describes the distribution of grain orientations (fabric). Fabric evolution is governed by the orientation mass balance, which depends on four distinct effects, interpreted as local rigid body rotation, grain rotation, rotation recrystallization (polygonization) and grain boundary migration (migration recrystallization), respectively. It is proven that the flow law of the CAFFE model is truly anisotropic despite the collinearity between the stress deviator and stretching tensors.     

Measurement of local size and velocity probability density distributions in two-phase suspension flows by Laser-Doppler technique

A technique is presented for the simultaneous measurement of the local number and velocity probability densities of a dilute two-phase suspension which has a distribution of particle sizes and a predominate direction of flow orientation such as in the cases of pipe and boundary-layer flows. It is shown that by a suitable scheme of discrimination on the amplitude as well as the residence time and frequency of the individual Laser-Doppler bursts, one can obtain the statistics on the size number density distribution and, for each size range, velocity distribution of the particulate phase together with the velocity probability distribution of the fluid phase.Results have been obtained for experiments conducted on a laminar uniform flow and a turbulent shear flow of a dilute glass particle-water suspension having a particle size distribution. Calibration needed for the scheme was accomplished by analyzing particle size and number density distribution data obtained from a Coulter particle sizing counter on a sample taken with an isokinetic probe.     

紊动流场中悬浮颗粒分布的随机理论

通过分析固体颗粒在紊动流场中的随机运动,建立了二维流场中垂直于时均流动的方向上颗粒随机位移的概率密度分布函数所满足的方程。由该方程解出的分布函数在一定条件下即相当于颗粒浓度分布函数。运用这一方法研究了[1]、[2]中报道的壁面附近颗粒浓度降低的现象。     

Construction of three-dimensional images of flow structure via particle tracking techniques

Construction of three-dimensional images of flow structure, based on the quantitative velocity field, is assessed for cases where experimental data are obtained using particle tracking technique. The experimental data are in the form of contiguous planes of particle images. These contiguous data planes are assumed to correspond to successive spatial realizations in steady flow, or to phase-referenced realizations in an unsteady flow.Given the particle images on contiguous planes, the in-plane velocity fields are determined. Then, the out-of-plane velocity field is obtained using a spectral interpolation method. Application of this method allows, in principle, construction of the three-dimensional vorticity field and the streamline patterns.A critical assessment is made of the uncertainties arising from the in-plane interpolation of the velocity field obtained from particle tracking and the evaluation of the out-of-plane velocity component. The consequences of such uncertainties on the reconstructed vorticity distributions and streamline patterns are addressed for two basic types of vortex flows: a columnar vortex, for which the streamlines are not closed and are spatially periodic in the streamwise direction; and for a spherical (Hill's) vortex exhibiting closed streamline patterns, and no spatial periodicity.