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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 523 records · Page 29

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Effect of out-of-roundness on the performance of a diesel engine connecting-rod bearing

In this paper, the dynamic performance of the Ruston and Hornsby VEB diesel engine connecting-rod bearing with circular and out-of-round profiles is analyzed. The effect of cavitation is considered by using a cavitation algorithm, which mimics JFO boundary conditions. The effect of mass inertia is accounted for by solving coupled nonlinear equations of motion. The journal profiles considered are circular, elliptical, semi-elliptical, and three lobe epicycloid. The predicted journal trajectory and other performance parameters for one complete load cycle are presented for all of the out-of-round profiles and are also compared with the predictions for the circular bearing.

Vijayaraghavan, D.↗

NASA NTRS1988

Stability analysis and control of tethered satellites

Issues of stability and stabilization of the Shuttle/Tethered Satellite System (TSS) are considered. The basic, nonlinear equations of motion of the TSS are derived using the system Lagrangian. A tension control strategy in feedback form is introduced to stabilize the TSS during station-keeping. In addition, the stability properties of a constant angle control method for deployment and retrieval are investigated using Liapunov stability theory.

Liaw, D.-C.↗

NASA NTRS1990

Bifurcations and sensitivity in parametric nonlinear programming

The parametric nonlinear programming problem is that of determining the behavior of solution(s) as a parameter or vector of parameters alpha belonging to R(sup r) varies over a region of interest for the problem: Minimize over x the set f(x, alpha):h(x, alpha) = 0, g(x, alpha) is greater than or equal to 0, where f:R(sup (n+r)) approaches R, h:R(sup (n+r)) approaches R(sup q) and g:R(sup (n+r)) approaches R(sup p) are assumed to be at least twice continuously differentiable. Some of these parameters may be fixed but not known precisely and others may be varied to enhance the performance of the system. In both cases a fundamentally important problem in the investigation of global sensitivity of the system is to determine the stability boundaries of the regions in parameter space which define regions of qualitatively similar solutions. The objective is to explain how numerical continuation and bifurcation techniques can be used to investigate the parametric nonlinear programming problem in a global sense. Thus, first the problem is converted to a closed system of parameterized nonlinear equations whose solution set contains all local minimizers of the original problem. This system, which will be represented as F(z,alpha) = O, will include all Karush-Kuhn-Tucker and Fritz John points, both feasible and infeasible solutions, and relative minima, maxima, and saddle points of the problem. The local existence and uniqueness of a solution path (z(alpha), alpha) of this system as well as the solution type persist as long as a singularity in the Jacobian D(sub z)F(z,alpha) is not encountered. Thus the nonsingularity of this Jacobian is characterized in terms of conditions on the problem itself. Then, a class of efficient predictor-corrector continuation procedures for tracing solution paths of the system F(z,alpha) = O which are tailored specifically to the parametric programming problem are described. Finally, these procedures and the obtained information are illustrated within the context of design optimization.

Lundberg, Bruce N.↗

NASA NTRS1993

Canonical forms of multidimensional steady inviscid flows

Canonical forms and canonical variables for inviscid flow problems are derived. In these forms the components of the system governed by different types of operators (elliptic and hyperbolic) are separated. Both the incompressible and compressible cases are analyzed, and their similarities and differences are discussed. The canonical forms obtained are block upper triangular operator form in which the elliptic and non-elliptic parts reside in different blocks. The full nonlinear equations are treated without using any linearization process. This form enables a better analysis of the equations as well as better numerical treatment. These forms are the analog of the decomposition of the one dimensional Euler equations into characteristic directions and Riemann invariants.

Taasan, Shlomo↗

NASA NTRS1994

A parametric sensitivity study for single-stage-to-orbit hypersonic vehicles using trajectory optimization

The class of hypersonic vehicle configurations with single stage-to-orbit (SSTO) capability reflect highly integrated airframe and propulsion systems. These designs are also known to exhibit a large degree of interaction between the airframe and engine dynamics. Consequently, even simplified hypersonic models are characterized by tightly coupled nonlinear equations of motion. In addition, hypersonic SSTO vehicles present a major system design challenge; the vehicle's overall mission performance is a function of its subsystem efficiencies including structural, aerodynamic, propulsive, and operational. Further, all subsystem efficiencies are interrelated, hence, independent optimization of the subsystems is not likely to lead to an optimum design. Thus, it is desired to know the effect of various subsystem efficiencies on overall mission performance. For the purposes of this analysis, mission performance will be measured in terms of the payload weight inserted into orbit. In this report, a trajectory optimization problem is formulated for a generic hypersonic lifting body for a specified orbit-injection mission. A solution method is outlined, and results are detailed for the generic vehicle, referred to as the baseline model. After evaluating the performance of the baseline model, a sensitivity study is presented to determine the effect of various subsystem efficiencies on mission performance. This consists of performing a parametric analysis of the basic design parameters, generating a matrix of configurations, and determining the mission performance of each configuration. Also, the performance loss due to constraining the total head load experienced by the vehicle is evaluated. The key results from this analysis include the formulation of the sizing problem for this vehicle class using trajectory optimization, characteristics of the optimal trajectories, and the subsystem design sensitivities.

Lovell, T. Alan↗

NASA NTRS1994

A Procedure for 3-D Contact Stress Analysis of Spiral Bevel Gears

Contact stress distribution of spiral bevel gears using nonlinear finite element static analysis is presented. Procedures have been developed to solve the nonlinear equations that identify the gear and pinion surface coordinates based on the kinematics of the cutting process and orientate the pinion and the gear in space to mesh with each other. Contact is simulated by connecting GAP elements along the intersection of a line from each pinion point (parallel to the normal at the contact point) with the gear surface. A three dimensional model with four gear teeth and three pinion teeth is used to determine the contact stresses at two different contact positions in a spiral bevel gearset. A summary of the elliptical contact stress distribution is given. This information will be helpful to helicopter and aircraft transmission designers who need to minimize weight of the transmission and maximize reliability.

Kumar, A.↗

NASA NTRS1994

Computation of iodine species concentrations in water

During an evaluation of the use of iodine as a water disinfectant and the development of methods for measuring various iodine species in water onboard Space Freedom, it became necessary to compute the concentration of the various species based on equilibrium principles alone. Of particular concern was the case when various amounts of iodine, iodide, strong acid, and strong base are added to water. Such solutions can be used to evaluate the performance of various monitoring methods being considered. The authors of this paper present an overview of aqueous iodine chemistry, a set of nonlinear equations which can be used to model the above case, and a computer program for solving this system of equations using the Newton-Raphson method. The program was validated by comparing results over a range of concentrations and pH values with those previously presented by Gottardi for a given pH. Use of this program indicated that there are multiple roots to many cases and selecting an appropriate initial guess is important. Comparison of program results with laboratory results for the case when only iodine is added to water indicates the program gives high pH values for the iodine concentrations normally used for water disinfection. Extending the model to include the effects of iodate formation results in the computer pH values being closer to those observed, but the model with iodate does not agree well for the case in which base is added in addition to iodine to raise the pH. Potential explanations include failure to obtain equilibrium conditions in the lab, inaccuracies in published values for the equilibrium constants, and inadequate model of iodine chemistry and/or the lack of adequate analytical methods for measuring the various iodine species in water.

Schultz, John R.↗

NASA NTRS1994

A combined algorithm for minimum time slewing of flexible spacecraft

The use of Pontryagin's Maximum Principle for the large-angle slewing of large flexible structures usually results in the so-called two-point boundary-value problem (TPBVP), in which many requirements (e.g., minimum time, small flexible amplitude, and limited control powers, etc.) must be satisfied simultaneously. The successful solution of this problem depends largely on the use of an efficient numerical computational algorithm. There are many candidate algorithms available for his problem (e.g., quasilinearization, gradient, and shooting, etc.). In this paper, a proposed algorithm, which combines the quasilinearization method with a time shortening technique and a shooting method, is applied to the minimum-time, three-dimensional, and large-angle maneuver of flexible spacecraft, particularly the orbiting Spacecraft Control Laboratory Experiment (SCOLE) configuration. Theoretically, the nonlinear TPBVP can be solved only through the shooting method to find the 'exact' switching times for the bang-bang controls. However, computationally, a suitable guess for the missing initial costates is crucial because the convergence range of the unknown initial costates is usually narrow, especially for systems with high dimensions and when a multi-bang-bang control strategy is needed. On the other hand, the problems of near minimum time attitude maneuver of general rigid spacecraft and fast slewing of flexible spacecraft have been examined by the authors through a numerical approach based on the quasilinearization algorithm with a time shortening technique. Computational results have demonstrated its broad convergence range and insensitivity to initial costate choices. Consequently, a combined approach is naturally suggested here to solve the minimum time slewing problem. That is, in the computational process, the quasilinearization method is used first to obtain a near minimum time solution. Then, the acquired converged initial costates from the quasilinearization approach are transformed (tailored) to and used as the initial costate guess for starting the shooting method. Finally, the shooting method takes over the remaining calculations until the minimum-time solution converges. The nonlinear equations of motion of the SCOLE are formulated by using Lagrange's equations, with the mast modeled as a continuous beam subject to three-dimensional deformations. The numerical results will be presented and some related computational issues will also be discussed.

Bainum, P. M.↗

NASA NTRS1994

Hierarchically Parallelized Constrained Nonlinear Solvers with Automated Substructuring

This paper develops a parallelizable multilevel multiple constrained nonlinear equation solver. The substructuring process is automated to yield appropriately balanced partitioning of each succeeding level. Due to the generality of the procedure,_sequential, as well as partially and fully parallel environments can be handled. This includes both single and multiprocessor assignment per individual partition. Several benchmark examples are presented. These illustrate the robustness of the procedure as well as its capability to yield significant reductions in memory utilization and calculational effort due both to updating and inversion.

Padovan, Joe↗

NASA NTRS1992

Aeroelastic modeling of composite rotor blades with straight and swept tips

This paper presents an analytical study of the aeroelastic behavior of composite rotor blades with straight and swept tips. The blade is modeled by beam type finite elements. A single finite element is used to model the swept tip. The nonlinear equations of motion for the FEM are derived using Hamilton's principle and based on a moderate deflection theory and accounts for: arbitrary cross-sectional shape, pretwist, generally anisotropic material behavior, transverse shears and out-of-plane warping. Numerical results illustrating the effects of tip sweep, anhedral and composite ply orientation on blade aeroelastic behavior are presented. It is shown that composite ply orientation has a substantial effect on blade stability. At low thrust conditions, certain ply orientations can cause instability in the lag mode. The flap-torsion coupling associated with tip sweep can also induce aeroelastic instability in the blade. This instability can be removed by appropriate ply orientation in the composite construction. These results illustrate the inherent potential for aeroelastic tailoring present in composite rotor blades with swept tips, which still remains to be exploited in the design process.

Yuan, Kuo-An↗

NASA NTRS1995

Aeroelasticity and structural optimization of composite helicopter rotor blades with swept tips

This report describes the development of an aeroelastic analysis capability for composite helicopter rotor blades with straight and swept tips, and its application to the simulation of helicopter vibration reduction through structural optimization. A new aeroelastic model is developed in this study which is suitable for composite rotor blades with swept tips in hover and in forward flight. The hingeless blade is modeled by beam type finite elements. A single finite element is used to model the swept tip. Arbitrary cross-sectional shape, generally anisotropic material behavior, transverse shears and out-of-plane warping are included in the blade model. The nonlinear equations of motion, derived using Hamilton's principle, are based on a moderate deflection theory. Composite blade cross-sectbnal properties are calculated by a separate linear, two-dimensional cross section analysis. The aerodynamic loads are obtained from quasi-steady, incompressible aerodynamics, based on an implicit formulation. The trim and steady state blade aeroelastic response are solved in a fully coupled manner. In forward flight, where the blade equations of motion are periodic, the coupled trim-aeroelastic response solution is obtained from the harmonic balance method. Subsequently, the periodic system is linearized about the steady state response, and its stability is determined from Floquet theory.

Yuan, K. A.↗

NASA NTRS1995

A description of the solar wind-magnetosphere coupling based on nonlinear filters

A nonlinear filtering method is introduced for the study of the solar wind -- magnetosphere coupling and related to earlier linear techniques. The filters are derived from the magnetospheric state, a representation of the magnetospheric conditions in terms of a few global variables, here the auroral electrojet indices. The filters also couple to the input, a representation of the solar wind variables, here the rectified electric field. Filter-based iterative prediction of the indices has been obtained for up to 20 hours. The prediction is stable with respect to perturbations in the initial magnetospheric state; these decrease exponentially at the rate of 30/min. The performance of the method is examined for a wide range of parameters and is superior to that of other linear and nonlinear techniques. In the magnetospheric state representation the coupling is modeled as a small number of nonlinear equations under a time-dependent input.

Vassiliadis, D.↗

NASA NTRS1996

Land Surface Temperature Measurements form EOS MODIS Data

We have developed a physics-based land-surface temperature (LST) algorithm for simultaneously retrieving surface band-averaged emissivities and temperatures from day/night pairs of MODIS (Moderate Resolution Imaging Spectroradiometer) data in seven thermal infrared bands. The set of 14 nonlinear equations in the algorithm is solved with the statistical regression method and the least-squares fit method. This new LST algorithm was tested with simulated MODIS data for 80 sets of band-averaged emissivities calculated from published spectral data of terrestrial materials in wide ranges of atmospheric and surface temperature conditions. Comprehensive sensitivity and error analysis has been made to evaluate the performance of the new LST algorithm and its dependence on variations in surface emissivity and temperature, upon atmospheric conditions, as well as the noise-equivalent temperature difference (NE(Delta)T) and calibration accuracy specifications of the MODIS instrument. In cases with a systematic calibration error of 0.5%, the standard deviations of errors in retrieved surface daytime and nighttime temperatures fall between 0.4-0.5 K over a wide range of surface temperatures for mid-latitude summer conditions. The standard deviations of errors in retrieved emissivities in bands 31 and 32 (in the 10-12.5 micrometer IR spectral window region) are 0.009, and the maximum error in retrieved LST values falls between 2-3 K. Several issues related to the day/night LST algorithm (uncertainties in the day/night registration and in surface emissivity changes caused by dew occurrence, and the cloud cover) have been investigated. The LST algorithms have been validated with MODIS Airborne Simulator (MAS) dada and ground-based measurement data in two field campaigns conducted in Railroad Valley playa, NV in 1995 and 1996. The MODIS LST version 1 software has been delivered.

Wan, Zhengming↗

NASA NTRS1961

The Effect of Solar Radiation Pressure on the Motion of an Artificial Satellite

The effects of solar radiation pressure on the motion of an artificial satellite are obtained, including the effects of the intermittent acceleration which results from the eclipsing of the satellite by the earth. Vectorial methods have been utilized to obtain the nonlinear equations describing the motion, and the method of Kryloff-Bogoliuboff has been applied in their solution.

Bryant, Robert W.↗

NASA NTRS2000

Estimation of Handling Qualities Parameters of the Tu-144 Supersonic Transport Aircraft from Flight Test Data

Low order equivalent system (LOES) models for the Tu-144 supersonic transport aircraft were identified from flight test data. The mathematical models were given in terms of transfer functions with a time delay by the military standard MIL-STD-1797A, "Flying Qualities of Piloted Aircraft," and the handling qualities were predicted from the estimated transfer function coefficients. The coefficients and the time delay in the transfer functions were estimated using a nonlinear equation error formulation in the frequency domain. Flight test data from pitch, roll, and yaw frequency sweeps at various flight conditions were used for parameter estimation. Flight test results are presented in terms of the estimated parameter values, their standard errors, and output fits in the time domain. Data from doublet maneuvers at the same flight conditions were used to assess the predictive capabilities of the identified models. The identified transfer function models fit the measured data well and demonstrated good prediction capabilities. The Tu-144 was predicted to be between level 2 and 3 for all longitudinal maneuvers and level I for all lateral maneuvers. High estimates of the equivalent time delay in the transfer function model caused the poor longitudinal rating.

Curry, Timothy J.↗

NASA NTRS2000

Modeling and Analysis of the Reverse Water Gas Shift Process for In-Situ Propellant Production

This report focuses on the development of mathematical models and simulation tools developed for the Reverse Water Gas Shift (RWGS) process. This process is a candidate technology for oxygen production on Mars under the In-Situ Propellant Production (ISPP) project. An analysis of the RWGS process was performed using a material balance for the system. The material balance is very complex due to the downstream separations and subsequent recycle inherent with the process. A numerical simulation was developed for the RWGS process to provide a tool for analysis and optimization of experimental hardware, which will be constructed later this year at Kennedy Space Center (KSC). Attempts to solve the material balance for the system, which can be defined by 27 nonlinear equations, initially failed. A convergence scheme was developed which led to successful solution of the material balance, however the simplified equations used for the gas separation membrane were found insufficient. Additional more rigorous models were successfully developed and solved for the membrane separation. Sample results from these models are included in this report, with recommendations for experimental work needed for model validation.

Whitlow, Jonathan E.↗

NASA NTRS1995

Evolution of Vapor Bubbles Nucleation Sites in Low Gravity

When liquid is expelled by a vapor bubble growing at a nucleation site on a superheated surface, a thin microlayer underneath the bubble is left behind. It is evaporated from the free microlayer surface that provides for bubble growth. The average thickness of the microlayer determining the evaporation rate increases with time if the latter does not exceed a threshold value associated with the burn-out crisis. The bubble is described as a spherical segment with its flattened part adjoining the microlayer. This introduces two independent variables - the radius of the spherical part of the bubble surface and the polar angle that defines the relative area of the flattened part. They are to be found out from a set of two strongly nonlinear equations resulting from mass and momentum conservation laws. The first one depends on both microlayer thickness and nonmonotonously changing bubble base area. The second involves two major factors favoring bubble detachment - the buoyancy and a force due to the initial momentum of vapor input into the bubble. The former force depends on gravity whereas the latter one does not. It is why the limiting regimes of bubble evolution that correspond to normal or moderately reduced gravity and to microgravity feature drastically different properties. In the first case, the buoyancy dominates and the bubble evolves in such a manner as to become a full sphere at a moment that can be viewed as that of detachment. The detachment volume grows as gravity decreases. In the second case, the buoyancy is negligible and the bubble stays near the surface, while its volume continues to increase for a sufficiently long time. The findings are discussed in connection with experimental data obtained under different gravity conditions, some unpublished experiments being included. They help to understand why the pool boiling heat transfer coefficient frequently increases as gravity falls down and eventually vanishes.

Buyevich, Yu A.↗

NASA NTRS2003

Integrated Power and Attitude Control for a Spacecraft with Flywheels and Control Moment Gyroscopes

A law is designed for simultaneous control of the orientation of an Earth-pointing spacecraft, the energy stored by counter-rotating flywheels, and the angular momentum of the flywheels and control moment gyroscopes used together as all integrated set of actuators for attitude control. General. nonlinear equations of motion are presented in vector-dyadic form, and used to obtain approximate expressions which are then linearized in preparation for design of control laws that include feedback of flywheel kinetic energy error as it means of compensating for damping exerted by rotor bearings. Two flywheel 'steering laws' are developed such that torque commanded by all attitude control law is achieved while energy is stored or discharged at the required rate. Using the International Space Station as an example, numerical simulations are performed to demonstrate control about a torque equilibrium attitude and illustrate the benefits of kinetic energy error feedback.

Roithmayr, Carlos M.↗