Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational geometry”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

National Combustion Code, a Multidisciplinary Combustor Design System, Will Be Transferred to the Commercial Sector

The NASA Lewis Research Center and Flow Parametrics will enter into an agreement to commercialize the National Combustion Code (NCC). This multidisciplinary combustor design system utilizes computer-aided design (CAD) tools for geometry creation, advanced mesh generators for creating solid model representations, a common framework for fluid flow and structural analyses, modern postprocessing tools, and parallel processing. This integrated system can facilitate and enhance various phases of the design and analysis process.

Steele, Gynelle C.↗
NASA NTRS2004

Overcoming Geometry-Induced Stiffness with IMplicit-Explicit (IMEX) Runge-Kutta Algorithms on Unstructured Grids with Applications to CEM, CFD, and CAA

My goal is to develop and implement efficient, accurate, and robust Implicit-Explicit Runge-Kutta (IMEX RK) methods [9] for overcoming geometry-induced stiffness with applications to computational electromagnetics (CEM), computational fluid dynamics (CFD) and computational aeroacoustics (CAA). IMEX algorithms solve the non-stiff portions of the domain using explicit methods, and isolate and solve the more expensive stiff portions using implicit methods. Current algorithms in CEM can only simulate purely harmonic (up to lOGHz plane wave) EM scattering by fighter aircraft, which are assumed to be pure metallic shells, and cannot handle the inclusion of coatings, penetration into and radiation out of the aircraft. Efficient MEX RK methods could potentially increase current CEM capabilities by 1-2 orders of magnitude, allowing scientists and engineers to attack more challenging and realistic problems.

Kanevsky, Alex↗
NASA NTRS2006

A Summary of the Experimental Results for a Generic Tractor-Trailer in the Ames Research Center 7- by 10-Foot and 12-Foot Wind Tunnels

Experimental measurements of a generic tractor-trailer were obtained in two wind tunnels at Ames Research Center. After a preliminary study at atmospheric conditions in the 7- by 10-Foot Wind Tunnel, additional testing was conducted at Reynolds numbers corresponding to full-scale highway speeds in the 12-Foot Pressure Wind Tunnel. To facilitate computational modeling, the 1:8-scale geometry, designated the Generic Conventional Model, included a simplified underbody and omitted many small-scale details. The measurements included overall and component forces and moments, static and dynamic surface pressures, and three-component particle image velocimetry. This summary report highlights the effects of numerous drag reduction concepts and provides details of the model installation in both wind tunnels. To provide a basis for comparison, the wind-averaged drag coefficient was tabulated for all configurations tested. Relative to the baseline configuration representative of a modern class-8 tractor-trailer, the most effective concepts were the trailer base flaps and trailer belly box providing a drag-coefficient reduction of 0.0855 and 0.0494, respectively. Trailer side skirts were less effective yielding a drag reduction of 0.0260. The database of this experimental effort is publicly available for further analysis.

Storms, Bruce L.↗
NASA NTRS2010

Propellant Sloshing Parameter Extraction from CFD Analysis

Propellant slosh is a potential source of disturbance critical to the stability of space vehicle. The sloshing dynamics is typically represented by a mechanical model of spring mass damper. This mechanical model is then included in the equation of motion of the entire vehicle for Guidance, Navigation and Control analysis. The typical parameters required by the mechanical model include natural frequency of the sloshing, sloshing mass, sloshing mass center coordinates, and critical damping coefficient. During the 1960 s US space program, these parameters were either computed from analytical solution for simple geometry or by experimental testing for the sub-scaled configurations. The purpose of this work is to demonstrate the soundness of a CFD approach in modeling the detailed fluid dynamics of tank sloshing and the excellent accuracy in extracting mechanical properties for different tank configurations and at different fill levels. The validation studies included straight cylinder against analytical solution, and sub-scaled Centaur LOX and LH2 tanks with and without baffles against experimental results. This effort shows that CFD technology can provide accurate mechanical parameters for any tank configuration, and is especially valuable to the future design of propellant tanks, as there is no previous experimental data available for the same size and configuration.

Yang, H. Q.↗
NASA NTRS2011

Linear Approximation to Optimal Control Allocation for Rocket Nozzles with Elliptical Constraints

In this paper we present a straightforward technique for assessing and realizing the maximum control moment effectiveness for a launch vehicle with multiple constrained rocket nozzles, where elliptical deflection limits in gimbal axes are expressed as an ensemble of independent quadratic constraints. A direct method of determining an approximating ellipsoid that inscribes the set of attainable angular accelerations is derived. In the case of a parameterized linear generalized inverse, the geometry of the attainable set is computationally expensive to obtain but can be approximated to a high degree of accuracy with the proposed method. A linear inverse can then be optimized to maximize the volume of the true attainable set by maximizing the volume of the approximating ellipsoid. The use of a linear inverse does not preclude the use of linear methods for stability analysis and control design, preferred in practice for assessing the stability characteristics of the inertial and servoelastic coupling appearing in large boosters. The present techniques are demonstrated via application to the control allocation scheme for a concept heavy-lift launch vehicle.

Orr, Jeb S.↗
NASA NTRS2013

Advances in SPICE Support of Planetary Science

SPICE is the de facto international standard for determining the geometric conditions-parameters such as altitude, lighting angles, and LAT/LON coverage of an instrument footprint-pertaining to scientific observations acquired by instruments on board robotic spacecraft. This system, comprised of data and allied software, is used for planning science observations and for analyzing the data returned from those observations. Use of SPICE is not a NASA requirement but is recommended by NASA's Planetary Data System and by the International Planetary Data Alliance. Owing in part to its reliability, stability, portability and user support, the use of SPICE has spread to many national space agencies, including those of the U.S., Europe (ESA), Japan, Russia and India. SPICE has been in use since the Magellan mission to Venus and so has many well-known capabilities. But the NAIF Team responsible for implementing SPICE continues to add new features; this presentation describes a number of these.

computational capability↗
NASA NTRS2022

Multi-Fidelity CFD Simulation of the Multirotor Test Bed

The high-fidelity and reduced order CFD method's validation and analysis are presented for the Multirotor Test Bed. The reduced-order method uses a rotor-disk model in the OVERFLOW CFD code and is an important aspect of this study. A rotor-trim option loosely coupling the Comprehensive code CAMRAD II with the OVERFLOW CFD code is included in the simulation that captures rotor blade aerodynamics, blade and rotor performance, rotor-rotor, and rotor-body interactions. Detailed comparison and validation were performed for both higher-fidelity full rotor unsteady flows and reduced-order rotor-disk methods. The tunnel geometry is included in the detailed computational effort. Computations of the low-Reynolds number flows were performed for a different number of in-plane rotors as well as for a single rotor modeled inside the tunnel.

ARMD↗
NASA NTRSDate not supplied

Computational Analysis of a Quiet Single-Main Rotor Helicopter for Air Taxi Operations

A computational study has been done on NASA’s quiet single-main rotor helicopter concept for urban air mobility using high-fidelity computational fluid dynamics, rotorcraft comprehensive analysis tools, and computational aeroacoustics. High-order accurate schemes, dual-time stepping, and the delayed detached-eddy simulation model have been employed. A loose-coupling approach between the flow solver and the rotorcraft comprehensive code is implemented to include vehicle trim and blade motions. Acoustic simulations based on the Ffowcs-Williams and Hawkings equations have been performed to compute the rotor noise. Different blade geometries and tip speeds are analyzed, and performance and acoustics results are compared. The vehicle has been simulated in hover and cruise, with flight conditions representative of an air taxi mission. The quiet single-main rotor helicopter is one of the conceptual designs intended to focus and guide NASA’s research activities in support of aircraft development for vertical take-off and landing air taxi operations.

ARMD↗
DOE OSTI2022

From PeleC to PeleACC, to PeleC++

PeleC is an Exascale Computing Project application for simulating compressible combustion in complex geometries. It has been built on top of the popular AMReX library. In the beginning of the Exascale Computing Project, PeleC was focused on KNL. It uses a mixture of C++, C, and kernels written in Fortran to obtain performance by focusing on vectorization. Recently we have taken two approaches in deciding PeleC's future for obtaining performance on exascale GPU machines. In the first programming model, we decorated the Fortran kernels with OpenACC directives. This expedited our ability to run at large scales on Summit's GPUs, where we achieved a significant speedup over the CPUs on Summit. The second programming model involved rewriting the Fortran kernels in C++ and using AMReX's Kokkos-like lambda abstractions for running on the GPU. This resulted in similar speedups on Summit's GPUs over merely utilizing the CPUs. Both approaches involved AMReX's management of memory transfers between the device and host. In this work, we compare and contrast the benefits and pitfalls to both programming approaches regarding performance, performance portability, and productivity. We also discuss advantages we have found in taking the time to modernize our code and why have chosen a specific pathway to prepare our code for the future DOE exascale machines.

exascale computing↗
DOE OSTI2022

Development of Numerical Model of Metal Foam with PCM for the Estimation of Effective Thermal Conductivity

Global warming due to climate change is a threat to humankind. Nuclear energy is one of the promising solutions to reduce fossil fuel usage. Nuclear energy can handle the base load, compensating for the volatility of renewable energy. If nuclear energy could achieve load following capability, the combination with renewable energy would be more suitable. Thermal energy storage (TES) is one of the options for enabling load following of nuclear reactors. The TES makes it possible to store surplus nuclear thermal energy and release it later as needed. In Idaho National Laboratory (INL), a new concept of latent heat TES integrated with high-temperature heat pipe has been proposed and is under development, which is called Heat pipe-Integrated Thermal Battery (HITB). HITB exchanges thermal energy between the reactor system and TES via heat pipe. The heat transferred to TES medium, made of phase change material (PCM), stores energy as sensible heat and/or latent heat. As PCM typically has poor thermal conductivity, however, various heat transfer enhancement techniques are required to achieve a rapid charging cycle. There are many techniques to enhance the heat transfer ability of TES medium such as disk, fin, and metal foam. Among them, metal foam is an appropriate option to enhance the heat transfer because it maximizes the heat transfer area through metal wicks. Metal foam is a lightweight metal structure that has a high porosity of over 0.9. The typical materials for metal foam are Aluminum, Copper, Nickel, and Silicon Carbide (SiC). Metal foam not only enhances heat transfer via conduction but also increases contact surface area. In the HITB design , the metal foam is being considered as one of the options to enhance the heat transfer of TES medium (PCM) [1]. To predict the enhanced thermal performance of TES, one should properly estimate the effective thermal conductivity of metal foam combined with PCM material or calculate heat transfer in distributed model. There are many experimental works that provides effective thermal conductivity of metal foam with various PCM [2,3]. Also, many theoretical models were developed based on the unit cell model of metal foam [4,5]. With a distributed model, on the other hand, detail heat transfer characteristics between metal foam and PCM material can be analyzed considering the geometry or buoyancy effect. However, due to the complex geometry of metal foam pores, the computational cost for three-dimensional modeling highly increases. Therefore, if metal foam structure can be modeled in simple and repetitive design, the computational cost would decrease Among the various metal foam models [2], lattice model is one of the simple and extendable design. The porosity and pores per inch (PPI) can be characterized by the size and spatial distance of lattice structure. If the three-dimensional metal foam model consists of lattice structure could properly estimate the heat transfer, which is characterized by effective thermal conductivity, it would be a good option to assess the thermal performance of metal foam with PCM. In this study, a three-dimensional numerical model was developed to simulate conductive heat transfer between metal foam and PCM. The three-dimensional lattice structure of square pillars was selected as a basic structure of the metal foam. The calculation result was characterized by the effective thermal conductivity of the whole domain. A sensitivity study was conducted for mesh size, domain size, and PPI to check whether the calculation result gives a converged result or not. Lastly, the effective thermal conductivity from the lattice model was compared with existing experimental data to validate the model result

25 ENERGY STORAGE↗
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↗
DOE OSTI2023

Metamodels for Rapid Analysis of Large Sets of Building Designs for Robotic Constructability: Technology Demonstration Using the NASA 3D Printed Mars Habitat Challenge

Disruptive robotic construction technologies such as additive deposition of cementitious materials like concrete (or "3D concrete printing") require the synchronous operation of multiple pieces of equipment in the production setup. In such an environment, it is crucial to simulate the robotic motions (for toolpath clashes) and the cementitious material behavior (for toolpath failures) to ensure fail-proof constructability of the envisioned building geometry. However, toolpath clash detection requires 4D simulations of the production setup, which are computationally graphics intensive, whereas toolpath failure detection requires actual 3D printing of test parts from the geometry to identify areas prone to failure while 3D printing, which is physically tedious. Both these processes, being computationally and physically intensive, have largely curtailed designers from simulating and exploring large sets of design options with varying geometries and toolpath configurations. To overcome this and allow designers to explore large sets of design possibilities, this paper proposes two novel computational metamodels capable of performing robotic toolpath clash detection and failure detection with significantly reduced times than the earlier approaches. The developed metamodels were used to rapidly simulate large sets of building design options for robotic constructability in the NASA 3D-Printed Mars Habitat Challenge.

clash detection↗
DOE OSTI2015

A generalized reconstruction framework for unconventional PET systems

Purpose: Quantitative estimation of the radionuclide activity concentration in positron emission tomography (PET) requires precise modeling of PET physics. The authors are focused on designing unconventional PET geometries for specific applications. This work reports the creation of a generalized reconstruction framework, capable of reconstructing tomographic PET data for systems that use right cuboidal detector elements positioned at arbitrary geometry using a regular Cartesian grid of image voxels. Methods: The authors report on a variety of design choices and optimization for the creation of the generalized framework. The image reconstruction algorithm is maximum likelihood-expectation–maximization. System geometry can be specified using a simple script. Given the geometry, a symmetry seeking algorithm finds existing symmetry in the geometry with respect to the image grid to improve the memory usage/speed. Normalization is approached from a geometry independent perspective. The system matrix is computed using the Siddon’s algorithm and subcrystal approach. The program is parallelized through open multiprocessing and message passing interface libraries. A wide variety of systems can be modeled using the framework. This is made possible by modeling the underlying physics and data correction, while generalizing the geometry dependent features. Results: Application of the framework for three novel PET systems, each designed for a specific application, is presented to demonstrate the robustness of the framework in modeling PET systems of unconventional geometry. Three PET systems of unconventional geometry are studied. (1) Virtual-pinhole half-ring insert integrated into Biograph-40: although the insert device improves image quality over conventional whole-body scanner, the image quality varies depending on the position of the insert and the object. (2) Virtual-pinhole flat-panel insert integrated into Biograph-40: preliminary results from an investigation into a modular flat-panel insert are presented. (3) Plant PET system: a reconfigurable PET system for imaging plants, with resolution of greater than 3.3 mm, is shown. Using the automated symmetry seeking algorithm, the authors achieved a compression ratio of the storage and memory requirement by a factor of approximately 50 for the half-ring and flat-panel systems. For plant PET system, the compression ratio is approximately five. The ratio depends on the level of symmetry that exists in different geometries. Conclusions: This work brings the field closer to arbitrary geometry reconstruction. A generalized reconstruction framework can be used to validate multiple hypotheses and the effort required to investigate each system is reduced. Memory usage/speed can be improved with certain optimizations.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2002

Space Radiation Transport Methods Development

Improved spacecraft shield design requires early entry of radiation constraints into the design process to maximize performance and minimize costs. As a result, we have been investigating high-speed computational procedures to allow shield analysis from the preliminary design concepts to the final design. In particular, we will discuss the progress towards a full three-dimensional and computationally efficient deterministic code for which the current HZETRN evaluates the lowest order asymptotic term. HZETRN is the first deterministic solution to the Boltzmann equation allowing field mapping within the International Space Station (ISS) in tens of minutes using standard Finite Element Method (FEM) geometry common to engineering design practice enabling development of integrated multidisciplinary design optimization methods. A single ray trace in ISS FEM geometry requires 14 milliseconds and severely limits application of Monte Carlo methods to such engineering models. A potential means of improving the Monte Carlo efficiency in coupling to spacecraft geometry is given in terms of reconfigurable computing and could be utilized in the final design as verification of the deterministic method optimized design.

Wilson, J. W.↗
NASA NTRS2010

Evolution of Geometric Sensitivity Derivatives from Computer Aided Design Models

The generation of design parameter sensitivity derivatives is required for gradient-based optimization. Such sensitivity derivatives are elusive at best when working with geometry defined within the solid modeling context of Computer-Aided Design (CAD) systems. Solid modeling CAD systems are often proprietary and always complex, thereby necessitating ad hoc procedures to infer parameter sensitivity. A new perspective is presented that makes direct use of the hierarchical associativity of CAD features to trace their evolution and thereby track design parameter sensitivity. In contrast to ad hoc methods, this method provides a more concise procedure following the model design intent and determining the sensitivity of CAD geometry directly to its respective defining parameters.

Jones, William T.↗
DOE OSTI2022

Guanosine Dianions Hydrated by One to Four Water Molecules

Intermolecular interactions as those present in molecule…water complexes may profoundly influence the physico-chemical properties of molecules. Here, we carried out an experimentally-computational study on the doubly-deprotonated guanosine monophosphate…water clusters, [dGMP-2H] 2- •nH 2 O (n=1–4), using a combination of negative anion photoelectron spectroscopy (NIPES) with molecular dynamics (MD) and quantum chemical (QM) calculations. Successive addition of water molecules to [dGMP-2H] 2- increases the experimental adiabatic detachment (ADE) and vertical detachment energy (VDE) by 0.5 to 0.1 eV, depending on the cluster size. In order to choose the representative conformations, we combined MD simulations with a clustering procedure to identify low energy geometries for which ADEs and VDEs were computed at the CAM-B3LYP/6-31++G(d,p) level. Our results demonstrate that the assumed approach lead to the sound geometries and energetics of the studied micro-solvates since the calculated ADEs and VDEs are in pretty good agreement with the experimental characteristics. The evolution of hydrogen bonding with cluster size indicates possibility of the occurrence of proton transfer for clusters comprising a larger number of water molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Guidelines for effective radiation transport for cable SGEMP modeling

This report describes experiences gained in performing radiation transport computations with the SCEPTRE radiation transport code for System Generated ElectroMagnetic Pulse (SGEMP) applications. SCEPTRE is a complex code requiring a fairly sophisticated user to run the code effectively, so this report provides guidance for analysts interested in performing these types of calculations. One challenge in modeling coupled photon/electron transport for SGEMP is to provide a spatial mesh that is sufficiently resolved to accurately model surface charge emission and charge deposition near material interfaces. The method that has been most commonly used to date to compute cable SGEMP typically requires a sub-micron mesh size near material interfaces, which may be difficult for meshing software to provide for complex geometries. We present here an alternative method for computing cable SGEMP that appears to substantially relax this requirement. The report also investigates the effect of refining the energy mesh and increasing the order of the angular approximation to provide some guidance on determining reasonable parameters for the energy/angular approximation needed for x-ray environments. Conclusions for γ-ray environments may be quite different and will be treated in a subsequent report. In the course of the energy-mesh refinement studies, a bug in the cross-section generation software was discovered that may cause underprediction of the result by as much as an order of magnitude for the test problem studied here, when the electron energy group widths are much smaller than those for the photons. Results will be presented and compared using cross sections generated before and after the fix. We also describe adjoint modeling, which provides sensitivity of the total charge drive to the source energy and angle of incidence, which is quite useful for comparing the effect of changing the source environment and for determining most stressing angle of incidence and source energy. This report focusses on cable SGEMP applications, but many of the conclusions will be directly applicable for box Internal ElectroMagnetic Pulse (IEMP) modeling as well.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS1977

Ray trajectories in a torus: An application of MACSYMA to a complex numerical computation

The study of ray trajectories of plasma waves in a torodial geometry using MACSYMA is an example of how symbolic, numerical, and graphical facilities can be used in concert to accomplish a complex computational goal. Computational features of this study which are of particular significance include: the derivation of code (i.e. writing functions to generate program fragments), the use of array functions to simplify the specification of a numerical iteration scheme, and the graphical presentation of the results. Mathematically, this study originates in the solution of a linear inhomogeneous partial differential equation in 3 dimensions by the method of characteristics. While it is possible to describe this equation compactly by using vector notation, and by specifying the spatial variation of the coefficients in terms of intermediate parameters, the transformation of the equation into a form amenable to solution is very tedious. A MACSYMA program is presented for obtaining description of the rf field structure excited by a waveguide located at the edge of a toroidal plasma confinement device.

Kulp, J. L.↗