Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear problems”

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,621 records · Page 90

Difference methods for problems with different time scales

The use of difference methods for weakly nonlinear systems of ordinary differential equations with rapidly oscillating solutions is considered, and a general approach which depends only on the smoothness of the coefficients and the nonlinearities is developed. In particular, one is led to a strategy which is suitable for the detection and resolution of turning points in such systems. A computational example is presented.

Scheid, R. E., Jr.↗

Nonlinear force-free magnetic fields: Calculation and applicatin to astrophysics

The problem concerned in this work is that of calculating magnetic field configurations in which the Lorentz force (vector)j x (vector)B is everywhere zero, subject to specified boundary conditions. The magnetic field is represented in terms of Clebsch variables in the form (vector)B = del alpha x del beta. These variables are constant on any field line. The most appropriate choice of boundary conditions is to specify the values of alpha and beta on the bounding surface. It is proposed that the field lines move in the direction of local Lorentz force and relax towards a force-free field configuration. This concept leads to an iteration procedure for modifying the variables alpha and beta that tends asymptotically towards the force-free state. This method is first applied to a simple problem in two rectangular dimensions; the calculation shows that the convergence of magnetic field energy to a minimum state (force-free) is close to exponential. This method is then applied to study some astrophysical force-free magnetic fields, such as the structures and evolution of magnetic fields of rotating sunspots and accretion disks. The implication of the results, as related to the mechanisms of solar flares, extragalactic radio sources and radio jets, are discussed.

Yang, Wei-Hong↗

Aircraft design for mission performance using nonlinear multiobjective optimization methods

A new technique which converts a constrained optimization problem to an unconstrained one where conflicting figures of merit may be simultaneously considered was combined with a complex mission analysis system. The method is compared with existing single and multiobjective optimization methods. A primary benefit from this new method for multiobjective optimization is the elimination of separate optimizations for each objective, which is required by some optimization methods. A typical wide body transport aircraft is used for the comparative studies.

Dovi, Augustine R.↗

A Scalable Reduced‐Order Model for the Steady Navier–Stokes Equations

Scaling up new scientific technologies from laboratory to industry often involves demonstrating performance on a larger scale. Computer simulations can accelerate design and predictions in the deployment process, though traditional numerical methods are computationally intractable even for intermediate pilot plant scales. Recently, the component reduced order modeling method has been developed to tackle this challenge by combining projection reduced order modeling and discontinuous Galerkin domain decomposition. However, while many scientific or engineering applications involve nonlinear physics, this method has only been demonstrated for various linear systems. In this work, the component reduced order modeling method is extended to steady Navier–Stokes flow, with application to general nonlinear physics in view. The large‐scale, global domain is decomposed into a combination of small‐scale unit component. Linear subspaces for flow velocity and pressure are identified via proper orthogonal decomposition over sample snapshots collected from each small‐scale unit component. Velocity bases are augmented with a pressure supremizer to satisfy the inf–sup condition for stable pressure prediction. Two different nonlinear reduced order modeling methods are employed and compared for efficient evaluation of nonlinear advection: A third‐order tensor projection operator and the empirical quadrature procedure. The proposed method is demonstrated on the flow over arrays of five different unit objects, achieving a 23‐fold speedup with less than 4% relative error in domains up to 256 times larger than the unit components. Furthermore, a numerical experiment with the pressure supremizer strongly indicates the need for a supremizer for stable pressure prediction. A comparison between the tensorial approach and the empirical quadrature procedure revealed a slight advantage of the empirical quadrature procedure. The framework is compared with an alternating Schwarz‐based reduced‐order approach, demonstrating improved efficiency and robustness for the DG‐based global solver while retaining flexibility for sub‐scale iterative solvers. The method is further extended to a coupled advection–diffusion and Navier–Stokes system, illustrating its applicability to multi‐physics problems and its potential for more general, inter‐coupled nonlinear systems.

42 ENGINEERING↗

Small Body GN&C Research Report: A Robust Model Predictive Control Algorithm with Guaranteed Resolvability

A robustly stabilizing MPC (model predictive control) algorithm for uncertain nonlinear systems is developed that guarantees the resolvability of the associated finite-horizon optimal control problem in a receding-horizon implementation. The control consists of two components; (i) feedforward, and (ii) feedback part. Feed-forward control is obtained by online solution of a finite-horizon optimal control problem for the nominal system dynamics. The feedback control policy is designed off-line based on a bound on the uncertainty in the system model. The entire controller is shown to be robustly stabilizing with a region of attraction composed of initial states for which the finite-horizon optimal control problem is feasible. The controller design for this algorithm is demonstrated on a class of systems with uncertain nonlinear terms that have norm-bounded derivatives, and derivatives in polytopes. An illustrative numerical example is also provided.

model predictive control algorithm (MPC)↗

Neural operators for stochastic modeling of nonlinear structural system response to natural hazards

Traditionally, neural networks have been employed to learn the mapping between finite-dimensional Euclidean spaces. However, recent research has opened up new horizons, focusing on the utilization of deep neural networks to learn operators capable of mapping infinite-dimensional function spaces. Here, in this work, we employ two state-of-the-art neural operators, the deep operator network (DeepONet) and the Fourier neural operator (FNO) for the prediction of the nonlinear time history response of structural systems exposed to natural hazards, such as earthquakes and windstorms. Specifically, we propose two architectures, a self-adaptive FNO and a fast Fourier transform-based DeepONet (DeepFNOnet), where we employ a FNO beyond the DeepONet to learn the discrepancy between the ground truth and the solution predicted by the DeepONet. To demonstrate the efficiency and applicability of the architectures, two problems are considered. In the first, we use the proposed model to predict the seismic nonlinear dynamic response of a six-story shear building subject to stochastic ground motions. In the second problem, we employ the operators to predict the wind-induced nonlinear dynamic response of a high-rise building while explicitly accounting for the stochastic nature of the wind excitation. In both cases, the trained metamodels achieve high accuracy while being orders of magnitude faster than their corresponding high-fidelity models.

DeepONet↗

The dynamic aspects of thermo-elasto-viscoplastic snap-through and creep buckling phenomena

Use of a mathematical model and solution methodology, to examine dynamic buckling and dynamic postbuckling behavior of shallow arches and spherical caps made of a realistic material and undergoing non-isothermal, elasto-viscoplastic deformation was examined. Thus, geometric as well as material type nonlinearities of higher order are included in this analysis. The dynamic stability problem is studied under impulsive loading and suddenly applied loading with loads of constant magnitude and infinite duration. A finite element model was derived directly from the incrementally formulated nonlinear shell equations, by using a tensor-oriented procedure. As an example of the results, the time history of the midspan displacement of a damped shallow circular arch is presented.

Riff, R.↗

Theoretical aspects of sounder antenna operation

The fundamental problem of calculating the radiation from an electric dipole immersed in a magnetoplasma is one of solving Maxwell's equations coupled to the self-consistent field-particle equations describing the plasma. For very small applied antenna voltages, the plasma dynamics may be linearized, and the plasma reaction currents computed self-consistently in terms of the voltage on the antenna. Thus, the short dipole can be idealized as a point dipole of moment P sub o = C sub a, V sub o el, where C sub a is the antenna element-to-plasma capacitance for a sheath of dimensions roughly several plasma Debye lengths and el is the dipole element's length. In actual practice, it is at this point that the first difficulty with the antenna theory arises. Even for small signals, the sheath around such an antenna will exhibit a complex impedance. For larger signals, nonlinear effects appear. Attemps are made to solve this problem.

Fredricks, R. W.↗

Implicit Total Variation Diminishing (TVD) schemes for steady-state calculations

The novel implicit and unconditionally stable, high resolution Total Variation Diminishing (TVD) scheme whose application to steady state calculations is presently examined is a member of a one-parameter family of implicit, second-order accurate systems developed by Harten (1983) for the computation of weak solutions for one-dimensional hyperbolic conservation laws. The scheme will not generate spurious oscillations for a nonlinear scalar equation and a constant coefficient system. Numerical experiments for a quasi-one-dimensional nozzle problem show that the experimentally determined stability limit correlates exactly with the theoretical stability limit for the nonlinear scalar hyberbolic conservation laws.

Yee, H. C.↗

Postbuckling analysis using a general purpose code

A new capability for solving postbuckling problems in shell structures is described. The matrix theory to adapt Newton's method to nonlinear finite element shell analysis is outlined first. The matrix theory is directed at writing consistent linear algebratic equations for problems where the tangent stiffness matrix is singular or nearly singular. The matrix theory suggests a change of variables as part of the usual iterative procedure in Newton's method. The change of variables is shown to be feasible for introduction into the algorithm programmed in general purpose codes for finite element analysis of structures. Numerical results from a new option that has been programmed in an existing general purpose code are presented. The analysis of shell structures for collapse and for branching at bifurcation loads is illustrated by the numerical examples.

Thurston, G. A.↗

A disturbance based control/structure design algorithm

Some authors take a classical approach to the simultaneous structure/control optimization by attempting to simultaneously minimize the weighted sum of the total mass and a quadratic form, subject to all of the structural and control constraints. Here, the optimization will be based on the dynamic response of a structure to an external unknown stochastic disturbance environment. Such a response to excitation approach is common to both the structural and control design phases, and hence represents a more natural control/structure optimization strategy than relying on artificial and vague control penalties. The design objective is to find the structure and controller of minimum mass such that all the prescribed constraints are satisfied. Two alternative solution algorithms are presented which have been applied to this problem. Each algorithm handles the optimization strategy and the imposition of the nonlinear constraints in a different manner. Two controller methodologies, and their effect on the solution algorithm, will be considered. These are full state feedback and direct output feedback, although the problem formulation is not restricted solely to these forms of controller. In fact, although full state feedback is a popular choice among researchers in this field (for reasons that will become apparent), its practical application is severely limited. The controller/structure interaction is inserted by the imposition of appropriate closed-loop constraints, such as closed-loop output response and control effort constraints. Numerical results will be obtained for a representative flexible structure model to illustrate the effectiveness of the solution algorithms.

Mclaren, Mark D.↗

On the onset of three-dimensionality and time-dependence in the Goertler vortex problem

The instability of large amplitude Goertler vortices in a growing boundary layer is discussed in the fully nonlinear regime. It is shown that a three-dimensional breakdown to a flow with wavy vortex boundaries similar to that which occurs in the Taylor vortex problem takes place. However, the instability is confined to the thin shear layers which have been shown to trap the region of vortex activity. The disturbance eigenfunctions decay exponentially away from the center of these layers so that the upper and lower shear layers can support independent modes of instability. The structure of the instability, in particular its location and speed of downstream propagation, is found to be entirely consistent with recent experimental results. Furthermore, it is shown that the upper and lower layers support wavy vortex instabilities with quite different frequencies. This result is again consistent with the available experimental observations.

Hall, Philip↗

Accelerating the solution of the S N equations with highly anisotropic scattering using the Fokker-Planck approximation

The discrete ordinates method can model forward-peaked transport problems accurately. However, convergence of discrete ordinates solution can become arbitrarily slow upon use of standard iterative procedures like source iteration. Standard zeroth and first moment-based acceleration methods like nonlinear diffusion acceleration and diffusion synthetic acceleration are ineffective in accelerating such problems because these methods do not correct higher order Legendre-moments of angular flux. In this paper, we explore the idea of using the Fokker-Planck approximation as a preconditioner to accelerate the convergence of solution of forward-peaked transport problems. Finally, we observe that such preconditioning can successfully reduce the iteration count by up to four orders of magnitude, and solver runtime by up to three orders of magnitude when compared to diffusion synthetic acceleration.

42 ENGINEERING↗

A Deformed Geometry Synthesis Technique for Determining Stacking and Cryogenically Induced Preloads for the Space Launch System

The Space Launch System (SLS) stacking and Core Stage (CS) fueling induce significant preloads that contribute to the liftoff pad separation “twang”. To accurately capture this, an approach is required that can replicate the physics of all SLS physical stacking steps, CS cryogenic shrinkage, associated geometric nonlinearities, and the transient behavior and decay of the preloads with changing boundary conditions as the vehicle separates from the pad. The Deformed Geometry Synthesis (DGS) approach presented here satisfies the above requirements. DGS determines induced preloads by modeling components in their deformed geometry states and then enforcing compatibility by closing the resulting “deadbands”. DGS seamlessly integrates into the multibody modal synthesis framework and does not require the use of artificial external loads to enforce preloads or post-processing steps to remove their influence. Since DGS iterates to solve for the deformed state inclusive of geometric nonlinearities, running linearized parametrics to exercise different potential orientations of ball jointed struts that connect the CS to Boosters for cryo-shrinkage analyses is entirely avoided. Relative to the transient behavior and decay of stacking and cryo-induced preloads with SLS liftoff pad separation, this is an area of considerable interest to the SLS program. To capture this in the most accurate way possible, DGS algorithms are designed to work with Henkel-Mar nonlinear pad separation algorithms which operate on the separating longitudinal and lateral degrees of freedom (DoFs) between the vehicle and the pad. As the separating DoFs release, in whatever manner as dictated by the interface geometries, interface loads and interface flexibilities as well as the external loading on the vehicle, the subject preloads generate a complex twang/decay time-trace as dictated by the physics of the problem. This paper presents DGS numerical verification against the closed-form solution for Timoshenko’s 3 ball-jointed strut preload problem. This problem is then extended by the authors to the geometric nonlinear case where DGS is compared to the Newton-Raphson solution of the nonlinear equations. Next, DGS is utilized to solve the SLS stacking and cryogenic shrinkage coupled loads analyses. Finally, Henkel-Mar pad separation simulations are executed that isolate the impact of the induced preloads’ twang and decay characteristics.

Space Launch System↗

An inverse dynamics approach to trajectory optimization and guidance for an aerospace plane

The optimal ascent problem for an aerospace planes is formulated as an optimal inverse dynamic problem. Both minimum-fuel and minimax type of performance indices are considered. Some important features of the optimal trajectory and controls are used to construct a nonlinear feedback midcourse controller, which not only greatly simplifies the difficult constrained optimization problem and yields improved solutions, but is also suited for onboard implementation. Robust ascent guidance is obtained by using combination of feedback compensation and onboard generation of control through the inverse dynamics approach. Accurate orbital insertion can be achieved with near-optimal control of the rocket through inverse dynamics even in the presence of disturbances.

Lu, Ping↗

Cardiovascular dynamics during space sickness and deconditioning

We are currently funded by NASA for the project, 'Cardiovascular Dynamics During Space Sickness and Deconditioning.' NASA has given priority to the investigation of two problems encountered in the long-term space flights currently being planned: (1) space motion sickness; and (2) cardiovascular deconditioning. We have proposed to use spectral and nonlinear dynamical analysis of heart rate data to quantify the presence of these problems and to evaluate countermeasures against them.

Goldberger, Ary L.↗

A finite element solution of transonic flow

The use of finite elements is explored in a field in which its use has previously not been deemed very feasible, that of transonic flow. The specific problem chosen is that of steady small-disturbance transonic flow. The nonlinear equations are formulated with an artificial viscosity term added to yield the proper domain of dependence and directional bias in supersonic regions and across imbedded shock waves. Justification is given for the problem and means of solution chosen, and the potential advantages of the finite element procedure over standard finite difference procedures are discussed. Several possible improvements on the method as presently derived are stated. Computational mesh requirements and certain mesh variations are described. Some results equivalent to finite difference calculations are given as a sample solution.

Tatum, K. E.↗

Constructing space difference schemes which satisfy a cell entropy inequality

A numerical methodology for solving convection problems is presented, using finite difference schemes which satisfy the second law of thermodynamics on a cell-by-cell basis in addition to the usual conservation laws. It is shown that satisfaction of a cell entropy inequality is sufficient, in some cases, to guarantee nonlinear stability. Some details are given for several one-dimensional problems, including the quasi-one-dimensional Euler equations applied to flow in a nozzle.

Merriam, Marshal L.↗