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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Natural convection in steady solidification - Finite element analysis of a two-phase Rayleigh-Benard problem

Galerkin finite-element approximations and Newton's method for solving free boundary problems are combined with computer-implemented techniques from nonlinear perturbation analysis to study solidification problems with natural convection in the melt. The Newton method gives rapid convergence to steady state velocity, temperature and pressure fields and melt-solid interface shapes, and forms the basis for algebraic methods for detecting multiple steady flows and assessing their stability. The power of this combination is demonstrated for a two-phase Rayleigh-Benard problem composed of melt and solid in a veritical cylinder with the thermal boundary conditions arranged so that a static melt with a flat melt-solid interface is always a solution. Multiple cellular flows bifurcating from the static state are detected and followed as Rayleigh number is varied. Changing the boundary conditions to approach those appropriate for the vertical Bridgman solidification system causes imperfections that eliminate the static state. The flow structure in the Bridgman system is related to those for the Rayleigh-Benard system by a continuous evolution of the boundary conditions.

Chang, C. J.↗
NASA NTRS1982

Combined four-wall interference assessment in two-dimensional airfoil tests

Two different procedures are examined for combining the correction method developed by Barnwell and Sewall for the effects of the sidewall boundary layer in two-dimensional wind-tunnel tests with the assessment and correction method due to Kemp for the effects of the upper and lower tunnel walls. One procedure utilizes the similarity transformation defined by Sewall to eliminate velocity perturbations induced by the sidewall boundary layer from consideration in the assessment and correction method as sources of error. The other procedure combines perturbations from all four walls before assessment and correction. In general, the two procedures yield different corrections to Mach number and angle of attack, either of which can be considered valid. Pressure distributions corrected for higher order interference effects are provided by either procedure.

Kemp, W. B., Jr.↗
NASA NTRS2013

Monte Carlo Bayesian Inference on a Statistical Model of Sub-gridcolumn Moisture Variability Using High-resolution Cloud Observations : Sensitivity Tests and Results - Part II

Part I presented a Monte Carlo Bayesian method for constraining a complex statistical model of GCM sub-gridcolumn moisture variability using high-resolution MODIS cloud data, thereby permitting large-scale model parameter estimation and cloud data assimilation. This part performs some basic testing of this new approach, verifying that it does indeed significantly reduce mean and standard deviation biases with respect to the assimilated MODIS cloud optical depth, brightness temperature and cloud top pressure, and that it also improves the simulated rotational-Ramman scattering cloud optical centroid pressure (OCP) against independent (non-assimilated) retrievals from the OMI instrument. Of particular interest, the Monte Carlo method does show skill in the especially difficult case where the background state is clear but cloudy observations exist. In traditional linearized data assimilation methods, a subsaturated background cannot produce clouds via any infinitesimal equilibrium perturbation, but the Monte Carlo approach allows finite jumps into regions of non-zero cloud probability. In the example provided, the method is able to restore marine stratocumulus near the Californian coast where the background state has a clear swath. This paper also examines a number of algorithmic and physical sensitivities of the new method and provides guidance for its cost-effective implementation. One obvious difficulty for the method, and other cloud data assimilation methods as well, is the lack of information content in the cloud observables on cloud vertical structure, beyond cloud top pressure and optical thickness, thus necessitating strong dependence on the background vertical moisture structure. It is found that a simple flow-dependent correlation modification due to Riishojgaard (1998) provides some help in this respect, by better honoring inversion structures in the background state.

Monte Carlo Bayesian inference↗
NASA NTRS1992

The mass of (1) Ceres from perturbations on (348) May

The most promising ground-based technique for determining the mass of a minor planet is the observation of the perturbations it induces in the motion of another minor planet. This method requires careful observation of both minor planets over extended periods of time. The mass of (1) Ceres has been determined from the perturbations on (348) May, which made three close approaches to Ceres at intervals of 46 years between 1891 and 1984. The motion of May is clearly influenced by Ceres, and by using different test masses for Ceres, a search was made to determine the mass of Ceres that minimizes the residuals in the observations of May.

Williams, Gareth V.↗
NASA NTRS1986

Dynamics of a discrete auroral arc

Porcupine Flight 4 data were used to determine the field-aligned currents associated with a southward moving discrete auroral arc in the postmidnight sector. Three different methods were used for determining the field-aligned current which should give identical results if the arcs are quasi-stationary and no parallel electric field exists between the payload and the dynamo region of the ionosphere. As long as the rocket is above the arc, the three methods agree. The integral of precipitating electron flux, the local magnetic field perturbations, and the divergence of the horizontal Pedersen current all indicate an upward current of 5 + or - 3 microamperes/sq m. Immediately north of the arc a strong downward current of about 10-20 microamperes/sq m is detected. The magnitude, however, is not well known because the rocket's velocity relative to the arc cannot be clearly established. Further north of the southward moving arc, the two methods that can be applied (magnetic field perturbations and divergence of the horizontal Pedersen current) yield contradictory results not only about the magnitude of the current but also about the direction of the current. It is suggested that this discrepancy is due to time-dependent electric field.

Bruening, K.↗
NASA NTRS1987

The use of atmospheric measurements to constrain model predictions of ozone change from chlorine perturbations

Atmospheric photochemistry models have been used to predict the sensitivity of the ozone layer to various perturbations. These same models also predict concentrations of chemical species in the present day atmosphere which can be compared to observations. Model results for both present day values and sensitivity to perturbation depend upon input data for reaction rates, photodissociation rates, and boundary conditions. A method of combining the results of a Monte Carlo uncertainty analysis with the existing set of present atmospheric species measurements is developed. The method is used to examine the range of values for the sensitivity of ozone to chlorine perturbations that is possible within the currently accepted ranges for input data. It is found that model runs which predict ozone column losses much greater than 10 percent as a result of present fluorocarbon fluxes produce concentrations and column amounts in the present atmosphere which are inconsistent with the measurements for ClO, HCl, NO, NO2, and HNO3.

Douglass, Anne R.↗
NASA NTRS2007

Scattering of Low Energy Electrons and Positrons from Hydrogenic Systems and Applications

While the electron scatters from the target, the target no longer stays in its original form. One of the first methods to take into account the distortion of the target at low incident energies is the method of polarized orbitals. In this method the wave function for the process is written using the first-order perturbation theory and the equation for the scattering function is derived from the Schradinger equation. This method has been very successful in calculating the phase shifts and therefore the cross sections at various energies. The total wave function can be used to calculate photoionization cross sections. The disadvantage of this approach is that the method is not variational and therefore does not provide bounds on the phase shifts. These difficulties can be overcome by using the Feshbach projection operator formalism. This approach has been employed for the scattering of electrons and positrons from targets. Results of various calculations will be discussed.

Bhatia, Anand K.↗
NASA NTRS2003

The Updated Umkehr Ozone Retrieval Algorithm and its Validation against Satellite Data

Improvements to the Umkehr ozone profile retrieval algorithm have been developed and are now being evaluated. The updated algorithm is able to simulate observations more accurately and provides data output that is easier to analyze. Among the new diagnostic capabilities that the updated algorithm provides is the averaging kernel (AK) method. The AK approach allows studying how the algorithm responds when a small perturbation is made in a particular layer of the atmosphere [Rodgers 1976, 1990]. We will use the AK method to define precisely what Umkehr should measure given a set of profiles measured by other platforms. This method allows us to compare trends and offsets in data more accurately than it has been done in the past. The updated Umkehr retrievals will be validated against SAGE II ozone profiles as well as SSBUV ozone profile data. We will discuss possible reasons for offset between data and differences in derived ozone profile trends. Considerable variability of the ozone profile within the 10-degree latitude envelope creates noise in the SAGE matching dataset and makes comparisons difficult. To eliminate this problem, the SAGE and Umkehr data had been previously de-seasonalized by subtracting the latitude/season dependent ozone climatology. However, the remaining noise in the ozone residuals was still considerably high for trend analysis and was attributed to longitude variability of SAGE sampling. The new ozone climatology (Labow, NASA) that has longitude dependent ozone variability will be used to minimize contribution of sampling noise in comparisons of satellite and ground station. The comparison of zenith-sky radiances (Umkehr N-value measurements) synthesized for a given set of SAGE profiles will be used to determine whether SAGE-derived N-values agree with the Umkehr-measured N-values. The instrumental effects will be discussed. Both the Umkehr data and SAGE II measurements will be analyzed for their information about ozone variability and loss and recovery rates at the mid- and upper (40 km) levels. The updated long-term Umkehr dataset can be used to provide high quality information for identifying signs of ozone recovery. The long Umkehr historical record can provide additional information for separating the dynamic and chemical mechanisms of depletion, and can help the community better understand climate change effects.

Petropavlovskikh, I.↗
NASA NTRS1981

Development of an efficient procedure for calculating the aerodynamic effects of planform variation

Numerical procedures to compute gradients in aerodynamic loading due to planform shape changes using panel method codes were studied. Two procedures were investigated: one computed the aerodynamic perturbation directly; the other computed the aerodynamic loading on the perturbed planform and on the base planform and then differenced these values to obtain the perturbation in loading. It is indicated that computing the perturbed values directly can not be done satisfactorily without proper aerodynamic representation of the pressure singularity at the leading edge of a thin wing. For the alternative procedure, a technique was developed which saves most of the time-consuming computations from a panel method calculation for the base planform. Using this procedure the perturbed loading can be calculated in about one-tenth the time of that for the base solution.

Mercer, J. E.↗
NASA NTRSDate not supplied

A Numerical Method for Computing the State Transition Matrix Using Poincare Integral Invariants

The Poincare integral invariants describe the volumes of sets in Hamiltonian phase space. We use these invariants to derive a new numerical procedure for obtaining the state transition matrix (STM), which can be applied to both conservative and nonconservative systems. The method is analogous to a finite difference approximation of the STM, where perturbed states are numerically propagated along with the reference trajectory. We discuss the mathematical similarities between this new STM and existing methods, show numerical results for orbital motion and uncertainty propagation, and discuss new insights afforded by the Hamiltonian properties of phase flow.

state transition matrix↗
NASA NTRS1974

A unified approach for the application of general perturbation theories to the artificial satellite problem

An Encke-type method is developed as well as a variation of parameters method, both of which use a non-Keplerian reference orbit. The regularized time is used in the numerical integration and the optimum value of n is found for each type of orbit investigated. Accuracy and computation time comparisons are made with a classical Cowell method. It should be noted that the element formulations developed and tested were found to be exceptionally numerically stable in the sense that it was possible to achieve numerical consistency order of very long integration periods - a property not available with the Cowell formulation - and therefore these formulations may be helpful for high precision calculations.

Alfriend, K. T.↗
NASA NTRS1989

Recovering the full velocity and density fields from large-scale redshift-distance samples

A new method for extracting the large-scale three-dimensional velocity and mass density fields from measurements of the radial peculiar velocities is presented. Galaxies are assumed to trace the velocity field rather than the mass. The key assumption made is that the Lagrangian velocity field has negligible vorticity, as might be expected from perturbations that grew by gravitational instability. By applying the method to cosmological N-body simulations, it is demonstrated that it accurately reconstructs the velocity field. This technique promises a direct determination of the mass density field and the initial conditions for the formation of large-scale structure from galaxy peculiar velocity surveys.

Bertschinger, Edmund↗
NASA NTRS1973

Solution of the stochastic control problem in unbounded domains.

Bellman's dynamic programming equation for the optimal index and control law for stochastic control problems is a parabolic or elliptic partial differential equation frequently defined in an unbounded domain. Existing methods of solution require bounded domain approximations, the application of singular perturbation techniques or Monte Carlo simulation procedures. In this paper, using the fact that Poisson impulse noise tends to a Gaussian process under certain limiting conditions, a method which achieves an arbitrarily good approximate solution to the stochastic control problem is given. The method uses the two iterative techniques of successive approximation and quasi-linearization and is inherently more efficient than existing methods of solution.

Robinson, P.↗
NASA NTRS1976

A lifting surface theory for the analysis of nonplanar lifting systems

A new nonlinear, nonplanar lifting surface theory is presented. The method is regarded as a lifting surface theory in that the effects of wing thickness are neglected, but none of the usual small perturbation assumptions inherent in most other lifting surface theories are made. The method represents nonplanar lifting systems by distributed vorticity, including the leading edge singular behavior characteristic of thin wings. The method is well suited to the computation of induced drag of nonplanar systems because leading edge suction is calculated from the leading edge singularity. The method has been used to compute the induced drag benefit of winglets (vortex diffusers), and the agreement with NASA experimental data is excellent.

Goldhammer, M. I.↗
NASA NTRS1975

Steady and Oscillatory, Subsonic and Supersonic, Aerodynamic Pressure and Generalized Forces for Complex Aircraft Configurations and Applications to Flutter

A general method for analyzing aerodynamic flows around complex configurations is presented. By applying the Green function method, a linear integral equation relating the unknown, small perturbation potential on the surface of the body, to the known downwash is obtained. The surfaces of the aircraft, wake and diaphragm (if necessary) are divided into small quadrilateral elements which are approximated with hyperboloidal surfaces. The potential and its normal derivative are assumed to be constant within each element. This yields a set of linear algebraic equations and the coefficients are evaluated analytically. By using Gaussian elimination method, equations are solved for the potentials at the centroids of elements. The pressure coefficient is evaluated by the finite different method; the lift and moment coefficients are evaluated by numerical integration. Numerical results are presented, and applications to flutter are also included.

Chen, L. T.↗