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 91 records · Page 5

State-constrained booster trajectory solutions via finite elements and shooting

This paper presents an extension of a FEM formulation based on variational principles. A general formulation for handling internal boundary conditions and discontinuities in the state equations is presented, and the general formulation is modified for optimal control problems subject to state-variable inequality constraints. Solutions which only touch the state constraint and solutions which have a boundary arc of finite length are considered. Suitable shape and test functions are chosen for a FEM discretization. All element quadrature (equivalent to one-point Gaussian quadrature over each element) may be done in closed form. The final form of the algebraic equations is then derived. A simple state-constrained problem is solved. Then, for a practical application of the use of the FEM formulation, a launch vehicle subject to a dynamic pressure constraint (a first-order state inequality constraint) is solved. The results presented for the launch-vehicle trajectory have some interesting features, including a touch-point solution.

Bless, Robert R.↗

Lightweight porous plastic plaque

The porosity and platability of various materials were investigated to determine a suitable substrate for nickel-plated electrodes. Immersion, ultrasonics, and flow-through plating techniques were tried using nonproprietary formulations, and proprietary phosphide and boride baths. Modifications to the selected material include variations in formulation and treatment, carbon loading to increase conductivity, and the incorporation of a grid. Problems to be solved relate to determining conductivities and porosities as a function of amount of nickel plated on the plastics; loading; charge and discharge curves of electrodes at different current densities; cell performance; and long-term degradation of electrodes.

Reid, M.↗

Computation of solar perturbations with Poisson series

Description of a project for computing first-order perturbations of natural or artificial satellites by integrating the equations of motion on a computer with automatic Poisson series expansions. A basic feature of the method of solution is that the classical variation-of-parameters formulation is used rather than rectangular coordinates. However, the variation-of-parameters formulation uses the three rectangular components of the disturbing force rather than the classical disturbing function, so that there is no problem in expanding the disturbing function in series. Another characteristic of the variation-of-parameters formulation employed is that six rather unusual variables are used in order to avoid singularities at the zero eccentricity and zero (or 90 deg) inclination. The integration process starts by assuming that all the orbit elements present on the right-hand sides of the equations of motion are constants. These right-hand sides are then simple Poisson series which can be obtained with the use of the Bessel expansions of the two-body problem in conjunction with certain interation methods. These Poisson series can then be integrated term by term, and a first-order solution is obtained.

Broucke, R.↗

A variational justification of the assumed natural strain formulation of finite elements

The objective is to study the assumed natural strain (ANS) formulation of finite elements from a variational standpoint. The study is based on two hybrid extensions of the Reissner-type functional that uses strains and displacements as independent fields. One of the forms is a genuine variational principle that contains an independent boundary traction field, whereas the other one represents a restricted variational principle. Two procedures for element level elimination of the strain field are discussed, and one of them is shown to be equivalent to the inclusion of incompatible displacement modes. Also, the 4-node C(exp 0) plate bending quadrilateral element is used to illustrate applications of this theory.

Militello, Carmelo↗

Formulation of cosmic-ray solar daily variation and its seasonal variation, produced from generalized stationary anisotropy of solar origin

In previous papers, a formulation was presented of cosmic ray daily variations produced from solar anisotropies stationary through a year, and also of their annual (or seasonal) modulation caused by the annual variation of the rotation axis of the Earth relative to that of the Sun. These anisotropies are symmetric for an arbitrary rotation around an axis. From observations of the tri-diurnal variation, it has been suggested that solar anisotropies also contain some axis-asymmetric term of the third order with respect to the IMF-axis. This suggestion has recently found support in a theoretical study by Munakata and Nagashima. According to their results, the terms of axis-asymmetry with respect to IMF-axis appear also in the 2nd order anisotropy, together with some different kinds of axis-symmetric terms. The contribution of these anisotropies to the daily variation is different from that of those discussed previously. The above mentioned formulation is extended to a case of a generalized anisotropy.

Tatsuoka, R.↗

Approximating the Helium Wavefunction in Positronium-Helium Scattering

In the Kohn variational treatment of the positronium- hydrogen scattering problem the scattering wave function is approximated by an expansion in some appropriate basis set, but the target and projectile wave functions are known exactly. In the positronium-helium case, however, a difficulty immediately arises in that the wave function of the helium target atom is not known exactly, and there are several ways to deal with the associated eigenvalue in formulating the variational scattering equations to be solved. In this work we will use the Kohn variational principle in the static exchange approximation to d e t e e the zero-energy scattering length for the Ps-He system, using a suite of approximate target functions. The results we obtain will be compared with each other and with corresponding values found by other approximation techniques.

DiRienzi, Joseph↗

Variational Asymptotic Homogenization of Finitely Deformed Viscoelastic-Viscoplastic Composites

The objective of this paper is to develop a constitutive model for finitely deformed viscoelastic-viscoplastic materials and a micromechanics approach to homogenizing composites consisting of such materials. The development of the constitutive model involves establishing a thermodynamic framework based on finite strain theory, developing a viscoelasticity and a viscoplasticity model based on the thermodynamic framework, developing a radial return algorithm based on a classic framework, and deriving a closed-from incremental constitutive relation in the spatial configuration. The development of the micromechanics approach involves pulling-back the above constitutive relation to the material configuration, formulating a variational statement with the resulting constitutive relation, discretizing variational statement in a finite-dimensional space, and solving the discretized variational statement using an Euler–Newton predictor–corrector method. The constitutive model is calibrated via monotonic uniaxial compression tests on a polymer, and the calibrated model is validated by comparing its predictions with the cyclic test data. It is shown capable of characterizing viscoelasticity, viscoplasticity, and complex loading paths. The micromechanics approach’s capabilities are demonstrated through homogenizing a unidirectional fiber-reinforced composite, subjected to uniaxial, biaxial, and shear loading, at different strain rates. It is demonstrated to be capable of handling rate dependence and complex loading paths. The present framework can be further improved by implementing more sophisticated viscoelasticity and viscoplasticity models in future work.

Finite element analysis↗

Halley comet rendezvous with a SEPS vehicle

An analysis of the performance of a Solar Electric Propulsion System (SEPS) vehicle rendezvousing with Halley's comet just prior to its Frebruary 1986 perihelion is described. A calculus of variations mathematical formulation is used to maximize Halley arrival mass while giving effect to the influence of solar array size, launch date, arrival date, and insertion hyperbolic excess velocity. Numerical sensitivity relief, thrust system modeling, trajectory characteristics and ion engine operating conditions are discussed and illustrated. Results indicate a rendezvous is feasible with a minimal advance in solar cell and ion engine technology.

Burrows, R. R.↗

A survey of the theory of the Earth's rotation

The theory of the Earth's rotation and the geophysical phenomena affecting it is examined. First principles are reviewed and the problem of polar motion and UT1 variations is formulated in considerable generality and detail. The effects of Earth deformations and the solid Earth tides are analyzed.

Cannon, W. H.↗

Comparison of Euler and full potential marching techniques for flows over complex configurations

Two recently developed aerodynamic prediction techniques based on the steady full potential equation and the unsteady Euler equations have been applied to a variety of three-dimensional supersonic flow problems exhibiting embedded subsonic regions. Both techniques utilize planar Gauss-Seidel relaxation in the marching direction and approximate factorization in the cross-flow plane. A conservative switching scheme and flux bias technique are employed in the full potential method to transition from the supersonic marching procedure to a subsonic relaxation algorithm and vice versa. A new unified approach with finite volume, high accuracy (up to third order) Total Variation Diminishing formulation (based on Roe's scheme) is used in the Euler solver. In the supersonic regions of the flow an 'infinitely large' time step is employed, and a finite time step is applied in the subsonic regions of the flow to reach the steady-state as a time-asymptote. Numerical solutions are obtained for a number of complex configurations, including: (1) an elliptic waverider, (2) a realistic fighter configuration, (3) the Space Shuttle, and (4) a Shuttle-like configuration. Both the Full Potential and Euler numerical results are in good agreement with available experimental data.

Szema, K. Y.↗

Euler solver for three-dimensional supersonic flows with subsonic pockets

A new finite-difference scheme has been developed to solve efficiently the unsteady Euler equations for three-dimensional inviscid supersonic flows with subsonic pockets. The technique utilizes planar Gauss-Seidel relaxation in the marching direction and approximate factorization n the crossflow plane. An 'infinitely large' time step is used in parts of the flowfield where the component of velocity in the marching direction is supersonic - here the Gauss-Seidel sweeps are restricted to the forward direction only, and the procedure reduces to simple space-marching; a finite time step is used in parts of the flowfield where the marching component of velocity is subsonic - here, backward and forward Gauss-Seidel sweeps are employed to allow for upstream and downstream propagation of signals, and a time-asymptotic steady state is obtained. The discretization formulas are based on finite-volume implementation of high accuracy (up to third-order) total variation diminishing formulations. Numerical solutions are obtained for an analytically defined forebody, a realistic fighter configuration, and the Space Shuttle. The results are in very good agreement with available experimental data and numerical solutions of the full-potential equation.

Chakravarthy, Sukumar R.↗

Variational finite element-tensor formulation for the large deflection random vibration of composite plates

A numerical integration routine is derived from a set of unified single step integration algorithms using a weighted satisfaction of the equilibrium equation governing the large deflection random response of laminated composite plates. The equilibrium equation is derived using a constant matrix large deflection finite-element formulation. In-plane inertia terms are considered in the formulation, however, rotary inertia terms are assumed negligible. Probability density, spectral density and autocorrelation functions of the maximum displacement and strain responses are presented for three acoustic excitation levels. Classical thin plate boundary conditions and pseudo white noise excitation are used in this investigation.

Robinson, Jay H.↗

Comparison of Models for Ball Bearing Dynamic Capacity and Life

Generalized formulations for dynamic capacity and life of ball bearings, based on the models introduced by Lundberg and Palmgren and Zaretsky, have been developed and implemented in the bearing dynamics computer code, ADORE. Unlike the original Lundberg-Palmgren dynamic capacity equation, where the elastic properties are part of the life constant, the generalized formulations permit variation of elastic properties of the interacting materials. The newly updated Lundberg-Palmgren model allows prediction of life as a function of elastic properties. For elastic properties similar to those of AISI 52100 bearing steel, both the original and updated Lundberg-Palmgren models provide identical results. A comparison between the Lundberg-Palmgren and the Zaretsky models shows that at relatively light loads the Zaretsky model predicts a much higher life than the Lundberg-Palmgren model. As the load increases, the Zaretsky model provides a much faster drop off in life. This is because the Zaretsky model is much more sensitive to load than the Lundberg-Palmgren model. The generalized implementation where all model parameters can be varied provides an effective tool for future model validation and enhancement in bearing life prediction capabilities.

Rolling-element Fatigue↗

Nonlinear Homogenization of Finitely Deformed Viscoelastic-Viscoplastic Composites Using Mechanics of Structure Genome

The objective of this paper is to develop a micromechanics approach to homogenizing finitely deformed viscoelastic-viscoplastic composites using the mechanics of structure genome. The incremental constitutive relation for glassy polymers, formulated in the spatial configuration, is implemented in the present approach.This involves (1) pulling-back the constitutive model to the material configuration and (2)choosing the deformation gradient tensor and the first Piola–Kirchhoff stress tensor as the strain and the stress measures during homogenization, respectively. An Euler–Newton predictor–corrector method is developed for homogenization. Each step involves formulating a variational statement using the mechanics of structure genome, discretizing the statement in a finite-dimensional space, and solving the problem using an Euler/multilevel Newton method. The present approach is demonstrated by homogenizing fiber- and particle-reinforced composites undergoing uniaxial, biaxial, or shear deformation, at different stain rates.

Multi-scale modeling, High Strain Composites, Visc↗

Computing Sensitivities in Evolutionary Systems: A Real-time Reduced Order Modeling Strategy

We present a new methodology for computing sensitivities in evolutionary systems using a model-driven low-rank approximation. To this end, we formulate a variational principle that seeks to minimize the distance between the time derivative of the reduced approximation and sensitivity dynamics. The first order optimality condition of the variational principle leads to a system of closed form evolution equations for an orthonormal basis and corresponding sensitivity coefficients. This approach allows for the computation of sensitivities with respect to a large number of parameters in an accurate and tractable manner by extracting correlations between different sensitivities on the fly. The presented method requires solving forward evolution equations, sidestepping the restrictions imposed by the forward/backward workflow of adjoint sensitivities. For example, the presented method, unlike the adjoint equation, does not impose any input/output load and can be used in applications in which real-time sensitivities are of interest. We demonstrate the utility of the method for three test cases: (1) computing sensitivity with respect to model parameters in the Rössler system, (2) computing sensitivity with respect to an infinite-dimensional forcing parameter in the chaotic Kuramoto--Sivashinsky equation, and (3) computing sensitivity with respect to reaction parameters for species transport in a turbulent reacting flow.

Reduced order model↗

Analysis of thermal stresses and metal movement during welding. I - Analytical study

This is the first part of a study of thermal stresses and metal movement during welding. This part discusses analysis of two-dimensional thermal stresses and metal movement during bead-on-plate and butt welding. A finite-element formulation has been derived, based on the variational principle. The formulation includes temperature dependence of material properties as well as the yield criterion.

Muraki, T.↗