Extended irreversible thermodynamics and its relation with other continuum approaches

We present an introduction to extended irreversible thermodynamics (EIT) as applied to polymer solutions in the presence of shear flow and of diffusion flux. We discuss with special attention the definition of chemical potential in non-equilibrium situations and its use in the analysis of shear-induced phase transitions. In the second part, we compare EIT with other contemporary continuum approaches: theories with internal variables, the GENERIC approach, and the matrix model. All these theories share an emphasis on the relations between dynamics and thermodynamics at a deeper level than in the classical theory, but each of them has some peculiar advantage in the analysis of some specific aspects of physical problems.     

Shear localisation in interfacial particle layers and its influence on Lissajous-plots

Interfacial rheological measurements often show in their nonlinear Lissajous-plots rhombus or saddle-like shapes indicating complex local deformation behaviour. A strong interacting silica particle and an almost not interacting clay particle were studied in respect to their interfacial rheological properties. Large amplitude oscillation shear measurements were performed with a bicone geometry and combined with optical measurements, from which particle tracks were calculated. A correlation was found between the appearance of shear localisation and Lissajous-plot shapes. Silica particles showed shear localisation at the bicone edge and rhombic plateaus in the Lissajous-plot, while the shear localisation for the clay particles was observed at the cup’s wall as saddle-like shaped Lissajous-plots     

Non-local continuum model and its numerical simulation for localization problem

A non-local continuum model for strain-softening simply taking plastic strain or damage variable as a non-local variable is derived by using the additive decomposition principle of finite deformation gradient. At the same time, variational equations, their finite element formulations and numerical convoluted integration algorithm of the model in current configuration usually called co-moving coordinate system are given. Stability and convergence of the model are proven by means of the weak convergence theorem of general function and the convoluted integration theory. Mathematical and physical properties of the characteristic size for material or structure are accounted for within the context of a statistical weighted or kernel function, and way is investigated. Numerical simulation shows that this model is suitable for to analyzing deformation localization problems. The project supported by the National Natural Science Foundation of China (No. 19632030).     

Accurate 3D viscous incompressible flow calculations with the FEM

A second-order-accurate (in both time and space) formulation is developed and implemented for solution of the three-dimensional incompressible Navier–Stokes equations involving high-Reynolds-number flows past complex configurations. For stabilization, only a fourth-order-accurate artificial dissipation term in the momentum equations is used. The finite element method (FEM) with an explicit time-marching scheme based on two-fractional-step integration is used for solution of the momentum equations. The element-by-element (EBE) technique is employed for solution of the auxiliary potential function equation in order to ease the memory requirements for matrix. The cubic cavity problem, the laminar flow past a sphere at various Reynolds numbers and the flow around the fuselage of a helicopter are successfully solved. Comparison of the results with the low-order solutions indicates that the flow details are depicted clearly even with coarse grids. © 1997 John Wiley & Sons, Ltd.     

曲线坐标系下平面二维浅水模型的修正与应用

借鉴有关弯道水流流速分布的研究成果,计入弯道环流引起的横向的动量交换,对曲线坐标系平面二维浅水方程做了修正;采用边界处正交的曲线网格生成技术,处理复杂的计算区域边界。采用修正后的模型对90°弯道水流进行了数值模拟,并与原模型的计算结果及实测资料进行了比较,结果表明该模型能够有效地模拟流线弯曲的复杂水流的水力特性。     

Converging Bounds for the Effective Shear Speed in 2D Phononic Crystals

Calculation of the effective quasistatic shear speed c in 2D solid phononic crystals is analyzed. The plane-wave expansion (PWE) and the monodromy-matrix (MM) methods are considered. For each method, the stepwise sequence of upper and lower bounds is obtained which monotonically converges to the exact value of c. It is proved that the two-sided MM bounds of c are tighter and their convergence to c is uniformly faster than that of the PWE bounds. Examples of the PWE and MM bounds of effective speed versus concentration of high-contrast inclusions are demonstrated.     

Characterization of the crossover from capillary invasion to viscous fingering to fracturing during drainage in a vertical 2D porous medium

We experimentally studied the displacement of a viscous wetting fluid (water) by an inviscid non-wetting fluid (air) injected at the bottom of a vertical Hele-Shaw cell filled with glass microbeads. In order to cover a wide parameter space, the permeability of the porous medium was varied by using different bead size ranges and diverse air flow rates were generated by means of a syringe pump. A LED light table was used to back illuminate the experimental cell, allowing a high speed camera to capture images of the drainage process at equal time intervals. The invasion occurred in intermittent bursts. Image processing of the bursts and fractal analysis showed successive transitions from capillary invasion to viscous fingering to fracturing during the same experiment, dependent on the medium permeability, the air injection flow rate, and the vertical position in the cell. The interplay between the capillary, viscous and gravity forces determines the nature of the invasion pattern and the transitions, from capillary invasion to viscous fingering with decreasing fluid pressure on one hand and from viscous fingering to fracturing with decreasing effective overburden pressure on the other hand.     

A numerical approach to model non-isothermal viscous flow through fibrous media with free surfaces

A numerical simulation is presented to predict the free surface and its interactions with heat transfer and cure for flow of a shear-thinning resin through the fibre preform the flow part of the simulation is based on the finite element/control volume method. Since the traditional control volume approach produces an error associated with a mass balance inconsistency, a new method which overcomes this issue is proposed, the element control volume method. The heat transfer and cure analysis in the simulation are based on the finite difference/control volume method. Since heat conduction is dominant in the through-thickness direction and most of the heat convection is in-plane, heat transfer and cure are solved in fully three-dimensional form. A simple concept of the boundary condition constant is introduced which models a realistic mould configuration with a heating element located at a distance behind the mould wall. The varying viscosity throughout the mould associated with the strain rate, temperature and degree of cure distribution may be accounted for in calculating the mould-filling pattern. This introduces a two-way coupling between momentum and energy transport in fibrous media during mould filling.     

Numerical investigation from rarefied flow to continuum by solving the Boltzmann model equation

Based on the Bhatnagar–Gross–Krook (BGK) Boltzmann model equation, the unified simplified velocity distribution function equation adapted to various flow regimes can be presented. The reduced velocity distribution functions and the discrete velocity ordinate method are developed and applied to remove the velocity space dependency of the distribution function, and then the distribution function equations will be cast into hyperbolic conservation laws form with non‐linear source terms. Based on the unsteady time‐splitting technique and the non‐oscillatory, containing no free parameters, and dissipative (NND) finite‐difference method, the gas kinetic finite‐difference second‐order scheme is constructed for the computation of the discrete velocity distribution functions. The discrete velocity numerical quadrature methods are developed to evaluate the macroscopic flow parameters at each point in the physical space. As a result, a unified simplified gas kinetic algorithm for the gas dynamical problems from various flow regimes is developed. To test the reliability of the present numerical method, the one‐dimensional shock‐tube problems and the flows past two‐dimensional circular cylinder with various Knudsen numbers are simulated. The computations of the related flows indicate that both high resolution of the flow fields and good qualitative agreement with the theoretical, DSMC and experimental results can be obtained. Copyright © 2003 John Wiley & Sons, Ltd.     

Simplest differential equation of stock price,its solution and relation to assumption of black-scholes model

1　TheSituationofMathematicalAnalysisofStockMarketandTwoAnalyticSystems　　Mostpublicationsunderthetitleofstockmarket,includinghundredsofauthoritativeguidetocontemporaryinvestmenttheories[1],belongtodatacollection ,summinguptheexperienceofmakingprofits,operatingskills,qualitativeanalysis,andtechnicalanalysisusingsimplemathematics,etc .Thework ,cannotbeconsidered ,ingeneral,asrulequantitativeinvestigationandisstillinadescribingstage.Theliteraturesonbehaviorofstockmarketinvestigatedbyanalytic…     

Use of B-spline surface to model large-deformation continuum plates: procedure and applications

The absolute nodal coordinate formulation (ANCF) has been used in the analysis of large deformation of flexible multibody systems that encompass belt drive, rotor blade, and cable applications. As demonstrated in the literature, the ANCF finite elements are ideal for isogeometric analysis. The purpose of this investigation is to establish a relationship between the B-splines, which are widely used in the geometric modeling, and the ANCF finite elements in order to construct continuum models of large-deformation geometries. This paper proposes a simplified approach to map the B-spline surfaces into ANCF thin plate elements. Matrix representation of the mapping process is established and examined through numerical examples successfully. The matrix representation of the mapping process is used because of its suitability of computer coding and to minimize the calculation time. The error estimation is carried out by analyzing the gap between the points of each ANCF element and the corresponding points of the portion of the B-spline surface. The Hausdorff distance is used to study the effect of the number of control points, the degree of interpolation, and the knot multiplicity on the mapped geometry. It is found that cubic interpolation is recommended for optimizing the accuracy of mapping the B-spline surface to ANCF thin plate elements. It is found that thin plate element in ANCF missing a number of basis functions which considered a source of error between the two surfaces, as well as it does not allow to converting the ANCF thin plate elements model to B-spline surface. In this investigation, an application example of modeling large-size wind turbine blade with uniform structure is illustrated. The use of the continuum plate elements in modeling flexible blades is more efficient because of the relative scale between the plate thickness and its length and width and the high flexibility of its structure. The numerical results are compared with the results of ANSYS code with a good agreement. The dynamic simulation for mapped surface model shows a numerical convergence, which ensures the ability of using the proposed approach for applications of dynamics for design and computer-aided design.     

Structure of kinematic and force fields in the Riemannian continuum model

In this paper, we consider a non-Euclidean continuum model for which the structure of defects in the material is characterized by an internal metric and scalar curvature. It is shown that the irrotational displacement field for points of this medium is composed of elastic displacements (in the absence of defects) and the field which characterizes the deviation of the internal geometry of the model from Euclidean geometry. The corresponding components of the internal stresses are the sum of elastic stresses and the self-equilibrated stresses determined by the scalar curvature. The exact solution for the vortex field of dislocations is constructed, and conditions of the existence of a nonzero stress field parametrized by a scalar curvature in the absence of external forces are formulated.     

A continuum traffic flow model with the consideration of coupling effect for two-lane freeways

A new higher-order continuum model is proposed by considering the coupling and lane changing effects of the vehicles on two adjacent lanes. A stability analysis of the proposed model provides the conditions that ensure its linear stability. Issues related to lane changing, shock waves and rarefaction waves, local clustering and phase transition are also investigated with numerical experiments. The simulation results show that the proposed model is capable of providing explanations to some particular traffic phenomena commonly observable in real traffic flows.     

Amplitude modulation and its relation to streamwise convection velocity

The paper is devoted to show the relation between the amplitude modulation and the convection velocity in a turbulent boundary layer. The analysis has been performed based upon velocity profiles measured with the hot-wire technique in a turbulent boundary layer with pressure gradient corresponding to turbomachinery conditions. The VITA detection method of small-scale structures was used in order to show the impact of the pressure gradient on phase averaged time-traces changes in the aspect of amplitude modulation. The physical significance of the large scales influence on the small scales, formerly known as the amplitude modulation, has been elucidated. It was demonstrated that the measure of amplitude modulation can be used as an indicator of the convection velocity of coherent structures. The convection velocity profile can be described in the inner region using the amplitude modulation term when the proper scale is applied. The convection velocity prediction in favourable and adverse pressure gradients shows that t the strongest impact of the pressure gradient occurs in the buffer layer. The stronger increase occurs in the buffer layer. It was also found that the skewness factor of phase-averaged VITA waveforms, is very similar in shape with the cross-product term of the streamwise skewness factor, which quantifies amplitude modulation. It indicates that the well-known modification of quadrant events in pressure gradient flows is a result of the amplitude modulation phenomenon.     

多自由度转子-轴承系统周期运动分岔及稳定性研究

建立考虑诸多因素的转子-轴承系统多自由度模型,将与Newmark结合的打靶法应用到多自由度转子-轴承系统的周期稳定性分析中。着重研究了转子-轴承系统失稳转速随系统偏心量、轴承间隙、润滑油动力粘度以及轴承长径比的变化规律,研究结果表明:提高系统偏心量、减小轴承间隙、增大润滑油动力粘度以及选择适当的轴承长径比均能提高转子-轴承系统的失稳转速;对于不同的参数值,系统表现出不同的分岔规律,系统发生半速涡动时表现为倍周期分岔或拟周期分岔,发生油膜振荡时则表现为拟周期分岔。     

Explicit formulations for evaluation of velocity gradients using boundary‐domain integral equations in 2D and 3D viscous flows

In this paper, explicit boundary‐domain integral equations for evaluating velocity gradients are derived from the basic velocity integral equations. A free term is produced in the new strongly singular integral equation, which is not included in recent formulations using the complex variable differentiation method (CVDM) to compute velocity gradients (Int. J. Numer. Meth. Fluids 2004; 45 :463–484; Int. J. Numer. Meth. Fluids 2005; 47 :19–43). The strongly singular domain integrals involved in the new integral equations are accurately evaluated using the radial integration method (RIM). Considerable computational time for evaluating integrals of velocity gradients can be saved by using present formulation than using CVDM. The formulation derived in this paper together with those presented in reference (Int. J. Numer. Meth. Fluids 2004; 45 :463–484) for 2D and in (Int. J. Numer. Meth. Fluids 2005; 47 :19–43) for 3D problems constitutes a complete boundary‐domain integral equation system for solving full Navier–Stokes equations using primitive variables. Three numerical examples for steady incompressible viscous flow are given to validate the derived formulations. Copyright © 2007 John Wiley & Sons, Ltd.     

Second law analysis of the 2D laminar flow of two-immiscible,incompressible viscous fluids in a channel

The laminar and fully developed flows of two immiscible fluids confined in a thin slit of constant wall heat fluxes are analyzed in terms of entropy generations due to irreversibility of forced convection heat transfer. The governing equations are analytically derived using expressions for velocity distributions. The derived equation for the dimensionless entropy generation number is used to interpret the relative importance of frictions to conduction by varying irreversibility distribution ratio ϕ. It is found that the minimum entropy generation takes place at the dimensionless half transverse distance (ξ) of 0.3 for values of ϕ higher than zero. The entropy generation near the plate increases more rapidly in fluid I than in fluid II as viscous dissipation effects becomes more important. The velocity profiles are found to be in agreement with the distributions of the dimensionless entropy generation number (N _S ) for two-immiscible incompressible flows in the slit.