Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solvers”

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,117 records · Page 62

Sonic Boom Ground Noise Minimization Via the Adjoint Method

A multidisciplinary design optimization methodology to directly minimize the ground-level noise generated by the sonic boom of a high altitude supersonic body is presented. A Cartesian Euler flow solver is coupled with an atmospheric propagation tool to create a ground-noise analysis capability for supersonic bodies. Adjoint formulations for both the flow solver and propagation tool are also coupled to compute sensitivities of shape variations in the body in a highly efficient manner. A gradient-based optimizer is then introduced to forge a valuable design capability. Output-based mesh adaptation that is driven directly by ground-level noise is employed to increase accuracy and provide error estimation. The design method is demonstrated first on a simple axisymmetric body with few design variables to evaluate the efficacy of the optimization scheme. Guided by the results of this initial case, the problem is then repeated with somewhat different design variables to further demonstrate the capabilities of the design method. Finally, the method is applied to a real-world problem by optimizing control surface deflection settings of a low-boom aircraft to minimize ground noise while maintaining trimmed, level flight.

ARMD↗

Content and Representation of Information Needed to Support Time-Constrained Problem Solving

NASA’s current mission-operations paradigm originated with Project Mercury and endured with minimum evolution through the Apollo Program, Space Shuttle Program, and ISS missions. At its foundation is a near-complete real-time dependence on a ground team to manage the combined state of the mission, vehicle, and crew. Utilizing many engineers and operators with broad and deep expertise; large, distributed datasets including extensive telemetry; and expansive analytical and computing power, this ground team has served as the safety net for crewed spaceflight missions over the past 60 years. This approach must change to address challenges associated with missions beyond low Earth orbit (BLEO), including infrequent resupply, reduced ability to evacuate, and delayed communications that prohibit real-time operational support. We anticipate that a necessary part of this change will be increased independence for the crew, as roles and responsibilities traditionally performed by ground teams move on board the vehicle. While many risks are associated with Earth-independent operations, one particular concern is ensuring that the crew will have adequate onboard support to perform urgent problem solving when communication with the ground is delayed or intermittent. A key resource that enables the ground team to respond to anomalies quickly and effectively is the extraordinary expertise and experience it possesses. It is comprised of 80+ experts on at any given time, with a combined 600+ years of system-specific experience across 22 unique console disciplines. A small crew will face the unprecedented challenge of independently responding to anomalies that have historically been handled by a team 20 times their size. Another important resource upon which the ground heavily relies to support procedure execution and anomaly response is data. The amount of telemetry data that each flight controller monitors is extensive. In addition, as the ground team works to further assess impacts, trouble shoot, identify workarounds, and oversee procedure execution, it accesses and synthesizes engineering and procedure information, as well as system build, test, and configuration documentation. It is not feasible nor useful to put all these data onboard as crews become more Earth independent. Each member of a small Mars mission small crew will have multiple roles beyond monitoring telemetry and data gathering, and multiple roles within anomaly resolution processes, thereby limiting their capacity for copious amounts of information. Moreover, while access is necessary, it alone is insufficient. Information will need to be compiled, refined, and represented appropriately to support the crew’s reduced attention and expertise. This work seeks to understand the content and representation of information needed to support time-constrained problem solving and decision making by the crew without real-time ground support. To build this understanding, we first surveyed the literature, focusing on how expert problem solvers construct and manipulate their mental models. Next, we interviewed expert problem solvers in spaceflight and analogous domains and surveyed industry solutions for data presentation. Finally, we analyzed current spaceflight operations by investigating flight controller anomaly resolution processes during ISS training simulations and real operational events. These methods led to creating a problem-solving framework that details common themes and features of attending to, assessing, analyzing, and acting on problems in complex, time-constrained domains. Using this framework and the results of our analysis, we identified conceptual data representations needed for crew-led problem-solving. Preliminary onboard user interface concepts to meet identified needs will be presented.

anomaly response↗

Application of a patched-grid algorithm to the F/A-18 forebody-leading-edge extension configuration

A patched-grid algorithm for the analysis of complex configurations with an implicit, upwind-biased Navier-Stokes solver is presented. Results from both a spatial-flux and a time-flux conservation approach to patching across zonal boundaries are presented. A generalized coordinate transformation with a biquadratic geometric element is used at the zonal interface in order to treat highly stretched viscous grids and arbitrarily shaped zonal boundaries. Applications are made to the F/A-18 forebody-Ieading-edge extension configuration at subsonic, high-alpha conditions. Computed surface flow patterns compare well with ground-based and flight-test results; the large effect of Reynolds number on the forebody flowfield is shown.

Leading edge extension↗

Rotational relaxation of molecular hydrogen at moderate temperatures

Using a coupled rotation-vibration-dissociation model the rotational relaxation times for molecular hydrogen as a function of final temperature (500-5000 K), in a hypothetical scenario of sudden compression, are computed. The theoretical model is based on a master equation solver. The bound-bound and bound-free transition rates have been computed using a quasiclassical trajectory method. A review of the available experimental data on the rotational relaxation of hydrogen is presented, with a critical overview of the method of measurements and data reduction, including the sources of errors. These experimental data are then compared with the computed results.

COUPLED ROTATION VIBRATION DIS↗

A Patched-Grid Algorithm for Complex Configurations Directed Towards the F/A-18 Aircraft

A patched-grid algorithm for the analysis of complex configurations with an implicit, upwind-biased Navier-Stokes solver is presented. Results from both a spatial-flux and a time-flux conservation approach to patching across zonal boundaries are presented. A generalized coordinate transformation with a biquadratic geometric element is used at the zonal interface in order to treat highly stretched viscous grids and arbitrarily-shaped zonal boundaries. Applications are made to the F-18 forebody-strake configuration at subsonic, high-alpha conditions. Computed surface flow patterns compare well with ground-based and flight-test results; the large effect of Reynolds number on the forebody flow-field is shown.

High alpha research vehicle↗

Large Scale Finite Element Modeling Using Scalable Parallel Processing

An iterative solver for use with finite element codes was developed for the Cray T3D massively parallel processor at the Jet Propulsion Laboratory. Finite element modeling is useful for simulating scattered or radiated electromagnetic fields from complex three-dimensional objects with geometry variations smaller than an electrical wavelength.

finite element modeling parallel processing iterat↗

Massively Parallel Solution of Poisson Equation on Coarse Grain MIMD Architectures

In this paper a new algorithm, designated as Fast Invariant Imbedding algorithm, for solution of Poisson equation on vector and massively parallel MIMD architectures is presented. This algorithm achieves the same optimal computational efficiency as other Fast Poisson solvers while offering a much better structure for vector and parallel implementation. Our implementation on the Intel Delta and Paragon shows that a speedup of over two orders of magnitude can be achieved even for moderate size problems.

parallel algorithms MIMD parallel architectures↗

Robust Path Planning and Feedback Design Under Stochastic Uncertainty

Autonomous vehicles require optimal path planning algorithms to achieve mission goals while avoiding obstacles and being robust to uncertainties. The uncertainties arise from exogenous disturbances, modeling errors, and sensor noise, which can be characterized via stochastic models. Previous work defined a notion of robustness in a stochastic setting by using the concept of chance constraints. This requires that mission constraint violation can occur with a probability less than a prescribed value.In this paper we describe a novel method for optimal chance constrained path planning with feedback design. The approach optimizes both the reference trajectory to be followed and the feedback controller used to reject uncertainty. Our method extends recent results in constrained control synthesis based on convex optimization to solve control problems with nonconvex constraints. This extension is essential for path planning problems, which inherently have nonconvex obstacle avoidance constraints. Unlike previous approaches to chance constrained path planning, the new approach optimizes the feedback gain as wellas the reference trajectory.The key idea is to couple a fast, nonconvex solver that does not take into account uncertainty, with existing robust approaches that apply only to convex feasible regions. By alternating between robust and nonrobust solutions, the new algorithm guarantees convergence to a global optimum. We apply the new method to an unmanned aircraft and show simulation results that demonstrate the efficacy of the approach.

autonomuys vehicles↗

Viscel: A general purpose computer program for analysis of linear viscoelastic structures, volume 2

The VISCEL program is a general purpose computer program developed for equilibrium analysis of linear viscoelastic structures. The program is written in FORTRAN 5 language to operate on the Univac 1108 computer under the EXEC 8 operating system. The program, an extension of the linear equilibrium problem solver ELAS, is an updated and extended version of its earlier form written for the IBM 7094 computer. Finite element matrix displacement approach coupled with the synchronized material property concept, utilizing incremental time steps, was adopted for the solution presented. The step-by-step procedure involves solution of recursive equations in the time domain, which takes into account the memory of material properties. Incremental and accumulative displacements and stresses are obtained at the end of each time step. In order to minimize the extent of computations resulting from accumulative effects of material memory, the program provides an option which enables the employment of constant time steps in the logarithmic scale. Program documentation is presented.

Gupta, K. K.↗

VISCEL: A general-purpose computer program for analysis of linear viscoelastic structures (user's manual), volume 1

This program, an extension of the linear equilibrium problem solver ELAS, is an updated and extended version of its earlier form (written in FORTRAN 2 for the IBM 7094 computer). A synchronized material property concept utilizing incremental time steps and the finite element matrix displacement approach has been adopted for the current analysis. A special option enables employment of constant time steps in the logarithmic scale, thereby reducing computational efforts resulting from accumulative material memory effects. A wide variety of structures with elastic or viscoelastic material properties can be analyzed by VISCEL. The program is written in FORTRAN 5 language for the Univac 1108 computer operating under the EXEC 8 system. Dynamic storage allocation is automatically effected by the program, and the user may request up to 195K core memory in a 260K Univac 1108/EXEC 8 machine. The physical program VISCEL, consisting of about 7200 instructions, has four distinct links (segments), and the compiled program occupies a maximum of about 11700 words decimal of core storage.

Gupta, K. K.↗

QA4, a language for artificial intelligence.

Introduction of a language for problem solving and specifically robot planning, program verification, and synthesis and theorem proving. This language, called question-answerer 4 (QA4), embodies many features that have been found useful for constructing problem solvers but have to be programmed explicitly by the user of a conventional language. The most important features of QA4 are described, and examples are provided for most of the material introduced. Language features include backtracking, parallel processing, pattern matching, set manipulation, and pattern-triggered function activation. The language is most convenient for use in an interactive way and has extensive trace and edit facilities.

Derksen, J. A. C.↗

Rapid finite-difference computation of subsonic and transonic aerodynamic flows

Rapid iterative (or semidirect) computation methods are developed for the finite-difference solution of the nonlinear equations of subsonic and transonic aerodynamics. At each iteration, a fast, direct elliptic algorithm solves the entire computation field. In an application to subsonic flow over a lifting airfoil, the full nonlinear stream-function equation is solved. Finally, a direct Cauchy-Riemann solver is used for the nonlinear transonic small-disturbance equations for a biconvex airfoil. At M = 0.7, t/c = 0.1 (subcritical), three iterations on a 39 x 32 mesh (totaling 2.45 sec on an IBM 360/67 computer) obtain convergence within 0.1%. A slightly supercritical case requires seven iterations (6.75 sec) for convergence within 1%.

Martin, E. D.↗

A generalized-capacity-matrix technique for computing aerodynamic flows

A numerical generalized-capacity-matrix technique is developed for application to aerodynamic flow computations. This technique allows the very fast direct (noniterative) numerical elliptic solvers to be used in problems with arbitrary internal boundaries and with a wide class of boundary conditions, including numerical application of the Kutta condition on an airfoil without iteration. Accuracy, speed, and usefulness of the technique are demonstrated with linear problems for potential flows over airfoil shapes. The method's main advantages, however, can be exploited within iterative procedures for a variety of complex flow problems governed by systems of equations not necessarily elliptic or linear.

Martin, E. D.↗

A simplified package for calculating with splines

This package is designed to solve some of the elementary problems of spline interpretation and least squares fit. The subroutines fall into three basic categories. The first category involves computations with a given spline and/or knot sequence; the second category involves routines which calculate the coefficients of splines which perform a certain task (such as interpolation or least squares fit); and the last category is a banded equation solver for specific linear equations.

Smith, P. W.↗

SSME structural computer program development: BOPACE programmer manual

The programming manual for the two-dimensional Boeing Plastic Analysis Capability for Engines (BOPACE) program written in FORTRAN IV is presented. The BOPACE flow logic is summarized along with subroutines and input/output files. Emphasis is placed on linear equation solver routines written as an independent package.

Source record↗

A fast semidirect method for computing transonic aerodynamic flows

A fast, semidirect, iterative computational method, previously introduced for finite-difference solution of subsonic and slightly supercritical flow over airfoils, is extended both to apply to strongly supercritical conditions and to include full second-order accuracy in computing inviscid flows over airfoils. The nonlinear small-disturbance equations are solved iteratively by a direct, linear, elliptic solver. General, fully conservative, type-dependent difference equations are formulated, including parabolic- and shock-point transition operators that provide consistency with the integral conservation laws. These equations specialize to either first-order or to fully second-order-accurate equations. Various free parameters are evaluated for rapid convergence of the first-order scheme. Resulting pressure distributions and computing times are compared with the improved Murman-Cole line-relaxation method.

Martin, E. D.↗

Accelerated iterative calculation of transonic nacelle flowfields

A method is presented for the calculation of inviscid, supercritical flowfields about axisymmetric inlet cowls. A finite-difference calculation is performed in a simple, rectangular domain obtained from the nacelle geometry by a nearly-conformal mapping procedure. Type-dependent finite-differences are constructed using a coordinate-independent, 'rotated' differencing scheme. Methods of accelerating convergence of the iterative solution are demonstrated including a hybrid fast-Poisson-solver/relaxation scheme and an extrapolated relaxation procedure. Calculated pressure distributions are compared with experimental data for a variety of Mach numbers and mass-flow ratios, and show generally good agreement.

Caughey, D. A.↗

Advances in the application of fast semidirect computational methods in transonic flow

The paper uses finite-difference algorithms called 'fast direct elliptic solvers' within an iteration scheme for the rapid solution of the equations of inviscid transonic aerodynamics. The methods are called 'direct' (or semidirect) because the entire computational field is solved at once rather than in successive traverses over the field. These semidirect iterative methods have been limited here to the investigation of two-dimensional steady inviscid flow over airfoils in subsonic free stream.

Martin, E. D.↗