Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mesh Refinement”

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 451 records · Page 25

Finite element Euler computations in three dimensions

A two-step explicit FEM solution algorithm for the three-dimensional compressible Euler and Navier-Stokes equations based on unstructured triangular and tetrahedral grids is described and demonstrated. The method represents an extension and refinement of the algorithms presented by Loehner et al. (1984 and 1985), Peraire et al. (1987), and Morgan et al. (1987). The formulation and numerical implementation are outlined; the mesh generation, data structures, and adaptive remeshing are explained; and results for a two-dimensional airfoil, a three-dimensional engine air intake, a B747 in landing configuration, and a generic fighter aircraft are presented in extensive graphics and discussed in detail.

Peraire, Jaime↗

Grid Refinement Techniques for the 𝜸-Re 𝜽 𝒕 Transition Model in FUN3D

There has been an increased focus on the overall accuracy and grid convergence of Reynolds-averaged Navier-Stokes (RANS)-based transition models from the recent AIAA and NATO-AVT workshops. Even though satisfactory grid convergence could be achieved for simple two-dimensional flow configurations, it required mesh counts that are substantially larger than those used in typical applications. In this paper, we focus our efforts on understanding how the grid resolution and topology influences the accuracy and convergence of results by studying the Schubauer and Skramstad flat-plate configuration and the NLF-0416 airfoil at an angle of attack equal to five degrees using FUN3D, a second-order finite-volume code. By focusing on these cases, we can analyze both natural and separation-induced transition scenarios. Multiple grid refinement strategies are investigated. First, we determine the relative effectiveness of zonal streamwise refinement in the transition region within structured grids as an alternative to the costly option of globally uniform refinement of a baseline grid. We also complement this zonal technique by globally varying the wall-normal resolution keeping the streamwise resolution fixed. The zonal streamwise refinement can accurately model natural transition in a flat-plate boundary layer and separation-induced transition, but struggles to accurately model natural transition in airfoil flows. A series of unstructured prismatic grids that have similar node counts and viscous wall spacings as the structured hexahedral grids are also tested, and they do not achieve grid convergence until an extremely fine resolution. Last, we employ adjoint-based unstructured grid adaptation in FUN3D to natural and separation-induced transition on the NLF-0416 airfoil. The adjoint-based refinement process converges to the same solution as the baseline family of structured grids, but leads to smaller errors on coarser grids.

Grid Adaptation↗

Fast projection algorithm for unstructured meshes

Two projection operators are presented which employ a very efficient point-location algorithm. For a wide variety of practical problems, the projection cost was reduced by a substantial factor when compared with simpler procedures. Problems involving high levels of local refinement do not impede the point-location algorithm, and the extension of these concept to three dimensions has been implemented with no additional difficulty. Nodal interpolation is noted to be an excellent projection operator when cost is a criterion, while consistent-mass L2-projection furnishes the highest accuracy.

Jansen, Kenneth↗

Swift: A Widefield Imaging Fabry Perot for Sofia

Contract was to pursue feasibility studies of the SOFIA Widefield Imaging Fabry-Perot (SWIFT). SWIFT was proposed as a two color 18 to 40 microns imaging Fabry-Perot that utilized two Rockwell/Boeing 256 x 256 pixel Si:Sb BIBs as detective devices. The colors were to be split between 26 and 30 microns using a MgO dichroic. The resolution achieving devices were to be a pair of cryogenic fully tunable scanning Fabry-Perot interferometers (FPIs), two in each band. For high resolving powers, a third, fixed FPI is inserted into the beam. The FPI mirrors were to be made of free standing metal mesh. We also proposed to look into a long wavelength (40 to 210 microns) band during the feasibility study period. We produced a proposal to USRA, submitted in July 1997 that substantially refined our ideas. We decided the long wavelength science was compelling, so the baseline wavelength coverage for SWIFT was widened to 17 to 205 microns. Under typical operations, we proposed to simultaneously image in two bands: 22 to 38 microns, and 50 to 205 microns. The bands were to be split by a cold CaF2 dichroic. The short wavelength (SW) band was to employ a 256 x 256 pixel Boeing/Rockwell Si:Sb BIB array, and the long wavelength (LW) band was to employ a Goddard 6 x 32 (upgradable to 32 x 32) element array of monolithic silicon "pop-up" bolometers as detective devices. The two color capability doubled the data taking efficiency, and ensured "perfect" registration between the images obtained in each band. For the SW band, the beam was to be fully sampled (0.7" pixels, 1.4 in. beam) at 17 microns, and over sampled at longer wavelengths. Even so, SWIFT has a very large (3 ft x 3 ft) field of view. To match the SW and LW fields of view (initially in one dimension only, but in 2-dimensions with 32 x 32 upgrade), SWIFT was to under sample at 63 microns (5.6 in pixels, 5.2 in beam) resulting in a 0.56 x 3 in (upgrade to 3 in x 3 in) field of view. Each band has both Lo-Res (R triple bond lambda/Delta(lambda) approx. 30 to 200) and Hi-Res (R approx. 1000 up to 10(exp 5)) modes, that can be changed in less than a minute. Spectral resolution is obtained with four identical cryogenic, fully tunable, scanning FPIS, two in each band, plus a set of fixed order FPI that can be rotated into the beam for good spectral purity at the highest resolving powers. SWIFT was proposed to explore three primary topics of scientific interest: 1. Protostars, Young Stellar Objects, and Protoplanetary Disks. 2. The Galactic Center. 3. External Galaxies.

Stacey, Gordon J.↗

HLPW-5: Overview and Workshop Summary

The Fifth AIAA CFD High-Lift Prediction Workshop was held with the goal of assessing the numerical prediction capability of current-generation computational fluid dynamics (CFD) technology for swept, medium/high-aspect-ratio wings in high-lift configurations. A key aspect of this endeavor was the use of Technology Focus Groups (TFG), which included both mesh generation and flow-solver experts working together to accelerate advancements for their particular CFD methodology, by addressing key questions of importance prior to the workshop. The high-lift version of the NASA Common Research Model (CRM-HL) configuration was the focus of this workshop, and was used for three unique test cases. Wind-tunnel data were available for comparison for one of the test cases. Altogether, 365 datasets of CFD results were submitted by 47 teams, with 41 teams contributing to the multiple configurations of Case 1, 40 to Case 2, and 18 to Case 3. This paper provides a high-level summary of the results and conclusions from the workshop. As concluded from past workshops, fixed-grid Reynolds-averaged Navier-Stokes methods continued to be inaccurate and inconsistent for high-lift flows near maximum lift. However, application of mesh-adaptation technology helped to achieve improved consistency. Scale-resolving methods appeared most promising for predicting high-lift flow physics, particularly at maximum lift. Best practices for these methods were refined over the course of the workshop and new challenges were identified.

Adam M Clark↗

Topology Synthesis of Structures Using Parameter Relaxation and Geometric Refinement

Typically, structural topology optimization problems undergo relaxation of certain design parameters to allow the existence of intermediate variable optimum topologies. Relaxation permits the use of a variety of gradient-based search techniques and has been shown to guarantee the existence of optimal solutions and eliminate mesh dependencies. This Technical Publication (TP) will demonstrate the application of relaxation to a control point discretization of the design workspace for the structural topology optimization process. The control point parameterization with subdivision has been offered as an alternative to the traditional method of discretized finite element design domain. The principle of relaxation demonstrates the increased utility of the control point parameterization. One of the significant results of the relaxation process offered in this TP is that direct manufacturability of the optimized design will be maintained without the need for designer intervention or translation. In addition, it will be shown that relaxation of certain parameters may extend the range of problems that can be addressed; e.g., in permitting limited out-of-plane motion to be included in a path generation problem.

Hull, P. V.↗

Reduced Navier Stokes Relaxation Procedures for Internal Flows

In spite of significant advancement in the field of high speed computing, flow calculations involving complex geometries and/or flow behavior still require large amounts of CPU time and memory. In order to predict such flows without sacrificing grid convergence and accuracy, adaptive gridding techniques that provide optimal resolution are highly desirable. The present work combines multigrid techniques and domain decomposition concepts to provide local, solution adaptive, grid refinement. Several viscous compressible and incompressible, two and three-dimensional, flows with strong inviscid interaction and/or axial flow reversal, are considered with a segmented multigrid domain decomposition (SMGDD) procedure for which uniform meshes result in each domain. A pressure-based form of flux-vector splitting is applied to the Navier-Stokes equations, which are represented by an implicit lowest-order reduced Navier-Stokes (RNS) system and a purely diffusive, higher-order, deferred-corrector. A trapezoidal or box-like form of discretization insures that all mass conservation properties are satisfied at interfacial and outflow boundaries, even for this primitive-variable non-staggered grid computation. The SMGDD technique presented herein has previously been applied for incompressible two dimensional flows. The present work offers improvement in the gridding strategy, by allowing for disjoint subdomains that provide optimal resolution of disparate flow features. It also extends the SMGDD technique to three dimensional compressible flows. Laminar and turbulent flow in a backward facing step channel is considered; although the procedure is applicable to more severe geometries. The standard K-epsilon model is applied for turbulence closure. For Re greater than 400, differences between two-dimensional theory and experiment are resolved through a three dimensional simulation, which confirms the experimentally observed three dimensionality of the recirculation patterns on the upper and lower surfaces.

Rubin, Stanley G.↗

ISO Guest Observer Data Analysis and LWS Instrument Team Activities

This project was granted a no-cost extension prompted by the request of the major subcontractor, the Naval Research Laboratory, which had not yet completed its tasks. As of July 2002, they had made substantial progress. They have successfully fabricated a metal mesh grid on polyimide, and also successfully fabricated a 2-layer metal mesh infrared filter using stacks of these metal mesh grids on polyimide; the actual layering was done at SAO. Both warm and cold spectroscopic tests were done on these fabricated devices. The measurements were in good agreement with the theory, and also reasonable performance in absolute terms. NRL is now working on fabricating a 3-layer metal mesh infrared filter, and a prototype is expected in the next month. Testing should occur before the end of the fiscal year. Finally, NRL has preliminarily agreed to hire a new postdoctoral person to refine the modeling of the filters based on the new measurements. The person should arrive this fall. NRL has a new Fourier Transform Spectrometer which will be delivered in the next month, and which will be used to facilitate the testing which has up to now been done in collaboration with NASA Goddard Space Flight Space Center.

Smith, Howard↗

Implementation of Active Thermal Control (ATC) for the Soil Moisture Active and Passive (SMAP) Radiometer

NASA's Earth Observing Soil Moisture Active and Passive (SMAP) Mission is scheduled to launch in November 2014 into a 685 kilometer near-polar, sun-synchronous orbit. SMAP will provide comprehensive global mapping measurements of soil moisture and freeze/thaw state in order to enhance understanding of the processes that link the water, energy, and carbon cycles. The primary objectives of SMAP are to improve worldwide weather and flood forecasting, enhance climate prediction, and refine drought and agriculture monitoring during its three year mission. The SMAP instrument architecture incorporates an L-band radar and an L-band radiometer which share a common feed horn and parabolic mesh reflector. The instrument rotates about the nadir axis at approximately 15 revolutions per minute, thereby providing a conically scanning wide swath antenna beam that is capable of achieving global coverage within three days. In order to make the necessary precise surface emission measurements from space, the electronics and hardware associated with the radiometer must meet tight short-term (instantaneous and orbital) and long-term (monthly and mission) thermal stabilities. Maintaining these tight thermal stabilities is quite challenging because the sensitive electronics are located on a fast spinning platform that can either be in full sunlight or total eclipse, thus exposing them to a highly transient environment. A passive design approach was first adopted early in the design cycle as a low-cost solution. With careful thermal design efforts to cocoon and protect all sensitive components, all stability requirements were met passively. Active thermal control (ATC) was later added after the instrument Preliminary Design Review (PDR) to mitigate the threat of undetected gain glitches, not for thermal-stability reasons. Gain glitches are common problems with radiometers during missions, and one simple way to avoid gain glitches is to use the in-flight set point programmability that ATC affords to operate the radiometer component away from the problematic temperature zone. A simple ThermXL model (10 nodes) was developed to exercise quick trade studies among various proposed control algorithms: Modified P control vs. PI control. The ThermXL results were then compared with the detailed Thermal Desktop (TD) model for corroboration. Once done, the simple ThermXL model was used to evaluate parameter effects such as temperature digitization, heater size and gain margin, time step, and voltage variation of power supply on the ATC performance. A Modified P control algorithm was implemented into the instrument flight electronics based on the ThermXL results. The thermal short-term stability margin decreased by 10 percent with ATC and a wide temperature error band (plus or minus 0.1 degrees Centigrade) compared to the original passive thermal design. However, a tighter temperature error band (plus or minus 0.1 degrees Centigrade) increased the thermal short-term stability margin by a factor of three over the passive thermal design. The current ATC design provides robust thermal control, tighter stability, and greater in-flight flexibility even though its implementation was prompted by non-thermal performance concerns.

Mikhaylov, Rebecca↗

A finite-volume Euler solver for computing rotary-wing aerodynamics on unstructured meshes

An unstructured-grid solver for the unsteady Euler equations has been developed for predicting the aerodynamics of helicopter rotor blades. This flow solver is a finite-volume scheme that computes flow quantities at the vertices of the mesh. Special treatments are used for the flux differencing and boundary conditions in order to compute rotary-wing flowfields, and these are detailed in the paper. The unstructured-grid solver permits adaptive grid refinement in order to improve the resolution of flow features such as shocks, rotor wakes and acoustic waves. These capabilities are demonstrated in the paper. Example calculations are presented for two hovering rotors. In both cases, adaptive-grid refinement is used to resolve high gradients near the rotor surface and also to capture the vortical regions in the rotor wake. The computed results show good agreement with experimental results for surface airloads and wake geometry.

Strawn, Roger C.↗

Implementation of a Smeared Crack Band Model in a Micromechanics Framework

The smeared crack band theory is implemented within the generalized method of cells and high-fidelity generalized method of cells micromechanics models to capture progressive failure within the constituents of a composite material while retaining objectivity with respect to the size of the discretization elements used in the model. An repeating unit cell containing 13 randomly arranged fibers is modeled and subjected to a combination of transverse tension/compression and transverse shear loading. The implementation is verified against experimental data (where available), and an equivalent finite element model utilizing the same implementation of the crack band theory. To evaluate the performance of the crack band theory within a repeating unit cell that is more amenable to a multiscale implementation, a single fiber is modeled with generalized method of cells and high-fidelity generalized method of cells using a relatively coarse subcell mesh which is subjected to the same loading scenarios as the multiple fiber repeating unit cell. The generalized method of cells and high-fidelity generalized method of cells models are validated against a very refined finite element model.

Pineda, Evan J.↗

Progressive Failure of a Unidirectional Fiber-Reinforced Composite Using the Method of Cells: Discretization Objective Computational Results

The smeared crack band theory is implemented within the generalized method of cells and high-fidelity generalized method of cells micromechanics models to capture progressive failure within the constituents of a composite material while retaining objectivity with respect to the size of the discretization elements used in the model. An repeating unit cell containing 13 randomly arranged fibers is modeled and subjected to a combination of transverse tension/compression and transverse shear loading. The implementation is verified against experimental data (where available), and an equivalent finite element model utilizing the same implementation of the crack band theory. To evaluate the performance of the crack band theory within a repeating unit cell that is more amenable to a multiscale implementation, a single fiber is modeled with generalized method of cells and high-fidelity generalized method of cells using a relatively coarse subcell mesh which is subjected to the same loading scenarios as the multiple fiber repeating unit cell. The generalized method of cells and high-fidelity generalized method of cells models are validated against a very refined finite element model.

Pineda, Evan J.↗

The 15 meter hoop-column antenna dynamics: Test and analysis

A 15 meter model of the hoop-column antenna concept has been vibration tested for model characterization and analytical model verification. Linear finite element analysis predicted the global vibration frequencies accurately. Good agreement between analysis and test data was obtained only after the analytical model was refined using static test data. As structures become more flexible, structural properties determined from static data become more accurate and should be used to update analytical models. Global vibration modes are not significantly affected by the surface mesh which permits simplified analytical models to be used for prediction of global behavior. These reduced models are believed sufficient for preliminary design and controls simulations where only global behavior is desired. The mesh modes were highly damped due to the knit mesh used for the reflector surface. These modes were also highly coupled and very difficult to measure in the laboratory. The inability to fully characterize the antenna mesh modes in the laboratory indicates robust methods for active surface vibration suppression will be needed. Fortunately, the surface mesh exhibits high passive damping which should be beneficial to active control systems.

Belvin, W. Keith↗

Gateway Element and Payload Materials Outgassing Analyses: HALO, HERMES, and ERSA

Gateway was intended to be humanity’s first space station around the Moon, but its development has been paused as the National Aeronautics and Space Administration (NASA) shifts focus to achieving the United States’ National Space Policy goals. Instead of an orbiting lunar outpost, NASA will now pursue the development of a lunar surface base to support a sustained human presence on the Moon. Before the program’s pause, Gateway’s Induced Environments team worked to ensure payloads and elements (i.e., modules) complied with induced environment requirements. Methods developed and insights gained from this work will have applicability to NASA’s Moon Base and the potential repurposing of Gateway elements and payloads, as well as to induced environments modeling for future space stations. The Gateway program’s induced environment included molecular contamination, electric thruster plume sputter and redeposition, and lunar dust transfer from the Human Landing System (HLS). Primary sources of external molecular contamination included materials outgassing, chemical thruster plume contamination, and vacuum venting. The focus of this paper will be on element- and payload-level materials outgassing analyses performed for Gateway Configuration 1, extending the previously-developed framework for Gateway system-level external molecular contamination modeling. Gateway Configuration 1 consisted of the Power and Propulsion Element (PPE) and the Habitation and Logistics Outpost (HALO). It also included payloads like the European Radiation Sensor Array (ERSA) attached to PPE and the Heliophysics Environmental and Radiation Measurement Experiment Suite (HERMES) attached to HALO. The element- and payload-level analyses to be introduced in this paper for HALO, HERMES, and ERSA enabled high-fidelity descriptions of Gateway’s external molecular contamination environment. Approaches to geometric modeling, meshing, outgassing rate assignment, molecular transport modeling, and analysis methodology will be presented. Element and payload contaminant deposition onto sensitive Gateway receiver surfaces will be summarized and results compared to induced environment requirements. While these results incorporate refinements made over the course of the program, they were not intended to be final. Therefore, modeling assumptions and inputs, potential improvements, and lessons-learned will be documented to inform future work on Moon Base, repurposed elements and payloads, and other space stations.

Gateway↗

Large HAWT wake measurement and analysis

From the theoretical fluid dynamics point of view, the wake region of a large horizontal-axis wind turbine has been defined and described, and numerical models of wake behavior have been developed. Wind tunnel studies of single turbine wakes and turbine array wakes have been used to verify the theory and further refine the numerical models. However, the effects of scaling, rotor solidity, and topography on wake behavior are questions that remain unanswered. In the wind tunnel studies, turbines were represented by anything from scaled models to tea strainers or wire mesh disks whose solidity was equivalent to that of a typical wind turbine. The scale factor compensation for the difference in Reynolds number between the scale model and an actual turbine is complex, and not typically accounted for. Though it is wise to study the simpler case of wakes in flat topography, which can be easily duplicated in the wind tunnel, current indications are that wind turbine farm development is actually occurring in somewhat more complex terrain. Empirical wake studies using large horizontal-axis wind turbines have not been thoroughly composited, and, therefore, the results have not been applied to the well-developed theory of wake structure. The measurement programs have made use of both in situ sensor systems, such as instrumented towers, and remote sensors, such as kites and tethered, balloonborne anemometers. We present a concise overview of the work that has been performed, including our own, which is based on the philosophy that the MOD-2 turbines are probably their own best detector of both the momentum deficit and the induced turbulence effect downwind. Only the momentum deficit aspects of the wake/machine interactions have been addressed. Both turbine power output deficits and wind energy deficits as measured by the onsite meteorological towers have been analyzed from a composite data set. The analysis has also evidenced certain topographic influences on the operation of spatially diverse wind turbines.

Miller, A. H.↗

Adjoint-Based, Three-Dimensional Error Prediction and Grid Adaptation

Engineering computational fluid dynamics (CFD) analysis and design applications focus on output functions (e.g., lift, drag). Errors in these output functions are generally unknown and conservatively accurate solutions may be computed. Computable error estimates can offer the possibility to minimize computational work for a prescribed error tolerance. Such an estimate can be computed by solving the flow equations and the linear adjoint problem for the functional of interest. The computational mesh can be modified to minimize the uncertainty of a computed error estimate. This robust mesh-adaptation procedure automatically terminates when the simulation is within a user specified error tolerance. This procedure for estimating and adapting to error in a functional is demonstrated for three-dimensional Euler problems. An adaptive mesh procedure that links to a Computer Aided Design (CAD) surface representation is demonstrated for wing, wing-body, and extruded high lift airfoil configurations. The error estimation and adaptation procedure yielded corrected functions that are as accurate as functions calculated on uniformly refined grids with ten times as many grid points.

Park, Michael A.↗

Multidimensional Unstructured-Grid Liquid Rocket Engine Nozzle Performance and Heat Transfer Analysis

The objective of this study is to conduct a unified computational analysis for computing design parameters such as axial thrust, convective and radiative wall heat fluxes for regeneratively cooled liquid rocket engine nozzles, so as to develop a computational strategy for computing those parameters through parametric investigations. The computational methodology is based on a multidimensional, finite-volume, turbulent, chemically reacting, radiating, unstructured-grid, and pressure-based formulation, with grid refinement capabilities. Systematic parametric studies on effects of wall boundary conditions, combustion chemistry, radiation coupling, computational cell shape, and grid refinement were performed and assessed. Under the computational framework of this study, it is found that the computed axial thrust performance, flow features, and wall heat fluxes compared well with those of available data and calculations, using a strategy of structured-grid dominated mesh, finite-rate chemistry, and cooled wall boundary condition.

Wang, Ten-See↗

Preparation and Analysis of the Geometry Models used in the 1st AIAA Geometry and Mesh Generation Workshop

The NASA High-Lift Common Research Model (HL-CRM) was the subject model chosen for the AIAA Geometry and Mesh Generation Workshop I (GMGW-1) and High-Lift Prediction Workshop III (HLPW-3). This paper describes how geometry models of the HLCRM were prepared for use in the workshops and reviews the analysis of their construction that was provided to workshop participants. The refinements made to the HL-CRM geometry model immediately after GMGW-1 are also presented.

Taylor, Nigel J.↗