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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

A method for phenomenological analysis of ecological data.

The experimental meaning of the phenomenological differential equations for a competing population is reviewed. It is concluded that it is virtually impossible to construct the differential equations precise enough for studying stability. We consider instead a method of phenomenological analysis which can be applied to a set of population curves. We suggest an ecological index calculated from the population curves, which indicates a group property of the entire system. As a function of time, the index is presumably insensitive to Volterra type fluctuations. A marked increase of the index's value however indicates a marked change of the environmental conditions. It is not easy to deduce the group property from the population curves alone, because a change in population is in general due to the superposition of external disturbances and Volterra fluctuations.

Huang, H.-W.↗
NASA NTRS1982

Spline approximations for nonlinear hereditary control systems

A sline-based approximation scheme is discussed for optimal control problems governed by nonlinear nonautonomous delay differential equations. The approximating framework reduces the original control problem to a sequence of optimization problems governed by ordinary differential equations. Convergence proofs, which appeal directly to dissipative-type estimates for the underlying nonlinear operator, are given and numerical findings are summarized.

Daniel, P. L.↗
NASA NTRS1983

A new approach to optimization for aerodynamic applications

A new approach to solving optimization problems that involve nonlinear partial differential equations is presented. The approach eliminates the need for an inner-outer iterative procedure, solving the partial differential equation only once, thereby reducing the cost of computation to an extent which would allow its use as a practical tool in optimization problems. The approach is tested on a single design parameter problem through the use of a specially developed scheme. The results indicate comparable convergence properties for the present iterative process and the standard iterative scheme. The presented ideas are also applicable to multidesign parameter problems.

Rizk, M. H.↗
NASA NTRS1988

Switch-on shocks

It is argued that parallel propagating Alfven waves steepen to second order in the wave amplitude, rather than first order, as is the case for oblique MHD waves. An analytical solution for the structure of weak switch-on shocks is derived. The Rankine-Hugoniot relations for finite amplitude switch-on shocks are obtained in a compact form. The differential equation for resistive switch-on shocks is solved, and the resistive critical Mach number is obtained. The differential equation for resistive switch-on shocks with finite ion inertia is identical to the one obtained in the pure resistive case. It is argued that efficient electron heat conduction can significantly extend the range of upstream parameters for which switch-on shocks are expected.

Kennel, C. F.↗
NASA NTRS1975

Almost periodic solutions to difference equations

The theory of Massera and Schaeffer relating the existence of unique almost periodic solutions of an inhomogeneous linear equation to an exponential dichotomy for the homogeneous equation was completely extended to discretizations by a strongly stable difference scheme. In addition it is shown that the almost periodic sequence solution will converge to the differential equation solution. The preceding theory was applied to a class of exponentially stable partial differential equations to which one can apply the Hille-Yoshida theorem. It is possible to prove the existence of unique almost periodic solutions of the inhomogeneous equation (which can be approximated by almost periodic sequences) which are the solutions to appropriate discretizations. Two methods of discretizations are discussed: the strongly stable scheme and the Lax-Wendroff scheme.

Bayliss, A.↗
NASA NTRS1977

Some MACSYMA program for solving difference equations

A set of MACSYMA programs are described for finding closed form solutions to linear recurrence relations in equations having either constant or variable coefficients. In the homogenous case, a polymonial equation is obtained and the solution to the recurrence relation can be written as a linear combination of the roots of the polynomial. Exponential generating functions are used to solve variable coefficient relations. Taking successive derivatives and using the recurrence relation, an ordinary differential equation is obtained. Expanding the solution to the differential equation in a Taylor series, shows that the nth term of the series is the solution to the recurrence relation. For second order recurrences, a check is made for those that can be solved in terms of Bessel functions.

Ivie, J.↗
NASA NTRS1980

An analysis of turbulent diffusion flame in axisymmetric jet

The kinetic theory of turbulent flow was employed to study the mixing limited combustion of hydrogen in axisymmetric jets. The integro-differential equations in two spatial and three velocity coordinates describing the combustion were reduced to a set of hyperbolic partial differential equations in the two spatial coordinates by a binodal approximation. The MacCormick's finite difference method was then employed for solution. The flame length was longer than that predicted by the flame-sheet analysis, and was found to be in general agreement with a recent experimental result. Increase of the turbulence energy and scale resulted in an enhancement of the combustion rate and, hence, in a shorter flame length. Details of the numerical method as well as of the physical findings are discussed.

Chung, P. M.↗
NASA NTRS1990

Theory of biaxial graded-index optical fiber

The problem of wave propagation in a biaxial graded-index fiber with circular symmetry is considered. The problem is formulated in terms of four first-order differential equations for the tangential components of the electric and magnetic fields. A general solution method for solving systems of differential equations is presented. This solution method is then used to solve the system of equations for a particular example of a biaxial graded-index fiber. Numerical results for the propagation constant in the fiber are also given.

Kawalko, Stephen F.↗
NASA NTRS1991

A method for the analysis of biaxial graded-index optical fibers

The problem of wave propagation is a biaxial graded-index fiber with circular symmetry is considered. The problem is formulated in terms of four first-order differential equations for the tangential components of the electric and magnetic fields. A general solution method for solving systems of differential equations is presented. This solution method is then used to solve the system of equations for a particular example of a biaxial graded-index fiber. Numerical results for the propagation constant in the fiber are given.

Kawalko, Stephen F.↗
NASA NTRS1993

Inclusion of transverse shear deformation in the exact buckling and vibration analysis of composite plate assemblies

The problem considered is the development of the necessary plate stiffnesses for use in the general purpose program VICONOPT for buckling and vibration of composite plate assemblies. The required stiffnesses include the effects of transverse shear deformation and are for sinusoidal response along the plate length as required in VICONOPT. The method is based on the exact solution of the plate differential equations for a composite laminate having fully populated A, B, and D stiffness matrices which leads to an ordinary differential equation of tenth order.

Anderson, Melvin S.↗
NASA NTRS1976

Lagrangian density for collisional plasma

For the purpose of deriving appropriate Lagrangians for plasma equations that include effects of energy loss, the paper examines the inverse problem of the calculus of variations for systems of first- and second-order quasi-linear partial differential equations. This results in convenient forms of the sufficient conditions under which the given differential equations are Euler-Lagrange equations of a Lagrangian. These conditions are then applied to determine the necessary transformation that converts equations, apparently not already in it, into Euler-Lagrange form. The appropriate Lagrangian for a warm collisional plasma is obtained, and the Lagrangian is derived for a resistive transmission line.

Peng, Y.-K. M.↗
NASA NTRS2009

Stability of Detached Solidification

Bridgman crystal growth can be conducted in the so-called "detached" solidification regime, where the growing crystal is detached from the crucible wall. A small gap between the growing crystal and the crucible wall, of the order of 100 micrometers or less, can be maintained during the process. A meniscus is formed at the bottom of the melt between the crystal and crucible wall. Under proper conditions, growth can proceed without collapsing the meniscus. The meniscus shape plays a key role in stabilizing the process. Thermal and other process parameters can also affect the geometrical steady-state stability conditions of solidification. The dynamic stability theory of the shaped crystal growth process has been developed by Tatarchenko. It consists of finding a simplified autonomous set of differential equations for the radius, height, and possibly other process parameters. The problem then reduces to analyzing a system of first order linear differential equations for stability. Here we apply a modified version of this theory for a particular case of detached solidification. Approximate analytical formulas as well as accurate numerical values for the capillary stability coefficients are presented. They display an unexpected singularity as a function of pressure differential. A novel approach to study the thermal field effects on the crystal shape stability has been proposed. In essence, it rectifies the unphysical assumption of the model that utilizes a perturbation of the crystal radius along the axis as being instantaneous. It consists of introducing time delay effects into the mathematical description and leads, in general, to stability over a broader parameter range. We believe that this novel treatment can be advantageously implemented in stability analyses of other crystal growth techniques such as Czochralski and float zone methods.

Mazuruk, K.↗
NASA NTRS2004

Variable Step Integration Coupled with the Method of Characteristics Solution for Water-Hammer Analysis, A Case Study

One-dimensional water-hammer modeling involves the solution of two coupled non-linear hyperbolic partial differential equations (PDEs). These equations result from applying the principles of conservation of mass and momentum to flow through a pipe, and usually the assumption that the speed at which pressure waves propagate through the pipe is constant. In order to solve these equations for the interested quantities (i.e. pressures and flow rates), they must first be converted to a system of ordinary differential equations (ODEs) by either approximating the spatial derivative terms with numerical techniques or using the Method of Characteristics (MOC). The MOC approach is ideal in that no numerical approximation errors are introduced in converting the original system of PDEs into an equivalent system of ODEs. Unfortunately this resulting system of ODEs is bound by a time step constraint so that when integrating the equations the solution can only be obtained at fixed time intervals. If the fluid system to be modeled also contains dynamic components (i.e. components that are best modeled by a system of ODEs), it may be necessary to take extremely small time steps during certain points of the model simulation in order to achieve stability and/or accuracy in the solution. Coupled together, the fixed time step constraint invoked by the MOC, and the occasional need for extremely small time steps in order to obtain stability and/or accuracy, can greatly increase simulation run times. As one solution to this problem, a method for combining variable step integration (VSI) algorithms with the MOC was developed for modeling water-hammer in systems with highly dynamic components. A case study is presented in which reverse flow through a dual-flapper check valve introduces a water-hammer event. The predicted pressure responses upstream of the check-valve are compared with test data.

Turpin, Jason B.↗
NASA NTRS1994

The Dynamics of Some Iterative Implicit Schemes

The global asymptotic nonlinear behavior of some standard iterative procedures in solving nonlinear systems of algebraic equations arising from four implicit linear multistep methods (LMMs) in discretizing 2 x 2 systems of first-order autonomous nonlinear ordinary differential equations is analyzed using the theory of dynamical systems. With the aid of parallel Connection Machines (CM-2 and CM-5), the associated bifurcation diagrams as a function of the time step, and the complex behavior of the associated 'numerical basins of attraction' of these iterative implicit schemes are revealed and compared. Studies showed that all of the four implicit LMMs exhibit a drastic distortion and segmentation but less shrinkage of the basin of attraction of the true solution than standard explicit methods. The numerical basins of attraction of a noniterative implicit procedure mimic more closely the basins of attraction of the differential equations than the iterative implicit procedures for the four implicit LMMs.

Yee, H. C.↗
NASA NTRS1981

An analytic solution to the classical two-body problem with drag

The two body problem was analyzed with a specific drag model. The model treats drag as a force proportional to the vector velocity and inversely proportional to the distance to the center of attraction. The solution is expressed in terms of known functions and is of a simple and compact form. The time of flight is expressed as a quadrature in the true anomaly. The results are: (1) development of a vector differential equation which allows analysis of an infinite number of gravitational and drag models; and (2) obtaining the solution of a linear differential equation using the inverse method of laplace transforms.

Mittleman, D.↗
NASA NTRS1984

Singular asymptotic expansions in nonlinear rotordynamics

During hot firing ground testing of the Space Shuttle's Main Engine, vibrations of the liquid oxygen pump occur at frequencies which cannot be explained by the linear Jeffcott model of the rotor. The model becomes nonlinear after accounting for deadband, side forces, and rubbing. Two phenomena present in the numerical solutions of the differential equations are unexpected periodic orbits of the rotor and tracking of the nonlinear frequency. A multiple scale asymptotic expansion of the differential equations is used to give an analytic explanation of these characteristics.

Day, W. B.↗
NASA NTRS1985

Singular asymptotic expansions in nonlinear rotordynamics

During hot firing ground testing of the Space shuttle's Main Engine, vibrations of the liquid oxygen pump occur at frequencies which cannot be explained by the linear Jeffcott model of the rotor. The model becomes nonlinear after accounting for deadband, side forces, and rubbing. Two phenomena present in the numerical solutions of the differential equations are unexpected periodic orbits of the rotor and tracking of the nonlinear frequency. A multiple scale asymptotic expansion of the differential equations is used to give an analytic explanation of these characteristics.

Day, W. B.↗