Search NASA⌕ Search

SEARCH · Search NASA

Results for “mesh data structures”

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 379 records · Page 21

Geomechanical Simulation of Big Hill Strategic Petroleum Reserve (Model Calibration)

A finite element numerical analysis model has been constructed that consists of a realistic mesh capturing the geometries of Big Hill (BH) Strategic Petroleum Reserve (SPR) site using the multi-mechanism deformation (M-D) salt constitutive model and including data taken daily of the wellhead pressure and level of the oil-brine interface. The salt creep rate is not uniform in the salt dome, and creep test data for BH salt is limited. Therefore, a model calibration is necessary to simulate the geomechanical behavior of the salt dome. Cavern volumetric closures of SPR caverns calculated from sonar survey reports are used for the field baseline measurement. The structure factor, A 2 , and transient strain limit factor, K o , in the M-D constitutive model are used for model calibration. An A 2 value obtained experimentally from the BH salt and K o value of WIPP salt are used as the baseline values. To adjust the magnitude of A 2 and K o , multiplication factors A2F and KOF are defined, respectively. The A2F and KOF values of the salt dome and salt drawdown layer of elements surrounding each SPR cavern have been determined through a number of back fitting analyses. The trendlines of the predictions and sonar data match up well for BH 101, 103, 104, 106, 110, 111, and 113. The prediction curves are close to the sonar data for BH 102 and 114. However, the prediction curves for BH 105, 107, 108, 109, and 112 are not close to the sonar data. An inconsistency was found in the sonar data, i.e. the sonar measurements of cavern volumes increase with time, during some periods for BH 101, 104, 106, 107, and 112. A follow-up report in 2019 will provide a resolution for these issues.

02 PETROLEUM↗

CAPRI (Computational Analysis PRogramming Interface): A Solid Modeling Based Infra-Structure for Engineering Analysis and Design Simulations

CAPRI is a CAD-vendor neutral application programming interface designed for the construction of analysis and design systems. By allowing access to the geometry from within all modules (grid generators, solvers and post-processors) such tasks as meshing on the actual surfaces, node enrichment by solvers and defining which mesh faces are boundaries (for the solver and visualization system) become simpler. The overall reliance on file 'standards' is minimized. This 'Geometry Centric' approach makes multi-physics (multi-disciplinary) analysis codes much easier to build. By using the shared (coupled) surface as the foundation, CAPRI provides a single call to interpolate grid-node based data from the surface discretization in one volume to another. Finally, design systems are possible where the results can be brought back into the CAD system (and therefore manufactured) because all geometry construction and modification are performed using the CAD system's geometry kernel.

Haimes, Robert↗

Volumetric Rendering on Wavelet-Based Adaptive Grid

Numerical modeling of physical phenomena frequently involves processes across a wide range of spatial and temporal scales. In the last two decades, the advancements in wavelet-based numerical methodologies to solve partial differential equations, combined with the unique properties of wavelet analysis to resolve localized structures of the solution on dynamically adaptive computational meshes, make it feasible to perform large-scale numerical simulations of a variety of physical systems on a dynamically adaptive computational mesh that changes both in space and time. Volumetric visualization of the solution is an essential part of scientific computing, yet the existing volumetric visualization techniques do not take full advantage of multi-resolution wavelet analysis and are not fully tailored for visualization of a compressed solution on the wavelet-based adaptive computational mesh. Our objective is to explore the alternatives for the visualization of time-dependent data on space-time varying adaptive mesh using volume rendering while capitalizing on the available sparse data representation. Two alternative formulations are explored. The first one is based on volumetric ray casting of multi-scale datasets in wavelet space. Rather than working with the wavelets at the finest possible resolution, a partial inverse wavelet transform is performed as a preprocessing step to obtain scaling functions on a uniform grid at a user-prescribed resolution. As a result, a solution in physical space is represented by a superposition of scaling functions on a coarse regular grid and wavelets on an adaptive mesh. An efficient and accurate ray casting algorithm is based just on these coarse scaling functions. Additional details are added during the ray tracing by taking an appropriate number of wavelets into account based on support overlap with the interpolation point, wavelet coefficient magnitude, and other characteristics, such as opacity accumulation (front to back ordering) and deviation from frontal viewing direction. The second approach is based on complementing of wavelet-based adaptive mesh to the traditional Adaptive Mesh Refinement (AMR) mesh. Both algorithms are illustrated and compared to the existing volume visualization software for Rayleigh-Benard thermal convection and electron density data sets in terms of rendering time and visual quality for different data compression of both wavelet-based and AMR adaptive meshes.

Vezolainen, Alexei V.↗

Technology requirements for large flexible space structures

Research, test, and demonstration experiments necessary for establishing a data base that will permit construction of large, lightweight flexible space structures meeting on-orbit pointing and surface precesion criteria are discussed. Attention is focused on the wrap-rib proof-of-concept antenna structures developed from technology used on the ATS-6 satellite. The target structure will be up to 150 m in diameter or smaller, operate at RF levels, be amenable to packaging for carriage in the Shuttle bay, be capable of being ground-tested, and permit on-orbit deployment and retraction. Graphite/epoxy has been chosen as the antenna ribs material, and the antenna mesh will be gold-plated Mo wire. A 55-m diam reflector was built as proof-of-concept with ground-test capability. Tests will proceed on components, a model, the entire structure, and in-flight. An analytical model has been formulated to characterize the antenna's thermal behavior. The flight test of the 55-m prototype in-orbit offers the chance to validate the analytical model and characterize the control, mechanical, and thermal characteristics of the antenna configuration.

Wada, B. K.↗

CFD Evaluation of Aviation Fuels for Commercial Supersonics Technology

An overview is provided of a CFD analysis of a Lean-Direct Injection (LDI) combustion concept for NASA’s Commercial Supersonic Transport (CST) program. The National Combustion Code (OpenNCC) was used to perform two-phase reacting flow computations with three different aviation fuels for UTRC’s Axially Controlled Stoichiometry (ACS) combustor at CST cruise conditions. All computations were performed with a consistent approach of mesh-generation, spray and chemical-kinetics modeling for the combustor. Emissions (EICO, EINOx) and performance characteristics were predicted with OpenNCC, to evaluate the impact of replacing Jet-A with RP2 and Iso-Paraffinic Kerosene (IPK) for supersonic combustor applications. The CFD predictions of flame structure for IPK fuel were significantly different in comparison to those of Jet-A and RP2 fuel. However, the CFD analysis predicted very similar EINOx emissions for Jet-A, RP-2 and IPK fuels. The predicted EINOx emissions for all three fuels also compared very well with measured experimental data.

gas turbine combustion↗

MESHMAKER (MM) V1.5

MeshMaker v1.5 is a code that describes the system geometry and discretizes the domain in problems of flow and transport through porous and fractured media that are simulated using the TOUGH+ [Moridis and Pruess, 2014] or TOUGH2 [Pruess et al., 1999; 2012] families of codes. It is a significantly modified and drastically enhanced version of an earlier simpler facility that was embedded in the TOUGH2 codes [Pruess et al., 1999; 2012], from which it could not be separated. The code (MeshMaker.f90) is a stand-alone product written in FORTRAN 95/2003, is written according to the tenets of Object-Oriented Programming, has a modular structure and can perform a number of mesh generation and processing operations. It can generate two-dimensional radially symmetric (r,z) meshes, and one-, two-, and three-dimensional rectilinear (Cartesian) grids in (x,y,z). The code generates the file MESH, which includes all the elements and connections that describe the discretized simulation domain and conforming to the requirements of the TOUGH+ and TOUGH2 codes. Multiple-porosity processing for simulation of flow in naturally fractured reservoirs can be invoked by means of a keyword MINC, which stands for Multiple INteracting Continua. The MINC process operates on the data of the primary (porous medium) mesh as provided on disk file MESH, and generates a secondary mesh containing fracture and matrix elements with identical data formats on file MINC.

MORIDIS, GEORGE↗

Selected computations of transonic cavity flows

An efficient diagonal scheme implemented in an overset mesh framework has permitted the analysis of geometrically complex cavity flows via the Reynolds averaged Navier-Stokes equations. Use of rapid hyperbolic and algebraic grid methods has allowed simple specification of critical turbulent regions with an algebraic turbulence model. Comparisons between numerical and experimental results are made in two dimensions for the following problems: a backward-facing step; a resonating cavity; and two quieted cavity configurations. In three-dimensions the flow about three early concepts of the stratospheric Observatory For Infrared Astronomy (SOFIA) are compared to wind-tunnel data. Shedding frequencies of resolved shear layer structures are compared against experiment for the quieted cavities. The results demonstrate the progress of computational assessment of configuration safety and performance.

Atwood, Christopher A.↗

Establishing Best Practices for X-57 Maxwell CFD Database Generation

The X-57 Maxwell is NASA’s latest electric airplane concept that has been simulated for aerodynamic performance using the structured overset and unstructured grid solvers within the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework as well as the unstructured polyhedral grid solver in Star-CCM+ for code-to-code comparison. In order to validate the predictions, comparisons were made between the CFD solutions and experimental data collected in the 12-foot Low-Speed Wind Tunnel at NASA Langley Research Center. The simulations are in preparation for the development of a comprehensive aerodynamic database which will assess aircraft performance at a variety of conditions. The findings from these simulations will establish the best practices for mesh resolution, numerical discretization, and turbulence modeling to be used for this database. Preliminary database results have shown that best-practices learned from the initial validation simulations will potentially reduce error in X-57 aerodynamic loads and moments relative to experiment by up to 14%.

Duensing, Jared C.↗

Complex Fault Geometry of the 2020 M ww 6.5 Monte Cristo Range, Nevada, Earthquake Sequence

On 15 May 2020 an M ww 6.5 earthquake occurred beneath the Monte Cristo Range in the Mina Deflection region of western Nevada. Rapid deployment of eight temporary seismic stations enabled detailed analysis of its productive and slowly decaying aftershock sequence (p=0.8), which included ~18,000 autodetected events in 3.5 months. Double-difference, waveform-based relative relocation of 16,714 earthquakes reveals a complex network of faults, many of which cross the inferred 35-km-long east–northeast-striking, left-lateral mainshock rupture. Seismicity aligns with left-lateral, right-lateral, and normal mechanism moment tensors of 128 of the largest earthquakes. The mainshock occurred near the middle of the aftershock zone at the intersection of two distinct zones of seismicity. In the western section, numerous subparallel, shallow, north-northeast-striking faults form a broad flower-structure-like fault mesh that coalesces at depth into a near-vertical, left-lateral fault. We infer the near-vertical fault to be a region of significant slip in the mainshock and an eastward extension of the left-lateral Candelaria fault. Near the mainshock hypocenter, seismicity occurs on a northeast-striking, west-dipping structure that extends north from the eastern Columbus Salt Marsh normal fault. Together, these two intersecting structures bound the Columbus Salt Marsh tectonic basin. East of this intersection and the mainshock hypocenter, seismicity occurs in a narrow, near-vertical, east-northeast-striking fault zone through to its eastern terminus. At the eastern end, the aftershock zone broadens and extends northwest toward the southern extension of the northwest-striking, right-lateral Petrified Springs fault system. The eastern section hosts significantly fewer aftershocks than the western section, but has more moment release. We infer that shallow aftershocks throughout the system highlight fault-fracture meshes that connect mapped fault systems at depth. Here, comparing earthquake data with surface ruptures and a simple geodetic fault model sheds light on the complexity of this recent M 6.5 Walker Lane earthquake.

58 GEOSCIENCES↗

Aeroelastic analysis of wings using the Euler equations with a deforming mesh

Modifications to the CFL3D three dimensional unsteady Euler/Navier-Stokes code for the aeroelastic analysis of wings are described. The modifications involve including a deforming mesh capability which can move the mesh to continuously conform to the instantaneous shape of the aeroelastically deforming wing, and including the structural equations of motion for their simultaneous time-integration with the governing flow equations. Calculations were performed using the Euler equations to verify the modifications to the code and as a first step toward aeroelastic analysis using the Navier-Stokes equations. Results are presented for the NACA 0012 airfoil and a 45 deg sweptback wing to demonstrate applications of CFL3D for generalized force computations and aeroelastic analysis. Comparisons are made with published Euler results for the NACA 0012 airfoil and with experimental flutter data for the 45 deg sweptback wing to assess the accuracy of the present capability. These comparisons show good agreement and, thus, the CFL3D code may be used with confidence for aeroelastic analysis of wings.

Robinson, Brian A.↗

Aeroelastic analysis of wings using the Euler equations with a deforming mesh

Modifications to the CFL3D three-dimensional unsteady Euler/Navier-Stokes code for the aeroelastic analysis of wings are described. The modifications involve including a deforming mesh capability which can move the mesh to continuously conform to the instantaneous shape of the aeroelastically deforming wing, and including the structural equations of motion for their simultaneous time-integration with the governing flow equations. Calculations were performed using the Euler equations to verify the modifications to the code and as a first-step toward aeroelastic analysis using the Navier-Stokes equations. Results are presented for the NACA 0012 airfoil and a 45 deg sweptback wing to demonstrate applications of CFL3D for generalized force computations and aeroelastic analysis. Comparisons are made with published Euler results for the NACA 0012 airfoil and with experimental flutter data for the 45 deg sweptback wing to assess the accuracy of the present capability. These comparisons show good agreement and, thus, the CFL3D code may be used with confidence for aeroelastic analysis of wings. The paper describes the modifications that were made to the code and presents results and comparisons which assess the capability.

Robinson, Brian A.↗

Construction of Hexahedral Elements Mesh Capturing Realistic Geometries of Bayou Choctaw SPR Site

The three-dimensional finite element mesh capturing realistic geometries of Bayou Choctaw site has been constructed using the sonar and seismic survey data obtained from the field. The mesh is consisting of hexahedral elements because the salt constitutive model is coded using hexahedral elements. Various ideas and techniques to construct finite element mesh capturing artificially and naturally formed geometries are provided. The techniques to reduce the number of elements as much as possible to save on computer run time with maintaining the computational accuracy is also introduced. The steps and methodologies could be applied to construct the meshes of Big Hill, Bryan Mound, and West Hackberry strategic petroleum reserve sites. The methodology could be applied to the complicated shape masses for not only various civil and geological structures but also biological applications such as artificial limbs.

42 ENGINEERING↗

Learning Distribution Grid Topologies: A Tutorial

Unveiling feeder topologies from data is of paramount importance to advance situational awareness and proper utilization of smart resources in power distribution grids. This tutorial summarizes, contrasts, and establishes useful links between recent works on topology identification and detection schemes that have been proposed for power distribution grids. The primary focus is to highlight methods that overcome the limited availability of measurement devices in distribution grids, while enhancing topology estimates using conservation laws of power-flow physics and structural properties of feeders. Grid data from phasor measurement units or smart meters can be collected either passively in the traditional way, or actively, upon actuating grid resources and measuring the feeder's voltage response. Analytical claims on feeder identifiability and detectability are reviewed under disparate meter placement scenarios. Such topology learning claims can be attained exactly or approximately so via algorithmic solutions with various levels of computational complexity, ranging from least-squares fits to convex optimization problems, and from polynomial-time searches over graphs to mixed-integer programs. Although the emphasis is on radial single-phase feeders, extensions to meshed and/or multiphase circuits are sometimes possible and discussed. Here this tutorial aspires to provide researchers and engineers with knowledge of the current state-of-the-art in tractable distribution grid learning and insights into future directions of work.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Analytical and Experimental Vibration Analysis of a Faulty Gear System

A comprehensive analytical procedure was developed for predicting faults in gear transmission systems under normal operating conditions. A gear tooth fault model is developed to simulate the effects of pitting and wear on the vibration signal under normal operating conditions. The model uses changes in the gear mesh stiffness to simulate the effects of gear tooth faults. The overall dynamics of the gear transmission system is evaluated by coupling the dynamics of each individual gear-rotor system through gear mesh forces generated between each gear-rotor system and the bearing forces generated between the rotor and the gearbox structure. The predicted results were compared with experimental results obtained from a spiral bevel gear fatigue test rig at NASA Lewis Research Center. The Wigner-Ville distribution (WVD) was used to give a comprehensive comparison of the predicted and experimental results. The WVD method applied to the experimental results were also compared to other fault detection techniques to verify the WVD's ability to detect the pitting damage, and to determine its relative performance. Overall results show good correlation between the experimental vibration data of the damaged test gear and the predicted vibration from the model with simulated gear tooth pitting damage. Results also verified that the WVD method can successfully detect and locate gear tooth wear and pitting damage.

Choy, F. K.↗

A Database of CFD-Based Buffet Forcing Functions for Artemis I Structural Response Evaluation

Time-accurate FUN3D simulations are utilized to estimate buffet-induced unsteady forces experienced by the Space Launch System during the Artemis I (AR01) flight. A set of FUN3D simulations was developed that employed time-accurate mesh translations to simulate the changing velocity and attitude based on the AR01 best estimated trajectory. In these simulations, referred to as accelerating-flow simulations, the freestream Mach number increased from 0.80 to 1.92. Additional time-accurate simulations were obtained at constant freestream Mach number equal to 0.95, 1.18, 1.70, thereby simulating stationary conditions experienced by the flow in a wind-tunnel. On the grounds of favorable comparisons between simulated and flight-measured environments, the FUN3D-based surface pressures were utilized to develop a buffet forcing function (BFF) database. This BFF database was analyzed to characterize the spatial distribution and frequency content of the buffet forces during transonic and supersonic portions of the AR01 flight. The region of interest is located downstream of the forward attachment (FA) hardware between the core stage and the solid rocket boosters where vortex shedding off the FA protuberance produces significant unsteadiness. The analysis reveals that, at transonic and supersonic conditions, buffet forces that are based on constant freestream Mach number data are a good approximation of those based on accelerating-flow simulations.

transonic buffet↗

A Database of CFD-Based Buffet Forcing Functions for Artemis I Structural Response Evaluation

Time-accurate FUN3D simulations are utilized to estimate buffet-induced unsteady forces experienced by the Space Launch System during the Artemis I (AR01) flight. A set of FUN3D simulations was developed that employed time-accurate mesh translations to simulate the changing velocity and attitude based on the AR01 best estimated trajectory. In these simulations, referred to as accelerating-flow simulations, the freestream Mach number increased from 0.80 to 1.92. Additional time-accurate simulations were obtained at constant freestream Mach number equal to 0.95, 1.18, 1.70, thereby simulating stationary conditions experienced by the flow in a wind-tunnel. On the grounds of favorable comparisons between simulated and flight-measured environments, the FUN3D-based surface pressures were utilized to develop a buffet forcing function (BFF) database. This BFF database was analyzed to characterize the spatial distribution and frequency content of the buffet forces during transonic and supersonic portions of the AR01 flight. The region of interest is located downstream of the forward attachment (FA) hardware between the core stage and the solid rocket boosters where vortex shedding off the FA protuberance produces significant unsteadiness. The analysis reveals that, at transonic and supersonic conditions, buffet forces that are based on constant freestream Mach number data are a good approximation of those based on accelerating-flow simulations.

FUN3D↗

Reduction of rotor-turbulence interaction noise in static fan noise testing

Three effective inlet cleanup methods are described. A flared reverse cone inlet is used to eliminate wakes from the fan casings and/or probe supports. Boundary layer suction is employed ahead of the fan rotor and on the outer flare of the cone to reduce the boundary layer turbulence and remove any residual wakes. To reduce the midstream turbulent intensity and length scales, a turbulence control structure is used which is constructed with both a layer of honeycomb and a fine mesh screen. The effects of these cleanup methods are quantified by measuring the far-field noise in an anechoic chamber, using a high-speed 20-in.-diam fan of the current high passband type. The changes in the turbulent field impinging on this rotor are quantified by mapping the streamwise and transverse turbulent properties (spectra, intensity and length scale) with crossed hot film probes. It is convincingly shown that it is possible to clean up the inlet flow of a static fan noise test facility to a point where the static acoustic data simulate flight data.

Kantola, R. A.↗

The NASA/industry Design Analysis Methods for Vibrations (DAMVIBS) program: Boeing Helicopters airframe finite element modeling

Mathematical models based on the finite element method of structural analysis, as embodied in the NASTRAN computer code, are routinely used by the helicopter industry to calculate airframe static internal loads used for sizing structural members. Historically, less reliance has been placed on the vibration predictions based on these models. Beginning in the early 1980's NASA's Langley Research Center initiated an industry wide program with the objective of engendering the needed trust in vibration predictions using these models and establishing a body of modeling guides which would enable confident future prediction of airframe vibration as part of the regular design process. Emphasis in this paper is placed on the successful modeling of the Army/Boeing CH-47D which showed reasonable correlation with test data. A principal finding indicates that improved dynamic analysis requires greater attention to detail and perhaps a finer mesh, especially the mass distribution, than the usual stress model. Post program modeling efforts show improved correlation placing key modal frequencies in the b/rev range with 4 percent of the test frequencies.

Gabel, R.↗