A proposed hybrid-computer method for functional optimization
Random perturbation procedure for control variable functions optimization using hybrid computer method for analog solutions and digital storage
Korn, G. A.↗
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.
A harmonic L4 orbit for the very restricted four-body problem
Harmonic L-4 orbit for very restricted four-body problem determined by method of general perturbations using Chebyshev series
EARTH-MOON SYSTEM↗
Study of diatomic molecules. 2: Intensities
The theory of perturbations, giving the diatomic effective Hamiltonian, is used for calculating actual molecular wave functions and intensity factors involved in transitions between states arising from Hund's coupling cases a,b, intermediate a-b, and c tendency. The Herman and Wallis corrections are derived, without any knowledge of the analytical expressions of the wave functions, and generalized to transitions between electronic states with whatever symmetry and multiplicity. A general method for studying perturbed intensities is presented using primarily modern spectroscopic numerical approaches. The method is used in the study of the ScO optical emission spectrum.
Femenias, J. L.↗
A parallel vectorized implementation of triple excitations in CCSD(T) - Application to the binding energies of the AlH3, AlH2F, AlHF2 and AlF3 dimers
An efficient method for various noniterative estimates of connected triple excitations in coupled-cluster theory is outlined and related to a similar expression occurring in Moller-Plesset perturbation theory. The method is highly vectorized and capable of utilizing multiple processors on a shared-memory machine, leading to computational rates in excess of one billion floating-point operations per second on four processors of a CRAY Y-MP. Using the new procedure, the binding energies of the D(2h) diborane-type dimers of AlH3, AlH2F, AlHF2, and AlF3 have been determined to be 32, 40, 20, and 47 kcal/mol, respectively. For Al2F6, the correlation procedure includes 232 molecular orbitals and over 1.5 x 10 to the 6th single and double coupled-cluster amplitudes, effectively accounting for over 2 x 10 to the 9th connected triple excitations.
Rendell, Alistair P.↗
Perturbations in rectangular coordinates by iteration.
Orbits computation by Picard successive approximations method, discussing iterative numerical perturbation techniques
Broucke, R.↗
Asymptotic continuation method for trajectory optimization
A continuation method is applied to a singular perturbation parameter to obtain a new numerical method for computing optimal trajectories. This method allows one to use simply calculated reduced order approximations as starting solutions and continue the perturbation parameter until the optimal full-order solution is obtained. The method does not require the calculation of higher order correction terms nor does it require the perturbation parameter to be small - thus, it has potentially superior convergence properties compared to conventional asymptotic expansions when the perturbation parameter is large. A simple trajectory optimization problem is considered to illustrate the method.
Washburn, R. B., Jr.↗
An explicit Runge-Kutta method for turbulent reacting flow calculations
The paper presents a numerical method for the solution of the conservation equations governing steady, reacting, turbulent viscous flow in two-dimensional geometries, in both Cartesian and axisymmetric coordinates. These equations are written in Favre-averaged form and closed with a first order model. A two-equation K-epsilon model, where low Reynolds number and compressibility effects are included, and a modified eddy-break up model are used to simulate fluid mechanics turbulence, chemistry and turbulence-combustion interaction. The solution is obtained by using a pseudo-unsteady method with improved perturbation propagation properties. The equations are discretized in space by using a finite volume formulation. An explicit multi-stage dissipative Runge-Kutta algorithm is then used to advance the flow equations in the pseudo-time. The method is applied to the computation of both diffusion and premixed turbulent reacting flows. The computed temperature distributions compare favorably with experimental data.
Boretti, A. A.↗
A method for comparing trajectories in optimum linear perturbation guidance schemes.
Comparison of trajectories in optimum linear perturbation guidance
Powers, W. F.↗
Advances in Land Data Assimilation at the NASA Goddard Space Flight Center
Research in land surface data assimilation has grown rapidly over the last decade. In this presentation we provide a brief overview of key research contributions by the NASA Goddard Space Flight Center (GSFC). The GSFC contributions to land assimilation primarily include the continued development and application of the Land Information System (US) and the ensemble Kalman filter (EnKF). In particular, we have developed a method to generate perturbation fields that are correlated in space, time, and across variables and that permit the flexible modeling of errors in land surface models and observations, along with an adaptive filtering approach that estimates observation and model error input parameters. A percentile-based scaling method that addresses soil moisture biases in model and observational estimates opened the path to the successful application of land data assimilation to satellite retrievals of surface soil moisture. Assimilation of AMSR-E surface soil moisture retrievals into the NASA Catchment model provided superior surface and root zone assimilation products (when validated against in situ measurements and compared to the model estimates or satellite observations alone). The multi-model capabilities of US were used to investigate the role of subsurface physics in the assimilation of surface soil moisture observations. Results indicate that the potential of surface soil moisture assimilation to improve root zone information is higher when the surface to root zone coupling is stronger. Building on this experience, GSFC leads the development of the Level 4 Surface and Root-Zone Soil Moisture (L4_SM) product for the planned NASA Soil-Moisture-Active-Passive (SMAP) mission. A key milestone was the design and execution of an Observing System Simulation Experiment that quantified the contribution of soil moisture retrievals to land data assimilation products as a function of retrieval and land model skill and yielded an estimate of the error budget for the SMAP L4_SM product. Terrestrial water storage observations from GRACE satellite system were also successfully assimilated into the NASA Catchment model and provided improved estimates of groundwater variability when compared to the model estimates alone. Moreover, satellite-based land surface temperature (LST) observations from the ISCCP archive were assimilated using a bias estimation module that was specifically designed for LST assimilation. As with soil moisture, LST assimilation provides modest yet statistically significant improvements when compared to the model or satellite observations alone. To achieve the improvement, however, the LST assimilation algorithm must be adapted to the specific formulation of LST in the land model. An improved method for the assimilation of snow cover observations was also developed. Finally, the coupling of LIS to the mesoscale Weather Research and Forecasting (WRF) model enabled investigations into how the sensitivity of land-atmosphere interactions to the specific choice of planetary boundary layer scheme and land surface model varies across surface moisture regimes, and how it can be quantified and evaluated against observations. The on-going development and integration of land assimilation modules into the Land Information System will enable the use of GSFC software with a variety of land models and make it accessible to the research community.
Reichle, Rolf↗
A numerical method for acoustic oscillations in tubes
A numerical method to obtain the neutral curve for the onset of acoustic oscillations in a helium-filled tube is described. Such oscillations can cause a serious heat loss in the plumbing associated with liquid helium dewars. The problem is modelled by a second-order, ordinary differential eigenvalue problem for the pressure perturbation. The numerical method to find the eigenvalues and track the resulting points along the neutral curve is tailored to this problem. The results show that a tube with a uniform temperature gradient along it is much more stable than one where the temperature suddenly jumps from the cold to the hot value in the middle of the tube.
Gary, John M.↗
Initiating Molecular Growth in the Interstellar Medium via Dimeric Complexes of Observed Ions and Molecules
A feasible initiation step for particle growth in the interstellar medium (ISM) is simulated by means of ab quantum chemistry methods. The systems studied are dimer ions formed by pairing nitrogen containing small molecules known to exist in the ISM with ions of unsaturated hydrocarbons or vice versa. Complexation energies, structures of ensuing complexes and electronic excitation spectra of the encounter complexes are estimated using various quantum chemistry methods. Moller-Plesset perturbation theory (MP2, Z-averaged perturbation theory (ZAP2), coupled cluster singles and doubles with perturbative triples corrections (CCSD(T)), and density functional theory (DFT) methods (B3LYP) were employed along with the correlation consistent cc-pVTZ and aug-cc-pVTZ basis sets. Two types of complexes are predicted. One type of complex has electrostatic binding with moderate (7-20 kcal per mol) binding energies, that are nonetheless significantly stronger than typical van der Waals interactions between molecules of this size. The other type of complex develops strong covalent bonds between the fragments. Cyclic isomers of the nitrogen containing complexes are produced very easily by ion-molecule reactions. Some of these complexes show intense ultraviolet visible spectra for electronic transitions with large oscillator strengths at the B3LYP, omegaB97, and equations of motion coupled cluster (EOM-CCSD) levels. The open shell nitrogen containing carbonaceous complexes especially exhibit a large oscillator strength electronic transition in the visible region of the electromagnetic spectrum.
Bera, Partha P.↗
Nonlinear dynamics of drops and bubbles and chaotic phenomena
Nonlinear phenomena associated with the dynamics of free drops and bubbles are investigated analytically, numerically and experimentally. Although newly developed levitation and measurement techniques have been implemented, the full experimental validation of theoretical predictions has been hindered by interfering artifacts associated with levitation in the Earth gravitational field. The low gravity environment of orbital space flight has been shown to provide a more quiescent environment which can be utilized to better match the idealized theoretical conditions. The research effort described in this paper is a closely coupled collaboration between predictive and guiding theoretical activities and a unique experimental program involving the ultrasonic and electrostatic levitation of single droplets and bubbles. The goal is to develop and to validate methods based on nonlinear dynamics for the understanding of the large amplitude oscillatory response of single drops and bubbles to both isotropic and asymmetric pressure stimuli. The first specific area on interest has been the resonant coupling between volume and shape oscillatory modes isolated gas or vapor bubbles in a liquid host. The result of multiple time-scale asymptotic treatment, combined with domain perturbation and bifurcation methods, has been the prediction of resonant and near-resonant coupling between volume and shape modes leading to stable as well as chaotic oscillations. Experimental investigations of the large amplitude shape oscillation modes of centimeter-size single bubbles trapped in water at 1 G and under reduced hydrostatic pressure, have suggested the possibility of a low gravity experiment to study the direct coupling between these low frequency shape modes and the volume pulsation, sound-radiating mode. The second subject of interest has involved numerical modeling, using the boundary integral method, of the large amplitude shape oscillations of charged and uncharged drops in the presence of a static or time-varying electric field. Theoretically predicted non linearity in the resonant frequency of the fundamental quadrupole mode has been verified by the accompanying experimental studies. Additional phenomena such as hysteresis in the frequency response of ultrasoncially levitated droplets in the presence of a time varying electric field, and mode coupling in the oscillations of ultrasonically modulated droplets, have also been uncovered. One of the results of this ground-based research has been the identification and characterization of phenomena strictly associated with the influence of the gravitational field. This has also allowed us to identify the specific requirements for potential microgravity investigations yielding new information not obtainable on Earth.
Trinh, Eugene H.↗
Methods of determining electron concentrations in the magnetosphere from nose whistlers
Whistler propagation in the magnetosphere was studied in detail to find accurate and economical means of determining the path latitude and the electron concentration along the path from whistler parameters of nose frequency and travel time at the nose. Longitudinal propagation in field aligned whistler ducts of cold plasma was assumed, and the earth's magnetic field was approximated by a centered dipole. The effects of whistler propagation in the earth-ionosphere waveguide and through the conjugate ionospheres were treated as small perturbations. Several alternative methods are described so that the most economical method may be chosen depending on the desired accuracy and the availability of a computer or a calculator.
Park, C. G.↗
Guidance, flight mechanics and trajectory optimization. Volume 6 - The N-body problem and special perturbation techniques
Analytical formulations and numerical integration methods for many body problem and special perturbative techniques
Allione, M. S.↗
Stability of nonuniform rotor blades in hover using a mixed formulation
A mixed formulation for calculating static equilibrium and stability eigenvalues of nonuniform rotor blades in hover is presented. The static equilibrium equations are nonlinear and are solved by an accurate and efficient collocation method. The linearized perturbation equations are solved by a one step, second order integration scheme. The numerical results correlate very well with published results from a nearly identical stability analysis based on a displacement formulation. Slight differences in the results are traced to terms in the equations that relate moments to derivatives of rotations. With the present ordering scheme, in which terms of the order of squares of rotations are neglected with respect to unity, it is not possible to achieve completely equivalent models based on mixed and displacement formulations. The one step methods reveal that a second order Taylor expansion is necessary to achieve good convergence for nonuniform rotating blades. Numerical results for a hypothetical nonuniform blade, including the nonlinear static equilibrium solution, were obtained with no more effort or computer time than that required for a uniform blade.
Stephens, W. B.↗
Meso-beta scale numerical simulation studies of terrain-induced jet streak mass and momentum perturbations
An in-depth analysis of observed gravity waves and their relationship to precipitation bands over the Montana mesonetwork during the 11-12 July 1981 CCOPE case study indicated two episodes of coherent waves. While geostrophic adjustment, shearing instability, and terrain were all implicated separately or in combination as possible wave generation mechanisms, the lack of upper-air data within the wave genesis region made it difficult to define the genesis processes from observations alone. The first part of this paper, 3D Numerical Modeling Studies of Terrain-Induced Mass/Momentum Perturbations, employs a mesoscale numerical model to help diagnose the intricate early wave generation mechanisms during the first observed gravity wave episode. The meso-beta scale numerical model is used to study various simulations of the role of multiple geostrophic adjustment processes in focusing a region for gravity wave genesis. The second part of this paper, Linear Theory and Theoretical Modeling, investigates the response of non-resting rotating homogeneous and continuously stratified Boussinesq models of the terrestrial atmosphere to temporally impulsive and uniformly propagating three-dimensional localized zonal momentum sources representative of midlatitude jet streaks. The methods of linear perturbation theory applied to the potential vorticity (PV) and wave field equations are used to study the geostrophic adjustment dynamics. The total zonal and meridional wind perturbations are separated into geostrophic and ageostrophic components in order to define and follow the evolution of both the primary and secondary mesocirculations accompanying midlatitude jetogenesis forced by geostrophic adjustment processes. This problem is addressed to help fill the gap in understanding the dynamics and structure of mesoscale inertia-gravity waves forced by geostrophic adjustment processes in simple two-dimensional quiescent current systems and those produced by mesoscale numerical models simulating the orographic and diabatic perturbation of three-dimensional quasi-geostrophically balanced synoptic scale jet streaks associated with complex baroclinic severe storm producing environments.
Lin, Yuh-Lang↗
Application of MACSYMA to first order perturbation theory in celestial mechanics
The application of MACSYMA to general first order perturbation theory in celestial mechanics is explored. Methods of derivation of small variations in the Keplerian orbital elements are developed. As an example of the methods, the small general relativistic perturbations on the two-body Newtonian motion, resulting from the rotation of the central body, are developed in detail.
Anderson, J. D.↗
Solution of the minimum time-to-climb problem by matched asymptotic expansions
Application of singular perturbation techniques to trajectory optimization problems of flight mechanics is discussed. The method of matched asymptotic expansions is used to obtain an approximate solution to the aircraft minimum time-to-climb problem. Outer, boundary-layer, and composite solutions are obtained to zeroth and first orders. A stability criterion is derived for the zeroth-order boundary-layer solutions (the theory requires a form of boundary-layer stability). A numerical example is considered for which it is shown that the stability criterion is satisfied and a useful numerical solution is obtained. The zeroth-order solution proves to be a poor approximation, but the first-order solution gives a good approximation for both the trajectory and the minimum time-to-climb. The computational cost of the singular perturbation solution is considerably less than that of a steepest descent solution. Thus singular perturbation methods appear to be promising for the solution of optimal control problems.
Ardema, M. D.↗