Search NASA⌕ Search

SEARCH · Search NASA

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

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

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

Quantum mechanical streamlines. I - Square potential barrier

Exact numerical calculations are made for scattering of quantum mechanical particles hitting a square two-dimensional potential barrier (an exact analog of the Goos-Haenchen optical experiments). Quantum mechanical streamlines are plotted and found to be smooth and continuous, to have continuous first derivatives even through the classical forbidden region, and to form quantized vortices around each of the nodal points. A comparison is made between the present numerical calculations and the stationary wave approximation, and good agreement is found between both the Goos-Haenchen shifts and the reflection coefficients. The time-independent Schroedinger equation for real wavefunctions is reduced to solving a nonlinear first-order partial differential equation, leading to a generalization of the Prager-Hirschfelder perturbation scheme. Implications of the hydrodynamical formulation of quantum mechanics are discussed, and cases are cited where quantum and classical mechanical motions are identical.

Hirschfelder, J. O.↗
NASA NTRS1980

An alternative approach to the numerical simulation of steady inviscid flow

A numerical procedure for the efficient simulation of steady inviscid flow is described and its utility is demonstrated. The method is uniformly valid for application in the subsonic, transonic and supersonic flow regimes. It does not rely on the introduction of additional assumptions beyond those necessary to obtain the Euler equations from the Navier-Stokes equations, nor does it make use of a time-asymptotic solution of the unsteady equations of motion. Application of the herein-defined surrogate equation technique allows the formulation of stable, fully-conservative, type-dependent finite difference equations for use in obtaining numerical solutions to systems of first-order partial differential equations, such as the steady-state Euler equations or their various approximations. Computational results are presented for the full Euler equations used to simulate rotational subsonic flow and for the transonic small disturbance equations. For the latter case, a computational efficiency greater than that obtained by means of the standard perturbation potential approach is indicated.

Johnson, G. M.↗
NASA NTRS1989

Nonlinear strain-displacement relations and flexible multibody dynamics

Dynamics of chains of flexible bodies undergoing large rigid body motions, but small elastic deflections are considered. The role of nonlinear strain-displacement relations in the development of the motion equations correct to first order in elastic deflections is investigated. The general form of these equations linearized only in the small elastic deflections is presented, and the relative significance of various nonlinear terms is studied both analytically and through the use of the numerical simulations. Numerical simulations are performed for a two link chain constrained to move in the plane, subject to hinge torques. Each link is modeled as a thin beam. Slew maneuver simulation results are compared for models with and without properly modeled kinematics of deformation. The goal of this case study is to quantify the importance of the terms in the equations of motion which arise from the inclusion of nonlinear strain-displacement relations. It is concluded that unless the consistently linearized equations in elastic deflections and speeds are available and necessary, the inconsistently (prematurely) linearized equations should be replaced in all cases by ruthlessly linearized equations: equations in which all nonlinear terms involving the elastic deflections and speeds are ignored.

Padilla, Carlos E.↗
NASA NTRS1974

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.↗
NASA NTRS1972

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.↗
NASA NTRS1984

Second-order-accurate spatial differencing for the transonic small-disturbance equation

Current methods for calculating transonic flows with the small-disturbance potential equation are typically only first-order accurate in the supersonic regions of the flow. However, calculations using the full-potential equation show significant improvements in accuracy when second-order methods are used instead of first-order methods. In this paper, algorithms for the small-disturbance equations are presented, for both steady and unsteady flows, with spatial differencing that is second-order accurate in both the subsonic and supersonic regions of the flow. These algorithms are stable, simple extensions of implicit monotone algorithms; the former are only first-order accurate spatially in the supersonic regions. Several calculations in one and two dimensions have been made, and the first- and second-order calculations are compared. In one dimension, the calculations include the four types of shock-wave motion. In two dimensions, the comparisons include steady flow over a Korn airfoil and unsteady flow over a pitching airfoil. All the comparisons in two dimensions show that the new methods are just as stable numerically as the old methods, but improve the accuracy of the solutions. The algorithm improvements can be implemented in present computer codes by making minor coding modifications.

Goorjian, P. M.↗
NASA NTRS1990

Stabilization and control of distributed systems with time-dependent spatial domains

This paper considers the problem of the stabilization and control of distributed systems with time-dependent spatial domains. The evolution of the spatial domains with time is described by a finite-dimensional system of ordinary differential equations, while the distributed systems are described by first-order or second-order linear evolution equations defined on appropriate Hilbert spaces. First, results pertaining to the existence and uniqueness of solutions of the system equations are presented. Then, various optimal control and stabilization problems are considered. The paper concludes with some examples which illustrate the application of the main results.

Wang, P. K. C.↗
NASA NTRS2002

High-Order Semi-Discrete Central-Upwind Schemes for Multi-Dimensional Hamilton-Jacobi Equations

We present the first fifth order, semi-discrete central upwind method for approximating solutions of multi-dimensional Hamilton-Jacobi equations. Unlike most of the commonly used high order upwind schemes, our scheme is formulated as a Godunov-type scheme. The scheme is based on the fluxes of Kurganov-Tadmor and Kurganov-Tadmor-Petrova, and is derived for an arbitrary number of space dimensions. A theorem establishing the monotonicity of these fluxes is provided. The spacial discretization is based on a weighted essentially non-oscillatory reconstruction of the derivative. The accuracy and stability properties of our scheme are demonstrated in a variety of examples. A comparison between our method and other fifth-order schemes for Hamilton-Jacobi equations shows that our method exhibits smaller errors without any increase in the complexity of the computations.

Bryson, Steve↗
NASA NTRS1986

Optimal control of aeroassisted plane change maneuver using feedback expansions

A guidance law for an aeroassisted plane change maneuver is developed by an asymptotic expansion technique using a small parameter which essentially represents the ratio of the inertial forces to the atmospheric forces. This guidance law minimizes the energy loss while meeting terminal constraints on the altitude, flight path angle, and heading angle. By neglecting the inertial forces, the resulting optimization problem is integrable and can be determined in closed form. This zeroth-order solution is the first term in an asymptotic series solution of the Hamilton-Jacobi-Bellman equation. The remaining terms are determined from the solution of a first-order, linear partial differential equation whose solution requires only quadrature integration. Our initial results in using this guidance scheme are encouraging.

Mishne, D.↗
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 NTRS1994

Numerical modeling of a cryogenic fluid within a fuel tank

The computational method developed to study the cryogenic fluid characteristics inside a fuel tank in a hypersonic aircraft is presented. The model simulates a rapid draining of the tank by modeling the ullage vapor and the cryogenic liquid with a moving interface. A mathematical transformation was developed and applied to the Navier-Stokes equations to account for the moving interface. The formulation of the numerical method is a transient hybrid explicit-implicit technique where the pressure term in the momentum equations is approximated to first order in time by combining the continuity equation with an ideal equation of state.

Greer, Donald S.↗
NASA NTRS1975

The geometric approach to sets of ordinary differential equations and Hamiltonian dynamics

The calculus of differential forms is used to discuss the local integration theory of a general set of autonomous first order ordinary differential equations. Geometrically, such a set is a vector field V in the space of dependent variables. Integration consists of seeking associated geometric structures invariant along V: scalar fields, forms, vectors, and integrals over subspaces. It is shown that to any field V can be associated a Hamiltonian structure of forms if, when dealing with an odd number of dependent variables, an arbitrary equation of constraint is also added. Families of integral invariants are an immediate consequence. Poisson brackets are isomorphic to Lie products of associated CT-generating vector fields. Hamilton's variational principle follows from the fact that the maximal regular integral manifolds of a closed set of forms must include the characteristics of the set.

Estabrook, F. B.↗
NASA NTRS2011

A First-Order Radiative Transfer Model for Microwave Radiometry of Forest Canopies at L-Band

In this study, a first-order radiative transfer (RT) model is developed to more accurately account for vegetation canopy scattering by modifying the basic Tau-Omega model (the zero-order RT solution). In order to optimally utilize microwave radiometric data in soil moisture (SM) retrievals over vegetated landscapes, a quantitative understanding of the relationship between scattering mechanisms within vegetation canopies and the microwave brightness temperature is desirable. The first-order RT model is used to investigate this relationship and to perform a physical analysis of the scattered and emitted radiation from vegetated terrain. This model is based on an iterative solution (successive orders of scattering) of the RT equations up to the first order. This formulation adds a new scattering term to the . model. The additional term represents emission by particles (vegetation components) in the vegetation layer and emission by the ground that is scattered once by particles in the layer. The model is tested against 1.4-GHz brightness temperature measurements acquired over deciduous trees by a truck-mounted microwave instrument system called ComRAD in 2007. The model predictions are in good agreement with the data, and they give quantitative understanding for the influence of first-order scattering within the canopy on the brightness temperature. The model results show that the scattering term is significant for trees and modifications are necessary to the . model when applied to dense vegetation. Numerical simulations also indicate that the scattering term has a negligible dependence on SM and is mainly a function of the incidence angle and polarization of the microwave observation. Index Terms Emission,microwave radiometry, scattering, soil, vegetation.

Kurum, Mehmet↗
NASA NTRS2010

A First-Order Radiative Transfer Model for Microwave Radiometry of Forest Canopies at L-Band

In this study, a new first-order radiative transfer (RT) model is developed to more accurately account for vegetation canopy scattering by modifying the basic r-co model (the zero-order RT solution). In order to optimally utilize microwave radiometric data in soil moisture (SM) retrievals over moderately to densely vegetated landscapes, a quantitative understanding of the relationship between scattering mechanisms within vegetation canopies and the microwave brightness temperature is desirable. A first-order RT model is used to investigate this relationship and to perform a physical analysis of the scattered and emitted radiation from vegetated terrain. The new model is based on an iterative solution (successive orders of scattering) of the RT equations up to the first order. This formulation adds a new scattering term to the i-w model. The additional term represents emission by particles (vegetation components) in the vegetation layer and emission by the ground that is scattered once by particles in the layer. The new model is tested against 1.4 GHz brightness temperature measurements acquired over deciduous trees by a truck-mounted microwave instrument system called ComRAD in 2007. The model predictions are in good agreement with the data and they give quantitative understanding for the influence of first-order scattering within the canopy on the brightness temperature. The model results show that the scattering term is significant for trees and modifications are necessary to the T-w model when applied to dense vegetation. Numerical simulations also indicate that the scattering term has a negligible dependence on SM and is mainly a function of the angle and polarization of the microwave observation.

Kurum, Mehmet↗
NASA NTRS2007

Application of the Finite Element Method in Atomic and Molecular Physics

The finite element method (FEM) is a numerical algorithm for solving second order differential equations. It has been successfully used to solve many problems in atomic and molecular physics, including bound state and scattering calculations. To illustrate the diversity of the method, we present here details of two applications. First, we calculate the non-adiabatic dipole polarizability of Hi by directly solving the first and second order equations of perturbation theory with FEM. In the second application, we calculate the scattering amplitude for e-H scattering (without partial wave analysis) by reducing the Schrodinger equation to set of integro-differential equations, which are then solved with FEM.

Shertzer, Janine↗
NASA NTRS1985

Quasiperiodically forced damped pendula and Schroedinger equations with quasiperiodic potentials - Implications of their equivalence

Certain first-order nonlinear ordinary differential equations exemplified by strongly damped, quasiperiodically driven pendula and Josephson junctions are isomorphic to Schroedinger equations with quasiperiodic potentials. The implications of this equivalence are discussed. In particular, it is shown that the transition to Anderson localization in the Schroedinger problem corresponds to the occurrence of a novel type of strange attractor in the pendulum problem. This transition should be experimentally observable in the frequency spectrum of the pendulum of Josephson junction.

Bondeson, A.↗