Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “discretization”

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 955 records · Page 53

Discrete time learning control in nonlinear systems

In this paper digital learning control methods are developed primarily for use in single-input, single-output nonlinear dynamic systems. Conditions for convergence of the basic form of learning control based on integral control concepts are given, and shown to be satisfied by a large class of nonlinear problems. It is shown that it is not the gross nonlinearities of the differential equations that matter in the convergence, but rather the much smaller nonlinearities that can manifest themselves during the short time interval of one sample time. New algorithms are developed that eliminate restrictions on the size of the learning gain, and on knowledge of the appropriate sign of the learning gain, for convergence to zero error in tracking a feasible desired output trajectory. It is shown that one of the new algorithms can give guaranteed convergence in the presence of actuator saturation constraints, and indicate when the requested trajectory is beyond the actuator capabilities.

Longman, Richard W.↗

Quantized discrete space oscillators

A quasi-canonical sequence of finite dimensional quantizations was found which has canonical quantization as its limit. In order to demonstrate its practical utility and its numerical convergence, this formalism is applied to the eigenvalue and 'eigenfunction' problem of several harmonic and anharmonic oscillators.

Uzes, C. A.↗

A low-altitude mechanism for mesoscale dynamics, structure, and current filamentation in the discrete aurora

The 2D nonlinear evolution of the ionization-driven adiabatic auroral arc instability is studied. We find: (1) the adiabatic auroral arc instability can fully develop on time scales of tens to hundreds of seconds and on spatial scales of tens to hundreds of kilometers; (2) the evolution of this instability leads to nonlinear 'hook-shaped' conductivity structures: (3) this instability can lead to parallel current filamentation over a wide range of scale sizes; and (4) the k-spectra of the density, electric field, and parallel current develop into inverse power laws in agreement with satellite observations. Comparison with mesoscale auroral phenomenology and current filamentation structures is made.

Keskinen, M. J.↗

An approximately factored incremental strategy for calculating consistent discrete aerodynamic sensitivity derivatives

An incremental strategy is presented for iteratively solving very large systems of linear equations, which are associated with aerodynamic sensitivity derivatives for advanced CFD codes. It is shown that the left-hand side matrix operator and the well-known factorization algorithm used to solve the nonlinear flow equations can also be used to efficiently solve the linear sensitivity equations. Two airfoil problems are considered as an example: subsonic low Reynolds number laminar flow and transonic high Reynolds number turbulent flow.

Korivi, V. M.↗

Effects of discrete-source outgassing on motions of periodic comets and discontinuous orbital anomalies

Isoline maps of the Style II nongravitational parameters A1 and A2 of periodic comets are presented for a baseline model as functions of the parameters of the inertially fixed spin vector and the source's location-dependent thrust angle. The parameters A1 and A2 for a collection of sources are equal to the sum of their values for the individual contributors, thus depending on the source distribution. The model is consistent with even the most extreme among detected perturbations, including instances of rapid temporal changes in A2 and the discontinuous orbital anomalies, which are interpreted as products of a sudden redistribution of the transferred momentum due to episodic activation of short-lived sources. It is argued that the lifetimes of active regions vary considerably from case to case, and forced precession of the nucleus strongly influences the comet's activity pattern.

Sekanina, Zdenek↗

A comparative study of turbulence decay using Navier-Stokes and a discrete particle simulation

A comparative study of the two dimensional temporal decay of an initial turbulent state of flow is presented using a direct Navier-Stokes simulation and a particle method, ranging from the near continuum to more rarefied regimes. Various topics related to matching the initial conditions between the two simulations are considered. The determination of the initial velocity distribution function in the particle method was found to play an important role in the comparison. This distribution was first developed by matching the initial Navier-Stokes state of stress, but was found to be inadequate beyond the near continuum regime. An alternative approach of using the Lees two-sided Maxwellian to match the initial strain-rate is discussed. Results of the comparison of the temporal decay of mean kinetic energy are presented for a range of Knudsen numbers. As expected, good agreement was observed for the near continuum regime, but the differences found for the more rarefied conditions were unexpectedly small.

Goswami, A.↗

Investigation for improving Global Positioning System (GPS) orbits using a discrete sequential estimator and stochastic models of selected physical processes

GEODYNII is a conventional batch least-squares differential corrector computer program with deterministic models of the physical environment. Conventional algorithms were used to process differenced phase and pseudorange data to determine eight-day Global Positioning system (GPS) orbits with several meter accuracy. However, random physical processes drive the errors whose magnitudes prevent improving the GPS orbit accuracy. To improve the orbit accuracy, these random processes should be modeled stochastically. The conventional batch least-squares algorithm cannot accommodate stochastic models, only a stochastic estimation algorithm is suitable, such as a sequential filter/smoother. Also, GEODYNII cannot currently model the correlation among data values. Differenced pseudorange, and especially differenced phase, are precise data types that can be used to improve the GPS orbit precision. To overcome these limitations and improve the accuracy of GPS orbits computed using GEODYNII, we proposed to develop a sequential stochastic filter/smoother processor by using GEODYNII as a type of trajectory preprocessor. Our proposed processor is now completed. It contains a correlated double difference range processing capability, first order Gauss Markov models for the solar radiation pressure scale coefficient and y-bias acceleration, and a random walk model for the tropospheric refraction correction. The development approach was to interface the standard GEODYNII output files (measurement partials and variationals) with software modules containing the stochastic estimator, the stochastic models, and a double differenced phase range processing routine. Thus, no modifications to the original GEODYNII software were required. A schematic of the development is shown. The observational data are edited in the preprocessor and the data are passed to GEODYNII as one of its standard data types. A reference orbit is determined using GEODYNII as a batch least-squares processor and the GEODYNII measurement partial (FTN90) and variational (FTN80, V-matrix) files are generated. These two files along with a control statement file and a satellite identification and mass file are passed to the filter/smoother to estimate time-varying parameter states at each epoch, improved satellite initial elements, and improved estimates of constant parameters.

Goad, Clyde C.↗

Wall interference and boundary simulation in a transonic wind tunnel with a discretely slotted test section

A computational simulation of a transonic wind tunnel test section with longitudinally slotted walls is developed and described herein. The nonlinear slot model includes dynamic pressure effects and a plenum pressure constraint, and each slot is treated individually. The solution is performed using a finite-difference method that solves an extended transonic small disturbance equation. The walls serve as the outer boundary conditions in the relaxation technique, and an interaction procedure is used at the slotted walls. Measured boundary pressures are not required to establish the wall conditions but are currently used to assess the accuracy of the simulation. This method can also calculate a free-air solution as well as solutions that employ the classical homogeneous wall conditions. The simulation is used to examine two commercial transport aircraft models at a supercritical Mach number for zero-lift and cruise conditions. Good agreement between measured and calculated wall pressures is obtained for the model geometries and flow conditions examined herein. Some localized disagreement is noted, which is attributed to improper simulation of viscous effects in the slots.

Al-Saadi, Jassim A.↗

Second order distorted born approximation for backscattering from a layer of discrete random medium

In recent years there has been increasing interest in scattering and depolarization characteristics of the vegetation canopies. Scattering models applied to the microwave remote sensing of vegetation canopies showed that multiple scattering effects can be important in simulating the backscattering coefficients correctly. In particular, in most applications, the cross-polarized backscattering coefficients are often underestimated by single scattering models. Recently, there have been concerted efforts to include the second order terms in the radiative transfer models of vegetation canopies in order to account for multiple scattering within the canopy. The coherent wave theory approach is extended to include multiple scattering effects to predict the coherent and incoherent backscattering contributions from a layer of vegetation canopy. The problem is initially formulated in terms of the exact equation for the correlation function of the field, i.e., the Bethe-Salpeter equation. Using fractional volume as a small parameter, a Foldy type approximation is made to obtain a more manageable correlation equation. This equation is iterated to obtain first and second order solutions. The iteration procedure assumes the variance of the field fluctuations are small compared to the coherent intensity. This assumption proved to be particularly successful in computing backscattering coefficients. First and second order backscattering coefficients are calculated from the iterants of the correlation equation. It is shown that the first order coefficients are the same as the distorted Born results used previously by the authors. These results contained enhancement terms in the direct-reflected contributions. The important contributions to second order backscattering are examined and interpreted in terms of scattering diagrams. Examples of situations in which second order backscattering coefficients are important are given.

Lang, Roger H.↗

The theoretical accuracy of Runge-Kutta time discretizations for the initial boundary value problem: A careful study of the boundary error

The conventional method of imposing time dependent boundary conditions for Runge-Kutta (RK) time advancement reduces the formal accuracy of the space-time method to first order locally, and second order globally, independently of the spatial operator. This counter intuitive result is analyzed in this paper. Two methods of eliminating this problem are proposed for the linear constant coefficient case: (1) impose the exact boundary condition only at the end of the complete RK cycle, (2) impose consistent intermediate boundary conditions derived from the physical boundary condition and its derivatives. The first method, while retaining the RK accuracy in all cases, results in a scheme with much reduced CFL condition, rendering the RK scheme less attractive. The second method retains the same allowable time step as the periodic problem. However it is a general remedy only for the linear case. For non-linear hyperbolic equations the second method is effective only for for RK schemes of third order accuracy or less. Numerical studies are presented to verify the efficacy of each approach.

Carpenter, Mark H.↗

Discrete crack growth analysis methodology for through cracks in pressurized fuselage structures

A methodology for simulating the growth of long through cracks in the skin of pressurized aircraft fuselage structures is described. Crack trajectories are allowed to be arbitrary and are computed as part of the simulation. The interaction between the mechanical loads acting on the superstructure and the local structural response near the crack tips is accounted for by employing a hierarchical modeling strategy. The structural response for each cracked configuration is obtained using a geometrically nonlinear shell finite element analysis procedure. Four stress intensity factors, two for membrane behavior and two for bending using Kirchhoff plate theory, are computed using an extension of the modified crack closure integral method. Crack trajectories are determined by applying the maximum tangential stress criterion. Crack growth results in localized mesh deletion, and the deletion regions are remeshed automatically using a newly developed all-quadrilateral meshing algorithm. The effectiveness of the methodology and its applicability to performing practical analyses of realistic structures is demonstrated by simulating curvilinear crack growth in a fuselage panel that is representative of a typical narrow-body aircraft. The predicted crack trajectory and fatigue life compare well with measurements of these same quantities from a full-scale pressurized panel test.

Potyondy, David O.↗

Time correlations between low and high energy gamma rays from discrete sources

Activities covered the following areas: (1) continuing analysis of the Cygnus Experiment data on the shadowing of cosmic rays by the moon and sun, which led to a direct confirmation of the angular resolution of the CYGNUS EAS array; and (2) development of analysis methods for the daily search overlapping with EGRET targets. To date, no steady emission of ultrahigh energy (UHE) gamma rays from any source has been detected by the Cygnus Experiment, but some evidence for sporadic emission had been found. Upper limits on steady fluxes from 49 sources in the northern hemisphere have been published. In addition, a daily search of 51 possible sources over the interval April 1986 to June 1992 found no evidence for emission. From these source lists, four candidates were selected for comparison with EGRET data.

Ellsworth, R. W.↗