Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “order”

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,171 records · Page 65

A High-Order Direct Solver for Helmholtz Equations with Neumann Boundary Conditions

In this study, a compact finite-difference discretization is first developed for Helmholtz equations on rectangular domains. Special treatments are then introduced for Neumann and Neumann-Dirichlet boundary conditions to achieve accuracy and separability. Finally, a Fast Fourier Transform (FFT) based technique is used to yield a fast direct solver. Analytical and experimental results show this newly proposed solver is comparable to the conventional second-order elliptic solver when accuracy is not a primary concern, and is significantly faster than that of the conventional solver if a highly accurate solution is required. In addition, this newly proposed fourth order Helmholtz solver is parallel in nature. It is readily available for parallel and distributed computers. The compact scheme introduced in this study is likely extendible for sixth-order accurate algorithms and for more general elliptic equations.

Sun, Xian-He↗

Formation Energy of Disordered Alloys form the Energetics of Ordered Compounds

The cluster expansion technique is used in combination with first-principles calculations of the total energy of ordered compounds to study the energetics of the disordered phase for a number of alloy systems. The effect of short-range order, as seen in the energy differences between an alloy with the configuration corresponding to finite temperature and the perfectly random one is studied. The role of relaxation due to large size mismatch is discussed in terms of an effective cluster volume approximation. Very good agreement with measured energies of formation is obtained when short-range order and relaxation effects are taken into account.

Amador, Carlos↗

Higher-Order Compact Schemes for Numerical Simulation of Incompressible Flows

A higher order accurate numerical procedure has been developed for solving incompressible Navier-Stokes equations for 2D or 3D fluid flow problems. It is based on low-storage Runge-Kutta schemes for temporal discretization and fourth and sixth order compact finite-difference schemes for spatial discretization. The particular difficulty of satisfying the divergence-free velocity field required in incompressible fluid flow is resolved by solving a Poisson equation for pressure. It is demonstrated that for consistent global accuracy, it is necessary to employ the same order of accuracy in the discretization of the Poisson equation. Special care is also required to achieve the formal temporal accuracy of the Runge-Kutta schemes. The accuracy of the present procedure is demonstrated by application to several pertinent benchmark problems.

Wilson, Robert V.↗

The Use of Higher-Order Difference Methods in Beam Vibration Analysis

Simple and higher-order difference methods for the solution for the natural frequencies of vibration of a uniform beam are compared. The same basic higher-order method is used throughout for the interior cells, but three different methods of boundary-condition representation are given. Tables and graphs of the error in mode frequencies, as compared with a continuous beam, are given for the various methods as a function of the number of cells. It is concluded that higher-order methods improve accuracy for a given number of cells, with essentially no change in the quantity of computing equipment required.

Greenwood, Donald T.↗

A Non-Dissipative Staggered Fourth-Order Accurate Explicit Finite Difference Scheme for the Time-Domain Maxwell's Equations

We consider a divergence-free non-dissipative fourth-order explicit staggered finite difference scheme for the hyperbolic Maxwell's equations. Special one-sided difference operators are derived in order to implement the scheme near metal boundaries and dielectric interfaces. Numerical results show the scheme is long-time stable, and is fourth-order convergent over complex domains that include dielectric interfaces and perfectly conducting surfaces. We also examine the scheme's behavior near metal surfaces that are not aligned with the grid axes, and compare its accuracy to that obtained by the Yee scheme.

Yefet, Amir↗

Molecular Static Third-Order Polarizabilities of Carbon-Cage Fullerene and Their Correlation with Three Geometric Properties: Symmetry, Aromaticity, and Size

The static third-order polarizabilities (gamma) of C60, C70, five isomers of C78 and two isomers of C84 were analyzed in terms of three properties, from a geometric point of view: symmetry, aromaticity and size. The polarizability values were based on the finite field approximation using a semiempirical Hamiltonian (AM1) and applied to molecular structures obtained from density functional theory calculations. Symmetry was characterized by the molecular group order. The selection of 6-member rings as aromatic was determined from an analysis of bond lengths. Maximum interatomic distance and surface area were the parameters considered with respect to size. Based on triple linear regression analysis, it was found that the static linear polarizability (alpha) and gamma in these molecules respond differently to geometrical properties: alpha depends almost exclusively on surface area while gamma is affected by a combination of number of aromatic rings, length and group order, in decreasing importance. In the case of alpha, valence electron contributions provide the same information as all-electron estimates. For gamma, the best correlation coefficients are obtained when all-electron estimates are used and when the dependent parameter is ln(gamma) instead of gamma.

Moore, C. E.↗

Efficient Simulation of Wing Modal Response: Application of 2nd Order Shape Sensitivities and Neural Networks

At the preliminary design stage of a wing structure, an efficient simulation, one needing little computation but yielding adequately accurate results for various response quantities, is essential in the search of optimal design in a vast design space. In the present paper, methods of using sensitivities up to 2nd order, and direct application of neural networks are explored. The example problem is how to decide the natural frequencies of a wing given the shape variables of the structure. It is shown that when sensitivities cannot be obtained analytically, the finite difference approach is usually more reliable than a semi-analytical approach provided an appropriate step size is used. The use of second order sensitivities is proved of being able to yield much better results than the case where only the first order sensitivities are used. When neural networks are trained to relate the wing natural frequencies to the shape variables, a negligible computation effort is needed to accurately determine the natural frequencies of a new design.

Kapania, Rakesh K.↗

A Two Colorable Fourth Order Compact Difference Scheme and Parallel Iterative Solution of the 3D Convection Diffusion Equation

A new fourth order compact difference scheme for the three dimensional convection diffusion equation with variable coefficients is presented. The novelty of this new difference scheme is that it Only requires 15 grid points and that it can be decoupled with two colors. The entire computational grid can be updated in two parallel subsweeps with the Gauss-Seidel type iterative method. This is compared with the known 19 point fourth order compact differenCe scheme which requires four colors to decouple the computational grid. Numerical results, with multigrid methods implemented on a shared memory parallel computer, are presented to compare the 15 point and the 19 point fourth order compact schemes.

Zhang, Jun↗

Higher-Order Theory for Functionally Graded Materials

This paper presents the full generalization of the Cartesian coordinate-based higher-order theory for functionally graded materials developed by the authors during the past several years. This theory circumvents the problematic use of the standard micromechanical approach, based on the concept of a representative volume element, commonly employed in the analysis of functionally graded composites by explicitly coupling the local (microstructural) and global (macrostructural) responses. The theoretical framework is based on volumetric averaging of the various field quantities, together with imposition of boundary and interfacial conditions in an average sense between the subvolumes used to characterize the composite's functionally graded microstructure. The generalization outlined herein involves extension of the theoretical framework to enable the analysis of materials characterized by spatially variable microstructures in three directions. Specialization of the generalized theoretical framework to previously published versions of the higher-order theory for materials functionally graded in one and two directions is demonstrated. In the applications part of the paper we summarize the major findings obtained with the one-directional and two-directional versions of the higher-order theory. The results illustrate both the fundamental issues related to the influence of microstructure on microscopic and macroscopic quantities governing the response of composites and the technologically important applications. A major issue addressed herein is the applicability of the classical homogenization schemes in the analysis of functionally graded materials. The technologically important applications illustrate the utility of functionally graded microstructures in tailoring the response of structural components in a variety of applications involving uniform and gradient thermomechanical loading.

Aboudi, Jacob↗

Radiation Heat Transfer Between Diffuse-Gray Surfaces Using Higher Order Finite Elements

This paper presents recent work on developing methods for analyzing radiation heat transfer between diffuse-gray surfaces using p-version finite elements. The work was motivated by a thermal analysis of a High Speed Civil Transport (HSCT) wing structure which showed the importance of radiation heat transfer throughout the structure. The analysis also showed that refining the finite element mesh to accurately capture the temperature distribution on the internal structure led to very large meshes with unacceptably long execution times. Traditional methods for calculating surface-to-surface radiation are based on assumptions that are not appropriate for p-version finite elements. Two methods for determining internal radiation heat transfer are developed for one and two-dimensional p-version finite elements. In the first method, higher-order elements are divided into a number of sub-elements. Traditional methods are used to determine radiation heat flux along each sub-element and then mapped back to the parent element. In the second method, the radiation heat transfer equations are numerically integrated over the higher-order element. Comparisons with analytical solutions show that the integration scheme is generally more accurate than the sub-element method. Comparison to results from traditional finite elements shows that significant reduction in the number of elements in the mesh is possible using higher-order (p-version) finite elements.

Gould, Dana C.↗

Low Order Equivalent System Identification for the Tu-144LL Supersonic Transport Aircraft

Low order equivalent system models were identified from flight test data for the Tu-144LL supersonic transport aircraft. Flight test maneuvers were executed by Russian and American test Pilots flying the aircraft from Zhukovsky airfield outside Moscow, Russia. Flight tests included longitudinal and lateral / directional maneuvers at supersonic cruise flight conditions. Piloted frequency sweeps and multi-step maneuvers were used to generate data for closed loop low order equivalent system modeling. Model parameters were estimated using a flexible. high accuracy Fourier transform and an equation error / output error (EE/OE) formulation in the frequency domain. Results were compared to parameter estimates obtained using spectral estimation and subsequent least squares fit to frequency response data in Bode plots. Modeling results from the two methods agreed well for both a frequency sweep and multiple concatenated multi-step maneuvers. For a single multi-step maneuvers the EE/OE method gave a better model fit with improved prediction capability. A summary of closed loop low order equivalent system identification results for the Tu-144LL, including estimated parameters, standard errors, and flying qualities level predictions, were computed and tabulated.

Morelli, Eugene A.↗

Identification of Low Order Equivalent System Models From Flight Test Data

Identification of low order equivalent system dynamic models from flight test data was studied. Inputs were pilot control deflections, and outputs were aircraft responses, so the models characterized the total aircraft response including bare airframe and flight control system. Theoretical investigations were conducted and related to results found in the literature. Low order equivalent system modeling techniques using output error and equation error parameter estimation in the frequency domain were developed and validated on simulation data. It was found that some common difficulties encountered in identifying closed loop low order equivalent system models from flight test data could be overcome using the developed techniques. Implications for data requirements and experiment design were discussed. The developed methods were demonstrated using realistic simulation cases, then applied to closed loop flight test data from the NASA F-18 High Alpha Research Vehicle.

Morelli, Eugene A.↗

Low Order Equivalent System for the Tu-144LL Supersonic Transport Aircraft

Low order equivalent system models were identified from flight test data for the Tu- 144LL supersonic transport aircraft. Flight test maneuvers were executed by Russian and American test pilots flying the aircraft from Zhukovsky airfield outside Moscow, Russia. Flight tests included longitudinal and lateral/directional maneuvers at supersonic cruise flight conditions. Piloted frequency sweeps and multi-step maneuvers were used to Generate data for p closed loop low order equivalent system modeling Model parameters were estimated using a flexible, high accuracy Fourier transform and an equation error output error (EE/OE) formulation in the frequency domain. Results were compared to parameter estimates obtained using spectral estimation and subsequent least squares fit to frequency response data in Bode plots. Modeling results from the two methods a-reed well for both a frequency sweep and multiple concatenated multi-step maneuvers. For a single multi-step maneuver, the EE/OE method gave a better model fit with improved prediction capability. A summary of closed loop low order equivalent system identification results for the Tu-144LL. including estimated parameters, standard errors, and flying qualities level predictions, were computed and tabulated.

Morelli, Eugene A.↗

HOTCFGM-2D: A Coupled Higher-Order Theory for Cylindrical Structural Components with Bi-Directionally Components with Bi-Directionally Graded Microstructures

The objective of this two-year project was to develop and deliver to the NASA-Glenn Research Center a two-dimensional higher-order theory, and related computer codes, for the analysis and design of cylindrical functionally graded materials/structural components for use in advanced aircraft engines (e.g., combustor linings, rotor disks, heat shields, brisk blades). To satisfy this objective, two-dimensional version of the higher-order theory, HOTCFGM-2D, and four computer codes based on this theory, for the analysis and design of structural components functionally graded in the radial and circumferential directions were developed in the cylindrical coordinate system r-Theta-z. This version of the higher-order theory is a significant generalization of the one-dimensional theory, HOTCFGM-1D, developed during the FY97 for the analysis and design of cylindrical structural components with radially graded microstructures. The generalized theory is applicable to thin multi-phased composite shells/cylinders subjected to steady-state thermomechanical, transient thermal and inertial loading applied uniformly along the axial direction such that the overall deformation is characterized by a constant average axial strain. The reinforcement phases are uniformly distributed in the axial direction, and arbitrarily distributed in the radial and circumferential direction, thereby allowing functional grading of the internal reinforcement in the r-Theta plane. The four computer codes fgmc3dq.cylindrical.f, fgmp3dq.cylindrical.f, fgmgvips3dq.cylindrical.f, and fgmc3dq.cylindrical.transient.f are research-oriented codes for investigating the effect of functionally graded architectures, as well as the properties of the multi-phase reinforcement, in thin shells subjected to thermomechanical and inertial loading, on the internal temperature, stress and (inelastic) strain fields. The reinforcement distribution in the radial and circumferential directions is specified by the user. The thermal and inelastic properties of the individual phases can vary with temperature. The inelastic phases are presently modeled by the power-law creep model generalized to multi-directional loading (within fgmc3dq.cylindrical.f and fgmc3dq.cylindrical.transient.f for steady-state and transient thermal loading, respectively), and incremental plasticity and GVIPS unified viscoplasticity theories (within the steady-state loading versions fgmp3dq.cylindrical.f and fgmgvips3dq.cylindrical.f).

Pindera, Marek-Jerzy↗

Reduced-Order Modeling: New Approaches for Computational Physics

In this paper, we review the development of new reduced-order modeling techniques and discuss their applicability to various problems in computational physics. Emphasis is given to methods ba'sed on Volterra series representations and the proper orthogonal decomposition. Results are reported for different nonlinear systems to provide clear examples of the construction and use of reduced-order models, particularly in the multi-disciplinary field of computational aeroelasticity. Unsteady aerodynamic and aeroelastic behaviors of two- dimensional and three-dimensional geometries are described. Large increases in computational efficiency are obtained through the use of reduced-order models, thereby justifying the initial computational expense of constructing these models and inotivatim,- their use for multi-disciplinary design analysis.

Beran, Philip S.↗

Validation of a High-Order Prefactored Compact Scheme on Nonlinear Flows with Complex Geometries

Three benchmark problems are solved using a sixth-order prefactored compact scheme employing an explicit 10th-order filter with optimized fourth-order Runge-Kutta time stepping. The problems solved are the following: (1) propagation of sound waves through a transonic nozzle; (2) shock-sound interaction; and (3) single airfoil gust response. In the first two problems, the spatial accuracy of the scheme is tested on a stretched grid, and the effectiveness of boundary conditions is shown. The solution stability and accuracy near a shock discontinuity is shown as well. Also, 1-D nonlinear characteristic boundary conditions will be evaluated. In the third problem, a nonlinear Euler solver will be used that solves the equations in generalized curvilinear coordinates using the chain rule transformation. This work, continuing earlier work on flat-plate cascades and Joukowski airfoils, will focus mainly on the effect of the grid and boundary conditions on the accuracy of the solution. The grids were generated using a commercially available grid generator, GridPro/az3000.

Hixon, Ray↗

New Third-Order Moments for the CBL

Turbulent convection is an inherently non-local phenomenon and a primary condition for a successful treatment of the CBL (convective boundary layer) is a reliable model of non-locality. In the dynamic equations governing the convective flux, the turbulent kinetic energy, etc., non-locality is represented by the third-order moments, TOMs. Since the simplest form, the so-called down gradient approximation (DGA), severely underestimates the TOMs (up to an order of magnitude), a more physical model is needed. In 1994, an analytical model was presented which was derived directly from the dynamical equations for the TOMs. It considerably improved the DGA but was a bit cumbersome to use and, more importantly, it was based on the quasi-normal (QN) approximation for the fourth-order moments. Here, we present a new analytic expression for the TOMs which is structurally simpler than the 1994 expression and which avoids the QN approximation. The resulting fit to the LES data is superior to that of the 1994 model.

Canuto, V. M.↗

Technique for Very High Order Nonlinear Simulation and Validation

Finding the sources of sound in large nonlinear fields via direct simulation currently requires excessive computational cost. This paper describes a simple technique for efficiently solving the multidimensional nonlinear Euler equations that significantly reduces this cost and demonstrates a useful approach for validating high order nonlinear methods. Up to 15th order accuracy in space and time methods were compared and it is shown that an algorithm with a fixed design accuracy approaches its maximal utility and then its usefulness exponentially decays unless higher accuracy is used. It is concluded that at least a 7th order method is required to efficiently propagate a harmonic wave using the nonlinear Euler equations to a distance of 5 wavelengths while maintaining an overall error tolerance that is low enough to capture both the mean flow and the acoustics.

Dyson, Rodger W.↗