Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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 433 records · Page 24

Structural response to discrete and continuous gusts of an airplane having wing bending flexibility and a correlation of calculated and flight results

An analysis is made of the structural response to gusts of an airplane having the degrees of freedom of vertical motion and wing bending flexibility and basic parameters are established. A convenient and accurate numerical solution of the response equations is developed for the case of discrete-gust encounter, an exact solution is made for the simpler case of continuous-sinusoidal-gust encounter, and the procedure is outlined for treating the more realistic condition of continuous random atmospheric turbulence, based on the methods of generalized harmonic analysis. Correlation studies between flight and calculated results are then given to evaluate the influence of wing bending flexibility on the structural response to gusts of two twin-engine transports and one four-engine bomber. It is shown that calculated results obtained by means of a discrete-gust approach reveal the general nature of the flexibility effects and lead to qualitative correlation with flight results. In contrast, calculations by means of the continuous-turbulence approach show good quantitative correlation with flight results and indicate a much greater degree of resolution of the flexibility effects.

Houbolt, John C↗

Continued Development of Expert System Tools for NPSS Engine Diagnostics

The objectives of this grant were to work with previously developed NPSS (Numerical Propulsion System Simulation) tools and enhance their functionality; explore similar AI systems; and work with the High Performance Computing Communication (HPCC) K-12 program. Activities for this reporting period are briefly summarized and a paper addressing the implementation, monitoring and zooming in a distributed jet engine simulation is included as an attachment.

Lewandowski, Henry↗

Maneuver simulations of flexible spacecraft by solving TPBVP

The optimal control of large angle rapid maneuvers and vibrations of a Shuttle mast reflector system is considered. The nonlinear equations of motion are formulated by using Lagrange's formula, with the mast modeled as a continuous beam. The nonlinear terms in the equations come from the coupling between the angular velocities, the modal coordinates, and the modal rates. Pontryagin's Maximum Principle is applied to the slewing problem, to derive the necessary conditions for the optimal controls, which are bounded by given saturation levels. The resulting two point boundary value problem (TPBVP) is then solved by using the quasilinearization algorithm and the method of particular solutions. In the numerical simulations, the structural parameters and the control limits from the Spacecraft Control Lab Experiment (SCOLE) are used. In the 2-D case, only the motion in the plane of an Earth orbit or the single axis slewing motion is discussed. In the 3-D slewing, the mast is modeled as a continuous beam subjected to 3-D deformations. The numerical results for both the linearized system and the nonlinear system are presented to compare the differences in their time response.

Bainum, Peter M.↗

3-D Flow Field Diagnostics and Validation Studies using Stereoscopic Tracking Velocimetry

The measurement of 3-D three-component velocity fields is of great importance in both ground and space experiments for understanding materials processing and fluid physics. Here, we present the investigation results of stereoscopic tracking velocimetry (STV) for measuring 3-D velocity fields. The effort includes diagnostic technology development, experimental velocity measurement, and comparison with analytical and numerical computation. The advantages of STV stems from the system simplicity for building compact hardware and in software efficiency for continual near-real-time process monitoring. It also has illumination flexibility for observing volumetric flow fields from arbitrary directions. STV is based on stereoscopic CCD observations of particles seeded in a flow. Neural networks are used for data analysis. The developed diagnostic tool is tested with a simple directional solidification apparatus using Succinonitrile. The 3-D velocity field in the liquid phase is measured and compared with results from detailed numerical computations. Our theoretical, numerical, and experimental effort has shown STV to be a viable candidate for reliably quantifying the 3-D flow field in materials processing and fluids experiments.

Cha, Soyoung Stephen↗

Three-Dimensional Stereoscopic Tracking Velocimetry and Experimental/Numerical Comparison of Directional Solidification

Measurement of three-dimensional (3-D) three-component velocity fields is of great importance in both ground and space experiments for understanding materials processing and fluid physics. The experiments in these fields most likely inhibit the application of conventional planar probes for observing 3-D phenomena. Here, we present the investigation results of stereoscopic tracking velocimetry (STV) for measuring 3-D velocity fields, which include diagnostic technology development, experimental velocity measurement, and comparison with analytical and numerical computation. STV is advantageous in system simplicity for building compact hardware and in software efficiency for continual near-real-time monitoring. It has great freedom in illuminating and observing volumetric fields from arbitrary directions. STV is based on stereoscopic observation of particles-Seeded in a flow by CCD sensors. In the approach, part of the individual particle images that provide data points is likely to be lost or cause errors when their images overlap and crisscross each other especially under a high particle density. In order to maximize the valid recovery of data points, neural networks are implemented for these two important processes. For the step of particle overlap decomposition, the back propagation neural network is utilized because of its ability in pattern recognition with pertinent particle image feature parameters. For the step of particle tracking, the Hopfield neural network is employed to find appropriate particle tracks based on global optimization. Our investigation indicates that the neural networks are very efficient and useful for stereoscopically tracking particles. As an initial assessment of the diagnostic technology performance, laminar water jets with and without pulsation are measured. The jet tip velocity profiles are in good agreement with analytical predictions. Finally, for testing in material processing applications, a simple directional solidification apparatus is built for experimenting with a metal analog of succinonitrile. Its 3-D velocity field at the liquid phase is then measured to be compared with those from numerical computation. Our theoretical, numerical, and experimental investigations have proven STV to be a viable candidate for reliably measuring 3-D flow velocities. With current activities are focused on further improving the processing efficiency, overall accuracy, and automation, the eventual efforts of broad experimental applications and concurrent numerical modeling validation will be vital to many areas in fluid flow and materials processing.

Lee, David↗

Quantum mechanical streamlines. I - Square potential barrier

Exact numerical calculations are made for scattering of quantum mechanical particles hitting a square two-dimensional potential barrier (an exact analog of the Goos-Haenchen optical experiments). Quantum mechanical streamlines are plotted and found to be smooth and continuous, to have continuous first derivatives even through the classical forbidden region, and to form quantized vortices around each of the nodal points. A comparison is made between the present numerical calculations and the stationary wave approximation, and good agreement is found between both the Goos-Haenchen shifts and the reflection coefficients. The time-independent Schroedinger equation for real wavefunctions is reduced to solving a nonlinear first-order partial differential equation, leading to a generalization of the Prager-Hirschfelder perturbation scheme. Implications of the hydrodynamical formulation of quantum mechanics are discussed, and cases are cited where quantum and classical mechanical motions are identical.

Hirschfelder, J. O.↗

The effects of attention sharing in a dynamic dual task environment

There are numerous examples of cases where the human operator is confronted with several tasks occurring simultaneously and continuously in time. The current study is an investigation into the nature of attention sharing between two continuous tasks with independent input-output modes. Eleven subjects were tested using a zero order compensatory control task with three levels of difficulty (input bandwidth) for each subject. As a secondary task on half of the trials, the subjects were also required to verbally shadow a random auditory input. Results from an extensive time and frequency domain analysis of the data are presented and discussed. The evidence supports a single channel model for continuous dual-task control.

Cliff, R. C.↗

A plate with a crack, stiffened by a partially debonded stringer

The effect of a partially debonded infinite stringer on the stress intensity factor at the crack tip is investigated. It is assumed that the stringer is bonded to an isotropic plate through an adhesive, along a line perpendicular to the crack. Due to the high stress concentration around the crack tip and on the adhesive, partial debonding is assumed to develop between the stringer and the plate. The integral equation, obtained from the continuity of displacements along the line where the stringer and the plate are bonded, is solved numerically to give the shear stress distribution between the plate and the stringer and the stress intensity factors at both crack tips.

Arin, K.↗

Numerical solution of the Navier-Stokes equations for arbitrary two-dimensional multi-element airfoils

Abstracts are presented on a method of numerical solution of the Navier-Stokes equation for the flow about arbitrary airfoils, using a numerically generated curvilinear coordinate system having a coordinate line coincident with the body contour. Results of continuing research are reported and include: application of the Navier-Stokes solution in the vorticity-stream function formulation to a number of single airfoils at Reynolds numbers up to 2000; programming of the Navier-Stokes solution for multiple airfoils in the primitive variable formulation; testing of the potential flow solution of multiple bodies; and development of a generalized coordinate system program.

Thompson, J. F.↗

Micromechanics of slip bands on a free surface

A micromechanics analysis for the formation and propagation of slip bands on the free surface of a polycrystal under monotonic loading is presented. For the growth of slip bands, the analysis satisfies the conditions of both equilibrium and displacement continuity, as well as the relation between slip and the resolved shear stress throughout the polycrystal. Numerical calculations show how the microstress field causes the concentration of plastic deformation in discrete sliding bands and give results which are in good qualitative agreement with known slip band observations on aluminum single crystals.

Lin, S. R.↗

A computer simulation of the afterburning processes occurring within solid rocket motor plumes in the troposphere

As part of a continuing study of the environmental effects of solid rocket motor (SRM) operations in the troposphere, a numerical model was used to simulate the afterburning processes occurring in solid rocket motor plumes and to predict the quantities of potentially harmful chemical species which are created. The calculations include the effects of finite-rate chemistry and turbulent mixing. It is found that the amount of NO produced is much less than the amount of HCl present in the plume, that chlorine will appear predominantly in the form of HCl although some molecular chlorine is present, and that combustion is complete as is evident from the predominance of carbon dioxide over carbon monoxide.

Gomberg, R. I.↗

Finite difference solutions to shocked acoustic waves

The MacCormack, Lambda and split flux finite differencing schemes are used to solve a one dimensional acoustics problem. Two duct configurations were considered, a uniform duct and a converging-diverging nozzle. Asymptotic solutions for these two ducts are compared with the numerical solutions. When the acoustic amplitude and frequency are sufficiently high the acoustic signal shocks. This condition leads to a deterioration of the numerical solutions since viscous terms may be required if the shock is to be resolved. A continuous uniform duct solution is considered to demonstrate how the viscous terms modify the solution. These results are then compared with a shocked solution with and without viscous terms. Generally it is found that the most accurate solutions are those obtained using the minimum possible viscosity coefficients. All of the schemes considered give results accurate enough for acoustic power calculations with no one scheme performing significantly better than the others.

Walkington, N. J.↗

Resonant absorption of phase-mixed Alfven surface waves in ideal and resistive magnetohydrodynamics - Initial-value problem

An initial-value approach to the study of Alfven surface waves, in which the linearized MHD equations are solved numerically in time and space, is considered. A disturbance of a particular wavelength is excited in a nonuniform plasma and its temporal evolution is simulated with numerical solutions of the ideal and resistive MHD equations in slab geometry. This procedure generates a continuous spectrum and therefore simulates phase-mixing and subsequent resonant absorption. Energy in the phase-mixed surface waves resistively dissipates with the absorption time and width scaling as resistivity to the -1/3 and 1/6 powers, respectively.

Steinolfson, R. S.↗

An embedding method for the steady Euler equations

Certain difficulties arise in connection with the numerical solution of a direct finite difference representation of the steady Euler equations. Johnson (1979, 1981, 1982) has, therefore, proposed a surrogate-equation technique, in which the first-order steady Euler equations are embedded in a certain second-order system of equations. The present paper is concerned with the theoretical justification for such an embedding approach. For the numerical solution of the two-dimensional steady Euler equations, it is shown that, under a continuity restriction, it is possible to solve a second-order embedded system together with appropriate additional boundary conditions. The result indicates that a more direct and potentially more efficient approach to the steady solutions exists than the alternative of solving the unsteady equations.

Chang, S.-H.↗

On the eigenvalue and eigenvector derivatives of a non-defective matrix

A novel approach is introduced to address the problem of existence of differentiable eigenvectors for a nondefective matrix which may have repeated eigenvalues. The existence of eigenvector derivatives for a unique set of continuous eigenvectors corresponding to a repeated eigenvalue is rigorously established for nondefective and analytic matrices. A numerically implementable method is then developed to compute the differentiable eigenvectors associated with repeated eigenvalues. The solutions of eigenvalue and eigenvector derivatives for repeated eigenvalues are then derived. An example is given to illustrate the validity of formulations developed in this paper.

Juang, Jer-Nan↗

Computational fluid dynamic control

A general technique is presented for modeling fluid, or gas, dynamic systems specifically for the development of control systems. The numerical methods which are generally used in computational fluid dynamics are borrowed to create either continuous-time or discrete-time models of the particular fluid system. The resulting equations can be either left in a nonlinear form, or easily linearized about an operating point. As there are typically very many states in these systems, the usual linear model reduction methods can be used on them to allow a low-order controller to be designed. A simple example is given which typifies many internal flow control problems. The resulting control is termed computational fluid dynamic control.

Hartley, Tom T.↗

A comparison of artificial compressibility and fractional step methods for incompressible flow computations

We have applied and compared the efficiency and accuracy of two commonly used numerical methods for the solution of Navier-Stokes equations. The artificial compressibility method augments the continuity equation with a transient pressure term and allows one to solve the modified equations as a coupled system. Due to its implicit nature, one can have the luxury of taking a large temporal integration step at the expense of higher memory requirement and larger operation counts per step. Meanwhile, the fractional step method splits the Navier-Stokes equations into a sequence of differential operators and integrates them in multiple steps. The memory requirement and operation count per time step are low, however, the restriction on the size of time marching step is more severe. To explore the strengths and weaknesses of these two methods, we used them for the computation of a two-dimensional driven cavity flow with Reynolds number of 100 and 1000, respectively. Three grid sizes, 41 x 41, 81 x 81, and 161 x 161 were used. The computations were considered after the L2-norm of the change of the dependent variables in two consecutive time steps has fallen below 10(exp -5).

Chan, Daniel C.↗

Transient scattering from two-dimensional dielectric cylinders of arbitrary shape

In this work, the problem of transient scattering by arbitrary shaped two-dimensional dielectric cylinders is solved using the marching-on-in-time (MOT) technique. The dielectric problem is approached via the surface equivalence principle. A pair of coupled integral equations are derived by enforcing the continuity of the electric and magnetic fields which are solved by using the method of moments. Numerical results are presented for two cross sections, viz. a circle and a square, and compared with inverse discrete Fourier transform (IDFT) techniques. In each case, good agreement is obtained with the IDFT solution.

Vechinski, Douglas A.↗