SEARCH · Engineering Papers

Results for “Voronoi tessellation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

33 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Electrical conductivity, ion pairing, and ion self-diffusion in aqueous NaCl solutions at elevated temperatures and pressures

We have performed classical molecular dynamics (MD) simulations of aqueous sodium chloride (NaCl) solutions from 298 to 674 K at 200 bars to understand the influence of ion pairing and ion self-diffusion on electrical conductivity in high-temperature/high-pressure salt solutions. Conductivity data obtained from the MD simulation highlight an apparent anomaly, namely, a conductivity maximum as temperature increases along an isobar, which has been also observed in experimental studies. By examining both velocity autocorrelation and cross-correlation terms of the Green-Kubo integral, we quantitatively demonstrate that the conductivity anomaly arises mainly from a competition between the single-ion self-diffusion and the contact ion pair formation. The velocity autocorrelation function in conjunction with structural analysis suggests that diffusive motion of ions is suppressed at high temperatures due to the persistence of an inner hydration shell. Here, the contribution of velocity cross-correlation functions between oppositely charged ions becomes significant at the onset of the conductivity decrease. Structural analysis based on Voronoi tessellation and pair correlation functions indicates that the fraction of contact ion pairs increases as temperature increases. Spatial decomposition of the electrical conductivity also indicates that the formation of contact ion pairs significantly decreases the electrical conductivity compared to Nernst-Einstein conductivity, but the contribution of distant opposite charges cannot be ignored except at the highest temperature due to unscreened long-range interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Detecting depinning and nonequilibrium transitions with unsupervised machine learning

Using numerical simulations of a model disk system, we demonstrate that a machine learning generated order-parameter-like measure can detect depinning transitions and different dynamic flow phases in systems driven far from equilibrium. Furthermore, we specifically consider monodisperse passive disks with short range interactions undergoing a depinning phase transition when driven over quenched disorder. The machine learning derived order-parameter-like measure identifies the depinning transition as well as different dynamical regimes, such as the transition from a flowing liquid to a phase separated liquid-solid state that is not readily distinguished with traditional measures such as velocity-force curves or Voronoi tessellation. The order-parameter-like measure also shows markedly distinct behavior in the limit of high density where jamming effects occur. Our results should be general to the broad class of particle-based systems that exhibit depinning transitions and nonequilibrium phase transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Amesh v1.0

Amesh generates discretized grids for numerical modeling of flow and transport problems in which the formulation is based on the integral finite difference method (IFDM). For example, the output of Amesh can be used directly as input to the TOUGH2 or TOUGH3 numerical simulators. Amesh can generate 1D, 2D, or 3D numerical grids for a given set of locations, i.e., the centers of each discrete sub-domain. In the 2D areal plane Voronoi tessellation is used.

Fuller, Peter↗

DOE OSTI2022

MPAS-Seaice (v1.0.0): sea-ice dynamics on unstructured Voronoi meshes

Abstract. We present MPAS-Seaice, a sea-ice model which uses the Model for Prediction Across Scales (MPAS) framework and spherical centroidal Voronoi tessellation (SCVT) unstructured meshes. As well as SCVT meshes, MPAS-Seaice can run on the traditional quadrilateral grids used by sea-ice models such as CICE. The MPAS-Seaice velocity solver uses the elastic–viscous–plastic (EVP) rheology and the variational discretization of the internal stress divergence operator used by CICE, but adapted for the polygonal cells of MPAS meshes, or alternatively an integral (“finite-volume”) formulation of the stress divergence operator. An incremental remapping advection scheme is used for mass and tracer transport. We validate these formulations with idealized test cases, both planar and on the sphere. The variational scheme displays lower errors than the finite-volume formulation for the strain rate operator but higher errors for the stress divergence operator. The variational stress divergence operator displays increased errors around the pentagonal cells of a quasi-uniform mesh, which is ameliorated with an alternate formulation for the operator. MPAS-Seaice shares the sophisticated column physics and biogeochemistry of CICE and when used with quadrilateral meshes can reproduce the results of CICE. We have used global simulations with realistic forcing to validate MPAS-Seaice against similar simulations with CICE and against observations. We find very similar results compared to CICE, with differences explained by minor differences in implementation such as with interpolation between the primary and dual meshes at coastlines. We have assessed the computational performance of the model, which, because it is unstructured, runs with 70 % of the throughput of CICE for a comparison quadrilateral simulation. The SCVT meshes used by MPAS-Seaice allow removal of equatorial model cells and flexibility in domain decomposition, improving model performance. MPAS-Seaice is the current sea-ice component of the Energy Exascale Earth System Model (E3SM).

58 GEOSCIENCES↗

NASA NTRS1991

The evolution of void-filled cosmological structures

1D, 2D, and 3D simulations are used here to investigate the salient features in the evolution of void-filled cosmological structures in universes with arbitrary values of Omega. It is found that the growth of a void as a function of time decreases significantly at the time corresponding to Omega = 0.5. In models constructed in 2D and 3D, suitable initial conditions lead to cellular structure with faceted voids similar to those observed in redshift surveys. Matter compressed to planes flows more rapidly toward condensations at the intersections than would be expected for spherical infall. The peculiar streaming velocities for void diameters of 5000 km/s should be observable. The simulations provide a more physical basis and dynamics for the bubbly and Voronois tesselation models used to derive statistical properties of cellular large-scale structure.

Regos, Eniko↗

NASA NTRS2000

Assessment of Internal Fabric of Particulate Materials

Particle arrangement and distribution within a soil matrix has long been recognized as having significant influence on the mechanical behavior of cohesionless soils. It is well known that two soil specimens having the same grain type (e.g., quartz, feldspar, etc.), same grain size distribution and relative density (or void ratio) can display completely different mechanical behavior. Because of the different fabric configurations in the otherwise similar specimens, they are likely to have different mechanical properties such as stress-strain response, anisotropy, dilatancy, etc. Soil Fabric is defined as the arrangement of particles, particle groups and associated pore space. In the literature, fabric analysis techniques are mainly classified as destructive (e.g., specimen stabilization, thin-sectioning, and microscopy), and nondestructive techniques (e.g., magnetic resonance imaging, ultrasonic testing, x-ray radiography, and computed tomography). Quantifying the void ratio and its distribution is the main parameter used to describe the fabric of particulate materials. This paper presents a comprehensive literature review of fabric analysis techniques applied to particulate materials. In addition, the results of a comprehensive investigation to quantify void ratio of sand specimens will be presented and discussed. The sand used in the experiments in a natural, uniform rounded to sub-rounded silica sand known as F-75 banding sand with mean particle size of 0.22 mm. Uniform specimens and specimens subjected to different axial-strain levels tested under triaxial and biaxial conditions are examined to evaluate void ratio evolution and its distribution using destructive thin-sectioning and nondestructive Computed Tomography (CT) techniques. Details of a new innovative polygon generation technique called Voronoi tessellation used to quantify void ratio of microscopic images of sand grains will be presented and compared to classical Oda's method. Finally, frequency distribution of void ratio calculated using CT analysis will be presented and discussed.

Alshibi, Khalid A.↗

NASA NTRS2005

Modeling for Ultrasonic Health Monitoring of Foams with Embedded Sensors

In this report analytical and numerical methods are proposed to estimate the effective elastic properties of regular and random open-cell foams. The methods are based on the principle of minimum energy and on structural beam models. The analytical solutions are obtained using symbolic processing software. The microstructure of the random foam is simulated using Voronoi tessellation together with a rate-dependent random close-packing algorithm. The statistics of the geometrical properties of random foams corresponding to different packing fractions have been studied. The effects of the packing fraction on elastic properties of the foams have been investigated by decomposing the compliance into bending and axial compliance components. It is shown that the bending compliance increases and the axial compliance decreases when the packing fraction increases. Keywords: Foam; Elastic properties; Finite element; Randomness

Wang, L.↗

NASA NTRS2019

Parallel High-Order Anisotropic Meshing Using Discrete Metric Tensors

This paper presents a metric-aligned meshing algorithm that relies on the Lp-Centroidal Voronoi Tesselation approach. A prototype of this algorithm was first presented at the Scitech conference of 2018 and this work is an extension to that paper. At the end of the previously presented work, a set of problems were mentioned which we are trying to address in this paper. First, we show a significant improvement in code performance since we were limited to present relatively benign (analytical) test cases. Second, we demonstrate here that we are able to rely on discrete metric data that is delivered by a Computational Fluid Dynamics (CFD) solver. Third, we demonstrate how to generate high-order curved elements that are aligned with the underlying discrete metric field.

Ekelschot, Dirk↗

DOE OSTI2020

Finding wombling boundaries in LHC data with Voronoi and Delaunay tessellations

We address the problem of finding a wombling boundary in point data generated by a general Poisson point process, a specific example of which is an LHC event sample distributed in the phase space of a final state signature, with the wombling boundary created by some new physics. We discuss the use of Voronoi and Delaunay tessellations of the point data for estimating the local gradients and investigate methods for sharpening the boundaries by reducing the statistical noise. The outcome from traditional wombling algorithms is a set of boundary cell candidates with relatively large gradients, whose spatial properties must then be scrutinized in order to construct the boundary and evaluate its significance. Here we propose an alternative approach where we simultaneously form and evaluate the significance of all possible boundaries in terms of the total gradient flux. We illustrate our method with several toy examples of both straight and curved boundaries with varying amounts of signal present in the data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Alternatives to Contour Visualizations for Power Systems Data

Electrical grids are geographical and topological structures whose voltage states are challenging to represent accurately and efficiently for visual analysis. The current common practice is to use colored contour maps, yet these can misrepresent the data. We examine the suitability of four alternative visualization methods for depicting voltage data in a geographically dense distribution system - Voronoi polygons, H3 tessellations, S2 tessellations, and a network-weighted contour map. We find that Voronoi tessellations and network-weighted contour maps more accurately represent the statistical distribution of the data than regular contour maps.

MATHEMATICS AND COMPUTING,POWER TRANSMISSION AND D↗

DOE OSTI2023

Alternatives to Contour Visualizations for Power Systems Data

Electrical grids are geographical and topological structures whose voltage states are challenging to represent accurately and efficiently for visual analysis. The current common practice is to use colored contour maps, yet these can misrepresent the data. We examine the suitability of four alternative visualization methods for depicting voltage data in a geographically dense distribution system-Voronoi polygons, H3 tessellations, S2 tessellations, and a network-weighted contour map. We find that Voronoi tessellations and network-weighted contour maps more accurately represent the statistical distribution of the data than regular contour maps.

MATHEMATICS AND COMPUTING↗

DOE OSTI2024

T RI M E ++: Multi-threaded triangular meshing in two dimensions

We present T RI M E ++, a multi-threaded software library designed for generating two-dimensional meshes for intricate geometric shapes using the Delaunay triangulation. Multi-threaded parallel computing is implemented throughout the meshing procedure, making it suitable for fast generation of large-scale meshes. Three iterative meshing algorithms are implemented: the DistMesh algorithm, the centroidal Voronoi diagram meshing, and a hybrid of the two. We compare the performance of the three meshing methods in T RI M E ++, and show that the hybrid method retains the advantages of the other two. The software library achieves significant parallel speedup when generating large-scale meshes containing between 10 4 to 10 7 points. T RI M E ++ can handle complicated geometries and generates adaptive meshes of high quality.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

ABC -Auxetics: An Implicit Design Approach for Negative Poisson's Ratio Materials

A novel methodology is introduced for designing auxetic (negative Poisson's ratio) structures based on topological principles and is demonstrated by investigating a new class of auxetics based on two-dimensional (2D) textile weave patterns. Conventional methodology for designing auxetic materials typically involves determining a single deformable block (a unit cell) of material whose shape results in auxetic behavior. Consequently, patterning such a unit cell in a 2D (or 3D) domain results in a larger structure that exhibits overall auxetic behavior. Such an approach naturally relies on some prior intuition and experience regarding which unit cells may be auxetic. Second, tuning the properties of the resulting structures is typically limited to parametric variations of the geometry of a specific type of unit cell. Thus, most of the currently known auxetic structures belong to a selected few classes of unit cell geometries that are explicitly defined in accordance with a specified topological (i.e., grid structure). Herein, a new class of auxetic structures is demonstrated that, while periodic, can be generated implicitly, i.e., without reference to a specific unit cell design. The approach leverages weave-based parameters (A–B–C), resulting in a rich design space for auxetics that is previously unexplored.

36 MATERIALS SCIENCE↗

NASA NTRS2017

JIGSAW-GEO (1.0): Locally Orthogonal Staggered Unstructured Grid Generation for General Circulation Modelling on the Sphere

An algorithm for the generation of non-uniform, locally orthogonal staggered unstructured spheroidal grids is described. This technique is designed to generate very high-quality staggered VoronoiDelaunay meshes appropriate for general circulation modelling on the sphere, including applications to atmospheric simulation, ocean-modelling and numerical weather prediction. Using a recently developed Frontal-Delaunay refinement technique, a method for the construction of high-quality unstructured spheroidal Delaunay triangulations is introduced. A locally orthogonal polygonal grid, derived from the associated Voronoi diagram, is computed as the staggered dual. It is shown that use of the Frontal-Delaunay refinement technique allows for the generation of very high-quality unstructured triangulations, satisfying a priori bounds on element size and shape. Grid quality is further improved through the application of hill-climbing-type optimisation techniques. Overall, the algorithm is shown to produce grids with very high element quality and smooth grading characteristics, while imposing relatively low computational expense. A selection of uniform and non-uniform spheroidal grids appropriate for high-resolution, multi-scale general circulation modelling are presented. These grids are shown to satisfy the geometric constraints associated with contemporary unstructured C-grid-type finite-volume models, including the Model for Prediction Across Scales (MPAS-O). The use of user-defined mesh-spacing functions to generate smoothly graded, non-uniform grids for multi-resolution-type studies is discussed in detail.

geophysical fluid dynamics↗

NASA NTRS2018

Generalised Primal-Dual Grids for Unstructured Co-Volume Schemes

The generation of high-quality staggered unstructured grids is considered, leading to the development of a new optimisation-based strategy designed to construct weighted 'Regular-Power' tessellations appropriate for co-volume type numerical discretisation techniques. This new framework aims to extend the conventional Delaunay-Voronoi primal-dual structure; seeking to assemble generalised orthogonal tessellations with enhanced geometric quality. The construction of these grids is motivated by the desire to improve the performance and accuracy of numerical methods based on unstructured co-volume type schemes, including various staggered grid techniques for the simulation of fluid dynamics and hyperbolic transport. In this study, a new hybrid optimisation strategy is proposed; seeking to optimise the geometry, topology and weights associated with general, two-dimensional Regular-Power tessellations using a combination of gradient-ascent and energy-based techniques. The performance of this new method is tested experimentally, with a range of complex, multi-resolution primal-dual grids generated for various coastal and regional ocean modelling applications.

power diagrams↗