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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Full 3D visualization tool-kit for Monte Carlo and deterministic transport codes

We propose a package of tools capable of translating the geometric inputs and outputs of many Monte Carlo and deterministic radiation transport codes into open source file formats. These tools are aimed at bridging the gap between trusted, widely-used radiation analysis codes and very powerful, more recent and commonly used visualization software, thus supporting the design process and helping with shielding optimization. Three main lines of development were followed: mesh-based analysis of Monte Carlo codes, mesh-based analysis of deterministic codes and Monte Carlo surface meshing. The developed kit is considered a powerful and cost-effective tool in the computer-aided design for radiation transport code users of the nuclear world, and in particular in the fields of core design and radiation analysis. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1987

The NESSUS finite element code

The objective of this development is to provide a new analysis tool which integrates the structural modeling versatility of a modern finite element code with the latest advances in the area of probabilistic modeling and structural reliability. Version 2.0 of the NESSUS finite element code was released last February, and is currently being exercised on a set of problems which are representative of typical Space Shuttle Main Engine (SSME) applications. NESSUS 2.0 allows linear elastostatic and eigenvalue analysis of structures with uncertain geometry, material properties and boundary conditions, which are subjected to a random mechanical and thermal loading environment. The NESSUS finite element code is a key component in a broader software system consisting of five major modules. NESSUS/EXPERT is an expert system under development at Southwest Research Institute, with the objective of centralizing all component-specific knowledge useful for conducting probabilistic analysis of typical Space Shuttle Main Engine (SSME) components. NESSUS/FEM contains the finite element code used for the structural analysis and parameter sensitivity evaluation of these components. The task of parametrizing a finite element mesh in terms of the random variables present is facilitated with the use of the probabilistic data preprocessor in NESSUS/PRE. An external database file is used for managing the bulk of the data generated by NESSUS/FEM.

Dias, J. B.↗
DOE OSTI2022

A phase-field model for quasi-dynamic nucleation, growth, and propagation of rate-and-state faults

Despite its critical role in the study of earthquake processes, numerical simulation of the entire stages of fault rupture remains a formidable task. The main challenges in simulating a fault rupture process include the complex evolution of fault geometry, frictional contact, and off-fault damage over a wide range of spatial and temporal scales. Here, in this paper, we develop a phase-field model for quasi-dynamic fault nucleation, growth, and propagation, which features two standout advantages: (i) it does not require any sophisticated algorithms to represent fault geometry and its evolution; and (ii) it allows for modeling fault nucleation, propagation, and off-fault damage processes with a single formulation. Built on a recently developed phase-field framework for shear fractures with frictional contact, the proposed formulation incorporates rate- and state-dependent friction, radiation damping, and their impacts on fault mechanics and off-fault damage. We show that the numerical results of the phase-field model are consistent with those obtained from well-verified approaches that model the fault as a surface of discontinuity, without suffering from the mesh convergence issue in the existing continuous approaches to fault rupture (e.g., the stress glut method). Further, through numerical examples of fault propagation in various settings, we demonstrate that the phase-field approach may open new opportunities for investigating complex earthquake processes that have remained overly challenging for the existing numerical methods.

58 GEOSCIENCES↗
DOE OSTI2012

Development Status of PROTEUS-MOC

The PROTEUS-MOC solver is intended to resolve the outstanding computational challenges that the PROTEUS-SN is limited by and bridge the gap between PROTEUS-SN and existing homogenized structured grid solvers. The primary push behind all of this work is to carry out coupled thermal-mechanical, fluid dynamics, and neutronics physics problems. The PROTEUS-MOC development has taken advantage of the combined experience with PROTEUS-SN and DeCART with the goal of producing a tool which combines the best features of both tools. Taking the DeCART tool as the model, it merges a subgroup cross section interpolation technique with a 2D method of characteristics approach and a 1D SN formulation. DeCART has been specifically demonstrated to be accurate on PWR and VHTR problems, but it has displayed several issues with regard to accuracy on thin axial mesh sizes which are attributable to the 2D1D scheme it is built upon. There are also notable problems with performance and scalability all of which we hope to resolve in PROTEUS-MOC. PROTEUS-MOC is at its core, based upon the success achieved with MOCFE which demonstrated excellent parallelization with respect to angle and good parallelization with respect to space. More recent efforts have indicated parallelization with respect to energy is also feasible. The work completed yielded a novel methodology where the axial basis functions are treated separately from the radial ones, but the geometry is extruded. This eliminates the need to carry out and store 3D ray tracing data and simultaneously allows the coarse axial node sizes typical with DeCART. To rapidly build a prototype, a standard power iteration algorithm was constructed using an unpreconditioned GMRES solution algorithm for solving the flux vector in each fission source iteration. The new PROTEUS-MOC solver was deployed on two test problems: Takeda and a heterogeneous cylindrical heavy water system. The Takeda benchmark problem is a standard homogeneous test problem which was previously found to be too difficult for the 3D MOCFE code to solve. Using the PROTEUS-MOC we were able to show space-angle convergence although full convergence requires further improvements in performance. The second benchmark problem was constructed more for future work than this year, but the results generated are promising indicating that it will be very useful in highlighting the weaknesses and strengths of the PROTEUS-MOC tool. While we are far from a production tool, the results included in this report give us confidence that we are on the right path to constructing an interceding code between the plethora of structured grid homogenized assembly solvers and the supercomputer consuming PROTEUS-SN code. With time, the creation of a suitable preconditioner along with typical optimization should reduce the performance shortcomings that are observable today.

Marin-Lafleche, A.↗
DOE OSTI2022

June 2022 Permeability Measurements in P-Tunnel, Nevada National Security Site

This document describes gas permeability measurements made in June 2022 in P-Tunnel at the Nevada National Security Site (NNSS). Four boreholes (AC-1, GI-5, GI-6, and HF-1) were tested with a dual packer system. The pressure versus flow rate data that were collected are analyzed using a numerical simulation method that captures compressible gas flow in partially saturated porous media using a close approximation to the geometry of the test setup. The first part of this document begins with data collection and a description of the packer assembly followed by details on data acquisition, measurement instrumentation, and equipment followed by a description of pre-test system verification. Next, the implementation of testing in the field is described, including how specific locations were chosen for testing. The second part of this document describes the numerical simulations used to back-out permeability from the pressure versus flow rate data. First, the numerical model is described, including assumptions. Next the details of the mesh that represents the experiment are described. A method for generating type curves in pressure versus flow rate space is presented and finally the data are mapped to permeability. We conclude with a discussion of the results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A novel approach for tetrahedral-element-based finite element simulations of anisotropic hyperelastic intervertebral disc behavior

Intervertebral discs are microstructurally complex spinal tissues that add greatly to the flexibility and mechanical strength of the human spine. Attempting to provide an adjustable basis for capturing a wide range of mechanical characteristics and to better address known challenges of numerical modeling of the disc, we present a robust finite-element-based model formulation for spinal segments in a hyperelastic framework using tetrahedral elements. We evaluate the model stability and accuracy using numerical simulations, with particular attention to the degenerated intervertebral discs and their likely skewed and narrowed geometry. To this end, 1) annulus fibrosus is modeled as a fiber-reinforced Mooney-Rivlin type solid for numerical analysis. 2) An adaptive state-variable dependent explicit time step is proposed and utilized here as a computationally efficient alternative to theoretical estimates. 3) Tetrahedral-element-based FE models for spinal segments under various loading conditions are evaluated for their use in robust numerical simulations. For flexion, extension, lateral bending, and axial rotation load cases, numerical simulations reveal that a suitable framework based on tetrahedral elements can provide greater stability and flexibility concerning geometrical meshing over commonly employed hexahedral-element-based ones for representation and study of spinal segments in various stages of degeneration.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1995

Slot spiral antenna modeling using hybrid/mixed FE-BI technique

This report is concerned with the numerical simulation of the printed slot spiral antenna (SSA) using the finite element - boundary integral (FE-BI) analysis. It has been reported that the FE-BI technique is suited for modeling microstrip patch antennas of any shape, printed on layered planar structure or cylindrical platform and fed with a coaxial cable or a microstripline network underneath the radiating elements. However, direct application of the hybrid FE-BI technique to thin slot spirals requires excessive sampling rates to accurately simulate the geometry. To alleviate the meshing/modeling difficulties encountered with cavity-backed slot antennas, we describe a mixed finite element-boundary integral formulation. As in the past, the boundary integral is used to describe the radiation of the slot above the cavity. However the cavity is now modeled using a suitable mix of edge and node elements. The latter are used only at the aperture of the thin slot so that the nodes follow the center line of the slot. In this manner, regular size elements can be used regardless of the slot's width and any meshing restrictions are substantially relaxed. The proposed mixed element FE-BI formulation introduces three different computational regions and as expected this complicates the generation of the discrete system.

Gong, Jian↗
NASA NTRS2005

Face Gear Drive With Helical Involute Pinion: Geometry, Generation by a Shaper and a Worm, Avoidance of Singularities and Stress Analysis

A new type of face-gear drive with intersected axes of rotation formed by a helical involute pinion and conjugated face-gear has been investigated. Generation of face-gears by a shaper free of undercutting and pointing has been investigated. A new method of grinding or cutting of face-gears by a worm of special shape has been developed. A computerized design procedure has been developed to avoid undercutting and pointing by a shaper or by a generating worm. Also, a method to determine the limitations of the helix angle magnitude has been developed. The method provides a localization of the bearing contact to reduce the shift of bearing contact caused by misalignment. The analytical method provides a simulation of the meshing and contact of misaligned gear drives. An automatic mesh generation method has been developed and used to conduct a 3D contact stress analysis of several teeth. The theory developed is illustrated with several examples.

Litvin, Faydor L.↗
DOE OSTI2015

Mechanical behavior simulation of MEMS-based cantilever beam using COMSOL multiphysics

This paper presents the studies of mechanical behavior of MEMS cantilever beam made of poly-silicon material, using the coupling of three application modes (plane strain, electrostatics and the moving mesh) of COMSOL Multi-physics software. The cantilevers playing a key role in Micro Electro-Mechanical Systems (MEMS) devices (switches, resonators, etc) working under potential shock. This is why they require actuation under predetermined conditions, such as electrostatic force or inertial force. In this paper, we present mechanical behavior of a cantilever actuated by an electrostatic force. In addition to the simplification of calculations, the weight of the cantilever was not taken into account. Different parameters like beam displacement, electrostatics force and stress over the beam have been calculated by finite element method after having defining the geometry, the material of the cantilever model (fixed at one of ends but is free to move otherwise) and his operational space.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

SU-E-T-82: Comparison of Several Lumbar Intervertebral Fusion Titanium Cages with Respect to Their Backscattering Properties

Purpose: Investigating the backscatter dose factor with regards to structure and geometry of the surface material. Methods: The titanium cages used for this study representing both prototypes and well established products are made of a laser-sintered titanium alloy (AditusV GmbH, Berlin, Germany). A set of four radiochromic EBT3 films was used in a stacked geometry to measure the range and the magnitude of the expected surface dose enhancement due to the in comparison to water increased secondary electron release from the material. The measurement geometry and the small thickness of radiochromic EBT3 film allowed the dose measurement at distances of 0.1 mm, 0.9 mm, 1.7 mm and 2.5 mm from the probe surfaces. Water reference measurements were taken under equal conditions, in order to allow the calculation of the relative dose enhancement at the surface of a probe. Measurements were performed within a water phantom. An Epson Expression 10000 XL flatbed scanner was used for digitization. Results: Sintered titanium showed a dose enhancement factor of 1.22 at the surface of the material. The factor can be reduced to less than 1.10 by utilizing mesh structures. In both cases, the dose enhancement factor decreased to less than 1.03 at a distance of 1.7mm indicating the low energy of scattered electrons. Conclusion: Backscattering of titanium cages should be considered in treatment planning, especially when the cages are located close to organs at risk. While mesh structures were introduced to improve bone fusion with the implant structure, the potentially harmful surface dose enhancement is significantly reduced.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1993

Micromechanics of crenulated fibers in carbon/carbon composites

The influence of crenulated noncircular fibers on the micromechanical stress states due to a transverse strain and to a temperature change in carbon/carbon composites is examined using the finite element method. Stresses at the interface of both fully bonded and fully disbonded fibers having two crenulation amplitudes and with two fiber volume fractions are presented. In each case, these interface stresses are compared to stresses at the interface of circular fibers which have the same degree of disbond and fiber volume fraction and are under the same loading conditions. For the disbonded cases, deformed meshes showing locations of fiber/matrix contact are also included. In addition to the interface stress states, selected composite properties are also computed and compared in each case examined. Interest in studying noncircular fibers stems from a desire to increase the transverse properties of carbon/carbon by introducing a mechanical interlocking between the fiber and the matrix. Results presented here indicate that this interlocking does in fact occur. Evidence from the interface stress data suggests, however, that any possible advantage of this interlocking may be outweighed by the disadvantage of stress concentrations which arise at the interface due to the crenulated geometry of the fibers.

Carapella, E. E.↗
DOE OSTI2016

A 4DCT imaging-based breathing lung model with relative hysteresis

To reproduce realistic airway motion and airflow, the authors developed a deforming lung computational fluid dynamics (CFD) model based on four-dimensional (4D, space and time) dynamic computed tomography (CT) images. A total of 13 time points within controlled tidal volume respiration were used to account for realistic and irregular lung motion in human volunteers. Because of the irregular motion of 4DCT-based airways, we identified an optimal interpolation method for airway surface deformation during respiration, and implemented a computational solid mechanics-based moving mesh algorithm to produce smooth deforming airway mesh. In addition, we developed physiologically realistic airflow boundary conditions for both models based on multiple images and a single image. Furthermore, we examined simplified models based on one or two dynamic or static images. By comparing these simplified models with the model based on 13 dynamic images, we investigated the effects of relative hysteresis of lung structure with respect to lung volume, lung deformation, and imaging methods, i.e., dynamic vs. static scans, on CFD-predicted pressure drop. The effect of imaging method on pressure drop was 24 percentage points due to the differences in airflow distribution and airway geometry. - Highlights: • We developed a breathing human lung CFD model based on 4D-dynamic CT images. • The 4DCT-based breathing lung model is able to capture lung relative hysteresis. • A new boundary condition for lung model based on one static CT image was proposed. • The difference between lung models based on 4D and static CT images was quantified.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Gaussian integral method for void fraction

Here, a novel method, the Gaussian Integral Method (GIM), is presented for calculating void fractions in Computational Fluid Dynamics–Discrete Element Method (CFD-DEM) simulations. GIM is versatile and applicable to various grid types, including structured and unstructured polyhedral meshes, without requiring special boundary treatments. An optimization technique is introduced to make GIM independent of grid resolution and type. The method is validated against experimental data from a fluidized bed, demonstrating that GIM produces realistic simulations closely resembling experimental observations. Additionally, unstructured polyhedral grids using GIM outperform structured grids of equivalent resolution, yielding results more aligned with experimental data. The gradient of the void fraction is computed in the CFD solver and utilized in the DEM solver for precise estimation at particle locations. Overall, GIM provides an effective solution for void fraction calculations in particulate media simulations with complex geometries, enhancing the accuracy and applicability of CFD-DEM simulations for industrial processes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1992

Improved inhomogeneous finite elements for fabric reinforced composite mechanics analysis

There is a need to do routine stress/failure analysis of fabric reinforced composite microstructures to provide additional confidence in critical applications and guide materials development. Conventional methods of 3-D stress analysis are time consuming to set up, run and interpret. A need exists for simpler methods of modeling these structures and analyzing the models. The principal difficulty is the discrete element mesh generation problem. Inhomogeneous finite elements are worth investigating for application to these problems because they eliminate the mesh generation problem. However, there are penalties associated with these elements. Their convergence rates can be slow compared to homogeneous elements. Also, there is no accepted method for obtaining detailed stresses in the constituent materials of each element. This paper shows that the convergence rate can be significantly improved by a simple device which substitutes homogeneous elements for the inhomogeneous ones. The device is shown to work well in simple one and two dimensional problems. However, demonstration of the application to more complex two and three dimensional problems remains to be done. Work is also progressing toward more realistic fabric microstructural geometries.

Foye, R. L.↗
NASA NTRS1986

Transition and mixing in axisymmetric jets and vortex rings

A class of impulsively started, axisymmetric, laminar jets produced by a time dependent joint source of momentum are considered. These jets are different flows, each initially at rest in an unbounded fluid. The study is conducted at three levels of detail. First, a generalized set of analytic creeping flow solutions are derived with a method of flow classification. Second, from this set, three specific creeping flow solutions are studied in detail: the vortex ring, the round jet, and the ramp jet. This study involves derivation of vorticity, stream function, entrainment diagrams, and evolution of time lines through computer animation. From entrainment diagrams, critical points are derived and analyzed. The flow geometry is dictated by the properties and location of critical points which undergo bifurcation and topological transformation (a form of transition) with changing Reynolds number. Transition Reynolds numbers were calculated. A state space trajectory was derived describing the topological behavior of these critical points. This state space derivation yielded three states of motion which are universal for all axisymmetric jets. Third, the axisymmetric round jet is solved numerically using the unsteady laminar Navier Stokes equations. These equations were shown to be self similar for the round jet. Numerical calculations were performed up to a Reynolds number of 30 for a 60x60 point mesh. Animations generated from numerical solution showed each of the three states of motion for the round jet, including the Re = 30 case.

Allen, G. A., Jr.↗
NASA NTRS1994

Effects of a Forward-swept Front Rotor on the Flowfield of a Counterrotation Propeller

The effects of a forward-swept front rotor on the flowfield of a counterrotation model propeller at takeoff conditions at zero degree angle of attack are studied by solving the unsteady three-dimensional Euler equations. The configuration considered is an uneven blade count counterrotation model with twelve forward-swept blades on the fore rotor and ten aft-swept blades on the aft rotor. The flowfield is compared with that of a reference aft-swept counterrotation geometry and Laser Doppler Velocimeter (LDV) measurements. At the operating conditions considered, the forward-swept blade experiences a higher tip loading and produces a stronger tip vortex compared to the aft-swept blade, consistent with the LDV and acoustic measurements. Neither the solution nor the LDV data indicated the formation of a leading edge vortex. The predicted radial distribution of the circumferentially averaged axial velocity at the measurement station agreed very closely with LDV data, while crossflow velocities showed poor agreement. The discrepancy between prediction and LDV data of tangential and radial velocities is due in part to the insufficient mesh resolution in the region between the rotors and in the tip region to track the tip vortex. The vortex is diffused by the time it arrives at the measurement station. The uneven blade count configuration requires the solution to be carried out for six blade passages of the fore rotor and five passages of the aft rotor, thus making grid refinement prohibitive.

Nallasamy, M.↗
NASA NTRS2009

Discrete Adjoint-Based Design Optimization of Unsteady Turbulent Flows on Dynamic Unstructured Grids

An adjoint-based methodology for design optimization of unsteady turbulent flows on dynamic unstructured grids is described. The implementation relies on an existing unsteady three-dimensional unstructured grid solver capable of dynamic mesh simulations and discrete adjoint capabilities previously developed for steady flows. The discrete equations for the primal and adjoint systems are presented for the backward-difference family of time-integration schemes on both static and dynamic grids. The consistency of sensitivity derivatives is established via comparisons with complex-variable computations. The current work is believed to be the first verified implementation of an adjoint-based optimization methodology for the true time-dependent formulation of the Navier-Stokes equations in a practical computational code. Large-scale shape optimizations are demonstrated for turbulent flows over a tiltrotor geometry and a simulated aeroelastic motion of a fighter jet.

Nielsen, Eric J.↗
NASA NTRS1998

Effects of regenerator geometry on pulse tube refrigerator performance

This paper gives results of the cooling performance of a double-inlet pulse tube refrigerator using various regenerators. The same pulse tube was used for all the experiments and measured 4.76 mm in diameter and 46.2 mm in length. A commercial linear compressor with a swept volume of 4 cm3 was used in these experiments. The operating conditions were held constant at a mean pressure of 2.0 MPa and a frequency of 54 Hz. Using finite difference software called REGEN3.1, developed at NIST, and recent experiment results, we optimized a series of regenerators based on dimensions, materials and screen packing. The values used for calculating the thermal conduction through stacked screens by REGEN3.1 were based on recent experimental results from NIST. The regenerator tubes were designed using 316 stainless steel and titanium materials. The regenerator matrices investigated were 400-mesh and 500-mesh stainless steel screen. The valve settings for both orifices were adjusted to minimize the no-load temperature for all regenerators. A cooling capacity curve from 0 to 3 W was also determined. The performance of the pulse tube refrigerator using the different regenerators is discussed. The experimental results from the various regenerators are evaluated and compared with their corresponding numerically calculated coefficient of performance (COP) and regenerator design as determined by REGEN3.1.

NASA Discipline Life Sciences Technologies↗