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,117 records · Page 62

End effects of a bifilar magnetic wiggler

The magnetic field at the entrance of a wiggler magnet in a free-electron laser (FEL) determines the character of the electron orbits and influences the FEL gain. To obtain the desired adiabatic entrance conditions, it is necessary to investigate different methods of connecting leads to the wiggler as well as methods of gradually reducing the wiggler field. The field of a current-driven, finite-length bifilar helical wiggler magnet is derived and is used to investigate the magnetic field end perturbations. For a limited range of wiggler parameters, flaring the wiggler end produces an adiabatic entrance region. Alternatively, gradually terminating the wiggler current at the wiggler entrance is shown to produce an adiabatic entrance for a large range of parameters.

Fajans, J.↗

Design sensitivity analysis with Applicon IFAD using the adjoint variable method

A numerical method is presented to implement structural design sensitivity analysis using the versatility and convenience of existing finite element structural analysis program and the theoretical foundation in structural design sensitivity analysis. Conventional design variables, such as thickness and cross-sectional areas, are considered. Structural performance functionals considered include compliance, displacement, and stress. It is shown that calculations can be carried out outside existing finite element codes, using postprocessing data only. That is, design sensitivity analysis software does not have to be imbedded in an existing finite element code. The finite element structural analysis program used in the implementation presented is IFAD. Feasibility of the method is shown through analysis of several problems, including built-up structures. Accurate design sensitivity results are obtained without the uncertainty of numerical accuracy associated with selection of a finite difference perturbation.

Frederick, Marjorie C.↗

Implementation of design sensitivity analysis with existing finite element codes

A numerical method is presented to implement structural design sensitivity analysis theory, using the versatility and convenience of existing finite element structural analysis programs. Design variables such as thickness and cross-sectional areas of components of individual members and built-up structures are considered. Structural performance functionals considered include displacement and stress. The method is also applicable for eigenvalue problem design sensitivity analysis. It is shown that calculations can be carried out outside existing finite element codes, using postprocessing data only. Thus, design sensitivity analysis software does not have to be embedded in an existing finite element code. Feasability of the method is shown through analysis of several problems, including a built-up structure. Accurate design sensitivity results are obtained without the uncertainty of numerical accuracy associated with selection of finite difference perturbations.

Choi, Kyung K.↗

Magnetic reversal spurts: Rain gauges for comet showers

Abrupt increases in the rate of magnetic reversals (magnetic reversal spurts) were first studied by many others. They hypothesized that spurts result from increased turbulence in the earth's core dynamo during episodes of intense bolide bombardment of the earth. Mechanisms for creating episodes of intense bombardment of the earth involve gravitational perturbation of the Oort cloud of comets, either by a hidden planet, a solar companion, or massive matter in the galactic plane. Herein, the time variation in reversal rate is analyzed using methods of statistical density estimation. A smooth, continuous estimate of reversal rate is obtained using an adaptive kernel method, in which the kernel width is adjusted as a function of reversal rate. The estimates near the ends of the data series (at 165 my ago and the present) are obtained by extending the data by reflection. The results show that the reversal spurts are not associated demonstrably with extinctions or well-dated impacts. If the spurts do record episodes of intense bombardment of the earth, then the mass extinctions do not, in general, occur at times of impacts. Furthermore, the large impact craters seen are not obviously related to the spurts, suggesting that the craters may have been caused by bolides of a different nature and with a different temporal pattern. However, the most simple explanation seems to be that the spurts do not record comet showers, either because the recording mechanism suggested by Muller and Morris is not effective or because comet showers are not triggered in the ways considered by Hut et al.

Lutz, T. M.↗

Direct numerical simulations of a spatially developing plane wake

In the present paper, direct numerical methods by which to simulate the spatially developing free shear flows in the transitional region are described and the numerical results of a spatially developing plane wake are presented. The incompressible time-dependent Navier-Stokes equations were solved using Pade finite difference approximations in the streamwise direction, a mapped pseudospectral Fourier method in the cross-stream direction, and a third-order compact Runge-Kutta scheme for time advancement. The unstable modes of the Orr-Sommerfeld equations were used to perturb the inlet of the wake. Statistical analyses were performed and some numerical results were compared with experimental measurements. When only the fundamental mode is forced, the energy spectra show amplification of the fundamental and its higher harmonics. In this case, unperturbed alternate vortices develop in the saturation region of the wake. The phase jitter around the fundamental frequency plays a critical role in generating vortices of random shape and spacing. Large- and small-scale distortions of the fundamental structure are observed. Pairing of vortices of the same sign is observed, as well as vortex coupling of vortices of the opposite sign.

Maekawa, Hiroshi↗

Bursting into space: alterations of sympathetic control by space travel

AIM: Astronauts return to Earth with reduced red cell masses and hypovolaemia. Not surprisingly, when they stand, their heart rates may speed inordinately, their blood pressures may fall, and some may experience frank syncope. We studied autonomic function in six male astronauts (average +/- SEM age: 40 +/- 2 years) before, during, and after the 16-day Neurolab space shuttle mission. METHOD: We recorded electrocardiograms, finger photoplethysmographic arterial pressures, respiration, peroneal nerve muscle sympathetic activity, plasma noradrenaline and noradrenaline kinetics, and cardiac output, and we calculated stroke volume and total peripheral resistance. We perturbed autonomic function before and during spaceflight with graded Valsalva manoeuvres and lower body suction, and before and after the mission with passive upright tilt. RESULTS: In-flight baseline sympathetic nerve activity was increased above pre-flight levels (by 10-33%) in three subjects, in whom noradrenaline spillover and clearance also were increased. Valsalva straining provoked greater reductions of arterial pressure, and proportionally greater sympathetic responses in space than on Earth. Lower body suction elicited greater increases of sympathetic nerve activity, plasma noradrenaline, and noradrenaline spillover in space than on Earth. After the Neurolab mission, left ventricular stroke volume was lower and heart rate was higher during tilt, than before spaceflight. No astronaut experienced orthostatic hypotension or pre-syncope during 10 min of post-flight tilting. CONCLUSION: We conclude that baseline sympathetic outflow, however measured, is higher in space than on earth, and that augmented sympathetic nerve responses to Valsalva straining, lower body suction, and post-flight upright tilt represent normal adjustments to greater haemodynamic stresses associated with hypovolaemia.

NASA Program Biomedical Research and Countermeasur↗

Aerosol-Induced Changes of Convective Cloud Anvils Produce Strong Climate Warming

The effect of aerosol on clouds poses one of the largest uncertainties in estimating the anthropogenic contribution to climate change. Small human-induced perturbations to cloud characteristics via aerosol pathways can create a change in the top-of-atmosphere radiative forcing of hundreds of Wm(exp-2) . Here we focus on links between aerosol and deep convective clouds of the Atlantic and Pacific Intertropical Convergence Zones, noting that the aerosol environment in each region is entirely different. The tops of these vertically developed clouds consisting of mostly ice can reach high levels of the atmosphere, overshooting the lower stratosphere and reaching altitudes greater than 16 km. We show a link between aerosol, clouds and the free atmosphere wind profile that can change the magnitude and sign of the overall climate radiative forcing. We find that increased aerosol loading is associated with taller cloud towers and anvils. The taller clouds reach levels of enhanced wind speeds that act to spread and thin the anvi1 clouds, increasing areal coverage and decreasing cloud optical depth. The radiative effect of this transition is to create a positive radiative forcing (warming) at top-of-atmosphere. Furthermore we introduce the cloud optical depth (r), cloud height (Z) forcing space and show that underestimation of radiative forcing is likely to occur in cases of non homogenous clouds. Specifically, the mean radiative forcing of towers and anvils in the same scene can be several times greater than simply calculating the forcing from the mean cloud optical depth in the scene. Limitations of the method are discussed, alternative sources of aerosol loading are tested and meteorological variance is restricted, but the trend of taller clouds; increased and thinner anvils associated with increased aerosol loading remains robust through all the different tests and perturbations.

Koren, I.↗

Two substorm studies of relations between westward electric fields in the outer plasmasphere, auroral activity, and geomagnetic perturbations

Temporal variations of the westward component of the magnetospheric convection electric field in the outer plasmasphere were compared to auroral activity near L = 7, and to variations in the geomagnetic field at middle and high latitudes. The substorms occurred on July 29, 1965 near 0530 UT and on August 20, 1965 near 0730 UT. The results on westward electric field E(w) were obtained by the whistler method using data from Eights, Antarctica (L is approximately 4). All sky camera records were obtained from Byrd, Antarctica, (L is approximately 7), located within about 1 hour of Eights in magnetic local time. It was found that E(w) within the outer plasmasphere increased rapidly to substorm levels about the time of auroral expansion at nearby longitudes. This behavior is shown to differ from results on E(w) from balloons, which show E(w) reaching enhanced levels prior to the expansion. A close temporal relation was found between the rapid, substorm associated increases in E(w) and a well known type of nightside geomagnetic perturbation. Particularly well defined was the correlation of E(w) rise and a large deviation of the D component at middle latitudes.

Carpenter, D. L.↗

Cislunar Trajectory Design and Maneuver Autonomy for NASA's Moon to Mars Architecture

NASA’s Moon to Mars architecture is an ambitious roadmap of manned cislunar and deep space exploration. The extensive amount of orbital assets required will place a significant burden on ground-based resources, such as communication networks and operations facilities. Spacecraft autonomy is essential for maintaining a vast number of complex missions beyond Earth orbit. To achieve full autonomy, spacecraft must be able to employ methods of robust maneuver design without an explicit dependence on commands sent from the ground. This level of autonomy is needed not only for stationkeeping, but also for outbound transfers. To address the need of spacecraft maneuver design autonomy, this work investigates the use of neural networks (NNs) in a supervised learning environment. A supervised learning approach for NNs allows for a curated training data set, consisting exclusively of perturbations applied to a desired mission concept of operations (ConOps). The proposed approach allows humans on the ground to design a specific mission ConOps before flight, then employ NNs to fly the mission robustly and autonomously. This investigation numerically tests maneuver autonomy in four highly sensitive regions of flight: orbit raising, translunar injection burns, powered lunar flybys, and invariant manifold insertion burns. These straining cases are contextualized by testing them in a demonstration mission, targeting an Earth-Moon L3 orbit. The study first establishes feasibility by automating impulsive burn maneuvers. However, some guidance algorithms will need more intensive commands, such as inertial pointing and angular rates. To validate this method, NN maneuver autonomy is applied to a finite burn model of the demonstration mission. The use of sequential, mission specific maneuvers provide an appropriate testbed to demonstrate the robustness of a NN trained on feasible perturbed states. Moreover, these scenarios provide preliminary proof-of-concept for fully autonomous missions that execute maneuvers without dependence upon explicit command uplinks. As a result, the technological advancement proposed in this work may significantly ease the strain on ground-based mission operations. This would enable complex and autonomous mission execution in cislunar and deep space regimes, filling a technology gap required to support future manned missions.

NASA↗

Study of different approximation in the calculation of g tensors - H2/+/

The theory of the g tensor in one-electron systems is briefly reviewed and calculations are performed in several ways for the hydrogen molecular ion in order to test approximations which must be made for larger systems. Approximate ground state wavefunctions are determined variationally. The first-order wavefunction with respect to the orbit-field perturbation is calculated, and the second-order g tensor is determined. The results of the various approximate calculations are compared and discussed. It is found that the linear combination of atomic orbitals method is rather poor and that two center integrals cannot be neglected.

De Montgolfier, P.↗

Flare loop radiative hydrodynamics. III - Nonlocal radiative transfer effects

The study has three goals. The first is to demonstrate that processes exist whose intrinsic nonlocal nature cannot be represented by local approximations. The second is to elucidate the physical nature and origins of these nonlocal processes. The third is to suggest that the methods and results described here may prove useful in constructing semiempirical models of the chromosphere by means more efficient than trial and error. Matrices are computed that describe the effect of a temperature perturbation at an arbitrary point in the loop on density, hydrogen ionized fraction, total radiative loss rate, and radiative loss rate of selected hydrogen lines and continua at all other points. It is found that the dominant nonlocal radiative transfer effects can be separated into flux divergence coefficient effects and upper level population effects. The former are most important when the perturbation takes place in a region of significant opacity. Upper level population effects arise in both optically thick and thin regions in response to nonlocal density, ionization, and interlocking effects.

Canfield, R. C.↗

Nonlinear systems approach to control system design

Consider some of the control system design methods for plants with nonlinear dynamics. If the nonlinearity is weak relative to the size of the operating region, then the linear methods apply directly. Fixed-gain design may be feasible even for significant nonlinearities. It may be possible to find a single gain which provides adequate control of the linear models at several perturbation points. If the nonlinearity is restricted to a sector, that fact may be used to obtain a fixed-gain controller. Otherwise, a gain may have to be associated with each perturbation point Pi. A gain schedule K(p(v)) is obtained by connecting the perturbation points by a function, say p(v), of the scheduling parameter v (i.e., speed). When the scheduling parameter must be multidimensional, this approach is difficult; the objective is to develop an easier procedure.

Meyer, G.↗

Rotordynamic coefficients for labyrinth seals calculated by means of a finite difference technique

The compressible, turbulent, time dependent and three dimensional flow in a labyrinth seal can be described by the Navier-Stokes equations in conjunction with a turbulence model. Additionally, equations for mass and energy conservation and an equation of state are required. To solve these equations, a perturbation analysis is performed yielding zeroth order equations for centric shaft position and first order equations describing the flow field for small motions around the seal center. For numerical solution a finite difference method is applied to the zeroth and first order equations resulting in leakage and dynamic seal coefficients respectively.

Nordmann, R.↗

Cascade Optimization Strategy for Aircraft and Air-Breathing Propulsion System Concepts

Design optimization for subsonic and supersonic aircraft and for air-breathing propulsion engine concepts has been accomplished by soft-coupling the Flight Optimization System (FLOPS) and the NASA Engine Performance Program analyzer (NEPP), to the NASA Lewis multidisciplinary optimization tool COMETBOARDS. Aircraft and engine design problems, with their associated constraints and design variables, were cast as nonlinear optimization problems with aircraft weight and engine thrust as the respective merit functions. Because of the diversity of constraint types and the overall distortion of the design space, the most reliable single optimization algorithm available in COMETBOARDS could not produce a satisfactory feasible optimum solution. Some of COMETBOARDS' unique features, which include a cascade strategy, variable and constraint formulations, and scaling devised especially for difficult multidisciplinary applications, successfully optimized the performance of both aircraft and engines. The cascade method has two principal steps: In the first, the solution initiates from a user-specified design and optimizer, in the second, the optimum design obtained in the first step with some random perturbation is used to begin the next specified optimizer. The second step is repeated for a specified sequence of optimizers or until a successful solution of the problem is achieved. A successful solution should satisfy the specified convergence criteria and have several active constraints but no violated constraints. The cascade strategy available in the combined COMETBOARDS, FLOPS, and NEPP design tool converges to the same global optimum solution even when it starts from different design points. This reliable and robust design tool eliminates manual intervention in the design of aircraft and of air-breathing propulsion engines where it eases the cycle analysis procedures. The combined code is also much easier to use, which is an added benefit. This paper describes COMETBOARDS and its cascade strategy and illustrates the capability of the combined design tool through the optimization of a subsonic aircraft and a high-bypass-turbofan wave-rotor-topped engine.

Patnaik, Surya N.↗

A Method for Generating Reduced-Order Linear Models of Multidimensional Supersonic Inlets

Simulation of high speed propulsion systems may be divided into two categories, nonlinear and linear. The nonlinear simulations are usually based on multidimensional computational fluid dynamics (CFD) methodologies and tend to provide high resolution results that show the fine detail of the flow. Consequently, these simulations are large, numerically intensive, and run much slower than real-time. ne linear simulations are usually based on large lumping techniques that are linearized about a steady-state operating condition. These simplistic models often run at or near real-time but do not always capture the detailed dynamics of the plant. Under a grant sponsored by the NASA Lewis Research Center, Cleveland, Ohio, a new method has been developed that can be used to generate improved linear models for control design from multidimensional steady-state CFD results. This CFD-based linear modeling technique provides a small perturbation model that can be used for control applications and real-time simulations. It is important to note the utility of the modeling procedure; all that is needed to obtain a linear model of the propulsion system is the geometry and steady-state operating conditions from a multidimensional CFD simulation or experiment. This research represents a beginning step in establishing a bridge between the controls discipline and the CFD discipline so that the control engineer is able to effectively use multidimensional CFD results in control system design and analysis.

Chicatelli, Amy↗

Design and Optimization of Low-thrust Orbit Transfers Using Q-law and Evolutionary Algorithms

Future space missions will depend more on low-thrust propulsion (such as ion engines) thanks to its high specific impulse. Yet, the design of low-thrust trajectories is complex and challenging. Third-body perturbations often dominate the thrust, and a significant change to the orbit requires a long duration of thrust. In order to guide the early design phases, we have developed an efficient and efficacious method to obtain approximate propellant and flight-time requirements (i.e., the Pareto front) for orbit transfers. A search for the Pareto-optimal trajectories is done in two levels: optimal thrust angles and locations are determined by Q-law, while the Q-law is optimized with two evolutionary algorithms: a genetic algorithm and a simulated-annealing-related algorithm. The examples considered are several types of orbit transfers around the Earth and the asteroid Vesta.

low-thrust orbit transfer↗

A (201)Hg+ Comagnetometer for (199)Hg+ Trapped Ion Space Atomic Clocks

A method has been developed for unambiguously measuring the exact magnetic field experienced by trapped mercury ions contained within an atomic clock intended for space applications. In general, atomic clocks are insensitive to external perturbations that would change the frequency at which the clocks operate. On a space platform, these perturbative effects can be much larger than they would be on the ground, especially in dealing with the magnetic field environment. The solution is to use a different isotope of mercury held within the same trap as the clock isotope. The magnetic field can be very accurately measured with a magnetic-field-sensitive atomic transition in the added isotope. Further, this measurement can be made simultaneously with normal clock operation, thereby not degrading clock performance. Instead of using a conventional magnetometer to measure ambient fields, which would necessarily be placed some distance away from the clock atoms, first order field-sensitive atomic transition frequency changes in the atoms themselves determine the variations in the magnetic field. As a result, all ambiguity over the exact field value experienced by the atoms is removed. Atoms used in atomic clocks always have an atomic transition (often referred to as the clock transition) that is sensitive to magnetic fields only in second order, and usually have one or more transitions that are first-order field sensitive. For operating parameters used in the (199)Hg(+) clock, the latter can be five orders of magnitude or more sensitive to field fluctuations than the clock transition, thereby providing an unambiguous probe of the magnetic field strength.

Burt, Eric A.↗

Assessing the Tangent Linear Behaviour of Common Tracer Transport Schemes and Their Use in a Linearised Atmospheric General Circulation Model

The linearity of a selection of common advection schemes is tested and examined with a view to their use in the tangent linear and adjoint versions of an atmospheric general circulation model. The schemes are tested within a simple offline one-dimensional periodic domain as well as using a simplified and complete configuration of the linearised version of NASA's Goddard Earth Observing System version 5 (GEOS-5). All schemes which prevent the development of negative values and preserve the shape of the solution are confirmed to have nonlinear behaviour. The piecewise parabolic method (PPM) with certain flux limiters, including that used by default in GEOS-5, is found to support linear growth near the shocks. This property can cause the rapid development of unrealistically large perturbations within the tangent linear and adjoint models. It is shown that these schemes with flux limiters should not be used within the linearised version of a transport scheme. The results from tests using GEOS-5 show that the current default scheme (a version of PPM) is not suitable for the tangent linear and adjoint model, and that using a linear third-order scheme for the linearised model produces better behaviour. Using the third-order scheme for the linearised model improves the correlations between the linear and non-linear perturbation trajectories for cloud liquid water and cloud liquid ice in GEOS-5.

Tangent Linear Model↗