A multithreaded parallel upwind sweep algorithm for the S_(N) transport equations discretized with discontinuous finite elements

The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can overcome the deficiencies of conventionally structured meshes in complex geometry modeling. A multithreaded parallel upwind sweep algorithm for S_(N) transport was proposed to achieve a more accurate geometric description and improve the computational efficiency. The spatial variables were discretized using the standard discontinuous Galerkin finite-element method. The angular flux transmission between neighboring meshes was handled using an upwind scheme. In addition, a combination of a mesh transport sweep and angular iterations was realized using a multithreaded parallel technique. The algorithm was implemented in the 2D/3D S_(N) transport code ThorSNIPE, and numerical evaluations were conducted using three typical benchmark problems:IAEA, Kobayashi-3i, and VENUS-3. These numerical results indicate that the multithreaded parallel upwind sweep algorithm can achieve high computational efficiency. ThorSNIPE, with a multithreaded parallel upwind sweep algorithm, has good reliability, stability, and high efficiency, making it suitable for complex shielding calculations.

Robust Falk-Neilan Finite Elements for the Reissner-Mindlin Plate

The family of Falk-Neilan P_(k)finite elements,combined with the Argyris P_(k+1)finite elements,solves the Reissner-Mindlin plate equation quasi-optimally and locking-free,on triangular meshes.The method is truly conforming or consistent in the sense that no projection/reduction is introduced.Theoretical proof and numerical confirmation are presented.

Use of the finite elements method for modelling the dynamic load impact on hydraulic leg equipped with gas accumulator

Procedures of preparation of numerical analysis,consisting in a simulation of cooperation of three different media: steel,liquid and gas undergoes dynamic load were discussed.Modelling of the initial static load of the mechanical system was presented.By using the MSC.Software products the following exemplary computer simulations were made: dynamic load impact on the hydraulic leg as well as effectiveness of the hydraulic leg protection against overload with help of gas accumulator.

FINITE ELEMENTS FOR THE LUBRICATED CONTACT BETWEEN SOLIDS IN METAL FORMING REOCESSES

In this paper, the lubrication problem in nmmerical sindation of forming processes is considered.After enumerationg the difficulties encountered when trying to solve such a problem with the finite element method, a generalization of the formaulation of Liu[4-6] for the thin film hydrodynamic lubrication re- gime is presented. This method is then aplied to a strip rolling simulation,using the Arbitrary La- grangian eulerian (ALE) formalism.

Multilayers Plates Buckling by Mixed Finite Elements

It is presented a buckling analysis ofmultilayered plates by the finite elements method with mixed unknowns （displacement and rotations, and transverse interlaminar stresses）. In the mixed model, each layer is analyzed as a single plate, where the continuity of displacement is ensured by Lagrange multipliers which represent static variables. This procedure is easy to be implemented from the computational point of view. A methodology to solve the eigen problem is presented based on the inverse iteration method. The model is verified successfully with results obtained by other authors.

Elastodynamic Infinite Elements with Modified Bessel Shape Functions for Dynamic Soil-Structure Interaction

The paper is devoted to formulations of decay and mapped elastodynamic infinite elements, based on modified Bessel shape functions. These elements are appropriate for Soil-Structure Interaction （SSI） problems, solved in time or frequency domain and can be treated as a new form of the recently proposed Elastodynamic Infinite Elements with United Shape Functions （EIEUSF） infinite elements. The formulation of 2D Horizontal type Infinite Elements （HIE） is demonstrated here, but by similar techniques 2D Vertical （VIE） and 2D Comer （CIE） Infinite Elements can also be formulated. Using elastodynamic infinite elements is the easier and appropriate way to achieve an adequate simulation including basic aspects of Soil-Structure Interaction. Continuity along the artificial boundary （the line between finite and infinite elements） is discussed as well and the application of the proposed elastodynamic infinite elements in the Finite Element Method （FEM） is explained in brief. Finally, a numerical example shows the computational efficiency of the proposed infinite elements.

Application of Finite Elements Method for Structural Analysis in a Coffee Harvester

Stress concentration and large displacements are usual problems in the components of the structure of agricultural machinery such harvesters coffee, and that finite element method (FEM) can be a tool to minimize its effects. The goal of this paper is to get results of stresses and displacements of a coffee harvester structure by using FEM for static simulation. The main parts of the coffee harvester analyzed were: engine frame, body right and left sides, front and rear end, main beam, coffee reservoir, wheels and fuel tank. Two different design concepts of a coffee harvester machine were analyzed (structure with rear wheels aligned and misaligned) and the results were compared. It was observed that the model with rear wheels misaligned showed maximum displacement lower than the model with rear wheels aligned. Although higher stress was found in the rear wheels misaligned, it was observed that average stresses for the misaligned wheels design were lower in most structural components analyzed. Based on FEM results, the coffee harvester machine with misaligned rear wheels was built and subjected to operational tests without showing any structural failure.

Computing of the Anchor by the Method of Three-Dimension Point-Radiate Infinite Elements

On the basis of the one-dimension infinite element theory, the coordinate translation and shape function of 3D point-radiate 8-node and 4-node infinite elements are derived. They are coupled with 20-node and 8-node finite elements to compute the compression distortion of the prestressed anchorage segment. The results indicate that when the prestressed force acts on the anchorage head and segment, the stresses and the displacements in the rock around the anchorage head and segment concentrate on the zone center with the anchor axis, and they decrease with exponential forms. Therefore,the stresses and the displacement spindles are formed. The calculating results of the infinite element are close to the theoretical results. This indicates the method is right. This article introduces a new way to study the mechanism of prestressed anchors. The obtained results have an important role in the research of the anchor mechanism and engineering application.

Unified analysis for stabilized methods of low-order mixed finite elements for stationary Navier-Stokes equations

A unified analysis is presented for the stabilized methods including the pres- sure projection method and the pressure gradient local projection method of conforming and nonconforming low-order mixed finite elements for the stationary Navier-Stokes equa- tions. The existence and uniqueness of the solution and the optimal error estimates are proved.

Analysis of the geometrical dependence of auxetic behavior in reentrant structures by finite elements

Materials with a negative Poisson＇s ratio（PR）are called auxetics;they are characterized by expansion/contraction when tensioned/compressed.Given this counterintuitive behavior,they present very particular characteristics and mechanical behavior.Geometrical models have been developed to justify and artificiall reproduce such materials＇ auxetic behavior.The focus of this study is the exploration of a reentrant model by analyzing the variation in the PR of reentrant structures as a function of geometrical and base material parameters.It is shown that,even in the presence of protruding ribs,there may not be auxetic behavior,and this depends on the geometry of each reentrant structure.Values determined for these parameters can be helpful as approximate reference data in the design and fabrication of auxetic lattices using reentrant geometries.

Finite Elements Based on Deslauriers-Dubuc Wavelets for Wave Propagation Problems

This paper presents the formulation of finite elements based on Deslauriers-Dubuc interpolating scaling functions, also known as Interpolets, for their use in wave propagation modeling. Unlike other wavelet families like Daubechies, Interpolets possess rational filter coefficients, are smooth, symmetric and therefore more suitable for use in numerical methods. Expressions for stiffness and mass matrices are developed based on connection coefficients, which are inner products of basis functions and their derivatives. An example in 1-D was formulated using Central Difference and Newmark schemes for time differentiation. Encouraging results were obtained even for large time steps. Results obtained in 2-D are compared with the standard Finite Difference Method for validation.

Ultraconvergence for averaging discontinuous finite elements and its applications in Hamiltonian system

This paper discusses the k-degree averaging discontinuous finite element solution for the initial value problem of ordinary differential equations. When k is even, the averaging numerical flux （the average of left and right limits for the discontinuous finite element at nodes） has the optimal-order ultraconvergence 2k ＋ 2. For nanlinear Hamiltonian systems （e.g., SchrSdinger equation and Kepler system） with momentum conservation, the discontinuous finite element methods preserve momentum at nodes. These properties are confirmed by numerical experiments.

MIXED FINITE ELEMENTS FOR PARABOLICINTEGRO-DIFFERENTIAL EQUATIONS

In this paper, we study mixed finite elements for parabolic integro-differential equations, and introduce a kind of nonclassical mixed projection, its optimal L-2 and h(-s) estimates are obtained. We define semi-discrete and full-discrete mixed finite elements for the equations, and obtain the optimal L-2 error estimates.

Nonconforming finite elements for the equation of planar elasticity

Two new locking-free nonconforming finite elements for the pure displacement planar elasticity problem are presented. Convergence rates of the elements are uniformly optimal with respect to A. The energy norm and L2 norm errors are proved to be O（h2） and O（h3）, respectively. Numerical tests confirm the theoretical analysis.

Analysis of Unbalances Forces Using Methods of Identification and Finite Elements

The study of damage in rotating machineries is of fundamental interest in the fields of machine and structure design. A rotating system, supported by bearings and under some dynamic conditions, can generate a variety of problems that are encountered in many different types of rotating machines. One of these problems is the unbalance due to non-homogeneous mass distribution along the shaft. One of the techniques which are widespread today is the identification of parameters and excitation forces that may well followed by monitoring the evolution and change of possible variations of these parameters. Although several methods for the identification of unbalance excitation force are available in the literature, none of them can be considered unrestricted to be applied for all rotating systems. In this study, two methodologies to identify unknown excitations, such as unbalance, have been proposed. This project refers to the analysis of unbalanced forces from displacement parameters and speed by using methods of identification by Fourier series and Legendre polynomials together with the finite element method, state observers in reasons of the problem of absence of signs of rotational displacement, bandpass filter were used to noise suppression of the data collected from the experimental part, Quasi-Newton method to minimize a function in which the bearing stiffness and its damping are unknowns, and also the experimental verification of the methodology, using for this system owned by a rotary mechanical vibrations of the Department of Mechanical Engineering of Faculty of Engineering, campus of llha Solteira.

Determination of Natural Frequencies of Multilayered Plates by Mixed Finite Elements

Natural frequencies for multilayer plates are calculated by mixed finite element method. The main object of this paper is to use the mixed model for multilayer plates, analyzing each layer as an isolated plate, where the continuity of displacements is achieved by Lagrange multipliers （representing static variables）. This procedure allows us to work with any model for single plate （so as to ensure the proper behavior of each layer）, and the complexity of the multilayer system is avoided by ensuring the condition of displacements by the Lagrange multipliers （static variables）. The plate is discretized by finite element modeling based on a primary hybrid model, where the domain is divided by quadrilateral, both for the displacement field and static variables. This mixed element for plates was implemented and several examples of vibrations have been verified successfully by the results obtained by other methods in the literature.

Development of a Methodology for Determination and Analysis of Thermal Displacements of Machine Tools Using Finite Elements Method and Artificial Neural Network

In the processes of manufacturing, MT （machine tools） plays an important role in the manufacture of work pieces with complex and high dimensional and geometric accuracy. Much of the errors of a machine tool are those which are thermally induced which are from internal and external heat sources acting on the machine. In this paper, a methodology for determining and analyzing the thermal deformation of machine tools using FEM （finite element method） and ANN （artificial neural networks） is presented. After modeling the machine using FEM is defined the location of the heat sources, it is possible to obtain the temperature gradient and the corresponding thermal deformation at predetermined periods. Results obtained with simulations using the software NX.7.5 showed that this methodology is an effective tool in determining the thermal deformation of the machine, correlating the temperature reading at strategic points with volumetric deformation at the tool tip. Therefore, the thermal analysis of the errors in the pair tool part can be established. After training and validation process, the network will be able to make the prediction of thermal errors just stating the temperature values of specific points of each heat source, providing a way for compensation of thermally induced errors.

Contribution to the Full 3D Finite Element Modelling of a Hybrid Stepping Motor with and without Current in the Coils

The paper presents our contribution to the full 3D finite element modelling of a hybrid stepping motor using COMSOL Multiphysics software. This type of four-phase motor has a permanent magnet interposed between the two identical and coaxial half stators. The calculation of the field with or without current in the windings (respectively with or without permanent magnet) is done using a mixed formulation with strong coupling. In addition, the local high saturation of the ferromagnetic material and the radial and axial components of the magnetic flux are taken into account. The results obtained make it possible to clearly observe, as a function of the intensity of the bus current or the remanent induction, the saturation zones, the lines, the orientations and the magnetic flux densities. 3D finite element modelling provide more accurate numerical data on the magnetic field through multiphysics analysis. This analysis considers the actual operating conditions and leads to the design of an optimized machine structure, with or without current in the windings and/or permanent magnet.

A stable implicit nodal integration-based particle finite element method(N-PFEM)for modelling saturated soil dynamics

In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.

Constructing Order Two Superconvergent WG Finite Elements on Rectangular Meshes

In this paper,we introduce a stabilizer free weak Galerkin(SFWG)finite element method for second order elliptic problems on rectangular meshes.With a special weak Gradient space,an order two superconvergence for the SFWG finite element solution is obtained,in both L 2 and H1 norms.A local post-process lifts such a Pk weak Galerkin solution to an optimal order Pk+2 solution.The numerical results confirm the theory.