Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Polynomial models”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Analysis of aircraft tires via semianalytic finite elements

A computational procedure is presented for the geometrically nonlinear analysis of aircraft tires. The tire was modeled by using a two-dimensional laminated anisotropic shell theory with the effects of variation in material and geometric parameters included. The four key elements of the procedure are: (1) semianalytic finite elements in which the shell variables are represented by Fourier series in the circumferential direction and piecewise polynomials in the meridional direction; (2) a mixed formulation with the fundamental unknowns consisting of strain parameters, stress-resultant parameters, and generalized displacements; (3) multilevel operator splitting to effect successive simplifications, and to uncouple the equations associated with different Fourier harmonics; and (4) multilevel iterative procedures and reduction techniques to generate the response of the shell.

Noor, Ahmed K.↗
NASA NTRSDate not supplied

Jig Twist Optimization of Mach 0.8 Transonic Truss-Braced Wing Aircraft

This paper presents a jig twist optimization study of Mach 0.8 Transonic Truss-Braced Wing (TTBW) aircraft using FUN3D. The solver has been validated against wind tunnel test data. An aeroelastic model deformation correction method is developed for the wind tunnel model to improve the validation. The flight condition corresponding to Mach 0.8 and design lift coefficient 0.695 is selected for the TTBW aircraft jig twist optimization to reduce the drag coefficient. The static aeroelastic simulation module in FUN3D is used for the aeroelastic simulation. Three different orders of the Chebyshev polynomials are used as shape function to represent the variations of the jig twist distribution along the wing span. The preliminary optimization results show that the Mach 0.8 TTBW aircraft with the optimized jig twist achieves a drag reduction of about 5.7 counts at the design condition when 5 𝑡ℎ order Chebyshev polynomials are used as shape function.

TTBW↗
NASA NTRS1982

Line shape of the non-thermal 6300 A O/1D/ emission

The two-population model of Schmitt, Abreu and Hays (1981) is used to calculate the line shape of the atomic oxygen metastable state, nonthermal O(1D) 6300 A emission, in order to simulate observations made from a space platform at different zenith angles and altitudes. The Addition theorem, for spherical harmonics of a Legendre polynomial expansion of the nonthermal population distribution function, is used to obtain nonthermal line shapes observed at zenith angles other than the local vertical one.

Schmitt, G. A.↗
NASA NTRS1995

Estimating errors in least-squares fitting

While least-squares fitting procedures are commonly used in data analysis and are extensively discussed in the literature devoted to this subject, the proper assessment of errors resulting from such fits has received relatively little attention. The present work considers statistical errors in the fitted parameters, as well as in the values of the fitted function itself, resulting from random errors in the data. Expressions are derived for the standard error of the fit, as a function of the independent variable, for the general nonlinear and linear fitting problems. Additionally, closed-form expressions are derived for some examples commonly encountered in the scientific and engineering fields, namely ordinary polynomial and Gaussian fitting functions. These results have direct application to the assessment of the antenna gain and system temperature characteristics, in addition to a broad range of problems in data analysis. The effects of the nature of the data and the choice of fitting function on the ability to accurately model the system under study are discussed, and some general rules are deduced to assist workers intent on maximizing the amount of information obtained form a given set of measurements.

Richter, P. H.↗
NASA NTRS1982

The optimization of convergence for Chebyshev polynomial methods in an unbounded domain

Grosch and Orszag (1977) have performed a numerical analysis of the problem of solving differential equations in a semiinfinite or infinite domain using Chebyshev polynomials. The principal limitation of the conducted study was that it was entirely empirical. Various differential equations were solved in different ways and the numbers were compared. The present investigation has the objective to extend the studies conducted by Grosch and Orszag by deriving asymptotic approximations to the Chebyshev coefficients of simple model functions. This approach makes it possible to conduct more systematic comparisons of different methods, extend the range of comparisons, and, perhaps most important, give simple analytic formulas for choosing the optimum domain size or mapping parameter L for various situations.

Boyd, J. P.↗
NASA NTRS1993

Random interactions in higher order neural networks

Recurrent networks of polynomial threshold elements with random symmetric interactions are studied. Precise asymptotic estimates are derived for the expected number of fixed points as a function of the margin of stability. In particular, it is shown that there is a critical range of margins of stability (depending on the degree of polynomial interaction) such that the expected number of fixed points with margins below the critical range grows exponentially with the number of nodes in the network, while the expected number of fixed points with margins above the critical range decreases exponentially with the number of nodes in the network. The random energy model is also briefly examined and links with higher order neural networks and higher order spin glass models made explicit.

Baldi, Pierre↗
NASA NTRS1996

Element Library for Three-Dimensional Stress Analysis by the Integrated Force Method

The Integrated Force Method, a recently developed method for analyzing structures, is extended in this paper to three-dimensional structural analysis. First, a general formulation is developed to generate the stress interpolation matrix in terms of complete polynomials of the required order. The formulation is based on definitions of the stress tensor components in term of stress functions. The stress functions are written as complete polynomials and substituted into expressions for stress components. Then elimination of the dependent coefficients leaves the stress components expressed as complete polynomials whose coefficients are defined as generalized independent forces. Such derived components of the stress tensor identically satisfy homogenous Navier equations of equilibrium. The resulting element matrices are invariant with respect to coordinate transformation and are free of spurious zero-energy modes. The formulation provides a rational way to calculate the exact number of independent forces necessary to arrive at an approximation of the required order for complete polynomials. The influence of reducing the number of independent forces on the accuracy of the response is also analyzed. The stress fields derived are used to develop a comprehensive finite element library for three-dimensional structural analysis by the Integrated Force Method. Both tetrahedral- and hexahedral-shaped elements capable of modeling arbitrary geometric configurations are developed. A number of examples with known analytical solutions are solved by using the developments presented herein. The results are in good agreement with the analytical solutions. The responses obtained with the Integrated Force Method are also compared with those generated by the standard displacement method. In most cases, the performance of the Integrated Force Method is better overall.

Kaljevic, Igor↗
NASA NTRS1989

Mathematical models for film sensitivity measurements

The quality of the pictorial record developed from photographic material depends on the composition of the material, the procedures used in its exposure and processing, as well as, the nature of the physical and temporal environment extant during the creation of the record. By holding many of the variables fixed, the effect of the two environmental parameters, namely temperature and aging are examined on the characteristic curve of a given film. Polynomial and exponential functions are evaluated as empirical formulas for the characteristic curve, and the sensitivity of derived coefficients to the selected parameters is assessed.

Hammond, Ernest C., Jr.↗
NASA NTRS1990

Modeling and synthesis of multicomputer interconnection networks

The type of interconnection network employed has a profound effect on the performance of a multicomputer and multiprocessor design. Adequate models are needed to aid in the design and development of interconnection networks. A novel modeling approach using statistical and optimization techniques is described. This method represents an attempt to compare diverse interconnection network designs in a way that allows not only the best of existing designs to be identified but to suggest other, perhaps hybrid, networks that may offer better performance. Stepwise linear regression is used to develop a polynomial surface representation of performance in a (k+1) space with a total of k quantitative and qualitative independent variables describing graph-theoretic characteristics such as size, average degree, diameter, radius, girth, node-connectivity, edge-connectivity, minimum dominating set size, and maximum number of prime node and edge cutsets. Dependent variables used to measure performance are average message delay and the ratio of message completion rate to network connection cost. Response Surface Methodology (RSM) optimizes a response variable from a polynomial function of several independent variables. Steepest ascent path may also be used to approach optimum points.

Standley, Hilda M.↗
NASA NTRS2023

Chemical Thermodynamics and the Mathematical Integration of Reaction Kinetics

Key advances in the development of numerical methods for non-reacting compressible flows have been enabled by translating physical requirements into concrete numerical guidelines, such as the satisfaction of entropy inequalities for shock-capturing techniques [Lax, Contributions to Nonlinear Functional Analysis (1971) 603-634]. In the present work, we present nonlinear numerical analysis tools that draw from Chemical Thermodynamics , the branch of Nonequilibrium Thermodynamics that deals with chemical reactions. Through Gibbs formalism, chemical thermodynamics provides a well-known theoretical expression for the chemical equilibrium constant of a reaction in terms of reduced chemical potentials. A less-known, yet extremely valuable result, due to [Krambeck, Arch. Ration. Mech. Anal. , 38 (1970) 317], states that when this expression is implemented, mass-action kinetic models are consistent with the dynamical prescriptions of the 2nd law of thermodynamics. For fixed-temperature ordinary differential equations modeling constant-volume reacting gas mixtures, this leads to a decreasing Helmholtz free energy. If the temperature is allowed to vary in accordance with conservation of energy (1st law), this leads to the statement of increasing entropy. These nonlinear prescriptions can, and should be, used to further develop temporal integration techniques for reaction kinetics. We demonstrate that Krambeck's result holds even when the equilibrium constants are approximated from data. We prove this result by constructing the implicit free energy and the implicit entropy inherent to a given approximation. This is first done for a 5-species, 17-reaction model problem for air. With this structure established, elements of discrete entropy-stability theory [Tadmor, Acta Numer. , 12 (2003) 451] are leveraged to examine the consistency of time-integration schemes with these prescriptions. Using chemical potentials, one can compute the respective contributions of the kinetics model and of the temporal scheme to free energy/entropy variations. We introduce a nonlinear-stable version of the Discontinuous-Galerkin (DG) scheme in time which shows robustness improvements over the standard linearly-stable version. Most notably, the maximum timestep that can be resolved with the nonlinearly-stable variant tends to grow with polynomial order, in contrast to the linearly-stable variant. We generalize our constructions to arbitrary systems of reversible chemical reactions, ultimately showing that the compressible reacting Euler system admits the opposite of the implicitly constructed thermodynamic entropy as a mathematical entropy . This lays important theoretical foundations towards robust scheme development [Harten, J. Comput. Phys. 49 (1983) 151-164].

STMD↗
NASA NTRS2005

Fitting Nonlinear Curves by use of Optimization Techniques

MULTIVAR is a FORTRAN 77 computer program that fits one of the members of a set of six multivariable mathematical models (five of which are nonlinear) to a multivariable set of data. The inputs to MULTIVAR include the data for the independent and dependent variables plus the user s choice of one of the models, one of the three optimization engines, and convergence criteria. By use of the chosen optimization engine, MULTIVAR finds values for the parameters of the chosen model so as to minimize the sum of squares of the residuals. One of the optimization engines implements a routine, developed in 1982, that utilizes the Broydon-Fletcher-Goldfarb-Shanno (BFGS) variable-metric method for unconstrained minimization in conjunction with a one-dimensional search technique that finds the minimum of an unconstrained function by polynomial interpolation and extrapolation without first finding bounds on the solution. The second optimization engine is a faster and more robust commercially available code, denoted Design Optimization Tool, that also uses the BFGS method. The third optimization engine is a robust and relatively fast routine that implements the Levenberg-Marquardt algorithm.

Hill, Scott A.↗
NASA NTRSDate not supplied

Preparing Quantum Many-Body Scar States on Quantum Computers

Quantum many-body scar states are highly excited eigenstates of many-body systems that exhibit atypical entanglement and correlation properties relative to other eigenstates at the same energy density. Scar states also give rise to infinitely long-lived coherent dynamics when the system is prepared in a special initial state having finite overlap with them. Many models with exact scar states have been constructed, but the fate of scarred eigenstates and dynamics when these models are perturbed is difficult to study with classical computational techniques. In this work, we propose state preparation protocols for individual scar states in a particular model, as well as superpositions of them that give rise to coherent dynamics. For superpositions of scar states, we propose both a linear depth unitary and a finite-depth nonunitary state preparation protocol, the latter of which uses measurement and postselection to reduce the circuit depth. For individual scarred eigenstates, we propose a circuit architecture with polynomial depth. We also provide proof of principle implementations of these protocols on superconducting quantum hardware.

Quantum Computing↗
NASA NTRS1974

Theoretical predictions of latitude dependencies in the solar wind

Results are presented which were obtained with the Winge-Coleman model for theoretical predictions of latitudinal dependencies in the solar wind. A first-order expansion is described which allows analysis of first-order latitudinal variations in the coronal boundary conditions and results in a second-order partial differential equation for the perturbation stream function. Latitudinal dependencies are analytically separated out in the form of Legendre polynomials and their derivative, and are reduced to the solution of radial differential equations. This analysis is shown to supply an estimate of how large the coronal variation in latitude must be to produce an 11 km/sec/deg gradient in the radial velocity of the solar wind, assuming steady-state processes.

Winge, C. R., Jr.↗
NASA NTRS1996

Development of Finite Elements for Two-Dimensional Structural Analysis Using the Integrated Force Method

The Integrated Force Method has been developed in recent years for the analysis of structural mechanics problems. This method treats all independent internal forces as unknown variables that can be calculated by simultaneously imposing equations of equilibrium and compatibility conditions. In this paper a finite element library for analyzing two-dimensional problems by the Integrated Force Method is presented. Triangular- and quadrilateral-shaped elements capable of modeling arbitrary domain configurations are presented. The element equilibrium and flexibility matrices are derived by discretizing the expressions for potential and complementary energies, respectively. The displacement and stress fields within the finite elements are independently approximated. The displacement field is interpolated as it is in the standard displacement method, and the stress field is approximated by using complete polynomials of the correct order. A procedure that uses the definitions of stress components in terms of an Airy stress function is developed to derive the stress interpolation polynomials. Such derived stress fields identically satisfy the equations of equilibrium. Moreover, the resulting element matrices are insensitive to the orientation of local coordinate systems. A method is devised to calculate the number of rigid body modes, and the present elements are shown to be free of spurious zero-energy modes. A number of example problems are solved by using the present library, and the results are compared with corresponding analytical solutions and with results from the standard displacement finite element method. The Integrated Force Method not only gives results that agree well with analytical and displacement method results but also outperforms the displacement method in stress calculations.

Kaljevic, Igor↗
NASA NTRS2009

Formal methods for modeling and analysis of hybrid systems

A technique based on the use of a quantifier elimination decision procedure for real closed fields and simple theorem proving to construct a series of successively finer qualitative abstractions of hybrid automata is taught. The resulting abstractions are always discrete transition systems which can then be used by any traditional analysis tool. The constructed abstractions are conservative and can be used to establish safety properties of the original system. The technique works on linear and non-linear polynomial hybrid systems: the guards on discrete transitions and the continuous flows in all modes can be specified using arbitrary polynomial expressions over the continuous variables. An exemplar tool in the SAL environment built over the theorem prover PVS is detailed. The technique scales well to large and complex hybrid systems.

Tiwari, Ashish↗
NASA NTRS1986

Plates and shells containing a surface crack under general loading conditions

The severity of the underlying assumptions of the line-spring model (LSM) are such that verification with three-dimensional solutions is necessary. Such comparisons show that the model is quite accurate, and therefore, its use in extensive parameter studies is justified. Investigations into the endpoint behavior of the line-spring model have led to important conclusions about the ability of the model to predict stresses in front of the crack tip. An important application of the LSM was to solve the contact plate bending problem. Here the flexibility of the model to allow for any crack shape is exploited. The use of displacement quantities as unknowns in the formulation of the problem leads to strongly singular integral equations, rather than singular integral equations which result from using displacement derivatives. The collocation method of solving the integral equations was found to be better and more convenient than the quadrature technique. Orthogonal polynomials should be used as fitting functions when using the LSM as opposed to simpler functions such as power series.

Joseph, Paul F.↗
NASA NTRS1976

Polyatomic molecule vibrations

Polyatomic molecule vibrations are analyzed as harmonic vibrations along normal coordinates. The energy eigenvalues are found for linear and nonlinear symmetric triatomic molecules for valence bond models of the potential function with arbitrary coupling coefficients; such models can usually be fitted to observed energy levels with reasonably good accuracy. Approximate normal coordinates for the H2O molecule are discussed. Degenerate vibrational modes such as occur in CO2 are analyzed and expressions for Fermi resonance between close-lying states of the same symmetry are developed. The bending modes of linear triatomic molecules are expressed in terms of Laguerre polynomials in cylindrical coordinates as well as in terms of Hermite polynomials in Cartesian coordinates. The effects of large-amplitude bending such as occur in the C3 molecule are analyzed, along with anharmonic effects, which split the usually degenerate bending mode energy levels. Finally, the vibrational frequencies, degeneracies, and symmetry properties of XY3, X2Y2, and XY4 type molecules are discussed.

Source record↗
NASA NTRS1988

Krylov subspace methods on supercomputers

A short survey of recent research on Krylov subspace methods with emphasis on implementation on vector and parallel computers is presented. Conjugate gradient methods have proven very useful on traditional scalar computers, and their popularity is likely to increase as three-dimensional models gain importance. A conservative approach to derive effective iterative techniques for supercomputers has been to find efficient parallel/vector implementations of the standard algorithms. The main source of difficulty in the incomplete factorization preconditionings is in the solution of the triangular systems at each step. A few approaches consisting of implementing efficient forward and backward triangular solutions are described in detail. Polynomial preconditioning as an alternative to standard incomplete factorization techniques is also discussed. Another efficient approach is to reorder the equations so as to improve the structure of the matrix to achieve better parallelism or vectorization. An overview of these and other ideas and their effectiveness or potential for different types of architectures is given.

Saad, Youcef↗