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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Construction of Hexahedral Elements Mesh Capturing Realistic Geometries of Bayou Choctaw SPR Site

The three-dimensional finite element mesh capturing realistic geometries of Bayou Choctaw site has been constructed using the sonar and seismic survey data obtained from the field. The mesh is consisting of hexahedral elements because the salt constitutive model is coded using hexahedral elements. Various ideas and techniques to construct finite element mesh capturing artificially and naturally formed geometries are provided. The techniques to reduce the number of elements as much as possible to save on computer run time with maintaining the computational accuracy is also introduced. The steps and methodologies could be applied to construct the meshes of Big Hill, Bryan Mound, and West Hackberry strategic petroleum reserve sites. The methodology could be applied to the complicated shape masses for not only various civil and geological structures but also biological applications such as artificial limbs.

42 ENGINEERING↗
NASA NTRSDate not supplied

HLPW-4/GMGW-3: Mesh Adaptation for RANS Technology Focus Group Workshop Summary

The AIAA workshop series helps guide the Computational Fluid Dynamics (CFD) community through publicly available test cases and impartial evaluation to establish the state-of-the-art in applied CFD. A summary of the Mesh Adaptation Technical Focus Group submissions for the joint Fourth CFD High Lift Prediction and Third Geometry & Mesh Generation Workshop is provided. The Common Research Model High-Lift variant (CRM-HL) is the target of this investigation. A 2D high-lift airfoil solution verification case was extracted from the 3D CRM-HL. Low variation is shown between 2D solutions with mesh adaptation to control estimated solution interpolation and output errors. As expected, variation for these mesh adaptation techniques is larger for the 3D CRM-HL than the 2D airfoil. More importantly, the 3D mesh adapted results have a lower variation than expert-crafted meshes of equal or larger size. Suspected multiple solutions to Reynolds-averaged Navier-Stokes with the Spalart-Allmaras turbulence model appear to resolve into a common solution with mesh resolution.

HLPW-4↗
DOE OSTI2020

Spectral Element Applications in Complex Nuclear Reactor Geometries: Tet-to-Hex Meshing

The spectral element code Nek5000 is an open-source, higher-order computational fluid dynamics code developed at Argonne National Laboratory. It is designed to solve incompressible Navier-Stokes equations, but it also has a low-Mach-number approximation feature available. Large eddy simulation is approached by explicit filtering of the velocity field (and other fields) to mimic the effect of dissipation due to the unresolved scale. The computational domain is decomposed into second-order hexahedral elements that conform to the boundaries. However, generating a high-quality pure-hexahedral mesh can be challenging for some problems. For simple geometries, traditional blocking methods can be used to decompose the domain into smaller blocks to generate a so-called structural mesh. A structural mesh can maintain good orthogonality but can have a highly skewed mesh to conform to the geometry, as well as unnecessary refinement in the far field. Moreover, for geometries with relative complexity, blocking the geometry becomes impossible. To address these issues, we adopted a let-to-hex strategy to generate a pure hexahedral mesh for Nek5000. First, we generate a pure tetrahedral mesh for the geometry; then we divide one tetrahedral element into four hexahedral elements. A pure tetrahedral mesh could be easily generated for complex geometries by using many current meshing codes. In this paper, we use the commercial codes ANSYS meshing and ANSYSICEM to generate the pure tetrahedral mesh and then convert it to a pure hexahedral mesh. Boundary layers are extruded in ANSYSICEM to maintain near-wall resolution.

ANSYS ICEM↗
DOE OSTI2025

Advancing the STS Neutron Moderator Design with an Automated Optimization Workflow and Unstructured Mesh Modeling

With the Second Target Station approaching its final design phase, a detailed neutronics evaluation of its critical components is necessary. Optimizing the dimensions of the two cold-source moderators that are at the heart of this facility presents a multi-objective optimization problem for which an accurate geometric description is crucial. We have applied a fully automated optimization workflow in which a detailed unstructured mesh geometry is automatically generated with Attila4MC, starting from a parametrized CREO geometry followed by preprocessing with SpaceClaim. With this geometry, a MCNP run is performed to calculate the brightness metrics, which are subsequently provided to the optimization algorithm in DAKOTA that provides new parameters and drives the optimization loop until convergence. In this paper, we show the results of the analysis that are used for the final design of the cylindrical and tube moderator. The optimization simulations provide a refinement to and confirmation of the conclusions of the previous design iteration. Additional to the optimization, a sensitivity study is performed to study the effect of minor geometry changes, which is important for the final engineering design. In conclusion, with these studies, we demonstrate that the automated workflow and high-fidelity unstructured mesh modeling are efficient tools for a thorough design evaluation.

DAKOTA↗
NASA NTRSDate not supplied

Aerothermodynamic Analysis of a Flexible Thermal Protection System under Reentry Loads

The Carryall Block 1 reentry vehicle being developed by Outpost Space utilizes a strut supported semirigid deployable heatshield. This consists of a flexible thermal protection system, a heat-resistant fabric stack, stretched over actuated spars. The advantages of a deployable heatshield include reduced heat loading and earlier deceleration in the trajectory. However, the nature of the flexible thermal protection system necessitates considering the loaded shape of the heat shield. The flexible thermal protection system will deflect under reentry loads leading to areas of higher heating rates as well as a reduced axial coefficient as compared to the nominal shape. The Carryall Block 1 is analyzed using NASA’s FUN3D and DPLR CFD solvers with a deflected shape based on the catenary equations. The aerodynamic results are found to be within a percent for both solvers and both structured and unstructured mesh types. Fluid Structure Interaction (FSI) analysis is currently a work in progress, using file I/O to communicate between FUN3D and LS-DYNA, a commercial nonlinear structural solver. Challenges in deforming the geometry, mesh, and initial results are presented here.

thermal protection system↗
DOE OSTI2018

High Voltage Insulation and Gas Absorption of Polymers in High Pressure Argon and Xenon Gases

This work advances the development of the high-pressure xenon gas Time projection Chamber (TPC) concept in the search for neutrino-less double beta decay in Xenon136, a promising isotopic candidate. In particular, this work advances the use of electroluminescence (EL) as a nearly noiseless gain mechanism that preserves energy resolution in measurement of the visible energy in the decay—the two emitted electrons. Energy resolution of ~1% FWHM at the Q-value of the decay, 2458 kEV, allows complete separation of the much more copious 2-neutrino mode. The goal is to realize the EL mechanism in a very large parallel plane, nearly optically transparent mesh geometry, on the order of one meter diameter. Additional effort by collaborators at Argonne National Laboratory provide design solutions for the electrostatic boundary conditions of the active volume of the TPC. Solutions for the large mesh have been found in either electroformed grids or a more conventional woven mesh made from stainless steel wire. Posts to maintain approximate planarity are needed to minimize electrostatic deflection. The work has shown that these posts are feasible with certain plastics that can support the high electric field with out discharges.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Three-dimensional analysis using polygons in particle and heavy ion transport code system (PHITS) - 25504

Recently a function to treat tetrahedral-mesh geometry was introduced in Particle and Heavy Ion Transport code System (PHITS). Owing to this function, three-dimensional analysis using polygons became possible. The work flow of the three-dimensional analysis using a particular polygon object was demonstrated, converting the data into tetrahedral mesh with the tetrahedral-mesh generator TetGen, computing the radiation transport with PHITS, and finally analyzing with three-dimensional data analysis application ParaView. (author)

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2021

Serpent neutronics model of Wendelstein 7-X for 14.1 MeV neutrons

In this work, a Serpent 2 neutronics model of the Wendelstein 7-X (W7-X) stellarator is prepared, and an response function for the Scintillating-Fibre neutron detector (SciFi) is calculated using the model. The neutronics model includes the simplified geometry for the key components of the stellarator itself as well as the torus hall. The objective of the model is to assess the 14.1 MeV neutron flux from deuteron-triton fusions in W7-X, where the neutrons are modelled only until they have slowed down to 1 MeV energy. The key messages of this article are: demonstration of unstructured mesh geometry usage for stellarators, W7-X in particular; technical documentation of the model and first insights in fast neutron behaviour in W7-X, especially related to the SciFi: the model indicates that the superconducting coils are the strongest scatterers and block neutrons from large parts of the plasma. The back-scattering from e.g. massive steel support structures is found to be small. The SciFi will detect neutrons from an extended plasma volume in contrast to having an effective line-of-sight.

14 MeV↗
DOE OSTI2016

TU-H-CAMPUS-IeP2-03: Development of 3D Printed Coronary Phantoms for In-Vitro CT-FFR Validation Using Data from 320- Detector Row Coronary CT Angiography

Purpose: To validate Computed Tomography Fractional Flow Reserve (CT-FFR) measurements with accurate 3D printed coronary phantoms. Methods: DICOM data from four phases in two patients imaged with a standard 320 × 0.5mm coronary CT acquisition (70–80% cardiac cycle) underwent semi-automated segmentation using a research workstation. Both patients had a >50% stenosis from the clinical image interpretation. Each volume was saved as a Stereo Lithographic (STL) file with 250 micron resolution. The 3D geometries were qualitatively assessed; the best of the four phases was 3D printed using a Stratasys Eden260V printer in Tango+, a rubber-like material that roughly emulates mechanical properties of human vasculature. We connected the model to a programmable pump and measured the pressure drop using pressure sensors embedded proximal and distal to the arterial stenosis. Next, the STL files used for the 3D printed models were uploaded in the ANSYS meshing tool (ICEM CFD 16.1). A standard meshing process was applied and the meshed geometry was directly imported in the ANSYS Fluent for Computational Flow Dynamics simulations. The CFD simulations were used to calculate the CT-FFR and compared to the bench top FFR measured in the 3D printed phantoms. Results: FFR-CT measurements and phantoms were completed in within an hour after the segmentation. Patient 1 had a 60% stenosis that resulted in a CT-FFR of 0.68. The second case had a 50% stenosis and a CT-FFR of 0.75. The average bench top FFR measurements were 0.72 and 0.80, respectively. Conclusion: This pilot investigation demonstrated the use of a bench-top coronary model for CT-FFR validation. The measurements and the CFD simulations agreed within 6%. Project supported by Support: Toshiba America Medical Systems Corp.and NIH grant R01-EB002873. Project supported by Toshiba America Medical Systems Corp.and partial support from NIH grant R01-EB002873.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

EBS Task Force: Task 9/FEBEX Modeling Final Report: Thermo-Hydrological Modeling with PFLOTRAN

This report outlines Sandia National Laboratories modeling studies applied to Stage 1 and Stage 2 of the Full-scale Engineered Barriers Experiment in Crystalline Host Rock (FEBEX) in situ test for the SKB EBS Task Force Task 9. The FEBEX test was a full-scale test conducted over ~18 years at the Grimsel, Switzerland Underground Research Laboratory (URL) managed by NAGRA. It involved emplacing simulated waste packages, in the form of welded cylindrical heaters, inside a tunnel in crystalline granitic rock and surrounded by a bentonite barrier and cement plug. Sensors emplaced within the bentonite monitored the wetting-up, heating, and drying out of the bentonite barrier, and the large resulting data set provides an excellent opportunity for validation of multiphysics Thermal-Hydrological (TH), Thermal-Hydrologic-Chemical (THC), and Thermal-Hydrological-Mechanical (THM) modeling approaches for underground nuclear waste storage and the performance of engineered bentonite barriers. The present status of the EBS Task Force is finalizing Task 9, which follows years of modeling studies of the FEBEX test, by many notable modeling teams (Gens et al., 2009; Sanchez et al. 2010; 2012; Samper et al., 2018). These modeling studies generally use two-dimensional axisymmetric meshes, ignoring threedimensional effects, gravity and asymmetric wetting and dry out of the bentonite engineered barrier. This study investigates these effects with use of the PFLOTRAN THC code with massively parallel computational methods in modeling FEBEX Stage 1 and Stage 2 results. The PFLOTRAN numerical code is an open source, state-of-the-art, massively parallel subsurface flow and reactive transport code operating in a high-performance computing environment (Hammond et al., 2014). Section 2 describes the applied partial differential equations describing mass, momentum and energy balance used in this study, considerations derived by assuming phase equilibrium between gas and liquid phases, constitutive equations for granite, cement plug, and bentonite domains, and specific approaches for use inthe PFLOTRAN code. Section 3 describes the geometry, meshing, and model set-up. Section 4 describes modeling results, Section 5 compares modeling results to field testing data, and Section 6 gives conclusions. The Appendix provides detailed information required by the EBSTask Force for final reporting.

42 ENGINEERING↗
DOE OSTI2016

Framework and algorithms for illustrative visualizations of time-varying flows on unstructured meshes

Photo- and physically realistic techniques are often insufficient for visualization of fluid flow simulations, especially for 3D and time-varying studies. Substantial research effort has been dedicated to the development of non-photorealistic and illustration-inspired visualization techniques for compact and intuitive presentation of such complex datasets. However, a great deal of work has been reproduced in this field, as many research groups have developed specialized visualization software. Additionally, interoperability between illustrative visualization software is limited due to diverse processing and rendering architectures employed in different studies. In this investigation, a framework for illustrative visualization is proposed, and implemented in MarmotViz, a ParaView plug-in, enabling its use on a variety of computing platforms with various data file formats and mesh geometries. Region-of-interest identification and feature-tracking algorithms incorporated into this tool are described. Implementations of multiple illustrative effect algorithms are also presented to demonstrate the use and flexibility of this framework. Here, by providing an integrated framework for illustrative visualization of CFD data, MarmotViz can serve as a valuable asset for the interpretation of simulations of ever-growing scale.

42 ENGINEERING↗
DOE OSTI2025

Coupled beam-target-moderator optimization for the Second Target Station

To support the design of the Second Target Station, that aims to provide the world’s highest peak brightness of cold neutrons, studies that optimize the dimensions of the target and moderators are invaluable. In this work, we investigate the influence of the target dimensions and beam profile on the performance and optimal size of the moderators. We perform optimization runs with detailed MCNP6.2 simulations using high-fidelity unstructured mesh geometries generated from parametrized CAD models. We demonstrate that small changes in target height do not influence the moderator performance if the beam dimensions are chosen adequately. We quantify the effect of beam footprint and target width on the moderator performance and show that the optimal moderator dimensions are insensitive to limited changes in target and beam profile.

Dakota↗
DOE OSTI2025

Solving sparse finite element problems on neuromorphic hardware

The finite element method (FEM) is one of the most important and ubiquitous numerical methods for solving partial differential equations (PDEs) on computers for scientific and engineering discovery. Applying the FEM to larger and more detailed scientific models has driven advances in high-performance computing for decades. Here we demonstrate that scalable spiking neuromorphic hardware can directly implement the FEM by constructing a spiking neural network that solves the large, sparse, linear systems of equations at the core of the FEM. We show that for the Poisson equation, a fundamental PDE in science and engineering, our neural circuit achieves meaningful levels of numerical accuracy and close to ideal scaling on modern, inherently parallel and energy-efficient neuromorphic hardware, specifically Intel’s Loihi 2 neuromorphic platform. We illustrate extensions to irregular mesh geometries in both two and three dimensions as well as other PDEs such as linear elasticity. Our spiking neural network is constructed from a recurrent network model of the brain’s motor cortex and, in contrast to black-box deep artificial neural network-based methods for PDEs, directly translates the well-understood and trusted mathematics of the FEM to a natively spiking neuromorphic algorithm.

Applied mathematics↗
DOE OSTI2015

Soft, curved electrode systems capable of integration on the auricle as a persistent brain–computer interface

Some recent advances in electrodes for noninvasive recording of electroencephalograms expand opportunities collecting such data for diagnosis of neurological disorders and brain–computer interfaces. Existing technologies, but, cannot be used effectively in continuous, uninterrupted modes for more than a few days due to irritation and irreversible degradation in the electrical and mechanical properties of the skin interface. We introduce a soft, foldable collection of electrodes in open, fractal mesh geometries that can mount directly and chronically on the complex surface topology of the auricle and the mastoid, to provide high-fidelity and long-term capture of electroencephalograms in ways that avoid any significant thermal, electrical, or mechanical loading of the skin. Experimental and computational studies establish the fundamental aspects of the bending and stretching mechanics that enable this type of intimate integration on the highly irregular and textured surfaces of the auricle. Furthermore, cell level tests and thermal imaging studies establish the biocompatibility and wearability of such systems, with examples of high-quality measurements over periods of 2 wk with devices that remain mounted throughout daily activities including vigorous exercise, swimming, sleeping, and bathing. Demonstrations include a text speller with a steady-state visually evoked potential-based brain–computer interface and elicitation of an event-related potential (P300 wave).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

SEFOR Core I-E

This repository contains the VTB model developed for the Southwest Experimental Fast Oxide Reactor (SEFOR) core configuration I-E using the NEAMS tools. The MOOSE Reactor Module was employed to set up the geometry mesh, MC2-3 was applied to generate multigroup cross sections at various temperatures, and Griffin was used to calculate the k-eff values of Core I-E across these temperature conditions. To validate the deterministic results from Griffin, reference solutions were obtained with the Monte Carlo code Shift. The reactivity feedback was then derived and compared against SEFOR experimental data from isothermal tests in order to assess the performance of Griffin. Reproducing these results requires licensed access to several code systems, including Shift, Griffin, and MC2-3. Instructions for obtaining and accessing these codes can be found at the following locations: https://www.ornl.gov/scale/releases for SCALE/Shift, https://mooseframework.inl.gov/ncrc/applications/ncrc_root_griffin.html for Griffin, and https://www.anl.gov/nse/mc23 for MC2-3.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

Supo Thermal Model Development II

This report describes the continuation of the Computational Fluid Dynamics (CFD) model of the Supo cooling system described in the report, Supo Thermal Model Development1, by Cynthia Buechler. The goal for this report is to estimate the natural convection heat transfer coefficient (HTC) of the system using the CFD results and to compare those results to remaining past operational data. Also, the correlation for determining radiolytic gas bubble size is reevaluated using the larger simulation sample size. The background, solution vessel geometry, mesh, material properties, and boundary conditions are developed in the same manner as the previous report. Although, the material properties and boundary conditions are determined using the appropriate experiment results for each individual power level.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗