Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equivalence Method”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Practical Methodology for the Inclusion of Nonlinear Slosh Damping in the Stability Analysis of Liquid-Propelled Space Vehicles

A method is presented that allows for the prediction of the magnitude of limit cycles due to adverse control-slosh interaction in liquid propelled space vehicles using non-linear slosh damping. Such a method is an alternative to the industry practice of assuming linear damping and relying on: mechanical slosh baffles to achieve desired stability margins; accepting minimal slosh stability margins; or time domain non-linear analysis to accept time periods of poor stability. Sinusoidal input describing functional analysis is used to develop a relationship between the non-linear slosh damping and an equivalent linear damping at a given slosh amplitude. In addition, a more accurate analytical prediction of the danger zone for slosh mass locations in a vehicle under proportional and derivative attitude control is presented. This method is used in the control-slosh stability analysis of the NASA Space Launch System.

Ottander, John A.↗
NASA NTRS2012

Development of the GPM Observatory Thermal Vacuum Test Model

A software-based thermal modeling process was documented for generating the thermal panel settings necessary to simulate worst-case on-orbit flight environments in an observatory-level thermal vacuum test setup. The method for creating such a thermal model involved four major steps: (1) determining the major thermal zones for test as indicated by the major dissipating components on the spacecraft, then mapping the major heat flows between these components; (2) finding the flight equivalent sink temperatures for these test thermal zones; (3) determining the thermal test ground support equipment (GSE) design and initial thermal panel settings based on the equivalent sink temperatures; and (4) adjusting the panel settings in the test model to match heat flows and temperatures with the flight model. The observatory test thermal model developed from this process allows quick predictions of the performance of the thermal vacuum test design. In this work, the method described above was applied to the Global Precipitation Measurement (GPM) core observatory spacecraft, a joint project between NASA and the Japanese Aerospace Exploration Agency (JAXA) which is currently being integrated at NASA Goddard Space Flight Center for launch in Early 2014. From preliminary results, the thermal test model generated from this process shows that the heat flows and temperatures match fairly well with the flight thermal model, indicating that the test model can simulate fairly accurately the conditions on-orbit. However, further analysis is needed to determine the best test configuration possible to validate the GPM thermal design before the start of environmental testing later this year. Also, while this analysis method has been applied solely to GPM, it should be emphasized that the same process can be applied to any mission to develop an effective test setup and panel settings which accurately simulate on-orbit thermal environments.

Yang, Kan↗
NASA NTRS1959

Investigation of the Effects of Low Reynolds Number Operation on the Performance of a Single-Stage Turbine with a Downstream Stator

High-altitude turbojet performance is adversely affected by the effects of low air density. This performance loss is evaluated as a Reynolds number effect, which represents the increased significance of high fluid viscous forces in relation to dynamic fluid forces as the Reynolds number is decreased. An analytical and experimental investigation of the effects of low Reynolds number operation on a single-stage, high-work-output turbine with a downstream stator was carried out at Reynolds numbers of 182,500, 39,600, and 23,000, based on average rotor-design flow conditions. At low Reynolds numbers and turbulent flow conditions, increased viscous losses caused decreased effective flow area, and thus decreased weight flow, torque, and over-all efficiency at a given equivalent speed and pressure ratio. Decreasing the Reynolds number from 182,500 to 23,000 at design equivalent speed resulted in a 5.00-point loss in peak over-all turbine efficiency for both theory and experiment. The choking equivalent weight flow decreased 2.30 percent for these conditions. Limiting loading work output was reached at design equivalent speed for all three Reynolds numbers. The value of limiting loading work output at design speed decreased 4.00 percent as Reynolds number was decreased from 182,500 to 23,000. A theoretical performance-prediction method using basic boundary-layer relations gave good agreement with experimental results over most of the performance range at a given Reynolds number if the experimental and analytical design operating conditions were carefully matched at the highest Reynolds number with regard to design performance parameters. High viscous losses in the inlet stator and rotor prevented the attainment of design equivalent work output at the lowest Reynolds number of 23,000.

Forrette, Robert E.↗
NASA NTRS1947

Distribution of wave drag and lift in the vicinity of wing tips at supersonic speeds

The point-source-distribution method of calculating the aerodynamic coefficients of thin wings at supersonic speeds was extended to include the effect of the region between the wing boundary and the foremost Mach wave from the wing leading edge. The effect of this region on the surface velocity potential has been determined by an equivalent function, which is evaluated over a portion of the wing surface. In this manner, the effect of angles of attack and yaw as well as the asymmetry of top and bottom wing surfaces may be calculated. As examples of the method, the pressure distribution on a thin plate wing of rectangular plan form as well as the lift and the drag coefficients as a function of Mach number, angle of attack, and aspect ratio are calculated. The equations for the surface velocity potential of several other plan forms are also included. (author)

AIRFOILS, FLAT-PLATE↗
NASA NTRS1984

An objective method for forecasting tropical cyclone motion using Nimbus and NOAA-2 infrared measurements

A statistical method for predicting tropical cyclone positions 12, 24, 48 and 72 hours after initial observation is described. Satellite measurements were at the 11.0 micron channel and at equivalent blackbody temperatures which give representations of the structure of the cyclones, as well as the circulation features surrounding them. In all, 197 individual measurements for blackbody temperature were made. Algorithms prepared from a single satellite image included 14 climatological and persistence type variables, and are based on eigenvector preprocessing and a series of independent screening tests. The independent tests showed that the average error made by the algorithms from the single satellite observation were comparable to the 48 hr Joint Typhoon Warning Center (JTWC) forecast and were about 10 percent better for 72 hr forecasts. Forecasts using only climatological persistence variables were found to be 20 percent worse than JTWC for 24 hr forecasts and 10 percent worse for 48 and 72 hr forecasts. When both satellite and nonsatellite variables were included, it was found that the method was comparable to JTWC for the 24 and 48 hr forecasts and about 25 percent better for the 72 hr forecasts.

Hunter, H. E.↗
NASA NTRS1988

Prediction of high temperature metal matrix composite ply properties

The application of the finite element method (superelement technique) in conjunction with basic concepts from mechanics of materials theory is demonstrated to predict the thermomechanical behavior of high temperature metal matrix composites (HTMMC). The simulated behavior is used as a basis to establish characteristic properties of a unidirectional composite idealized an as equivalent homogeneous material. The ply properties predicted include: thermal properties (thermal conductivities and thermal expansion coefficients) and mechanical properties (moduli and Poisson's ratio). These properties are compared with those predicted by a simplified, analytical composite micromechanics model. The predictive capabilities of the finite element method and the simplified model are illustrated through the simulation of the thermomechanical behavior of a P100-graphite/copper unidirectional composite at room temperature and near matrix melting temperature. The advantage of the finite element analysis approach is its ability to more precisely represent the composite local geometry and hence capture the subtle effects that are dependent on this. The closed form micromechanics model does a good job at representing the average behavior of the constituents to predict composite behavior.

Caruso, J. J.↗
NASA NTRS1976

Numerical Airfoil Optimization Using a Reduced Number of Design Coordinates

A method is presented for numerical airfoil optimization whereby a reduced number of design coordinates are used to define the airfoil shape. The approach is to define the airfoil as a linear combination of shapes. These basic shapes may be analytically or numerically defined, allowing the designer to use his insight to propose candidate designs. The design problem becomes one of determining the participation of each such function in defining the optimum airfoil. Examples are presented for two-dimensional airfoil design and are compared with previous results based on a polynomial representation of the airfoil shape. Four existing NACA airfoils are used as basic shapes. Solutions equivalent to previous results are achieved with a factor of more than 3 improvements in efficiency, while superior designs are demonstrated with an efficiency greater than 2 over previous methods. With this shape definition, the optimization process is shown to exploit the simplifying assumptions in the inviscid aerodynamic analysis used here, thus demonstrating the need to use more advanced aerodynamics for airfoil optimization.

Vanderplaats, G. N.↗
DOE OSTI2021

Fast multiscale contrast independent preconditioners for linear elastic topology optimization problems

The goal of this work is to present a fast and viable approach for the numerical solution of the high-contrast state problems arising in topology optimization. The optimization process is iterative, and the gradients are obtained by an adjoint analysis, which requires the numerical solution of large high-contrast linear elastic problems with features spanning several length scales. The size of the discretized problems forces the utilization of iterative linear solvers with solution time dependent on the quality of the preconditioner. The lack of clear separation between the scales, as well as the high-contrast, imposes severe challenges on the standard preconditioning techniques. Thus, here we propose new methods for the high-contrast elasticity equation with performance independent of the high-contrast and the multi-scale structure of the elasticity problem. The solvers are based on two-levels domain decomposition techniques with a carefully constructed coarse level to deal with the high-contrast and multi-scale nature of the problem. The construction utilizes spectral equivalence between scalar diffusion and each displacement block of the elasticity problems and, in contrast to previous solutions proposed in the literature, is able to select the appropriate dimension of the coarse space automatically. The new methods inherit the advantages of domain decomposition techniques, such as easy parallelization and scalability. Finally, the presented numerical experiments demonstrate the excellent performance of the proposed methods.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

Analysis of microstrip dipoles and slots transversely coupled to a microstrip line using the FDTD method

Printed dipole elements and their complement, linear slots, are elementary radiators that have found use in low-profile antenna arrays. Low-profile antenna arrays, in addition to their small size and low weight characteristics, offer the potential advantage of low-cost, high-volume production with easy integration with active integrated circuit components. The design of such arrays requires that the radiation and impedance characteristics of the radiating elements be known. The FDTD (Finite-Difference Time-Domain) method is a general, straight-forward implementation of Maxwell's equations and offers a relatively simple way of analyzing both printed dipole and slot elements. Investigated in this work is the application of the FDTD method to the analysis of printed dipole and slot elements transversely coupled to an infinite transmission line in a multilayered configuration. Such dipole and slot elements may be used in dipole and slot series-fed-type linear arrays, where element offsets and interelement line lengths are used to obtain the desired amplitude distribution and beam direction, respectively. The design of such arrays is achieved using transmission line theory with equivalent circuit models for the radiating elements. In an equivalent circuit model, the dipole represents a shunt impedance to the transmission line, where the impedance is a function of dipole offset, length, and width. Similarly, the slot represents a series impedance to the transmission line. The FDTD method is applied to single dipole and slot elements transversely coupled to an infinite microstrip line using a fixed rectangular grid with Mur's second order absorbing boundary conditions. Frequency-dependent circuit and scattering parameters are obtained by saving desired time-domain quantities and using the Fourier transform. A Gaussian pulse excitation is applied to the microstrip transmission line, where the resulting reflected signal due to the presence of the radiating element is used to determine the equivalent element impedance.

Tulintseff, A. N.↗
DOE OSTI2022

COG Beta-Effective Benchmarks

COG 11.3 calculations with the “new method” using ENDF/B-VIII.0 data are shown to produce results with |C/E – 1| < 3σ in all cases, whereas the other methods produce a few minor outliers. Furthermore, the “new method” using ENDF/B-VIII.0 data is observed to produce “over” and “under” results consistent with 50/50 expectations whereas the other methods appear to be slightly over or under predicting. JEFF-3.3 results are slightly larger overall than ENDF/B-VIII.0 results. In all but a few cases, it is also observed that $F_d/F_t ≥ (vF)_d/(vF)_t ≥ 1 – (k_{prompt}/k_{eff})$ and all estimators produce equivalent results for the most thermal systems. In those other cases, this relation holds within the given statistics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

A Cooperative Game Theory-based Approach to Under-frequency Load Shedding Control

This paper proposes a cooperative game theory-based under-frequency load shedding (UFLS) approach for frequency stability and control in power systems. UFLS is a crucial factor for frequency stability and control especially in power grids with high penetration of renewable energy sources and restructured power systems. Conventional UFLS methods, most of which are off-line, usually shed fixed amounts of predetermined loads based on a predetermined schedule which can lead to over or under curtailment of load. This paper presents a co-operative game theory-based two-stage strategy to effectively and precisely determine locations and amounts of loads to be shed for UFLS control. In the first stage, the total amount of loads to be shed, also referred to as deficit in generation or the disturbance power, is computed using the initial rate of change of frequency (ROCOF) referred to the equivalent inertial center. In the second stage, the Shapley value, one of the solution concepts of cooperative game theory, is used to determine load shedding amounts and locations. The proposed method is implemented on the reduced 9-bus 3-machine Western Electricity Coordinating Council (WECC) system and simulated on Real-time Digital Simulators (RTDS). The results show that the proposed UFLS approach can effectively return the system to normal state after disturbances.

42 ENGINEERING↗
DOE OSTI2022

On the S-matrix of Ising field theory in two dimensions

We explore the analytic structure of the non-perturbative S-matrix in arguably the simplest family of massive non-integrable quantum field theories: the Ising field theory (IFT) in two dimensions, which may be viewed as the Ising CFT deformed by its two relevant operators, or equivalently, the scaling limit of the Ising model in a magnetic field. Our strategy is that of collider physics: we employ Hamiltonian truncation method (TFFSA) to extract the scattering phase of the lightest particles in the elastic regime, and combine it with S-matrix bootstrap methods based on unitarity and analyticity assumptions to determine the analytic continuation of the 2 → 2 S-matrix element to the complex s-plane. Focusing primarily on the “high temperature” regime in which the IFT interpolates between that of a weakly coupled massive fermion and the E 8 affine Toda theory, we will numerically determine 3-particle amplitudes, follow the evolution of poles and certain resonances of the S-matrix, and exclude the possibility of unknown wide resonances up to reasonably high energies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Statistical Uncertainty in Paleoclimate Proxy Reconstructions

A quantitative analysis of any environment older than the instrumental record relies on proxies. Uncertainties associated with proxy reconstructions are often underestimated, which can lead to artificial conflict between different proxies, and between data and models. In this paper, using ordinary least squares linear regression as a common example, we describe a simple, robust and generalizable method for quantifying uncertainty in proxy reconstructions. We highlight the primary controls on the magnitude of uncertainty, and compare this simple estimate to equivalent estimates from Bayesian, nonparametric and fiducial statistical frameworks. We discuss when it may be possible to reduce uncertainties, and conclude that the unexplained variance in the calibration must always feature in the uncertainty in the reconstruction. This directs future research toward explaining as much of the variance in the calibration data as possible. We also advocate for a “data-forward” approach, that clearly decouples the presentation of proxy data from plausible environmental inferences.

58 GEOSCIENCES↗
DOE OSTI2022

Super‐Droplet Method to Simulate Lagrangian Microphysics of Nuclear Fallout in a Homogeneous Cloud

Abstract Nuclear detonations produce hazardous local and global particles or fallout. Predicting fallout size, chemical components, and location is necessary to inform officials and determine immediate guidance for the public. However, existing nuclear detonation fallout models prescribe the particle size distributions based on limited observations. In this work, we apply the super‐droplet method, which is a numerical modeling technique developed for cloud microphysics, to simulate size distributions of particles in a mushroom cloud formed post‐detonation of a nuclear device. We model fallout formation and evolution with homogeneous nucleation and condensation of a single species and a Monte Carlo coagulation algorithm. We verify the numerical methods representing coagulation and condensation processes against analytical test problems. Additionally, we explore several scenarios for the integral system mass and yield in equivalent kilotons (kt) of TNT (trinitrotoluene). The fallout size distribution median diameter d pg follows a scaling law based on the integral system mass m v 0 kg and yield Y kt: nm. We test the effect of cloud turbulence, enhanced nucleation and growth, and vapor volatility with a sensitivity study. The range in median diameter predictions for simulations of historical tests performed over the Pacific encompass the measurements of particles sampled from the cloud caps. Predicted median particle size ranges up to 217, 123, 86, and 35 nm for historical tests with yields of 0.2, 0.7, 2, and 10 Mt, respectively. This work can be expanded in many different directions to build a more predictive model for fallout formation.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

InAs sidewall tunnel diodes enabled by surface states

Negative differential resistance (NDR), where the device current decreases with increasing bias voltage, is a representative phenomenon where quantum mechanics induces counterintuitive physical behavior and offers promising applications such as high-frequency oscillators, amplifiers, and multilevel logic circuits. While the NDR behavior has been extensively studied in various materials and devices, the role of surface properties in NDR, particularly in InAs-based diodes, remains underexplored. In this work, we report the observation of NDR in vertically structured InAs p + n diodes that exhibit a peak-valley current ratio of ∼6, which is suitably high for applications. Circumference-normalized current–voltage characterization revealed that the NDR originates from band-to-band tunneling between the valence band of p + -InAs and the conduction band of an n + -InAs surface, where the n + surface is due to surface states on otherwise n − -InAs. In addition, by comparing devices with various surface passivation methods (without intentional passivation, benzocyclobutene polymer, and silicon nitride), we found that the surface termination significantly affects the NDR characteristics. We present an equivalent circuit model to explain the observed device behavior. These findings offer insights into surface-enabled NDR phenomena and present new knobs for engineering NDR devices.

Circuit theorems↗
DOE OSTI2025

Quantifying Climate Change Effects of Bioenergy and BECCS: Critical Considerations and Guidance on Methodology

Bioenergy is a critical element in many national and international climate change mitigation efforts, including as a carbon dioxide removal strategy combined with the capture and durable geological storage of flue gas emissions (BECCS). However, divergent results on the effectiveness of bioenergy as a climate change mitigation measure are reported in the scientific literature. Climate impacts of bioenergy depend on case-specific factors, primarily biophysical features of the biomass production system, and the design and efficiency of conversion and capture processes. Estimates of climate impacts are also strongly affected by methodological choices and assumptions, and much of the divergence between studies derives from differences in the assumed alternate use of the land or feedstock, the alternate energy source and the system boundaries applied. We present a methodology to support robust estimates of the climate change effects of bioenergy systems, updating the standard methodology developed by the International Energy Agency's Technology Collaboration Program on Bioenergy. We provide guidance on the key choices including the reference land use and energy system that bioenergy is assumed to displace, spatial and temporal system boundaries, co-product handling, climate forcers considered, metrics applied and time horizon of impact assessment. Researchers should consider the whole bioenergy system including all life cycle stages, and choose system boundaries, reference systems and treatment of co-products that are consistent with the intended application of the results. The assessment should be normalised to a functional unit that can be compared with other systems delivering an equivalent quantity of the same function. All significant climate forcers should be included, and climate effects should be quantified using appropriate impact assessment methods that distinguish the impact of time. Consistency in methodology and interpretation will facilitate comparison between studies of different bioenergy systems.

09 BIOMASS FUELS↗
DOE OSTI2023

Expanding-box Quasilinear Model of the Solar Wind

The expanding-box model of the solar wind has been adopted in the literature within the context of magnetohydrodynamics, hybrid, and full particle-in-cell simulations to investigate the dynamic evolution of the solar wind. The present paper extends such a method to the framework of self-consistent quasilinear kinetic theory. It is shown that the expanding-box quasilinear methodology is largely equivalent to the inhomogeneous steady-state quasilinear model discussed earlier in the literature, but a distinction regarding the description of wave dynamics between the two approaches is also found. The expanding-box quasilinear formalism is further extended to include the effects of a spiraling solar-wind magnetic field as well as collisional age effects. The present finding shows that the expanding-box quasilinear approach and the steady-state global-kinetic models may be employed interchangeably in order to address other more complex problems associated with the solar-wind dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗