Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear equation systems”

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 1,477 records · Page 82

Global organization of tectonic deformation on Venus

The geographic organization of surface deformation on Venus as on Earth is a key to understanding the global tectonic system. To date we have mapped the distribution of three unambiguous tectonic land forms on Venus: (1) linear foldbelts analogous to those at plate margins of the Earth; (2) linear rift zones, analogous to continental rifts on the Earth; and (3) distributed plains deformation in the form of wrinkle ridges and extensional faults and fractures. The linear foldbelts are the dominant structural style in the Northern Hemisphere; ninety percent of the planet's foldbelts lie above the equator. In contrast, compressive deformation in the Southern Hemisphere is dominated by two large, sweeping patterns of wrinkle ridges. The two hemispheres are divided by an equatorial region that is largely covered by rift zones and several large tessera blocks. A tectonic model of generally poleward convergence of the Northern Hemisphere explains the distribution of foldbelts and rift zones. In our model, a northern hemispherical plate (or system of plates) moves poleward and deforms along discrete, predominately longitudinal bands. We recognize four types of foldbelts based on their relationships to other large-scale tectonic features on Venus. There are foldbelts that lie within the low plains, foldbelts associated with coronae, novae and chasmata, foldbelts that lie at the margins of poly-deformed tessera plateaus, and the folded mountain belts around Lakshmi Planum. We see a geometric increase in the area of fold belts when normalized to percent area at a given latitude. This increase is consistent with our model of poleward convergence. Also, the orientations of most foldbelts are either approximately north-south or parallel to lines of latitude in the northern hemisphere. This observation is also consistent with the model in that the longitudinal bands are the result of the decreasing area of the sphere as the plate moves poleward and the latitudinal belts are the direct result of poleward compression. The trends of wrinkle ridges have been mapped over the planet and several large, sweeping patterns evidently reflect long-wavelength topography. Using wrinkle ridges as paleostress indicators, we have developed local and regional stress trajectory maps.

Bilotti, Frank↗

Learning physics-based reduced-order models from data using nonlinear manifolds

Here we present a novel method for learning reduced-order models of dynamical systems using nonlinear manifolds. First, we learn the manifold by identifying nonlinear structure in the data through a general representation learning problem. The proposed approach is driven by embeddings of low-order polynomial form. A projection onto the nonlinear manifold reveals the algebraic structure of the reduced-space system that governs the problem of interest. The matrix operators of the reduced-order model are then inferred from the data using operator inference. Numerical experiments on a number of nonlinear problems demonstrate the generalizability of the methodology and the increase in accuracy that can be obtained over reduced-order modeling methods that employ a linear subspace approximation.

97 MATHEMATICS AND COMPUTING↗

Control Design for a Marine Hydrokinetic Cycloturbine Vehicle

Abstract Marine Hydrokinetic (MHK) cycloturbines generate sustainable power by exploiting tidal currents. By powering the turbines and using pitching foils for control, a vehicle comprised of MHK cycloturbines also has the ability to station keep and maneuver. The vehicle consists of four counter-rotating cycloturbines, with hydrofoils oriented perpendicular to the flow in a paddlewheel configuration. Lift and drag generated from these foils sum together to produce thrust. An experimentally tuned simulation model that solves the six-degrees-of-freedom rigid body equations of motion for the MHK vehicle subject to hydrodynamic, hydrostatic, and propulsive forces is used to aid the design of vehicle controllers. Global feedback controllers are initially designed by applying classical control methods to an approximate linear model of the system dynamics. A higher performing nonlinear controller is designed using the nonlinear dynamic inversion (NDI) method. NDI accounts for the nonlinearities of the MHK system and therefore is suitable for a wide range of operating conditions. The response of the classical and NDI controllers to speed, depth, roll, pitch, and yaw commands are evaluated and compared in simulation. The classical controller outperforms the NDI controller for small amplitude maneuvers, although the degradation with NDI is minor. However, in the nonlinear operating regime the NDI controller outperforms the classical controller and the classical controller exhibits instability.

Automation & Control Systems↗

Power-spectral-density relationship for retarded differential equations

The power spectral density (PSD) relationship between input and output of a set of linear differential-difference equations of the retarded type with real constant coefficients and delays is discussed. The form of the PSD relationship is identical with that applicable to unretarded equations. Since the PSD relationship is useful if and only if the system described by the equations is stable, the stability must be determined before applying the PSD relationship. Since it is sometimes difficult to determine the stability of retarded equations, such equations are often approximated by simpler forms. It is pointed out that some common approximations can lead to erroneous conclusions regarding the stability of a system and, therefore, to the possibility of obtaining PSD results which are not valid.

Barker, L. K.↗

Scaled equation of state parameters for gases in the critical region

In the light of recent theoretical developments, the paper presents an accurate characterization of anomalous thermodynamic behavior of xenon, helium 4, helium 3, carbon dioxide, steam and oxygen in the critical region. This behavior is associated with long range fluctuations in the system and the physical properties depend primarily on a single variable, namely, the correlation length. A description of the thermodynamic behavior of fluids in terms of scaling laws is formulated, and the two successfully used scaled equations of state (NBS equation and Linear Model parametric equation) are compared. Methods for fitting both equations to experimental equation of state data are developed and formulated, and the optimum fit for each of the two scaled equations of the above gases are presented and the results are compared. By extending the experimental data for the above one-component fluids to partially miscible binary liquids, superfluid liquid helium, ferromagnets and solids exhibiting order-disorder transitions, the principle of universality is concluded. Finally by using this principle, the critical regions for nine additional fluids are described.

Sengers, J. M. H. L.↗

Space construction base control system

Several approaches for an attitude control system are studied and developed for a large space construction base that is structurally flexible. Digital simulations were obtained using the following techniques: (1) the multivariable Nyquist array method combined with closed loop pole allocation, (2) the linear quadratic regulator method. Equations for the three-axis simulation using the multilevel control method were generated and are presented. Several alternate control approaches are also described. A technique is demonstrated for obtaining the dynamic structural properties of a vehicle which is constructed of two or more submodules of known dynamic characteristics.

Kaczynski, R. F.↗

Computational attributes of the integral form of the equation of transfer

Difficulties can arise in radiative and neutron transport calculations when a highly anisotropic scattering phase function is present. In the presence of anisotropy, currently used numerical solutions are based on the integro-differential form of the linearized Boltzmann transport equation. This paper, departs from classical thought and presents an alternative numerical approach based on application of the integral form of the transport equation. Use of the integral formalism facilitates the following steps: a reduction in dimensionality of the system prior to discretization, the use of symbolic manipulation to augment the computational procedure, and the direct determination of key physical quantities which are derivable through the various Legendre moments of the intensity. The approach is developed in the context of radiative heat transfer in a plane-parallel geometry, and results are presented and compared with existing benchmark solutions. Encouraging results are presented to illustrate the potential of the integral formalism for computation. The integral formalism appears to possess several computational attributes which are well-suited to radiative and neutron transport calculations.

Frankel, J. I.↗

Recent developments in high order K-exact reconstruction on unstructured meshes

This paper presents recent improvements in high-order K-exact reconstruction on unstructured meshes. The new reconstruction procedures are incorporated into a basic upwind finite-volume scheme suitable for solving scalar advection-diffusion equations as well as the Euler and Navier-Stokes equations. Numerical calculations are performed comparing the present method with lower order accurate reconstruction procedures (piecewise constant and piecewise linear) and various competing technologies such as the fluctuation splitting method of Roe (1987) and Deconinck et al. (1992) and a system-variant of the streamline diffusion Petrov-Galerkin method developed by Hansbo (1991) and Hansbo and Johnson (1991). Five test problems are used in the numerical comparisons: scalar circular advection, transonic and supersonic Euler flow, laminar boundary-layer flow, and general compressible Navier-Stokes flow.

Barth, Timothy J.↗

Spectral response of a bilinear oscillator

An approximate analytical procedure is presented to estimate the response power spectral density of a randomly excited spring/mass/damper system having a bilinear spring. The approximate expression for the response spectrum is developed by representing the nonlinear oscillator as a linear system having a natural frequency that depends on the envelope of the random response. This approximate representation of the system leads to estimates of the response spectrum that agree extremely well with those obtained by direct numerical simulation of the governing equation.

Miles, R. N.↗

A Model Predictive Control to Improve Grid Resilience

The following article details a model predictive control (MPC) to improve grid resilience when faced with variable generation resources. This topic is of significant interest to utility power systems where distributed intermittent energy sources will increase significantly and be relied on for electric grid ancillary services. Previous work on MPCs has focused on narrowly targeted control applications such as improving electric vehicle (EV) charging infrastructure or reducing the cost of integrating Energy Storage Systems (ESSs) into the grid. In contrast, this article develops a comprehensive treatment of the construction of an MPC tailored to electric grids and then applies it integration of intermittent energy resources. To accomplish this, the following article includes a description of a reduced order model (ROM) of an electric power grid based on a circuit model, an optimization formulation that describes the MPC, a collocation method for solving linear time-dependent differential algebraic equations (DAEs) that result from the ROM, and an overall strategy for iteratively refining the behavior of the MPC. Next, the algorithm is validated using two separate numerical experiments. First, the algorithm is compared to an existing MPC code and the results are verified by a numerically precise simulation. It is shown that this algorithm produces a control comparable to existing algorithms and the behavior of the control carefully respects the bounds specified. Second, the MPC is applied to a small nine bus system that contains a mix of turbine-spinning-machine-based and intermittent generation in order to demonstrate the algorithm’s utility for resource planning and control of intermittent resources. This study demonstrates how the MPC can be tuned to change the behavior of the control, which can then assist with the integration of intermittent resources into the grid. The emphasis throughout the paper is to provide systematic treatment of the topic and produce a novel nonlinear control compatible design framework applicable to electric grids and the control of variable resources. This differs from the more targeted application-based focus in most presentations.

microgrid↗

Linear and nonlinear dynamic analysis of redundant load path bearingless rotor systems

The goal of this research is to develop the transfer matrix method to treat nonlinear autonomous boundary value problems with multiple branches. The application is the complete nonlinear aeroelastic analysis of multiple-branched rotor blades. Once the development is complete, it can be incorporated into the existing transfer matrix analyses. There are several difficulties to be overcome in reaching this objective. The conventional transfer matrix method is limited in that it is applicable only to linear branch chain-like structures, but consideration of multiple branch modeling is important for bearingless rotors. Also, hingeless and bearingless rotor blade dynamic characteristics (particularly their aeroelasticity problems) are inherently nonlinear. The nonlinear equations of motion and the multiple-branched boundary value problem are treated together using a direct transfer matrix method. First, the formulation is applied to a nonlinear single-branch blade to validate the nonlinear portion of the formulation. The nonlinear system of equations is iteratively solved using a form of Newton-Raphson iteration scheme developed for differential equations of continuous systems. The formulation is then applied to determine the nonlinear steady state trim and aeroelastic stability of a rotor blade in hover with two branches at the root. A comprehensive computer program is developed and is used to obtain numerical results for the (1) free vibration, (2) nonlinearly deformed steady state, (3) free vibration about the nonlinearly deformed steady state, and (4) aeroelastic stability tasks. The numerical results obtained by the present method agree with results from other methods.

Murthy, V. R.↗

Ballistic limit regression analysis for Space Station Freedom meteoroid and space debris protection system

Relationships defining the ballistic limit of Space Station Freedom's (SSF) dual wall protection systems have been determined. These functions were regressed from empirical data found in Marshall Space Flight Center's (MSFC) Hypervelocity Impact Testing Summary (HITS) for the velocity range between three and seven kilometers per second. A stepwise linear least squares regression was used to determine the coefficients of several expressions that define a ballistic limit surface. Using statistical significance indicators and graphical comparisons to other limit curves, a final set of expressions is recommended for potential use in Probability of No Critical Flaw (PNCF) calculations for Space Station. The three equations listed below represent the mean curves for normal, 45 degree, and 65 degree obliquity ballistic limits, respectively, for a dual wall protection system consisting of a thin 6061-T6 aluminum bumper spaced 4.0 inches from a .125 inches thick 2219-T87 rear wall with multiple layer thermal insulation installed between the two walls. Normal obliquity is d(sub c) = 1.0514 v(exp 0.2983 t(sub 1)(exp 0.5228). Forty-five degree obliquity is d(sub c) = 0.8591 v(exp 0.0428) t(sub 1)(exp 0.2063). Sixty-five degree obliquity is d(sub c) = 0.2824 v(exp 0.1986) t(sub 1)(exp -0.3874). Plots of these curves are provided. A sensitivity study on the effects of using these new equations in the probability of no critical flaw analysis indicated a negligible increase in the performance of the dual wall protection system for SSF over the current baseline. The magnitude of the increase was 0.17 percent over 25 years on the MB-7 configuration run with the Bumper II program code.

Jolly, William H.↗

A numerical solution of the Navier-Stokes equations for supercritical fluid thermodynamic analysis

An explicit numerical solution of the compressible Navier-Stokes equations is applied to the thermodynamic analysis of supercritical oxygen in the Apollo cryogenic storage system. The wave character is retained in the conservation equations which are written in the basic fluid variables for a two-dimensional Cartesian coordinate system. Control-volume cells are employed to simplify imposition of boundary conditions and to ensure strict observance of local and global conservation principles. Non-linear real-gas thermodynamic properties responsible for the pressure collapse phenomonon in supercritical fluids are represented by tabular and empirical functions relating pressure and temperature to density and internal energy. Wall boundary conditions are adjusted at one cell face to emit a prescribed mass flowrate. Scaling principles are invoked to achieve acceptable computer execution times for very low Mach number convection problems. Detailed simulations of thermal stratification and fluid mixing occurring under low acceleration in the Apollo 12 supercritical oxygen tank are presented which model the pressure decay associated with de-stratification induced by an ordinary vehicle maneuver and heater cycle operation.

Heinmiller, P. J.↗

Stability of a steady, large amplitude whistler wave.

Study of the behavior of weak electrostatic waves in a collisionless magnetoplasma supporting a steady large amplitude whistler wave. All waves are assumed to propagate parallel to a uniform background magnetic field B sub zero. In the presence of the whistler wave fields each particle executes an oscillatory motion parallel to B sub zero, in addition to a translation along B sub zero and transverse motions. This oscillation causes the Landau resonance to be replaced by a series of new resonances between particles and the electrostatic modes. A distribution function for the perturbed plasma is constructed by solving the Vlasov equation, linearized in the electrostatic wave amplitudes. A dispersion relation is obtained and solved approximately for the growth/damping rate of the perturbations. Growing electrostatic modes are found to be approximately uncoupled. Trapped particles have a strong influence on the stability of the system.

Palmadesso, P. J.↗

Stability and Relative Stability of Linear Systems with Many Constant Time Delays

A method of determining the stability of linear systems with many constant time delays is developed. This technique, an extension of the tau-decomposition method, is used to examine not only the stability but also the relative stability of retarded systems with many delays and a class of neutral equations with one delay. Analytical equations are derived for partitioning the delay space of a retarded system with two time delays. The stability of the system in each of the regions defined by the partitioning curves in the parameter plane is determined using the extended tau-decomposition method. In addition, relative stability boundaries are defined using the extended tau-decompositon method in association with parameter plane techniques. Several applications of the extended tau-decomposition method are presented and compared with stability results obtained from other analyses. In all cases the results obtained using the method outlined herein coincide with and extend those of previous investigations. The extended tau-decomposition method applied to systems with time delays requires less computational effort and yields more complete stability analyses than previous techniques.

Barker, Larry Keith↗

Comparative study of flare control laws

A digital 3-D automatic control law was developed to achieve an optimal transition of a B-737 aircraft between various initial glid slope conditions and the desired final touchdown condition. A discrete, time-invariant, optimal, closed-loop control law presented for a linear regulator problem, was extended to include a system being acted upon by a constant disturbance. Two forms of control laws were derived to solve this problem. One method utilized the feedback of integral states defined appropriately and augmented with the original system equations. The second method formulated the problem as a control variable constraint, and the control variables were augmented with the original system. The control variable constraint control law yielded a better performance compared to feedback control law for the integral states chosen.

Nadkarni, A. A.↗

General Rotorcraft Aeromechanical Stability Program (GRASP): Theory manual

The general rotorcraft aeromechanical stability program (GRASP) was developed to calculate aeroelastic stability for rotorcraft in hovering flight, vertical flight, and ground contact conditions. GRASP is described in terms of its capabilities and its philosophy of modeling. The equations of motion that govern the physical system are described, as well as the analytical approximations used to derive them. The equations include the kinematical equation, the element equations, and the constraint equations. In addition, the solution procedures used by GRASP are described. GRASP is capable of treating the nonlinear static and linearized dynamic behavior of structures represented by arbitrary collections of rigid-body and beam elements. These elements may be connected in an arbitrary fashion, and are permitted to have large relative motions. The main limitation of this analysis is that periodic coefficient effects are not treated, restricting rotorcraft flight conditions to hover, axial flight, and ground contact. Instead of following the methods employed in other rotorcraft programs. GRASP is designed to be a hybrid of the finite-element method and the multibody methods used in spacecraft analysis. GRASP differs from traditional finite-element programs by allowing multiple levels of substructure in which the substructures can move and/or rotate relative to others with no small-angle approximations. This capability facilitates the modeling of rotorcraft structures, including the rotating/nonrotating interface and the details of the blade/root kinematics for various types. GRASP differs from traditional multibody programs by considering aeroelastic effects, including inflow dynamics (simple unsteady aerodynamics) and nonlinear aerodynamic coefficients.

Hodges, Dewey H.↗

Can a variable alpha induce limit cycle behavior and exponential luminosity decay in transient soft x ray sources?

There has been, recently, a revival of the stability problem of accretion disks. Much of this renewed interest is due to recent observational data on transient soft X-ray novae, which are low-mass X-ray binaries. It is widely believed that nonsteady mass transfer from the secondary onto the compact primary, through an accretion disk, is the reason for the observed spectacular events in the form of often repetitive outbursts, with recurrence times ranging from 1 to 60 yr and duration time on the scale of months. Though not having reached yet a consensus about the nature of the mechanism that regulates the mass transfer, the disk thermal instability model seems to be favored by the fact that the rise in the hard X-ray luminosity is prior to the rise in the soft X-ray luminosity, while the mass transfer instability model seems to be hindered by the fact that the luminosity during quiescence is unable to trigger the thermal instability. However, it should be stressed that, remarkably, the X-ray light curves of these X-ray novae all show overall exponential decays, a feature quite difficult to reproduce in the framework of the viscous disk model, which yields powerlike luminosity decay. Taking into account this observational constraint, we have studied the temporal evolution of perturbations in the accretion rate, under the assumption that alpha is radial and parameter dependent. The chosen dependence is such that the model can reproduce limit cycle behavior (the system is locally unstable but globally stable). However, the kind of dependence we are looking for in alpha does not allow us to use the usual Shakura and Sunyaev procedure in the sense that we no longer can obtain a linearized continuity equation without explicit dependence on the accretion rate. This is so because now we cannot eliminate the accretion rate by using the angular momentum conservation equation.

Meirellesfilho, C.↗