Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bilinear control”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

GRID3D-v2: An updated version of the GRID2D/3D computer program for generating grid systems in complex-shaped three-dimensional spatial domains

In order to generate good quality systems for complicated three-dimensional spatial domains, the grid-generation method used must be able to exert rather precise controls over grid-point distributions. Several techniques are presented that enhance control of grid-point distribution for a class of algebraic grid-generation methods known as the two-, four-, and six-boundary methods. These techniques include variable stretching functions from bilinear interpolation, interpolating functions based on tension splines, and normalized K-factors. The techniques developed in this study were incorporated into a new version of GRID3D called GRID3D-v2. The usefulness of GRID3D-v2 was demonstrated by using it to generate a three-dimensional grid system in the coolent passage of a radial turbine blade with serpentine channels and pin fins.

Steinthorsson, E.↗
NASA NTRS1995

Optimal guidance law development for an advanced launch system

The objective of this research effort was to develop a real-time guidance approach for launch vehicles ascent to orbit injection. Various analytical approaches combined with a variety of model order and model complexity reduction have been investigated. Singular perturbation methods were first attempted and found to be unsatisfactory. The second approach based on regular perturbation analysis was subsequently investigated. It also fails because the aerodynamic effects (ignored in the zero order solution) are too large to be treated as perturbations. Therefore, the study demonstrates that perturbation methods alone (both regular and singular perturbations) are inadequate for use in developing a guidance algorithm for the atmospheric flight phase of a launch vehicle. During a second phase of the research effort, a hybrid analytic/numerical approach was developed and evaluated. The approach combines the numerical methods of collocation and the analytical method of regular perturbations. The concept of choosing intelligent interpolating functions is also introduced. Regular perturbation analysis allows the use of a crude representation for the collocation solution, and intelligent interpolating functions further reduce the number of elements without sacrificing the approximation accuracy. As a result, the combined method forms a powerful tool for solving real-time optimal control problems. Details of the approach are illustrated in a fourth order nonlinear example. The hybrid approach is then applied to the launch vehicle problem. The collocation solution is derived from a bilinear tangent steering law, and results in a guidance solution for the entire flight regime that includes both atmospheric and exoatmospheric flight phases.

Calise, Anthony J.↗
NASA NTRS1972

On the algebraic structure of bilinear systems.

It is shown that a particular bilinear model is both quite general and easy to work with. A basic structure theory is developed with the aid of previous results. Some preliminary ideas are discussed together with the system interconnection, the canonical form, questions of controllability, aspects of observability, and equivalent realizations. It is pointed out that in actually determining equivalent realizations for systems and in the classification of systems, the results available in the study of Lie algebras are of fundamental importance.

Brockett, R. W.↗
NASA NTRS1989

Sliding control of pointing and tracking with operator spline estimation

It is shown how a variable structure control technique could be implemented to achieve precise pointing and good tracking of a deformable structure subject to fast slewing maneuvers. The correction torque that has to be applied to the structure is based on estimates of upper bounds on the model errors. For a rapid rotation of the deformable structure, the elastic response can be modeled by oscillators driven by angular acceleration, and where stiffness and damping coefficients are also angular velocity and acceleration dependent. By transforming this slew-driven elastic dynamics into bilinear form (be regarding the vector made up of the angular velocity, squared angular velocity and angular acceleration components, which appear in the coefficients as the input to the deformation dynamics), an operator spline can be constructed, that gives a low order estimate of the induced disturbance. Moreover, a worst case error bound between the estimated deformation and the unknown exact deformation is also generated, which can be used where required in the sliding control correction.

Dwyer, Thomas A. W., III↗
NASA NTRS2009

Improved Fabrication of Ceramic Matrix Composite/Foam Core Integrated Structures

The use of hybridized carbon/silicon carbide (C/SiC) fabric to reinforce ceramic matrix composite face sheets and the integration of such face sheets with a foam core creates a sandwich structure capable of withstanding high-heatflux environments (150 W/cm2) in which the core provides a temperature drop of 1,000 C between the surface and the back face without cracking or delamination of the structure. The composite face sheet exhibits a bilinear response, which results from the SiC matrix not being cracked on fabrication. In addition, the structure exhibits damage tolerance under impact with projectiles, showing no penetration to the back face sheet. These attributes make the composite ideal for leading edge structures and control surfaces in aerospace vehicles, as well as for acreage thermal protection systems and in high-temperature, lightweight stiffened structures. By tailoring the coefficient of thermal expansion (CTE) of a carbon fiber containing ceramic matrix composite (CMC) face sheet to match that of a ceramic foam core, the face sheet and the core can be integrally fabricated without any delamination. Carbon and SiC are woven together in the reinforcing fabric. Integral densification of the CMC and the foam core is accomplished with chemical vapor deposition, eliminating the need for bond-line adhesive. This means there is no need to separately fabricate the core and the face sheet, or to bond the two elements together, risking edge delamination during use. Fibers of two or more types are woven together on a loom. The carbon and ceramic fibers are pulled into the same pick location during the weaving process. Tow spacing may be varied to accommodate the increased volume of the combined fiber tows while maintaining a target fiber volume fraction in the composite. Foam pore size, strut thickness, and ratio of face sheet to core thickness can be used to tailor thermal and mechanical properties. The anticipated CTE for the hybridized composite is managed by the choice of constituents, varying fiber tow sizes and constituent part ratios. This structural concept provides high strength and stiffness at low density 1.06 g/cm3 in panels tested. Varieties of face sheet constructions are possible, including variations in fiber type and weave geometry. The integrated structures possible with this composite could eliminate the need for non-load-bearing thermal protection systems on top of a structural component. The back sheet can readily be integrated to substructures through the incorporation of ribs. This would eliminate weight and cost for aerospace missions.

Hurwitz, Frances I.↗
NASA NTRS2009

Inert Welding/Brazing Gas Filters and Dryers

The use of hybridized carbon/silicon carbide (C/SiC) fabric to reinforce ceramic matrix composite face sheets and the integration of such face sheets with a foam core creates a sandwich structure capable of withstanding high-heat-flux environments (150 W/sq cm) in which the core provides a temperature drop of 1,000 C between the surface and the back face without cracking or delamination of the structure. The composite face sheet exhibits a bilinear response, which results from the SiC matrix not being cracked on fabrication. In addition, the structure exhibits damage tolerance under impact with projectiles, showing no penetration to the back face sheet. These attributes make the composite ideal for leading-edge structures and control surfaces in aerospace vehicles, as well as for acreage thermal protection systems and in high-temperature, lightweight stiffened structures. By tailoring the coefficient of thermal expansion (CTE) of a carbon fiber containing ceramic matrix composite (CMC) face sheet to match that of a ceramic foam core, the face sheet and the core can be integrally fabricated without any delamination. Carbon and SiC are woven together in the reinforcing fabric. Integral densification of the CMC and the foam core is accomplished with chemical vapor deposition, eliminating the need for bond-line adhesive. This means there is no need to separately fabricate the core and the face sheet, or to bond the two elements together, risking edge delamination during use. Fibers of two or more types are woven together on a loom. The carbon and ceramic fibers are pulled into the same "pick" location during the weaving process. Tow spacing may be varied to accommodate the increased volume of the combined fiber tows while maintaining a target fiber volume fraction in the composite. Foam pore size, strut thickness, and ratio of face sheet to core thickness can be used to tailor thermal and mechanical properties. The anticipated CTE for the hybridized composite is managed by the choice of constituents, varying fiber tow sizes and constituent part ratios. This structural concept provides high strength and stiffness at low density 1.06 g/cu cm in panels tested. Varieties of face sheet constructions are possible, including variations in fiber type and weave geometry. The integrated structures possible with this composite could eliminate the need for non-load-bearing thermal protection systems on top of a structural component. The back sheet can readily be integrated to substructures through the incorporation of ribs. This would eliminate weight and cost for aerospace missions.

Goudy, Jerry↗
NASA NTRS1991

A methodology for using nonlinear aerodynamics in aeroservoelastic analysis and design

A methodology is presented for using the Volterra-Wiener theory of nonlinear systems in aeroservoelastic (ASE) analyses and design. The theory is applied to the development of nonlinear aerodynamic response models that can be defined in state-space form and are, therefore, appropriate for use in modern control theory. The theory relies on the identification of nonlinear kernels that can be used to predict the response of a nonlinear system due to an arbitrary input. A numerical kernel identification technique, based on unit impulse responses, is presented and applied to a simple bilinear, single-input single-output (SISO) system. The linear kernel (unit impulse response) and the nonlinear second-order kernel of the system are numerically-identified and compared with the exact, analytically-defined and linear and second-order kernels. This kernel identification technique is then applied to the CAP-TSD (Computational Aeroelasticity Program-Transonic Small Disturbance) code for identification of the linear and second-order kernels of a NACA64A010 rectangular wing undergoing pitch at M = 0.5, M = 8.5 (transonic), and M = 0.93 (transonic). Results presented demonstrate the feasibility of this approach for use with nonlinear, unsteady aerodynamic responses.

Silva, Walter A.↗
NASA NTRS1991

A methodology for using nonlinear aerodynamics in aeroservoelastic analysis and design

A methodology is presented for using the Volterra-Wiener theory of nonlinear systems in aeroservoelastic (ASE) analyses and design. The theory is applied to the development of nonlinear aerodynamic response models that can be defined in state-space form and are, therefore, appropriate for use in modern control theory. The theory relies on the identification of nonlinear kernels that can be used to predict the response of a nonlinear system due to an arbitrary input. A numerical kernel identification technique, based on unit impulse responses, is presented and applied to a simple bilinear, single-input-single-output system. The linear kernel (unit impulse response) and the nonlinear second-order kernel of the system are numerically-identified and compared with the exact, analytically-defined linear and second-order kernels. This kernel identification technique is then applied to the CAP-TSD code for identification of the linear and second-order kernels of a NACA64A010 rectangular wing undergoing pitch at M = 0.5, M = 0.85 (transonic), and M = 0.93 (transonic). Results presented demonstrate the feasibility of this approach for use with nonlinear, unsteady aerodynamic responses.

Silva, Walter A.↗
DOE OSTI2020

Cavitating Acoustics Modeling in Sierra/SD

Motivated by the need to simulate the effects of underwater explosion on ship structures, we develop a new cavitating acoustics formulation. The proposed approach is consistent with existing methods where the cavitation phenomenon is captured with a bilinear constitutive law. However, the new formulation is in terms of velocity potential, as opposed to the existing displacement-potential and pressure formulations. Also unique to the proposed formulation is a new generalized time-stepping procedure specific to cavitating acoustics, which has the ability to introduce numerical damping to control frothing. Numerical examples of varying complexity are presented to illustrate the effectiveness of the proposed approach and the ability to use velocity potential as a primary field variable for cavitating acoustics simulations.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2023

Update on flavor diagonal nucleon charges

This talk provides an update on the calculation of matrix elements of flavor diagonal axial, scalar and tensor quark bilinear operators between the nucleon ground state. The simulations are done using Wilson-clover fermions on a sea of eight 2+1+1-flavor HISQ ensembles generated by the MILC collaboration. We discuss the signal in the sum of the connected and disconnected contributions for the up, down and strange quarks, control over fits to remove excited state contamination, and the simultaneous chiral-continuum fit used to extract the charges.

Park, Sungwoo↗
NASA NTRS1993

Non-oscillatory and non-diffusive solution of convection problems by the iteratively reweighted least-squares finite element method

A comparative description is presented for the least-squares FEM (LSFEM) for 2D steady-state pure convection problems. In addition to exhibiting better control of the streamline derivative than the streamline upwinding Petrov-Galerkin method, numerical convergence rates are obtained which show the LSFEM to be virtually optimal. The LSFEM is used as a framework for an iteratively reweighted LSFEM yielding nonoscillatory and nondiffusive solutions for problems with contact discontinuities; this method is shown to convect contact discontinuities without error when using triangular and bilinear elements.

Jiang, Bo-Nan↗
NASA NTRS1981

Programmer's manual for MMLE3, a general FORTRAN program for maximum likelihood parameter estimation

The MMLE3 is a maximum likelihood parameter estimation program capable of handling general bilinear dynamic equations of arbitrary order with measurement noise and/or state noise (process noise). The basic MMLE3 program is quite general and, therefore, applicable to a wide variety of problems. The basic program can interact with a set of user written problem specific routines to simplify the use of the program on specific systems. A set of user routines for the aircraft stability and control derivative estimation problem is provided with the program. The implementation of the program on specific computer systems is discussed. The structure of the program is diagrammed, and the function and operation of individual routines is described. Complete listings and reference maps of the routines are included on microfiche as a supplement. Four test cases are discussed; listings of the input cards and program output for the test cases are included on microfiche as a supplement.

Maine, R. E.↗
NASA NTRS1998

Interior Noise Reduction by Adaptive Feedback Vibration Control

The objective of this project is to investigate the possible use of adaptive digital filtering techniques in simultaneous, multiple-mode identification of the modal parameters of a vibrating structure in real-time. It is intended that the results obtained from this project will be used for state estimation needed in adaptive structural acoustics control. The work done in this project is basically an extension of the work on real-time single mode identification, which was performed successfully using a digital signal processor (DSP) at NASA, Langley. Initially, in this investigation the single mode identification work was duplicated on a different processor, namely the Texas Instruments TMS32OC40 DSP. The system identification results for the single mode case were very good. Then an algorithm for simultaneous two mode identification was developed and tested using analytical simulation. When it successfully performed the expected tasks, it was implemented in real-time on the DSP system to identify the first two modes of vibration of a cantilever aluminum beam. The results of the simultaneous two mode case were good but some problems were identified related to frequency warping and spurious mode identification. The frequency warping problem was found to be due to the bilinear transformation used in the algorithm to convert the system transfer function from the continuous-time domain to the discrete-time domain. An alternative approach was developed to rectify the problem. The spurious mode identification problem was found to be associated with high sampling rates. Noise in the signal is suspected to be the cause of this problem but further investigation will be needed to clarify the cause. For simultaneous identification of more than two modes, it was found that theoretically an adaptive digital filter can be designed to identify the required number of modes, but the algebra became very complex which made it impossible to implement in the DSP system used in this study. The on-line identification algorithm developed in this research will be useful in constructing a state estimator for feedback vibration control.

Lim, Tae W.↗
NASA NTRS1980

User's manual for MMLE3, a general FORTRAN program for maximum likelihood parameter estimation

A user's manual for the FORTRAN IV computer program MMLE3 is described. It is a maximum likelihood parameter estimation program capable of handling general bilinear dynamic equations of arbitrary order with measurement noise and/or state noise (process noise). The theory and use of the program is described. The basic MMLE3 program is quite general and, therefore, applicable to a wide variety of problems. The basic program can interact with a set of user written problem specific routines to simplify the use of the program on specific systems. A set of user routines for the aircraft stability and control derivative estimation problem is provided with the program.

Maine, R. E.↗
DOE OSTI2025

Effective field theory for radiative corrections to charged-current processes. II. Axial-vector coupling

Here we discuss the hadronic structure-dependent radiative corrections to the axial-vector coupling that controls single-nucleon weak charged-current processes—commonly denoted by 𝑔 𝐴 . We match the Standard Model at the GeV scale onto chiral perturbation theory at next-to-leading order in the one-nucleon sector, in the presence of electromagnetic and weak interactions. As a result, we provide a representation for the corrections to 𝑔 𝐴 in terms of infrared finite convolutions of simple kernels with the single-nucleon matrix elements of time-ordered products of two and three quark bilinears (vector, axial-vector, and pseudoscalar). We discuss strategies to determine the required nonperturbative input from data, lattice-QCD (+QED), and possibly hadronic models. This work paves the way for a precise comparison of the values of the ratio 𝑔 𝐴 ⁡/𝑔 𝑉 extracted from experiment and from lattice QCD, which constrain physics beyond the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

On the Stability of Power Transmission Systems Under Persistent Inverter Attacks: A Bi-Linear Matrix Approach

We investigate the stability and robustness properties of a power transmission system under persistent deceiving attacks on inverter-interfaced energy resources. The attacks can corrupt the damping coefficients in the inverters' controllers and measurements of the frequency at the points of coupling. Leveraging tools from hybrid dynamical systems theory, we characterize a broad family of persistent (and not necessarily periodic) attacks acting on the inverters, under which the stability properties of the transmission system can be shown to not be compromised. To address potentially conservative conditions identified through conventional bounding techniques, sufficient conditions on the average activation time of the attacks are identified via Lyapunov theory, as well as the formulation and solution of a class of bilinear matrix inequalities (BMI). The results are obtained for constant and slowly time-varying loads via input-to-state stability (ISS) tools. Numerical simulations on the IEEE 39-bus test system are also presented.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1994

The Flux-integral Method for Multidimensional Convection and Diffusion

The flux-integral method is a procedure for constructing an explicit, single-step, forward-in-time, conservative, control volume update of the unsteady, multidimensional convection-diffusion equation. The convective plus diffusive flux at each face of a control-volume cell is estimated by integrating the transported variable and its face-normal derivative over the volume swept out by the convecting velocity field. This yields a unique description of the fluxes, whereas other conservative methods rely on nonunique, arbitrary pseudoflux-difference splitting procedures. The accuracy of the resulting scheme depends on the form of the subcell interpolation assumed, given cell-average data. Cellwise constant behavior results in a (very artificially diffusive) first-order convection scheme. Second-order convection-diffusion schemes correspond to cellwise linear (or bilinear) subcell interpolation. Cellwise quadratic subcell interpolants generate a highly accurate convection-diffusion scheme with excellent phase accuracy. Under constant-coefficient conditions, this is a uniformly third-order polynomial interpolation algorithm (UTOPIA).

Leonard, B. P.↗
NASA NTRS1974

Estimation and detection of signals in multiplicative noise

We define a class of detection-estimation problems on matrix Lie groups in which the observation noise is multiplicative in nature. By examining the differential versions of the hypotheses, which are bilinear, we are able to derive the relevant likelihood ratio formula and the associated optimal estimation equations for the signal given the observations and the assumption that the signal is present. These estimation equations are of interest in their own right, in that they represent a finite-dimensional optimal solution to a nonlinear estimation problem and consist of a Kalman-Bucy filter along with the on-line computation of the solution of the associated Riccati equation, which is driven by the observations. The usefulness of these results is illustrated via an example concerning the detection of an actuator failure in a rigid-body rotational control system.

Willsky, A. S.↗