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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1979

Energy conserving norms for the solution of hyperbolic systems of partial differential equations

The hyperbolic system of partial differential equations with a real constant square coefficient matrix A is considered. The problem of finding an energy conserving norm for the solution of the system is reduced to the problem of characterizing those matrices appearing in the boundary conditions which satisfy two specific matrix equations. Necessary and sufficient conditions on the coefficient matrix A and the matrices appearing in boundary conditions are derived for an energy conserving norm. The conditions serve as criteria on a given system which determine whether or not the solution will have its energy conserved in some norm. Examples of specific systems and boundary conditions are also provided.

Gunzburger, M. D.↗

NASA NTRS1979

Sequential triangulation of orbital photography

The feasibility of structuring the satellite photogrammetric triangulation as an iterative Extended Kalman estimation algorithm is demonstrated. Comparative numerical results of the sequential against batch estimation algorithm are presented. Difficulty of accurately modeling of the attitude motion is overcome by utilizing the on-board angular rate measurements. Solutions of the differential equations and the evaluation of state transition matrix are carried out numerically.

Rajan, M.↗

NASA NTRS2010

Nonlinear Dynamics of Cantilever-Sample Interactions in Atomic Force Microscopy

The interaction of the cantilever tip of an atomic force microscope (AFM) with the sample surface is obtained by treating the cantilever and sample as independent systems coupled by a nonlinear force acting between the cantilever tip and a volume element of the sample surface. The volume element is subjected to a restoring force from the remainder of the sample that provides dynamical equilibrium for the combined systems. The model accounts for the positions on the cantilever of the cantilever tip, laser probe, and excitation force (if any) via a basis set of set of orthogonal functions that may be generalized to account for arbitrary cantilever shapes. The basis set is extended to include nonlinear cantilever modes. The model leads to a pair of coupled nonlinear differential equations that are solved analytically using a matrix iteration procedure. The effects of oscillatory excitation forces applied either to the cantilever or to the sample surface (or to both) are obtained from the solution set and applied to the to the assessment of phase and amplitude signals generated by various acoustic-atomic force microscope (A-AFM) modalities. The influence of bistable cantilever modes of on AFM signal generation is discussed. The effects on the cantilever-sample surface dynamics of subsurface features embedded in the sample that are perturbed by surface-generated oscillatory excitation forces and carried to the cantilever via wave propagation are accounted by the Bolef-Miller propagating wave model. Expressions pertaining to signal generation and image contrast in A-AFM are obtained and applied to amplitude modulation (intermittent contact) atomic force microscopy and resonant difference-frequency atomic force ultrasonic microscopy (RDF-AFUM). The influence of phase accumulation in A-AFM on image contrast is discussed, as is the effect of hard contact and maximum nonlinearity regimes of A-AFM operation.

Cantrell, John H.↗

DOE OSTI2025

Parametric matrix models

We present a general class of machine learning algorithms called parametric matrix models. In contrast with most existing machine learning models that imitate the biology of neurons, parametric matrix models use matrix equations that emulate physical systems. Similar to how physics problems are usually solved, parametric matrix models learn the governing equations that lead to the desired outputs. Parametric matrix models can be efficiently trained from empirical data, and the equations may use algebraic, differential, or integral relations. While originally designed for scientific computing, we prove that parametric matrix models are universal function approximators that can be applied to general machine learning problems. After introducing the underlying theory, we apply parametric matrix models to a series of different challenges that show their performance for a wide range of problems. For all the challenges tested here, parametric matrix models produce accurate results within an efficient and interpretable computational framework that allows for input feature extrapolation.

Computational science↗

NASA NTRS1968

Molecular collisions. VIII.

Molecular collisions, formulating description as set of differential equations and developing eigenvectors and eigenvalues of Hermitian matrix associated with coupling

Curtiss, C. F.↗

NASA NTRS1983

An analytical investigation of shape control of large space structures by applied temperatures

An analytical procedure for the static shape control of flexible space structures subjected to thermal distortions is developed which is based on prescribing temperatures in control elements having much higher coefficients of thermal expansion than the main structure. The temperatures at the control elements are defined so as to minimize the overall thermal distortion of the structure from its ideal shape, and a matrix equation is obtained which can be solved for the set of optimum control temperatures. A formulation of the procedure for continuous structures governed by differential equations and a formulation for discrete (finite element modeled) structures governed by matrix equations are presented. The equations from the continuous formulation are employed for the shape control of a simple beam distorted by nonuniform heating, and the discrete formulation is applied in a general purpose finite-element structural analysis computer program for the shape control of a 750 m radiometer antenna reflector dish subjected to orbital heating. A reduction in thermal distortion by a factor of nearly 50 was obtained with the use of only seven control elements. Results for four different sets of control locations for the antenna are presented in which reductions in distortion of up to a factor of four were obtained.

Haftka, R. T.↗

NASA NTRS1984

Acoustooptic linear algebra processors - Architectures, algorithms, and applications

Architectures, algorithms, and applications for systolic processors are described with attention to the realization of parallel algorithms on various optical systolic array processors. Systolic processors for matrices with special structure and matrices of general structure, and the realization of matrix-vector, matrix-matrix, and triple-matrix products and such architectures are described. Parallel algorithms for direct and indirect solutions to systems of linear algebraic equations and their implementation on optical systolic processors are detailed with attention to the pipelining and flow of data and operations. Parallel algorithms and their optical realization for LU and QR matrix decomposition are specifically detailed. These represent the fundamental operations necessary in the implementation of least squares, eigenvalue, and SVD solutions. Specific applications (e.g., the solution of partial differential equations, adaptive noise cancellation, and optimal control) are described to typify the use of matrix processors in modern advanced signal processing.

Casasent, D.↗

NASA NTRS1977

Comparison of some optimal control methods for the design of turbine blades

This paper attempts a comparative study of some numerical methods for the optimal control design of turbine blades whose vibration characteristics are approximated by Timoshenko beam idealizations with shear and incorporating simple boundary conditions. The blade was synthesized using the following methods: (1) conjugate gradient minimization of the system Hamiltonian in function space incorporating penalty function transformations, (2) projection operator methods in a function space which includes the frequencies of vibration and the control function, (3) epsilon-technique penalty function transformation resulting in a highly nonlinear programming problem, (4) finite difference discretization of the state equations again resulting in a nonlinear program, (5) second variation methods with complex state differential equations to include damping effects resulting in systems of inhomogeneous matrix Riccatti equations some of which are stiff, (6) quasi-linear methods based on iterative linearization of the state and adjoint equation. The paper includes a discussion of some substantial computational difficulties encountered in the implementation of these techniques together with a resume of work presently in progress using a differential dynamic programming approach.

Desilva, B. M. E.↗

NASA NTRS1972

Transient, compressible heat and mass transfer in porous media using the strongly implicit iteration procedure.

Coupled nonlinear partial differential equations describing heat and mass transfer in a porous matrix are solved in finite difference form with the aid of a new iterative technique (the strongly implicit procedure). Example numerical results demonstrate the characteristics of heat and mass transport in a porous matrix such as a charring ablator. It is emphasized that multidimensional flow must be considered when predicting the thermal response of a porous material subjected to nonuniform boundary conditions.

Curry, D. M.↗

NASA NTRS1979

Buckling of rotating beams

The stability of a beam subjected to compressive centrifugal forces arising from steady rotation about an axis which does not pass through the clamped end of the beam is analyzed to determine the critical rotational speeds for buckling in the inplane and out-of-plane directions. The differential equations of motion are solved numerically using an integrating matrix method in combination with an eigenanalysis to determine the eigenvalues from which stability is assessed. The results clarify several differences which have been identified in the literature relating to the proper behavior of the critical rotational speed for buckling as the radius of rotation of the clamped end of the beam is reduced.

White, W. F., Jr.↗

NASA NTRS1980

On time discretizations for spectral methods

New methods are introduced for the time integration of the Fourier and Chebyshev methods of solution for dynamic differential equations. These methods are unconditionally stable, even though no matrix inversions are required. Time steps are chosen by accuracy requirements alone. For the Fourier method both leapfrog and Runge-Kutta methods are considered. For the Chebyshev method only Runge-Kutta schemes are tested. Numerical calculations are presented to verify the analytic results. Applications to the shallow water equations are presented.

Gottlieb, D.↗

NASA NTRS1986

Calculations of the stability of some axisymmetric flows proposed as a model of vortex breakdown

The term vortex breakdown refers to the abrupt and drastic changes of structure that can sometimes occur in swirling flows. It was conjectured that the bubble type of breakdown can be viewed as an axisymmetric wave traveling upstream in a primarily columnar vortex flow. In this scenario the wave's upstream progress is impeded only when it reaches a critical amplitude and it loses stability to some nonaxisymmetric disturbance. The stability of some axisymmetric wavy flows to three dimensional disturbances, viewing the amplitude of the wave as a bifurcation parameter is examined. The stability of a set of related columnar vortex flows, constructed by taking the two dimensional flow at a single axial location and extending it throughout the domain without variation, is investigated. The method used will be to expand the perturbation velocity in a series of divergence free vectors which ensures that the continuity equation for the incompressible fluid is satisfied exactly by the computed velocity field. Projections of the stability equation onto the space of inviscid vector fields eliminated the pressure term from the equation and reduces the differential eigen problem to a generalized matrix eigen problem. Results are presented both for the one dimensional, columnar vortex flows and also for the wavy bubble flow.

Mhuiris, N. M. G.↗

NASA NTRS1991

The coupling of radiative transfer to quasi 1-D flows with thermochemical nonequilibrium

Quasi-one-dimensional nonequilibrium nozzle flows with coupled radiative transfer are considered. The strongly coupled formulation of radiation and flowfield leads to a governing set of integro-differential equations. A fully implicit numerical method using the full matrix inversion or block iteration methods is presented to solve these equations. The nonequilibrium gas model consists of two chemical species, molecular and atomic nitrogen. The thermodynamic state of the gas is described by two temperatures, translational-rotational and vibrational, and the thermal radiation is assumed to be governed by the vibrational temperature. In radiative transfer, gases are assumed to be absorbing and emitting, and a detailed spectral dependency of the absorption coefficient is prescribed for a non-gray gas. The numerical solutions of strongly radiating nonequilibrium flows are presented for both gray and non-gray gases.

Gokcen, Tahir↗

DOE OSTI2025

Greedy emulators for nuclear two-body scattering

Applications of reduced basis method emulators are increasing in low-energy nuclear physics because they enable fast and accurate sampling of high-fidelity calculations, enabling robust uncertainty quantification. Here, in this paper, we develop, implement, and test two model-driven emulators based on the (Petrov-)Galerkin projection using the prototypical test case of two-body scattering with the Minnesota potential and a more realistic local chiral potential. The high-fidelity scattering equations are solved with the matrix Numerov method, a reformulation of the popular Numerov recurrence relation for solving special second-order differential equations as a linear system of coupled equations. A novel error estimator based on reduced-space residuals is applied to an active learning approach (a greedy algorithm) to choosing training samples (“snapshots”) for the emulator and contrasted with a proper orthogonal decomposition (POD) approach. Both approaches allow for computationally efficient offline-online decompositions, but the greedy approach requires many fewer snapshot calculations. These developments set the groundwork for emulating scattering observables based on chiral nucleon-nucleon and three-nucleon interactions and optical models, where computational speed-ups are necessary for Bayesian uncertainty quantification. Our emulators and error estimators are widely applicable to linear systems.

Bayesian methods↗

NASA NTRS1996

Fast Multigrid Techniques in Total Variation-Based Image Reconstruction

Existing multigrid techniques are used to effect an efficient method for reconstructing an image from noisy, blurred data. Total Variation minimization yields a nonlinear integro-differential equation which, when discretized using cell-centered finite differences, yields a full matrix equation. A fixed point iteration is applied with the intermediate matrix equations solved via a preconditioned conjugate gradient method which utilizes multi-level quadrature (due to Brandt and Lubrecht) to apply the integral operator and a multigrid scheme (due to Ewing and Shen) to invert the differential operator. With effective preconditioning, the method presented seems to require Omicron(n) operations. Numerical results are given for a two-dimensional example.

Oman, Mary Ellen↗

NASA NTRS1989

An integrating matrix formulation for buckling of rotating beams including the effects of concentrated masses

The integrating matrix technique of computational mechanics is extended to include the effects of concentrated masses. The stability of a flexible rotating beam with discrete masses is analyzed to determine the critical rotational speeds for buckling in the inplane and out-of-plane directions. In this problem, the beam is subjected to compressive centrifugal forces arising from steady rotation about an axis which does not pass through the clamped end of the beam. To determine the eigenvalues from which stability is assessed, the differential equations of motion are solved numerically by combining the extended integrating matrix method with an eigenanalysis. Stability boundaries for a discrete mass representation of a uniform beam are shown to asymptotically approach the stability boundaries for the corresponding continuous mass beam as the number of concentrated masses is increased. An error in the literature is also noted for the discrete mass problem concerning the behavior of the critical rotational speed for inplane buckling as the radius of rotation of the clamped end of the beam is reduced.

Lakin, William D.↗

NASA NTRS1972

Precision state and filter weighting matrix extrapolation, revision 3

The matrix extrapolation routine presented is a coded algorithm for the numerical solution of modified forms of the basic differential equations which are satisfied by the geocentric state vector of a spacecraft's center of mass and by the filter-weighting matrix. The equations are described in detail and the various input and output variables designated. The capability is thus established to extrapolate any spacecraft's geocentric state vector, either backwards or forwards in time, through a force field consisting of the earth's primary central-force gravitational attraction, and a superimposed perturbing acceleration.

Robertson, W. M.↗

NASA NTRS1973

Mixed finite-difference scheme for free vibration analysis of noncircular cylinders

A mixed finite-difference scheme is presented for the free-vibration analysis of simply supported closed noncircular cylindrical shells. The problem is formulated in terms of eight first-order differential equations in the circumferential coordinate which possess a symmetric coefficient matrix and are free of the derivatives of the elastic and geometric characteristics of the shell. In the finite-difference discretization, two interlacing grids are used for the different fundamental unknowns in such a way as to avoid averaging in the difference-quotient expressions used for the first derivative. The resulting finite-difference equations are symmetric. The inverse-power method is used for obtaining the eigenvalues and eigenvectors.

Noor, A. K.↗