Multi-Scale Meshes for Finite Element and Finite Volume Methods: Ative Device and Guided-Wave Modeling
Explore the source record for details and available documents.
electromagnetic field↗
SEARCH · Engineering Papers
Results for “Mesh Adaptation”
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.
Using the GeoFEST Faulted Region Simulation System
GeoFEST (the Geophysical Finite Element Simulation Tool) simulates stress evolution, fault slip and plastic/elastic processes in realistic materials, and so is suitable for earthquake cycle studies in regions such as Southern California. Many new capabilities and means of access for GeoFEST are now supported. New abilities include MPI-based cluster parallel computing using automatic PYRAMID/Parmetis-based mesh partitioning, automatic mesh generation for layered media with rectangular faults, and results visualization that is integrated with remote sensing data. The parallel GeoFEST application has been successfully run on over a half-dozen computers, including Intel Xeon clusters, Itanium II and Altix machines, and the Apple G5 cluster. It is not separately optimized for different machines, but relies on good domain partitioning for load-balance and low communication, and careful writing of the parallel diagonally preconditioned conjugate gradient solver to keep communication overhead low. Demonstrated thousand-step solutions for over a million finite elements on 64 processors require under three hours, and scaling tests show high efficiency when using more than (order of) 4000 elements per processor. The source code and documentation for GeoFEST is available at no cost from Open Channel Foundation. In addition GeoFEST may be used through a browser-based portal environment available to approved users. That environment includes semi-automated geometry creation and mesh generation tools, GeoFEST, and RIVA-based visualization tools that include the ability to generate a flyover animation showing deformations and topography. Work is in progress to support simulation of a region with several faults using 16 million elements, using a strain energy metric to adapt the mesh to faithfully represent the solution in a region of widely varying strain.
Geophyical Finite Element Simulation Tool (GeoFEST↗
Adaptive computations of flow around a delta wing with vortex breakdown
An adaptive unstructured mesh solution method for the three-dimensional Euler equations was used to simulate the flow around a sharp edged delta wing. Emphasis was on the breakdown of the leading edge vortex at high angle of attack. Large values of entropy, which indicate vortical regions of the flow, specified the region in which adaptation was performed. The aerodynamic normal force coefficients show excellent agreement with wind tunnel data measured by Jarrah, and demonstrate the importance of adaptation in obtaining an accurate solution. The pitching moment coefficient and the location of vortex breakdown are compared with experimental data measured by Hummel and Srinivasan, showing good agreement in cases in which vortex breakdown is located over the wing.
Modiano, David L.↗
Physics-Based Modeling and Simulation of Self-Reacting Friction Stir Welding Using Computational Fluid Dynamics
A physics-based model was developed to simulate the behavior of material in a self-reacting friction stir welding (SR-FSW) process for the joining of metals. This steady-state model builds upon fundamental computational fluid dynamic (CFD) principles within Ansys Fluent to solve the discretized equations. The effective viscosity is calculated using a viscoplastic model using a Sheppard-Wright formulation of flow stress. Numerous advancements have been made in the incorporated physics including (1) temperature-dependent material properties; (2) locally adaptable flow and thermal boundary conditions; and (3) adapting material properties in nugget in response to microstructural changes. Simulation strategies to accelerate computation and improve numerical stability include adapting the mesh refinement and solver relaxation factors during simulation. The result is a highly robust and computationally efficient model capable of providing the material flow and temperature history across the domain. As material history determines the local microstructure and ultimately weld strength, an accurate and detailed physics-based model has the potential to accelerate SR-FSW process development. The model is highly adaptable to changes in process parameters, tool design, or alloy.
Process Modeling↗
Implicit solvers for unstructured meshes
Implicit methods for unstructured mesh computations are developed and tested. The approximate system which arises from the Newton-linearization of the nonlinear evolution operator is solved by using the preconditioned generalized minimum residual technique. These different preconditioners are investigated: the incomplete LU factorization (ILU), block diagonal factorization, and the symmetric successive over-relaxation (SSOR). The preconditioners have been optimized to have good vectorization properties. The various methods are compared over a wide range of problems. Ordering of the unknowns, which affects the convergence of these sparse matrix iterative methods, is also investigated. Results are presented for inviscid and turbulent viscous calculations on single and multielement airfoil configurations using globally and adaptively generated meshes.
Venkatakrishnan, V.↗
Implicit solvers for unstructured meshes
Implicit methods were developed and tested for unstructured mesh computations. The approximate system which arises from the Newton linearization of the nonlinear evolution operator is solved by using the preconditioned GMRES (Generalized Minimum Residual) technique. Three different preconditioners were studied, namely, the incomplete LU factorization (ILU), block diagonal factorization, and the symmetric successive over relaxation (SSOR). The preconditioners were optimized to have good vectorization properties. SSOR and ILU were also studied as iterative schemes. The various methods are compared over a wide range of problems. Ordering of the unknowns, which affects the convergence of these sparse matrix iterative methods, is also studied. Results are presented for inviscid and turbulent viscous calculations on single and multielement airfoil configurations using globally and adaptively generated meshes.
Venkatakrishnan, V.↗
Conference on Advanced Technology for Future Space Systems, Hampton, Va., May 8-10, 1979, Technical Papers
Propulsion systems for spacecraft, satellite communications technology, the design of large light-weight erectable structures for assembly in space, electronics and information processing for spacecraft, and self-diagnostic, fault-tolerant controls based on high memory and processing capabilities are discussed. Topics of the papers include the design of large delta wings for earth-to-orbit transports, dual-fuel propulsion units, magnetoplasmadynamic thrusters, heating rates on blunt-nosed bodies at various angles of attack, remote manipulators for space assembly tasks, solar electric propulsion for planetary missions, deployable space platforms with multiple payloads, the design of large offset-fed antennas, a nonlinear stress-strain relationship for metallic meshes, and adaptive sensors for spacecraft.
Source record↗
Dynamic Load Balancing for Adaptive Computations on Distributed-Memory Machines
Dynamic load balancing is central to adaptive mesh-based computations on large-scale parallel computers. The principal investigator has investigated various issues on the dynamic load balancing problem under NASA JOVE and JAG rants. The major accomplishments of the project are two graph partitioning algorithms and a load balancing framework. The S-HARP dynamic graph partitioner is known to be the fastest among the known dynamic graph partitioners to date. It can partition a graph of over 100,000 vertices in 0.25 seconds on a 64- processor Cray T3E distributed-memory multiprocessor while maintaining the scalability of over 16-fold speedup. Other known and widely used dynamic graph partitioners take over a second or two while giving low scalability of a few fold speedup on 64 processors. These results have been published in journals and peer-reviewed flagship conferences.
Source record↗
Design of Unstructured Adaptive (UA) NAS Parallel Benchmark Featuring Irregular, Dynamic Memory Accesses
We describe the design of a new method for the measurement of the performance of modern computer systems when solving scientific problems featuring irregular, dynamic memory accesses. The method involves the solution of a stylized heat transfer problem on an unstructured, adaptive grid. A Spectral Element Method (SEM) with an adaptive, nonconforming mesh is selected to discretize the transport equation. The relatively high order of the SEM lowers the fraction of wall clock time spent on inter-processor communication, which eases the load balancing task and allows us to concentrate on the memory accesses. The benchmark is designed to be three-dimensional. Parallelization and load balance issues of a reference implementation will be described in detail in future reports.
Feng, Hui-Yu↗
Three Dimensional Multi-Fluid Simulation of Comet Halley
The Cometary Atmosphere Simulator (CASIM) is designed to simulate in 2-D and 3-D the complex interaction between the cometary atmosphere and the hypersonic solar wind using a multi-fluid approach. Our simulator is based on the solution of multi-fluid equations using an efficient adaptive Cartesian mesh. It is designed to use the capabilities of highly parallel super-cluster computers.
Benna, M.↗
The common component architecture (CCA) applied to sequential and parallel computational electromagnetic applications
Explore the source record for details and available documents.
parallel computing software engineering componenet↗
Kestrel Results at Liftoff Conditions for a Space Launch System Configuration in Proximity to the Launch Tower
Aerodynamic data books for Space Launch System vehicles require databases for the integrated forces and moments and section loads during liftoff and transition to the ascent phase of flight. While the force and moment database can be generated from wind tunnel results, computational analyses are necessary to provide the extensive surface information required to generate proper lineloads. Of the two flight regimes, the liftoff problem is the more costly and complex situation to simulate, as it requires modeling of the vehicle in proximity to the launch tower. The effects of massive separation on the leeward pressure fields of both the tower and vehicle are not well captured with RANS methods, necessitating the use of more advanced methods, such as Delayed Detached Eddy Simulation, in conjunction with computational grids sufficiently refined to resolve the wakes. Details on the computational setup for employing the Kestrel flow solver to address the liftoff problem are presented. The methodology involves the use of independent unstructured near-body grids for the vehicle and the tower, overset by a solution adaptive Cartesian off-body grid. Results from the simulations are compared to experimental results from a test in the NASA Langley Research Center 14- by 22-Foot Subsonic Tunnel.
Computational fluid dynamics↗
Pterodactyl: Aerodynamic and Aerothermal Modeling for a Symmetric Deployable Earth Entry Vehicle with Flaps
NASA’s Pterodactyl project has investigated a deployable atmospheric entry vehicle integrated with a flap control system that provides precision targeting during reentry. The control system consists of eight flaps mounted at the edge of a heatshield that can deflect in and out of the flow. An aerodynamic and aerothermodynamic analysis process was developed for preliminary vehicle design and aerodynamic database generation using engineering and CFD tools with varying levels of fidelity. The objective of this analysis was to: 1) understand the inherent aerodynamics, 2) provide an aerodynamics database for stability and control analysis, and 3) provide qualitative and quantitative aeroheating analysis for Thermal Protection System modeling of the flaps. A high fidelity Euler code, Cart3D, was used to resolve complex flow features such as secondary shocks and shock impingement. A wide range of supersonic (Mach 2) and hypersonic Mach numbers (up to Mach 40) were tested with an Earth atmosphere model. A process was developed to utilize the adaptive volume mesh generator utility of the perfect gas model of Cart3D to create meshes with high cell efficiency and numerical stability prior to running a 2nd order accurate solution with a real gas model. High-fidelity aerothermal CFD simulations were performed using US3D to further improve the aerothermal analysis on hypersonic flow around a complex entry vehicle shape including viscosity, chemical reactions of air species, vibrational energy, and catalytic surface reactions. In this process, it was found that the flap control system provided multi-axis control that can be utilized for entry precision targeting. Additionally, the Pterodactyl vehicle can achieve up to a trim L/D of 0.2. Finally, the increased fidelity of the aerothermal heating environments revealed that the shear stress contributes to increasing the heating on the flaps.
Deployable Entry Vehicle↗
Pterodactyl: Aerodynamic and Aerothermal Modeling for a Symmetric Deployable Earth Entry Vehicle with Flaps
NASA’s Pterodactyl project has investigated a deployable atmospheric entry vehicle integrated with a flap control system that provides precision targeting during reentry. The control system consists of eight flaps mounted at the edge of a heatshield that can deflect in and out of the flow. An aerodynamic and aerothermodynamic analysis process was developed for preliminary vehicle design and aerodynamic database generation using engineering and CFD tools with varying levels of fidelity. The objective of this analysis was to: 1) understand the inherent aerodynamics, 2) provide an aerodynamics database for stability and control analysis, and 3) provide qualitative and quantitative aeroheating analysis for Thermal Protection System modeling of the flaps. A high fidelity Euler code, Cart3D, was used to resolve complex flow features such as secondary shocks and shock impingement. A wide range of supersonic (Mach 2) and hypersonic Mach numbers (up to Mach 40) were tested with an Earth atmosphere model. A process was developed to utilize the adaptive volume mesh generator utility of the perfect gas model of Cart3D to create meshes with high cell efficiency and numerical stability prior to running a 2nd order accurate solution with a real gas model. High-fidelity aerothermal CFD simulations were performed using US3D to further improve the aerothermal analysis on hypersonic flow around a complex entry vehicle shape including viscosity, chemical reactions of air species, vibrational energy, and catalytic surface reactions. In this process, it was found that the flap control system provided multi-axis control that can be utilized for entry precision targeting. Additionally, the Pterodactyl vehicle can achieve up to a trim L/D of 0.2. Finally, the increased fidelity of the aerothermal heating environments revealed that the shear stress contributes to increasing the heating on the flaps.
Deployable Entry Vehicle↗
The use of cluster computer systems for NASA/JPL applications
In this paper we outline the results of using three generations of cluster machines at JPL for a select set of applications drawn from a larger development effort.
cluster↗
Cartesian based grid generation/adaptive mesh refinement
Grid adaptation has recently received attention in the computational fluid dynamics (CFD) community as a means to capture the salient features of a flowfield by either moving grid points of a structured or by adding cells in an unstructured manner. An approach based on a background cartesian mesh is investigated from which the geometry is 'cut' out of the mesh. Once the mesh is obtained, a solution on this coarse grid is found, that indicates which cells need to be refined. This process of refining/solving continues until the flow is grid refined in terms of a user specified global parameter (such as drag coefficient etc.). The advantages of this approach are twofold: the generation of the base grid is independent of the topology of the bodies or surfaces around/through which the flow is to be computed, and the resulting grid (in uncut regions) is highly isotropic, so that the truncation error is low. The flow solver (which, along with the grid generation is still under development) uses a completely unstructured data base, and is a finite volume, upwinding scheme. Current and future work will address generating Navier-Stokes suitable grids by using locally aligned and normal face/cell refining. The attached plot shows a simple grid about two turbine blades.
Coirier, William J.↗
Self-Avoiding Walks over Adaptive Triangular Grids
In this paper, we present a new approach to constructing a "self-avoiding" walk through a triangular mesh. Unlike the popular approach of visiting mesh elements using space-filling curves which is based on a geometric embedding, our approach is combinatorial in the sense that it uses the mesh connectivity only. We present an algorithm for constructing a self-avoiding walk which can be applied to any unstructured triangular mesh. The complexity of the algorithm is O(n x log(n)), where n is the number of triangles in the mesh. We show that for hierarchical adaptive meshes, the algorithm can be easily parallelized by taking advantage of the regularity of the refinement rules. The proposed approach should be very useful in the run-time partitioning and load balancing of adaptive unstructured grids.
Heber, Gerd↗
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↗