Search NASA⌕ Search

SEARCH · Search NASA

Results for “polyhedral model”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

A quantitative approach to the topology of large-scale structure

A quantitative measure of the topology of large-scale structure: the genus of density contours in a smoothed density distribution, is described and applied. For random phase (Gaussian) density fields, the mean genus per unit volume exhibits a universal dependence on threshold density, with a normalizing factor that can be calculated from the power spectrum. If large-scale structure formed from the gravitational instability of small-amplitude density fluctuations, the topology observed today on suitable scales should follow the topology in the initial conditions. The technique is illustrated by applying it to simulations of galaxy clustering in a flat universe dominated by cold dark matter. The technique is also applied to a volume-limited sample of the CfA redshift survey and to a model in which galaxies reside on the surfaces of polyhedral 'bubbles'. The topology of the evolved mass distribution and 'biased' galaxy distribution in the cold dark matter models closely matches the topology of the density fluctuations in the initial conditions. The topology of the observational sample is consistent with the random phase, cold dark matter model.

Gott, J. Richard, III↗
DOE OSTI2020

An Integrated Multiscale Experimental-Numerical Analysis on Reconsolidation of Salt-Clay Mixture for Disposal of Heat-Generating Waste (Final NEUP Technical Report)

The overall purpose of this research is to improve understanding of THMC coupling effect on the reconsolidation of granular (or crushed) salt-clay mixture used for seal systems of shafts and drifts in salt repositories. This proposed work is partially motivated by the recent work on the Waste Isolation Pilot Plant (WIPP) that shows the promising sealing capability of clay-salt mixture compared to crushed salt. In particular, primary emphasis is to develop a fully integrated multiscale experiment-numerical study to determine and explain what leads to the superior sealing ability of the clay-salt mixture. These research activities are designed to seek further understanding of (1) why clay additives may enhance the fluid trapping and (2) whether this flow barrier effect may prevail under different combinations of temperature, confining pressure, deviatoric stress and other foreseeable environmental factors. If successful, this enhanced flow trapping ability of the seal provides significant improvement to the seal and repository performance and therefore make the repository safer in the long-term. The experiment component includes microstructural investigation and macroscopic tests on a reconsolidated salt-clay mixture. In the microstructural study, the goal is to (1) characterize microscopic distributions of distinct phases (e.g., clay, salt crystal boundaries, trapped brine, and pore) to examine the connectivity of the pore network inside the salt-clay mixture with different amounts of clay additive and moisture content and (2) analyze multiscale imaging data to reconstruct the polycrystalline microstructures for numerical simulations. Meanwhile, macroscopic tests are performed to analyze how clay alters the failure/creep mechanisms in the salt-clay mixture. Microscopic and macroscopic experimental observations will both be used to calibrate and validate a multiscale model that explicitly simulates the capillary and multiphase flow in the connected pores and the deformation due to the presence of intra-crystalline brine at the pore scale via a new polyhedral discrete element–lattice Boltzmann method (DEM-LBM) coupling model. The pore-scale simulations are homogenized via an upscaling procedure that converts pore-scale information (e.g. force exerted on grain boundary, sliding, pressure-solution) to continuum measures (e.g. Cauchy stress, Darcy’s flow) at each integration point in the macroscopic multiphase TMHC model. This multiscale scheme will allow coupling be- tween high-fidelity simulations of brine-salt-clay interaction and the macroscopic TMHC model. The multiscale model helps the understanding of how the trapped brine inclusion affects the pressure-solution mechanism with the presence of clay and moisture. This work brings new insight into the sealing capacity of salt-clay mixture under elevated temperature over a long period of time - a key to evaluating the potential of salt-clay mixture usage for salt repositories.

42 ENGINEERING↗
DOE OSTI2022

Hydrogen bond arrangements in (H 2 O) 20, 24, 28 clathrate hydrate cages: Optimization and many-body analysis

Here we provide a detailed study of hydrogen bonding arrangements, relative stability, residual entropy, and an analysis of the many-body effects in the (H 2 O) 20 (D-cage), (H 2 O) 24 (T-cage), and (H 2 O) 28 (H-cage) hollow cages making up structures I (sI) and II (sII) of clathrate hydrate lattices. Based on the enumeration of the possible hydrogen bonding networks for a fixed oxygen atom scaffold, the residual entropy (S 0 ) of these three gas phase cages was estimated at 0.754 82, 0.754 44, and 0.754 17 · Nk b , where N is the number of molecules and k b is Boltzmann’s constant. A previously identified descriptor of enhanced stability based on the relative arrangement and connectivity of nearest-neighbor fragments on the polyhedral water cluster [strong-weak-effective-bond model] also applies to the larger hollow cages. The three cages contain a maximum of 7, 9, and 11 such preferable arrangements of trans nearest dimer pairs with one “free” OH bond on the donor molecule (t1d dimers). The Many-Body Expansion (MBE) up to the 4-body suggests that the many-body terms vary nearly linearly with the cluster binding energy. Using a hierarchical approach of screening the relative stability of networks starting from optimizations with the TIP4P, TTM2.1-F, and MB-pol classical potentials, subsequently refining at more accurate levels of electronic structure theory (DFT and MP2), and finally correcting for zero-point energy, we were able to identify a group of four low-lying isomers of the (H 2 O) 24 T-cage, two of which are antisymmetric and the other two form a pair of antipode configurations.

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

Polyhedral Relaxations for Optimal Pump Scheduling of Potable Water Distribution Networks

The classic pump scheduling or optimal water flow (OWF) problem for water distribution networks (WDNs) minimizes the cost of power consumption for a given WDN over a fixed time horizon. In its exact form, the OWF is a computationally challenging mixed-integer nonlinear program (MINLP). It is complicated by nonlinear equality constraints that model network physics, discrete variables that model operational controls, and intertemporal constraints that model changes to storage devices. To address the computational challenges of the OWF, this paper develops tight polyhedral relaxations of the original MINLP, derives novel valid inequalities (or cuts) using duality theory, and implements novel optimization-based bound tightening and cut generation procedures. The efficacy of each new method is rigorously evaluated by measuring empirical improvements in OWF primal and dual bounds over 45 literature instances. The evaluation suggests that our relaxation improvements, model strengthening techniques, and a thoughtfully selected polyhedral relaxation partitioning scheme can substantially improve OWF primal and dual bounds, especially when compared with similar relaxation-based techniques that do not leverage these new methods.

bound tightening↗
DOE OSTI2022

Multi-service battery energy storage system optimization and control

Battery energy storage systems (BESS) have become fundamental part of modern power systems due to their capability to provide multiple grid services. As the renewable penetration increases, BESS procurement is also expected to increase where it is envisioned to play a systematic and strategical role in power systems planning and operation. Hence, in this paper we present a multiple grid service procurement and operation for BESS - ranging from energy arbitrage, reserve/regulation services, power factor correction, and demand management. The proposed framework considers an optimal multi-temporal dimension, designed to be operable for both planning and real-time operation. Moreover, non-linearity inherent to BESS services and uncertainty associated to market forecasts variables are addressed using techniques such as polyhedral norms and robust optimization approaches. Here, the developed model is tested using a utility-scaled BESS and the obtained results show the effectiveness of the systematic BESS multi-service planning and operation approach.

25 ENERGY STORAGE↗
DOE OSTI2022

Atomistic Simulations of the Elastic Compression of Platinum Nanoparticles

Abstract The elastic behavior of nanoparticles depends strongly on particle shape, size, and crystallographic orientation. Many prior investigations have characterized the elastic modulus of nanoscale particles using experiments or simulations; however their reported values vary widely depending on the methods for measurement and calculation. To understand these discrepancies, we used classical molecular dynamics simulation to model the compression of platinum nanoparticles with two different polyhedral shapes and a range of sizes from 4 to 20 nm, loaded in two different crystal orientations. Multiple standard methods were used to calculate the elastic modulus from stress-vs-strain data for each nanoparticle. The magnitudes and particle-size dependence of the resulting moduli varied with calculation method and, even for larger nanoparticles where bulk-like behavior may be expected, the effective elastic modulus depended strongly on shape and orientation. Analysis of per-atom stress distributions indicated that the shape- and orientation-dependence arise due to stress triaxiality and inhomogeneity across the particle. When the effective elastic modulus was recalculated using a representative volume element in the center of a large nanoparticle, the elastic modulus had the expected value for each orientation and was shape independent. It is only for single-digit nanoparticles that meaningful differences emerged, where even the very center of the particle had a lower modulus due to the effect of the surface. These findings provide better understanding of the elastic properties of nanoparticles and disentangle geometric contributions (such as stress triaxiality and spatial inhomogeneity) from true changes in elastic properties of the nanoscale material.

36 MATERIALS SCIENCE↗
NASA NTRS2017

UAV Inspection of Electrical Transmission Infrastructure with Path Conformance Autonomy and Lidar-Based Geofences NASA Report on UTM Reference Mission Flights at Southern Company Flights November 2016

Flights at low altitudes in close proximity to electrical transmission infrastructure present serious navigational challenges: GPS and radio communication quality is variable and yet tight position control is needed to measure defects while avoiding collisions with ground structures. To advance unmanned aerial vehicle (UAV) navigation technology while accomplishing a task with economic and societal benefit, a high voltage electrical infrastructure inspection reference mission was designed. An integrated air-ground platform was developed for this mission and tested in two days of experimental flights to determine whether navigational augmentation was needed to successfully conduct a controlled inspection experiment. The airborne component of the platform was a multirotor UAV built from commercial off-the-shelf hardware and software, and the ground component was a commercial laptop running open source software. A compact ultraviolet sensor mounted on the UAV can locate 'hot spots' (potential failure points in the electric grid), so long as the UAV flight path adequately samples the airspace near the power grid structures. To improve navigation, the platform was supplemented with two navigation technologies: lidar-to-polyhedron preflight processing for obstacle demarcation and inspection distance planning, and trajectory management software to enforce inspection standoff distance. Both navigation technologies were essential to obtaining useful results from the hot spot sensor in this obstacle-rich, low-altitude airspace. Because the electrical grid extends into crowded airspaces, the UAV position was tracked with NASA unmanned aerial system traffic management (UTM) technology. The following results were obtained: (1) Inspection of high-voltage electrical transmission infrastructure to locate 'hot spots' of ultraviolet emission requires navigation methods that are not broadly available and are not needed at higher altitude flights above ground structures. (2) The sensing capability of a novel airborne UV detector was verified with a standard ground-based instrument. Flights with this sensor showed that UAV measurement operations and recording methods are viable. With improved sensor range, UAVs equipped with compact UV sensors could serve as the detection elements in a self-diagnosing power grid. (3) Simplification of rich lidar maps to polyhedral obstacle maps reduces data volume by orders of magnitude, so that computation with the resultant maps in real time is possible. This enables real-time obstacle avoidance autonomy. Stable navigation may be feasible in the GPS-deprived environment near transmission lines by a UAV that senses ground structures and compares them to these simplified maps. (4) A new, formally verified path conformance software system that runs onboard a UAV was demonstrated in flight for the first time. It successfully maneuvered the aircraft after a sudden lateral perturbation that models a gust of wind, and processed lidar-derived polyhedral obstacle maps in real time. (5) Tracking of the UAV in the national airspace using the NASA UTM technology was a key safety component of this reference mission, since the flights were conducted beneath the landing approach to a heavily used runway. Comparison to autopilot tracking showed that UTM tracking accurately records the UAV position throughout the flight path.

Moore, Andrew J.↗
NASA NTRS1992

Divergence preserving discrete surface integral methods for Maxwell's curl equations using non-orthogonal unstructured grids

Several new discrete surface integral (DSI) methods for solving Maxwell's equations in the time-domain are presented. These methods, which allow the use of general nonorthogonal mixed-polyhedral unstructured grids, are direct generalizations of the canonical staggered-grid finite difference method. These methods are conservative in that they locally preserve divergence or charge. Employing mixed polyhedral cells, (hexahedral, tetrahedral, etc.) these methods allow more accurate modeling of non-rectangular structures and objects because the traditional stair-stepped boundary approximations associated with the orthogonal grid based finite difference methods can be avoided. Numerical results demonstrating the accuracy of these new methods are presented.

Madsen, Niel K.↗
NASA NTRS2014

Coupled Fluid-Structure Interaction Analysis of Solid Rocket Motor with Flexible Inhibitors

Flexible inhibitors are generally used in solid rocket motors (SRMs) as a means to control the burning of propellant. Vortices generated by the flow of propellant around the flexible inhibitors have been identified as a driving source of instabilities that can lead to thrust oscillations in launch vehicles. Potential coupling between the SRM thrust oscillations and structural vibration modes is an important risk factor in launch vehicle design. As a means to predict and better understand these phenomena, a multidisciplinary simulation capability that couples the NASA production CFD code, Loci/CHEM, with CFDRC's structural finite element code, CoBi, has been developed. This capability is crucial to the development of NASA's new space launch system (SLS). This paper summarizes the efforts in applying the coupled software to demonstrate and investigate fluid-structure interaction (FSI) phenomena between pressure waves and flexible inhibitors inside reusable solid rocket motors (RSRMs). The features of the fluid and structural solvers are described in detail, and the coupling methodology and interfacial continuity requirements are then presented in a general Eulerian-Lagrangian framework. The simulations presented herein utilize production level CFD with hybrid RANS/LES turbulence modeling and grid resolution in excess of 80 million cells. The fluid domain in the SRM is discretized using a general mixed polyhedral unstructured mesh, while full 3D shell elements are utilized in the structural domain for the flexible inhibitors. Verifications against analytical solutions for a structural model under a steady uniform pressure condition and under dynamic modal analysis show excellent agreement in terms of displacement distribution and eigenmode frequencies. The preliminary coupled results indicate that due to acoustic coupling, the dynamics of one of the more flexible inhibitors shift from its first modal frequency to the first acoustic frequency of the solid rocket motor. This insight could have profound implications for SRM and flexible inhibitor designs for current and future launch vehicles including SLS.

Yang, H. Q.↗
NASA NTRS1987

On autonomous terrain model acquistion by a mobile robot

The following problem is considered: A point robot is placed in a terrain populated by an unknown number of polyhedral obstacles of varied sizes and locations in two/three dimensions. The robot is equipped with a sensor capable of detecting all the obstacle vertices and edges that are visible from the present location of the robot. The robot is required to autonomously navigate and build the complete terrain model using the sensor information. It is established that the necessary number of scanning operations needed for complete terrain model acquisition by any algorithm that is based on scan from vertices strategy is given by the summation of i = 1 (sup n) N(O sub i)-n and summation of i = 1 (sup n) N(O sub i)-2n in two- and three-dimensional terrains respectively, where O = (O sub 1, O sub 2,....O sub n) set of the obstacles in the terrain, and N(O sub i) is the number of vertices of the obstacle O sub i.

Rao, N. S. V.↗
NASA NTRS2023

Evaluation of Voronoi Meshes for Large Eddy Simulations of High Lift Aerodynamics

Numerical sensitivity to 3 different Voronoi seeding methods is investigated for Large Eddy Simulations (LES). A second order accurate, non-dissipative finite volume discretization is used to systematically investigate the effects of different polyhedral Voronoi mesh types using a sequence of three canonical problems with increasing complexity. First, inviscid isentropic vortex propagation is studied to demonstrate the substantial reduction of errors for rhombic dodecahedron and truncated octahedral cell types over Cartesian hexagonal cells of identical spacing. Furthermore, the reduction of errors induced at cell-size transitions (grid-coarsening interfaces) due to Lloyd smoothing iterations is quantified. It is shown that by utilizing an appropriate viscous flux discretization, a constant coefficient subgrid scale model is sufficient for non-linear stability at the 2:1 cell-size transitions on polyhedral grids, although some further error reduction does occur when smoothing is utilized. Next, forced homogeneous isotropic turbulence at an asymptotically large Reynolds number is studied to demonstrate the non-dissipative character of the inviscid flux discretization, and the non-linear robustness and accuracy offered by the viscous flux discretization using a subgrid scale model for all three Voronoi grid types. Finally,Wall-Modeled Large Eddy Simulations (WMLES) are performed to study the high-lift aerodynamics on the McDonnell Douglas 30P30N multi-element airfoil at two distinct grid levels and for two distinct Voronoi cell types. The formulation is shown to predict the aerodynamic loading with high accuracy at all angles of attack when sufficient resolution is reached, and the hexagonal prism grid topology, while computationally more expensive, has higher effective resolution compared to the Cartesian grid topology with the same spacing.

TTT↗
NASA NTRS2023

Evaluation of Voronoi Meshes for Large Eddy Simulations of High Lift Aerodynamics

Numerical sensitivity to 3 different Voronoi seeding methods is investigated for Large Eddy Simulations (LES). A second order accurate, non-dissipative finite volume discretization is used to systematically investigate the effects of different polyhedral Voronoi mesh types using a sequence of three canonical problems with increasing complexity. First, inviscid isentropic vortex propagation is studied to demonstrate the substantial reduction of errors for rhombic dodecahedron and truncated octahedral cell types over Cartesian hexagonal cells of identical spacing. Furthermore, the reduction of errors induced at cell-size transitions (grid-coarsening interfaces) due to Lloyd smoothing iterations is quantified. It is shown that by utilizing an appropriate viscous flux discretization, a constant coefficient subgrid scale model is sufficient for non-linear stability at the 2:1 cell-size transitions on polyhedral grids, although some further error reduction does occur when smoothing is utilized. Next, forced homogeneous isotropic turbulence at an asymptotically large Reynolds number is studied to demonstrate the non-dissipative character of the inviscid flux discretization, and the non-linear robustness and accuracy offered by the viscous flux discretization using a subgrid scale model for all three Voronoi grid types. Finally, Wall-Modeled Large Eddy Simulations (WMLES) are performed to study the high-lift aerodynamics on the McDonnell Douglas 30P30N multi-element airfoil at two distinct grid levels and for two distinct Voronoi cell types. The formulation is shown to predict the aerodynamic loading with high accuracy at all angles of attack when sufficient resolution is reached, and the hexagonal prism grid topology, while computationally more expensive, has higher effective resolution compared to the Cartesian grid topology with the same spacing.

TTT↗
NASA NTRS1987

Real-time model-based vision system for object acquisition and tracking

A machine vision system is described which is designed to acquire and track polyhedral objects moving and rotating in space by means of two or more cameras, programmable image-processing hardware, and a general-purpose computer for high-level functions. The image-processing hardware is capable of performing a large variety of operations on images and on image-like arrays of data. Acquisition utilizes image locations and velocities of the features extracted by the image-processing hardware to determine the three-dimensional position, orientation, velocity, and angular velocity of the object. Tracking correlates edges detected in the current image with edge locations predicted from an internal model of the object and its motion, continually updating velocity information to predict where edges should appear in future frames. With some 10 frames processed per second, real-time tracking is possible.

Wilcox, Brian↗
DOE OSTI2019

Analytical Modeling of Cache Behavior for Affine Programs

Modern optimizing compiler design integrates program transformation strategies aimed at reducing data movement to/from main memory by exploiting the data cache hierarchy. But in reality, the effect of these transformations on the actual cache miss count is ignored in compilers due to the lack of precise compile-time models of misses for hierarchical caches. The current state of practice for cache miss analysis is based on accurate simulation. However, simulation requires time proportional to the dataset/problem size as well as the number of distinct cache configurations of interest to be evaluated. This paper takes a fundamentally different approach, made possible by focusing on polyhedral programs with static control flow: instead of relying on costly simulation, this paper develops the first closed-form solution for exact modeling of misses in a set-associative cache hierarchy. It enabling simulation-less selection of program transformations at compile-time to optimize cache misses. A push-button tool implementing the approach is developed and used for extensive validation of the complete framework.

affine programs, cache modeling, set associativity↗
DOE OSTI2020

Phase Diagram and Structure Map of Binary Nanoparticle Superlattices from a Lennard-Jones Model

A first-principles prediction of the binary nanoparticle phase diagram assembled by solvent evaporation has eluded theoretical approaches. In this paper, we show that a binary system interacting through the Lennard-Jones (LJ) potential contains all experimental phases in which nanoparticles are effectively described as quasi hard spheres. We report a phase diagram consisting of 53 equilibrium phases, whose stability is quite insensitive to the microscopic details of the potentials, thus giving rise to some type of universality. Furthermore, we show that binary lattices may be understood as consisting of certain particle clusters, i.e., motifs, that provide a generalization of the four conventional Frank–Kasper polyhedral units. Our results show that metastable phases share the very same motifs as equilibrium phases. Finally, we discuss the connection with packing models, phase diagrams with repulsive potentials, and the prediction of likely experimental superlattices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1988

Model-based orientation-independent 3-D machine vision techniques

Orientation-dependent techniques for the identification of a three-dimensional object by a machine vision system are represented in parts. In the first part, the data consist of intensity images of polyhedral objects obtained by a single camera, while in the second part, the data consist of range images of curved objects obtained by a laser scanner. In both cases, the attributed graphic representation of the object surface is used to drive the respective algorithm. In this representation, a graph node represents a surface patch and a link represents the adjacency between two patches. The attributes assigned to nodes are moment invariants of the corresponding face for polyhedral objects. For range images, the Gaussian curvature is used as a segmentation criterion for providing symbolic shape attributes. Identification is achieved by an efficient graph-matching algorithm used to match the graph obtained from the data to a subgraph of one of the model graphs stored in the commputer memory.

De Figueiredo, R. J. P.↗
DOE OSTI2025

Improvements to MOOSE user workflow through polyhedral elements, automation, and concise physics syntax

The MOOSE framework is a foundational capability used by the NEAMS program to create over 15 different simulation tools for advanced nuclear reactors. Due to MOOSE's broad use, improvements to the framework in support of modeling and simulation goals are critical to the program. Such improvements can take many forms, including optimization, improved user experience, streamlined application programming interfaces (APIs), parallelism, and new capabilities. The work described in this report was conducted in direct support of NEAMS tools and includes: addition of support for polyhedral elements, incorporation of mesh smoothers for mesh repair, integration of the Physics and ActionComponents systems, expansion of the Convergence system, and exploration of automated input file generation. These five areas of development are fundamental capabilities that will be leveraged by many NEAMS applications.

97 - MATHEMATICS AND COMPUTING↗
NASA NTRS2018

Analysis of High-Speed Aerodynamics of a Swept Wing with Seamless Flaps

Computational fluid dynamics (CFD) analysis was conducted to evaluate the high-speed aerodynamics of a Gulfstream GIII airplane (Gulfstream Aerospace Corporation, Savannah, Georgia) swept wing modified with an experimental seamless, compliant flap called the Adaptive Compliant Trailing Edge (ACTE) flap for airworthiness. The high-speed aerodynamics of the modified ACTE wing were analyzed at a Mach number of 0.85 and an altitude of 40,000 ft. A polyhedral finite-volume unstructured Reynolds-averaged Navier-Stokes CFD code was used with a two-equation turbulence model. Wing aerodynamics sensitivities to lift coefficients, airplane engine operation, vortex generators (VGs), and ACTE flap deflections were examined. High-speed flow over the clean GIII wing without VGs was found to be complex with outboard shock and boundary-layer interactions causing flow separation over the upper wing surface upstream of the aileron, which could result in buffet and controllability problems at higher lift coefficient values. The VGs were found to improve the high-speed aerodynamics of the clean wing, as expected; however, ACTE flap deflections of as little as 1 to 2 deg were found to negate all of the aerodynamic benefits provided by the VGs at high speeds.

wings↗