Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational methods”

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 271 records · Page 15

Thermally-induced, geometrically nonlinear response of symmetrically laminated composite plates

This paper discusses the thermally-induced geometrically nonlinear response of symmetrically laminated composite plates. The plate response is due to a temperature increase that is uniform in the plane of the plate but has a slight gradient through the thickness. The case of a completely uniform temperature increase but with an initial out-of-plane imperfection in the plate is also considered. Because they are closely allied problems, thermal buckling and postbuckling are discussed. Using variational methods in conjunction with a Rayleigh-Ritz formulation, these responses are investigated for two laminates, a (+/- 45/0/90)s and a (+/- 45/02)s, under two different simple support conditions.

Meyers, C. A.↗

DUKSUP - A high thrust trajectory optimization code

Designing missions on expendable launch vehicles (ELV's) includes determining launch vehicle performance capabilities and trajectory characteristics over the range of mission requirements for suitable launch periods. This analysis depends on mathematically modeling both the launch vehicle and the mission requirements. Generally the result is a mathematical model that is described by an objective function to be optimized subject to an assortment of algebraic and dynamic constraints. About 30 years ago, engineers at Lewis Research Center (LeRC) undertook the task of creating a software code to solve 3D versions of such problems and based it on the Calculus of Variations/Optimal Control Theory. One of these codes, DUKSUP, has been in use at LeRC for nearly 25 years, during which time it has played an important role in a large number of studies and actual missions. Currently it is being used by about 12 analysts in the Center's Advanced Space Analysis Office (ASAO) to do mission design, feasibility studies, corroboration of contractor data and planning studies for the Space Exploration Initiative (SEI). Today, it is one of the few ELV mission analysis production codes based on variational methods. With future ELV missions in mind, ASAO is presently creating a new code to upgrade DUKSUP's capabilities.

Balkanyi, Leslie R.↗

Thermal buckling of symmetrically laminated composite plates

This paper discusses an investigation into the thermal buckling of symmetrically laminated composite plates. In this study thermal buckling is investigated for laminates under two different simple support conditions, fixed and sliding. These laminates are subjected to the conditions of a uniform temperature change, and a temperature change varying linearly along the length of the plate. The effects of the principal material axes not being aligned with the edges of the plate are also investigated. The buckling response is studied using variational methods, specifically the Trefftz criterion. A Rayleigh-Ritz formulation is used to obtain numerical results from the formulations for the prebuckling response and the buckling response.

Meyers, C. A.↗

Ellipsoidal figures of equilibrium - Compressible models

The results of Chandrasekhar (1969) are generalized to polytropes, using a formalism based on ellipsoidal energy variational principle to construct approximate stellar equilibrium solutions and study their stability. After reviewing the energy variational method and describing the approach, several equivalent stability conditions are established and secular vs. dynamical instabilities are discussed. Then, the equilibrium structure equations are derived for isolated, rotating polytropes, and axisymmetric configurations (compressible Maclaurin spheroids) are considered. Particular attention is given to triaxial configurations, either in a state of uniform rotation (generalizing the classical Jacobi ellipsoids) or with internal fluid motions of uniform vorticity (the compressible analogues of Riemann-S ellipsoids) and to the stability of these single star configurations. The compressible generalizations of the Roche and Roche-Riemann problems for a polytrope in orbit about a point-mass companion are solved, and the generalized Darwin problem for two identical polytropes in a binary is considered.

Lai, Dong↗

Electronic and total energy properties of ternary and quaternary semiconductor compounds, alloys and superlattices

This proposal was mainly concerned with the theoretical study of semiconductor compounds, alloys, and superlattices of interest for photovoltaic applications. In the last year (1991) a study was devoted to metal/graphite bonding in relation to use of graphite fiber reinforcement of Cu for high thermal conductivity applications. The main research topics addressed during the full period of the grant are briefly described: studies of the In-Ga-As ternary system; band-offsets at common anion and InAs/GaSb/AlSb heterojunctions; alloy theory (cluster variation method); and Cu/graphite bonding. Most of the work was described more extensively in previous yearly reports and renewal applications and in publications. The last topic is described more fully in a separate report attached. A list of publications resulting directly from this grant or from other grants but related to this work and of conference presentations is given at the end.

Lambrecht, Walter R. L.↗

Mapping ERS-1 wind fields over north west Atlantic using a variational objective analysis

A variational method is implemented to produce five day mean gridded ERS-1 analyzed wind fields in the north west Atlantic with the aim of providing for wind forcing of basin scale ocean models. The method consists of minimizing a cost functional, designed to measure misfits to prescribed weighted constraints which express a smoothed behavior and the proximity to input data vectors and curl. The weights are empirically determined by comparison with independent ship and buoy data over four five day periods. Root mean square differences between analyzed winds and independent data are thus further decreased: they range from 0.8 up to 1.8 m/s. Some problems present in the initial data remain; they are principally due to incomplete data coverage (instrumental problems), and possibly unresolved ambiguities. The resulting curl fields are smoothed and show coherent patterns. A comparison with the European Center for Medium range Weather Forecasting (ECMWF) analysis is encouraging.

Siefridt, L.↗

Plasma-screening effects on the electron-impact excitation of hydrogenic ions in dense plasmas

Plasma-screening effects are investigated on electron-impact excitation of hydrogenic ions in dense plasmas. Scaled cross sections Z(exp 4) sigma for 1s yields 2s and 1s yields 2p are obtained for a Debye-Hueckel model of the screened Coulomb interaction. Ground and excited bound wave functions are modified in the screened Coulomb potential (Debye-Hueckel model) using the Ritz variation method. The resulting atomic wave functions and their eigenenergies agree well with the numerical and high-order perturbation theory calculations for the interesting domain of the Debye length not less than 10. The Born approximation is used to describe the continuum states of the projectile electron. Plasma screening effects on the atomic electrons cannot be neglected in the high-density cases. Including these effects, the cross sections are appreciably increased for 1s yields 2s transitions and decreased for 1s yields 2p transitions.

Jung, Young-Dae↗

An analysis of nearfield normal mode amplitude anomalies of the Landers earthquake

The 1992 Landers earthquake (M(sub w) = 7.3) occurred in the middle of the TERRAscope network. Long-period Rayleigh waves recorded at the TERRAscope stations (delta less than or = 3 deg) after traveling around the Earth show large amplitude anomalies, one order of magnitude larger than spherical Earth predictions up to a period of about 600 s. The ground motions over the epicentral region at and after the arrival of R4-5 are in phase at all stations. These observations are inconsistent with the nearly vertical strike slip mechanism of the Landers earthquake. Synthetic seismograms for a rotating, elliptic, and laterally heterogeneous Earth model calculated by the variational method agree well with the observed waveforms. Calculations for various 3D Earth models demonstrate that the amplitudes are very sensitive to the large scale aspherical structure in the crust and the mantle. The anomalies for modes shorter than 300 s period can be explained by lateral heterogeneity shallower than the upper mantle. Rotation of the Earth and lower mantle heterogeneity are required to explain mode amplitudes at longer periods. Current whole mantle seismic tomographic models can fully explain the observed amplitudes longer than 300 s. To assess the effect of the high order lateral heterogeneity in the mantle more precise estimate of the crustal correction is required.

Watada, Shingo↗

Finite Element Analysis of Geodesically Stiffened Cylindrical Composite Shells Using a Layerwise Theory

Layerwise finite element analyses of geodesically stiffened cylindrical shells are presented. The layerwise laminate theory of Reddy (LWTR) is developed and adapted to circular cylindrical shells. The Ritz variational method is used to develop an analytical approach for studying the buckling of simply supported geodesically stiffened shells with discrete stiffeners. This method utilizes a Lagrange multiplier technique to attach the stiffeners to the shell. The development of the layerwise shells couples a one-dimensional finite element through the thickness with a Navier solution that satisfies the boundary conditions. The buckling results from the Ritz discrete analytical method are compared with smeared buckling results and with NASA Testbed finite element results. The development of layerwise shell and beam finite elements is presented and these elements are used to perform the displacement field, stress, and first-ply failure analyses. The layerwise shell elements are used to model the shell skin and the layerwise beam elements are used to model the stiffeners. This arrangement allows the beam stiffeners to be assembled directly into the global stiffness matrix. A series of analytical studies are made to compare the response of geodesically stiffened shells as a function of loading, shell geometry, shell radii, shell laminate thickness, stiffener height, and geometric nonlinearity. Comparisons of the structural response of geodesically stiffened shells, axial and ring stiffened shells, and unstiffened shells are provided. In addition, interlaminar stress results near the stiffener intersection are presented. First-ply failure analyses for geodesically stiffened shells utilizing the Tsai-Wu failure criterion are presented for a few selected cases.

Gerhard, Craig Steven↗

An MHD variational principle that admits reconnection

The variational approach of Pfirsch and Sudan's averaged magnetohydrodynamics (MHD) to the stability of a line-tied current layer is summarized. The effect of line-tying on current sheets that might arise in line-tied magnetic flux tubes by estimating the growth rates of a resistive instability using a variational method. The results show that this method provides a potentially new technique to gauge the stability of nearly ideal magnetohydrodynamic systems. The primary implication for the stability of solar coronal structures is that tearing modes are probably constant at work removing magnetic shear from the solar corona.

Rilee, M. L.↗

Positive Anharmonicities: The Oxonide Anion as an Example

An accurate ab initio quartic force field for the ozonide anion has been determined at the singles and doubles coupled-cluster level of theory that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), using the augmented valence triple-zeta correlation consistent one-particle basis set. Convergence of the harmonic frequencies with respect to the one-particle basis set has been demonstrated by computing quadratic force fields at the CCSD(T) level using augmented valence double-zeta and quadruple-zeta basis sets. Fundamental vibrational frequencies have been determined via second-order ro-vibrational perturbation theory and also using exact variational methods. Agreement is very good, and they both show that the antisymmetric O-O stretch, v 3, possesses a positive anharmonic correction (that is, the fundamental frequency is predicted to be higher in energy than the harmonic frequency). Comparison of the O_3 and O3 quartic force fields shows that the positive anharmonic correction is the result of a particularly large F3333 symmetry internal coordinate force constant. The reasonableness of this force constant has been tested by computing both F33 and F3333 at the CCSD, CCSD(T), and CASPT2 levels of theory. A discussion of known positive anharmonicities for stretching modes is presented.

Lee, Timothy J.↗

Structural Modeling of a Five-Meter Thin Film Inflatable Antenna/Concentrator

Inflatable structures have been the subject of renewed interest in recent years for space applications such as communications antennas, solar thermal propulsion, and space solar power. A major advantage of using inflatable structures in space is their extremely light weight. An obvious second advantage is on-orbit deployability and related space savings in the launch configuration. A recent technology demonstrator flight for inflatable structures was the Inflatable Antenna Experiment (IAE) that was deployed on orbit from the Shuttle Orbiter. Although difficulty was encountered in the inflation/deployment phase, the flight was successful overall and provided valuable experience in the use of such structures. Several papers on static structural analysis of inflated cylinders have been written, describing different techniques such as linear shell theory, and nonlinear and variational methods, but very little work had been done in dynamics of inflatable structures until recent years. In 1988 Leonard indicated that elastic beam bending modes could be utilized in approximating lower-order frequencies of inflatable beams. Main, et al. wrote a very significant 1995 paper describing results of modal tests of inflated cantilever beams and the determination of effective material properties. Changes in material properties for different pressures were also discussed, and the beam model was used in a more complex structure. The paper demonstrated that conventional finite element analysis packages could be very useful in the analysis of complex inflatable structures. The purposes of this paper are to discuss the methodology for dynamically characterizing a large 5-meter thin film inflatable reflector, and to discuss the test arrangement and results. Nonlinear finite element modal results are compared to modal test data. The work is significant and of considerable interest to researchers because of 1) the large size of the structure, making it useful for scaling studies, and 2) application of commercially available finite element software for modeling pressurized thin-film structures.

Smalley, Kurt B.↗

Snap-Through of Unsymmetric Laminates Using Piezocomposite Actuators

The paper discusses the concept of using a piezoceramic actuator bonded to one side of a two-layer unsymmetric cross-ply [0/90]T laminate to provide the moments necessary to snap the laminate from one stable equilibrium shape to another. This concept could be applied to the morphing of structures. A model of this concept, which is based on the Rayleigh-Ritz technique and the use of energy and variational methods, is developed. The experimental phase of the study is discussed, including the measurement of the voltage level needed to snap the laminate. The voltage measurements and shapes are compared with predictions of the models and the agreement between measurements and the predictions are reasonable, both qualitatively and quantitatively. Suggestions for future activities are presented.

Schultz, Marc R.↗

Electron-Molecule Col1isions: Quantitative Approaches, and the Legacy of Aaron Temkin

This article, on electron-molecule collisions, is dedicated to the legacy of my good friend and sometime collaborator, Aaron Temkin on his retirement from the NASA-Goddard Space Flight Center after many years of work at the highest intellectual level in the theoretical treatment of electron-atom and electron-molecule scattering. Aaron's contributions to the manner in which we think about electron-molecule collisions is clear to all of us who have worked in this field. I doubt that the great progress that has occurred in the computational treatment of such complex collision problems could have happened without these contributions. For a brief historical account, see the discussion of Temkin's contribution to electron-molecule scattering in the first article of this volume by Dr. A. K. Bhatia. In this article, I will concentrate on the application of the so called, non-adiabatic R-matrix theory, to vibrational excitation and dissociative attachment, although I will also present some results applying the Linear Algebraic and Kohn-Variational methods to vibrational excitation. As a starting point for almost all computationally effective approaches to electron-molecule collisions, is the fixed nuclei approximation. That is, one recognizes, just as one does with molecular bound states, that there is a separation of electronic(fast) and nuclear(s1ow) degrees of freedom. This separation makes it possible to "freeze" the nuclei in space, calculate the collision parameters for the frozen molecule and then, somehow to add back the vibrations and rotations. The manner in which this is done, depends on the details of the collision problem. It is the work of Aaron and a number of other researchers that has provided the guidance necessary to resolve these issues.

Schneider, B.I.↗

Testing the Reliability of Cluster Mass Indicators with a Systematics Limited Dataset

We present the mass X-ray observable scaling relationships for clusters of galaxies using the XMM-Newton cluster catalog of Snowden et al. Our results are roughly consistent with previous observational and theoretical work, with one major exception. We find 2-3 times the scatter around the best fit mass scaling relationships as expected from cluster simulations or seen in other observational studies. We suggest that this is a consequence of using hydrostatic mass, as opposed to virial mass, and is due to the explicit dependence of the hydrostatic mass on the gradients of the temperature and gas density profiles. We find a larger range of slope in the cluster temperature profiles at radii 500 than previous observational studies. Additionally, we find only a weak dependence of the gas mass fraction on cluster mass, consistent with a constant. Our average gas mass fraction results also argue for a closer study of the systematic errors due to instrumental calibration and modeling method variations between analyses. We suggest that a more careful study of the differences between various observational results and with cluster simulations is needed to understand sources of bias and scatter in cosmological studies of galaxy clusters.

Juett, Adrienne M.↗

Optical-Quality Thin Polymer Membranes

A method of fabricating both curved and flat thin polymer membranes of optical quality has been developed. The method was originally intended to enable the fabrication of lightweight membrane imaging and interferometric optics, possibly with apertures multiple meters wide, for use in scientific instruments that would operate in outer space. The method may also be applicable to the fabrication of lightweight membrane optics for terrestrial use. The method involves flow-casting of a soluble polymer with mechanical and environmental controls that provide nearly ideal conditions for the formation of a membrane. The preferred environmental conditions and other details of the process depend on the choice of polymer and substrate material and on the shape and size of the membrane to be cast. Once the polymer has dried to a membrane, it is cured with convective heating, then released. Membranes with root-mean-square surface roughnesses of less than 10.5 can be produced routinely by this method. Variations in the thicknesses of the membranes have ranged from 1/3 wavelength down to as little as 1/20 wavelength (at a wavelength of 633 nm). Membranes fabricated thus far have had diameters up to 0.5 m, and there appears to be no major obstacle to scaling up to multiple-meter diameters.

Moore, James↗

Implementation of Topology Optimization Using OpenMDAO

Topology optimization is one of the widely known branches among the structural optimization, and it distinguishes itself by a large design domain and versatility. It can determine the optimal design out of an infinite number of configurations, thereby drawing interest from both industry and academia with regards to its applicability to additive manufacturing. However, its implementation is often a daunting task for engineers in practice. One of the issues is the programming effort required whenever the implementation requires any changes, ranging from subtle tweaks to drastic changes in the problem definition, and the derivative computation must be correspondingly updated after any change. The implementation, therefore, is not only time-consuming but also repetitive and susceptible to human-induced errors. In this regard, topology optimization implementations stand to benefit from reusability, ease of restructuring, and modularity. In this work, we propose OpenMDAO, a computational framework for multidisciplinary design optimization (MDO), as a generic platform for topology optimization. Two widely used topology optimization techniquesâ€"density-based and level-setâ€"are implemented as a demonstration. These techniques are implemented in a decomposed manner, with the aid of the modular architecture of OpenMDAO as well as state-of-the-art numerical methods. Variations on the density-based topology optimization approach are shown to demonstrate the modularity and automation for derivative computation that OpenMDAO provides.

Chung, Hayoung↗

Quantum simulation of $\phi^4$ theories in qudit systems

We discuss the implementation of quantum algorithms for lattice $\Phi^4$ theory on circuit quantum electrodynamics (cQED) system. The field is represented on qudits in a discretized field amplitude basis. The main advantage of qudit systems is that its multi-level characteristic allows the field interaction to be implemented only with diagonal single-qudit gates. Considering the set of universal gates formed by the single-qudit phase gate and the displacement gate, we address initial state preparation and single-qudit gate synthesis with variational methods.

Kurkcuoglu, Doga Murat↗