SEARCH · Engineering Papers
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.
Second order equations in terms of auxiliary variables similar to potential and stream functions are obtained by applying a weighted least squares formulation to a first order system. The additional boundary conditions which are necessary to solve the higher order equations are determined and numerical results are presented for the Cauchy-Riemann equations.
Three second-order accurate numerical methods - MacCormack's method, Lambda scheme and Gabutti scheme - are introduced to solve the quasi-linear, hyperbolic partial differential equations describing transient flows in closed conduits. The details of these methods and the treatment of boundary conditions are presented and the results computed by using these methods for a typical piping system are compared. It is shown that for the same accuracy, second-order methods require considerably lesser number of computational nodes and computer time as compared to those required by the first-order methods.
Results from the application of a second generation low-order panel method to a number of powerplant installations are presented. The method, which retains the advantages of ease of use and low operating cost of earlier panel methods while retaining the aerodynamic modelling rigor of higher-order methods was used to predict the flow in and around typical nacelles and aircraft for a wide range of conditions. The ability of the second generation low-order analysis to calculate internal flows without the leakage problems associated with earlier programs is demonstrated and correlation between calculated and measured surface pressures and flow behavior, including predicting streamlines and locating regions of separated flow, is presented.
A quartic spline, S(4,2), is proposed which overcomes some of the difficulties associated with the use of splines S(5,3) and S(3,1) and provides fourth-order accurate results with relatively few grid points. The accuracy of spline S(4,2) is comparable to or better than that of the fourth-order box scheme and the compact differencing scheme. The use of spline S(4,2) is suggested as a possible way of obtaining fourth-order accurate solutions to Navier-Stokes equations.
The development of a 2 deg latitude by 2.5 deg longitude version of the GLAS Fourth Order General Circulation Model has necessitated the specification of a correspondingly fine horizontal resolution orographic data set to act as a lower boundary condition for the GCM. The S1 topography attempts to relate the mean areally averaged height of such terrain while the Q3 topography emphasizes the effect of the flow of the tallest peaks in such regions by specifying an exaggerated significant height. In order to assess the relative merits of these two orographic data sets. A series of tests was initiated to evaluate the comparative weather forecasting and climate simulation skills of the GLAS Fourth Order GCM when the S1 and Q3 topographies are applied to the lower boundary. Three sets of Northern Hemispheric winter forecasts at the 2 deg x 2.5 deg resolution indicate that the Q3 topography gives slightly better results over the Northern Hemisphere than does the S1 topography. It is apparent that more forecasts will be necessary at both horizontal resolutions before any definitive conclusions are made regarding the choice of an optimal orographic data set.
The accuracy and stability of the second order upwind differencing scheme was investigated. The solution algorithm employed is based on a coupled solution of the nonlinear finite difference equations by the multigrid technique. Calculations have been made of the driven cavity flow for several Reynolds numbers and finite difference grids. In comparison with the hybrid differencing, the second order upwind differencing is somewhat more accurate but it is not monotonically accurate with mesh refinement. Also, the convergence of the solution algorithm deteriorates with the use of the second order upwind differencing.
Sufficient excitation conditions are derived for three adaptive observation schemes for the case in which the order of the observer is less than that of the plant being observed. The reduced-order excitation requirement, found to be identical for all three observers, is that the input contain more than (n+N-1)/2 frequencies, where n is the observer order, and N is the dimension of the portion of the plant which is both controllable and observable.
A higher-order shear deformation theory is used to demonstrate the natural frequencies and buckling loads of elastic plates. The theory accounts for parabolic distribution of the transverse shear strains through the thickness of the plate and rotary inertia. Exact solutions of simply supported plates are obtained and the results are compared with the exact solutions of three-dimensional elasticity theory, the first-order shear deformation theory, and the classical plate theory. The present theory predicts the frequencies and buckling loads more accurately when compared to the first-order and classical plate theories.
In the strong fluctuation theory for a bounded layer of random discrete scatterers, the second moments of the fields in the second-order distorted Born approximation are obtained for copolarized and cross-polarized fields. The backscattering cross sections per unit area are calculated by including the mutual coherence of the fields due to the coincidental ray paths, and that due to the opposite ray paths, corresponding to the ladder and cross terms in the Feynman diagramatic representation. It is proved that the contributions from ladder and cross terms for the copolarized backscattering cross sections are the same, while the contributions for the cross-polarized backscattering cross sections are of the same order. The bistatic scattering coefficients in the second-order approximation for both the ladder and cross terms are also obtained. The contributions from the cross terms explain the enhancement in the backscattering direction.
The behavior of higher order elastic properties, which are much more sensitive to material state than are second order properties, has been studied for steel alloys AISI 1016, 1045, 1095, and 8620 by measuring the stress derivative of the acoustic natural velocity to determine the stress acoustic constants (SAC's). Results of these tests show a 20 percent linear variation of SAC's with carbon content as well as even larger variations with prestrain (plastic deformation). The use of higher order elastic characterization permits quantitative evaluation of solids and may prove useful in studies of fatigue and fracture.
Second order equations in terms of auxiliary variables similar to potential and stream functions are obtained by applying a weighted least squares formulation to a first order system. The additional boundary conditions which are necessary to solve the higher order equations are determined and numerical results are presented for the Cauchy-Riemann equations.
The classical approach to the investigation of small-scale intermittency in turbulence is based on the higher-order derivative correlations such as skewness and flatness factors. In the study of the small scales, numerial simulations can provide more detail than experiments. In the present paper, a variety of velocity- and scalar-derivative correlations are calculated over a range of Reynolds numbers. Particular attention is given to third- and fourth-order correlations, taking into account also some fifth- and sixth-order correlations to allow comparisons with the phenomenological models. The governing equations are the incompresssible Navier-Stokes equation for the velocity and transport equation for a passive scalar. Two numerical codes are used for the simulations presented. Attention is given to details regarding the numerical method used, forcing, simulation parameters, spectra and skewnesses, and graphics.
The physical properties of the commonly used second-order closure models are examined theoretically for rotating turbulent flows. Comparisons are made with results which are a rigorous consequence of the Navier-Stokes equations for the problem of a fully-developed turbulent channel flow in a rapidly rotating framework. It is demonstrated that all existing second-order closures yield spurious physical results for this test problem of rotating channel flow. In fact, the results obtained are shown to be substantially more unphysical than those obtained from the simpler K-epsilon and K-l models. Modifications in the basic structure of these second-order closure models are proposed which can alleviate this problem.
The standard point vortex method has recently been shown to be of high order of accuracy for problems on the whole plane, when using a uniform initial subdivision for assigning the vorticity to the points. If obstacles are present in the flow, this high order deteriorates to first or second order. New vortex methods are introduced which are of arbitrary accuracy (under regularity assumptions) regardless of the presence of bodies and the uniformity of the initial subdivision.
Efficient solution of the equations from spectral discretizations is essential if the high-order accuracy of these methods is to be realized. Direct solution of these equations is rarely feasible, thus iterative techniques are required. A preconditioning scheme for first-order Chebyshev collocation operators is proposed herein, in which the central finite difference mesh is finer than the collocation mesh. Details of the proper techniques for transferring information between the meshes are given here, and the scheme is analyzed by examination of the eigenvalue spectra of the preconditioned operators. The effect of artificial viscosity required in the inversion of the finite difference operator is examined. A second preconditioning scheme, involving a high-order upwind finite difference operator of the van Leer type is also analyzed to provide a comparison with the present scheme. Finally, the performance of the present scheme is verified by application to several test problems.
Asymptotic solutions of order k to the nth are developed for the reduced wave equation. Here k is a dimensionless wave number and n is the arbitrary order of the approximation. These approximations are an extension of geometric acoustics theory, and provide corrections to that theory in the form of multiplicative functions which satisfy parabolic partial differential equations. These corrections account for the diffraction effects caused by variation of the field normal to the ray path and the interaction of these transverse variations with the variation of the field along the ray. The theory is applied to the example of radiation from a piston, and it is demonstrated that the higher order approximations are more accurate for decreasing values of k.
The work includes an investigation of the applicability of the nucleation theory to second and higher order thermodynamic transitions in the Ehrenfest sense, and a number of significant conclusions relevant to first order transitions, as well. The underlying theoretical method consisted of expanding the Gibbs' free energy in a Maclarin or Taylor series and then using fundamental thermodynamic determinable quantities, and interpreting the results. Work was performed on the existence and interpretation of an interfacial energy between phases in a second order transition in addition to an investigation of the solid-liquid interfacial energy for various polymers. Extensive considerations were devoted to various aspects of a particular polymer, polyvinylidene fluoride (PVDF or PVF2), including an experimetal investigation of the effects of an applied electric field on the morphology of melt crystallization and on the nucleation and growth of polarized domains.
A third-order shear deformation theory of laminated composite plates and shells is developed, the Navier solutions are derived, and its finite element models are developed. The theory allows parabolic description of the transverse shear stresses, and therefore the shear correction factors of the usual shear deformation theory are not required in the present theory. The theory also accounts for the von Karman nonlinear strains. Closed-form solutions of the theory for rectangular cross-ply and angle-ply plates and cross-ply shells are developed. The finite element model is based on independent approximations of the displacements and bending moments (i.e., mixed finite element model), and therefore, only C sup o -approximation is required. The finite element model is used to analyze cross-ply and angle-ply laminated plates and shells for bending and natural vibration. Many of the numerical results presented here should serve as references for future investigations. Three major conclusions resulted from the research: First, for thick laminates, shear deformation theories predict deflections, stresses and vibration frequencies significantly different from those predicted by classical theories. Second, even for thin laminates, shear deformation effects are significant in dynamic and geometrically nonlinear analyses. Third, the present third-order theory is more accurate compared to the classical and firt-order theories in predicting static and dynamic response of laminated plates and shells made of high-modulus composite materials.