Search NASA⌕ Search

SEARCH · Search NASA

Results for “Code Robustness”

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 109 records · Page 6

Concept Demonstration of Dopant Selective Reactive Etching (DSRIE) in Silicon Carbide

Accurate quantification of combustor pressure dynamics for the primary purpose of experimental validation of computational fluid dynamics (CFD) codes requires the use of robust, reliable and sensitive pressure sensors that can resolve sub--pound-per-square-inch pressure levels in high temperature environments (i.e., combustor). The state of the art microfabricated piezoresistive silicon carbide (SiC) pressure sensors that we have developed are capable of operating reliably at 600 degrees Centigrade. This technology was used in support of the ARMD ISRP-ERA (NASA's Aeronautics Research Mission Directorate, Integrated System Research Project - Environmentally Responsible Aviation) program to quantify combustor thermoacoustic instabilities. The results showed that while the SiC pressure sensors survived the high temperature and measured instabilities, the diaphragm (force collector) was not thin enough to be sensitive in resolving sub-pound-per-square-inch pressures; 30 meters is the thinnest diaphragm achievable with conventional reactive ion etching (RIE) processes. Therefore, this precludes its use for sub-pound-per-square-inch pressure measurement with high fidelity. In order to effectively resolve sub-pound-per-square-inch pressures, a thinner more sensitive diaphragm (10 meters) is needed. To achieve this would require a new and innovative fabrication process technique.

static↗

Parallelized Carbon Oxidation and Surface Recession Model in Direct-Simulation Monte Carlo Code, SPARTA

Ablative thermal protection system (TPS) materials for spacecraft are composites that often consist of a carbon-based reinforcement and a polymeric matrix. During Earth re-entry, they endure high-temperature oxidation and surface recession. Oxidation is an important mechanism for ablation, sometimes leading to the weakening, spallation, or failure of the oxidized fibers at the surface and in the char layer. However, more details are required including accurate material properties of the fiber microstructure, whether the fibers recede homogenously or localized at pits, and the role of pyrolysis outgassing in order to evaluate the role oxidation plays in the degradation and failure mechanisms of these materials. In this work, we demonstrate a parallelizable oxidation-driven ablation model developed for detailed, large-scale simulations in the DSMC code SPARTA. We also develop robust algorithms which handle the conservation of the surface state after an ablation step. Finally, we verify our model for both simple and more complex chemistry as well as microstructures with reaction statistics, oxidation depth calculations, and detailed visualization of oxidation-induced material recession at the microscale.

V Arias↗

A computer code for multiphase all-speed transient flows in complex geometries. MAST version 1.0

The operation of the MAST code, which computes transient solutions to the multiphase flow equations applicable to all-speed flows, is described. Two-phase flows are formulated based on the Eulerian-Lagrange scheme in which the continuous phase is described by the Navier-Stokes equation (or Reynolds equations for turbulent flows). Dispersed phase is formulated by a Lagrangian tracking scheme. The numerical solution algorithms utilized for fluid flows is a newly developed pressure-implicit algorithm based on the operator-splitting technique in generalized nonorthogonal coordinates. This operator split allows separate operation on each of the variable fields to handle pressure-velocity coupling. The obtained pressure correction equation has the hyperbolic nature and is effective for Mach numbers ranging from the incompressible limit to supersonic flow regimes. The present code adopts a nonstaggered grid arrangement; thus, the velocity components and other dependent variables are collocated at the same grid. A sequence of benchmark-quality problems, including incompressible, subsonic, transonic, supersonic, gas-droplet two-phase flows, as well as spray-combustion problems, were performed to demonstrate the robustness and accuracy of the present code.

Chen, C. P.↗

A numerical code for a three-dimensional magnetospheric MHD equilibrium model

Two dimensional and three dimensional MHD equilibrium models were begun for Earth's magnetosphere. The original proposal was motivated by realizing that global, purely data based models of Earth's magnetosphere are inadequate for studying the underlying plasma physical principles according to which the magnetosphere evolves on the quasi-static convection time scale. Complex numerical grid generation schemes were established for a 3-D Poisson solver, and a robust Grad-Shafranov solver was coded for high beta MHD equilibria. Thus, the effects were calculated of both the magnetopause geometry and boundary conditions on the magnetotail current distribution.

Voigt, G.-H.↗

Lewis inverse design code (LINDES): Users manual

The method of complex characteristics and hodograph transformation for the design of shockless airfoils was introduced by Bauer, Garabedian, and Korn and has been extended by the author to design subcritical and supercritical cascades with high solidities and large inlet angles. This new capability was achieved by introducing a new conformal mapping of the hodograph domain onto an ellipse and expanding the solution in terms of Chebyshev polynomials. A new computer code, the NASA Lewis inverse design code, was developed based on this idea. This new design code is an efficient method for the design of airfoils in cascade. In particular, the design of subcritical cascades of airfoils is a very fast, robust, and versatile process. The inverse design code can be made to interact with a turbulent boundary layer calculation to obtain airfoils with no separated flows at the design condition. This report is intended to serve as a users manual for this design code. Material previously reported by the author is included here for completeness and quick access to the user. The manual contains a description of the method followed by a discussion of the design procedure and examples. The input parameters necessary to run the code are then described and their default values given. Output listings corresponding to six different blade shapes designed with the code are given, as well as the necessary input data to reproduce the computer runs. The examples have been chosen to show that a wide range of applications can be covered with the code, ranging from supercritical propeller sections to wind tunnel turning vanes that can operate with a large inlet flow angle range.

Sanz, Jose M.↗

Simulation of patch and slot antennas using FEM with prismatic elements and investigations of artificial absorber mesh termination schemes

Year 1 progress can be characterized with four major achievements which are crucial toward the development of robust, easy to use antenna analysis code on doubly conformal platforms. (1) A new FEM code was developed using prismatic meshes. This code is based on a new edge based distorted prism and is particularly attractive for growing meshes associated with printed slot and patch antennas on doubly conformal platforms. It is anticipated that this technology will lead to interactive, simple to use codes for a large class of antenna geometries. Moreover, the codes can be expanded to include modeling of the circuit characteristics. An attached report describes the theory and validation of the new prismatic code using reference calculations and measured data collected at the NASA Langley facilities. The agreement between the measured and calculated data is impressive even for the coated patch configuration. (2) A scheme was developed for improved feed modeling in the context of FEM. A new approach based on the voltage continuity condition was devised and successfully tested in modeling coax cables and aperture fed antennas. An important aspect of this new feed modeling approach is the ability to completely separate the feed and antenna mesh regions. In this manner, different elements can be used in each of the regions leading to substantially improved accuracy and meshing simplicity. (3) A most important development this year has been the introduction of the perfectly matched interface (PMI) layer for truncating finite element meshes. So far the robust boundary integral method has been used for truncating the finite element meshes. However, this approach is not suitable for antennas on nonplanar platforms. The PMI layer is a lossy anisotropic absorber with zero reflection at its interface. (4) We were able to interface our antenna code FEMA_CYL (for antennas on cylindrical platforms) with a standard high frequency code. This interface was achieved by first generating equivalent magnetic currents across the antenna aperture using the FEM code. These currents were employed as the sources in the high frequency code.

Gong, J.↗

A hybrid algorithm for transonic airfoil and wing design

The present method for the design of transonic airfoils and wings employs a predictor/corrector approach in which an analysis code calculates the flowfield for an initial geometry, then modifies it on the basis of the difference between calculated and target pressures. This allows the design method to be straightforwardly coupled with any existing analysis code, as presently undertaken with several two- and three-dimensional potential flow codes. The results obtained indicate that the method is robust and accurate, even in the cases of airfoils with strongly supercritical flow and shocks. The design codes are noted to require computational resources typical of current pure-inverse methods.

Campbell, Richard L.↗

Vectorized schemes for conical potential flow using the artificial density method

A method is developed to determine solutions to the full-potential equation for steady supersonic conical flow using the artificial density method. Various update schemes used generally for transonic potential solutions are investigated. The schemes are compared for speed and robustness. All versions of the computer code have been vectorized and are currently running on the CYBER-203 computer. The update schemes are vectorized, where possible, either fully (explicit schemes) or partially (implicit schemes). Since each version of the code differs only by the update scheme and elements other than the update scheme are completely vectorizable, comparisons of computational effort and convergence rate among schemes are a measure of the specific scheme's performance. Results are presented for circular and elliptical cones at angle of attack for subcritical and supercritical crossflows.

Bradley, P. F.↗

Finite element simulation of temperature dependent free surface flows

The method of Engelman and Sani (1984) for a finite-element simulation of incompressible surface flows with a free and/or moving fluid interface, such as encountered in crystal growth and coating and polymer technology, is extended to temperature-dependent flows, including the effect of temperature-dependent surface tension. The basic algorithm of Saito and Scriven (1981) and Ruschak (1980) has been generalized and implemented in a robust and versatile finite-element code that can be employed with relative ease for the simulation of free-surface problems in complex geometries. As a result, the costly dependence on the Newton-Raphson algorithm has been eliminated by replacing it with a quasi-Newton iterative method, which nearly retains the superior convergence properties of the Newton-Raphson method.

Engelman, M. S.↗

Structural Modeling of a Five-Meter Thin Film Inflatable Antenna/Concentrator With Rigidized Support Struts

Dynamic characterization of a non-rigidized thin film inflatable antenna/solar concentrator structure with rigidized composite support struts is described in detail. A two-step finite element modeling approach in MSC/NASTRAN is utilized, consisting of: (1) a nonlinear static pressurization procedure used to obtain the updated stiffness matrix, and (2) a modal "restart" eigen solution that uses the modified stiffness matrix. Unique problems encountered in modeling of this large 5-m lightweight inflatable are identified, including considerable difficulty in obtaining convergence in the nonlinear pressurization solution. It was found that the extremely thin polyimide film material (.001 in or I mil) presents tremendous problems in obtaining a converged solution when internal pressure loading is applied. It was concluded that the ratios of film thickness to other geometric dimensions such as torus cross-sectional and ring diameter and lenticular diameter are the critical parameters for convergence of the pressurization procedure. Comparison of finite element predictions for frequency and mode shapes with experimental results indicated reasonable agreement considering the complexity of the structure, the film-to-air interaction, and the nonlinear material properties of the film. It was also concluded that analysis should be done using different finite element to codes to determine if a more robust and stable solution can be obtained.

Smalley, Kurt B.↗

Carbon Nanostructure Examined by Lattice Fringe Analysis of High Resolution Transmission Electron Microscopy Images

The dimensions of graphitic layer planes directly affect the reactivity of soot towards oxidation and growth. Quantification of graphitic structure could be used to develop and test correlations between the soot nanostructure and its reactivity. Based upon transmission electron microscopy images, this paper provides a demonstration of the robustness of a fringe image analysis code for determining the level of graphitic structure within nanoscale carbon, i.e. soot. Results, in the form of histograms of graphitic layer plane lengths, are compared to their determination through Raman analysis.

VanderWal, Randy L.↗

Carbon Nanostructure Examined by Lattice Fringe Analysis of High Resolution Transmission Electron Microscopy Images

The dimensions of graphitic layer planes directly affect the reactivity of soot towards oxidation and growth. Quantification of graphitic structure could be used to develop and test correlations between the soot nanostructure and its reactivity. Based upon transmission electron microscopy images, this paper provides a demonstration of the robustness of a fringe image analysis code for determining the level of graphitic structure within nanoscale carbon, i.e., soot. Results, in the form of histograms of graphitic layer plane lengths, are compared to their determination through Raman analysis.

VanderWal, Randy L.↗

Jet Noise Modeling for Suppressed and Unsuppressed Aircraft in Simulated Flight

This document describes the development of further extensions and improvements to the jet noise model developed by Modern Technologies Corporation (MTC) for the National Aeronautics and Space Administration (NASA). The noise component extraction and correlation approach, first used successfully by MTC in developing a noise prediction model for two-dimensional mixer ejector (2DME) nozzles under the High Speed Research (HSR) Program, has been applied to dual-stream nozzles, then extended and improved in earlier tasks under this contract. Under Task 6, the coannular jet noise model was formulated and calibrated with limited scale model data, mainly at high bypass ratio, including a limited-range prediction of the effects of mixing-enhancement nozzle-exit chevrons on jet noise. Under Task 9 this model was extended to a wider range of conditions, particularly those appropriate for a Supersonic Business Jet, with an improvement in simulated flight effects modeling and generalization of the suppressor model. In the present task further comparisons are made over a still wider range of conditions from more test facilities. The model is also further generalized to cover single-stream nozzles of otherwise similar configuration. So the evolution of this prediction/analysis/correlation approach has been in a sense backward, from the complex to the simple; but from this approach a very robust capability is emerging. Also from these studies, some observations emerge relative to theoretical considerations. The purpose of this task is to develop an analytical, semi-empirical jet noise prediction method applicable to takeoff, sideline and approach noise of subsonic and supersonic cruise aircraft over a wide size range. The product of this task is an even more consistent and robust model for the Footprint/Radius (FOOTPR) code than even the Task 9 model. The model is validated for a wider range of cases and statistically quantified for the various reference facilities. The possible role of facility effects will thus be documented. Although the comparisons that can be accomplished within the limited resources of this task are not comprehensive, they provide a broad enough sampling to enable NASA to make an informed decision on how much further effort should be expended on such comparisons. The improved finalized model is incorporated into the FOOTPR code. MTC has also supported the adaptation of this code for incorporation in NASA s Aircraft Noise Prediction Program (ANOPP).

Stone, James R.↗

TRANAIR - A computer code for transonic analyses of arbitrary configurations

Attention is given to a new approach to solving full potential equations about arbitrary configurations. Numerical algorithms from such fields as finite elements, preconditioned Krylov subspace methods, discrete Fourier analysis, and integral equations are combined to take advantage of the size and speed of current and emerging supercomputers. On the basis of this appraoch, a robust, efficient and easy to use computer code referred to as TRANAIR has been developed for transonic analysis of complex geometries.

Samant, S. S.↗

Tool for Generation of MAC/GMC Representative Unit Cell for CMC/PMC Analysis

This document describes a recently developed analysis tool that enhances the resident capabilities of the Micromechanics Analysis Code with the Generalized Method of Cells (MAC/GMC) 4.0. This tool is especially useful in analyzing ceramic matrix composites (CMCs), where higher fidelity with improved accuracy of local response is needed. The tool, however, can be used for analyzing polymer matrix composites (PMCs) as well. MAC/GMC 4.0 is a composite material and laminate analysis software developed at NASA Glenn Research Center. The software package has been built around the concept of the generalized method of cells (GMC). The computer code is developed with a user friendly framework, along with a library of local inelastic, damage, and failure models. Further, application of simulated thermomechanical loading, generation of output results, and selection of architectures to represent the composite material have been automated to increase the user friendliness, as well as to make it more robust in terms of input preparation and code execution. Finally, classical lamination theory has been implemented within the software, wherein GMC is used to model the composite material response of each ply. Thus, the full range of GMC composite material capabilities is available for analysis of arbitrary laminate configurations as well. The primary focus of the current effort is to provide a graphical user interface (GUI) capability that generates a number of different user-defined repeating unit cells (RUCs). In addition, the code has provisions for generation of a MAC/GMC-compatible input text file that can be merged with any MAC/GMC input file tailored to analyze composite materials. Although the primary intention was to address the three different constituents and phases that are usually present in CMCs-namely, fibers, matrix, and interphase-it can be easily modified to address two-phase polymer matrix composite (PMC) materials where an interphase is absent. Currently, the tool capability includes generation of RUCs for square packing, hexagonal packing, and random fiber packing as well as RUCs based on actual composite micrographs. All these options have the fibers modeled as having a circular cross-sectional area. In addition, a simplified version of RUC is provided where the fibers are treated as having a square cross section and are distributed randomly. This RUC facilitates a speedy analysis using the higher fidelity version of GMC known as HFGMC. The first four mentioned options above support uniform subcell discretization. The last one has variable subcell sizes due to the primary intention of keeping the RUC size to a minimum to gain the speed ups using the higher fidelity version of MAC. The code is implemented within the MATLAB (The Mathworks, Inc., Natick, MA) developmental framework; however, a standalone application that does not need a priori MATLAB installation is also created with the aid of the MATLAB compiler.

Materials Engineering↗

A Compilation of MATLAB Scripts and Functions for MACGMC Analyses

The primary aim of the current effort is to provide scripts that automate many of the repetitive pre- and post-processing tasks associated with composite materials analyses using the Micromechanics Analysis Code with the Generalized Method of Cells. This document consists of a compilation of hundreds of scripts that were developed in MATLAB (The Mathworks, Inc., Natick, MA) programming language and consolidated into 16 MATLAB functions. (MACGMC). MACGMC is a composite material and laminate analysis software code developed at NASA Glenn Research Center. The software package has been built around the generalized method of cells (GMC) family of micromechanics theories. The computer code is developed with a user-friendly framework, along with a library of local inelastic, damage, and failure models. Further, application of simulated thermo-mechanical loading, generation of output results, and selection of architectures to represent the composite material have been automated to increase the user friendliness, as well as to make it more robust in terms of input preparation and code execution. Finally, classical lamination theory has been implemented within the software, wherein GMC is used to model the composite material response of each ply. Thus, the full range of GMC composite material capabilities is available for analysis of arbitrary laminate configurations as well. The pre-processing tasks include generation of a multitude of different repeating unit cells (RUCs) for CMCs and PMCs, visualization of RUCs from MACGMC input and output files and generation of the RUC section of a MACGMC input file. The post-processing tasks include visualization of the predicted composite response, such as local stress and strain contours, damage initiation and progression, stress-strain behavior, and fatigue response. In addition to the above, several miscellaneous scripts have been developed that can be used to perform repeated Monte-Carlo simulations to enable probabilistic simulations with minimal manual intervention. This document is formatted to provide MATLAB source files and descriptions of how to utilize them. It is assumed that the user has a basic understanding of how MATLAB scripts work and some MATLAB programming experience.

Materials Engineering↗

Modification and Validation of Conceptual Design Aerodynamic Prediction Method HASC95 With VTXCHN

A conceptual/preliminary design level subsonic aerodynamic prediction code HASC (High Angle of Attack Stability and Control) has been improved in several areas, validated, and documented. The improved code includes improved methodologies for increased accuracy and robustness, and simplified input/output files. An engineering method called VTXCHN (Vortex Chine) for predicting nose vortex shedding from circular and non-circular forebodies with sharp chine edges has been improved and integrated into the HASC code. This report contains a summary of modifications, description of the code, user's guide, and validation of HASC. Appendices include discussion of a new HASC utility code, listings of sample input and output files, and a discussion of the application of HASC to buffet analysis.

Aerodynamics prediction↗

Flow of GE90 Turbofan Engine Simulated

The objective of this task was to create and validate a three-dimensional model of the GE90 turbofan engine (General Electric) using the APNASA (average passage) flow code. This was a joint effort between GE Aircraft Engines and the NASA Lewis Research Center. The goal was to perform an aerodynamic analysis of the engine primary flow path, in under 24 hours of CPU time, on a parallel distributed workstation system. Enhancements were made to the APNASA Navier-Stokes code to make it faster and more robust and to allow for the analysis of more arbitrary geometry. The resulting simulation exploited the use of parallel computations by using two levels of parallelism, with extremely high efficiency.The primary flow path of the GE90 turbofan consists of a nacelle and inlet, 49 blade rows of turbomachinery, and an exhaust nozzle. Secondary flows entering and exiting the primary flow path-such as bleed, purge, and cooling flows-were modeled macroscopically as source terms to accurately simulate the engine. The information on these source terms came from detailed descriptions of the cooling flow and from thermodynamic cycle system simulations. These provided boundary condition data to the three-dimensional analysis. A simplified combustor was used to feed boundary conditions to the turbomachinery. Flow simulations of the fan, high-pressure compressor, and high- and low-pressure turbines were completed with the APNASA code.

Veres, Joseph P.↗