Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical partial 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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Calculation of Scattering Amplitude Without Partial Analysis: Inclusion of Exchange - II

There was a method for calculating the whole scattering amplitude, f(Omega(sub k)), directly. The idea was to calculate the complete wave function Psi numerically, and use it in an integral expression for f, which can be reduced to a 2 dimensional quadrature. The original application was for e-H scattering without exchange. There the Schrodinger reduces a 2-d partial differential equation (pde), which was solved using the finite element method (FEM). Here we extend the method to the exchange approximation. The S.E. can be reduced to a pair of coupled pde's, which are again solved by the FEM. The formal expression for f(Omega(sub k)) consists two integrals, f+/- = f(sub d) +/- f(sub e); f(sub d) is formally the same integral as the no-exchange f. We have also succeeded in reducing f(sub e) to a 2-d integral. Results will be presented at the meeting.

Temkin, Aaron↗
NASA NTRS1985

Unsteady Heat Transfer Due to Time-dependent Free Stream Velocity

The present work was undertaken in order to study the unsteady combined convection from a horizontal circular cylinder to a transverse flow. A coordinate perturbation method is used to transform the governing set of partial differential equations into a system of ordinary differential equations. The free stream time-dependent velocity was assumed to be sinusoidal and the boundary layer response due to both low as well as high frequencies of oscillation will be studied. Currently numerical solutions are being obtained for the distribution of the unsteady Nusselt number and the friction factor.

Gorla, R. S. R.↗
NASA NTRS1987

Implicit finite difference methods on composite grids

Techniques for eliminating time lags in the implicit finite-difference solution of partial differential equations are investigated analytically, with a focus on transient fluid dynamics problems on overlapping multicomponent grids. The fundamental principles of the approach are explained, and the method is shown to be applicable to both rectangular and curvilinear grids. Numerical results for sample problems are compared with exact solutions in graphs, and good agreement is demonstrated.

Mastin, C. Wayne↗
NASA NTRS1989

Essentially nonoscillatory spectral Fourier methods for shock wave calculations

An essentially nonoscillatory spectral Fourier method for the solution of hyperbolic partial differential equations is presented. The method is based on adding a nonsmooth function to the trigonometric polynomials which are the usual basis functions for the Fourier method. The high accuracy away from the shock is enhanced by using filters. Numerical results confirm that essentially no oscillations develop in the solution.

Cai, Wei↗
NASA NTRS2010

Algorithm for Stabilizing a POD-Based Dynamical System

This algorithm provides a new way to improve the accuracy and asymptotic behavior of a low-dimensional system based on the proper orthogonal decomposition (POD). Given a data set representing the evolution of a system of partial differential equations (PDEs), such as the Navier-Stokes equations for incompressible flow, one may obtain a low-dimensional model in the form of ordinary differential equations (ODEs) that should model the dynamics of the flow. Temporal sampling of the direct numerical simulation of the PDEs produces a spatial time series. The POD extracts the temporal and spatial eigenfunctions of this data set. Truncated to retain only the most energetic modes followed by Galerkin projection of these modes onto the PDEs obtains a dynamical system of ordinary differential equations for the time-dependent behavior of the flow. In practice, the steps leading to this system of ODEs entail numerically computing first-order derivatives of the mean data field and the eigenfunctions, and the computation of many inner products. This is far from a perfect process, and often results in the lack of long-term stability of the system and incorrect asymptotic behavior of the model. This algorithm describes a new stabilization method that utilizes the temporal eigenfunctions to derive correction terms for the coefficients of the dynamical system to significantly reduce these errors.

Kalb, Virginia L.↗
NASA NTRS1999

Continuum Modeling of the Dynamics of Externally Injection-locked Coupled Oscillator Arrays

Mutually injection-locked arrays of electronic oscillators provide a novel means of controlling the aperture phase of a phased-array antenna, thus achieving the advantages of spatial power combining while retaining the ability to steer the radiated beam. In a number of design concepts, one or more of the oscillators are injection locked to a signal from an external master-oscillator. The behavior of such a system has been analyzed by numerical solution of a system of nonlinear differential equations which, due to its complexity, yields limited insight into the relationship between the injection signals and the aperture phase. In this paper, we develop a continuum model, which results in a single partial differential equation for the aperture phase as a function of time. Solution of the equation is effected by means of the Laplace transformation and yields detailed information concerning the dynamics of the array under the influence of the external injection signals.

injection locked↗
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 NTRS1976

Diffraction of a shock wave by a compression corner. I - Regular reflection

The unsteady, two-dimensional flowfield resulting from the interaction of a moving planar shock wave with a compression corner is determined using a second-order, discontinuity-fitting, finite-difference approach. The time-dependent Euler equations are transformed to normalize the distance between the body and peripheral shock and to include the existing self-similar property of the flow. The resulting set of partial differential equations in conservation-law form is then solved in a time-dependent fashion using MacCormack's scheme. The vortical singularity, which lies on the body surface, and the single reflected shock are both treated as discontinuities in the numerical procedure. The results of the numerical simulation compare quite favorably with existing experimental interferograms and yield better flowfield resolution than previous first-order, shock-capturing, numerical solutions.

Kutler, P.↗
NASA NTRS1980

Iterative methods based upon residual averaging

Iterative methods for solving boundary value problems for systems of nonlinear partial differential equations are discussed. The methods involve subtracting an average of residuals from one approximation in order to arrive at a subsequent approximation. Two abstract methods in Hilbert space are given and application of these methods to quasilinear systems to give numerical schemes for such problems is demonstrated. Potential theoretic matters related to the iteration schemes are discussed.

Neuberger, J. W.↗
NASA NTRS1991

On the anomaly of velocity-pressure decoupling in collocated mesh solutions

The use of various pressure correction algorithms originally developed for fully staggered meshes can yield a velocity-pressure decoupled solution for collocated meshes. The mechanism that causes velocity-pressure decoupling is identified. It is shown that the use of a partial differential equation for the incremental pressure eliminates such a mechanism and yields a velocity-pressure coupled solution. Example flows considered are a three dimensional lid-driven cavity flow and a laminar flow through a 90 deg bend square duct. Numerical results obtained using the collocated mesh are in good agreement with the measured data and other numerical results.

Kim, Sang-Wook↗
NASA NTRS1982

On current aspects of finite element computational fluid mechanics for turbulent flows

A set of nonlinear partial differential equations suitable for the description of a class of turbulent three-dimensional flow fields in select geometries is identified. On the basis of the concept of enforcing a penalty constraint to ensure accurate accounting of ordering effects, a finite element numerical solution algorithm is established for the equation set and the theoretical aspects of accuracy, convergence and stability are identified and quantized. Hypermatrix constructions are used to formulate the reduction of the computational aspects of the theory to practice. The robustness of the algorithm, and the computer program embodiment, have been verified for pertinent flow configurations.

Baker, A. J.↗
NASA NTRS1989

Modeling and control of flexible space stations (slew maneuvers)

Large orbiting space structures are expected to experience mechanical vibrations arising from several disturbing forces such as those induced by shuttle takeoff or docking and crew movements. The problem is considered of modeling and control of large space structures subject to these and other disturbing forces. The system consists of a (rigid) massive body, which may play the role of experimental modules located at the center of the space station and flexible configurations, consisting of several beams, forming the space structure. A complete dynamic model of the system was developed using Hamilton's principle. This model consists of radial equations describing the translational motion of the central body, rotational equations describing the attitude motions of the body and several beam equations governing the vibration of the flexible members (platform) including appropriate boundary conditions. In summary, the dynamics of the space structure is governed by a complex system of interconnected partial and ordinary differential equations. Using Lyapunov's approach the asymptotic stability of the space structure is investigated. For asymptotic stability of the rest state (nominal trajectory), feedback controls are suggested. In the investigation, stability of the slewing maneuvers is also considered. Several numerical results are presented for illustration of the impact of coupling and the effectiveness of the stabilizing controls. Some insight is provided into the complexity of modeling, analysis and stabilization of actual space structures.

Ahmed, N. U.↗
NASA NTRS1994

Chaotic dynamics of corotating magnetospheric convection

The corotating plasma convection system of the Jovian magnetosphere is analyzed. The macroscopic (mhd) model introduced by Summers and Mu, (1992) that incorporates the effects of microdiffusion is extended by including previously neglected density effects. We reduce the governing partial differential equations to a third-order ordinary differential system by the Galerkin technique of mode truncation. We carry out such a severe truncation partly in the interests of tractability, and leave open the question of the efficacy of adding additional modes. Exhaustive numerical integrations are carried out to calculate the long-term solutions, and we discover that a rich array of plasma motions is possible, dependent on the value of the height-integrated ionospheric Pederson conductivity Sigma. If Sigma is less than a certain critical value Sigma(sub c), then plasma motion can be expected to be chaotic (or periodic), while if Sigma is greater than Sigma(sub c), then steady state convection is expected. In the former case, whether the plasma motion is chaotic or periodic (and, if periodic, the magnitude of the period) can be very sensitive to the value of Sigma. The value of Sigma(sub c), which is a function of a parameter q that occurs in the assumed form of the stationary radial profile (varies as L(exp -q) of the plasma mass per unit magnetic flux, lies well within the accepted range of values of Sigma for Jupiter, i.e. Sigma greater than or equal to 0.1 mho and less than or equal to 10 mho.

Summers, Danny↗
NASA NTRS2018

Implementation and Use of a Convergence Test for Iterative Wind Tunnel Balance Load Predictions

A modified and extended version of a convergence test for wind tunnel strain-gage balance load iterations was implemented. The test uses an upper bound of the Lipschitz constant to assess convergence characteristics of balance load predictions if the Iterative Method is applied. Convergence is expected within the use envelope of the balance whenever this upper bound is less than the threshold of one. It is explained in great detail how the convergence test can be applied to the two load iteration equation types that are currently being used in the aerospace testing community. In addition, the application of the test to balances with bi-directional output characteristics is discussed. It is also shown how numerical differentiation can be used in order to obtain partial derivatives that are needed for the calculation of the Lipschitz constant. Finally, machine calibration data of NASA's MC60E six-component force balance is selected to demonstrate both implementation and use of the convergence test.

calibration analysis↗
NASA NTRS1988

Observability of discretized partial differential equations

It is shown that complete observability of the discrete model used to assimilate data from a linear partial differential equation (PDE) system is necessary and sufficient for asymptotic stability of the data assimilation process. The observability theory for discrete systems is reviewed and applied to obtain simple observability tests for discretized constant-coefficient PDEs. Examples are used to show how numerical dispersion can result in discrete dynamics with multiple eigenvalues, thereby detracting from observability.

Cohn, Stephen E.↗
NASA NTRS1997

Frequency Response Calculations of Input Characteristics of Cavity-Backed Aperture Antennas Using AWE with Hybrid FEM/MoM Technique

Application of Asymptotic Waveform Evaluation (AWE) is presented in conjunction with a hybrid Finite Element Method (FEM)/Method of Moments (MoM) technique to calculate the input characteristics of cavity-backed aperture antennas over a frequency range. The hybrid FEM/MoM technique is used to form an integro-partial-differential equation to compute the electric field distribution of the cavity-backed aperture antenna. The electric field, thus obtained, is expanded in a Taylor series around the frequency of interest. The coefficients of 'Taylor series (called 'moments') are obtained using the frequency derivatives of the integro-partial-differential Equation formed by the hybrid FEM/MoM technique. Using the moments, the electric field in the cavity is obtained over a frequency range. Using the electric field at different frequencies, the input characteristics of the antenna are obtained over a wide frequency band. Numerical results for an open coaxial line, probe fed cavity, and cavity-backed microstrip patch antennas are presented. Good agreement between AWE and the exact solution over the frequency range is observed.

Reddy, C. J.↗
NASA NTRS1998

Application of Model Based Parameter Estimation for Fast Frequency Response Calculations of Input Characteristics of Cavity-Backed Aperture Antennas Using Hybrid FEM/MoM Technique

Model Based Parameter Estimation (MBPE) is presented in conjunction with the hybrid Finite Element Method (FEM)/Method of Moments (MoM) technique for fast computation of the input characteristics of cavity-backed aperture antennas over a frequency range. The hybrid FENI/MoM technique is used to form an integro-partial- differential equation to compute the electric field distribution of a cavity-backed aperture antenna. In MBPE, the electric field is expanded in a rational function of two polynomials. The coefficients of the rational function are obtained using the frequency derivatives of the integro-partial-differential equation formed by the hybrid FEM/ MoM technique. Using the rational function approximation, the electric field is obtained over a frequency range. Using the electric field at different frequencies, the input characteristics of the antenna are obtained over a wide frequency range. Numerical results for an open coaxial line, probe-fed coaxial cavity and cavity-backed microstrip patch antennas are presented. Good agreement between MBPE and the solutions over individual frequencies is observed.

Reddy C. J.↗