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 289 records · Page 16

Efficient real space formalism for hybrid density functionals

We present an efficient real space formalism for hybrid exchange-correlation functionals in generalized Kohn–Sham density functional theory (DFT). In particular, we develop an efficient representation for any function of the real space finite-difference Laplacian matrix by leveraging its Kronecker product structure, thereby enabling the time to solution of associated linear systems to be highly competitive with the fast Fourier transform scheme while not imposing any restrictions on the boundary conditions. We implement this formalism for both the unscreened and range-separated variants of hybrid functionals. We verify its accuracy and efficiency through comparisons with established planewave codes for isolated as well as bulk systems. In particular, we demonstrate up to an order-of-magnitude speedup in time to solution for the real space method. We also apply the framework to study the structure of liquid water using ab initio molecular dynamics, where we find good agreement with the literature. Overall, the current formalism provides an avenue for efficient real-space DFT calculations with hybrid density functionals.

Chemistry↗

New constraints on warm dark matter from the Lyman- α forest power spectrum

The forest of Lyman-α absorption lines detected in the spectra of distant quasars encodes information on the nature and properties of dark matter and the thermodynamics of diffuse baryonic material. Its main observable—the 1D flux power spectrum (FPS)—should exhibit a suppression on small scales and an enhancement on large scales in warm dark matter (WDM) cosmologies compared to standard Λ⁢CDM. Here, we present an unprecedented suite of 1080 high-resolution cosmological hydrodynamical simulations run with the graphics processing unit-accelerated code cholla to study the evolution of the Lyman-α forest under a wide range of physically motivated gas thermal histories along with different free-streaming lengths of WDM thermal relics in the early Universe. A statistical comparison of synthetic data with the forest FPS measured down to the smallest velocity scales ever probed at redshifts 4.0≲z≲5.2 [E. Boera et al., Revealing reionization with the thermal history of the intergalactic medium: New constraints from the Ly⁢α flux power spectrum, Astrophys. J. 872, 101 (2019)] yields a lower-limit m WDM >3.1 keV (95% C.L.) for the WDM particle mass and constrains the amplitude and spectrum of the photoheating and photoionizing background produced by star-forming galaxies and active galactic nuclei at these redshifts. Interestingly, our Bayesian inference analysis appears to weakly favor WDM models with a peak likelihood value at the thermal relic mass of m WDM =4.5 keV. In conclusion, we find that the suppression of the FPS from free-streaming saturates at k≳0.1 s km -1 because of peculiar velocity smearing, and this saturated suppression combined with a slightly lower gas temperature provides a moderately better fit to the observed small-scale FPS for WDM cosmologies.

79 ASTRONOMY AND ASTROPHYSICS↗

Aeroelastic stability analysis of the AD-1 manned oblique-wing aircraft

The AD-1 manned flight test program was conducted to evaluate the stability, control and handling characteristics of oblique wing aircraft. The results of the aeroelastic stability analysis are presented for both the wing alone and the wing with ailerons. A comparison was made between the results obtained using the traditional k-method of flutter analysis and the results using the PK or British method of flutter analysis. Studies were performed using the latest version of the NASTRAN computer code as well as the PASS/FLUT program.

Rutkowski, M. J.↗

Radiative Heat Transport and Optical Characterization of High Temperature Molten Salts

The goal of this project was to advance the field of molten fluoride radiative heat transfer (RHT), through optical property measurements, heat transfer loop construction and operation, and computational fluid dynamics simulations. The following outcomes have been accomplished. A comprehensive review has been conducted to compile all available RHT-relevant optical absorption data for molten fluorides. The various measurement methods have also been compiled and critiqued. Thermal emission measurement methods have been refined and documented. A unique thermal emissivity measurement facility has been constructed and used to measure the emissivity of flowing FLiBe-exposed SS316 samples. A slight increase has been observed in the emissivity compared to unexposed SS316. A dual reflectance-emission setup was partially demonstrated for the purpose of molten fluoride absorption measurement. Preliminary infrared edge data for FLiNaK was obtained. A visible near-infrared (NIR) measurement setup has been constructed and will be used to measure absorption of corrosion products in FLiBe. A new multi-band thermal radiation method using blackbody-weighted spectral bands for the absorption-coefficient was deployed to investigate the thermal performance of the UW Natural Circulation FLiBe Loop under various power levels and operating conditions. The results showed an increasing and decreasing trend in the overall heat transfer with optical thickness with a maximum at optical thickness τ D ≈Ο(1). A forced convection FLiBe loop was constructed for the purpose of studying corrosion and thermal hydraulic phenomena. A vertical cantilever pump is used to achieve up to 5 gpm and 5 m/s flow velocity in hot and cold test sections. A methodology has been established for modeling radiative heat transfer (RHT) in molten salt coolants for the purpose of quantifying the contribution of thermal radiation in conditions relevant to FHR or MSR systems. The comparison between the results obtained for forced and mixed convection showed that in the studied ranges of optical thickness, the maximum effect of RHT on overall heat transfer was 60% improvement observed at z/D = 250 for forced convection; whereas for mixed convection, that value was higher at 170% observed at z/D = 30. The highest effect of RHT on the overall energy distribution was observed at τ D ≈3.4. Negligible effect (<20%) was observed above τ D =18 and below τ D ≈0.002. A CFD code-to-code comparison was conducted using COMSOL Multiphysics at UW-Madison and STAR-CCM at MIT, for a forced and mixed convection laminar flow for a grey medium in a vertical heated using both COMSOL Multiphysics at UW-Madison and STAR-CCM at MIT. The relative error was the largest at τ D =40 for the radiative heat flux, bulk temperature, and Nusselt number. On the other hand, the relative error was the least at τ D =0.002 for the conductive heat flux, bulk temperature, and Nusselt number.

Anderson, Mark↗

Radiation from advanced solid rocket motor plumes

The overall objective of this study was to develop an understanding of solid rocket motor (SRM) plumes in sufficient detail to accurately explain the majority of plume radiation test data. Improved flowfield and radiation analysis codes were developed to accurately and efficiently account for all the factors which effect radiation heating from rocket plumes. These codes were verified by comparing predicted plume behavior with measured NASA/MSFC ASRM test data. Upon conducting a thorough review of the current state-of-the-art of SRM plume flowfield and radiation prediction methodology and the pertinent data base, the following analyses were developed for future design use. The NOZZRAD code was developed for preliminary base heating design and Al2O3 particle optical property data evaluation using a generalized two-flux solution to the radiative transfer equation. The IDARAD code was developed for rapid evaluation of plume radiation effects using the spherical harmonics method of differential approximation to the radiative transfer equation. The FDNS CFD code with fully coupled Euler-Lagrange particle tracking was validated by comparison to predictions made with the industry standard RAMP code for SRM nozzle flowfield analysis. The FDNS code provides the ability to analyze not only rocket nozzle flow, but also axisymmetric and three-dimensional plume flowfields with state-of-the-art CFD methodology. Procedures for conducting meaningful thermo-vision camera studies were developed.

Farmer, Richard C.↗

A general panel sizing computer code and its application to composite structural panels

A computer code for obtaining the dimensions of optimum (least mass) stiffened composite structural panels is described. The procedure, which is based on nonlinear mathematical programming and a rigorous buckling analysis, is applicable to general cross sections under general loading conditions causing buckling. A simplified method of accounting for bow-type imperfections is also included. Design studies in the form of structural efficiency charts for axial compression loading are made with the code for blade and hat stiffened panels. The effects on panel mass of imperfections, material strength limitations, and panel stiffness requirements are also examined. Comparisons with previously published experimental data show that accounting for imperfections improves correlation between theory and experiment.

Anderson, M. S.↗

Early flight test experience with Cockpit Displayed Traffic Information (CDTI)

Coded symbology, based on the results of early human factors studies, was displayed on the electronic horizontal situation indicator and flight tested on an advanced research aircraft in order to subject the coded traffic symbology to a realistic flight environment and to assess its value by means of a direct comparison with simple, uncoded traffic symbology. The tests consisted of 28 curved, decelerating approaches, flown by research-pilot flight crews. The traffic scenarios involved both conflict-free and blunder situations. Subjective pilot commentary was obtained through the use of a questionnaire and extensive pilot debriefing sessions. The results of these debriefing sessions group conveniently under either of two categories: display factors or task performance. A major item under the display factor category was the problem of display clutter. The primary contributors to clutter were the use of large map-scale factors, the use of traffic data blocks, and the presentation of more than a few aircraft. In terms of task performance, the cockpit displayed traffic information was found to provide excellent overall situation awareness.

Abbott, T. S.↗

SATCOM antenna siting study on a P-3C using the NEC-BSC V3.1

The location of a UHF SATCOM antenna on a P-3C aircraft is studied using the NEC-Basic Scattering Code V3.1 (NEC-BSC3). The NEC-BSC3 is a computer code based on the uniform theory of diffraction. The code is first validated for this application using scale model measurements. In general, the comparisons are good except in 10 degree regions near the nose and tail of the aircraft. Patterns for various antenna locations are analyzed to achieve a prescripted performance.

Bensman, D.↗

Experimental and numerical study of swept ramp injection into a supersonic flowfield

Time-averaged measurements of pressure, temperature, velocity, and injectant mole fraction are presented using the planar laser-induced iodine fluorescence (PLIIF) techniques in the complex 3D compressible flowfield around a swept ramp fuel injector. Within the range of thermodynamic conditions present in the test case studied, the technique's accuracy is estimated to be 4 percent for pressure, 4 percent for temperature, 4 percent for velocity, and 3 percent for injectant mole fraction. Comparisons to numerical simulations using the SPARK 3D Navier-Stokes computer code with an algebraic turbulence model are made at the centerplane of the flowfield as well as on three crossflow planes downstream of the injector, and good agreement is found. A weak asymmetry in the incoming flowfield appears to be amplified by boundary layer separation occuring when the the ramp generated shock reflects off the tunnel walls. In the near field of the injector, laminar calculations show better agreement to the measurements than turbulent calculations.

Donohue, James M.↗

Numerical Analysis of the Trailblazer Inlet Flowfield for Hypersonic Mach Numbers

A study of the Trailblazer vehicle inlet was conducted using the Global Air Sampling Program (GASP) code for flight Mach numbers ranging from 4-12. Both perfect gas and finite rate chemical analysis were performed with the intention of making detailed comparisons between the two results. Inlet performance was assessed using total pressure recovery and kinetic energy efficiency. These assessments were based upon a one-dimensional stream-thrust-average of the axisymmetric flowfield. Flow visualization utilized to examine the detailed shock structures internal to this mixed-compression inlet. Kinetic energy efficiency appeared to be the least sensitive to differences between the perfect gas and finite rate chemistry results. Total pressure recovery appeared to be the most sensitive discriminator between the perfect gas and finite rate chemistry results for flight Mach numbers above Mach 6. Adiabatic wall temperature was consistently overpredicted by the perfect gas model for flight Mach numbers above Mach 4. The predicted shock structures were noticeably different for Mach numbers from 6-12. At Mach 4, the perfect gas and finite rate chemistry models collapse to the same result.

Steffen, C. J., Jr.↗

Software for Collaborative Engineering of Launch Rockets

The Rocket Evaluation and Cost Integration for Propulsion and Engineering software enables collaborative computing with automated exchange of information in the design and analysis of launch rockets and other complex systems. RECIPE can interact with and incorporate a variety of programs, including legacy codes, that model aspects of a system from the perspectives of different technological disciplines (e.g., aerodynamics, structures, propulsion, trajectory, aeroheating, controls, and operations) and that are used by different engineers on different computers running different operating systems. RECIPE consists mainly of (1) ISCRM a file-transfer subprogram that makes it possible for legacy codes executed in their original operating systems on their original computers to exchange data and (2) CONES an easy-to-use filewrapper subprogram that enables the integration of legacy codes. RECIPE provides a tightly integrated conceptual framework that emphasizes connectivity among the programs used by the collaborators, linking these programs in a manner that provides some configuration control while facilitating collaborative engineering tradeoff studies, including design to cost studies. In comparison with prior collaborative-engineering schemes, one based on the use of RECIPE enables fewer engineers to do more in less time.

Stanley, Thomas Troy↗

Modeling and Prediction of the Noise from Non-Axisymmetric Jets

The new source model was combined with the original sound propagation model developed for rectangular jets to produce a new version of the rectangular jet noise prediction code. This code was validated using a set of rectangular nozzles whose geometries were specified by NASA. Nozzles of aspect ratios two, four and eight were studied at jet exit Mach numbers of 0.5, 0.7 and 0.9, for a total of nine cases. Reynolds-averaged Navier-Stokes solutions for these jets were provided to the contactor for use as input to the code. Quantitative comparisons of the predicted azimuthal and polar directivity of the acoustic spectrum were made with experimental data provided by NASA. The results of these comparisons, along with a documentation of the propagation and source models, were reported in a journal article publication (Ref. 4). The complete set of computer codes and computational modules that make up the prediction scheme, along with a user's guide describing their use and example test cases, was provided to NASA as a deliverable of this task. The use of conformal mapping, along with simplified modeling of the mean flow field, for noise propagation modeling was explored for other nozzle geometries, to support the task milestone of developing methods which are applicable to other geometries and flow conditions of interest to NASA. A model to represent twin round jets using this approach was formulated and implemented. A general approach to solving the equations governing sound propagation in a locally parallel nonaxisymmetric jet was developed and implemented, in aid of the tasks and milestones charged with selecting more exact numerical methods for modeling sound propagation, and developing methods that have application to other nozzle geometries. The method is based on expansion of both the mean-flowdependent coefficients in the governing equation and the Green's function in series of orthogonal functions. The method was coded and tested on two analytically prescribed mean flows which were meant to represent noise reduction concepts being considered by NASA. Testing (Ref. 5) showed that the method was feasible for the types of mean flows of interest in jet noise applications. Subsequently, this method was further developed to allow use of mean flow profiles obtained from a Reynolds-averaged Navier-Stokes (RANS) solution of the flow. Preliminary testing of the generalized code was among the last tasks completed under this contract. The stringent noise-reduction goals of NASA's Fundamental Aeronautics Program suggest that, in addition to potentially complex exhaust nozzle geometries, next generation aircraft will also involve tighter integration of the engine with the airframe. Therefore, noise generated and propagated by jet flows in the vicinity of solid surfaces is expected to be quite significant, and reduced-order noise prediction tools will be needed that can deal with such geometries. One important source of noise is that generated by the interaction of a turbulent jet with the edge of a solid surface (edge noise). Such noise is generated, for example, by the passing of the engine exhaust over a shielding surface, such as a wing. Work under this task supported an effort to develop a RANS-based prediction code for edge noise based on an extension of the classical Rapid Distortion Theory (RDT) to transversely sheared base flows (Refs. 6 and 7). The RDT-based theoretical analysis was applied to the generic problem of a turbulent jet interacting with the trailing edge of a flat plate. A code was written to evaluate the formula derived for the spectrum of the noise produced by this interaction and results were compared with data taken at NASA Glenn for a variety of jet/plate configurations and flow conditions (Ref. 8). A longer-term goal of this task was to work toward the development of a high-fidelity model of sound propagation in spatially developing non-axisymmetric jets using direct numerical methods for solving the relevant equations. Working with NASA Glenn Acoustics Branch personnel, numerical methods and boundary conditions appropriate for use in a high-resolution calculation of the full equations governing sound propagation in a steady base flow were identified. Computer codes were then written (by NASA) and tested (by OAI) for an increasingly complex set of flow conditions to validate the methods. The NASA-supplied codes were ported to the High-End Computing resources of the NASA Advanced Supercomputing facility for testing and validation against analytical (where possible) and independent numerical solutions. The cases which were completed during the course of this contract were solutions of the two-dimensional linearized Euler equations with no mean flow, a uniform mean flow and a nonuniform mean flow representative of a parallel flow jet.

Leib, Stewart J.↗

Multi-Fidelity CFD Simulation of the Multirotor Test Bed

The high-fidelity and reduced order CFD method's validation and analysis are presented for the Multirotor Test Bed. The reduced-order method uses a rotor-disk model in the OVERFLOW CFD code and is an important aspect of this study. A rotor-trim option loosely coupling the Comprehensive code CAMRAD II with the OVERFLOW CFD code is included in the simulation that captures rotor blade aerodynamics, blade and rotor performance, rotor-rotor, and rotor-body interactions. Detailed comparison and validation were performed for both higher-fidelity full rotor unsteady flows and reduced-order rotor-disk methods. The tunnel geometry is included in the detailed computational effort. Computations of the low-Reynolds number flows were performed for a different number of in-plane rotors as well as for a single rotor modeled inside the tunnel.

ARMD↗

Anisotropic angular scattering models of elastic electron-neutral collisions for Monte Carlo plasma simulations

Many laboratory and industrial plasma applications require accurate modeling techniques to understand the interplay between microscopic and macroscopic processes. A prime example of this interplay is how particle and Monte Carlo (MC) simulation codes describe angular scattering of electrons following elastic scattering events. The forward peaked nature of high energy electron elastic scattering is relatively trivial to accurately describe in plasma simulations. However, for lower energy collisions, which produce near isotropic or backward peaked differential cross sections, there is not a strong consensus among the plasma modeling community on how to best describe these angular scattering trends. Here, in this study, we propose a systematic method to approximate the aforementioned non-trivial angular scattering behavior with a formula that can be readily implemented in particle-in-cell (PIC) and/or MC plasma simulation codes. The present approach is specifically applied to fusion relevant atomic hydrogen and helium, as well as for molecular hydrogen, and results are also applicable to the atomic isotopes and homonuclear molecular isotopologues of these species. Comparisons between the present angular distribution function and benchmark scattering data were used to validate the proposed models. In addition, two-term Boltzmann calculations and PIC direct simulation MC simulations revealed that the proposed angular distribution function is accurate, agreeing very well with benchmark convergent close-coupling scattering calculations, and electron transport measurements. These studies confirmed that the present angular distribution function model can be utilized without the need of renormalization to the momentum transfer cross section (as opposed to using the elastic scattering integrated cross section), which has been suggested by several studies in order to correct for deficient angular scattering models, and to agree with transport measurements. Hence, the present anisotropic angular scattering model can be utilized to accurately model the momentum transfer as well as the electron trajectories of elastic collisions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Using NASA's Worldview to Convey Scientific Finding

NASA Worldview is an interactive interface for browsing full-resolution, global satellite imagery. The application can show single or multiple layers of satellite imagery with over 700 layers available. Each layer has options such as opacity and color palettes which can be adjusted to differentiate layers or underlying data values. The timeline feature allows end-users to study historical events: changes over time can be visualized using the animation tool and exported to animated GIFs. Users can even compare two separate days using the A|B comparison feature. Every time the state of the application changes, it is saved in the URL making it easy to share findings with others.The application is accessible online for anyone to use at any time and the source code is available on GitHub.com. The Worldview source code was designed to be customizable allowing end users to turn on/off features and brand the application for their specific use case. Worldview has code contributors and end-users from all over the world who use the official tool and create personal instances to test and convey their own data.From a team of scientists studying weather over the Philippines, to a group in France studying high environmental risks, Worldview's features and portability allows end-users all over the World to gain valuable earth science insights in a multitude of ways. This presentation will emphasis how end-users are using Worldview to convey their findings.

Rice, Zachary↗

CFD and Experimental Data Comparisons for Conventional and AFC-Enabled CRM High-Lift Configurations

Numerical simulations have been performed for a simplified high-lift (SHL) version of the Common Research Model (CRM) configuration, where the Fowler flaps of a representative conventional high-lift (CRM-HL) configuration are replaced by a set of simple hinged flaps. These hinged flaps are equipped with integrated modular active flow control (AFC) cartridges on the suction surface of the flap shoulder, and the resulting geometry is known as the CRM-SHL-AFC configuration. The main objective is to make use of AFC devices on the CRM-SHL-Configuration to produce the aerodynamic performance (lift) comparable to that of the CRM-HL configuration over a large angle of attack range encompassing maximum lift conditions. For comparison purposes, computations are also per-formed for the CRM-HL configuration. In the current paper, PowerFLOWR©, a CFD code based on the Lattice Boltzmann method (LBM), is used to simulate the entire flow field associated with the CRM-SHL-AFC configuration equipped with a combination of two different types of AFC actuators. The transonic version of the PowerFLOWR© code that has been validated for high-speed flows is used to simulate the flow field generated by the high-momentum actuators required to mitigate reversed flow regions on the suction surfaces of the main wing and the flap. This study is focused on the AFC systems and actuator arrangements that had emerged based on the parametric studies conducted previously at the nominal landing condition. Comparisons of the numerical solutions for lift and surface pressures are presented here with the experimental data, demonstrating the usefulness of CFD for predicting the flow field and lift characteristics of AFC-enabled high-lift configurations over a broad angle of attack range. The numerical solutions predict the expected trends in aerodynamic forces with angle of attack variation.

Veer N Vatsa↗

Status of QCD precision predictions for Drell–Yan rapidity distributions

We compute differential distributions for Drell–Yan processes at the LHC and the Tevatron colliders at next-to-next-to-leading order in perturbative QCD, including fiducial cuts on the decay leptons in the final state. The comparison of predictions obtained with four different codes shows excellent agreement, once linear power corrections from the fiducial cuts are included in those codes that rely on phase-space slicing subtraction schemes. For Z-boson production we perform a detailed study of the symmetric cuts on the transverse momenta of the decay leptons. Predictions at fixed order in perturbative QCD for those symmetric cuts, typically imposed in experiments, suffer from an instability. We show how this can be remedied by an all-order resummation of the fiducial transverse momentum spectrum, and we comment on the choice of cuts for future experimental analyses.

Alekhin, S. [University of Hamburg (Germany)]↗

Calculated viscous and scale effects on transonic aeroelasticity

A viscous-inviscid interactive coupling method is used for the computation of unsteady transonic flows. A lagextrainment integral boundary layer method is used with a transonic small disturbance potential code to compute the transonic aeroelastic response for two wing flutter models. By varying the modeled length scale, viscous effects may be studied as the Reynolds number per reference chordlength varies. Appropriate variation of modeled frequencies and generalized masses then allows comparison of responses for varying scales or Reynolds number. Two wing planforms are studied: one a four percent thick swept wing and the other a typical business jet wing. Calculations for both wings show limit cycle oscillations at transonic speeds in the vicinity of minimum flutter speed indices.

Edwards, John W.↗