Numerical solution of the turbulent flow of an incompressible fluid in curved planar and axially symmetrical channels

The governing equations for axially symmetric flow, where the Reynolds stresses are expressed by scalar turbulent viscosity, are the Reynolds equations. The turbulence model k, ? is used in the well-known form for fully developed turbulent flow.The numerical method, a continuation of the MAC system¹, is adapted so that even for high Reynolds cell numbers precision (δx²) can be achieved for the steady flow. Irregular cells join the rectangular network on the curved surface. Von Neumann's stability condition of the linearised numerical system is investigated. Special problems concerning the numerical solution of the turbulence model equations are stated and a special procedure is worked out to ensure that the fields k, ? do not converge to physically meaningless values. The program for the computer is universal in that the boundary problems can be assigned by input data.As an example, an axially symmetrical diffuser with an area ratio of widening 1.40 is computed. Fields of velocity and pressure at the wall as well as fields v_T and k are assessed. The results are compared with an experiment. The conclusion is that this method is suitable for the problems mentioned in this study as well as for unsteady flow.     

Method of reverberation ray matrix for static analysis of planar framed structures composed of anisotropic Timoshenko beam members

Based on the method of reverberation ray matrix(MRRM), a reverberation matrix for planar framed structures composed of anisotropic Timoshenko(T) beam members containing completely hinged joints is developed for static analysis of such structures.In the MRRM for dynamic analysis, amplitudes of arriving and departing waves for joints are chosen as unknown quantities. However, for the present case of static analysis, displacements and rotational angles at the ends of each beam member are directly considered as unknown quantities. The expressions for stiffness matrices for anisotropic beam members are developed. A corresponding reverberation matrix is derived analytically for exact and unified determination on the displacements and internal forces at both ends of each member and arbitrary cross sectional locations in the structure. Numerical examples are given and compared with the finite element method(FEM) results to validate the present model. The characteristic parameter analysis is performed to demonstrate accuracy of the present model with the T beam theory in contrast with errors in the usual model based on the Euler-Bernoulli(EB) beam theory. The resulting reverberation matrix can be used for exact calculation of anisotropic framed structures as well as for parameter analysis of geometrical and material properties of the framed structures.     

Homogenization of two-phase elasto-plastic composite materials and structures: Study of tangent operators,cyclic plasticity and numerical algorithms

We develop homogenization schemes and numerical algorithms for two-phase elasto-plastic composite materials and structures. A Hill-type incremental formulation enables the simulation of unloading and cyclic loadings. It also allows to handle any rate-independent model for each phase. We study the crucial issue of tangent operators: elasto-plastic (or “continuum”) versus algorithmic (or “consistent”), and anisotropic versus isotropic. We apply two methods of extraction of isotropic tangent moduli. We compare mathematically the stiffnesses of various tangent operators. All rate equations are discretized in time using implicit integration. We implemented two homogenization schemes: Mori–Tanaka and a double inclusion model, and two plasticity models: classical J₂ plasticity and Chaboche’s model with non-linear kinematic and isotropic hardenings. We consider composites with different properties and present several discriminating numerical simulations. In many cases, the results are validated against finite element (FE) or experimental data. We integrated our homogenization code into the FE program ABAQUS using a user material interface UMAT. A two-scale procedure allows to compute realistic structures made of non-linear composite materials within reasonable CPU time and memory usage; examples are shown.     

Stability and non-linear second-order elastic analyses of beam and framed structures with semi-rigid connections using the cross method

The main objective of this publication is to present an extended version of the Moment Distribution Method (MDM) for the stability and non-linear second-order analysis of indeterminate beams and framed structures made of beam-columns of symmetrical cross-section including the combined effects of shear and bending deformations, axial loads, and semi-rigid connections. The proposed method along each member has the following advantages: (1) it can be utilized in the first- and second-order analyses (including buckling analysis) of indeterminate beams and framed structures made of beam-columns with rigid, semi-rigid, and simple end connections; (2) the effects of semi-rigid connections are condensed into the bending stiffness and fixed-end moments without introducing additional degrees of freedom and equations of equilibrium; and (3) it is accurate, powerful, practical, versatile, and an excellent teaching tool. Analytical studies indicate that shear deformations, semi-rigid connections, and axial loads increase the lateral deflections and affect the internal moments and reactions of continuous beams and framed structures. These effects must be taken into account particularly in slender structures and when they are made of beam or columns with relatively low effective shear areas (like laced columns, columns with batten plates or with perforated cover plates, and columns with open webs) or with low shear stiffness (like short columns made of laminated composites with low shear modulus G when compared to their elastic modulus E) making the shear stiffness GA_s of the same order of magnitude as EI/L². These effects become even more significant when the external supports are not perfectly clamped. Three comprehensive examples are included that show the effectiveness of the proposed method.     

A powerful force-based approach for the limit analysis of three-dimensional frames

For the estimation of the strength of a structure, one could avoid detailed elastoplastic analysis and resort, instead, to direct limit analysis methods that are formulated within linear programming. This work describes the application of the force method to the limit analysis of three-dimensional frames. For the limit analysis of a framed structure, the force method, being an equilibrium-based approach, is better suited than the displacement method and results, generally, to faster solutions. Nevertheless, the latter has been used mostly, since it has a better potential for automation. The difficulty for the direct computerization of the force method is to automatically pick up the structure’s redundant forces. Graph theory concepts may be used to accomplish this task, and a numerical procedure was proposed for the optimal plastic design of plane frames. An analogous approach is developed herein for the limit analysis of space frames which is computationally more cumbersome than the limit analysis of plane frames. The proposed procedure results in hypersparse matrices, and in conjunction with the kinematic upper bound linear program which is solved by a sparse solver, the proposed method appears computationally very efficient. It is also proved that it is much more effective than any displacement-based formulation. The robustness and efficiency of the approach are testified by numerical examples for grillages and multi-storey frames that are included.     

Air flow aerodynamic on a wall-mounted hemisphere for various turbulent boundary layers

Air flow field around a surface-mounted hemisphere of a fixed height for two different turbulent boundary layers (thin and thick) are investigated experimentally and numerically. Flow measurements are performed in a wind tunnel using hot-wire anemometer and streamwise component of velocity fluctuation are calculated using a special developed program of the hardware system. Mean surface pressure coefficients and velocity field for the same hemisphere are determined by the numerical simulation. Turbulent flow field and intensity are measured for two types of boundary layers and compared at various sections in both streamwise and spanwise directions. Numerical scheme based on finite volume and SIMPLE algorithm is used to treat pressure and velocity coupling. Studies are performed for Reynolds number, Re_H = 32,000. Based on the numerical simulation using RNG k–ε turbulence model, flow pathlines, separation region and recirculation area are determined for the two types of turbulent boundary layer flows and complex flow field and recirculation regions are identified and presented graphically.     

An isotropic-plasticity-based constitutive model for martensitic reorientation and shape-memory effect in shape-memory alloys

In this work, we develop an isotropic-plasticity-based constitutive model for initially martensitic shape-memory alloys (SMA) which exhibit martensitic reorientation and the shape-memory effect. The constitutive model is then implemented in the [Abaqus reference manuals. 2006. Providence, R.I.] finite-element program by writing a user-material subroutine. The results from the constitutive model and numerical procedure are then compared to representative physical experiments conducted on polycrystalline rod and sheet Ti–Ni. The constitutive model and the numerical simulations are able to reproduce the stress–strain responses from these physical experiments to good accuracy. Finally, two different boundary value problems utilizing the one-way shape-memory effect are studied: (a) the deformation of an arterial stent, and (b) a micro-clamper. We show that our constitutive model can be used to model the response of the aforementioned boundary value examples.     

Modeling and simulation of an acoustic well stimulation method

This paper presents a mathematical model and a numerical procedure to simulate an acoustic well stimulation (AWS) method for enhancing the permeability of the rock formation surrounding oil and gas wells. The AWS method considered herein aims to exploit the well-known permeability-enhancing effect of mechanical vibrations in acoustically porous materials, by transmitting time-harmonic sound waves from a sound source device—placed inside the well—to the well perforations made into the formation. The efficiency of the AWS is assessed by quantifying the amount of acoustic energy transmitted from the source device to the rock formation in terms of the emission frequency and the well configuration. A simple methodology to find optimal emission frequencies for a given well configuration is presented. The proposed model is based on the Helmholtz equation, a sound-hard boundary condition at the casing, and an impedance boundary condition that effectively accounts for the porous solid–fluid interaction at the interface between the rock formation and the well perforations. Exact non-reflecting boundary conditions derived from Dirichlet-to-Neumann maps are utilized to truncate the circular cylindrical waveguides considered in the model. The resulting boundary value problem is then numerically solved by means of the finite element method. A variety of numerical examples are presented in order to demonstrate the effectiveness of the proposed procedure for finding optimal emission frequencies.     

Investigating effects of various base restraints on the nonlinear inelastic static and seismic responses of steel frames

This paper deals with effects of various base restraints on the nonlinear inelastic static and seismic response of plane and space steel frames. The inelastic behavior is captured by a plastic fiber beam-column method, in which the beam-column member is monitored by integration points along the member length, and the cross-section is meshed into several sub-sections. The second-order effects are considered through the use of stability functions and the geometric stiffness matrix. The effect of shear deformation is also taken into account. The column-base restraint is simulated by using a multi-spring connection element developed by authors. The independent hardening model is employed for performing hysteretic loops of rotational springs under seismic loadings, whereas mathematical models are adopted for representing moment-rotation curves of those springs under static loadings. The accuracy and efficiency of the proposed program are compared with an experimental test, numerical examples of previous studies or SAP2000 commercial software. The results show significant differences of steel frames with various base restraints. Nonlinear semi-rigid base restraints perform strongly damping ability due to energy dissipation through moment-rotation hysteretic loops.     

大坝混凝土三维细观力学数值模型研究

在细观结构层次上将大坝混凝土作为骨料、固化水泥砂浆及其粘结界面组成的复合材料,建立了大坝混凝土三维细观力学数值模型。该模型既能够反映混凝土及其细观各相材料在荷载作用下的损伤演化过程,又考虑了动载作用的应变率强化效应。给出了该数值模型求解方法,并编制出能够在普通PC机上运行的串行程序。加载过程既可按荷载控制又可按位移控制。同时,为了减少求解自由度应用了分尺度方法以使最小骨料和固化水泥砂浆混合后其力学性能与一种复合介质等效。通过混凝土湿筛和三级配试件的静、动(冲击)弯拉数值计算验证了本文计算方法和程序正确有效。另外,在串行程序的基础上,优化了刚度矩阵的存储方式,采用双门槛不完全Cholesky分解(ICT)预条件的共轭梯度法(CG),完成了能够在Sun Fire 6800服务器实现并行计算的并行程序改造,从而大大提高了计算效率。     

Numerical simulation of vortex-induced vibration of a circular cylinder at low mass-damping using RANS code

Fundamental research on vortex-induced vibration (VIV) of a circular cylinder is still needed to build more rational VIV analysis tools for slender marine structures. Numerical results are presented for the response of an elastically mounted rigid cylinder at low mass damping constrained to oscillate transversely to a free stream. A two-dimensional Reynolds-averaged Navier–Stokes (RANS) code equipped with the SST k−ω turbulence model is applied for the numerical calculations. The numerical results are compared in detail with recent experimental and computational work. The Reynolds-averaging procedure erases the random disturbances in the vortex shedding process, so that the comparison between experimental data and the numerical results obtained by RANS codes may reveal some random characteristics of the VIV response. How random disturbance affects the observation in the experiments is discussed in this paper and the issues influencing the appearance of the upper branch in experiments are especially investigated. The absence of the upper branch in RANS simulations is explained in depth on account of discrepancies, which exist between experiments and RANS simulations. In addition, the formation of the 2P vortex shedding mode and its transition through the lock-in region are well reproduced in this investigation.     

Optimal truss design based on an algorithm using optimality criteria

A computational scheme is presented for the calculation of the optimal design of trusses. Constraints on the design variables (the cross-sectional areas) are considered. Linearly elastic behavior is assumed, and optimality criteria are derived, based on strain energy considerations. As in mathematical programming techniques, the optimum is approached through a sequence of designs, each differing slightly from its predecessor. The design changes to be made at each stage of the procedure are determined by application of the optimablity criteria. The formulation is sufficiently general to permit the solution of the problem of predicting both optimal member size and member layout-given the loads and the location of the joints. The procedure is illustrated with a number of numerical examples.     

Designing wing airfoils with active flow controls

A method of designing wing airfoils in separationless flow with suction of a portion of the external flow and reactive jet injection from the rear of the body (the total pressure and the density in the jet are different from those in the freestream) within the framework of the ideal incompressible fluid model is proposed. It is shown that this method of active flow control makes it possible considerably to increase the airfoil lift as compared with the same airfoil with no suction or injection. Examples of the design of such airfoils in separationless flow are presented. The reliability of the results obtained is confirmed by a numerical experiment using the Fluent program.     

A new beam element for analyzing geometrical and physical nonlinearity

Based on Timoshenko＇s beam theory and Vlasov＇s thin-walled member theory, a new model of spatial thin-walled beam element is developed for analyzing geometrical and physical nonlinearity, which incorporates an interior node and independent interpolations of bending angles and warp and takes diversified factors into consideration, such as traverse shear deformation, torsional shear deformation and their coupling, coupling of flexure and torsion, and the second shear stress. The geometrical nonlinear strain is formulated in updated Lagarange （UL） and the corresponding stiffness matrix is derived. The perfectly plastic model is used to account for physical nonlinearity, and the yield rule of von Mises and incremental relationship of Prandtle-Reuss are adopted. Elastoplastic stiffness matrix is obtained by numerical integration based on the finite segment method, and a finite element program is compiled. Numerical examples manifest that the proposed model is accurate and feasible in the analysis of thin-walled structures.     

Determination of the non-linear parameter (mobility factor) of the Giesekus constitutive model using LAOS procedure

The paper introduces a novel procedure to determine the non-linear parameter of the Giesekus model, in relation to the characterization of the non-linear oscillatory shear regime of viscoelastic polymer solutions based on polyacrylamide. Instead of using the shear-thinning viscosity as the representative non-linear effect, the third harmonic in the Fourier spectrum of the shear stress response signal is considered for computing the mobility factor. The fluid is subjected to large amplitude oscillatory shear (LAOS) and its response is recorded. Deviations of this signal from the sinusoidal form are specific to each material and gives both qualitative and quantitative measures of the non-linearity. By fitting the material response with the corresponding numerical solutions of the n-modes Giesekus constitutive relation, one can extract the values of the non-linear α_i-parameters that describe the fluid rheology. It is demonstrated that this procedure, which can be successfully applied to semi-concentrated polymer solutions, provides better results than the classical viscosity-fit method.     

Numerical investigation of the slope discontinuities in large deformation finite element formulations

In many multibody system applications, the system components are made of structural elements that can have different orientations, leading to slope discontinuities. In this paper, a numerical investigation of a new procedure that can be used to model structures with slope discontinuities in the finite element absolute nodal coordinate formulation (ANCF) is presented. This procedure can be applied to model slope discontinuities in the case of commutative rotations of gradient deficient elements that are used for modeling thin beam and plate structures. An important special case to which the proposed procedure can be applied is the case of all planar gradient deficient ANCF finite elements. The use of the proposed method leads to a constant orthogonal element transformation that describes an arbitrary initial configuration. As a consequence, one obtains, in the case of large commutative rotations and large deformations, a constant mass matrix for structures which have complex geometry. The procedure used in this investigation to model slope discontinuities requires the use of the concept of the intermediate finite element coordinate system. For each finite element, a new set of gradient coordinates that define, at the discontinuity node, the element deformation with respect to the intermediate element coordinate system is introduced. These new gradient coordinates are assumed to be equal for the two finite elements at the point of intersection. That is, the change of the gradients of two elements at the intersection point from their respective intermediate initial reference configuration is assumed to be the same. This procedure leads to a set of linear algebraic equations that define the orthogonal transformation matrix for the finite element. Numerical examples are presented in order to demonstrate the use of the proposed procedure for modeling slope discontinuities.     

Improving separated-flow predictions using an anisotropy-capturing subgrid-scale model

The major conclusion of this paper is that resolution requirements for large-eddy simulation (LES) of flow separation and reattachment can be significantly reduced using the anisotropy-capturing explicit algebraic subgrid-scale (SGS) stress model (EASSM) of Marstorp et al. (J. Fluid Mech., vol. 639, 2009, pp. 403–432), instead of the conventional isotropic dynamic eddy-viscosity model (DEVM). LES of flow separation in a channel with streamwise periodic hill-shaped constrictions and spanwise homogeneity is performed at coarse resolutions for which it is observed that flow separation cannot be predicted without a SGS model and cannot be correctly predicted by the DEVM, while reasonable predictions are obtained with the EASSM. It is shown that the lower resolution requirements by the EASSM, compared to the DEVM, is not only due its nonlinear formulation, but also due to the better formulation of its eddy-viscosity part. The improvements obtained with the EASSM have previously been demonstrated using higher-order numerical solvers for channel flows. In this study, it is observed that these improvements still remain using a low-order code with significant inherent numerical dissipation.     

Adaptive and anisotropic mesh strategy for thin shell problems. Case of inhibited parabolic shells

The aim of this paper is to show the reliability of an adaptive and anisotropic mesh procedure for thin shell problems. We consider singular perturbation problems only for parabolic shells whose behavior is described by the Koiter model. The corresponding system of equations, which depends on the relative thickness ε of the shell, is elliptic except at the limit for ε = 0 where it is parabolic. In a first part of this paper, we study theoretically the phenomena of internal layers appearing during the singular perturbation process, when the loading is somewhat singular. These layers have very different structures either they are along or across the asymptotic lines of the middle surface of the shell. In a second part, numerical computations are performed using a finite element software coupled with an adaptive anisotropic mesh generator. This technique enables to approach accurately the singularities and the layers predicted by the theory especially for very small values of the thickness. The efficiency of such a procedure in comparison with uniform meshes is put in a prominent position.     

Numerical Simulation of Non-premixed Turbulent Combustion Using the Eddy Dissipation Concept and Comparing with the Steady Laminar Flamelet Model

Numerical simulations of the Sandia flame CHNa and the Sydney bluff-body stabilized flame HM1E are reported and the results are compared to available experimental data. The numerical method is based on compressible URANS formulations which were implemented recently in the OpenFOAM toolbox. In this study, the calculations are carried out using the conventional compressible URANS approach and a standard k- ?? turbulence model. The Eddy Dissipation Concept with a detailed chemistry approach is used for the turbulence-chemistry interaction. The syngas (CO/H₂) chemistry diluted by 30 % nitrogen in the Sandia flame CHNa and CH₄/H₂ combustion in the Sydney flame HM1E are described by the full GRI-3.0 mechanism. A robust implicit Runge-Kutta method (RADAU5) is used for integrating stiff ordinary differential equations to calculate the reaction rates. The radiation is treated by the P1-approximation model. Both target flames are predicted with the Steady Laminar Flamelet model using the commercial code ANSYS FLUENT as well. In general, there is good agreement between present simulations and measurements for both flames, which indicates that the proposed numerical method is suitable for this type of combustion, provides acceptable accuracy and is ready for further combustion application development.