Engineering PapersSearch

SEARCH · Engineering Papers

Results for “CHARACTERISTIC EQUATION”

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 19 records

A high order accurate finite element algorithm for high Reynolds number flow prediction

A Galerkin-weighted residuals formulation is employed to establish an implicit finite element solution algorithm for generally nonlinear initial-boundary value problems. Solution accuracy, and convergence rate with discretization refinement, are quantized in several error norms, by a systematic study of numerical solutions to several nonlinear parabolic and a hyperbolic partial differential equation characteristic of the equations governing fluid flows. Solutions are generated using selective linear, quadratic and cubic basis functions. Richardson extrapolation is employed to generate a higher-order accurate solution to facilitate isolation of truncation error in all norms. Extension of the mathematical theory underlying accuracy and convergence concepts for linear elliptic equations is predicted for equations characteristic of laminar and turbulent fluid flows at nonmodest Reynolds number. The nondiagonal initial-value matrix structure introduced by the finite element theory is determined intrinsic to improved solution accuracy and convergence. A factored Jacobian iteration algorithm is derived and evaluated to yield a consequential reduction in both computer storage and execution CPU requirements while retaining solution accuracy.

Baker, A. J.

Invariant poles feedback control of flexible highly variable spacecraft.

Description of a technique for single-axis control of a model of a highly flexible space station. Active damping of lower frequency flexibility modes is employed. In the control technique, referred to as invariant poles feedback control, feedback gains are adjusted so that the closed-loop system characteristic equation is matched to that of a reference model. Hence closed-loop system poles will not move; they will be invariant (provided that bending frequencies and parameters can be identified accurately). This is accomplished by obtaining the system characteristic equation in closed form; equating respective coefficients between terms of like powers in s in the system and reference model characteristic equations; and solving for the feedback gains. The feedback gains are explicit functions of system plant parameters and the coefficients of the reference model characteristic equation, and are easily programmed for the digital computer.

Mendel, J. M.

Invariant poles feedback control of flexible, highly variable spacecraft.

This paper describes a technique for single-axis control of a model of a highly flexible Space Station. Active damping of lower frequency flexibility modes is employed. In the control technique, referred to as invariant poles feedback control (IPFC), feedback gains are adjusted so that the closed-loop system's characteristic equation is matched to that of a reference model; hence, closed-loop system's poles will not move - they will be invariant (provided bending frequencies and parameters can be identified accurately). This is accomplished by obtaining the system's characteristic equation in closed form; equating respective coefficients between terms of like powers in s in the system and reference model characteristic equations; and, solving for the feedback gains. The feedback gains are explicit functions of system plant parameters and the coefficients of the reference model's characteristic equation, and are easily programmed for the digital computer.

Mendel, J. M.

Two bonded half planes with a crack going through the interface

The plane problem of two bonded elastic half planes containing a finite crack perpendicular to and going through the interface is considered. The problem is formulated as a system of singular integral equations with generalized Cauchy kernels. Even though the system has three irregular points, it is shown that the unknown functions are algebraically related at the irregular point on the interface and the integral equations can be solved by a method developed previously. The system of integral equations is shown to yield the same characteristic equation as that for two bonded quarter planes in the general case of the through crack, and the characteristic equation for a crack tip terminating at the interface in the special case. The numerical results given in the paper include the stress intensity factors at the crack tips, the normal and shear components of the stress intensity factors at the singular point on the interface, and the crack surface displacements.

Erdogan, F.

Anisotropic elliptic optical fibers

The exact characteristic equation for an anisotropic elliptic optical fiber is obtained for odd and even hybrid modes in terms of infinite determinants utilizing Mathieu and modified Mathieu functions. A simplified characteristic equation is obtained by applying the weakly guiding approximation such that the difference in the refractive indices of the core and the cladding is small. The simplified characteristic equation is used to compute the normalized guide wavelength for an elliptical fiber. When the anisotropic parameter is equal to unity, the results are compared with the previous research and they are in close agreement. For a fixed value normalized cross-section area or major axis, the normalized guide wavelength lambda/lambda(sub 0) for an anisotropic elliptic fiber is small for the larger value of anisotropy. This condition indicates that more energy is carried inside of the fiber. However, the geometry and anisotropy of the fiber have a smaller effect when the normalized cross-section area is very small or very large.

Kang, Soon Ahm

The three dimensional motion and stability of a rotating space station-cable-counterweight configuration

The three dimensional equations of motion for a cable connected space station-counterweight system are developed using a Lagrangian formulation. The system model employed allows for cable and end body damping and restoring effects. The equations are then linearized about the equilibrium motion and nondimensionalized. To first degree, the out-of-plane equations uncouple from the in-plane equations. From the general in-plane characteristic equation, necessary conditions for stability are obtained. The Routh-Hurwitz necessary and sufficient conditions for stability are derived for the general out-of-plane characteristic equation. Special cases of the in-plane and out-of-plane equations are then examined for stability. Time constants for the least damped mode are obtained for a range of system parameters by numerical examination of the roots of the in-plane and out-of-plane characteristic polynomials.

Bainum, P. M.

Topological Control of Triply Periodic Minimal Surfaces for Thermal Design and Advanced Manufacturing: A Gyroid Case Study

Recently, there has been a heightened interest in using triply periodic minimal surfaces (TPMSs) in the design of compact process engineering components. The benefits of high surface area per unit volume, modular form, and inherent periodicity provide a holistic self-supporting network and flow-conducive features. Applications of importance include thermal power management, biomimetic scaffolds and structures, and feasibility of advanced manufacturing. This study presents a novel approach to the manipulation of the characteristic Schwarz-G, or gyroid TPMS, for thermal design in the context of advanced manufacturing. The study presents relationships between design parameters and resulting surface area as a target response using the characteristic equation of a gyroid. Through parametric control, the characteristic equation is manipulated to produce a 20-fold increase in achievable area over a baseline design characteristic of 25.4 mm through controlled combinations of design parameters. A second relationship is presented as a function of the maximum area achieved and manipulated design parameters. Through the analysis, the study presents a framework to identify and maximize the achievable area of TPMSs for advanced manufacturing and thermal management applications.

gyroid

Generalized characteristics method for elastic wave propagation problems

Characteristic equations are derived in generalized curvilinear coordinates. Linear elastic, isotropic, and homogeneous constitutive equations have been used in the derivation. The generalized characteristic equations readily lend themselves to any requirements of space dimension and geometry. A simple boundary value problem is solved to indicate the applicability of these equations.

Ziv, M.

Unsteady density-current equations for highly curved terrain

New nonlinear partial differential equations containing terrain curvature and its rate of change are derived that describe the flow of an atmospheric density current. Unlike the classical hydraulic-type equations for density currents, the new equations are valid for two-dimensional, gradually varied flow over highly curved terrain, hence suitable for computing unsteady (or steady) flows over arbitrary mountain/valley profiles. The model assumes the atmosphere above the density current exerts a known arbitrary variable pressure upon the unknown interface. Later this is specialized to the varying hydrostatic pressure of the atmosphere above. The new equations yield the variable velocity distribution, the interface position, and the pressure distribution that contains a centrifugal component, often significantly larger than its hydrostatic component. These partial differential equations are hyperbolic, and the characteristic equations and characteristic directions are derived. Using these to form a characteristic mesh, a hypothetical unsteady curved-flow problem is calculated, not based upon observed data, merely as an example to illustrate the simplicity of their application to unsteady flows over mountains.

Sivakumaran, N. S.

The Physics Imposed on a Streaming Operator by Spherical Transport Problems

The streaming operator, which generates a displacement of a particle on a straight line at a constant speed in transport theory, is derived algebraically from a spherical coordinate formulation of Newton’s second law. This derivation leads to an operator that has more partial derivatives than a Cartesian coordinate formulation of the operator. The additional partial derivatives, which are with respect to the normalized velocity variables of a particle, take into account the intrinsic curvature of a ball. Moreover, these partial derivatives mitigate ray effects, which arise when a finite number of normalized velocities (also called directions or discrete ordinates) are used to simulate a continuous S 2 sphere of directions, by rotating the polar axis of the S 2 sphere into the radial direction of the coordinate system. As a consequence of this rotation, the number of actual discrete ordinates is greatly amplified to an enormous number of effective discrete ordinates by a multiplier that is equal to the number of patches that partitions a spherical surface. In addition to the derivation of the streaming operator, we provide in closed form a solution to the system of characteristic equations that is equivalent to the streaming operator. Furthermore, the solution to the system of characteristic equations enables the construction of an integral operator that is the inverse to the streaming operator. Examples in which ray effects are immensely mitigated by spherical coordinates are presented.

integral operator

A comparative study of nonreflecting far-field boundary condition procedures for unsteady transonic flow computation

Various nonreflecting far-field boundary condition procedures are compared by implementing them in the computer code LTRAN2. This code solves the implicit finite-difference representation of the small-disturbance equations for transonic flows about airfoils. The first- and second-approximate nonreflecting conditions, as proposed by Engquist and Majda, are compared with the condition derived from the full-characteristic equation. The far-field boundary conditions and the description of the algorithm for implementing these conditions in LTRAN2 are discussed. Various cases are computed and compared with results from the older, more conventional procedures. One concludes that the full-characteristic equation produces the most effective results, thus allowing the far-field boundary to be located closer to the airfoil; this decreases the computer time required to obtain the solution because fewer mesh points are required.

Kwak, D.

Transverse vibration and buckling of a cantilevered beam with tip body under constant axial base acceleration

The planar transverse bending behavior of a uniform cantilevered beam with rigid tip body subject to constant axial base acceleration was analyzed. The beam is inextensible and capable of small elastic transverse bending deformations only. Two classes of tip bodies are recognized: (1) mass centers located along the beam tip tangent line; and (2) mass centers with arbitrary offset towards the beam attachment point. The steady state response is studied for the beam end condition cases: free, tip mass, tip body with restricted mass center offset, and tip body with arbitrary mass center offset. The first three cases constitute classical Euler buckling problems, and the characteristic equation for the critical loads/accelerations are determined. For the last case a unique steady state solution exists. The free vibration response is examined for the two classes of tip body. The characteristic equation, eigenfunctions and their orthogonality properties are obtained for the case of restricted mass center offset. The vibration problem is nonhomogeneous for the case of arbitrary mass center offset. The exact solution is obtained as a sum of the steady state solution and a superposition of simple harmonic motions.

Storch, J.

Design implementation in model-reference adaptive systems

The derivation of an approximate error characteristic equation describing the transient system error response is given, along with a procedure for selecting adaptive gain parameters so as to relate to the transient error response. A detailed example of the application and implementation of these methods for a space shuttle type vehicle is included. An extension of the characteristic equation technique is used to provide an estimate of the magnitude of the maximum system error and an estimate of the time of occurrence of this maximum after a plant parameter disturbance. Techniques for relaxing certain stability requirements and the conditions under which this can be done and still guarantee asymptotic stability of the system error are discussed. Such conditions are possible because the Lyapunov methods used in the stability derivation allow for overconstraining a problem in the process of insuring stability.

Boland, J. S., III

A feedback control model for network flow with multiple pure time delays

A control model describing a network flow hindered by multiple pure time (or transport) delays is formulated. Feedbacks connect each desired output with a single control sector situated at the origin. The dynamic formulation invokes the use of differential difference equations. This causes the characteristic equation of the model to consist of transcendental functions instead of a common algebraic polynomial. A general graphical criterion is developed to evaluate the stability of such a problem. A digital computer simulation confirms the validity of such criterion. An optimal decision making process with multiple delays is presented.

Press, J.