SEARCH · Engineering Papers
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.
This paper describes dynamic trade-off evalation (DTE), a new technique that has been developed to improve the performance of real-time problem solving systems. The DTE technique is most suitable for environments in which the requirement for meeting time constraints is of equal or greater importance to that of providing optimally intelligent solutions.
Helicopter flowfields are highly unsteady, nonlinear and three-dimensional. In forward flight and in hover, the rotor blades interact with the tip vortex and wake sheet developed by either itself or the other blades. This interaction, known as blade-vortex interactions (BVI), results in unsteady loading of the blades and can cause a distinctive acoustic signature. Accurate and cost-effective computational fluid dynamic solutions that capture blade-vortex interactions can help rotor designers and engineers to predict rotor performance and to develop designs for low acoustic signature. Such a predictive method must preserve a blade's shed vortex for several blade revolutions before being dissipated. A number of researchers have explored the requirements for this task. This paper will outline some new capabilities that have been added to the NASA Ames' OVERFLOW code to improve its overall accuracy for both vortex capturing and unsteady flows. To highlight these improvements, a number of case studies will be presented. These case studies consist of free convection of a 2-dimensional vortex, dynamically pitching 2-D airfoil including light-stall, and a full 3-D unsteady viscous solution of a helicopter rotor in forward flight In this study both central and upwind difference schemes are modified to be more accurate. Central difference scheme is chosen for this simulation because the flowfield is not dominated by strong shocks. The feature of shock-vortex interaction in such a flow is less important than the dominant blade-vortex interaction. The scheme is second-order accurate in time and solves the thin-layer Navier-Stokes equations in fully-implicit manner at each time-step. The spatial accuracy is either second and fourth-order central difference or third-order upwind difference using Roe-flux and MUSCLE scheme. This paper will highlight and demonstrate the methods for several sample cases and for a helicopter rotor. Preliminary computations on a rotor were performed by using this method and are in the process of documentation.
Modeling and simulation are important to understand laser operation and to optimize and design device functions. Numerical simulation of VCSEL (Vertical Cavity Surface Emitting Lasers) has been largely based on solving time-independent Helmholtz equation or time dependent coupled mode equations. There are various advantages for choosing these approaches. However, the disadvantages are also apparent. The former cannot handle dynamical mode competition seen in VCSELs, while the latter assumes a given type and number of modes a priori. Furthermore, the microscopic physics of heterstructures and electron-hole plasma is very often represented by a few parameters such as linear gain coefficients and the linewidth enhancement factor. These are over simplification of space and frequency (wavelength) dependent gain and refractive index functions. When the space-time dynamical operation of VCSELs becomes important, these simple approximations become questionable. In this paper, we apply a recently developed model for edge-emitting lasers to a gain guided VCSEL for space-time domain simulation. This model takes into account the actual nonlinear dependence of gain and refractive index on frequency and carrier density within the frame work of the effective Bloch equations. The corresponding partial differential equations are solved directly by finite difference methods. Laser behavior with increasing pumping current is investigated in detail. Special attention is paid to the dynamical competition of the transverse modes.
A novel conversion algorithm is presented that combines the fidelity of indirect optimization methods with the generality of direct methods to more easily solve time-optimal, finite-burn pseudo-rendezvous problems. An algorithm is described that converts a set of multiple-impulses, representing the entirety or a portion of a high- or low-thrust maneuver, to an exact time optimal finite-burn trajectory for a thrust limited, constant exhaust velocity spacecraft. A pseudo-rendezvous problem is one that yields a solution whose final time, position and velocity state is equal to that of the original post-impulsive trajectory. An iterative adjoint-control transformation is used to initialize the optimal control two-point boundary value problem. Examples are shown for both high and low-thrust non-coplanar Earth orbit transfers, as well as a low-thrust Hohmann-type Earth-Mars transfer.
Two computational methods to determine an open loop time optimal control sequence for a simple single spool turbojet engine are described by a set of nonlinear differential equations. Both methods are modifications of widely accepted algorithms which can solve fixed time unconstrained optimal control problems with a free right end. Constrained problems to be considered have fixed right ends and free time. Dynamic programming is defined on a standard problem and it yields a successive approximation solution to the time optimal problem of interest. A feedback control law is obtained and it is then used to determine the corresponding open loop control sequence. The Fletcher-Reeves conjugate gradient method has been selected for adaptation to solve a nonlinear optimal control problem with state variable and control constraints.
This paper describes an intelligent user interface that is currently under development. The interface serves as a front end for real-time spacecraft monitoring software. The software operates under circumstances in which neither an intelligent human without automated assistance nor an automated system without intelligence are sufficiently effective. The user interface is supported by dynamic trade-off evaluation (DTE), a new technique that has been developed to automate general types of performance trade-offs in real-time problem solving systems. In this application, DTE is used to perform dynamic input data management for the purpose of determining which input data should be monitored in time constrained situations and how to present the monitoring information to a human analyst who has the ultimate responsibility for the spacecraft. This application has demonstrated that DTE can be used to dynamically vary the data that is monitored, making it possible to detect and correctly analyze all anomalous data by examining only a subset of the total input data. In carefully structured experimental evaluations that use real spacecraft data and real decision making, DTE provides the ability to handle a three-fold increase in input data (in real-time) without loss of performance and to intelligently present the information to a mission analyst.
A general method for computing laminar and turbulent boundary layers for two-dimensional time-dependent flows is presented. The method uses an eddy-viscosity formulation to model the Reynolds shear-stress term and a very efficient numerical method to solve the governing equations. The model was applied to steady two-dimensional and three-dimensional flows and was shown to give good results. A discussion of the numerical method and the results obtained by the present method for both laminar and turbulent flows are discussed. Based on these results, the method is efficient and suitable for solving time-dependent laminar and turbulent boundary layers.
A theoretical overview of the response coefficient method for heat conduction transients with time-dependent input forcing functions is presented with a number of illustrative applications. The method may be the most convenient and economical if the same problem is to be solved many times with different input-time histories or if the solution time is relatively long. The method is applicable to a wide variety of problems, including irregular geometries, position-dependent boundary conditions, position-dependent physical properties, and nonperiodic irregular input histories. Nonuniform internal energy generation rates within the structure can also be handled by the method. The area of interest is long-time solutions, in which initial condition is unimportant, and not the early transient period. The method can be applied to one dimensional problems in cartesian, cylindrical, and spherical coordinates as well as to two dimensional problems in cartesian and cylindrical coordinates.
A code has been developed for simulating high Reynolds number transonic flow fields of arbitrary configuration. An explicit finite-difference method with time splitting is used to solve the time-dependent equations for compressible turbulent flow. A nonorthogonal computational mesh of arbitrary configuration facilitates the description of the flow field. The code is applied to simulate the flow over an 18 percent thick circular-arc biconvex airfoil at zero angle of attack and free-stream Mach number of 0.775. A simple mixing-length model is used to describe the turbulence and chord Reynolds numbers of 1, 2, 4, and 10 million are considered. The solution describes in sufficient detail both the shock-induced and trailing-edge separation regions, and provides the profile and friction drag.
A code has been developed for simulating high Reynolds number transonic flow fields of arbitrary configuration. An explicit finite-difference method with time splitting is used to solve the time-dependent equations for compressible turbulent flow. A nonorthogonal computational mesh of arbitrary configuration facilitates the description of the flow field. The code is applied to simulate the flow over a two-dimensional 18 percent thick circular-arc biconvex airfoil at zero angle of attack for several different Reynolds numbers and a free-stream Mach number of 0.775.
A laboratory design facility for digital microprocessor implementation of linear-quadratic-Gaussian feedback compensators is described. Outputs from user interactive programs for solving infinite time horizon LQ regulator and Kalman filter problems were conditioned for implementation on the laboratory microcomputer system. The software consisted of two parts: an offline high-level program for solving the LQ Ricatti equations and generating associated feedback and filter gains and a cross compiler/macro assembler which generates object code for the target microprocessor system. A PDP 11/70 with a UNIX operating system was used for all high level program and data management, and the target microprocessor system is an Intel MDS (8080-based processor). Application to the control of a two dimensional inverted pendulum is presented and issues in expanding the design/prototyping system to other target machine architectures are discussed.
The paper describes a laboratory design facility for digital microprocessor implementation of Linear-Quadratic-Gaussian feedback compensators. Outputs from user interactive programs for solving infinite time horizon LQ regulator and Kalman filter problems are conditioned for implementation on a laboratory microcomputer system. The software consists of two parts: (1) an off-line high-level program for solving the LQ Ricatti equations and generating associated feedback and filter gains, and (2) a cross compiler/macro assembler which generates object code for the target microprocessor system. Application to the control of a two dimensional inverted pendulum and expanding the design/prototyping system to other target machine architectures are discussed.
The fuel optimal spin axis reorientation problem for large and small maneuver times is solved using a new numerical technique. The approach, which employs the adaptive grid bisection search, is shown to be useful in determining the fuel optimal impulsive control strategy, particularly when rapid maneuvers are desired. For large enough maneuver times, the method produces the two-impulse reorientation maneuver. For even larger maneuver times, the existence of multiple optimal pulsing schemes is demonstrated.
This paper presents a linear Integer Programming model for managing air traffic flow in the United States. The decision variables in the model are departure delays and predeparture reroutes of aircraft whose trajectories are predicted to cross weather-impacted regions of the National Airspace System. The model assigns delays to a set of flights while ensuring their trajectories are free of any conflicts with weather. In a deterministic setting, there is no airborne holding due to unexpected weather incursion in a flight s path. The model is applied to solve a large-scale traffic flow management problem with realistic weather data and flight schedules. Experimental results indicate that allowing rerouting can reduce departure delays by nearly 57%, but it is associated with an increase in total airborne time due to longer routes flown by aircraft. The computation times to solve this problem were significantly lower than those reported in the earlier studies.
The gust response of a coupled hingeless rotor-fuselage system is studied in both hover and forward flight. Each rotor blade undergoes flap bending, lag bending, and torsional deflections. The blades are discretized into beam elements, each with fifteen nodal degrees of freedom. The fuselage is permitted three translational (vertical, longitudinal, and lateral) and two rotational (pitch and roll) degrees of freedom. The formulation considers a three-dimensional gust field wherein each gust velocity component can have an arbitrary variation in space and time. Aerodynamic loads are obtained using quasisteady strip theory. Wake-induced effects are introduced through dynamic inflow modeling. Dynamic stall and reverse flow effects are also included. Equations governing the rotor-fuselage gust response are linearized about the vehicle propulsive trim state and the blade steady-state defected positions, and solved by time integration. The effects of several parameters on the helicopter gust response are illustrated, including dynamic inflow, lag stiffness, forward speed, gust profile, gust penetration rate, and gust velocity direction.
Approximation method for solving singe-degree-of- freedom differential roll equation with time dependent coefficient - ballistic trajectory anallysis
We are interested in illumination conditions and the temperature distribution within the upper two meters of regolith near the lunar poles. Here, areas exist receiving almost constant illumination near areas in permanent shadow, which were identified as potential exploration sites for future missions. For our study a numerical simulation of the illumination and thermal environment for lunar near-polar regions is needed. Our study is based on high-resolution, twenty meters per pixel and 400 x 400 km large polar Digital Terrain Models (DTMs), which were derived from Lunar Orbiter Laser Altimeter (LOLA) data. Illumination conditions were simulated by synthetically illuminating the LOLA DTMs using the horizon method considering the Sun as an extended source. We model polar illumination for the central 50 x 50 km subset and use it as an input at each time-step (2 h) to evaluate the heating of the lunar surface and subsequent conduction in the sub-surface. At surface level we balance the incoming insolation with the subsurface conduction and radiation into space, whereas in the sub-surface we consider conduction with an additional constant radiogenic heat source at the bottom of our two-meter layer. Density is modeled as depth-dependent, the specific heat parameter as temperature-dependent and the thermal conductivity as depth- and temperature-dependent. We implemented a fully implicit finite-volume method in space and backward Euler scheme in time to solve the one-dimensional heat equation at each pixel in our 50 x 50 km DTM. Due to the non-linear dependencies of the parameters mentioned above, Newton's method is employed as the non-linear solver together with the Gauss-Seidel method as the iterative linear solver in each Newton iteration. The software is written in OpenCL and runs in parallel on the GPU cores, which allows for fast computation of large areas and long time scales.
Multidisciplinary design optimization (MDO) for large-scale engineering problems poses many challenges (e.g., the design of an efficient concurrent paradigm for global optimization based on disciplinary analyses, expensive computations over vast data sets, etc.) This work focuses on the application of distributed schemes for massively parallel architectures to MDO problems, as a tool for reducing computation time and solving larger problems. The specific problem considered here is configuration optimization of a high speed civil transport (HSCT), and the efficient parallelization of the embedded paradigm for reasonable design space identification. Two distributed dynamic load balancing techniques (random polling and global round robin with message combining) and two necessary termination detection schemes (global task count and token passing) were implemented and evaluated in terms of effectiveness and scalability to large problem sizes and a thousand processors. The effect of certain parameters on execution time was also inspected. Empirical results demonstrated stable performance and effectiveness for all schemes, and the parametric study showed that the selected algorithmic parameters have a negligible effect on performance.