Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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.

At least 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
NASA NTRS1966

Matrix methods

Partial differential equations and matrix method for obtaining equations of motion in orbit

MOTION EQUATION↗
NASA NTRS1977

General instability analysis

HOLBOAT computer program provides instability analysis of inhomogeneous, anisotropic, right-circular cylinder or segment under combined loading. Program is based on Kirchhoff-Love hypothesis, general anisotropic constitutive equations, and Flugge's differential equations of equilibrium.

Bianca, C. J.↗
NASA NTRS1976

Stability of neutral equations with constant time delays

A method was developed for determining the stability of a scalar neutral equation with constant coefficients and constant time delays. A neutral equation is basically a differential equation in which the highest derivative appears both with and without a time delay. Time delays may appear also in the lower derivatives or the independent variable itself. The method is easily implemented, and an illustrative example is presented.

Barker, L. K.↗
NASA NTRS1976

A three-dimensional, time-dependent model of Mobile Bay

A three-dimensional, time-variant mathematical model for momentum and mass transport in estuaries was developed and its solution implemented on a digital computer. The mathematical model is based on state and conservation equations applied to turbulent flow of a two-component, incompressible fluid having a free surface. Thus, bouyancy effects caused by density differences between the fresh and salt water, inertia from thare river and tidal currents, and differences in hydrostatic head are taken into account. The conservation equations, which are partial differential equations, are solved numerically by an explicit, one-step finite difference scheme and the solutions displayed numerically and graphically. To test the validity of the model, a specific estuary for which scaled model and experimental field data are available, Mobile Bay, was simulated. Comparisons of velocity, salinity and water level data show that the model is valid and a viable means of simulating the hydrodynamics and mass transport in non-idealized estuaries.

Pitts, F. H.↗
NASA NTRS1984

Fracture, failure, and fragmentation

A general approach is required for describing matter of behavior when the failure is likely to involve growth and coalescence of a large number of fractures. Failures of this kind appear frequently in rapid dynamic processes, particularly in the formation of spall fragments. An approach to formulating constitutive relations that accounts for the opening, shear and growth of an ensemble of cracks is discussed. The approach accounts for plastic flow accompanying fragmentation. The resulting constitutive relations were incorporated into a Lagrangian computer program. A theoretical approach to coalescence is described. The simplest formulation uses a linear Liouville equation, with crack growth limited by the mean free path of cracks, assumed constant. This approach allows for an anisotropic distribution of cracks. An alternative approach in which the decrease of the mean free path with increasing crack size is accounted for, but the crack distribution is assumed isotropic is described. A reduction of the governing Liouville equation to an ordinary differential equation of third order is possible, and the result can be used to determine how mean free path decreases with increasing crack size.

Dienes, J. K.↗
NASA NTRS1987

A computational method for the determination of the response of a linear system

A number of issues are discussed relating to the analysis of a linear damped oscillator equation with a forcing term where the right side of the equation, g(t), is known only at fixed, equal time intervals. (A particular class of such equations has application to civil earthquake engineering.) It is noted that the equation is not a differential equation because the condition of uniqueness does not obtain. A general computational method is presented for calculating x(t), based on an unconditionally stable finite difference technique. The work generalizes the results given in the recent paper of Ly (1984).

Mickens, R. E.↗
NASA NTRS1999

Transverse Mode Dynamics of VCSELs Through Space-Time Domain Simulation

Modeling and simulation are important to understand laser operation and to optimize and design device functions. Numerical simulation of VCSEL (Vertical Cavity Surface Emitting Lasers) has been largely based on solving time-independent Helmholtz equation or time dependent coupled mode equations. There are various advantages for choosing these approaches. However, the disadvantages are also apparent. The former cannot handle dynamical mode competition seen in VCSELs, while the latter assumes a given type and number of modes a priori. Furthermore, the microscopic physics of heterstructures and electron-hole plasma is very often represented by a few parameters such as linear gain coefficients and the linewidth enhancement factor. These are over simplification of space and frequency (wavelength) dependent gain and refractive index functions. When the space-time dynamical operation of VCSELs becomes important, these simple approximations become questionable. In this paper, we apply a recently developed model for edge-emitting lasers to a gain guided VCSEL for space-time domain simulation. This model takes into account the actual nonlinear dependence of gain and refractive index on frequency and carrier density within the frame work of the effective Bloch equations. The corresponding partial differential equations are solved directly by finite difference methods. Laser behavior with increasing pumping current is investigated in detail. Special attention is paid to the dynamical competition of the transverse modes.

Goorjian, Peter M.↗
NASA NTRS2008

An Initial Investigation of the Effects of Turbulence Models on the Convergence of the RK/Implicit Scheme

A three-stage Runge-Kutta (RK) scheme with multigrid and an implicit preconditioner has been shown to be an effective solver for the fluid dynamic equations. This scheme has been applied to both the compressible and essentially incompressible Reynolds-averaged Navier-Stokes (RANS) equations using the algebraic turbulence model of Baldwin and Lomax (BL). In this paper we focus on the convergence of the RK/implicit scheme when the effects of turbulence are represented by either the Spalart-Allmaras model or the Wilcox k-! model, which are frequently used models in practical fluid dynamic applications. Convergence behavior of the scheme with these turbulence models and the BL model are directly compared. For this initial investigation we solve the flow equations and the partial differential equations of the turbulence models indirectly coupled. With this approach we examine the convergence behavior of each system. Both point and line symmetric Gauss-Seidel are considered for approximating the inverse of the implicit operator of the flow solver. To solve the turbulence equations we use a diagonally dominant alternating direction implicit (DDADI) scheme. Computational results are presented for three airfoil flow cases and comparisons are made with experimental data. We demonstrate that the two-dimensional RANS equations and transport-type equations for turbulence modeling can be efficiently solved with an indirectly coupled algorithm that uses the RK/implicit scheme for the flow equations.

Swanson, R. C.↗
NASA NTRS2012

Colonel Blotto Games and Lancaster's Equations: A Novel Military Modeling Combination

Military strategists face a difficult task when engaged in a battle against an adversarial force. They have to predict both what tactics their opponent will employ and the outcomes of any resultant conflicts in order to make the best decision about their actions. Game theory has been the dominant technique used by analysts to investigate the possible actions that an enemy will employ. Traditional game theory can be augmented by use of Lanchester equations, a set of differential equations used to determine the outcome of a conflict. This paper demonstrates a novel combination of game theory and Lanchester equations using Colonel Blotto games. Colonel Blotto games, which are one of the oldest applications of game theory to the military domain, look at the allocation of troops and resources when fighting across multiple areas of operation. This paper demonstrates that employing Lanchester equations within a game overcomes some of practical problems faced when applying game theory.

Collins, Andrew J.↗
NASA NTRS1977

Derivation of an integral equation for transonic flows

It is noted that the nonlinear partial differential equation for the perturbation velocity potential and boundary conditions describing steady inviscid compressible transonic flow past a thin two-dimensional airfoil can be transformed into a singular integrodifferential equation and that differentiation of the latter yields an integral equation. Two forms of this integral equation currently exist: one for the singularity that is enclosed in an infinitely long strip of vanishing thickness and the other for the singularity that is enclosed in a vanishing circle. In the present article, a more general integral equation is derived by enclosing the singularity in a vanishing rectangular cavity of arbitrary aspect ratio. The two existing forms of this equation are deduced as special cases distinguished by the respective values for the aspect ratio (infinity for the first form and unity for the second).

Ogana, W.↗
NASA NTRS1977

Symbolic computation of recurrence equations for the Chebyshev series solution of linear ODE's

If a linear ordinary differential equation with polynomial coefficients is converted into integrated form then the formal substitution of a Chebyshev series leads to recurrence equations defining the Chebyshev coefficients of the solution function. An explicit formula is presented for the polynomial coefficients of the integrated form in terms of the polynomial coefficients of the differential form. The symmetries arising from multiplication and integration of Chebyshev polynomials are exploited in deriving a general recurrence equation from which can be derived all of the linear equations defining the Chebyshev coefficients. Procedures for deriving the general recurrence equation are specified in a precise algorithmic notation suitable for translation into any of the languages for symbolic computation. The method is algebraic and it can therefore be applied to differential equations containing indeterminates.

Geddes, K. O.↗
DOE OSTI2021

Implicit-Explicit Multirate Infinitesimal GARK Methods

This work focuses on the development of a new class of high-order accurate methods for multirate time integration of systems of ordinary differential equations. Unlike other recent work in this area, the proposed methods support mixed implicit-explicit (IMEX) treatment of the slow time scale. In addition to allowing this slow time scale flexibility, the proposed methods utilize a so-called infinitesimal formulation for the fast time scale through definition of a sequence of modified “fast" initial-value problems that may be solved using any viable algorithm. We name the proposed class as implicit-explicit multirate infinitesimal generalized-structure additive Runge--Kutta (IMEX-MRI-GARK) methods. In addition to defining these methods, we prove that they may be viewed as specific instances of GARK methods and derive a set of order conditions on the IMEX-MRI-GARK coefficients to guarantee both third and fourth order accuracy for the overall multirate method. Additionally, we provide three specific IMEX-MRI-GARK methods, two of order three and one of order four. We conclude with numerical simulations on two multirate test problems, demonstrating the methods' predicted convergence rates and comparing their efficiency against both legacy IMEX multirate schemes and recent third and fourth order implicit MRI-GARK methods.

97 MATHEMATICS AND COMPUTING↗