SEARCH · Engineering Papers

Results for “mesh optimization”

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 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Adaptive Mesh Refinement Simulations for Turbulent Reacting Flow

With the increased availability of exascale computing hardware, detailed simulations of realistic devices can be performed at practically relevant time and length scales. Insights into the multiscale driving mechanisms in compressible reacting flow systems with complex geometry, such as combustors, can be used for design optimization and technology improvements. However, to effectively perform these simulations, advanced numerical algorithms must be used to maintain solution accuracy without incurring undue computational costs. PeleC, part of the Pele suite of codes, leverages block-structured adaptive mesh refinement (AMR) through the AMReX library to capture fine-scale flow features in compressible reacting flows. In this talk, we discuss recent improvements to the numerical algorithms, particularly in regard to describing flows at complex boundary structures, and PeleC's performance on exascale computing hardware. We will demonstrate that PeleC is well-suited for modern, extreme-scale, heterogenous compute platforms.

combustion↗

NASA NTRS1998

Performance Analysis and Portability of the PLUM Load Balancing System

The ability to dynamically adapt an unstructured mesh is a powerful tool for solving computational problems with evolving physical features; however, an efficient parallel implementation is rather difficult. To address this problem, we have developed PLUM, an automatic portable framework for performing adaptive numerical computations in a message-passing environment. PLUM requires that all data be globally redistributed after each mesh adaption to achieve load balance. We present an algorithm for minimizing this remapping overhead by guaranteeing an optimal processor reassignment. We also show that the data redistribution cost can be significantly reduced by applying our heuristic processor reassignment algorithm to the default mapping of the parallel partitioner. Portability is examined by comparing performance on a SP2, an Origin2000, and a T3E. Results show that PLUM can be successfully ported to different platforms without any code modifications.

Oliker, Leonid↗

NASA NTRS2004

Improved Thermoplastic/Iron-Particle Transformer Cores

A method of fabricating improved transformer cores from composites of thermoplastic matrices and iron-particles has been invented. Relative to commercially available laminated-iron-alloy transformer cores, the cores fabricated by this method weigh less and are less expensive. Relative to prior polymer-matrix/ iron-particle composite-material transformer cores, the cores fabricated by this method can be made mechanically stronger and more magnetically permeable. In addition, whereas some prior cores have exhibited significant eddy-current losses, the cores fabricated by this method exhibit very small eddy-current losses. The cores made by this method can be expected to be attractive for use in diverse applications, including high-signal-to-noise transformers, stepping motors, and high-frequency ignition coils. The present method is a product of an experimental study of the relationships among fabrication conditions, final densities of iron particles, and mechanical and electromagnetic properties of fabricated cores. Among the fabrication conditions investigated were molding pressures (83, 104, and 131 MPa), and molding temperatures (250, 300, and 350 C). Each block of core material was made by uniaxial-compression molding, at the applicable pressure/temperature combination, of a mixture of 2 weight percent of LaRC (or equivalent high-temperature soluble thermoplastic adhesive) with 98 weight percent of approximately spherical iron particles having diameters in the micron range. Each molded block was cut into square cross-section rods that were used as core specimens in mechanical and electromagnetic tests. Some of the core specimens were annealed at 900 C and cooled slowly before testing. For comparison, a low-carbon-steel core was also tested. The results of the tests showed that density, hardness, and rupture strength generally increased with molding pressure and temperature, though the correlation was rather weak. The weakness of the correlation was attributed to the pores in the specimens. The maximum relative permeabilities of cores made without annealing ranged from 30 to 110, while those of cores made with annealing ranged from 900 to 1,400. However, the greater permeabilities of the annealed specimens were not associated with noticeably greater densities. The major practical result of the investigation was the discovery of an optimum distribution of iron-particle sizes: It was found that eddy-current losses in the molded cores were minimized by using 100 mesh (corresponding to particles with diameters less than or equal to 100 m) iron particles. The effect of optimization of particle sizes on eddy-current losses is depicted in the figure.

Wincheski, Russell A.↗

NASA NTRS2014

Active Vibration Control for Helicopter Interior Noise Reduction Using Power Minimization

This report describes work performed by United Technologies Research Center (UTRC) for NASA Langley Research Center (LaRC) under Contract NNL11AA06C. The objective of this program is to develop technology to reduce helicopter interior noise resulting from multiple gear meshing frequencies. A novel active vibration control approach called Minimum Actuation Power (MAP) is developed. MAP is an optimal control strategy that minimizes the total input power into a structure by monitoring and varying the input power of controlling sources. MAP control was implemented without explicit knowledge of the phasing and magnitude of the excitation sources by driving the real part of the input power from the controlling sources to zero. It is shown that this occurs when the total mechanical input power from the excitation and controlling sources is a minimum. MAP theory is developed for multiple excitation sources with arbitrary relative phasing for single or multiple discrete frequencies and controlled by a single or multiple controlling sources. Simulations and experimental results demonstrate the feasibility of MAP for structural vibration reduction of a realistic rotorcraft interior structure. MAP control resulted in significant average global vibration reduction of a single frequency and multiple frequency excitations with one controlling actuator. Simulations also demonstrate the potential effectiveness of the observed vibration reductions on interior radiated noise.

Mendoza, J.↗

DOE OSTI2017

CAD-Integrated Direct-to-HPC Simulation for Design and Optimization

Explore the source record for details and available documents.

97 MATHEMATICS AND COMPUTING↗

NASA NTRSDate not supplied

A Design Framework for Thick and Thin Tow-Steered Composites Using Mechanics of Structure Genome

The design of tow-steered composites is one of the most popular and promising topics under the class of variable stiffness structures. In this work, the authors propose a new design framework and tool for the optimization of tow-steered composites. Mechanics of structure genome (MSG) provides accurate computation of plate/shell section properties including thick and highly curved structures. The computed section properties can be directly used in commercial computer-aided engineering tools. Tow-steered composites are modeled through a design framework wherein tow paths are parameterized in a general way and projected onto a finite element mesh, and MSG calculates the local shell/plate properties. The result is a complete workflow from design parameter input to structural performance evaluation. An open-source optimization software package is used to enable constrained optimization of structural stiffness for the broad tow-steered composite design space. Numerical examples are provided to demonstrate the capabilities of this tool and the promise and potential of tow-steered designs that optimize structural performance while satisfying manufacturing constraints.

Su Tian↗

NASA NTRSDate not supplied

A Design Framework for Thick and Thin Tow-Steered Composites Using Mechanics of Structure Genome

Su Tian↗

NASA NTRS2015

Acoustically Tailored Composite Rotorcraft Fuselage Panels

A rotorcraft roof sandwich panel has been redesigned to optimize sound power transmission loss (TL) and minimize structure-borne sound for frequencies between 1 and 4 kHz where gear meshing noise from the transmission has the most impact on speech intelligibility. The roof section, framed by a grid of ribs, was originally constructed of a single honeycomb core/composite face sheet panel. The original panel has coincidence frequencies near 700 Hz, leading to poor TL across the frequency range of 1 to 4 kHz. To quiet the panel, the cross section was split into two thinner sandwich subpanels separated by an air gap. The air gap was sized to target the fundamental mass-spring-mass resonance of the double panel system to less than 500 Hz. The panels were designed to withstand structural loading from normal rotorcraft operation, as well as 'man-on-the-roof' static loads experienced during maintenance operations. Thin layers of VHB 9469 viscoelastomer from 3M were also included in the face sheet ply layups, increasing panel damping loss factors from about 0.01 to 0.05. Measurements in the NASA SALT facility show the optimized panel provides 6-11 dB of acoustic transmission loss improvement, and 6-15 dB of structure-borne sound reduction at critical rotorcraft transmission tonal frequencies. Analytic panel TL theory simulates the measured performance quite well. Detailed finite element/boundary element modeling of the baseline panel simulates TL slightly more accurately, and also simulates structure-borne sound well.

Hambric, Stephen↗

NASA NTRS2016

Soil Moisture Active Passive Mission Mass Properties Optimization and Analysis

The Soil Moisture Active Passive (SMAP) mission creates global soil moisture measurements every 2-3 days by spinning a lightweight, ~6 meter diameter deployable mesh reflector. This represents the largest spinning antenna ever to be used on a NASA satellite. The paper discusses the mass properties challenges that were overcome in implementing the mission, including the definition of a combined imbalance term called the Effective Product of Inertia, the selection of an optimum Spun Instrument focal length and edge offset that fostered a reduced mass properties sensitivity, the inclusion of late game Spun Instrument configuration flexibility, the use of balance mass locations that leveraged the large moment arm provided by the reflector, and a subassembly measurement campaign that reduced mass properties uncertainties on testable components.

Hendry, Morgan L., IV↗

DOE OSTI2026

A meshing framework for digital twins for extrusion based additive manufacturing

Additive manufacturing (AM) allows for manufacturing of complex three-dimensional geometries not typically realizable with standard manufacturing practices. The internal microstructure of AM components has a significant impact on mechanical, vibrational, and shock properties and permits richer design space when this is controllable. Due to complex interactions of internal geometry of an extrusion-based AM component, it is common practice to assume homogeneous behavior or to perform characterization testing on specific toolpath configurations. To avoid testing or material waste, it is necessary to develop a consistently accurate numerical simulation framework with relevant boundary value problems that can handle the complicated geometry of internal material microstructure present in AM components. Herein, a framework is proposed to directly create computational meshes suitable for finite element analysis (FEA) of the fine-scale features generated from extrusion-based AM tool paths to maintain a strong process–structure–property-performance linkage. This mesh can be manually or automatically analyzed using standard FEA simulations such as quasi-static preloading or modal analysis. The framework allows an in-silico assessment of a target AM geometry where fine-scale features greatly impact quantities of design interest such as in soft elastomeric lattices where toolpath infill can greatly influence the self-contact of a structure in compression, which we use as a motivating exemplar. This approach greatly reduces both time and resource waste present in traditional build and test design cycles for non-intuitive design spaces, and acts as a tool for use in the production of a key component of a digital twin, a mesh suitable for finite element analysis. In conclusion, it also further allows for the exploration of toolpath infill to optimize component properties beyond simple linear properties such as density and stiffness.

Additive manufacturing↗

NASA NTRS1987

Concurrent and vectorized mixed time, explicit nonlinear structural dynamics algorithms

A nonlinear structural dynamics program with an element library that exploits parallel processing is described. The aim is to exploit scheduling-allocation so that parallel processing and vectorization can effectively be treated in a general purpose program with explicit time integration and different time steps in different parts of the mesh. The program uses an element group scheme, which, as a by-product, also provides an automatic scheme for assigning different time steps to different parts of the mesh. The program has been tested on the Alliant FX/8; it shows a fivefold improvement in speed over compiler optimization.

Belytschko, Ted↗

NASA NTRS1997

RIACS

Topics considered include: high-performance computing; cognitive and perceptual prostheses (computational aids designed to leverage human abilities); autonomous systems. Also included: development of a 3D unstructured grid code based on a finite volume formulation and applied to the Navier-stokes equations; Cartesian grid methods for complex geometry; multigrid methods for solving elliptic problems on unstructured grids; algebraic non-overlapping domain decomposition methods for compressible fluid flow problems on unstructured meshes; numerical methods for the compressible navier-stokes equations with application to aerodynamic flows; research in aerodynamic shape optimization; S-HARP: a parallel dynamic spectral partitioner; numerical schemes for the Hamilton-Jacobi and level set equations on triangulated domains; application of high-order shock capturing schemes to direct simulation of turbulence; multicast technology; network testbeds; supercomputer consolidation project.

Oliger, Joseph↗

NASA NTRS2009

Discrete Adjoint-Based Design Optimization of Unsteady Turbulent Flows on Dynamic Unstructured Grids

An adjoint-based methodology for design optimization of unsteady turbulent flows on dynamic unstructured grids is described. The implementation relies on an existing unsteady three-dimensional unstructured grid solver capable of dynamic mesh simulations and discrete adjoint capabilities previously developed for steady flows. The discrete equations for the primal and adjoint systems are presented for the backward-difference family of time-integration schemes on both static and dynamic grids. The consistency of sensitivity derivatives is established via comparisons with complex-variable computations. The current work is believed to be the first verified implementation of an adjoint-based optimization methodology for the true time-dependent formulation of the Navier-Stokes equations in a practical computational code. Large-scale shape optimizations are demonstrated for turbulent flows over a tiltrotor geometry and a simulated aeroelastic motion of a fighter jet.

Nielsen, Eric J.↗

DOE OSTI2017

Optimisation of dynamic reactive power sources using mesh adaptive direct search

Dynamic reactive power sources can be used to effectively mitigate the fault-induced delayed voltage recovery and transient voltage instability issues. When many var sources need to be installed at planned locations, optimisation of their sizes is a complicated non-linear optimisation problem due to its non-convexity and the dependence of the constraint on time-series trajectories of post-fault voltages. Solving this optimal sizing problem needs to utilise both a non-linear optimisation solver and a power system differential-algebraic equation solver. This study proposes a new approach for solving this problem under both a single contingency and multiple contingencies by using the mesh adaptive direct search algorithm interfaced with a power system simulator. The proposed approach is validated by case studies on a North American Eastern Interconnection model to optimise the sizes of planned STATCOMs against critical contingencies.

Huang, Weihong↗

NASA NTRS1976

The multigrid method: Fast relaxation

A multi-level grid method was studied as a possible means of accelerating convergence in relaxation calculations for transonic flows. The method employs a hierarchy of grids, ranging from very coarse (e.g. 4 x 2 mesh cells) to fine (e.g. 64 x 32); the coarser grids are used to diminish the magnitude of the smooth part of the residuals, hopefully with far less total work than would be required with optimal iterations on the finest grid. To date the method was applied quite successfully to the solution of the transonic small-disturbance equation for the velocity potential in conservation form. Nonlifting transonic flow past a parabolic arc airfoil is the example studied, with meshes of both constant and variable step size.

South, J. C., Jr.↗

NASA NTRS1977

Application of a multi-level grid method to transonic flow calculations

A multi-level grid method has been studied as a possible means of accelerating convergence in relaxation calculations for transonic flows. The method employs a hierarchy of grids, ranging from very coarse (e.g., 8 x 2 mesh cells) to fine (e.g., 128 x 32); the coarser grids are used to diminish the magnitude of the smooth part of the residuals, hopefully with far less total work than would be required with, say, optimal SLOR iterations on the finest grid. The method was applied to the solution of the transonic small-disturbance equation for the velocity potential in the conservation form. Nonlifting transonic flow past a parabolic-arc airfoil is the example studied, with meshes of both constant and variable step size.

South, J. C., Jr.↗

DOE OSTI2025

A Block-Structured Adaptive Mesh Framework to Solve Radiation Transfer Equation in Irregular Embedded Geometries

Radiation transport arises in various scientific, industrial, and medical fields, and understanding its effect in applications is needed to make accurate predictions, safety assessments and performance optimizations. Solving the Radiation Transport Equation (RTE) is challenging due to its integro-differential nature, which involves both differential and integral terms. The differential term describes the change in radiation intensity due to absorption and emission, while the integral term accounts for scattering. The accurate modeling of radiation is further complicated in many applications due to the complex, irregular geometries. Various methods exist for solving the RTE, including the zonal, Monte Carlo, spherical harmonics, discrete ordinates, and finite volume methods. Traditional mesh-based approaches, which rely on structured or unstructured meshes, struggle with irregular geometries due to: a) the difficulty of conforming structured grids to irregular domains, b) challenges in enforcing boundary conditions correctly, and c) the additional computational cost of unstructured mesh methods. This work presents a second-order accurate method for solving the RTE in irregular geometries. The radiation intensity is discretized using the finite-volume method in both spatial and angular directions on regular Cartesian grid blocks. Leveraging the block-structured adaptive mesh refinement (AMR) framework provided by AMReX, our method refines the grid locally to reduce spatial discretization error, ensuring a converged numerical solution while minimizing computational costs elsewhere. A two-stage deferred correction approach is employed: First, a first-order discretization on grid blocks is solved using an algebraic multigrid method in HYPRE. Second, a correction term is applied explicitly to achieve second-order accuracy. The correction term is calculated by approximating the radiation flux on cell faces using a Total Variation Diminishing (TVD) scheme. This approach ensures quick convergence of the multigrid method while preserving higher-order accuracy of the numerical solution. Irregular geometries are resolved as embedded boundaries (EB), resulting in both cut cells and regular cells. In cut cells, we modify the fluxes using face fractions and incorporate additional contributions from EB boundary conditions. To ensure higher-order convergence near the EB interface, the correction term is modified by interpolating the radiation intensity to fictitious ghost points. The implementation takes advantage of modern supercomputers by leveraging AMReX’sMPI/X parallelization strategy where X can be MPI or a GPU accelerator including CUDA, HIP and DPC++. We validate our solver using classical test cases, both with and without EB, demonstrating accuracy and efficiency. Additionally, we analyze the impact of adaptive mesh refinement on solution accuracy and computational cost, highlighting the advantages of our approach for high-resolution radiation transport simulations.

computational fluid dynamics (CFD)↗

DOE OSTI2018

Three dimensional optimization for subsonic axial turbines operating at high unsteady inlet Mach number

In this paper, a three dimensional optimization strategy is presented to improve the performance of a subsonic axial turbine exposed to high Mach number unsteady inflows encountered in pressure gain combustors. The optimization aims at altering the stator hub and shroud end wells. A multi-objective differential evolution optimizer was used to enhance specific work and efficiency. The optimization routine consisted of a multi-step approach where in a first step an optimization was carried out where the turbine was subjected to a steady high Mach inflow. The turbine stage was parametrized via an in-house turbine design tool based on Bezier curves. Seven design parameters were selected to contour the stator end walls and rotor hub end wall. The parametrized turbine stage was meshed and the steady Reynolds-Averaged Navier-Stokes equations were solved via the commercial solver CFD++ and calculation time of each stage was around 10 hours. Results of the steady optimization identified three major design parameters of stator end walls, which were taken as inputs for the next phase unsteady optimization. For the unsteady simulations, a large inlet fluctuation was imposed and the turbine was modeled with the mixing plane approach. Total calculation time of each individual was 48 hours. An efficiency improvement of 10% from the steady optimization and 13% from the unsteady optimization was found compared to the baseline case. Finally, full unsteady simulations with two stator and three rotor passages were performed on selected steady and unsteady optimized configurations.

03 NATURAL GAS↗