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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

A conservative type-dependent full potential method for the treatment of supersonic flows with embedded subsonic regions

A nonlinear method based on the full potential equation in conservation form, cast in an arbitrary coordinate system, has been developed to treat predominantly supersonic flows with embedded subsonic regions. This type of flow field occurs frequently near the fuselage-canopy junction area and wing leading edge regions for a moderately swept fighter configuration. The method uses the theory of characteristics to accurately monitor the type-dependent flow field. A conservative switching scheme is developed to transition from the supersonic marching algorithm to a subsonic relaxation procedure, and vice versa. An implicit approximate factorization scheme is employed to solve the finite-differenced equation. Results are shown for a few configurations, including a wing-body-wake realistic fighter model having embedded subsonic regions.

Shankar, V.↗
NASA NTRS2005

Incorporation of a Chemical Equilibrium Equation of State into LOCI-Chem

Renewed interest in development of advanced high-speed transport, reentry vehicles and propulsion systems has led to a resurgence of research into high speed aerodynamics. As this flow regime is typically dominated by hot reacting gaseous flow, efficient models for the characteristic chemical activity are necessary for accurate and cost effective analysis and design of aerodynamic vehicles that transit this regime. The LOCI-Chem code recently developed by Ed Luke at Mississippi State University for NASA/MSFC and used by NASA/MSFC and SSC represents an important step in providing an accurate, efficient computational tool for the simulation of reacting flows through the use of finite-rate kinetics [3]. Finite rate chemistry however, requires the solution of an additional N-1 species mass conservation equations with source terms involving reaction kinetics that are not fully understood. In the equilibrium limit, where the reaction rates approach infinity, these equations become very stiff. Through the use of the assumption of local chemical equilibrium the set of governing equations is reduced back to the usual gas dynamic equations, and thus requires less computation, while still allowing for the inclusion of reacting flow phenomenology. The incorporation of a chemical equilibrium equation of state module into the LOCI-Chem code was the primary objective of the current research. The major goals of the project were: (1) the development of a chemical equilibrium composition solver, and (2) the incorporation of chemical equilibrium solver into LOCI-Chem. Due to time and resource constraints, code optimization was not considered unless it was important to the proper functioning of the code.

Cox, Carey F.↗
NASA NTRS1993

CFD validation of subsonic turbulent planar shear layers

The primary objective of the present study is to assess the limitations and capabilities of RPLUS in predicting the entrainment, mixing, and burning characteristics of a high subsonic planar shear layer; these computations are performed in conjunction with an experiment being conducted at NASA Lewis Research Center. Turbulence is modeled by a two-equation k-epsilon closure modified for compressibility effects. Combustion is modeled by two different models: one in which a finite-rate laminar chemistry model is used for H2-O2 oxidation, and the other based on a composition joint pdf approach to account for the turbulence-chemistry interaction effects. Between upwind and centered differencing schemes, minimal differences are observed for the flow properties of a nonreacting shear layer case studied. In comparison with experimental data, computed growth rates are underpredicted in both reacting and nonreacting cases of a shear layer, but more severely for the reacting flows. However, the qualitative agreement between computation and experiment is reasonably good for the other observable characteristics, and indicative of potential for further improvements.

Lai, H. T.↗
NASA NTRS2015

A Persistent Feature of Multiple Scattering of Waves in the Time-Domain: A Tutorial

The equations for frequency-domain multiple scattering are derived for a scalar or electromagnetic plane wave incident on a collection of particles at known positions, and in the time-domain for a plane wave pulse incident on the same collection of particles. The calculation is carried out for five different combinations of wave types and particle types of increasing geometrical complexity. The results are used to illustrate and discuss a number of physical and mathematical characteristics of multiple scattering in the frequency- and time-domains. We argue that frequency-domain multiple scattering is a purely mathematical construct since there is no temporal sequencing information in the frequency-domain equations and since the multi-particle path information can be dispelled by writing the equations in another mathematical form. However, multiple scattering becomes a definite physical phenomenon in the time-domain when the collection of particles is illuminated by an appropriately short localized pulse.

particles↗
NASA NTRS1977

Flow study in the cross sectional planes of a turbine scroll

A numerical study of the nonviscous flow characteristics in the cross-sectional planes of a radial inflow turbine scroll is presented. The velocity potential is used in the formulation to determine the flow velocity in these planes resulting from the continuous mass discharge. The effect of the through flow velocity is simulated by a continuous distribution of source/sink in the cross-section. A special iterative procedure is devised to handle the solution of the resulting Poisson's differential equation with Neumann boundary conditions in a domain with generally curved boundaries. The analysis is used to determine the effects of the radius of curvature, the location of the scroll section and its geometry on the flow characteristics in the turbine scroll.

Hamed, A.↗
NASA NTRS1989

Light scattering characteristics of various aerosol types derived from multiple wavelength lidar observations

The present study demonstrates the potential of a multiple-wavelength lidar for discriminating between several aerosol types on the basis of the wavelength dependence of the aerosol backscatter coefficient. The two-component lidar equation was solved under the assumption of similarity in the derived profiles of backscatter coefficients for each wavelength. It is shown that a three-wavelength lidar system operating at 300, 600, and 1064nm can provide unique information for discriminating between various aerosol types (continental, maritime, Saharan-dust, stratospheric aerosols in a tropopause fold event, and tropical forest aerosols). Mie calculations were made using in situ aerosol data and aerosol models to compare with the lidar results. The disagreement between the theoretical and empirical results in some cases was substantial. These differences may be partly due to uncertainties in the lidar data analysis and aerosol characteristics and also due to the conventional assumption of aerosol sphericity for the aerosol Mie calculations.

Sasano, Yasuhiro↗
DOE OSTI2021

Mixing Cell Model: A One-Dimensional Numerical Model for Assessment of Water Flow and Contaminant Transport in the Unsaturated Zone

This report describes the Mixing Cell Model code, a one dimensional model for water flow and solute transport in the unsaturated zone under steady state or transient flow conditions. The model is based on the principles and assumptions underlying mixing-cell model formulations. The unsaturated zone is discretized into a series of independent mixing cells. Each cell may have unique hydrologic, lithologic, and sorptive properties. Ordinary differential equations describe the material (water and solute) balance within each cell. Water-flow equations are derived from the continuity equation, assuming that unit gradient conditions exist at all times in each cell. Pressure gradients are considered implicitly through model discretization. Unsaturated hydraulic conductivity and moisture contents are determined by the material specific moisture-characteristic curves. Solute-transport processes include explicit treatment of advective processes, first order chain decay, and linear sorption reactions. Dispersion is addressed through implicit and explicit dispersion. Implicit dispersion is an inherent feature of all mixing-cell models and originates from the formulation of the problem in terms of mass balance around fully mixed volume elements. Expressions are provided that relate implicit dispersion to the physical dispersion of the system. Two FORTRAN codes were developed to solve the water flow and solute-transport equations: (1) the Mixing Cell Model for Flow (MCMF) solves transient water-flow problems and (2) the Mixing Cell Model for Transport (MCMT) solves the solute-transport problem. The transient water-flow problem is typically solved first by estimating the water flux through each cell in the model domain as a function of time using the MCMF code. These data are stored in either ASCII or binary files that are later read by the solute transport code (MCMT). Code output includes solute pore water concentrations, water and solute inventories in each cell and at each specified output time, and water and solute fluxes through each cell and specified output time. Computer run times for coupled transient water flow and solute transport were typically several seconds on a 2 GHz Intel Pentium IV desktop computer. The model was benchmarked against analytical solutions and finite element approximations to the partial differential equations (PDE) describing unsaturated flow and transport. Differences between the maximum solute flux estimated by the mixing cell model and the PDE models were typically less than two percent. This revision includes an option for a fixed concentration lower boundary condition for diffusive fluxes for versions 020321 and later.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1998

A Technique for Assimilating GOES-Derived Land Surface Products into Regional Models to Improve the Representation of Land Surface Forcing

As the parameterizations of surface energy budgets in regional models have become more complete physically, models have the potential to be much more realistic in simulations of coupling between surface radiation, hydrology, and surface energy transfer. Realizing the importance of properly specifying the surface energy budget, many institutions are using land-surface models to represent the lower boundary forcing associated with biophysical processes and soil hydrology. However, the added degrees of freedom due to inclusion of such land-surface schemes require the specification of additional parameters within the model system such as vegetative resistances, green vegetation fraction, leaf area index, soil physical and hydraulic characteristics, stream flow, runoff, and the vertical distribution of soil moisture. A technique has been developed for assimilating GOES-IR skin temperature tendencies into the surface energy budget equation of a mesoscale model so that the simulated rate of temperature change closely agrees with the satellite observations. A critical assumption of the technique is that the availability of moisture (either from the soil or vegetation) is the least known term in the model's surface energy budget. Therefore, the simulated latent heat flux, which is a function of surface moisture availability, is adjusted based upon differences between the modeled and satellite-observed skin temperature tendencies. An advantage of this technique is that satellite temperature tendencies are assimilated in an energetically consistent manner that avoids energy imbalances and surface stability problems that arise from direct assimilation of surface shelter temperatures. The fact that the rate of change of the satellite skin temperature is used rather than the absolute temperature means that sensor calibration is not as critical. An advantage of this technique for short-range forecasts is that it does not require a complex land-surface formulation within the atmospheric model. As a result, the need to specify poorly known soil and vegetative characteristics is eliminated. The GOES assimilation technique has been incorporated into the PSU/NCAR MM5. Results will be presented to demonstrate the ability of the assimilation scheme to improve short-term simulations of near-surface air temperature and mixing ratio during the warm season for several selected cases during June of 1995 and June 1997 which exhibit a variety of atmospheric and land-surface conditions. In addition, the simulations produced by the assimilation technique will be compared with those produced by the MM5/BATS coupled model initialized with soil water content derived from a soil hydrology model.

Lapenta, William M.↗
DOE OSTI2022

Physics-Informed Sparse Gaussian Process for Probabilistic Stability Analysis of Large-Scale Power System with Dynamic PVs and Loads

This work proposes a physics-informed sparse Gaussian process (SGP) for probabilistic stability assessment of large-scale power systems in the presence of uncertain dynamic PVs and loads. The differential and algebraic equations considering uncertainties from dynamic PVs and loads are reformulated to a nonlinear mapping relationship that allows the application of SGP. Thanks to the nonparametric characteristic of Gaussian process, the proposed framework does not require distributions of uncertain inputs and this distinguishes it from existing approaches. As the original Gaussian process is not scalable to large-scale systems with high dimensional uncertain inputs, this paper develops the SGP with a stochastic variational inference technique. It leads to approximately two orders of complex reduction. A data pre-processing step is also introduced to tackle the coexistence of stable and unstable cases by sample clustering and constructing separate SGPs. The probabilistic transient stability index is analyzed to assess system stability under different uncertain dynamics loads and PVs. Comparisons are performed with the sampling-based, the polynomial chaos expansion-based, and traditional Gaussian process-based methods on the modified IEEE 118-bus and Texas 2000-bus systems under various scenarios, including different levels of uncertainties and the existence of nonlinear correlations among dynamic PVs. The impacts of data quality and quantity issues are also investigated. It is shown that the proposed SGP achieves significantly improved computational efficiency while maintaining high accuracy with a limited number of data.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Helmholtz vibrations in bowed strings

We report for almost 160 years, it has been known that Helmholtz oscillations, unique to vibrating strings in bowed instruments (violin, cello, etc.), have two distinct regimes: “slip” and “stick.” During the slip regime, the force at the bow-string interaction is attributed to friction between the sliding bow hair and the vibrating string, with a friction coefficient that decreases with increasing relative velocity. Yet the hair-string interaction during the stick regime is less understood. We propose that the interaction force during the stick regime is proportional to the product of the longitudinal acoustic impedance of the bow hair to the relative bow-string velocity. We validate this hypothesis by solving the string's differential equation of motion, including an enhanced formulation to avoid parasitic high-frequency oscillations. This physical model enables us to analyze, in real time, the characteristics of the Helmholtz oscillations, including the string shape, excitation of harmonics, Schelleng ripples, and string energy, showing that the bowed string gains energy during the stick regime and loses energy during the slip regime.

36 MATERIALS SCIENCE↗
NASA NTRS1972

Effect of support flexibilty and damping on the dynamic response of a single mass flexible rotor in elastic bearings

A steady state analysis of the shaft and the bearing housing motion was made by assuming synchronous precession of the system. The conditions under which the support system would act as a dynamic vibration absorber at the rotor critical speed were studied; plots of the rotor and support amplitudes, phase angles, and forces transmitted were evaluated by the computer, and the performance curves were automatically plotted by a CalComp plotter unit. Curves are presented on the optimization of the support housing characteristics to attenuate the rotor unbalance response over the entire rotor speed range. The complete transient motion including rotor unbalance was examined by integrating the equations of motion numerically using a modified fourth order Runge-Kutta procedure, and the resulting whirl orbits were plotted by the CalComp plotter unit. The results of the transient analysis are discussed with regards to the design optimization procedure derived from the steady-state analysis.

Kirk, R. G.↗
NASA NTRS1972

Optimal Quasi-steady Plasma Thruster system characteristics.

The overall characteristics of a generalized Quasi-steady Plasma Thruster (QPT) system consisting of thruster head, power conditioning network, propellant supply subsystem are studied. Energy balance equations for the system are coupled with component mass relationships in order to determine overall system mass and performance. Power supply power levels varying from 100 to 10,000 watts with thruster power levels ranging from 300 kw to 30 Mw employing argon as the propellant are considered. The manner in which overall system mass, average thrust, and burn time vary as a function power supply power level, quasi-steady power level, and pulse time are studied. Results indicate the existence of optimum pulse times when system mass is employed as an optimization criterion.

Ludwig, D. E.↗
NASA NTRS1976

Effects of mean flow on duct mode optimum suppression rates

The nature of the solution to the convected acoustic wave equation and associated boundary conditions for rectangular ducts containing uniform mean flow is examined in terms of the complex mapping between the wall admittance and characteristic mode eigenvalues. It is shown that the Cremer optimum suppression criteria must be modified to account for the effects of flow below certain critical values of the nondimensional frequency parameter of duct height divided by sound wavelength. The implications of these results on the design of low frequency suppressors are considered.

Kraft, R. E.↗
NASA NTRS1978

Measurements of S II optical emission from the thermal plasma of Jupiter

Photometric spectra of the forbidden S II emission lines at 6716 and 6731 A originating in the Jovian magnetosphere are examined. The summed apparent emission rate in the two sulfur lines is plotted against the central-meridian longitude of Jupiter in System III. It is found that the total sulfur brightness exhibited a distinct minimum when the Jovian magnetic dipole was in the plane defined by the earth-Jupiter vector and the Jovian rotational axis. The observations are shown to be qualitatively consistent with an annular emitting region at the magnetic equator with a thickness of about 1 Jupiter radius and an inner radius of 4 to 5 Jupiter radii. The observations also imply a (logarithmic) characteristic electron density of approximately 3.7 per cu cm and a (logarithmic) average column abundance of about 11.5 per sq cm for singly ionized sulfur.

Brown, R. A.↗
NASA NTRS1980

Flutter analysis of an airplane with multiple structural nonlinearities in the control system

Experience has shown that the flutter prediction process for airplanes can be greatly affected by strong concentrated nonlinearities which may be localized in the linking elements of the control mechanism, in the pivot joints of variable-sweep-wing systems, and in the connecting points between wing and pylon-mounted external stores. The principle of equivalent linearization offers an efficent possibility for solving the related nonlinear flutter equations in the frequency domain as a complement to the well-known time domain procedures. Taking as an example an airplane with nonlinear control characteristics, it is demonstrated how the equivalent linearization approach can be extended to rather complicated systems with multiple sets of strongly interacting, concentrated nonlinearities.

Breitbach, E. J.↗
NASA NTRS1980

A diagonal form of an implicit approximate-factorization algorithm with application to a two dimensional inlet

A modification of an implicit approximate-factorization finite-difference algorithm applied to the two dimensional Euler and Navier-Stokes equations in general curvilinear coordinates is presented for supersonic free stream flow about and through inlets. The modification transforms the coupled system of equations into an uncoupled diagonal form which requires less computation work. For steady-state applications the resulting diagonal algorithm retains the stability and accuracy characteristics of the original algorithm. Solutions are given for inviscid and laminar flow about a two dimensional wedge inlet configuration. Comparisons are made between computed results and exact theory.

Chaussee, D. S.↗
NASA NTRS1981

A diagonal form of an implicit approximate-factorization algorithm

A modification of an implicit approximate-factorization finite-difference algorithm applied to partial differential equations is presented. This algorithm is applied to the two- and three-dimensional Euler equations in general curvilinear coordinates. The modification transforms the coupled system of equations into an uncoupled diagonal form that requires less computational work. For steady-state applications, the resulting diagonal algorithm retains the stability and accuracy characteristics of the original algorithm. The diagonal algorithm reduces the storage requirement of the implicit solution process and therefore has an important effect on the application of implicit finite-difference schemes to vector processors. Results are presented for realistic two-dimensional transonic flow fields about airfoils. Computation costs are reduced to 24-34%.

Pulliam, T. H.↗
NASA NTRS1981

A high-resolution study of the gas flow in barred spirals

The considered investigation is concerned with conditions in that type of spiral galaxy which is known as a 'barred spiral'. The lack of star formation and the smoothness of the dark dust lanes may indicate that the bar region is pervaded by a more homogeneous interstellar medium with a lower viscosity than that characteristic of the outer regions. The effects of such a reduced viscosity are investigated on the basis of high-resolution time-dependent gasdynamical equations. Calculations modeling the gas flow in a given bar potential are presented. The calculations were performed on grids with various mesh sizes and consequently had different resolutions and numerical viscosities. A reduction in the numerical diffusion and viscosity was accompanied by changes in the flow pattern, including the appearance of a strong shock in the bar region and the occurrence of postshock outflow.

Van Albada, G. D.↗