SEARCH · Engineering Papers

Results for “Code of Accounts”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 631 records · Page 35

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Boundary layer simulator (BLIMPJ) improvement

Improvements were made to advance the current Boundary Layer Integral Matrix Procedure - Version J (BLIMPJ) containing previously modeled simplified calculation methods accounting for surface roughness, relaminarization, and thick boundary layer effects. These improvements have potential applications in the design of the future Orbit Transfer Vehicles (OTV) engines. The coded wall roughness module was validated successfully against measured data from the LeRC 60-15 degree half-angle conical nozzle. The relaminarization model was modified to include wall cooling and roughness based on a set of test data from the same nozzle. The model was compared against typical data sets for validation in different flow regimes. The thrust loss calculation for thick boundary layer effects was computerized where the iterations between TDK and BLIMPJ proceeded automatically until convergence was achieved. An example is given for a typical OTV nozzle with an area ratio 1293.

Praharaj, S. C.↗

NASA NTRS1989

Local stiffener and skin pocket buckling prediction by special PASCO modeling technique: Correlation to test data

Waffle panels are often used on fuselage structures such as that of the Space Shuttle. The waffle panel design is an efficient design for carrying biaxial, in-plane, and shear loads. The WAFFLE program was designed for application on waffle panels. The Panel Analysis and Sizing Code (PASCO) program was designed for analyzing and sizing uniaxially stiffened panels. The application of the PASCO program in conjunction with the WAFFLE program is discussed to account for both the fillet radius and the presence of stiffness in both directions. The results of the tests are used to verify that these adjustments are valid and necessary if accurate analysis of the waffle panel is to be achieved.

Yin, Dah N.↗

NASA NTRS1996

Prediction Capability for Transonic Nozzle Afterbodies Improved

Under the High Speed Research (HSR) Program, the NASA Lewis Research Center, the NASA Langley Research Center, and the aerospace industry have been developing exhaust nozzle concepts for a future High Speed Civil Transport (HSCT). For such an aircraft to be both economically viable and environmentally acceptable, the exhaust nozzles must combine highly efficient operation throughout the flight envelope with low noise levels at takeoff. A team composed of researchers from Lewis, Langley, and McDonnell Douglas created a two-dimensional data base using computational fluid dynamics. First, each group validated their computational fluid dynamics code against existing experimental data to gain confidence in the code's results. Next, the organizations analyzed several configurations, and the results were used to create the new data base. Then, the analysis methods were updated to account for the effects of two-dimensional nozzles. This will result in a more accurate analysis of the propulsion system for the High Speed Civil Transport.

Source record↗

NASA NTRSDate not supplied

Numerical Simulation of Lean Blowout of Alternative Fuels in 7-element Lean Direct Injector

This research presents the result of numerically simulating 7-element swirl-venturi Lean Direct Injector (SV-LDI) lean blowout (LBO) experiments conducted at NASA Glenn Research Center in May of 2019. After simulating a cold flow case to confirm the pressure drop agrees well with the experiment, additional cases with two different fuels (an average jet fuel and a Gevo alcohol-to-jet fuel) from the National Jet Fuels Combustion Program (NJFCP) were computed to numerically determine the LBO condition. The procedure to approach the LBO follows the method used in the experiment where the air mass flow rate is gradually increased while the fuel supply is maintained. Transient history of global heat release rate as a function of air flow rate is presented, as well as temperature contours at different conditions to give a visual representation of the flame state. The Open National Combustion Code (OpenNCC) used in this research adopted reduced HyChem (Hybrid Chemistry) models along with k-LES turbulence model and a Lagrangian spray model that takes into account droplet internal temperature distribution affected by the shear force on the droplet surface. The transport equations and chemical reaction terms are integrated together to enhance conservation of chemical species that are especially important in the near LBO conditions. After showing the computed range of LBO agrees well with the experimental measurements, time averaged solutions of both fuels at their initial condition and at their limiting condition just before LBO are compared in detail to facilitate the understanding of LBO mechanism.

combustion↗

NASA NTRS2023

Partition Thickness Considerations for Additively Manufactured Acoustic Liners

Three types of uniform-depth liners are evaluated to explore the effects of partition thickness on the surface impedance achieved with additively manufactured liners. A transmission line code is used to predict the effects of sound transmission through empty chambers and wire mesh facesheets, and is combined with the Motsinger and Kraft model to account for sound transmission through perforated facesheets. The inclusion of partitions causes a blockage effect, i.e., a portion of the surface is ‘blocked’ (nonactive surface). For liners with no facesheet, these blockage effects are incorporated simply by accounting for the change in cross sectional area between the individual chamber and the unit cell that includes half of a partition thickness around the perimeter of this chamber. Comparison of impedances predicted in this manner with data acquired in the NASA Langley Normal Incidence Tube (NIT) confirms the efficacy of this modeling approach. Two approaches are considered to account for the inclusion of a wiremesh or perforated plate facesheet on these same cores. The first accounts for the effects of partition thickness before including the transfer impedance across the facesheet, while the second reverses these steps. A comparison of data acquired in the NASA NIT with modeled impedances suggests the first approach is best when the facesheet is a wire mesh, but the second approach is preferred when a perforated plate facesheet is used. It is hypothesized that this is due to the fact that a lumped element model is used to compute the transfer impedance across a wire mesh, while the corresponding transfer impedance across a perforated sheet explicitly incorporates the efffects of the liner core. However, comparisons of data acquired in the NASA Langley Grazing Flow Impedance Tube with modeled impedances seem slightly better when the first approach is employed with a perforated plate facesheet. Thus, the effects of partition thickness on liners with perforated sheets subjected to grazing incidence sound require further review

acoustic↗

NASA NTRS1997

Coupling of Sph and Finite Element Codes for Multi-Layer Orbital Debris Shield Design

Particle-based hydrodynamics models offer distinct advantages over Eulerian and Lagrangian hydrocodes in particular shock physics applications. Particle models are designed to avoid the mesh distortion and state variable diffusion problems which can hinder the effective use of Lagrangian and Eulerian codes respectively. However conventional particle-in-cell and smooth particle hydrodynamics methods employ particles which are actually moving interpolation points. A new particle-based modeling methodology, termed Hamiltonian particle hydrodynamics, was developed by Fahrenthold and Koo (1997) to provide an alternative, fully Lagrangian, energy-based approach to shock physics simulations. This alternative formulation avoids the tensile and boundary instabilities associated with standard smooth particle hydrodynamics formulations and the diffusive grid- to-particle mapping schemes characteristic of particle-in-cell methods. In the work described herein, the method of Fahrenthold and Koo has been extended, by coupling the aforementioned hydrodynamic particle model to a hexahedral finite element based description of the continuum dynamics. The resulting continuum model retains all of the features (including general contact-impact effects) of Hamiltonian particle hydrodynamics, while in addition accounting for tensile strength, plasticity, and damage effects important in the simulation of hypervelocity impact on orbital debris shielding. A three dimensional, vectorized, and autotasked implementation of the extended particle method described here has been coded for application to orbital debris shielding design. Source code for the pre-processor (PREP), analysis code (EXOS), post-processor (POST), and rezoner (ZONE), have been delivered separately, along with a User's Guide describing installation and application of the software.

Fahrenthold, Eric P.↗

DOE OSTI2019

Core polarization effects: Oscillator strengths, transition probabilities and radiative lifetimes of levels in Bi IV

In this work a systematic analysis of Bi IV atomic properties is presented. Ab initio relativistic Hartree–Fock calculations in an intermediate coupling (IC) formalism and taking into account Core Polarization Effects (CPE) are used. We use for the IC calculations the standard method of least square fitting of experimental energy levels by means of computer codes from Cowan which was modified in order to include CPE. Our calculations have been made from 17 configurations of Bi IV: 5d{sup 10}6s{sup 2}, 5d{sup 10}6p{sup 2}, 5d{sup 10}6sns (n= 7–10) and 5d{sup 10}6snd (n= 6–9) for even parity and 5d{sup 10}6snp (n=6,7), 5d{sup 10}6snf (n=5,6), 5d{sup 9}6s{sup 2}6p and 5d{sup 9}6s{sup 2}nf (n=5,6) for odd parity. Transition probabilities and oscillator strengths for 100 lines and some radiative lifetimes of Bi IV have been calculated. Besides a comparison of our results with the scarce experimental values available in the bibliography are presented. The presence of spectral lines of bismuth in stellar atmospheres has been reported in different stars. This is the motivation by this work can be very interest in astrophysical area. Recently the Bi IV is taking importance because it is considered as the main dopant in the prospective used for the design and synthesis of new scintillators.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2019

Spatially dependent embedded self-shielding method for nonuniform temperature distribution

The embedded self-shielding method (ESSM) is simple and robust, and applicable to any geometrical configuration to account for resonance self-shielding by performing fixed source transport calculations. ESSM is very similar to the subgroup method, which has been successful in the MPACT neutronics code developed by CASL. The subgroup method introduces errors in a least-squares fitting process that is used to generate subgroup quadrature sets and cannot be applied to thermal energy groups, as ESSM introduces error for spatially dependent self-shielded cross sections inside the fuel. This study focused on integrating the conventional ESSM with an improved Stoker-Weiss method to obtain spatially dependent resonance self-shielded cross sections for a nonuniform fuel temperature distribution. Here, the proposed method would serve as a substitute for the subgroup method in high fidelity simulations using CASL’s VERA. Benchmark results show that this new, spatially dependent ESSM (SDESSM) is an appropriate candidate to replace the subgroup method.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Probabilistic effects of porosity and chemical kinetics on the shock initiation of an octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) based explosive

In this work, we investigate the effects of porosity and chemical reactivity on the shock initiation of an HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine)-based polymer bonded explosive due to thin flyer impact using two-dimensional mesoscale simulations. Microstructural samples are computationally generated, and the calculated shock response is compared to experiments. The sample dimensions are 0.1×1.05 mm 2 , and each sample accounts for grain scale heterogeneities, including randomized grain morphology and pore spacing. Simulations are performed using the Lawrence Livermore National Laboratory (LLNL) arbitrary Lagrangian–Eulerian multi-physics code ALE3D. Chemistry and equation of states are handled with LLNL's thermochemical code Cheetah. We employ a modified Arrhenius-based chemical kinetics model that scales the temperature exponent fitting parameter with flyer shock conditions. This new model is calibrated to the mean response of experimental shock initiation thresholds for three separate Kapton flyer thicknesses. A simplified shock initiation criterion is proposed, which classifies the macroscale initiation response based on cross-sectional temperature profiles of the shocked sample. This criterion is used to further investigate the probabilistic response of the sample to variations in the mean porosity and chemical kinetics rate, and we find an inverse correlation between the shock pressure and the confidence interval of the experimental threshold velocity. Furthermore, the methodology presented in this work for quantifying the probabilistic nature of shock initiation is material agnostic. These mesoscale studies may be used to inform continuum scale models and expedite the development of future, tailored energetic formulations.

36 MATERIALS SCIENCE↗

DOE OSTI2022

RAFT (Response Amplitudes of Floating Turbines) [SWR-22-40]

RAFT is a Python code for frequency-domain analysis of floating wind turbines. RAFT works in the frequency domain to provide efficient computations of a floating wind system's response spectra accounting for platform hydrodynamics, quasi-static mooring reactions, rotor aerodynamics, and linearized control. It can be used to calculate response amplitude operators (RAO's), power spectral densities, mean/static properties, and response metrics based on mean and root-mean-square values of each output. RAFT serves as the Level 1 model in the WEIS (Wind Energy with Integrated Servo-control) toolkit, and features an OpenMDAO wrapper.

Hall, Matthew↗

NASA NTRS1990

Code For Analysis Of Wing-And-Flap Systems

SUBAERF2 program developed to provide for analysis of subsonic aerodynamics and design by iteration of low-speed wing/flap systems. Based on linearized-theory lifting-surface solution but also accounts for some nonlinear characteristics. Effects of leading- and trailing-edge flaps included. New and improved version of programs described in LAR-13116 and LAR-12987, and replaces both of them. Written in FORTRAN V.

Darden, Christine M.↗

NASA NTRSDate not supplied

Modeling Stellar Jitter and its Effects on the Stellar Spectrum for More Precise Radial Velocity Measurements

Detection of Earth-size exoplanets is a technological and data analysis challenge. Future progress in Earth-mass exoplanet detection is expected from the development of extreme precision radial velocity (EPRV) measurements. The required precision of these measurements necessitates developing a new physics-based data analysis methodology to discriminate planetary signals from host-star-related effects, taking stellar variability and instrumental uncertainties into account. In this work, we investigate and quantify stellar disturbances of the planet-hosting solar-type star HD121504 from 3D radiative modeling obtained with the StellarBox code. The model was used to determine the statistical properties of the turbulent plasma and to obtain synthetic spectroscopic observations for several Fe I lines at different locations on the stellar disk to mimic high-resolution spectroscopic observations. Presented are the initial results of our effort to characterize stellar jitter properties as well as future objectives.

SMD↗

NASA NTRSDate not supplied

Coupling CFD and Material Response for Analysis of Mars Entry

In computing the response of an ablating thermal protection system during atmospheric entry, the aerothermal environment and material response are generally computed separately with a blowing correction term in the material response model to account for the blowing of char and pyrolysis gases. In this work, we apply a coupled approach in which pyrolysis blowing gases, computed in the PATO material response code, are used with a blowing boundary condition in the DPLR hypersonic CFD code. This leads to an iterative method in which blowing products from PATO are input into DPLR to update surface heating estimates. The full iterative method, with the addition of radiative heating estimates using the NEQAIR radiation solver, is shown in Fig. 1. The method is demonstrated on a sphere case with the environment and material properties based on the Mars Science Laboratory entry. Future work includes utilizing this method in computing full 3D material response during the Mars 2020 entry.

Ablation↗

NASA NTRS1999

Computation of Supersonic Jet Mixing Noise Using PARC Code With a kappa-epsilon Turbulence Model

A number of modifications have been proposed in order to improve the jet noise prediction capabilities of the MGB code. This code which was developed at General Electric, employees the concept of acoustic analogy for the prediction of turbulent mixing noise. The source convection and also refraction of sound due to the shrouding effect of the mean flow are accounted for by incorporating the high frequency solution to Lilley's equation for cylindrical jets (Balsa and Mani). The broadband shock-associated noise is estimated using Harper-Bourne and Fisher's shock noise theory. The proposed modifications are aimed at improving the aerodynamic predictions (source/spectrum computations) and allowing for the non- axisymmetric effects in the jet plume and nozzle geometry (sound/flow interaction). In addition, recent advances in shock noise prediction as proposed by Tam can be employed to predict the shock-associated noise as an addition to the jet mixing noise when the flow is not perfectly expanded. Here we concentrate on the aerodynamic predictions using the PARC code with a k-E turbulence model and the ensuing turbulent mixing noise. The geometry under consideration is an axisymmetric convergent-divergent nozzle at its design operating conditions. Aerodynamic and acoustic computations are compared with data as well as predictions due to the original MGB model using Reichardt's aerodynamic theory.

Khavaran, A.↗

NASA NTRS1985

CELCAP: A Computer Model for Cogeneration System Analysis

A description of the CELCAP cogeneration analysis program is presented. A detailed description of the methodology used by the Naval Civil Engineering Laboratory in developing the CELCAP code and the procedures for analyzing cogeneration systems for a given user are given. The four engines modeled in CELCAP are: gas turbine with exhaust heat boiler, diesel engine with waste heat boiler, single automatic-extraction steam turbine, and back-pressure steam turbine. Both the design point and part-load performances are taken into account in the engine models. The load model describes how the hourly electric and steam demand of the user is represented by 24 hourly profiles. The economic model describes how the annual and life-cycle operating costs that include the costs of fuel, purchased electricity, and operation and maintenance of engines and boilers are calculated. The CELCAP code structure and principal functions of the code are described to how the various components of the code are related to each other. Three examples of the application of the CELCAP code are given to illustrate the versatility of the code. The examples shown represent cases of system selection, system modification, and system optimization.

Source record↗

DOE OSTI2024

The influence of cloud cover on the reliability of satellite-based solar resource data

Satellite-based solar resource data are often developed and validated by using binary cloudiness categories: clear sky or overcast cloudy sky. To investigate the reliability of solar resource data in partially cloudy conditions, we estimate cloud fraction using two distinct algorithms: a physical retrieval model using surface observed global horizontal irradiance (GHI) and direct normal irradiance (DNI) and a temporal average of cloud mask data estimated by the observed DNI. Our analysis reveals a significant presence of scattered clouds, broken clouds, and mismatches between satellite- and surface-based cloud data at 17 surface sites across the contiguous United States, though confidently clear and cloudy conditions collectively account for more than 70 % of the data. Solar radiation is computed using the National Solar Radiation Database (NSRDB) algorithm and validated using surface observations. Here, our findings suggest that, in the presence of scattered clouds, NSRDB data for clear-sky conditions can be subject to significant overestimation. In cloudy-sky conditions classified by satellite data, DNI computed by the Fast All-sky Radiation Model for Solar applications with DNI (FARMS-DNI) can be underestimated when limited clouds are detected by surface observations. The bias observed in several cloudiness categories indicates that the NSRDB is exceptionally accurate in confidently clear conditions. However, clear-sky conditions with scattered clouds and mismatched cloud data contribute significantly to the overall uncertainties in the NSRDB. Therefore, future improvements in solar resource data should involve development and implementation of satellite-derived cloud fraction and should consider a novel radiative transfer model accounting for amplified cloud reflection. The evaluation within cloudiness categories also provides a physical rationale for the superior performance of FARMS-DNI compared to the Direct Insolation Simulation Code (DISC) in both cloudy-sky and all-sky conditions.

14 SOLAR ENERGY↗

DOE OSTI2025

High-Fidelity CFD Simulation of Mixed Convection and Forced Convection in a Pebble Bed Test Reactor Core

The Hermes low-power [35-MW(thermal)] reactor will be built and operated by Kairos Power LLC (KP) to demonstrate its fluoride salt-cooled high-temperature reactor (FHR) technology. In the KP FHR, the reactor core is composed of randomly packed pebbles with TRISO fuel particles inside with FLiBe flow upward through the core acting as a coolant. Previous numerical and experimental studies have been limited to either a small-size bed or to a lack of detailed measurements for heat transfer. Here, to address the lack of high-fidelity heat transfer data in a real-size FHR core, in this study, we simulated a pebble bed core with 34 374 pebbles randomly packed, similar to the Hermes reactor's size. The core radius was 14 times that of the pebble diameter, while the core height was 45 times. In this work, we were particularly interested in a mixed convection regime, where buoyancy is important. Therefore, we performed several large-eddy simulations at different Reynolds numbers (160 to 1000) with gravitational force included. The spectral element computational fluid dynamics code NekRS with graphics processing unit acceleration was used for this study. The low-Mach number approximation was applied to address property changes in the FLiBe and to account for buoyancy. A pure hexahedral mesh with 60 million elements was generated by the Voronoi cell method. At the polynomial order of 5, the total degrees of freedom was 7.5 billion. The developed case in this work is the first of its kind in terms of size and complexity. The local numerical data across the domain were obtained and compared with empirical correlations. After examining the data, we found the following conclusions. For pressure drop, the Reger correlation predicted less than a 5% error. On the other hand, for heat transfer, the Wakao correlation outperformed the others. Based on our findings, we recommend the use of the Wakao correlation for the Nusselt number calculation, and for pressure drop, the KTA (Kerntechnischer Ausschuss) correclation, among the available experimental correlations. In conclusion, the Reger direct numerical simulation-driven correlation for pressure drops should also be considered, given its best agreement with our calculations.

Mixed Convection↗

NASA NTRS1989

Flux-split algorithms for the multi-dimensional Euler equations with real gases

Upwind algorithms are developed for the numerical solution of the multidimensional Euler equations for real gases. Flux-splitting methods are derived which account for a general equation of state. Approximations to the state equation based on physical arguments result in simplified algorithms which may be implemented into existing perfect-gas codes. Applications of the method to several high-Mach-number high-temperature flows are presented for two and three space dimensions.

Grossman, B.↗