SEARCH · Search NASA

Results for “Complete linearization method”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS2004

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

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

Murphy, Patrick C.↗

DOE OSTI2020

Hole–hole Tamm–Dancoff-approximated density functional theory: A highly efficient electronic structure method incorporating dynamic and static correlation

The study of photochemical reaction dynamics requires accurate as well as computationally efficient electronic structure methods for the ground and excited states. While time-dependent density functional theory (TDDFT) is not able to capture static correlation, complete active space self-consistent field methods neglect much of the dynamic correlation. Hence, inexpensive methods that encompass both static and dynamic electron correlation effects are of high interest. Here, we revisit hole–hole Tamm–Dancoff approximated (hh-TDA) density functional theory for this purpose. The hh-TDA method is the hole–hole counterpart to the more established particle–particle TDA (pp-TDA) method, both of which are derived from the particle–particle random phase approximation (pp-RPA). In hh-TDA, the N-electron electronic states are obtained through double annihilations starting from a doubly anionic (N+2 electron) reference state. In this way, hh-TDA treats ground and excited states on equal footing, thus allowing for conical intersections to be correctly described. Furthermore, the treatment of dynamic correlation is introduced through the use of commonly employed density functional approximations to the exchange-correlation potential. Additionally, we show that hh-TDA is a promising candidate to efficiently treat the photochemistry of organic and biochemical systems that involve several low-lying excited states—particularly those with both low-lying ππ* and nπ* states where inclusion of dynamic correlation is essential to describe the relative energetics. In contrast to the existing literature on pp-TDA and pp-RPA, we employ a functional-dependent choice for the response kernel in pp- and hh-TDA, which closely resembles the response kernels occurring in linear response and collinear spin-flip TDDFT.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Multimode two-dimensional vibronic spectroscopy. I. Orientational response and polarization-selectivity

Two-dimensional Electronic–Vibrational (2D EV) spectroscopy and two-dimensional Vibrational–Electronic (2D VE) spectroscopy are among the newest additions to the coherent multidimensional spectroscopy toolbox, and they are directly sensitive to vibronic couplings. In this first of two papers, the complete orientational response functions are developed for a model system consisting of two coupled anharmonic oscillators and two electronic states in order to simulate polarization-selective 2D EV and 2D VE spectra with arbitrary combinations of linearly polarized electric fields. Here, we propose analytical methods to isolate desired signals within complicated spectra and to extract the relative orientation between vibrational and vibronic dipole moments of the model system using combinations of polarization-selective 2D EV and 2D VE spectral features. Time-dependent peak amplitudes of coherence peaks are also discussed as means for isolating desired signals within the time-domain. This paper serves as a field guide for using polarization-selective 2D EV and 2D VE spectroscopies to map coupled vibronic coordinates on the molecular frame.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1976

A general unified approach to modelling switching-converter power stages

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

Middlebrook, R. D.↗

NASA NTRS1978

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

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

Ruo, S. Y.↗

NASA NTRS1983

Seismic velocity structure of the lunar mantle

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

Nakamura, Y.↗

NASA NTRS1987

A simple method to model truss-beams as equivalent continua

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

Barton, O., Jr.↗

NASA NTRS2007

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

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

Horta, Lucas G.↗

NASA NTRS2009

Precise Stabilization of the Optical Frequency of WGMRs

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

Savchenkov, Anatoliy↗

DOE OSTI2022

The bottom-up EFT: complete UV resonances of the SMEFT operators

The standard model effective field theory (SMEFT) provides systematic parameterization of all possible new physics above the electroweak scale. According to the amplitude-operator correspondence, an effective operator can be decomposed into a linear combination of several j-basis operators, which correspond to local amplitudes carrying certain spin and gauge quantum numbers in a particular scattering channel. Based on the Poincare and gauge symmetries of scattering amplitude, we construct the j-basis using the Casimir method for both the Lorentz and gauge sectors. The quantum numbers of the j-basis operators fix the quantum numbers of any intermediate state in the corresponding amplitudes, such as a UV resonance. This can be re-interpreted as the j-basis/UV correspondence, thus obtaining the j-bases in all partitions of fields for an operator amounts to finding all of its UV origins at tree level, constituting the central part of the bottom-up EFT framework. Applying the j-basis analysis to SMEFT, we obtain a complete list of possible tree-level UV origins of the effective operators at the dimension 5, 6, 7, and all the bosonic operators at dimension 8.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2019

A primer on resurgent transseries and their asymptotics

The computation of observables in general interacting theories, be them quantum mechanical, field, gauge or string theories, is a non-trivial problem which in many cases can only be addressed by resorting to perturbative methods. In most physically interesting problems these perturbative expansions result in asymptotic series with zero radius of convergence. These asymptotic series then require the use of resurgence and transseries in order for the associated observables to become nonperturbatively well-defined. Resurgence encodes the complete large-order asymptotic behaviour of the coefficients from a perturbative expansion, generically in terms of (multi) instanton sectors and for each problem in terms of its Stokes constants. Some observables arise from linear problems, and have a finite number of instanton sectors and associated Stokes constants; some other observables arise from nonlinear problems, and have an infinite number of instanton sectors and Stokes constants. By means of two very explicit examples, and with emphasis on a pedagogical style of presentation, here this work aims at serving as a primer on the aforementioned resurgent, large-order asymptotics of general perturbative expansions. This includes discussions of transseries, Stokes phenomena, generalized steepest-descent methods, Borel transforms, nonlinear resonance, and alien calculus. Furthermore, resurgent properties of transseries – usually described mathematically via alien calculus – are recast in equivalent physical languages: either a “statistical mechanical” language, as motions in chains and lattices; or a “conformal field theoretical” language, with underlying Virasoro-like algebraic structures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Equation-of-Motion Coupled-Cluster Cumulant Green’s Function for Excited States and X-ray Spectra

Green’s function methods within many-body perturbation theory provide a general framework for treating electronic correlations in excited states and spectra. Conventional methods using the Dyson equation or the cumulant expansion are typically based on the GW self-energy approximation. In order to extend this approximation in molecular systems, a non-perturbative real-time coupled-cluster cumulant Green’s function approach has been introduced, where the cumulant is obtained as the solution to a set of coupled first order, non-linear differential equations. This approach naturally includes non-linear corrections to conventional cumulant Green’s function techniques where the cumulant is linear in the GW self-energy. The method yields the spectral function associated with the core Green’s function, which is directly related to the x-ray photoemission spectra (XPS) of molecular systems. The approach also yields very good results for binding energies and satellite excitations. The x-ray absorption spectrum (XAS) is then obtained as a convolution of the core spectral function and an effective one-body absorption spectrum. Here this approach is extended to include the full coupled-cluster-singles (CCS) core Green’s function by including the complete form of the non-linear contributions to the cumulant as well as all single, double, and triple cluster excitations in the CC amplitude equations. This approach naturally builds in orthogonality and shake-up effects analogous to those in the Mahan-Noizeres-de Dominicis edge singularity corrections that enhance the XAS near the edge. Themethod is illustrated for the XPS and XAS of NH 3 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Response of Graphite to Dynamic Loading and Hypervelocity Jet Impacts

The compressive strengths of three varieties of high purity graphite, PCEA, NBG-18, and NBG-25, as well as the depth of penetration of small-scale charges into these materials was experimentally determined. These grades are similar in density, ranging from 1.80 – 1.85 g/cc, and nominal apparent porosity, ranging from 18% to 20%, but provide a wide range in maximum grain or particle size from 10s to 1000s of µm. Two very different manufacturing methods are also represented; PCEA is extruded while NBG-18 and NBG-25 are iso-molded. The quasistatic and dynamic strengths of each grade were determined on a load frame and split-Hopkinson pressure bar, respectively. The depth of penetration (DOP) of two small-scale shaped charges, the Teledyne RP-1 and RP-4, was determined against graphite. The global response of the RP-4 impacts was markedly different as the PCEA samples remained intact while all the NBG-25 samples split into 2 or 3 pieces after the jet penetration had completed. However, for all tests, the trusted DOPs fell within 2 cm. Preliminary hydrocode modeling of the penetration events used existing models that were not designed for graphite. The results can be tuned to reasonably reproduce the DOP, but the wound channel geometry is not reproduced well. A model designed for graphite would need to represent graphite’s non-linear and energy dissipation characteristics.

36 MATERIALS SCIENCE↗

NASA NTRS2002

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

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

Wells, Scott R.↗

NASA NTRS2007

Analytic Guidance for the First Entry in a Skip Atmospheric Entry

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

Garcia-Llama, Eduardo↗

NASA NTRS1978

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

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

Manson, S. S.↗

NASA NTRS1993

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

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

Tamma, Kumar K.↗

DOE OSTI2020

Substation Secondary Asset Health Monitoring and Management System (SSHM)

Electric Power Group, LLC (EPG) was awarded DOE OE0000850 to design, develop, and demonstrate a real-time software application for substation secondary equipment health monitoring and management at host utility American Electric Power (AEP), a cost share partner. This project addresses the need for monitoring substation equipment health and providing operators with tools to identify and take pre-emptive action to avoid catastrophic equipment failure. By monitoring synchrophasor data in real time, data anomalies that indicate potential asset failure can be detected and alerts can be sent to operators in time to take corrective actions. EPG developed two data driven algorithms to identify abnormal equipment signature patterns as well as a method using substation linear state estimation (SLSE). The software has been deployed on hardened PC’s and tested and validated for cost share partner AEP’s substations - 138kV and 765 kV. The SSHM software has been accepted by AEP and final demonstration of DOE - OE0000850 for AEP and DOE was successfully completed on March 17th, 2020. EPG developed the Substation Secondary Asset Health Monitoring (SSHM) Platform to analyze equipment failure signatures in PMU data and alert substation personnel when pre-emptive inspection, repairs and other actions are warranted to prevent potential catastrophic failures. This is the first system that automatically monitors the health of substation secondary assets and alerts users to emerging failures using synchrophasor measurements. SSHM uses high resolution PMU data from PTs, CTs, and CCVTs to analyze equipment signatures and identify anomalies that are precursor indicators of potential equipment failure. When an anomaly is detected, there is a likelihood of potential failure of monitored equipment, the system alerts the user with visual alarms including alarm trend charts and indicator lights on a oneline diagram of the substation.

24 POWER TRANSMISSION AND DISTRIBUTION↗