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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,261 records · Page 70

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Transient Resonances in EMRIs: A Scalar Toy Model

We present calculations of the scalar self-force for a scalar point-charge following different $r\theta$-resonant geodesics in Kerr spacetime. Using this scalar model, we present two methods for constructing the self-force for resonant sources using two different phase parameterizations of the point-particle motion. We relate these two methods via a simple shifting relation that holds for both the scalar and gravitational self-forces. We discuss how this shifting relation allows us to compare self-force results both at and away from resonances. Finally, we provide new evidence that conservative scalar perturbations can drive the leading-order evolution of the orbital dynamics during resonances (i.e., the orbit-averaged rate of change of the Carter constant depends on contributions from the conservative self-force). This evidence conflicts with a commonly assumed conjecture that these integrated conservative contributions should vanish during both resonant and non-resonant motion.

Zachary Nasipak↗

NASA NTRS1972

Singular perturbations in the state regulator problem

Most of the results of singular perturbation theory have been concerned with initial value problems whereas optimal control problems are of two-point boundary value type. The portions of this theory applicable to the open loop state regulator problem are reviewed. For obtaining approximate solutions to the state regulator problem the method of matched asymptotic expansions is employed. This method has been developed in connection with certain fluid mechanics problems and is applicable to nonlinear as well as linear problems. It has been found in the past to be advantageous not to formulate this method generally but to apply it to each individual problem and this approach is adopted here. A general recipe for the method is given and its application is illustrated by using the method to obtain an approximate solution to a simple, specific state regulator problem.

Ardema, M. D.↗

NASA NTRS1976

Matrix perturbation for structural dynamic analysis

The objective of the effort was to investigate methodologies to reduce the cost of evaluating changes in dynamic loads when small modifications are made in the structure. A matrix perturbation technique has been developed to calculate the dynamic responses of a structural system that has been modified from the original design. The calculation is based on the results of the original design and the assumption that the structural modification is small. The advantage of the method is an update of the dynamic response due to design changes without performing an entire analysis. This procedure can be used in a design load analysis cycle in which the structural design is subject to frequent changes. A sample problem is given to demonstrate the validity of the technique.

Chen, J. C.↗

NASA NTRS1982

Solution of Von Karman's plate equation with perturbation and series summation

A nonliner series summation technique was combined with Chien's small parameter perturbation technique to solve the problem of the finite deflection of a simply supported circular plate subjected to a uniform transverse load. Problems involving finite deflections of plates with simply supported or free edges were regarded as almost intractable. The great advantage of incorporating the summation method into the solution procedure is demonstrated.

Watson, L. G.↗

NASA NTRS1990

The role of secondary resonances in the orbital history of Miranda

A simple 'perturbed pendulum' model is developed in order to ascertain the origin and dynamics of the secondary resonances in the tidal evolution of Miranda and Umbriel through the 1:3 mean motion commensurability. Long-term evolution is addressed on the basis of the computation-accelerating method of algebraic mappings. It is found that the present 4.34-deg value of Miranda's orbital inclination is probably a result of capture into the primary resonance, as well as subsequent capture into the 3/1 secondary resonance. The damping of an eccentricity as high as 0.035 by tidal dissipation in the satellite may account for the surface features of Miranda.

Malhotra, Renu↗

NASA NTRS1992

An extensive ab initio study of the structures, vibrational spectra, quadratic force fields, and relative energetics of three isomers of Cl2O2

Quantum mechanical computational methods are employed for an ab initio investigation of: (1) the molecular properties of the lowest isomers of the ClO dimer; and (2) predicted molecular and thermochemical properties. Techniques employed include electron correlation and particularly singles and doubles coupled-cluster (CCSD) theory with or without perturbational estimates of the effects of connected triple excitations. The isomers ClOClO and ClClO2 are found to have higher energies than the ClOOCl isomer, and the theoretical vibrational frequencies of the isomers are well correlated with experimental data. Experimental values of the heat of formation for the isomers are also compared with calculations based on an isodesmic reaction with Cl2O, H2O, and HOOH.

Lee, Timothy J.↗

NASA NTRS1996

Near-Optimal Operation of Dual-Fuel Launch Vehicles

A near-optimal guidance law for the ascent trajectory from earth surface to earth orbit of a fully reusable single-stage-to-orbit pure rocket launch vehicle is derived. Of interest are both the optimal operation of the propulsion system and the optimal flight path. A methodology is developed to investigate the optimal throttle switching of dual-fuel engines. The method is based on selecting propulsion system modes and parameters that maximize a certain performance function. This function is derived from consideration of the energy-state model of the aircraft equations of motion. Because the density of liquid hydrogen is relatively low, the sensitivity of perturbations in volume need to be taken into consideration as well as weight sensitivity. The cost functional is a weighted sum of fuel mass and volume; the weighting factor is chosen to minimize vehicle empty weight for a given payload mass and volume in orbit.

Ardema, M. D.↗

NASA NTRS2016

C/NOFS Remote Sensing of Ionospheric Reflectance

Alfvn waves play critical roles in the electrodynamic coupling of plasmas at magnetically conjugate regions in near-Earth space. Associated electric (E*) and magnetic (dec B*) field perturbations sampled by sensors on satellites in low-Earth orbits are generally super positions of incident and reflected waves. However, lack of knowledge about ionospheric reflection coefficients (alpha) hinders understanding of generator outputs and load absorption of Alfvn wave energies. Here we demonstrate a new method for estimating using satellite measurements of ambient E* and dec B* then apply it to a case in which the Communication Navigation Outage Forecasting System (CNOFS) satellite flew conjugate to the field of view of a 630.0 nm all-sky imager at El Leoncito, Argentina, while medium-scale traveling ionosphere disturbances were detected in its field of view. In regions of relatively large amplitudes of E* and B*,calculated values ranged between 0.67 and 0.88. This implies that due to impedance mismatches, the generator ionosphere puts out significantly more electromagnetic energy than the load can absorb. Our analysis also uncovered caveats concerning the methods range of applicability in regions of low E* and B*. The method can be validated in future satellite-based auroral studies where energetic particle precipitation fluxes can be used to make independent estimates of alpha.

Burke, W. J.↗

NASA NTRS2017

Low-Thrust Trajectory Optimization with Simplified SQP Algorithm

The problem of low-thrust trajectory optimization in highly perturbed dynamics is a stressing case for many optimization tools. Highly nonlinear dynamics and continuous thrust are each, separately, non-trivial problems in the field of optimal control, and when combined, the problem is even more difficult. This paper de-scribes a fast, robust method to design a trajectory in the CRTBP (circular restricted three body problem), beginning with no or very little knowledge of the system. The approach is inspired by the SQP (sequential quadratic programming) algorithm, in which a general nonlinear programming problem is solved via a sequence of quadratic problems. A few key simplifications make the algorithm presented fast and robust to initial guess: a quadratic cost function, neglecting the line search step when the solution is known to be far away, judicious use of end-point constraints, and mesh refinement on multiple shooting with fixed-step integration.In comparison to the traditional approach of plugging the problem into a “black-box” NLP solver, the methods shown converge even when given no knowledge of the solution at all. It was found that the only piece of information that the user needs to provide is a rough guess for the time of flight, as the transfer time guess will dictate which set of local solutions the algorithm could converge on. This robustness to initial guess is a compelling feature, as three-body orbit transfers are challenging to design with intuition alone. Of course, if a high-quality initial guess is available, the methods shown are still valid.We have shown that endpoints can be efficiently constrained to lie on 3-body repeating orbits, and that time of flight can be optimized as well. When optimizing the endpoints, we must make a trade between converging quickly on sub-optimal endpoints or converging more slowly on end-points that are arbitrarily close to optimal. It is easy for the mission design engineer to adjust this trade based on the problem at hand.The biggest limitation to the algorithm at this point is that multi-revolution transfers (greater than 2 revolutions) do not work nearly as well. This restriction comes in because the relationship between node 1 and node N becomes increasingly nonlinear as the angular distance grows. Trans-fers with more than about 1.5 complete revolutions generally require the line search to improve convergence. Future work includes: Comparison of this algorithm with other established tools; improvements to how multiple-revolution transfers are handled; parallelization of the Jacobian computation; in-creased efficiency for the line search; and optimization of many more trajectories between a variety of 3-body orbits.

Parrish, Nathan L.↗

NASA NTRS2015

Methods and Apparatuses for Signaling with Geometric Constellations in a Raleigh Fading Channel

Communication systems are described that use signal constellations, which have unequally spaced (i.e., `geometrically` shaped) points. In many embodiments, the communication systems use specific geometric constellations that are capacity optimized at a specific SNR (signal to noise ratio). In addition, ranges within which the constellation points of a capacity optimized constellation can be perturbed and are still likely to achieve a given percentage of the optimal capacity increase compared to a constellation that maximizes d (sub min) (i.e. minimum distance between constellations) are also described. Capacity measures that are used in the selection of the location of constellation points include, but are not limited to, parallel decode (PD) capacity and joint capacity.

Barsoum, Maged F.↗

NASA NTRSDate not supplied

Metal Hydroxide and Oxyhydroxide Gas Database

In a high temperature combustion situation such as a gas turbine or rocket engine, combustion products include water vapor. This water vapor reacts with oxides in the hot-gas-path to form volatile hydroxides and oxyhydroxides. The oxides may be structural oxide components or protective oxides on alloys or non-oxide ceramics. In either case, this enhanced volatilization leads to degradation of components. Reliable thermochemical data for many metal hydroxides and oxyhydroxides has not been available. We review a database put together at NASA over many years based on experimental and ab initio computations for thermochemical data of Al, Si, Ti, Cr, Mn, Y, Zr, La, Ta, Gd, and Yb hydroxides and oxyhydroxides. Experimental methods for measuring thermodynamic quantities of metal hydroxides and oxyhydroxides are limited since the required high temperatures and high pressures preclude many conventional methods. Nonetheless, transpiration has proven a valuable tool in our laboratories for the Si-O-H, Cr-O-H, and Ti-O-H systems. Data are analyzed with second and third law methods. These data can be leveraged with quantum chemistry to provide reliable free energy functions for the third law calculation and comparison to the derived enthalpies of formation. Both wave function and DFT methods can be utilized for deriving the thermodynamic functions for the various metal hydroxides and oxyhydroxides. We have found that, in general, reasonable molecular geometries and spectroscopic data can be derived with a relatively low level of theory, i.e. DFT with the B3LYP functional. Accurate calculations of the enthalpy of formation require higher levels of theory, such as the Coupled Cluster, Singles, Doubles, and perturbative Triples CCSD(T) approach with complete basis sets (CBS). Also, corrections for hindered rotations and anharmonicity are necessary. For some symmetries, the anharmonic correction has proven to be a challenge and we have found that this is only treated properly by certain codes. The database lists the fHo(298), So(298), and Cp in the form Cp = A + BT + CT2 + D/T + E/T2, where A, B, C, D, and E are the fitting constants. Thus, it can be easily used in any of the common computational thermochemistry codes. To date, the database contains 51 gaseous species of the formula MOx(OH)y.

Metal Hydroxides↗

NASA NTRS2017

Methods and Apparatuses for Signaling with Geometric Constellations in a Raleigh Fading Channel

Communication systems are described that use signal constellations, which have unequally spaced (i.e. `geometrically` shaped) points. In many embodiments, the communication systems use specific geometric constellations that are capacity optimized at a specific SNR, over the Raleigh fading channel. In addition, ranges within which the constellation points of a capacity optimized constellation can be perturbed and are still likely to achieve a given percentage of the optimal capacity increase compared to a constellation that maximizes d.sub.min, are also described. Capacity measures that are used in the selection of the location of constellation points include, but are not limited to, parallel decode (PD) capacity and joint capacity.

Barsoum, Maged F.↗

NASA NTRS1981

A family of solution algorithms for nonlinear structural analysis based on relaxation equations

A family of hierarchical algorithms for nonlinear structural equations are presented. The algorithms are based on the Davidenko-Branin type homotopy and shown to yield consistent hierarchical perturbation equations. The algorithms appear to be particularly suitable to problems involving bifurcation and limit point calculations. An important by-product of the algorithms is that it provides a systematic and economical means for computing the stepsize at each iteration stage when a Newton-like method is employed to solve the systems of equations. Some sample problems are provided to illustrate the characteristics of the algorithms.

Park, K. C.↗

NASA NTRS1988

Linear statistical analysis of maneuver optimization techniques

It is shown how a method for maneuver design during the Jovian satellite tour phase of the Galileo mission requires a full-scale redesign of the tour at each maneuver. Software was developed to provide insight into the statistical properties of this process; it is based on the use of linear perturbations to provide trajectory propagation and a least squares representation of the optimization. The result is a Monte Carlo program, making it possible to evaluate the proposed strategy.

Maize, Earl H.↗

NASA NTRS1989

Wentzel-Kramers-Brillouin method in the Bargmann representation

It is demonstrated that the Bargmann representation of quantum mechanics is ideally suited for semiclassical analysis, using as an example the WKB method applied to the bound-state problem in a single well of one degree of freedom. For the harmonic oscillator, this WKB method trivially gives the exact eigenfunctions in addition to the exact eigenvalues. For an anharmonic well, a self-consistent variational choice of the representation greatly improves the accuracy of the semiclassical ground state. Also, a simple change of scale illuminates the relationship of semiclassical versus linear perturbative expansions, allowing a variety of multidimensional extensions.

Voros, A.↗

NASA NTRS1992

A numerical study of two-dimensional vortex shedding from rectangular cylinders

An efficient time-marching, non-iterative calculation method is used to analyze time-dependent flows around rectangular cylinders. The turbulent flow in the wake region of a square section cylinder is analyzed using an anisotropic k-epsilon model. Initiation and subsequent development of the vortex shedding phenomenon is naturally captured once a perturbation is introduced in the flow. Transient calculations using standard eddy-viscosity and an anisotropic k-epsilon model averaged over an integral number of cycles to get the fluctuating energy (organized and turbulent) are compared with experimental data. It is shown that the anisotropic k-epsilon model resolves the anisotropy of the Reynolds stresses and gives mean energy distribution closer to the experiment than the standard k-epsilon model.

Hadid, A. H.↗

NASA NTRS1996

Investigation of the Interaction of External Disturbances with Roughened Flat Plate Boundary Layer

The interaction of an external disturbance with a laminar boundary layer over a flat plate with distributed roughness is investigated using combined experimental and numerical methods. The experiment is modeled with an unsteady boundary layer code using second order backward differencing. The simulation includes the second order scattering from roughness elements at and near the first streamwise station of predicted boundary layer instability. A comparison of experimental measurements of the boundary layer perturbation due to the wake from a vibrating ribbon with the computed first order forced boundary layer perturbation showed excellent agreement. Second order roughness induced eigenfunctions from boundary layer theory are examined and compared with other forms of excitation

Davis, Sanford S.↗

NASA NTRS2012

Helioseismology of a Realistic Magnetoconvective Sunspot Simulation

We compare helioseismic travel-time shifts measured from a realistic magnetoconvective sunspot simulation using both helioseismic holography and time-distance helioseismology, and measured from real sunspots observed with the Helioseismic and Magnetic Imager instrument on board the Solar Dynamics Observatory and the Michelson Doppler Imager instrument on board the Solar and Heliospheric Observatory. We find remarkable similarities in the travel-time shifts measured between the methodologies applied and between the simulated and real sunspots. Forward modeling of the travel-time shifts using either Born or ray approximation kernels and the sound-speed perturbations present in the simulation indicates major disagreements with the measured travel-time shifts. These findings do not substantially change with the application of a correction for the reduction of wave amplitudes in the simulated and real sunspots. Overall, our findings demonstrate the need for new methods for inferring the subsurface structure of sunspots through helioseismic inversions.

Sun-interior↗