Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “code comparisons”

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 253 records · Page 14

Application of CFD to sonic boom near and mid flow-field prediction

A 3-D parabolized Navier-Stokes (PNS) code was used to calculate the supersonic overpressures from three different geometries at near- and mid-flow fields. Wind tunnel data is used for code validation. Comparison of the computed results with different grid refinements is shown. It is observed that a large number of grid points is needed to resolve the tail shock/expansion fan interaction. Therefore, an adaptive grid approach is employed to calculate the flow field. The agreement between the numerical results and the wind tunnel data confirms that computational fluid dynamics can be applied to the problem of sonic boom prediction.

Cheung, Samson H.↗

High-order ENO schemes applied to two- and three-dimensional compressible flow

High order essentially non-oscillatory (ENO) finite difference schemes are applied to the 2-D and 3-D compressible Euler and Navier-Stokes equations. Practical issues, such as vectorization, efficiency of coding, cost comparison with other numerical methods, and accuracy degeneracy effects, are discussed. Numerical examples are provided which are representative of computational problems of current interest in transition and turbulence physics. These require both nonoscillatory shock capturing and high resolution for detailed structures in the smooth regions and demonstrate the advantage of ENO schemes.

Shu, Chi-Wang↗

Application of CFD to sonic boom near and mid flow-field prediction

A three-dimensional parabolized Navier-Stokes (PNS) code has been used to calculate the supersonic overpressures from three different geometries at near- and mid-flow fields. Wind-tunnel data is used for code validation. Comparison of the computed results with different grid refinements is shown in this paper. It is observed that a large number of grid points is needed to resolve the tail shock/expansion fan interaction. Therefore, an adaptive grid approach is employed to calculate the flow field. The agreement between the numerical results and the wind-tunnel data confirms that computational fluid dynamics can be applied to the problem of sonic boom prediction.

Cheung, Samson H.↗

Global synchronization algorithms for the Intel iPSC/860

In a distributed memory multicomputer that has no global clock, global processor synchronization can only be achieved through software. Global synchronization algorithms are used in tridiagonal systems solvers, CFD codes, sequence comparison algorithms, and sorting algorithms. They are also useful for event simulation, debugging, and for solving mutual exclusion problems. For the Intel iPSC/860 in particular, global synchronization can be used to ensure the most effective use of the communication network for operations such as the shift, where each processor in a one-dimensional array or ring concurrently sends a message to its right (or left) neighbor. Three global synchronization algorithms are considered for the iPSC/860: the gysnc() primitive provided by Intel, the PICL primitive sync0(), and a new recursive doubling synchronization (RDS) algorithm. The performance of these algorithms is compared to the performance predicted by communication models of both the long and forced message protocols. Measurements of the cost of shift operations preceded by global synchronization show that the RDS algorithm always synchronizes the nodes more precisely and costs only slightly more than the other two algorithms.

Seidel, Steven R.↗

Seals Flow Code Development 1993

Seals Workshop of 1993 code releases include SPIRALI for spiral grooved cylindrical and face seal configurations; IFACE for face seals with pockets, steps, tapers, turbulence, and cavitation; GFACE for gas face seals with 'lift pad' configurations; and SCISEAL, a CFD code for research and design of seals of cylindrical configuration. GUI (graphical user interface) and code usage was discussed with hands on usage of the codes, discussions, comparisons, and industry feedback. Other highlights for the Seals Workshop-93 include environmental and customer driven seal requirements; 'what's coming'; and brush seal developments including flow visualization, numerical analysis, bench testing, T-700 engine testing, tribological pairing and ceramic configurations, and cryogenic and hot gas facility brush seal results. Also discussed are seals for hypersonic engines and dynamic results for spiral groove and smooth annular seals.

Liang, Anita D.↗

A Generalized Fluid System Simulation Program to Model Flow Distribution in Fluid Networks

This paper describes a general purpose computer program for analyzing steady state and transient flow in a complex network. The program is capable of modeling phase changes, compressibility, mixture thermodynamics and external body forces such as gravity and centrifugal. The program's preprocessor allows the user to interactively develop a fluid network simulation consisting of nodes and branches. Mass, energy and specie conservation equations are solved at the nodes; the momentum conservation equations are solved in the branches. The program contains subroutines for computing "real fluid" thermodynamic and thermophysical properties for 33 fluids. The fluids are: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutane, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride and ammonia. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. Seventeen different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include: pipe flow, flow through a restriction, non-circular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct, labyrinth seal, parallel plates, common fittings and valves, pump characteristics, pump power, valve with a given loss coefficient, and a Joule-Thompson device. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. This paper also illustrates the application and verification of the code by comparison with Hardy Cross method for steady state flow and analytical solution for unsteady flow.

Majumdar, Alok↗

Mach Number Effects on Turbine Blade Transition Length Prediction

The effect of a Mach number correction on a model for predicting the length of transition was investigated. The transition length decreases as the turbulent spot production rate increases. Much of the data for predicting the spot production rate comes from low speed flow experiments. Recent data and analysis showed that the spot production rate is affected by Mach number. The degree of agreement between analysis and data for turbine blade heat transfer without film cooling is strongly dependent of accurately predicting the length of transition. Consequently, turbine blade heat transfer data sets were used to validate a transition length turbulence model. A method for modifying models for the length of transition to account for Mach number effects is presented. The modification was made to two transition length models. The modified models were incorporated into the two-dimensional Navier-Stokes code, RVCQ3D. Comparisons were made between predicted and measured midspan surface heat transfer for stator and rotor turbine blades. The results showed that accounting for Mach number effects significantly improved the agreement with the experimental data.

Boyle, R. J.↗

Follow-on Low Noise Fan Aerodynamic Study

The focus of the project was to investigate the effects of turbulence models on the prediction of rotor wake structures. The Advanced Ducted Propfan Analysis (ADPAC) code was modified through the incorporation of the Spalart-Allmaras one-equation turbulence model. Suitable test cases were solved numerically using ADPAC employing the Spalart-Allmaras turbulence model and another prediction code for comparison. A near-wall spacing study was also completed to determine the adequate spacing of the first computational cell off the wall. Solutions were also collected using two versions of the algebraic Baldwin-Lomax turbulence model in ADPAC. The effects of the turbulence model on the rotor wake definition was examined by obtaining ADPAC solutions for the Low Noise Fan rotor-only steady-flow case using the standard algebraic Baldwin-Lomax turbulence model, a modified version of the Baldwin-Lomax turbulence model and the one-equation Spalart-Allmaras turbulence model. The results from the three different turbulence modeling techniques were compared with each other and the available experimental data. These results include overall rotor performance, spanwise exit profiles, and contours of axial velocity taken along constant axial locations and along blade-to-blade surfaces. Wake characterizations were also performed on the experimental and ADPAC predicted results including the definition of a wake correlation function. Correlations were evaluated for wake width and wake depth. Similarity profiles of the wake shape were also compared between all numerical solutions and experimental data.

Heidegger, Nathan J.↗

CFD Predictions of Sonic-Boom Characteristics for Unmodified and Modified SR-71 Configurations

Shaped sonic-boom signatures refer to signatures that look something other than the typical N-waves. Shaped sonic-boom signatures such as "flat-top," "ramp-type," or "hybrid-type" waveforms have been shown to reduce the subjective loudness without requiring reductions in overpressure peaks. The shaping of sonic-boom signatures requires increasing the shock rise time and changes in frequency spectra. So far, a flat-top waveform was shown to be achievable in wind tunnels; however, the influence of long propagation distance and real atmosphere on shaped signatures should be addressed using flight tests. Two different approaches have been proposed for sonic-boom minimization flight tests. The first approach, proposed by Eagle Aerospace, is for a flight test using a modified BQM-34 "FIREBEE" remotely piloted vehicle. The 30-foot long FIREBEE has a steady state flight condition at the Mach number and altitude of interest, and it can be recovered by helicopter from the water. As an alternative approach, a modified SR-71 vehicle has been proposed by the McDonnell Douglas Corporation. Benefits of the SR-71 include its variable geometry supersonic inlets, small cockpit bulge, higher Mach number capabilities, slender design, and longer length (105 foot). The present investigation addresses the sonic-boom analysis for the second vehicle.The objective of the current investigation is to assess the feasibility of a modified SR-71 configuration, with McDonnell Douglas-designed fuselage modifications, intended to produce shaped sonic-boom signatures on the ground. The present study describes the use of a higher-order computational fluid dynamics (CFD) method to predict the sonic-boom characteristics for both unmodified and modified SR-71 configurations. An Euler unstructured grid methodology is used to predict the near-field, three-dimensional pressure patterns generated by both SR-71 models. The computed near-field pressure signatures are extrapolated to specified distances below the aircraft down to impingement on the ground using the code MDBOOM. Comparisons of the near-field pressure signatures with available flight-test data are presented in the current paper.

Fouladi, Kamran↗

Comparison of the NASCAP/GEO, POLAR, SEE Charging Handbook, and NASCAP-2K.1 Spacecraft Charging Codes

The NASA Charging Analyzer Program (NASCAP) spacecraft charging software developed by Maxwell Technologies has been widely used for the past fifteen to twenty years in satellite design and investigation of spacecraft charging related anomalies. Individual versions of the NASCAP software are available for use in low inclination, low Earth orbit environments (NASCAP[LEO) and geostationary orbit environments (NASCAP/GEO). In addition, the Potentials of Large objects in the Auroral Region (POLAR) code is available for use in LEO polar orbit environments. NASCAP/GEO and POLAR were both written in the 1980's using algorithms appropriate for the computers of the time. They solve the Poisson-Vlasov system for currents and densities assuming limited speed and memory of computer systems standard for the day. In addition, use of the charging models requires individual input files that are not readily transported into the various codes to facilitate comparison of results by the user.

Neergaard, Linda E.↗

Preliminary Structural Sizing and Alternative Material Trade Study of CEV Crew Module

This paper presents the results of a preliminary structural sizing and alternate material trade study for NASA s Crew Exploration Vehicle (CEV) Crew Module (CM). This critical CEV component will house the astronauts during ascent, docking with the International Space Station, reentry, and landing. The alternate material design study considers three materials beyond the standard metallic (aluminum alloy) design that resulted from an earlier NASA Smart Buyer Team analysis. These materials are graphite/epoxy composite laminates, discontinuously reinforced SiC/Al (DRA) composites, and a novel integrated panel material/concept known as WebCore. Using the HyperSizer (Collier Research and Development Corporation) structural sizing software and NASTRAN finite element analysis code, a comparison is made among these materials for the three composite CM concepts considered by the 2006 NASA Engineering and Safety Center Composite Crew Module project.

Bednarcyk, Brett A.↗

Assessment of Radiated Fan Noise Prediction Capabilities Using Static Engine Test Data

This paper describes further assessment of the CDUCT-LaRC code via comparison with static engine test data. In an effort to improve confidence in the use of CDUCT-LaRC for liner optimization studies addressing realistic three-dimensional geometries, inlet radiated fan noise predictions were performed at 54% and 87% engine speed settings. Predictions were then compared with far-field measurements to assess the approach and implementation. The particular configurations were chosen to exercise the three-dimensional capability of CDUCT-LaRC and it s applicability to realistic configurations and conditions. At the 54% engine speed setting, the predictions capture the general directivity and acoustic treatment effects quite well. Comparisons of the predicted and measured directivity at the 87% power setting were more problematic. This was likely due in part to the difficulties in source specification and possibly the nonlinear nature of buzz-saw tones at this engine operating condition. Overall, the approach captured the basic trends and provided a conservative estimate of liner effects from which relative performance metrics could be inferred.

Nark, Douglas M.↗

Atomic Data and Spectral Line Intensities for NI XVII

Electron impact collision strengths, energy levels, oscillator strengths, and spontaneous radiative decay rates are calculated for Ni XVII. We include in the calculations the 23 lowest configurations, corresponding to 159 fine-structure levels: 3l3l', 3l4l0'' , and 3s5l0''' , with l,l' = s,p,d, l'' = s,p,d, f, and l''' = s,p,d. Collision strengths are calculated at five incident energies for all transitions at varying energies above the threshold of each transition. One additional energy, very close to the threshold of each transition, has also been included. Calculations have been carried out using the Flexible Atomic Code in the distorted wave approximation. Additional calculations have been performed with the University College London suite of codes for comparison. Excitation rate coefficients are calculated as a function of electron temperature by assuming a Maxwellian electron velocity distribution. Using the excitation rate coefficients and the radiative transition rates of the present work, statistical equilibrium equations for level populations are solved at electron densities covering the range of 10(exp 8) - 10(exp 14) / cubic cm and at an electron temperature of logT(sub e)e(K) = 6.5, corresponding to the maximum abundance of Ni XVII. Spectral line intensities are calculated, and their diagnostic relevance is discussed. This dataset will be made available in the next version of the CHIANTI database

Bhatia, A. K.↗

Development and Evaluation of an Order-N Formulation for Multi-Flexible Body Space Systems

This paper presents development of a generic recursive Order-N algorithm for systems with rigid and flexible bodies, in tree or closed-loop topology, with N being the number of bodies of the system. Simulation results are presented for several test cases to verify and evaluate the performance of the code compared to an existing efficient dense mass matrix-based code. The comparison brought out situations where Order-N or mass matrix-based algorithms could be useful.

Multibody dynamics↗

Numerical and Experimental Study of Transport Phenomena in Directional Solidification of Succinonitrile

A numerical and experimental study of the growth of succinonitrile (SCN) using a horizontal Bridginan furnace and transparent glass ampoule was conducted. Two experiments were considered: one in which the temperature profile was fixed relative to the ampoule (no-growth case); and a second in which the thermal profile was translated at a constant rate (steady growth case). Measured temperature profiles on the outer surface of the ampoule were used as thermal boundary conditions for the modelling. The apparent heat capacity formulation combined with the variable viscositymeth~ was used to model the phase change in SeN. Both 2-D and 3-D models were studied and numerical solutions obtained using the commercial finite element code, FIDAP1. Comparison of the numerical results to experimental data showed excellent agreement. The complex 3-D shallow-cavity flow in the melt, differences between 2-D and 3-D models, effects of natural convection on the thermal gradient and shape of the solid/liquid interface, and the sensitivity of simulations to specific assumptions, are also discussed.

de Groh, Henry C., III↗

Output-Based Adaptive Meshing Applied to Space Launch System Booster Separation Analysis

This paper presents details of Computational Fluid Dynamic (CFD) simulations of the Space Launch System during solid-rocket booster separation using the Cart3D inviscid code with comparisons to Overflow viscous CFD results and a wind tunnel test performed at NASA Langley Research Center's Unitary PlanWind Tunnel. The Space Launch System (SLS) launch vehicle includes two solid-rocket boosters that burn out before the primary core stage and thus must be discarded during the ascent trajectory. The main challenges for creating an aerodynamic database for this separation event are the large number of basis variables (including orientation of the core, relative position and orientation of the boosters, and rocket thrust levels) and the complex flow caused by the booster separation motors. The solid-rocket boosters are modified from their form when used with the Space Shuttle Launch Vehicle, which has a rich flight history. However, the differences between the SLS core and the Space Shuttle External Tank result in the boosters separating with much narrower clearances, and so reducing aerodynamic uncertainty is necessary to clear the integrated system for flight. This paper discusses an approach that has been developed to analyze about 6000 wind tunnel simulations and 5000 flight vehicle simulations using Cart3D in adaptive-meshing mode. In addition, a discussion is presented of Overflow viscous CFD runs used for uncertainty quantification. Finally, the article presents lessons learned and improvements that will be implemented in future separation databases.

Booster↗

Modeling of Multiscale Solar Dynamics for Understanding Drivers of Space Weather

Understanding the solar dynamics is critical for improving our capabilities to forecast the evolution of space weather conditions. We take advantage of currently available computational capabilities to model solar dynamics from the deep interior to the corona and investigate mechanisms that may drive space weather conditions. The simulations are performed using the 3D radiative MHD code StellarBox. Comparison of synthetic spectroscopic observables obtained from numerical simulations and actual observations allows us to uncover physical processes associated with observed phenomena. To facilitate a transition from modeling short-term physical phenomena to developing a reliable forecast-oriented model, we suggest using the data assimilation approach. It allows us to cross-analyze dynamo model solutions and observations and to consider possible uncertainties and errors. In this presentation, we briefly summarize current multi-scale modeling capabilities and results and discuss ongoing developments to build a reliable physics-based forecast-oriented model of solar activity.

Irina N Kitiashvili↗

Modeling of Multiscale Solar Dynamics for Understanding Drivers of Space Weather

Understanding the solar dynamics is critical for improving our capabilities to forecast the evolution of space weather conditions. We take advantage of currently available computational capabilities to model solar dynamics from the deep interior to the corona and investigate mechanisms that may drive space weather conditions. The simulations are performed using the 3D radiative MHD code StellarBox. Comparison of synthetic spectroscopic observables obtained from numerical simulations and actual observations allows us to uncover physical processes associated with observed phenomena. To facilitate a transition from modeling short-term physical phenomena to developing a reliable forecast-oriented model, we suggest using the data assimilation approach. It allows us to cross-analyze dynamo model solutions and observations and to consider possible uncertainties and errors. In this presentation, we briefly summarize current multi-scale modeling capabilities and results and discuss ongoing developments to build a reliable physics-based forecast-oriented model of solar activity

Irina N Kitiashvili↗