Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical calculations: variational”

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 271 records · Page 15

Hot gas environment around STOVL aircraft in ground proximity. II - Numerical study

Ingestion of hot exhaust gases by the engines of STOVL aircraft has been an important research problem for several years. The hot gas environment around STOVL aircraft is three-dimensional and turbulent. In this study, the Navier-Stokes equations governing the hot gas ingestion flowfield are solved by an efficient finite-difference calculation procedure. The complete geometry including the head wind and the fuselage is simulated. Four demonstration calculations with variations in the height of the fuselage and the head wind velocity are presented. It is shown that the calculation procedure efficiently provides a solution to the governing equations and produces realistic descriptions of the flow and temperature fields.

Tafti, D. K.↗

Non-LTE profiles of the Al I autoionization lines

A non-LTE formulation is given for the transfer of radiation in the autoionizing lines of neutral aluminum at 1932 and 1936 A through both the Bilderberg and Harvard-Smithsonian model atmospheres. Numerical solutions for the common source function of these lines and their theoretical line profiles are calculated and compared with the corresponding LTE profiles. The results show that the non-LTE profiles provide a better match with the observations. They also indicate that the continuous opacity of the standard solar models should be increased in this wavelength region if the center-limb variations of observed and theoretical profiles of these lines are to be in reasonable agreement.

Finn, G. D.↗

Flap-lag-torsion aeroelastic stability of a circulation control rotor in forward flight

The aeroelastic stability of a circulation control rotor blade undergoing three degrees of motion (flap, lag, and torsion) is investigated in forward flight. Quasi-steady strip theory is used to evaluate the aerodynamics forces; and the airfoil characteristics are from data tables. The propulsive and the auxiliary power trims are calculated from vehicle and rotor equilibrium equations through the numerical integration of element forces in azimuth as well as in radial directions. The nonlinear time dependent periodic blade response is calculated using an iterative procedure based on Floquet theory. The periodic perturbation equations are solved for stability using Floquet transition matrix theory. The effects of several parameters on blade stability are examined, including advance ratio, collective pitch, thrust level, shaft tilt, structural stiffnesses variation, and propulsive and auxiliary power trims.

Chopra, Inderjit↗

Critical evaluation of Jet-A spray combustion using propane chemical kinetics in gas turbine combustion simulated by KIVA-II

Numerical solutions of the Jet-A spray combustion were obtained by means of the KIVA-II computer code after Jet-A properties were added to the 12 chemical species the program had initially contained. Three different reaction mechanism models are considered. The first model consists of 131 reactions and 45 species; it is evaluated by comparing calculated ignition delay times with available shock tube data, and it is used in the evaluation of the other two simplified models. The simplified mechanisms consider 45 reactions and 27 species and 5 reactions and 12 species, respectively. In the prediction of pollutants NOx and CO, the full mechanism of 131 reactions is considered to be more reliable. The numerical results indicate that the variation of the maximum flame temperature is within 20 percent as compared with that of the full mechanism of 131 reactions. The chemical compositions of major components such as C3H8, H2O, O2, CO2, and N2 are of the same order of magnitude. However, the concentrations of pollutants are quite different.

Nguyen, H. L.↗

Evaluation of F/A-18A HARV inlet flow analysis with flight data

The F/A-18A aircraft has experienced engine stalls at high angles-of-attack and yaw flight conditions which were outside of its flight envelope. Future aircraft may be designed to operate routinely in this flight regime. Therefore, it is essential that an understanding of the inlet flow field at these flight conditions be obtained. Due to the complex interactions of the fuselage and inlet flow fields, a study of the flow within the inlet must also include external effects. Full Navier-Stokes (FNS) calculations on the F/A-18A High Alpha Research Vehicle (HARV) inlet for several angles-of-attack with sideslip and free stream Mach numbers have been obtained. The predicted forebody/fuselage surface static pressures agreed well with flight data. The surface static pressures along the inlet lip are in good agreement with the numerical predictions. The major departure in agreement is along the bottom of the lip at 30 deg and 60 deg angle-of-attack where a possible streamwise flow separation is not being predicted by the code. The circumferential pressure distributions at the engine face are in very good agreement with the numerical results. The variation in surface static pressure in the circumferential direction is very small with the exception of 60 angle-of-attack. Although the simulation does not include the effect of the engine, it appears that this omission has a second order effect on the circumferential pressure distribution. An examination of the unsteady flight test data base has shown that the secondary vortex migrates a significant distance with time. In fact, the extent of this migration increases with angle-of-attack with increasing levels of distortion. The effects of the engine on this vortex movement is unknown. This implies that the level of flow unsteadiness increases with increasing distortion. Since the computational results represent an asymptotic solution driven by steady boundary conditions, these numerical results may represent an arbitrary point in time. A comparison of the predicted total pressure contours with flight data indicates that the numerical results are within the excursion range of the unsteady data which is the best the calculations can attain unless an unsteady simulation is performed.

Smith, C. Frederic↗

Venus. I - Carbon monoxide distribution and molecular-line searches

An observational program to study variations of the vertical distribution of CO in the Venus atmosphere is presented. Measurements of the J = 0 - 1 absorption line at 2.6 mm wavelength are reported for two phase angles in 1977, one near eastern elongation (February) and the other near inferior conjunction (April). The two spectra are significantly different, with the April absorption line being narrower and deeper. The results of numerical inversion calculations show that the CO mixing ratio increases by a factor of approximately 100 between 78 and 100 km and that the CO abundance above approximately 100 km is greatest on the night-side hemisphere. These conclusions are in qualitative agreement with theoretical models. In addition to the CO observations, a search for other molecules was made to provide further information on the composition of the Venus middle atmosphere. The J = 0 - 1 transition of (C-13)O was detected and upper limits were derived for nine other molecules.

Wilson, W. J.↗

Planetary Accretion, Oxygen Isotopes and the Central Limit Theorem

The accumulation of presolar dust into increasingly larger aggregates (CAIs and Chondrules, Asteroids, Planets) should result in a very drastic reduction in the numerical spread in oxygen isotopic composition between bodies of similar size, in accord with the Central Limit Theorem. Observed variations in oxygen isotopic composition are many orders of magnitude larger than would be predicted by a simple, random accumulation model that begins in a well-mixed nebula - no matter which size-scale objects are used as the beginning or end points of the calculation. This discrepancy implies either that some as yet unspecified process acted on the solids in the Solar Nebula to increase the spread in oxygen isotopic composition during each and every stage of accumulation or that the nebula was heterogeneous and maintained this heterogeneity throughout most of nebular history. Large-scale nebular heterogeneity would have significant consequences for many areas of cosmochemistry, including the application of some well-known isotopic systems to the dating of nebular events or the prediction of bulk compositions of planetary bodies on the basis of a uniform cosmic abundance.

Nuth, Joseph A., III↗

A Look at the Truths and Misconceptions of the Variational Quantum Eigensolver and the Implications of Overparameterization

In this work, we investigate loss landscapes of the variational quantum eigensolver (VQE) by quantifying the number of local minima through empirical analyses. We focus on minimal models in chemistry and physics so that we can do a complete analysis using more computationally expensive tools. We employ Hessian eigenvalue calculations and the nudged elastic band algorithm to characterize these landscapes. Our results expand upon the existing literature by highlighting the optimization challenges faced by VQE. We find that, as the number of parameters in our ansatz increases, the number of basins increases while the corresponding loss function values converge toward the global minimum value. This observation implies that overparameterization may lead to an ``effective convexity'' in VQE loss landscapes, a phenomenon supported by theoretical and numerical work in classical machine learning.

quantum computing↗

Numerical calculation of boundary-induced interference in slotted or perforated wind tunnels including viscous effects in slots

A numerical method is presented for calculating the incompressible boundary-induced interference in wind tunnels of rectangular cross section with slotted or perforated walls. The method includes a wall representation which is capable of satisfying a generalized homogeneous boundary condition including the effects of viscosity within the slots. The effects of viscosity in the slots are found to be very significant. The method allows for a variation in the boundary conditions along the tunnel walls. The model can be any configuration and can be located anywhere in the test section. The interference can be computed at any point in the test section.

Keller, J. D.↗

A time variable model of Earth's albedo

A time variable model of Earth's albedo was prepared for use in climate studies and as an aid to the interpretation of satellite Earth radiation budget data. The features of the model include: a 10 deg latitude 10 deg longitude grid for numerical integration, surface albedo specified at 1 month intervals, calculation of zenith angle effect for surface albedo and of the additional effect of the atmosphere on the albedo. Percent cloud cover is specified for 29 different climatological cloud type regions at 8 times of the day for 12 months of the year. Cloud albedos were specified for each of the cloud climatological types. Diurnal and monthly variations of this model are described and results are compared with a model which is based on satellite measurements. A computer program was also written for use in studying the sampling effects in satellite radiation budget measurements. An example of the results of calculations with this program are compared with a previous study of the sampling effects. This program for satellite orbit characteristics is to be combined with the time-variable albedo model for further study of the sampling problem.

Bartman, F. L.↗

The disturbance flow field produced by an evolving vortex

The flow field of a vortex in a viscous shear flow is found by constructing a uniformly valid asymptotic expansion consisting of an inner solution field represented, to lowest order, by a two dimensional, nonliner, inviscid Stuart vortex and an outer solution field represented, to lowest order, by either a two dimensional parallel or self similar viscous flow. The technique involves scaling both the transverse and streamwise coordinates in the vicinity of the vortex as well as allowing for a slow variation of the outer viscous flow. Criteria are established for both the size of the vortical structure and proximity to the boundary surfaces. The composite solution is a consistent mathematical picture of the flow field at a fixed streamwise location as the vortical structure evolves past this point. Such a formulation is also useful in the specification of boundary or initial conditions in numerical fluid dynamic calculations, where an inconsistent setting of these conditions leads to spurious results for rather long computation times.

Gatski, T. B.↗

Acoustic radiation from lifting airfoils in compressible subsonic flow

The far field acoustic radiation from a lifting airfoil in a three-dimensional gust is studied. The acoustic pressure is calculated using the Kirchhoff method, instead of using the classical acoustic analogy approach due to Lighthill. The pressure on the Kirchhoff surface is calculated using an existing numerical solution of the unsteady flow field. The far field acoustic pressure is calculated in terms of these values using Kirchhoff's formula. The method is validated against existing semi-analytical results for a flat plate. The method is then used to study the problem of an airfoil in a harmonic three-dimensional gust, for a wide range of Mach numbers. The effect of variation of the airfoil thickness and angle of attack on the acoustic far field is studied. The changes in the mechanism of sound generation and propagation due to the presence of steady loading and nonuniform mean flow are also studied.

Atassi, Hafiz M.↗

Acoustic radiation from lifting airfoils in compressible subsonic flow

The far field acoustic radiation from a lifting airfoil in a three-dimensional gust is studied. The acoustic pressure is calculated using the Kirchhoff method, instead of using the classical acoustic analogy approach due to Lighthill. The pressure on the Kirchhoff surface is calculated using an existing numerical solution of the unsteady flow field. The far field acoustic pressure is calculated in terms of these values using Kirchhoff's formula. The method is validated against existing semi-analytical results for a flat plate. The method is then used to study the problem of an airfoil in a harmonic three-dimensional gust, for a wide range of Mach numbers. The effect of variation of the airfoil thickness and angle of attack on the acoustic far field is studied. The changes in the mechanism of sound generation and propagation due to the presence of steady loading and non-uniform mean flow are also studied.

Atassi, Hafiz M.↗

Three-dimensional elastic-plastic finite-element analyses of constraint variations in cracked bodies

Three-dimensional elastic-plastic (small-strain) finite-element analyses were used to study the stresses, deformations, and constraint variations around a straight-through crack in finite-thickness plates for an elastic-perfectly plastic material under monotonic and cyclic loading. Middle-crack tension specimens were analyzed for thicknesses ranging from 1.25 to 20 mm with various crack lengths. Three local constraint parameters, related to the normal, tangential, and hydrostatic stresses, showed similar variations along the crack front for a given thickness and applied stress level. Numerical analyses indicated that cyclic stress history and crack growth reduced the local constraint parameters in the interior of a plate, especially at high applied stress levels. A global constraint factor alpha(sub g) was defined to simulate three-dimensional effects in two-dimensional crack analyses. The global constraint factor was calculated as an average through-the-thickness value over the crack-front plastic region. Values of alpha(sub g) were found to be nearly independent of crack length and were related to the stress-intensity factor for a given thickness.

Newman, J. C., Jr.↗

An asymptotic theory of supersonic propeller noise

A theory for predicting the noise field of a propeller with a realistic blade geometry is presented. The theory, which utilizes a large blade count approximation, provides an efficient formula for predicting the radiation of sound from all three sources of propeller noise. Comparisons with full numerical integration indicate that the noise levels predicted by this formula are quite accurate. Calculations based on this method also show that the radiation from the Lighthill quadrupole source is rather substantial when compared with thickness and loading noise for high speed propellers. A preliminary application of the theory to the problem of the sensitivity of the peak noise levels generated by a supersonic propeller to the variations in its tip helical Mach number has produced a trend that is in qualitative agreement with the experimental observations.

Envia, Edmane↗

Effects of wing flexibility and variable air lift upon wing bending moment during landing impacts of a small seaplane

A smooth-water-landing investigation was conducted with a small seaplane to obtain experimental wing-bending-moment time histories together with time histories of the various parameters necessary for the prediction of wing bending moments during hydrodynamic forcing functions. The experimental results were compared with calculated results which include inertia-load effects and the effects of air-load variation during impact. The responses of the fundamental mode were calculated with the use of the measured hydrodynamic forcing functions. From these responses, the wing bending moments due to the hydrodynamic load were calculated according to the procedure given in R.M. No. 2221. The comparison of the time histories of the experimental and calculated wing bending moments showed good agreement both in phase relationship of the oscillations and in numerical values.

Merten, Kenneth F↗

Numerical optimization techniques for bound circulation distribution for minimum induced drag of nonplanar wings: Basic formulations

A theoretical method was developed for determining the optimum span load distribution for minimum induced drag for subsonic nonplanar configurations. The undistorted wing wake is assumed to have piecewise linear variation of shed vortex sheet strength, resulting in a quadratic variation of bound circulation and span load. The optimum loading is obtained either through a direct technique, whereby derivatives of the drag expression are calculated analytically in terms of the unknown wake vortex sheet strengths. Both techniques agree well with each other and with available exact solutions for minimum induced drag.

Kuhlman, J.↗

Numerical calculation of finite amplitude effects in unstable laminar boundary layers

A numerical study of one stage of the transition region from laminar to turbulent flow of the boundary layer over a flat plate is performed. Benney-Lin theory is applied to the laminar profile, which is closer to that of the experiment, i.e., the Blasius profile. Several laminar profiles are studied, all belonging to the Falkner-Skan family. The flow is found to depend principally upon the disturbance spanwise wavenumber, beta, whose variation results, in some instances, in a system of counter-rotating vortices, one on top of the other, for the mean secondary flow.

Antar, B. N.↗