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,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Application of a full potential method for computation of three-dimensional supersonic flows

A nonlinear aerodynamic analysis technique based on the full potential equation in conservative form has been modified to permit treatment of supersonic flows with embedded subsonic regions (typically near the fuselage-canopy juncture and the wing leading edge). Solution procedures for the equations do not require any specific form of geometry or physical grid system. This results in the capability to analyze easily very complex geometries provided the posed problem lies within the isentropic restrictions of the full potential theory. Characteristic signal propagation theory is used to monitor the type dependent flow and a conservative switching scheme is employed to transition from the supersonic marching algorithm to a subsonic relaxation procedure and vice versa. An implicit approximate factorization scheme is used to solve the finite-difference equations. These modifications now permit analysis of fully three-dimensional flowfields including the interference effects due to lifting surface wakes. Improved grid generation capability allows analysis of complete complex aircraft geometries (fuselage, wing, tail, wing wake, and tail wake). Results are presented showing very good correlations with experimental surface pressure data and aerodynamic force data at both design and off-design operating points. Configurations examined include several waverider concepts, an arrow wing-body with wake, an advanced tactical fighter concept, and a fighter forebody-canard configuration.

Jones, K. M.↗
DOE OSTI2024

Constraints on initial baryon stopping and equation of state from directed flow

Our investigation focuses on the rapidity-dependent directed flow, v 1 (y), of identified hadrons in Au+Au collisions across a broad range of √S NN from 7.7 to 200 GeV. Employing a (3+1)-dimensional hybrid framework, our study successfully reproduces the characteristic features of the measured v 1 (y) for both mesons and baryons across the considered beam energies. Notably, our analysis reveals the constraining power of baryonic v 1 (y) on the initial baryon stopping mechanism. Together with mesonic v 1 (y), the directed flow serves as a crucial tool for probing the equation of state governing dense nuclear matter at finite chemical potentials.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1979

An analytical technique for predicting the characteristics of a flexible wing equipped with an active flutter-suppression system and comparison with wind-tunnel data

An analytical technique for predicting the performance of an active flutter-suppression system is presented. This technique is based on the use of an interpolating function to approximate the unsteady aerodynamics. The resulting equations are formulated in terms of linear, ordinary differential equations with constant coefficients. This technique is then applied to an aeroelastic model wing equipped with an active flutter-suppression system. Comparisons between wind-tunnel data and analysis are presented for the wing both with and without active flutter suppression. Results indicate that the wing flutter characteristics without flutter suppression can be predicted very well but that a more adequate model of wind-tunnel turbulence is required when the active flutter-suppression system is used.

Abel, I.↗
NASA NTRS1952

Transonic flow past a wedge profile with detached bow wave

A theoretical study has been made of the aerodynamic characteristics at zero angle of attack of a thin, doubly symmetrical double-wedge profile in the range of supersonic flight speed in which the bow wave is detached. The analysis utilizes the equations of the transonic small-disturbance theory and involves no assumptions beyond those implicit in this theory. The mixed flow about the front half of the profile is calculated by relaxation solution of boundary conditions along the shock polar and sonic line. The purely subsonic flow about the rear of the profile is found by means of the method of characteristics specialized to the transonic small-disturbance theory. Complete calculations were made for four values of the transonic similarity parameter. These were found sufficient to bridge the gap between the previous results of Guderley and Yoshihara at a Mach number of 1 and the results which are readily obtained when the bow wave is attached and the flow is completely supersonic.

Vincenti, Walter G↗
NASA NTRS2016

Optimal Growth in Hypersonic Boundary Layers

The linear form of the parabolized linear stability equations is used in a variational approach to extend the previous body of results for the optimal, nonmodal disturbance growth in boundary-layer flows. This paper investigates the optimal growth characteristics in the hypersonic Mach number regime without any high-enthalpy effects. The influence of wall cooling is studied, with particular emphasis on the role of the initial disturbance location and the value of the spanwise wave number that leads to the maximum energy growth up to a specified location. Unlike previous predictions that used a basic state obtained from a self-similar solution to the boundary-layer equations, mean flow solutions based on the full Navier-Stokes equations are used in select cases to help account for the viscous- inviscid interaction near the leading edge of the plate and for the weak shock wave emanating from that region. Using the full Navier-Stokes mean flow is shown to result in further reduction with Mach number in the magnitude of optimal growth relative to the predictions based on the self-similar approximation to the base flow.

Navier Stokes equations↗
NASA NTRS1992

A computational examination of directional stability for smooth and chined forebodies at high-alpha

Computational Fluid Dynamics (CFD) has been used to study aircraft forebody flowfields at low-speed, angle-of-attack conditions with sideslip. The purpose is to define forebody geometries which provide good directional stability characteristics under these conditions. The flows over the experimentally investigated F-5A forebody and chine type configuration, previously computed by the authors, were recomputed with better grid topology and resolution. The results were obtained using a modified version of CFL3D (developed at NASA Langley) to solve either the Euler equations or the Reynolds equations employing the Baldwin-Lomax turbulence model with the Degani-Schiff modification to account for massive crossflow separation. Based on the results, it is concluded that current CFD methods can be used to investigate the aerodynamic characteristics of forebodies to achieve desirable high angle-of-attack characteristics. An analytically defined generic forebody model is described, and a parametric study of various forebody shapes was then conducted to determine which shapes promote a positive contribution to directional stability at high angle-of-attack. An unconventional approach for presenting the results is used to illustrate how the positive contribution arises. Based on the results of this initial parametric study, some guidelines for aerodynamic design to promote positive directional stability are presented.

Ravi, Ramakrishnan↗
DOE OSTI2021

A scalable exponential-DG approach for nonlinear conservation laws: With application to Burger and Euler equations

In this work, we propose an Exponential DG framework for partial differential equations. We decompose 7 governing equations into linear and nonlinear parts to which we apply the discontinuous Galerkin 8 (DG) spatial discretization. In particular, we construct the linear part using Jacobian that effectively 9 capture stiff characteristics in the system. The former is integrated analytically, whereas the latter 10 is approximated. This approach i) is stable with a large Courant number (Cr > 1); ii) supports 11 high-order solutions both in time and space; iii) is computationally favorable compared to IMEX 12 DG methods with no preconditioner; iv) becomes comparable to explicit RKDG methods on uniform 13 mesh and beneficial on non-uniform grid for Euler equations; v) is scalable in a modern massively 14 parallel computing architecture due to its explicit nature of exponential time integrators and com15 pact communication stencil of DG method. Numerical results demonstrate the performance of our 16 proposed methods through various examples. We also discuss the stability and convergence analysis 17 for our exponential DG scheme in the context of Burgers equation.

42 ENGINEERING↗
DOE OSTI2020

Beam impedance and heating analysis of the diagnostic stripline

In this work the heat transfer analysis has been performed for the stripline kicker, design of which is used at many accelerator facilities as a part of the transverse bunch-by-bunch feed-back system, including NSLS-II. The experimental data of the electrode temperature have been collected from a special diagnostic stripline designed with six infrared view ports. The temperature analysis is based on the power loss calculations and the loss factor computed by the 3-D electromagnetic simulation code GdfidL. To estimate the power loss, we analyze the real part of the longitudinal impedance and the loss factor for two different regimes. We use the ANSYS finite element code to correlate the temperature with the power loss. Lambertson and Shafer formalisms of the longitudinal and transverse beam impedances are discussed for the simplified stripline geometries. The required characteristic impedances and the geometric factors are determined analytically for the simplified geometries and compared with the results for the designed stripline geometry using the 2-D POISSON code by solving the Laplace’s equation in two dimensions. Lambertson’s equations are compared with numerical results. We introduce the frequency range, within which we can apply the analytical approach.

43 PARTICLE ACCELERATORS↗
NASA NTRS1980

Analytical and experimental study of mean flow and turbulence characteristics inside the passages of an axial flow inducer

The inviscid and viscid effects existing within the passages of a three bladed axial flow inducer operating at a flow coefficient of 0.065 are investigated. The blade static pressure and blade limiting streamline angle distributions were determined and the three components of mean velocity, turbulence intensities, and turbulence stresses were measured at locations inside the inducer blade passage utilizing a rotating three sensor hotwire probe. Applicable equations were derived for the hotwire data reduction analysis and solved numerically to obtain the appropriate flow parameters. The three dimensional inviscid flow in the inducer was predicted by numerically solving the exact equations of motion, and the three dimensional viscid flow was predicted by incorporating the dominant viscous terms into the exact equations. The analytical results are compared with the experimental measurements and design values where appropriate. Radial velocities are found to be of the same order as axial velocities within the inducer passage, confirming the highly three dimensional characteristic of inducer flow. Total relative velocity distribution indicate a substantial velocity deficiency near the tip at mid-passage which expands significantly as the flow proceeds toward the inducer trailing edge. High turbulence intensities and turbulence stresses are concentrated within this core region. Considerable wake diffusion occurs immediately downstream of the inducer trailing edge to decay this loss core. Evidence of boundary layer interactions, blade blockage effects, radially inward flows, annulus wall effects, and backflows are all found to exist within the long, narrow passages of the inducer.

Gorton, C. A.↗
NASA NTRS1992

Development of iterative techniques for the solution of unsteady compressible viscous flows

The development of efficient iterative solution methods for the numerical solution of two- and three-dimensional compressible Navier-Stokes equations is discussed. Iterative time marching methods have several advantages over classical multi-step explicit time marching schemes, and non-iterative implicit time marching schemes. Iterative schemes have better stability characteristics than non-iterative explicit and implicit schemes. In this work, another approach based on the classical conjugate gradient method, known as the Generalized Minimum Residual (GMRES) algorithm is investigated. The GMRES algorithm has been used in the past by a number of researchers for solving steady viscous and inviscid flow problems. Here, we investigate the suitability of this algorithm for solving the system of non-linear equations that arise in unsteady Navier-Stokes solvers at each time step.

Sankar, Lakshmi N.↗
NASA NTRS1991

Simulating boundary layer transition with low-Reynolds-number k-epsilon turbulence models. I - An evaluation of prediction characteristics. II - An approach to improving the predictions

The capability of two k-epsilon low-Reynolds number (LRN) turbulence models, those of Jones and Launder (1972) and Lam and Bremhorst (1981), to predict transition in external boundary-layer flows subject to free-stream turbulence is analyzed. Both models correctly predict the basic qualitative aspects of boundary-layer transition with free stream turbulence, but for calculations started at low values of certain defined Reynolds numbers, the transition is generally predicted at unrealistically early locations. Also, the methods predict transition lengths significantly shorter than those found experimentally. An approach to overcoming these deficiencies without abandoning the basic LRN k-epsilon framework is developed. This approach limits the production term in the turbulent kinetic energy equation and is based on a simple stability criterion. It is correlated to the free-stream turbulence value. The modification is shown to improve the qualitative and quantitative characteristics of the transition predictions.

Schmidt, R. C.↗
NASA NTRS1975

Propagation of acoustic modes in the transitional ionosphere

Continuum theory is unsatisfactory for describing wave propagation in the transitional ionosphere, since the mean free path is of the same order as the characteristic length (i.e., Knudsen number in this region is of the order of unity). Therefore a transition model in which isotropic electrons behave as a fluid while the ions are governed by the kinetic equation modified by ion-electron collision effects is proposed to study acoustic wave propagation in the upper atmosphere in the altitude region from 500 to 2000 km. The results show that the dissipation of acoustic waves by electrons is through the viscous and thermal conduction effects, and the dissipation of acoustic waves by ions is through Landau damping. Comparisons of the characteristics of acoustic mode propagation and wave dissipation mechanisms in the collision-dominated, transitional from collisional to collisionless, and collisionless media are also discussed.

Hung, R. J.↗
DOE OSTI2025

Molecular fluctuations inhibit intermittency in compressible turbulence

In the standard picture of fully developed turbulence, highly intermittent hydrodynamic fields are nonlinearly coupled across scales, where local energy cascades from large scales into dissipative vortices and large density gradients. Microscopically, however, constituent fluid molecules are in constant thermal (Brownian) motion, but the role of molecular fluctuations in large-scale turbulence is largely unknown, and with rare exceptions, it has historically been considered irrelevant at scales larger than the molecular mean free path. Recent theoretical and computational investigations have shown that molecular fluctuations can impact energy cascade at Kolmogorov length scales. Here, we show that molecular fluctuations not only modify energy spectrum at wavelengths larger than the Kolmogorov length in compressible turbulence, but also significantly inhibit spatio-temporal intermittency across the entire dissipation range. Using large-scale direct numerical simulations of computational fluctuating hydrodynamics, we demonstrate that the extreme intermittency characteristic of turbulence models is replaced by nearly Gaussian statistics in the dissipation range. These results demonstrate that the compressible Navier–Stokes equations should be augmented with molecular fluctuations to accurately predict turbulence statistics across the dissipation range. Our findings have significant consequences for turbulence modelling in applications such as astrophysics, reactive flows and hypersonic aerodynamics, where dissipation-range turbulence is approximated by closure models.

compressible turbulence↗
DOE OSTI2025

Non-Hermitian quantum mechanics approach for extracting and emulating continuum physics based on bound-state-like calculations: Detailed description

Here, this work applies a reduced basis method to study the continuum physics of a finite quantum system—either few or many-body. Specifically, I develop reduced-order models, or emulators, for the underlying inhomogeneous Schrödinger equation and train the emulators against the equation's bound-state-like solutions at complex energies. The emulators rapidly and accurately interpolate and extrapolate the matrix elements of the Hamiltonian resolvent operator (Green's function) across a parameter space that includes both complex energy and other real-valued physical inputs in the Schrödinger equation. The spectra, discretized and compressed as the result of emulation, and the associated resolvent matrix elements (or amplitudes), have the defining characteristics of non-Hermitian quantum mechanics calculations, featuring complex eigenenergies with negative imaginary parts and branch cuts moved below the real axis in the complex energy plane. Therefore, one now has a method that extracts continuum physics from bound-state-like calculations and emulates those extractions in the input parameter space. Building on a prior Letter [Zhang, Phys. Rev. Lett. 135, 242501 (2025)], this article provides the full theoretical details, a comprehensive analysis of the method's performance, and a brief discussion of how it can be coupled with existing continuum approaches to perform emulations in their input parameter spaces.

ab initio calculations↗
DOE OSTI2023

Numerical Methods to Evaluate Hyperelastic Transducers: Hexagonal Distributed Embedded Energy Converters

Hexagonal distributed embedded energy converters, also known as hexDEECs, are centimeter-scale energy transducers that leverage variable capacitance to generate electricity when their hyperelastic structure is dynamically deformed. To better understand, characterize, and optimize hexDEEC designs, a series of numerical methods and techniques were developed to model the hyperelastic mechanics of hexDEECs, electrostatic properties, and electricity generation characteristics. The numerical methods developed for the hyperelastic structural analysis were corroborated by empirical results from another study, and the models and equations for capacitance, electrostatic forces, and electrical potential energy were derived from fundamental electrostatic equations. These methods and techniques were implemented within the STAR-CCM+ multiphysics software Version 2020.3 (15.06.008) environment. Results from this analysis revealed methodologies and techniques necessary to model the energy converters, which will enable future exploration and optimization of more specific designs and corresponding applications.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1977

Application of the TRPLT1 element to large amplitude free vibrations of plates

A finite element formulation is developed for analyzing large amplitude free flexural vibrations of thin plates in NASTRAN. Stress distributions in the plate, in addition to deflection shapes and nonlinear frequencies are determined. Linearized equations of motion governing large amplitude oscillations of plates and a linearized geometrical stiffness matrix are presented. The solution procedure and convergence characteristics are discussed. The quasi-linear geometrical stiffness matrix for an eighteen degree-of-freedom higher order triangular plate element is evaluated by using a seven-point numerical integration. Nonlinear frequencies for square, rectangular, rhombic, and isosceles triangular plates, with edges simply supported or clamped, are compared with earlier solutions. The present formulation is found to give results entirely adequate for engineering purposes.

Mei, C.↗
NASA NTRS1979

Random harmonic analysis program, L221 (TEV156). Volume 2: Supplemental system design and maintenenace document

Volume 2 of a two volume document is presented. A computer program, L222 (TEV 156), available for execution on the CDC 6600 computer is described. The program is capable of calculating steady-state solutions for linear second-order differential equations due to sinusoidal forcing functions. From this, steady-state solutions, generalized coordinates, and load frequency responses may be determined. Statistical characteristics of loads for the forcing function spectral shape may also be calculated using random harmonic analysis techniques. The particular field of application of the program is the analysis of airplane response and loads due to continuous random air turbulence.

Graham, M. L.↗
NASA NTRS1980

The stratospheric sulfate aerosol layer - Processes, models, observations, and simulations

After briefly reviewing the observational data on the stratospheric sulfate aerosol layer, the chemical and physical processes that are likely to fix the properties of the layer are discussed. We present appropriate continuity equations for aerosol particles, and show how to solve the equations on a digital computer. Simulations of the unperturbed aerosol layer by various published models are discussed and the sensitivity of layer characteristics to variations in several aerosol model parameters is studied. We discuss model applications to anthropogenic pollution problems and demonstrate that moderate levels of aerospace activity (supersonic transport and Space Shuttle operations) will probably have only a negligible effect on global climate. Finally, we evaluate the possible climatic effect of a ten-fold increase in the atmospheric abundance of carbonyl sulfide.

Whitten, R. C.↗