Search NASA⌕ Search

SEARCH · Search NASA

Results for “Triangular meshes”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

Adaptive remeshing for compressible flow computations

The present, quality-enhancing adaptive-mesh procedure for two-dimensional Euler equation steady state solutions is implemented by means of linear triangular elements and an explicit time-stepping scheme, in conjunction with a finite element solution algorithm. The meshes thus generated typically take the form of stretched elements in the vicinity of one-dimensional flow features; a considerable variation in element size may thereby emerge which allows the desired high-quality solutions to be obtained with commensurately high computational efficiency.

Peraire, J.↗

TokaMaker: An open-source time-dependent Grad-Shafranov tool for the design and modeling of axisymmetric fusion devices

Here, in this paper, we present a new static and time-dependent MagnetoHydroDynamic (MHD) equilibrium code, TokaMaker, for axisymmetric configurations of magnetized plasmas, based on the well-known Grad-Shafranov equation. This code utilizes finite element methods on an unstructured triangular grid to enable capturing accurate machine geometry and simple mesh generation from engineering-like descriptions of present and future devices. The new code is designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. A detailed description of the numerical methods of the code, including a novel formulation of the boundary conditions for free-boundary equilibria, and validation of the implementation of those methods using both analytic test cases and cross-code validation is shown. Results show expected convergence across tested polynomial degree for analytic and cross-code test cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A mixed volume grid approach for the Euler and Navier-Stokes equations

An approach for solving the compressible Euler and Navier-Stokes equations upon meshes composed of nearly arbitrary polyhedra is described. Each polyhedron is constructed from an arbitrary number of triangular and quadrilateral face elements, allowing the unified treatment of tetrahedral, prismatic, pyramidal, and hexahedral cells, as well the general cut cells produced by Cartesian mesh approaches. The basics behind the numerical approach and the resulting data structures are described. The accuracy of the mixed volume grid approach is assessed by performing a grid refinement study upon a series of hexahedral, tetrahedral, prismatic, and Cartesian meshes for an analytic inviscid problem. A series of laminar validation cases are made, comparing the results upon differing grid topologies to each other, to theory, and experimental data. A computation upon a prismatic/tetrahedral mesh is made simulating the laminar flow over a wall/cylinder combination.

Coirier, William J.↗

Adaptive unstructured meshing for thermal stress analysis of built-up structures

An adaptive unstructured meshing technique for mechanical and thermal stress analysis of built-up structures has been developed. A triangular membrane finite element and a new plate bending element are evaluated on a panel with a circular cutout and a frame stiffened panel. The adaptive unstructured meshing technique, without a priori knowledge of the solution to the problem, generates clustered elements only where needed. An improved solution accuracy is obtained at a reduced problem size and analysis computational time as compared to the results produced by the standard finite element procedure.

Dechaumphai, Pramote↗

Adaptive unstructured meshing for thermal stress analysis of built-up structures

An adaptive unstructured meshing technique for mechanical and thermal stress analysis of built-up structures has been developed. A triangular membrane finite element and a new plate bending element are evaluated on a panel with a circular cutout and a frame stiffened panel. The adaptive unstructured meshing technique, without a priori knowledge of the solution to the problem, generates clustered elements only where needed. An improved solution accuracy is obtained at a reduced problem size and analysis computational time as compared to the results produced by the standard finite element procedure.

Dechaumphai, Pramote↗

An adaptive remeshing method for finite element thermal analysis

A finite element remeshing approach that makes use of quadrilateral and triangular elements is described. The approach uses the solution on a previous mesh to create a new mesh. Meshes are completely unstructured with highly refined elements in regions of steep gradients and larger elements where gradients are smaller. Studies of convergence rates for heat conduction problems with exact solutions show that for problems with highly localized solution variations, the remeshing approach gives smaller solution errors with fewer unknowns than refinement of uniform, structured meshes.

Thornton, Earl A.↗

Adaptive remeshing method for finite-element thermal analysis

A finite-element remeshing approach that makes use of quadrilateral and triangular elements is described. The approach uses the solution on a previous mesh to create a new mesh. Meshes are completely unstructured with highly refined elements in regions of steep gradients and larger elements where gradients are smaller. Studies of convergence rates for heat conduction problems with exact solutions show that for problems with highly localized solution variations, the remeshing approach gives smaller solution errors with fewer unknowns than refinement of uniform, structured meshes.

Thornton, Earl A.↗

Accuracy of Gradient Reconstruction on Grids with High Aspect Ratio

Gradient approximation methods commonly used in unstructured-grid finite-volume schemes intended for solutions of high Reynolds number flow equations are studied comprehensively. The accuracy of gradients within cells and within faces is evaluated systematically for both node-centered and cell-centered formulations. Computational and analytical evaluations are made on a series of high-aspect-ratio grids with different primal elements, including quadrilateral, triangular, and mixed element grids, with and without random perturbations to the mesh. Both rectangular and cylindrical geometries are considered; the latter serves to study the effects of geometric curvature. The study shows that the accuracy of gradient reconstruction on high-aspect-ratio grids is determined by a combination of the grid and the solution. The contributors to the error are identified and approaches to reduce errors are given, including the addition of higher-order terms in the direction of larger mesh spacing. A parameter GAMMA characterizing accuracy on curved high-aspect-ratio grids is discussed and an approximate-mapped-least-square method using a commonly-available distance function is presented; the method provides accurate gradient reconstruction on general grids. The study is intended to be a reference guide accompanying the construction of accurate and efficient methods for high Reynolds number applications

Thomas, James↗

NASTRAN finite element idealization study

The investigation of the effects of variations of mesh refinement and mesh pattern were conducted using a basic rectangular mesh pattern. When employing the constant strain TRMEM element, the basic rectangular pattern was subdivided into triangles. This subdivision employs two different triangular patterns and allows results to be obtained which demonstrate the effect of modelling bias. Errors in tip deflection, direct stress, and shearing stress as a function of mesh size and element aspect ratio were obtained as well as mid-span stress distributions. All problems were solved on an IBM 360/95 computer using MacNeal-Schwendler Version MSC-38 Rigid Format-1 of the NASTRAN computer program. While NASTRAN uses double precision arithmetic for the solution of the global equations for displacements, subsequent computations to obtain element stresses are carried out in single precision. This suggests that some improvement in stress recovery might be expected when using the higher precision DCD machines.

Case, W. R.↗

Self-Avoiding Walks Over Adaptive Triangular Grids

Space-filling curves is a popular approach based on a geometric embedding for linearizing computational meshes. We present a new O(n log n) combinatorial algorithm for constructing a self avoiding walk through a two dimensional mesh containing n triangles. We show that for hierarchical adaptive meshes, the algorithm can be locally adapted and easily parallelized by taking advantage of the regularity of the refinement rules. The proposed approach should be very useful in the runtime partitioning and load balancing of adaptive unstructured grids.

Heber, Gerd↗

New plate and shell elements for NASTRAN

A new higher order triangular plate-bending finite element is presented which possesses high accuracy for practical mesh subdivisions and which uses only translations and rotations as grid point degrees of freedom. The element has 18 degrees of freedom, the transverse displacement and two rotations at the vertices and mid-side grid points of the triangle. The transverse displacement within the element is approximated by a quintic polynomial; the bending strains thus vary cubically within the element. Transverse shear flexibility is taken into account in the stiffness formulation. Two examples of static and dynamic analysis are included to show the behavior of the element.

Narayanaswami, R.↗

A Solar Sail Shape Modeling Approach for Attitude Control Design and Analysis

Solar sails operating in the space environment experience deformations in sail shape that result in relatively large disturbance torques which dictate the required performance of the spacecraft attitude control and momentum management systems. These deformations are driven by thermal loads on the booms (due to uneven solar heating), manufacturing and assembly tolerances, and variations in membrane tension. The Solar Cruiser spacecraft utilizes a four-quadrant sail design with four 30-meter length booms and four triangular sail membranes, creating a square sail structure of >1600 m2. Medium-fidelity mesh models were developed based on a characteristic deformed shape. A series of parametric studies were conducted using this shape paradigm to determine worst-case deformed sail shapes which produce bounding disturbance torques. A large database of shapes was produced, and the forces and moments induced by each individual shape were calculated using a Rios-Reyes reduced order generalized sail model. Two were selected as reference worst-case shapes for the Solar Cruiser mission: one which produced the highest pitch/yaw root-sum-squared (RSS) torque, and one which produced the highest roll torque. The results showed that the worst-case shapes at high solar incidence angles induce significantly higher (2-10x) disturbance torques than an ideal, flat-plate sail. Even with considerable safety margins, assuming an ideal sail is unlikely to sufficiently bound the disturbances, which is critical when designing the attitude control system and sizing actuators. Accurate sail shape modeling methodologies should therefore be employed on future solar sail missions.

Solar Sail↗

A Solar Sail Shape Modeling Approach for Attitude Control Design and Analysis

Solar sails operating in the space environment experience deformations in sail shape that result in relatively large disturbance torques which dictate the required performance of the spacecraft attitude control and momentum management systems. These deformations are driven by thermal loads on the booms (due to uneven solar heating), manufacturing and assembly tolerances, and variations in membrane tension. The Solar Cruiser spacecraft utilizes a four-quadrant sail design with four 30-meter length booms and four triangular sail membranes, creating a square sail structure of >1600 m2. Medium-fidelity mesh models were developed based on a characteristic deformed shape. A series of parametric studies were conducted using this shape paradigm to determine worst-case deformed sail shapes which produce bounding disturbance torques. A large database of shapes was produced, and the forces and moments induced by each individual shape were calculated using a Rios-Reyes reduced order generalized sail model [1]. Two were selected as reference worst-case shapes for the Solar Cruiser mission: one which produced the highest pitch/yaw root-sum-squared (RSS) torque, and one which produced the highest roll torque. The results showed that the worst-case shapes at high solar incidence angles induce significantly higher (2-10x) disturbance torques than an ideal, flat-plate sail. Even with considerable safety margins, assuming an ideal sail is unlikely to sufficiently bound the disturbances, which is critical when designing the attitude control system and sizing actuators. Accurate sail shape modeling methodologies should therefore be employed on future solar sail missions. References [1] L. Rios-Reyes and D. J. Scheeres, “Generalized Model for Solar Sails,” Journal of Spacecraft and Rockets, Vol. 42, Jan. 2005, pp. 182–185, 10.2514/1.9054.

Solar Sail↗

A Novel Four-Node Quadrilateral Smoothing Element for Stress Enhancement and Error Estimation

A four-node, quadrilateral smoothing element is developed based upon a penalized-discrete-least-squares variational formulation. The smoothing methodology recovers C1-continuous stresses, thus enabling effective a posteriori error estimation and automatic adaptive mesh refinement. The element formulation is originated with a five-node macro-element configuration consisting of four triangular anisoparametric smoothing elements in a cross-diagonal pattern. This element pattern enables a convenient closed-form solution for the degrees of freedom of the interior node, resulting from enforcing explicitly a set of natural edge-wise penalty constraints. The degree-of-freedom reduction scheme leads to a very efficient formulation of a four-node quadrilateral smoothing element without any compromise in robustness and accuracy of the smoothing analysis. The application examples include stress recovery and error estimation in adaptive mesh refinement solutions for an elasticity problem and an aerospace structural component.

Tessler, A.↗

Analytical integration of the tractions induced by non-singular dislocations on an arbitrary shaped triangular quadratic element

Here, an analytical model is proposed to evaluate the nodal force induced by a segment of dislocation upon an arbitrary shaped triangular element. This calculation is required in hybrid methods that associate dislocation dynamics to boundary or finite element to solve simultaneously the evolution of large ensembles of dislocations with complex boundary conditions. Nodal forces are defined as the triple integration of the unbalanced traction field induced by a straight dislocation upon the surface of the element. Following our previous approach (Queyreau et al 2014 Modelling Simul. Mater. Sci. Eng. 22 035004) on a simpler geometry and in the case of linear isotropic elasticity, triple integrals are solved by sequences of integration by parts that exhibit recurrence relations. The traction field is defined and finite everywhere even at the core of dislocations, thanks to the use of the non-singular stress expression formulated by Cai et al (2006 J. Mech. Phys. Solids 54 561–587). The nodal force expressions can be used when considering both a single convolution or double convolution of the Green's function with the core distribution. A solution is also proposed for the case of a semi-infinite segment through the study of the asymptotic behavior of the analytical expressions. The proposed approach is exact and very computationally efficient. The choice of an arbitrary shaped triangular element and quadratic shape functions allow the consideration of complex geometries and comply with automatic meshing procedures. These analytical expressions could also be employed to estimate dislocation interactions with interfaces.

42 ENGINEERING↗

A structured and unstructured remeshing method for high speed flows

An adaptive remeshing method using both triangular and quadrilateral elements suitable for high speed flows is presented. For inviscid flows the method generates completely unstructured meshes. For viscous flows the boundary layer edge is identified adaptively and a structured mesh is generated in the boundary layer, and an unstructured mesh is generated in the inviscid region. Examples of inviscid and viscous mesh adaptations for high speed flows are presented. A comparison is made between first order and higher order finite element algorithms when used in association with the remeshing method.

Vemaganti, Gururaja R.↗

Finite element analysis of high speed compressible flows using mesh refinement/movement procedures

An adaptive mesh refinement procedure for analyzing high-speed inviscid and viscous compressible flows is described. The adaptation procedure which uses both quadrilateral and triangular elements is implemented with an explicit finite element formulation. Elements in regions of strong and weak gradients are refined or coarsened based on inviscid and viscous indicators. Nodal locations are also adaptively moved to better resolve flow features. The effectiveness of the finite element procedure is demonstrated by modeling flows with complex shock structure and viscous-inviscid interactions. Numerical results are compared with experimental data.

Ramakrishnan, R.↗

Adaptive finite element analysis of hypersonic laminar flows for aerothermal load predictions

The use of an adaptive mesh refinement procedure for analyzing hypersonic laminar flows with application to aerothermal load predictions is described. The adaptation procedure, which uses both quadrilateral and triangular elements, is implemented with the multistep Galerkin-Runge-Kutta scheme. Elements that lie in regions of strong gradients are refined based on indicators to obtain better definition of flow features. The effectiveness of the adaptive procedure is demonstrated by modeling Mach 11.7 flow over a 15-deg ramp. Numerical results are compared with predictions of strong interaction theories and experimental data. Surface quantities such as heating rates and pressure loads, critical for the effective design of high-speed vehicles, are found to be in good agreement with experimental values.

Ramakrishnan, R.↗