Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Complex geometry”

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 55 records · Page 3

Analysis of GPU Data Access Patterns on Complex Geometries for the D3Q19 Lattice Boltzmann Algorithm

GPU performance of the lattice Boltzmann method (LBM) depends heavily on memory access patterns. When implemented with GPUs on complex domains, typically, geometric data is accessed indirectly and lattice data is accessed lexicographically. Although there are a variety of other options, no study has examined the relative efficacy between them. Here, we examine a suite of memory access schemes via empirical testing and performance modeling. We find strong evidence that semi-direct is often better suited than the more common indirect addressing, providing increased computational speed and reducing memory consumption. For the layout, we find that the Collected Structure of Arrays (CSoA) and bundling layouts outperform the common Structure of Array layout; on V100 and P100 devices, CSoA consistently outperforms bundling, however the relationship is more complicated on K40 devices. When compared to state-of-the-art practices, our recommendations lead to speedups of 10–40 percent and reduce memory consumption up to 17 percent. Using performance modeling and computational experimentation, we determine the mechanisms behind the accelerations. We demonstrate that our results hold across multiple GPUs on two leadership class systems, and present the first near-optimal strong results for LBM with arterial geometries run on GPUs.

42 ENGINEERING↗

An isoparametric spectral element method for solution of the Navier-Stokes equations in complex geometry

High-order (p-type) finite element methods combine the geometric flexibility of standard low-order finite element (or finite volume) techniques with the rapid convergence properties of spectral methods. Various p-type schemes have been proposed for elliptic problems. However, an application of p-type methods to more complex equations, in particular the passive scalar or Navier-Stokes equations, is complicated by the nature of the equations and by efficiency considerations as regards the solution procedure. Patera (1984) has proposed a p-type method for the Navier-Stokes equations, taking into account the arising additional complications. The present study is concerned with an isoparametric spectral element discretization for the solution of flow problems in fully general, curvy geometries. Attention is given to the isoparametric formulation for elliptic equations and a time-splitting scheme for the Navier-Stokes equations.

Korczak, K. Z.↗

Generation of unstructured grids and Euler solutions for complex geometries

Algorithms are described for the generation and adaptation of unstructured grids in two and three dimensions, as well as Euler solvers for unstructured grids. The main purpose is to demonstrate how unstructured grids may be employed advantageously for the economic simulation of both geometrically as well as physically complex flow fields.

Loehner, Rainald↗

Counter-Current Flow Limitation Studies in Complex Geometries Utilizing Interface Capturing Simulations Coupled with PID Flow Rate Controller

In nuclear thermal-hydraulic studies, counter-current flow limitation (CCFL) typically refers to steam rising at a fast rate such that it prevents coolant from draining down within a confined channel. CCFL is a crucial issue in nuclear reactor safety analysis. This study investigates CCFL in debris bed channels using high-resolution interface-capturing simulations. A novel proportional-integral-derivative flow rate controller is developed to efficiently achieve the CCFL conditions. Verification studies confirm that CCFL occurs under the same conditions with or without the controller, demonstrating that PID control ensures accurate prediction. Three debris bed channel geometries were examined: a cylindrical channel, a channel with small obstacles, and a channel with large obstacles. Results show that obstacles significantly impact flow behavior, interfacial shear, wall shear, and pressure gradients required for CCFL. Furthermore, the comparison with experimental data confirmed that simulations incorporating geometric complexities align more closely with experimental CCFL conditions. A pressure gradient correlation was also developed for CCFL prediction.

Counter-current flow limitation↗

Detection of Micro-Leaks Through Complex Geometries Under Mechanical Load and at Cryogenic Temperature

Polymer Matrix Composite (PMC) hydrogen tanks have been proposed as an enabling technology for reducing the weight of Single-Stage-to-Orbit reusable launch vehicles where structural mass has a large impact on vehicle performance. A key development issue of these lightweight structures is the leakage of hydrogen through the composite material. The rate of hydrogen leakage can be a function of the material used, method of 6 fabrication used to manufacture the tank, mechanical load the tank must react, internal damage-state of the material, and the temperatures at which the tank must operate. A method for measuring leakage through a geometrically complex structure at cryogenic temperature and under mechanical load was developed, calibrated and used to measure hydrogen leakage through complex X-33 liquid-hydrogen tank structure sections.

Rivers, H. Kevin↗

A stochastic scan strategy for grain structure control in complex geometries using electron beam powder bed fusion

Spatial control of microstructure within a three-dimensional component has been a dream of materials scientists for centuries. However, limitations in traditional manufacturing processes prevent detailed control over the distribution of microstructures in a single part. Here, we demonstrate the ability to control grain structure and crystallographic texture during metal additive manufacturing for arbitrary cross-sections of a practical size, with profound implications for the design and optimization of next-generation products. The key to this advance is a new geometry agnostic scan path algorithm that manipulates the spatial distribution of solidification conditions. Utilizing a fundamental understanding of solidification dynamics and a model of the heat transfer during processing, we have designed this algorithm to manipulate the natural competition between epitaxial dendrite growth and grain nucleation. With this algorithm, we successfully controlled the grain structure of Ni-based superalloy IN718 in the shape of the Mona Lisa.

36 MATERIALS SCIENCE↗

Internal 2-D Surface Temperature Measurements for Large Complex Geometries

For 2D-temperature monitoring applications, a variant of EIT (Electrical Impedance Tomography) is evaluated computationally in this work. Literature examples of poor sensor performance in the center of the 2D domains away from the side electrodes motivated this study which seeks to overcome some of the previously noted shortcomings. In particular, the use of ‘sensing skins’ with novel tailored baseline conductivities were examined using the EIDORS package for EIT. It was found that the best approach for detecting a hot spot depends on several factors such as the current injection (stimulation) patterns, the measurement patterns, and the reconstruction algorithms. For a well-performing combination of these factors, tailored baseline conductivities were assessed and compared to the baseline uniform conductivity. It was discovered that for some EIT applications, a tailored distribution needs to be smooth and that sudden changes in the conductivity gradients should be avoided. Still, the benefits in terms of improved EIT performance were small for conditions for which the EIT measurements had been ‘optimized’ for the uniform baseline case. Within the limited scope of this study, only two specific cases showed benefits from tailored distributions. For one case, a smooth tailored distribution with increased baseline conductivity in the center provided a better separation of two centrally located hot spots. For another case, a smooth tailored distribution with reduced conductivity in the center provided better estimates of the magnitudes of two hot spots near the center of the sensing skin.

47 OTHER INSTRUMENTATION↗

Spectral methods for problems in complex geometries

Techniques that permit the efficient application of spectral methods to solve problems in nearly arbitrary geometries are presented. These methods were found to be viable alternatives to finite difference and finite element processes. The spectral methods applied are extensions of the standard techniques of separation of variables to the solution of arbitrarily complicated problems.

Orszag, S. A.↗

An interactive approach to surface-fitting complex geometries for flowfield applications

Numerical flowfield methods require a geometry subprogram which can calculate body coordinates, slopes, and radii of curvature for typical aircraft and spacecraft configurations. The objective of this paper is to develop a new surface-fitting technique which addresses two major problems with existing geometry packages: computer storage requirements and the time required of the user for the initial set-up of the geometry model. In the present method, coordinates of cross sections are fit in a least-squares sense using segments of general conic sections. After fitting each cross section, the next step is to blend the cross-sectional curve-fits in the longitudinal direction using general conics to fit specific meridional half-planes. For the initial setup of the geometry model, an interactive, completely menu-driven computer code has been developed to allow the user to make modifications to the initial fit for a given cross section or meridional cut. Graphic displays are provided to assist the user in the visualization of the effect of each modification. The completed model may be viewed from any angle using the code's three-dimensional graphics package. Geometry results for the modeling of the Space Shuttle and a proposed Aeroassist Flight Experiment (AFE) geometry are presented, in addition to calculated heat-transfer rates based on these models.

Dejarnette, Fred R.↗

A computer code for multiphase all-speed transient flows in complex geometries. MAST version 1.0

The operation of the MAST code, which computes transient solutions to the multiphase flow equations applicable to all-speed flows, is described. Two-phase flows are formulated based on the Eulerian-Lagrange scheme in which the continuous phase is described by the Navier-Stokes equation (or Reynolds equations for turbulent flows). Dispersed phase is formulated by a Lagrangian tracking scheme. The numerical solution algorithms utilized for fluid flows is a newly developed pressure-implicit algorithm based on the operator-splitting technique in generalized nonorthogonal coordinates. This operator split allows separate operation on each of the variable fields to handle pressure-velocity coupling. The obtained pressure correction equation has the hyperbolic nature and is effective for Mach numbers ranging from the incompressible limit to supersonic flow regimes. The present code adopts a nonstaggered grid arrangement; thus, the velocity components and other dependent variables are collocated at the same grid. A sequence of benchmark-quality problems, including incompressible, subsonic, transonic, supersonic, gas-droplet two-phase flows, as well as spray-combustion problems, were performed to demonstrate the robustness and accuracy of the present code.

Chen, C. P.↗

A Numerical Method for Direct Simulation of Turbulence in Complex Geometries

The ultimate goal of this work is to study the flow inside a channel with riblets on one of the two walls. The method has been tested for two dimensional flows in the presence of bodies with a geometrical singularity and for three dimensional flows inside domains described by Cartesian coordinates. The results have been compared with previous numerical simulations and with experimental results. The cases considered are: (1) the growth of Orr-Sommerfield waves in plane Poiseulle flow; (2) the flow over a backward facing step; (3) the flow past a wedge; and (4) the flow inside a narrow channel. Finally, the case of a channel with two large riblets on a wall has been simulated. In this case a limited number of grid points is sufficient, and in spite of the slow convergence for the pressure solver, one is able to obtain solutions with a reasonable amount of computer time. At present solutions with very fine grids in all three directions can not be obtained, due to lack of a fast pressure solver for general curvilinear coordinates.

P Orlandi↗