6,478 research outputs found
Geometry Modeling for Unstructured Mesh Adaptation
The quantification and control of discretization error is critical to obtaining reliable simulation results. Adaptive mesh techniques have the potential to automate discretization error control, but have made limited impact on production analysis workflow. Recent progress has matured a number of independent implementations of flow solvers, error estimation methods, and anisotropic mesh adaptation mechanics. However, the poor integration of initial mesh generation and adaptive mesh mechanics to typical sources of geometry has hindered adoption of adaptive mesh techniques, where these geometries are often created in Mechanical Computer- Aided Design (MCAD) systems. The difficulty of this coupling is compounded by two factors: the inherent complexity of the model (e.g., large range of scales, bodies in proximity, details not required for analysis) and unintended geometry construction artifacts (e.g., translation, uneven parameterization, degeneracy, self-intersection, sliver faces, gaps, large tolerances be- tween topological elements, local high curvature to enforce continuity). Manual preparation of geometry is commonly employed to enable fixed-grid and adaptive-grid workflows by reducing the severity and negative impacts of these construction artifacts, but manual process interaction inhibits workflow automation. Techniques to permit the use of complex geometry models and reduce the impact of geometry construction artifacts on unstructured grid workflows are models from the AIAA Sonic Boom and High Lift Prediction are shown to demonstrate the utility of the current approach
Boundary-Conforming Free-Surface Flow Computations: Interface Tracking for Linear, Higher-Order and Isogeometric Finite Elements
The simulation of certain flow problems requires a means for modeling a free
fluid surface; examples being viscoelastic die swell or fluid sloshing in
tanks. In a finite-element context, this type of problem can, among many other
options, be dealt with using an interface-tracking approach with the
Deforming-Spatial-Domain/Stabilized-Space-Time (DSD/SST) formulation. A
difficult issue that is connected with this type of approach is the
determination of a suitable coupling mechanism between the fluid velocity at
the boundary and the displacement of the boundary mesh nodes. In order to avoid
large mesh distortions, one goal is to keep the nodal movements as small as
possible; but of course still compliant with the no-penetration boundary
condition. Standard displacement techniques are full velocity, velocity in a
specific coordinate direction, and velocity in normal direction. In this work,
we investigate how the interface-tracking approach can be combined with
isogeometric analysis for the spatial discretization. If NURBS basis functions
of sufficient order are used for both the geometry and the solution, both a
continuous normal vector as well as the velocity are available on the entire
boundary. This circumstance allows the weak imposition of the no-penetration
boundary condition. We compare this option with an alternative that relies on
strong imposition at discrete points. Furthermore, we examine several coupling
methods between the fluid equations, boundary conditions, and equations for the
adjustment of interior control point positions.Comment: 20 pages, 16 figure
Solid/FEM integration at SNLA
The effort at Sandia National Labs. on the methodologies and techniques being used to generate strict hexahedral finite element meshes from a solid model is described. The functionality of the modeler is used to decompose the solid into a set of nonintersecting meshable finite element primitives. The description of the decomposition is exported, via a Boundary Representative format, to the meshing program which uses the information for complete finite element model specification. Particular features of the program are discussed in some detail along with future plans for development which includes automation of the decomposition using artificial intelligence techniques
Generalizing the advancing front method to composite surfaces in the context of meshing constraints topology
International audienceBeing able to automatically mesh composite geometry is an important issue in the context of CAD-FEA integration. In some specific contexts of this integration, such as using virtual topology or meshing constraints topology (MCT), it is even a key requirement. In this paper, we present a new approach to automatic mesh generation over composite geometry. The proposed mesh generation approach is based on a generalization of the advancing front method (AFM) over curved surfaces. The adaptation of the AFM to composite faces (composed of multiple boundary representation (B-Rep) faces) involves the computation of complex paths along these B-Rep faces, on which progression of the advancing front is based. Each mesh segment or mesh triangle generated through this progression on composite geometry is likely to lie on multiple B-Rep faces and consequently, it is likely to be associated with a composite definition across multiple parametric spaces. Collision tests between new front segments and existing mesh elements also require specific and significant adaptations of the AFM, since a given front segment is also likely to lie on multiple B-Rep faces. This new mesh generation approach is presented in the context of MCT, which requires being able to handle composite geometry along with non-manifold boundary configurations, such as edges and vertices lying in the interior domain of B-Rep faces
Mesh adaptation for high-order flow simulations
Mesh adaptation has only been considered for high-order flow simulations in recent years and many techniques are still to be made more robust and efficient with curvilinear meshes required by these high-order methods. This thesis covers the developments made to improve the mesh generation and adaptation capabilities of the open-source spectral/hp element framework Nektar++ and its dedicated mesh utility NekMesh.
This thesis first covers the generation of quality initial meshes typically required before an iterative adaptation procedure can be used. For optimal performance of the spectral/hp element method, quadrilateral and hexahedral meshes are preferred and two methods are presented to achieve this, either entirely or partially. The first method, inspired from cross field methods, solves a Laplace problem to obtain a guiding field from which a valid two-dimensional quadrilateral block decomposition can be automatically obtained. In turn, naturally curved meshes are generated. The second method takes advantage of the medial axis to generate structured partitions in the boundary layer region of three-dimensional domains. The method proves to be robust in generating hybrid high-order meshes with boundary layer aligned prismatic elements near boundaries and tetrahedral elements elsewhere.
The thesis goes on to explore the adaptation of high-order meshes for the simulation of flows using a spectral/hp element formulation. First a new approach to moving meshes, referred to here as r-adaptation, based on a variational framework, is described. This new r-adaptation module is then enhanced by p-adaptation for the simulation of compressible inviscid flows with shocks. Where the flow is smooth, p-adaptation is used to benefit from the spectral convergence of the spectral/hp element methods. Where the flow is discontinuous, e.g. at shock waves, r-adaptation clusters nodes together to better capture these field discontinuities. The benefits of this dual, rp-adaptation approach are demonstrated through two-dimensional benchmark test cases.Open Acces
Isogeometric Analysis and Harmonic Stator-Rotor Coupling for Simulating Electric Machines
This work proposes Isogeometric Analysis as an alternative to classical
finite elements for simulating electric machines. Through the spline-based
Isogeometric discretization it is possible to parametrize the circular arcs
exactly, thereby avoiding any geometrical error in the representation of the
air gap where a high accuracy is mandatory. To increase the generality of the
method, and to allow rotation, the rotor and the stator computational domains
are constructed independently as multipatch entities. The two subdomains are
then coupled using harmonic basis functions at the interface which gives rise
to a saddle-point problem. The properties of Isogeometric Analysis combined
with harmonic stator-rotor coupling are presented. The results and performance
of the new approach are compared to the ones for a classical finite element
method using a permanent magnet synchronous machine as an example
Exploiting the Design Freedom of RM
This paper details how Rapid Manufacturing (RM) can overcome the restrictions imposed by the
inherent process limitations of conventional manufacturing techniques and become the enabling
technology in fabricating optimal products. A new design methodology capable of exploiting
RM’s increased design freedom is therefore needed. Inspired by natural world structures of trees
and bones, a multi-objective, genetic algorithm based topology optimisation approach is
presented. This combines multiple unit cell structures and varying volume fractions to create a
heterogeneous part structure which exhibits a uniform stress distribution.Mechanical Engineerin
Merging enriched Finite Element triangle meshes for fast prototyping of alternate solutions in the context of industrial maintenance
A new approach to the merging of Finite Element (FE) triangle meshes is proposed. Not only it takes into account the geometric aspects, but it also considers the way the semantic information possibly associated to the groups of entities (nodes, faces) can be maintained. Such high level modification capabilities are of major importance in all the engineering activities requiring fast modifications of meshes without going back to the CAD model. This is especially true in the context of industrial maintenance where the engineers often have to solve critical problems in very short time. Indeed, in this case, the product is already designed, the CAD models are not necessarily available and the FE models might be tuned. Thus, the product behaviour has to be studied and improved during its exploitation while prototyping directly several alternate solutions. Such a framework also finds interest in the preliminary design phases where alternative solutions have to be simulated. The algorithm first removes the intersecting faces in an n-ring neighbourhood so that the filling of the created holes produces triangles whose sizes smoothly evolve according to the possibly heterogeneous sizes of the surrounding triagles. The holefilling algorithm is driven by an aspect ratio factor which ensures that the produced triangulation fits well the FE requirements. It is also constrained by the boundaries of the groups of entities gathering together the simulation semantic. The filled areas are then deformed to blend smoothly with the surroundings meshes
- …