Search results

The purpose of this paper is to apply the variational multi-scale framework to the finite element approximation of the compressible Navier-Stokes equations written in conservation form. Even though this formulation is relatively well known, some particular features that have been applied with great success in other flow problems are incorporated.

The orthogonal subgrid scales, the non-linear tracking of these subscales, and their time evolution are applied. Moreover, a systematic way to design the matrix of algorithmic parameters from the perspective of a Fourier analysis is given, and the adjoint of the non-linear operator including the volumetric part of the convective term is defined. Because the subgrid stabilization method works in the streamline direction, an anisotropic shock capturing method that keeps the diffusion unaltered in the direction of the streamlines, but modifies the crosswind diffusion is implemented. The artificial shock capturing diffusivity is calculated by using the orthogonal projection onto the finite element space of the gradient of the solution, instead of the common residual definition. Temporal derivatives are integrated in an explicit fashion.

Subsonic and supersonic numerical experiments show that including the orthogonal, dynamic, and the non-linear subscales improve the accuracy of the compressible formulation. The non-linearity introduced by the anisotropic shock capturing method has less effect in the convergence behavior to the steady state.

A complete investigation of the stabilized formulation of the compressible problem is addressed.

Error indicators are introduced as part of a simple computational procedure for improving the accuracy of the finite element solutions for plate and shell problems. The procedure is based on using an initial (coarse) grid and a refined (enriched) grid, and approximating the solution for the refined grid by a linear combination of a few global approximation vectors (or modes) which are generated by solving two uncoupled sets of equations in the coarse grid unknowns and the additional degrees of freedom of the refined grid. The global approximation vectors serve as error indicators since they provide quantitative pointwise information about the sensitivity of the different response quantities to the approximation used. The three key elements of the computational procedure are: (a) use of mixed finite element models with discontinuous stress resultants at the element interfaces; (b) operator splitting, or additive decomposition of the finite element arrays for the refined grid into the sum of the coarse grid arrays and correction terms (representing the refined grid contributions); and (c) application of a reduction method through successive use of the finite element method and the classical Bubnov—Galerkin technique. The finite element method is first used to generate a few global approximation vectors (or modes). Then the amplitudes of these modes are computed by using the Bubnov—Galerkin technique. The similarities between the proposed computational procedure and a preconditioned conjugate gradient (PCG) technique are identified and are exploited to generate from the PCG technique pointwise error indicators. The effectiveness of the proposed procedure is demonstrated by means of two numerical examples of an isotropic toroidal shell and a laminated anisotropic cylindrical panel.

In modeling an aircraft wing, structural idealizations are often employed in hand calculations to simplify the structural analysis. In real applications of structural design, analysis and optimization, finite element methods are used because of the complexity of the geometry, combined and complex loading conditions. The purpose of this paper is to give a comprehensive study on the effect of using different structural idealizations on the design, analysis and optimization of thin walled semi-monocoque wing structures in the preliminary design phase.

In the design part of the paper, wing structures are designed by employing two different structural idealizations that are typically used in the preliminary design phase. In the structural analysis part, finite element analysis of one of the designed wing configurations is performed using six different one and two dimensional finite element pairs which are typically used to model the sub-elements of semi-monocoque wing structures. Finally in the optimization part, wing structure is optimized for minimum weight by using finite element models which have the same six different finite element pairs used in the analysis phase.

Based on the results presented in the paper, it is concluded that with the simplified methods, preliminary sizing of the wing configurations can be performed with enough confidence as long as the simplified method based designs are also optimized iteratively, which is what is practiced in the design phase of this study.

This research aims at investigating the effect of using different one and two dimensional element pairs on the final analyzed and optimized configurations of the wing structure, and conclusions are inferred with regard to the sensitivity of the optimized wing configurations with respect to the choice of different element types in the finite element model.

A finite element method for the analysis of combined radiative and conductive heat transport in a finite axisymmetric configuration is presented. The appropriate integro‐differential governing equations for a grey and non‐scattering medium with grey and diffuse walls are developed and solved for several model problems. We consider axisymmetric, cylindrical geometries with top and bottom boundaries of arbitrary convex shape. The method is accurate for media of any optical thickness and is capable of handling a wide array of axisymmetric geometries and boundary conditions. Several techniques are presented to reduce computational overhead, such as employing a Swartz‐Wendroff approximation and cut‐off criteria for evaluating radiation integrals. The method is successfully tested against several cases from the literature and is applied to some additional example problems to demonstrate its versatility. Solution of a free‐boundary, combined‐mode heat transfer problem representing the solidification of a semitransparent material, the Bridgman growth of an yttrium aluminium garnet (YAG) crystal, demonstrates the utility of this method for analysis of a complex materials processing system. The method is suitable for application to other research areas, such as the study of glass processing and the design of combustion furnace systems.

Presents a procedure for obtaining an improved finite element solution of boundary problems by estimating the principle of exact displacement method in the finite element technique. The displacement field is approximated by two types of functions: the shape functions satisfying the homogeneous differential equilibrium equation and the full clamping element functions as a particular solution of the differential equation between the nodes. The full clamping functions represent the solution of the full clamping state on finite elements. An improved numerical solution of displacements, strains, stresses and internal forces, not only at nodes but over the whole finite element, is obtained without an increase of the global basis, because the shape functions are orthogonal with the full clamping functions. This principle is generally applicable to different finite elements. The contribution of introducing two types of functions based on the principle of the exact displacement method is demonstrated in the solution procedure of frame structures and thin plates.

The boundary element method is a useful method for the analysis of field problems involving unbounded regions. Therefore, the method can be used advantageously in combination with the finite element method. This is sometimes called a combination method and it is suitable as a picture‐frame technique. Although this technique attains good accuracy, the matrix of the discretized equation is not banded, since it is a dense matrix. In this paper, we propose an infinite boundary element which divides the unbounded region radially. By the use of this element, the bandwidth of the discretized system matrix does not increase beyond that of the finite element region and its original matrix structure is maintained. The infinite boundary element can also be applied to homogeneous unbounded field problems, for which the Green's function of the mirror image is difficult to use. To illustrate the validity of the proposed technique, some numerical calculations are demonstrated and the results are compared with those of the usual combination method and the method using the hybrid‐type infinite element.

A multi‐variable finite element algorithm based on the generalized Galerkin method is presented. This algorithm takes advantages from the hybrid finite element method and the quasi‐conforming element method. Therefore, it has more flexibility and makes the hybrid and quasi‐conforming concepts easily applied to various physical problems. The formulation of a hybrid quadrilateral finite element with a unified singularity for both thermal and stress analyses has been derived using this algorithm. Some numerical examples are included to illustrate advantages of the hybrid algorithm in both stress and thermal analyses. It is also convenient to employ the algorithm in finite displacement problems.

We investigate the application of a least squares finite element method for the solution of fluid flow problems. The least squares finite element method is based on the minimization of the L2 norm of the equation residuals. Upon discretization, the formulation results in a symmetric, positive definite matrix system which enables efficient iterative solvers to be used. The other motivations behind the development of least squares finite element methods are the applicability of higher order elements and the possibility of using the norm associated to the least squares functional for error estimation. For steady incompressible flows, we develop a method employing linear and quadratic triangular elements and compare their respective accuracy. For steady compressible flows, an implicit conservative least squares scheme which can capture shocks without the addition of artificial viscosity is proposed. A refinement strategy based upon the use of the least squares residuals is developed and several numerical examples are used to illustrate the capabilities of the method when implemented on unstructured triangular meshes.

The paper presents a numerical version of a method which makes it possible to calculate any quasi‐stationary two‐dimensional field in non‐ferromagnetic conductors of limited cross‐sections for unknown boundary conditions on the conductor surfaces. The principle of the method depends on the application of a finite‐element network within the current region and the method of variable division within the no‐current region. The conditions of the continuity of the vector potential and the continuity of the tangential component of the magnetic induction vector are considered in the form of functions defined by the Bubnov‐Galerkin method. Since the separation of the variables is used in the region without current, it is easy to satisfy the conditions at infinity for different types of input. The shapes of the conductors in the current region are, on the other hand, illustrated by discretization by means of the finite element method. An example of the application of the method for two‐dimensional analysis of power losses in two thick solid circular conductors placed in a sinusoidally varying transverse field is presented. Plots of power losses were drawn based on numerical computation.

In finite element analysis of pile driving, the nodes of the finite element mesh are the most important locations for output stresses. Especially at the pile‐soil interface, it is essential to obtain accurate nodal stresses. Several global and local stress smoothing methods available in the literature were reviewed and examined. Global methods are found to be computationally expensive, so results obtained from several local stress smoothing methods are compared. It is shown that accurate nodal stresses can be obtained by approximating the stress distribution inside four‐element patches by a polynomial with order equal to the order of the shape functions. Equally good results can be obtained by approximating the stress distribution inside each element by a bilinear surface. When a method taking into account both equilibrium and boundary conditions was applied, a set of ill‐conditioned matrices was produced for the four‐element patches. Such methods are therefore not recommended.