Search NASA⌕ Search

SEARCH · Search NASA

Results for “code comparison study”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Experimental studies of 90 deg corner cascades in the National Full-Scale Aerodynamic Complex

Two 90 deg, fixed geometry, turning vane cascades were replaced in the 40- by 80-foot wind tunnel circuit of the National Full-Scale Aerodynamic Complex as a part of a major facility upgrade. Prior to the replacement of these vanes, small-scale, two-dimensional tests were conducted to assess the aerodynamic performance of the original, 1930's technology vane geometry and a new vane geometry which was developed using a modern cascade design code. The small scale tests were conducted in the 1/10th-Scale Cascade Test Facility and included measurements of total pressure, static pressure, and flow angularity distributions in the wake of the cascade and static pressure distributions along the vane profiles. Test results and a detailed comparison of the two vanes designs are presented. Based on these results, the new design was chosen and subsequently installed in the wind tunnel. During the integrated systems test of the upgraded facility, two of the 52 vanes in the first cascade were instrumented near midspan. The full-scale test results are presented, these show the cascade pressure loss coefficient to be slightly larger than obtained in the small-scale tests. Elements of the three-dimensional flow field upstream of the cascade, which effect can impact cascade performance, are discussed and recommendations for future research are given.

Meyn, Larry A.↗
DOE OSTI2024

High- and Mid-Fidelity Modeling Comparison for a Floating Marine Turbine System: Preprint

There is a lack of suitable numerical tools, particularly open-source tools, that can be used for designing and optimizing marine turbine systems. The National Renewable Energy Laboratory (NREL) has added features to their widely used mid-fidelity wind turbine modeling code, OpenFAST, to enable modeling of axial flow marine turbines. This necessitated the addition of several physical effects relevant to marine turbines that can be neglected for wind turbines. These include buoyant loads, added mass and inertia loads, wave-current superposition, and changes to the coordinate systems. This updated version of OpenFAST allows for the modeling of both fixed and floating marine turbine systems at a speed comparable to real time. While efficient for long simulations, large sets of load cases, and design studies, mid-fidelity codes cannot capture all of the potentially important physical phenomenon impacting marine turbine systems. High-fidelity computational fluid dynamics (CFD) simulations can capture more flow effects with fewer assumptions and provide detailed body pressure mapping and flow-field information. It is important to compare predictions between mid-fidelity and high-fidelity codes, both to verify the models and to understand the limitations. A floating marine turbine system designed by NREL was modeled both with OpenFAST and with the commercial CFD code, STARCCM+. The CFD model used a 3-D unsteady Reynolds-averaged Navier-Stokes solver for a volume-of-fluid numerical wave and current tank. The blade-resolved simulations used the sliding-interface technique for the spinning rotor and an overset grid to accommodate the rigid-body motion of the floating system. The mooring system was modelled with a custom coupling of the CFD solver with the open-source code, MoorDyn. This improves upon the existing quasi-static catenary solver in STARCCM+, which lacks seabed contact or line-to-line connections. Spatial and temporal convergence studies were conducted. The simulation results for a combined current and wave condition are compared between OpenFAST and CFD, highlighting the capabilities of the mid-fidelity code and identifying the areas where a high-fidelity approach is needed.

CFD↗
NASA NTRS2006

Comparison of Damage Path Predictions for Composite Laminates by Explicit and Standard Finite Element Analysis Tools

Splitting, ultimate failure load and the damage path in center notched composite specimens subjected to in-plane tension loading are predicted using progressive failure analysis methodology. A 2-D Hashin-Rotem failure criterion is used in determining intra-laminar fiber and matrix failures. This progressive failure methodology has been implemented in the Abaqus/Explicit and Abaqus/Standard finite element codes through user written subroutines "VUMAT" and "USDFLD" respectively. A 2-D finite element model is used for predicting the intra-laminar damages. Analysis results obtained from the Abaqus/Explicit and Abaqus/Standard code show good agreement with experimental results. The importance of modeling delamination in progressive failure analysis methodology is recognized for future studies. The use of an explicit integration dynamics code for simple specimen geometry and static loading establishes a foundation for future analyses where complex loading and nonlinear dynamic interactions of damage and structure will necessitate it.

Bogert, Philip B.↗
DOE OSTI2022

High-yield recombinant bacterial expression of 13 C-, 15 N-labeled, serine-16 phosphorylated, murine amelogenin using a modified third generation genetic code expansion protocol

Amelogenin constitutes ~ 90% of the enamel matrix in the secretory stage of amelogenesis, a still poorly understood process that results in the formation of the hardest and most mineralized tissue in vertebrates – enamel. Most biophysical research with amelogenin uses recombinant protein expressed in Escherichia coli. In addition to providing copious amounts of protein, recombinant expression allows 13 C- and 15 N-labelling for detailed structural studies using NMR spectroscopy. However, native amelogenin is phosphorylated at one position, Ser-16 in murine amelogenin, and a lingering question overshadowing decades of enamel research with recombinant amelogenin is the effect of phosphorylation on its biophysical properties. Using a modified genetic code expansion protocol we have expressed and purified uniformly 13 C-, 15 N-labeled murine amelogenin (pS16M179) with ~95% of the protein being correctly phosphorylated. Homogeneous phosphorylation was achieved using commercially available, enriched, 13 C-, 15 N-labeled media and protein expression was induced with isopropyl β-D-1-thiogalactopyranoside at 310 K. Phosphoserine incorporation was verified from one-dimensional 31P NMR spectra, comparison of 1 H- 15 N HSQC spectra, Phos-tag SDS PAGE, and mass spectrometry. Phosphorus-31 NMR spectra for pS16M179 under conditions known to trigger amelogenin self-assembly into nanospheres confirm nanosphere models with buried N-termini. Lambda phosphatase treatment of these nanospheres results in the dephosphorylation of pS16M179, confirming that smaller oligomers and monomers with exposed N-termini are in equilibrium with nanospheres. Such 13 C-, 15 N-labeling of amelogenin with accurately encoded phosphoserine incorporation will accelerate biomineralization research to understand amelogenesis and stimulate the expanded use of genetic code expansion protocols to introduce phosphorylated amino acids into proteins.

31P-NMR↗
NASA NTRS1981

The seasonal-cycle climate model

The seasonal cycle run which will become the control run for the comparison with runs utilizing codes and parameterizations developed by outside investigators is discussed. The climate model currently exists in two parallel versions: one running on the Amdahl and the other running on the CYBER 203. These two versions are as nearly identical as machine capability and the requirement for high speed performance will allow. Developmental changes are made on the Amdahl/CMS version for ease of testing and rapidity of turnaround. The changes are subsequently incorporated into the CYBER 203 version using vectorization techniques where speed improvement can be realized. The 400 day seasonal cycle run serves as a control run for both medium and long range climate forecasts alsensitivity studies.

Marx, L.↗
NASA NTRS1985

Analysis of turbulent underexpanded jets. II - Shock noise features using SCIPVIS

SCIPVIS, the computational model discussed by Dash et al. (1985), is assessed in predicting the complicated flow structure associated with shock-containing plumes. In addition, the analysis in this study examines this code's applicability as a basic part of a program for estimating broadband shock noise radiation. The results of this study show that excellent agreement exists between predicted and measured static pressure distributions for both underexpanded and overexpanded flow cases considered. Of the three turbulence closure models incorporated in the SCIPVIS code, the kW model of Spalding produces the most uniform agreement with measurement. The k-epsilon-2 model of Launder consistently overestimates plume spreading for supersonic jets with exit Mach numbers in the 1-2 range. Dash's (1983) k-epsilon-2-cc, compressibility-corrected version of Launder's model underestimates plume spreading. Good qualitative agreement was also obtained between the measured longitudinal turbulence intensity and that predicted by the code for the same trial case. Comparison of measured and predicted broadband shock noise spectrum peak values were found to be in excellent agreement. This utilized a variant of the Harper-Bourne and Fisher (1973) phase-array model: the effective shock spacing was reinterpreted as the value of the end of the plume potential core, determined herein by the SCIPVIS code.

Seiner, J. M.↗
NASA NTRSDate not supplied

Using the Diurnal Variability in GeoNEX TOA Reflectances for Earth Monitoring

Observations from the third-generation geostationary satellite instruments (GOES 16/17 ABI, Himawari 8/9 AHI, and etc.) have spatial resolution and spectral band configurations comparable to flagship LEO sensors (e.g., MODIS/VIIRS). More importantly, these data are acquired at very high temporal resolution, faithfully recording the variations of the full disk of Earth at every 5-10 minutes. They thus provide unique information about Earth’s atmosphere and surface. In order to explore the unique information content of geostationary data, this study systematically analyzes the diurnal variability in the GeoNEX L1G TOA reflectance products and compares them to simulated results by state-of-the-art radiative transfer codes. Our results show that • The smoothness of the TOA reflectance diurnal cycle provides a convenient and reliable way to identify stable atmospheric conditions and filter out passing clouds/shadows. • The diurnal variability of the blue band (0.47µm) reflectance is regulated mainly by atmospheric optical conditions over a majority of land cover types. As such, the diurnal variability of the blue band data allows us to retrieve AOD without invoking the use of spectral band ratios (SRC) as in previous algorithms. • In comparison, the diurnal variability of the short-wave infrared band (2.2µm) BRFs is mainly regulated by surface reflectance and the sun-target-satellite geometry. This information allows us to test and, if suitable, retrieve surface BRDF parameters. • Spectral band ratios, especially those between the 2.2µm and 0.47µm bands, are not constant but vary by locations and sun-target-satellite geometries. Our analysis clearly demonstrates that the information provided in high-frequent geostationary observations is unique and complementary to LEO sensors. Therefore, a synergy of GEO and LEO (and other) sensors has the great potential to improve existing remote sensing models and algorithms for better Earth monitoring.

Diurnal Variability↗
NASA NTRS1992

Numerical studies of the fluid and optical fields associated with complex cavity flows

Numerical solutions for the flowfield about several cavity configurations have been computed using the Reynolds averaged Navier-Stokes equations. Comparisons between numerical and experimental results are made in two dimensions for free shear layers and a rectangular cavity, and in three dimensions for the transonic aero-window problem of the Stratospheric Observatory for Infrared Astronomy (SOFIA). Results show that dominant acoustic frequencies and magnitudes of the self excited resonant cavity flows compare well with the experiment. In addition, solution sensitivity to artificial dissipation and grid resolution levels are determined. Optical path distortion due to the flow field is modelled geometrically and is found to match the experiment. The fluid field was computed using a diagonalized scheme within an overset mesh framework. An existing code, OVERFLOW, was utilized with the additions of characteristic boundary condition and output routines required for reduction of the unsteady data. The newly developed code is directly applicable to a generalized three dimensional structured grid zone. Details are provided in a paper included in Appendix A.

Atwood, Christopher A.↗
NASA NTRS1996

Euler Technology Assessment for Preliminary Aircraft Design: Compressibility Predictions by Employing the Cartesian Unstructured Grid SPLITFLOW Code

The objective of the second phase of the Euler Technology Assessment program was to evaluate the ability of Euler computational fluid dynamics codes to predict compressible flow effects over a generic fighter wind tunnel model. This portion of the study was conducted by Lockheed Martin Tactical Aircraft Systems, using an in-house Cartesian-grid code called SPLITFLOW. The Cartesian grid technique offers several advantages, including ease of volume grid generation and reduced number of cells compared to other grid schemes. SPLITFLOW also includes grid adaption of the volume grid during the solution to resolve high-gradient regions. The SPLITFLOW code predictions of configuration forces and moments are shown to be adequate for preliminary design, including predictions of sideslip effects and the effects of geometry variations at low and high angles-of-attack. The transonic pressure prediction capabilities of SPLITFLOW are shown to be improved over subsonic comparisons. The time required to generate the results from initial surface data is on the order of several hours, including grid generation, which is compatible with the needs of the design environment.

Finley, Dennis B.↗
NASA NTRS1997

Efficient Gradient-Based Shape Optimization Methodology Using Inviscid/Viscous CFD

The formerly developed preconditioned-biconjugate-gradient (PBCG) solvers for the analysis and the sensitivity equations had resulted in very large error reductions per iteration; quadratic convergence was achieved whenever the solution entered the domain of attraction to the root. Its memory requirement was also lower as compared to a direct inversion solver. However, this memory requirement was high enough to preclude the realistic, high grid-density design of a practical 3D geometry. This limitation served as the impetus to the first-year activity (March 9, 1995 to March 8, 1996). Therefore, the major activity for this period was the development of the low-memory methodology for the discrete-sensitivity-based shape optimization. This was accomplished by solving all the resulting sets of equations using an alternating-direction-implicit (ADI) approach. The results indicated that shape optimization problems which required large numbers of grid points could be resolved with a gradient-based approach. Therefore, to better utilize the computational resources, it was recommended that a number of coarse grid cases, using the PBCG method, should initially be conducted to better define the optimization problem and the design space, and obtain an improved initial shape. Subsequently, a fine grid shape optimization, which necessitates using the ADI method, should be conducted to accurately obtain the final optimized shape. The other activity during this period was the interaction with the members of the Aerodynamic and Aeroacoustic Methods Branch of Langley Research Center during one stage of their investigation to develop an adjoint-variable sensitivity method using the viscous flow equations. This method had algorithmic similarities to the variational sensitivity methods and the control-theory approach. However, unlike the prior studies, it was considered for the three-dimensional, viscous flow equations. The major accomplishment in the second period of this project (March 9, 1996 to March 8, 1997) was the extension of the shape optimization methodology for the Thin-Layer Navier-Stokes equations. Both the Euler-based and the TLNS-based analyses compared with the analyses obtained using the CFL3D code. The sensitivities, again from both levels of the flow equations, also compared very well with the finite-differenced sensitivities. A fairly large set of shape optimization cases were conducted to study a number of issues previously not well understood. The testbed for these cases was the shaping of an arrow wing in Mach 2.4 flow. All the final shapes, obtained either from a coarse-grid-based or a fine-grid-based optimization, using either a Euler-based or a TLNS-based analysis, were all re-analyzed using a fine-grid, TLNS solution for their function evaluations. This allowed for a more fair comparison of their relative merits. From the aerodynamic performance standpoint, the fine-grid TLNS-based optimization produced the best shape, and the fine-grid Euler-based optimization produced the lowest cruise efficiency.

Baysal, Oktay↗
NASA NTRS2003

RSRM and ETM03 Internal Flow Simulations and Comparisons

ETM03 (Engineering Test Motor-03) is an extended length RSRM (Reusable Solid Rocket Motor) designed to increase motor performance and create more severe internal environments compared with the standard four-segment RSRM motor configuration. This is achieved primarily through three unique design features. First is the incorporation of an additional RSRM center segment, second is a slight increase in throat diameter, and third is the use of an Extended Aft Exit Cone (EAEC). As a result of these design features, parameters such as web time, action time, head end pressure, web time average pressure, maximum thrust, mass flow rate, centerline Mach number, pressure and thrust integrals have all increased compared with nominal RSRM values. In some cases these increases are substantial. The primary objective of the ETM03 test program is to provide a platform for RSRM component margin testing. Test results will not only provide direct data concerning component performance under more adverse conditions, but serve as a second design data point for developing, validating and enhancing component analytical modeling techniques. To help component designers assess how the changes in motor environment will affect performance, internal flow simulations for both the nominal RSRM and ETM03 motor designs were completed to obtain comparisons of aero-thermal boundary conditions and system loads. Full geometries for both motors were characterized with two-dimensional axi-symmetric models at burn times of 1, 20, 54, 67 and 80-seconds. A sixth set considered burn times of 110 and 117-seconds for RSRM and ETM03, respectively. The simulations were performed using the computational fluid dynamics (CFD) commercial code FLUENT (trademark). Of particular interest were any differences between the two motor environments that could lead to a significant increase in system loads, or in internal insulation and/or nozzle component charring and erosion in ETM03 compared with RSRM. Based on these comparative analyses conducted in this study, the objective of ETM03 will be achieved by providing a more adverse operating environment for motor components than the nominal RSRM environment. For example: Higher chamber pressure drop in ETM03 than in RSRM; higher centerline Mach numbers approaching the nozzle in ETM03 than in RSRM; higher heat transfer rates for the internal insulation and nozzle components in ETM03 than in RSRM; and higher levels of droplet impingement and slag accumulation in ETM03 than in the RSRM.

Ahmad, R. A.↗
NASA NTRS1997

Analysis of Emission Spectra from Arc-jet Shock Layer Flows

This paper reports computational comparisons with experimental studies of a nonequilibrium blunt body shock layer in a high enthalpy arc-jet wind tunnel at NASA Ames Research Center. The primary objective of this work is to investigate the existence of a thermochemical equilibrium region in the shock layer. The existence of such an equilibrium region is of interest for following reasons: (1) to understand the equilibration process behind the shock in an arc-jet flow environment; (2) to interpret measured surface heat transfer data for purpose of determining surface catalytic efficiency, and (3) to determine the total enthalpy from the spectroscopic measurements. The paper will present an analysis of the experimental data obtained in the arc-jet wind tunnel. Experimental data includes measurements of emission spectra of radiation emanating from a shock layer formed in front of a 6-inch flat-faced cylinder. The measurements, obtained using a two dimensional CCD camera mounted on a spectrograph, provide spatially resolved spectra along the stagnation streamline of the model. Computational analysis includes simulation of nonequilibrium flow in the arc-jet facility (flow in the conical nozzle and shock layer in front of a flat-faced cylinder) using 2-D/axisymmetric Navier-Stokes codes and prediction of the radiation spectra from the axisymmetric flowfield using NEQAIR radiation code. Various line-of-sight averaged flow properties such as vibrational and rotational temperatures, species number densities within the shock layer are deduced from the experimental spectra. Comparison of the computed and experimental line-of-sight averaged flow properties provides assessment of thermochemical equilibration processes in an arc-jet shock layer.

Gokcen, T.↗
DOE OSTI2022

A scalable compressible volume of fluid solver using a stratified flow model

Here, in this paper, a Volume of Fluid method which uses a Stratified Flow model for flux calculation is proposed to solve a compressible multiphase flow in a way that has high parallel efficiency. The solver is part of an open-source computing tool to study fuel entrainment and combustion in the application of paraffin-based hybrid rocket motors. This paper focuses on the compressible Volume of Fluid solver for the liquid and gaseous phases in this environment. This solver uses perfect gas and stiffened gas models for equation of state, with the ability to easily add other models as desired. Additionally a comparison of a Riemann solver versus an AUSM+up scheme is conducted, showing improvement with the later in the solution of multiphase flows. This solver successfully reproduces shock tube cases with the expected accuracy, as well as simulating two-dimensional shear and gravity-driven flows. Results also show that the code is able to utilize a Stratified Flow model to evaluate shear flow without the need for interface reconstruction or gradient calculations. Finally, it is shown that the solver has near-ideal scaling under strong scaling tests as well as good performance in static scaling, giving improved performance over current options and implying future advancement of high-performance computing efficiency for multiphase flow solvers.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

Sidewall boundary-layer effects in two-dimensional airfoil testing

Theoretical studies required to evaluate and validate the streamlined wall test section of the 0.3 m transonic Cryogenic Tunnel were initiated. The various aspects that are being considered presently are deviations of the real wall shape from the true streamline shape, wall adjustment strategy, and the influence of the sidewall boundary layers. Since the top and bottom walls are supported at a finite number of jack points, the true wall shape will be different from the ideal streamline shape. This is determined by calculating the structural shape for cases for which the exact streamline shape can be calculated. For the structural shape calculations, the MSCNASTRAN code is being used. To start with, comparisons are made using simple singularities for model representation and also for a flat plate at angle of attack.

Murthy, A. V.↗
NASA NTRS2015

Time-Accurate Unsteady Pressure Loads Simulated for the Space Launch System at Wind Tunnel Conditions

A transonic flow field about a Space Launch System (SLS) configuration was simulated with the Fully Unstructured Three-Dimensional (FUN3D) computational fluid dynamics (CFD) code at wind tunnel conditions. Unsteady, time-accurate computations were performed using second-order Delayed Detached Eddy Simulation (DDES) for up to 1.5 physical seconds. The surface pressure time history was collected at 619 locations, 169 of which matched locations on a 2.5 percent wind tunnel model that was tested in the 11 ft. x 11 ft. test section of the NASA Ames Research Center's Unitary Plan Wind Tunnel. Comparisons between computation and experiment showed that the peak surface pressure RMS level occurs behind the forward attach hardware, and good agreement for frequency and power was obtained in this region. Computational domain, grid resolution, and time step sensitivity studies were performed. These included an investigation of pseudo-time sub-iteration convergence. Using these sensitivity studies and experimental data comparisons, a set of best practices to date have been established for FUN3D simulations for SLS launch vehicle analysis. To the author's knowledge, this is the first time DDES has been used in a systematic approach and establish simulation time needed, to analyze unsteady pressure loads on a space launch vehicle such as the NASA SLS.

Alter, Stephen J.↗
DOE OSTI2025

Numerical simulation of coupled THM behaviour of full-scale EBS in backfilled experimental gallery in the Horonobe URL

Bentonite-based engineered barrier system (EBS) is a key component of many repository designs for the geological disposal of high-level radioactive waste. Given the complexity and interaction of the phenomena affecting the barrier system, coupled thermo-hydro-mechanical (THM) numerical analyses are a potentially useful tool for a better understanding of their behaviour. In this context, a Task (the Horonobe EBS experiment) was undertaken to study, using numerical analyses, the thermo-hydro-mechanical (and thermo-hydro) interactions in bentonite based engineered barriers within the international cooperative project DECOVALEX 2023. One full-scale in-situ experiment and four laboratory experiments, largely complementary, were selected for modelling. The Horonobe EBS experiment is a temperature-controlled non-isothermal experiment combined with artificial groundwater injection. The Horonobe EBS experiment consists of the heating and cooling phases. Six research teams performed the THM or TH (depended on research team approach) numerical analyses using a variety of computer codes, formulations and constitutive laws. Finally, for each experiment, the basic features of the analyses are described and the comparison between calculations and laboratory experiments and field observations are presented and discussed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1974

Large deflection elastic-plastic dynamic response of stiffened shells of revolution

This paper presents the formulation and check-out problems for a computer code DYNAPLAS, which analyzes the large deflection elastic-plastic dynamic response of stiffened shells of revolution. The formulation for spacial discretization is by the finite element method with finite differences being used for the evaluation of the pseudo forces due to material and geometric nonlinearities. Time integration is by the Houbolt method or central differences. The stiffeners may be due to concentrated or distributed eccentric rings and spring supports at arbitrary angles around the circumference of the elements. Check-out problems include the comparison of solutions from DYNAPLAS with experimental and other computer solutions for rings and conical and cylindrical shells. A hypothetical submarine including stiffeners and missile tube is studied under a combination of hydrostatic and dynamically applied asymmetrical pressure loadings.

Stricklin, J. A.↗
NASA NTRS1994

Validation of Vortex-Lattice Method for loads on wings in lift-generated wakes

A study is described that evaluates the accuracy of vortex-lattice methods when they are used to compute the loads induced on aircraft as they encounter lift-generated wakes. The evaluation is accomplished by use of measurements made in the 80- by 120-foot wind tunnel of the lift, rolling-moment, and downwash in the wake of three configurations of a model of a subsonic transport aircraft. The downwash measurements are used as input for a vortex-lattice code in order to compute the lift and rolling moment induced on wings that have a span of 0.186, 0.510, or 1.022 times the span of the wake-generating model. Comparison of the computed results with the measured lift and rolling moment distributions are used to determine the accuracy of the vortex-lattice code. It was found that the vortex-lattice method is very reliable as long as the span of the encountering of following wing is less than about 0.2 of the generator span. As the span of the following wing increases above 0.2, the vortex-lattice method continues to correctly predict the trends and nature of the induced loads, but it overpredicts the magnitude of the loads by increasing amounts. The increase in deviation of the computed from the measured loads with size of the following wing is attributed to the increase in distortion of the structure of the vortex wake as it approaches and passes the larger following wings.

Rossow, J.↗