Search NASA⌕ Search

SEARCH · Search NASA

Results for “meshing framework”

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

PLUM: Parallel Load Balancing for Adaptive Unstructured Meshes

Mesh adaption is a powerful tool for efficient unstructured-grid computations but causes load imbalance among processors on a parallel machine. We present a novel method called PLUM to dynamically balance the processor workloads with a global view. This paper presents the implementation and integration of all major components within our dynamic load balancing strategy for adaptive grid calculations. Mesh adaption, repartitioning, processor assignment, and remapping are critical components of the framework that must be accomplished rapidly and efficiently so as not to cause a significant overhead to the numerical simulation. A data redistribution model is also presented that predicts the remapping cost on the SP2. This model is required to determine whether the gain from a balanced workload distribution offsets the cost of data movement. Results presented in this paper demonstrate that PLUM is an effective dynamic load balancing strategy which remains viable on a large number of processors.

Oliker, Leonid↗

Parallel Load Balancing for Adaptive Unstructured Meshes

Mesh adaption is a powerful tool for efficient unstructured-grid computations but causes load imbalance among processors on a parallel machine. We describe a novel method to dynamically balance the processor workloads with a global view. Mesh question, repartitioning, processor assignment, and remapping are critical components of the framework that must be accomplished rapidly and efficiently so as not to cause a significant overhead to the numerical simulation. A data redistribution model will also be presented that predicts the remapping cost. This model is required to determine whether the gain from a balanced workload distribution offsets the cost of data movement. Results presented will demonstrate that this is an effective dynamic load balancing strategy which remains viable on a large number of processors.

Biswas, Rupak↗

A Numerical Investigation of Parachute Deployment in Supersonic Flow

A validated computational fluid-structure interaction method applied to supersonic parachute inflation is extended to consider a unique, parallel self-contact algorithm, porosity interface conditions on the canopy, improved flow field sampling procedures for obtaining quality loading on the parachute canopy, adaptive mesh refinement, and improved treatment of thin geometries in an immersed boundary framework. These extensions are discussed in detail and demonstrated individually on test problems. Finally, the developments are brought together for demonstration on a sub-scale MSL parachute geometry.

Boustani, Jonathan↗

Fan Noise Predictions of the NASA Source Diagnostic Test using Unsteady Simulations with LAVA Part II - Tonal and Broadband Noise Assessment

The NASA Source Diagnostic Test (SDT) campaign experimental data is used for validation of a sliding mesh technique recently implemented within the Launch, Ascent, and Vehicle Aerodynamics computational framework for time-accurate simulation of rotating fans. The far-field acoustics are analyzed in this work, building upon the aerodynamic validation studies previously published in Part I. Two modeling approaches are explored: the unsteady Reynolds-averaged Navier Stokes (uRANS) with the negative Spalart-Allmaras (SA-neg) turbulence model closure, and a hybrid Reynolds-averaged Navier Stokes/large-eddy simulation (RANS/LES) paradigm employing a zonal detached-eddy simulation (ZDES) closure with enhanced shielding protection. Two convective flux scheme approaches with different dissipation properties are also explored with ZDES. The Ffowcs-Williams and Hawkings (FW-H) permeable surface approach is used for propagation of the near-field acoustics to the far-field microphone locations. This work analyzes the low-speed (approach) condition, characterized by a fan rotation speed of 7808 rotations-per-minute (RPM). Three different grid levels ranging between 200 million and 1.1 billion grid points are considered. Results show good prediction of the broadband noise levels at sideline angles ranging between 70° and 110°. The forward arc observers show an under-prediction of the overall sound pressure levels (OASPL) even at the fine grid level. A breakdown of the inlet and exhaust contributions reveals a steep drop-off in the broadband noise levels past a blade-passing frequency (BPF) of 1.5, potentially caused by a lack of resolved small-scale turbulent fluctuations in the interstage region. Past 110° the OASPL are over-predicted by up to 10 dB due to an over-prediction of the low-frequency broadband noise levels in the aft arc. The detuning of BPF2 caused by small deviations in the blade stagger angle around the wheel is captured, and a corresponding decrease in the sound power level for this cut-on tone is observed.

AATT↗

Parallel Tetrahedral Mesh Adaptation with Dynamic Load Balancing

The ability to dynamically adapt an unstructured grid is a powerful tool for efficiently solving computational problems with evolving physical features. In this paper, we report on our experience parallelizing an edge-based adaptation scheme, called 3D_TAG. using message passing. Results show excellent speedup when a realistic helicopter rotor mesh is randomly refined. However. performance deteriorates when the mesh is refined using a solution-based error indicator since mesh adaptation for practical problems occurs in a localized region., creating a severe load imbalance. To address this problem, we have developed PLUM, a global dynamic load balancing framework for adaptive numerical computations. Even though PLUM primarily balances processor workloads for the solution phase, it reduces the load imbalance problem within mesh adaptation by repartitioning the mesh after targeting edges for refinement but before the actual subdivision. This dramatically improves the performance of parallel 3D_TAG since refinement occurs in a more load balanced fashion. We also present optimal and heuristic algorithms that, when applied to the default mapping of a parallel repartitioner, significantly reduce the data redistribution overhead. Finally, portability is examined by comparing performance on three state-of-the-art parallel machines.

Oliker, Leonid↗

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↗

Load Balancing Sequences of Unstructured Adaptive Grids

Mesh adaption is a powerful tool for efficient unstructured grid computations but causes load imbalance on multiprocessor systems. To address this problem, we have developed PLUM, an automatic portable framework for performing adaptive large-scale numerical computations in a message-passing environment. This paper makes several important additions to our previous work. First, a new remapping cost model is presented and empirically validated on an SP2. Next, our load balancing strategy is applied to sequences of dynamically adapted unstructured grids. Results indicate that our framework is effective on many processors for both steady and unsteady problems with several levels of adaption. Additionally, we demonstrate that a coarse starting mesh produces high quality load balancing, at a fraction of the cost required for a fine initial mesh. Finally, we show that the data remapping overhead can be significantly reduced by applying our heuristic processor reassignment algorithm.

Biswas, Rupak↗

A Cartesian grid approach with hierarchical refinement for compressible flows

Many numerical studies of flows that involve complex geometries are limited by the difficulties in generating suitable grids. We present a Cartesian boundary scheme for two-dimensional, compressible flows that is unfettered by the need to generate a computational grid and so it may be used, routinely, even for the most awkward of geometries. In essence, an arbitrary-shaped body is allowed to blank out some region of a background Cartesian mesh and the resultant cut-cells are singled out for special treatment. This is done within a finite-volume framework and so, in principle, any explicit flux-based integration scheme can take advantage of this method for enforcing solid boundary conditions. For best effect, the present Cartesian boundary scheme has been combined with a sophisticated, local mesh refinement scheme, and a number of examples are shown in order to demonstrate the efficacy of the combined algorithm for simulations of shock interaction phenomena.

Quirk, James J.↗

Tangle-Free Mesh Motion for Ablation Simulations

Problems involving mesh motion-which should not be mistakenly associated with moving mesh methods, a class of adaptive mesh redistribution techniques-are of critical importance in numerical simulations of the thermal response of melting and ablative materials. Ablation is the process by which material vaporizes or otherwise erodes due to strong heating. Accurate modeling of such materials is of the utmost importance in design of passive thermal protection systems ("heatshields") for spacecraft, the layer of the vehicle that ensures survival of crew and craft during re-entry. In an explicit mesh motion approach, a complete thermal solve is first performed. Afterwards, the thermal response is used to determine surface recession rates. These values are then used to generate boundary conditions for an a posteriori correction designed to update the location of the mesh nodes. Most often, linear elastic or biharmonic equations are used to model this material response, traditionally in a finite element framework so that complex geometries can be simulated. A simple scheme for moving the boundary nodes involves receding along the surface normals. However, for all but the simplest problem geometries, evolution in time following such a scheme will eventually bring the mesh to intersect and "tangle" with itself, inducing failure. This presentation demonstrates a comprehensive and sophisticated scheme that analyzes the local geometry of each node with help from user-provided clues to eliminate the tangle and enable simulations on a wide-class of difficult problem geometries. The method developed is demonstrated for linear elastic equations but is general enough that it may be adapted to other modeling equations. The presentation will explicate the inner workings of the tangle-free mesh motion algorithm for both two and three-dimensional meshes. It will show abstract examples of the method's success, including a verification problem that demonstrates its accuracy and correctness. The focus of the presentation will be on the algorithm; specifics on how the techniques may be used in spacecraft design will be not discussed.

Droba, Justin↗

Dust Particle Aeroheating Calculations for Mars Entry Hypersonic Flows

The purpose of the current research effort is to predict particle trajectories in a hypersonic flowfield while accounting for all relevant particle-fluid and particle-particle interactions. The resulting particle solver is titled DUst Simulation & Tracking (DUST) and works in conjunction with the US3D CFD flow solver. The key elements underpinning the current work are an efficient mesh-localiztion algorithm, time-driven hard-sphere model, point-to-point MPI framework, coarse-graining using computational parcels, and high-order Adams-Bashforth time stepping. The coupled US3D-DUST framework will be applied to conduct a multi-physics examination of dust-laden flows around the the Schiarparelli capsule.

Dusty Flows↗

Dust Particle Aeroheating Calculations for Mars Entry Hypersonic Flows

The purpose of the current research effort is to predict particle trajectories in a hypersonic flowfield while accounting for all relevant particle-fluid and particle-particle interactions. The resulting particle solver is titled DUst Simulation & Tracking (DUST) and works in conjunction with the US3D CFD flow solver. The key elements underpinning the current work are an efficient mesh-localiztion algorithm, time-driven hard-sphere model, point-to-point MPI framework, coarse-graining using computational parcels, and high-order Adams-Bashforth time stepping. The coupled US3D-DUST framework will be applied to conduct a multi-physics examination of dust-laden flows around the the Mars 2020 capsule.

Dusty Flows↗

A Modular Conjugate Heat Transfer Optimization Framework for Thermal Management of Electric Aircraft

Conjugate heat transfer (CHT) analysis and optimization is a powerful method for improving thermal management, as it simultaneously resolves the temperature distribution in both fluid and solid domains. This paper presents a modular, discrete adjoint-based CHT optimization capability integrated within the OpenMDAO/MPhys framework. A unique feature of the proposed framework is its flexibility to extend to multidisciplinary optimization, including aero-structural-thermal applications. The fluid domain is modeled using a finite-volume Computational Fluid Dynamics (CFD) solver, and the solid domain with a conduction heat transfer solver. A mixed Neumann-Dirichlet boundary condition is developed to enable full submersion of the solid geometry within the fluid domain, while ensuring consistent temperature and heat flux coupling at the CHT interface. Gradient-based optimization is performed; the gradients are efficiently computed using the discrete adjoint solvers implemented in DAFoam. To demonstrate the method, this paper considers two cases related to electric aircraft thermal management: a U-bend heat exchanger and an actively cooled battery pack. The U-bend case aims to minimize pressure loss while maximizing heat flux by changing the pipe geometry. The optimized design reduces pressure loss by 52.7% and increases total heat flux by 2.3%. In the battery pack case, a 3-by-3 cell configuration is cooled by ambient airflow, with constant heat generation prescribed in the cells. The battery casing shape serves as the design variable, and the objective function is a weighted sum of pressure loss and pack weight, subject to a maximum temperature constraint. The optimized design achieves a 44.6% reduction in pressure loss and a 1.5% reduction in weight, while satisfying the thermal constraint. To ensure the reliability of the optimized designs, this study validates coarse-mesh, steady-state predictions against fine-mesh unsteady simulations, demonstrating consistency within acceptable errors. This work demonstrates the potential of the developed framework to enable rapid, high-fidelity design of thermal management systems for electric aircraft.

heat transfer↗

Modeling Progressive Failure of Bonded Joints Using a Single Joint Finite Element

Enhanced finite elements are elements with an embedded analytical solution which can capture detailed local fields, enabling more efficient, mesh-independent finite element analysis. In the present study, an enhanced finite element is applied to generate a general framework capable of modeling an array of joint types. The joint field equations are derived using the principle of minimum potential energy, and the resulting solutions for the displacement fields are used to generate shape functions and a stiffness matrix for a single joint finite element. This single finite element thus captures the detailed stress and strain fields within the bonded joint, but it can function within a broader structural finite element model. The costs associated with a fine mesh of the joint can thus be avoided while still obtaining a detailed solution for the joint. Additionally, the capability to model non-linear adhesive constitutive behavior has been included within the method, and progressive failure of the adhesive can be modeled by using a strain-based failure criteria and re-sizing the joint as the adhesive fails. Results of the model compare favorably with experimental and finite element results.

Stapleton, Scott E.↗

Predicting Quadcopter Drone Noise Using the Lattice Boltzmann Method

The market for new vertical takeoff and landing vehicles, including autonomous urban air taxis and drones for applications such as package delivery, imaging, and surveillance, is growing rapidly. However, aerodynamic noise continues to be the biggest roadblock to community acceptance and adoption. To predict the aerodynamic noise generated by an isolated quadcopter drone, derived from from first principles, we used the Lattice Boltzmann flow solver within NASA’s Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework. The solver’s computational efficiency, and the complete absence of labor-intensive manual volume mesh generation in the workflow, are key to making routine aeroacoustic analysis of urban air taxis and drones from first principles possible.

Cadieux, Francois↗

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↗

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↗

Image Segmentation, Registration, Compression, and Matching

A novel computational framework was developed of a 2D affine invariant matching exploiting a parameter space. Named as affine invariant parameter space (AIPS), the technique can be applied to many image-processing and computer-vision problems, including image registration, template matching, and object tracking from image sequence. The AIPS is formed by the parameters in an affine combination of a set of feature points in the image plane. In cases where the entire image can be assumed to have undergone a single affine transformation, the new AIPS match metric and matching framework becomes very effective (compared with the state-of-the-art methods at the time of this reporting). No knowledge about scaling or any other transformation parameters need to be known a priori to apply the AIPS framework. An automated suite of software tools has been created to provide accurate image segmentation (for data cleaning) and high-quality 2D image and 3D surface registration (for fusing multi-resolution terrain, image, and map data). These tools are capable of supporting existing GIS toolkits already in the marketplace, and will also be usable in a stand-alone fashion. The toolkit applies novel algorithmic approaches for image segmentation, feature extraction, and registration of 2D imagery and 3D surface data, which supports first-pass, batched, fully automatic feature extraction (for segmentation), and registration. A hierarchical and adaptive approach is taken for achieving automatic feature extraction, segmentation, and registration. Surface registration is the process of aligning two (or more) data sets to a common coordinate system, during which the transformation between their different coordinate systems is determined. Also developed here are a novel, volumetric surface modeling and compression technique that provide both quality-guaranteed mesh surface approximations and compaction of the model sizes by efficiently coding the geometry and connectivity/topology components of the generated models. The highly efficient triangular mesh compression compacts the connectivity information at the rate of 1.5-4 bits per vertex (on average for triangle meshes), while reducing the 3D geometry by 40-50 percent. Finally, taking into consideration the characteristics of 3D terrain data, and using the innovative, regularized binary decomposition mesh modeling, a multistage, pattern-drive modeling, and compression technique has been developed to provide an effective framework for compressing digital elevation model (DEM) surfaces, high-resolution aerial imagery, and other types of NASA data.

Yadegar, Jacob↗

Impact of Load Balancing on Unstructured Adaptive Grid Computations for Distributed-Memory Multiprocessors

The computational requirements for an adaptive solution of unsteady problems change as the simulation progresses. This causes workload imbalance among processors on a parallel machine which, in turn, requires significant data movement at runtime. We present a new dynamic load-balancing framework, called JOVE, that balances the workload across all processors with a global view. Whenever the computational mesh is adapted, JOVE is activated to eliminate the load imbalance. JOVE has been implemented on an IBM SP2 distributed-memory machine in MPI for portability. Experimental results for two model meshes demonstrate that mesh adaption with load balancing gives more than a sixfold improvement over one without load balancing. We also show that JOVE gives a 24-fold speedup on 64 processors compared to sequential execution.

Biswas, Rupak↗