Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mesh data structure”

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 199 records · Page 11

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

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

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

An investigation of turbulence models

The accuracy to which a turbulent boundary layer or wake can be predicted numerically depends on the validity of the turbulence closure model used. The modeling of turbulence physics is one of the most difficult problems in computational fluid dynamics (CFD). In fact, it is one of the pacing factors in the development of CFD. In general, there are three main approaches to the description of trubulence physics. First is turbulence modeling in which the Reynolds averaged Navier-Stokes equations are used and some closure approximation is made for the Reynolds stresses. A second approach to turbulence is large eddy simulation (LES) in which the computational mesh is taken to be fine enough that the large scale structure of the turbulence can be calculated directly. An empirical assumption must be made for the small scale sub-grid turbulence. The third approach is direct simulation. In this technique the Navier-Strokes equations are solved directly on a mesh which if fine enough to resolve the smallest length scale of the turbulence. The Reynolds averaged equations are not used and no closure assumption is required. These last two approaches require extensive computer resources and as such are not engineering tools. The purpose of the work was to investigate the various engineering turbulence models for accuracy and ease of programming. This involved comparison of the models with each other and with experimental data.

Doria, Michael L.↗

An Examination of Parameters Affecting Large Eddy Simulations of Flow Past a Square Cylinder

Separated flow over a bluff body is analyzed via large eddy simulations. The turbulent flow around a square cylinder features a variety of complex flow phenomena such as highly unsteady vortical structures, reverse flow in the near wall region, and wake turbulence. The formation of spanwise vortices is often times artificially suppressed in computations by either insufficient depth or a coarse spanwise resolution. As the resolution is refined and the domain extended, the artificial turbulent energy exchange between spanwise and streamwise turbulence is eliminated within the wake region. A parametric study is performed highlighting the effects of spanwise vortices where the spanwise computational domain's resolution and depth are varied. For Re=22,000, the mean and turbulent statistics computed from the numerical large eddy simulations (NLES) are in good agreement with experimental data. Von-Karman shedding is observed in the wake of the cylinder. Mesh independence is illustrated by comparing a mesh resolution of 2 million to 16 million. Sensitivities to time stepping were minimized and sampling frequency sensitivities were nonpresent. While increasing the spanwise depth and resolution can be costly, this practice was found to be necessary to eliminating the artificial turbulent energy exchange.

Mankbadi, M. R.↗

Structure of the Turbulent Separated Flow Around a Stalled Airfoil

Hot-wire measurements were made in the boundary layer, the separated region, and the near wake for flow past a NACA 4412 airfoil at maximum lift. The Reynolds number based on chord was 1,500,000. Special care was taken to achieve a two-dimensional mean flow. Data were obtained at several thousand locations in the flow field. These data include intermittency, two components of mean velocity, and mean values for three double, four triple, and five quadruple products of two velocity fluctuations. No information was obtained about the third (spanwise) velocity component. Smoothing and interpolating routines were used to determine intermittency, two components of mean velocity, and mean values of three double, four triple, and five quadruple products of two velocity fluctuations on a fine rectangular mesh aligned with the airfoil chord. The data are presented in contour plots, in three-dimensional plots, and in tabular form. The format used to store the experimental data in digital form is described and a computer program which illustrates how this data can be accessed is presented.

Wadcock, A. J.↗

Hypersonic Flow Computations on Unstructured Meshes

A method for computing inviscid hypersonic flow over complex configurations using unstructured meshes is presented. The unstructured grid solver uses an edge{based finite{volume formulation. Fluxes are computed using a flux vector splitting scheme that is capable of representing constant enthalpy solutions. Second{order accuracy in smooth flow regions is obtained by linearly reconstructing the solution, and stability near discontinuities is maintained by locally forcing the scheme to reduce to first-order accuracy. The implementation of the algorithm to parallel computers is described. Computations using the proposed method are presented for a sphere-cone configuration at Mach numbers of 5.25 and 10.6, and a complex hypersonic re-entry vehicle at Mach numbers of 4.5 and 9.8. Results are compared to experimental data and computations made with established structured grid methods. The use of the solver as a screening tool for rapid aerodynamic assessment of proposed vehicles is described.

Bibb, K. L.↗

Comparison of High-Order and Low-Order Methods for Large-Eddy Simulation of a Compressible Shear Layer

The objective of this work is to compare a high-order solver with a low-order solver for performing Large-Eddy Simulations (LES) of a compressible mixing layer. The high-order method is the Wave-Resolving LES (WRLES) solver employing a Dispersion Relation Preserving (DRP) scheme. The low-order solver is the Wind-US code, which employs the second-order Roe Physical scheme. Both solvers are used to perform LES of the turbulent mixing between two supersonic streams at a convective Mach number of 0.46. The high-order and low-order methods are evaluated at two different levels of grid resolution. For a fine grid resolution, the low-order method produces a very similar solution to the highorder method. At this fine resolution the effects of numerical scheme, subgrid scale modeling, and filtering were found to be negligible. Both methods predict turbulent stresses that are in reasonable agreement with experimental data. However, when the grid resolution is coarsened, the difference between the two solvers becomes apparent. The low-order method deviates from experimental results when the resolution is no longer adequate. The high-order DRP solution shows minimal grid dependence. The effects of subgrid scale modeling and spatial filtering were found to be negligible at both resolutions. For the high-order solver on the fine mesh, a parametric study of the spanwise width was conducted to determine its effect on solution accuracy. An insufficient spanwise width was found to impose an artificial spanwise mode and limit the resolved spanwise modes. We estimate that the spanwise depth needs to be 2.5 times larger than the largest coherent structures to capture the largest spanwise mode and accurately predict turbulent mixing.

large eddy simulations↗

Comparison of High-Order and Low-Order Methods for Large-Eddy Simulation of a Compressible Shear Layer

The objective of this work is to compare a high-order solver with a low-order solver for performing large-eddy simulations (LES) of a compressible mixing layer. The high-order method is the Wave-Resolving LES (WRLES) solver employing a Dispersion Relation Preserving (DRP) scheme. The low-order solver is the Wind-US code, which employs the second-order Roe Physical scheme. Both solvers are used to perform LES of the turbulent mixing between two supersonic streams at a convective Mach number of 0.46. The high-order and low-order methods are evaluated at two different levels of grid resolution. For a fine grid resolution, the low-order method produces a very similar solution to the high-order method. At this fine resolution the effects of numerical scheme, subgrid scale modeling, and filtering were found to be negligible. Both methods predict turbulent stresses that are in reasonable agreement with experimental data. However, when the grid resolution is coarsened, the difference between the two solvers becomes apparent. The low-order method deviates from experimental results when the resolution is no longer adequate. The high-order DRP solution shows minimal grid dependence. The effects of subgrid scale modeling and spatial filtering were found to be negligible at both resolutions. For the high-order solver on the fine mesh, a parametric study of the spanwise width was conducted to determine its effect on solution accuracy. An insufficient spanwise width was found to impose an artificial spanwise mode and limit the resolved spanwise modes. We estimate that the spanwise depth needs to be 2.5 times larger than the largest coherent structures to capture the largest spanwise mode and accurately predict turbulent mixing.

large eddy simulations↗

A point implicit unstructured grid solver for the Euler and Navier-Stokes equations

An upwind finite element technique that uses cell centered quantities and implicit and/or explicit time marching has been developed for computing hypersonic laminar viscous flows using adaptive unstructured triangular grids. A structured grid of quadrilaterals is laid out near the body surface. For inviscid flows the method is stable at Courant numbers of over 100,000. A first order basic scheme and a higher order flux corrected transport (FCT) scheme have been implemented. This technique has been applied to the problem of predicting type III and IV shock wave interactions on a cylinder, with a view of simulating the pressure and heating rate augmentation caused by an impinging shock on the leading edge of a cowl lip of an engine inlet. The predictions of wall pressure and heating rates compare very well with experimental data. The flow features are very distinctly captured with a sequence of adaptively generated grids. The adaptive mesh generator and the upwind Navier-Stokes solver are combined in a set of programs called LARCNESS, an acronym for Langley Adaptive Remeshing Code and Navier-Stokes Solver.

Thareja, Rajiv R.↗

Composite Fuselage Technology

Work was conducted over a ten-year period to address the central issue of damage in primary load-bearing aircraft composite structure, specifically fuselage structure. This included the three facets of damage resistance, damage tolerance, and damage arrest. Experimental, analytical, and numerical work was conducted in order to identify and better understand the mechanisms that control the structural behavior of fuselage structures in their response to the three aspects of damage. Furthermore, work was done to develop straightforward design methodologies that can be employed by structural designers in preliminary design stages to make intelligent choices concerning the material, layup, and structural configurations so that a more efficient structure with structural integrity can be designed and built. Considerable progress was made towards achieving these goals via this work. In regard to damage tolerance considerations, the following were identified as important effects: composite layup and associated orthotropy/structural anisotropy, specifics of initial local damage mechanisms, role of longitudinal versus hoop stress, and large deformation and associated geometric nonlinearity. Means were established to account for effects of radius and for the nonlinear response. In particular, nondimensional parameters were identified to characterize the importance of nonlinearity in the response of pressurized cylinders. This led to the establishment of a iso-nonlinear-error plot for reference in structural design. Finally, in the case of damage tolerance, the general approach of the original methodology to predict the failure pressure involving extending basic plate failure data by accounting for the local stress intensification was accomplished for the general case by accounting for the mechanisms noted by utilizing the capability of the STAGS finite element code and numerically calculating the local stress intensification for the particular configuration to be considered. For the issue of damage arrest, placement of and configuration of stiffeners (including stiffener curvature), and magnitude and orientation of principal strains due to local bending were found to be key considerations. Means were established to account for stiffener effectiveness quantitatively based on radius, slit size, stiffener curvature' and relative bending stifffiesses involved. Geometric nonlinearity was also found to play an - 24 - important role here. Furthermore, it was determined that damage propagation is controlled by different mechanisms (hoop stress versus flapping stress and the associated factors involved in each) depending upon the direction of damage propagation. This latter item results in an inability to scale these phenomena in one test due to the different factors involved. Finally, the importance of shell curvature and associated instability in response to transverse loading including impact were found to be important considerations in damage resistance. A technique, involving asymmetric meshing of a finite element mesh, was developed to predict this behavior and showed excellent correlation with experimental results. Further details of these ten years of work are presented herein with references made to the fourteen documents produced during this work where full details can be found. Implications of this work are discussed and recommendations made. Although it is clear that there is more work to be done to fully understand composite fuselage technology and specifically the overall issue of damage in primary load-bearing composite structures, important understanding and capability has been extended via this work.

Lagace, Paul A.↗

A comparison of cosmological hydrodynamic codes

We present a detailed comparison of the simulation results of various hydrodynamic codes. Starting with identical initial conditions based on the cold dark matter scenario for the growth of structure, with parameters h = 0.5 Omega = Omega(sub b) = 1, and sigma(sub 8) = 1, we integrate from redshift z = 20 to z = O to determine the physical state within a representative volume of size L(exp 3) where L = 64 h(exp -1) Mpc. Five indenpendent codes are compared: three of them Eulerian mesh-based and two variants of the smooth particle hydrodynamics 'SPH' Lagrangian approach. The Eulerian codes were run at N(exp 3) = (32(exp 3), 64(exp 3), 128(exp 3), and 256(exp 3)) cells, the SPH codes at N(exp 3) = 32(exp 3) and 64(exp 3) particles. Results were then rebinned to a 16(exp 3) grid with the exception that the rebinned data should converge, by all techniques, to a common and correct result as N approaches infinity. We find that global averages of various physical quantities do, as expected, tend to converge in the rebinned model, but that uncertainites in even primitive quantities such as (T), (rho(exp 2))(exp 1/2) persists at the 3%-17% level achieve comparable and satisfactory accuracy for comparable computer time in their treatment of the high-density, high-temeprature regions as measured in the rebinned data; the variance among the five codes (at highest resolution) for the mean temperature (as weighted by rho(exp 2) is only 4.5%. Examined at high resolution we suspect that the density resolution is better in the SPH codes and the thermal accuracy in low-density regions better in the Eulerian codes. In the low-density, low-temperature regions the SPH codes have poor accuracy due to statiscal effects, and the Jameson code gives the temperatures which are too high, due to overuse of artificial viscosity in these high Mach number regions. Overall the comparison allows us to better estimate errors; it points to ways of improving this current generation ofhydrodynamic codes and of suiting their use to problems which exploit their best individual features.

Kang, Hyesung↗