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.

SourceResearch connectionFrom yearThrough year

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

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS1993

A nonreflecting upper boundary condition for anelastic nonhydrostatic mesoscale gravity-wave models

A sponge layer is formulated to prevent spurious reflection of vertically propagating quasi-stationary gravity waves at the upper boundary of a two-dimensional numerical anelastic nonhydrostatic model. The sponge layer includes damping of both Newtonian-cooling type and Rayleigh-friction type, whose coefficients are determined in such a way that the reflectivity of wave energy at the bottom of the layer is zero. Unlike the formulations in earlier studies, our formulation includes the effects of vertical discretization, vertical mean density variation, and nonhydrostaticity. This sponge formulation is found effective in suppressing false downward reflection of waves for various types of quasi-stationary forcing.

Kim, Young-Joon↗

NASA NTRS1992

Improved assumed-stress hybrid shell element with drilling degrees of freedom for linear stress, buckling, and free vibration analyses

An improved four-node quadrilateral assumed-stress hybrid shell element with drilling degrees of freedom is presented. The formulation is based on Hellinger-Reissner variational principle and the shape functions are formulated directly for the four-node element. The element has 12 membrane degrees of freedom and 12 bending degrees of freedom. It has nine independent stress parameters to describe the membrane stress resultant field and 13 independent stress parameters to describe the moment and transverse shear stress resultant field. The formulation encompasses linear stress, linear buckling, and linear free vibration problems. The element is validated with standard tests cases and is shown to be robust. Numerical results are presented for linear stress, buckling, and free vibration analyses.

Rengarajan, Govind↗

NASA NTRS1974

Calculation of precision satellite orbits with nonsingular elements /VOP formulation/

Review of some results obtained in an effort to develop efficient, high-precision trajectory computation processes for artificial satellites by optimum selection of the form of the equations of motion of the satellite and the numerical integration method. In particular, the matching of a Gaussian variation-of-parameter (VOP) formulation is considered which is expressed in terms of equinoctial orbital elements and partially decouples the motion of the orbital frame from motion within the orbital frame. The performance of the resulting orbit generators is then compared with the popular classical Cowell/Gauss-Jackson formulation/integrator pair for two distinctly different orbit types - namely, the orbit of the ATS satellite at near-geosynchronous conditions and the near-circular orbit of the GEOS-C satellite at 1000 km.

Velez, C. E.↗

NASA NTRS1992

A variational method for finite element stress recovery: Applications in one-dimension

It is well-known that stresses (and strains) calculated by a displacement-based finite element analysis are generally not as accurate as the displacements. In addition, the calculated stress field is typically discontinuous at element interfaces. Because the stresses are typically of more interest than the displacements, several procedures have been proposed to obtain a smooth stress field, given the finite element stresses, and to improve the accuracy. Hinton and Irons introduced global least squares smoothing of discrete data defined on a plane using a finite element formulation. Tessler and co-workers recently developed a conceptually similar formulation for smoothing of two-dimensional data based on a discrete least square approximation with a penalty constraint. The penalty constraint results in a stress field which is C(exp 1)-continuous, a result not previously obtained. The approach requires additional, 'smoothing' finite element analysis and for their two-dimensional application, they used a conforming C(exp 0)-continuous triangular finite element based on a conforming plate element. This paper presents the results of a detailed investigation into the application of Tessler's smoothing procedure to the smoothing of finite element stresses from one-dimensional problems. Although the one-dimensional formulation has some practical applicability, such as in truss, beam, axisymmetric mechanics, and one-dimensional heat conduction, the primary motivation for developing the one-dimensional smoothing case is to explore the characteristics of the general smoothing strategy. In particular, it is used to describe the behavior of the method and to explore the suitability of criteria proposed for the smoothing analysis. Prior to presenting numerical results, the variational formulation of the smoothing strategy is presented and a criterion for the smoothing analysis is described.

Riggs, H. Ronald↗

NASA NTRS1994

Advanced data assimilation in strongly nonlinear dynamical systems

Advanced data assimilation methods are applied to simple but highly nonlinear problems. The dynamical systems studied here are the stochastically forced double well and the Lorenz model. In both systems, linear approximation of the dynamics about the critical points near which regime transitions occur is not always sufficient to track their occurrence or nonoccurrence. Straightforward application of the extended Kalman filter yields mixed results. The ability of the extended Kalman filter to track transitions of the double-well system from one stable critical point to the other depends on the frequency and accuracy of the observations relative to the mean-square amplitude of the stochastic forcing. The ability of the filter to track the chaotic trajectories of the Lorenz model is limited to short times, as is the ability of strong-constraint variational methods. Examples are given to illustrate the difficulties involved, and qualitative explanations for these difficulties are provided. Three generalizations of the extended Kalman filter are described. The first is based on inspection of the innovation sequence, that is, the successive differences between observations and forecasts; it works very well for the double-well problem. The second, an extension to fourth-order moments, yields excellent results for the Lorenz model but will be unwieldy when applied to models with high-dimensional state spaces. A third, more practical method--based on an empirical statistical model derived from a Monte Carlo simulation--is formulated, and shown to work very well. Weak-constraint methods can be made to perform satisfactorily in the context of these simple models, but such methods do not seem to generalize easily to practical models of the atmosphere and ocean. In particular, it is shown that the equations derived in the weak variational formulation are difficult to solve conveniently for large systems.

Miller, Robert N.↗

NASA NTRS1987

Equation-of-state spinning fluids in the Einstein-Cartan theory

The relativistic fluid equations may be completed in two physically distinct methods. One method assumes the mass, rho, (or particle number) is conserved, while the other method assumes an equation of state of the form P = P(rho). A variational principle for the mass conservation method both with and without an intrinsic spin for the fluid was constructed earlier (Ray and Smalley, 1982 and 1983). A variational principle for the fluid described by an equation of state both with and without spin is formulated. In all cases the variational principle is set in the Einstein-Cartan metric-torsion U4 geometry. The results for general relativity follow as a special case.

Ray, John R.↗

NASA NTRS1971

A Stability Analysis of Cylindrical Panels Using a Finite Element Formulation

A cylindrical finite element suitable for the linear stability analysis of cylindrical shells is developed. Energy principles and variational methods lead to a problem formulation which lends itself to physical interpretations of the governing matrices of the finite element. By properly grouping the terms which result from taking the second variation of the potential energy of the element, it is possible to identify three distinct types of matrices. The first matrix is the conventional stiffness matrix; the second is an initial stress stiffness matrix; and the third is an initial displacement stiffness matrix. With the assumption of linearity, the buckling problem is stated in terms of the classical linear real eigenvalue equation. This problem formulation was programmed on the CDC 6600 series computer. The computer program is used to analyze the buckling of a variety of structures. Columns, arches, flat plates and curved panels with and without cutouts are considered.

Snyder, R. E.↗

NASA NTRS1972

Void growth in an elastic-plastic medium.

The uniaxial deformation of an elastic-plastic medium containing a doubly periodic square array of circular cylindrical voids is studied under plane-strain conditions. Both the effects of geometrical nonlinearities resulting from large deformation and physical nonlinearities arising from plastic material behavior are included in formulating the problem. A variational principle is used as the basis for implementing a finite-element solution. Results are obtained for the change in void shape and size under increasing overall strain, the overall tensile behavior of the material with voids, and the development of the plastic zone about a void.

Needleman, A.↗

NASA NTRS1974

Vehicle control for fuel optimization

The problem of predicting the minimum fuel trajectory for a six degree of freedom vehicle which has a motion characterized by the first order differential equations of translational and rotational dynamics is considered. The thrust direction and center of gravity of the vehicle are assumed to be fixed with respect to the vehicle. Thrust magnitude and the control moment are used as control variables and appear linearly in the equation of motion. Pontryagin's Maximum Principle is used to solve the variational problem. With this formulation, the extremal controls are bang-bang with the exception of the singular case.

Dannenberg, K. D.↗

NASA NTRS1974

A dual extremum principle for a population equation

A dual extremum principle for the Verhulst-Pearl population equation is constructed using a complementary variational technique. The dual formulation utilizes a minimum principle recently developed by Leitmann to convert the functional optimization problem into a parameter optimization problem.

Chan, W. L.↗

NASA NTRS1976

A new formulation for anisotropic radiative transfer problems. I - Solution with a variational technique

The equations of radiative transfer in anisotropically scattering media are reformulated as linear operator equations in a single independent variable. The resulting equations are suitable for solution by a variety of standard mathematical techniques. The operators appearing in the resulting equations are in general nonsymmetric; however, it is shown that every bounded linear operator equation can be embedded in a symmetric linear operator equation and a variational solution can be obtained in a straightforward way. For purposes of demonstration, a Rayleigh-Ritz variational method is applied to three problems involving simple phase functions. It is to be noted that the variational technique demonstrated is of general applicability and permits simple solutions for a wide range of otherwise difficult mathematical problems in physics.

Cheyney, H., III↗

NASA NTRS1976

Optimal perturbation models for averaged orbit generation

Averaging techniques applied to the variation of parameters (VOP) formulation of the equations of motion are being investigated as methods for long-term prediction of artificial satellite orbits. Analytically averaged equations were compared with numerically averaged equations with respect to accuracy and efficiency for computation of zonal and nonresonant third-body perturbations. Numerically averaged equations were also evaluated for computation of long-period effects from resonant third-body, tesseral harmonic, and atmospheric drag perturbations. Guidelines will be presented for application of averaged VOP equations to a broad class of orbits.

Long, A. C.↗

NASA NTRS1979

Antenna feedhorn software upgrade

The HYBRIDHORN computer program was developed to serve as an item of general purpose antenna feedhorn design and analysis software. The formulation contains a small flare angle approximation which is subject to question for designs such as the S- and X-band feedhorn. Additionally, the original formulation did not allow azimuthal variation indexes other than unity. The HYBRIDHORN program was upgraded to correct both of these deficiencies. A large flare angle formulation was found. In the upgrade, all of the major program elements were converted to Univac 1108 compatible structured FORTRAN (SFTRAN) for ease of software maintenance. The small and large angle formulations are described and sample numerical results are presented.

Potter, P. D.↗

NASA NTRS1981

On Hybrid and mixed finite element methods

Three versions of the assumed stress hybrid model in finite element methods and the corresponding variational principles for the formulation are presented. Examples of rank deficiency for stiffness matrices by the hybrid stress model are given and their corresponding kinematic deformation modes are identified. A discussion of the derivation of general semi-Loof elements for plates and shells by the hybrid stress method is given. It is shown that the equilibrium model by Fraeijs de Veubeke can be derived by the approach of the hybrid stress model as a special case of semi-Loof elements.

Pian, T. H. H.↗

NASA NTRS1986

A higher order theory of laminated composite cylindrical shells

A new higher order theory has been proposed for the analysis of composite cylindrical shells. The formulation allows for arbitrary variation of inplane displacements. Governing equations are presented in the form of a hierarchy of sets of partial differential equations. Each set describes the shell behavior to a certain degree of approximation. The natural frequencies of simply-supported isotropic and laminated shells and stresses in a ring loaded composite shell have been determined to various orders of approximation and compared with three dimensional solutions. These numerical studies indicate the improvements achievable in estimating the natural frequencies and the interlaminar shear stresses in laminated composite cylinders.

Krishna Murthy, A. V.↗

NASA NTRS1993

Transverse cracking and stiffness reduction in composite laminates

A study of transverse cracking mechanism in composite laminates is presented using a singular hybrid finite element model. The model provides the global structural response as well as the precise local crack-tip stress fields. An elasticity basis for the problem is established by employing Lekhnitskii's complex variable potentials and method of eigenfunction expansion. Stress singularities associated with the transverse crack are obtained by decomposing the deformation into the symmetric and antisymmetric modes and proper boundary conditions. A singular hybrid element is thereby formulated based on the variational principle of a modified hybrid functional to incorporate local crack singularities. Axial stiffness reduction due to transverse cracking is studied. The results are shown to be in very good agreement with the existing experimental data. Comparison with simple shear lag analysis is also given. The effects of stress intensity factors and strain energy density on the increase of crack density are analyzed. The results reveal that the parameters approach definite limits when crack densities are saturated, an evidence of the existence of characteristic damage state.

Yuan, F. G.↗

NASA NTRS1993

Hybrid Analytical Technique for Nonlinear Vibration Analysis of Thin-Walled Beams

A two-step hybrid analytical technique is presented for the nonlinear vibration analysis of thin-walled beams. The first step involves the generation of various-order perturbation functions using the Linstedt-Poincare perturbation technique. The second step consists of using the perturbation functions as coordinate (or approximation) functions and then computing both the amplitudes of these functions and the nonlinear frequency of vibration via a direct variational procedure. The analytical formulation is based on a form of the geometrically nonlinear beam theory with the effects of in-plane inertia, rotatory inertia, and transverse shear deformation included. The effectiveness of the proposed technique is demonstrated by means of a numerical example of thin-walled beam with a doubly symmetric I-section. The solutions obtained using a single-spatial mode were compared with those obtained using multiple-spatial modes. The standard of comparison was taken to be the frequencies obtained by the direct integration/fast Fourier transform (FFT) technique. The nonlinear frequencies obtained by the hybrid technique were shown to converge to the corresponding ones obtained by the direct integration/fast Fourier transform (FFT) technique well beyond the range of applicability of the perturbation technique. The frequencies and total strain energy of the beam were overestimated by using a single-spatial mode.

Noor, Ahmed K.↗

NASA NTRS2009

Stability of Dosage Forms in the Pharmaceutical Payload Aboard Space Missions

Efficacious pharmaceuticals with adequate shelf lives are essential for successful space medical operations. Stability of pharmaceuticals, therefore, is of paramount importance for assuring the health and wellness of astronauts on future space exploration missions. Unique physical and environmental factors of space missions may contribute to the instability of pharmaceuticals, e.g., radiation, humidity and temperature variations. Degradation of pharmaceutical formulations can result in inadequate efficacy and/or untoward toxic effects, which could compromise astronaut safety and health. Methods: Four identical pharmaceutical payload kits containing 31 medications in different dosage forms (liquid, tablet, capsule, ointment and suppository) were transported to the International Space Station aboard the Space Shuttle (STS-121). One of the 4 kits was stored on the Shuttle and the other 3 were stored on the International Space Station (ISS) for return to Earth at 6-month interval aboard a pre-designated Shuttle flight for each kit. The kit stored on the Shuttle was returned to Earth aboard STS-121 and 2 kits from ISS were returned on STS 117 and STS-122. Results: Analysis of standard physical and chemical parameters of degradation was completed for pharmaceuticals returned by STS-121 after14 days, STS - 117 after11 months and STS 122 after 19 months storage aboard ISS. Analysis of all flight samples along with ground-based matching controls was completed and results were compiled. Conclusion: Evaluation of results from the shuttle (1) and ISS increments (2) indicate that the number of formulations degraded in space increased with duration of storage in space and was higher in space compared to their ground-based counterparts. Rate of degradation for some of the formulations tested was faster in space than on Earth. Additionally, some of the formulations included in the medical kits were unstable, more so in space than on the ground. These results indicate that the space flight environment may adversely affect the shelf life of pharmaceuticals aboard space missions.

Du, Brian J.↗