Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear problems”

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,081 records · Page 60

Fixed-trim re-entry guidance analysis

The terminal guidance problem for a fixed-trim re-entry body is formulated with the objective of synthesizing a closed-loop steering law. A transformation of variables and subsequent linearization of the motion, with the sight-line to the target as a reference, reduces the order of the state system for the guidance problem. The reduced order system, although nonlinear and time-varying, is simple enough to lend itself to synthesis of a class of guidance laws. A generalization of the feedforward device of classical control theory is successfully employed for compensation of roll autopilot lags. The proposed steering law exhibits superior miss-distance performance in a computational comparison with existing fixed-trim guidance laws.

Gracey, C.↗

Improved finite element methodology for integrated thermal structural analysis

An integrated thermal-structural finite element approach for efficient coupling of thermal and structural analysis is presented. New thermal finite elements which yield exact nodal and element temperatures for one dimensional linear steady state heat transfer problems are developed. A nodeless variable formulation is used to establish improved thermal finite elements for one dimensional nonlinear transient and two dimensional linear transient heat transfer problems. The thermal finite elements provide detailed temperature distributions without using additional element nodes and permit a common discretization with lower order congruent structural finite elements. The accuracy of the integrated approach is evaluated by comparisons with analytical solutions and conventional finite element thermal structural analyses for a number of academic and more realistic problems. Results indicate that the approach provides a significant improvement in the accuracy and efficiency of thermal stress analysis for structures with complex temperature distributions.

Dechaumphai, P.↗

Nonlinear Green's function method for unsteady transonic flows

Advantages to employing Green's function in describing unsteady three-dimensional transonic flows are explored. The development of the function for application to linear subsonic and supersonic unsteady aerodynamics is reviewed. It is shown that unique solutions are possible for external flows, with all functional expressions being defined in Prandtl-Glauert space. The development of methods of using the Green's function for transonic flows is traced, noting the necessity of including the effects of significant nonlinear terms. The steady-state problem is considered to demonstrate the shock-capturing ability of the method and the usefulness of the function in the incompressible, subsonic, transonic, and supersonic areas of potential unsteady three-dimensional flows around complex configurations. Computational time is asserted to be an order of magnitude less than with finite difference methods.

Tseng, K.↗

Improved finite element methodology for integrated thermal structural analysis

An integrated thermal-structural finite element approach for efficient coupling of thermal and structural analyses is presented. New thermal finite elements which yield exact nodal and element temperature for one dimensional linear steady state heat transfer problems are developed. A nodeless variable formulation is used to establish improved thermal finite elements for one dimensional nonlinear transient and two dimensional linear transient heat transfer problems. The thermal finite elements provide detailed temperature distributions without using additional element nodes and permit a common discretization with lower order congruent structural finite elements. The accuracy of the integrated approach is evaluated by comparisons with analytical solutions and conventional finite element thermal-structural analyses for a number of academic and more realistic problems. Results indicate that the approach provides a significant improvement in the accuracy and efficiency of thermal stress analysis for structures with complex temperature distributions.

Dechaumphai, P.↗

Numerical computations on one-dimensional inverse scattering problems

An approximate method to determine the index of refraction of a dielectric obstacle is presented. For simplicity one dimensional models of electromagnetic scattering are treated. The governing equations yield a second order boundary value problem, in which the index of refraction appears as a functional parameter. The availability of reflection coefficients yield two additional boundary conditions. The index of refraction by a k-th order spline which can be written as a linear combination of B-splines is approximated. For N distinct reflection coefficients, the resulting N boundary value problems yield a system of N nonlinear equations in N unknowns which are the coefficients of the B-splines.

Dunn, M. H.↗

Finite-element methods for steady solidification problems

Four Galerkin finite-element methods are tested for solving the free-boundary problem that describes steady solidification. The formulations differ in the solution method used to account for the unknown shape of the melt/solid interface, in the interphase condition (either balance of heat flux or equilibrium of temperature) distinguished for locating the interface, and in the technique used for solving the systems of algebraic equations that result from the finite-element approximations. Methods that use the melting point isotherm to locate the melt/solid interface are found more accurate and efficient than formulations based on the interfacial energy balance. Solution by a Galerkin-Newton algorithm of the free-boundary problem transformed to a fixed domain is most efficient when the field problem in each phase is made nonlinear by including radiation from the melt and solid to the surroundings.

Ettouney, H. M.↗

Numerical computations on one-dimensional inverse scattering problems

An approximate method to determine the index of refraction of a dielectric obstacle is presented. For simplicity one dimensional models of electromagnetic scattering are treated. The governing equations yield a second order boundary value problem, in which the index of refraction appears as a functional parameter. The availability of reflection coefficients yield two additional boundary conditions. The index of refraction by a k-th order spline which can be written as a linear combination of B-splines is approximated. For N distinct reflection coefficients, the resulting N boundary value problems yield a system of N nonlinear equations in N unknowns which are the coefficients of the B-splines.

Dunn, M. H.↗

Interactive application of quadratic expansion of chi-square statistic to nonlinear curve fitting

This report contains a detailed theoretical description of an all-purpose, interactive curve-fitting routine that is based on P. R. Bevington's description of the quadratic expansion of the Chi-Square statistic. The method is implemented in the associated interactive, graphics-based computer program. Taylor's expansion of Chi-Square is first introduced, and justifications for retaining only the first term are presented. From the expansion, a set of n simultaneous linear equations is derived, then solved by matrix algebra. A brief description of the code is presented along with a limited number of changes that are required to customize the program of a particular task. To evaluate the performance of the method and the goodness of nonlinear curve fitting, two typical engineering problems are examined and the graphical and tabular output of each is discussed. A complete listing of the entire package is included as an appendix.

Badavi, F. F.↗

Modeling of multi-rotor torsional vibrations in rotating machinery using substructuring

The application of FEM modeling techniques to the analysis of torsional vibrations in complex rotating systems is described and demonstrated, summarizing results reported by Soares (1985). A substructuring approach is used for determination of torsional natural frequencies and resonant-mode shapes, steady-state frequency-sweep analysis, identification of dynamically unstable speed ranges, and characterization of transient linear and nonlinear systems. Results for several sample problems are presented in diagrams, graphs, and tables. STORV, a computer code based on this approach, is in use as a preliminary design tool for drive-train torsional analysis in the High Altitude Wind Tunnel at NASA Lewis.

Soares, Fola R.↗

An approach to probabilistic finite element analysis using a mixed-iterative formulation

An efficient algorithm for computing the response sensitivity of finite element problems based on a mixed-iterative formulation is proposed. This method does not involve explicit differentiation of the tangent stiffness array and can be used with formulations for which a consistent tangent stiffness is not readily available. The method has been successfully applied to probabilistic finite element analysis of problems using the proposed mixed formulation, and this exercise has provided valuable insights regarding the extension of the method to a more general class of problems to include material and geometric nonlinearities.

Dias, J. B.↗

Analytic theory for the selection of Saffman-Taylor fingers in the presence of thin film effects

An analytic theory is presented for the width selection of Saffman-Taylor fingers in the presence of thin film effect. In the limit of small capillary number Ca and small gap to width ratio epsilon, such that epsilon much less than Ca much less than 1, it is found that fingers with relative width lambda less than 1/2 are possible such that lambda squared (1-lambda)/(1-2 lambda) = k epsilon/Ca(sup 3/2), where the positive constant k depends on the branch of solution and equals 2.776 for the first branch. A fully nonlinear analysis is necessary in this problem even to obtain the correct scaling law. It is also shown how in principle, the selection rule for arbitrary Ca can be obtained.

Tanveer, S.↗

Optimal fabrication processes for unidirectional metal-matrix composites: A computational simulation

A method is proposed for optimizing the fabrication process of unidirectional metal matrix composites. The temperature and pressure histories are optimized such that the residual microstresses of the composite at the end of the fabrication process are minimized and the material integrity throughout the process is ensured. The response of the composite during the fabrication is simulated based on a nonlinear micromechanics theory. The optimal fabrication problem is formulated and solved with non-linear programming. Application cases regarding the optimization of the fabrication cool-down phases of unidirectional ultra-high modulus graphite/copper and silicon carbide/titanium composites are presented.

Saravanos, D. A.↗

Plasma density enhancements created by the ionization of the Earth's upper atmosphere by artificial electron beams

Analytical calculations and experimental observations relating to the interaction with the Earth's upper atmosphere of electron beams emitted from low altitude spacecraft are presented. The problem is described by two coupled nonlinear differential equations in the up-going (along a magnetic field line) and down-going differential energy flux. The equations are solved numerically, using the MSIS atmospheric model and the IRI ionospheric model. The results form the model compare well with recent observations from the CHARGE 2 sounding rocket experiment. Two aspects of the beam-neutral atmosphere interaction are discussed. First, the limits on the electron beam current that can be emitted from a spacecraft without substantial spacecraft charging are investigated. This is important because the charging of the spacecraft to positive potentials limits the current and the escape energy of the beam electrons and thereby limits the ionization of the neutral atmosphere. As an example, we find from CHARGE 2 observations and from the model calculations that below about 180 km, secondary electrons generated through the ionization of the neutral atmosphere by 1 to 10 keV electron beams from sounding rockets, completely balance the beam current, thereby allowing the emission of very high beam currents. Second, the amount of plasma production in the beam-streak is discussed. Results are shown for selected values of the beam energy, spacecraft velocity, and spacecraft altitude.

Neubert, Torsten↗

A technique for locating function roots and for satisfying equality constraints in optimization

A new technique for locating simultaneous roots of a set of functions is described. The technique is based on the property of the Kreisselmeier-Steinhauser function which descends to a minimum at each root location. It is shown that the ensuing algorithm may be merged into any nonlinear programming method for solving optimization problems with equality constraints.

Sobieszczanski-Sobieski, Jaroslaw↗

Analytic theory for the selection of Saffman-Taylor fingers in the presence of thin film effects

The present analytic theory for the width selection of Saffman-Taylor (1958) fingers in the presence of the thin film effect establishes that, in the limit of a small capillary number and a small gap-to-width ratio, fingers whose relative width is smaller than 1/2 are possible. It is established that a fully nonlinear analysis is required for this problem in order to obtain even the correct (and rather preliminary) scaling law. The way in which the selection rule for arbitrary small capillary number is obtainable is also presented.

Tanveer, S.↗

A two-time-scale autopilot for high-performance aircraft

A two-time-scale autopilot is proposed for the Aircraft Controls Design Challenge problem. This control law uses a nonlinear aircraft model constructed from the given vehicle simulation. The vehicle model is partitioned into slow translational dynamics and fast rotational dynamics. Feedback linearization is then employed to synthesize control laws for these two-time scales. Due to the nature of the synthesis, the control law is suitable for automatic trajectory following, and also for pilot control.

Menon, P. K. A.↗

Projection methods for the numerical solution of Markov chain models

Projection methods for computing stationary probability distributions for Markov chain models are presented. A general projection method is a method which seeks an approximation from a subspace of small dimension to the original problem. Thus, the original matrix problem of size N is approximated by one of dimension m, typically much smaller than N. A particularly successful class of methods based on this principle is that of Krylov subspace methods which utilize subspaces of the form span(v,av,...,A(exp m-1)v). These methods are effective in solving linear systems and eigenvalue problems (Lanczos, Arnoldi,...) as well as nonlinear equations. They can be combined with more traditional iterative methods such as successive overrelaxation, symmetric successive overrelaxation, or with incomplete factorization methods to enhance convergence.

Saad, Youcef↗

Periodic trim solutions with HP-version finite elements in time

Finite Element in Time has been proven to be a powerful alternative solving strategy for the rotor craft trim problem. Additionally, Finite Element Method in Time has been developed in various versions like time-marching framework, Galerkin framework, Rayleigh-Ritz framework, and mixed formulation. Recently, this method was applied to the rotorcraft trim problem to obtain linearized solutions. The rotorcraft trim problem consists of trying to find a period solution for period-coefficient, differential equations subject to side constraints where certain force and momentum balance equations are forced to be equal to zero. There are free (or trim) parameters that are chosen to meet these side constraints. This project aims at expanding the application, in terms of the rotorcraft trim problem, from a linearized solution to nonlinear solution.

Source record↗