Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mesh Adaptivity”

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 1,153 records · Page 64

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.↗

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↗

A high-order finite-volume method for hyperbolic conservation laws on locally-refined grids

We present a fourth-order accurate finite-volume method for solving time-dependent hyperbolic systems of conservation laws on Cartesian grids with multiple levels of refinement. The underlying method is a generalization of that in [5] to nonlinear systems, and is based on using fourth-order accurate quadratures for computing fluxes on faces, combined with fourth-order accurate Runge?Kutta discretization in time. To interpolate boundary conditions at refinement boundaries, we interpolate in time in a manner consistent with the individual stages of the Runge-Kutta method, and interpolate in space by solving a least-squares problem over a neighborhood of each target cell for the coefficients of a cubic polynomial. The method also uses a variation on the extremum-preserving limiter in [8], as well as slope flattening and a fourth-order accurate artificial viscosity for strong shocks. We show that the resulting method is fourth-order accurate for smooth solutions, and is robust in the presence of complex combinations of shocks and smooth flows.

97 MATHEMATICS AND COMPUTING↗

Detailed Simulations of Shock-Bifurcation and Ignition of an Argon-diluted Hydrogen/Oxygen Mixture in a Shock Tube

Detailed simulations of the bifurcation and ignition of an Argon-diluted Hydrogen/Oxygen mixture in the two-stage weak ignition regime are performed. An adaptive mesh-refinement (AMR) technique is employed to resolve all relevant physical scales that are associated with the viscous boundary-layer, the reaction front, and the shock-wave. A high-order hybrid WENO/central-differencing method is used as spatial discretization scheme, and a detailed chemical mechanism is employed to describe the combustion of the H2/O2 mixture. The operating conditions considered in this study are p = 5 bar and T = 1100 K, and fall in the third explosion limit. The computations show that the mixing of the thermally stratified fluid, carrying different momentum and enthalpy, introduces inhomogeneities in the core-region behind the reflected shock. These inhomogeneities act as localized ignition kernels. During the induction period, these kernels slowly expand and eventually transition to a detonation wave that rapidly consumes the unburned mixture.In competition with this detonation wave are the presence of secondary ignition kernels that appear in the unreacted core-region between reflected shock and detonation wave.

Ihme, Matthias↗

Characterizing Biomass Feedstock Transport Properties Using State of the Art Imaging and Computational Techniques

The microstructure of lignocellulosic biomass determines heat and mass transfer during conversion processes. We present a novel method for characterizing the transport properties of biomass using advanced imaging and computational techniques. The microstructure of two woody feedstocks, red oak and Douglas fir, before and after pyrolysis, is revealed using X-ray computed tomography (XCT). Transport properties are calculated from the XCT images, and principal permeability tensors are calculated using an immersed boundary-based finite volume solver to model gas flow through the geometries. We observe that the permeabilities of native biomass are distinctly anisotropic, however, this anisotropy is greatly reduced after pyrolysis.

adaptive mesh refinement↗

Direct Numerical Simulation of Flame-Wall Interaction for Low-Carbon Gas Turbine Combustion

Green hydrogen (H2) and ammonia (NH3) are emerging as carbon free alternatives to hydrocarbon fuels. In gas turbines, flames interact with the combustor liner (wall), which affects the pollutant emissions, the burning efficiency, and the thermal load on the liner. We lack understanding of this important flame-wall interaction (FWI) for alternative fuels. FWI occurs at the scale of the flame thickness, and such scales are only resolved by direct numerical simulation. Here, we resolve the FWI for NH3 and H2 flames. Preliminary results of two distinct projects are presented: 2D laminar anchored V-flame, and 3D turbulent swirling flame. For the laminar flame, we show flame quenching at the wall and resulting leakage of NH3 and increased N2O, but decreased NO emissions. For the swirling flame, we show CO emissions, and comparison of flame shapes and location to experimental measurements. Ultimately we show that quenching and pollutant emissions strongly affected by fuel/air ratio, and that an inclined wall can stabilize very lean turbulent CH4/H2 flames, but prone to large CO emissions.

adaptive mesh refinement↗

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↗