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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

The Vertical Linear Fractional Initialization Problem

This paper presents a solution to the initialization problem for a system of linear fractional-order differential equations. The scalar problem is considered first, and solutions are obtained both generally and for a specific initialization. Next the vector fractional order differential equation is considered. In this case, the solution is obtained in the form of matrix F-functions. Some control implications of the vector case are discussed. The suggested method of problem solution is shown via an example.

Lorenzo, Carl F.↗
NASA NTRS2000

R-Function Relationships for Application in the Fractional Calculus

The F-function, and its generalization the R-function, are of fundamental importance in the fractional calculus. It has been shown that the solution of the fundamental linear fractional differential equation may be expressed in terms of these functions. These functions serve as generalizations of the exponential function in the solution of fractional differential equations. Because of this central role in the fractional calculus, this paper explores various intrarelationships of the R-function, which will be useful in further analysis. Relationships of the R-function to the common exponential function, e(t), and its fractional derivatives are shown. From the relationships developed, some important approximations are observed. Further, the inverse relationships of the exponential function, el, in terms of the R-function are developed. Also, some approximations for the R-function are developed.

Lorenzo, Carl F.↗
NASA NTRS2002

The Statistical Mechanics of Ideal Homogeneous Turbulence

Plasmas, such as those found in the space environment or in plasma confinement devices, are often modeled as electrically conducting fluids. When fluids and plasmas are energetically stirred, regions of highly nonlinear, chaotic behavior known as turbulence arise. Understanding the fundamental nature of turbulence is a long-standing theoretical challenge. The present work describes a statistical theory concerning a certain class of nonlinear, finite dimensional, dynamical models of turbulence. These models arise when the partial differential equations describing incompressible, ideal (i.e., nondissipative) homogeneous fluid and magnetofluid (i.e., plasma) turbulence are Fourier transformed into a very large set of ordinary differential equations. These equations define a divergenceless flow in a high-dimensional phase space, which allows for the existence of a Liouville theorem, guaranteeing a distribution function based on constants of the motion (integral invariants). The novelty of these particular dynamical systems is that there are integral invariants other than the energy, and that some of these invariants behave like pseudoscalars under two of the discrete symmetry transformations of physics, parity, and charge conjugation. In this work the 'rugged invariants' of ideal homogeneous turbulence are shown to be the only significant scalar and pseudoscalar invariants. The discovery that pseudoscalar invariants cause symmetries of the original equations to be dynamically broken and induce a nonergodic structure on the associated phase space is the primary result presented here. Applicability of this result to dissipative turbulence is also discussed.

Shebalin, John V.↗
NASA NTRS2003

Using Bond Graphs for Articulated, Flexible Multi-bodies, Sensors, Actuators, and Controllers with Application to the International Space Station

Conceptually, modeling of flexible, multi-body systems involves a formulation as a set of time-dependent partial differential equations. However, for practical, engineering purposes, this modeling is usually done using the method of Finite Elements, which approximates the set of partial differential equations, thus generalizing the approach to all continuous media. This research investigates the links between the Bond Graph method and the classical methods used to develop system models and advocates the Bond Graph Methodology and current bond graph tools as alternate approaches that will lead to a quick and precise understanding of a flexible multi-body system under automatic control. For long endurance, complex spacecraft, because of articulation and mission evolution the model of the physical system may change frequently. So a method of automatic generation and regeneration of system models that does not lead to implicit equations, as does the Lagrange equation approach, is desirable. The bond graph method has been shown to be amenable to automatic generation of equations with appropriate consideration of causality. Indeed human-interactive software now exists that automatically generates both symbolic and numeric system models and evaluates causality as the user develops the model, e.g. the CAMP-G software package. In this paper the CAMP-G package is used to generate a bond graph model of the International Space Station (ISS) at an early stage in its assembly, Zvezda. The ISS is an ideal example because it is a collection of bodies that are articulated, many of which are highly flexible. Also many reaction jets are used to control translation and attitude, and many electric motors are used to articulate appendages, which consist of photovoltaic arrays and composite assemblies. The Zvezda bond graph model is compared to an existing model, which was generated by the NASA Johnson Space Center during the Verification and Analysis Cycle of Zvezda.

Montgomery, Raymond C.↗
NASA NTRS1999

A Viscoelastic Hybrid Shell Finite Element

An elastic large displacement thick-shell hybrid finite element is modified to allow for the calculation of viscoelastic stresses. Internal strain variables are introduced at he element's stress nodes and are employed to construct a viscous material model. First order ordinary differential equations relate the internal strain variables to the corresponding elastic strains at the stress nodes. The viscous stresses are computed from the internal strain variables using viscous moduli which are a fraction of the elastic moduli. The energy dissipated by the action of the viscous stresses in included in the mixed variational functional. Nonlinear quasi-static viscous equilibrium equations are then obtained. Previously developed Taylor expansions of the equilibrium equations are modified to include the viscous terms. A predictor-corrector time marching solution algorithm is employed to solve the algebraic-differential equations. The viscous shell element is employed to numerically simulate a stair-step loading and unloading of an aircraft tire in contact with a frictionless surface.

Johnson, Arthur↗
DOE OSTI2022

Modeling Heat Transfer and Reaction Kinetics of Biomass in Pyrolysis Feeding Systems

Biomass feedstocks can be used to produce low carbon-intensity chemicals and liquid fuels by displacing fossil-fuel sources. Pyrolysis is the thermal decomposition of biomass in the absence of oxygen and can be used to produce crude biofuels. A typical feeding mechanism for pyrolysis reactors is by lock hopper followed by a horizontal auger feeder. Particle agglomeration and plugging of the auger screw are regularly occurring phenomena that is detrimental to the process since it results in the need for periodic inspection and cleaning. The underlying mechanisms leading to agglomeration and plugging are unclear, but a working hypothesis is that heat from the reactor raises the temperature of the biomass in the feeder to the point at which preliminary decomposition reactions occur and produce "sticky" products. Using principles of heat and mass transfer, we derived one-dimensional differential equations for heat flow from the pyrolysis reactor through the auger feeder. Temperature profiles were solved from the system of one-dimensional differential equations using the shooting method. Subsequently, we used the temperature profiles as a one-way coupling to a reaction model for the biomass in the feeder. Kinetic models were adopted from the literature and used to predict the formation of intermediates, products, char, extractives, and metaplastics as well as phase changes along the auger feeder. Species were identified that correlated with softening of the biomass, causing it to become sticky and plug the feeder. We will discuss the details of this work, our findings, and the direction of future work that includes fully coupled transport and reaction-kinetics simulations.

BIOMASS FUELS↗
NASA NTRS1994

Global Asymptotic Behavior of 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 three models of 2 x 2 systems of first-order autonomous nonlinear ordinary differential equations (ODEs) is analyzed using the theory of dynamical systems. The iterative procedures include simple iteration and full and modified Newton iterations. The results are compared with standard Runge-Kutta explicit methods, a noniterative implicit procedure, and the Newton method of solving the steady part of the ODEs. Studies showed that aside from exhibiting spurious asymptotes, all of the four implicit LMMs can change the type and stability of the steady states of the differential equations (DEs). They also exhibit a drastic distortion but less shrinkage of the basin of attraction of the true solution than standard nonLMM explicit methods. The simple iteration procedure exhibits behavior which is similar to standard nonLMM explicit methods except that spurious steady-state numerical solutions cannot occur. The numerical basins of attraction of the noniterative implicit procedure mimic more closely the basins of attraction of the DEs and are more efficient than the three iterative implicit procedures for the four implicit LMMs. Contrary to popular belief, the initial data using the Newton method of solving the steady part of the DEs may not have to be close to the exact steady state for convergence. These results can be used as an explanation for possible causes and cures of slow convergence and nonconvergence of steady-state numerical solutions when using an implicit LMM time-dependent approach in computational fluid dynamics.

Yee, H. C.↗
NASA NTRS1991

Optimization of transistor design including large signal device/circuit interactions at extremely high frequencies (20-100+GHz)

Transistor design for extremely high frequency applications requires consideration of the interaction between the device and the circuit to which it is connected. Traditional analytical transistor models are to approximate at some of these frequencies and may not account for variations of dopants and semiconductor materials (especially some of the newer materials) within the device. Physically based models of device performance are required. These are based on coupled systems of partial differential equations and typically require 20 minutes of Cray computer time for a single AC operating point. A technique is presented to extract parameters from a few partial differential equation solutions for the device to create a nonlinear equivalent circuit model which runs in approximately 1 second of personal computer time. This nonlinear equivalent circuit model accurately replicates the contact current properties of the device as computed by the partial differential solver on which it is based. Using the nonlinear equivalent circuit model of the device, optimization of systems design can be performed based on device/circuit interactions.

Levy, Ralph↗
NASA NTRS1964

On an Approximation Theorem of Kupka and Smale

Simplified and generalized geometrical proof of Kupka and Smale approximation theorem concerning differential equations defined on closed, compact, infinitely differentiable manifold

THEOREM PROVING↗
NASA NTRS1982

The optimization of convergence for Chebyshev polynomial methods in an unbounded domain

Grosch and Orszag (1977) have performed a numerical analysis of the problem of solving differential equations in a semiinfinite or infinite domain using Chebyshev polynomials. The principal limitation of the conducted study was that it was entirely empirical. Various differential equations were solved in different ways and the numbers were compared. The present investigation has the objective to extend the studies conducted by Grosch and Orszag by deriving asymptotic approximations to the Chebyshev coefficients of simple model functions. This approach makes it possible to conduct more systematic comparisons of different methods, extend the range of comparisons, and, perhaps most important, give simple analytic formulas for choosing the optimum domain size or mapping parameter L for various situations.

Boyd, J. P.↗
NASA NTRS1983

Postbuckling of long orthotropic plates in combined shear and compression

The nonlinear large-deflection partial differential equations of von karman for orthotropic plates loaded in combined shear and compression are converted into a set of first-order nonlinear ordinary differential equations by assuming trigonometric functions in one direction. These equations are solved numerically using a two point boundary problem solver which makes use of Newton's method. Results are obtained which determine the postbuckling behavior of rectangular plates with loading up to about three times the buckling load. Both isotropic and orthotropic composite plates are considered. Results show that orthotropic plates may behave quite differently than isotropic plates and that in-plane boundary conditions are important for plates loaded in shear.

Stein, M.↗
NASA NTRS1986

Predicting the Dynamic Behavior of High-Speed Roller Bearings

Computer program developed for solving nonlinear differential equations of motion of high-speed, lightly loaded, cylindrical roller bearing. Program organized into three major parts: Main program, controls input and output; differential equation solver; and subroutines, used to evaluate vector function. Program has capability of performing either twodimensional or three-dimensional simulation.

Conry, T. F.↗
NASA NTRS1990

Mathematical modeling of the flow field and particle motion in a rotating bioreactor at unit gravity and microgravity

The biotechnology group at NASA Johnson Space Center is developing systems for culturing mammalian cells that stimulate some aspect of microgravity and provide a low shear environment for microgravity-based studies on suspension and anchorage dependent cells. The design of these vessels for culturing cells is based on the need to suspend cells and aggregates of cells and microcarrier beads continually in the culturing medium. The design must also provide sufficient circulation for adequate mass transfer of nutrients to the cells and minimize the total force on the cells. Forces, resulting from sources such as hydrodynamic fluid shear and collisions of cells and walls of the vessels, may damage delicate cells and degrade the formation of three dimensional structures. This study examines one particular design in both unit gravity and microgravity based on two concentric cylinders rotating in the same direction at different speeds to create a Couette flow between them. A numerical simulation for the flow field and the trajectories of particles in the vessel. The flow field for the circulation of the culturing medium is modeled by the Navier-Stokes equations. The forces on a particle are assumed to be drag from the fluid's circulation, buoyancy from the gravitational force and centrifugal force from the rotation of the vessel. The problem requires first solving the system of partial differential equations for the fluid flow by a finite difference method and then solving the system of ordinary differential equations for the trajectories by Gear's stiff method. Results of the study indicate that the trajectories in unit gravity and microgravity are very similar except for small spatial deviations on the fast time scale in unit gravity. The total force per unit cross sectional area on a particle in microgravity, however, is significantly smaller than the corresponding value in unit gravity, which is also smaller than anticipated. Hence, this study indicates that this design for a bioreactor with optimal rates of rotation can provide a good environment for culturing cells in microgravity with adequate circulation and minimal force on the cells.

Boyd, Ernest J.↗
DOE OSTI2025

A physical basis for cosmological correlators from cuts

Significant progress has been made in our understanding of the analytic structure of FRW wavefunction coefficients, facilitated by the development of efficient algorithms to derive the differential equations they satisfy. Moreover, recent findings indicate that the twisted cohomology of the associated hyperplane arrangement defining FRW integrals overestimates the number of integrals required to define differential equations for the wave-function coefficient. We demonstrate that the associated dual cohomology is automatically organized in a way that is ideal for understanding and exploiting the cut/residue structure of FRW integrals. Utilizing this understanding, we develop a systematic approach to organize compatible sequential residues, which dictates the physical subspace of FRW integrals for any n -site, ℓ-loop graph. In particular, the physical subspace of tree-level FRW wavefunction coefficients is populated by differential forms associated to cuts/residues that factorize the integrand of the wavefunction coefficient into only flat space amplitudes. After demonstrating the validity of our construction using intersection theory, we develop simple graphical rules for cut tubings that enumerate the space of physical cuts and, consequently, differential forms without any calculation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1983

Some applications of Lie groups in astrodynamics

Differential equations that arise in astrodynamics are examined from the standpoint of Lie group theory. A summary of the Lie method is given for first degree differential equations. The Kepler problem in Hamiltonian form is treated by this method. Extension of the Lie method to optimal trajectories is outlined.

Jackson, A. A.↗
NASA NTRS1995

Combining Comparison Functions and Finite Element Approximations in CFD

In a variety of potential flow applications, the modal element method has been shown to significantly reduce the numerical grid, employ a more precise grid termination boundary condition, and give theoretical insight to the flow physics. The method employs eigenfunctions to replace the numerical grid over significant portions of the flow field. Generally, a numerical grid is employed around obstacles with complex geometry while eigenfunctions are applied to regions in the flow field where the boundary conditions can easily be satisfied. To handle a wider class of computational fluid dynamics (CFD) problems, the present paper extends the modal element to include function approximations which do not satisfy the governing differential equation. To accomplish this task, a double modal series approximation and weighted residual constraints are developed to force the comparison functions to satisfy the governing differential equation and to interface properly with the finite element solution. As an example, the method is applied to the problem of potential flow in a channel with two-dimensional cylindrical like obstacles. The calculated flow fields are in excellent agreement with exact analytical solutions.

Baumeister, Kenneth J.↗