Search NASA⌕ Search

SEARCH · Search NASA

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 199 records · Page 11

Offset Compound Gear Drive

The Offset Compound Gear Drive is an in-line, discrete, two-speed device utilizing a special offset compound gear that has both an internal tooth configuration on the input end and external tooth configuration on the output end, thus allowing it to mesh in series, simultaneously, with both a smaller external tooth input gear and a larger internal tooth output gear. This unique geometry and offset axis permits the compound gear to mesh with the smaller diameter input gear and the larger diameter output gear, both of which are on the same central, or primary, centerline. This configuration results in a compact in-line reduction gear set consisting of fewer gears and bearings than a conventional planetary gear train. Switching between the two output ratios is accomplished through a main control clutch and sprag. Power flow to the above is transmitted through concentric power paths. Low-speed operation is accomplished in two meshes. For the purpose of illustrating the low-speed output operation, the following example pitch diameters are given. A 5.0 pitch diameter (PD) input gear to 7.50 PD (internal tooth) intermediate gear (0.667 reduction mesh), and a 7.50 PD (external tooth) intermediate gear to a 10.00 PD output gear (0.750 reduction mesh). Note that it is not required that the intermediate gears on the offset axis be of the same diameter. For this example, the resultant low-speed ratio is 2:1 (output speed = 0.500; product of stage one 0.667 reduction and stage two 0.750 stage reduction). The design is not restricted to the example pitch diameters, or output ratio. From the output gear, power is transmitted through a hollow drive shaft, which, in turn, drives a sprag during which time the main clutch is disengaged.

Stevens, Mark A.↗

Unstructured and adaptive mesh generation for high Reynolds number viscous flows

A method for generating and adaptively refining a highly stretched unstructured mesh suitable for the computation of high-Reynolds-number viscous flows about arbitrary two-dimensional geometries was developed. The method is based on the Delaunay triangulation of a predetermined set of points and employs a local mapping in order to achieve the high stretching rates required in the boundary-layer and wake regions. The initial mesh-point distribution is determined in a geometry-adaptive manner which clusters points in regions of high curvature and sharp corners. Adaptive mesh refinement is achieved by adding new points in regions of large flow gradients, and locally retriangulating; thus, obviating the need for global mesh regeneration. Initial and adapted meshes about complex multi-element airfoil geometries are shown and compressible flow solutions are computed on these meshes.

Mavriplis, Dimitri J.↗

Unstructured and adaptive mesh generation for high Reynolds number viscous flows

A method for generating and adaptively refining a highly stretched unstructured mesh suitable for the computation of high-Reynolds-number viscous flows about arbitrary two-dimensional geometries was developed. The method is based on the Delaunay triangulation of a predetermined set of points and employs a local mapping in order to achieve the high stretching rates required in the boundary-layer and wake regions. The initial mesh-point distribution is determined in a geometry-adaptive manner which clusters points in regions of high curvature and sharp corners. Adaptive mesh refinement is achieved by adding new points in regions of large flow gradients, and locally retriangulating; thus, obviating the need for global mesh regeneration. Initial and adapted meshes about complex multi-element airfoil geometries are shown and compressible flow solutions are computed on these meshes.

Mavriplis, D. J.↗

An unstructured body-of-revolution electromagnetic particle-in-cell algorithm with radial perfectly matched layers and dual polarizations

A novel electromagnetic particle-in-cell algorithm has been developed for fully kinetic plasma simulations on unstructured (irregular) meshes in complex body-of-revolution geometries. The algorithm, implemented in the BORPIC++ code, utilizes a set of field scalings and a coordinate mapping, reducing the Maxwell field problem in a cylindrical system to a Cartesian finite element Maxwell solver in the meridian plane. The latter obviates the cylindrical coordinate singularity in the symmetry axis. The choice of an unstructured finite element discretization enhances the geometrical flexibility of the BORPIC++ solver compared to the more traditional finite difference solvers. Symmetries in Maxwell’s equations are explored to decompose the problem into two dual polarization states with isomorphic representations that enable code reuse. The particle-in-cell scatter and gather steps preserve charge conservation at the discrete level. Our previous algorithm (BORPIC+) discretized the E and B field components of TE Φ and TM Φ polarizations on the finite element (primal) mesh. Here, we employ a new field-update scheme. Using the same finite element (primal) mesh, this scheme advances two sets of field components independently: (1) E and B of TE Φ polarized fields, (E z , E ρ , B Φ ) and (2) D and H of TM Φ polarized fields, (D Φ , H z , H ρ ). Since these field updates are not explicitly coupled, the new field solver obviates the coordinate singularity, which otherwise arises at the cylindrical symmetric axis, ρ = 0 when defining the discrete Hodge matrices (generalized finite element mass matrices). Here, a cylindrical perfectly matched layer is implemented as a boundary condition in the radial direction to simulate open space problems, with periodic boundary conditions in the axial direction. We investigate effects of charged particles moving next to the cylindrical perfectly matched layer. We model azimuthal currents arising from rotational motion of charged rings, which produce TMΦ polarized fields. Several numerical examples are provided to illustrate the first application of the algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Towards a NEAMS-based high-fidelity model of the MARVEL reactor

This report outlines the progress of Idaho National Laboratory in developing a high-fidelity and high-resolution model of the Microreactor Applications Research Validation and Evaluation reactor. The model was developed under the Nuclear Energy Advanced Modeling and Simulation microreactor application driver at Idaho National Laboratory. The overarching objective of this activity is the development of a high-fidelity multiphysics MARVEL model using NEAMS tools, and to verify and validate NEAMS tools against MARVEL reference simulation and experimental data, respectively. This is a unique opportunity to conduct multiphysics analysis on a soon-to-be-deployed microreactor. This multiphysics model developed under the NEAMS-funded INL microreactor application driver leverages three single-physics models coupled via the MOOSE’s MultiApp and Transfer systems. The latter systems enable in-memory data transfer between MOOSE-based and MOOSE-wrapped applications. The first single-physics model, that functions as main application, leverages Griffin to model the neutron transport in the core through the discontinuous finite element (DFEM) discrete ordinates solver (SN). Several optimization flags that were developed by the Griffin developer team were beta-tested to enhance the solver’s performance. These include the combined use of using_average_xs and update_averaged_xs_on that enable to avoid expensive on-the-fly cross sections evaluations at each linear iterations in favor of evaluations of the macroscopic cross sections at each Picard iteration. The second single-physics model uses BISON to handle solid heat transfer and asymptotic hydrogen redistribution analysis in the fuel. While the model returns consistent results for the temperature and hydrogen distribution in the fuel, a mismatch was noticed in the calculated temperature in the reflector due to the value of the gap conductance used in our model. Ongoing investigations are being performed to assess the origin of this discrepancy. Finally, the System Analysis Module (SAM) was used to model the flow of the sodium-potassium eutectic in the primary loop. A first verification was also performed showing good agreement in terms of mass flow rate and inlet temperature. All mesh files were generated using the MOOSE Reactor module, removing the need for external meshing tools. Notably, this workscope represents one of the initial applications of the MOOSE Reactor module for modeling highly irregular geometries. The use of the reactor module significantly streamlined the mesh generation process. The full multiphysics mode, that combines all the single physics models, was leveraged to conduct initial steady-state multiphysics simulations to compute power, and temperature distribution in the reactor. Initial testing was performed for transient simulations as well. In this case, the new checkpoint restart capability for eigenvalue calculations was tested showing the capability for streamlined restart of transient calculations. Future work will focus on improving the fidelity of the model by performing comprehensive code-to-code comparisons. For instance, the full-core Griffin neutronics model will be benchmarked against MCNP reference results, that were provided by the MARVEL design team. Additionally, the SAM T/H model will be verified against reference RELAP-5 results for selected accident scenarios. Besides code-to-code verification exercises, the model fidelity will be improved by replacing the single-channel SAM model with a more complex SAM-Pronghorn coupled model, in which the sub-channel capability is deployed to obtain radial temperature resolution in the coolant. This model will be developed in synergy with the NEAMS thermal hydraulics team.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Image-Based Computational Fluid Dynamics in Blood Vessel Models: Toward Developing a Prognostic Tool to Assess Cardiovascular Function Changes in Prolonged Space Flights

One of NASA's objectives is to be able to perform a complete, pre-flight, evaluation of cardiovascular changes in astronauts scheduled for prolonged space missions. Computational fluid dynamics (CFD) has shown promise as a method for estimating cardiovascular function during reduced gravity conditions. For this purpose, MRI can provide geometrical information, to reconstruct vessel geometries, and measure all spatial velocity components, providing location specific boundary conditions. The objective of this study was to investigate the reliability of MRI-based model reconstruction and measured boundary conditions for CFD simulations. An aortic arch model and a carotid bifurcation model were scanned in a 1.5T Siemens MRI scanner. Axial MRI acquisitions provided images for geometry reconstruction (slice thickness 3 and 5 mm; pixel size 1x1 and 0.5x0.5 square millimeters). Velocity acquisitions provided measured inlet boundary conditions and localized three-directional steady-flow velocity data (0.7-3.0 L/min). The vessel walls were isolated using NIH provided software (ImageJ) and lofted to form the geometric surface. Constructed and idealized geometries were imported into a commercial CFD code for meshing and simulation. Contour and vector plots of the velocity showed identical features between the MRI velocity data, the MRI-based CFD data, and the idealized-geometry CFD data, with less than 10% differences in the local velocity values. CFD results on models reconstructed from different MRI resolution settings showed insignificant differences (less than 5%). This study illustrated, quantitatively, that reliable CFD simulations can be performed with MRI reconstructed models and gives evidence that a future, subject-specific, computational evaluation of the cardiovascular system alteration during space travel is feasible.

Chatzimavroudis, George P.↗

Tangle-Free Finite Element Mesh Motion for Ablation Problems

Mesh motion is the process by which a computational domain is updated in time to reflect physical changes in the material the domain represents. Such a technique is needed in the study of the thermal response of ablative materials, which erode when strong heating is applied to the boundary. Traditionally, the thermal solver is coupled with a linear elastic or biharmonic system whose sole purpose is to update mesh node locations in response to altering boundary heating. Simple mesh motion algorithms rely on boundary surface normals. In such schemes, evolution in time will eventually cause the mesh to intersect and "tangle" with itself, causing failure. Furthermore, such schemes are greatly limited in the problems geometries on which they will be successful. This paper presents a comprehensive and sophisticated scheme that tailors the directions of motion based on context. By choosing directions for each node smartly, the inevitable tangle can be completely avoided and mesh motion on complex geometries can be modeled accurately.

Droba, Justin↗

Spatial modeling algorithms for reactions and transport in biological cells

Biological cells rely on precise spatiotemporal coordination of biochemical reactions to control their functions. Such cell signaling networks have been a common focus for mathematical models, but they remain challenging to simulate, particularly in realistic cell geometries. Here we present Spatial Modeling Algorithms for Reactions and Transport (SMART), a software package that takes in high-level user specifications about cell signaling networks and then assembles and solves the associated mathematical systems. SMART uses state-of-the-art finite element analysis, via the FEniCS Project software, to efficiently and accurately resolve cell signaling events over discretized cellular and subcellular geometries. We demonstrate its application to several different biological systems, including yes-associated protein (YAP)/PDZ-binding motif (TAZ) mechanotransduction, calcium signaling in neurons and cardiomyocytes, and ATP generation in mitochondria. Throughout, we utilize experimentally derived realistic cellular geometries represented by well-conditioned tetrahedral meshes. These scenarios demonstrate the applicability, flexibility, accuracy and efficiency of SMART across a range of temporal and spatial scales.

59 BASIC BIOLOGICAL SCIENCES↗

Recent developments in finite element analysis for transonic airfoils

The prediction of aerodynamic forces in the transonic regime generally requires a flow field calculation to solve the governing non-linear mixed elliptic-hyperbolic partial differential equations. Finite difference techniques were developed to the point that design and analysis application are routine, and continual improvements are being made by various research groups. The principal limitation in extending finite difference methods to complex three-dimensional geometries is the construction of a suitable mesh system. Finite element techniques are attractive since their application to other problems have permitted irregular mesh elements to be employed. The purpose of this paper is to review the recent developments in the application of finite element methods to transonic flow problems and to report some recent results.

Hafez, M. M.↗

Numerical simulation of the tip vortex off a low-aspect-ratio wing at transonic speed

The viscous transonic flow around a low aspect ratio wing was computed by an implicit, three dimensional, thin-layer Navier-Stokes solver. The grid around the geometry of interest is obtained numerically as a solution to a Dirichlet problem for the cube. A low aspect ratio wing with large sweep, twist, taper, and camber is the chosen geometry. The topology chosen to wrap the mesh around the wing with good tip resolution is a C-O type mesh. The flow around the wing was computed for a free stream Mach number of 0.82 at an angle of attack of 5 deg. At this Mach number, an oblique shock forms on the upper surface of the wing, and a tip vortex and three dimensional flow separation off the wind surface are observed. Particle path lines indicate that the three dimensional flow separation on the wing surface is part of the roots of the tip vortex formation. The lifting of the tip vortex before the wing trailing edge is observed by following the trajectory of particles release around the wing tip.

Mansour, N. N.↗

Stress analyses of B-52 pylon hooks

The NASTRAN finite element computer program was used in the two dimensional stress analysis of B-52 carrier aircraft pylon hooks: (1) old rear hook (which failed), (2) new rear hook (improved geometry), (3) new DAST rear hook (derated geometry), and (4) front hook. NASTRAN model meshes were generated by the aid of PATRAN-G computer program. Brittle limit loads for all the four hooks were established. The critical stress level calculated from NASTRAN agrees reasonably well with the values predicted from the fracture mechanics for the failed old rear hook.

Ko, W. L.↗

Implicit time-marching solution of the Navier-Stokes equations for thrust reversing and thrust vectoring nozzle flows

An implicit finite volume method is investigated for the solution of the compressible Navier-Stokes equations for flows within thrust reversing and thrust vectoring nozzles. Thrust reversing nozzles typically have sharp corners, and the rapid expansion and large turning angles near these corners are shown to cause unacceptable time step restrictions when conventional approximate factorization methods are used. In this investigation these limitations are overcome by using second-order upwind differencing and line Gauss-Siedel relaxation. This method is implemented with a zonal mesh so that flows through complex nozzle geometries may be efficiently calculated. Results are presented for five nozzle configurations including two with time varying geometries. Three cases are compared with available experimental data and the results are generally acceptable.

Imlay, S. T.↗

RCC Plug Repair Thermal Tools for Shuttle Mission Support

A thermal math model for the Space Shuttle Reinforced Carbon-Carbon (RCC) Plug Repair was developed to increase the confidence in the repair entry performance and provide a real-time mission support tool. The thermal response of the plug cover plate, local RCC, and metallic attach hardware can be assessed with this model for any location on the wing leading edge. The geometry and spatial location of the thermal mesh also matches the structural mesh which allows for the direct mapping of temperature loads and computation of the thermoelastic stresses. The thermal model was correlated to a full scale plug repair radiant test. To utilize the thermal model for flight analyses, accurate predictions of protuberance heating were required. Wind tunnel testing was performed at CUBRC to characterize the heat flux in both the radial and angular directions. Due to the complexity of the implementation of the protuberance heating, an intermediate program was developed to output the heating per nodal location for all OML surfaces in SINDA format. Three Design Reference Cases (DRC) were evaluated with the correlated plug thermal math model to bound the environments which the plug repair would potentially be used.

Rodriguez, Alvaro C.↗

A Numerical Investigation of Parachute Deployment in Supersonic Flow

A validated computational fluid-structure interaction method applied to supersonic parachute inflation is extended to consider a unique, parallel self-contact algorithm, porosity interface conditions on the canopy, improved flow field sampling procedures for obtaining quality loading on the parachute canopy, adaptive mesh refinement, and improved treatment of thin geometries in an immersed boundary framework. These extensions are discussed in detail and demonstrated individually on test problems. Finally, the developments are brought together for demonstration on a sub-scale MSL parachute geometry.

Boustani, Jonathan↗

A Concept for the Deployment of a Large Lunar Crater Radio Telescope Using Teams of Tethered Robots

Kilometer-scale craters on the far side of the Moon have unique potential as future locations for large radio telescopes, which can observe the universe at wavelengths and frequencies (> 10 m, < 30 MHz) not possible with conventional Earth or orbital-based approaches. Distinct advantages of building a Lunar Crater Radio Telescope (LCRT) on the far side include i) isolation from radio noise due to the Earth’s ionosphere, orbiting satellites, and the Sun, ii) days of uninterrupted dark/cold sky viewing during lunar night, and iii) terrain geometry naturally suited for constructing the largest mesh antenna structure in the Solar System. A key challenge to constructing LCRT on the Moon is related to the complexity of deploying a 1-km diameter antenna and hanging receiver within a lunar crater whose diameter, depth, and slope are 3-5 km, 1 km, and ~30 degrees respectively. In this paper, we first evaluate the trade space for deploying a large, complex structure within a crater, and then provide a more detailed concept evaluation of our favored approach, which employs coordinated teams of tethered rovers to extract and suspend a folded antenna from a lander at the base of a crater. NASA’s Jet Propulsion Laboratory in collaboration with California Institute of Technology (Caltech) have developed a novel robotic system for accessing extremely steep terrains; the Axel rover is a two-wheeled rugged terrain vehicle that is supported by an electro-mechanical tether that provides power, data, and tensile support from a top-side anchor location. Recently, a pair of Axel robots have been used in a DuAxel configuration that allows for four-wheel driving and repeated passive anchoring at different locations. The DuAxel system has unique advantages for deploying an LCRT antenna, including the ability to deploy from a lander near a crater, drive a distance to the crater rim to deploy an Axel, and later, retract the deployed Axel in order to sequentially lift up sections of the antenna. Our proposed concept involves delivering a packaged antenna and receiver to the bottom-center of a crater floor on a lander, then later sending a team of multiple DuAxel rovers to retrieve guide wires from the lander, which are pulled to the top of the crater. We explore this concept in detail and provide some initial quantitative analysis to demonstrate the feasibility of our system with respect to the spatial and mass properties of the antenna as juxtaposed to DuAxel capabilities. Finally, we outline next steps towards validating our concept on the way to a future lunar deployment opportunity.

Hallinan, Greg↗

Tangle-Free Mesh Motion for Ablation Simulations

Problems involving mesh motion-which should not be mistakenly associated with moving mesh methods, a class of adaptive mesh redistribution techniques-are of critical importance in numerical simulations of the thermal response of melting and ablative materials. Ablation is the process by which material vaporizes or otherwise erodes due to strong heating. Accurate modeling of such materials is of the utmost importance in design of passive thermal protection systems ("heatshields") for spacecraft, the layer of the vehicle that ensures survival of crew and craft during re-entry. In an explicit mesh motion approach, a complete thermal solve is first performed. Afterwards, the thermal response is used to determine surface recession rates. These values are then used to generate boundary conditions for an a posteriori correction designed to update the location of the mesh nodes. Most often, linear elastic or biharmonic equations are used to model this material response, traditionally in a finite element framework so that complex geometries can be simulated. A simple scheme for moving the boundary nodes involves receding along the surface normals. However, for all but the simplest problem geometries, evolution in time following such a scheme will eventually bring the mesh to intersect and "tangle" with itself, inducing failure. This presentation demonstrates a comprehensive and sophisticated scheme that analyzes the local geometry of each node with help from user-provided clues to eliminate the tangle and enable simulations on a wide-class of difficult problem geometries. The method developed is demonstrated for linear elastic equations but is general enough that it may be adapted to other modeling equations. The presentation will explicate the inner workings of the tangle-free mesh motion algorithm for both two and three-dimensional meshes. It will show abstract examples of the method's success, including a verification problem that demonstrates its accuracy and correctness. The focus of the presentation will be on the algorithm; specifics on how the techniques may be used in spacecraft design will be not discussed.

Droba, Justin↗

Using SpaceClaimTD Direct for Modeling Components with Complex Geometries for the Thermal Desktop-Based Advanced Stirling Radioisotope Generator Model

SpaceClaim/TD Direct has been used extensively in the development of the Advanced Stirling Radioisotope Generator (ASRG) thermal model. This paper outlines the workflow for that aspect of the task and includes proposed best practices and lessons learned. The ASRG thermal model was developed to predict component temperatures and power output and to provide insight into the prime contractor's thermal modeling efforts. The insulation blocks, heat collectors, and cold side adapter flanges (CSAFs) were modeled with this approach. The model was constructed using mostly TD finite difference (FD) surfaces/solids. However, some complex geometry could not be reproduced with TD primitives while maintaining the desired degree of geometric fidelity. Using SpaceClaim permitted the import of original CAD files and enabled the defeaturing/repair of those geometries. TD Direct (a SpaceClaim add-on from CRTech) adds features that allowed the "mark-up" of that geometry. These so-called "mark-ups" control how finite element (FE) meshes are to be generated through the "tagging" of features (e.g. edges, solids, surfaces). These tags represent parameters that include: submodels, material properties, material orienters, optical properties, and radiation analysis groups. TD aliases were used for most tags to allow analysis to be performed with a variety of parameter values. "Domain-tags" were also attached to individual and groups of surfaces and solids to allow them to be used later within TD to populate objects like, for example, heaters and contactors. These tools allow the user to make changes to the geometry in SpaceClaim and then easily synchronize the mesh in TD without having to redefine the objects each time as one would if using TDMesher. The use of SpaceClaim/TD Direct helps simplify the process for importing existing geometries and in the creation of high fidelity FE meshes to represent complex parts. It also saves time and effort in the subsequent analysis.

Power Systems↗

Using SpaceClaim/TD Direct for Modeling Components with Complex Geometries for the Thermal Desktop-Based Advanced Stirling Radioisotope Generator Model

SpaceClaim/TD Direct has been used extensively in the development of the Advanced Stirling Radioisotope Generator (ASRG) thermal model. This paper outlines the workflow for that aspect of the task and includes proposed best practices and lessons learned. The ASRG thermal model was developed to predict component temperatures and power output and to provide insight into the prime contractors thermal modeling efforts. The insulation blocks, heat collectors, and cold side adapter flanges (CSAFs) were modeled with this approach. The model was constructed using mostly TD finite difference (FD) surfaces solids. However, some complex geometry could not be reproduced with TD primitives while maintaining the desired degree of geometric fidelity. Using SpaceClaim permitted the import of original CAD files and enabled the defeaturing repair of those geometries. TD Direct (a SpaceClaim add-on from CRTech) adds features that allowed the mark-up of that geometry. These so-called mark-ups control how finite element (FE) meshes were generated and allowed the tagging of features (e.g. edges, solids, surfaces). These tags represent parameters that include: submodels, material properties, material orienters, optical properties, and radiation analysis groups. TD aliases were used for most tags to allow analysis to be performed with a variety of parameter values. Domain-tags were also attached to individual and groups of surfaces and solids to allow them to be used later within TD to populate objects like, for example, heaters and contactors. These tools allow the user to make changes to the geometry in SpaceClaim and then easily synchronize the mesh in TD without having to redefine these objects each time as one would if using TD Mesher.The use of SpaceClaim/TD Direct has helped simplify the process for importing existing geometries and in the creation of high fidelity FE meshes to represent complex parts. It has also saved time and effort in the subsequent analysis.

Power Systems↗