Search results

The purpose of this paper is to discuss about evaluating the integrals involving B-spline wavelet on the interval (BSWI), in wavelet finite element formulations, using Gauss Quadrature.

In the proposed scheme, background cells are placed over each BSWI element and Gauss quadrature rule is defined for each of these cells. The nodal discretization used for BSWI WFEM element is independent to the selection of number of background cells used for the integration process. During the analysis, background cells of various lengths are used for evaluating the integrals for various combination of order and resolution of BSWI scaling functions. Numerical examples based on one-dimensional (1D) and two-dimensional (2D) plane elasto-statics are solved. Problems on beams based on Euler Bernoulli and Timoshenko beam theory under different boundary conditions are also examined. The condition number and sparseness of the formulated stiffness matrices are analyzed.

It is found that to form a well-conditioned stiffness matrix, the support domain of every wavelet scaling function should possess sufficient number of integration points. The results are analyzed and validated against the existing analytical solutions. Numerical examples demonstrate that the accuracy of displacements and stresses is dependent on the size of the background cell and number of Gauss points considered per background cell during the analysis.

The current paper gives the details on implementation of Gauss Quadrature scheme, using a background cell-based approach, for evaluating the integrals involved in BSWI-based wavelet finite element method, which is missing in the existing literature.

Unscented transformation (UT) and point estimate method (PEM) are two efficient algorithms for probabilistic power flow (PPF) computation. This paper aims to show the relevance between UT and PEM and to derive a rule to determine the accuracy controlling parameters for UT method.

The authors derive the underlying sampling strategies of UT and PEM and check them in different probability spaces, where quadrature nodes are selected.

Gauss-type quadrature rule can be used to determine the accuracy controlling parameters of UT. If UT method and PEM select quadrature nodes in two probability spaces related by a linear transform, these two algorithms are equivalent.

It shows that UT method can be conveniently extended to (km + 1) scheme (k = 4; 6; : : : ) by Gauss-type quadrature rule.

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.

The purpose of this paper is to find an accurate, efficient and easy-to-implement point estimate method (PEM) for the statistical moments of random systems.

First, by the theoretical and numerical analysis, the approximate reference variables for the frequently used nine types of random variables are obtained; then by combining with the dimension-reduction method (DRM), a new method which consists of four sub-methods is proposed; and finally, several examples are investigated to verify the characteristics of the proposed method.

Two types of reference variables for the frequently used nine types of variables are proposed, and four sub-methods for estimating the moments of responses are presented by combining with the univariate and bivariate DRM.

In this paper, the number of nodes of one-dimensional integrals is determined subjectively and empirically; therefore, determining the number of nodes rationally is still a challenge.

Through the linear transformation, the optimal reference variables of random variables are presented, and a PEM based on the linear transformation is proposed which is efficient and easy to implement. By the numerical method, the quasi-optimal reference variables are given, which is the basis of the proposed PEM based on the quasi-optimal reference variables, together with high efficiency and ease of implementation.

A few earlier studies presented infeasible heatline trajectories for natural convection within annular domains involving an inner circular cylinder and outer square/circular enclosure. The purpose of this paper is to revisit and illustrate the correct heatline trajectories for various test cases.

Galerkin finite element based methodology and space adaptive grid have been used to simulate natural convective flows within the annular domains. The prediction of heatlines involves derivatives at the nodes, which are evaluated based on finite element basis functions and contributions from neighboring elements.

The heatlines in the earlier work indicate infeasible heat flow paths such as heat flow from one portion to the other of isothermal hot walls and heat flow across the adiabatic walls. Current results illustrate physically consistent heat flow paths involving perpendicularly emerging heatlines from hot to cold walls for conductive transport, long heat flow paths around the closed-loop heatline cells for convective transport and parallel layout of heatlines to the adiabatic walls. Results also demonstrate complex heatlines involving multiple flow vortices and complex flow structures.

Current work translates heatfunctions from energy flux vectors, which are determined by using basis sets. This work demonstrates the expected heatline trajectories for various scenarios involving conductive and convective heat transport within enclosures with an inner hot object as a first attempt, and the results are precursors for the understanding of energy flow estimates.

The boundary integral method (BIM) is very attractive to practicing engineers as it reduces the dimensionality of the problem by one, thereby making the procedure computationally inexpensive compared to its peers. The principal feature of this technique is the limitation of all its computations to only the boundaries of the domain. Although the procedure is well developed for the Laplace equation, the Poisson equation offers some computational challenges. Nevertheless, the literature provides a couple of solution methods. This paper revisits an alternate approach that has not gained much traction within the community. The purpose of this paper is to address the main bottleneck of that approach in an effort to popularize it and critically evaluate the errors introduced into the solution by that method.

The primary intent in the paper is to work on the particular solution of the Poisson equation by representing the source term through a Fourier series. The evaluation of the Fourier coefficients requires a rectangular domain even though the original domain can be of any arbitrary shape. The boundary conditions for the homogeneous solution gets modified by the projection of the particular solution on the original boundaries. The paper also develops a new Gauss quadrature procedure to compute the integrals appearing in the Fourier coefficients in case they cannot be analytically evaluated.

The current endeavor has developed two different representations of the source terms. A comprehensive set of benchmark exercises has successfully demonstrated the effectiveness of both the methods, especially the second one. A subsequent detailed analysis has identified the errors emanating from an inadequate number of boundary nodes and Fourier modes, a high difference in sizes between the particular solution and the original domains and the used Gauss quadrature integration procedures. Adequate mitigation procedures were successful in suppressing each of the above errors and in improving the solution accuracy to any desired level. A comparative study with the finite difference method revealed that the BIM was as accurate as the FDM but was computationally more efficient for problems of real-life scale. A later exercise minutely analyzed the heat transfer physics for a fin after validating the simulation results with the analytical solution that was separately derived. The final set of simulations demonstrated the applicability of the method to complicated geometries.

First, the newly developed Gauss quadrature integration procedure can efficiently compute the integrals during evaluation of the Fourier coefficients; the current literature lacks such a tool, thereby deterring researchers to adopt this category of methods. Second, to the best of the author’s knowledge, such a comprehensive error analysis of the solution method within the BIM framework for the Poisson equation does not currently exist in the literature. This particular exercise should go a long way in increasing the confidence of the research community to venture into this category of methods for the solution of the Poisson equation.

This paper aims to focus on mathematical model development issues, necessary for a better prediction of dynamic responses of articulated rotor helicopters.

The methodology is laid out based on model development for an articulated main rotor, using the theories of aeroelastisity, finite element and state‐space represented indicial‐based unsteady aerodynamics. The model is represented by a set of nonlinear partial differential equations for the main rotor within a state‐space representation for all other parts of helicopter dynamics. The coupled rotor and fuselage formulation enforces the use of numerical solution techniques for trim and linearization calculations. The mathematical model validation is carried out by comparing model responses against flight test data for a known configuration.

Improvements in dynamic prediction of both on‐axis and cross‐coupled responses of helicopter to pilot inputs are observed.

Further work is required for investigation of the unsteady aerodynamics, a state‐space representation, within various compatible dynamic inflow models to describe the helicopter response characteristics.

The results of this work support ongoing research on the development of highly accurate helicopter flight dynamic mathematical models. These models are used as engineering tools both for designing new aerial products such as modernized agile helicopters and optimization of the old version products at minimum time and expense.

Provides further information on the mathematical model development problems associated with advanced helicopter flight dynamics research.

The purpose of the paper is to study natural convection within porous square and triangular geometries (design 1: regular isosceles triangle, design 2: inverted isosceles triangle) subjected to discrete heating with various locations of double heaters along the vertical (square) or inclined (triangular) arms.

Galerkin finite element method is used to solve the governing equations for a wide range of modified Darcy number, Da_m = 10⁻⁵–10⁻² with various fluid saturated porous media, Pr_m = 0.015 and 7.2 at a modified Rayleigh number, Ra_m = 10⁶ involving the strategic placement of double heaters along the vertical or inclined arms (types 1-3). Adaptive mesh refinement is implemented based on the lengths of discrete heaters. Finite element based heat flow visualization via heatlines has been adopted to study heat distribution at various portions.

The strategic positioning of the double heaters (types 1-3) and the convective heatline vortices depict significant overall temperature elevation at both Da_m = 10⁻⁴ and 10⁻² compared to type 0 (single heater at each vertical or inclined arm). Types 2 and 3 are found to promote higher temperature uniformity and greater overall temperature elevation at Da_m = 10⁻². Overall, the triangular design 2 geometry is also found to be optimal in achieving greater temperature elevation for the porous media saturated with various fluids (Pr_m).

Multiple heaters (at each side [left or right] wall) result in enhanced temperature elevation compared to the single heater (at each side [left or right] wall). The results of the current work may be useful for the material processing, thermal storage and solar heating applications.

The heatline approach is used to visualize the heat flow involving double heaters along the side (left or right) arms (square and triangular geometries) during natural convection involving porous media. The heatlines depict the trajectories of heat flow that are essential for thermal management involving larger thermal elevation. The mixing cup or bulk average temperature values are obtained for all types of heating (types 0-3) involving all geometries, and overall temperature elevation is examined based on higher mixing cup temperature values.

The purpose of this paper is to evaluate some numerical integration strategies used in generalized (G)/extended finite element method (XFEM) to solve linear elastic fracture mechanics problems. A range of parameters are here analyzed, evidencing how the numerical integration error and the computational efficiency are improved when particularities from these examples are properly considered.

Numerical integration strategies were implemented in an existing computational environment that provides a finite element method and G/XFEM tools. The main parameters of the analysis are considered and the performance using such strategies is compared with standard integration results.

Known numerical integration strategies suitable for fracture mechanics analysis are studied and implemented. Results from different crack configurations are presented and discussed, highlighting the necessity of alternative integration techniques for problems with singularities and/or discontinuities.

This study presents a variety of fracture mechanics examples solved by G/XFEM in which the use of standard numerical integration with Gauss quadratures results in loss of precision. It is discussed the behaviour of subdivision of elements and mapping of integration points strategies for a range of meshes and cracks geometries, also featuring distorted elements and how they affect strain energy and stress intensity factors evaluation for both strategies.

For those who are working with the finite element method, one may somehow wonder why a n point Gauss‐Legendre quadrature formula can integrate a 2n—1 order polynomial exactly. The purpose of this note is to unravel the secret.