Search NASASearch

SEARCH · Search NASA

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 55 records · Page 3

Unstructured mesh generation and adaptivity

An overview of current unstructured mesh generation and adaptivity techniques is given. Basic building blocks taken from the field of computational geometry are first described. Various practical mesh generation techniques based on these algorithms are then constructed and illustrated with examples. Issues of adaptive meshing and stretched mesh generation for anisotropic problems are treated in subsequent sections. The presentation is organized in an education manner, for readers familiar with computational fluid dynamics, wishing to learn more about current unstructured mesh techniques.

Mavriplis, D. J.

GEODYN system support program, volume 4

The GEODYN Orbit Determination and Geodetic Parameter Estimation System consists of a set of computer programs designed to determine and analyze definitive satellite orbits and their associated geodetic and measurement parameters. This manual describes the Support Programs used by the GEODYN System. The mathematics and programming descriptions are detailed. The operational procedures of each program are presented. GEODYN ancillary analysis programs may be grouped into three different categories: (1) orbit comparison - DELTA (2) data analysis using reference orbits - GEORGE, and (3) pass geometry computations - GROUNDTRACK. All of the above three programs use one or more tapes written by the GEODYN program in either a data reduction or orbit generator run.

Mullins, N. E.

Low‐dimensional manifold learning for uncertainty quantification in complex multi‐scale stochastic systems

Broadly speaking, the goals of the project are to develop techniques to use manifold learning to develop reduced‐order and surrogate models for "hyper‐reduction" of very high‐dimensional complex multi‐scale systems. This is being achieved by employing a newly proposed form of manifold projection and learning that leverages recent advancements in computational geometry and data‐driven modeling. In particular, we are applying a manifold projection technique to project the solutions of very high‐dimensional systems onto the so‐called Grassmannmanifold, a Reimannian manifold comprised of orthonormal matrices. We then apply data‐driven machine learning techniques to classify the solutions on the manifold (e.g. clustering techniques) according to their proximity on the manifold and leverage a further nonlinear dimension reduction to organize the structured data on the manifold. Finally, we are developing novel techniques that enable us to directly interpolate the hyper‐reduced data such that we can predict the solution of the complex, high‐ dimensional system without need to call the full expensive computational model. Given their adherence to the underlying structure of the solution of the physical system, it is expected that these approximate solutions will be sufficiently constrained so as to (approximately) adhere to physical principles.

97 MATHEMATICS AND COMPUTING

Algorithm 1049: The Delaunay Density Diagnostic

Accurate approximation of a real-valued function depends on two aspects of the available data: the density of inputs within the domain of interest and the variation of the outputs over that domain. There are few methods for assessing whether the density of inputs is sufficient to identify the relevant variations in outputs—i.e., the “geometric scale” of the function—despite the fact that sampling density is closely tied to the success or failure of an approximation method. In this article, we introduce a general purpose, computational approach to detecting the geometric scale of real-valued functions over a fixed domain using a deterministic interpolation technique from computational geometry. The algorithm is intended to work on scalar data in moderate dimensions (2–10). Our algorithm is based on the observation that a sequence of piecewise linear interpolants will converge to a continuous function at a quadratic rate (in L 2 norm) if and only if the data are sampled densely enough to distinguish the feature from noise (assuming sufficiently regular sampling). We present numerical experiments demonstrating how our method can identify feature scale, estimate uncertainty in feature scale, and assess the sampling density for fixed (i.e., static) datasets of input–output pairs. Finally, we include analytical results in support of our numerical findings and have released lightweight code that can be adapted for use in a variety of data science settings.

97 MATHEMATICS AND COMPUTING

A Geometric Volume of Fluid-Based Multiphase Flow Solver Extension to the Reacting Flow Solver, PeleLM

A new algorithm is presented to simulate multiphase flows with surface tension in a pathway for spray combustion simulation. The algorithm combines capabilities from two open-source packages, including the interface reconstruction library (IRL), a library of computational geometry routines to enable the volume of fluid (VOF) method, and PeleLM, a solver for the reacting Navier-Stokes equations. Additionally, surface tension is implemented using the continuum surface force (CSF) model with an improved height function technique in the volume fraction field. Spurious errors in volume fraction arising from our combined strategy are corrected through a topology-based method that improves both numerical stability and accuracy. Multiple validation simulations are conducted, including (i) translations and rotations of Zalesak's disk, (ii) a stationary circular droplet with surface tension, (iii) an oscillating elliptical droplet, and (iv) three-dimensional deformation of a spherical droplet. Results indicate that the combined scheme retains the favorable properties of each of the component algorithms.

42 ENGINEERING

Radiofrequency sheath rectification on WEST: application of the sheath-equivalent dielectric layer technique in tokamak geometry *

Radiofrequency sheath rectification is a phenomenon relevant to the operation of Ion Cyclotron Range of Frequencies (ICRFs) actuators in tokamaks. Techniques to model the sheath rectification on 3D ICRF antenna geometries have only recently become available (Shiraiw et al 2023 Nucl. Fusion 63 026024; Beers et al 2021 Phys. Plasmas 28 093503). In this work, we apply the 'sheath-equivalent dielectric layer' technique, used previously only on linear devices (Beers et al 2021 Phys. Plasmas 28 103508), in tokamak geometry, computing rectified sheath potentials on the WEST ICRF antenna. Advancing the state of the art in sheath rectification modeling, we compute the sheath potentials not just on the limiters, but also on the Faraday Screen bars. The calculations show a peak rectified DC potential of 300 V on the limiters and 500 V on the Faraday screen. Assuming a typical sputtering yield curve, the RF sheath rectification increases the sputtering yield from the limiters by a factor of 2.6 w.r.t. the sputtering due to the non-rectified thermal sheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Assessment of Microstructure Prediction Capabilities for Powder Bed Fusion Stainless Steel 316

The Advanced Materials and Manufacturing Technologies program aims to accelerate the development, qualification, demonstration, and deployment of advanced materials and manufacturing technologies to enable reliable and economical nuclear energy. However, the characteristic process-structure-property relationships of additive manufacturing (AM) materials pose challenges for the qualification and certification of AM nuclear components. In particular, component-scale variations in microstructure and properties can be driven by localized changes in melt pool dynamics due to how process parameters interact with different part geometries. Computational modeling tools can play a crucial role in predicting and controlling this variability. This report presents final results on process modeling tools designed to predict microstructure variability in additively manufactured stainless steel 316 parts. It details the software packages and physical modeling approaches employed to simulate an AM component within an automated process modeling workflow. Results are demonstrated through comparisons between predicted microstructures and experimental measurements across various representative processing conditions. The report concludes by discussing identified challenges and future opportunities for connecting the developed simulation workflow with mechanics simulations for prediction of part performance.

36 MATERIALS SCIENCE

Finite element formulations for compressible flows

Researchers started their studies on the development and application of computational methods for compressible flows. Particular attention was given to proper numerical treatment of sharp layers occurring in such problems and to general mesh generation capabilities for intricate computational geometries. Mainly finite element methods enhanced with several state-of-the art techniques (such as the streamline-upwind/Petrov-Galerkin, discontinuity capturing, adaptive implicit-explicit, and trouped element-by-element approximate factorization schemes) were employed.

Tezduyar, Tayfun E.

An ab initio study of Fe(CO)n, n = 1,5, and Cr(CO)6

Ab initio calculations have been performed for Cr(CO)6 and Fe(CO)n, n = 1,5. Basis sets of better than double zeta quality are used, and correlation is included using the modified coupler-pair functional method. The computed geometries and force constants are in reasonable agreement with experiment. The sequential bond dissociation energies of CO from Fe(CO)5 are estimated to be: 39, 31, 25, 22, and greater than 5 kcal/mol. It is noted that the first bond dissociation energy is relative to the singlet ground state of Fe(CO)5 and the lowest singlet state of Fe(CO)4, whereas the second is relative to the ground triplet states of Fe(CO)4 and Fe(CO)3. In addition, the binding energy for Fe-CO would be modified to 18 kcal/mol if dissociation occurred to the Fe(5F) excited state asymptote. The CO binding energies for Fe and Cr are found to be in poorer agreement with experiment than those found in a previous study on Ni(CO)4. The origins of this difference are discussed.

Barnes, Leslie A.

Emerging CFD technologies and aerospace vehicle design

With the recent focus on the needs of design and applications CFD, research groups have begun to address the traditional bottlenecks of grid generation and surface modeling. Now, a host of emerging technologies promise to shortcut or dramatically simplify the simulation process. This paper discusses the current status of these emerging technologies. It will argue that some tools are already available which can have positive impact on portions of the design cycle. However, in most cases, these tools need to be integrated into specific engineering systems and process cycles to be used effectively. The rapidly maturing status of unstructured and Cartesian approaches for inviscid simulations makes suggests the possibility of highly automated Euler-boundary layer simulations with application to loads estimation and even preliminary design. Similarly, technology is available to link block structured mesh generation algorithms with topology libraries to avoid tedious re-meshing of topologically similar configurations. Work in algorithmic based auto-blocking suggests that domain decomposition and point placement operations in multi-block mesh generation may be properly posed as problems in Computational Geometry, and following this approach may lead to robust algorithmic processes for automatic mesh generation.

Aftosmis, Michael J.

On the Sensitivity of Piezoceramics and Piezopolymers in Structural Integrity Monitoring of Large Trusses

An analytical assessment has been made of the reliability of using integrated microactuators and sensors in the form of piezoceramics and piezopolymers as joint integrity monitors in trussed systems. The concept is first implemented for a simple structure which consists of two truss members with a 45 deg lift angle joined at the apex. A piezoceramic patch (or piezopolymer film) bonded on the surface of one of the members at a location near the joint is used as a collocated actuator/sensor. The overall structural dynamic response under an excitation was modeled by finite element method. Different degrees of nodal constraints at the joints representing various degrees of joint integrity are employed. The resulting dynamic response showed distinct responses for varying joint stiffnesses. Parallel experimental work on a truss model using a multichannel data acquisition system and a digital signal analyzer confirms the results from analysis. We further studied the sensitivity of the micro-sensors to the behavior of joints of large arch truss structure. Results obtained for large trusses with many degrees of freedom indicate optimum locations of sensors for which the dynamic response signatures are distinct and distinguishable for relatively small changes in joint integrity and/or structural geometry. Computations based on finite element modeling show that locating the single actuator/sensor at the joint corresponding to the first loss of static stability appear optimal. Hence, static stability analysis of complex trusses can give us a good indication of the optimum placement of sensors for maximum response. This observation is important if few distributed sensors and actuators are available for placement in constructed facilities made from large trusses with many degrees of freedom. As an extension of this work a dynamic response signature identification technique to monitor in-service degradation of joints is under development for application to the monitoring of the integrity of adhesive joints in composite structures.

Abatan, A. O.

A Comparison of the Accuracy of Different Functionals

The errors in the computed geometries, zero-point energies, and atomization energies of molecules containing only first and second row atoms are compared for several levels of theory, including Hartree-Fock, second-order Moller-Plesset perturbation theory (MP2), and density functional theory (DFT) using five different functionals, including two hybrid functionals. The 6-31G* and 6-311+G(3df,2p) basis set are used. Overall, the B3LYP hybrid functional yields the best results.

Bauschlicher, Charles W., Jr.

Faster Aerodynamic Simulation With Cart3D

A NASA-developed aerodynamic simulation tool is ensuring the safety of future space operations while providing designers and engineers with an automated, highly accurate computer simulation suite. Cart3D, co-winner of NASA's 2002 Software of the Year award, is the result of over 10 years of research and software development conducted by Michael Aftosmis and Dr. John Melton of Ames Research Center and Professor Marsha Berger of the Courant Institute at New York University. Cart3D offers a revolutionary approach to computational fluid dynamics (CFD), the computer simulation of how fluids and gases flow around an object of a particular design. By fusing technological advancements in diverse fields such as mineralogy, computer graphics, computational geometry, and fluid dynamics, the software provides a new industrial geometry processing and fluid analysis capability with unsurpassed automation and efficiency.

Source record

Geometric Reasoning for Automated Planning

An important aspect of mission planning for NASA s operation of the International Space Station is the allocation and management of space for supplies and equipment. The Stowage, Configuration Analysis, and Operations Planning teams collaborate to perform the bulk of that planning. A Geometric Reasoning Engine is developed in a way that can be shared by the teams to optimize item placement in the context of crew planning. The ISS crew spends (at the time of this writing) a third or more of their time moving supplies and equipment around. Better logistical support and optimized packing could make a significant impact on operational efficiency of the ISS. Currently, computational geometry and motion planning do not focus specifically on the optimized orientation and placement of 3D objects based on multiple distance and containment preferences and constraints. The software performs reasoning about the manipulation of 3D solid models in order to maximize an objective function based on distance. It optimizes for 3D orientation and placement. Spatial placement optimization is a general problem and can be applied to object packing or asset relocation.

Clement, Bradley J.

Simulation of Complex Cracking in Plain Weave C/SiC Composite under Biaxial Loading

Finite element analysis is performed on a mesh, based on computed geometry of a plain weave C/SiC composite with assumed internal stacking, to reveal the pattern of internal damage due to biaxial normal cyclic loading. The simulation encompasses intertow matrix cracking, matrix cracking inside the tows, and separation at the tow-intertow matrix and tow-tow interfaces. All these dissipative behaviors are represented by traction-separation cohesive laws. Not aimed at quantitatively predicting the overall stress-strain relation, the simulation, however, does not take the actual process of fiber debonding into account. The fiber tows are represented by a simple rule-of-mixture model where the reinforcing phase is a hypothetical one-dimensional material. Numerical results indicate that for the plain weave C/SiC composite, 1) matrix-crack initiation sites are primarily determined by large intertow matrix voids and interlayer tow-tow contacts, 2) the pattern of internal damage strongly depends on the loading path and initial stress, 3) compressive loading inflicts virtually no damage evolution. KEY WORDS: ceramic matrix composite, plain weave, cohesive model, brittle failure, smeared crack model, progressive damage, meso-mechanical analysis, finite element.

Cheng, Ron-Bin

An Overview of Modifications Applied to a Turbulence Response Analysis Method for Flexible Aircraft Configurations

A software program and associated methodology to study gust loading on aircraft exists for a classification of geometrically simplified flexible configurations. This program consists of a simple aircraft response model with two rigid and three flexible symmetric degrees of freedom and allows for the calculation of various airplane responses due to a discrete one-minus-cosine gust as well as continuous turbulence. Simplifications, assumptions, and opportunities for potential improvements pertaining to the existing software program are first identified, then a revised version of the original software tool is developed with improved methodology to include more complex geometries, additional excitation cases, and output data so as to provide a more useful and accurate tool for gust load analysis. Revisions are made in the categories of aircraft geometry, computation of aerodynamic forces and moments, and implementation of horizontal tail mode shapes. In order to improve the original software program to enhance usefulness, a wing control surface and horizontal tail control surface is added, an extended application of the discrete one-minus-cosine gust input is employed, a supplemental continuous turbulence spectrum is implemented, and a capability to animate the total vehicle deformation response to gust inputs in included. These revisions and enhancements are implemented and an analysis of the results is used to validate the modifications.

Funk, Christie J.

Dynamic Airspace Configuration

In air traffic management systems, airspace is partitioned into regions in part to distribute the tasks associated with managing air traffic among different systems and people. These regions, as well as the systems and people allocated to each, are changed dynamically so that air traffic can be safely and efficiently managed. It is expected that new air traffic control systems will enable greater flexibility in how airspace is partitioned and how resources are allocated to airspace regions. In this talk, I will begin by providing an overview of some previous work and open questions in Dynamic Airspace Configuration research, which is concerned with how to partition airspace and assign resources to regions of airspace. For example, I will introduce airspace partitioning algorithms based on clustering, integer programming optimization, and computational geometry. I will conclude by discussing the development of a tablet-based tool that is intended to help air traffic controller supervisors configure airspace and controllers in current operations.

airspace

Aircraft geometry verification with enhanced computer generated displays

A method for visual verification of aerodynamic geometries using computer generated, color shaded images is described. The mathematical models representing aircraft geometries are created for use in theoretical aerodynamic analyses and in computer aided manufacturing. The aerodynamic shapes are defined using parametric bi-cubic splined patches. This mathematical representation is then used as input to an algorithm that generates a color shaded image of the geometry. A discussion of the techniques used in the mathematical representation of the geometry and in the rendering of the color shaded display is presented. The results include examples of color shaded displays, which are contrasted with wire frame type displays. The examples also show the use of mapped surface pressures in terms of color shaded images of V/STOL fighter/attack aircraft and advanced turboprop aircraft.

Cozzolongo, J. V.