The prediction of spatially periodic flows using a finite‐volume model

A periodic boundary condition has been developed that can be used in conjunction with a specified flow rate to produce accurate results in spatially periodic geometries. This condition is useful in situations where the flow rate is known, or more importantly, in cases where the pressure gradient is not known a priori, such as in countercurrent flows. Using the present condition, the flow rate is imposed at the inlet in terms of a bulk velocity, but the velocity field evolves as part of the solution. The condition is formulated to be suitable for both fixed and moving periodic domains. For the case of a moving domain, a correction is introduced to account for changes in the instantaneous velocity through the periodic edges. Under periodic conditions, these corrections integrate to zero over a complete (temporal) period. The new periodic condition is shown to produce accurate results for flat and wavy‐walled channels under both induced flow and countercurrent conditions. Copyright © 2003 John Wiley & Sons, Ltd.     

Macromolecular deformation in periodic extensional flows

The response of a dumbbell model for dilute polymer solutions is examined for periodic sequences of homogeneous flows chosen to approximate the unsteady kinematics appropriate for flow through porous media. The evolution of particle shape depends on the kinematic history of the flow and on microrheological properties of the dumbbell such as a variable friction factor and a nonlinear spring. The effect of vorticity on macromolecular stretching is found to differ qualitatively from results for steady flows, and entrance effects that persist over timescales much larger than the intrinsic macromolecular relaxation time are observed.     

On the interpretation of data from Converging Flow Rheometers

A reappraisal of data obtained from a Converging Flow Rheometer (CFR) is presented, together with new results for a specific polymer solution.Particular emphasis is placed on the interpretation of the experimental pressure data in terms of a planar extensional viscosity. It is suggested that previous interpretations, while yielding viscosity levels that appear reasonable, nevertheless fail to give the qualitative behaviour that might be expected on the basis of predictions from well accepted constitutive models. This, in the authors' opinions, arises because certain fluids, i.e. those that are highly tension-thickening, cease to flow in accordance with the assumed kinematics at high flow rates.By adapting a recently proposed approximate analysis for flow through a contraction it is shown that better qualitative behaviour, for the planar extensional viscosity, can be obtained from the Converging Flow Rheometer.     

Periodic,almost periodic and chaotic motions of forced self-sustained oscillators

Steady motions of the Van der Pol oscillator and an oscillator with hysteresis are studied numerically in this paper. Some features of periodic, almost periodic and chaotic motions of forced self-sustained oscillators are investigated. This paper has been presented at the ICTAM XVI Lyngby.     

Theoretical and numerical investigation of HF elastic wave propagation in two-dimensional periodic beam lattices

The elastic wave propagation phenomena in two-dimensional periodic beam lattices are studied by using the Bloch wave transform. The numerical modeling is applied to the hexagonal and the rectangular beam lattices, in which, both the in-plane （with respect to the lattice plane） and out-of-plane waves are considered. The dispersion relations are obtained by calculating the Bloch eigenfrequencies and eigenmodes. The frequency bandgaps are observed and the influence of the elastic and geometric properties of the primitive cell on the bandgaps is studied. By analyzing the phase and the group velocities of the Bloch wave modes, the anisotropic behaviors and the dispersive characteristics of the hexagonal beam lattice with respect to the wave prop- agation are highlighted in high frequency domains. One im- portant result presented herein is the comparison between the first Bloch wave modes to the membrane and bend- ing/transverse shear wave modes of the classical equivalent homogenized orthotropic plate model of the hexagonal beam lattice. It is shown that, in low frequency ranges, the homog- enized plate model can correctly represent both the in-plane and out-of-plane dynamic behaviors of the beam lattice, its frequency validity domain can be precisely evaluated thanks to the Bloch modal analysis. As another important and original result, we have highlighted the existence of the retro- propagating Bloch wave modes with a negative group veloc- ity, and of the corresponding ＂retro-propagating＂ frequency bands.     

The possibility of measuring flow velocity using a low-intensity periodic pulse discharge

A direct method of measuring the flow velocity in a supersonic wind tunnel by creating a low-intensity periodic pulse discharge is proposed.     

Assessment of nonlinear strain measures for extensional and shearing flows of polymer melts

From stress-strain experiments in extensional and shearing flows, nonlinear strain measures and effective damping functions are derived for a polyisobutylene melt. The strain measures determined in planar extensional flow and in simple shear flow coincide. Experimental results are compared with predictions of two molecular theories, the Doi-Edwards model and the molecular stress function approach of Wagner and Schaeffer. Discrepancies between theories and experiment lead to a reconsideration of the classification of extensional flows. The symmetry of the flow field is identified and quantified as an important parameter influencing the strain measure, and a unifying strain measure for general extensional and shearing flows of polymer melts is presented.     

Numerical investigation of viscous fingering in Hele-Shaw cell with spatially periodic variation of depth

Viscous fingering in a modified Hele-Shaw cell is numerically investigated. The cell allows periodic variation of depth in the lateral direction. The wavenumber n of the depth perturbation has great influence on fingering patterns. For n = 1, the fingering pattern due to the interface instability remains the same as that in the conventional Hele-Shaw cell, while the depth variation causes the steady finger to be a little narrower. For n = 2, four different fingering patterns are captured, similar to the available experimental observations in a modified Hele-Shaw cell containing a centered step-like occlusion. It is found that new fingering patterns appear as n further increases, among which, two patterns with spatial oscillation along both edges of the finger are particularly interesting. One is a symmetric oscillatory finger for n = 3, and the other is an asymmetric one for n = 4. The influence of capillary number on fingering patterns is studied for n = 3 and 4. We find that spatial oscillation of the finger nearly ceases at moderate capillary numbers and occurs again as the capillary number increases further. Meanwhile, the wide finger shifts to the narrow one. It is accompanied by a sudden decrease in the finger width which otherwise decreases continuously as the capillary number increases. The wavenumber and the amplitude of depth perturbation have little effect on the finger width.     

Extension of a combined analytical/numerical initial value problem solver for unsteady periodic flow

Here we describe analytical and numerical modifications that extend the Differential Reduced Ejector/ mixer Analysis (DREA), a combined analytical/numerical, multiple species ejector/mixing code developed for preliminary design applications, to apply to periodic unsteady flow. An unsteady periodic flow modelling capability opens a range of pertinent simulation problems including pulse detonation engines (PDE), internal combustion engine ICE applications, mixing enhancement and more fundamental fluid dynamic unsteadiness, e.g. fan instability/vortex shedding problems. Although mapping between steady and periodic forms for a scalar equation is a classical problem in applied mathematics, we will show that extension to systems of equations and, moreover, problems with complex initial conditions are more challenging. Additionally, the inherent large gradient initial condition singularities that are characteristic of mixing flows and that have greatly influenced the DREA code formulation, place considerable limitations on the use of numerical solution methods. Fortunately, using the combined analytical–numerical form of the DREA formulation, a successful formulation is developed and described. Comparison of this method with experimental measurements for jet flows with excitation shows reasonable agreement with the simulation. Other flow fields are presented to demonstrate the capabilities of the model. As such, we demonstrate that unsteady periodic effects can be included within the simple, efficient, coarse grid DREA implementation that has been the original intent of the DREA development effort, namely, to provide a viable tool where more complex and expensive models are inappropriate. Copyright © 2002 John Wiley & Sons, Ltd.     

树脂基三维编织复合材料粘弹性性能的有限元分析

根据材料的细观结构,采用APDL语言分别建立了纤维束和三维编织复合材料两级单胞的参数化几何模型;推导了Prony级数表示的树脂粘弹性本构方程,对模型进行了组分材料参数设置;对纤维束单胞模型进行扫掠式网格划分,对三维编织复合材料单胞模型进行线-面-体式网格划分;对两级单胞模型均施加合理的边界条件,使单胞边界上的位移满足周期性和连续性。以有限元模型为基础,计算了三维编织复合材料的粘弹性能,并给出了材料粘弹性效应随工艺参数变化的规律。计算结果表明:三维编织复合材料编织方向的粘弹性效应随编织角的增大而增强,随纤维体积比的增大而减弱。该结果与已有实验结论一致。     

Implicit solution of time spectral method for periodic unsteady flows

The present paper investigates the implicit solution of time spectral model for periodic unsteady flows. In the time spectral model, the physical time derivative is approximated using spectral method. The robustness issues associated with implicit solution of time spectral model are analyzed and validated by numerical results. It is found that spectral approximation of the time derivative weakens the diagonal dominance property of the Jacobian matrix, resulting in the deterioration of stability and convergence speed. In this paper we propose to solve the coupled governing equations implicitly using multigrid preconditioned generalized minimal residual (GMRES) method, which demonstrates favorable convergence speed. Also it is demonstrated that the current method is insensitive to the variations of frequency and number of harmonics. Comparison of computation results with dual time step unsteady computation validates the high efficiency of the current method. Copyright © 2009 John Wiley & Sons, Ltd.     

Simulation of flow through porous anode in MFC at higher power density

Microbial fuel cell (MFC) is a novel environmental friendly energy device which has received great attention due to its technology for producing electricity directly from organic or inorganic matter by using bacteria as catalyst. To date, many experiments have been carried out to achieve the maximum power output with advective flow through porous anode to the cathode in the MFC. However, the precise mechanical mechanism of flow through anode and the quantified relationship between electrode spacing and MFC performance are not yet clearly understood. It has been found experimentally that the power output can be increased apparently at certain electrode spacing configuration. Based on these available experimental data, this paper investigates the effect of spacing between electrodes, the Darcy number of porous anode and the Reynolds number on the power production performance of MFC by using lattice Boltzmann method. The numerical simulation results present that the distance between electrodes significantly influences the flow velocity and residence time of the organic matter attached to the anode in the MFC. Moreover, it is found that the Darcy number of porous anode and the Reynolds number can regulate the output efficiency of MFC. These results perform better understanding of the complex phenomena of MFC and will be helpful to optimize MFC design.     

多孔格栅均匀化模型平压仿真分析

为了研究均匀化方法在一种多孔格栅结构中的应用,从格栅单胞尺度入手,建立了一种适用于有限元仿真分析的三维周期性边界条件.以ABAQUS作为分析平台,对周期性边界条件下的格栅单胞模型进行了平压仿真分析,并将仿真结果与文献实验结果对比,验证了该边界条件的可靠性.利用均匀化理论建立了格栅单胞力学平衡方程,得到了格栅均匀化模型....     

纺织复合材料和结构多尺度耦合的数值分析

研究了纺织复合材料和结构多尺度耦合的数值分析模型。建立了微、细观单胞,给出了纺织复合材料平均弹性常数的逐级分析方法,着重研究了由宏观结构、到细观纤维束、再到微观纤维三个尺度耦合的应力分析方案。对于常用的板壳状纺织复合材料结构,在面内载荷下,假设每层细观单胞的平均面内应变是一致的,在弯曲、横向剪切及扭曲等非面内载荷下,在内力等效条件下将沿厚度方向连续分布的宏观应力简化为阶梯状分布,忽略了每层细观单胞范围内宏观应力沿厚度方向的梯度变化,由此利用细观单胞模型实现宏观应力与细观应力之间的传递,再利用微观单胞可得到纤维尺度的微观应力。最后以一种三维机织复合材料为例,用上述多尺度耦合的模型逐级分析了材料的平均弹性常数,并沿相反方向,由宏观结构分析逐级计算出纤维束尺度和纤维尺度的细、微观应力的局部波动。     

A Lagrangian criterion of unsteady flow separation for two-dimensional periodic flows

The present paper proposes a Lagrangian criterion of unsteady flow separation for two-dimensional periodic flows based on the principle of weighted averaging zero skin-friction given by Haller(HALLER, G. Exact theory of unsteady separation for two-dimensional flows. Journal of Fluid Mechanics, 512, 257–311(2004)). By analyzing the distribution of the finite-time Lyapunov exponent(FTLE) along the no-slip wall, it can be found that the periodic separation takes place at the point of the zero FTLE.This new criterion is verified with an analytical solution of the separation bubble and a numerical simulation of lid-driven cavity flows.     

A three‐dimensional vortex particle‐in‐cell method for vortex motions in the vicinity of a wall

A new vortex particle‐in‐cell method for the simulation of three‐dimensional unsteady incompressible viscous flow is presented. The projection of the vortex strengths onto the mesh is based on volume interpolation. The convection of vorticity is treated as a Lagrangian move operation but one where the velocity of each particle is interpolated from an Eulerian mesh solution of velocity–Poisson equations. The change in vorticity due to diffusion is also computed on the Eulerian mesh and projected back to the particles. Where diffusive fluxes cause vorticity to enter a cell not already containing any particles new particles are created. The surface vorticity and the cancellation of tangential velocity at the plate are related by the Neumann conditions. The basic framework for implementation of the procedure is also introduced where the solution update comprises a sequence of two fractional steps. The method is applied to a problem where an unsteady boundary layer develops under the impact of a vortex ring and comparison is made with the experimental and numerical literature. Copyright © 2001 John Wiley & Sons, Ltd.