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 1,477 records · Page 82

Enhancing Application Performance Using Mini-Apps: Comparison of Hybrid Parallel Programming Paradigms

In many fields, real-world applications for High Performance Computing have already been developed. For these applications to stay up-to-date, new parallel strategies must be explored to yield the best performance; however, restructuring or modifying a real-world application may be daunting depending on the size of the code. In this case, a mini-app may be employed to quickly explore such options without modifying the entire code. In this work, several mini-apps have been created to enhance a real-world application performance, namely the VULCAN code for complex flow analysis developed at the NASA Langley Research Center. These mini-apps explore hybrid parallel programming paradigms with Message Passing Interface (MPI) for distributed memory access and either Shared MPI (SMPI) or OpenMP for shared memory accesses. Performance testing shows that MPI+SMPI yields the best execution performance, while requiring the largest number of code changes. A maximum speedup of 23 was measured for MPI+SMPI, but only 11 was measured for MPI+OpenMP.

Lawson, Gary↗

Computation of three-dimensional, rotational flow through turbomachinery blade rows for improved aerodynamic design studies

The performance of a three dimensional computer code developed for predicting the flowfield in stationary and rotating turbomachinery blade rows is described in this study. The four stage Runge-Kutta numerical integration scheme is used for solving the governing flow equations and yields solution to the full, three dimensional, unsteady Euler equations in cylindrical coordinates. This method is fully explicit and uses the finite volume, time marching procedure. In order to demonstrate the accuracy and efficiency of the code, steady solutions were obtained for several cascade geometries under widely varying flow conditions. Computed flowfield results are presented for a fully subsonic turbine stator and a low aspect ratio, transonic compressor rotor blade under maximum flow and peak efficiency design conditions. Comparisons with Laser Anemometer measurements and other numerical predictions are also provided to illustrate that the present method predicts important flow features with good accuracy and can be used for cost effective aerodynamic design studies.

Subramanian, S. V.↗

Numerical study of hydrogen-air supersonic combustion by using elliptic and parabolized equations

The two-dimensional Navier-Stokes and species continuity equations are used to investigate supersonic chemically reacting flow problems which are related to scramjet-engine configurations. A global two-step finite-rate chemistry model is employed to represent the hydrogen-air combustion in the flow. An algebraic turbulent model is adopted for turbulent flow calculations. The explicit unsplit MacCormack finite-difference algorithm is used to develop a computer program suitable for a vector processing computer. The computer program developed is then used to integrate the system of the governing equations in time until convergence is attained. The chemistry source terms in the species continuity equations are evaluated implicitly to alleviate stiffness associated with fast chemical reactions. The problems solved by the elliptic code are re-investigated by using a set of two-dimensional parabolized Navier-Stokes and species equations. A linearized fully-coupled fully-implicit finite difference algorithm is used to develop a second computer code which solves the governing equations by marching in spce rather than time, resulting in a considerable saving in computer resources. Results obtained by using the parabolized formulation are compared with the results obtained by using the fully-elliptic equations. The comparisons indicate fairly good agreement of the results of the two formulations.

Chitsomboon, T.↗

A case study for the real-time experimental evaluation of the VIPER microprocessor

An experiment to evaluate the applicability of the Verifiable Integrated Processor for Enhanced Reliability (VIPER) microprocessor to real time control is described. The VIPER microprocessor was invented by the Royal Signals and Radar Establishment (RSRE), U.K., and is an example of the use of formal mathematical methods for developing electronic digital systems with a high degree of assurance on the system design and implementation correctness. The experiment consisted of selecting a control law, writing the control law algorithm for the VIPER processor, and providing real time, dynamic inputs into the processor and monitoring the outputs. The control law selected and coded for the VIPER processor was the yaw damper function of an automatic landing program for a 737 aircraft. The mechanisms for interfacing the VIPER Single Board Computer to the VAX host are described. Results include run time experiences, performance evaluation, and comparison of VIPER and FORTRAN yaw damper algorithm output for accuracy estimation.

Carreno, Victor A.↗

Finite element implementation of state variable-based viscoplasticity models

The implementation of state variable-based viscoplasticity models is made in a general purpose finite element code for structural applications of metals deformed at elevated temperatures. Two constitutive models, Walker's and Robinson's models, are studied in conjunction with two implicit integration methods: the trapezoidal rule with Newton-Raphson iterations and an asymptotic integration algorithm. A comparison is made between the two integration methods, and the latter method appears to be computationally more appealing in terms of numerical accuracy and CPU time. However, in order to make the asymptotic algorithm robust, it is necessary to include a self adaptive scheme with subincremental step control and error checking of the Jacobian matrix at the integration points. Three examples are given to illustrate the numerical aspects of the integration methods tested.

Iskovitz, I.↗

Comparison of turbulence models for powered-lift flow fields

The development of a practical turbulence model for the complex flows with strong streamline curvature associated with powered-lift aircraft was investigated. The approach was to modify existing models that were compatible with the widely used kappa-epsilon model already available in an existing 3D Navier-Stokes code. The state-of-the-art of current turbulence models was examined. Several models that provide curvature corrections to the kappa-epsilon model were numerically investigated. An algebraic Reynolds stress model, which provides additional physics, was also examined. Test cases including an impinging jet, an impinging jet in a crossflow, and two adjacent impinging jets with the resultant upwash fountain were calculated. The first two test cases could be modeled reasonsbly well with some of the models. The upwash fountain could not be modeled with sufficient accuracy, which is consistent with the findings of other investigators.

Roberts, Donald W.↗

Finite element implementation of state variable-based viscoplasticity models

The implementation of state variable-based viscoplasticity models is made in a general purpose finite element code for structural applications of metals deformed at elevated temperatures. Two constitutive models, Walker's and Robinson's models, are studied in conjunction with two implicit integration methods: the trapezoidal rule with Newton-Raphson iterations and an asymptotic integration algorithm. A comparison is made between the two integration methods, and the latter method appears to be computationally more appealing in terms of numerical accuracy and CPU time. However, in order to make the asymptotic algorithm robust, it is necessary to include a self adaptive scheme with subincremental step control and error checking of the Jacobian matrix at the integration points. Three examples are given to illustrate the numerical aspects of the integration methods tested.

Chang, T. Y. P.↗

Summary of completed project

The primary goal of the work in this proposal has been the determination of the requirement for acceleration of particles, especially electrons, during solar flares. This is one of the important and yet unsolved problems in solar flares and other astronomical objects. This is done by detailed theoretical investigation of interaction of accelerated particles with ambient plasma including interaction with particles, waves, and large scale magnetic fields. Comparison with observed impulsive phase radiation such as X-rays, gamma-rays, microwaves, and type III radiation is used for constraining the models and for determination of the characteristics of the accelerated particles. Steady state and time-dependent numerical codes have been developed based on the Fokker-Planck method for study of the particle acceleration and transport during a flare. Using these, the correlation between characteristics of radiation at various wavelengths have been investigated and some constraints have been set on the model parameters. During the last three year period of this grant, dissertation research of three graduate students have been partially supported by this grant. Each of these theses has dealt with various observed properties of flare emission with particular emphasis on the interrelationship between radiation in at least two energy bands: X-rays and gamma-rays, X-rays and microwaves, and X-rays and type III radiation. This has allowed considerable tightening of constraints on models, and inclusion of wave particle interaction has allowed direct calculations of acceleration processes.

Petrosian, Vahe↗

Mars atmospheric dust properties: A synthesis of Mariner 9, Viking, and Phobos observations

We have modified a doubling-and-adding code to reanalyze the Mariner 9 IRIS spectra of Mars atmospheric dust as well as Viking IRTM EPF sequences in the 7, 9, and 20 micron channels. The code is capable of accurate emission/ absorption/scattering radiative transfer calculations over the 5-30 micron wavelength region for variable dust composition and particle size inputs, and incorporates both the Viking IRTM channel weightings and the Mariner 9 IRIS wavelength resolution for direct comparisons to these datasets. We adopt atmospheric temperature profiles according to the algorithm of Martin (1986) in the case of the Viking IRTM comparisons, and obtained Mariner 9 IRIS temperature retrievals from the 15 micron CO2 band for the case of the IRIS comparisons. We consider palagonite as the primary alternative to the montmorillonite composition of Mars atmospheric dust, based on several considerations. Palagonite absorbs in the ultraviolet and visible wavelength region due to its Fe content. Palagonite is also, in principal, consistent with the observed lack of clays on the Mars surface. Furthermore, palagonite does not display strong, structured absorption near 20 microns as does montmorillonite (in conflict with the IRIS observations). We propose that a palagonite composition with particle sizes roughly one-half that of the Toon et al. (1977) determination provide a much improved model to Mars atmospheric dust. Since palagonite is a common weathering product of terrrestrial basalts, it would not be unreasonable for palagonite to be a major surface component for Mars. The lack of even a minor component of Al-rich clays on the surface of Mars could be consistent with a palagonite composition for Mars dust if the conditions for basalt weathering on Mars were sufficiently anhydrous. Variations in palagonite composition could also lead to the inability of the modeled palagonite to fit the details of the 9 micron absorbtion indicated by the IRIS observations.

Clancy, R. T.↗

Prediction of Film Cooling on Gas Turbine Airfoils

A three-dimensional Navier-Stokes analysis tool has been developed in order to study the effect of film cooling on the flow and heat transfer characteristics of actual turbine airfoils. An existing code (Arnone et al., 1991) has been modified for the purpose. The code is an explicit, multigrid, cell-centered, finite volume code with an algebraic turbulence model. Eigenvalue scaled artificial dissipation and variable-coefficient implicit residual smoothing are used with a full-multigrid technique. Moreover, Mayle's transition criterion (Mayle, 1991) is used. The effects of film cooling have been incorporated into the code in the form of appropriate boundary conditions at the hole locations on the airfoil surface. Each hole exit is represented by several control volumes, thus providing an ability to study the effect of hole shape on the film-cooling characteristics. Comparison is fair with near mid-span experimental data for four and nine rows of cooling holes, five on the shower head, and two rows each on the pressure and suction surfaces. The computations, however, show a strong spanwise variation of the heat transfer coefficient on the airfoil surface, specially with shower-head cooling.

Garg, Vijay K.↗

Dynamic response of a pulsed Burke-Schumann diffusion flame

Turbulent flames are often envisioned as an ensemble of random vortices interacting with the combustion process. A better understanding of the vortex-flame interactions therefore would be useful in improving the modeling of turbulent diffusion flames. Substantial simplification may be made by investigating controlled interactions in a laminar flame, as opposed to random interactions in a turbulent flame. The general goals of the research project are to improve our understanding of (1) the influence of buoyancy on co-flow diffusion flames and (2) the effects of buoyancy on vortex-flame interactions in co-flow diffusion flames. As a first step toward objective (2), we conducted a joint experimental and numerical investigation of the vortex-flame interaction. Vortices were produced by mechanically pulsing the fuel flow at a low frequency, e.g., 10 Hz. Experiments were conducted using a nonflickering Burke-Schumann flame in both microgravity (mu-g) and normal gravity (1g) as a means of varying the buoyant force without modification of the pressure (i.e., density). The effects of buoyant convection may then be determined by a comparison of the mu-g and 1g results. The mu-g results may also reveal the important mechanisms which are masked or overwhelmed by buoyant convection in 1g. A numerical investigation was conducted using a validated, time-accurate numerical code to study the underlying physics during the flame interaction and to assist the interpretation of the experimental results.

Sheu, Jyh-Cherng↗

Hypersonic Laminar Viscous Flow Past Spinning Cones at Angle of Attack

Computational results are presented for hypersonic viscous flow past spinning sharp and blunt cones of angle of attack, obtained with a parabolic Navier-Stokes marching code. The code takes into account the asymmetries in the flowfield resulting from spinning motion and computes the asymmetric shock shape, cross-flow and streamwise shear, heat transfer, cross-flow separation, and vortex structure. The Magnus force and moments are also computed. Comparisons are made with other theoretical analyses based on boundary-layer and boundary-region equations, and an anomaly is discovered in the displacement thickness contribution to the Magnus force when compared with boundary-layer results.

Agarwal, Ramesh↗

Unsteady Analysis of Inlet-Compressor Acoustic Interactions Using Coupled 3-D and 1-D CFD Codes

It is well known that the dynamic response of a mixed compression supersonic inlet is very sensitive to the boundary condition imposed at the subsonic exit (engine face) of the inlet. In previous work, a 3-D computational fluid dynamics (CFD) inlet code (NPARC) was coupled at the engine face to a 3-D turbomachinery code (ADPAC) simulating an isolated rotor and the coupled simulation used to study the unsteady response of the inlet. The main problem with this approach is that the high fidelity turbomachinery simulation becomes prohibitively expensive as more stages are included in the simulation. In this paper, an alternative approach is explored, wherein the inlet code is coupled to a lesser fidelity 1-D transient compressor code (DYNTECC) which simulates the whole compressor. The specific application chosen for this evaluation is the collapsing bump experiment performed at the University of Cincinnati, wherein reflections of a large-amplitude acoustic pulse from a compressor were measured. The metrics for comparison are the pulse strength (time integral of the pulse amplitude) and wave form (shape). When the compressor is modeled by stage characteristics the computed strength is about ten percent greater than that for the experiment, but the wave shapes are in poor agreement. An alternate approach that uses a fixed rise in duct total pressure and temperature (so-called 'lossy' duct) to simulate a compressor gives good pulse shapes but the strength is about 30 percent low.

Suresh, A.↗

Analysis of Measured and Predicted Acoustics From an XV-15 Flight Test

Flight acoustic and vehicle state data from an XV-15 acoustic flight test are examined. Flight predictions using TiltRotor Aeroacoustic Code (TRAC) are performed for a level flight (repeated) and four descent conditions (including a blade vortex interaction (BVI)). The assumptions and procedures used for TRAC flight predictions as well as the variability in flight measurements, which are used for input and comparison to predictions, are investigated in detail. Differences were found in the measured vehicle airspeed, altitude, glideslope, and vehicle orientation (yaw, pitch, and roll angle) between each of the repeat runs. These differences violate some of the prediction assumptions and significantly impacted the resulting acoustic predictions. Multiple acoustic pulses, with a variable time between the pulses, were found in the measured acoustic time histories for the repeat runs. These differences could be attributed in part to the variability in vehicle orientation. Acoustic predictions that used the measured vehicle orientation for the repeat runs captured this multiple pulse variability. Thickness noise was found to be dominant on approach for all the cases, except the BVI condition. After the aircraft passed overhead, broadband noise and low frequency loading noise were dominant. The predicted LowSPL time histories compared well with measurement on approach to the array for the non-BVI conditions and poorly for the BVI condition. Accurate prediction of the lift share between the rotor and fuselage must be known in order to improve predictions. At a minimum, measurements of the rotor thrust and tip-path-plane angle are critical to further develop accurate flight acoustic prediction capabilities.

Boyd, D. Douglas, Jr.↗

CFD Results for an Axisymmetric Isentropic Relaxed Compression Inlet

The OVERFLOW code was used to calculate the flow field for a family of five relaxed compression inlets, which were part of a screening study to determine a configuration most suited to the application of microscale flow control technology as a replacement for bleed. Comparisons are made to experimental data collected for each of the inlets in the 1- by 1-Foot Supersonic Wind Tunnel at the NASA Glenn Research Center (GRC) to help determine the suitability of computational fluid dynamics (CFD) as a tool for future studies of these inlets with flow control devices. Effects on the wind tunnel results of the struts present in a high subsonic flow region accounted for most of the inconsistency between the results. Based on the level of agreement in the present study, it is expected that CFD can be used as a tool to aid in the design of a study of this class of inlets with flow control.

Hirt, Stefanie M.↗

Aerothermodynamics of Blunt Body Entry Vehicles

In this chapter, the aerothermodynamic phenomena of blunt body entry vehicles are discussed. Four topics will be considered that present challenges to current computational modeling techniques for blunt body environments: turbulent flow, non-equilibrium flow, rarefied flow, and radiation transport. Examples of comparisons between computational tools to ground and flight-test data will be presented in order to illustrate the challenges existing in the numerical modeling of each of these phenomena and to provide test cases for evaluation of Computational Fluid Dynamics (CFD) code predictions.

Hollis, Brian R.↗

The Omics of Stem Cell Mediated Regeneration: A Pilot Single Cell RNA-Seq Study of Mechanotransduction

Mechanical forces are potent modulators of stem cell based tissue regenerative mechanisms, inducing cell fate decisions and tissue specific commitment. A unique platform for investigating mechanotransduction is spaceflight, where microgravity and altered fluid mechanics provide a loading-null experimental condition. Seminal investigations of regenerative capacity in a wholly regenerative species, the newt model, and in a variety of totipotent and adult stem cell populations have demonstrated the detrimental effects of unloading on maintenance of stem cell based regeneration. Of particular interest is the observation that unloading interferes with the transition of stem cell pools from proliferative state to differentiation commitment. In this work we sought to test the hypothesis that gravity mechanotransduction regulates stem cell tissue regenerative processes by modulating stem cell proliferation and differentiation fates at specific cell cycle stages. To do this, clonally-derived ESCs were plated on a collagen matrix and expanded for 36 hours before re-plating on a non-adherent culture dish in the absence of leukemia inhibitory factor (LIF) to form spheroid aggregate EBs. After formation, the EBs were transferred to a collagen matrix coated culture dishes and given 4 days to allow implantation and outgrowth. In parallel, totipotent ESCs were plated 24 hours before mechanical stimulation on collagen matrix culture dishes in the presence of LIF to maintain totipotency and serve as un-differentiation committed controls. The EBs and ESCs were then subjected to either a 60 minute pulse of gravity (static loading) or 60 minutes of cyclic stretch (dynamic loading) mechanotransduction. Six hours post-stimulation, we used a 10X Genomics Single Cell controller to generate bar-coded single cell Illumina libraries and sequenced expressomes for 5,000 static loaded cells, representative of a change in gravity mechanotransduction, 5,000 dynamic loaded cells, representative of tissue loading associate with physiologic function, and 5,000 unstimulated 1g control cells. The comparison of these 3 libraries by cluster assignment based on like gene expression patterns show substantial alteration in cluster geometry due to mechanical loading. Specifically the mechanically loaded EB outgrowth cells to retain potency markers (PAX6, SOX2, CD34) and suppress early commitment markers (Dhh, VCAN, Igf1). Whereas the EBs cultured under the non-stimulated conditions display clear departure from the ESC expressome with lineage commitment markers upregulated and several tissue specific markers being expressed (BMP "early musculoskeletal development, Mesp1" early cardiovascular cell lineage). These markers are not seen in the mechano-stimulated cultures or the totipotent ESC cultures. Comparison of like clusters between our experimental conditions revealed an array of regenerative and stem cell genes are significantly mechano-regulated. Of particular importance CDKN1a/p21, a gene shown by previous investigation of our research team to be significantly upregulated in unloading, was suppressed in the static and dynamic loaded EBS. In addition to CDKN1a/p21 many genes related to cell cycle and transitory differentiation markers had elevated expression in the mechano-stimulated EBs, but surprisingly these trends were not observed in the ESC cultures. This study is the first of its kind investigating for mechano-signaling and mechano-regulated pathways, and has alre

Single cell RNA-sequecing↗

Comparisons of the polypeptide chains of globins.

Discussion of the amino acid differences and minimum base differences per codon due to mutations which took place during divergent evolution of vertebrates from a common ancestral gene. The ?random mutation model' of evolution is examined by comparing a carp alpha Hb chain and six mammalian alpha Hb chains in terms of the genetic code. The occurrence of recognizable three-base changes is analyzed and a summary of the distribution of changes in the hemoglobin and myoglobin chains is given for 148 sites.

Jukes, T. H.↗