PML absorbing boundary conditions for the linearized and nonlinear Euler equations in the case of oblique mean flow

For the case of uniform mean flow in an arbitrary direction, perfectly matched layer (PML) absorbing boundary conditions are presented for both the linearized and nonlinear Euler equations. Although linear perfectly matched side layers with an oblique mean flow have been studied in previous works, we propose in the present paper a construction of corner layer equations that are dynamically stable. Stability issues are investigated by examining the dispersion relations of linear waves supported by the corner layer equations. For increased efficiency, a pseudo mean flow is included in the derivation of the PML equations for the nonlinear case. Numerical examples are given to support the validity of the proposed equations. Specifically, the linear PML formulation is tested for the case of acoustic, vorticity, and entropy waves traveling with an oblique mean flow. The nonlinear formulation is tested with an isentropic vortex moving diagonally with a constant velocity. Copyright © 2008 John Wiley & Sons, Ltd.     

SIMILARITY SOLUTIONS OF BOUNDARY LAYER EQUATIONS FOR A SPECIAL NON-NEWTONIAN FLUID IN A SPECIAL COORDINATE SYSTME

Two dimensional equations of steady motion for third order fluids are expressed in a special coordinate system generated by the potential flow corresponding to an inviscid fluid. For the inviscid flow around an arbitrary object, the streamlines are the phicoordinates and velocity potential lines are psi-coordinates which form an orthogonal curvilinear set of coordinates. The outcome, boundary layer equations, is then shown to be independent of the body shape immersed into the flow. As a first approximation, assumption that second grade terms are negligible compared to viscous and third grade terms. Second grade terms spoil scaling transformation which is only transformation leading to similarity solutions for third grade fluid. By ~sing Lie group methods, infinitesimal generators of boundary layer equations are calculated. The equations are transformed into an ordinary differential system. Numerical solutions of outcoming nonlinear differential equations are found by using combination of a Runge-Kutta algorithm and shooting technique.     

Extension of domain‐free discretization method to simulate compressible flows over fixed and moving bodies

This paper is the first endeavour to present the local domain‐free discretization (DFD) method for the solution of compressible Navier–Stokes/Euler equations in conservative form. The discretization strategy of DFD is that for any complex geometry, there is no need to introduce coordinate transformation and the discrete form of governing equations at an interior point may involve some points outside the solution domain. The functional values at the exterior dependent points are updated at each time step to impose the wall boundary condition by the approximate form of solution near the boundary. Some points inside the solution domain are constructed for the approximate form of solution, and the flow variables at constructed points are evaluated by the linear interpolation on triangles. The numerical schemes used in DFD are the finite element Galerkin method for spatial discretization and the dual‐time scheme for temporal discretization. Some numerical results of compressible flows over fixed and moving bodies are presented to validate the local DFD method. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical solution of solidification in a square prism using an algebraic grid generation technique

The solidification of an infinitely long square prism was analyzed numerically. A front fixing technique along with an algebraic grid generation scheme was used, where the finite difference form of the energy equation is solved for the temperature distribution in the solid phase and the solid–liquid interface energy balance is integrated for the new position of the moving solidification front. Results are given for the moving solidification boundary with a circular phase change interface. An algebraic grid generation scheme was developed for two-dimensional domains, which generates grid points separated by equal distances in the physical domain. The current scheme also allows the implementation of a finer grid structure at desired locations in the domain. The method is based on fitting a constant arc length mesh in the two computational directions in the physical domain. The resulting simultaneous, nonlinear algebraic equations for the grid locations are solved using the Newton-Raphson method for a system of equations. The approach is used in a two-dimensional solidification problem, in which the liquid phase is initially at the melting temperature, solved by using a front-fixing approach. The difference of the current study lies in the fact that front fixing is applied to problems, where the solid–liquid interface is curved such that the position of the interface, when expressed in terms of one of the coordinates is a double valued function. This requires a coordinate transformation in both coordinate directions to transform the complex physical solidification domain to a Cartesian, square computational domain. Due to the motion of the solid–liquid interface in time, the computational grid structure is regenerated at every time step.     

比例边界坐标插值方法在谱元法中的应用——无穷域Euler方程的数值模拟

将比例边界坐标插值方法引入谱元法, 构成比例边界谱单元, 对无穷域Euler方程进行数值模拟.阐述了比例边界谱单元的基本使用方法以及基于比例边界谱元的Runge-Kutta间断Galerkin方法求解Euler方程的过程;计算了无穷域圆柱和NACA0012翼型绕流问题, 并与已有结果进行了比较, 显示了计算结果的正确性.用基于比例边界谱元的间断Galerkin方法求解无穷域Euler方程时, 最多只需将求解域划分为2个子域, 避免了一般谱方法将求解域划分为9个或者27个子域的麻烦. 比例边界谱单元为无穷域Euler方程的直接求解提供了一个可供参考的方法.     

Hamiltonian principle based stress singularity analysis near crack corners of multi-material junctions

This paper presents a new method for the stress singularity analysis near the crack corners of a multi-material junctions. The stress singularities near the crack corners of multi-dissimilar isotropic elastic material junctions are studied analytically in terms of the methods developed in Hamiltonian system. The governing equations of plane elasticity in a sectorial domain are derived in Hamiltonian form via variable substitution and variational principle respectively. Both of the methods of global state variable separation and symplectic eigenfunction expansion are used to find the analytical solution of the problem. The relationships among the state vectors in different material spaces are obtained by means of coordinate transformation and consistent conditions between the two adjacent domains. The expression of the original problem is thus changed into a new form where the solutions of symplectic generalized eigenvalues and eigenvectors are needed. The closed form of expressions is established for the stress singularity analysis near the corner with arbitrary vertex angles. Numerical results are presented with several chosen angles and multi-material constants. To show the potential of the new method proposed, a semi-analytical finite element is furthermore developed for the numerical analysis of crack problems.     

A velocity transformation method for the nonlinear dynamic simulation of railroad vehicle systems

The solution of the constrained multibody system equations of motion using the generalized coordinate partitioning method requires the identification of the dependent and independent coordinates. Using this approach, only the independent accelerations are integrated forward in time in order to determine the independent coordinates and velocities. Dependent coordinates are determined by solving the nonlinear constraint equations at the position level. If the constraint equations are highly nonlinear, numerical difficulties can be encountered or more Newton–Raphson iterations may be required in order to achieve convergence for the dependent variables. In this paper, a velocity transformation method is proposed for railroad vehicle systems in order to deal with the nonlinearity of the constraint equations when the vehicles negotiate curved tracks. In this formulation, two different sets of coordinates are simultaneously used. The first set is the absolute Cartesian coordinates which are widely used in general multibody system computer formulations. These coordinates lead to a simple form of the equations of motion which has a sparse matrix structure. The second set is the trajectory coordinates which are widely used in specialized railroad vehicle system formulations. The trajectory coordinates can be used to obtain simple formulations of the specified motion trajectory constraint equations in the case of railroad vehicle systems. While the equations of motion are formulated in terms of the absolute Cartesian coordinates, the trajectory accelerations are the ones which are integrated forward in time. The problems associated with the higher degree of differentiability required when the trajectory coordinates are used are discussed. Numerical examples are presented in order to examine the performance of the hybrid coordinate formulation proposed in this paper in the analysis of multibody railroad vehicle systems.     

On the equations of motion for accelerated frames

In this communication we present the equations of Euler generalized for the motion of a body in an accelerated reference frame using the generalized work-energy principle. The equivalence among the generalized Euler equation, the generalized Lagrange equation, and the generalized Kane equation are shown when applied to the motion of a body of a holonomic system that depend onn generalized coordinates. Therefore when the generalized coordinates can be reduced to two sets of independent coordinates, the generalized Euler equation can be split into two uncoupled equations that are not independent of each other.Universidade da Beira Interior, Covilhã, Portugal. Published in Prikladnaya Mekhanika, Vol. 31, No. 9, pp. 79–89, September, 1995.     

可折叠航天结构展开动力学分析

以欧拉参数为广义坐标（准坐标）,相对角速度和相对移动速度为广义速率,采用Ｋａｎｅ方程的Ｈｕｓｔｏｎ形式建立多体系统的运动力学方程。由伪上三角分解求约束Ｊａｃｏｂｉ矩阵的正交补阵,约简约束力,从而将运动方程由微分几何方程（ＤＡＥ）变为常微分方程（ＯＤＥ）,并由Ｇｅａｒ法对ＯＤＥ积分求出运动历程。最后给出一伸展臂数值分析算例。     

On the numerical solution of tracked vehicle dynamic equations

In this investigation, the solution of the nonlinear dynamic equations of the multibody tracked vehicle systems are obtained using different procedures. In the first technique, which is based on the augmented formulation that employes the absolute Cartesian coordinates and Lagrange multipliers, the generalized coordinate partitioning of the constraint Jacobian matrix is used to determine the independent coordinates and the associated independent differential equations. An iterative Newton-Raphson algorithm is used to solve the nonlinear constraint equations for the dependent variables. The numerical problems encountered when one set of independent coordinates is used during the simulation of large scale tracked vehicle systems are demonstrated and their relationship to the track dynamics is discussed. The second approach employed in this investigation is the velocity transformation technique. One of the versions of this technique is discussed in this paper and the numerical problems that arise from the use of inconsistent system of kinematic equations are reported. In the velocity transformation technique, the tracked vehicle system is assumed to consist of two kinematically decoupled subsystems; the first subsystem consists of the chassis, the rollers, the sprocket and the idler, while the second subsystem consists of the track which is represented as a closed kinematic chain that consists of rigid links connected by revolute joints. It is demonstrated that the use of one set of recursive equations leads to numerical difficulties because of the change in the track configuration. Singular configurations can be avoided by repeated changes in the recursive equations. The sensitivity of the predictor-corrector multistep numerical integration schemes to the method of formulating the state equations is demonstrated. The numerical results presented in this investigation are obtained using a planner tracked vehicle model that consists of fifty four rigid bodies.     

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.     

On general form of Navier-Stokes equations and implicit factored scheme

A general weak conservative form of Navier-Stokes equations expressed with respect to non-orthogonal Curvilinear coordinates and with primitive variables was obtained by using tensor analysis technique, where the contravariant and covariant velocity components were employed. Compared with the current coordinate transformation method, the established equations are concise and forthright, and they are more convenient to be used for solving problems in body-fitted curvilinear coordinate system. An implicit factored scheme for solving the equations is presented with detailed discussions in this paper. For n-dimensional flow the algorithm requires n-steps and for each step only a block tridiagonal matrix equation needs to be solved. It avoids inverting the matrix for large systems of equations and enhances the speed of arithmetic. In this study, the Beam- Warming’s implicit factored schceme is extended and developed in non-orthogonal curvilinear coordinate system.     

Absorbing Boundary Conditions

Absorbing boundary conditions for computational aeroacoustics (CAA) are reviewed. Commonly used absorbing zonal techniques, such as sponge layers and buffer zones, as well as perfectly matched layers (PML) are discussed. The basic ideas and central results of these methods are surveyed and summarized. Special attention will be given to the recently emerged PML technique and its application to CAA. Numerical examples are presented for PML in duct acoustics. A comparison of PML and non-PML absorbing boundary conditions will also be given.     

An immersed boundary method for simulation of inviscid compressible flows

In this paper, an immersed boundary method for simulating inviscid compressible flows governed by Euler equations is presented. All the mesh points are classified as interior computed points, immersed boundary points (interior points closest to the solid boundary), and exterior points that are blanked out of computation. The flow variables at an immersed boundary point are determined via the approximate form of solution in the direction normal to the wall boundary. The normal velocity is evaluated by applying the no‐penetration boundary condition, and therefore, the influence of solid wall in the inviscid flow is taken into account. The pressure is computed with the local simplified momentum equation, and the density and the tangential velocity are evaluated by using the constant‐entropy relation and the constant‐total‐enthalpy relation, respectively. With a local coordinate system, the present method has been extended easily to the three‐dimensional case. The present work is the first endeavor to extend the idea of hybrid Cartesian/immersed boundary approach to compressible inviscid flows. The tedious task of handling multi‐valued points can be eliminated, and the overshoot resulting from the extrapolation for the evaluation of flow variables at exterior points can also be avoided. In order to validate the present method, inviscid compressible flows over fixed and moving bodies have been simulated. All the obtained numerical results show good agreement with available data in the literature. Copyright © 2014 John Wiley & Sons, Ltd.     

Dynamics of a mass–spring–beam with 0:1:1 internal resonance using the analytical and continuation method

The nonlinear harmonic response of an autoparametric system comprised of a linear oscillator with a vertically attached flexural beam is investigated and the capability of the beam as a vibration absorber is assessed. A weak torsional spring is used for constraining the rotation of the beam giving rise to an almost non-flexural rotational mode with a low frequency. The system parameters are also tuned to enforce the zero-to-one-to-one internal resonance condition. The Lagrange’s formulation accompanied by the assumed-mode method is used to derive the discretized equations based on the order three nonlinear Euler–Bernoulli beam theory. An analytical solution is developed based on the method of multiple scales where the generalized coordinate corresponding to the non-flexural rotational mode is approximated by higher order perturbation expansion than the other coordinates, due to much larger contribution of the non-flexural rotation to the response. Comprehensive response and bifurcation analysis are performed using analytical and direct numerical solutions. The results are obtained for vertically-aligned and also initially inclined beams and various interesting behaviors are recognized for different non-dimensional system parameters. Different types of bifurcations such as the Pitch-fork, Hopf, Period-doubling and symmetry breaking bifurcations are observed in the solution of slow-flow equations and some of them are found to be beneficial for vibration absorption in a limited range of excitation amplitudes and frequencies.     

二维稳态辐射声场的光滑有限元-完美匹配层解法

针对外场声学有限元计算精度偏低的问题, 将光滑有限元技术引入到二维稳态辐射声场预测中, 提出了光滑有限元-完美匹配层解法. 该解法采用完美匹配层截断声场计算域, 并将其离散为等参四边形单元, 采用指数吸收函数实现完美匹配层内参数坐标和笛卡尔坐标的映射关系, 采用光滑声压梯度技术计算辐射声场刚度矩阵, 将形函数梯度的域内积分转换为形函数域边界积分. 某汽车二维声腔辐射声场的数值分析结果表明, 与标准有限元-完美匹配层相比, 光滑有限元-完美匹配层解法在完美匹配层内的声波吸收效果更好, 在计算域内的数值计算精度更高, 具有良好的工程应用前景.      

坐标变换法数值求解微通道行波电场电渗流

采用坐标变换法数值求解了耦合的Poisson-Nernst-Planck (PNP)方程和Navier-Stokes(NS)方程, 研究二维狭窄微通道行波电场电渗流数值解. 数值结果表明,坐标变换法能有效降低电渗流解数值解在双电层的高梯度, 有效改善数值解的收敛性和稳定性. 坐标变换的电渗流数值解和原始坐标下的数值解完全一致. 坐标变换后采用简单的网格也能得到和原始坐标下复杂网格相同的解. 给出了滑移边界的近似解与完整的PNP-NS数值解的比较. 在双电层厚度与微通道深度比值(λ/H)很小的情况下(相对深通道), 两者的解基本一致. 但在λ/H较大时(相对浅通道)滑移边界的解高于电渗流速度.