SEARCH · Search NASA

Results for “variable coefficients”

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 253 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS2006

Singular-Arc Time-Optimal Trajectory of Aircraft in Two-Dimensional Wind Field

This paper presents a study of a minimum time-to-climb trajectory analysis for aircraft flying in a two-dimensional altitude dependent wind field. The time optimal control problem possesses a singular control structure when the lift coefficient is taken as a control variable. A singular arc analysis is performed to obtain an optimal control solution on the singular arc. Using a time-scale separation with the flight path angle treated as a fast state, the dimensionality of the optimal control solution is reduced by eliminating the lift coefficient control. A further singular arc analysis is used to decompose the original optimal control solution into the flight path angle solution and a trajectory solution as a function of the airspeed and altitude. The optimal control solutions for the initial and final climb segments are computed using a shooting method with known starting values on the singular arc The numerical results of the shooting method show that the optimal flight path angle on the initial and final climb segments are constant. The analytical approach provides a rapid means for analyzing a time optimal trajectory for aircraft performance.

Nguyen, Nhan↗

NASA NTRS1989

Optimum interface properties for metal matrix composites

Due to the thermal expansion coefficient mismatch (CTE) between the fiber and the matrix, high residual sresses exist in metal matrix composite systems upon cool down from processing temperature to room temperature. An interface material can be placed between the fiber and the matrix to reduce the high tensile residual stresses in the matrix. A computer program was written to minimize the residual stress in the matrix subject to the interface material properties. The decision variables are the interface modulus, thickness and thermal expansion coefficient. The properties of the interface material are optimized such that the average distortion energy in the matrix and the interface is minimized. As a result, the only active variable is the thermal expansion coefficient. The optimum modulus of the interface is always the minimum allowable value and the interface thickness is always the maximum allowable value, independent of the fiber/matrix system. The optimum interface thermal expansion coefficient is always between the values of the fiber and the matrix. Using this analysis, a survey of materials was conducted for use as fiber coatings in some specific composite systems.

Ghosn, Louis J.↗

NASA NTRS2019

Characterization of Flat-Plate Heat Pipe Functionality for Fuel Cell Application

For many exothermic systems such as proton exchange membrane (PEM) fuel cells, cooling is necessary to maintain consistent and uniform operating temperatures which are generally in the range of 60 to 80°C. A potential solution is a heat pipe employing the principle of conduction in combination with a phase change. The purpose of this study was to fully characterize the operation of a heat pipe with unique geometry designed for fuel cell heat removal at operating temperatures up to 100°C and to assess the suitability of the resulting performance for cooling system integration. Test variables consist of cold sink thermal contact conductance coefficient improving material usage, heat rejection temperature, contact pressure and area, heat pipe orientation, operational hours, and variability between units. The test results showed that heat pipe performance improves through increased condenser and heat sink contact pressure, thermal contact coating usage, greater condenser area utilization, and colder heat sink temperature. Higher initial heat pipe temperature and greater age or operation time reduce the heat pipe capabilities. Placing the condenser lower than the evaporator is the only heat pipe orientation that affects evaporator plate temperature. Condensing fluid in the lower section reduces heat transfer capacity as that orientation requires liquid water to transport by wicking against gravity.

Smith, Phillip J.↗

DOE OSTI2026

Statistical Uncertainty of Inhalation Dose Coefficients in Consequence Management: Propagated Dose Uncertainty in ICRP 66 Human Respiratory Tract Model

Reference inhalation dose models rely on deterministic biokinetics and reference computational phantoms, limiting their applicability to the variability present in population-specific exposures encountered in emergency response scenarios. Here, this study introduces REDCAL, a Python-based computational framework developed to propagate uncertainty in inhalation dose coefficients using the International Commission on Radiological Protection (ICRP) Publication 66 Human Respiratory Tract Model. REDCAL integrates ICRP deposition and clearance models, systemic biokinetics, and governing physics principles, and leverages Sandia National Laboratories’ Dakota toolkit for uncertainty quantification via Latin Hypercube Sampling. REDCAL was validated against DCAL, with biokinetic retention results differing by less than 1% and effective dose coefficients by less than 2% across all tested radionuclides. Stochastic sampling introduced variability in dose coefficients, with geometric standard deviations (GSD) in committed effective dose coefficients (CEDC) ranging from 1.0 to 1.5, based on lognormal distribution fits. Analysis demonstrated that variations in the activity median aerodynamic diameter (AMAD) notably influenced the computed CEDC values. Smaller particles (<1 µm) increased doses by 20–30% due to deeper lung deposition and prolonged retention for alpha emitting radionuclides, such as 241 Am and 239 Pu. Radionuclides with fast clearance, such as 133 I, demonstrated a dose reduction exceeding 50%, as AMAD increased beyond 5 µm due to upper airway deposition and rapid mucociliary clearance. The greatest GSD among the radionuclides reported in this study was for 241 Am. In most cases, the largest GSDs in the CEDC were associated with larger particle sizes, an expected outcome, as ICRP Publication 66 defines GSD in particle size as a function of AMAD, resulting in an extended tail of the lognormal distribution. The findings support improved inhalation dose assessments and enhance consequence management strategies for the U.S. Federal Radiological Monitoring and Assessment Center by quantifying uncertainty in dose coefficients and strengthening decision-making for emergency response scenarios.

Biokinetic Modeling↗

NASA NTRS1980

Hyperstability in variable structure systems

A variable structure control law is a discontinuous feedback law whose coefficients switch on hypersurfaces defined in the state space. Feedback systems with variable structure control laws are often referred to as variable structure systems (VSS). The main feature of VSS is the sliding motion which can occur at the intersection of the switching surfaces. Demonstrating existence and reaching in the case of vector control in n-dimensional space is in general a problem in stability of nonlinear systems. The present investigation is concerned with the projection of motion on an m-dimensional subspace, taking into account aspects of the hyperstability concept in VSS design. The considered concepts are illustrated with the aid of an example.

Calise, A. J.↗

NASA NTRS2011

Broken Ergodicity in MHD Turbulence in a Spherical Domain

Broken ergodicity (BE) occurs in Fourier method numerical simulations of ideal, homogeneous, incompressible magnetohydrodynamic (MHD) turbulence. Although naive statistical theory predicts that Fourier coefficients of fluid velocity and magnetic field are zero-mean random variables, numerical simulations clearly show that low-wave-number coefficients have non-zero mean values that can be very large compared to the associated standard deviation. In other words, large-scale coherent structure (i.e., broken ergodicity) in homogeneous MHD turbulence can spontaneously grow out of random initial conditions. Eigenanalysis of the modal covariance matrices in the probability density functions of ideal statistical theory leads to a theoretical explanation of observed BE in homogeneous MHD turbulence. Since dissipation is minimal at the largest scales, BE is also relevant for resistive magnetofluids, as evidenced in numerical simulations. Here, we move beyond model magnetofluids confined by periodic boxes to examine BE in rotating magnetofluids in spherical domains using spherical harmonic expansions along with suitable boundary conditions. We present theoretical results for 3-D and 2-D spherical models and also present computational results from dynamical simulations of 2-D MHD turbulence on a rotating spherical surface. MHD turbulence on a 2-D sphere is affected by Coriolus forces, while MHD turbulence on a 2-D plane is not, so that 2-D spherical models are a useful (and simpler) intermediate stage on the path to understanding the much more complex 3-D spherical case.

Shebalin, John V.↗

NASA NTRS1983

Equivalent absorption coefficients generated from frequency probability distributions

A flexible and computationally accurate method of treating aerosol scattering in spectral regions in which gaseous absorption is important is described. In the method, line-by-line absorption coefficients are computed as a function of pressure, temperature, and absorber gas for the spectral region of interest. The coefficients are sorted into a probability distribution which is converted into a cumulative probability distribution, which in turn can be inverted due to its monotonic nature. The inverted distribution is a smooth curve giving the absorption coefficient as a function of an independent variable on the domain. The frequency integration of the radiative transfer equation can then be performed by a quadrature technique with values of the absorption coefficient determined from the inverted distribution curve. The method is illustrated by applying it to the 9.6 micron band of ozone.

Ackerman, T. P.↗

NASA NTRS1999

Optimal Frequency-Domain System Realization with Weighting

Several approaches are presented to identify an experimental system model directly from frequency response data. The formulation uses a matrix-fraction description as the model structure. Frequency weighting such as exponential weighting is introduced to solve a weighted least-squares problem to obtain the coefficient matrices for the matrix-fraction description. A multi-variable state-space model can then be formed using the coefficient matrices of the matrix-fraction description. Three different approaches are introduced to fine-tune the model using nonlinear programming methods to minimize the desired cost function. The first method uses an eigenvalue assignment technique to reassign a subset of system poles to improve the identified model. The second method deals with the model in the real Schur or modal form, reassigns a subset of system poles, and adjusts the columns (rows) of the input (output) influence matrix using a nonlinear optimizer. The third method also optimizes a subset of poles, but the input and output influence matrices are refined at every optimization step through least-squares procedures.

Juang, Jer-Nan↗

NASA NTRS1999

Optimal System Realization in Frequency Domain

Several approaches are presented to identify an experimental system model directly from frequency response data. The formulation begins with a matrix-fraction description as the model structure. Frequency weighting such as exponential weighting is introduced to solve a weighted least-squares problem to obtain the coefficient matrices for the matrix-fraction description. A multi-variable state-space model can then be formed using the coefficient matrices of the matrix-fraction description. An approach is introduced to fine-tune the model using non-linear programming methods to minimize the desired cost function. The method deals with the model in the real Schur or modal form and reassigns a subset of system poles using a nonlinear optimizer. At every optimization step, the input and output influence matrices are refined through least-squares procedures. The proposed approaches are used to identify an analytical model for a NASA testbed from experimental data.

Juang, Jer-Nan↗

NASA NTRS2019

Exploiting the Convergence of Evidence in Satellite Data for Advanced Weather Index Insurance Design

The goal of drought-related weather index insurance (WII) is to protect smallholder farmers against the risk of weather shocks and to increase their agricultural productivity. Estimates of precipitation and vegetation greenness are the two dominant satellite datasets. However, ignoring additional moisture- and energy-related processes that influence the response of vegetation to rainfall leads to an incomplete representation of the hydrologic cycle. This study evaluates the added value of considering multiple independent satellite-based variables to design, calibrate, and validate weather insurance indices on the African continent. The satellite data include two rainfall datasets, soil moisture, the evaporative stress index (ESI), and vegetation greenness. We limit artificial advantages by resampling all datasets to the same spatial (0.25°) and temporal (monthly) resolution, although datasets with a higher spatial resolution might have an added value, if considered as the single source of information for localized applications. A higher correlation coefficient between the moisture-focused variables and the normalized difference vegetation index (NDVI), an indicator for vegetation vigor, provides evidence for the datasets’ capability to capture agricultural drought conditions on the ground. The Climate Hazards Group Infrared Precipitation with Stations (CHIRPS) rainfall dataset, soil moisture, and ESI show higher correlations with the (lagged) NDVI in large parts of Africa, for different land covers and various climate zones, than the African Rainfall Climatology, version 2 (ARC2), rainfall dataset, which is often used in WII. A comparison to drought years as reported by farmers in Ethiopia, Senegal, and Zambia indicates a high “hit rate” of all satellite-derived anomalies regarding the detection of severe droughts but limitations regarding moderate drought events.

drought↗

NASA NTRS2002

Vibration Reduction of Helicopter Blade Using Variable Dampers: A Feasibility Study

In the report, the investigation of controlling helicopter-blade lead-lag vibration is described. Current practice of adding passive damping may be improved to handle large dynamic range of the blade with several peaks of vibration resonance. To minimize extra-large damping forces that may damage the control system of blade, passive dampers should have relatively small damping coefficients, which in turn limit the effectiveness. By providing variable damping, a much larger damping coefficient to suppress the vibration can be realized. If the damping force reaches the maximum allowed threshold, the damper will be automatically switched into the mode with smaller damping coefficient to maintain near-constant damping force. Furthermore, the proposed control system will also have a fail-safe feature to guarantee the basic performation of a typical passive damper. The proposed control strategy to avoid resonant regions in the frequency domain is to generate variable damping force in combination with the supporting stiffness to manipulate the restoring force and conservative energy of the controlled blade system. Two control algorithms are developed and verified by a prototype variable damper, a digital controller and corresponding algorithms. Primary experiments show good potentials for the proposed variable damper: about 66% and 82% reductions in displacement at 1/3 length and the root of the blade respectively.

Lee, George C.↗

NASA NTRS1999

Selection of Optical Glasses Using Buchdahl's Chromatic Coordinate

This investigation attempted to extend the method of reducing the size of glass catalogs to a global glass selection technique with the hope of guiding glass catalog offerings. Buchdahl's development of optical aberration coefficients included a transformation of the variable in the dispersion equation from wavelength to a chromatic coordinate omega defined as omega = (lambda - lambda(sub 0))/ 1 + 2.5(lambda - lambda(sub 0)) where lambda is the wavelength at which the wavelength is calculated and lambda(sub 0) is a base wavelength about which the expansion is performed. The advantage of this approach is that the dispersion equation may be written in terms of a simple power series and permits direct calculation of dispersion coefficients. While several promising examples were given, a systematic application of the technique to an entire glass catalog and analysis of the subsequent predictions was not performed. The goal of this work was to apply the technique in a systematic fashion to glasses in the Schoft catalog and assess the quality of the predictions.

Griffin, DeVon W.↗

DOE OSTI2025

Optimal Complementarity Analysis of Potential Floating Solar Co-Located With Existing Hydropower Assets Across the Contiguous United States

The U.S. is expected to double its rate of renewable capacity from 2024 to 2030. However, the stochastic nature of renewable energy poses challenges to the operation and reliability of our power grid. The combined generation from renewable energy sources, with dispatchable sources (such as hydropower) operating as a hybrid energy plant, could mitigate this variability. In this paper, the complementarity analysis of selected U.S. reservoirs with existing hydropower assets (EHAs) and potential floating photovoltaics (FPVs) is conducted for the continuous U.S. (CONUS). The optimal FPV capacity for each site is determined by minimizing the variability of the combined output, while adhering to the FPV potential. Our results indicate that over 50% of the analyzed reservoirs achieve a stability coefficient exceeding 0.5, leading to a less-variable output after optimization. Finally, we analyze the complementary hydro-FPV hybrid reservoirs by considering both the Pearson correlation coefficient and the stability coefficient on daily, monthly, and yearly scales. Summaries are included of locations of theoretical FPVs co-located with hydropower plants that exhibit high complementarity based on the selected metrics.

13 - HYDRO ENERGY↗

NASA NTRS2002

A Numerical Scheme for Ordinary Differential Equations Having Time Varying and Nonlinear Coefficients Based on the State Transition Matrix

A variable order method of integrating initial value ordinary differential equations that is based on the state transition matrix has been developed. The method has been evaluated for linear time variant and nonlinear systems of equations. While it is more complex than most other methods, it produces exact solutions at arbitrary time step size when the time variation of the system can be modeled exactly by a polynomial. Solutions to several nonlinear problems exhibiting chaotic behavior have been computed. Accuracy of the method has been demonstrated by comparison with an exact solution and with solutions obtained by established methods.

Bartels, Robert E.↗

NASA NTRS2004

Modeling Turbulent Combustion for Variable Prandtl and Schmidt Number

This report consists of two abstracts submitted for possible presentation at the AIAA Aerospace Science Meeting to be held in January 2005. Since the submittal of these abstracts we are continuing refinement of the model coefficients derived for the case of a variable Turbulent Prandtl number. The test cases being investigated are a Mach 9.2 flow over a degree ramp and a Mach 8.2 3-D calculation of crossing shocks. We have developed an axisymmetric code for treating axisymmetric flows. In addition the variable Schmidt number formulation was incorporated in the code and we are in the process of determining the model constants.

Hassan, H. A.↗

NASA NTRS2008

Vertical Jump Height is more Strongly Associated with Velocity and Work Performed Prior to Take-off

Vertical jump (VJ) height is commonly used as a measure of athletic capability in strength and power sports. Although VJ has been shown to be a predictor of athletic performance, it is not clear which kinetic ground reaction force (GRF) variables, such as peak force (PF), peak power (PP), peak velocity (PV), total work (TW) or impulse (Imp) are the best correlates. To determine which kinetic variables (PF, PP, PV, TW, and Imp) best correlate with VJ height. Twenty subjects (14 males, 6 females) performed three maximal countermovement VJs on a force platform (Advanced Mechanical Technology, Inc., Watertown, MA, USA). VJ jump height was calculated as the difference between standing reach and the highest reach point measured using a Vertec. PF, PP, PV, TW, and Imp were calculated using the vertical GRF data sampled at 1000 Hz from the lowest point in the countermovement through the concentric portion until take-off. GRF data were normalized to body mass measured using a standard scale (Detecto, Webb City, MO, USA). Correlation coefficients were computed between each GRF variable and VJ height using a Pearson correlation. VJ height (43.4 plus or minus 9.1 cm) was significantly correlated (p less than 0.001) with PF (998 plus or minus 321 N; r=0.51), PP (1997 plus or minus 772 W; r=0.69), PV (2.66 plus or minus 0.40 m (raised dot) s(sup -1); r=0.85), TW (259 plus or minus 93.0 kJ; r=0.82), and Imp (204 plus or minus 51.1 N(raised dot)s; r=0.67). Although all variables were correlated to VJ height, PV and TW were more strongly correlated to VJ height than PF, PP, and Imp. Therefore, since TW is equal to force times displacement, the relative displacement of the center of mass along with the forces applied during the upward movement of the jump are critical determinants of VJ height. PV and TW are key determinants of VJ height, and therefore successful training programs to increase VJ height should focus on rapid movement (PV) and TW by increasing power over time rather than focusing on PF alone.

Bentley, J. R.↗

NASA NTRS2005

Qualifying a Bonding Process for the Space Interferometry Mission

The Space Interferometry Mission consists of three parallel Michelson interferometers that will be capable of detecting extrasolar planets with a high degree of accuracy and precision. High levels of stability must be met in order to fulfill the scientific requirements of this mission. To attain successful measurements the coefficient of thermal expansion between optics and bonding material must be minimized without jeopardizing the integrity of the bonds. Optic-to-optic bonds have been analyzed to better understand variables such as the effects of the coefficient of thermal expansion differences between optics and bonding materials, and materials have been chosen for the project based on these analyses. A study was conducted to determine if a reliable, repeatable process for bonding by wicking adhesive could be obtained using a low-viscosity epoxy and ultra-low expansion glass. A process of creating a methodology of bonding fused silica optics with Z-6020 silane primer and Epo-Tek 301 epoxy will be discussed.

Joyce, Gretchen P.↗

NASA NTRS1992

Reduced basis technique for calculating sensitivity coefficients of nonlinear structural response

An efficient reduced basis technique is presented for calculating the sensitivity of nonlinear structural response to variations in the design variables. The structure is discretized by using two-field mixed finite element models. The vector of structural response and its sensitivity coefficients (derivatives with respect to design variables) are each expressed as a linear combination of a small number of basis (or global approximation) vectors. The Bubnov-Galerkin technique is then used to approximate each of the finite element equations governing the response and the sensitivity coefficients by a small number of algebraic equations in the amplitudes of these vectors. The path derivatives (derivatives of the response vector with respect to path parameters, e.g., load parameters) are used as basis vectors for approximating the response. A combination of the path derivatives and their derivatives with respect to the design variables is used for approximating the sensitivity coefficients. The potential of the proposed technique is discussed and its effectiveness is demonstrated by means of numerical examples of laminated composite plates subjected to mechanical and thermal loads.

Noor, Ahmed K.↗