Search NASA⌕ Search

SEARCH · Search NASA

Results for “Initial conditions problem”

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 307 records · Page 17

The terminal area automated path generation problem

The automated terminal area path generation problem in the advanced Air Traffic Control System (ATC), has been studied. Definitions, input, output and the interrelationships with other ATC functions have been discussed. Alternatives in modeling the problem have been identified. Problem formulations and solution techniques are presented. In particular, the solution of a minimum effort path stretching problem (path generation on a given schedule) has been carried out using the Newton-Raphson trajectory optimization method. Discussions are presented on the effect of different delivery time, aircraft entry position, initial guess on the boundary conditions, etc. Recommendations are made on real-world implementations.

Hsin, C.-C.↗

Senstitivity analysis of horizontal heat and vapor transfer coefficients for a cloud-topped marine boundary layer during cold-air outbreaks

The effects of external parameters on the surface heat and vapor fluxes into the marine atmospheric boundary layer (MABL) during cold-air outbreaks are investigated using the numerical model of Stage and Businger (1981a). These fluxes are nondimensionalized using the horizontal heat (g1) and vapor (g2) transfer coefficient method first suggested by Chou and Atlas (1982) and further formulated by Stage (1983a). In order to simplify the problem, the boundary layer is assumed to be well mixed and horizontally homogeneous, and to have linear shoreline soundings of equivalent potential temperature and mixing ratio. Modifications of initial surface flux estimates, time step limitation, and termination conditions are made to the MABL model to obtain accurate computations. The dependence of g1 and g2 in the cloud topped boundary layer on the external parameters (wind speed, divergence, sea surface temperature, radiative sky temperature, cloud top radiation cooling, and initial shoreline soundings of temperature, and mixing ratio) is studied by a sensitivity analysis, which shows that the uncertainties of horizontal transfer coefficients caused by changes in the parameters are reasonably small.

Chang, Y. V.↗

Goertler vortices in supersonic and hypersonic boundary layers

The problem of Goertler vortices in compressible boundary layers over concave walls is considered. At O(1) wavelengths, the instability is governed by parabolic partial differential equations that are solved numerically to determine the effect of various initial conditions on the development of Goertler vortex instability in compressible boundary layers. The results show that both the velocity and temperature fluctuations may lead to a Goertler vortex. The vortex growth rates determined from the present method are found to differ somewhat from those given by a normal mode solution. At both the supersonic and hypersonic Mach numbers, cooling has a small destabilizing effect. In addition, the most unstable disturbances shift toward lower wavelengths because of thinning of the boundary layer. The results also show that compressibility has a stabilizing effect on the Goertler instability, while the effect of an adverse pressure gradient is found to be destabilizing. The behavior of the Goertler vortex structure with Mach number is also examined. At hypersonic Mach numbers, vortices are located near the edge of the boundary layer for adiabatic wall conditions. However, the entire boundary layer is affected when the wall is cooled.

Spall, R. E.↗

Automated Knowledge Discovery From Simulators

A computational method, SimLearn, has been devised to facilitate efficient knowledge discovery from simulators. Simulators are complex computer programs used in science and engineering to model diverse phenomena such as fluid flow, gravitational interactions, coupled mechanical systems, and nuclear, chemical, and biological processes. SimLearn uses active-learning techniques to efficiently address the "landscape characterization problem." In particular, SimLearn tries to determine which regions in "input space" lead to a given output from the simulator, where "input space" refers to an abstraction of all the variables going into the simulator, e.g., initial conditions, parameters, and interaction equations. Landscape characterization can be viewed as an attempt to invert the forward mapping of the simulator and recover the inputs that produce a particular output. Given that a single simulation run can take days or weeks to complete even on a large computing cluster, SimLearn attempts to reduce costs by reducing the number of simulations needed to effect discoveries. Unlike conventional data-mining methods that are applied to static predefined datasets, SimLearn involves an iterative process in which a most informative dataset is constructed dynamically by using the simulator as an oracle. On each iteration, the algorithm models the knowledge it has gained through previous simulation trials and then chooses which simulation trials to run next. Running these trials through the simulator produces new data in the form of input-output pairs. The overall process is embodied in an algorithm that combines support vector machines (SVMs) with active learning. SVMs use learning from examples (the examples are the input-output pairs generated by running the simulator) and a principle called maximum margin to derive predictors that generalize well to new inputs. In SimLearn, the SVM plays the role of modeling the knowledge that has been gained through previous simulation trials. Active learning is used to determine which new input points would be most informative if their output were known. The selected input points are run through the simulator to generate new information that can be used to refine the SVM. The process is then repeated. SimLearn carefully balances exploration (semi-randomly searching around the input space) versus exploitation (using the current state of knowledge to conduct a tightly focused search). During each iteration, SimLearn uses not one, but an ensemble of SVMs. Each SVM in the ensemble is characterized by different hyper-parameters that control various aspects of the learned predictor - for example, whether the predictor is constrained to be very smooth (nearby points in input space lead to similar output predictions) or whether the predictor is allowed to be "bumpy." The various SVMs will have different preferences about which input points they would like to run through the simulator next. SimLearn includes a formal mechanism for balancing the ensemble SVM preferences so that a single choice can be made for the next set of trials.

Burl, Michael↗

VISUAL CONTROL OF RENDEZVOUS

Space rendezvous involving human occupants will be required in many phases of space flight. In order that rendezvous capability be concurrent with other phases of manned space flight to insure that such tasks as the rotation and addition of crews for manned space station, supply of stations and launching complexes for deep-space probes, and the emergency rescue of personnel may be carried out, it is desirable to determine the ability of a human pilot to control the rendezvous maneuver. The National Aeronautics and Space Administration has conducted a study program to determine vehicle requirements and control techniques for a number of proposed methods for performing piloted space rendezvous. The final step in each phase of the rendezvous study program has been a simulation of the man-machine combination to test its feasibility. Initially, an analytical study was conducted to determine a range of conditions that would exist at the beginning of the terminal phase of rendezvous (ref. 1). Piloted systems were investigated for rendezvous control using these initial conditions. One such study of piloted rendezvous in which all the required control information was displayed to the simulator pilot on instruments is presented in reference 2. While the control of the rendezvous maneuver by use of the instrument display was shown to be feasible, the sufficiently accurate measurement of the necessary information to operate the instrument was realized to be a difficult problem. Consideration was therefore given to the possibility of maneuvering range and range rate required for the control of the rendezvous from visual sightings made during the initial control application. To provide this capability, an analytical study was made. This analytical study, presented in reference 3, derived three techniques for obtaining range and range rate visually and conducted a simplified preliminary situation to test the three analytical techniques. The study reported herein is a continuation of the simulation work begun in reference 3, and is aimed at determining the ability of a pilot to use visual information to reduce or eliminate the need for some or all of the data that would normally be displayed on instruments.

Rendezvous↗

Estimation of distributed parameter systems

The estimation and control of large flexible space structures pose significant new control technology problems. One major problem area is the closed-loop stability of a structure that has been modeled with truncated modal dynamics. Stability problems arise due to the control design process beginning with a deficient model. This work provides the necessary conditions for the optimal estimation of infinite dimensional (partial-differential equation) systems. This approach can be particularly useful for initial control studies and for gaining considerable insight into what the optimal estimator for a truly infinite dimensional system should be. A detailed example of the estimation of the continuous shape of a string in tension is presented.

Schaechter, D. B.↗

CREKID: A computer code for transient, gas-phase combustion of kinetics

A new algorithm was developed for fast, automatic integration of chemical kinetic rate equations describing homogeneous, gas-phase combustion at constant pressure. Particular attention is paid to the distinguishing physical and computational characteristics of the induction, heat-release and equilibration regimes. The two-part predictor-corrector algorithm, based on an exponentially-fitted trapezoidal rule, includes filtering of ill-posed initial conditions, automatic selection of Newton-Jacobi or Newton iteration for convergence to achieve maximum computational efficiency while observing a prescribed error tolerance. The new algorithm was found to compare favorably with LSODE on two representative test problems drawn from combustion kinetics.

Pratt, D. T.↗

Energy and maximum norm estimates for nonlinear conservation laws

We have devised a technique that makes it possible to obtain energy estimates for initial-boundary value problems for nonlinear conservation laws. The two major tools to achieve the energy estimates are a certain splitting of the flux vector derivative f(u)(sub x), and a structural hypothesis, referred to as a cone condition, on the flux vector f(u). These hypotheses are fulfilled for many equations that occur in practice, such as the Euler equations of gas dynamics. It should be noted that the energy estimates are obtained without any assumptions on the gradient of the solution u. The results extend to weak solutions that are obtained as point wise limits of vanishing viscosity solutions. As a byproduct we obtain explicit expressions for the entropy function and the entropy flux of symmetrizable systems of conservation laws. Under certain circumstances the proposed technique can be applied repeatedly so as to yield estimates in the maximum norm.

Olsson, Pelle↗

Cyclogenesis in a saturated environment

The dynamics of baroclinic wave growth in a saturated environment is examined using linear and nonlinear models employing a parameterization of latent heat release that assumes all rising air is saturated, and saturation equivalent potential temperature is conserved on ascent. Piecewise potential vorticity (PV) diagnostics are used to interpret the results. When the stability to vertical displacements in saturated air is allowed to increase with height, as it must in an atmosphere with a constant, positive lapse rate of potential temperature, the growth rates of the most unstable modes of the Eady problem grow only marginally faster than the modes of the dry problem. The vertical variation of moist static stability produces a gradient of moist potential vorticity in the rising air, eliminating the short-wave cutoff present in the dry Eady problem. The destabilization of the short waves is shown to be associated with the interaction between surface potential temperature anomalies and diabatically generated lower-tropospheric potential vorticity anomalies. Nonlinear primitive equation simulations, starting from normal-mode initial conditions, show that while the dry wave grows at nearly the linear growth rate until maximum amplitude is reached, the moist wave grows significantly faster than the linear growth rate at finite amplitude. This enhanced growth is associated with the rapid amplification of a mesoscale PV anomaly generated by latent heat release at the warm front. The rapid amplification of the surface cyclone results from the superposition of the circulation associated with this misoscale PV anomaly upon the circulation associated with the surface and upper boundary potential temperature anomalies. Additional integrations with finite-amplitude initial conditions more typical of atmospheric conditions exhibit similar behavior. It is suggested that many of the rapid cyclogenesis events that occur as upper-troposheric PV anomalies cross the east coasts of continents may arise from the rapid generation of PV anomalies by condensational heating in the moist maritime lower troposphere.

Whitaker, Jeffrey S.↗

Comparison between Computational and Experimental Wall-Shear Stress and Pressure Measurements in a Supersonic Wind Tunnel

We study the performance of two selected RANS models (SST-$k$-$\omega$ and Baldwin-Lomax) to simulate the mean wall-shear and pressure data obtained from a distributed-roughness-related experimental campaign conducted at Ahmic Aerospace, LLC, located in Ohio, US. A structured compressible Navier-Stokes in-house solver, DPLR, is used to conduct the RANS simulations. As the final objective is to perform the studies under high-enthalpy flow conditions, such as inside the AHF facility at NASA-Ames Research Center (ARC), the current investigation focuses on the low-enthalpy cases using Ahmic's high-Reynolds-number supersonic wind tunnel as the initial step to simplify the flow physics. Furthermore, the wind tunnel simulations are conducted with perfectly-smooth walls to further simplify the problem and perform the initial testing of different solver-based RANS models. The Baldwin-Lomax model is observed to be more accurate than the SST model. However, when these RANS models and their roughness augmentation models are utilized to study their accuracy with respect to the experimental datasets, where the roughness effects are activated, both models perform well for a low roughness height case. However, as the roughness height increases, the SST model more accurately predicts the wall data (shear and pressure) than the Baldwin-Lomax model. Therefore, this present work has built a solid foundation for selecting the best RANS model to study the wall data under high-enthalpy flow conditions at the AHF facility soon.

Prakash Shrestha↗

Code for Multiblock CFD and Heat-Transfer Computations

The NASA Glenn Research Center General Multi-Block Navier-Stokes Convective Heat Transfer Code, Glenn-HT, has been used extensively to predict heat transfer and fluid flow for a variety of steady gas turbine engine problems. Recently, the Glenn-HT code has been completely rewritten in Fortran 90/95, a more object-oriented language that allows programmers to create code that is more modular and makes more efficient use of data structures. The new implementation takes full advantage of the capabilities of the Fortran 90/95 programming language. As a result, the Glenn-HT code now provides dynamic memory allocation, modular design, and unsteady flow capability. This allows for the heat-transfer analysis of a full turbine stage. The code has been demonstrated for an unsteady inflow condition, and gridding efforts have been initiated for a full turbine stage unsteady calculation. This analysis will be the first to simultaneously include the effects of rotation, blade interaction, film cooling, and tip clearance with recessed tip on turbine heat transfer and cooling performance. Future plans call for the application of the new Glenn-HT code to a range of gas turbine engine problems of current interest to the heat-transfer community. The new unsteady flow capability will allow researchers to predict the effect of unsteady flow phenomena upon the convective heat transfer of turbine blades and vanes. Work will also continue on the development of conjugate heat-transfer capability in the code, where simultaneous solution of convective and conductive heat-transfer domains is accomplished. Finally, advanced turbulence and fluid flow models and automatic gridding techniques are being developed that will be applied to the Glenn-HT code and solution process.

Fabian, John C.↗

Cosmological Inflation: A Personal Perspective

Approximately twenty five years ago a novel proposal was made to explain two of the outstanding cosmological conundrums, namely those of the Horizon Problem and the Flatness Problem of the Universe. These are the fact that widely separated parts of the sky that have never been in causal contact during the evolution of the Universe have apparently the same CMB temperature and the fact that the mean density of the Universe is very close to the critical one, i.e. very close to the density that separates the closed and open models. These coincidences implied that the corresponding initial condition of the Universe must have been set to exquisite accuracy. This novel proposal posted that at these very early times, the energy density of the Universe was dominated by a fluid which had the equation state attributed to the vacuum (i.e. dominated by tension rather than pressure) and that this led to an exponential expansion of the Universe which was "inflated" by many orders of magnitude of its original size. It was then shown that this "inflation" could provide a resolution of the above outstanding problems. The talk will cover the speaker's personal perspective and contributions to this idea and the subsequent developments over the following 25 years since its inception.

Kazanas, D.↗

Stability and Interaction of Coherent Structure in Supersonic Reactive Wakes

A theoretical formulation and analysis is presented for a study of the stability and interaction of coherent structure in reacting free shear layers. The physical problem under investigation is a premixed hydrogen-oxygen reacting shear layer in the wake of a thin flat plate. The coherent structure is modeled as a periodic disturbance and its stability is determined by the application of linearized hydrodynamic stability theory which results in a generalized eigenvalue problem for reactive flows. Detailed stability analysis of the reactive wake for neutral, symmetrical and antisymmetrical disturbance is presented. Reactive stability criteria is shown to be quite different from classical non-reactive stability. The interaction between the mean flow, coherent structure and fine-scale turbulence is theoretically formulated using the von-Kaman integral technique. Both time-averaging and conditional phase averaging are necessary to separate the three types of motion. The resulting integro-differential equations can then be solved subject to initial conditions with appropriate shape functions. In the laminar flow transition region of interest, the spatial interaction between the mean motion and coherent structure is calculated for both non-reactive and reactive conditions and compared with experimental data wherever available. The fine-scale turbulent motion determined by the application of integral analysis to the fluctuation equations. Since at present this turbulence model is still untested, turbulence is modeled in the interaction problem by a simple algebraic eddy viscosity model. The applicability of the integral turbulence model formulated here is studied parametrically by integrating these equations for the simple case of self-similar mean motion with assumed shape functions. The effect of the motion of the coherent structure is studied and very good agreement is obtained with previous experimental and theoretical works for non-reactive flow. For the reactive case, lack of experimental data made direct comparison difficult. It was determined that the growth rate of the disturbance amplitude is lower for reactive case. The results indicate that the reactive flow stability is in qualitative agreement with experimental observation.

Menon, Suresh↗

Weighted least squares stationary approximations to linear systems.

Investigation of the problem of replacing a certain time-varying linear system by a stationary one. Several quadratic criteria are proposed to aid in determining suitable candidate systems. One criterion for choosing the matrix B (in the stationary system) is initial-condition dependent, and another bounds the 'worst case' homogeneous system performance. Both of these criteria produce weighted least square fits.

Bierman, G. J.↗

Effects of control saturation on the command response of statically unstable aircraft

Hanson and Stengel (1981) have described the effects of saturating controls on the initial-condition response of statically unstable aircraft. In an analysis of the stability augmentation system, consideration was given to maximizing the region of stable response. The present investigation is concerned with an extension of the control saturation analysis to the problem of command response, taking into account the main problem of control design for satisfactory flying qualities (as perceived by the pilot). A model for examining the effects of control saturation on longitudinal stability and command response is developed and applied to an aircraft with a static instability. Three parameters prescribing longitudinal flying qualities are examined with respect to their ability to identify proper flying qualities in the presence of control saturation. The presented analysis provides also satisfactory guidelines for establishing command saturation/stability boundaries.

Hanson, G. D.↗

Eigenmode analysis of unsteady one-dimensional Euler equations

The initial boundary value problem describing the evolution of unsteady linearized perturbations of a steady, uniform subsonic flow is analyzed. The eigenmodes and eigenfrequencies of the system are derived and several examples are presented to illustrate the effect of different boundary conditions on the exponential decay rate of the eigenmodes. The resultant implications for the stability and convergence rates of finite difference computations are discussed.

Giles, M.↗

Topics in spectral methods

After detailing the construction of spectral approximations to time-dependent mixed initial boundary value problems, a study is conducted of differential equations of the form 'partial derivative of u/partial derivative of t = Lu + f', where for each t, u(t) belongs to a Hilbert space such that u satisfies homogeneous boundary conditions. For the sake of simplicity, it is assumed that L is an unbounded, time-independent linear operator. Attention is given to Fourier methods of both Galerkin and pseudospectral method types, the Galerkin method, the pseudospectral Chebyshev and Legendre methods, the error equation, hyperbolic partial differentiation equations, and time discretization and iterative methods.

Gottlieb, D.↗

Laser-heated rotating specimen autoignition test

Specimens of 440 C steel were rotated in a chamber pressurized with oxygen gas and heated with a 5-kW CO2 laser to determine the temperature required for autoignition to occur. Tests included exposures of static and rotating (25,000 rpm) specimens in oxygen pressurized to 5.51 MPa, and with focused laser fluences of more than 3.5 billion W/sq m. Specimen surface temperatures were monitored with a scanning infrared camera. Temperature measurement difficulties were experienced due to a problem with internal reflection inside the test chamber; however, posttest specimen examinations confirmed that surface melt (1371 C) was achieved in several tests. No sustained combustion was initiated in any rotating specimen. One static specimen was ignited. Results indicated that conditions necessary for autoignition of 440 C steel are more dependent on specimen geometry and available heat removal mechanisms. Sustained combustion occurred in the ignited static specimen with an estimated 130 C/sec cooling rate due to conduction. The rotating specimens could not sustain combustion due to a greater conductive/convective cooling rate of about 4000 C/sec and ejection of molten material. These results were applied to the Space Shuttle Main Engine (SSME) oxygen turbopump bearings to conclude that the LOX-cooled 440 C steel bearings cannot sustain combustion initiated by skidding friction.

Au, A. C.↗