Search NASA⌕ Search

SEARCH · Search NASA

Results for “Variational calculus”

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.

224 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

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.↗
DOE OSTI2017

Convective mixing in vertically-layered porous media: The linear regime and the onset of convection

In this paper, we study the effect of permeability heterogeneity on the stability of gravitationally unstable, transient, diffusive boundary layers in porous media. Permeability is taken to vary periodically in the horizontal plane normal to the direction of gravity. In contrast to the situation for vertical permeability variation, the horizontal perturbation structures are multimodal. We therefore use a two-dimensional quasi-steady eigenvalue analysis as well as a complementary initial value problem to investigate the stability behavior in the linear regime, until the onset of convection. We find that thick permeability layers enhance instability compared with thin layers when heterogeneity is increased. On the contrary, for thin layers the instability is weakened progressively with increasing heterogeneity to the extent that the corresponding homogeneous case is more unstable. For high levels of heterogeneity, we find that a small change in the permeability field results in large variations in the onset time of convection, similar to the instability event in the linear regime. However, this trend does not persist unconditionally because of the reorientation of vorticity pairs due to the interaction of evolving perturbation structures with heterogeneity. Consequently, an earlier onset of instability does not necessarily imply an earlier onset of convection. A resonant amplification of instability is observed within the linear regime when the dominant perturbation mode is equal to half the wavenumber of permeability variation. On the other hand, a substantial damping occurs when the perturbation mode is equal to the harmonic and sub-harmonic components of the permeability wavenumber. Finally, the phenomenon of such harmonic interactions influences both the onset of instability as well as the onset of convection.

36 MATERIALS SCIENCE↗
NASA NTRS2010

Modeling the Impact of Space Suit Components and Anthropometry on the Center of Mass of a Seated Crewmember

The Crew Impact Attenuation System (CIAS) is the energy-absorbing strut concept that dampens Orion Crew Exploration Vehicle (CEV) landing loads to levels sustainable by the crew. Significant COM variations across suited crew configurations would amplify the inertial effects of the pallet and potentially create unacceptable crew loading during launch and landing. The objective of this study was to obtain data needed for dynamic simulation models by quantifying the effects of posture, suit components, and the expected range of anthropometry on the COM of a seated individual. Several elements are required for the COM calculation of a suited human in a seated position: anthropometry, body segment mass, suit component mass, suit component location relative to the body, and joint angles defining the seated posture. Three-dimensional (3D) human body models, suit mass data, and vector calculus were utilized to compute the COM positions for 12 boundary manikins in two different seated postures. The analysis focused on two objectives: (1) quantify how much the wholebody COM varied from the smallest to largest subject and (2) quantify the effects of the suit components on the overall COM in each seat configuration. The location of the anterior-posterior COM varied across all boundary manikins by about 7 cm, and the vertical COM varied by approximately 9 to 10 cm. The mediolateral COM varied by 1.2 cm from the midline sagittal plane for both seat configurations. The suit components caused an anterior shift of the total COM by approximately 2 cm and a shift to the right along the mediolateral axis of 0.4 cm for both seat configurations. When the seat configuration was in the standard posture the suited vertical COM shifted inferiorly by as much as 1 cm, whereas in the CEV posture the vertical COM had no appreciable change. These general differences were due to the high proportion of suit mass located in the boots and lower legs and their corresponding distance from the body COM, as well as to the prevalence of suit components on the right side of the body.

Blackledge, Christopher↗
NASA NTRS1975

The geometric approach to sets of ordinary differential equations and Hamiltonian dynamics

The calculus of differential forms is used to discuss the local integration theory of a general set of autonomous first order ordinary differential equations. Geometrically, such a set is a vector field V in the space of dependent variables. Integration consists of seeking associated geometric structures invariant along V: scalar fields, forms, vectors, and integrals over subspaces. It is shown that to any field V can be associated a Hamiltonian structure of forms if, when dealing with an odd number of dependent variables, an arbitrary equation of constraint is also added. Families of integral invariants are an immediate consequence. Poisson brackets are isomorphic to Lie products of associated CT-generating vector fields. Hamilton's variational principle follows from the fact that the maximal regular integral manifolds of a closed set of forms must include the characteristics of the set.

Estabrook, F. B.↗
NASA NTRS1991

On the Power of Abstract Interpretation

Increasingly sophisticated applications of static analysis place increased burden on the reliability of the analysis techniques. Often, the failure of the analysis technique to detect some information my mean that the time or space complexity of the generated code would be altered. Thus, it is important to precisely characterize the power of static analysis techniques. We follow the approach of Selur et. al. who studied the power of strictness analysis techniques. Their result can be summarized by saying 'strictness analysis is perfect up to variations in constants.' In other words, strictness analysis is as good as it could be, short of actually distinguishing between concrete values. We use this approach to characterize a broad class of analysis techniques based on abstract interpretation including, but not limited to, strictness analysis. For the first-order case, we consider abstract interpretations where the abstract domain for data values is totally ordered. This condition is satisfied by Mycroft's strictness analysis that of Sekar et. al. and Wadler's analysis of list-strictness. For such abstract interpretations, we show that the analysis is complete in the sense that, short of actually distinguishing between concrete values with the same abstraction, it gives the best possible information. We further generalize these results to typed lambda calculus with pairs and higher-order functions. Note that products and function spaces over totally ordered domains are not totally ordered. In fact, the notion of completeness used in the first-order case fails if product domains or function spaces are added. We formulate a weaker notion of completeness based on observability of values. Two values (including pairs and functions) are considered indistinguishable if their observable components are indistinguishable. We show that abstract interpretation of typed lambda calculus programs is complete up to this notion of indistinguishability. We use denotationally-oriented arguments instead of the detailed operational arguments used by Selur et. al.. Hence, our proofs are much simpler. They should be useful for further future improvements.

Reddy, Uday S.↗
DOE OSTI2018

The Lorenzo–Hartley’s function for fractional calculus and its applications pertaining to fractional order modelling of anomalous relaxation in dielectrics

Considering the applications and a great significance of the study pertaining to generalized functions in applied sciences, the author investigated the Lorenzo–Hartley’s function which is a generalization of classical functions widely used in fractional calculus. The Laplace transform pairs are derived and the requirements for the Lorenzo–Hartley’s function G{sub ν,μ,δ}(a,c,t) to be completely monotonic (for t>0) are investigated. The author has shown the applications of this generalized function to describe the relaxation models, particularly in dielectrics. The Lorenzo–Hartley’s function G{sub ν,μ,δ}(a,c,t) of a real variable t is considered to investigate a computable mathematical framework for standard Debye and non-Debye relaxation processes in dielectric materials. The non-negative spectral distribution function is obtained for the corresponding response function. It is also demonstrated that the classical models like Cole–Cole (C–C), Davidson–Cole (D–C), and Havriliak–Negami (H–N) for non-Debye relaxation and standard Debye relaxation are the particular cases of the Lorenzo–Hartley’s function. Some of the study cases are also worked out to visualize the effects of variations of parameters on the response function and corresponding spectral distribution function. A generalized and unified fractional relaxation differential equation which governs response functions for classical dielectrics models pertaining to non-Debye relaxation (C–C, D–C, and H–N) is also established in the present investigation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Energy and momentum conservation in the Euler–Poincaré formulation of local Vlasov–Maxwell-type systems

The action principle by Low for the classic Vlasov–Maxwell system contains a mix of Eulerian and Lagrangian variables. This renders the Noether analysis of reparametrization symmetries inconvenient, especially since the well-known energy- and momentum-conservation laws for the system are expressed in terms of Eulerian variables only. While an Euler–Poincaré formulation of Vlasov–Maxwell-type systems, effectively starting with Low's action and using constrained variations for the Eulerian description of particle motion, has been known for a while Cendra et al, it is hard to come by a documented derivation of the related energy- and momentum-conservation laws in the spirit of the Euler–Poincaré machinery. To our knowledge only one such derivation exists in the literature so far, dealing with the so-called guiding-center Vlasov–Darwin system Sugama et al. The present exposition discusses a generic class of local Vlasov–Maxwell-type systems, with a conscious choice of adopting the language of differential geometry to exploit the Euler–Poincaré framework to its full extent. After reviewing the transition from a Lagrangian picture to an Eulerian one, we demonstrate how symmetries generated by isometries in space lead to conservation laws for linear- and angular-momentum density and how symmetry by time translation produces a conservation law for energy density. Here, we also discuss what happens if no symmetries exist. Finally, two explicit examples will be given—the classic Vlasov–Maxwell and the drift-kinetic Vlasov–Maxwell—and the results expressed in the language of regular vector calculus for familiarity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Fractional Laplacian spectral approach to turbulence in a dusty plasma monolayer

This work presents an analytical investigation of anomalous diffusion and turbulence in a dusty plasma monolayer, where energy transport across scales leads to the spontaneous formation of spatially disordered patterns. Many-body simulations of 10,000-particle dusty plasma monolayers are used to demonstrate how the global dynamics depend on the statistical properties of the dust assembly for realistic laboratory conditions. We find that disorder due to variations in the dust size distribution and charge-driven nonlocal interactions resulting in anomalous dust diffusion are key factors for the onset of instabilities. The resulting dynamics exhibit features of inertial turbulence over slightly more than half a decade of scales proportional or smaller than the Debye shielding length. These processes are examined analytically using a recently developed Fractional Laplacian Spectral (FLS) technique, which identifies the active energy channels as a function of scale, disorder concentration, and features of the nonlocal-interactions. The predictions from the theoretical (spectral) analysis demonstrate agreement with the results from the many-body (kinetic) simulations, thus providing a powerful tool for the study of active turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗