Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “characteristic equations”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Nonlinear dynamic simulation of single- and multi-spool core engines

In this paper a new computational method for accurate simulation of the nonlinear dynamic behavior of single- and multi-spool core engines, turbofan engines, and power generation gas turbine engines is presented. In order to perform the simulation, a modularly structured computer code has been developed which includes individual mathematical modules representing various engine components. The generic structure of the code enables the dynamic simulation of arbitrary engine configurations ranging from single-spool thrust generation to multi-spool thrust/power generation engines under adverse dynamic operating conditions. For precise simulation of turbine and compressor components, row-by-row calculation procedures were implemented that account for the specific turbine and compressor cascade and blade geometry and characteristics. The dynamic behavior of the subject engine is calculated by solving a number of systems of partial differential equations, which describe the unsteady behavior of the individual components. In order to ensure the capability, accuracy, robustness, and reliability of the code, comprehensive critical performance assessment and validation tests were performed. As representatives, three different transient cases with single- and multi-spool thrust and power generation engines were simulated. The transient cases range from operating with a prescribed fuel schedule, to extreme load changes, to generator and turbine shut down.

Schobeiri, T.↗
NASA NTRS1952

Matrix Method of Determining the Longitudinal-Stability Coefficients and Frequency Response of an Aircraft from Transient Flight Data

A matrix method is presented for determining the longitudinal-stability coefficients and frequency response of an aircraft from arbitrary maneuvers. The method is devised so that it can be applied to time-history measurements of combinations of such simple quantities as angle of attack, pitching velocity, load factor, elevator angle, and hinge moment to obtain the over-all coefficients. Although the method has been devised primarily for the evaluation of stability coefficients which are of primary interest in most aircraft loads and stability studies, it can be used also, with a simple additional computation, to determine the frequency-response characteristics. The entire procedure can be applied or extended to other problems which can be expressed by linear differential equations.

Donegan, James J↗
NASA NTRS1958

A Theoretical and Experimental Study of Planing Surfaces Including Effects of Cross Section and Plan Form

A summary is given of the background and present status of the pure-planing theory for rectangular flat plates and v-bottom surfaces. The equations reviewed are compared with experiment. In order to extend the range of available planing data, the principal planing characteristics for models having sharp bottom surfaces having constant angles of dead rise of 20 degrees and 40 degrees. Planing data were also obtained for flat-plate surfaces with very slightly rounded chines for which decreased lift and drag coefficients are obtained.

Shuford, Charles L , Jr↗
NASA NTRS1993

Proposed framework for thermomechanical life modeling of metal matrix composites

The framework of a mechanics of materials model is proposed for thermomechanical fatigue (TMF) life prediction of unidirectional, continuous-fiber metal matrix composites (MMC's). Axially loaded MMC test samples are analyzed as structural components whose fatigue lives are governed by local stress-strain conditions resulting from combined interactions of the matrix, interfacial layer, and fiber constituents. The metallic matrix is identified as the vehicle for tracking fatigue crack initiation and propagation. The proposed framework has three major elements. First, TMF flow and failure characteristics of in situ matrix material are approximated from tests of unreinforced matrix material, and matrix TMF life prediction equations are numerically calibrated. The macrocrack initiation fatigue life of the matrix material is divided into microcrack initiation and microcrack propagation phases. Second, the influencing factors created by the presence of fibers and interfaces are analyzed, characterized, and documented in equation form. Some of the influences act on the microcrack initiation portion of the matrix fatigue life, others on the microcrack propagation life, while some affect both. Influencing factors include coefficient of thermal expansion mismatch strains, residual (mean) stresses, multiaxial stress states, off-axis fibers, internal stress concentrations, multiple initiation sites, nonuniform fiber spacing, fiber debonding, interfacial layers and cracking, fractured fibers, fiber deflections of crack fronts, fiber bridging of matrix cracks, and internal oxidation along internal interfaces. Equations exist for some, but not all, of the currently identified influencing factors. The third element is the inclusion of overriding influences such as maximum tensile strain limits of brittle fibers that could cause local fractures and ensuing catastrophic failure of surrounding matrix material. Some experimental data exist for assessing the plausibility of the proposed framework.

Halford, Gary R.↗
NASA NTRS1993

Energy transfer and constrained simulations in isotropic turbulence

The defining characteristic of turbulent flows is their ability to dissipate energy, even in the limit of zero viscosity. The Euler equations, if constrained in such a way that the velocity derivatives remain bounded, conserve energy. But when they arise as the limit of the Navier-Stokes (NS) equations, when the Reynolds number goes to infinity, there is persuasive empirical evidence that the gradients become singular as just the right function of Re for the dissipation to remain non-zero and to approach a well defined limit. It is generally believed that this limiting value of the dissipation is a property of the Euler equations themselves, independent of the particular dissipative mechanism involved, and that it can be normalized with the large scale properties of the turbulent flow (e.g. the kinetic energy per unit volume u'(exp 2)/2, and the integral scale L) without reference to the Reynolds number or to other dissipative quantities. This is usually taken to imply that the low wave number end of the energy spectrum, far from the dissipative range, is also independent of the particular mechanism chosen to dispose of the energy transfer. In the following sections, we present some numerical experiments on the effect of substituting different dissipation models into the truncated Euler equations. We will see that the effect is mainly felt in the 'near dissipation' range of the energy spectrum, but that this range can be quite wide in some cases, contaminating a substantial range of wave numbers. In the process, we will develop a 'practical' approximation to the subgrid energy transfer in isotropic turbulence, and we will gain insight into the structure of the nonlinear interactions among turbulent scales of comparable size, and into the nature of energy backscatter. Some considerations on future research directions are offered at the end.

Jimenez, Javier↗
NASA NTRS1995

Numerical simulations of small amplitude acoustic wave propagation in a converging-diverging nozzle

Numerical simulations of the small amplitude acoustic wave propagation in a converging-diverging nozzle proposed in Category 5 are performed and presented here. The quasi 1-D unsteady flow equations in conservative form are discretized by the DRP scheme with the artificial damping developed by Tam and Webb. Characteristic boundary conditions are used at both the inlet and outlet of the nozzle. The effect of different numerical implementations of the subsonic inflow boundary conditions on the convergence of the solution to the steady state is studied. In the case of a subsonic outflow in which a shock is formed in the nozzle, the interaction between acoustics and the shock is also investigated.

Dong, Thomas Z.↗
NASA NTRS1959

Modified Matrix Method for Calculating Steady-State Span Loading on Flexible Wings in Subsonic Flight

A method is presented for shortening the computations required to determine the steady-state span loading on flexible wings in subsonic flight. The method makes use of tables of downwash factors to find the necessary aerodynamic-influence coefficients for the application of lifting-line theory. Explicit matrix equations of equilibrium are converted into a matrix power series with a finite number of terms by utilizing certain characteristic properties of matrices. The number of terms in the series is determined by a trial-and-error process dependent upon the required accuracy of the solution. Spanwise distributions of angle of attack, airload, shear, bending moment, and pitching moment are readily obtained as functions of qm(sub R) where q denotes the dynamic pressure and mR denotes the lift-curve slope of a rigid wing. This method is intended primarily to make it practical to solve steady-state aeroelastic problems on the ordinary manually operated desk calculators, but the method is also readily adaptable to automatic computing equipment.

Gainer, Patrick A.↗
NASA NTRS2002

Towers for Earth Launch

This report lists some characteristics of a hypothetical 15 kilometer tower for launching spacecraft, the advantages of launching from high altitude, and some equations pertaining to launch from a 15 kilometer tower.

Landis, Geoffrey A.↗
NASA NTRS2004

Computed Linear/Nonlinear Acoustic Response of a Cascade for Single/Multi Frequency Excitation

This paper examines mode generation and propagation characteristics of a 2-D cascade due to incident vortical disturbances using a time domain approach. Full nonlinear Euler equations are solved employing high order accurate spatial differencing and time marching techniques. The solutions show the generation and propagation of mode orders that are expected from theory. Single frequency excitations show linear response over a wide range of amplitudes. The response for multi-frequency excitations tend to become nonlinear due to interaction between frequencies and self interaction.

Nallasamy, M.↗
NASA NTRS2016

Numerical Simulation of a Complete Low-Speed Wind Tunnel Circuit

A numerical simulation of the complete circuit of the NASA Langley 14 x 22-ft low-speed wind tunnel is described. Inside the circuit, all turning vanes are modeled as well as the five flow control vanes downstream of the 1st corner. The fan drive system is modeled using an actuator disk for the fan blades coupled with the fan nacelle. All the surfaces are modeled as viscous walls except the turning vanes, which were modeled as inviscid surfaces. NASA Langley's TetrUSS unstructured grid software was used for grid generation and flow simulation. Two turbulence models were employed in the present study, namely, the one-equation Spalart-Allmaras model and the shear stress transport (SST) model of Menter. The paper shows the flow characteristics in the circuit and compares the results with experimental data where available.

Nayani, Sudheer N.↗
NASA NTRS2018

Metabrain for Embedded Cognition (MBEC)

This study presents the application of Hidden Markov Models (HMM) to determine specialized features without expert input. Specifically, the application of such a method for classification of high multi-path fading is targeted, for demonstrating the feasibility of such an approach. This is the first step in the development of a meta-brain for embedded cognition (M-BEC) suite that can be used to apply machine learning to various communication systems at NASA GRC. The project explores the concept of fading and how it affects communication systems in a negative way. Currently, supervised learning methods are used to study the effects of fading on space links. However, such models rely on expert features to make predictions as to the state of a link and whether fading is present. This project offers the possibility of having the HMM learn what characteristics are important and make predictions based on those characteristics. This project explores Hidden Markov Models, their theory and applications to various problems, as well as the underlying equations and assumptions. A preliminary result is presented and recommendations are made as to the use of such an approach for communications systems.

Propagation↗
NASA NTRS2022

A Generalized Dual-Frequency Ratio (DFR) Approach for Rain Retrievals

The dual-frequency ratio of radar reflectivity factors (DFR) has been shown to be a useful quantity as it is independent of the number concentration of the particle size distribution and primarily a function of the mass-weighted particle diameter, Dm. A drawback of DFR-related methods for rain estimation, however, is the non-unique relationship between Dm and DFR. At Ku- and Ka-band frequencies, two solutions for Dm exist when DFR is less than zero. This ambiguity generates multiple solutions for the range profiles of the particle size parameters. We investigate characteristics of these solutions for both the initial-value (forward) and final-value (backward) forms of the equations. To choose one of among many possible range profiles of Dm, number concentration and rain rate, R, independently measured path attenuations are used. For the backward approach the possibility exists of dispensing with externally measured path attenuations by achieving consistency between the input and output path attenuations. The methods are tested by means of a simulation based on disdrometer-measured raindrop size distributions and the results are compared with a simplified version of the operational R-Dm method.

Robert Meneghini↗
NASA NTRSDate not supplied

Separation and Transition on a Cone-Cylinder-Flare: Computational Investigations

Base flow computation and stability analysis were conducted for hypersonic flow over a cone-cylinder-flare (CCF) geometry for conditions that correspond to the experimental runs carried out in three wind tunnels. Owing to the presence of an attached boundary layer, a separation bubble induced by a shock-boundary layer interaction, and a reattachment region, the chosen flow configuration is physically rich. The complexity of this flowfield encompasses a combination of convective instabilities developing on the cone, global instabilities in the separation bubble, and shear-layer modes and streaks in the reattachment region. Thus, the selected CCF configuration provides the opportunity for a comprehensive comparison of the currently available methodologies for analyzing boundary layer instabilities. Various tools are used for the analysis, including global stability codes as well as convective instability analyses based on a local theory, a weakly non-parallel analysis, and tools that are applicable to strongly non-parallel flows. The paper presents a comparison of the convective instability characteristics based on different methodologies, such as linear stability theory (LST), the harmonic form of linearized Navier-Stokes equations (HLNSE), and resolvent analysis. The CCF configuration provided an effective framework for conducting a detailed cross-validation of this type, which had not yet been addressed in existing literature. This document is accompanied by a companion paper that is focused on the experimental aspects of the CCF configuration. Both papers are being presented in a dedicated session that highlights the research activities of the NATO STO Research Task Group AVT-346.

Boundary Layer Transition↗
NASA NTRSDate not supplied

Shock Capturing via Limiting for High-Order Methods including Discontinuous Galerkin

High-order methods, such as discontinuous Galerkin (DG), spectral, and flux reconstruction (FR), are prone to generating unwanted oscillations near shocks and discontinuities. Conventional limiting techniques, while effective in suppressing oscillations near shocks, often compromise accuracy near extrema, where the solution is only first-order accurate. This paper introduces a novel limiting technique for these high-order schemes, aimed at effectively managing shocks while preserving accuracy. The key idea is to expand the standard monotonicity limits to provide “room” near smooth extrema, ensuring that limiting has no effect and thus preserving accuracy. Near a discontinuity, these expanded limits effectively reduce to the original monotonicity limits, suppressing oscillations. Additional motivation is drawn from a formula for the derivative of Radau polynomials, which depicts the behavior of oscillations resulting from discontinuities. This behavior leads to a simplification by applying the limits to the sum of magnitudes of all modes, linear and higher degree. Unlike typical approaches, which rely on successful detection to activate limiting, our limiter depends continuously on the data, there by avoiding potential issues if detection fails. To reduce computing time, efficient criteria for detecting smooth regions where limiting is unnecessary are presented. Combined with detection, the continuous dependence on the data is lost, but the method is more economical. A notable characteristic of the entire process is its simplicity in both concept and implementation. Numerical tests for advection and Euler equations are conducted to demonstrate the effectiveness of the proposed method.

numerical methods↗
NASA NTRSDate not supplied

Shock Capturing via Limiting for High-Order Methods including Discontinuous Galerkin

High-order methods, such as discontinuous Galerkin (DG), spectral, and flux reconstruction (FR), are prone to generating unwanted oscillations near shocks and discontinuities. Conventional limiting techniques, while effective in suppressing oscillations near shocks, often compromise accuracy near extrema, where the solution is only first-order accurate. This paper introduces a novel limiting technique for these high-order schemes, aimed at effectively managing shocks while preserving accuracy. The key idea is to expand the standard monotonicity limits to provide “room” near smooth extrema, ensuring that limiting has no effect and thus preserving accuracy. Near a discontinuity, these expanded limits effectively reduce to the original monotonicity limits, suppressing oscillations. Additional motivation is drawn from a formula for the derivative of Radau polynomials, which depicts the behavior of oscillations resulting from discontinuities. This behavior leads to a simplification by applying the limits to the sum of magnitudes of all modes, linear and higher degree. Unlike typical approaches, which rely on successful detection to activate limiting, our limiter depends continuously on the data, there by avoiding potential issues if detection fails. To reduce computing time, efficient criteria for detecting smooth regions where limiting is unnecessary are presented. Combined with detection, the continuous dependence on the data is lost, but the method is more economical. A notable characteristic of the entire process is its simplicity in both concept and implementation. Numerical tests for advection and Euler equations are conducted to demonstrate the effectiveness of the proposed method.

numerical methods↗
DOE OSTI2021

Shifts in hydroclimatology of US megaregions in response to climate change

Most of the population and economic growth in the United States occurs in megaregions as the clustered metropolitan areas, whereas climate change may amplify negative impacts on water and natural resources. This study assesses shifts in regional hydroclimatology of fourteen US megaregions in response to climate change over the 21st century. Hydroclimatic projections were simulated using the Variable Infiltration Capacity (VIC) model driven by three downscaled climate models from the Multivariate Adaptive Constructed Analogs (MACA) dataset to cover driest to wettest future conditions in the conterminous United States (CONUS). Shifts in the regional hydroclimatolgy and basin characteristics of US megaregions were represented as a combination of changes in the aridity and evaporative indices using the Budyko framework and Fu's equation. Changes in the climate types of US megaregions were estimated using the Fine Gaussian Support Vector Machine (SVM) method. The results indicate that Los Angeles, San Diego, and San Francisco are more likely to experience less arid conditions with some shifts from Continental to Temperate climate type while the hydroclimatology of Houston may become drier with some shifts from Temperate to Continental climate type. Additionally, water yield is likely to decrease in Seattle. Change in the hydroclimatology of Denver and Phoenix highly depends on the selected climate model. However, the basin characteristics of Phoenix have the highest sensitivity to climate change. Overall, the hydroclimatic conditions of Los Angeles, San Diego, Phoenix, Denver, and Houston have the highest sensitivity to climate change. Understanding of future shifts in hydroclimatology of megaregions can help decision-makers to attenuate negative consequences by implementing appropriate adaptation strategies, particularly in the water-scare megaregions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Validation and parameterization of a novel physics-constrained neural dynamics model applied to turbulent fluid flow

We report, in fluid physics, data-driven models to enhance or accelerate time to solution are becoming increasingly popular for many application domains, such as alternatives to turbulence closures, system surrogates, or for new physics discovery. In the context of reduced order models of high-dimensional time-dependent fluid systems, machine learning methods grant the benefit of automated learning from data, but the burden of a model lies on its reduced-order representation of both the fluid state and physical dynamics. In this work, we build a physics-constrained, data-driven reduced order model for Navier–Stokes equations to approximate spatiotemporal fluid dynamics in the canonical case of isotropic turbulence in a triply periodic box. The model design choices mimic numerical and physical constraints by, for example, implicitly enforcing the incompressibility constraint and utilizing continuous neural ordinary differential equations for tracking the evolution of the governing differential equation. We demonstrate this technique on a three-dimensional, moderate Reynolds number turbulent fluid flow. In assessing the statistical quality and characteristics of the machine-learned model through rigorous diagnostic tests, we find that our model is capable of reconstructing the dynamics of the flow over large integral timescales, favoring accuracy at the larger length scales. More significantly, comprehensive diagnostics suggest that physically interpretable model parameters, corresponding to the representations of the fluid state and dynamics, have attributable and quantifiable impact on the quality of the model predictions and computational complexity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1977

Non-linear duct acoustics and its application to fan noise

Quite often the scalar Helmholtz equation is assumed to be the fundamental equation of duct acoustics and various relations are derived from it. An investigation is, therefore, conducted regarding the underlying assumptions leading to the Helmholtz equation. It is found that serious errors, even well below 186 dB, are possible if the assumption of linearity were to be made. In duct acoustics these errors are highest near the cutoff points and in the case of ducts with transonic flow. It is pointed out that serious errors could arise if a nonlinear problem is analyzed by linear methods. At times, an entire phenomenon remains unpredictable by linear methods. Attention is given to the characteristics of nonlinear systems, the phenomenon of multiple pure tones in the case of noise from the aircraft engine fans, the limitations of the conventional acoustics equation, duct shape and nonlinearity, and the numerical simulation of a fan with 8 blades.

Vaidya, P. G.↗