Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational topology”

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 577 records · Page 32

Topological or rotational non-Abelian gauge fields from Einstein-Skyrme holography

We report analytically known states at non-zero temperature which may serve as a powerful tool to reveal common topological and thermodynamic properties of systems ranging from the QCD phase diagram to topological phase transitions in condensed matter materials. In the holographically dual gravity theory, these are analytic solutions to a five-dimensional non-linear-sigma (Skyrme) model dynamically coupled to Einstein gravity. This theory is shown to be holographically dual to \( \mathcal{N} \) = 4 Super-Yang-Mills theory coupled to an SU(2)-current. All solutions are fully backreacted asymptotically Anti-de Sitter (AdS) black branes or holes. One family of global AdS black hole solutions contains non-Abelian gauge field configurations with positive integer Chern numbers and finite energy density. Larger Chern numbers increase the Hawking-Page transition temperature. In the holographically dual field theory this indicates a significant effect on the deconfinement phase transition. Black holes with one Hawking temperature can have distinct Chern numbers, potentially enabling topological transitions. A second family of analytic solutions, rotating black branes, is also provided. These rotating solutions induce states with propagating charge density waves in the dual field theory. We compute the Hawking temperature, entropy density, angular velocity and free energy for these black holes/branes. These correspond to thermodynamic data in the dual field theory. For these states the energy-momentum tensor, (non-)conserved current, and topological charge are interpreted.

79 ASTRONOMY AND ASTROPHYSICS↗

Computing a Stable Distance on Merge Trees

Distances on merge trees facilitate visual comparison of collections of scalar fields. Two desirable properties for these distances to exhibit are 1) the ability to discern between scalar fields which other, less complex topological summaries cannot and 2) to still be robust to perturbations in the dataset. The combination of these two properties, known respectively as stability and discriminativity, has led to theoretical distances which are either thought to be or shown to be computationally complex and thus their implementations have been scarce. In order to design similarity measures on merge trees which are computationally feasible for more complex merge trees, many researchers have elected to loosen the restrictions on at least one of these two properties. The question still remains, however, if there are practical situations where trading these desirable properties is necessary. Here we construct a distance between merge trees which is designed to retain both discriminativity and stability. While our approach can be expensive for large merge trees, we illustrate its use in a setting where the number of nodes is small. This setting can be made more practical since we also provide a proof that persistence simplification increases the outputted distance by at most half of the simplified value. As a result, we demonstrate our distance measure on applications in shape comparison and on detection of periodicity in the von Kármán vortex street.

97 MATHEMATICS AND COMPUTING↗

Lessons Learned from Numerical Simulations of the F-16XL Aircraft at Flight Conditions

Nine groups participating in the Cranked Arrow Wing Aerodynamics Project International (CAWAPI) project have contributed steady and unsteady viscous simulations of a full-scale, semi-span model of the F-16XL aircraft. Three different categories of flight Reynolds/Mach number combinations were computed and compared with flight-test measurements for the purpose of code validation and improved understanding of the flight physics. Steady-state simulations are done with several turbulence models of different complexity with no topology information required and which overcome Boussinesq-assumption problems in vortical flows. Detached-eddy simulation (DES) and its successor delayed detached-eddy simulation (DDES) have been used to compute the time accurate flow development. Common structured and unstructured grids as well as individually-adapted unstructured grids were used. Although discrepancies are observed in the comparisons, overall reasonable agreement is demonstrated for surface pressure distribution, local skin friction and boundary velocity profiles at subsonic speeds. The physical modeling, steady or unsteady, and the grid resolution both contribute to the discrepancies observed in the comparisons with flight data, but at this time it cannot be determined how much each part contributes to the whole. Overall it can be said that the technology readiness of CFD-simulation technology for the study of vehicle performance has matured since 2001 such that it can be used today with a reasonable level of confidence for complex configurations.

Rizzi, Arthur↗

What was Learned from Numerical Simulations of F-16XL (CAWAPI) at Flight Conditions

Nine groups participating in the CAWAPI project have contributed steady and unsteady viscous simulations of a full-scale, semi-span model of the F-16XL aircraft. Three different categories of flight Reynolds/Mach number combinations were computed and compared with flight-test measurements for the purpose of code validation and improved understanding of the flight physics. Steady-state simulations are done with several turbulence models of different complexity with no topology information required and which overcome Boussinesq-assumption problems in vortical flows. Detached-eddy simulation (DES) and its successor delayed detached-eddy simulation (DDES) have been used to compute the time accurate flow development. Common structured and unstructured grids as well as individually-adapted unstructured grids were used. Although discrepancies are observed in the comparisons, overall reasonable agreement is demonstrated for surface pressure distribution, local skin friction and boundary velocity profiles at subsonic speeds. The physical modeling, be it steady or unsteady flow, and the grid resolution both contribute to the discrepancies observed in the comparisons with flight data, but at this time it cannot be determined how much each part contributes to the whole. Overall it can be said that the technology readiness of CFD-simulation technology for the study of vehicle performance has matured since 2001 such that it can be used today with a reasonable level of confidence for complex configurations.

Rizzi, Arthur↗

Lagrangian–Eulerian multidensity topology optimization with the material point method

Abstract In this paper, a hybrid Lagrangian–Eulerian topology optimization (LETO) method is proposed to solve the elastic force equilibrium with the Material Point Method (MPM). LETO transfers density information from freely movable Lagrangian carrier particles to a fixed set of Eulerian quadrature points. This transfer is based on a smooth radial kernel involved in the compliance objective to avoid the artificial checkerboard pattern. The quadrature points act as MPM particles embedded in a lower‐resolution grid and enable a subcell multidensity resolution of intricate structures with a reduced computational cost. A quadrature‐level connectivity graph‐based method is adopted to avoid the artificial checkerboard issues commonly existing in multiresolution topology optimization methods. Numerical experiments are provided to demonstrate the efficacy of the proposed approach.

Li, Yue↗

Computer generation of symbolic network functions - A new theory and implementation.

A new method is presented for obtaining network functions in which some, none, or all of the network elements are represented by symbolic parameters (i.e., symbolic network functions). Unlike the topological tree enumeration or signal flow graph methods generally used to derive symbolic network functions, the proposed procedure employs fast, efficient, numerical-type algorithms to determine the contribution of those network branches that are not represented by symbolic parameters. A computer program called NAPPE (for Network Analysis Program using Parameter Extractions) and incorporating all of the concepts discussed has been written. Several examples illustrating the usefulness and efficiency of NAPPE are presented.

Alderson, G. E.↗

Topology optimization of additively manufactured fluidic components free of internal support structures

This report integrates projection-based approaches for implementing overhang constraints with fluid topology optimization to design additively manufactured fluidic components that do not require internal support structures. Internal support structures are challenging, potentially impossible to remove in complex fluid networks, yet would degrade fluidic performance if left within the part. This presents a major challenge in coupling topology optimization with additive manufacturing in the design of fluid components. The proposed approach leverages past work in overhang constraints for solid mechanics, including projection formulations and adjoint sensitivity analysis, and considers incompressible Navier–Stokes equations for laminar fluid flow. The approach is demonstrated on 2D and 3D fluid topology optimization problems to minimize pressure drop and constrain mass flow rates in pipes and manifolds. Resulting designs are crisp, logical, satisfy performance constraints, and satisfy overhang constraints, eliminating the need for internal support structures.

42 ENGINEERING↗

Minimum feature size control in level set topology optimization via density fields

A level set topology optimization approach that uses an auxiliary density field to nucleate holes during the optimization process and achieves minimum feature size control in optimized designs is explored. The level set field determines the solid-void interface and the density field describes the distribution of a fictitious porous material using the solid isotropic material with penalization. These fields are governed by two sets of independent optimization variables which are initially coupled using a penalty for hole nucleation. The strength of the density field penalization and projection is gradually increased during the optimization process to promote a 0-1 density distribution. In addition, a second penalty regulates the evolution of the density field in the void phase. The treatment of the density field combined with the second penalty mitigate the appearance of small design features. The minimum feature size of optimized designs is controlled by the radius of the linear filter applied to the density optimization variables. The structural response is predicted by the extended finite element method, the sensitivities by the adjoint method, and the optimization variables are updated by a gradient-based optimization algorithm. Numerical examples investigate the robustness of this approach with respect to algorithmic parameters and mesh refinement. The results show the applicability of the combined density level set topology optimization approach for both optimal hole nucleation and for minimum feature size control in 2D and 3D. This comes, however, at the cost of a more complex problem formulation and additional computational cost due to an increased number of optimization variables.

42 ENGINEERING↗

Parallelized solvers for heat conduction formulations

Based on multilevel partitioning, this paper develops a structural parallelizable solution methodology that enables a significant reduction in computational effort and memory requirements for very large scale linear and nonlinear steady and transient thermal (heat conduction) models. Due to the generality of the formulation of the scheme, both finite element and finite difference simulations can be treated. Diverse model topologies can thus be handled, including both simply and multiply connected (branched/perforated) geometries. To verify the methodology, analytical and numerical benchmark trends are verified in both sequential and parallel computer environments.

Padovan, Joe↗

AdS 3 gravity and random CFT

We compute the path integral of three-dimensional gravity with negative cosmological constant on spaces which are topologically a torus times an interval. These are Euclidean wormholes, which smoothly interpolate between two asymptotically Euclidean AdS 3 regions with torus boundary. From our results we obtain the spectral correlations between BTZ black hole microstates near threshold, as well as extract the spectral form factor at fixed momentum, which has linear growth in time with small fluctuations around it. The low-energy limit of these correlations is precisely that of a double-scaled random matrix ensemble with Virasoro symmetry. Our findings suggest that if pure three-dimensional gravity has a holographic dual, then the dual is an ensemble which generalizes random matrix theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

AdS 3 wormholes from a modular bootstrap

In recent work we computed the path integral of three-dimensional gravity with negative cosmological constant on spaces which are topologically a torus times an interval. Here we employ a modular bootstrap to show that the amplitude is completely fixed by consistency conditions and a few basic inputs from gravity. This bootstrap is notably for an ensemble of CFTs, rather than for a single instance. We also compare the 3d gravity result with the Narain ensemble. The former is well-approximated at low temperature by a random matrix theory ansatz, and we conjecture that this behavior is generic for an ensemble of CFTs at large central charge with a chaotic spectrum of heavy operators.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Fractal breakage of porous carbonate sand particles: Microstructures and mechanisms

In situ, three-dimensional (3D) characterizations of particle breakage in porous carbonate sands are presented, for the first time, with synchrotron-based micro computed tomography. Evolution of grain-scale characteristics are identified and quantified via elaborate image processing and topology analyses. The sequential 3D images reveal distinctly different fracture mechanisms for carbonate sands from silica sands. The angular shape of carbonate sand particles facilitates bending fracture, and particles with a lower sphericity and a higher porosity are more prone to break. 3D crack networks extracted from fractured particles imply considerable cleavage along initial pores. The fractal dimension of crack networks increases with external loading due to crack branching via cleavage. The resultant fragment size distribution also appears fractal and the fractal feature is valid down to the breakage limit of calcium carbonate. Crack propagation along the initial pores reduces the energy barrier for particle breakage and thus fracture strength of particles

Carbonate sands↗

Automatic Lane-Level Road Network Extraction from Aerial Imagery for Transportation Digital Twins

Accurate road networks are essential for credible traffic microsimulation and transportation digital twins, yet high-definition maps are often difficult to obtain due to limited availability, high cost, or proprietary restrictions. Some build networks from crowdsourced data, such as OpenStreetMap, but these sources often contain geometric and semantic inconsistencies. Others create networks manually, a process that is labor-intensive and difficult to scale. To address these limitations, this work presents an end-to-end pipeline that automatically extracts georeferenced, lane-level road networks from publicly available high-resolution satellite imagery and converts them into simulation-ready assets. The developed end-to-end pipeline has three primary modules: (1) A computer-vision-based module first detects directed lane geometries and intersection layouts. (2) A heuristic-based topology construction module then identifies approach and exit legs and establishes conflict-free lane-to-lane connections. (3) Finally, an automatic simulation-building module converts the extracted network into standard formats, e.g., OpenDRIVE, and generates routable SUMO networks. The framework supports both complete network construction from scratch and local-scale refinement of existing networks through lane-count correction, transition recovery, and geometric regularization. The proposed pipeline provides a practical pathway to generate traffic simulation networks from satellite imagery, significantly reducing manual reconstruction effort and enabling scalable, continuously updated transportation digital twins.

Guo, Hetian [University of Georgia, Athens] (ORCID↗

Computer analysis of general linear networks using digraphs.

Investigation of the application of digraphs in analyzing general electronic networks, and development of a computer program based on a particular digraph method developed by Chen. The Chen digraph method is a topological method for solution of networks and serves as a shortcut when hand calculations are required. The advantage offered by this method of analysis is that the results are in symbolic form. It is limited, however, by the size of network that may be handled. Usually hand calculations become too tedious for networks larger than about five nodes, depending on how many elements the network contains. Direct determinant expansion for a five-node network is a very tedious process also.

Mcclenahan, J. O.↗

State criminal justice telecommunications (STACOM). Volume 4: Network design software user's guide

A user's guide to the network design program is presented. The program is written in FORTRAN V and implemented on a UNIVAC 1108 computer under the EXEC-8 operating system which enables the user to construct least-cost network topologies for criminal justice digital telecommunications networks. A complete description of program features, inputs, processing logic, and outputs is presented, and a sample run and a program listing are included.

Lee, J. J.↗

Global simulations of the three-dimensional magnetosphere

Global three-dimensional computer simulations of the magnetosphere using a particle MHD code, reproduce the steady-state Dungey magnetospheric topology in three dimensions. The formation of a compression zone downstream of the tail neutral line that is probably bounded by wake shocks is observed. This compression zone changes its cross-section with distance downstream.

Leboeuf, J. N.↗

Automatic algebraic coordinate generation

A computer software system has been developed to automatically generate two-dimensional coordinates from algebraic transformations. For topologically complex regions, a smooth assembly of the transformations can be used to automatically produce a composite mesh where a general gridded format is retained. The algebraic mesh generation system consists of a collection of operator subroutines which are applied to an established data structure and which automatically perform the necessary parts of mesh construction from a sequence of multisurface transformations. The system operators are discussed, taking into account the data base, the order of application, direct surface generators, geometric surface operators, surface generators from existing surfaces, transverse operators, mesh operators, assembly operators, and data visualization operators. Attention is given to applications related to airfoils.

Eiseman, P. R.↗

Magnetohydrodynamic simulation of the coronal transient associated with the solar limb flare of 1980, June 29, 18:21 UT

The spatial and temporal behavior of excess and depleted density regions of the coronal transient accompanying the west limb solar flare of 18:21 UT on June 29, 1980, is modeled mathematically on the basis of SMM-X-ray-polychrometer data, and the results are compared to observations made with the radio spectrograph at Harvard Radio Astronomy Station at Fort Davis and with the NCAR/HAO Mark III K-coronameter at Mauna Loa. Input data for the model include T(max) = about 20 x 10 to the 6th K, n(max) = about 4 x 10 to the 11th/cu cm, and an assumed ejection velocity of 200 km/sec. Computations using an improved 2D nonplane MHD model are carried out for locally open and closed magnetic topologies. The spatially wide, large-amplitude, temporarily steepened MHD wave predicted by the model for both magnetic topologies is shown to agree well with the observations, except that the predicted density enhancement exceeded observed values by at least 50 percent for the closed field and by a factor of 3 for the open field. This discrepancy is seen as an indication that the mass emitting soft X-rays was confined in closed-field regions near the sun during the obervation period.

Wu, S. T.↗