Search NASA⌕ Search

SEARCH · Search NASA

Results for “System of 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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Group-kinetic theory of turbulence

The two phases are governed by two coupled systems of Navier-Stokes equations. The couplings are nonlinear. These equations describe the microdynamical state of turbulence, and are transformed into a master equation. By scaling, a kinetic hierarchy is generated in the form of groups, representing the spectral evolution, the diffusivity and the relaxation. The loss of memory in formulating the relaxation yields the closure. The network of sub-distributions that participates in the relaxation is simulated by a self-consistent porous medium, so that the average effect on the diffusivity is to make it approach equilibrium. The kinetic equation of turbulence is derived. The method of moments reverts it to the continuum. The equation of spectral evolution is obtained and the transport properties are calculated. In inertia turbulence, the Kolmogoroff law for weak coupling and the spectrum for the strong coupling are found. As the fluid analog, the nonlinear Schrodinger equation has a driving force in the form of emission of solitons by velocity fluctuations, and is used to describe the microdynamical state of turbulence. In order for the emission together with the modulation to participate in the transport processes, the non-homogeneous Schrodinger equation is transformed into a homogeneous master equation. By group-scaling, the master equation is decomposed into a system of transport equations, replacing the Bogoliubov system of equations of many-particle distributions. It is in the relaxation that the memory is lost when the ensemble of higher-order distributions is simulated by an effective porous medium. The closure is thus found. The kinetic equation is derived and transformed into the equation of spectral flow.

Tchen, C. M.↗
NASA NTRS1993

Optimal fixed-finite-dimensional compensator for Burgers' equation with unbounded input/output operators

The problem of using reduced order dynamic compensators to control a class of nonlinear parabolic distributed parameter systems was considered. Concentration was on a system with unbounded input and output operators governed by Burgers' equation. A linearized model was used to compute low-order-finite-dimensional control laws by minimizing certain energy functionals. Then these laws were applied to the nonlinear model. Standard approaches to this problem employ model/controller reduction techniques in conjunction with linear quadratic Gaussian (LQG) theory. The approach used is based on the finite dimensional Bernstein/Hyland optimal projection theory which yields a fixed-finite-order controller.

Burns, John A.↗
NASA NTRS1974

Nonlinear equations of motion for the elastic bending and torsion of twisted nonuniform rotor blades

The equations of motion are developed by two complementary methods, Hamilton's principle and the Newtonian method. The resulting equations are valid to second order for long, straight, slender, homogeneous, isotropic beams undergoing moderate displacements. The ordering scheme is based on the restriction that squares of the bending slopes, the torsion deformation, and the chord/radius and thickness/radius ratios are negligible with respect to unity. All remaining nonlinear terms are retained. The equations are valid for beams with mass centroid axis and area centroid (tension) axis offsets from the elastic axis, nonuniform mass and stiffness section properties, variable pretwist, and a small precone angle. The strain-displacement relations are developed from an exact transformation between the deformed and undeformed coordinate systems. These nonlinear relations form an important contribution to the final equations. Several nonlinear structural and inertial terms in the final equations are identified that can substantially influence the aeroelastic stability and response of hingeless helicopter rotor blades.

Hodges, D. H.↗
DOE OSTI2021

A Segregated Approach for Modeling the Electrochemistry in the 3-D Microstructure of Li-Ion Batteries and Its Acceleration Using Block Preconditioners

Abstract Battery performance is strongly correlated with electrode microstructure. Electrode materials for lithium-ion batteries have complex microstructure geometries that require millions of degrees of freedom to solve the electrochemical system at the microstructure scale. A fast-iterative solver with an appropriate preconditioner is then required to simulate large representative volume in a reasonable time. In this work, a finite element electrochemical model is developed to resolve the concentration and potential within the electrode active materials and the electrolyte domains at the microstructure scale, with an emphasis on numerical stability and scaling performances. The block Gauss-Seidel (BGS) numerical method is implemented because the system of equations within the electrodes is coupled only through the nonlinear Butler–Volmer equation, which governs the electrochemical reaction at the interface between the domains. The best solution strategy found in this work consists of splitting the system into two blocks—one for the concentration and one for the potential field—and then performing block generalized minimal residual preconditioned with algebraic multigrid, using the FEniCS and the Portable, Extensible Toolkit for Scientific Computation libraries. Significant improvements in terms of time to solution (six times faster) and memory usage (halving) are achieved compared with the MUltifrontal Massively Parallel sparse direct Solver. Additionally, BGS experiences decent strong parallel scaling within the electrode domains. Last, the system of equations is modified to specifically address numerical instability induced by electrolyte depletion, which is particularly valuable for simulating fast-charge scenarios relevant for automotive application.

25 ENERGY STORAGE↗
NASA NTRS1992

Invariant two-dimensional structures of vorticity field in nonviscous fluid

It is shown that a fluid flow which are invariant under translations along their vorticity field directions can be described in terms of a formally 2D system. Its nonviscous behavior of the scalar vorticity is governed by Hamiltonian dynamics, similar to the 2D line vortex case. Applying a statistical mechanical method, the equilibrium distribution of the vorticity field for such a general 2D system is shown to obey a nonlinear partial differential equation that is a generalization of the sinh-Poisson equation derived by Montgomery and Joyce (1973, 1974). Thus, if the domain containing the fluid is finite, similar negative temperature states may occur and the distribution of the vorticity intensity may become localized in space. The results are believed to have applications in studying large-scale partial and temporal properties of flows having small but nonzero viscosity. The analysis can probably be extended to other fluid systems with curved magnetic field lines.

Chen, Hudong↗
NASA NTRS1973

Continuous and discrete describing function analysis of the LST system

A describing function of the control moment gyros (CMG) frictional nonlinearity is derived using the analytic torque equation. Computer simulation of the simplified Large Space Telescope (LST) system with the analytic torque expression is discussed along with the transfer functions of the sampled-data LST system, and the discrete describing function of the GMC frictionality.

Kuo, B. C.↗
NASA NTRS1985

DEAN: A program for dynamic engine analysis

The Dynamic Engine Analysis program, DEAN, is a FORTRAN code implemented on the IBM/370 mainframe at NASA Lewis Research Center for digital simulation of turbofan engine dynamics. DEAN is an interactive program which allows the user to simulate engine subsystems as well as a full engine systems with relative ease. The nonlinear first order ordinary differential equations which define the engine model may be solved by one of four integration schemes, a second order Runge-Kutta, a fourth order Runge-Kutta, an Adams Predictor-Corrector, or Gear's method for still systems. The numerical data generated by the model equations are displayed at specified intervals between which the user may choose to modify various parameters affecting the model equations and transient execution. Following the transient run, versatile graphics capabilities allow close examination of the data. DEAN's modeling procedure and capabilities are demonstrated by generating a model of simple compressor rig.

Sadler, G. G.↗
NASA NTRS1985

DEAN - A program for Dynamic Engine ANalysis

The Dynamic Engine Analysis Program, DEAN, is a FORTRAN code implemented on the IBM/370 mainframe at NASA Lewis Research Center for digital simulation of turbofan engine dynamics. DEAN is an interactive program which allows the user to simulate engine subsystems as well as full engine systems with relative ease. The nonlinear first order ordinary differential equations which define the engine model may be solved by one of four integration schemes, a second order Runge-Kutta, a fourth order Runge-Kutta, an Adams Predictor-Corrector, or Gear's method for still systems. The numerical data generated by the model equations are displayed a specified intervals between which the user may choose to modify various parameters affecting the model equations and transient execution. Following the transient run, versatile graphics capabilities allow close examination of the data. DEAN's modeling procedure and capabilities are demonstrated by generating a model of simple compressor rig.

Sadler, G. G.↗
NASA NTRS1995

Nonlinear evolution of magnetic flux ropes. 2: Finite beta plasma

In this second paper on the evolution of magnetic flux ropes we study the effects of gas pressure. We assume that the energy transport is described by a polytropic relationship and reduce the set of ideal MHD equations to a single, second-order, nonlinear, ordinary differential equation for the evolution function. For this conservative system we obtain a first integral of motion. To analyze the possible motions, we use a mechanical analogue -- a one-dimensional, nonlinear oscillator. We find that the effective potential for such an oscillator depends on two parameters: the polytropic index gamma and a dimensionless quantity kappa the latter being a function of the plasma beta, the strength of the azimuthal magnetic field relative to the axial field of the flux rope, and gamma. Through a study of this effective potential we classify all possible modes of evolution of the system. In the main body of the paper, we focus on magnetic flux ropes whose field and gas pressure increase steadily towards the symmetry axis. In this case, for gamma greater than 1 and all values of kappa, only oscillations are possible. For gamma less than 1, however, both oscillations and expansion are allowed. For gamma less than 1 and kappa below a critical value, the energy of the nonlinear oscillator determines whether the flux rope will oscillate or expand to infinity. For gamma less than 1 and kappa above critical, however, only expansion occurs. Thus by increasing kappa while keeping gamma fixed (less than 1), a phase transition occurs at kappa = kappa(sub critical) and the oscillatory mode disappears. We illustrate the above theoretical considerations by the example of a flux rope of constant field line twist evolving self-similarly. For this example, we present the full numerical MHD solution. In an appendix to the paper we catalogue all possible evolutions when (1) either the magnetic field or (2) the gas pressure decreases monotonically toward the axis. We find that in these cases critical conditions can occur for gamma greater than 1. While in most cases the flux rope collapses, there are notable exceptions when, for certain ranges of kappa and gamma, collapse may be averted.

Osherovich, V. A.↗
DOE OSTI2022

Modeling the AC Power Flow Equations with Optimally Compact Neural Networks: Application to Unit Commitment

Nonlinear power flow constraints render a variety of power system optimization problems computationally intractable. Emerging research shows, however, that the nonlinear AC power flow equations can be successfully modeled using neural networks. These neural networks can be exactly transformed into mixed integer linear programs and embedded inside challenging optimization problems, thus replacing nonlinearities that are intractable for many applications with tractable piecewise linear approximations. Such approaches, though, suffer from an explosion of the number of binary variables needed to represent the neural network. Accordingly, this paper develops a technique for training an "optimally compact'' neural network, i.e., one that can represent the power flow equations with a sufficiently high degree of accuracy while still maintaining a tractable number of binary variables. We demonstrate the use of this neural network as an approximator of the nonlinear power flow equations by embedding it in the AC unit commitment problem, transforming the problem from a mixed integer nonlinear program into a more manageable mixed integer linear program. We use the 14-, 57-, and 89-bus networks as test cases and compare the AC-feasibility of commitment decisions resulting from the neural network, DC, and linearized power flow approximations. Our results show that the neural network model outperforms both the DC and linearized power flow approximations when embedded in the unit commitment problem. The neural network formulation most often selects a feasible unit commitment schedule, and furthermore, it only s

AC power flow↗
NASA NTRS1990

Algorithm for "Bang-Bang" Control Laws

Switching times updated to minimize errors in final positions. Algorithm computes "bang-bang" control laws for single- or multiple-input systems, both those describable by linear dynamical equations and those describable by nonlinear dynamical equations that are linear in control vectors. Algorithm needed because analytical solutions of "bang-bang" control problems intractable for all but simplest systems.

Wen, John Ting-Yung↗
NASA NTRS1991

Large Angle Transient Dynamics (LATDYN) demonstration problem manual

LATDYN is a computer code for modeling the Large Angle Transient Dynamics of structures. The objective in developing the code was to investigate new techniques for analyzing flexible deformation and control/structure interaction problems associated with large angular motions of spacecraft. Such motions may consist of pointing the entire spacecraft or articulation of individual components, events which occur frequently during construction, operation, and maintenance of large spacecraft. This type of analysis is beyond the routine capability of conventional analytical tools without simplifying assumptions. In some instances, the motion may be sufficiently slow and the spacecraft (or component) sufficiently rigid to simplify the analysis of dynamics and controls by making pseudo-static and/or rigid body assumptions. LATDYN introduces a new approach to the problem by combining finite element structural analysis, multi-body dynamics, and control system analysis in a single tool. It includes a new type of finite element that can deform and rotate through large angles at the same time, and which can be connected to other finite elements either rigidly or through mechanical joints. LATDYN also provides symbolic capabilities for modeling control systems which are interfaced directly with the finite element structural model. Thus, the nonlinear equations representing the structural model are integrated along with the equations representing sensors, processing, and controls as a coupled system.

Wu, Shih-Chin↗
NASA NTRS1991

Large Angle Transient Dynamics (LATDYN) user's manual

A computer code for modeling the large angle transient dynamics (LATDYN) of structures was developed to investigate techniques for analyzing flexible deformation and control/structure interaction problems associated with large angular motions of spacecraft. This type of analysis is beyond the routine capability of conventional analytical tools without simplifying assumptions. In some instances, the motion may be sufficiently slow and the spacecraft (or component) sufficiently rigid to simplify analyses of dynamics and controls by making pseudo-static and/or rigid body assumptions. The LATDYN introduces a new approach to the problem by combining finite element structural analysis, multi-body dynamics, and control system analysis in a single tool. It includes a type of finite element that can deform and rotate through large angles at the same time, and which can be connected to other finite elements either rigidly or through mechanical joints. The LATDYN also provides symbolic capabilities for modeling control systems which are interfaced directly with the finite element structural model. Thus, the nonlinear equations representing the structural model are integrated along with the equations representing sensors, processing, and controls as a coupled system.

Abrahamson, A. Louis↗
NASA NTRS1996

Validation Tests of TASS for Application to 3-D Vortex Simulations

Direct analytical solutions can be useful in validating the Core formulation of numerical systems. In this document an exact analytical solution to the nonlinear Navier-Stokes equation is compared to the numerical results from the three-dimensional Terminal Area Simulation System (TASS). This exact solution, of which the derivation is included, is for Beltrami type flow. Direct comparison of TASS to the analytical Beltrami solution is then used in evaluating the accuracy of TASS.

Switzer, George F.↗
NASA NTRS2010

Forces Associated with Nonlinear Nonholonomic Constraint Equations

A concise method has been formulated for identifying a set of forces needed to constrain the behavior of a mechanical system, modeled as a set of particles and rigid bodies, when it is subject to motion constraints described by nonholonomic equations that are inherently nonlinear in velocity. An expression in vector form is obtained for each force; a direction is determined, together with the point of application. This result is a consequence of expressing constraint equations in terms of dot products of vectors rather than in the usual way, which is entirely in terms of scalars and matrices. The constraint forces in vector form are used together with two new analytical approaches for deriving equations governing motion of a system subject to such constraints. If constraint forces are of interest they can be brought into evidence in explicit dynamical equations by employing the well-known nonholonomic partial velocities associated with Kane's method; if they are not of interest, equations can be formed instead with the aid of vectors introduced here as nonholonomic partial accelerations. When the analyst requires only the latter, smaller set of equations, they can be formed directly; it is not necessary to expend the labor to form the former, larger set first and subsequently perform matrix multiplications.

Roithmayr, Carlos M.↗
DOE OSTI2023

DRIPS: A framework for dimension reduction and interpolation in parameter space

Reduced-order models are often used to describe the behavior of complex systems, whose simulation with a full model is too expensive, or to extract salient features from the full model’s output. We introduce a new model-reduction framework DRIPS (dimension reduction and interpolation in parameter space) that combines the offline local model reduction with the online parameter interpolation of reduced-order bases (ROBs). The offline step of this framework relies on dynamic mode decomposition (DMD) to build a low-rank linear surrogate model, equipped with a local ROB, for quantities of interest derived from the training data generated by repeatedly solving the (nonlinear) high-fidelity model for multiple parameter points. The online step consists of the construction of a parametric reduced-order model for each target/test point in the parameter space, with the interpolation of ROBs done on a Grassman manifold and the interpolation of reduced-order operators done on a matrix manifold. The DMD component enables DRIPS to model (typically low-dimensional) quantities of interest directly, without having to access the (typically high-dimensional and possibly nonlinear) operators in a high-fidelity model that governs the dynamics of the underlying high-dimensional state variables, as required in projection-based reduced-order modeling. A series of numerical experiments suggests that DRIPS yields a model reduction, which is computationally more efficient than the commonly used projection-based proper orthogonal decomposition; it does so without requiring a prior knowledge of the governing equation for quantities of interest. Furthermore, for the nonlinear systems considered, DRIPS is more accurate than Gaussian-process interpolation (Kriging).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1977

Dynamic response and stability analysis of flexible, multibody systems

A general version of Lagrange's equations, including auxiliary nonholonomic, rheonomic conditions of constraint, is used in the dynamic simulation and stability analysis of interconnected flexible bodies. Modeling of the nonlinear flexible/rigid dynamic coupling effects, the interaction forces/torques, and the elastic deformation effects is discussed. A digital computer program is developed to obtain time-domain solution for the nonlinear response of systems represented as a collection of individual bodies, numerical linearization of system-governing equations, time-domain solution for the perturbation response about a nominal state, and a frequency-domain stability analysis corresponding to the linearization. The digital simulation code is employed to study the dynamic behavior of a typical satellite and a spacecraft with deployable experiment booms.

Bodley, C. S.↗