Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Perturbative 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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Unsteady lifting-line theory as a singular-perturbation problem

Unsteady lifting-line theory is developed for a wing of large aspect ratio oscillating at low frequency in inviscid incompressible flow. The wing is assumed to have a rigid chord but a flexible span. Use of the method of matched asymptotic expansions reduces the problem from a singular integral equation to quadrature. The pressure field and airloads, for a prescribed wing shape and motion, are obtained in closed form as expansions in inverse aspect ratio. A rigorous definition of unsteady induced downwash is also obtained. Numerical calculations are presented for an elliptic wing in pitch and heave; compared with numerical lifting-surface theory, computation time is reduced significantly. The present work also identifies and resolves errors in the unsteady lifting line theory of James (1975), and points out a limitation in that of Van Holten (1975, 1976, 1977).

Ahmadi, A.↗
NASA NTRS1989

Actuator placement by degree of controllability including the effect of actuator mass

Various concepts of the degree of controllability (DOC) have been previously developed for the purpose of choosing actuator locations in the control of large flexible spacecraft. A perturbation technique is presented to efficiently account for the actuator mass in the computation of the fuel-optimal, time-optimal, and energy-optimal DOC definitions. Methods are developed for computing each of the three DOCs, when the mass of the actuator is nonnegligible. The energy-optimal DOC is employed to find optimal actuator locations for a large spacecraft under study by the NASA Langley Research Center. The mass of the actuator is seen to have a significant effect on the optimal actuator locations and is shown to spread the actuator locations more evenly for large mass actuators. Results indicate that when the actuator mass is large when compared to the structural mass it should be included in selecting actuator locations.

Longman, R. W.↗
NASA NTRS1992

Dynamics and stability of a tethered centrifuge in low earth orbit

The three-dimensional attitude dynamics of a spaceborne tethered centrifuge for artificial gravity experiments in low earth orbit is analyzed using two different methods. First, the tethered centrifuge is modeled as a dumbbell with a straight viscoelastic tether, point tip-masses, and sophisticated environmental models such as nonspherical gravity, thermal perturbations, and a dynamic atmospheric model. The motion of the centrifuge during spin-up, de-spin, and steady-rotation is then simulated. Second, a continuum model of the tether is developed for analyzing the stability of lateral tether oscillations. Results indicate that the maximum fluctuation about the 1-g radial acceleration level is less than 0.001 g; the time required for spin-up and de-spin is less than one orbit; and lateral oscillations are stable for any practical values of the system parameters.

Quadrelli, B. M.↗
NASA NTRS1996

Implementation of Hybrid V-Cycle Multilevel Methods for Mixed Finite Element Systems with Penalty

The goal of this paper is the implementation of hybrid V-cycle hierarchical multilevel methods for the indefinite discrete systems which arise when a mixed finite element approximation is used to solve elliptic boundary value problems. By introducing a penalty parameter, the perturbed indefinite system can be reduced to a symmetric positive definite system containing the small penalty parameter for the velocity unknown alone. We stabilize the hierarchical spatial decomposition approach proposed by Cai, Goldstein, and Pasciak for the reduced system. We demonstrate that the relative condition number of the preconditioner is bounded uniformly with respect to the penalty parameter, the number of levels and possible jumps of the coefficients as long as they occur only across the edges of the coarsest elements.

Lai, Chen-Yao G.↗
NASA NTRS2015

Numerical Computation of a Continuous-thrust State Transition Matrix Incorporating Accurate Hardware and Ephemeris Models

A significant body of work exists showing that providing a nonlinear programming (NLP) solver with expressions for the problem constraint gradient substantially increases the speed of program execution and can also improve the robustness of convergence, especially for local optimizers. Calculation of these derivatives is often accomplished through the computation of spacecraft's state transition matrix (STM). If the two-body gravitational model is employed as is often done in the context of preliminary design, closed form expressions for these derivatives may be provided. If a high fidelity dynamics model, that might include perturbing forces such as the gravitational effect from multiple third bodies and solar radiation pressure is used then these STM's must be computed numerically. We present a method for the power hardward model and a full ephemeris model. An adaptive-step embedded eight order Dormand-Prince numerical integrator is discussed and a method for the computation of the time of flight derivatives in this framework is presented. The use of these numerically calculated derivatieves offer a substantial improvement over finite differencing in the context of a global optimizer. Specifically the inclusion of these STM's into the low thrust missiondesign tool chain in use at NASA Goddard Spaceflight Center allows for an increased preliminary mission design cadence.

Low Thrust↗
NASA NTRS2015

Plasma Perturbations in High-Speed Probing of Hall Thruster Discharge Chambers: Quantification and Mitigation

An experimental investigation is presented to quantify the effect of high-speed probing on the plasma parameters inside the discharge chamber of a 6-kW Hall thruster. Understanding the nature of these perturbations is of significant interest given the importance of accurate plasma measurements for characterizing thruster operation. An array of diagnostics including a high-speed camera and embedded wall probes is employed to examine in real time the changes in electron temperature and plasma potential induced by inserting a high-speed reciprocating Langmuir probe into the discharge chamber. It is found that the perturbations onset when the scanning probe is downstream of the electron temperature peak, and that along channel centerline, the perturbations are best characterized as a downstream shift of plasma parameters by 15-20% the length of the discharge chamber. A parametric study is performed to investigate techniques to mitigate the observed probe perturbations including varying probe speed, probe location, and operating conditions. It is found that the perturbations largely disappear when the thruster is operated at low power and low discharge voltage. The results of this mitigation study are discussed in the context of recommended methods for generating unperturbed measurements of the discharge chamber plasma.

Jorns, Benjamin A.↗
NASA NTRSDate not supplied

Using Anharmonic PAH Vibrational Spectra to Simulate Fully Anharmonic Cascade Emission Spectra

In recent years, we have developed vibrational second-order perturbation theory (VPT2) approaches that can be applied to very large molecules, such as polycyclic aromatic hydrocarbon (PAH) molecules. One of the main difficulties in applying standard VPT2 methods to PAHs is the large number of Fermi resonances that occur, especially in the C-H stretching region of the infrared spectrum. In order to properly account for these resonances, it is necessary to set up resonance polyads for each irreducible representation of the symmetry point group to which the molecule belongs, and for specific wavenumber regions, such as the C-H stretching region. Diagonalization of the polyads allows us to determine vibrational band positions, the intensity sharing due to the resonances, and hence there will be many more vibrational bands present than there are C-H bonds. Using these data, a library of fully anharmonic temperature-dependent spectra, with proper treatment of polyad resonances, can be determined and subsequently used to model the cascade emission spectra of PAH molecules – the type of spectra directly observed by astronomers. The theoretical approach will be described, with comparison to experiment where possible, and the latest fully anharmonic PAH cascade emission spectra will be discussed.

Timothy J Lee↗
NASA NTRS2009

A Critical Study of Agglomerated Multigrid Methods for Diffusion

Agglomerated multigrid techniques used in unstructured-grid methods are studied critically for a model problem representative of laminar diffusion in the incompressible limit. The studied target-grid discretizations and discretizations used on agglomerated grids are typical of current node-centered formulations. Agglomerated multigrid convergence rates are presented using a range of two- and three-dimensional randomly perturbed unstructured grids for simple geometries with isotropic and highly stretched grids. Two agglomeration techniques are used within an overall topology-preserving agglomeration framework. The results show that multigrid with an inconsistent coarse-grid scheme using only the edge terms (also referred to in the literature as a thin-layer formulation) provides considerable speedup over single-grid methods but its convergence deteriorates on finer grids. Multigrid with a Galerkin coarse-grid discretization using piecewise-constant prolongation and a heuristic correction factor is slower and also grid-dependent. In contrast, grid-independent convergence rates are demonstrated for multigrid with consistent coarse-grid discretizations. Actual cycle results are verified using quantitative analysis methods in which parts of the cycle are replaced by their idealized counterparts.

Thomas, James L.↗
NASA NTRS2004

ENSO Bred Vectors in Coupled Ocean-Atmosphere General Circulation Models

The breeding method has been implemented in the NASA Seasonal-to-Interannual Prediction Project (NSIPP) Coupled General Circulation Model (CGCM) with the goal of improving operational seasonal to interannual climate predictions through ensemble forecasting and data assimilation. The coupled instability as cap'tured by the breeding method is the first attempt to isolate the evolving ENSO instability and its corresponding global atmospheric response in a fully coupled ocean-atmosphere GCM. Our results show that the growth rate of the coupled bred vectors (BV) peaks at about 3 months before a background ENSO event. The dominant growing BV modes are reminiscent of the background ENSO anomalies and show a strong tropical response with wind/SST/thermocline interrelated in a manner similar to the background ENSO mode. They exhibit larger amplitudes in the eastern tropical Pacific, reflecting the natural dynamical sensitivity associated with the presence of the shallow thermocline. Moreover, the extratropical perturbations associated with these coupled BV modes reveal the variations related to the atmospheric teleconnection patterns associated with background ENSO variability, e.g. over the North Pacific and North America. A similar experiment was carried out with the NCEP/CFS03 CGCM. Comparisons between bred vectors from the NSIPP CGCM and NCEP/CFS03 CGCM demonstrate the robustness of the results. Our results strongly suggest that the breeding method can serve as a natural filter to identify the slowly varying, coupled instabilities in a coupled GCM, which can be used to construct ensemble perturbations for ensemble forecasts and to estimate the coupled background error covariance for coupled data assimilation.

Yang, S. C.↗
NASA NTRS2017

NIAC Swarm Flyby Gravimetry Phase II Report

This NASA Innovative Advanced Concept (NIAC) grant has enabled the research and development of a method for conducting small body gravimetry from a spacecraft, using relative measurements to a set of deployed test-masses. The test-masses are tracked from a host spacecraft, which dispenses them near to the small body's surface. Thanks to this close proximity, the probes' orbits can be highly perturbed, which yields useful gravimetric measurements. The most readily achievable approach for tracking the probes is to use an optical instrument on- board the spacecraft. The probes then need only be reflective to sunlight. This implementation, called optical gravimetry (OpGrav), has the fewest requirements for the host spacecraft and probes.The results of this study indicate that OpGrav is feasible and offers meaningful improvement over existing methods. Parametric studies suggest roughly an order of magnitude improvement in accuracy or asteroid accessibility (how small an asteroid one can measure) over Earth-based Doppler-only mass estimation. This exponentially expands the number of potential near-Earth objects that one could study, which has implications for planetary defense.As a sample mission, we evaluated OpGrav as an added instrument on a main- belt asteroid tour mission. In this case, simulations show that OpGrav would increase the number of asteroid mass estimates from 3 of 9 to 7 of 9. That is, OpGrav has sufficient sensitivity to offer utility in missions for which it is not explicitly designed for.We designed and fabricated a prototype hardware implementation for this concept called the Small-body In-situ Multi-probe Mass Estimation Experiment (SIMMEE). This hardware provides a basis for many inputs into the simulations and grounds the models with physical values. The primary design driver for the hardware is a long life, on the order of five years prior to operation, and a need for high pointing accuracy to enable flybys of the smallest objects.The next steps include further hardware testing and extension of the concept to rendezvous cases. We believe that this concept offers planetary scientists a new and relevant means of better understanding small-bodies.

Atchison, Justin A.↗
NASA NTRS2010

Method of Obtaining High Resolution Intrinsic Wire Boom Damping Parameters for Multi-Body Dynamics Simulations

The goal of NASA's Magnetospheric MultiScale (MMS) mission is to understand magnetic reconnection with sensor measurements from four spinning satellites flown in a tight tetrahedron formation. Four of the six electric field sensors on each satellite are located at the end of 60- meter wire booms to increase measurement sensitivity in the spin plane and to minimize motion coupling from perturbations on the main body. A propulsion burn however, might induce boom oscillations that could impact science measurements if oscillations do not damp to values on the order of 0.1 degree in a timely fashion. Large damping time constants could also adversely affect flight dynamics and attitude control performance. In this paper, we will discuss the implementation of a high resolution method for calculating the boom's intrinsic damping, which was used in multi-body dynamics simulations. In summary, experimental data was obtained with a scaled-down boom, which was suspended as a pendulum in vacuum. Optical techniques were designed to accurately measure the natural decay of angular position and subsequently, data processing algorithms resulted in excellent spatial and temporal resolutions. This method was repeated in a parametric study for various lengths, root tensions and vacuum levels. For all data sets, regression models for damping were applied, including: nonlinear viscous, frequency-independent hysteretic, coulomb and some combination of them. Our data analysis and dynamics models have shown that the intrinsic damping for the baseline boom is insufficient, thereby forcing project management to explore mitigation strategies.

Yew, Alvin G.↗
NASA NTRS1966

Vainshtein Method Applied to Electron Impact Excitation and Ionization of Atoms.

In most calculations on inelastic collision of electrons with atoms there has been a general lack of agreement with observation mainly due to the treatment of the interelectronic interaction as a perturbation. In a recent calculation by Vainshtein et al., full account is given to this interaction by treating the problem as a binary collision between the incident and atomic electrons with the motion of their center of mass in the Coulomb field of the nucleus, and marked improvement in comparison of theory and observation has resulted. Three approximations are applied in this calculation which will be described briefly.

COLLISION↗
NASA NTRS1989

Gradient methods for identification of distributed parameter systems

Parameter-dependence properties are developed in the context of a linear abstract Cauchy problem governed by a parameter-dependent operator. Of particular interest are problems in which the parameter induces an unbounded perturbation of the evolution operator. Conditions are stated under which the derivative of the state with respect to the parameter possesses certain smoothness properties. These properties lead to local convergence results for a gradient-based estimation algorithm based on quasi-linearization. Numerical results concerning a delay-differential equation indicate the effect of a forcing function on parameter sensitivity.

Brewer, Dennis W.↗
NASA NTRS1973

Gravitational perturbations of equatorial orbits.

Asymptotic solutions are developed for the motion of a geocentric satellite in the equatorial plane due to gravitational perturbations such as nonsphericity (especially oblateness) of the primary body. Axisymmetric potentials are considered. A class of transformations is developed and the equations of motion are solved by the method of generalized multiple scales. Further it is shown that the equations of motion can be transformed into the required form to within any specified degree of accuracy. The transformations form an Abelian group of infinite order which leaves the differential equations of motion invariant. Solutions are developed in terms of elementary functions instead of elliptic or other higher transcendental functions and are shown to agree with known results.

Ramnath, R.↗
NASA NTRS1979

Robust stability of linear systems: Some computational considerations

The cases of both additive and multiplicative perturbations were discussed and a number of relationships between the two cases were given. A number of computational aspects of the theory were also discussed, including a proposed new method for evaluating general transfer or frequency response matrices. The new method is numerically stable and efficient, requiring only operations to update for new values of the frequency parameter.

Laub, A. J.↗
NASA NTRS1985

Unsteady transonic flows - Introduction, current trends, applications

The computational treatment of unsteady transonic flows is discussed, reviewing the historical development and current techniques. The fundamental physical principles are outlined; the governing equations are introduced; three-dimensional linearized and two-dimensional linear-perturbation theories in frequency domain are described in detail; and consideration is given to frequency-domain FEMs and time-domain finite-difference and integral-equation methods. Extensive graphs and diagrams are included.

Yates, E. C., Jr.↗
NASA NTRS1989

Resonance in cylindrical-rectangular and wraparound microstrip structures

A rigorous analysis of the resonance frequency problem of both the cylindrical-rectangular and the wraparound microstrip structure is presented. The problem is formulated in terms of a set of vector integral equations. Using Galerkin's method to solve the integral equations, the complex resonance frequencies are studied with sinusoidal basis functions which incorporate the edge singularity. The complex resonance frequencies are computed using a perturbation approach. Modes suitable for resonator or antenna applications are investigated. The edge singularity of the patch current is shown to have no significant effect on the accuracy of the results. It is shown that the HE10 modes of the cylindrical-rectangular and wraparound patches are more appropriate for resonator applications. The HE01 and TE01 modes of the cylindrical-rectangular and wraparound patches, respectively, are efficient radiating modes.

Ali, Sami M.↗
NASA NTRS1994

Towards enhancing and delaying disturbances in free shear flows

The family of shear flows comprising the jet, wake, and the mixing layer are subjected to perturbations in an inviscid incompressible fluid. By modeling the basic mean flows as parallel with piecewise linear variations for the velocities, complete and general solutions to the linearized equations of motion can be obtained in closed form as functions of all space variables and time when posed as an initial value problem. The results show that there is a continuous as well as the discrete spectrum that is more familiar in stability theory and therefore there can be both algebraic and exponential growth of disturbances in time. These bases make it feasible to consider control of such flows. To this end, the possibility of enhancing the disturbances in the mixing layer and delaying the onset in the jet and wake is investigated. It is found that growth of perturbations can be delayed to a considerable degree for the jet and the wake but, by comparison, cannot be enhanced in the mixing layer. By using moving coordinates, a method for demonstrating the predominant early and long time behavior of disturbances in these flows is given for continuous velocity profiles. It is shown that the early time transients are always algebraic whereas the asymptotic limit is that of an exponential normal mode. Numerical treatment of the new governing equations confirm the conclusions reached by use of the piecewise linear basic models. Although not pursued here, feedback mechanisms designed for control of the flow could be devised using the results of this work.

Criminale, W. O.↗