SEARCH · Search NASA

Results for “Complete linearization method”

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.

SourceResearch connectionFrom yearThrough year

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

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Numerical Study of Three-dimensional Spatial Instability of a Supersonic Flat Plate Boundary Layer

The behavior of spatially growing three-dimensional waves in a supersonic boundary layer was studied numerically by solving the complete Navier-Stokes equations. Satisfactory comparison with linear parallel and non-parallel stability theories, and experiment are obtained when a small amplitude inflow disturbance is used. The three-dimensional unsteady Navier-Stokes equations are solved by a finite difference method which is fourth-order and second-order accurate in the convection and viscous terms respectively, and second-order accurate in time. Spanwise periodicity is assumed. The inflow disturbance is composed of eigenfunctions from linear stability theory. By increasing the amplitude of the inflow disturbance, nonlinear effects in the form of a relaxation type oscillation of the time signal of rho(u) are observed.

Maestrello, Lucio↗

NASA NTRS1991

On optimal infinite impulse response edge detection filters

The authors outline the design of an optimal, computationally efficient, infinite impulse response edge detection filter. The optimal filter is computed based on Canny's high signal to noise ratio, good localization criteria, and a criterion on the spurious response of the filter to noise. An expression for the width of the filter, which is appropriate for infinite-length filters, is incorporated directly in the expression for spurious responses. The three criteria are maximized using the variational method and nonlinear constrained optimization. The optimal filter parameters are tabulated for various values of the filter performance criteria. A complete methodology for implementing the optimal filter using approximating recursive digital filtering is presented. The approximating recursive digital filter is separable into two linear filters operating in two orthogonal directions. The implementation is very simple and computationally efficient, has a constant time of execution for different sizes of the operator, and is readily amenable to real-time hardware implementation.

Sarkar, Sudeep↗

NASA NTRS1996

Necessary and Sufficient Conditions for LTI Representations of Adaptive Systems with Sinusoidal Regressors

This paper establishes necessary and sufficient coneditions for an adaptive system with a sinusoidal regressor (ie., a regressor comprised of sinusoidal signals) to admit an exact linear time-invariant (LTI) representation. These conditions are important because a large number of adaptive systems used in practice have sinusoidal regressors, and the stability, convergence and robustness properties of systems having LTI representations can be completely analyzed by well-known methods.

adaptive↗

NASA NTRS1980

Stability of nonuniform rotor blades in hover using a mixed formulation

A mixed formulation for calculating static equilibrium and stability eigenvalues of nonuniform rotor blades in hover is presented. The static equilibrium equations are nonlinear and are solved by an accurate and efficient collocation method. The linearized perturbation equations are solved by a one step, second order integration scheme. The numerical results correlate very well with published results from a nearly identical stability analysis based on a displacement formulation. Slight differences in the results are traced to terms in the equations that relate moments to derivatives of rotations. With the present ordering scheme, in which terms of the order of squares of rotations are neglected with respect to unity, it is not possible to achieve completely equivalent models based on mixed and displacement formulations. The one step methods reveal that a second order Taylor expansion is necessary to achieve good convergence for nonuniform rotating blades. Numerical results for a hypothetical nonuniform blade, including the nonlinear static equilibrium solution, were obtained with no more effort or computer time than that required for a uniform blade.

Stephens, W. B.↗

NASA NTRS1976

An improved method for the prediction of completely three-dimensional aerodynamic load distributions of configurations with leading edge vortex separation

The application of a higher-order subsonic potential flow panel method to the solution of three-dimensional flow about wing and wing-body combinations with leading-edge vortex separation is presented. The governing equations are the linear flow differential equation and nonlinear boundary conditions which require that the flow be parallel to the wing and body surfaces and that the free vortex sheet, springing from the leading and trailing edges, be aligned with the local flow and support no pressure jump. The vortex core is modeled as a simple line vortex which receives vorticity from the free sheet through a connecting sheet. The Kutta condition is imposed on all appropriate edges of the wing. This set of nonlinear equations is solved by an iterative procedure. The Goethert rule accounts for compressibility. The method has been programmed for the CDC 6600. Delta wings, gothic wings, arrow wings, cambered wings, and wing with body have been analyzed. Initial studies involving variations of panel density, vortex sheet sizing, Jacobian update, and initial geometry demonstrate that the present method generally exhibits good convergence characteristics.

Rubbert, p. E.↗

NASA NTRS2004

Progressive Aerodynamic Model Identification From Dynamic Water Tunnel Test of the F-16XL Aircraft

Development of a general aerodynamic model that is adequate for predicting the forces and moments in the nonlinear and unsteady portions of the flight envelope has not been accomplished to a satisfactory degree. Predicting aerodynamic response during arbitrary motion of an aircraft over the complete flight envelope requires further development of the mathematical model and the associated methods for ground-based testing in order to allow identification of the model. In this study, a general nonlinear unsteady aerodynamic model is presented, followed by a summary of a linear modeling methodology that includes test and identification methods, and then a progressive series of steps suggesting a roadmap to develop a general nonlinear methodology that defines modeling, testing, and identification methods. Initial steps of the general methodology were applied to static and oscillatory test data to identify rolling-moment coefficient. Static measurements uncovered complicated dependencies of the aerodynamic coefficient on angle of attack and sideslip in the stall region making it difficult to find a simple analytical expression for the measurement data. In order to assess the effect of sideslip on the damping and unsteady terms, oscillatory tests in roll were conducted at different values of an initial offset in sideslip. Candidate runs for analyses were selected where higher order harmonics were required for the model and where in-phase and out-of-phase components varied with frequency. From these results it was found that only data in the angle-of-attack range of 35 degrees to 37.5 degrees met these requirements. From the limited results it was observed that the identified models fit the data well and both the damping-in-roll and the unsteady term gain are decreasing with increasing sideslip and motion amplitude. Limited similarity between parameter values in the nonlinear model and the linear model suggest that identifiability of parameters in both terms may be a problem. However, the proposed methodology can still be used with careful experiment design and carefully selected values of angle of attack, sideslip, amplitude, and frequency of the oscillatory data.

Murphy, Patrick C.↗

NASA NTRS1976

A general unified approach to modelling switching-converter power stages

A method for modeling switching-converter power stages is developed, whose starting point is the unified state-space representation of the switched networks and whose end result is either a complete state-space description or its equivalent small-signal low-frequency linear circuit model. A new canonical circuit model is proposed, whose fixed topology contains all the essential input-output and control properties of any dc-to-dc switching converter, regardless of its detailed configuration, and by which different converters can be characterized in the form of a table conveniently stored in a computer data bank to provide a useful tool for computer aided design and optimization. The new canonical circuit model predicts that, in general, switching action introduces both zeros and poles into the duty ratio to output transfer function in addition to those from the effective filter network.

Middlebrook, R. D.↗

NASA NTRS1978

Improved sonic-box computer program for calculating transonic aerodynamic loads on oscillating wings with thickness

A computer program was developed to account approximately for the effects of finite wing thickness in transonic potential flow over an oscillation wing of finite span. The program is based on the original sonic box computer program for planar wing which was extended to account for the effect of wing thickness. Computational efficiency and accuracy were improved and swept trailing edges were accounted for. Account for the nonuniform flow caused by finite thickness was made by application of the local linearization concept with appropriate coordinate transformation. A brief description of each computer routine and the applications of cubic spline and spline surface data fitting techniques used in the program are given, and the method of input was shown in detail. Sample calculations as well as a complete listing of the computer program listing are presented.

Ruo, S. Y.↗

NASA NTRS1983

Seismic velocity structure of the lunar mantle

The recently completed set of seismic arrival times from the Apollo lunar seismic network are inverted to estimate the average seismic velocities in three sections of the lunar mantle: two for the upper mantle and one for the middle mantle. The method used is a variation of the linearized least squares inversion where the inversion is accomplished in steps. The estimated average velocities in the upper mantle decrease from Vp = 7.74 km/sec and Vs = 4.49 km/sec in the section above 270-km depth to Vp = 7.46 km/sec and Vs = 4.25 km/sec in the section between 270- and 500-km depth, confirming the earlier finding of negative gradients based on seismic amplitude variations. The average velocities in the middle mantle between the depths of 500 km and 1000 km of Vp = 8.26 km/sec and Vs = 4.65 km/sec are significantly higher than those in the upper mantle, contradicting earlier estimates based on more limited data. The higher velocities may suggest initial melting of the moon down to at least 1000-km depth.

Nakamura, Y.↗

NASA NTRS1987

A simple method to model truss-beams as equivalent continua

A new method is presented for obtaining stiffness constants for structures modeled as three-dimensional linear elastic beams. The three-dimensional lattice is comprised of multiple two-dimensional lattices. Matrix manipulation is used to determine the contribution of each two-dimensional substructure to the complete three dimensional lattice.

Barton, O., Jr.↗

NASA NTRS2007

NASA Langley's Approach to the Sandia's Structural Dynamics Challenge Problem

The objective of this challenge is to develop a data-based probabilistic model of uncertainty to predict the behavior of subsystems (payloads) by themselves and while coupled to a primary (target) system. Although this type of analysis is routinely performed and representative of issues faced in real-world system design and integration, there are still several key technical challenges that must be addressed when analyzing uncertain interconnected systems. For example, one key technical challenge is related to the fact that there is limited data on target configurations. Moreover, it is typical to have multiple data sets from experiments conducted at the subsystem level, but often samples sizes are not sufficient to compute high confidence statistics. In this challenge problem additional constraints are placed as ground rules for the participants. One such rule is that mathematical models of the subsystem are limited to linear approximations of the nonlinear physics of the problem at hand. Also, participants are constrained to use these models and the multiple data sets to make predictions about the target system response under completely different input conditions. Our approach involved initially the screening of several different methods. Three of the ones considered are presented herein. The first one is based on the transformation of the modal data to an orthogonal space where the mean and covariance of the data are matched by the model. The other two approaches worked solutions in physical space where the uncertain parameter set is made of masses, stiffnesses and damping coefficients; one matches confidence intervals of low order moments of the statistics via optimization while the second one uses a Kernel density estimation approach. The paper will touch on all the approaches, lessons learned, validation 1 metrics and their comparison, data quantity restriction, and assumptions/limitations of each approach. Keywords: Probabilistic modeling, model validation, uncertainty quantification, kernel density

Horta, Lucas G.↗

NASA NTRS2009

Precise Stabilization of the Optical Frequency of WGMRs

Crystalline whispering gallery mode resonators (CWGMRs) made of crystals with axial symmetry have ordinary and extraordinary families of optical modes. These modes have substantially different thermo-refractive constants. This results in a very sharp dependence of differential detuning of optical frequency on effective temperature. This frequency difference compared with clock gives an error signal for precise compensation of the random fluctuations of optical frequency. Certain crystals, like MgF2, have turnover points where the thermo-refractive effect is completely nullified. An advantage for applications using WGMRs for frequency stabilization is in the possibility of manufacturing resonators out of practically any optically transparent crystal. It is known that there are crystals with negative and zero thermal expansion at some specific temperatures. Doping changes properties of the crystals and it is possible to create an optically transparent crystal with zero thermal expansion at room temperature. With this innovation s stabilization technique, the resultant WGMR will have absolute frequency stability The expansion of the resonator s body can be completely compensated for by nonlinear elements. This results in compensation of linear thermal expansion (see figure). In three-mode, the MgF2 resonator, if tuned at the turnover thermal point, can compensate for all types of random thermal-related frequency drift. Simplified dual-mode method is also available. This creates miniature optical resonators with good short- and long-term stability for passive secondary frequency ethalon and an active resonator for active secondary frequency standard (a narrowband laser with long-term stability).

Savchenkov, Anatoliy↗

NASA NTRS2002

Application of Sliding Mode Methods to the Design of Reconfigurable Flight Control Systems

Observer-based sliding mode control is investigated for application to aircraft reconfigurable flight control. A comprehensive overview of reconfigurable flight control is given, including, a review of the current state-of-the-art within the subdisciplines of fault detection, parameter identification, adaptive control schemes, and dynamic control allocation. Of the adaptive control methods reviewed, sliding mode control (SMC) appears very promising due its property of invariance to matched uncertainty. An overview of sliding mode control is given and its remarkable properties are demonstrated by example. Sliding mode methods, however, are difficult to implement because unmodeled parasitic dynamics cause immediate and severe instability. This presents a challenge for all practical applications with limited bandwidth actuators. One method to deal with parasitic dynamics is the use of an asymptotic observer in the feedback path. Observer-based SMC is investigated, and a method for selecting observer gains is offered. An additional method for shaping the feedback loop using a filter is also developed. It is shown that this SMC prefilter is equivalent to a form of model reference hedging. A complete design procedure is given which takes advantage of the sliding mode boundary layer to recast the SMC as a linear control law. Frequency domain loop shaping is then used to design the sliding manifold. Finally, three aircraft applications are demonstrated. An F-18/HARV is used to demonstrate a SISO pitch rate tracking controller. It is also used to demonstrate a MIMO lateral-directional roll rate tracking controller. The last application is a full linear six degree-of-freedom advanced tailless fighter model. The observer-based SMC is seen to provide excellent tracking with superior robustness to parameter changes and actuator failures.

Wells, Scott R.↗

NASA NTRS2007

Analytic Guidance for the First Entry in a Skip Atmospheric Entry

This paper presents an analytic method to generate a reference drag trajectory for the first entry portion of a skip atmospheric entry. The drag reference, expressed as a polynomial function of the velocity, will meet the conditions necessary to fit the requirements of the complete entry phase. The generic method proposed to generate the drag reference profile is further simplified by thinking of the drag and the velocity as density and cumulative distribution functions respectively. With this notion it will be shown that the reference drag profile can be obtained by solving a linear algebraic system of equations. The resulting drag profile is flown using the feedback linearization method of differential geometric control as guidance law with the error dynamics of a second order homogeneous equation in the form of a damped oscillator. This approach was first proposed as a revisited version of the Space Shuttle Orbiter entry guidance. However, this paper will show that it can be used to fly the first entry in a skip entry trajectory. In doing so, the gains in the error dynamics will be changed at a certain point along the trajectory to improve the tracking performance.

Garcia-Llama, Eduardo↗

NASA NTRS1978

Interpolation and extrapolation of creep rupture data by the minimum commitment method. III - Analysis of multiheats

An outline is presented of approaches for treating multiheats on the basis of the focal point convergence method. The method has been employed in the case of two highly characterized multiheats, including a 304 stainless steel studied in Japan and a low alloy carbon steel studied in England. The method makes use of the same functional form for all materials. Only the constants are varied for each multiheat. Completely computerized procedures are employed for the determination of the constants. Once the basic analysis has been performed, the representation of various members in the same system is achieved by adding linear expressions of log stress, changing only two constants in the equations to represent a selected heat.

Manson, S. S.↗

NASA NTRS1993

Conjugate gradient based projection - A new explicit methodology for frictional contact

With special attention towards the applicability to parallel computation or vectorization, a new and effective explicit approach for linear complementary formulations involving a conjugate gradient based projection methodology is proposed in this study for contact problems with Coulomb friction. The overall objectives are focussed towards providing an explicit methodology of computation for the complete contact problem with friction. In this regard, the primary idea for solving the linear complementary formulations stems from an established search direction which is projected to a feasible region determined by the non-negative constraint condition; this direction is then applied to the Fletcher-Reeves conjugate gradient method resulting in a powerful explicit methodology which possesses high accuracy, excellent convergence characteristics, fast computational speed and is relatively simple to implement for contact problems involving Coulomb friction.

Tamma, Kumar K.↗

NASA NTRS1987

High resolution upwind schemes for the three-dimensional incompressible Navier-Stokes equations

Based on flux-difference splitting, implicit high resolution schemes are constructed for efficient computations of steady-state solutions to the three-dimensional, incompressible Navier-Stokes equations in curvilinear coordinates. These schemes use first-order accurate Euler backward-time differencing and second-order central differencing for the viscous shear fluxes. Up to third-order accurate upwind differencing is achieved through a reconstruction of the solution from its cell averages. The reconstruction is accomplished by linear interpolation, where the node stencils are selected such that in regions of smooth solution the flow is highly resolved while spurious oscillations in regions of rapid changes in gradient are still suppressed. Fairly rapid convergence to steady-state solutions is attained with a completely vectorizable hybrid time-marching method. Flows around a sharp-edged delta wing are computed with the maximum accuracy of the upwind-differencing restricted to first-, second-, and third-order, to illustrate the effect of accuracy on the global and on the local vortical flow fields. The results are validated with experimental data.

Hartwich, PETER-M.↗

NASA NTRS1996

Simulation of Three-Dimensional Symmetric and Asymmetric Instabilities in Attachment-Line Boundary Layers

On a swept wing, contamination along the leading edge, Tollmien-Schlichting waves, stationary or traveling crossflow vortices, and/or Taylor-Gortler vortices can cause the catastrophic breakdown of laminar to turbulent flow, which leads to increased skin-friction drag for the aircraft. The discussion in this Note will be limited to disturbances which evolve along the attachment line (leading edge of swept wing). If the Reynolds number of the attachment-line boundary layer is greater than some critical value, then the complete wing is inevitably engulfed in turbulent flow. Essentially, there are two critical Reynolds number points that must be considered. The first is for small-amplitude disturbances, and the second is for bypass transition. The present study will use direct numerical simulations to validate a linear 2D-eigenvalue prediction method based on parabolized stability equations by Lin and Malik. This method is considered because it suggests that a number of symmetric and asymmetric modes exist and are stable or unstable on the attachment line depending on the Reynolds number. If validated, the approach would predict a number of modes which are linearly damped in the Reynolds number regime 100 to 245; however, these modes may grow nonlinearly and provide an explanation to this region.

Joslin, Ronald D.↗