A method for the numerical integration of coupled first-order differential equations with greatly different time constants.
Numerical integration of coupled first order ODE OF greatly differing time constant
DIFFERENTIAL EQUATION↗
SEARCH · Search NASA
Results for “Coupled first order differential equations”
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.
Solving the Cauchy-Riemann equations on parallel computers
Discussed is the implementation of a single algorithm on three parallel-vector computers. The algorithm is a relaxation scheme for the solution of the Cauchy-Riemann equations; a set of coupled first order partial differential equations. The computers were chosen so as to encompass a variety of architectures. They are: the MPP, and SIMD machine with 16K bit serial processors; FLEX/32, an MIMD machine with 20 processors; and CRAY/2, an MIMD machine with four vector processors. The machine architectures are briefly described. The implementation of the algorithm is discussed in relation to these architectures and measures of the performance on each machine are given. Simple performance models are used to describe the performance. These models highlight the bottlenecks and limiting factors for this algorithm on these architectures. Conclusions are presented.
Fatoohi, Raad A.↗
Singular perturbation equations for flexible satellites
Force equations of motion of the individual flexible elements of a satellite were obtained in a previous paper. Moment equations of motion of the composite bodies of a flexible satellite are to be developed using two sets of equations which form the basic system for any dynamic model of flexible satellites. This basic system consists of a set of N-coupled, nonlinear, ordinary, or partial differential equations, for a flexible satellite with n generalized, structural position coordinates. For single composite body satellites, N is equal to (n + 3); for dual-spin systems, N is equal to (n + 9). These equations involve time derivatives up to the second order. The study shows a method of avoiding this linearization by reducing the N equations to 3 or 9 nonlinear, coupled, first order, ordinary, differential equations involving only the angular velocities of the composite bodies. The solutions for these angular velocities lead to linear equations in the n generalized structural position coordinates, which can be solved by known methods.
Huang, T. C.↗
Nonlinear behavior of acoustic waves in combustion chambers
The nonlinear growth and limiting amplitude of acoustic waves in a combustion chamber are considered. A formal framework is provided within which practical problems can be treated with a minimum of effort and expense. The general conservation equations were expanded in two small parameters, one characterizing the mean flow field and one measuring the amplitude of oscillations, and then combined to yield a nonlinear inhomogeneous wave equation. The unsteady pressure and velocity fields were expressed as syntheses of the normal modes of the chamber, but with unknown time-varying amplitudes. This procedure yielded a representation of a general unsteady field as a system of coupled nonlinear oscillators. The system of nonlinear equations was treated by the method of averaging to produce a set of coupled nonlinear first order differential equations for the amplitudes and phases of the modes. The analysis is applicable to any combustion chamber. The most interesting applications are probably to solid rockets, liquid rockets, or thrust augmentors on jet engines.
Culick, F. E. C.↗
Nonlinear behavior of acoustic waves in combustion chambers. I, II
The general problem of the nonlinear growth and limiting amplitude of acoustic waves in a combustion chamber is treated in three parts: (1) the general conservation equations are expanded in two small parameters, and then combined to yield a nonlinear inhomogeneous wave equation, (2) the unsteady pressure and velocity fields are expressed as a synthesis of the normal modes of the chamber, but with unknown time-varying amplitudes, and (3) the system of nonlinear equations is treated by the method of averaging to produce a set of coupled nonlinear first order differential equations for the amplitudes and phases of the modes. This approximate analysis is applied to the investigation of the unstable motions in a solid propellant rocket engine and in a T burner.
Culick, F. E. C.↗
Comparison of Computed and Measured Performance of a Pulsed Inductive Thruster Operating on Argon Propellant
Pulsed inductive plasma accelerators are electrodeless space propulsion devices where a capacitor is charged to an initial voltage and then discharged through a coil as a high-current pulse that inductively couples energy into the propellant. The field produced by this pulse ionizes the propellant, producing a plasma near the face of the coil. Once a plasma is formed if can be accelerated and expelled at a high exhaust velocity by the Lorentz force arising from the interaction of an induced plasma current and the magnetic field. A recent review of the developmental history of planar-geometry pulsed inductive thrusters, where the coil take the shape of a flat spiral, can be found in Ref. [1]. Two concepts that have employed this geometry are the Pulsed Inductive Thruster (PIT)[2, 3] and the Faraday Accelerator with Radio-frequency Assisted Discharge (FARAD)[4]. There exists a 1-D pulsed inductive acceleration model that employs a set of circuit equations coupled to a one-dimensional momentum equation. The model was originally developed and used by Lovberg and Dailey[2, 3] and has since been nondimensionalized and used by Polzin et al.[5, 6] to define a set of scaling parameters and gain general insight into their effect on thruster performance. The circuit presented in Fig. 1 provides a description of the electrical coupling between the current flowing in the thruster I1 and the plasma current I2. Recently, the model was upgraded to include an equation governing the deposition of energy into various modes present in a pulsed inductive thruster system (acceleration, magnetic flux generation, resistive heating, etc.)[7]. An MHD description of the plasma energy density evolution was tailored to the thruster geometry by assuming only one-dimensional motion and averaging the plasma properties over the spatial dimensions of the current sheet to obtain an equation for the time-evolution of the total energy. The equation set governing the dynamics of the coupled electrodynamic-current sheet system is composed of first-order, coupled ordinary differential equations that can be easily solved numerically without having to resort to much more complex 2-D finite element plasma simulations.
Polzin, Kurt A.↗
Multiple scattering from finite inhomogeneous media
Utilizing the characteristic information concerning the apparent phase constant difference between the electric and magnetic fields propagating in an inhomogeneous medium, a theoretical analysis of the multiple scattering of electromagnetic waves in finite inhomogeneous media is presented. The solution is obtained by first approximating the coefficients of a pair of exact coupled first-order differential equations and then solving the equations by first-order iteration. The present first-order approximate solution with multiple scattering considerations is shown to be more accurate than the WKB solution. Methods to improve the accuracy of the first-order solution further are discussed. Application of the solution to slowly varying finite media with periodic properties demonstrates the validity of the solution. The same approach can be extended to frequencies in the optical region by retaining additional terms in the coefficients of the coupled differential equations.
Tang, C. C. H.↗
On the use of a finite difference method for solving anisotropic scattering problems
A new method of solving the radiative transfer equation is developed in which the scattering and absorption coefficients may have arbitrary variations with depth, and in which both internal (thermal) emission and incident radiation are allowed. Specular and diffuse reflection at both boundaries also is taken into account. The method begins by forming a paired set of coupled first-order differential equations for the symmetric and antisymmetric parts of the radiation field after writing the scattering integral as a numerical quadrature. These differential equations are broken into finite difference form, in which the symmetric and antisymmetric parts of the radiation field are found on alternate grid points. Numerical results for a number of test problems are shown, demonstrating that the method is very fast, that it returns specific intensities and fluxes that are accurate to at least a percent, and that it can be applied to optically thick problems.
Barkstrom, B. R.↗
Analysis of the single-blow transient testing technique for perforated plate heat exchangers
An analysis of the single-blow transient testing technique for heat exchanger surfaces is made for perforated plate and similar discontinuous surfaces. The model assumes that there is no temperature variation across each plate, and allows for axial conduction in spacers: the resulting axial fluid temperature profile is discontinuous, rather than the usual continuous profile. Numerical solutions are obtained to the resulting set of coupled first-order differential equations for a step change in inlet fluid temperature. Results are presented in tabular form, which allow the heat transfer coefficients to be calculated from test data using the maximum slope technique.
Rodriguez, Jose I.↗
Singular perturbation equations for flexible satellites
The dynamic model of a flexible satellite with n generalized structural position coordinates requires the solution of a set of N coupled nonlinear ordinary or partial differential equations. For single composite body satellites, N is equal to (n + 3). For dual-spin systems, N is equal to (n + 9). These equations usually involve time derivatives up to the second order. For large values of n, linearization of the system has so far been the only practicable way of solution. The present study shows a method of avoiding this linearization by reducing the N equations to three or nine nonlinear, coupled, first-order ordinary differential equations involving only the angular velocities of the composite bodies. The solutions for these angular velocities lead to linear equations in the n generalized structural position coordinates, which can then be solved by known methods.
Huang, T. C.↗
A theory of post-stall transients in multistage axial compression systems
A theory is presented for post stall transients in multistage axial compressors. The theory leads to a set of coupled first-order ordinary differential equations capable of describing the growth and possible decay of a rotating-stall cell during a compressor mass-flow transient. These changing flow features are shown to have a significant effect on the instantaneous compressor pumping characteristic during unsteady operation, and henace on the overall system behavior. It is also found from the theory that the ultimate mode of system response, stable rotating stall or surge, depends not only on the B parameter but also on other parameters, such as the compressor length-to-diameter ratio. Small values of this latter quantity tend to favor the occurrence of surge, as do large values of B. A limited parametric study is carried out to show the impact of the different system features on transient behavior. Based on analytical and numerical results, several specific topics are suggested for future research on post-stall transients.
Moore, F. K.↗
A minimum principle for smooth first-order distributed systems
Consideration of the problem of characterizing optimal controls for a class of distributed parameter systems. The system dynamics are characterized mathematically by a finite number of coupled partial differential equations involving first-order time and space derivatives of the state variables. Boundary conditions on the state are in the form of a finite number of algebraic relations between the state and boundary control variables. A few directions are outlined in which the results obtained may be extended further.
Johnson, T. L.↗
Optimal control of first order distributed systems
The problem of characterizing optimal controls for a class of distributed-parameter systems is considered. The system dynamics are characterized mathematically by a finite number of coupled partial differential equations involving first-order time and space derivatives of the state variables, which are constrained at the boundary by a finite number of algebraic relations. Multiple control inputs, extending over the entire spatial region occupied by the system ("distributed controls') are to be designed so that the response of the system is optimal. A major example involving boundary control of an unstable low-density plasma is developed from physical laws.
Johnson, T. L.↗
A finite difference method of solving anisotropic scattering problems
A new method of solving radiative transfer problems is described including a comparison of its speed with that of the doubling method, and a discussion of its accuracy and suitability for computations involving variable optical properties. The method uses a discretization in angle to produce a coupled set of first-order differential equations which are integrated between discrete depth points to produce a set of recursion relations for symmetric and anti-symmetric angular sums of the radiation field at alternate depth points. The formulation given here includes depth-dependent anisotropic scattering, absorption, and internal sources, and allows arbitrary combinations of specular and non-Lambertian diffuse reflection at either or both boundaries. Numerical tests of the method show that it can return accurate emergent intensities even for large optical depths. The method is also shown to conserve flux to machine accuracy in conservative atmospheres
Barkstrom, B. R.↗
Exponential integration algorithms applied to viscoplasticity
Four, linear, exponential, integration algorithms (two implicit, one explicit, and one predictor/corrector) are applied to a viscoplastic model to assess their capabilities. Viscoplasticity comprises a system of coupled, nonlinear, stiff, first order, ordinary differential equations which are a challenge to integrate by any means. Two of the algorithms (the predictor/corrector and one of the implicits) give outstanding results, even for very large time steps.
Freed, Alan D.↗
Updated Chemical Kinetics and Sensitivity Analysis Code
An updated version of the General Chemical Kinetics and Sensitivity Analysis (LSENS) computer code has become available. A prior version of LSENS was described in "Program Helps to Determine Chemical-Reaction Mechanisms" (LEW-15758), NASA Tech Briefs, Vol. 19, No. 5 (May 1995), page 66. To recapitulate: LSENS solves complex, homogeneous, gas-phase, chemical-kinetics problems (e.g., combustion of fuels) that are represented by sets of many coupled, nonlinear, first-order ordinary differential equations. LSENS has been designed for flexibility, convenience, and computational efficiency. The present version of LSENS incorporates mathematical models for (1) a static system; (2) steady, one-dimensional inviscid flow; (3) reaction behind an incident shock wave, including boundary layer correction; (4) a perfectly stirred reactor; and (5) a perfectly stirred reactor followed by a plug-flow reactor. In addition, LSENS can compute equilibrium properties for the following assigned states: enthalpy and pressure, temperature and pressure, internal energy and volume, and temperature and volume. For static and one-dimensional-flow problems, including those behind an incident shock wave and following a perfectly stirred reactor calculation, LSENS can compute sensitivity coefficients of dependent variables and their derivatives, with respect to the initial values of dependent variables and/or the rate-coefficient parameters of the chemical reactions.
Radhakrishnan, Krishnan↗
Implementation of a partitioned algorithm for simulation of large CSI problems
The implementation of a partitioned numerical algorithm for determining the dynamic response of coupled structure/controller/estimator finite-dimensional systems is reviewed. The partitioned approach leads to a set of coupled first and second-order linear differential equations which are numerically integrated with extrapolation and implicit step methods. The present software implementation, ACSIS, utilizes parallel processing techniques at various levels to optimize performance on a shared-memory concurrent/vector processing system. A general procedure for the design of controller and filter gains is also implemented, which utilizes the vibration characteristics of the structure to be solved. Also presented are: example problems; a user's guide to the software; the procedures and algorithm scripts; a stability analysis for the algorithm; and the source code for the parallel implementation.
Alvin, Kenneth F.↗
Viscous Effects in the Elastodynamics of Thick Beams
A viscoelastic higher-order thick beam finite element formulation is extended to include elastodynamic deformations. The material constitutive law is a special differential form of the Maxwell solid. In the constitutive model, the elastic strains and the conjugate viscous strains are coupled through a system of first- order ordinary differential equations. The total time-dependent stress is the superposition of its elastic and viscous components. The elastodynamic equations of motion are derived from the virtual work principle. Computational examples are carried out for a thick orthotropic cantilevered beam. A quasi-static relaxation problem is employed as a validation test for the elastodynamic algorithm. The elastodynamic code is demonstrated by analyzing the damped vibrations of the beam which is deformed and then released to freely vibrate.
Johnson, A. R.↗