Search NASA⌕ Search

SEARCH · Search NASA

Results for “Matrix differential 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 217 records · Page 12

Nonlinear Equations of Motion for Cantilever Rotor Blades in Hover with Pitch Link Flexibility, Twist, Precone, Droop, Sweep, Torque Offset, and Blade Root Offset

Nonlinear equations of motion for a cantilever rotor blade are derived for the hovering flight condition. The blade is assumed to have twist, precone, droop, sweep, torque offset and blade root offset, and the elastic axis and the axes of center of mass, tension, and aerodynamic center coincident at the quarter chord. The blade is cantilevered in bending, but has a torsional root spring to simulate pitch link flexibility. Aerodynamic forces acting on the blade are derived from strip theory based on quasi-steady two-dimensional airfoil theory. The equations are hybrid, consisting of one integro-differential equation for root torsion and three integro-partial differential equations for flatwise and chordwise bending and elastic torsion. The equations are specialized for a uniform blade and reduced to nonlinear ordinary differential equations by Galerkin's method. They are linearized for small perturbation motions about the equilibrium operating condition. Modal analysis leads to formulation of a standard eigenvalue problem where the elements of the stability matrix depend on the solution of the equilibrium equations. Two different forms of the root torsion equation are derived that yield virtually identical numerical results. This provides a reasonable check for the accuracy of the equations.

Hodges, D. H.↗

Integral-matrix procedure for boundary-layer problems

Program, BLIMP, provides fast, highly accurate solution to general class of gas-phase boundary layer flow problems encompassing broad range of boundary conditions. Program is capable of obtaining accurate and economical solutions to governing differential equations of momentum, energy, and species.

Gross, K. W.↗

Fast methods to numerically integrate the Reynolds equation for gas fluid films

The alternating direction implicit (ADI) method is adopted, modified, and applied to the Reynolds equation for thin, gas fluid films. An efficient code is developed to predict both the steady-state and dynamic performance of an aerodynamic journal bearing. An alternative approach is shown for hybrid journal gas bearings by using Liebmann's iterative solution (LIS) for elliptic partial differential equations. The results are compared with known design criteria from experimental data. The developed methods show good accuracy and very short computer running time in comparison with methods based on an inverting of a matrix. The computer codes need a small amount of memory and can be run on either personal computers or on mainframe systems.

Dimofte, Florin↗

General Solution for Theoretical Packet Data Loss Rate

Communications systems which transfer blocks ("frames") of data must use a marker ("frame synchronization pattern") for identifying where a block begins. A technique ("frame synchronization strategy") is used to locate the start of each frame and maintain synchronization as additional blocks are processed. A device which strips out the frame synchronization pattern [FSP] and provides an "end of frame" pulse is called a frame synchronizer. As clock and data errors are introduced into the system, the start-of-block marker becomes displaced and/or corrupted. The capability of the frame synchronizer to stay locked to the pattern under these conditions is a figure of merit for the frame synchronization strategy. It is important to select a strategy which will stay locked nearly all the time at bit error rates where the data is usable. ("Bit error rate" [BER] is the fraction of binary bits which are inverted by passage through a communication system.) The fraction of frames that are discarded because the frame synchronizer is not locked is called "Percent Data Loss" or "Packet Data Loss rate" [PDL]. A general approach for accurately predicting PDL given BER was developed in Theoretical Percent Data Loss Calculation and Measurement Accuracy, T. P. Kelly, LESC-30554, December 1992. Kelly gave a solution in terms of matrix equations, and only addressed "level" channel encoding. This paper goes on to give a closed-form polynomial solution for the most common class of frame synchronizer strategies, and will also address "mark" and "space" (differential) channel encoding, and burst error environments. The paper is divided into four sections and follows a logically ordered presentation, with results developed before they are evaluated. However, most readers will derive the greatest benefit from this paper by treating the results as reference material. The result developed for differential encoding can be extended to other applications (like block codes) where the probability is needed that a block contains only a certain number of errors.

Lansdowne, Chatwin↗

Estimation and detection of signals in multiplicative noise

We define a class of detection-estimation problems on matrix Lie groups in which the observation noise is multiplicative in nature. By examining the differential versions of the hypotheses, which are bilinear, we are able to derive the relevant likelihood ratio formula and the associated optimal estimation equations for the signal given the observations and the assumption that the signal is present. These estimation equations are of interest in their own right, in that they represent a finite-dimensional optimal solution to a nonlinear estimation problem and consist of a Kalman-Bucy filter along with the on-line computation of the solution of the associated Riccati equation, which is driven by the observations. The usefulness of these results is illustrated via an example concerning the detection of an actuator failure in a rigid-body rotational control system.

Willsky, A. S.↗

Enhancing photoionization rate calculations in low-temperature plasmas using spectral methods

Photoionization plays a central role in the development of streamer discharges and other non-equilibrium plasma phenomena. It creates seed electrons, which are essential for positive streamer propagation, allowing the ionization front to move forward. Because of this, accurate modeling of photoionization is very important for predicting streamer behavior and plasma evolution. The photoionization process in air (N 2 – O 2 mixture) is often described by the Zheleznyak model (1982). This model is usually solved through Helmholtz-type equations that approximate the Zheleznyak photoionization model (Zheleznyak et al. 1982) as Partial Differential Equations (PDEs). Conventional numerical methods, such as the Finite Difference Method (FDM) or Finite Volume Method (FVM), are widely used to solve these equations. Although they are prevalent, the computational cost of these methods due to their need for matrix operations and iterative solver is demanding. To address this challenge, this work develops a spectral solver based on the Fast Fourier Transform (FFT) combined with Discrete Cosine Transform (DCT) and Discrete Sine Transform (DST) to calculate the photoionization rate efficiently in an axisymmetric cylindrical domain. This method naturally satisfies the boundary conditions used in the model and converts the PDE into algebraic ones in spectral space. Thus, avoids the need for iterative matrix solvers. When compared with FDM results, it is demonstrated that the new solver not only maintains accuracy, but also reduces the computational cost, showing a performance increase of approximately 100 compared to FDM over a wide range of problem sizes. The method is parallelized using Message Passing Interface (MPI) and has been integrated into a fluid plasma model for streamer simulation. Here, this FFT-based approach provides a fast and reliable alternative for calculating photoionization in fluid models, helping large-scale plasma simulations run faster and efficiently, and allows higher-resolution simulation without extra computational cost.

Axisymmetric system↗

A generalized analysis of dual-polarization radar measurements of rain

In the study of meteorological phenomena by radar, it has been shown that dual-polarization radar measurements can yield useful data about precipitation structures within the radar beam. The techniques available for dual-polarized radar measurement are differential reflectivity, or Z(DR), linear depolarization ratio (LDR), and circular depolarization ratio (CDR). This paper presents a theoretical analysis of backscattering from a volume filled with raindrops using Rayleigh scattering theory. A generalized radar equation is given in matrix form, and the dependence of backscatter measurements in terms of a rainstorm model including raindrop shape, size, and canting angle is examined. Calculated results which show the dependence of Z(DR) on beam elevation angle and polarization angle are presented. Some results for LDR and CDR are presented, showing how these quantities depend on raindrop parameters.

Stapor, D. P.↗

ORACLS: A system for linear-quadratic-Gaussian control law design

A modern control theory design package (ORACLS) for constructing controllers and optimal filters for systems modeled by linear time-invariant differential or difference equations is described. Numerical linear-algebra procedures are used to implement the linear-quadratic-Gaussian (LQG) methodology of modern control theory. Algorithms are included for computing eigensystems of real matrices, the relative stability of a matrix, factored forms for nonnegative definite matrices, the solutions and least squares approximations to the solutions of certain linear matrix algebraic equations, the controllability properties of a linear time-invariant system, and the steady state covariance matrix of an open-loop stable system forced by white noise. Subroutines are provided for solving both the continuous and discrete optimal linear regulator problems with noise free measurements and the sampled-data optimal linear regulator problem. For measurement noise, duality theory and the optimal regulator algorithms are used to solve the continuous and discrete Kalman-Bucy filter problems. Subroutines are also included which give control laws causing the output of a system to track the output of a prescribed model.

Armstrong, E. S.↗

Vibrational and rotational transitions in low-energy electron-diatomic-molecule collisions. I - Close-coupling theory in the moving body-fixed frame. II - Hybrid theory and close-coupling theory: An /l subscript z-prime/-conserving close-coupling approximation

A detailed vibrational-rotational (V-R) close-coupling formulation of electron-diatomic-molecule scattering is developed in which the target molecular axis is chosen to be the z-axis and the resulting coupled differential equation is solved in the moving body-fixed frame throughout the entire interaction region. The coupled differential equation and asymptotic boundary conditions in the body-fixed frame are given for each parity, and procedures are outlined for evaluating V-R transition cross sections on the basis of the body-fixed transition and reactance matrix elements. Conditions are discussed for obtaining identical results from the space-fixed and body-fixed formulations in the case where a finite truncated basis set is used. The hybrid theory of Chandra and Temkin (1976) is then reformulated, relevant expressions and formulas for the simultaneous V-R transitions of the hybrid theory are obtained in the same forms as those of the V-R close-coupling theory, and distorted-wave Born-approximation expressions for the cross sections of the hybrid theory are presented. A close-coupling approximation that conserves the internuclear axis component of the incident electronic angular momentum (l subscript z-prime) is derived from the V-R close-coupling formulation in the moving body-fixed frame.

Choi, B. H.↗

High precision computing with charge domain devices and a pseudo-spectral method therefor

The present invention enhances the bit resolution of a CCD/CID MVM processor by storing each bit of each matrix element as a separate CCD charge packet. The bits of each input vector are separately multiplied by each bit of each matrix element in massive parallelism and the resulting products are combined appropriately to synthesize the correct product. In another aspect of the invention, such arrays are employed in a pseudo-spectral method of the invention, in which partial differential equations are solved by expressing each derivative analytically as matrices, and the state function is updated at each computation cycle by multiplying it by the matrices. The matrices are treated as synaptic arrays of a neural network and the state function vector elements are treated as neurons. In a further aspect of the invention, moving target detection is performed by driving the soliton equation with a vector of detector outputs. The neural architecture consists of two synaptic arrays corresponding to the two differential terms of the soliton-equation and an adder connected to the output thereof and to the output of the detector array to drive the soliton equation.

Barhen, Jacob↗

Robot arm dynamic model reduction for control

General methods are described by which the mathematical complexities of explicit and exact state equations of robot arms can be reduced to a simplified and compact state equation representation without introducing significant errors into the robot arm dynamic model. The model reduction methods are based on homogeneous coordinates and on the Langrangian algorithm for robot arm dynamics, and utilize matrix, vector and numeric analysis techniques. The derivation of differential vector representation of centripetal and Coriolis forces which has not yet been established in the literature is presented.

Bejczy, A. K.↗

Linear and nonlinear dynamic analysis of redundant load path bearingless rotor systems

The goal of this research is to develop the transfer matrix method to treat nonlinear autonomous boundary value problems with multiple branches. The application is the complete nonlinear aeroelastic analysis of multiple-branched rotor blades. Once the development is complete, it can be incorporated into the existing transfer matrix analyses. There are several difficulties to be overcome in reaching this objective. The conventional transfer matrix method is limited in that it is applicable only to linear branch chain-like structures, but consideration of multiple branch modeling is important for bearingless rotors. Also, hingeless and bearingless rotor blade dynamic characteristics (particularly their aeroelasticity problems) are inherently nonlinear. The nonlinear equations of motion and the multiple-branched boundary value problem are treated together using a direct transfer matrix method. First, the formulation is applied to a nonlinear single-branch blade to validate the nonlinear portion of the formulation. The nonlinear system of equations is iteratively solved using a form of Newton-Raphson iteration scheme developed for differential equations of continuous systems. The formulation is then applied to determine the nonlinear steady state trim and aeroelastic stability of a rotor blade in hover with two branches at the root. A comprehensive computer program is developed and is used to obtain numerical results for the (1) free vibration, (2) nonlinearly deformed steady state, (3) free vibration about the nonlinearly deformed steady state, and (4) aeroelastic stability tasks. The numerical results obtained by the present method agree with results from other methods.

Murthy, V. R.↗

Coilgun Acceleration Model Containing Multiple Interacting Coils

A coilgun operates by pulsing current through an axially-arranged series of independently-controlled coils inductively interacting with a small, electrically-conductive, azimuthally-symmetric projectile to accelerate it to high velocities. The electrical circuits are programmed to pulse current through the coils in such a way so as to impart further electromagnetic acceleration in each stage. A method is developed to calculate the mutual inductance between the coils and between each coil and the projectile. These terms are used to write a system of first-order ordinary differential equations governing the projectile velocity and the current flow in each coil. While the inclusion of the electromagnetic interactions between coils significantly complicates the equation set as more coil sets are included in the problem, casting the problem symbolically in mass matrix form permits solution using standard numerical Runge-Kutta techniques. Comparing a projectile with a single-turn to that comprised of nine-turns, the inductance of the former is much smaller, but this leads to a greater induced projectile current. The lower inductance and greater current appear to offset each other with little difference in the acceleration profile for the two cases. For the limited cases studied, coils with a discharge half-cycle equal to the time for a projectile to transit from one coil to the next yield increased efficiency.

Polzin, Kurt A.↗

Estimation and detection of signals in multiplicative noise

A class of detection-estimation problems on matrix Lie groups is defined in which the observation noise is multiplicative in nature. By examining the differential versions of the hypotheses, which are bilinear in nature, it is possible to derive the relevant likelihood ratio formula and the associated optimal estimation equations for the signal given the observations and the assumption that the signal is present. These estimation equations are of interest in their own right, in that they represent a finite dimensional optimal solution to a nonlinear estimation problem and can be viewed as consisting of a Kalman-Bucy filter along with the on-line computation of the solution of the associated Riccati equation, which is driven by the observations. The usefulness of these results is illustrated via an example concerning the detection of an actuator failure in a rigid body rotational control system.

Willsky, A. S.↗

On substructuring algorithms and solution techniques for the numerical approximation of partial differential equations

Substructuring methods are in common use in mechanics problems where typically the associated linear systems of algebraic equations are positive definite. Here these methods are extended to problems which lead to nonpositive definite, nonsymmetric matrices. The extension is based on an algorithm which carries out the block Gauss elimination procedure without the need for interchanges even when a pivot matrix is singular. Examples are provided wherein the method is used in connection with finite element solutions of the stationary Stokes equations and the Helmholtz equation, and dual methods for second-order elliptic equations.

Gunzburger, M. D.↗

A finite element algorithm for sound propagation in axisymmetric ducts containing compressible mean flow

The described investigation is concerned with he development of a finite element scheme which can be used in a study of the acoustics of aircraft-engine ducts. In the absence of suitable variational principles for acoustic fluctuations within an aircraft fan engine, an acoustic analysis must proceed directly from the differential equations which describe compressible flow. The derived equations cannot be solved algebraically. The numerical technique used for solving them makes use of a linear rectangular element of a type considered by Zienkiewicz (1971). Attention is given to aspects of element derivation, the global matrix assembly, the solution of the matrix equation, questions of acoustic attenuation, and illustrations of the potential of the current model in duct optimization.

Abrahamson, A. L.↗

A Walsh Function Module Users' Manual

The solution of partial differential equations (PDEs) with Walsh functions offers new opportunities to simulate many challenging problems in mathematical physics. The approach was developed to better simulate hypersonic flows with shocks on unstructured grids. It is unique in that integrals and derivatives are computed using simple matrix multiplication of series representations of functions without the need for divided differences. The product of any two Walsh functions is another Walsh function - a feature that radically changes an algorithm for solving PDEs. A FORTRAN module for supporting Walsh function simulations is documented. A FORTRAN code is also documented with options for solving time-dependent problems: an advection equation, a Burgers equation, and a Riemann problem. The sample problems demonstrate the usage of the Walsh function module including such features as operator overloading, Fast Walsh Transforms in multi-dimensions, and a Fast Walsh reciprocal.

Gnoffo, Peter A.↗

Implicit Formulation of Muscle Dynamics in OpenSim

Astronauts lose bone and muscle mass during spaceflight. Exercise countermeasure is the primary method for counteracting bone and muscle mass loss in space. New spacecraft exercise device concepts are currently being developed for the NASAs new crew exploration vehicle. The NASA Digital Astronaut Project (DAP) uses computational modeling to help determine if the new exercise devices will be effective as countermeasures. The NASA Digital Astronaut Project is developing the ability to utilize predictive simulation to provide insight into the change in kinematics and kinetics with a change in device and gravitational environment (1-g versus 0-g). For example, in space exercise the subject's body weight is applied in addition to the loads prescribed for musculoskeletal maintenance. How and where these loads are applied obviously directly impacts bone and tissue loads. Additionally, due to space vehicle structural requirements, exercise devices are often placed on vibration isolation systems. This changes the apparent impedance or stiffness of the device as seen by the user. Data collection under these conditions is often impractical and limited. Predictive modeling provides a means to have a virtual subject to test hypotheses. Predictive simulation provides a virtual subject for which we are able to perform studies such as sensitivity to device loading and vibration isolation without the need for laboratory kinematic or kinetic test data.Direct Collocation optimization provides an efficient means to perform task based optimization and predictive modeling. It is relatively straight forward to structure a physical exercise task in a Direct Collocation mathematical formulation: perform a motion such that you start at an initial pose, achieve a given amount of deflection i.e a squat, return to the initial pose, and minimize muscle activation cost. Direct Collocation is advantageous in that it does not require numerical integration to evaluate the objective function. Instead, the system dynamics are transformed to discrete time and the optimizer is constrained such that the solution is not considered to be a valid unless the dynamic equations are satisfied at all time points. The simulation and optimization are effectively done simultaneously. Due to the implicit integration, time steps can be more coarse than in a differential equation solver. In a gait scenario this means that that the model constraints and cost function are evaluated at 100 nodes in the gait cycle versus 10,000 integration steps in a variable-step forward dynamic simulation. Furthermore, no time is wasted on accurate simulations of movements that are far from the optimum. Constrained optimization algorithms require a Jacobian matrix that contains the partial derivatives of each of the dynamic constraints with respect to of each of the state and control variables at all time points. This is a large but sparse matrix. An implicit dynamics formulation requires computation of the dynamic residuals f as a function of the states x and their derivatives, and controls u:f(x, dxdt, u) 0If the dynamics of musculoskeletal system are formulated implicitly, the Jacobian elements are often available analytically, eliminating the need for numerical differentiation; this is obviously computationally advantageous. Additionally, implicit formulation of musculoskeletal dynamics do not suffer from singularities from low mass bodies, zero muscle activation, or other stiff system or

physical exercise↗