Discontinuous Galerkin finite element method applied to the coupled unsteady Stokes/Cahn-Hilliard equations

Two-phase flows driven by the interfacial dynamics are studied by tracking implicitly interfaces in the framework of the Cahn-Hilliard theory. The fluid dynamics is described by the Stokes equations with an additional source term in the momentum equation taking into account the capillary forces. A discontinuous Galerkin finite element method is used to solve the coupled Stokes/Cahn-Hilliard equations. The Cahn-Hilliard equation is treated as a system of two coupled equations corresponding to the advection-diffusion equation for the phase field and a nonlinear elliptic equation for the chemical potential. First, the variational formulation of the Cahn-Hilliard equation is presented. A numerical test is achieved showing the optimal order in error bounds. Second, the variational formulation in discontinuous Galerkin finite element approach of the Stokes equations is recalled, in which the same space of approximation is used for the velocity and the pressure with an adequate stabilization technique. The rates of convergence in space and time are evaluated leading to an optimal order in error bounds in space and a second order in time with a backward differentiation formula at the second order. Numerical tests devoted to two-phase flows are provided on ellipsoidal droplet retraction, on the capillary rising of a liquid in a tube, and on the wetting drop over a horizontal solid wall.     

刚-柔耦合动力学系统的建模理论研究

刚－柔耦合动力学系统的传统的混合坐标方法是零次近似方法,在建模过程中,直接套用的结构动力学的小变形假设,忽略了变形位移的高次耦合变形量．本文对柔性梁建立较零次近似更精确的高次耦合动力学模型,从连续介质力学理论出发,在变形位移中,计及横向位移引起的轴向缩短,导出变形位移的二次耦合量．用一致质量有限元方法对梁进行离散,基于Jourdain速度变分原理导出大范围运动为自由的柔性梁的刚－柔耦合动力学方程．计算了柔性重力摆的角速度和摆端点的横向变形,揭示零次近似模型和耦合模型的刚－柔耦合动力学性质的根本差异．     

Analytic evaluation of periodic responses of a forced nonlinear oscillator

The periodic responses of a strongly nonlinear, single-degree-of-freedom forced oscillator with weak excitation and damping are examined. The presented methodology is based on a regular perturbation expansion, whose first term is the solution of the unforced, and undamped nonlinear problem. Higher order approximations are computed by explicitly solving linear differential equations possessing a periodically varying coefficient. The general theory is used for studying the periodic steady state motions of the periodically forced system. Moreover, it is shown that the presented analysis can be used to analytically study the orbital stability of the identified steady state motions. The proposed method can also be used for studying periodic responses due to nonperiodic transient forces, provided that these responses are close to the O(1) periodic generating solution.     

A justification of nonlinear properly invariant plate theories

A single asymptotic derivation of three classical nonlinear plate theories is presented in a setting which preserves the frame-invariance properties of three-dimensional finite elasticity. By a successive scaling of the external loading on the three-dimensional body, the nonlinear membrane theory, the nonlinear inextensional theory and the von Kármán equations are derived as the leading-order terms in the asymptotic expansion of finite elasticity. The governing equations of the nonlinear inextensional theory are of particular interest where 1) plane-strain kinematics and plane-stress constitutive equations are derived simultaneously from the asymptotic analysis, 2) the theory can be phrased as a minimization problem over the space of isometric deformations of a surface, and 3) the local equilibrium equations are identical to those arising in the one-director Cosserat shell model. Furthermore, it can be concluded that with a regular, single-scale asymptotic expansion it is not possible to obtain a system of plate equations in which finite membrane strain and finite bending strain occur simultaneously in the leading-order term of an asymptotic analysis.     

准晶数学弹性力学和缺陷力学

对准晶数学弹性理论的基本概念和基本框架作了介绍,在此基础上分别针对目前已经发现的几类一维准晶、二维准晶和三维准晶讨论了其数学弹性的理论体系．为了求解准晶弹性的边值问题或初值一边值问题,还必须发展相应的方法论．物理工作者在研究准晶位错弹性问题中发展了Ｇｒｅｅｎ函数方法．针对一维与二维准晶弹性中几类问题提出了分解与叠加程序,这一程序的使用,使极其复杂的准晶弹性问题得到简化,进而引进位移函数或应力函数,把数目。庞大的准晶弹性基本方程化成一个或少数几个高阶偏微分方程,进一步使求解步骤大为简化．对三维立方准晶弹性也采用了类似步骤使求解过程大为简化．在以上化简的基础上,发展了准晶弹性的边值问题或初值一边值问题的复交函数方法和 Ｆｏｕｒｉｅｒ分析方法,求得了一系列准晶位错问题和裂纹问题的分析解（古典解）．在研究准晶弹性的边值问题古典解的同时,也讨论了同这些边值问题相对应的变分问题和广义解（弱解）以及这种弱解的数值方法──有限元法．在物理学家工作基础上开展的这些工作可以看作对经典数学弹性理论和方法、经典Ｖｏｌｔｅｒｒａ位错理论、普通结构材料断裂力学和经典有限元的某些发展．此外,还把一维六方准晶弹性动力学的结果与统计物理的某些     

Thermal buckling of functionally graded circular plates based on higher order shear deformation plate theory

In this research, thermal buckling of circular plates compose of functionally graded material (FGM) is considered. Equilibrium and stability equations of a FGM circular plate under thermal loads are derived, based on the higher order shear deformation plate theory (3rd order plate theory). Assuming that the material properties vary as a power form of the thickness coordinate variable z and using the variational method, the system of fundamental partial differential equations is established. A buckling analysis of a functionally graded circular plate (FGCP) under various types of thermal loads is carried out and the result are given in closed-form solutions. The results are compared with the critical buckling temperature obtained for FGCP based on first order (1st order plate theory) and classical plate theory (0 order plate theory) given in the literature. The study concludes that higher order shear deformation theory accurately predicts the behavior of FGCP, whereas the first order and classical plate theory overestimates buckling temperature.     

Fluid models for relativistic electron beams

A system of equations is provided that may be used in the study of relativistic charged particle beams. The equations are based upon the equations of the kinetic theory for first, second and third order moments and the system is closed by letting the third order moment depend on the lower order ones. The form of that dependence is formally equal to the explicit constitutive function given by extended thermodynamics. However, here the contributions to the third order moment can be classed as being different in order of magnitude, because there is a smallness parameter characterizing the small dispersion of the particle beam. The resulting system of equations is quite specific, it is fully covariant and it is equivalent to a symmetric hyperbolic system thus ensuring existence and uniqueness of solutions.     

Internal waves in a two‐layer system using fully nonlinear internal‐wave equations

In order to understand the nonlinear effect in a two‐layer system, fully nonlinear strongly dispersive internal‐wave equations, based on a variational principle, were proposed in this study. A simple iteration method was used to solve the internal‐wave equations in order to solve the equations stably. The applicability of the proposed numerical computation scheme was confirmed to agree with linear dispersion relation theoretically obtained from variational principle. The proposed computational scheme was also shown to reproduce internal waves including higher‐order nonlinear effect from the analysis of internal solitary waves in a two‐layer system. Furthermore, for the second‐order numerical analysis, the balance of nonlinearity and dispersion was found to be similar to the balance assumed in the KdV theory and the Boussinesq‐type equations. Copyright © 2009 John Wiley & Sons, Ltd.     

Reduction of dimensions in random,elastic soil medium

The paper deals with an application of the plane strain analysis in a stochastic three-dimensional soil medium. In a framework of random elasticity theory, the geostatical state of stresses and the problem of a unit force acting in a statistically homogeneous half-space are considered. Only the modulus of elasticity is considered to be random and is modelled as a three-dimensional (3-D) homogeneous random field. As the result of imposed constrains due to the plane strain assumption the additional body and surface forces are induced. In order to determine them, additional equations must be introduced. The equations in a form of constrain relations are proposed in this paper. These equations are also valid for a case of uniformly distributed external loading.First, the two-dimensional (2-D) problem and its reduction to the uni-axial strain state, for the gravity forces and uniform, unlimited surface loading is considered. Then, it is generalised into a 2-D schematization of the 3-D state. Next, the problem of a unit force acting in a statistically homogeneous half-space is considered. For a 3-D state of stress and strain the resulting stresses are compared with those for a 2-D state. These stresses for the multidimensional state of strain and stress are presented as a sum of two components. The first one reflects plane strain state stresses and is given in a form of a 3-D random field. This term allows for incorporating a spatial, 3-D soil variability into a two-dimensional analysis. The second component can be treated as a correction term and it represents the longitudinal influence of a 3-D analysis.Some numerical results are presented in this paper. The proposed method can be regarded as a framework for further research aiming at application to a variety of geotechnical problems, for which the plane strain state is assumed.     

Method of relaxed streamline upwinding for hyperbolic conservation laws

In this work a new finite element based Method of Relaxed Streamline Upwinding is proposed to solve hyperbolic conservation laws. Formulation of the proposed scheme is based on relaxation system which replaces hyperbolic conservation laws by semi-linear system with stiff source term also called as relaxation term. The advantage of the semi-linear system is that the nonlinearity in the convection term is pushed towards the source term on right hand side which can be handled with ease. Six symmetric discrete velocity models are introduced in two dimensions which symmetrically spread foot of the characteristics in all four quadrants thereby taking information symmetrically from all directions. Proposed scheme gives exact diffusion vectors which are very simple. Moreover, the formulation is easily extendable from scalar to vector conservation laws. Various test cases are solved for Burgers equation (with convex and non-convex flux functions), Euler equations and shallow water equations in one and two dimensions which demonstrate the robustness and accuracy of the proposed scheme. New test cases are proposed for Burgers equation, Euler and shallow water equations. Exact solution is given for two-dimensional Burgers test case which involves normal discontinuity and series of oblique discontinuities. Error analysis of the proposed scheme shows optimal convergence rate. Moreover, spectral stability analysis gives implicit expression of critical time step.     

瞬变流摩阻计算及摩阻对水力瞬变的影响

摩阻损失的精确计算对于长输管道的水力瞬变分析是很重要的一个环节、过去为减少计算工作量，采用固定摩阻系数和流量一阶近似法计算瞬交流摩阻，精度有限．本跟踪液流流态变化，采用变摩阻系数和流量二阶插值法计算瞬交流摩阻，提高了摩阻计算精度．同时，从理论上分析了管道摩阻对水力瞬变的影响，澄清了一些模糊认识．     

Higher Order Terms in Extended Thermodynamics

Extended Thermodynamics is an elegant and powerful material theory which yields, when applied, a symmetric hyperbolic system of evolution equations for the independent field variables. This system is obtained by imposing general physical principles, such as that of material objectivity and the entropy principle. With respect to thermodynamic equilibrium, the latter is carried through to second order. However, the following question arises: If this principle is imposed to even higher order terms, could it then yield restrictions also on the lower order terms, beyond those previously obtained in the literature? In this paper the entropy principle is applied to fourth order so that constitutive functions up to this order are obtained. In the process of computations many complicated equations involving only the lower order terms are deduced, which must be satisfied as identities. These equations, after long and tedious calculations, turn indeed out to be identically satisfied. This fact cannot be casual, thus confirming that this theory is trustworthy. Received October 28, 1997     

Center manifold and multivariable approximants applied to non-linear stability analysis

This paper presents a research devoted to the study of instability phenomena in non-linear model with a constant brake friction coefficient. This paper outlines the stability analysis and a procedure to reduce and simplify the non-linear system, in order to obtain limit cycle amplitudes. The center manifold approach, the multivariable approximants theory, and the alternate frequency/time domain (AFT) method are applied. Brake vibrations, and more specifically heavy trucks grabbing are concerned. The modelling introduces sprag-slip mechanism based on dynamic coupling due to buttressing. The non-linearity is expressed as a polynomial with quadratic and cubic terms. This model does not require the use of brake negative coefficient, in order to predict the instability phenomena. Finally, the center manifold approach, the multivariable approximants, and the AFT method are used in order to obtain equations for the limit cycle amplitudes. These methods allow the reduction of the number of equations of the original system in order to obtain a simplified system, without loosing the dynamics of the original system, as well as the contributions of non-linear terms. The goal is the validation of this procedure for a complex non-linear model by comparing results obtained by solving the full system and by using these methods. The brake friction coefficient is used as an unfolding parameter of the fundamental Hopf bifurcation point.     

Theoretical models of vegetable drying by convection

The results of experiments are often used to model empirical phenomena. However, the term model is applied in various meanings. A model is usually treated as an abstract formal structure that can replace a material system considered as original, in respect to the aim of modeling. Certain formal structures may be treated as theoretical models of empirical phenomena. On the other hand, a material system can also be referred to as a model of an abstract system, e.g., a set of equations or a hypothesis. Such a material system, if it is a distinct empirical interpretation of the language of a given theory, is then called a real model. Both kinds of models are applied in drying technology, but the second one is more inventive. The mathematical structures are treated as empirical formulae or as theoretical models properly derived from true or legitimated promises of a given theory. The advantages of some mathematical theoretical models of drying processes versus empirical formulae are discussed. The creation of new mathematical theoretical models of convection drying kinetics of some shrinking solids is presented and analyzed. One of the above models was also hypothetically suggested for modeling the drying of cut vegetables in a fluidized-bed. Despite its initial acceptance due to peer empirical justification on cut carrots and celery, it still requires further theoretical analysis. Other models indicated here are theoretical models of vegetable drying in a tunnel drier. These models are created by deduction from laws of heat and mass transfer theory and its basic equations. XI Polish Drying Symposium, Poznań, Poland, 13–16 September 2005.     

Three-dimensional free vibration analysis of rotating laminated conical shells: layerwise differential quadrature (LW-DQ) method

This paper focuses on the free vibration analysis of thick, rotating laminated composite conical shells with different boundary conditions based on the three-dimensional theory, using the layerwise differential quadrature method (LW-DQM). The equations of motion are derived applying the Hamilton’s principle. In order to accurately account for the thickness effects, the layerwise theory is used to discretize the equations of motion and the related boundary conditions through the thickness of the shells. Then, the equations of motion as well as the boundary condition equations are transformed into a set of algebraic equation applying the DQM in the meridional direction. This study demonstrates the applicability, accuracy, stability and the fast rate of convergence of the present method, for free vibration analyses of rotating thick laminated conical shells. The presented results are compared with those of other shell theories obtained using conventional methods and a special case where the angle of the conical shell approaches zero, that is, a cylindrical shell and excellent agreements are achieved.     

基于进气射流的压气机失速主动控制系统仿真

针对一轴流压气机，构造了利用进气射流对旋转失速进行主动控制的控制系统，发展了以Ｍｏｏｒｅ－Ｇｒｅｉｔｚｅｒ模型为基础的系统数学模型并进行了数值仿真研究。结果表明，该控制系统可有效地抑制压气机中旋转失速的发生。     

Geometry and differentiability requirements in multibody railroad vehicle dynamic formulations

The dynamic equations of multibody railroad vehicle systems can be formulated using different sets of generalized coordinates; examples of these sets of coordinates are the absolute Cartesian and trajectory coordinates. The absolute coordinate based formulations do not require introducing an intermediate track coordinate system since all the absolute coordinates are defined in the global system. On the other hand, when the trajectory coordinates are used, a track coordinate system that follows the motion of a body in the railroad vehicle system is introduced. This track coordinate system is defined by the track geometry and the distance traveled by the body along the track centerline. The configuration of the body with respect to the track coordinate system is defined using five coordinates; two translations and three Euler angles. In this paper, the formulations based on the absolute and trajectory coordinates are compared. It is shown that these two sets of coordinates require different degrees of differentiability and smoothness. When an elastic contact formulation is used to study the wheel/rail dynamic interaction, there are significant differences in the order of the derivatives required in both formulations. In fact, as demonstrated in this study, in the absence of a contact constraint formulation, higher order derivatives with respect to geometric parameters are still required when the equations are formulated using the trajectory coordinates. The formulation of the constraints used in the analysis of the wheel/rail contact is discussed and it is shown that when the absolute coordinates are used, only third order derivatives need to be evaluated. The relationship between the track frame used in railroad vehicle dynamics and the Frenet frame used in the theory of curves to describe the curve geometry is also discussed in this paper. Based on the analysis presented in this paper, the advantages and drawbacks of a hybrid method which employs both the absolute and trajectory coordinates and planar contact conditions in order to reduce the number of contact constraints and relax the differentiability requirements are discussed. In this method, the absolute coordinates are used to formulate the equations of motion of the railroad vehicle system. The absolute coordinate solution can be used to determine the trajectory coordinates and their time derivatives. Using the trajectory coordinates, the motion of the body in the vehicle with respect to the track coordinate system can be predicted and used in the formulation of the planar contact model.     

Multicomponent fluid flow analysis using a new set of conservation equations

In this work hydrodynamics of multicomponent ideal gas mixtures have been studied. Starting from the kinetic equations, the Eulerian approach is used to derive a new set of conservation equations for the multicomponent system where each component may have different velocity and kinetic temperature. The equations are based on the Grad's method of moment derived from the kinetic model in a relaxation time approximation (RTA). Based on this model which contains separate equation sets for each component of the system, a computer code has been developed for numerical computation of compressible flows of binary gas mixture in generalized curvilinear boundary conforming coordinates. Since these equations are similar to the Navier-Stokes equations for the single fluid systems, the same numerical methods are applied to these new equations. The Roe's numerical scheme is used to discretize the convective terms of governing fluid flow equations. The prepared algorithm and the computer code are capable of computing and presenting flow fields of each component of the system separately as well as the average flow field of the multicomponent gas system as a whole. Comparison of the present code results with those of a more common algorithm based on the mixture theory in a supersonic converging-diverging nozzle provides the validation of the present formulation. Afterwards, a more involved nozzle cooling problem with a binary ideal gas (helium-xenon) is chosen to compare the present results with those of the ordinary mixture theory. The present model provides the details of the flow fields of each component separately which is not available otherwise. It is also shown that the separate fluids treatment, such as the present study, is crucial when considering time scales on the order of (or shorter than) the intercollisions relaxation times.     

A two-fluid model for reactive dilute solid–liquid mixtures with phase changes

Based on the Eulerian spatial averaging theory and the Müller–Liu entropy principle, a two-fluid model for reactive dilute solid–liquid mixtures is presented. Initially, some averaging theorems and properties of average quantities are discussed and, then, averaged balance equations including interfacial source terms are postulated. Moreover, constitutive equations are proposed for a reactive dilute solid–liquid mixture, where the formation of the solid phase is due to a precipitation chemical reaction that involves ions dissolved in the liquid phase. To this end, principles of constitutive theory are used to propose linearized constitutive equations that account for diffusion, heat conduction, viscous and drag effects, and interfacial deformations. A particularity of the model is that the mass interfacial source term is regarded as an independent constitutive variable. The obtained results show that the inclusion of the mass interfacial source term into the set of independent constitutive variables permits to easily describe the phase changes associated with precipitation chemical reactions.     

A finite volume solver for 1D shallow‐water equations applied to an actual river

This paper describes the numerical solution of the 1D shallow‐water equations by a finite volume scheme based on the Roe solver. In the first part, the 1D shallow‐water equations are presented. These equations model the free‐surface flows in a river. This set of equations is widely used for applications: dam‐break waves, reservoir emptying, flooding, etc. The main feature of these equations is the presence of a non‐conservative term in the momentum equation in the case of an actual river. In order to apply schemes well adapted to conservative equations, this term is split in two terms: a conservative one which is kept on the left‐hand side of the equation of momentum and the non‐conservative part is introduced as a source term on the right‐hand side. In the second section, we describe the scheme based on a Roe Solver for the homogeneous problem. Next, the numerical treatment of the source term which is the essential point of the numerical modelisation is described. The source term is split in two components: one is upwinded and the other is treated according to a centred discretization. By using this method for the discretization of the source term, one gets the right behaviour for steady flow. Finally, in the last part, the problem of validation is tackled. Most of the numerical tests have been defined for a working group about dam‐break wave simulation. A real dam‐break wave simulation will be shown. Copyright © 2002 John Wiley & Sons, Ltd.