Search NASASearch

SEARCH · Search NASA

Results for “meshing”

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 163 records · Page 9

Exploring Unstructured Mesh Adaptation for Hybrid Reynolds-Averaged Navier–Stokes/Large Eddy Simulation

Mesh adaptation methods for the Reynolds-averaged Navier–Stokes (RANS) equations are rapidly maturing and beginning to impact the design of aerospace vehicles. RANS turbulence modeling improvements have slowed and may stagnate. Wall-modeled large eddy simulation (LES) and hybrid RANS/LES (HRLES) may provide an improved modeling capability but require specialized expertise to construct appropriate meshes and are considered be too ex-pensive for routine practical use. The realization of the CFD Vision 2030 Study includes improving geometry linkage, mesh generation/adaptation, and turbulence modeling/resolving methods for automated management of errors and uncertainties of physics-based, predictive modeling that can set the stage for ensuring a vehicle is in compliance with a regulation or specification (i.e., certification or qualification by analysis). An exploration of mesh adaptation for HRLES is performed to document synergies and challenges between mesh adaptation and HRLES. Vortex breakdown over a delta wing is examined to show the improvement of HRLES over RANS turbulence modeling approaches. A high lift configuration is shown to demonstrate complex geometry capability. Research and development opportunities are identified to advocate for continuing investments that may allow HRLES to enter routine practical use as a tool for aerospace vehicle analysis and design.

Michael A Park

Mesh Effects on Flow Solutions for a 2-D Multi-Element Airfoil Using Structured Overset Methods

Current best-practices for mesh generation in a structured overset framework are evaluated to determine their effects on the accuracy and convergence of the flow solution on a two-dimensional sectional cut of the High-Lift Common Research Model (HLCRM) wing. Mesh parameters are varied independently and correlated with force and moment coefficients from the flow solutions to identify the trends and sensitivities of the parameter variations, and the characteristics of a best practice mesh. Grid convergence and mesh adaption studies are also performed to determine the minimum mesh quality needed in order for solution-based adaption to reach an accurate flow solution.

T^3

Nearfield Anisotropic Mesh Adaptivity for the Third AIAA Sonic Boom Workshop

The Third AIAA Sonic Boom Workshop provides a unique opportunity to verify nearfield Computational Fluid Dynamics (CFD) tools. The workshop gathered nearfield CFD pressures from families of Mach-aligned, manually-tailored meshes on the C608 Low Boom Flight Test Demonstrator and the shock-plume interaction wind tunnel model. Here two classical adaptive strategies, multiscale and goal-oriented, are compared to the results obtained on tailored grids. The multiscale strategy is implemented independently in two separate toolsets. Details of the adapted mesh density are shown that resolve the complex interaction of boundary layers, shocks, expansions, and vortical structures. Mesh convergence of nearfield pressure signatures and their integral is shown. These detailed comparisons between adapted and manually-tailored meshes and independent implementations of mesh adaptation demonstrate the readiness of these methods for controlling the discretization error of nearfield sonic boom predication.

Julien Vanharen

An Automated Meshing Framework for Progressive Damage Analysis of Fabrics Using CompDam

A meshing and model construction framework for progressive damage analysis of fabric composites is proposed. A mesoscale idealization is used such that the fiber tows are considered homogeneous and transversely isotropic continua, and the matrix is a homogeneous isotropic continuum. By using this idealization, relatively mature damage modeling techniques developed for laminates of unidirectional plies are transferrable to fabrics. Cohesive interface elements and a continuum damage mechanics (CDM) model implemented in the software CompDam are used to account for the key damage modes including cracking at the tow/matrix interface, transverse matrix cracking in the tows, and fiber breakage in the tows. Matrix degradation is accounted for by plasticity. The modeling approach is demonstrated for a 3D orthogonal woven fabric. The main novelty of the proposed model construction framework is in its flexibility resulting from a formulation that is agnostic to the fabric architecture and implementation using mesh based operations. The algorithm does not make any assumptions about the tow sizes, shape, path, or proximity to other tows, thus it is applicable to a broad range of fabric architectures. The matrix region mesh surrounding the tows is constructed from domain boundary and tow surface meshes (i.e. through a series of mesh-based operations) such that the model domain can be arbitrarily shaped.

Andrew Cole Bergan

Learning Topological Operations on Meshes with Application to Block Decomposition of Polygons

We present a learning based framework for mesh quality improvement on unstructured triangular and quadrilateral meshes. Our model learns to improve mesh quality according to a prescribed objective function purely via self-play reinforcement learning with no prior heuristics. The actions performed on the mesh are standard local and global element operations. The goal is to minimize the deviation of the node degrees from their ideal values, which in the case of interior vertices leads to a minimization of irregular nodes.

97 MATHEMATICS AND COMPUTING

Computing material volume fractions on a superimposed mesh as applied to Monte Carlo particle transport simulations

Here, we present a newly implemented ray tracing algorithm in OpenMC for efficiently computing material volume fractions on superimposed meshes in complex geometries. By firing rays along each coordinate direction through the geometry, the approach accumulates track-length data in each mesh element, thereby determining the fractional composition of each material. Scaling studies on three different models—a random tetrahedra configuration, the Frascati Neutron Generator ITER dose rate benchmark, and a stellarator design—show excellent parallel performance, with nearly linear speedup on modern multi-threaded and distributed-memory systems. An analysis of the residual error relative to high-resolution reference solutions demonstrated that under optimal conditions it decreases as 1/R, where R is the number of rays fired, making it straightforward to achieve user-prescribed accuracy. This new functionality enables practical, mesh-based approaches for detailed nuclear analyses in production Monte Carlo workflows without resorting to expensive, fully conformal or unstructured meshing.

Monte Carlo

Application of mesh refinement to relativistic magnetic reconnection

During relativistic magnetic reconnection, antiparallel magnetic fields undergo a rapid change in topology, releasing a large amount of energy in the form of non-thermal particle acceleration. This work explores the application of mesh refinement to 2D reconnection simulations to efficiently model the inherent disparity in length-scales. We have systematically investigated the effects of mesh refinement and determined necessary modifications to the algorithm required to mitigate non-physical artifacts at the coarse–fine interface. We have used the ultrahigh-order pseudo-spectral analytical time-domain Maxwell solver to analyze how its use can mitigate the numerical dispersion that occurs with the finite-difference time-domain (or “Yee”) method. Absorbing layers are introduced at the coarse–fine interface to eliminate spurious effects that occur with mesh refinement. We also study how damping the electromagnetic fields and current density in the absorbing layer can help prevent the non-physical accumulation of charge and current density at the coarse–fine interface. Using a mesh refinement ratio of 8 for two-dimensional magnetic reconnection simulations, we obtained good agreement with the high-resolution baseline simulation, using only 36% of the macroparticles and 71% of the node-hours needed for the baseline. The methods presented here are especially applicable to 3D systems where higher memory savings are expected than in 2D, enabling comprehensive, computationally efficient 3D reconnection studies in the future.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Diffusion Limit–Preserving Lumped DFEMs on AMR Meshes

Here, we present sweep-compatible, novel upwinding recipes for the bilinear discontinuous (BLD) finite element method (FEM) that allows lumped BLD to be used on adaptive mesh refinement (AMR) meshes for thick transport applications without adding additional degrees of freedom at hanging nodes that exist on refinement boundaries. We analyze the properties of the upwinding and lumping that are needed for BLD to get the thick diffusion limit on such meshes, present results demonstrating locking with the wrong recipe, and present results showing error convergence and robustness properties for two diffusive problems on a variety of AMR meshes.

42 ENGINEERING

Improvements to the LAURA Mesh Adaptation Algorithm

This paper describes recent improvements to the LAURA code that removes the requirementof unbroken body-normal lines when using the automated mesh adaptation algorithm inLAURA responsible for mesh alignment with shock and adjusting near-wall mesh spacing. Fornon-simple geometry, restricting block topologies such that blocks span the domain from thesurface wall boundaries to domain inflow boundary can be difficult or impossible. Removing thisrequirement enables novel block topologies that simplify the meshing workflow for LAURA users.The new adaptation algorithm implementation in LAURA is described and several examples ofthe adaptation are presented to demonstrate the efficacy of the adaptation algorithm.

LAURA

Improvements to the LAURA Mesh Adaptation Algorithm

This paper describes recent improvements to the LAURA code that removes the requirementof unbroken body-normal lines when using the automated mesh adaptation algorithm inLAURA responsible for mesh alignment with shock and adjusting near-wall mesh spacing. Fornon-simple geometry, restricting block topologies such that blocks span the domain from thesurface wall boundaries to domain inflow boundary can be difficult or impossible. Removing thisrequirement enables novel block topologies that simplify the meshing workflow for LAURA users.The new adaptation algorithm implementation in LAURA is described and several examples ofthe adaptation are presented to demonstrate the efficacy of the adaptation algorithm.

LAURA

Thermomechanics coupling to Monte Carlo particle transport on unstructured mesh geometries using Cardinal

Geometry deformation due to thermal expansion influences neutron transport in many systems. Studying this phenomenon involves coupling models for neutronics, thermal hydraulics, and solid mechanics. To enable high fidelity modeling of these coupled physics, new capabilities were introduced in Cardinal, coupling OpenMC Monte Carlo particle transport models with MOOSE thermomechanical physics on unstructured moving-mesh geometries. In this work, we present a fully open-source capability leveraging on-the-fly mesh skinning to automatically regenerate OpenMC geometry, which allows multiphysics feedback from temperature, density, and geometry changes. The new capability is verified using an analytic benchmark slab problem, which couples S 2 neutron transport with thermal conduction, convective boundary conditions, Doppler-broadened cross sections, and nonlinear thermal expansion effects along the heated slab. Cardinal reproduces the analytic solutions for the neutron flux, heating, k eff , and temperature with demonstrated convergence in various error terms including mesh resolution and cross section temperature library spacing. For the nominal benchmark conditions and with a fine mesh, maximum relative errors for neutron flux, temperature, and heating are lower than 1%, while errors in integral quantities such as k eff and slab length are within 1 pcm and 48 µm, respectively. This work (i) presents a new numerical approach to thermomechanics coupling with OpenMC models, (ii) is the first (to our knowledge) to utilize a mechanical partial differential equation (PDE) solution to solve the (Griesheimer and Kooreman, 2022) analytic benchmark, and (iii) develops this verified capability within an open-source package.

97 - MATHEMATICS AND COMPUTING

Instantaneous mesh load factor ( K γ ) measurements in a wind turbine gearbox using fiber-optic strain sensors

The mesh load factor, K γ , describes how loads are shared between planet gears and has become one of the key design challenges in modern wind turbine gearboxes. Planet load sharing directly impacts tooth root stresses, a critical driver of torque density and gearbox reliability. Experimental evaluation of K γ is typically performed from sun gear tooth root strain gauge measurements, which are complex. Furthermore, such measurements can only provide an average value of load sharing. The present study describes an alternative method to evaluate the mesh load factor in wind turbine gearboxes based on fiber-optic strain sensors installed on the outer surface of the fixed ring gear. We present the results of an extensive measurement campaign to evaluate this novel sensing solution installed on the input planetary stage of a 2-MW wind turbine gearbox at the National Renewable Energy Laboratory's Flatirons Campus (Colorado, USA). The number of strain sensors on the ring gear was selected as an integer multiple of the number of planets, which has enabled an instantaneous evaluation of the mesh load factor. The effect of operating conditions on the planet load-sharing behavior of the gearbox has been investigated. The mesh load factor measured for operating conditions close to rated was below 1.05, well below IEC 61400-4 standard requirements.

17 WIND ENERGY

automesh: Automatic mesh generation in Rust

automesh is an open-source Rust software program that uses a segmentation, typically generated from a 3D image stack, to create a finite element mesh, composed either of hexahedral (volumetric) or triangular (isosurface) elements. automesh converts between segmentation formats (.npy, .spn) and mesh formats (.exo, .inp, .mesh, .stl, .vtk). automesh can defeature voxel domains, apply Laplacian and Taubin smoothing, and output mesh quality metrics. automesh uses an internal octree for fast performance.

Hovey, Chad Brian [Sandia National Laboratories (S

Grain2Mesh: Mesh Generation for Grain-Scale Nonlinear Elasticity Modeling

The nonlinear hysteretic behavior of rocks under cyclic loading is a crucial area of study in geomechanics. The macroscopic response of a variety of materials has been found to be contingent upon the behavior of the micro-scale structure. This project aims to develop a functional and maintainable software package for generating a multi-phase numerical mesh and accompanying simulation files for finite element modeling used in computational mechanics solvers. Meshes generated from images often lack key preprocessing that reduces noise and prevents mesh element distortion that can increase computational cost. By incorporating user feedback throughout, grain2mesh ensures a high-fidelity mesh that can be used to model grain-scale interactions such as shearing, crack propagation, and interfacial material contrast. Scientific applications of this software include material fracturing, stress-strain analysis for natural and engineered materials, and nonlinear meso-scale analysis.

54 ENVIRONMENTAL SCIENCES

A mesh gradient technique for numerical optimization

A class of successive-improvement optimization methods in which directions of descent are defined in the state space along each trial trajectory are considered. The given problem is first decomposed into two discrete levels by imposing mesh points. Level 1 consists of running optimal subarcs between each successive pair of mesh points. For normal systems, these optimal two-point boundary value problems can be solved by following a routine prescription if the mesh spacing is sufficiently close. A spacing criterion is given. Under appropriate conditions, the criterion value depends only on the coordinates of the mesh points, and its gradient with respect to those coordinates may be defined by interpreting the adjoint variables as partial derivatives of the criterion value function. In level 2, the gradient data is used to generate improvement steps or search directions in the state space which satisfy the boundary values and constraints of the given problem.

Willis, E. A., Jr.

Some mesh generation requirements and methods

Discretized solution algorithms, which find solutions of field equations in a two or three dimensional field, generally use meshes which are fitted to the field boundary to allow convenient formulation of boundary conditions there. A mesh is defined to be the image of a rectangular grid in computational space under a mesh mapping which maps computational space into physical space. It is not necessary that all of computational space be mapped onto the region of interest in physical space. Parts of it can be excised to give a better fit to the boundary. Many different excisions can be made to fit a single boundary; the choice depends on the mesh arrangement desired in the field.

Dickson, L. J.

Measurement of losses of mesh membrane material for reflector applications with an S-band radiometer

The spatial resolution, the frequency of revisits, and the accuracy requirements for the detection of Earth surface parameters from space force the satellite designers to consider large space structures for microwaves. Some Earth surface parameters are only detectable with radiometers. Because the large reflectors should be much lighter than solid and deployable reflectors, mesh membrane material is considered to be used as the reflector. It is essential to determine the emissivity of the reflecting material even when close to zero to estimate the accuracy of the radiometer measurement. An existing S-band radiometer was used to determine the loss of a 1 sq m pretensioned mesh first in a field experiment and later in a well shielded laboratory set-up. The models for retrieving the emissivity of the mesh are described in detail and the equivalent losses for different mesh positions were calculated.

Blume, H. J. C.

Application of closed-form solutions to a mesh point field in silicon solar cells

A computer simulation method is discussed that provides for equivalent simulation accuracy, but that exhibits significantly lower CPU running time per bias point compared to other techniques. This new method is applied to a mesh point field as is customary in numerical integration (NI) techniques. The assumption of a linear approximation for the dependent variable, which is typically used in the finite difference and finite element NI methods, is not required. Instead, the set of device transport equations is applied to, and the closed-form solutions obtained for, each mesh point. The mesh point field is generated so that the coefficients in the set of transport equations exhibit small changes between adjacent mesh points. Application of this method to high-efficiency silicon solar cells is described; and the method by which Auger recombination, ambipolar considerations, built-in and induced electric fields, bandgap narrowing, carrier confinement, and carrier diffusivities are treated. Bandgap narrowing has been investigated using Fermi-Dirac statistics, and these results show that bandgap narrowing is more pronounced and that it is temperature-dependent in contrast to the results based on Boltzmann statistics.

Lamorte, M. F.