Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geometry meshing”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Generation of gear tooth surfaces by application of CNC machines

This study will demonstrate the importance of application of computer numerically controlled (CNC) machines in generation of gear tooth surfaces with new topology. This topology decreases gear vibration and will extend the gear capacity and service life. A preliminary investigation by a tooth contact analysis (TCA) program has shown that gear tooth surfaces in line contact (for instance, involute helical gears with parallel axes, worm gear drives with cylindrical worms, etc.) are very sensitive to angular errors of misalignment that cause edge contact and an unfavorable shape of transmission errors and vibration. The new topology of gear tooth surfaces is based on the localization of bearing contact, and the synthesis of a predesigned parabolic function of transmission errors that is able to absorb a piecewise linear function of transmission errors caused by gear misalignment. The report will describe the following topics: description of kinematics of CNC machines with six degrees of freedom that can be applied for generation of gear tooth surfaces with new topology. A new method for grinding of gear tooth surfaces by a cone surface or surface of revolution based on application of CNC machines is described. This method provides an optimal approximation of the ground surface to the given one. This method is especially beneficial when undeveloped ruled surfaces are to be ground. Execution of motions of the CNC machine is also described. The solution to this problem can be applied as well for the transfer of machine tool settings from a conventional generator to the CNC machine. The developed theory required the derivation of a modified equation of meshing based on application of the concept of space curves, space curves represented on surfaces, geodesic curvature, surface torsion, etc. Condensed information on these topics of differential geometry is provided as well.

Litvin, F. L.↗
NASA NTRS1995

The National Grid Project: A system overview

The National Grid Project (NGP) is a comprehensive numerical grid generation software system that is being developed at the National Science Foundation (NSF) Engineering Research Center (ERC) for Computational Field Simulation (CFS) at Mississippi State University (MSU). NGP is supported by a coalition of U.S. industries and federal laboratories. The objective of the NGP is to significantly decrease the amount of time it takes to generate a numerical grid for complex geometries and to increase the quality of these grids to enable computational field simulations for applications in industry. A geometric configuration can be discretized into grids (or meshes) that have two fundamental forms: structured and unstructured. Structured grids are formed by intersecting curvilinear coordinate lines and are composed of quadrilateral (2D) and hexahedral (3D) logically rectangular cells. The connectivity of a structured grid provides for trivial identification of neighboring points by incrementing coordinate indices. Unstructured grids are composed of cells of any shape (commonly triangles, quadrilaterals, tetrahedra and hexahedra), but do not have trivial identification of neighbors by incrementing an index. For unstructured grids, a set of points and an associated connectivity table is generated to define unstructured cell shapes and neighboring points. Hybrid grids are a combination of structured grids and unstructured grids. Chimera (overset) grids are intersecting or overlapping structured grids. The NGP system currently provides a user interface that integrates both 2D and 3D structured and unstructured grid generation, a solid modeling topology data management system, an internal Computer Aided Design (CAD) system based on Non-Uniform Rational B-Splines (NURBS), a journaling language, and a grid/solution visualization system.

Gaither, Adam↗
NASA NTRS1999

Computer Aided Modeling to Determine the Effectiveness of Resistive Exercises as Countermeasures for Bone Mineral Density Loss

Due to the loss of gravitational loading, astronauts have a tendency to lose bone mineral density in their lumbar spine and lower extremities on orbit. NASA requires astronauts to perform exercises during space flight to help reduce the amount of demineralization. To test these exercises on earth, 17 week bed rest studies are conducted that consist of specific diet and exercise regimes. Developing a finite element model of these exercises will help to quantify the stress distribution imposed by of each of these exercises. To help develop this model, MRI images are acquired from individuals participating in the bed rest studies. The MRIs can be used to create a subject specific model of each individual for testing. The MRIs are processed in the Magnetic Resonance Imaging Data Transfer System program to develop a three-dimensional finite element model of the femur for evaluation. Modifications were made to the MRIDTS that simplified the model creation process. These modifications made it possible to construct two separate models of different portions of a bone simultaneously and then later connect them manually. This helped alleviate the warping problem associated with the drastic changes in geometry found in some body parts, such as the joints. The code was also modified to incorporate material properties of various bone components into the model. Interior meshing was also incorporated into the program to allow for both the cortical shell and the entire bone to be modeled. A prototype model of the right femur of an adult female is being constructed and tested to determine the feasibility of finite element analysis as a tool for evaluating exercise effectiveness. The model is being run through the ANSYS finite element program on the Alabama Super Computer Network. After the model is validated, models of bedrest subjects can be generated to investigate exercise countermeasures.

Murphy, Benjamin M.↗
NASA NTRS2000

An Application of the Quadrature-Free Discontinuous Galerkin Method

The process of generating a block-structured mesh with the smoothness required for high-accuracy schemes is still a time-consuming process often measured in weeks or months. Unstructured grids about complex geometries are more easily generated, and for this reason, methods using unstructured grids have gained favor for aerodynamic analyses. The discontinuous Galerkin (DG) method is a compact finite-element projection method that provides a practical framework for the development of a high-order method using unstructured grids. Higher-order accuracy is obtained by representing the solution as a high-degree polynomial whose time evolution is governed by a local Galerkin projection. The traditional implementation of the discontinuous Galerkin uses quadrature for the evaluation of the integral projections and is prohibitively expensive. Atkins and Shu introduced the quadrature-free formulation in which the integrals are evaluated a-priori and exactly for a similarity element. The approach has been demonstrated to possess the accuracy required for acoustics even in cases where the grid is not smooth. Other issues such as boundary conditions and the treatment of non-linear fluxes have also been studied in earlier work This paper describes the application of the quadrature-free discontinuous Galerkin method to a two-dimensional shear layer problem. First, a brief description of the method is given. Next, the problem is described and the solution is presented. Finally, the resources required to perform the calculations are given.

Lockard, David P.↗
NASA NTRS2010

A New Approach for Constructing Highly Stable High Order CESE Schemes

A new approach is devised to construct high order CESE schemes which would avoid the common shortcomings of traditional high order schemes including: (a) susceptibility to computational instabilities; (b) computational inefficiency due to their local implicit nature (i.e., at each mesh points, need to solve a system of linear/nonlinear equations involving all the mesh variables associated with this mesh point); (c) use of large and elaborate stencils which complicates boundary treatments and also makes efficient parallel computing much harder; (d) difficulties in applications involving complex geometries; and (e) use of problem-specific techniques which are needed to overcome stability problems but often cause undesirable side effects. In fact it will be shown that, with the aid of a conceptual leap, one can build from a given 2nd-order CESE scheme its 4th-, 6th-, 8th-,... order versions which have the same stencil and same stability conditions of the 2nd-order scheme, and also retain all other advantages of the latter scheme. A sketch of multidimensional extensions will also be provided.

Chang, Sin-Chung↗
NASA NTRS2015

Evaluation of the Anisotropic Radiative Conductivity of a Low-Density Carbon Fiber Material from Realistic Microscale Imaging

The radiative heat transfer inside a low-density carbon fiber insulator is analyzed using a three-dimensional direct simulation model. A robust procedure is presented for the numerical calculation of the geometric configuration factor to compute the radiative energy exchange processes among the small discretized surface areas of the fibrous material. The methodology is applied to a polygonal mesh of a fibrous insulator obtained from three-dimensional microscale imaging of the real material. The anisotropic values of the radiative conductivity are calculated for that geometry. The results yield both directional and thermal dependence of the radiative conductivity.

Evaluation↗
DOE OSTI2012

Modelling the graphite fracture mechanisms

In order to define a design criterion for graphite components, it is important to identify the physical phenomena responsible for the graphite fracture, to include them in a more effective modelling. In a first step, a large panel of experiments have been realised in order to build up an important database; results of tensile tests, 3 and 4 point bending tests on smooth and notched specimens have been analysed and have demonstrated an important geometry related effects on the behavior up to fracture. Then, first simulations with an elastic or an elastoplastic bilinear constitutive law have not made it possible to simulate the experimental fracture stress variations with the specimen geometry, the fracture mechanisms of the graphite being at the microstructural scale. That is the reason why a specific F.E. model of the graphite structure has been developed in which every graphite grain has been meshed independently, the crack initiation along the basal plane of the particles as well as the crack propagation and coalescence have been modelled too. This specific model has been used to test two different approaches for fracture initiation: a critical stress criterion and two criteria of fracture mechanic type. They are all based on crystallographic considerations as a global critical stress criterion gave unsatisfactory results. The criteria of fracture mechanic type being extremely unstable and unable to represent the graphite global behaviour up to the final collapse, the critical stress criterion has been preferred to predict the results of the large range of available experiments, on both smooth and notched specimens. In so doing, the experimental observations have been correctly simulated: the geometry related effects on the experimental fracture stress dispersion, the specimen volume effects on the macroscopic fracture stress and the crack propagation at a constant stress intensity factor. In addition, the parameters of the criterion have been related to experimental observations: the local crack initiation stress of 8 MPa corresponds to the non-linearity apparition on the global behavior observed experimentally and the the maximal critical stress defined for the particle of 30 MPa is equivalent to the fracture stress of notched specimens. This innovative combination of crack modelling and a local crystallographic critical stress criterion made it possible to understand that cleavage initiation and propagation in the graphite microstructure was driven by a mean critical stress criterion. (authors)

36 MATERIALS SCIENCE↗
DOE OSTI2023

Preserving Superconvergence of Spectral Elements for Curved Domains via $h$ and $p$-Geometric Refinement

Spectral element methods (SEM), extensions of finite element methods (FEM), have emerged as significant techniques for solving partial differential equations in physics and engineering. SEM can potentially deliver superior accuracy due to the potential superconvergence in nodal solutions for well-shaped tensor-product elements. However, the accuracy of SEM often degrades in complex geometries due to geometric inaccuracies near curved boundaries and the loss of superconvergence with simplicial or non-tensor-product elements. To overcome the first issue, we propose using h- and p-geometric refinement, which refines the mesh near high-curvature regions and increases the degree of geometric basis functions, respectively. We show that when using mixed-element meshes with tensor-product elements in the interior of the domain, curvature-based geometric refinement near boundaries can improve the accuracy of the interior elements by reducing pollution errors and preserving the superconvergence in nodal solutions. To address the second issue, we introduce a post-processing technique using the adaptive extended stencil finite element method (AES-FEM) to recover the accuracy near the curved boundaries. The combination of curvature-based geometric refinement and accurate post-processing offers an effective and easier-to-implement alternative to methods reliant on exact geometries. We demonstrate our techniques by solving the convection-diffusion equation in 2D and show up to two orders of magnitude of improvement in the solution accuracy, even when the elements are poorly shaped near boundaries.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Reducing Ohmic Resistances in Membrane Capacitive Deionization Using Micropatterned Ion-Exchange Membranes, Ionomer Infiltrated Electrodes, and Ionomer-Coated Nylon Meshes

Membrane capacitive deionization (MCDI) is an emerging water desalination platform that is compact, electrified, and does not require high-pressure piping. Herein, highly conductive poly(phenylene alkylene) ion-exchange membranes (IEMs) are micropatterned with different surface geometries for MCDI. The micropatterned membranes increase the interfacial area with the liquid stream leading to a 700 mV reduction in cell voltage when operating at constant current (2 mA cm -2 ; 2000 ppm NaCl feed) while improving the energy normalized adsorbed salt (ENAS) value by 1.4 times. Combining the micropatterned poly(phenylene alkylene) IEMs with poly(phenylene alkylene) ionomer-filled electrodes reduces the cell voltage by 1000 mV and improves the ENAS values by 2.3 times relative to the base case. This reduction in cell voltage allows for higher current density operation (i.e., 3–4 mA cm -2 ) . The reduction in cell voltage is ascribed to the ameliorating ohmic resistances related to ion transport at the membrane-process stream interface and in the carbon cloth electrode. Finally, porous ionic conductors are implemented into the spacer channel with flat and micropatterned IEM configurations and ionomer infiltrated electrodes. For the configuration with flat IEMs, the porous ionic conductor improves ENAS values across the current density regime (2–4 mA cm -2 ), while for micropatterned IEMs it gets improved only at 4 mA cm -2 .

42 ENGINEERING↗
NASA NTRS2000

Complicated Geometry

Using NASA SBIR funding, CFD Research Corporation has developed CFD-GEOM, an extension of traditional computer-aided drawing (CAD) software. CFD-GEOM can provide modeling and interactivity of computational fluid dynamics (CFD) latest field, mesh-generation and allows for quick and easy updating of a grid in response to changes in the CAD model.

Source record↗
DOE OSTI2012

Analysis of the pool critical assembly pressure vessel benchmark using pentran

The internationally circulated Pool Critical Assembly (PCA) Pressure Vessel Benchmark was analyzed using the PENTRAN Parallel Sn code system for the geometry, material, and source specifications as described in the PCA Benchmark documentation. This research focused on utilizing the BUGLE-96 cross section library and accompanying reaction rates, while examining both adaptive differencing on a coarse mesh basis as well as Directional Theta Weighted Sn differencing in order to compare the calculated PENTRAN results to measured data. The results show good comparison with the measured data as well as to the calculated results provided from TORT for the BUGLE-96 cross sections and reaction rates, which suggests PENTRAN is a viable and reliable code system for calculation of light water reactor neutron shielding and dosimetry calculations. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1995

Heat transfer predictions for two turbine nozzle geometries at high Reynolds and Mach numbers

Predictions of turbine vane and endwall heat transfer and pressure distributions are compared with experimental measurements for two vane geometries. The differences in geometries were due to differences in the hub profile, and both geometries were derived from the design of a high rim speed turbine (HRST). The experiments were conducted in the Isentropic Light Piston Facility (ILPF) at Pyestock at a Reynolds number of 5.3 x 10(exp 6), a Mach number of 1.2, and a wall-to-gas temperature ratio of 0.66. Predictions are given for two different steady-state three-dimensional Navier-Stokes computational analyses. C-type meshes were used, and algebraic models were employed to calculate the turbulent eddy viscosity. The effects of different turbulence modeling assumptions on the predicted results are examined. Comparisons are also given between predicted and measured total pressure distributions behind the vane. The combination of realistic engine geometries and flow conditions proved to be quite demanding in terms of the convergence of the CFD solutions. An appropriate method of grid generation, which resulted in consistently converged CFD solutions, was identified.

Boyle, R. J.↗
DOE OSTI2021

CFD Simulation of Helium Flow Loop Test Section

A helium flow loop is being assembled at Oak Ridge National Laboratory to analyze heat transfer enhancement for systems such as blanket and divertor components. To efficiently identify optimum geometries for heat transfer enhancement in these applications, simulation work is performed to optimize test section designs that are built and tested in the helium flow loop that operates at 4 MPa and a mass flow rate of 100 g/s. Different ribbed geometries that examine rib shape, rib height, rib orientation, rib spacing, and three dimensional orientation are modeled and simulated in STAR-CCM+ to compare their ability to remove heat and mitigate pressure drop. Following the simulations, models are selected and manufactured for the helium flow loop tests. Simulations initially focus on a hydrodynamic study to determine the appropriate mesh and physics models and then add a heat flux to analyze the heat transfer abilities of the models. The simulations are run in steady state and use a Reynolds-averaged Navier-Stokes k-ε turbulence model. The helium is modeled as an ideal gas. The simulation explores models of geometries that enhance the heat transfer and decrease pressure drop with an overall goal of increasing fluid collision with the wall. Enhanced geometries are simulated to select appropriate designs for manufacturing, and preliminary experimental results are used to validate the simulations. Furthermore, the factors that are being analyzed in the comparison between the experimental and the simulated results include matching thermocouple temperatures, pressure drop, roughness, and fluid velocity.

42 ENGINEERING↗
DOE OSTI2016

Arbitrary-Order Conservative and Consistent Remapping and a Theory of Linear Maps: Part II

The focus on this series of articles is on the generation of accurate, conservative, consistent, and (optionally) monotone linear offline maps. This paper is the second in the series. It extends on the first part by describing four examples of 2D linear maps that can be constructed in accordance with the theory of the earlier work. The focus is again on spherical geometry, although these techniques can be readily extended to arbitrary manifolds. The four maps include conservative, consistent, and (optionally) monotone linear maps (i) between two finite-volume meshes, (ii) from finite-volume to finite-element meshes using a projection-type approach, (iii) from finite-volume to finite-element meshes using volumetric integration, and (iv) between two finite-element meshes. Arbitrary order of accuracy is supported for each of the described nonmonotone maps.

58 GEOSCIENCES↗
DOE OSTI2025

CaptureFOAM

CaptureFOAM provides utilities for fast and accurate direct air capture (DAC) contactor simulations, based on the OpenFOAM framework. At the core, CaptureFOAM is a suite of mass transfer and chemistry models built upon the OpenFOAM thin film model, which allows coarser mesh resolution and less constrictive time step requirements when compared to a traditional volume of fluid approach. This allows CaptureFOAM to conduct device-scale simulations for complex geometries including commercial packing materials and novel contactor shapes, matching experimental results in a wide variety of applications.

Kincaid, Kellis [Oak Ridge National Laboratory (OR↗
NASA NTRS1987

Composite grid and finite-volume LU implicit scheme for turbine flow analysis

A composite grid was generated in an attempt to improve grid quality for a typical turbine blade with large camber in terms of mesh control, smoothness, and orthogonality. This composite grid consists of the C grid (or O grid) in the immediate vicinity of the blade and the H grid in the upstream region and in the middle of the blade passage between the C grids. It provides a good boundary layer resolution around the leading edge region for viscous calculation, has orthogonality at the blade surface and slope continuity at the C-H (or O-H) interface, and has flexibility in controlling the mesh distribution in the upstream region without using excessive grid points. This composite grid eliminates the undesirable qualities of a single grid when generated for a typical turbine geometry. A finite-volume lower-upper (LU) implicit scheme can be used in solving for the turbine flows on the composite grid. This grid has a special grid node that is connected to more than four neighboring nodes in two dimensions and to more than six nodes in three dimensions. But the finite-volume approach poses no problem at the special point because each interior cell has only four neighboring cells in two dimensions and only six cells in three dimensions. The finite-volume LU implicit scheme was demonstrated to be robust and efficient for both external and internal flows in a broad flow regime.

Choo, Yung K.↗
NASA NTRS1987

Composite grid and finite-volume LU implicit scheme for turbine flow analysis

A composite grid was generated in an attempt to improve grid quality for a typical turbine blade with large camber in terms of mesh control, smoothness, and orthogonality. This composite grid consists of the C grid (or O grid) in the immediate vicinity of the blade and the H grid in the upstream region and in the middle of the blade passage between the C grids. It provides a good boundary layer resolution around the leading edge region for viscous calculation, has orthogonality at the blade surface and slope continuity at the C-H (or O-H) interface, and has flexibility in controlling the mesh distribution in the upstream region without using excessive grid points. This composite grid eliminates the undesirable qualities of a single grid when generated for a typical turbine geometry. A finite-volume lower-upper (LU) implicit schemes can be used in solving for the turbine flows on the composite grid. This grid has a special grid node that is connected to more than four neighboring nodes in two dimensions and to more than six nodes in three dimensions. But the finite-volume approach poses no problem at the special point because each interior cell has only four neighboring cells in two dimensions and only six cells in three dimensions. The finite-volume LU implicit scheme was demonstrated to be robust and efficient for both external and internal flows in a broad flow regime.

Choo, Yung K.↗
NASA NTRS1993

BLASIM - A computational tool to assess ice impact damage on engine blades

A portable computer code called BLASIM is developed at NASA LeRC to assess the ice impact damage on aircraft engine blades. In addition to the ice impact analyses, the code is also capable of carrying out static, dynamic, resonance margin and flutter analyses. The blade can be solid, hollow, superhybrid or composite material. An optional preprocessor (input generator) is also developed to generate input to the code through interactive process. The blade geometry can be defined either by a series of airfoils at discrete input stations or by a finite element grid. The code employs a coarse fixed finite element mesh with triangular plate finite elements and has quick turnaround time. The ice piece is modeled as an equivalent spherical object and has the velocity opposite to that of the aircraft with direction parallel to the engine axis. For the local impact damage assessment, the impact force is considered as a distributed load acting over a region around the impact point and the average radial strain of the finite elements along the leading edge is taken as a measure of the local damage. To estimate the damage at the blade root, the impact is considered to be an impulse and a combined stress failure criteria is employed. Parametric studies for local and root ice impact damage, and post-impact dynamics are discussed for solid and composite blades.

Reddy, E. S.↗