Search NASA⌕ Search

SEARCH · Search NASA

Results for “Variational formulation”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Template Matching Used for Small Body Optical Navigation with Poorly Detailed Objects

Object and template matching becomes difficult when an image lacks detail. This is particularly worrisome when typical matching techniques, cross-correlation, log-polar mapping, and key point matching fail. Work herein describes a formulation that identifies objects of interest, estimates the affine transformation between a template object and scene using Principal Component Analysis (PCA), and provides a fit value for the objects and template incorporating Hu's Moments. The algorithm presented is tested on synthetic images and images obtained from the OSIRIS-REx mission while the spacecraft was approaching its target, Bennu. Results for the current formulation show that, with the presence of large-scale variations and rotation, the fitting scheme performs well when compared with other techniques.

Lyzhoft, Joshua R.↗
NASA NTRS2024

Optimal Control Allocation for Distributed Electric Propulsion in a Series/Parallel Partial Hybrid Powertrain

The SUbsonic Single Aft eNgine (SUSAN) Electrofan is a NASA concept jet transport aircraft with a 2040 entry-into-service date. It utilizes electrified aircraft propulsion (EAP) to enable propulsive and aerodynamic benefits to reduce fuel usage, emissions, and cost. The powertrain consists of a single thrust producing, boundary layer-ingesting (BLI) turbofan gas turbine engine (GTE) with generators driving a series/parallel partial hybrid EAP system. The architecture includes 16 underwing contrarotating BLI fans, eight on each side, in a mailslot configuration. The 16 fans run on power extracted from the GTE through four 5 MW motor/generators connected to the Low-Pressure Spool, and a single 1 MW motor/generator on the High-Pressure Spool. The distributed fans can be used by the flight control to augment or replace the rudder function. At top of climb, the power extracted from the GTE for the fans is boosted by batteries. The design provides redundancy, and the capacity for boost means that the fans are designed to be able to provide additional thrust when necessary. These features can be leveraged in case of a fan or generator failure. This paper sets up the optimal control problem of setpoint determination for individual fans in the distributed propulsion system, accounting for electrical string efficiencies, saturations, and failures. The solution minimizes power consumption while maintaining thrust and torque on the airframe for maneuvering. Additionally, thrust that would have been lost due to temporary fan speed or power saturation is optimally redistributed to maintain overall desired thrust and torque on the aircraft. The power extraction range constraints derive from the gas turbine engine design and the small amount of variation allowed for the engine to maintain operability. The problem formulation allows the number and location of fan failures for which the thrust and torque can be maintained to be investigated, which has implications for certification. Simulations of a coordinated turn utilizing the distributed electric propulsion for yaw rate control under different failure scenarios demonstrate the robustness of the powertrain design to failures and help define its limitations.

Distributed Electric Propulsion↗
NASA NTRS1991

Joint inversion of shear wave travel time residuals and geoid and depth anomalies for long-wavelength variations in upper mantle temperature and composition along the Mid-Atlantic Ridge

Measurements were carried out for SS-S differential travel time residuals for nearly 500 paths crossing the northern Mid-Atlantic Ridge, assuming that the residuals are dominated by contributions from the upper mantle near the surface bounce point of the reflected phase SS. Results indicate that the SS-S travel time residuals decrease linearly with square root of age, to an age of 80-100 Ma, in general agreement with the plate cooling model. A joint inversion was formulated of travel time residuals and geoid and bathymetric anomalies for lateral variation in the upper mantle temperature and composition. The preferred inversion solutions were found to have variations in upper mantle temperature along the Mid-Atlantic Ridge of about 100 K. It was calculated that, for a constant bulk composition, such a temperature variation would produce about a 7-km variation in crustal thickness, larger than is generally observed.

Sheehan, Anne F.↗
NASA NTRS1991

Integrated flight/propulsion control design for a STOVL aircraft using H-infinity control design techniques

Results are presented from an application of H(infinity) control design methodology to a centralized integrated flight/propulsion control (IFPC) system design for a supersonic short take-off and vertical landing (STOVL) fighter aircraft in transition flight. The emphasis is on formulating the H(infinity) control design problem such that the resulting controller provides robustness to modeling uncertainties and model parameter variations with flight condition. Experience gained from a preliminary H(infinity)=based IFPC design study performed earlier is used as the base to formulate the robust H(infinity) control design problem and improve the previous design. Detailed evaluation results are presented for a reduced-order controller obtained from the improved H(infinity) control design showing that the control design meets the specified nominal performance objectives as well as provides stability robustness for variations in plant system dynamics with chnages in aircraft trim speed within the transition flight envelope.

Garg, Sanjay↗
NASA NTRS2004

An Interface Damage Model for the Simulation of Delamination Under Variable-Mode Ratio in Composite Materials

A thermodynamically consistent damage model for the simulation of progressive delamination under variable mode ratio is presented. The model is formulated in the context of the Damage Mechanics (DM). The constitutive equations that result from the variation of the free energy with damage are used to model the initiation and propagation of delamination. A new delamination initiation criterion is developed to assure that the formulation can account for changes in the loading mode in a thermodynamically consistent way. Interfacial penetration of two adjacent layers after complete decohesion is prevented by the formulation of the free energy. The model is implemented into the commercial finite element code ABAQUS by means of a user-written decohesion element. Finally, the numerical predictions given by the model are compared with experimental results.

Turon, Albert↗
NASA NTRS1982

Alternative ways for formulation of hybrid stress elements

An element stiffness matrix can be derived by the conventional potential energy principle and, indirectly, also by generalized variational principles, such as the Hu-Washizu principle and the Hellinger-Reissner principle. The present investigation has the objective to show an approach which is concerned with the formulation of incompatible elements for solid continuum and for plate bending problems by the Hellinger-Reissner principle. It is found that the resulting scheme is equivalent to that considered by Tong (1982) for the construction of hybrid stress elements. In Tong's scheme the inversion of a large flexibility matrix can be avoided. It is concluded that the introduction of additional internal displacement modes in mixed finite element formulations by the Hellinger-Reissner principle and the Hu-Washizu principle can lead to element stiffness matrices which are equivalent to the assumed stress hybrid method.

Pian, T. H. H.↗
NASA NTRS1991

Integrated flight/propulsion control design for a STOVL aircraft using H-infinity control design techniques

Results are presented from an application of H-infinity control design methodology to a centralized integrated flight propulsion control (IFPC) system design for a supersonic Short Takeoff and Vertical Landing (STOVL) fighter aircraft in transition flight. The emphasis is on formulating the H-infinity control design problem such that the resulting controller provides robustness to modeling uncertainties and model parameter variations with flight condition. Experience gained from a preliminary H-infinity based IFPC design study performed earlier is used as the basis to formulate the robust H-infinity control design problem and improve upon the previous design. Detailed evaluation results are presented for a reduced order controller obtained from the improved H-infinity control design showing that the control design meets the specified nominal performance objectives as well as provides stability robustness for variations in plant system dynamics with changes in aircraft trim speed within the transition flight envelope. A controller scheduling technique which accounts for changes in plant control effectiveness with variation in trim conditions is developed and off design model performance results are presented.

Garg, Sanjay↗
NASA NTRS1990

A 4-node assumed-stress hybrid shell element with rotational degrees of freedom

An assumed-stress hybrid/mixed 4-node quadrilateral shell element is introduced that alleviates most of the deficiencies associated with such elements. The formulation of the element is based on the assumed-stress hybrid/mixed method using the Hellinger-Reissner variational principle. The membrane part of the element has 12 degrees of freedom including rotational or drilling degrees of freedom at the nodes. The bending part of the element also has 12 degrees of freedom. The bending part of the element uses the Reissner-Mindlin plate theory which takes into account the transverse shear contributions. The element formulation is derived from an 8-node isoparametric element. This process is accomplished by assuming quadratic variations for both in-plane and out-of-plane displacement fields and linear variations for both in-plane and out-of-plane rotation fields along the edges of the element. In addition, the degrees of freedom at midside nodes are approximated in terms of the degrees of freedom at corner nodes. During this process the rotational degrees of freedom at the corner nodes enter into the formulation of the element. The stress field are expressed in the element natural-coordinate system such that the element remains invariant with respect to node numbering.

Aminpour, Mohammad A.↗
NASA NTRS1974

McLaughlin and Mars

McLaughlin formulated a novel explanation for the Martian albedo markings and their variations. This explanation, the volcanic-aeolian hypothesis, is one of the very few prespacecraft views of planet-wide phenomena on Mars which have stood the test of time. The distribution of albedo markings is considered along with secular changes on Mars and seasonal changes in region Syrtis Major. Mariner 9 has demonstrated that aeolian transport is the dominant factor in determining the distribution of albedo markings on Mars.

Veverka, J.↗
NASA NTRS1983

A mixed finite element for the analysis of laminated plates

A new mixed shear-flexible finite element based on the Hellinger-Reissner's variational principle is developed. The element is constructed using a mixed formulation of the shear deformation theory of laminated composite plates, and consists of three displacements, two shear rotations, and three moments as the independent degrees of freedom. The numerical convergence and accuracy characteristics of the element are investigated for bending of laminated anisotropic composite plates. The element is relatively simple to construct and has better accuracy and convergence features when compared to other conventional finite elements.

Putcha, N. S.↗
NASA NTRS1986

Implicit multigrid algorithms for the three-dimensional flux split Euler equations

The full approximation scheme multigrid method is applied to several implicit flux-split algorithms for solving the three-dimensional Euler equations in a body fitted coordinate system. Each uses a variation of approximate factorization and is implemented in a finite volume formulation. The algorithms are all vectorizable with little or no scalar computations required. The flux vectors are split into upwind components using both the splittings of Steger-Warming and Van Leer. Results comparing pressure distributions with experimental data using both splitting types are shown. The stability and smoothing rate of each of the schemes are examined using a Fourier analysis of the complete system of equations. Results are presented for three-dimensional subsonic, transonic, and supersonic flows which demonstrate substantially improved convergence rates with the multigrid algorithm. The influence of using both a V-cycle and a W-cycle on the convergence is examined. Using the multigrid method on both subsonic and transonic wing calculations, the final lift coefficient is obtained to within 0.1 percent of its final value in a few as 15 cycles for a mesh with over 210,000 points. A spectral radius of 0.89 is achieved for both subsonic and transonic flow over the ONERA M6 wing while a spectral radius of 0.83 is obtained for supersonic flow over an analytically defined forebody. Results compared with experiment for all cases show good agreement.

Anderson, W. K.↗
NASA NTRS1988

Three-dimensional multigrid algorithms for the flux-split Euler equations

The Full Approximation Scheme (FAS) multigrid method is applied to several implicit flux-split algorithms for solving the three-dimensional Euler equations in a body fitted coordinate system. Each of the splitting algorithms uses a variation of approximate factorization and is implemented in a finite volume formulation. The algorithms are all vectorizable with little or no scalar computation required. The flux vectors are split into upwind components using both the splittings of Steger-Warming and Van Leer. The stability and smoothing rate of each of the schemes are examined using a Fourier analysis of the complete system of equations. Results are presented for three-dimensional subsonic, transonic, and supersonic flows which demonstrate substantially improved convergence rates with the multigrid algorithm. The influence of using both a V-cycle and a W-cycle on the convergence is examined.

Anderson, W. Kyle↗
NASA NTRS1989

A weak Hamiltonian finite element method for optimal control problems

A temporal finite element method based on a mixed form of the Hamiltonian weak principle is developed for dynamics and optimal control problems. The mixed form of Hamilton's weak principle contains both displacements and momenta as primary variables that are expanded in terms of nodal values and simple polynomial shape functions. Unlike other forms of Hamilton's principle, however, time derivatives of the momenta and displacements do not appear therein; instead, only the virtual momenta and virtual displacements are differentiated with respect to time. Based on the duality that is observed to exist between the mixed form of Hamilton's weak principle and variational principles governing classical optimal control problems, a temporal finite element formulation of the latter can be developed in a rather straightforward manner. Several well-known problems in dynamics and optimal control are illustrated. The example dynamics problem involves a time-marching problem. As optimal control examples, elementary trajectory optimization problems are treated.

Hodges, Dewey H.↗
NASA NTRS1991

Weak Hamiltonian finite element method for optimal control problems

A temporal finite element method based on a mixed form of the Hamiltonian weak principle is developed for dynamics and optimal control problems. The mixed form of Hamilton's weak principle contains both displacements and momenta as primary variables that are expanded in terms of nodal values and simple polynomial shape functions. Unlike other forms of Hamilton's principle, however, time derivatives of the momenta and displacements do not appear therein; instead, only the virtual momenta and virtual displacements are differentiated with respect to time. Based on the duality that is observed to exist between the mixed form of Hamilton's weak principle and variational principles governing classical optimal control problems, a temporal finite element formulation of the latter can be developed in a rather straightforward manner. Several well-known problems in dynamics and optimal control are illustrated. The example dynamics problem involves a time-marching problem. As optimal control examples, elementary trajectory optimization problems are treated.

Hodges, Dewey H.↗
NASA NTRS1990

A weak Hamiltonian finite element method for optimal control problems

A temporal finite element method based on a mixed form of the Hamiltonian weak principle is developed for dynamics and optimal control problems. The mixed form of Hamilton's weak principle contains both displacements and momenta as primary variables that are expanded in terms of nodal values and simple polynomial shape functions. Unlike other forms of Hamilton's principle, however, time derivatives of the momenta and displacements do not appear therein; instead, only the virtual momenta and virtual displacements are differentiated with respect to time. Based on the duality that is observed to exist between the mixed form of Hamilton's weak principle and variational principles governing classical optimal control problems, a temporal finite element formulation of the latter can be developed in a rather straightforward manner. Several well-known problems in dynamics and optimal control are illustrated. The example dynamics problem involves a time-marching problem. As optimal control examples, elementary trajectory optimization problems are treated.

Hodges, Dewey H.↗
NASA NTRS1991

Dynamic mechanical properties of N-phenylnadimide modified PMR polyimide composites

Temperature-frequency dependence of alpha, beta, and gamma transitions was determined using a Rheometrics dynamic spectrometer on a series of unidirectional Celion 6000/N-phenylnadimide (PN) modified PMR polyimide composites. The objective was to see if any correlations exist between crosslinked network structure and dynamic mechanical properties. Variation in crosslinked network structures was achieved by altering the polyimide formulation through addition of various quantities of PN into the standard PMR-15 composition. As a control, PMR-15 composite system exhibited well-defined alpha, beta, and gamma transitions in the regions of 360, 100, and -120 C, respectively. Their activation energies were estimated to be 232, 60, and 14 kcal/mole, respectively. Increasing the amount of PN concentration caused lowering of the activation energies of the three relaxations, a decrease of the glass transition temperature, and increasing intensities of the three damping peaks, compared to the control PMR-15 counterpart. These dynamic mechanical responses were in agreement with formation of a more flexible copolymer from PN and PMR-15 prepolymer.

Pater, Ruth H.↗
NASA NTRS1994

Analysis and design of numerical schemes for gas dynamics 1: Artificial diffusion, upwind biasing, limiters and their effect on accuracy and multigrid convergence

The theory of non-oscillatory scalar schemes is developed in this paper in terms of the local extremum diminishing (LED) principle that maxima should not increase and minima should not decrease. This principle can be used for multi-dimensional problems on both structured and unstructured meshes, while it is equivalent to the total variation diminishing (TVD) principle for one-dimensional problems. A new formulation of symmetric limited positive (SLIP) schemes is presented, which can be generalized to produce schemes with arbitrary high order of accuracy in regions where the solution contains no extrema, and which can also be implemented on multi-dimensional unstructured meshes. Systems of equations lead to waves traveling with distinct speeds and possibly in opposite directions. Alternative treatments using characteristic splitting and scalar diffusive fluxes are examined, together with modification of the scalar diffusion through the addition of pressure differences to the momentum equations to produce full upwinding in supersonic flow. This convective upwind and split pressure (CUSP) scheme exhibits very rapid convergence in multigrid calculations of transonic flow, and provides excellent shock resolution at very high Mach numbers.

Jameson, Antony↗
NASA NTRS1999

Development of a Three-Dimensional PSE Code for Compressible Flows: Stability of Three-Dimensional Compressible Boundary Layers

A program is developed to investigate the linear stability of three-dimensional compressible boundary layer flows over bodies of revolutions. The problem is formulated as a two dimensional (2D) eigenvalue problem incorporating the meanflow variations in the normal and azimuthal directions. Normal mode solutions are sought in the whole plane rather than in a line normal to the wall as is done in the classical one dimensional (1D) stability theory. The stability characteristics of a supersonic boundary layer over a sharp cone with 50 half-angle at 2 degrees angle of attack is investigated. The 1D eigenvalue computations showed that the most amplified disturbances occur around x(sub 2) = 90 degrees and the azimuthal mode number for the most amplified disturbances range between m = -30 to -40. The frequencies of the most amplified waves are smaller in the middle region where the crossflow dominates the instability than the most amplified frequencies near the windward and leeward planes. The 2D eigenvalue computations showed that due to the variations in the azimuthal direction, the eigenmodes are clustered into isolated confined regions. For some eigenvalues, the eigenfunctions are clustered in two regions. Due to the nonparallel effect in the azimuthal direction, the eigenmodes are clustered into isolated confined regions. For some eigenvalues, the eigenfunctions are clustered in two regions. Due to the nonparallel effect in the azimuthal direction, the most amplified disturbances are shifted to 120 degrees compared to 90 degrees for the parallel theory. It is also observed that the nonparallel amplification rates are smaller than that is obtained from the parallel theory.

Balakumar, P.↗