Search NASA⌕ Search

SEARCH · Search NASA

Results for “backward”

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 181 records · Page 10

A numerical study of the effects of reverse sweep on a scramjet inlet performance

A comparative numerical study of performance parameters of a similar and an opposite sweep sidewall compression inlet is made. The focus of the study is the investigation of the impact of alternate backward-forward sweep on the compression sidewalls as opposed to back-ward sweep on all the sidewalls. Two equivalent scramjet inlet configurations are designed for this purpose. These inlets have the same wetted areas of compression and expansion and same height and width; but in one inlet all the compression surface are swept back (similar sweep inlet) whereas in the other inlet, alternate surfaces are swept backward and forward (opposite sweep inlet). The cowl closure in both cases begins at the start of the throat region. A three-dimensional Navier-Stokes code is used to calculate the flow through these inlets. Results of these calculations are used to compare the two designs for their performance and flow quality. Effects of boundary- layer ingestion on the performance and overall flow features are also investigated.

Kumar, Ajay↗

Global pressure relaxation for laminar two-dimensional internal flow

This study extends the reduced Navier-Stokes (RNS) global pressure relaxation procedure developed by Rubin and co-workers for external flow to internal flow applications. The streamwise pressure gradient is split into a backward-differenced or initial value component, as in boundary layer marching, and a forward-differenced or boundary value component that represents the elliptic downstream effects. The streamwise convection terms are upwind-differenced and all other streamwise derivatives are backward-differenced. A standard boundary layer marching technique imbedded in a conventional line relaxation technique is obtained. For compressible flow the pressure iteration determines the interior flow interaction as well as the inlet mass flux that is consistent with the outflow pressure boundary condition. Results have been computed for incompressible flow in both rectangular and curved channels, and for subsonic compressible flow in the simulation of an aerofoil in a wind tunnel. Converged solutions were obtained over a range of Reynolds numbers generating small to moderately large separation bubbles.

Rosenbaum, D.↗

Development of a two photon/laser induced fluorescence technique for the detection of atmospheric OH radicals

The development of a new mid-IR laser source was the primary goal. Backward propagating stimulated D2 Raman frequency down conversion of a commercially available 1.06 micron Nd:YAG laser was shown to generate an efficient source of 1.56 micron radiation with near diffraction limited beam quality. The efficient generation of a 2.9 micron laser source was also achieved using backward propagating CH4 Raman frequency down conversion of the 1.56 micron pump. Slightly higher efficiencies were obtained for frequency down conversion of the 1.06 micron Nd:YAG using the H2 Raman shift yielding a near diffraction limited source in the 200 mJ range at 1.9 micron. Similar conversion efficiencies are anticipated as a result of extending the wavelength coverage of recently available Ti:sapphire pulse laser to not only cover the 740 to 860 nm fundamental wavelength range but also the .95 to 1.15 and 1.06 to 1.33 micron range using D2 and H2, respectively. The anticipated sensitivity of a TP-LIF OH sensor using this mid-IR source would give signal limited detection of 1.4 x 10(exp 5) OH/cu cm under boundary layer conditions and 5.5 x 10(exp 4) OH/cu cm under free troposphere sampling conditions for a five minute signal integration period. This level of performance coupled with the techniques non-perturbing nature and freedom from both interferences and background would allow reliable tropospheric OH measurement to be obtained under virtually any ambient condition of current interest, including interstitial and sampling.

Bradshaw, John↗

Instabilities and large structures in reattaching boundary layers

The underlying premise in this work is that large structures in turbulent shear flows are not universal and similar large structures in different flows owe their origins to basically similar instabilities. The turbulent boundary-layer flow over a backward-facing step and also the reattaching transitional boundary layer behind a large spanwise rod have been studied experimentally in parallel primarily to understand the origin of the curious large structures observed in the far downstream region (after 39 step heights) of the former flowfield. Simultaneous flow visualization and hot-wire anemometry and also plain-flow visualization have been carried out. In the backward-facing step, the three-dimensional mixing layer developing after detachment persists even after reattachment. The flow does not recover fully back to a regular turbulent boundary layer even far downstream of reattachment. The turbulence-producing instability of a turbulent boundary layer continues to compete with the wall-bounded mixing layer instability, giving rise to the curious schismatic arrays of large structures.

Bandyopadhyay, Promode R.↗

Experimental analysis of vortex trapping techniques

The use of suction to trap a vortex is investigated experimentally. Vortex trapping using a finite sink and cross-flow obstructions is shown to be viable in the case of low-speed flow over a flat plate airfoil. A dual-fence geometry is found to be more efficient than a backward-facing step due to the greater restriction of vortex movement in the streamwise direction. The backward-facing step geometry provides for an increase in the sink strength required for trapping, while the more efficient dual-fence geometry permits the use of a lower suction force for trapping.

Riddle, Todd W.↗

Planning Assembly Of Large Truss Structures In Outer Space

Report dicusses developmental algorithm used in systematic planning of sequences of operations in which large truss structures assembled in outer space. Assembly sequence represented by directed graph called "assembly graph", in which each arc represents joining of two parts or subassemblies. Algorithm generates assembly graph, working backward from state of complete assembly to initial state, in which all parts disassembled. Working backward more efficient than working forward because it avoids intermediate dead ends.

De Mello, Luiz S. Homem↗

Asymptotic integration algorithms for first-order ODEs with application to viscoplasticity

When constructing an algorithm for the numerical integration of a differential equation, one must first convert the known ordinary differential equation (ODE), which is defined at a point, into an ordinary difference equation (O(delta)E), which is defined over an interval. Asymptotic, generalized, midpoint, and trapezoidal, O(delta)E algorithms are derived for a nonlinear first order ODE written in the form of a linear ODE. The asymptotic forward (typically underdamped) and backward (typically overdamped) integrators bound these midpoint and trapezoidal integrators, which tend to cancel out unwanted numerical damping by averaging, in some sense, the forward and backward integrations. Viscoplasticity presents itself as a system of nonlinear, coupled first-ordered ODE's that are mathematically stiff, and therefore, difficult to numerically integrate. They are an excellent application for the asymptotic integrators. Considering a general viscoplastic structure, it is demonstrated that one can either integrate the viscoplastic stresses or their associated eigenstrains.

Freed, Alan D.↗

A fast, uncoupled, compressible, two-dimensional, unsteady boundary layer algorithm with separation for engine inlets

A finite difference boundary layer algorithm was developed to model viscous effects when an inviscid core flow solution is given. This algorithm solved each boundary layer equation separately, then iterated to find a solution. Solving the boundary layer equations sequentially was 2.4 to 4.0 times faster than solving the boundary layer equations simultaneously. This algorithm used a modified Baldwin-Lomax turbulence model, a weighted average of forward and backward differencing of the pressure gradient, and a backward sweep of the pressure. With these modifications, the boundary layer algorithm was able to model flows with and without separation. The number of grid points used in the boundary layer algorithm affected the stability of the algorithm as well as the accuracy of the predictions of friction coefficients and momentum thicknesses. Results of this boundary layer algorithm compared well with experimental observations of friction coefficients and momentum thicknesses. In addition, when used interactively with an inviscid flow algorithm, this boundary layer algorithm corrected for viscous effects to give a good match with experimental observations for pressures in a supersonic inlet.

Roach, Robert L.↗

Digital control algorithms for microgravity isolation systems

New digital control algorithms were developed to achieve the desired acceleration transmissibility function. The attractive electromagnets have been taken as actuators. The relative displacement and the acceleration of the mass were used as feedback signals. Two approaches were developed to find that controller transfer function in Z-domain, which yields the desired transmissibility at each frequency. In the first approach, the controller transfer function is obtained by assuming that the desired transmissibility is known in Z-domain. Since the desired transmissibility H sub d(S) = 1/(tauS+1)(exp 2) is given in S-domain, the first task is to obtain the desired transmissibility in Z-domain. There are three methods to perform this task: bilinear transformation, and backward and forward rectangular rules. The bilinear transformation and backward rectangular rule lead to improper controller transfer functions, which are physically not realizable. The forward rectangular rule does lead to a physically realizable controller. However, this controller is found to be marginally stable because of a pole at Z=1. In order to eliminate this pole, a hybrid control structure is proposed. Here the control input is composed of two parts: analog and digital. The analog input simply represents the velocity (or the integral of acceleration) feedback; and the digital controller which uses only relative displacement signal, is then obtained to achieve the desired closed-loop transfer function. The stability analysis indicates that the controller transfer function is stable for typical values of sampling period. In the second approach, the aforementioned hybrid control structure is again used. First, an analog controller transfer function corresponding to relative displacement feedback is obtained to achieve the transmissibility as 1/(tauS+1)(exp 2). Then the transfer function for the digital control input is obtained by discretizing this analog controller transfer function via bilinear transformation. The stability of the resulting Z-domain closed loop system is analyzed. Also, the frequency response of the Z-domain closed-loop transfer function is determined to evaluate the performance of the control system.

Sinha, Alok↗

Sufficient conditions for a local minimum of the Bolza problem with multiple terminal point constraints

Sufficient conditions for a weak relative minimum are derived for a form of the Bolza problem of variational calculus. The derivation ties together first-order and second-order conditions in a unified consistent manner, addresses controllability considerations in some detail, and clarifies some small inconsistencies in earlier work. The resulting second-order conditions of optimality involve the integration of fewer backward-sweep matrix elements than the standard conditions in the literature for problems with an unspecified final time. As a result, the backward-sweep matrices have the same general structure and dimensionality whether the final time is specified or free, with the terminal values of two of three sweep matrices being more complicated in the latter case.

Wood, Lincoln J.↗

A fast, uncoupled, compressible, two-dimensional, unsteady boundary layer algorithm with separation for engine inlets

A finite difference boundary layer algorithm was developed to model viscous effects when an inviscid core flow solution is given. This algorithm solved each boundary layer equation separately, then iterated to find a solution. Solving the boundary layer equations sequentially was 2.4 to 4.0 times faster than solving the boundary layer equations simultaneously. This algorithm used a modified Baldwin-Lomax turbulence model, a weighted average of forward and backward differencing of the pressure gradient, and a backward sweep of the pressure. With these modifications, the boundary layer algorithm was able to model flows with and without separation. The number of grid points used in the boundary layer algorithm affected the stability of the algorithm affected the stability of the algorithm as well as the accuracy of the predictions of friction coefficients and momentum thicknesses. Results of this boundary layer algorithm compared well with experimental observations of friction coefficients and momentum thicknesses. In addition, when used interactively with an inviscid flow algorithm, this boundary layer algorithm corrected for viscous effects to give a good match with experimental observations for pressures in a supersonic inlet.

Roach, Robert L.↗

Diabatic dynamic initialization

A generalized dynamical adjustment procedure was applied to a diabatic model to produce balanced initial conditions. Namely, backward adiabatic model integration is followed by forward diabatic model integration, with a high frequency (low pass) filter in the form of the Euler-backward time differencing scheme being applied throughout the whole integration. As a result of the application of such a diabatic dynamic initialization procedure within the Goddard Laboratory for Atmospheres (GLA) 4-D data assimilation system, the following properties of forecasts from initialized fields are achieved right from the beginning of the usual forecast integration: (1) the forecast tendencies (and fields) are free of any noise due to imbalance in initial conditions; (2) the shocks related to an initial imbalance between model physics and dynamics, and especially the substantial initial imbalance of precipitation and evaporation fields, or the initial spin-up effect, are practically removed. Diabatic dynamic initialization has been compared with implicit nonlinear normal mode initialization, and found to be superior in removing the initial spin-up effect, and also in improving the tropical structure. The diabatic dynamic initialization procedure was successfully tested for the GLA system with the use of all conventional data and the GLA satellite data retrievals. It allows a smooth data insertion without any shocks or imbalances, which is highly desirable for efficient functioning of 4-D data assimilation systems. The developed initialization procedure is computationally efficient and in principle easily applicable to different forecast models.

Fox-Rabinovitz, Michael S.↗

Diabatic dynamic initialization

A version of the dynamical adjustment procedure is generalized in order to make the procedure applicable to a diabatic model. In this new initialization procedure, backward adiabatic model integration is followed by forward diabatic model integration, with a high-frequency (low-pass) time filter in the form of the Euler backward time-differencing scheme applied throughout the integration. The new procedure was applied to the Goddard Laboratory for Atmospheres 4D data-assimilation system. It was found that, right from the beginning of the forecast integration, the forecast tendencies (and fields) were free of any noise due to imbalance in initial conditions, and the shocks related to an initial imbalance between model physics and dynamics, or the initial spinup effect, were practically removed.

Fox-Rabinovitz, Michael S.↗

Large eddy simulation of a boundary layer with concave streamwise curvature

One of the most exciting recent developments in the field of large eddy simulation (LES) is the dynamic subgrid-scale model. The dynamic model concept is a general procedure for evaluating model constants by sampling a band of the smallest scales actually resolved in the simulation. To date, the procedure has been used primarily in conjunction with the Smagorinsky model. The dynamic procedure has the advantage that the value of the model constant need not be specified a priori, but rather is calculated as a function of space and time as the simulation progresses. This feature makes the dynamic model especially attractive for flows in complex geometries where it is difficult or impossible to calibrate model constants. The dynamic model was highly successful in benchmark tests involving homogeneous and channel flows. Having demonstrated the potential of the dynamic model in these simple flows, the overall direction of the LES effort at CTR shifted toward an evaluation of the model in more complex situations. The current test cases are basic engineering-type flows for which Reynolds averaged approaches were unable to model the turbulence to within engineering accuracy. Flows currently under investigation include a backward-facing step, wake behind a circular cylinder, airfoil at high angles of attack, separated flow in a diffuser, and boundary layer over a concave surface. Preliminary results from the backward-facing step and cylinder wake simulations are encouraging. Progress on the LES of a boundary layer on a concave surface is discussed. Although the geometry of a concave wall is not very complex, the boundary layer that develops on its surface is difficult to model due to the presence of streamwise Taylor-Gortler vortices. These vortices arise as a result of a centrifugal instability associated with the convex curvature.

Lund, Thomas S.↗

New techniques for reversible compression of multispectral data

While spatial correlations are adequately exploited by standard lossless image compression techniques, little success has been attained in exploiting spectral correlations when dealing with multispectral image data. In this paper, we present some new lossless image compression techniques that capture spectral correlations as well as spatial correlation in a simple and elegant manner. The schemes are based on the notion of a prediction tree, which defines a non-causal prediction model for an image. We present a backward adaptive technique and a forward adaptive technique. We then give a computationally efficient way of approximating the backward adaptive technique. The approximation gives good results and is extremely easy to compute. Simulation results show that for high spectral resolution images, significant savings can be made by using spectral correlations in addition to spatial correlations. Furthermore, the increase in complexity incurred in order to make these gains is minimal.

Memon, N. D.↗

Simulation study of the role of ion kinetics in low-frequency wave train evolution

The evolution of uniform, parallel propagating, low-frequency (less than or approx. ion cyclotron) wave trains is followed with a one-dimensional hybrid numerical code with fluid electrons and particle ions. We show that moderate amplitude (delta B/B less than 1/2) wave trains give instabilities and saturated states which differ completely from pure fluid evolution. This is most clearly seen when beta greater than 1 and instability exists for wavenumbers both below and above the wavenumber of an inital, left-handed wave train or pump wave. For corresponding parameters a fluid theory gives only a narrow range of instability above the pump wavenumber where decay and beat instabilities can occur. in simulations wave energy inverse cascades to smaller wavenumbers and into a greater number of forward than backward going waves. In fluids energy by decay goes mostly to backward ones of smaller wavenumber, and energy by beat goes mostly to forward ones of larger wavenumber. Neither fluid instability explains simulation results. The instability is saturated by thermalizing ions and sometimes exciting small wavenumber electrostatic or acoustic modes. In contrast, saturation in fluids first occurs by generating the harmonics of the growing linear modes. Harmonic generation is mostly absent in simulations. Simulations are carried out to long times and mostly reach a limit beyond which no further significant evolution can occur. Application to Alfvenic fluctuations in the solar wind is discussed.

Vasquez, Bernard J.↗

Beat, modulational, and decay instabilities of a circularly polarized Alfven wave

A circulary polarized low-frequency electomagnetic pump wave propoagating along an ambient magnetic field is known to be unstable to the growth of several parallel-propagating parametric instabilties. If ion-cyclotron effects are retained in a two-fluid description, the dispersion relation is a sixth-order polynomial. We present a series of new analytical approximations to this dispersion relation. We emphasize new results for the beat instability that occurs as an interaction of the forward prpagating upper sideband with the backward propagating lower sideband. The nature of the beat instabitlity depends on beta = (v(sub sound)/v(sub A)(exp 2) and on the sense of polarization of the pump wave. The beat and decay instabilities can occur together if the pump is left-handed (i.e., ion resonant) and if beta is less than or approximately 1, but they cannot occur together if the pump is right-handed. For a left-handed pump the beat mode is the only instability if beta is greater than or approximately 1. If the pump is right-handed and beta is greater than or approximately 1, then the beat instability exists only when the pump amplitude exceeds a threshold value, and the beat will be the only instability if the pump amplitude is large enough to stabilize the modulational instability. If the pump is left-handed and beta is less than or approximately 1, then the beat mode is stabilized when the pump amplitude becomes sufficiently large. The beat instability primarily produces a forward propagating transverse wave in the upper sideband. Thus if beta is greater than or approximately 1, the instabilities considered here do not produce the backward propagating waves which are thought to affect turbulence and the evolution of cross helicity in the solar wind. New analytical results are presented also for the decay and modulational instabilites when beta is approximately equal to 1.

Hollweg, Joseph V.↗

Parametric instabilities of parallel-propagating Alfven waves: Some analytical results

We consider the stability of a circularly polarized Alfven wave (the pump wave) which propagates parallel to the ambient magnetic field. Only parallel-propagating perturbations are considered, and we ignore dispersive effects due to the ion cyclotron frequency. The dissipationless MHD equations are used throughout; thus possibibly important effects arising from Landau and transit time damping are omitted. We derive a series of analytical approximations to the dispersion relation using A = (Delta B/B(sub O))(exp 2) as a small expansion parameter; Delta B is the pump amplitude, and B(sub O) is the ambient magnetic field strength. We find that the plasma beta (the square of the ratio of the sound speed to the Alfven speed) plays a crucial role in determining the behavior of the parametric instabilities of the pump. If 0 less than beta less than 1 we find the familiar result that the pump decays into a forward propagating sound wave and a backward propagating Alfven wave with maximum growth rate gamma(sub max) varies A(sup 1/2), but beta cannot be too close to 0 or to 1. If beta approx. 1, we find gamma(sub max) varies A(sup 3/4), if beta greater than 1, we find gamma(sub max) varies A(sup 3/2), while if beta approx. 0, we obtain gamma(sub max) varies A(sup 1/3); moreover, if beta approx. 0 there is a nearly purely growing instability. In constrast to the familiar decay instability, for which the backward propagating Alfven wave has lower frequency and wavenumber than the pump, we find that if beta greater than or approx. equal to 1 the instability is really a beat instability which is dominated by a transverse wave which is forward propagating and has frequency and wavenumber which are nearly twice the pump values. Only the decay instability for 0 less than beta less than 1 can be regarded as producing two recognizable normal modes, namely, a sound wave and an Alfven wave. We discuss how the different characteristics of the instabilities may affect the evolution of Alfven waves in the solar wind. However, for a solar wind in which beta approx. 1 the growth times of the instabilities are probably too long for these instabilities to have an appreciable effect inside 1 AU.

Jayanti, V.↗