SEARCH · Engineering Papers

Results for “Curved meshes”

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 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Vector diffraction analysis of reflector antennas with mesh surfaces

Reflector antennas with mesh surfaces are used extensively in many satellite and ground antenna systems. A strip-aperture modeling of commonly used mesh surfaces is presented which provides considerable versatility in characterizing the mesh cells. The mesh transmission coefficients are constructed using a Floquet-modal expansion in conjuction with two dominant aperture modes. To account for the mesh local coordinates, the Eulerian angle transformation is invoked to obtain the total induced current on the curved reflector surface. General formulas are presented to show how the solid surface induced current is modified due to the transmission through the mesh. The effects of a variety of mesh configurations on both the co-polar and cross-polar patterns of reflector antennas are studied by numerically evaluating the vector diffraction integral using the Jacobi-Bessel expansion. For some special cases, a comparison is made with the results of the commonly used wire-grid formulation. Many of the numerical data are tailored to the dimensions of a conceptually designed mesh deployable offset reflector of the land mobile satellite system (LMSS).

Rahmat-Samii, Y.↗

DOE OSTI2022

An extrapolation method to remove spurious stress concentration in micromechanical analyses of composites using pixel-based meshes

It is easier to generate a finite element mesh of complicated geometries with voxel-based meshes compared to conformal meshes. However, the local stresses in voxel-based meshes are higher owing to artificial stress concentration caused due to jagged lines that approximate a curve. Here, we propose a method that accounts for this anomaly which would lead to accurate prediction of maximum von Mises stress in voxel-based meshes. The method is based on the variation of reciprocal of stresses and associated numerical gradient. The prediction of the proposed method is compared against the results of a conformal mesh by considering four two-dimensional example problems. The accuracy of the proposed method is verified using different mesh densities, material properties, and integration schemes. The proposed method was also implemented on a conformal mesh to verify that it filters out just the spurious stress concentration and not the realistic stress concentration.

36 MATERIALS SCIENCE↗

NASA NTRSDate not supplied

Development of a Combined Cohesive and Virtual Crack-Closure Technique Approach to Represent R-Curves

Resistance curve (R-curve) effects due to fiber bridging, crack migration, and other blunting mechanisms are common in composite laminates. These mechanisms can dramatically increase the resistance to crack propagation but can be challenging to simulate. Delamination with R-curves can be analyzed using cohesive zone modeling (CZM) or the virtual crack closure technique (VCCT). Large fracture process zones can be simulated with CZM, but they require highly refined meshes. Coarser meshes can be used with VCCT, but this method is only applicable to small fracture process zones with R-curves defined as functions of position. Therefore, a technique with the computational efficiency of VCCT and the natural ability of cohesive elements to represent large fracture process zones is desirable. An approach is proposed that starts as a CZM in which the cohesive traction separation law (TSL) is separated into high-strength (HS) and low-strength (LS) components. The HS part, which is responsible for the mesh requirements of the CZM analysis, is replaced by VCCT. The combined method is evaluated by analyzing the response of double cantilevered beam specimens. Two specimen layup configurations, [0/90/90/0] 3s and [0/90/90/0] 9s , are evaluated. The R-curve response from each configuration is determined and applied to the constitutive properties of models built with CZM and with VCCT. The results and computational efficiency of the CZM, VCCT, and combined approaches are compared. The results are indicative that a combined cohesive/VCCT approach can enable progressive failure analyses to retain the computational efficiency of VCCT with the ability of the cohesive elements to capture R-curve effects.

Virtual Crack-Closure Technique↗

NASA NTRS1980

Development of glass encapsulation techniques for terrestrial photovoltaic arrays

Two parallel development programs for terrestrial solar cell module encapsulation and fabrication are reviewed, including the status of electrostatic bonding as an encapsulation technique. Current designs of electrostatically bonded modules are discussed, fabrication of which is now routine. The design of a high-efficiency module is presented, for which performance reliability features are discussed. This design is compatible with later generation changes such as the introduction of encapsulation by electrostatic bonding. Application of wire mesh contacts by electrostatic bonding has resulted in I-V curve fill factors of 0.74. Cell designs that would allow existing Pyrex glass to be used in electrostatically bonded modules are presented.

Younger, P. R.↗

NASA NTRS1995

SEEK: A FORTRAN optimization program using a feasible directions gradient search

This report describes the use of computer program 'SEEK' which works in conjunction with two user-written subroutines and an input data file to perform an optimization procedure on a user's problem. The optimization method uses a modified feasible directions gradient technique. SEEK is written in ANSI standard Fortran 77, has an object size of about 46K bytes, and can be used on a personal computer running DOS. This report describes the use of the program and discusses the optimizing method. The program use is illustrated with four example problems: a bushing design, a helical coil spring design, a gear mesh design, and a two-parameter Weibull life-reliability curve fit.

Savage, M.↗

DOE OSTI2026

Spatially Accelerated Winding Numbers for Curved Geometry

The generalized winding number (GWN) is a scalar field that supports robust containment queries on curved geometry, including non-watertight, overlapping, and nested boundary representations. While queries can be easily parallelized over samples, direct evaluation on parametric curves and surfaces remains costly for large and complex models. Fast, state-of-the-art GWN approaches leverage a spatial index to approximate the GWN, typically coupled with a Taylor expansion which approximates the GWN contribution for far clusters of geometric primitives. However, such methods operate only on discrete inputs such as triangle meshes and point clouds, and would introduce containment errors near boundaries if applied to curved input. We extend support for fast GWN evaluation over arbitrary collections of NURBS curves in 2D and trimmed NURBS patches in 3D via a Bounding Volume Hierarchy that stores efficiently precomputed moment data in the hierarchy nodes. When querying the hierarchy, approximations for far clusters are used alongside direct evaluation for nearby NURBS primitives, achieving sub-linear complexity while preserving the geometric features in the vicinity of the query point. Central to our performance improvements is an adaptive subdivision strategy for NURBS primitives during a preprocessing phase, creating better spatial partitions while retaining the same accuracy for containment decisions as a direct evaluation. We demonstrate the performance and accuracy of our approach across a large collection of 2D and 3D datasets.

Computer science↗

NASA NTRS1996

Generating Voronoi Diagrams for Curved Shapes with Divide-and-Conquer

Voronoi diagrams have been used in many practical applications including solid modeling, Numerical Control machining, Finite Element mesh generation, etc. Investigating the properties of these diagrams is an active research topic with many fruitful results. Computational techniques for Voronoi diagrams, however, have been concentrated on low order elements, for points, lines, polygons, quadratic curves and surfaces. This paper presents a divide-and-conquer scheme computing the diagrams for planar shapes bounded by closed curves commonly encountered in Computer Aided Design. These curves include analytical curves and splines. The algorithm first divides the boundary curve into sections, delimited by curve points with the minimal curvatures. Bisector branches for these sections are then generated and merged recursively to obtain the final diagram. A detailed example in the paper shows the steps of the generation and merging of the bisectors. A simple analysis of the complexity of the algorithm is also presented.

Chou, Jin J.↗

DOE OSTI2022

Hybrid Curve Fitting For Reducing Motion Commands In Object Construction

Existing slicing software for additive manufacturing typically requires a triangulated mesh as its input. Triangulated meshes are approximate representations of exact CAD models. Despite the loss in dimensional accuracy, triangulated meshes are used because they are computationally easier to cross-section and offset than the exact geometry in CAD format. When a triangulated object is prepared, the resulting machine instructions include only linear motion commands. Numerous modern motion controllers can move in arc and spline motions; however, the absence of slicing software that supports curvature prevents these commands from being leveraged. To address this limitation, this paper presents a method for the hybrid reconstruction of arcs and splines as a post-processing step to traditional slicing. This method can greatly reduce the number of motion commands required to construct an object by printing smooth curved surfaces. This concise representation of tool-pathing allows for more even extrusion and is computed without a major impact on slicing time.

Wade, Charles↗

DOE OSTI2022

June 2022 Permeability Measurements in P-Tunnel, Nevada National Security Site

This document describes gas permeability measurements made in June 2022 in P-Tunnel at the Nevada National Security Site (NNSS). Four boreholes (AC-1, GI-5, GI-6, and HF-1) were tested with a dual packer system. The pressure versus flow rate data that were collected are analyzed using a numerical simulation method that captures compressible gas flow in partially saturated porous media using a close approximation to the geometry of the test setup. The first part of this document begins with data collection and a description of the packer assembly followed by details on data acquisition, measurement instrumentation, and equipment followed by a description of pre-test system verification. Next, the implementation of testing in the field is described, including how specific locations were chosen for testing. The second part of this document describes the numerical simulations used to back-out permeability from the pressure versus flow rate data. First, the numerical model is described, including assumptions. Next the details of the mesh that represents the experiment are described. A method for generating type curves in pressure versus flow rate space is presented and finally the data are mapped to permeability. We conclude with a discussion of the results.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1980

Generation of boundary and boundary-layer fitting grids

The details of extended physical processes, such as the gas dynamic flow over an airfoil, the reactive flow through a combustor, or the electric field in a multi-contact transistor, are understood by solving the differential equations of a mathematical model of the process. The accuracy of finite difference methods for the numerical solution of the equations is increased if the underlying mesh fits the region boundaries and is closely spaced in regions where the solution is rapidly varying. Automatic methods for producing a satisfactorily adjusted mesh were developed for one dimensional problems. In one simple, effective scheme of this kind the unknown function and the distribution of mesh modes are found simultaneously, the nodes being placed so that they correspond to points uniformly spaced on the solution curve. In a two dimensional generalization, the nodes correspond to points equally spaced on the solution surface in two directions that are as nearly orthogonal as possible. Examples of such meshes are shown.

Ablow, C. M.↗

DOE OSTI2025

Replace Human Intelligence with Fast and Smart Geometric Reasoning and Graph Neural Network to Accelerate Next Gen ModSim Workflows

We present an agent-guided approach to CAD geometry decomposition that automates hex/hybrid meshing with graph neural networks (GNNs) to accelerate next-generation ModSim workflows. Our end-to-end pipeline (i) reduces 3D boundary-representation (B-Rep) models to a 2D chordal axis skeleton (CAT) and then to a 1D bipartite graph of surface and curve nodes, (ii) assigns per node labels as Cubit® WebCut actions, (iii) trains a multi-action GNN under supervised learning, and (iv) predicts five surface-node and three curve-node actions on out-of-distribution test geometries. Each graph node carries geometric, topological, and meshing attributes drawn from the B-Rep “skin” and CAT “skeleton,” with two-way mappings across 3D↔2D↔1D representations to maintain traceability back to 3D CAD. The supervised learning model exhibits stable convergence of the binary cross-entropy loss and achieves 98.7% accuracy on unseen lattice models. To operationalize decision-making, we rank predicted commands by geometric significance and prototyped the agent-guided workflow through the Cubit® Meshing PowerTool GUI. As a stretch goal, we explore reinforcement learning (RL) to reduce or remove label requirements and to learn policies for action sequences that maximize total reward (e.g., size of hex-meshable regions and resulting hex mesh quality). When all-hex meshing is not feasible, the agent assists in producing hybrid meshes—prioritizing hex in critical regions and transitioning to tetrahedral elements (tets) elsewhere—maintaining fidelity while ensuring robustness. The overarching objective is to replace manual, heuristics-based decomposition with data-driven, reproducible automation, cutting meshing turnaround time by orders of magnitude. We anticipate direct impact on simulation workflows through intelligent, scalable decomposition of complex CAD models into hex-meshable subdomains.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Closing the Loop between In Situ Stress Complexity and EGS Fracture Complexity

42 ENGINEERING↗

DOE OSTI2025

Global Gyrokinetic Simulations of Isotope Effects under Ambipolar Electric Fields and Advances Toward Whole-Volume Modeling

We review global gyrokinetic simulation studies on plasma transport in the Large Helical Device using XGC-S. XGC-S is an extended version of X-point Gyrokinetic Code for stellarators and has been progressively verified throughout the code development process. Verification tests of neoclassical transport successfully demonstrate the generation of an ambipolar electric field due to ripple-trapped particles. We perform quasi-linear analyses of the ion temperature gradient mode under the influence of the ambipolar electric field. The results reveal that the ambipolar electric field and the heavy hydrogen component in mixed isotope plasmas can lead to the favorable isotope effect observed in recent deuterium experiments. We also present recent efforts in code development toward whole-volume simulations, including the helical divertor region. A mesh generation scheme based on field-line tracing and the construction of curved surfaces perpendicular to the magnetic field would be promising for global field calculations in the whole-volume simulations.

Basic Plasma Phenomena and Gas Discharges↗

NASA NTRS1983

Multi-anode microchannel arrays - New detectors for imaging and spectroscopy in space

Consideration is given to the construction and operation of multi-anode microchannel array detector systems having formats as large as 256 x 1024 pixels. Such arrays are being developed for imaging and spectroscopy at soft X-ray, ultraviolet and visible wavelengths from balloons, sounding rockets and space probes. Both discrete-anode and coincidence-anode arrays are described. Two types of photocathode structures are evaluated: an opaque photocathode deposited directly on the curved-channel MCP and an activated cathode deposited on a proximity-focused mesh. Future work will include sensitivity optimization in the different wavelength regions and the development of detector tubes with semitransparent proximity-focused photocathodes.

Timothy, J. G.↗

NASA NTRS1999

Self-Avoiding Walks Over Adaptive Triangular Grids

Space-filling curves is a popular approach based on a geometric embedding for linearizing computational meshes. We present a new O(n log n) combinatorial algorithm for constructing a self avoiding walk through a two dimensional mesh containing n triangles. We show that for hierarchical adaptive meshes, the algorithm can be locally adapted and easily parallelized by taking advantage of the regularity of the refinement rules. The proposed approach should be very useful in the runtime partitioning and load balancing of adaptive unstructured grids.

Heber, Gerd↗

NASA NTRS1960

Measurement of Screen-Size Effects on Intensity, Scale, and Spectrum of Turbulence in a Free Subsonic Jet

The results are reported of hot-wire anemometer measurements of the fluctuating longitudinal component of the turbulent velocities in the mean flow downstream of screens in an air jet. These measurements have been analyzed by well-established techniques to give the influence of tile screen mesh size on the turbulent intensity, scale, and the power-spectral-density. The results show a linear dependence of the intensity upon the screen mesh size for locations within the central core of the air jet. The spectral-density curves show that the screens redistribute the turbulent energy from the low frequencies (<1000 cps) to the high frequencies (>1000 cps). The effects of the screens are overwhelmed in the mixing region of the jet flow by the turbulence levels existing there. The large pressure drops occurring across the screens reduce the velocity of the jet as compared to the jet without screens by approximately one-third for the velocity and range of mesh sizes investigated and reported in this report. The turbulence scale is a linear function of distance from the nozzle exit and is somewhat greater than comparable jets without screens.

Howard, Charles D.↗

DOE OSTI2023

Impact of artificial topological changes on flow and transport through fractured media due to mesh resolution

Abstract We performed a set of numerical simulations to characterize the interplay of fracture network topology, upscaling, and mesh refinement on flow and transport properties in fractured porous media. We generated a set of generic three-dimensional discrete fracture networks at various densities, where the radii of the fractures were sampled from a truncated power-law distribution, and whose parameters were loosely based on field site characterizations. We also considered five network densities, which were defined using a dimensionless version of density based on percolation theory. Once the networks were generated, we upscaled them into a single continuum model using the upscaled discrete fracture matrix model presented by Sweeney et al. (2019). We considered steady, isothermal pressure-driven flow through each domain and then simulated conservative, decaying, and adsorbing tracers using a pulse injection into the domain. For each simulation, we calculated the effective permeability and solute breakthrough curves as quantities of interest to compare between network realizations. We found that selecting a mesh resolution such that the global topology of the upscaled mesh matches the fracture network is essential. If the upscaled mesh has a connected pathway of fracture (higher permeability) cells but the fracture network does not, then the estimates for effective permeability and solute breakthrough will be incorrect. False connections cannot be eliminated entirely, but they can be managed by choosing appropriate mesh resolution and refinement for a given network. Adopting octree meshing to obtain sufficient levels of refinement leads to fewer computational cells (up to a 90% reduction in overall cell count) when compared to using a uniform resolution grid and can result in a more accurate continuum representation of the true fracture network.

58 GEOSCIENCES↗

NASA NTRS2001

A Parallel Cartesian Approach for External Aerodynamics of Vehicles with Complex Geometry

This workshop paper presents the current status in the development of a new approach for the solution of the Euler equations on Cartesian meshes with embedded boundaries in three dimensions on distributed and shared memory architectures. The approach uses adaptively refined Cartesian hexahedra to fill the computational domain. Where these cells intersect the geometry, they are cut by the boundary into arbitrarily shaped polyhedra which receive special treatment by the solver. The presentation documents a newly developed multilevel upwind solver based on a flexible domain-decomposition strategy. One novel aspect of the work is its use of space-filling curves (SFC) for memory efficient on-the-fly parallelization, dynamic re-partitioning and automatic coarse mesh generation. Within each subdomain the approach employs a variety reordering techniques so that relevant data are on the same page in memory permitting high-performance on cache-based processors. Details of the on-the-fly SFC based partitioning are presented as are construction rules for the automatic coarse mesh generation. After describing the approach, the paper uses model problems and 3- D configurations to both verify and validate the solver. The model problems demonstrate that second-order accuracy is maintained despite the presence of the irregular cut-cells in the mesh. In addition, it examines both parallel efficiency and convergence behavior. These investigations demonstrate a parallel speed-up in excess of 28 on 32 processors of an SGI Origin 2000 system and confirm that mesh partitioning has no effect on convergence behavior.

Aftosmis, M. J.↗