Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear equations”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

A generalized procedure for constructing an upwind-based TVD scheme

A generalized formulation for constructing second- and higher-order accurate TVD (total variation diminishing) schemes is presented. A given scheme is made TVD by limiting antidiffusive flux differences with some nonlinear functions, so-called limiters. The general idea of the formulation and its mathematical proof of Harten's TVD conditions is shown by applying the Lax-Wendroff method to a scalar nonlinear equation and constant-coefficient system of conservation laws. For the system of equations, several definitions are derived for the argument used in the limiter function and present their performance to numerical experiments. Then the formulation is formally extended to the nonlinear system of equations. It is demonstrated that use of the present procedure allows easy conversion of existing central or upwind, and second- or higher-order differencing schemes so as to preserve monotonicity and to yield physically admissible solutions. The formulation is simple mathematically as well as numerically; neither matrix-vector multiplication nor Riemann solver is required. Roughly twice as much computational effort is needed as compared to conventional scheme. Although the notion of TVD is based on the initial value problem, application to the steady Euler equations of the formulation is also made. Numerical examples including various ranges of problems show both time- and spatial-accuracy in comparison with exact solutions.

Liou, Meng-Sing↗
NASA NTRS1991

Maneuver simulations of flexible spacecraft by solving TPBVP

The optimal control of large angle rapid maneuvers and vibrations of a Shuttle mast reflector system is considered. The nonlinear equations of motion are formulated by using Lagrange's formula, with the mast modeled as a continuous beam. The nonlinear terms in the equations come from the coupling between the angular velocities, the modal coordinates, and the modal rates. Pontryagin's Maximum Principle is applied to the slewing problem, to derive the necessary conditions for the optimal controls, which are bounded by given saturation levels. The resulting two point boundary value problem (TPBVP) is then solved by using the quasilinearization algorithm and the method of particular solutions. In the numerical simulations, the structural parameters and the control limits from the Spacecraft Control Lab Experiment (SCOLE) are used. In the 2-D case, only the motion in the plane of an Earth orbit or the single axis slewing motion is discussed. In the 3-D slewing, the mast is modeled as a continuous beam subjected to 3-D deformations. The numerical results for both the linearized system and the nonlinear system are presented to compare the differences in their time response.

Bainum, Peter M.↗
NASA NTRS1976

Acoustic radiation from sources convected at transonic speeds

Each supersonic jet contains a long segment of transonic flow, and it is known through acoustic experiments that the strongest noise production region is embedded in this segment. The linearized convected wave equation breaks down in the transonic regime and a nonlinear equation must be used. Amplitude, directivity, and the transonic similarity rule for noise radiation can be established through proper matching of the displacement amplitude between the source and the wave field. This paper is a first step towards establishing a general theory of transonic aeroacoustics.

Pao, S. P.↗
NASA NTRS1985

GMRES acceleration of computational fluid dynamics codes

The generalized minimal residual algorithm (GMRES) is a conjugate-gradient like method that applies directly to nonsymmetric linear systems of equations. In this paper, GMRES is modified to handle nonlinear equations characteristic of computational fluid dynamics. Attention is devoted to the concept of preconditioning and the role it plays in assuring rapid convergence. A formulation is developed that allows GMRES to be preconditioned by the solution procedures already built into existing computer codes. Examples are provided that demonstrate the ability of GMRES to greatly improve the robustness and rate of convergence of current state-of-the-art fluid dynamics codes. Theoretical aspects of GMRES are presented that explain why it works. Finally, the advantage GMRES enjoys over related methods such as conjugate gradients are discussed.

Wigton, L. B.↗
NASA NTRS1986

The role of computerized symbolic manipulation in rotorcraft dynamics analysis

The potential role of symbolic manipulation programs in development and solution of the governing equations for rotorcraft dynamics problems is discussed and illustrated. Nonlinear equations of motion for a helicopter rotor blade represented by a rotating beam are developed making use of the computerized symbolic manipulation program MACSYMA. The use of computerized symbolic manipulation allows the analyst to concentrate on more meaningful tasks, such as establishment of physical assumptions, without being sidetracked by the tedious and trivial details of the algebraic manipulations. Furthermore, the resulting equations can be produced, if necessary, in a format suitable for numerical solution. A perturbation-type solution for the resulting dynamical equations is shown to be possible with a combination of symbolic manipulation and standard numerical techniques. This should ultimately lead to a greater physical understanding of the behavior of the solution than is possible with purely numerical techniques. The perturbation analysis of the flapping motion of a rigid rotor blade in forward flight is presented, for illustrative purposes, via computerized symbolic manipulation with a method that bypasses Floquet theory.

Crespo Da Silva, Marcelo R. M.↗
NASA NTRSDate not supplied

Wall-Modeled Large Eddy Simulation Method for Unstructured-Grid Navier-Stokes Solvers

This paper reports on the implementation and assessment of a Wall-Modeled Large-Eddy Simulation (WMLES) methodology in an unstructured-grid, node-centered flow solver, FUN3D that is developed and supported at the NASA Langley Research Center. Finite-volume (FV) and finite-element (FE) discretization schemes considered in the study provide formal second-order spatial accuracy. Large-Eddy Simulations (LES) resolve large-scale turbulent-flow features and filter out small-scale effects using the Vreman subgrid-scale model. At solid-wall boundaries, a shear-stress model is employed to provide a proper boundary-flux closure. The nonlinear equations are integrated in time using either an optimized backward difference formula or an implicit multistage Runge-Kutta temporal scheme. The implicit equations at each time step are solved by strong nonlinear iteration schemes. WMLES demonstrations are shown for two high-lift configurations, namely, the McDonnell Douglas 30P30N multielement airfoil and a NASA High-Lift Common Research Model. Results show that the WMLES approaches implemented in the FV and FE discretization methods produce consistent solutions and are capable of capturing key aerodynamic characteristics and flow structures for high-lift configurations at a wide range of angles of attack including maximum-lift conditions. In the 30P30N example, correct trends in the variations of integrated aerodynamic forces and moments, surface pressure distributions, and boundary-layer profiles are captured as the Reynolds number is increased.

CFD; turbulence modeling; High-Lift flow simulatio↗
NASA NTRS1985

A Lagrange-D'Alembert formulation of the equations of motion of a helicopter carrying an externally suspended load

The exact nonlinear equations of motion are derived for a helicopter with an extenal load suspended by fore and aft, rigid-link cables. Lagrange's form of D'Alembert's principle is used. Ten degrees of freedom are necessary to represent the motion of this system in an inertial reference frame: six for the helicopter relative to inertial space and four for the load relative to the helicopter.

Weber, J. M.↗
NASA NTRS1990

Parallel homotopy curve tracking on a hypercube

An investigation is conducted to find good parallel algorithms for solving systems of nonlinear equations using probability-one homotopy methods. Particular attention is paid to algorithms for the hypercube. Methods for one of the most computationally expensive steps of the homotopy approach, the computation of the kernel of the Jacobian matrix of the homotopy map, are studied. General nonlinear systems of equations with small and dense Jacobian matrices are considered, however, polynomial systems are not, since their structure leads to different strategies for parallelism. The mathematics behind the homotopy algorithm is summarized and the use of orthogonal factorizations is discussed. Parallel algorithms for orthogonal factorizations and triangular system solving are described. Computational results are presented and discussed.

Chakraborty, A.↗
NASA NTRS1993

The barotropic normal modes in certain shear flows and the traveling waves in the atmosphere

It is shown analytically and numerically that in certain shear flows the linearized nondivergent barotropic vorticity equation has a limited number of neutral normal modes. The latitudinal structures of these shear flows can be expressed as polynomials of the sine of latitude. The first few such shear flows resemble the gross features of the zonal winds in the atmosphere of the earth at different times and altitudes. The spatial structures of the neutral normal modes in these shear flows are spherical harmonics, and, as a consequence, these modes are also the exact solutions of the fully nonlinear equation because the nonlinear interaction term vanishes identically. The spatial structures of the observed 5-, 4-, 2-, and 16-day free traveling waves in the atmosphere are often identified with the spherical harmonics with indices of (m, n) = ( 1, 2), (2, 3), (3, 3), and ( 1, 4), which are known previously as the neutral normal modes of the nondivergent barotropic vorticity equation in a motionless background state. Our results could explain why these free traveling waves can survive the shearing effects of zonal flows that are far different from rest because these spherical harmonics are also normal modes in certain shear flows that resemble the observations of the atmosphere.

Chen, Ping↗
NASA NTRS1978

Liquid jet impingement normal to a disk in zero gravity

An experimental and analytical investigation was conducted to determine the free surface shapes of circular jets impinging normal to sharp-edged disks in zero gravity. Experiments conducted in a zero gravity drop tower yielded three distinct flow patterns which were classified in terms of the relative effects of surface tension and inertial forces. An order of magnitude analysis was conducted indicating regions where viscous forces were not significant when computing free surface shapes. The free surface analysis was simplified by transforming the governing potential flow equations and boundary conditions into the inverse plane. The resulting nonlinear equations were solved numerically and comparisons were made with the experimental data for the inertia dominated regime.

Labus, T. L.↗
NASA NTRS1974

The application of MINIQUASI to thermal program boundary and initial value problems

The feasibility of applying the solution techniques of Miniquasi to the set of equations which govern a thermoregulatory model is investigated. For solving nonlinear equations and/or boundary conditions, a Taylor Series expansion is required for linearization of both equations and boundary conditions. The solutions are iterative and in each iteration, a problem like the linear case is solved. It is shown that Miniquasi cannot be applied to the thermoregulatory model as originally planned.

Source record↗
NASA NTRS1992

A quadratic-tensor model algorithm for nonlinear least-squares problems with linear constraints

A new algorithm for solving nonlinear least-squares and nonlinear equation problems is proposed which is based on approximating the nonlinear functions using the quadratic-tensor model by Schnabel and Frank. The algorithm uses a trust region defined by a box containing the current values of the unknowns. The algorithm is found to be effective for problems with linear constraints and dense Jacobian matrices.

Hanson, R. J.↗
NASA NTRS2020

Application of Quasi-Static Modal Analysis to an Orion Multi-Purpose Crew Vehicle Test

Bolted structural joints often exhibit load-dependent stiffness and energy dissipation that leads to nonlinear, amplitude dependent frequency and damping in the structure. As an alternative to direct integration of the nonlinear equations of motion, quasi-static modal analysis (QSMA) determines the dependence of frequency and damping on response amplitude using loading behavior from nonlinear static analyses. QSMA has previously been demonstrated to substantially reduce computational cost and maintain accuracy relative to full nonlinear dynamic simulation. This work explores the applicability of QSMA to a complex, large-scale aerospace structure. QSMA is employed to analyze a nonlinear model of test hardware developed to support the Orion Multi-Purpose Crew Vehicle program, which exhibited nonlinear behavior during dynamic testing at flight-like load levels. In addition to the extraction of amplitude-dependent frequency and damping curves, a Bouc-Wen hysteresis model was used in conjunction with the quasi-static results to develop nonlinear, uncoupled, time-domain modal equations of motion for the structure. Excellent agreement was observed between the reduced and full-order nonlinear models, encouraging future employment of QSMA to support accurate and efficient model reduction of structures with bolted joint nonlinearities.

QSMA↗
NASA NTRS1988

The nonlinear dynamics of a spacecraft coupled to the vibration of a contained fluid

The dynamics of a linear spacecraft mode coupled to a nonlinear low gravity slosh of a fluid in a cylindrical tank is investigated. Coupled, nonlinear equations of motion for the fluid-spacecraft dynamics are derived through an assumed mode Lagrangian method. Unlike linear fluid slosh models, this nonlinear slosh model retains two fundamental slosh modes and three secondary modes. An approximate perturbation solution of the equations of motion indicates that the nonlinear coupled system response involves fluid-spacecraft modal resonances not predicted by either a linear, or a nonlinear, uncoupled slosh analysis. Experimental results substantiate the analytical predictions.

Peterson, Lee D.↗
NASA NTRS1994

Multilevel algorithms for nonlinear optimization

Multidisciplinary design optimization (MDO) gives rise to nonlinear optimization problems characterized by a large number of constraints that naturally occur in blocks. We propose a class of multilevel optimization methods motivated by the structure and number of constraints and by the expense of the derivative computations for MDO. The algorithms are an extension to the nonlinear programming problem of the successful class of local Brown-Brent algorithms for nonlinear equations. Our extensions allow the user to partition constraints into arbitrary blocks to fit the application, and they separately process each block and the objective function, restricted to certain subspaces. The methods use trust regions as a globalization strategy, and they have been shown to be globally convergent under reasonable assumptions. The multilevel algorithms can be applied to all classes of MDO formulations. Multilevel algorithms for solving nonlinear systems of equations are a special case of the multilevel optimization methods. In this case, they can be viewed as a trust-region globalization of the Brown-Brent class.

Alexandrov, Natalia↗