SEARCH · Engineering Papers

Results for “CALCULATION METHODS”

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 271 records · Page 15

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Impacts of Biofuel Blending on MCCI Ignition Delay with Review of Methods for Defining Cycle-by-Cycle Ignition Points from Noisy Cylinder Pressure Data

Conventional diesel combustion, also known as Mixing-Controlled Compression Ignition (MCCI), is expected to be the primary power source for medium- and heavy-duty vehicles for decades to come. Displacing petroleum-based ultra-low-sulfur diesel (ULSD) as much as possible with low-net-carbon biofuels will become necessary to help mitigate effects on climate change. Neat biofuels may have difficulty meeting current diesel fuel standards but blends of 30% biofuel in ULSD show potential as ‘drop-in’ fuels. These blends must not make significant changes to the combustion phasing of the MCCI process if they are to be used interchangeably with neat ULSD. An important aspect of MCCI phasing is the ignition delay (ID), i.e. the time between the start of fuel injection and the initial premixed autoignition that initiates the MCCI process. Bench experiments can evaluate the expected ignition delay of a fuel via cetane number (CN) or alternative methods such as the indicated cetane number (ICN); however, neither CN nor ICN correlate perfectly with the ignition delay measured in actual engine experiments. Furthermore, there is no standardized methodology on how to quantify MCCI ignition delay from engine cylinder pressure measurements, creating difficulties in cross-study comparison. In this study, several engine ignition delay calculation methods are evaluated for robustness in deriving ignition delay on both a cycle-averaged and cycle-to-cycle basis. Eight biofuel blends with varying ICN, oxygen concentration and other fuel properties were used to study the different methods. This yields a thorough analysis of how certain biofuel blends affect ignition delay and the entire MCCI process, as well as a thorough evaluation of the differences between the ID calculation methods. Many of these methods are equally valid, but the choice of method has a significant impact on the resulting ID, which must be carefully considered when evaluating results across multiple studies.

47 OTHER INSTRUMENTATION↗

DOE OSTI2021

Impacts of Biofuel Blending on MCCI Ignition Delay with Review of Methods for Defining Cycle-by-Cycle Ignition Points from Noisy Cylinder Pressure Data: Preprint

47 OTHER INSTRUMENTATION↗

DOE OSTI2023

A Better Method to Calculate Fuel Burnup in Pebble Bed Reactors Using Machine Learning

Burnup measurement is an important step in material control and accountancy (MC&A) at nuclear reactors, and may be done by examining gamma spectra of fuel samples. Traditional approaches rely on known correlations to specific photopeaks (e.g. 137 Cs) and operate via a standard linear regression method. However, the quality of these regression methods is limited even in the best case, and is significantly poorer at short fuel cool-down times, due to the elevated radiation background by short life-time isotopes, and self-shielding effect of the fuel. For practical operation of pebble bed reactors (PBRs), quick measurements (in minutes) and short cooling times (in hours) are required from a safety and security perspective. We investigated the efficacy and performance of machine learning (ML) methods to predict the burnup of the pebble fuel from full gamma spectra (rather than specific discrete photopeaks) and found a full-spectrum ML approach to far outperform baseline regression predictions in all measurement and cooling conditions - including in operational-like measurement conditions. We also performed model and data ablation experiments to determine the relative performance impact of our ML methods' capacity to model data nonlinearities and the inherent additional information in full spectra. Applying our ML methods, we found a number of surprising results, including improved accuracy at shorter fuel cooling times (the opposite of the norm), remarkable robustness to spectrum compression (via rebinning), and competitive burnup predictions even when using background signal only (i.e. explicitly omitting known isotope photopeaks).

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

NASA NTRS2004

New Reduced Two-Time Step Method for Calculating Combustion and Emission Rates of Jet-A and Methane Fuel With and Without Water Injection

A simplified kinetic scheme for Jet-A, and methane fuels with water injection was developed to be used in numerical combustion codes, such as the National Combustor Code (NCC) or even simple FORTRAN codes that are being developed at Glenn. The two time step method is either an initial time averaged value (step one) or an instantaneous value (step two). The switch is based on the water concentration in moles/cc of 1x10(exp -20). The results presented here results in a correlation that gives the chemical kinetic time as two separate functions. This two step method is used as opposed to a one step time averaged method previously developed to determine the chemical kinetic time with increased accuracy. The first time averaged step is used at the initial times for smaller water concentrations. This gives the average chemical kinetic time as a function of initial overall fuel air ratio, initial water to fuel mass ratio, temperature, and pressure. The second instantaneous step, to be used with higher water concentrations, gives the chemical kinetic time as a function of instantaneous fuel and water mole concentration, pressure and temperature (T4). The simple correlations would then be compared to the turbulent mixing times to determine the limiting properties of the reaction. The NASA Glenn GLSENS kinetics code calculates the reaction rates and rate constants for each species in a kinetic scheme for finite kinetic rates. These reaction rates were then used to calculate the necessary chemical kinetic times. Chemical kinetic time equations for fuel, carbon monoxide and NOx were obtained for Jet-A fuel and methane with and without water injection to water mass loadings of 2/1 water to fuel. A similar correlation was also developed using data from NASA's Chemical Equilibrium Applications (CEA) code to determine the equilibrium concentrations of carbon monoxide and nitrogen oxide as functions of overall equivalence ratio, water to fuel mass ratio, pressure and temperature (T3). The temperature of the gas entering the turbine (T4) was also correlated as a function of the initial combustor temperature (T3), equivalence ratio, water to fuel mass ratio, and pressure.

Molnar, Melissa↗

DOE OSTI2024

An introduction to Spent Nuclear Fuel decay heat for Light Water Reactors: a review from the NEA WPNCS

This paper summarized the efforts performed to understand decay heat estimation from existing spent nuclear fuel (SNF), under the auspices of the Working Party on Nuclear Criticality Safety (WPNCS) of the OECD Nuclear Energy Agency. Needs for precise estimations are related to safety, cost, and optimization of SNF handling, storage, and repository. The physical origins of decay heat (a more correct denomination would be decay power) are then introduced, to identify its main contributors (fission products and actinides) and time-dependent evolution. Due to limited absolute prediction capabilities, experimental information is crucial; measurement facilities and methods are then presented, highlighting both their relevance and our need for maintaining the unique current full-scale facility and developing new ones. The third part of this report is dedicated to the computational aspect of the decay heat estimation: calculation methods, codes, and validation. Different approaches and implementations currently exist for these three aspects, directly impacting our capabilities to predict decay heat and to inform decision-makers. Finally, recommendations from the expert community are proposed, potentially guiding future experimental and computational developments. One of the most important outcomes of this work is the consensus among participants on the need to reduce biases and uncertainties for the estimated SNF decay heat. If it is agreed that uncertainties (being one standard deviation) are on average small (less than a few percent), they still substantially impact various applications when one needs to consider up to three standard deviations, thus covering more than 95% of cases. The second main finding is the need of new decay heat measurements and validation for cases corresponding to more modern fuel characteristics: higher initial enrichment, higher average burnup, as well as shorter and longer cooling time. Similar needs exist for fuel types without public experimental data, such as MOX, VVER, or CANDU fuels. A third outcome is related to SNF assemblies for which no direct validation can be performed, representing the vast majority of cases (due to the large number of SNF assemblies currently stored, or too short or too long cooling periods of interest). A few solutions are possible, depending on the application. For the final repository, systematic measurements of quantities related to decay heat can be performed, such as neutron or gamma emission. This would provide indications of the SNF decay heat at the time of encapsulation. For other applications (short- or long-term cooling), the community would benefit from applying consistent and accepted recommendations on calculation methods, for both decay heat and uncertainties. This would improve the understanding of the results and make comparisons easier.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS1986

Acoustic-Liner Admittance in a Duct

Method calculates admittance from easily obtainable values. New method for calculating acoustic-liner admittance in rectangular duct with grazing flow based on finite-element discretization of acoustic field and reposing of unknown admittance value as linear eigenvalue problem on admittance value. Problem solved by Gaussian elimination. Unlike existing methods, present method extendable to mean flows with two-dimensional boundary layers as well. In presence of shear, results of method compared well with results of Runge-Kutta integration technique.

Watson, W. R.↗

NASA NTRS1989

Calculation of turbulence-driven secondary motion in ducts with arbitrary cross section

Calculation methods for turbulent duct flows are generalized for ducts with arbitrary cross-sections. The irregular physical geometry is transformed into a regular one in computational space, and the flow equations are solved with a finite-volume numerical procedure. The turbulent stresses are calculated with an algebraic stress model derived by simplifying model transport equations for the individual Reynolds stresses. Two variants of such a model are considered. These procedures enable the prediction of both the turbulence-driven secondary flow and the anisotropy of the Reynolds stresses, in contrast to some of the earlier calculation methods. Model predictions are compared to experimental data for developed flow in triangular duct, trapezoidal duct and a rod-bundle geometry. The correct trends are predicted, and the quantitative agreement is mostly fair. The simpler variant of the algebraic stress model procured better agreement with the measured data.

Demuren, A. O.↗

NASA NTRS1990

Calculation of turbulence-driven secondary motion in ducts with arbitrary cross-section

Demuren, A. O.↗

NASA NTRS2005

Simplified Two-Time Step Method for Calculating Combustion Rates and Nitrogen Oxide Emissions for Hydrogen/Air and Hydorgen/Oxygen

A simplified single rate expression for hydrogen combustion and nitrogen oxide production was developed. Detailed kinetics are predicted for the chemical kinetic times using the complete chemical mechanism over the entire operating space. These times are then correlated to the reactor conditions using an exponential fit. Simple first order reaction expressions are then used to find the conversion in the reactor. The method uses a two-time step kinetic scheme. The first time averaged step is used at the initial times with smaller water concentrations. This gives the average chemical kinetic time as a function of initial overall fuel air ratio, temperature, and pressure. The second instantaneous step is used at higher water concentrations (> 1 x 10(exp -20) moles/cc) in the mixture which gives the chemical kinetic time as a function of the instantaneous fuel and water mole concentrations, pressure and temperature (T4). The simple correlations are then compared to the turbulent mixing times to determine the limiting properties of the reaction. The NASA Glenn GLSENS kinetics code calculates the reaction rates and rate constants for each species in a kinetic scheme for finite kinetic rates. These reaction rates are used to calculate the necessary chemical kinetic times. This time is regressed over the complete initial conditions using the Excel regression routine. Chemical kinetic time equations for H2 and NOx are obtained for H2/air fuel and for the H2/O2. A similar correlation is also developed using data from NASA s Chemical Equilibrium Applications (CEA) code to determine the equilibrium temperature (T4) as a function of overall fuel/air ratio, pressure and initial temperature (T3). High values of the regression coefficient R2 are obtained.

Molnar, Melissa↗

DOE OSTI2023

Loop current fluctuations and quantum critical transport

We study electrical transport at quantum critical points (QCPs) associated with loop current ordering in a metal, focusing specifically on models of the “Hertz-Millis” type. At the infrared (IR) fixed point and in the absence of disorder, the simplest such models have infinite DC conductivity and zero incoherent conductivity at nonzero frequencies. However, we find that a particular deformation, involving N N species of bosons and fermions with random couplings in flavor space, admits a finite incoherent, frequency-dependent conductivity at the IR fixed point, \sigma(\omega>0)\sim\omega^{-2/z} σ ( ω > 0 ) ∼ ω − 2 / z , where z z is the boson dynamical exponent. Leveraging the non-perturbative structure of quantum anomalies, we develop a powerful calculational method for transport. The resulting "anomaly-assisted large N N expansion" allows us to extract the conductivity systematically. Although our results imply that such random-flavor models are problematic as a description of the physical N = 1 N = 1 system, they serve to illustrate some general conditions for quantum critical transport as well as the anomaly-assisted calculational methods. In addition, we revisit an old result that irrelevant operators generate a frequency-dependent conductivity, \sigma(\omega>0) \sim \omega^{-2(z-2)/z} σ ( ω > 0 ) ∼ ω − 2 ( z − 2 ) / z , in problems of this kind. We show explicitly, within the scope of the original calculation, that this result does not hold for any order parameter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2005

Summary of Simplified Two Time Step Method for Calculating Combustion Rates and Nitrogen Oxide Emissions for Hydrogen/Air and Hydrogen/Oxygen

A simplified single rate expression for hydrogen combustion and nitrogen oxide production was developed. Detailed kinetics are predicted for the chemical kinetic times using the complete chemical mechanism over the entire operating space. These times are then correlated to the reactor conditions using an exponential fit. Simple first order reaction expressions are then used to find the conversion in the reactor. The method uses a two time step kinetic scheme. The first time averaged step is used at the initial times with smaller water concentrations. This gives the average chemical kinetic time as a function of initial overall fuel air ratio, temperature, and pressure. The second instantaneous step is used at higher water concentrations (greater than l x 10(exp -20)) moles per cc) in the mixture which gives the chemical kinetic time as a function of the instantaneous fuel and water mole concentrations, pressure and temperature (T(sub 4)). The simple correlations are then compared to the turbulent mixing times to determine the limiting properties of the reaction. The NASA Glenn GLSENS kinetics code calculates the reaction rates and rate constants for each species in a kinetic scheme for finite kinetic rates. These reaction rates are used to calculate the necessary chemical kinetic times. This time is regressed over the complete initial conditions using the Excel regression routine. Chemical kinetic time equations for H2 and NOx are obtained for H2/Air fuel and for H2/O2. A similar correlation is also developed using data from NASA's Chemical Equilibrium Applications (CEA) code to determine the equilibrium temperature (T(sub 4)) as a function of overall fuel/air ratio, pressure and initial temperature (T(sub 3)). High values of the regression coefficient R squared are obtained.

Marek, C. John↗

NASA NTRS1987

A method of calculating the total flow from a given sea surface topography

Using a simple dynamical model of a wind-driven ocean circulation of the Stommel type, and an analytical basis developed to objectively analyze the sea surface height residuals from an altimeter and, in the process, to determine the total flow instead of just the near surface geostrophic component associated with the given sea surface topography. The method is based on first deriving the solution to the forced problem for a given wind stress required to develop a hypothetical true or perfect data field and to establishing the basis for the objective analysis. The stream function and the surface height field for the forced problem are developed in terms of certain characteristic functions with the same expansion coefficients for both fields. These characteristic functions are simply the solutions for a homogeneous elliptic equation for the stream function and the solutions of an inhomogeneous balance equation for the height field. For the objective analysis, using a sample of randomly selected height values from the true data field, the height field characteristic functions are used to fit the given topography in a least squares sense. The resulting expansion coefficients then permit the synthesis of the total flow field via the stream function characteristic modes and the solution is perfectly well behaved even along the equator. The method of solution is easily adaptable to realistic ocean basis by straight forward numerical methods. The analytical basis of the theory and the results for an ideal rectangular basin on a beta plane are described.

Rao, Desiraju B.↗

DOE OSTI2020

Impact of boundary layer simulation on predicting radioactive pollutant dispersion: a case study for HANARO Research Reactor using the WRF-MMIF-CALPUFF modeling system

Wind plays an important role in cases of unexpected radioactive pollutant dispersion, deciding distribution and concentration of the leaked substance. The accurate prediction of wind has been challenging in numerical weather prediction models, especially near the surface because of the complex interaction between turbulent flow and topographic effect. As such, in this study, we investigated the characteristics of atmospheric dispersion of radioactive material (i.e. 137 Cs) according to the simulated boundary layer around the HANARO research nuclear reactor in Korea using the Weather Research and Forecasting (WRF)-Mesoscale Model Interface (MMIF)-California Puff (CALPUFF) model system. We examined the impacts of orographic drag on wind field, stability calculation methods, and planetary boundary layer parameterizations on the dispersion of radioactive material under a radioactive leaking scenario. We found that inclusion of the orographic drag effect in the WRF model improved the wind prediction most significantly over the complex terrain area, leading the model system to estimate the radioactive concentration near the reactor more conservatively. We also emphasized the importance of the stability calculation method and employing the skillful boundary layer parameterization to ensure more accurate low atmospheric conditions, in order to simulate more feasible spatial distribution of the radioactive dispersion in leaking scenarios.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1977

Approximate method for calculating free vibrations of a large-wind-turbine tower structure

A set of ordinary differential equations were derived for a simplified structural dynamic lumped-mass model of a typical large-wind-turbine tower structure. Dunkerley's equation was used to arrive at a solution for the fundamental natural frequencies of the tower in bending and torsion. The ERDA-NASA 100-kW wind turbine tower structure was modeled, and the fundamental frequencies were determined by the simplified method described. The approximate fundamental natural frequencies for the tower agree within 18 percent with test data and predictions analyzed.

Das, S. C.↗

NASA NTRSDate not supplied

Relating Land Cover Change to Runoff Distribution Using NASA Earth Observations in Riley County, Kansas

In recent years, Riley County, Kansas experienced high levels of flooding potentially due to changes in land use and impervious surface. To analyze the increase in floods and runoff patterns over time, this study contrasted two methods of curve number (CN) calculation, which estimated runoff from rainfall events, from 2006-2020. The conventional calculation method used a look-up table and tracked land use/land cover change. This method resulted in static tables from various inputs such as soil type, surface imperviousness, slope, and crop phenology to estimate runoff versus infiltration. These tables allowed for phenology specific curve numbers and accounted for various farming techniques but required many inputs and lacked flexibility to include seasonal variations in vegetation. In contrast, the dynamic method is a more modern and simple calculation employing normalized difference vegetation index (NDVI) data from Landsat 5, Landsat 7, and Landsat 8 satellites compiled over the rainy season each year to calculate curve numbers based off seasonal vegetation. This dynamic method allowed for more precise analysis of runoff variability because NDVI is a more accurate, real-time measurement for land cover and it could be calculated with greater temporal detail. The results will allow decision makers to make more informed decisions on resiliency strategies to address future flooding.

Trista Brophy↗

NASA NTRS2019

An Analytical Method to Calculate Effective Elastic Properties of General Multifunctional Honeycomb Cores in Sandwich Composites

Sandwich composite structures are ideal configurations in which to incorporate additional functionality beyond load-carrying capabilities. The inner core-walls can be layered to incorporate other functions such as power storage for a battery. In this work we investigate an assemblage of analytical tools to compute effective properties that allow complex layered core architectures to be homogenized into a single continuum layer. This provides a great increase in computational efficiency to numerically simulate the structural response of multifunctional sandwich structures under applied loads. We present a coupled analytical method including an extensive numerical verification of the accuracy of this method.

Sather, Erik↗

DOE OSTI2022

Centrifuge modeling and numerical analysis on lateral performance of mono-bucket foundation for offshore wind turbines

The lateral loading performance of the mono-bucket foundation for offshore wind turbines is investigated in this study through a series of centrifuge tests and numerical analysis. Six mono-bucket models are considered to assess the influences of aspect ratio under the ultimate condition. The FE model is validated against the centrifuge test results, and the failure mechanism of the foundation is revealed. The parametric study and sensitivity analysis of the soil parameters are performed in the numerical simulations. A modified theoretical calculation method is proposed to estimate the ultimate lateral capacity of the suction bucket foundation by solving the limit equilibrium equations. A correction term is added in the calculation using the “m” method. The position of the rotation center is calculated, demonstrating a better accuracy compared to the previously reported geometric method. The modified calculation method is verified against the field test, the centrifuge test, and the laboratory test. The method is proved to be applicable in estimating the lateral performance of the mono-bucket foundation in a general form. Finally, this study aims to provide well-documented data and design references for practical engineers.

17 WIND ENERGY↗

NASA NTRS1990

Boundary-integral method for calculating aerodynamic sensitivities with illustration for lifting-surface theory

The technique of implicit differentiation has been used in combination with linearized lifting-surface theory to derive analytical expressions for aerodynamic sensitivities (i.e., rates of change of lifting pressures with respect to general changes in aircraft geometry, including planform variations) for steady or oscillating planar or nonplanar lifting surfaces in subsonic, sonic, or supersonic flow. The geometric perturbation is defined in terms of a single variable, and the user need only provide simple expressions or similar means for defining the continuous or discontinuous global or local perturbation of interest. Example expressions are given for perturbations of the sweep, taper, and aspect ratio of a wing with trapezoidal semispan planform. The present process appears to be readily adaptable to more general surface-panel methods.

Yates, E. Carson, Jr.↗