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,369 records · Page 76

Turbulent MHD transport coefficients - An attempt at self-consistency

In this paper, some multiple scale perturbation calculations of turbulent MHD transport coefficients begun in earlier papers are first completed. These generalize 'alpha effect' calculations by treating the velocity field and magnetic field on the same footing. Then the problem of rendering such calculations self-consistent is addressed, generalizing an eddy-viscosity hypothesis similar to that of Heisenberg for the Navier-Stokes case. The method also borrows from Kraichnan's direct interaction approximation. The output is a set of integral equations relating the spectra and the turbulent transport coefficients. Previous 'alpha effect' and 'beta effect' coefficients emerge as limiting cases. A treatment of the inertial range can also be given, consistent with a -5/3 energy spectrum power law. In the Navier-Stokes limit, a value of 1.72 is extracted for the Kolmogorov constant. Further applications to MHD are possible.

Chen, H.↗

Enhanced viscous flow drag reduction using acoustic excitation

Large eddy break up devices (LEBUs) constitute a promising method of obtaining drag reduction in a turbulent boundary layer. Enhancement of the LEBU effectiveness by exciting its trailing edge with acoustic waves phase locked to the large scale structure influencing the momentum transfer to the wall is sought. An initial estimate of the required sound pressure level for an effective pulse was obtained by considering the magnitude of the pressure perturbations at the near wake of a thin plate in inviscid flow. Detailed skin friction measurments were obtained in the flow region downstream of a LEBU excited with acoustic waves. The data are compared with skin friction measurements of a simply manipulated flow, without acoustic excitation and with a plain flow configuration. The properties and the scales of motion in the flow regime downstream of the acoustically excited LEBU are studied. A parametric study based upon the characteristics of the acoustic input was pursued in addition to the careful mapping of the drag reduction phenomenon within the acoustically manipulated boundary layer. This study of boundary layer manipulation has lead to improved skin friction drag reduction and further understanding of the turbulent boundary layer.

Nagel, R. T.↗

Incorporation of a priori gravity field information in satellite orbit determination using bin parameters

A method to determine satellite orbits using tracking data and a priori gravitational field is described. The a priori constraint on the orbit dynamics is determined by the covariance matrix of the spherical harmonic coefficients for the gravity model, so that the optimal combination of the measurements and gravitational field is achieved. A set of bin parameters is introduced to represent the perturbation of the gravitational field on the position of the satellite orbit. The covariance matrix of a conventional gravity model is transformed into that for the bin parameters by the variational partial derivatives. The covariance matrices of the bin parameters and the epoch state are combined to form the covariance matrix of the satellite positions at the measurement times. The combined matrix is used as the a priori information to estimate the satellite positions with measurements.

Wu, Jiun-Tsong↗

Accurate ab Initio Quartic Force Fields, Vibrational Frequencies, and Heats of Formation for FCN, FNC, ClCN, and ClNC

The XCN and XNC (X = F, Cl) isomers have been investigated using the CCSD(T) method in conjunction with correlation consistent basis sets. Equilibrium geometries, harmonic frequencies, anharmonic constants, fundamental frequencies, and heats of formation have been evaluated. Agreement with experiment for the fundamental frequencies is very good, even for nu(sub 2), for CICN, which is subject to a strong Fermi resonance with 2nu(sub 3). It is also shown that a second-order perturbation theory approach to solving the nuclear Schroedinger equation gives results in excellent agreement with essentially exact variational calculations. This is true even for nu(sub 2) of ClCN, provided that near-singular terms are eliminated from the perturbation theory formulas and the appropriate Fermi interaction energy matrix is then diagonalized. A band at 615/cm, tentatively assigned as the Cl-N stretch in ClNC in matrix isolation experiments, is shown not to be due to ClNC. Accurate atomization energies are determined and are used to evaluate accurate heats of formation (3.1 +/- 1.5, 33.2 +/- 1.5, 72.6 +/- 1.5, and 75.9 +/- 1.5 kcal/mol for FCN, ClCN, FNC, and ClNC, respectively). It is expected that the theoretical heats of formation for FCN, FNC, and ClNC are the most accurate available.

Lee, Timothy J.↗

Fractal mechanisms in the electrophysiology of the heart

The mathematical concept of fractals provides insights into complex anatomic branching structures that lack a characteristic (single) length scale, and certain complex physiologic processes, such as heart rate regulation, that lack a single time scale. Heart rate control is perturbed by alterations in neuro-autonomic function in a number of important clinical syndromes, including sudden cardiac death, congestive failure, cocaine intoxication, fetal distress, space sickness and physiologic aging. These conditions are associated with a loss of the normal fractal complexity of interbeat interval dynamics. Such changes, which may not be detectable using conventional statistics, can be quantified using new methods derived from "chaos theory.".

NASA Discipline Cardiopulmonary↗

Mean Element Propagations Using Numerical Averaging

The long-term evolution characteristics (and stability) of an orbit are best characterized using a mean element propagation of the perturbed two body variational equations of motion. The averaging process eliminates short period terms leaving only secular and long period effects. In this study, a non-traditional approach is taken that averages the variational equations using adaptive numerical techniques and then numerically integrating the resulting EOMs. Doing this avoids the Fourier series expansions and truncations required by the traditional analytic methods. The resultant numerical techniques can be easily adapted to propagations at most solar system bodies.

orbit propagations↗

Approximate Green's function methods for HZE transport in multilayered materials

A nonperturbative analytic solution of the high charge and energy (HZE) Green's function is used to implement a computer code for laboratory ion beam transport in multilayered materials. The code is established to operate on the Langley nuclear fragmentation model used in engineering applications. Computational procedures are established to generate linear energy transfer (LET) distributions for a specified ion beam and target for comparison with experimental measurements. The code was found to be highly efficient and compared well with the perturbation approximation.

Wilson, John W.↗

A synoptic analysis of the first AVE-SESAME '79 period

Key features of a severe convection observed during April 10-11, 1979 as part of the Atmospheric Variability Experiment-Severe Environmental Storms and Mesoscale Experiment (AVE-SESAME) are examined. Three-hourly rawinsonde readings from 23 stations were taken, and vertical motion and divergence parameters are considered. The data were converted into a 127 km grid at the surface, and at 50 mb intervals from 900 mb to 100 mb by an objective analysis scheme, while a kinematic method was used to compute vertical motion. A weak upper tropospheric short wave trough was found to propagate from New Mexico into the Texas panhandle, while a jet maximum propagated eastward. The development of a strong wind maximum over Oklahoma and Kansas was associated with a rapid increase in upper-level divergence and the development of a small-scale pressure perturbation in the Texas panhandle, as well as a low-level jet and convergence, which led to rapid changes over the Red River Valley, where stability was decreased.

Moore, J. T.↗

Electronic structures of HgTe and CdTe surfaces and HgTe/CdTe interfaces

A Green's-function method has been used to study the surface and interface electronic structures of the II-VI compounds HgTe and CdTe. Localized surface and resonance states near the cation-terminated (100) surface of CdTe and the anion-terminated surface of HgTe have been found for the ideal surfaces. The energies and strengths of these surface states are altered by surface perturbations. The bulk states near the surface are drastically modified by the creation of the surface, but the band gaps remain unchanged. Numerical evaluation of the local densities of states at the Gamma and J points shows that, at the (100) interface of HgTe/CdTe, the previously observed surface states are no longer present. However, in the interface region, bulk states of one material penetrate some distance into the other material.

Schick, J. T.↗

Suppressing Transients In Digital Phase-Locked Loops

Loop of arbitrary order starts in steady-state lock. Method for initializing variables of digital phase-locked loop reduces or eliminates transients in phase and frequency typically occurring during acquisition of lock on signal or when changes made in values of loop-filter parameters called "loop constants". Enables direct acquisition by third-order loop without prior acquisition by second-order loop of greater bandwidth, and eliminates those perturbations in phase and frequency lock occurring when loop constants changed by arbitrarily large amounts.

Thomas, J. B.↗

Dynamic Stability and Gravitational Balancing of Multiple Extended Bodies

Feasibility of a non-invasive compensation scheme was analyzed for precise positioning of a massive extended body in free fall using gravitational forces influenced by surrounding source masses in close proximity. The N-body problem of classical mechanics is a paradigm used to gain insight into the physics of the equivalent N-body problem subject to control forces. The analysis addressed how a number of control masses move around the proof mass so that the proof mass position can be accurately and remotely compensated when exogenous disturbances are acting on it, while its sensitivity to gravitational waves remains unaffected. Past methods to correct the dynamics of the proof mass have considered active electrostatic or capacitive methods, but the possibility of stray capacitances on the surfaces of the proof mass have prompted the investigation of other alternatives, such as the method presented in this paper. While more rigorous analyses of the problem should be carried out, the data show that, by means of a combined feedback and feed-forward control approach, the control masses succeeded in driving the proof mass along the specified trajectory, which implies that the proof mass can, in principle, be balanced via gravitational forces only while external perturbations are acting on it. This concept involves the dynamic stability of a group of massive objects interacting gravitationally under active control, and can apply to drag-free control of spacecraft during missions, to successor gravitational wave space borne sensors, or to any application requiring flying objects to be precisely controlled in position and attitude relative to another body via gravitational interactions only.

Quadrelli, Marco↗

Verifying Implementation of the Dryden Turbulence Model and MIL-F-8785 Gust Gradient

Turbulence modeling in human-in-the-loop simulation is important to assessing aircraft handling qualities and pilot performance and to provide additional realism for pilot training. In the simulation community, the Dryden turbulence spectra is a popular choice for modeling the linear turbulent gusts because its rational form is efficiently reproduced by passing white noise through linear filters. The MIL-F-8785 gust gradients similarly use additional linear filters to model the gradient of the turbulent gust over the wing, and it represents the gust gradients as perturbations to the air-relative rotational rates. The Cockpit Motion Facility at NASA Langley Research Center (LaRC) models continuous random turbulence using the Dryden one-dimensional spectra and MIL-F-8785 gust gradient. The facility recently reviewed and updated its verification of these models as part of an initiative to improve motion cueing under turbulence. This exercise introduced improved methods for verifying the turbulence models and led to rediscovery of model assumptions that informed improvements to implementation.

Madden, Michael M.↗

Early Information Parameter-Set Analysis for Satellite Close Approaches using Machine Learning

Understanding orbital mechanics is essential in space flight and navigation applications, and leveraging modern force models for flight path projection remains an important aspect in space mission design and operation. However, force models do not capture all the dynamics or perturbations in the space environment and thus are subject to errors in predicting the state vectors. The further out the predicted miss distance between spacecraft is from the time of closest approach (TCA), the larger the propagated errors in the predicted miss distance at TCA is. The dependency on these force models for spacecraft flight state prediction calls for a more reliable method that can quantify, or even reduce, these propagated errors. With recent advances in the field artificial intelligence, specifically in machine and deep learning algorithms, a model that implements these approaches can improve on the modern force model approach. The goal for this work is to provide an early-information decision-making threshold, in order to prioritize risk assessment implementation, given the ongoing increase of space objects. In analyzing the relationship of several parameters from conjunction data messages(CDMs) and solar information, early information becomes viable in miss distance prediction with unsupervised learning techniques, which learn the parameters that are linked together with miss distance and probability of collision (Pc) variables. Another approach implemented for identifying relationships within CDMs is supervised learning, in which a shallow neural network binary classifier learns to distinguish events with Pc values¡108. These parameters detected in the unsupervised process are then applied to a regression neural network, which predicts the miss distance at TCA for a specific event within a given uncertainty bound. For the regression neural network, a Long Short Term Memory (LSTM) neural network is implemented, which yields memory about each time step in an event. Using an LSTM network, the model learns to predict miss distance within 0.2km of the value measured at TCA. Although there is a limited amount of "close miss" data to train a network, the network learns to associate parameters, like large energy dissipation rates with the secondary object, with an elevated Pc

Brianna I. Robertson↗

Singular perturbation equations for flexible satellites

The dynamic model of a flexible satellite with n generalized structural position coordinates requires the solution of a set of N coupled nonlinear ordinary or partial differential equations. For single composite body satellites, N is equal to (n + 3). For dual-spin systems, N is equal to (n + 9). These equations usually involve time derivatives up to the second order. For large values of n, linearization of the system has so far been the only practicable way of solution. The present study shows a method of avoiding this linearization by reducing the N equations to three or nine nonlinear, coupled, first-order ordinary differential equations involving only the angular velocities of the composite bodies. The solutions for these angular velocities lead to linear equations in the n generalized structural position coordinates, which can then be solved by known methods.

Huang, T. C.↗

Imp 6 measurements in the distant polar cusp during substorms

High time-resolution data from the magnetic field, plasma, energetic particle, and VLF wave experiments performed aboard Imp 6 in a study of the distant dayside cusp during substorms are described. The cusp was studied when its location was slightly equatorward of its normal location and the geomagnetic dipole was tilted in the appropriate direction. The data support both reconnection and diffusion as methods of particle entry to the magnetosphere. The evidence (1) indicates an acceleration process to explain enhancements of 400 to 600-km/sec protons above their magnetosheath intensities, and (2) suggests convection of field lines over the polar cap as a means of explaining the lack of low-energy protons near the low-latitude boundary of the cusp. Magnetic field fluctuations, a perturbation vector, ion cyclotron waves, and an abrupt change in the intensity of both whistler waves and electrostatic waves are characterized.

Fairfield, D. H.↗

Computation of asymmetric supersonic flows around cones at large incidence

The Schiff-Steger parabolized Navier-Stokes (PNS) code has been modified to allow computation of conical flowfields around cones at high incidence. The improved algorithm of Degani and Schiff has been incorporated with the PNS code. This algorithm adds the cross derivative and circumferential viscous terms to the original PNS code and modifies the algebraic eddy viscosity turbulence model to take into account regions of so called cross-flow separation. Assuming the flowfield is conical (but not necessarily symmetric) a marching stepback procedure is used: the solution is marched one step downstream using improved PNS code and the flow variables are then scaled to place the solution back to the original station. The process is repeated until no change in the flow variables is observed with further marching. The flow variables are then constant along rays of the flowfield. The experiments obtained by Bannik and Nebbeling were chosen as a test case. In these experiments a cone of 7.5 deg. half angle at Mach number 2.94 and Reynolds number 1.372 x 10(7) was tested up 34 deg. angle of attack. At high angle of attack nonconical asymmetric leeward side vortex patterns were observed. In the first set of computations, using an earlier obtained solution of the above cone for angle of attack of 22.6 deg. and at station x=0.5 as a starting solution, the angle of attack was gradually increased up to 34 deg. During this procedure the grid was carfully adjusted to capture the bow shock. A stable, converged symmetric solution was obtained. Since the numerical code converged to a symmetric solution which is not the physical one, the stability was tested by a random perturbation at each point. The possible effect of surface roughness or non perfect body shape was also investigated. It was concluded that although the assumption of conical viscous flows can be very useful for certain cases, it can not be used for the present case. Thus the second part of the investigation attempted to obtain a marching (in space) solution with the PNS method using the conical solution as initial data. Finally, the solution of the full Navier-Stokes equations was carried out.

Degani, David↗

An investigation into critical aspects of a new form of low energy lunar transfer, the Belbruno-Miller trajectories

This paper documents an integrated Belbruno-Miller (B-M) trajectory and its corresponding injection period. The B-M trajectories use Weak Stability Boundaries (WSB) resulting from four-body perturbative dynamics between the earth, moon, sun, and spacecraft to significantly reduce maneuver requirements for lunar transfer. It is determined that the presented nominal trajectory has a viable 4-day injection period considering maneuver constraints. The nominal B-M trajectory produces saving varying from 150 m/s to 222 m/s over traditional Hohmann means thus pointing to the practical utility of the WSB method. Energy savings are extracted at the price of extended flight time of approximately 6-months. Usefulness of this procedure was recently dramatized by the Japanese spacecraft Hiten when it arrived a the moon on October 2, 1991 following its entry into B-M trajectory on April 25, 1991.

Krish, V.↗

Theoretical study of the transonic lift of a double-wedge profile with detached bow wave

A theoretical study is described of the aerodynamic characteristics at small angle of attack of a thin, double-wedge profile in the range of supersonic flight speed in which the bow wave is detached. The analysis is carried out within the framework of the transonic (nonlinear) small-disturbance theory, and the effects of angle of attack are regarded as a small perturbation on the flow previously calculated at zero angle. The mixed flow about the front half of the profile is calculated by relaxation solution of a suitably defined boundary-value problem for transonic small-disturbance equation in the hodograph plane (i.e., the Tricomi equation). The purely supersonic flow about the rear half is found by an extension of the usual numerical method of characteristics. Analytical results are also obtained, within the framework of the same theory, for the range of speed in which the bow wave is attached and the flow is completely supersonic.

Vincenti, Walter G↗