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 595 records · Page 33

Parameter adaptive estimation of random processes

This paper is concerned with the parameter adaptive least squares estimation of random processes. The main result is a general representation theorem for the conditional expectation of a random variable on a product probability space. Using this theorem along with the general likelihood ratio expression, the least squares estimate of the process is found in terms of the parameter conditioned estimates. The stochastic differential for the a posteriori probability and the stochastic differential equation for the a posteriori density are found by using simple stochastic calculus on the representations obtained. The results are specialized to the case when the parameter has a discrete distribution. The results can be used to construct an implementable recursive estimator for certain types of nonlinear filtering problems. This is illustrated by some simple examples.

Caglayan, A. K.↗

Generation of linear dynamic models from a digital nonlinear simulation

The results and methodology used to derive linear models from a nonlinear simulation are presented. It is shown that averaged positive and negative perturbations in the state variables can reduce numerical errors in finite difference, partial derivative approximations and, in the control inputs, can better approximate the system response in both directions about the operating point. Both explicit and implicit formulations are addressed. Linear models are derived for the F 100 engine, and comparisons of transients are made with the nonlinear simulation. The problem of startup transients in the nonlinear simulation in making these comparisons is addressed. Also, reduction of the linear models is investigated using the modal and normal techniques. Reduced-order models of the F 100 are derived and compared with the full-state models.

Daniele, C. J.↗

Research on an augmented Lagrangian penalty function algorithm for nonlinear programming

The augmented Lagrangian (ALAG) Penalty Function Algorithm for optimizing nonlinear mathematical models is discussed. The mathematical models of interest are deterministic in nature and finite dimensional optimization is assumed. A detailed review of penalty function techniques in general and the ALAG technique in particular is presented. Numerical experiments are conducted utilizing a number of nonlinear optimization problems to identify an efficient ALAG Penalty Function Technique for computer implementation.

Frair, L.↗

Nonlinear singularly perturbed optimal control problems with singular arcs

Singular perturbation techniques are studied for dealing with singular arc problems by analyzing a relatively low-order but otherwise general system. This system encompasses many flight mechanic problems including Goddard's problem and a version of the minimum time-to-climb problem. Boundary layer solutions are constructed which are stable and reach the outer solution in a finite time. A uniformly valid composite solution is then formed from the reduced and boundary layer solutions. The value of the approximate solution is that it is relatively easy to obtain and does not involve singular arcs. To illustrate the utility of the results, the technique is used to obtain an approximate solution of a simplified version of the aircraft minimum time-to-climb problem.

Ardema, M. D.↗

Profiles of electrified drops and bubbles

Axisymmetric equilibrium shapes of conducting drops and bubbles, (1) pendant or sessile on one face of a circular parallel-plate capacitor or (2) free and surface-charged, are found by solving simultaneously the free boundary problem consisting of the augmented Young-Laplace equation for surface shape and the Laplace equation for electrostatic field, given the surface potential. The problem is nonlinear and the method is a finite element algorithm employing Newton iteration, a modified frontal solver, and triangular as well as quadrilateral tessellations of the domain exterior to the drop in order to facilitate refined analysis of sharply curved drop tips seen in experiments. The stability limit predicted by this computer-aided theoretical analysis agrees well with experiments.

Basaran, O. A.↗

Utilization of normal mode initial conditions for detecting errors in the dynamics part of primitive equation global models

When a global atmospheric basic state has constant angular velocity and its temperature varies with altitude only, there exist normal mode solutions to the linearized global primitive equations. The use of these normal modes, which have known behavior in time, is superior to the use of the Rossby-Haurwitz wave as initial conditions for detecting errors in the dynamics part of primitive equation global models. With these initial conditions, integration through only one time step is sufficient to detect many formulation and coding errors. Other tests are still required for detecting problems of nonlinear instability and conservation of integral properties, however.

Chao, W. C.↗

Nonlinear optimization with linear constraints using a projection method

Nonlinear optimization problems that are encountered in science and industry are examined. A method of projecting the gradient vector onto a set of linear contraints is developed, and a program that uses this method is presented. The algorithm that generates this projection matrix is based on the Gram-Schmidt method and overcomes some of the objections to the Rosen projection method.

Fox, T.↗

Calculation of exact vibration modes for plane grillages by the dynamic stiffness method

A dynamic stiffness method is developed for the calculation of the exact modal parameters for plane grillages which consist of straight and uniform beams with coincident elastic and inertial axes. Elementary bending-torsion beam theory is utilized, and bending translation is restricted to one direction. The exact bending-torsion dynamic stiffness matrix is obtained for a straight and uniform beam element with coincident elastic and inertial axes. The element stiffness matrices are assembled using the standard procedure of the static stiffness method to form the dynamic stiffness matrix of the complete grillage. The exact natural frequencies, mode shapes, and generalized masses of the grillage are then calculated by solving a nonlinear eigenvalue problem based on the dynamic stiffness matrix. The exact modal solutions for an example grillage are calculated and compared with the approximate solutions obtained by using the finite element method.

Hallauer, W. L., Jr.↗

A refined finite element for vibration analysis of twisted blades based on beam theory

A finite element method of discretizing beam segments of pretwisted rotating blades is presented. Employing the matrix displacement method, stiffness and mass properties are developed from basic mechanics of a pretwisted beam theory. By introducing the proper displacement functions, the effect of rotor blade rotational motion on the stiffness matrix is obtained systematically from the kinetic energy expression. Comparing with other beam elements the derivation of this element is more fundamental. This allows one to apply the same approach to more complicated problems including nonlinear effects or complex dynamic motions. Illustrative examples are given comparing numerical results with available data and other numerical solutions from rotating and nonrotating force fields. These examples show that accurate prediction of vibration frequencies for pretwisted blades can be obtained by employing a quite modest number of degrees of freedom.

Sisto, F.↗

Robust penalty method for structural synthesis

The Sequential Unconstrained Minimization Technique (SUMT) offers an easy way of solving nonlinearly constrained problems. However, this algorithm frequently suffers from the need to minimize an ill-conditioned penalty function. An ill-conditioned minimization problem can be solved very effectively by posing the problem as one of integrating a system of stiff differential equations utilizing concepts from singular perturbation theory. This paper evaluates the robustness and the reliability of such a singular perturbation based SUMT algorithm on two different problems of structural optimization of widely separated scales. The report concludes that whereas conventional SUMT can be bogged down by frequent ill-conditioning, especially in large scale problems, the singular perturbation SUMT has no such difficulty in converging to very accurate solutions.

Kamat, M. P.↗

Hybrid perturbation/Bubnov-Galerkin technique for nonlinear thermal analysis

A two step hybrid analysis technique to predict the nonlinear steady state temperature distribution in structures and solids is presented. The technique is based on the regular perturbation expansion and the classical Bubnov-Galerkin approximation. The functions are obtained by using the regular perturbation method. These functions are selected as coordinate functions and the classical Bubnov-Galerkin technique is used to compute their amplitudes. The potential of the proposed hybrid technique for the solution of nonlinear thermal problems is discussed. The effectiveness of this technique is demonstrated by the effects of conduction, convection, and radiation modes of heat transfer. It is indicated that the hybrid technique overcomes the two major drawbacks of the classical techniques: (1) the requirement of using a small parameter in the regular perturbation method; and (2) the arbitrariness in the choice of the coordinate functions in the Bubnov-Galerkin technique. The proposed technique extends the range of applicability of the regular perturbation method and enhances the effectiveness of the Bubnov-Galerkin technique.

Noor, A. K.↗

The CMC/3DPNS computer program for prediction of three-dimension, subsonic, turbulent aerodynamic juncture region flow. Volume 2: Users' manual

The CMC fluid mechanics program system was developed to transmit the theoretical solution of finite element numerical solution methodology, applied to nonlinear field problems into a versatile computer code for comprehensive flow field analysis. Data procedures for the CMC 3 dimensional Parabolic Navier-Stokes (PNS) algorithm are presented. General data procedures a juncture corner flow standard test case data deck is described. A listing of the data deck and an explanation of grid generation methodology are presented. Tabulations of all commands and variables available to the user are described. These are in alphabetical order with cross reference numbers which refer to storage addresses.

Manhardt, P. D.↗

Theoretical analyses of baroclinic flows

Completed and ongoing research activities are discussed briefly, including a three-dimensional, linear stability analysis of the baroclinic Hadley cell and a numerical model of the baroclinic flow between two rotating concentric spheres. This model simulates axisymmetric flow in the Atmospheric General Circulation Experiment configuration. A computer code designed to solve the strongly nonlinear stability problem for the Eady basic state is mentioned.

Antar, B. N.↗

Equivalent source modeling of the core magnetic field using magsat data

Experiments are carried out on fitting the main field using different numbers of equivalent sources arranged in equal area at fixed radii at and inside the core-mantle boundary. In fixing the radius for a given series of runs, the convergence problems that result from the extreme nonlinearity of the problem when dipole positions are allowed to vary are avoided. Results are presented from a comparison between this approach and the standard spherical harmonic approach for modeling the main field in terms of accuracy and computational efficiency. The modeling of the main field with an equivalent dipole representation is found to be comparable to the standard spherical harmonic approach in accuracy. The 32 deg dipole density (42 dipoles) corresponds approximately to an eleventh degree/order spherical harmonic expansion (143 parameters), whereas the 21 dipole density (92 dipoles) corresponds to approximately a seventeenth degree and order expansion (323 parameters). It is pointed out that fixing the dipole positions results in rapid convergence of the dipole solutions for single-epoch models.

Mayhew, M. A.↗

Dynamic effects of combustion

The dynamic effects of combustion are due to the evolution of exothermic energy and its deposition in the compressible medium where the process takes place. The paper examines the dynamics of combustion phenomena, including ignition, turbulent flame propagation (inflammation), explosion, and detonation, with emphasis on their exothermic characteristics. Ignition and explosion are treated as problems of nonlinear mechanics, and their dynamic behavior is described in terms of phase space models and cinematographic laser shear interferograms. The results of a numerical random vortex model of turbulent flame propagation are confirmed in a combustion tunnel experiment, where it was observed that a fresh mixture of burnt and unburnt gases can sustain combustion with a relatively small expenditure of overall mass flow, due to the increasing specific volume of burnt gases inside the flame front. An isentropic pressure wave is found to precede the accelerating flame in the process of detonation, and components of this presssure wave are shown to propagate at local sonic velocities.

Oppenheim, A. K.↗

Differential eigenvalue problems in which the parameter appears nonlinearly

Several methods are examined for determining the eigenvalues of a system of equations in which the parameter appears nonlinearly. The equations are the result of the discretization of differential eigenvalue problems using a finite Chebyshev series. Two global methods are considered which determine the spectrum of eigenvalues without an initial estimate. A local iteration scheme with cubic convergence is presented. Calculations are performed for a model second order differential problem and the Orr-Sommerfeld problem for plane Poiseuille flow.

Bridges, T. J.↗

Three-Dimensional, Subsonic, Turbulent Juncture Region Flow

CMC3DPNS predicts three-dimensional, subsonic, turbulent aerodynamic juncture region flow. Versatile program for comprehensive flow-field analysis applies finite-element methodology to nonlinear field problems. Order-ofmagnitude analysis of subsonic, three-dimensional, steady time-averaged NavierStokes equations for semibounded aerodynamic juncture geometries yields para bolic Navier-Stokes equations. CMC3DPNS solves these equations to predict juncture-region flow.

Baker, A. J.↗

The CMC/3DPNS computer program for prediction of three-dimensional, subsonic, turbulent aerodynamic juncture region flow. Volume 3: Programmers' manual

The CMC fluid mechanics program system was developed to transmit the theoretical evolution of finite element numerical solution methodology, applied to nonlinear field problems into a versatile computer code for comprehensive flow field analysis. A detailed view of the code from the standpoint of a computer programmer's use is presented. A system macroflow chart and detailed flow charts of several routines necessary to interact with a theoretican/user to modify the operation of this program are presented. All subroutines and details of usage, primarily for input and output routines are described. Integer and real scalars and a cross reference list denoting subroutine usage for these scalars are outlined. Entry points in dynamic storage vector IZ; the lengths of each vector accompanying the scalar definitions are described. A listing of the routines peculiar to the standard test case and a listing of the input deck and printout for this case are included.

Orzechowski, J. A.↗