Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix approximation”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Accurate Measurement of Bone Density with QCT

The objective of this study was to determine the accuracy of bone density measurement with a new OCT technology. A phantom was fabricated using two materials, a water-equivalent compound and hydroxyapatite (HA), combined in precise proportions (QRM GrnbH, Germany). The phantom was designed to have the approximate physical size and range in bone density as a human calcaneus, with regions of 0, 50, 100, 200, 400, and 800 mg/cc HA. The phantom was scanned at 80, 120 and 140 KVp with a GE CT/i HiSpeed Advantage scanner. A ring of highly attenuating material (polyvinyl chloride or teflon) was slipped over the phantom to alter the image by introducing non-axi-symmetric beam hardening. Images were corrected with a new OCT technology using an estimate of the effective X-ray beam spectrum to eliminate beam hardening artifacts. The algorithm computes the volume fraction of HA and water-equivalent matrix in each voxel. We found excellent agreement between expected and computed HA volume fractions. Results were insensitive to beam hardening ring material, HA concentration, and scan voltage settings. Data from all 3 voltages with a best fit linear regression are displays.

Cleek, Tammy M.↗
NASA NTRS2004

Hexagonal Pixels and Indexing Scheme for Binary Images

A scheme for resampling binaryimage data from a rectangular grid to a regular hexagonal grid and an associated tree-structured pixel-indexing scheme keyed to the level of resolution have been devised. This scheme could be utilized in conjunction with appropriate image-data-processing algorithms to enable automated retrieval and/or recognition of images. For some purposes, this scheme is superior to a prior scheme that relies on rectangular pixels: one example of such a purpose is recognition of fingerprints, which can be approximated more closely by use of line segments along hexagonal axes than by line segments along rectangular axes. This scheme could also be combined with algorithms for query-image-based retrieval of images via the Internet. A binary image on a rectangular grid is generated by raster scanning or by sampling on a stationary grid of rectangular pixels. In either case, each pixel (each cell in the rectangular grid) is denoted as either bright or dark, depending on whether the light level in the pixel is above or below a prescribed threshold. The binary data on such an image are stored in a matrix form that lends itself readily to searches of line segments aligned with either or both of the perpendicular coordinate axes. The first step in resampling onto a regular hexagonal grid is to make the resolution of the hexagonal grid fine enough to capture all the binaryimage detail from the rectangular grid. In practice, this amounts to choosing a hexagonal-cell width equal to or less than a third of the rectangular- cell width. Once the data have been resampled onto the hexagonal grid, the image can readily be checked for line segments aligned with the hexagonal coordinate axes, which typically lie at angles of 30deg, 90deg, and 150deg with respect to say, the horizontal rectangular coordinate axis. Optionally, one can then rotate the rectangular image by 90deg, then again sample onto the hexagonal grid and check for line segments at angles of 0deg, 60deg, and 120deg to the original horizontal coordinate axis. The net result is that one has checked for line segments at angular intervals of 30deg. For even finer angular resolution, one could, for example, then rotate the rectangular-grid image +/-45deg before sampling to perform checking for line segments at angular intervals of 15deg.

Johnson, Gordon G.↗
NASA NTRS1994

Microstructurally based variations on the dwell fatgue life of titanium alloy IMI 834

An experimental study was undertaken to determine the role of microstructure on the fatigue life reduction observed in titanium alloy IMI 834 under dwell loading conditions. The wave forms compared were a trapezoid with 15 and 30 second hold times at the maximum test load and a baseline, 10 Hertz, haversine. The stress ratio for both loading wave forms was 0.10. The fatigue loading of each specimen was conducted in a vacuum within a scanning electron microscope chamber which minimized the possibility that the laboratory environment would adversely affect the material behavior. Two microstructural conditions were investigated in the experimental program. The first involved standard 'disk' material with equiaxed alpha in a transformed beta matrix. The second material was cut from the same disk forging as the first but was heat treated to obtain a martensitic alpha prime microstructure. Tensile tests were performed prior to the onset of the fatigue loading portion of the study, and it was determined that the yield strengths of the specimens from both material conditions were within ten percent. The maximum fatigue loads were chosen to be 72 percent of the average yield strength for both materials as determined from the tensile tests. It was found that the cycles to failure from the 10 Hertz loading wave form were reduced by a factor of approximately five when the loading was changed to the trapezoidal wave form for the standard 'disk' material. The fatigue life reduction for the martensitic structure under identical test conditions was approximately 1.75. The improvement observed with the martensitic structure also was accompanied by an increase in overall fatigue life for the wave forms tested. This paper will review the results and conclusions of this effort.

Thomsen, Mark L.↗
NASA NTRS1989

The effect of hydrogen on the deformation behavior of a single crystal nickel-base superalloy

The effect of hydrogen on the tensile deformation behavior of PWA 1480 is presented. Tensile tests were interrupted at different plastic strain levels to observe the development of the dislocation structure. Transmission electron microscopy (TEM) foils were cut perpendicular to the tensile axis to allow the deformation of both phases to be simultaneously observed as well as parallel to zone axes (III) to show the superdislocations on their slip planes. Similar to other nickel-base superalloys, hydrogen was detrimental to the room temperature tensile properties of PWA 1480. There was little effect on strength, however the material was severely embrittled. Even without hydrogen, the elongation-to-failure was only approximately 3 percent. The tensile fracture surface was made up primarily of ductile voids with regions of cleavage fracture. These cleavage facets are the eutectic (gamma') in the microstructure. It was shown by quantitative fractography that hydrogen embrittles the eutectic (gamma') and causes the crack path to seek out and fracture through the eutectic (gamma'). There was two to three times the amount of cleavage on the fracture surface of the hydrogen-charged samples than on the surface of the uncharged samples. The effect of hydrogen can also be seen in the dislocation structure. There is a marked tendency for dislocation trapping in the gamma matrix with and without hydrogen at all plastic strain levels. Without hydrogen there is a high dislocation density in the gamma matrix leading to strain exhaustion in this region and failure through the matrix. The dislocation structure at failure with hydrogen is slightly different. The TEM foils cut parallel to zone axes (III) showed dislocations wrapping around gamma precipitates. Zone axes (001) foils show that there is a lower dislocation density in the gamma matrix which can be linked to the effects of hydrogen on the fracture behavior. The primary activity in the gamma precipitates is in the form of superlattice intrinsic stacking faults (SISFs). These faults have also been reported in other ordered alloys and superalloys.

Walston, W. S.↗
NASA NTRS2007

Development of a Tunable Electromechanical Acoustic Liner for Engine Nacelles

This report describes the development of a tunable electromechanical Helmholtz resonator (EMHR) for engine nacelles using smart materials technology. This effort addresses both near-term and long-term goals for tunable electromechanical acoustic liner technology for the Quiet Aircraft Technology (QAT) Program. Analytical models, i.e. lumped element model (LEM) and transfer matrix (TM) representation of the EMHR, have been developed to predict the acoustic behavior of the EMHR. The models have been implemented in a MATLAB program and used to compare with measurement results. Moreover, the prediction performance of models is further improved with the aid of parameter extraction of the piezoelectric backplate. The EMHR has been experimentally investigated using standard two-microphone method (TMM). The measurement results validated both the LEM and TM models of the EMHR. Good agreement between predicted and measured impedance is obtained. Short- and open circuit loads define the limits of the tuning range using resistive and capacitive loads. There is approximately a 9% tuning limit under these conditions for the non-optimized resonator configuration studied. Inductive shunt loads result in a 3 degree-of-freedom DOF) system and an enhanced tuning range of over 20% that is not restricted by the short- and open-circuit limits. Damping coefficient ' measurements for piezoelectric backplates in a vacuum chamber are also performed and indicate that the damping is dominated by the structural damping losses, such as compliant boundaries, and other intrinsic loss mechanisms. Based on models of the EMHR, a Pareto optimization design of the EMHR has been performed for the EMHR with non-inductive loads. The EMHR with non-inductive loads is a 2DOF system with two resonant fiequencies. The tuning ranges of the two resonant frequencies of the EMHR with non-inductive loads cannot be optimized simultaneously; a trade-off (i.e., a Pareto solution) must be reached. The Pareto solution provides the information for a designer that shows how design trade-offs can be used to satisfy specific design requirements. The optimization design of the EMHR with inductive loads aims at optimal tuning of these three resonant fiequencies. The results indicate that it is possible to keep the acoustic reactance of the resonator close to a constant over a given frequency range. An effort to mimic the second layer of the NASA 2DOF liner using a piezoelectric composite diaphragm has been made. The optimal acoustic reactance of the second layer of the NASA 2DOF liner is achieved using a thin PVDF composite diaphragm, but matching the acoustic resistance requires further investigation. Acoustic energy harvesting is achieved by connecting the EMHR to an energy reclamation circuit that converts the ac voltage signal across the piezoceramic to a conditioned dc signal. Energy harvesting experiment yields 16 m W continuous power for an incident SPL of 153 dB. Such a level is sufficient to power a variety of low power electronic devices. Finally, technology transfer has been achieved by converting the original NASA ZKTL FORTRAN code to a MATLAB code while incorporating the models of the EMHR. Initial studies indicate that the EMHR is a promising technology that may enable lowpower, light weight, tunable engine nacelle liners. This technology, however, is very immature, and additional developments are required. Recommendations for future work include testing of sample EMHR liner designs in NASA Langley s normal incidence dual-waveguide and the grazing-incidence flow facility to evaluating both the impedance characteristics as well as the energy reclamation abilities. Additional design work is required for more complex tuning circuits with greater performance. Poor electromechanical coupling limited the electromechanical tuning capabilities of the proof of concept EMHR. Different materials than those studies and perhaps novel composite material systems may dramatically improvehe electromechanical coupling. Such improvements are essential to improved mimicking of existing double layer liners.

Liu, Fei↗
NASA NTRS2009

New insights into the origin of magnetite crystals in ALH84001 carbonate disks

Martian meteorite ALH84001 preserves evidence of interaction with aqueous fluids while on Mars in the form of microscopic carbonate disks which are believed to have precipitated approximately 3.9 Ga ago at beginning of the Noachian epoch. Intimately associated within and throughout these carbonate disks are nanocrystal magnetites (Fe3O4) with unusual chemical and physical properties, whose origins have become the source of considerable debate. One group of hypotheses argues that these Fe3O4 are the product of partial thermal decomposition of the host carbonate. Alternatively, the origins of Fe3O4 and carbonate may be unrelated; that is, from the perspective of the carbonate the magnetite is allochthonous. We have sought to resolve between these hypotheses through the detailed characterized of the compositional and structural relationships of the carbonate disks and associated magnetites with the orthopyroxene matrix in which they are embedded. We focus this discussion on the composition of ALH84001 magnetites and then compare these observations with those from experimental thermal decomposition studies of sideritic carbonates under a range of plausible geological heating scenarios.

Thomas-Keprta, K. L.↗
NASA NTRS2004

Aluminum-Alloy-Matrix/Alumina-Reinforcement Composites

Isotropic composites of aluminum-alloy matrices reinforced with particulate alumina have been developed as lightweight, high-specific-strength, less-expensive alternatives to nickel-base and ferrous superalloys. These composites feature a specific gravity of about 3.45 grams per cubic centimeter and specific strengths of about 200 MPa/(grams per cubic centimeter). The room-temperature tensile strength is 100 ksi (689 MPa) and stiffness is 30 Msi (206 GPa). At 500 F (260 C), these composites have shown 80 percent retention in strength and 95 percent retention in stiffness. These materials also have excellent fatigue tolerance and tribological properties. They can be fabricated in net (or nearly net) sizes and shapes to make housings, pistons, valves, and ducts in turbomachinery, and to make structural components of such diverse systems as diesel engines, automotive brake systems, and power-generation, mining, and oil-drilling equipment. Separately, incorporation of these metal matrix composites within aluminum gravity castings for localized reinforcement has been demonstrated. A composite part of this type can be fabricated in a pressure infiltration casting process. The process begins with the placement of a mold with alumina particulate preform of net or nearly net size and shape in a crucible in a vacuum furnace. A charge of the alloy is placed in the crucible with the preform. The interior of the furnace is evacuated, then the furnace heaters are turned on to heat the alloy above its liquidus temperature. Next, the interior of the furnace is filled with argon gas at a pressure about 900 psi (approximately equal to 6.2 MPa) to force the molten alloy to infiltrate the preform. Once infiltrated, the entire contents of the crucible can be allowed to cool in place, and the composite part recovered from the mold.

Kashalikar, Uday↗
NASA NTRS2012

Wide Range SET Pulse Measurement

A method for measuring a wide range of SET pulses is demonstrated. Use of dynamic logic, faster than ordinary CMOS, allows capture of short pulses. A weighted binning of SET lengths allows measurement of a wide range of pulse lengths with compact circuitry. A pulse-length-conservative pulse combiner tree routes SETs from combinational logic to the measurement circuit, allowing SET measurements in circuits that cannot easily be arranged in long chains. The method is applied to add-multiplex combinational logic, and to an array of NFET routing switches, at .35 micron. Pulses are captured in a chain of Domino Logic AND gates. Propagation through the chain is frozen on the trailing edge by dropping low the second "enable" input to the AND gates. Capacitive loading is increased in the latter stages to create an approximately logarithmic weighted binning, so that a broad range of pulse lengths can be captured with a 10 stage capture chain. Simulations show pulses can be captured which are 1/5th the length of those typically captured with leading edge triggered latch methods, and less than the length of those captured with a trailing edge latch method. After capture, the pulse pattern is transferred to an SEU protected shift register for readout. 64 instances of each of two types of logic are used as targets. One is a full adder with a 4 to 1 mux on its inputs. The other is a 4 x 4 NFET routing matrix. The outputs are passed through buffered XNOR comparators to identify pulses, which are merged in a buffered not-nand (OR) tree designed to avoid pulse absorption as much as possible. The output from each of the two test circuits are input into separate pulse measurement circuits. Test inputs were provided so that the circuit could be bench tested and calibrated. A third SET measurement circuit with no inputs was used to judge the contribution from direct hits on the measurement circuit. Heavy ions were used with an LET range from 12 to 176. At LET of 21 and below, the very small number of SETs were not significantly higher in the test over the control circuits. At higher LET the test circuit SETs are one or two orders of magnitude greater than for the control circuit. The NFET circuit produces more and slightly longer SETs as expected. But the differences do not appear to be significant enough to modify strategies now used to avoid capture of SETs in chips such as FPGAs. Complete data and graphs will be in the full paper / presentation. In the summary figure below left, NOCL is the reference circuit without any input, and number of stages triggered is plotted. Simulation at right shows the smallest pulse captured (stage 2) at about 300 ps. Our conclusion is that the method is promising, but that improvements in the merge network are desirable before applying in a deep submicron process

Shuler, Robert L.↗
NASA NTRS1988

Improving the accuracy of central difference schemes

General difference approximations to the fluid dynamic equations require an artificial viscosity in order to converge to a steady state. This artificial viscosity serves two purposes. One is to suppress high frequency noise which is not damped by the central differences. The second purpose is to introduce an entropy-like condition so that shocks can be captured. These viscosities need a coefficient to measure the amount of viscosity to be added. In the standard scheme, a scalar coefficient is used based on the spectral radius of the Jacobian of the convective flux. However, this can add too much viscosity to the slower waves. Hence, it is suggested that a matrix viscosity be used. This gives an appropriate viscosity for each wave component. With this matrix valued coefficient, the central difference scheme becomes closer to upwind biased methods.

Turkel, Eli↗
NASA NTRS1993

Supersonic flow visualization of a nacelle in close proximity to a simulated wing

A flow visualization study was made in the 9 x 9 inch supersonic wind tunnel at Wichita State University to examine shock and boundary layer flow interaction for a nacelle in close proximity to the lower surface of a simulated wing. The test matrix included variations of angle of attack from -2 degrees to +4 degrees, nacelle-wing gap from 0.5 to 3-nacelle inlet diameter (0.12 inch), and Reynolds number based on nacelle length (1.164 inch) from 1.16 x 10(exp 6) to 1.45 x 10(exp 6) at a nominal Mach number of 2. Schlieren pictures of wing and nacelle flowfield were recorded by a video camera during each tunnel run. Results show that the nacelle inlet shock wave remains attached to the inlet lip and its impingement does not significantly affect the wing boundary layer. At the nacelle trailing edge location, the wing boundary layer thickness is approximately one nacelle inlet diameter at alpha = 0 degrees and it decreases with increase of angle of attack.

Biber, Kasim↗
NASA NTRS1997

Microphysical Retrievals Over Stratiform Rain Using Measurements from an Airborne Dual-Wavelength Radar-Radiometer

The need to understand the complementarity of the radar and radiometer is important not only to the Tropical Rain Measuring Mission (TRMM) program but to a growing number of multi-instrumented airborne experiment that combine single or dual-frequency radars with multichannel radiometers. The method of analysis used in this study begins with the derivation of dual-wavelength radar equations for the estimation of a two-parameter drop size distribution (DSD). Defining a "storm model" as the set of parameters that characterize snow density, cloud water, water vapor, and features of the melting layer, then to each storm model there will usually correspond a set of range-profiled drop size distributions that are approximate solutions of the radar equations. To test these solutions, a radiative transfer model is used to compute the brightness temperatures for the radiometric frequencies of interest. A storm model or class of storm models is considered optimum if it provides the best reproduction of the radar and radiometer measurements. Tests of the method are made for stratiform rain using simulated storm models as well as measured airborne data. Preliminary results show that the best correspondence between the measured and estimated radar profiles usually can be obtained by using a moderate snow density (0.1-0.2 g/cu cm), the Maxwell-Garnett mixing formula for partially melted hydrometeors (water matrix with snow inclusions), and low to moderate values of the integrated cloud liquid water (less than 1 kg/sq m). The storm-model parameters that yield the best reproductions of the measured radar reflectivity factors also provide brightness temperatures at 10 GHz that agree well with the measurements. On the other hand, the correspondence between the measured and modeled values usually worsens in going to the higher frequency channels at 19 and 34 GHz. In searching for possible reasons for the discrepancies, It is found that changes in the DSD parameter Mu, the radar constants, or the path-integrated attenuation can affect the high frequency channels significantly. In particular, parameters that cause only modest increases in the median mass diameter of the snow, and which have a minor effect on the radar returns or the low frequency brightness temperature, can produce a strong cooling of the 34 GHz brightness temperature.

Meneghini, Robert↗
NASA NTRS2004

Integrated-Optic Oxygen Sensors

Compact optical oxygen sensors with self-calibration capabilities are undergoing development. A sensor of this type features a single-chip, integrated-optic design implemented by photolithographic fabrication of optical waveguides in a photosensitive porous glass. The porosity serves as both a matrix for retention of an oxygen-sensitive fluorescent indicator chemical and a medium for diffusion of oxygen to the chemical from the ambient air to be monitored. Each sensor includes at least one such waveguide exposed to the atmosphere and at least one covered with metal for isolation from the atmosphere. The covered one serves as a reference channel. In operation, the concentration of oxygen is deduced from the intensity and lifetime of the fluorescence in the exposed channel, with the help of calibration data acquired via the reference channel. Because the sensory chemical is placed directly in and throughout the cross section of the light path, approximately 99 percent of the light in the waveguide is available for interaction with the chemical, in contradistinction to only about 1 percent of the light in an optical sensor that utilizes evanescentwave coupling. Hence, a sensor of this type is significantly more sensitive.

Mendoza, Edgar A.↗
NASA NTRS1996

Methodology for Sensitivity Analysis, Approximate Analysis, and Design Optimization in CFD for Multidisciplinary Applications

An incremental iterative formulation together with the well-known spatially split approximate-factorization algorithm, is presented for solving the large, sparse systems of linear equations that are associated with aerodynamic sensitivity analysis. This formulation is also known as the 'delta' or 'correction' form. For the smaller two dimensional problems, a direct method can be applied to solve these linear equations in either the standard or the incremental form, in which case the two are equivalent. However, iterative methods are needed for larger two-dimensional and three dimensional applications because direct methods require more computer memory than is currently available. Iterative methods for solving these equations in the standard form are generally unsatisfactory due to an ill-conditioned coefficient matrix; this problem is overcome when these equations are cast in the incremental form. The methodology is successfully implemented and tested using an upwind cell-centered finite-volume formulation applied in two dimensions to the thin-layer Navier-Stokes equations for external flow over an airfoil. In three dimensions this methodology is demonstrated with a marching-solution algorithm for the Euler equations to calculate supersonic flow over the High-Speed Civil Transport configuration (HSCT 24E). The sensitivity derivatives obtained with the incremental iterative method from a marching Euler code are used in a design-improvement study of the HSCT configuration that involves thickness. camber, and planform design variables.

Taylor, Arthur C., III↗
NASA NTRSDate not supplied

Towards Formalization of Advanced Linear Algebra with Applications to Dynamical Systems using PVS

Linear Algebra is essential for numerous aerospace problems of interest. Formal reasoning about hybrid systems that contain variables modeled by differential equations rely on concepts from Linear Algebra such as eigenvalues, matrix decompositions, and matrix valued functions. For example, the long-term dynamics of a system of differential equations depend on the stability/instability of its equilibrium points, which often reduces to an eigenvalue problem. This talk will embark on a quest to formalize theorems and results about eigenvalues and eigenvectors using PVS. We shall start our journey with 2 x 2 complex matrices, where we will apply our PVS code to a simple example of a dynamical system. Since it can be difficult or impossible to give simple expressions of eigenvalues for larger matrices (i.e. 5 x 5 or higher), we then move towards specifying the power method for verified computation of eigenvalue approximations in PVS. This effort requires development of multivariate complex arithmetic. At the end of the day, having such additions to the PVS NASA libraries will help move towards the use of formal methods to verify concepts of control theory and system level verification.

Linear Algebra↗
NASA NTRS2012

Characterization and Analysis of Triaxially Braided Polymer Composites under Static and Impact Loads

In order to design impact resistant aerospace components made of triaxially-braided polymer matrix composite materials, a need exists to have reliable impact simulation methods and a detailed understanding of the material behavior. Traditional test methods and specimen designs have yielded unrealistic material property data due to material mechanisms such as edge damage. To overcome these deficiencies, various alternative testing geometries such as notched flat coupons have been examined to alleviate difficulties observed with standard test methods. The results from the coupon level tests have been used to characterize and validate a macro level finite element-based model which can be used to simulate the mechanical and impact response of the braided composites. In the analytical model, the triaxial braid unit cell is approximated by using four parallel laminated composites, each with a different fiber layup, which roughly simulates the braid architecture. In the analysis, each of these laminated composites is modeled as a shell element. Currently, each shell element is considered to be a smeared homogeneous material. Simplified micromechanics techniques and lamination theory are used to determine the equivalent stiffness properties of each shell element, and results from the coupon level tests on the braided composite are used to back out the strength properties of each shell element. Recent improvements to the model include the incorporation of strain rate effects into the model. Simulations of ballistic impact tests have been carried out to investigate and verify the analysis approach.

Goldberg, Robert K.↗
NASA NTRS1976

Two-electron bond-orbital model, 2

The two-electron bond-orbital model of tetrahedrally-coordinated solids is generalized and its application extended. All intrabond matrix elements entering the formalism are explicitly retained, including the direct overlap S between the anion and cation sp3 hybrid wavefunctions. Complete analytic results are obtained for the six two-electron eigenvalues and eigenstates of the anion-cation bond in terms of S, one-electron parameters V2 and V3, and two-electron correlation parameters V4, V5 and V6. Refined formulas for the dielectric constant and the nuclear exchange and pseudodipolar coefficients, as well as new expressions for the valence electron density, polarity of the bond and the cohesive energy, are then derived. The theory gives a good account of the experimentally observed trends in all properties considered and approximate quantitative agreement is achieved for the pseudodipolar coefficient.

Huang, C.↗
NASA NTRS1992

Closed-loop transfer recovery with observer-based controllers. I - Analysis. II - Design

A detailed study is presented of three fundamental issues related to the problem of closed-loop transfer (CLT) recovery. The first issues concerns what can and cannot be achieved for a given system and for an arbitrary target CLT function (TCLTF). The second issue involves developing necessary and/or sufficient conditions for a TCLTF to be recoverable either exactly or approximately. The third issue involves the necessary and/or sufficient conditions on a given system such that it has at least one recoverable TCLTF. The results of the analysis identify some fundamental limitations of the given system as a consequence of its structural properties which enables designers to appreciate at the outset different design limitations incurred in the synthesis of output-feedback controllers. Then, the actual design of full-order or reduced-order observer-based controllers is addressed which will achieve as close as possibly the desired TCLTF. Three design methods are considered: (1) the ATEA method, (2) a method that minimizes the H2-norm of a recovery matrix, and (3) a method that minimizes the respective H(infinity) norm. The relative merits of the methods are discussed.

Chen, Ben M.↗
NASA NTRS1995

Galactic oscillator symmetry

Riemann ellipsoids model rotating galaxies when the galactic velocity field is a linear function of the Cartesian coordinates of the galactic masses. In nuclear physics, the kinetic energy in the linear velocity field approximation is known as the collective kinetic energy. But, the linear approximation neglects intrinsic degrees of freedom associated with nonlinear velocity fields. To remove this limitation, the theory of symplectic dynamical symmetry is developed for classical systems. A classical phase space for a self-gravitating symplectic system is a co-adjoint orbit of the noncompact group SP(3,R). The degenerate co-adjoint orbit is the 12 dimensional homogeneous space Sp(3,R)/U(3), where the maximal compact subgroup U(3) is the symmetry group of the harmonic oscillator. The Hamiltonian equations of motion on each orbit form a Lax system X = (X,F), where X and F are elements of the symplectic Lie algebra. The elements of the matrix X are the generators of the symplectic Lie algebra, viz., the one-body collective quadratic functions of the positions and momenta of the galactic masses. The matrix F is composed from the self-gravitating potential energy, the angular velocity, and the hydostatic pressure. Solutions to the hamiltonian dynamical system on Sp(3,R)/U(3) are given by symplectic isospectral deformations. The Casimirs of Sp(3,R), equal to the traces of powers of X, are conserved quantities.

Rosensteel, George↗