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.

At least 19 records

All sources NASA NTRS DOE OSTI

NASA NTRS2000

Quantifying Void Ratio in Granular Materials Using Voronoi Tessellation

Voronoi technique was used to calculate the local void ratio distribution of granular materials. It was implemented in an application-oriented image processing and analysis algorithm capable of extracting object edges, separating adjacent particles, obtaining the centroid of each particle, generating Voronoi polygons, and calculating the local void ratio. Details of the algorithm capabilities and features are presented. Verification calculations included performing manual digitization of synthetic images using Oda's method and Voronoi polygon system. The developed algorithm yielded very accurate measurements of the local void ratio distribution. Voronoi tessellation has the advantage, compared to Oda's method, of offering a well-defined polygon generation criterion that can be implemented in an algorithm to automatically calculate local void ratio of particulate materials.

Alshibli, Khalid A.↗

NASA NTRS1991

Quasi-periodic structures in the large-scale galaxy distribution and three-dimensional Voronoi tessellation

It has been suggested that the void and wall structure associated with the large-scale galaxy distribution might be qualitatively, or perhaps even physically, modeled by a Voronoi tessellation, and that such structure might account for the surprisingly regular, sharp peaks in the galaxy redshift distributions obtained from 'pencil beam' surveys. Taking cell wall crossings by random line segments to correspond to such redshift peaks, an exact expression is derived for the distribution of spacings of these intersections in a three-dimensional Voronoi tessellation. This result verifies that the spacings are nonrandom and quasi-periodic, qualitatively resembling the observed pattern, even though the cell wall structure is generated from randomly placed seeds. Finally, moments of the spacing distribution are used to show that apparently periodic samples, similar to those recently reported, represent only one to two sigma fluctuations in a Voronoi tesselation.

Ikeuchi, Satoru↗

DOE OSTI2020

Centroidal Voronoi tessellation based methods for optimal rain gauge location prediction

With more satellite and model precipitation data becoming available, new analytical methods are needed that can take advantage of emerging data patterns to make well informed predictions in many hydrological applications. Herein we propose a new strategy where we extract precipitation variability patterns and use correlation map to build the resulting density map that serves as an input to centroidal Voronoi tessellation construction that optimizes placement of precipitation gauges. We provide results of numerical experiments based on the data from the Alto-Adige region in Northern Italy and Oklahoma and compare them against actual gauge locations. This method provides an automated way for choosing new gauge locations and can be generalized to include physical constraints and to tackle other types of resource allocation problems.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2018

Robust Metric-Aligned Quad-Dominant Meshing Using L(sub p) Centroidal Voronoi Tessellation

We introduce a meshing algorithm that can be used to both generate and adapt meshes for bounded domains in an anisotropic manner. This is particularly beneficial when anisotropic flow features like shock waves or contact discontinuities are present in the computational domain. The algorithm presented in this paper is based upon meshing under the imposed Riemannian metric tensor, which controls the orientation and size of the mesh elements. In this way there is no need for user intervention to recognize these features. We demonstrate that the method indeed aligns the elements with the underlying metric and produces right-angled simplices that can be recombined into quadrilateral elements. The aim is to eventually incorporate this meshing strategy in the monolithic high-order spectral element solver that is currently being developed at NASA Ames. This paper has two main contributions: First, we demonstrate that we can generate quad-dominant metric-aligned meshes for bounded domains using a generalized form of L(sub p)-Centroidal Voronoi Tessellation (L(sub p)-CVT). Unlike previous works, we do not rely on a background mesh and discretize the bounded domain in a hierarchical way by first discretizing the boundaries and then the volume using the underlying metric. Second, we present an alternative for clipping the Voronoi cells on the boundary, which is common practice in CVT-based meshing algorithms, by reconstructing the Voronoi cells using the defined metric field. In this way we avoid the geometrical complexity of the clipping procedure and we show that we evaluate the energy and its gradient correctly. We show that the reconstruction of the computational domain is consistent with the Lloyds’ algorithm that is used to compute the L(sub p)-CVT.

Centroidal↗

NASA NTRS2006

Improving Neural Network Generalization Ability Using Outlier Analysis and Voronoi Tessellation

The data used in this study was obtained from the Sloan Digital Sky Survey (SDSS), which provides astronomers with what is currently the most extensive mapping of the universe, covering 25% of the sky and cataloging the spectral properties (e.g., luminosity, color, surface temperature) of over 100 million celestial objects. Images generated by the SDSS are collected through 5 filters named u, g, r, l, and z that have respective wavelengths of 3540, 4750, 6222, 7632, and 9049 A. By measuring the photometric redshifts of the aforementioned wavelengths of a galaxy, astronomers can ascertain the extent to which galaxy is receding from which the distance to the galaxy can be calculated. Data collected for a small select group of galaxies (approximately 30, 000) contains accurate measurements of the galaxies' redshifts, in addition to measurements of their spectral properties. The above dataset containing both redshift measurements as well as spectral properties of the selected galaxies served as the training set for the purposes of this study; the data set containing only the spectra properties of a separate group of galaxies served as the test set.

Ho, Michelle↗

DOE OSTI2021

The Area Localized Coupled Model for Analytical Mean Flow Prediction in Arbitrary Wind Farm Geometries

This work introduces the area localized coupled (ALC) model, which extends the applicability of approaches that couple classical wake superposition models and atmospheric boundary layer models to wind farms with arbitrary layouts. Coupling wake and top–down boundary layer models is particularly challenging since the latter requires averaging over planform areas associated with turbine-specific regions of the flow that need to be specified. The ALC model uses Voronoi tessellation to define this local area around each turbine. A top–down description of a developing internal boundary layer is then applied over Voronoi cells upstream of each turbine to estimate the local mean velocity profile. Coupling between the velocity at hub-height based on this localized top–down model and a wake model is achieved by enforcing a minimum least-square-error in mean velocity in each cell. The wake model in the present implementation takes into account variations in wind farm inflow velocity and represents the wake profile behind each turbine as a super-Gaussian function that smoothly transitions between a top-hat shape in the region immediately following the turbine to a Gaussian profile downstream. Detailed comparisons to large-eddy simulation (LES) data from two different wind farms demonstrate the efficacy of the model in accurately predicting both wind farm power output and local turbine hub-height velocity for different wind farm geometries. These validations using data generated from two different LES codes demonstrate the model's versatility with respect to capturing results from different simulation setups and wind farm configurations.

49 EE - Wind and Water Power Program - Wind (EE-4W↗

DOE OSTI2022

Voronoi cell analysis: The shapes of particle systems

Many physical systems can be studied as collections of particles embedded in space, often evolving in time. Natural questions arise concerning how to characterize these arrangements—are they ordered or disordered? If they are ordered, how are they ordered and what kinds of defects do they possess? Voronoi tessellations, originally introduced to study problems in pure mathematics, have become a powerful and versatile tool for analyzing countless problems in pure and applied physics. Here we explain the basics of Voronoi tessellations and the shapes that they produce and describe how they can be used to characterize many physical systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

An extension to V ORO ++ for multithreaded computation of Voronoi cells

V ORO ++ is a software library written in C++ for computing the Voronoi tessellation, a technique in computational geometry that is widely used for analyzing systems of particles. V ORO ++ was released in 2009 and is based on computing the Voronoi cell for each particle individually. Here, we take advantage of modern computer hardware, and extend the original serial version to allow for multithreaded computation of Voronoi cells via the OpenMP application programming interface. We test the performance of the code, and demonstrate that it can achieve parallel efficiencies greater than 95% in many cases. Further, the multithreaded extension follows standard OpenMP programming paradigms, allowing it to be incorporated into other programs. We provide an example of this using the VoroTop software library, performing a multithreaded Voronoi cell topology analysis of up to 102.4 million particles.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Meshless discretization of the discrete-ordinates transport equation with integration based on Voronoi cells

The time-dependent, gray, linear radiation transport equation is discretized using the meshless local Petrov-Galerkin method with reproducing kernels. The integration is performed using a Voronoi tessellation, which creates a partition of unity that only depends on the position and extent of the kernels. The resolution of the integration automatically follows the particles and requires no manual adjustment. The discretization includes streamline-upwind Petrov-Galerkin stabilization to prevent oscillations and improve numerical conditioning. The angular quadrature is selectively refineable to increase angular resolution in chosen directions. The time discretization is done using backward Euler. The transport solve for each direction and the solve for the scattering source are both done using Krylov iterative methods. The results indicate first-order convergence in time and second-order convergence in space for linear reproducing kernels.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Explainable machine learning reveals that local structural motifs encode the thermodynamic state across the CuZr metallic glass-forming range

Metallic glasses derive their properties from the statistics of local atomic motifs rather than from long-range order, yet a quantitative, chemistry-specific link between motif populations and the underlying glassy state has remained elusive. In this work we combine large-scale molecular dynamics, Voronoi tessellation, deep neural networks, and SHapley Additive exPlanations (SHAP) to identify which local structural motifs define the glassy state of Cu—Zr metallic glasses. A dataset of 17,180 atomistic configurations spanning ten compositions (Cu 20 Zr 80 –Cu 80 Zr 20 ) and four quench rates (10 9 –10 12 K/s) is used to train a feed-forward neural network that regresses temperature across the 50–2000 K liquid–supercooled–glass range, achieving a mean absolute error of 19.89 K and R 2 = 0.9974, confirming that the local structural state is faithfully encoded in motif-level structure. SHAP analysis then reveals that a tightly coupled near-icosahedral family of motifs (coordination numbers (CN) 11–13, including the full icosahedron 001200 and its single-atom-perturbation sibling 10930) collectively encodes the thermodynamic state of the system across the full glass-forming range. The CN = 11–13 ordered members carry negative SHAP values at high populations, tracking the most deeply-quenched configurations, while 10930 shows the reversed signature consistent with its role as a soft-spot host whose population shrinks as the icosahedral network deepens. The analysis demonstrates that explainable machine learning can isolate the minimal motif vocabulary defining the glassy state and recovers the near-icosahedral building blocks previously identified by data-driven analyses of Cu—Zr. The approach provides a general, chemistry-specific route for characterizing the structural state of disordered materials.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Excess free volume and structural properties of inert gas condensation synthesized nanoparticles based CuZr nanoglasses

Abstract Nanoglass (NG) as a new structure-tunable material has been investigated using both experiments and computational modeling. Experimentally, inert gas condensation (IGC) is commonly employed to prepare metallic glass (MG) nanoparticles that are consolidated using cold compression to generate an NG. In computational modeling, various methods have been used to generate NGs. However, due to the high computational cost involved, heretofore modeling investigations have not followed the experimental synthesis route. In this work, we use molecular dynamics simulations to generate an NG model by consolidating IGC-prepared Cu 64 Zr 36 nanoparticles following a workflow similar to that of experiments. The resulting structure is compared with those of NGs produced following two alternative procedures previously used: direct generation employing Voronoi tessellation and consolidation of spherical nanoparticles carved from an MG sample. We focus on the characterization of the excess free volume and the Voronoi polyhedral statistics in order to identify and quantify contrasting features of the glass-glass interfaces in the three NG samples prepared using distinct methods. Results indicate that glass-glass interfaces in IGC-based NGs are thicker and display higher structural contrast with their parent MG structure. Nanoparticle-based methods display excess free volume exceeding 4%, in agreement with experiments. IGC-prepared nanoparticles, which display Cu segregation to their surfaces, generate the highest glass-glass interface excess free volume levels and the largest relative interface volume with excess free volume higher than 3%. Voronoi polyhedral analysis indicates a sharp drop in the full icosahedral motif fraction in the glass-glass interfaces in nanoparticle-based NG as compared to their parent MG.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Multi-Sensor Optimal Motion Planning for Radiological Contamination Surveys by Using Prediction-Difference Maps

Distributed and networked mobile sensor platforms using unmanned aerial and/or ground vehicles to survey areas of interest offer a safer and more efficient method for radiological contamination mapping; however, most applications rely on uniformly sweeping of the area in a raster-type motion without utilizing the information available in a dynamic sense. We have developed a fully autonomous optimal motion planning procedure for networks with two or more mobile sensors. The procedure utilizes well-established concepts of Gaussian processes in combination with control laws based on centroidal Voronoi tessellations to achieve optimal next-iteration sensor movements. A new method of informing optimal motion planning is proposed, whereby the absolute difference between the prior and current full-map prediction, referred to as the prediction-difference map, is used as the spatial density function within each Voronoi cell, providing immediate and iterative feedback for dynamic use of available information. The Gaussian process regression model used to estimate the contamination in unvisited locations also provides prediction uncertainties, and can be used as a quantitative metric to assess the confidence in the calculated contamination map; these estimates and prediction uncertainties are unavailable for standard uniform survey routines as they can only produce maps in the vicinity of observed locations. We present through simulation the achievable performance gains from using this new method by directly comparing to a uniform survey method. Results show that using the prediction-difference maps to inform motion planning procedures offers a faster rate of producing an accurate and convergent map relative to a uniform survey route.

47 OTHER INSTRUMENTATION↗

DOE OSTI2021

Global field reconstruction from sparse sensors with Veronoi tessellation-assisted deep learning

Achieving accurate and robust global situational awareness of a complex time-evolving field from a limited number of sensors has been a longstanding challenge. This reconstruction problem is especially difficult when sensors are sparsely positioned in a seemingly random or unorganized manner, which is often encountered in a range of scientific and engineering problems. Moreover, these sensors can be in motion and can become online or offline over time. The key leverage in addressing this scientific issue is the wealth of data accumulated from the sensors. As a solution to this problem, we propose a data-driven spatial field recovery technique founded on a structured grid-based deep-learning approach for arbitrary positioned sensors of any numbers. It should be noted that the naïve use of machine learning becomes prohibitively expensive for global field reconstruction and is furthermore not adaptable to an arbitrary number of sensors. In the present work, we consider the use of Voronoi tessellation to obtain a structured-grid representation from sensor locations enabling the computationally tractable use of convolutional neural networks. One of the central features of the present method is its compatibility with deep-learning based super-resolution reconstruction techniques for structured sensor data that are established for image processing. The proposed reconstruction technique is demonstrated for unsteady wake flow, geophysical data, and three-dimensional turbulence. The current framework is able to handle an arbitrary number of moving sensors, and thereby overcomes a major limitation with existing reconstruction methods. The presented technique opens a new pathway towards the practical use of neural networks for real-time global field estimation.

Fukami, Kai↗

DOE OSTI2021

Comparing numerical accuracy and stability for different horizontal discretizations in MPAS-Ocean

This manuscript investigates the effectiveness of two possible horizontal discretizations for the global ocean model MPAS-Ocean, both applied to Spherical Centroidal Voronoi Tessellations (SCVTs). The first discretization is TRiSK, a C-grid, finite-volume method, that possesses many desirable mimetic properties, but has a low order accuracy. The second discretization was introduced for the first time by Peixoto (2016), and consists of modifications to the TRiSK scheme designed to achieve at least first-order accuracy in the L ∞ norm, with the loss of some of the mimetic properties. Tests on shallow-water and primitive-equation models show that the scheme due to Peixoto is indeed more accurate, but presents stability issues with respect to TRiSK. Here, TRiSK is indeed found to be often more stable in time and more robust with respect to errors in the geometric properties of the grid.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Meso-scale modelling of compressive fracture in concrete with irregularly shaped aggregates

This paper presents a meso-scale modelling framework to investigate the fracture process in concrete subjected to uniaxial and biaxial compression accounting for its mesostructural characteristics. 3D mesostructure of concrete consisting of coarse aggregates, mortar and interfacial transition zone between them was developed using an in-house code based on the Voronoi tessellation and splining method, which enables to generate the realistic-look aggregates with controllable structural features such as content, location, size and shape. Based on the generated 3D mesostructure, the concrete damage plasticity approach was employed to simulate the compressive fracture behaviour of concrete in terms of crack morphology and stress-strain response against the shape parameters of aggregate. Results indicate that the shape of aggregate has a negligible effect on compressive strength of concrete, which is highly associated with the random location and size distribution of aggregate. The aggregate irregularity has a significant influence on crack initiation and growth of concrete.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Effect of lattice orientation on compressive properties of selective laser sintered nylon lattice coupons

As AM lattices become more popular in medical devices, it is important to consider how lattice design parameters affect the mechanical integrity and performance of a device. This research investigated the effect of lattice orientation on the compressive mechanical response of five common lattice geometries: Hexagonal Honeycomb (Hex), Diamond, Voronoi Tessellation Method (VTM), sheet-based Gyroid, and Face-Centered-Cubic (FCC). Samples were tested at two relative densities and eight orientations, printed in Nylon 12 (PA2200) on an EOS P396. The mechanical response was compared to simulated finite element analysis (FEA) response for Hex and FCC lattices. Hex lattices displayed significant orientation dependance. Diamond lattices demonstrated some orientation dependance when at a lower relative density. Gyroid, FCC, and VTM all displayed minimal orientation dependance. Increasing lattice relative density appeared to reduce anisotropy among tested orientations of all lattice geometries except Hex. FEA simulations were able to produce a trend response that matched the mechanical response trends based on the investigated orientations. Furthermore, these results support the need for multi axial lattice mechanical test considerations. It also bolsters the need for lattice FEA models to have a sufficient comparator when identifying “worst-case” lattice test scenarios.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Achieving geometric accuracy in FFT-based micromechanical models using conformal grid

Owing to its efficiency, simplicity and robustness, the FFT-based method has become the standard for computation of mechanical fields in a heterogeneous periodic unit cell. One of the main disadvantages of the FFT-based method is the inaccurate representation of the initial microstructure on a regular grid of voxels, which can be alleviated through the use of distorted initial grids. Here, in this paper, a method for generation of distorted initial grids conforming to the microstructural features (e.g. straight/curved boundaries) is proposed. The method determines the positions of the grid nodes in the initial configuration by solving a system of springs connecting the nodes. Microstructures consisting of layers, Voronoi tessellation and circular/spherical inclusions are considered, and mechanical fields simulated using the FFT-based method. It is found that distorted initial grids, conforming to the microstructural features, lead to more accurate mechanical fields in comparison to the corresponding non-distorted initial grid solution. The effect of initial grid distortion on the convergence of the FFT-based method is analyzed and discussed.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Clustering of inertial particles in turbulent flow through a porous unit cell

We report direct numerical simulation is used to investigate effects of turbulent flow in the confined geometry of a face-centred cubic porous unit cell on the transport, clustering and deposition of fine particles at different Stokes numbers (St = 0.01, 0.1, 0.5, 1, 2) and at a pore Reynolds number of 500. Particles are advanced using one-way coupling and the collision of particles with pore walls is modelled as perfectly elastic with specular reflection. Tools for studying inertial particle dynamics and clustering developed for homogeneous flows are adapted to take into account the embedded, curved geometry of the pore walls. The pattern and dynamics of clustering are investigated using the volume change of Voronoi tesselation in time to analyse the divergence and convergence of the particles. Similar to the case of homogeneous, isotropic turbulence, the cluster formation is present at large volumes, while cluster destruction is prominent at small volumes and these effects are amplified with the Stokes number. However, unlike homogeneous, isotropic turbulence, the formation of a large number of very small volumes was observed at all Stokes numbers and attributed to the collision of particles with the pore wall. Multiscale wavelet analysis of the particle number density indicates that the peak of the energy density spectrum, representative of enhanced particle clustering, shifts towards larger scales with an increase in the Stokes number. Scale-dependent skewness and flatness quantify the intermittent void and cluster distribution, with cluster formation observed at small scales for all Stokes numbers, and void regions at large scales for large Stokes numbers.

42 ENGINEERING↗