Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation codes”

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 127 records · Page 7

Verification and Validation of Elastodynamic Simulation Software for Aerospace Research

Physics-based simulation of nondestructive evaluation (NDE) inspection can help to advance the inspectability and reliability of mechanical systems. However, NDE simulations applicable to non-idealized mechanical components often require large compute domains and long run times. This has prompted development of custom NDE simulation software tailored to high performance computing (HPC) hardware. Verification and validation (V&V) is an integral part of developing this software to ensure implementations are robust and applicable to inspection problems, producing tools and simulations suitable for computational NDE research. This presentation addresses factors common to V&V of several elastodynamic simulation codes applicable to ultrasonic NDE. Examples are drawn from in-house simulation software at NASA Langley Research Center, ranging from ensuring reliability in a 1D heterogeneous media wave equation solver to the V&V needs of 3D cluster-parallel elastodynamic software. Factors specific to a research environment are addressed, where individual simulation results can be as relevant as the software product itself. Distinct facets of V&V are discussed including testing to establish software reliability, employing systematic approaches for consistency with fundamental conservation laws, establishing the numerical stability of algorithms, and demonstrating concurrence with empirical data. This talk also addresses V&V practices for small groups of researchers. This includes establishing resources (e.g. time and personnel) for V&V during project planning to mitigate and control the risk of setbacks. Similarly, we identify ways for individual researchers to use V&V during simulation software development itself to both speed up the development process and reduce incurred technical debt.

NDE↗

Highlights of Project Accomplishments (2000-2004): Redshift Dependence of the Interaction-Activity Connection

We completed this project by developing the new galaxy-galaxy collision simulation code that we originally proposed. We included star formation heuristically, along with gas recycling and energetic feedback into the interstellar medium of the galaxies. We ran several test simulations. And we finally validated the code by running models to emulate the double-ring galaxy AM 0644-741. Most of the exotic and unique features of this collisional ring galaxy were matched by our models, thereby validating our algorithm, which may now be used by us and by the community for additional dynamic studies of galaxies. A paper has been submitted to the Astrophysical Journal describing our algorithm and the results of the AM 0644-741 modeling.

Borne, Kirk D.↗

Numerical simulation of turbulent jet noise, part 2

Results on the numerical simulation of jet flow fields were used to study the radiated sound field, and in addition, to extend and test the capabilities of the turbulent jet simulation codes. The principal result of the investigation was the computation of the radiated sound field from a turbulent jet. In addition, the computer codes were extended to account for the effects of compressibility and eddy viscosity, and the treatment of the nonlinear terms of the Navier-Stokes equations was modified so that they can be computed in a semi-implicit way. A summary of the flow model and a description of the numerical methods used for its solution are presented. Calculations of the radiated sound field are reported. In addition, the extensions that were made to the fundamental dynamical codes are described. Finally, the current state-of-the-art for computer simulation of turbulent jet noise is summarized.

Metcalfe, R. W.↗

Development of UMARC (University of Maryland Advanced Rotorcraft Code)

The University of Maryland Advanced Rotorcraft Code (UMARC) is a user-friendly, FEM-based comprehensive helicopter rotor simulation code of high numerical robustness and computational efficiency. UMARC formulates the rotor-fuselage equations using Hamilton's principle, and are discretized using finite elements in space and time. The FEM formulation allows the code to analyze a wide variety of rotor designs. Dynamic inflow modeling is used for unsteady wake inflow computations. Predicted stability, response, and blade-load data are validated with experimental data for several configurations, including representative articulated, hingeless, and bearingless rotors.

Bir, Gunjit↗

Shock Position Control for Mode Transition in a Turbine Based Combined Cycle Engine Inlet Model

A dual flow-path inlet for a turbine based combined cycle (TBCC) propulsion system is to be tested in order to evaluate methodologies for performing a controlled inlet mode transition. Prior to experimental testing, simulation models are used to test, debug, and validate potential control algorithms which are designed to maintain shock position during inlet disturbances. One simulation package being used for testing is the High Mach Transient Engine Cycle Code simulation, known as HiTECC. This paper discusses the development of a mode transition schedule for the HiTECC simulation that is analogous to the development of inlet performance maps. Inlet performance maps, derived through experimental means, describe the performance and operability of the inlet as the splitter closes, switching power production from the turbine engine to the Dual Mode Scram Jet. With knowledge of the operability and performance tradeoffs, a closed loop system can be designed to optimize the performance of the inlet. This paper demonstrates the design of the closed loop control system and benefit with the implementation of a Proportional-Integral controller, an H-Infinity based controller, and a disturbance observer based controller; all of which avoid inlet unstart during a mode transition with a simulated disturbance that would lead to inlet unstart without closed loop control.

Csank, Jeffrey T.↗

Advanced Shock Position Control for Mode Transition in a Turbine Based Combined Cycle Engine Inlet Model

A dual flow-path inlet system is being tested to evaluate methodologies for a Turbine Based Combined Cycle (TBCC) propulsion system to perform a controlled inlet mode transition. Prior to experimental testing, simulation models are used to test, debug, and validate potential control algorithms. One simulation package being used for testing is the High Mach Transient Engine Cycle Code simulation, known as HiTECC. This paper discusses the closed loop control system, which utilizes a shock location sensor to improve inlet performance and operability. Even though the shock location feedback has a coarse resolution, the feedback allows for a reduction in steady state error and, in some cases, better performance than with previous proposed pressure ratio based methods. This paper demonstrates the design and benefit with the implementation of a proportional-integral controller, an H-Infinity based controller, and a disturbance observer based controller.

Csank, Jeffrey T.↗

Validation of NASA Thermal Ice Protection Computer Codes: The Validation of Antice - Part 3

An experimental program was generated by the Icing Technology Branch at NASA Glenn Research Center to validate two ice protection simulation codes: (1) LEWICE/Thermal for transient electrothermal de-icing and anti-icing simulations, and (2) ANTICE for steady state hot gas and electrothermal anti-icing simulations. An electrothermal ice protection system was designed and constructed integral to a 36 inch chord NACA0012 airfoil. The model was fully instrumented with thermo-couples, RTD'S, and heat flux gages. Tests were conducted at several icing environmental conditions during a two week period at the NASA Glenn Icing Research Tunnel. Experimental results of running-wet and evaporative cases were compared to the ANTICE computer code predictions and are presented in this paper.

Al-Khalil, Kamel M.↗

Implementation and Characterization of Three-Dimensional Particle-in-Cell Codes on Multiple-Instruction-Multiple-Data Massively Parallel Supercomputers

A three-dimensional electrostatic particle-in-cell (PIC) plasma simulation code has been developed on coarse-grain distributed-memory massively parallel computers with message passing communications. Our implementation is the generalization to three-dimensions of the general concurrent particle-in-cell (GCPIC) algorithm. In the GCPIC algorithm, the particle computation is divided among the processors using a domain decomposition of the simulation domain. In a three-dimensional simulation, the domain can be partitioned into one-, two-, or three-dimensional subdomains ("slabs," "rods," or "cubes") and we investigate the efficiency of the parallel implementation of the push for all three choices. The present implementation runs on the Intel Touchstone Delta machine at Caltech; a multiple-instruction-multiple-data (MIMD) parallel computer with 512 nodes. We find that the parallel efficiency of the push is very high, with the ratio of communication to computation time in the range 0.3%-10.0%. The highest efficiency (> 99%) occurs for a large, scaled problem with 64(sup 3) particles per processing node (approximately 134 million particles of 512 nodes) which has a push time of about 250 ns per particle per time step. We have also developed expressions for the timing of the code which are a function of both code parameters (number of grid points, particles, etc.) and machine-dependent parameters (effective FLOP rate, and the effective interprocessor bandwidths for the communication of particles and grid points). These expressions can be used to estimate the performance of scaled problems--including those with inhomogeneous plasmas--to other parallel machines once the machine-dependent parameters are known.

Lyster, P. M.↗

MOD-0 wind turbine dynamics test correlations

The behavior of the teetered, downwind, free yaw, MOD-0 wind turbine, as represented by NASA dynamic test data, was used to support confidence in the Hamilton Standard computer code simulations. Trim position, performance at trim, and teeter response as predicted by the computer codes were compared to test results. Using the computer codes, other possible configurations for MOD-0 were investigated. Several new test configurations are recommended for exploring free yaw behavior. It is shown that eliminating rotor tilt and optimizing cooling and blade twist can contribute to good free yaw behavior and stability. The effects of rotor teeter, teeter gravity balance, inflow and other physical and operating parameters were also investigated.

Brooks, B. M.↗

Self-consistent simulation of plasma interactions with secondary-emitting insulators

A cylindrical particle-in-cell (PIC) plasma simulation code applicable to plasma densities encountered in low Earth orbit (LEO) is described. The simulated geometries include that of a plain disk and a disk surrounded by a dielectric. Both configurations are mounted upon a ground plate in contact with a plasma environment. Techniques allowing simulation of dielectric charging using PIC time scales are discussed. Current versus voltage characteristic curves are calculated and the results are compared to experimental data.

Brandon, S. T.↗

The role of tip clearance in high-speed fan stall

A numerical experiment has been carried out to define the near-stall casing endwall flowfield of a high-speed fan rotor. The experiment used a simulation code incorporating a simple clearance model, whose calibration is presented. The results of the simulation show that the interaction of the tip leakage vortex and the in-pasage shock plays a major role in determining the fan flow range. More specifically, the computations imply that it is the area increase of this vortex as it passes through the in-passage shock, which is the source of the blockage associated with stall. In addition, for fans of this type, it is the clearance over the forward portion of the fan blade which controls the flow processes leading to stall.

Adamczyk, J. J.↗

Thermal-Flow Code for Modeling Gas Dynamics and Heat Transfer in Space Shuttle Solid Rocket Motor Joints

A new, thermal-flow simulation code, called SFLOW. has been developed to model the gas dynamics, heat transfer, as well as O-ring and flow path erosion inside the space shuttle solid rocket motor joints by combining SINDA/Glo, a commercial thermal analyzer. and SHARPO, a general-purpose CFD code developed at Thiokol Propulsion. SHARP was modified so that friction, heat transfer, mass addition, as well as minor losses in one-dimensional flow can be taken into account. The pressure, temperature and velocity of the combustion gas in the leak paths are calculated in SHARP by solving the time-dependent Navier-Stokes equations while the heat conduction in the solid is modeled by SINDA/G. The two codes are coupled by the heat flux at the solid-gas interface. A few test cases are presented and the results from SFLOW agree very well with the exact solutions or experimental data. These cases include Fanno flow where friction is important, Rayleigh flow where heat transfer between gas and solid is important, flow with mass addition due to the erosion of the solid wall, a transient volume venting process, as well as some transient one-dimensional flows with analytical solutions. In addition, SFLOW is applied to model the RSRM nozzle joint 4 subscale hot-flow tests and the predicted pressures, temperatures (both gas and solid), and O-ring erosions agree well with the experimental data. It was also found that the heat transfer between gas and solid has a major effect on the pressures and temperatures of the fill bottles in the RSRM nozzle joint 4 configuration No. 8 test.

Wang, Qunzhen↗

Three-Dimensional Simulation of Traveling-Wave Tube Cold-Test Characteristics Using CST MICROWAVE STUDIO

The electromagnetic field simulation software package CST MICROWAVE STUDIO (MWS) was used to compute the cold-test parameters - frequency-phase dispersion, on-axis impedance, and attenuation - for a traveling-wave tube (TWT) slow-wave circuit. The results were compared to experimental data, as well as to results from MAFIA, another three-dimensional simulation code from CST currently used at the NASA Glenn Research Center (GRC). The strong agreement between cold-test parameters simulated with MWS and those measured experimentally demonstrates the potential of this code to reduce the time and cost of TWT development.

Chevalier, Christine T.↗

Computer simulation of auroral kilometric radiation

An investigation is carried out of the linear amplification and nonlinear saturation of electromagnetic waves associated with two components of electrons consisting of a cold Maxwellian background and a weakly relativistic population of electrons possessing a loss-cone distribution. The goal is to understand the mechanism responsible for the production of auroral kilometric radiation (AKR). A relativistic 1-2/2 dimensional electromagnetic particle simulation code is used in studying an initial value problem modeling the AKR source region. The simulation makes it possible to follow the growth of all modes of radiation, including the ordinary mode (O-mode), slow extraordinary mode (Z-mode), whistler mode and fast extraordinary mode (X-mode), past the point of saturation.

Wagner, J. S.↗

Dynamics of Unusual Debris Flows on Martian Sand Dunes

Gullies that dissect sand dunes in Russell impact crater often display debris flow-like deposits in their distal reaches. The possible range of both the rheological properties and the flow rates are estimated using a numerical simulation code of a Bingham plastic flow to help explain the formation of these features. Our simulated results are best explained by a rapid debris flow. For example, a debris flow with the viscosity of 10(exp 2) Pa s and the yield strength of 10(exp 2) Pa can form the observed deposits with a flow rate of 0.5 cu m/s sustained over several minutes and total discharged water volume on the order of hundreds of cubic meters, which may be produced by melting a surface layer of interstitial ice within the dune deposits to several centimeters depth.

Miyamoto, Hideaki↗

Imposing a Lagrangian Particle Framework on an Eulerian Hydrodynamics Infrastructure in Flash

In many astrophysical simulations, both Eulerian and Lagrangian quantities are of interest. For example, in a galaxy cluster merger simulation, the intracluster gas can have Eulerian discretization, while dark matter can be modeled using particles. FLASH, a component-based scientific simulation code, superimposes a Lagrangian framework atop an adaptive mesh refinement Eulerian framework to enable such simulations. The discretization of the field variables is Eulerian, while the Lagrangian entities occur in many different forms including tracer particles, massive particles, charged particles in particle-in-cell mode, and Lagrangian markers to model fluid structure interactions. These widely varying roles for Lagrangian entities are possible because of the highly modular, flexible, and extensible architecture of the Lagrangian framework. In this paper, we describe the Lagrangian framework in FLASH in the context of two very different applications, Type Ia supernovae and galaxy cluster mergers, which use the Lagrangian entities in fundamentally different ways.

Dubey, A.↗

High Resolution Aerospace Applications using the NASA Columbia Supercomputer

This paper focuses on the parallel performance of two high-performance aerodynamic simulation packages on the newly installed NASA Columbia supercomputer. These packages include both a high-fidelity, unstructured, Reynolds-averaged Navier-Stokes solver, and a fully-automated inviscid flow package for cut-cell Cartesian grids. The complementary combination of these two simulation codes enables high-fidelity characterization of aerospace vehicle design performance over the entire flight envelope through extensive parametric analysis and detailed simulation of critical regions of the flight envelope. Both packages. are industrial-level codes designed for complex geometry and incorpor.ats. CuStomized multigrid solution algorithms. The performance of these codes on Columbia is examined using both MPI and OpenMP and using both the NUMAlink and InfiniBand interconnect fabrics. Numerical results demonstrate good scalability on up to 2016 CPUs using the NUMAIink4 interconnect, with measured computational rates in the vicinity of 3 TFLOP/s, while InfiniBand showed some performance degradation at high CPU counts, particularly with multigrid. Nonetheless, the results are encouraging enough to indicate that larger test cases using combined MPI/OpenMP communication should scale well on even more processors.

Mavriplis, Dimitri J.↗

An Iterative Soft-Decision Decoding Algorithm

This paper presents a new minimum-weight trellis-based soft-decision iterative decoding algorithm for binary linear block codes. Simulation results for the RM(64,22), EBCH(64,24), RM(64,42) and EBCH(64,45) codes show that the proposed decoding algorithm achieves practically (or near) optimal error performance with significant reduction in decoding computational complexity. The average number of search iterations is also small even for low signal-to-noise ratio.

Lin, Shu↗