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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1971

An analysis and comparison of several trajectory optimization methods

The sensitivities of the convergence characteristics of the methods to initially assumed parameters and trial solution, convergence times, computer logic, and storage requirements are discussed. Numerical comparison of the convergence characteristics is made by considering a minimum time, low thrust, Earth-Mars transfer trajectory. A modified quasi-linearization method reduces convergence time by approximately 70% when compared with the generalized Newton-Raphson method and allows the terminal boundary to be specified by a general function of the problem variables. A uniquely specified and easily determined, time dependent weighting matrix for the gradient techniques accelerates the shaping of the optimal control program and improves the convergence characteristics during the terminal iterations. Convergence envelopes, indicating how sensitive the convergence characteristics are to initially assumed parameters, are plotted for the perturbation and quasi-linearization methods. Several iteration schemes are proposed which increase the size of the convergence envelopes and decrease the sensitivity of the method to initially assumed parameters.

Lewallen, J. M.↗
NASA NTRS1990

Singular perturbation analysis of AOTV-related trajectory optimization problems

The problem of real time guidance and optimal control of Aeroassisted Orbit Transfer Vehicles (AOTV's) was addressed using singular perturbation theory as an underlying method of analysis. Trajectories were optimized with the objective of minimum energy expenditure in the atmospheric phase of the maneuver. Two major problem areas were addressed: optimal reentry, and synergetic plane change with aeroglide. For the reentry problem, several reduced order models were analyzed with the objective of optimal changes in heading with minimum energy loss. It was demonstrated that a further model order reduction to a single state model is possible through the application of singular perturbation theory. The optimal solution for the reduced problem defines an optimal altitude profile dependent on the current energy level of the vehicle. A separate boundary layer analysis is used to account for altitude and flight path angle dynamics, and to obtain lift and bank angle control solutions. By considering alternative approximations to solve the boundary layer problem, three guidance laws were derived, each having an analytic feedback form. The guidance laws were evaluated using a Maneuvering Reentry Research Vehicle model and all three laws were found to be near optimal. For the problem of synergetic plane change with aeroglide, a difficult terminal boundary layer control problem arises which to date is found to be analytically intractable. Thus a predictive/corrective solution was developed to satisfy the terminal constraints on altitude and flight path angle. A composite guidance solution was obtained by combining the optimal reentry solution with the predictive/corrective guidance method. Numerical comparisons with the corresponding optimal trajectory solutions show that the resulting performance is very close to optimal. An attempt was made to obtain numerically optimized trajectories for the case where heating rate is constrained. A first order state variable inequality constraint was imposed on the full order AOTV point mass equations of motion, using a simple aerodynamic heating rate model.

Calise, Anthony J.↗
NASA NTRSDate not supplied

Tropical Tropospheric Ozone Trends (1990 to 2022): A Re-evaluation Based on SHADOZ and IAGOS Profiles and TOMS/OMI Columns

Changes in tropical tropospheric ozone (TTO) are of importance because this region spans roughly a third of the Earth and portions of it are experiencing variability in trends of ozone precursors (CO, NO x , CH 4 and nonmethane hydrocarbons) associated with economic growth and fires. In addition to ozone changes affecting radiative forcing, tropical ozone is an important source of the OH radical and thus, the oxidizing capacity of the planet (Thompson, 1992). Recent studies examining TTO trends satellite and in-situ observations over the past ~25 years include: Thompson et al., JGR, 2021; Gaudel et al., ACP, 2023; Stauffer et al., ACP, 2023. The results show considerable regional and seasonal variability in TTO trends and sensitivity to data selection, frequency, and statistical method used. The satellite data vary most widely in method, time period and reliability. Here we revisit trends for the 1990-2022 period with the best-characterized buv-based satellite products that span that period (derived from TOMS and OMI, Ziemke et al., 2019). The satellite-based trends are compared to trends based on in-situ data from the Southern Hemisphere Additional Ozonesondes (SHADOZ) network (1998-2022), measurements from selected pre-SHADOZ and IAGOS commercial aircraft data (Gaudel et al., 2023). Among sensitivities examined are the dependence of ozone trends on start and end years, impacts of ENSO events and the COVID-19 perturbation to emissions. Two statistical methods are used, quantile regression (QR) and multiple linear regression (MLR). Trends of total TTO, ozone segments in the boundary layer (to ~700 hPa), and free troposphere (700-300 hPa) are compared.

ozone, OMI, tropospheric ozone, SHADOZ↗
NASA NTRS1980

Numerical solutions of transonic flows by parametric differentiation and integral equation techniques

The paper discusses results of an exploratory study of the advantages obtained by combining integral equation methods with the method of parametric differentiation in the treatment of transonic flow problems. In the proposed method, the nonlinear unsteady transonic flow equation for small perturbations is transformed into a linear equation by the use of the method of parametric differentiation. The linear equation is split into a pair of weakly coupled partial differential equations by writing the transformed perturbation potential as the sum of a steady component and an unsteady component. The solution of the steady equation as an integral equation is based on Ogana's treatment (1978). As a test case, the formulation is applied to predict the steady transonic flow over a nonlifting parabolic-arc airfoil.

Sivaneri, N. T.↗
NASA NTRS2003

Prospects for an Improved Measurement of Experimental Limit on G-dot

The orbital motion of an ultra-drag-free satellite, such as the large test body of the SEE (Satellite Energy Exchange) satellite, known as the "Shepherd," may possibly provide the best test for time variation of the gravitational constant G at the level of parts in 10(exp 14). Scarcely anything could be more significant scientifically than the incontestable discovery that a fundamental "constant" of Nature is not constant. A finding of non-zero (G-dot)/G would clearly mark the boundaries where general relativity is valid, and specify the onset of new physics. The requirements for measuring G-dot at the level proposed by SEE will require great care in treating perturbation forces. In the present paper we concentrate on the methods for dealing with the gravitational field due to possible large manufacturing defects in the SEE observatory. We find that, with adequate modeling of the perturbation forces and cancellation methods, the effective time-averaged acceleration on the SEE Shepherd will be approx. 10(exp -18) g (10(exp -17) m/sq s).

Sanders, Alvin J.↗
NASA NTRS1989

On equivalence of discrete-discrete and continuum-discrete design sensitivity analysis

Developments in design sensitivity analysis (DSA) method have been made using two fundamentally different approaches as shown. In the first approach, a discretized structural finite element model is used to carry out DSA. There are three different methods in the discrete DSA approach: finite difference, semi-analytical, and analytical methods. The finite difference method is a popular one due to its simplicity, but a serious shortcoming of the method is the uncertainty in the choice of a perturbation step size of design variables. In the semi-analytical method, the derivatives of stiffness matrix is computed by finite differences, whereas in the analytical method, the derivatives are obtained analytically. For the shape design variable, computation of analytical derivative of stiffness matrix is quite costly. Because of this, the semi-analytical method is a popular choice in discrete shape DSA approach. However, recently, Barthelemy and Haftka presented that the semi-analytical method can have serious accuracy problems for shape design variables in structures modeled by beam, plate, truss, frame, and solid elements. They found that accuracy problems occur even for a simple cantilever beam. In the second approach, a continuum model of the structure is used to carry out DSA.

Choi, Kyung K.↗
NASA NTRS1999

Time-of-Arrival Lightning Location Retrieval Using an Oblate Spheroidal Earth Model

The problem of retrieving lightning ground strike location on an oblate spheroidal Earth using a network of 4 or more time-of-arrival sensors is considered. A recently developed analytic method for obtaining such retrievals on a spherical Earth surface is perturbed resulting in an iterative procedure to get correction terms. The perturbation procedure consists of applying a vector Newton's method to eqs. relating the distances from the lightning location to each sensor along a geodesic and the times of arrival of the wave produced by the lightning source at each sensor.

Solakiewicz, Richard↗
NASA NTRS2014

Can Eccentric Debris Disks Be Long-lived? A First Numerical Investigation and Application to Zeta(exp 2) Reticuli

Context. Imaging of debris disks has found evidence for both eccentric and offset disks. One hypothesis is that they provide evidence for massive perturbers, for example, planets or binary companions, which sculpt the observed structures. One such disk was recently observed in the far-IR by the Herschel Space Observatory around Zeta2 Reticuli. In contrast with previously reported systems, the disk is significantly eccentric, and the system is several Gyr old. Aims. We aim to investigate the long-term evolution of eccentric structures in debris disks caused by a perturber on an eccentric orbit around the star. We hypothesise that the observed eccentric disk around Zeta2 Reticuli might be evidence of such a scenario. If so, we are able to constrain the mass and orbit of a potential perturber, either a giant planet or a binary companion. Methods. Analytical techniques were used to predict the effects of a perturber on a debris disk. Numerical N-body simulations were used to verify these results and further investigate the observable structures that may be produced by eccentric perturbers. The long-term evolution of the disk geometry was examined, with particular application to the Zeta2 Reticuli system. In addition, synthetic images of the disk were produced for direct comparison with Herschel observations. Results. We show that an eccentric companion can produce both the observed offsets and eccentric disks. These effects are not immediate, and we characterise the timescale required for the disk to develop to an eccentric state (and any spirals to vanish). For Zeta2 Reticuli, we derive limits on the mass and orbit of the companion required to produce the observations. Synthetic images show that the pattern observed around Zeta2 Reticuli can be produced by an eccentric disk seen close to edge-on, and allow us to bring additional constraints on the disk parameters of our model (disk flux and extent). Conclusions. We conclude that eccentric planets or stellar companions can induce long-lived eccentric structures in debris disks. Observations of such eccentric structures thus provide potential evidence of the presence of such a companion in a planetary system. We considered the specific example of Zeta2 Reticuli, whose observed eccentric disk can be explained by a distant companion (at tens of AU) on an eccentric orbit (ep greater than approx. 0.3).

planetary systems↗
NASA NTRS1974

Development of solution techniques for nonlinear structural analysis

Nonlinear structural solution methods in the current research literature are classified according to order of the solution scheme, and it is shown that the analytical tools for these methods are uniformly derivable by perturbation techniques. A new perturbation formulation is developed for treating an arbitrary nonlinear material, in terms of a finite-difference generated stress-strain expansion. Nonlinear geometric effects are included in an explicit manner by appropriate definition of an applicable strain tensor. A new finite-element pilot computer program PANES (Program for Analysis of Nonlinear Equilibrium and Stability) is presented for treatment of problems involving material and geometric nonlinearities, as well as certain forms on nonconservative loading.

Vos, R. G.↗
NASA NTRS1976

Identification of natural frequencies and modal damping ratios of aerospace structures from response data

An analysis of one and multidegree of freedom systems with classical damping is presented. Definition and minimization of error functions for each system are discussed. Systems with classical and nonclassical normal modes are studied, and results for first order perturbation are given. An alternative method of matching power spectral densities is provided, and numerical results are reviewed.

Michalopoulos, C. D.↗
NASA NTRS1979

Test problems for inviscid transonic flow

Solving of test problems with the TRANDES program is discussed. This method utilizes the full, inviscid, perturbation potential flow equation in a Cartesian grid system that is stretched to infinity. This equation is represented by a nonconservative system of finite difference equations that includes at supersonic points a rotated difference scheme and is solved by column relaxation. The solution usually starts from a zero perturbation potential on a very coarse grid (typically 13 by 7) followed by several grid halvings until a final solution is obtained on a fine grid (97 by 49).

Carlson, L. A.↗
NASA NTRS1981

The supersonic fan engine - An advanced concept in supersonic cruise propulsion

Engine performance and mission studies were conducted for a novel turbofan engine concept incorporating a supersonic through-flow fan, and comparisons were made with two supersonic transport (SST) engine concepts of equivalent thrust and technological sophistication. It was found that in the case of an SST with a cruise speed of Mach 2.32, the through-flow fan engine may yield ranges 10 to 20% greater than the two alternatives considered. The engine has a conventional core, with the supersonic fan being driven by a concentric low-pressure turbine that is uncoupled with the single, high pressure turbine/compressor core spool. Among the topics discussed are the methods of analysis employed and perturbation studies concerning supersonic fan adiabatic efficiency, fan discharge characteristics and propulsion system weight.

Franciscus, L. C.↗
NASA NTRS1989

Computing induced velocity perturbations due to a helicopter fuselage in a free stream

The velocity field of a representative helicopter fuselage in a free stream is computed. Perturbation velocities due to the fuselage are computed in a plan above the location of the helicopter rotor (rotor removed). The velocity perturbations computed by a source-panel model of the fuselage are compared with experimental measurements taken with a laser velocimeter. Three paneled fuselage models are studied: fuselage shape, fuselage shape with hub shape, and a body of revolution. The velocity perturbations computed for both fuselage shape models agree well with the measured velocity field except in the close vicinity of the rotor hub. In the hub region, without knowing the extent of separation, modeling of the effective source shape is difficult. The effects of the fuselage perturbations are not well-predicted with a simplified ellipsoid fuselage. The velocity perturbations due to the fuselage at the plane of the measurements have magnitudes of less than 8 percent of free-stream velocity. The velocity perturbations computed by the panel method are tabulated for the same locations at which previously reported rotor-inflow velocity measurements were made.

Berry, John D.↗
NASA NTRS1998

A Nonlinear Physics-Based Optimal Control Method for Magnetostrictive Actuators

This paper addresses the development of a nonlinear optimal control methodology for magnetostrictive actuators. At moderate to high drive levels, the output from these actuators is highly nonlinear and contains significant magnetic and magnetomechanical hysteresis. These dynamics must be accommodated by models and control laws to utilize the full capabilities of the actuators. A characterization based upon ferromagnetic mean field theory provides a model which accurately quantifies both transient and steady state actuator dynamics under a variety of operating conditions. The control method consists of a linear perturbation feedback law used in combination with an optimal open loop nonlinear control. The nonlinear control incorporates the hysteresis and nonlinearities inherent to the transducer and can be computed offline. The feedback control is constructed through linearization of the perturbed system about the optimal system and is efficient for online implementation. As demonstrated through numerical examples, the combined hybrid control is robust and can be readily implemented in linear PDE-based structural models.

Smith, Ralph C.↗
NASA NTRS1995

A Coupled-Cluster Study of the Molecular Structure, Vibrational Spectrum, and Heats of Formation of XONO2 (X=H, F, Cl)

The equilibrium structures, harmonic vibrational frequencies, dipole moments, and IR intensities of nitric acid, fluorine nitrate, and chlorine nitrate have been investigated by using the singles and doubles coupled cluster method that also includes a perturbational estimate of the effects of connected triple excitations, CCSD(T). A standard triple-zeta double-polarized basis set was utilized. The equilibrium geometries and vibrational spectra of HONO2 and ClONO2 are shown to be in excellent agreement with the available experimental data. The ab initio vibrational spectrum of FONO2 is also shown to be in excellent agreement with experiment. Unlike the FOOF and FON molecules, but similar to the cis- and trans-FONO molecules, FONO2 is shown to possess normal bond distances. The bonding in FONO2 is shown to be more similar to that in ClONO2 than that in HONO2, although there are still significant differences, especially in the partial atomic charges as deduced from Mulliken populations. This causes FONO2 to possess almost no dipole moment, which is very different to both HONO2 and ClONO2. By using large atomic natural orbital basis sets, CCSD(T) energies are computed for four isodesmic reactions in order to determine an accurate heat of formation for FONO2. Our best estimate for DELTA H(raised circle) (sub f,298) is 3.1 plus or minus 2.0 kcal/mol, indicating that the F-ONO2 bond energy is 31.3 kcal/mol.

Lee, Timothy J.↗
NASA NTRS1998

Adaptation in a rotating artificial gravity environment

The centripetal force generated by a rotating space vehicle is a potential source of artificial gravity. Minimizing the cost of such a vehicle dictates using the smallest radius and highest rotation rate possible, but head movements made at high rotation rates generate disorienting, nauseogenic cross-coupled semicircular canal stimulation. Early studies suggested 3 or 4 rpm as the highest rate at which humans could adapt to this vestibular stimulus. These studies neglected the concomitant Coriolis force actions on the head/neck system. We assessed non-vestibular Coriolis effects by measuring arm and leg movements made in the center of a rotating room turning at 10 rpm and found that movement endpoints and trajectories are initially deviated; however, subjects readily adapt with 10-20 additional movements, even without seeing their errors. Equilibrium point theories of motor control errantly predict that Coriolis forces will not cause movement endpoint errors so that subjects will not have to adapt their reaching movements during rotation. Adaptation of movement trajectory acquired during Coriolis force perturbations of one arm transfers to the unexposed arm but there is no intermanual transfer of endpoint adaptation indicating that neuromotor representations of movement endpoint and trajectory are separable and can adapt independently, also contradictory to equilibrium point theories. Touching a surface at the end of reaching movements is required for complete endpoint adaptation in darkness but trajectory adapts completely with or without terminal contact. We have also made the first kinematic measurements of unconstrained head movements during rotation, these movements show rapid adaptation to Coriolis force perturbations. Our results point to methods for achieving full compensation for rotation up to 10 rpm. Copyright 1998 Published by Elsevier Science B.V.

Non-NASA Center↗