SEARCH · Engineering Papers

Results for “sum of squares”

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 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Gain weighted eigenspace assignment

This report presents the development of the gain weighted eigenspace assignment methodology. This provides a designer with a systematic methodology for trading off eigenvector placement versus gain magnitudes, while still maintaining desired closed-loop eigenvalue locations. This is accomplished by forming a cost function composed of a scalar measure of error between desired and achievable eigenvectors and a scalar measure of gain magnitude, determining analytical expressions for the gradients, and solving for the optimal solution by numerical iteration. For this development the scalar measure of gain magnitude is chosen to be a weighted sum of the squares of all the individual elements of the feedback gain matrix. An example is presented to demonstrate the method. In this example, solutions yielding achievable eigenvectors close to the desired eigenvectors are obtained with significant reductions in gain magnitude compared to a solution obtained using a previously developed eigenspace (eigenstructure) assignment method.

Davidson, John B.↗

NASA NTRS1994

Active minimization of energy density in three-dimensional enclosures

The objective of this research project is to further investigate and develop a novel approach for actively controlling the sound field in enclosures. Typically the acoustic field in an enclosure has been controlled by minimizing the sum of the squared pressures from several microphones distributed throughout the enclosure. The approach being investigated in this project involves minimizing the acoustic energy density at the sensor locations, rather than the squared pressure. Previous research in a simple one-dimensional enclosure showed that improved global attenuation of the acoustic field is often obtained by minimizing the energy density, rather than the pressure. The current project builds on the previous research by extending the method of controlling the acoustic energy density to three-dimensional enclosures. The results will establish if improved control can still be expected in a more general enclosure. Pending successful results, the method could be applied to control problems such as attenuating the acoustic noise in an aircraft fuselage, an automobile cabin, or other general enclosures. The research project was set up as a two-year project designed to achieve both numerical and experimental results. The primary focus of the first year of research (now being completed) was on the analytical/numerical modeling of the method of controlling the acoustic energy density. During the second year, the research focuses on experimental verification of the approach and extending our understanding of the method.

Sommerfeldt, Scott D.↗

NASA NTRS1993

New dynamic variables for rotating spacecraft

This paper introduces two new seven-parameter representations for spacecraft attitude dynamics modeling. The seven parameters are the three components of the total system angular momentum in the spacecraft body frame; the three components of the angular momentum in the inertial reference frame; and an angle variable. These obey a single constraint as do parameterizations that include a quaternion; in this case the constraint is the equality of the sum of the squares of the angular momentum components in the two frames. The two representations are nonsingular if the system angular momentum is non-zero and obeys certain orientation constraints. The new parameterizations of the attitude matrix, the equations of motion, and the relation of the solution of these equations to Euler angles for torque-free motion are developed and analyzed. The superiority of the new parameterizations for numerical integration is shown in a specific example.

Markley, F. Landis↗

NASA NTRS1996

Active minimization of energy density in three-dimensional enclosures

The objective of this study was to further investigate and develop a novel approach for actively controlling the sound field in enclosures that is based on the acoustic energy density. Typically the acoustic field in an enclosure has been controlled by minimizing the sum of the squared pressures from several microphones distributed throughout the enclosure. The approach investigated in this study involved minimizing the acoustic energy density at the sensor locations, rather than the squared pressure. Research previous to this study in a simple one-dimensional enclosure showed that improved global attenuation of the acoustic field is often obtained by minimizing the energy density, rather than the pressure. The current study built on the previous research by extending the method of controlling the acoustic energy density to three-dimensional enclosures. The study was intended to help establish if improved control can still be expected in a more general enclosure. The study was designed to be both analytical/numerical and experimental in nature.

Sommerfeldt, Scott D.↗

NASA NTRS1996

Dual Brushless Resolver Rate Sensor

This invention relates to dual analog angular rate sensors which are implemented without the use of mechanical brushes. A resolver rate sensor which includes two brushless resolvers which are mechanically coupled to the same output shaft is provided with inputs which are provided to each resolver by providing the first resolver with a DC input and the second resolver with an AC sinusoidal input. A trigonometric identity in which the sum of the squares of the sin and cosine components equal one is used to advantage in providing a sensor of increased accuracy. The first resolver may have a fixed or variable DC input to permit dynamic adjustment of resolver sensitivity thus permitting a wide range of coverage. Novelty and advantages of the invention reside in the excitation of a resolver with a DC signal and in the utilization of two resolvers and the trigonometric identity of cos(exp 2)(theta) + sin(exp 2)(theta) = 1 to provide an accurate rate sensor which is sensitive to direction and accurate through zero rate.

Howard, David E.↗

NASA NTRS1997

Dual Brushless Resolver Rate Sensor

A resolver rate sensor is disclosed in which dual brushless resolvers are mechanically coupled to the same output shaft. Diverse inputs are provided to each resolver by providing the first resolver with a DC input and the second resolver with an AC sinusoidal input. A trigonometric identity in which the sum of the squares of the sin and cosine components equal one is used to advantage in providing a sensor of increased accuracy. The first resolver may have a fixed or variable DC input to permit dynamic adjustment of resolver sensitivity thus permitting a wide range of coverage. In one embodiment of the invention the outputs of the first resolver are directly inputted into two separate multipliers and the outputs of the second resolver are inputted into the two separate multipliers, after being demodulated in a pair of demodulator circuits. The multiplied signals are then added in an adder circuit to provide a directional sensitive output. In another embodiment the outputs from the first resolver is modulated in separate modulator circuits and the output from the modulator circuits are used to excite the second resolver. The outputs from the second resolver are demodulated in separate demodulator circuit and added in an adder circuit to provide a direction sensitive rate output.

Howard, David E.↗

NASA NTRS1996

Adjustment of Sensor Locations During Thermal Property Parameter Estimation

The temperature dependent thermal properties of a material may be evaluated from transient temperature histories using nonlinear parameter estimation techniques. The usual approach is to minimize the sum of the squared errors between measured and calculated temperatures at specific locations in the body. Temperature measurements are usually made with thermocouples and it is customary to take thermocouple locations as known and fixed during parameter estimation computations. In fact, thermocouple locations are never known exactly. Location errors on the order of the thermocouple wire diameter are intrinsic to most common instrumentation procedures (e.g., inserting a thermocouple into a drilled hole) and additional errors can be expected for delicate materials, difficult installations, large thermocouple beads, etc.. Thermocouple location errors are especially significant when estimating thermal properties of low diffusively materials which can sustain large temperature gradients during testing. In the present work, a parameter estimation formulation is presented which allows for the direct inclusion of thermocouple positions into the primary parameter estimation procedure. It is straightforward to set bounds on thermocouple locations which exclude non-physical locations and are consistent with installation tolerances. Furthermore, bounds may be tightened to an extent consistent with any independent verification of thermocouple location, such as x-raying, and so the procedure is entirely consonant with experimental information. A mathematical outline of the procedure is given and its implementation is illustrated through numerical examples characteristic of light-weight, high-temperature ceramic insulation during transient heating. The efficacy and the errors associated with the procedure are discussed.

Milos, Frank S.↗

NASA NTRS1999

Method of Conjugate Radii for Solving Linear and Nonlinear Systems

This paper describes a method to solve a system of N linear equations in N steps. A quadratic form is developed involving the sum of the squares of the residuals of the equations. Equating the quadratic form to a constant yields a surface which is an ellipsoid. For different constants, a family of similar ellipsoids can be generated. Starting at an arbitrary point an orthogonal basis is constructed and the center of the family of similar ellipsoids is found in this basis by a sequence of projections. The coordinates of the center in this basis are the solution of linear system of equations. A quadratic form in N variables requires N projections. That is, the current method is an exact method. It is shown that the sequence of projections is equivalent to a special case of the Gram-Schmidt orthogonalization process. The current method enjoys an advantage not shared by the classic Method of Conjugate Gradients. The current method can be extended to nonlinear systems without modification. For nonlinear equations the Method of Conjugate Gradients has to be augmented with a line-search procedure. Results for linear and nonlinear problems are presented.

Nachtsheim, Philip R.↗

NASA NTRS2006

Equivalent Plate Analysis of Aircraft Wing with Discrete Source Damage

An equivalent plate procedure is developed to provide a computationally efficient means of matching the stiffness and frequencies of flight vehicle wing structures for prescribed loading conditions. First, the equivalent plate is used to match the stiffness of a stiffened panel without damage and the stiffness of a stiffened panel with damage. For both stiffened panels, the equivalent plate models reproduce the deformation of a corresponding detailed model exactly for the given loading conditions. Once the stiffness was matched, the equivalent plate models were then used to predict the frequencies of the panels. Two analytical procedures using the lumped-mass matrix were used to match the first five frequencies of the corresponding detailed model. In both the procedures, the lumped-mass matrix for the equivalent plate is constructed by multiplying the diagonal terms of the consistent-mass matrix by a proportionality constant. In the first procedure, the proportionality constant is selected such that the total mass of the equivalent plate is equal to that of the detailed model. In the second method, the proportionality constant is selected to minimize the sum of the squares of the errors in a set of pre-selected frequencies between the equivalent plate model and the detailed model. The equivalent plate models reproduced the fundamental first frequency accurately in both the methods. It is observed that changing only the mass distribution in the equivalent plate model did not provide enough flexibility to match all of the frequencies.

Krishnamurthy, T.↗

NASA NTRS2005

High Frequency QRS ECG Accurately Detects Cardiomyopathy

High frequency (HF, 150-250 Hz) analysis over the entire QRS interval of the ECG is more sensitive than conventional ECG for detecting myocardial ischemia. However, the accuracy of HF QRS ECG for detecting cardiomyopathy is unknown. We obtained simultaneous resting conventional and HF QRS 12-lead ECGs in 66 patients with cardiomyopathy (EF = 23.2 plus or minus 6.l%, mean plus or minus SD) and in 66 age- and gender-matched healthy controls using PC-based ECG software recently developed at NASA. The single most accurate ECG parameter for detecting cardiomyopathy was an HF QRS morphological score that takes into consideration the total number and severity of reduced amplitude zones (RAZs) present plus the clustering of RAZs together in contiguous leads. This RAZ score had an area under the receiver operator curve (ROC) of 0.91, and was 88% sensitive, 82% specific and 85% accurate for identifying cardiomyopathy at optimum score cut-off of 140 points. Although conventional ECG parameters such as the QRS and QTc intervals were also significantly longer in patients than controls (P less than 0.001, BBBs excluded), these conventional parameters were less accurate (area under the ROC = 0.77 and 0.77, respectively) than HF QRS morphological parameters for identifying underlying cardiomyopathy. The total amplitude of the HF QRS complexes, as measured by summed root mean square voltages (RMSVs), also differed between patients and controls (33.8 plus or minus 11.5 vs. 41.5 plus or minus 13.6 mV, respectively, P less than 0.003), but this parameter was even less accurate in distinguishing the two groups (area under ROC = 0.67) than the HF QRS morphologic and conventional ECG parameters. Diagnostic accuracy was optimal (86%) when the RAZ score from the HF QRS ECG and the QTc interval from the conventional ECG were used simultaneously with cut-offs of greater than or equal to 40 points and greater than or equal to 445 ms, respectively. In conclusion 12-lead HF QRS ECG employing RAZ scoring is a simple, accurate and inexpensive screening technique for cardiomyopathy. Although HF QRS ECG is highly sensitive for cardiomyopathy, its specificity may be compromised in patients with cardiac pathologies other than cardiomyopathy, such as uncomplicated coronary artery disease or multiple coronary disease risk factors. Further studies are required to determine whether HF QRS might be useful for monitoring cardiomyopathy severity or the efficacy of therapy in a longitudinal fashion.

Schlegel, Todd T.↗

NASA NTRS2013

Estimating the Properties of Hard X-Ray Solar Flares by Constraining Model Parameters

We wish to better constrain the properties of solar flares by exploring how parameterized models of solar flares interact with uncertainty estimation methods. We compare four different methods of calculating uncertainty estimates in fitting parameterized models to Ramaty High Energy Solar Spectroscopic Imager X-ray spectra, considering only statistical sources of error. Three of the four methods are based on estimating the scale-size of the minimum in a hypersurface formed by the weighted sum of the squares of the differences between the model fit and the data as a function of the fit parameters, and are implemented as commonly practiced. The fourth method is also based on the difference between the data and the model, but instead uses Bayesian data analysis and Markov chain Monte Carlo (MCMC) techniques to calculate an uncertainty estimate. Two flare spectra are modeled: one from the Geostationary Operational Environmental Satellite X1.3 class flare of 2005 January 19, and the other from the X4.8 flare of 2002 July 23.We find that the four methods give approximately the same uncertainty estimates for the 2005 January 19 spectral fit parameters, but lead to very different uncertainty estimates for the 2002 July 23 spectral fit. This is because each method implements different analyses of the hypersurface, yielding method-dependent results that can differ greatly depending on the shape of the hypersurface. The hypersurface arising from the 2005 January 19 analysis is consistent with a normal distribution; therefore, the assumptions behind the three non- Bayesian uncertainty estimation methods are satisfied and similar estimates are found. The 2002 July 23 analysis shows that the hypersurface is not consistent with a normal distribution, indicating that the assumptions behind the three non-Bayesian uncertainty estimation methods are not satisfied, leading to differing estimates of the uncertainty. We find that the shape of the hypersurface is crucial in understanding the output from each uncertainty estimation technique, and that a crucial factor determining the shape of hypersurface is the location of the low-energy cutoff relative to energies where the thermal emission dominates. The Bayesian/MCMC approach also allows us to provide detailed information on probable values of the low-energy cutoff, Ec, a crucial parameter in defining the energy content of the flare-accelerated electrons. We show that for the 2002 July 23 flare data, there is a 95% probability that Ec lies below approximately 40 keV, and a 68% probability that it lies in the range 7-36 keV. Further, the low-energy cutoff is more likely to be in the range 25-35 keV than in any other 10 keV wide energy range. The low-energy cutoff for the 2005 January 19 flare is more tightly constrained to 107 +/- 4 keV with 68% probability.

X-rays↗

DOE OSTI2023

An Assessment of Persistent Acoustic Monitoring of a Nuclear Reactor during Full Power Generation

Persistent low-frequency (<180 Hz) acoustic detection took place within the boundaries of Oak Ridge National Laboratory to monitor full power operations of the High Flux Isotope Reactor. Three acoustic sensors were installed at distances of 69, 101, and 914 m from the northeast corner of the cooling towers to monitor and assess four reactor power generation cycles. Features were extracted from power spectral density calculations where data were collected during reactor on and off operations. Diverse spectral features were present during full reactor power, including a 21.4 Hz fundamental frequency and ascending harmonics. Using bandpass filters, these related frequencies were isolated and summed, and the root mean square energy was calculated. The method of isolating and summing characteristic features provided a significant improvement in identifying acoustic behavior related to reactor power when the raw signals were obscured by noise.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1982

The Apollo 15 regolith - Chemical modeling and mare/highland mixing

The chemical compositions of thirty-one A-15 soils from nine sampling stations were studied. Mixing models were calculated using five components: 'representative A-15 mare basalt', anorthosite, low-K Fra Mauro basalt, green glass, and medium-K KREEP (Laul and Papike, 1980). A linear-least-squares mixing model (Boynton et al., 1975) was used with sum check and chi-square tests of the fit of components to soil composition. The soils show significant compositional variability, in some instances over relatively small distances within a sampling station. The models indicate that green-glass abundance does not vary systematically with location and that a significant KREEP component is found in all Apennine Front soils. The models also show the soils from the rim of Hadley Rille to be the most mare-enriched, and soils from Station 2 to be the most highland-enriched. The regolith from Station 8-LM is anomalously high in KREEP component relative to the other mare soils, possibly because of ray material.

Walker, R. J.↗

DOE OSTI2020

Wave function analysis with a maximum flow algorithm

An efficient algorithm for computing the maximum-flow path in a network is applied to the identification of the dominant configuration state functions (CSFs) in a graphically contracted function (GCF), configuration interaction, wave function. The flow network is a space of spin-adapted CSFs represented by a Shavitt graph, wherein the nodes correspond to orbital occupations and spin quantum numbers. The graph nodes are connected by arcs, and an arc density is defined as sums of the associated squared CSF coefficients. A max-min approach determines an upper bound to the maximum possible incoming flow for each graph node. A backtracking step generates a candidate walk and is followed by a limited search of alternative branching paths for the dominant CSF. The arc density contributions are removed from the graph, and the algorithm is reapplied to the updated graph. This list of generated walks can be partitioned in order to guarantee that the dominant CSFs have been identified. All of the steps in this algorithm are computationally efficient and do not depend on the potentially large dimension of the underlying linear CSF expansion space. An analysis of low-lying valence states of C-2 illustrates the method.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1972

Simulation of random envelope processes.

Efficient and practical methods of simulating stationary and nonstationary random envelope processes are presented. The stationary envelope processes are simulated by using the fast Fourier transform while the nonstationary envelope processes are simulated as the square root of the sum of a series of cosine functions and a series of sine functions with random phase angles. Typical applications of the envelope simulation are the simulations of peaks and troughs which play an important role in the analyses of the first excursion probability, fatigue and crack propagation. In particular, applications to the crack propagation under random loadings are demonstrated in detail.

Yang, J.-N.↗

NASA NTRS1981

A method for obtaining reduced-order control laws for high-order systems using optimization techniques

A method of synthesizing reduced-order optimal feedback control laws for a high-order system is developed. A nonlinear programming algorithm is employed to search for the control law design variables that minimize a performance index defined by a weighted sum of mean-square steady-state responses and control inputs. An analogy with the linear quadractic Gaussian solution is utilized to select a set of design variables and their initial values. To improve the stability margins of the system, an input-noise adjustment procedure is used in the design algorithm. The method is applied to the synthesis of an active flutter-suppression control law for a wind tunnel model of an aeroelastic wing. The reduced-order controller is compared with the corresponding full-order controller and found to provide nearly optimal performance. The performance of the present method appeared to be superior to that of two other control law order-reduction methods. It is concluded that by using the present algorithm, nearly optimal low-order control laws with good stability margins can be synthesized.

Mukhopadhyay, V.↗

NASA NTRS1985

Optimal generalized multistep integration formulae for real-time digital simulation

The problem of discretizing a dynamical system for real-time digital simulation is considered. Treating the system and its simulation as stochastic processes leads to a statistical characterization of simulator fidelity. A plant discretization procedure based on an efficient matrix generalization of explicit linear multistep discrete integration formulae is introduced, which minimizes a weighted sum of the mean squared steady-state and transient error between the system and simulator outputs.

Moerder, D. D.↗

NASA NTRS1985

Theoretical and experimental investigations of sensor location for optimal aeroelastic system state estimation

One of the major concerns in the design of an active control system is obtaining the information needed for effective feedback. This involves the combination of sensing and estimation. A sensor location index is defined as the weighted sum of the mean square estimation errors in which the sensor locations can be regarded as estimator design parameters. The design goal is to choose these locations to minimize the sensor location index. The choice of the number of sensors is a tradeoff between the estimation quality based upon the same performance index and the total costs of installing and maintaining extra sensors. An experimental study for choosing the sensor location was conducted on an aeroelastic system. The system modeling which includes the unsteady aerodynamics model developed by Stephen Rock was improved. Experimental results verify the trend of the theoretical predictions of the sensor location index for different sensor locations at various wind speeds.

Liu, G.↗