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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Potential flow analysis of glaze ice accretions on an airfoil

The results of an analytical/experimental study of the flow fields about an airfoil with leading edge glaze ice accretion shapes are presented. Tests were conducted in the Icing Research Tunnel to measure surface pressure distributions and boundary layer separation reattachment characteristics on a general aviation wing section to which was affixed wooden ice shapes which approximated typical glaze ice accretions. Comparisons were made with predicted pressure distributions using current airfoil analysis codes as well as the Bristow mixed analysis/design airfoil panel code. The Bristow code was also used to predict the separation reattachment dividing streamline by inputting the appropriate experimental surface pressure distribution.

Zaguli, R. J.↗
NASA NTRS1992

High spatial resolution measurements in a single stage ram accelerator

High spatial resolution experimental tube wall pressure measurements of ram accelerator gas dynamic phenomena are presented in this paper. The ram accelerator is a ramjet-in-tube device which operates in a manner similar to that of a conventional ramjet. The projectile resembles the centerbody of a ramjet and travels supersonically through a tube filled with a combustible gaseous mixture, with the tube acting as the outer cowling. Pressure data are recorded as the projectile passes by sensors mounted in the tube wall at various locations along the tube. Utilization of special highly instrumented sections of tube has allowed the recording of gas dynamic phenomena with high resolution. High spatial resolution tube wall pressure data from the three regimes of propulsion studied to date (subdetonative, transdetonative, and superdetonative) in a single stage gas mixture are presented and reveal the three-dimensional character of the flow field induced by projectile fins and the canting of the fins and the canting of the projectile body relative to the tube wall. Also presented for comparison to the experimental data are calculations made with an inviscid, three-dimensional CFD code. The knowledge gained from these experiments and simulations is useful in understanding the underlying nature of ram accelerator propulsive regimes, as well as assisting in the validation of three-dimensional CFD coded which model unsteady, chemically reactive flows.

Hinkey, J. B.↗
NASA NTRS2005

Laminar Heating Validation of the OVERFLOW Code

OVERFLOW, a structured finite difference code, was applied to the solution of hypersonic laminar flow over several configurations assuming perfect gas chemistry. By testing OVERFLOW's capabilities over several configurations encompassing a variety of flow physics a validated laminar heating was produced. Configurations tested were a flat plate at 0 degrees incidence, a sphere, a compression ramp, and the X-38 re-entry vehicle. This variety of test cases shows the ability of the code to predict boundary layer flow, stagnation heating, laminar separation with re-attachment heating, and complex flow over a three-dimensional body. In addition, grid resolutions studies were done to give recommendations for the correct number of off-body points to be applied to generic problems and for wall-spacing values to capture heat transfer and skin friction. Numerical results show good comparison to the test data for all the configurations.

Lillard, Randolph P.↗
NASA NTRS2012

Longwave Band-by-band Cloud Radiative Effect and its Application in GCM Evaluation

The cloud radiative effect (CRE) of each longwave (LW) absorption band of a GCM fs radiation code is uniquely valuable for GCM evaluation because (1) comparing band-by-band CRE avoids the compensating biases in the broadband CRE comparison and (2) the fractional contribution of each band to the LW broadband CRE (f(sub CRE)) is sensitive to cloud top height but largely insensitive to cloud fraction, presenting thus a diagnostic metric to separate the two macroscopic properties of clouds. Recent studies led by the first author have established methods to derive such band ]by ]band quantities from collocated AIRS and CERES observations. We present here a study that compares the observed band-by-band CRE over the tropical oceans with those simulated by three different atmospheric GCMs (GFDL AM2, NASA GEOS-5, and CCCma CanAM4) forced by observed SST. The models agree with observation on the annual ]mean LW broadband CRE over the tropical oceans within +/-1W/sq m. However, the differences among these three GCMs in some bands can be as large as or even larger than +/-1W/sq m. Observed seasonal cycles of f(sub CRE) in major bands are shown to be consistent with the seasonal cycle of cloud top pressure for both the amplitude and the phase. However, while the three simulated seasonal cycles of f(sub CRE) agree with observations on the phase, the amplitudes are underestimated. Simulated interannual anomalies from GFDL AM2 and CCCma CanAM4 are in phase with observed anomalies. The spatial distribution of f(sub CRE) highlights the discrepancies between models and observation over the low-cloud regions and the compensating biases from different bands.

Huang, Xianglei↗
DOE OSTI2022

Predicting deviations in ex-core detector response in Krsko NPP during an earthquake with Monte-Carlo neutron transport methods

To study the response of ex-core detectors during an earthquake, an approach using Monte Carlo code MCNP was developed. As a first step to predict the response of the ex-core detectors a detailed MCNP model of the Krsko reactor core was created. A detailed MCNP core model was used to calculate 3D power distributions inside the core. Calculated power distributions were verified by comparison to the CORD-2 calculations, which is currently used for core design calculation verification of the Krsko nuclear power plant. Due to the outside the core position of the neutron detectors, the hybrid code ADVANTG is used to generate weight windows to speed up neutron transport outside the reactor core. To be able to use ADVANTG, fixed neutron source had to be reconstructed from the criticality core calculation. It was found that the power range ex-core detector movement has a negligible effect on the value of the thermal neutron flux at the active part of the detector. However, the 5 cm radial movement of the intermediate-range detector leads to 7 % - 8 % change in thermal neutron flux within the active intermediate-range detector region. The analysis continued with the evaluation of the effect of core barrel movement on the ex-core detector response. It was determined that the 2 mm core barrel radial oscillation can lead to ∼ 4% change in thermal neutron flux within the active detector region. It was evaluated that reactor core movement inside the baffle can lead up to ∼ 1.9% and ∼ 1.3% deviation in power-range and intermediate-range detector signal respectively. Analysis showed that the mechanical movement of ex-core neutron detectors could not explain the fluctuations in the ex-core detector signal. However, combined core barrel and reactor core inside baffle oscillations could be a probable reason for the observed fluctuations in the ex-core detector signal during an earthquake. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1996

Secondary particle contribution to LET spectra on LDEF

Four experiments utilizing passive detectors (P0006, P0004, A0015, M0004) were flown on LDEF to study the radiation environment. These experiments have been summarized in a companion paper (Benton et al., 1996). One of the experimental goals was to measure LET spectra at different locations and shielding depths with plastic nuclear track detectors (PNTD). It was found that the LET spectra extended well above the LET cutoff imposed by the geomagnetic field on GCR particle penetration into LEO. The high LET particles detected were mostly short-range (range < 2000 m), indicating that they were secondaries produced locally within the PNTD. The presence of these high LET particle fluences is important for the determination of dose equivalent because of the high Quality Factors (Q) involved. A relatively small fraction of particle fluence can contribute a large fraction of dose equivalent. Short-range, inelastic secondary particles produced by trapped protons in the South Atlantic Anomaly (SAA) were found to be a major contributor to the LET spectra above 100 keV/micrometer. The LET spectra were found to extend beyond the approximately 137 keV/micrometer relativistic GCR Fe peak to over 1000 keV/micrometer. The high LET tail of the LET spectra was measured in CR-39 and polycarbonate PNTDs using different techniques. GCR made a relatively modest contribution to the LET spectra as compared to the contributions from short-range secondary particles and stopping protons. LET spectra intercomparisons were made between LDEF measurements and exposures to 154 MeV accelerated proton beams. The similarities support the role of nuclear interactions by trapped protons as the major source of secondary particles in the PNTDs. Also techniques were employed to reduce the range cutoff for detection of the short-range secondaries to approximately 1 micrometer, so that essentially all secondary particles were included in the LET spectra. This has allowed a more realistic assessment of secondary contribution to dose equivalent. Comparisons of measured and calculated LET spectra have been made that demonstrate the need for more accurate modeling of secondary particles in radiation transport codes. Comparisons include preliminary calculations in which attempts have been made to include secondary particles.

long duration↗
DOE OSTI2025

Convergence study of wakefield simulations with GdfidL and ECHO3D

The interaction of charged particle beams with vacuum chamber components gives rise to electromagnetic wakefields, whose frequency-domain representation is known as beam coupling impedance. Geometric impedance arising from discontinuities and transitions in the vacuum chamber is the focus of this study. Minimizing this impedance is essential to mitigate adverse collective effects in modern storage rings operating with high-intensity particle beams. Accurate and reliable impedance simulations is a key factor of the vacuum chamber design. This paper presents the results of a convergence study of two widely used electromagnetic solvers, GdfidL and ECHO3D, applied to key vacuum-chamber components of the National Synchrotron Light Source II (NSLS-II) storage ring. Detailed comparisons are performed for several geometries, including flange absorbers, RF bellows, button-type beam position monitors, and an in-vacuum undulator (IVU). The results show notable differences in convergence and computational efficiency between the two codes. While GdfidL provides highly resolved results and serves as a common benchmark tool, ECHO3D yields consistent results with coarser meshes, significantly reducing simulation time and memory demands. Simulations with a full-geometry IVU model demonstrate that simplified taper-transition models can miss important impedance contributions. In conclusion, these findings provide practical guidelines for efficient and accurate impedance modeling to optimize design of vacuum chamber components for accelerators.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Code validation for the simulation of supersonic viscous flow about the F-16XL

The viewgraphs and discussion on code validation for the simulation of supersonic viscous flow about the F-16XL are provided. Because of the large potential gains related to laminar flow on the swept wings of supersonic aircraft, interest in the applications of laminar flow control (LFC) techniques in the supersonic regime has increased. A supersonic laminar flow control (SLFC) technology program is currently underway within NASA. The objective of this program is to develop the data base and design methods that are critical to the development of laminar flow control technology for application to supersonic transport aircraft design. Towards this end, the program integrates computational investigations underway at NASA Ames-Moffett and NASA Langley with flight-test investigations being conducted on the F-16XL at the NASA Ames-Dryden Research Facility in cooperation with Rockwell International. The computational goal at NASA Ames-Moffett is to integrate a thin-layer Reynolds averaged Navier-Stokes flow solver with a stability analysis code. The flow solver would provide boundary layer profiles to the stability analysis code which in turn would predict transition on the F-16XL wing. To utilize the stability analysis codes, reliable boundary layer data is necessary at off-design cases. Previously, much of the prediction of boundary layer transition has been accomplished through the coupling of boundary layer codes with stability theory. However, boundary layer codes may have difficulties at high Reynolds numbers, of the order of 100 million, and with the current complex geometry in question. Therefore, a reliable code which solves the thin-layer Reynolds averaged Navier-Stokes equations is needed. Two objectives are discussed, the first in greater depth. The first objective is method verification, via comparisons of computations with experiment, of the reliability and robustness of the code. To successfully implement LFC techniques to the F-16XL wing, the flow about the leading edge must be maintained as laminar flow. Therefore, the second objective is to focus on a series of numerical simulations with different values of angle of attack, alpha, and Reynolds numbers. The purpose of the simulations is to study their effects on the two main factors which precipitate transition to turbulence at leading edges of highly swept wings (e.g., 'spanwise contamination' and 'crossflow instability').

Flores, Jolen↗
NASA NTRS1997

High Resolution Far-Infrared Studies

We have been obtained high-resolution data (20 ft at 50 microns and 30 ft at 100 microns) on the KAO using Paul Harvey's 2 x 10 element photometer in both scanning and nodding modes. The practical flux limit for scanning is about 100 Jy. For fainter sources, a nodding (beam-switching) mode, which spend more time on the source, is used. This technique has been used successfully on objects as faint as 10 Jy; the 1 sigma noise for a 1 hour integration is about 1 Jy. Although not as sensitive as space-based instruments, the higher spatial resolution afforded by the KAO is essential in studying the far-infrared emission associated with young stars; in several cases we have been able to distinguish emission from multiple sources which were blended in the IRAS beam. In addition, comparison of fluxes in the KAO beam to those in the much larger IRAS beam provides information on the extended low-level emission arising from the surrounding region. We have developed a number of codes for producing model intensity distributions.

Mundy, Lee G.↗
NASA NTRS2005

Computation of Sound Propagation by Boundary Element Method

This report documents the development of a Boundary Element Method (BEM) code for the computation of sound propagation in uniform mean flows. The basic formulation and implementation follow the standard BEM methodology; the convective wave equation and the boundary conditions on the surfaces of the bodies in the flow are formulated into an integral equation and the method of collocation is used to discretize this equation into a matrix equation to be solved numerically. New features discussed here include the formulation of the additional terms due to the effects of the mean flow and the treatment of the numerical singularities in the implementation by the method of collocation. The effects of mean flows introduce terms in the integral equation that contain the gradients of the unknown, which is undesirable if the gradients are treated as additional unknowns, greatly increasing the sizes of the matrix equation, or if numerical differentiation is used to approximate the gradients, introducing numerical error in the computation. It is shown that these terms can be reformulated in terms of the unknown itself, making the integral equation very similar to the case without mean flows and simple for numerical implementation. To avoid asymptotic analysis in the treatment of numerical singularities in the method of collocation, as is conventionally done, we perform the surface integrations in the integral equation by using sub-triangles so that the field point never coincide with the evaluation points on the surfaces. This simplifies the formulation and greatly facilitates the implementation. To validate the method and the code, three canonic problems are studied. They are respectively the sound scattering by a sphere, the sound reflection by a plate in uniform mean flows and the sound propagation over a hump of irregular shape in uniform flows. The first two have analytical solutions and the third is solved by the method of Computational Aeroacoustics (CAA), all of which are used to compare the BEM solutions. The comparisons show very good agreements and validate the accuracy of the BEM approach implemented here.

Guo, Yueping↗
NASA NTRSDate not supplied

Transition prediction in hypersonic regime on complex geometries with RANS-based models

In the near future, RANS computations will continue to play a significant role in the design of hypersonic vehicles with complex geometries. Therefore, it is imperative to continue testing, benchmarking, and refining the RANS models. In the present study, we evaluate RANS-like, transport equations-based models for predicting laminar-turbulent transition over a full-scale scale model of the BOLT flight configuration that was tested in the CUBRC LENS-II wind tunnel test facility. Based on the availability of the experimental results, comparisons are made between computations from several pre-existing transition models and Computational Fluid Dynamics (CFD) codes, with an emphasis on using the same computational meshes and flow conditions for all computations. The analysis covers the sensitivity of the transition predictions to the input parameters for five different transition models (four mainstream models designed for low-speed flows and a uniquely high-speed model), grid resolution, and the details of model implementation across three different flow solvers. The results show that the phenomenological models can describe significant aspects of the measured transition front. However, a number of additional improvements are required before these models can offer more reliable estimates of transition in high-speed flows.

CFD↗
NASA NTRSDate not supplied

Transition Prediction in Hypersonic Regime on Complex Geometries With Rans-Based Models

In the near future, RANS computations will continue to play a significant role in the design of hypersonic vehicles with complex geometries. Therefore, it is imperative to continue testing, benchmarking, and refining the RANS models. In the present study, we evaluate RANS-like, transport equations-based models for predicting laminar-turbulent transition over a full-scale scale model of the BOLT flight configuration that was tested in the CUBRC LENS-II wind tunnel test facility. Based on the availability of the experimental results, comparisons are made between computations from several pre-existing transition models and Computational Fluid Dynamics (CFD) codes, with an emphasis on using the same computational meshes and flow conditions for all computations. The analysis covers the sensitivity of the transition predictions to the input parameters for five different transition models (four mainstream models designed for low-speed flows and a uniquely high-speed model), grid resolution, and the details of model implementation across three different flow solvers. The results show that the phenomenological models can describe significant aspects of the measured transition front. However, a number of additional improvements are required before these models can offer more reliable estimates of transition in high-speed flows.

CFD↗
NASA NTRS1981

Ground effect hover characteristics of a large-scale twin tilt-nacelle V/STOL model

This paper is a summary of an analysis of the ground-effect characteristics of a large-scale twin-engine, tilt-nacelle V/STOL model. The analysis considers data from the flow field beneath the full-scale model, as well as small-scale model test data, and makes comparisons with jet-ground interactions predicted by a computer code. The data from the large-scale test comprise ground-plane surface temperatures, static pressure distribution and wall-jet total-pressure profiles, fuselage undersurface static pressures, and model forces and moments. The results indicate that the near-field flow is more complex than is indicated by either the small-scale uniform jet studies or the computer predictions. The far-field flow characteristics do show some similarity for these three cases.

Dudley, M. R.↗
NASA NTRS1992

Optimal design of geodesically stiffened composite cylindrical shells

An optimization system based on general-purpose finite element code CSM Testbed and optimization program ADS is described. The system can be used to obtain minimum-mass designs of composite shell structures with complex stiffening arrangements. Ply thicknesses, ply orientations, and stiffener heights can be used as design variables. Buckling, displacement, and material failure constraints can be imposed on the design. The system is used to conduct a preliminary design study of geodesically stiffened shells. For comparison purposes, optimal designs of unstiffened shells, and ring and longitudinal stringer stiffened shells are also studied. Trends in the design of geodesically stiffened shells are identified. Features that enhance the capabilities and efficiency of the design system are described.

Gendron, G.↗
NASA NTRS1993

Analytical and experimental studies of a short compact subsonic diffuser for a two-dimensional supersonic inlet

An experimental study of a two-dimensional supersonic inlet with a short compact subsonic diffuser, length to exit diameter (dl/d) ratio of 1.25, was conducted to investigate the impact of the short diffuser on inlet performance at low speeds and to assess the diffuser subsonic performance for a simulated diffuser flow corresponding to high-speed inlet conditions near the design flight Mach number of 2.2. For the low-speed testing, a drooped lip was employed to improve the inlet performance at a high angle of attack. For the simulated high-speed testing, air was blown through slots or discrete nozzles as an active boundary-layer control. The results from the low-speed performance test were compared with the results from a previous test program on the same inlet with a long subsonic diffuser (dl/d = 4.5). The comparison indicates that inlet recovery was not affected by the use of the short diffuser for either the baseline (no droop) or the drooped cowl lip configuration. However, the inlet baseline distortion for the short diffuser configuration was substantially higher than for the long diffuser. A comparison of the two configurations with a 70 deg drooped lip showed no significant difference in distortion. For the portion of the experimental program in which diffuser conditions for high-speed flight were simulated, diffuser-induced flow separation occurred. This separation was predicted from an analytical study that used the Hess potential flow panel method and the Herring two-dimensional boundary-layer analysis computer codes. The flow separated mainly on the diffuser ramp. Subsequent tests in which boundary-control systems were utilized showed that blowing with either slots or discrete nozzles could suppress the flow separation in the short subsonic diffuser, thereby substantially improving the diffuser performance.

Iek, Chanthy↗
NASA NTRS1998

Development of Curved-Plate Elements for the Exact Buckling Analysis of Composite Plate Assemblies Including Transverse Shear Effects

The analytical formulation of curved-plate non-linear equilibrium equations that include transverse-shear-deformation effects is presented. A unified set of non-linear strains that contains terms from both physical and tensorial strain measures is used. Using several simplifying assumptions, linearized, stability equations are derived that describe the response of the plate just after bifurcation buckling occurs. These equations are then modified to allow the plate reference surface to be located a distance z(c), from the centroid surface which is convenient for modeling stiffened-plate assemblies. The implementation of the new theory into the VICONOPT buckling and vibration analysis and optimum design program code is described. Either classical plate theory (CPT) or first-order shear-deformation plate theory (SDPT) may be selected in VICONOPT. Comparisons of numerical results for several example problems with different loading states are made. Results from the new curved-plate analysis compare well with closed-form solution results and with results from known example problems in the literature. Finally, a design-optimization study of two different cylindrical shells subject to uniform axial compression is presented.

McGowan, David M.↗
NASA NTRS2014

Steady and Unsteady Nozzle Simulations Using the Conservation Element and Solution Element Method

This paper presents results from computational fluid dynamic (CFD) simulations of a three-stream plug nozzle. Time-accurate, Euler, quasi-1D and 2D-axisymmetric simulations were performed as part of an effort to provide a CFD-based approach to modeling nozzle dynamics. The CFD code used for the simulations is based on the space-time Conservation Element and Solution Element (CESE) method. Steady-state results were validated using the Wind-US code and a code utilizing the MacCormack method while the unsteady results were partially validated via an aeroacoustic benchmark problem. The CESE steady-state flow field solutions showed excellent agreement with solutions derived from the other methods and codes while preliminary unsteady results for the three-stream plug nozzle are also shown. Additionally, a study was performed to explore the sensitivity of gross thrust computations to the control surface definition. The results showed that most of the sensitivity while computing the gross thrust is attributed to the control surface stencil resolution and choice of stencil end points and not to the control surface definition itself.Finally, comparisons between the quasi-1D and 2D-axisymetric solutions were performed in order to gain insight on whether a quasi-1D solution can capture the steady and unsteady nozzle phenomena without the cost of a 2D-axisymmetric simulation. Initial results show that while the quasi-1D solutions are similar to the 2D-axisymmetric solutions, the inability of the quasi-1D simulations to predict two dimensional phenomena limits its accuracy.

Nozzle gross thrust calculations↗
DOE OSTI2022

Fast Reactor Physics Model Verification Studies using ARC and PyARC Workflows

PyARC was recently developed at Argonne National Laboratory to automate many of the tasks required in the ARC (Argonne Reactor Computation) fast reactor simulation workflow, from input file generation, code execution, data transfer between ARC codes, and output postprocessing. PyARC will likely be the path forward to train new users of the ARC codes with the goal of wide adoption by the national laboratories, academia, and industry. In particular, for the ANL-JAEA collaboration under the Civil Nuclear Working Group (CNWG) project agreement NE-01, PyARC will be used to model the Joyo and EBR-II reactors for comparisons with measured data and calculated results from JAEA (Task 3: Fast Reactor Fuel and Core). As an additional avenue for verification and validation, this report investigates the use of PyARC towards a variety of existing ARC-based reactor models, in order to understand its efficacy in replicating the behavior of base ARC codes and better understand any limitations within modeling realistic fast reactor problems. To this end, PyARC was used to model the Joyo MKI, RBEC Benchmark-M, PRISM Mod-B, and EBR-II Run 138B cores, and its results were compared to those from existing ARC-based models. It was found that for hexagonal-based geometries PyARC was able to replicate the behavior of ARC codes to within 10 pcm for small reactor cores, and ~150pcm difference in eigenvalue for larger cores. These discrepancies are attributed primarily to differences in local mesh refinement options between ARC and PyARC, which currently cannot be resolved with PyARC’s latest version (1.6.0). In some of these cases, PyARC was used to model steady-state problems with initial core compositions originating from a prior REBUS depletion calculation. While PyARC was not designed to support such steady-state calculations, workarounds were applied to replicate the behavior of ARC-based calculations as closely as possible. Thus, these results demonstrate the wide extent to which they can be applied to fast reactor problems while still providing immense benefit to the user in terms of automating and standardizing common routines within the fast reactor analysis workflow. This study concluded that PyARC will be suitable for modeling the steady-state conditions of the EBR-II and Joyo fast reactors as part of the CNWG project agreement.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗