Search results
1 – 10 of 22Ò. À. Bèg, H.S. Takhar and V.M. Soundalgekar
Numerical results generated by a highly efficient finite‐difference method (originated by Keller for aerodynamical flows at the California Institute of Technology in 1970), and a…
Abstract
Numerical results generated by a highly efficient finite‐difference method (originated by Keller for aerodynamical flows at the California Institute of Technology in 1970), and a robust double shooting Runge‐Kutta‐Merson scheme are presented for the boundary layer equations representing the convection flow of a viscous incompressible fluid past a hot vertical flat plate embedded in a non‐Darcy porous medium. Viscous dissipation due to mechanical work is included in the temperature field equation. The computations for both solution techniques are compared at the leading edge (ξ = 0.0) and found to be in excellent agreement. The effects of the viscous heating parameter (Ec), thermal conductivity ratio (λ) and a Darcy porous parameter (Re/GrDa) on the fluid velocities, temperatures, local shear stress and wall heat transfer rate are discussed with applications to geothermal and industrial flows.
Details
Keywords
This article describes a method for mesh generation, suitable for applications of the finite‐element method, which proceeds fully automatically from a geometric model of the…
Abstract
This article describes a method for mesh generation, suitable for applications of the finite‐element method, which proceeds fully automatically from a geometric model of the object provided by a CAD‐system. It first generates a coarse mesh which is then adapted to fit the finite‐element problem. A resulting system of equations can be solved by a Gaussian‐type matrix method with as few computations as are necessary for a well‐banded matrix, but without the need for node or element numbering.
B.M. Nicolaï and J. De Baerdemaeker
Derives a first order perturbation algorithm for the computation of mean values and (co‐) variances of the transient temperature field in conduction heated materials with random…
Abstract
Derives a first order perturbation algorithm for the computation of mean values and (co‐) variances of the transient temperature field in conduction heated materials with random field parameters. Considers both linear as well as non‐linear heat conduction problems. The algorithm is advantageous in terms of computer time compared to the Monte Carlo method. The computer time can further be reduced by appropriate transformation of the random vectors resulting from the discretization of the random fields.
Details
Keywords
Charbel Farhat and Edward Wilson
Computational algorithms for finite element dynamic analysis of large‐scale structural problems that exploit both concurrent and parallel features of multiple instruction multiple…
Abstract
Computational algorithms for finite element dynamic analysis of large‐scale structural problems that exploit both concurrent and parallel features of multiple instruction multiple data streams computers are presented. A new computer program architecture is used in which large finite element domains are automatically divided into subdomains. The number of subdomains generated is equal to the number of available processors. The spatial solution is obtained using a basis of orthogonal vectors. The temporal solution is computed exactly. Discussion is focused on the concurrent generation of global Ritz vectors. Examples run on a hypercube multiprocessor confirm the potential of the proposed scheme.
C.A.M. Oliveira and R.A.F. Martins
The Semiloof shell element stiffness and mass matrices are analysed. Various integration rules for the stiffness matrix are used, and the influence of these rules on the existence…
Abstract
The Semiloof shell element stiffness and mass matrices are analysed. Various integration rules for the stiffness matrix are used, and the influence of these rules on the existence of mechanisms and on the element spectra is studied. Some methods for lumping the mass matrix are attempted with special reference to a method imposing a given behaviour of the spectra of eigenvalues.
Eligiusz W. Postek, Roland W. Lewis and David T. Gethin
This paper sets out to present developments of a numerical model of squeeze casting process.
Abstract
Purpose
This paper sets out to present developments of a numerical model of squeeze casting process.
Design/methodology/approach
The entire process is modelled using the finite element method. The mould filling, associated thermal and thermomechanical equations are discretized using the Galerkin method. The front in the filling analysis is followed using volume of fluid method and the advection equation is discretized using the Taylor Galerkin method. The coupling between mould filling and the thermal problem is achieved by solving the thermal equation explicitly at the end of each time step of the Navier Stokes and advection equations, which allows one to consider the actual position of the front of the filling material. The thermomechanical problem is defined as elasto‐visco‐plastic described in a Lagrangian frame and is solved in the staggered mode. A parallel version of the thermomechanical program is presented. A microstructural solidification model is applied.
Findings
During mould filling a quasi‐static Arbitrary Lagrangian Eulerian (ALE) is applied and the resulting temperatures distribution is used as the initial condition for the cooling phase. During mould filling the applied pressure can be used as a control for steering the distribution of the solidified fractions.
Practical implications
The presented model can be used in engineering practice. The industrial examples are shown.
Originality/value
The quasi‐static ALE approach was found to be applicable to model the industrial SQC processes. It was found that the staggered scheme of the solution of the thermomechanical problem could parallelize using a multifrontal parallel solver.
Details
Keywords
Youcef Boutora and Noureddine Takorabet
This paper aims to propose a novel direct method for indefinite algebraic linear systems. It is well adapted for sparse linear systems, such as those of two-dimensional (2-D…
Abstract
Purpose
This paper aims to propose a novel direct method for indefinite algebraic linear systems. It is well adapted for sparse linear systems, such as those of two-dimensional (2-D) finite elements problems, especially for coupled systems.
Design/methodology/approach
The proposed method is developed on an example of an indefinite symmetric matrix. The algorithm of the method is given next, and a comparison between the numbers of operations required by the method and the Cholesky method is also given. Finally, an application on a magnetostatic problem for classical methods (Gauss and Cholesky) shows the relative efficiency of the proposed method.
Findings
The proposed method can be used advantageously for 2-D finite elements in stepping methods without using a block decomposition of matrices.
Research limitations/implications
This method is advantageous for direct linear solving for 2-D problems, but it is not recommended at this time for three-dimensional problems.
Originality/value
The proposed method is the first direct solver for algebraic linear systems proposed since more than a half century. It is not limited for symmetric positive systems such as many of direct and iterative methods.
Details
Keywords
Kaveh Salmalian, Ali Alijani and Habib Ramezannejad Azarboni
The purpose of this study is to investigate the post-buckling analysis of functionally graded columns by using three analytical, approximate and numerical methods. A pre-defined…
Abstract
Purpose
The purpose of this study is to investigate the post-buckling analysis of functionally graded columns by using three analytical, approximate and numerical methods. A pre-defined function as an initial assumption for the post-buckling path is introduced to solve the differential equation. The finite difference method is used to approximate the lateral deflection of the column based on the differential equation. Moreover, the finite element method is used to derive the tangent stiffness matrix of the column.
Design/methodology/approach
The non-linear buckling analysis of functionally graded materials is carried out by using three analytical, finite difference and finite element methods. The elastic deformation and Euler-Bernoulli beam theory are considered to establish the constitutive and kinematics relations, respectively. The governing differential equation of the post-buckling problem is derived through the energy method and the calculus variation.
Findings
An incremental iterative solution and the perturbation of the displacement vector at the critical buckling point are performed to determine the post-buckling path. The convergence of the finite element results and the effects of geometric and material characteristics on the post-buckling path are investigated.
Originality/value
The key point of the research is to compare three methods and to detect error sources by considering the derivation process of relations. This comparison shows that a non-incremental solution in the analytical and finite difference methods and an initial assumption in the analytical method lead to an error in results. However, the post-buckling path in the finite element method is traced by the updated tangent stiffness matrix in each load step without any initial limitation.
Details
Keywords
W. Wagner and P. Wriggers
The practical behaviour of problems exhibiting bifurcation with secondary branches cannot be studied in general by using standard path‐following methods such as arc‐length…
Abstract
The practical behaviour of problems exhibiting bifurcation with secondary branches cannot be studied in general by using standard path‐following methods such as arc‐length schemes. Special algorithms have to be employed for the detection of bifurcation and limit points and furthermore for branch‐switching. Simple methods for this purpose are given by inspection of the determinant of the tangent stiffness matrix or the calculation of the current stiffness parameter. Near stability points, the associated eigenvalue problem has to be solved in order to calculate the number of existing branches. The associated eigenvectors are used for a perturbation of the solution at bifurcation points. This perturbation is performed by adding the scaled eigenvector to the deformed configuration in an appropriate way. Several examples of beam and shell problems show the performance of the method.
This paper offers various state‐space representations in the context of applications of the system control theory to dynamic economic systems and examines interrelationships…
Abstract
This paper offers various state‐space representations in the context of applications of the system control theory to dynamic economic systems and examines interrelationships between the alternative representations in both economics literature and system control engineering literature. In particular, some characteristics of various state‐space forms are assessed with respect to the structural properties of each form, thereby demonstrating the relative advantages and disadvantages of different realization methods presented in this paper.