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 955 records · Page 53

The nonlinear dynamics of a spacecraft coupled to the vibration of a contained fluid

The dynamics of a linear spacecraft mode coupled to a nonlinear low gravity slosh of a fluid in a cylindrical tank is investigated. Coupled, nonlinear equations of motion for the fluid-spacecraft dynamics are derived through an assumed mode Lagrangian method. Unlike linear fluid slosh models, this nonlinear slosh model retains two fundamental slosh modes and three secondary modes. An approximate perturbation solution of the equations of motion indicates that the nonlinear coupled system response involves fluid-spacecraft modal resonances not predicted by either a linear, or a nonlinear, uncoupled slosh analysis. Experimental results substantiate the analytical predictions.

Peterson, Lee D.↗

Improved calculation of optimum design sensitivity

Optimum design parameter sensitivity analysis has become an important topic in recent years. The principal reasons for obtaining optimum design sensitivity information with respect to various problem parameters are (l) to predict revised optimum designs, associated with specified perturbations of the problem parameters, without re-optimizing the problem and, (2) to provide sensitivity information in multilevel optimization strategies. Methods for calculating these derivatives were proposed by several authors. The most important drawback in estimating parameter sensitivities by the current methods is that they do not allow for changes in the active constraint set. Changes in the active constraint set produce discontinuities in the optimum design sensitivities, and therefore, a correct formulation has to be cast in terms of directional derivatives. Optimum design parameter sensitivity analysis is addressed in terms of directional derivatives so as to include possible discontinuities and it also presents a method for estimating optimum design sensitivities using finite differences.

Sepulveda, Abdon E.↗

Nonlinear damping model for flexible structures

The study of nonlinear damping problem of flexible structures is addressed. Both passive and active damping, both finite dimensional and infinite dimensional models are studied. In the first part, the spectral density and the correlation function of a single DOF nonlinear damping model is investigated. A formula for the spectral density is established with O(Gamma(sub 2)) accuracy based upon Fokker-Planck technique and perturbation. The spectral density depends upon certain first order statistics which could be obtained if the stationary density is known. A method is proposed to find the approximate stationary density explicitly. In the second part, the spectral density of a multi-DOF nonlinear damping model is investigated. In the third part, energy type nonlinear damping model in an infinite dimensional setting is studied.

Zang, Weijian↗

Dependence of the Element Patterns of HYDROSTAR on Mutual Coupling

HYDROSTAR is a hybrid synthesis radiometer, intended for spaceborne applications, which employs a real aperture (waveguide stick antenna) for resolution along track and employs aperture synthesis to obtain resolution across track. This L-band system is an extension of the successful aircraft prototype, the Electronically Steered Thinned Array Radiometer (ESTAR). A proof-of-concept, full size system was constructed (45 wavelength in the synthesis direction) and subjected to extensive antenna pattern and associated microwave component measurements in an outdoor antenna test facility. HYDROSTAR employs a thinned array of 16 elements each 5.8 inches long in the along track dimension and spanning 9.5 inches across track. Each element is a narrow-wall shunt slot array with 36 slots. The polarization is linear (along-track). In the across track dimension, the antennas are deployed in a minimum redundancy array which has 90 independent baselines spaced in integer multiples of half a wavelength. The closest spacings used are for the first three elements at each end, which are spaced by only one-half wavelength. This study was intended to assess how closely each of 16 stick element patterns compare with their nominal values (as an individual, isolated radiator), when installed in their intended composite thinned array configuration. Extensive pattern measurements of the 16 elements that constitute the HYDROSTAR antenna subsystem were conducted to observe their relative features along the synthesis plane. Data was also collected for the mutual coupling between pairs of selected antennas. All antenna patterns had features different from that of an isolated element indicating some level of interaction among neighboring radiators. Those elements which had nearest neighbors at least five wavelengths away were located near the middle of the array. Their radiation patterns displayed sonic small, symmetric ripple across their full azimuth range. The patterns of elements that lie within 2.5 wavelengths of their neighbors showed stronger and asymmetric features. These are believed to be caused by mutual coupling among these structures. Evidence for this was seen when an antenna position was displaced by 0.05 wavelengths, Its pattern and those of its near neighbors were seen to change. Displacement within the plane of the array were observed to have different effects than displacements out-of-plane. A program of data analysis and theoretical development is in progress to provide a physical interpretation of the properties of these antenna patterns and to develop methods which can optimize the performance of this synthetic aperture imaging system. This includes compensation for pattern asymmetries and element position perturbation.

Weissman, David E.↗

Orbit Determination Accuracy for Comets on Earth-Impacting Trajectories

The results presented show the level of orbit determination accuracy obtainable for long-period comets discovered approximately one year before collision with Earth. Preliminary orbits are determined from simulated observations using Gauss' method. Additional measurements are incorporated to improve the solution through the use of a Kalman filter, and include non-gravitational perturbations due to outgassing. Comparisons between observatories in several different circular heliocentric orbits show that observatories in orbits with radii less than 1 AU result in increased orbit determination accuracy for short tracking durations due to increased parallax per unit time. However, an observatory at 1 AU will perform similarly if the tracking duration is increased, and accuracy is significantly improved if additional observatories are positioned at the Sun-Earth Lagrange points L3, L4, or L5. A single observatory at 1 AU capable of both optical and range measurements yields the highest orbit determination accuracy in the shortest amount of time when compared to other systems of observatories.

Kay-Bunnell, Linda↗

A higher order panel method for general analysis and design applications in subsonic flow

A higher-order panel method is described for numerical solution of boundary-value problems relating to steady inviscid irrotational incompressible subsonic fluid flow in a domain. Both Neumann and Dirichlet boundary conditions are treated; two types of auxiliary conditions are used to remove the degrees of freedom that arise from specifying only the derivative of the perturbation velocity potential. Four general network types and two expansions of the induced potential kernel are employed in the numerical solution. Some results are presented which illustrate the modeling options and numerical characteristics of the method.

Johnson, F. T.↗

Interchange Reconnection Dynamics in A Solar Coronal Pseudo-Streamer

Context. The generation of the slow solar wind remains an open problem in heliophysics. One of the current theories among those aimed at explaining the injection of coronal plasma in the interplanetary medium is based on interchange reconnection. It assumes that the exchange of magnetic connectivity between closed and open fields allows the injection of coronal plasma in the interplanetary medium to travel along the newly reconnected open field. However, the exact mechanism underlying this effect is still poorly understood. Aims. Our objective is to study this scenario in a particular magnetic structure of the solar corona: a pseudo-streamer. This topological structure lies at the interface between open and closed magnetic field and is thought to be involved in the generation of the slow solar wind. Methods. We performed innovative 3D magnetohydrodynamic (MHD) simulations of the solar corona with a pseudo-streamer, using the Adaptively Refined MHD Solver (ARMS). By perturbing the quasi-steady ambient state with a simple photospheric, large-scale velocity flow, we were able to generate a complex dynamics of the open-and-closed boundary of the pseudo-streamer. We studied the evolution of the connectivity of numerous field lines to understand its precise dynamics. Results. We witnessed different scenarios of opening of the magnetic field initially closed under the pseudo-streamer: one-step interchange reconnection dynamics, along with more complex scenarios, including a coupling between pseudo-streamer and helmet streamer, as well as back-and-forth reconnections between open and closed connectivity domains. Finally, our analysis revealed large-scale motions of a newly opened magnetic field high in the corona that may be explained by slipping reconnection. Conclusions. By introducing a new analysis method for the magnetic connectivity evolution based on distinct closed-field domains, this study provides an understanding of the precise dynamics underway during the opening of a closed field, which enables the injection of closed-field, coronal plasma in the interplanetary medium. Further studies shall provide synthetic observations for these diverse outgoing flows, which could be measured by Parker Solar Probe and Solar Orbiter.

magnetic fields↗

Verifiable Adaptive Control with Analytical Stability Margins by Optimal Control Modification

This paper presents a verifiable model-reference adaptive control method based on an optimal control formulation for linear uncertain systems. A predictor model is formulated to enable a parameter estimation of the system parametric uncertainty. The adaptation is based on both the tracking error and predictor error. Using a singular perturbation argument, it can be shown that the closed-loop system tends to a linear time invariant model asymptotically under an assumption of fast adaptation. A stability margin analysis is given to estimate a lower bound of the time delay margin using a matrix measure method. Using this analytical method, the free design parameter n of the optimal control modification adaptive law can be determined to meet a specification of stability margin for verification purposes.

Nguyen, Nhan T.↗

Probabilistic Structures Analysis Methods (PSAM) for select space propulsion system components

The basic formulation for probabilistic finite element analysis is described and demonstrated on a few sample problems. This formulation is based on iterative perturbation that uses the factorized stiffness on the unperturbed system as the iteration preconditioner for obtaining the solution to the perturbed problem. This approach eliminates the need to compute, store and manipulate explicit partial derivatives of the element matrices and force vector, which not only reduces memory usage considerably, but also greatly simplifies the coding and validation tasks. All aspects for the proposed formulation were combined in a demonstration problem using a simplified model of a curved turbine blade discretized with 48 shell elements, and having random pressure and temperature fields with partial correlation, random uniform thickness, and random stiffness at the root.

Source record↗

Efficient algorithms for use in probabilistic finite element analysis

This paper investigates the use of Fast Probability Integration (FPI) algorithms in a Finite Element environment. A method allowing the representation of correlated fields in terms of a vector of uncorrelated transformed variables, based on the spectral decomposition of the variance-covariance matrix is developed. The response of the deterministic model corresponding to selected perturbations of these uncorrelated variables is then obtained via a Newton-type iterative scheme. The results of the perturbed problems are used to construct a local representation of the model's behavior in the neighborhood of the deterministic state, which the FPI algorithm will use to estimate the reliability of the system. Although the proposed strategy has thus far only been applied to linear elastostatics, the extension of the method to a broader class of problems appears to be feasible.

Dias, J. B.↗

On the simultaneous improvement of a satellite orbit and determination of sea surface topography using altimeter data

A method is presented in satellite altimetry that attempts to simultaneously determine the geoid and sea surface toography with minimum wavelengths of about 500 km and to reduce the radial orbit errors caused by geopotential uncertainties. The modeling of these errors is made using the linearized Lagrangian perturbation theory. Observation equations are developed using sea surface heights and crossover discrepancies as observables. A minimum variance solution with prior information can then provide estimates of parametrs representing the sea surface topography and corrections to the orbit. The potential of the method is demonstrated in a solution where simulated geopotential errors and the Levitus sea surface topography are used to generate the observables for a Seasat 3 day arc. The simulation results suggest that the method can be used to efficiently process real altimeter data.

Engelis, Theodossios↗

High speed flow past wings

The analytical solution to the transonic small perturbation equation which describes steady compressible flow past finite wings at subsonic speeds can be expressed as a nonlinear integral equation with the perturbation velocity potential as the unknown function. This known formulation is substituted by a system of nonlinear algebraic equations to which various methods are applicable for its solution. Due to the presence of mathematical discontinuities in the flow solutions, however, a main computational difficulty was to ensure uniqueness of the solutions when local velocities on the wing exceeded the speed of sound. For continuous solutions this was achieved by embedding the algebraic system in an one-parameter operator homotopy in order to apply the method of parametric differentiation. The solution to the initial system of equations appears then as a solution to a Cauchy problem where the initial condition is related to the accompanying incompressible flow solution. In using this technique, however, a continuous dependence of the solution development on the initial data is lost when the solution reaches the minimum bifurcation point. A steepest descent iteration technique was therefore, added to the computational scheme for the calculation of discontinuous flow solutions. Results for purely subsonic flows and supersonic flows with and without compression shocks are given and compared with other available theoretical solutions.

Norstrud, H.↗

Transmission of sound through high subsonic flows in non-uniform ducts

The transmission of sound through a near-sonic throat in a variable area duct is studied using the method of matched asymptotic expansions applied to a quasi-one-dimensional model. Nonlinear interactions between the mean flow and the acoustic perturbation in the throat region are shown to cause significant effects on the sound field upstream and downstream of the throat. Numerical results are presented for several special applications of the general nonlinear theory.

Myers, M. K.↗

A model for helicopter guidance on spiral trajectories

A point mass model is developed for helicopter guidance on spiral trajectories. A fully coupled set of state equations is developed and perturbation equations suitable for 3-D and 4-D guidance are derived and shown to be amenable to conventional state variable feedback methods. Control variables are chosen to be the magnitude and orientation of the net rotor thrust. Using these variables reference controls for nonlevel accelerating trajectories are easily determined. The effects of constant wind are shown to require significant feedforward correction to some of the reference controls and to the time. Although not easily measured themselves, the controls variables chosen are shown to be easily related to the physical variables available in the cockpit.

Mendenhall, S.↗

Evaluation of 3-D graphics software: A case study

An efficient 3-D geometry graphics software package which is suitable for advanced design studies was developed. The advanced design system is called GRADE--Graphics for Advanced Design. Efficiency and ease of use are gained by sacrificing flexibility in surface representation. The immediate options were either to continue development of GRADE or to acquire a commercially available system which would replace or complement GRADE. Test cases which would reveal the ability of each system to satisfy the requirements were developed. A scoring method which adequately captured the relative capabilities of the three systems was presented. While more complex multi-attribute decision methods could be used, the selected method provides all the needed information without being so complex that it is difficult to understand. If the value factors are modestly perturbed, system Z is a clear winner based on its overall capabilities. System Z is superior in two vital areas: surfacing and ease of interface with application programs.

Lores, M. E.↗

Neutral Atmosphere Properties Determining D-region Electron Densities

The increasing discoveries of various manifestations of meteorological control of the D region ionization and the growth of techniques for its measurement provide a challenge to meteorologists to test their insight into middle atmosphere processes with the physical interpretation of D layer phenomena. Models for ion production due to photoionization of minor atmospheric nitric oxide by quasi-monochromatic solar Lyman-alpha radiation are presented. A ground based measuring technique using low frequency radio reflection heights is briefly described and an approach to the interpretation of data acquired by this method is discussed. It is shown that D region electron density variations can provide an efficient diagnostic tool for the detection of perturbations of the circulation state of the middle atmosphere.

Taubenheim, J.↗

Orbital stability of the unseen solar companion linked to periodic extinction events

Evidence from three-dimensional numerical modelling is presented that only cometary orbits with a limited range in inclination with respect to the galactic plane are formally stable for the length of time required to cause periodic extinction events. The calculations were done using Cowell's method employing a fourth-order Runge-Kutta integration scheme in an inertial reference frame in orbit about the Galaxy. Tidal perturbations in the radial direction due to the Galaxy and the Coriolis forces are included. The vertical component of the gravitational field of the galactic disk is superimposed on these forces. The results indicate that orbits for Nemesis that are inclined at more than 30 deg to the galactic plane are not allowed and suggests that the search for Nemesis should be concentrated toward the plane of the Galaxy. Perturbations by passing stars or molecular clouds may make even the low-inclination orbits unstable.

Torbett, M. V.↗

A review of hot wire anemometry in transonic flows

The present paper provides a review of hot wire anemometry for compressible flows, giving particular attention to the transonic flow problem. It is pointed out that the first and most important definitive work in hot wire anemometry for compressible flows was reported by Kovasznay (1953). The existence of three independent fluctuating modes in compressible flows for small perturbations was found, taking into account the vorticity mode, the entropy mode, and the sound-wave mode. A review of Kovasznays' method for supersonic flows is also presented, and advances reported by Markovin (1956) are examined. Attention is given to experiments conducted by Horstman and Rose (1977), a general solution to the hot wire problem at transonic conditions sought by Stainback et al. (1983), and some apparent problems.

Stainback, P. C.↗