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,747 records · Page 97
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

An algorithm for retrieving water vapor profiles in clear and cloudy atmospheres from 183 GHz radiometric measurements - Simulation studies

The latent heat represented by atmospheric water vapor is extremely important to the energetics of the earth system. Future satellites (NOAA and DMSP) will carry microwave radiometers designed to measure the profile of water vapor globally. The problem of retrieving water vapor from the measurements is highly nonlinear even in clear atmospheres and the addition of clouds only makes it more so. In this paper, an algorithm with several novel features, which will retrieve water vapor profiles from microwave radiometric measurements even in the presence of clouds, is developed. Simulations with this algorithm show a vertical resolution on the order of 3 km and that clouds are well handled in many, but not all, circumstances. The most surprising result is that clouds can actually improve the vertical resolution of the retrieval.

Wilheit, Thomas T.↗
NASA NTRS1990

Computation of convective flow with gravity modulation in rectangular cavities

In this work, a computational study is presented for the investigation of gravity modulation (g-jitter) effects in thermally driven cavity flows at terrestrial and microgravity environments. The two-dimensional, time-dependent Navier-Stokes equations are numerically integrated by a time-split method using direct matrix solvers. Computations at terrestrial gravity are utilized to assess the effects of adiabatic side-wall boundary conditions as well as the full nonlinearity of the governing equations on the sinusoidally forced Benard problem studied by Gresho and Sani. The low-g calculations focus on the establishment of critical frequency ranges and consider the effects of modulation direction and randomness. The applicability of linear analysis in the excitable frequency range at low g is also discussed.

Biringen, S.↗
NASA NTRS2008

Multiplexed Predictive Control of a Large Commercial Turbofan Engine

Model predictive control is a strategy well-suited to handle the highly complex, nonlinear, uncertain, and constrained dynamics involved in aircraft engine control problems. However, it has thus far been infeasible to implement model predictive control in engine control applications, because of the combination of model complexity and the time allotted for the control update calculation. In this paper, a multiplexed implementation is proposed that dramatically reduces the computational burden of the quadratic programming optimization that must be solved online as part of the model-predictive-control algorithm. Actuator updates are calculated sequentially and cyclically in a multiplexed implementation, as opposed to the simultaneous optimization taking place in conventional model predictive control. Theoretical aspects are discussed based on a nominal model, and actual computational savings are demonstrated using a realistic commercial engine model.

Richter, hanz↗
NASA NTRS2015

Entropy Stable Staggered Grid Spectral Collocation for the Burgers' and Compressible Navier-Stokes Equations

Staggered grid, entropy stable discontinuous spectral collocation operators of any order are developed for Burgers' and the compressible Navier-Stokes equations on unstructured hexahedral elements. This generalization of previous entropy stable spectral collocation work [1, 2], extends the applicable set of points from tensor product, Legendre-Gauss-Lobatto (LGL) to a combination of tensor product Legendre-Gauss (LG) and LGL points. The new semi-discrete operators discretely conserve mass, momentum, energy and satisfy a mathematical entropy inequality for both Burgers' and the compressible Navier-Stokes equations in three spatial dimensions. They are valid for smooth as well as discontinuous flows. The staggered LG and conventional LGL point formulations are compared on several challenging test problems. The staggered LG operators are significantly more accurate, although more costly to implement. The LG and LGL operators exhibit similar robustness, as is demonstrated using test problems known to be problematic for operators that lack a nonlinearly stability proof for the compressible Navier-Stokes equations (e.g., discontinuous Galerkin, spectral difference, or flux reconstruction operators).

Carpenter, Mark H.↗
NASA NTRS2019

Modeling and simulation challenges for softgoods in space deployable structures

This paper examines modeling and simulation challenges for deployable spacecraft structures with softgood components. Previous studies highlighted the numerical challenges due to the mixture of high stiffness and low stiffness contact in this class of structures. This paper examines the use of the commercial LS-DYNA nonlinear finite element code on a set of relevant benchmark problems.

Mobrem, Mehran↗
NASA NTRS1976

Stability of flow of a thermoviscoelastic fluid between rotating coaxial circular cylinders

The stability problem of thermoviscoelastic fluid flow between rotating coaxial cylinders is investigated using nonlinear thermoviscoelastic constitutive equations due to Eringen and Koh. The velocity field is found to be identical with that of the classical viscous case and the case of the viscoelastic fluid, but the temperature and pressure fields are found to be different. By imposing some physically reasonable mechanical and geometrical restrictions on the flow, and by a suitable mathematical analysis, the problem is reduced to a characteristic value problem. The resulting problem is solved and stability criteria are obtained in terms of critical Taylor numbers. In general, it is found that thermoviscoelastic fluids are more stable than classical viscous fluids and viscoinelastic fluids under similar conditions.

Ghandour, N. N.↗
NASA NTRS1990

A parallel algorithm for nonlinear convection-diffusion equations

A parallel algorithm for the efficient solution of nonlinear time-dependent convection-diffusion equations with small parameter on the diffusion term is presented. The method is based on a physically motivated domain decomposition that is dictated by singular perturbation analysis. The analysis is used to determine regions where certain reduced equations may be solved in place of the full equation. The method is suitable for the solution of problems arising in the simulation of fluid dynamics. Experimental results for a nonlinear equation in two-dimensions are presented.

Scroggs, Jeffrey S.↗
DOE OSTI2025

Observation of Joule–Thomson photon-gas expansion

In recent years, a self-consistent optical thermodynamic framework has emerged that offers a systematic methodology to understand, harness and exploit the complex collective dynamics of multimode nonlinear systems. These developments now allow consideration of a series of longstanding problems in optics, including the prospect of funnelling the entire power flowing in a multimode system into its ground state, for which no methodology currently exists. Here, we demonstrate an all-optical Joule-Thomson expansion process mediated by photon-photon interactions whereby the temperature of the optical gas drops abruptly to zero. Our experiments in various configurations of coupled multicore nonlinear waveguide arrangements illustrate how light undergoing expansion-induced cooling can be channelled from arbitrary input states into the fundamental mode with near-unity efficiency. We show that the stability of the post-expansion state is ensured through an irreversible process of energy conversion. The all-optical thermodynamic phenomena explored in this study may enable innovative techniques where various uncorrelated but identical sources are merged into a unified spatially coherent state, offering a route for direct beam combining.

Kirsch, Marco S. (ORCID:0000000342798228)↗
NASA NTRS1991

Hierarchic Extensions in the Static and Dynamic Analysis of Elastic Beams

Approximate solutions of static and dynamic beam problems by the p-version of the finite element method are investigated. Within a hierarchy of engineering beam idealizations, rigorous formulations of the strain and kinetic energies for straight and circular beam elements are presented. These formulations include rotating coordinate system effects and geometric nonlinearities to allow for the evaluation of vertical axis wind turbines, the motivating problem for this research. Hierarchic finite element spaces, based on extensions of the polynomial orders used to approximate the displacement variables, are constructed. The developed models are implemented into a general purpose computer program for evaluation. Quality control procedures are examined for a diverse set of sample problems. These procedures include estimating discretization errors in energy norm and natural frequencies, performing static and dynamic equilibrium checks, observing convergence for qualities of interest, and comparison with more exacting theories and experimental data. It is demonstrated that p-extensions produce exponential rates of convergence in the approximation of strain energy and natural frequencies for the class of problems investigated.

Watson, Robert A.↗
NASA NTRS2007

Comparing Consider-Covariance Analysis with Sigma-Point Consider Filter and Linear-Theory Consider Filter Formulations

Recent literature in applied estimation theory reflects growing interest in the sigma-point (also called unscented ) formulation for optimal sequential state estimation, often describing performance comparisons with extended Kalman filters as applied to specific dynamical problems [c.f. 1, 2, 3]. Favorable attributes of sigma-point filters are described as including a lower expected error for nonlinear even non-differentiable dynamical systems, and a straightforward formulation not requiring derivation or implementation of any partial derivative Jacobian matrices. These attributes are particularly attractive, e.g. in terms of enabling simplified code architecture and streamlined testing, in the formulation of estimators for nonlinear spaceflight mechanics systems, such as filter software onboard deep-space robotic spacecraft. As presented in [4], the Sigma-Point Consider Filter (SPCF) algorithm extends the sigma-point filter algorithm to the problem of consider covariance analysis. Considering parameters in a dynamical system, while estimating its state, provides an upper bound on the estimated state covariance, which is viewed as a conservative approach to designing estimators for problems of general guidance, navigation and control. This is because, whether a parameter in the system model is observable or not, error in the knowledge of the value of a non-estimated parameter will increase the actual uncertainty of the estimated state of the system beyond the level formally indicated by the covariance of an estimator that neglects errors or uncertainty in that parameter. The equations for SPCF covariance evolution are obtained in a fashion similar to the derivation approach taken with standard (i.e. linearized or extended) consider parameterized Kalman filters (c.f. [5]). While in [4] the SPCF and linear-theory consider filter (LTCF) were applied to an illustrative linear dynamics/linear measurement problem, in the present work examines the SPCF as applied to nonlinear sequential consider covariance analysis, i.e. in the presence of nonlinear dynamics and nonlinear measurements. A simple SPCF for orbit determination, exemplifying an algorithm hosted in the guidance, navigation and control (GN&C) computer processor of a hypothetical robotic spacecraft, was implemented, and compared with an identically-parameterized (standard) extended, consider-parameterized Kalman filter. The onboard filtering scenario examined is a hypothetical spacecraft orbit about a small natural body with imperfectly-known mass. The formulations, relative complexities, and performances of the filters are compared and discussed.

Lisano, Michael E.↗
NASA NTRS2002

Nonlinear Viscoelastic Mechanics of Cross-linked Rubbers

The paper develops a general theory for finite rubber viscoelasticity, and specifies it in the form, convenient for solving problems important for rubber, tire and space industries. Based on the quasi-linear approach of non-equilibrium thermodynamics, a general nonlinear theory has been developed for arbitrary nonisothermal deformations of viscoelastic solids. In this theory, the constitutive equations are presented as the sum of known equilibrium (rubber elastic) and non-equilibrium (liquid polymer viscoelastic) terms. These equations are then simplified using several modeling arguments. Stability constraints for the proposed constitutive equations are also discussed. It is shown that only strong ellipticity criteria are applicable for assessing stability of the equations governing viscoelastic solids.

Freed, Alan D.↗
NASA NTRS2003

Unified Models of Turbulence and Nonlinear Wave Evolution in the Extended Solar Corona and Solar Wind

The PI (Cranmer) and Co-I (A. van Ballegooijen) made significant progress toward the goal of building a "unified model" of the dominant physical processes responsible for the acceleration of the solar wind. The approach outlined in the original proposal comprised two complementary pieces: (1) to further investigate individual physical processes under realistic coronal and solar wind conditions, and (2) to extract the dominant physical effects from simulations and apply them to a one-dimensional and time-independent model of plasma heating and acceleration. The accomplishments in the report period are thus divided into these two categories: 1a. Focused Study of Kinetic MHD Turbulence. We have developed a model of magnetohydrodynamic (MHD) turbulence in the extended solar corona that contains the effects of collisionless dissipation and anisotropic particle heating. A turbulent cascade is one possible way of generating small-scale fluctuations (easy to dissipate/heat) from a pre-existing population of low-frequency Alfven waves (difficult to dissipate/heat). We modeled the cascade as a combination of advection and diffusion in wavenumber space. The dominant spectral transfer occurs in the direction perpendicular to the background magnetic field. As expected from earlier models, this leads to a highly anisotropic fluctuation spectrum with a rapidly decaying tail in the parallel wavenumber direction. The wave power that decays to high enough frequencies to become ion cyclotron resonant depends on the relative strengths of advection and diffusion in the cascade. For the most realistic values of these parameters, though, there is insufficient power to heat protons and heavy ions. The dominant oblique waves undergo Landau damping, which implies strong parallel electron heating. We thus investigated the nonlinear evolution of the electron velocity distributions (VDFs) into parallel beams and discrete phase-space holes (similar to those seen in the terrestrial magnetosphere) which are an alternate means of heating protons via stochastic interactions similar to particle-particle collisions. 1b. Focused Study of the Multi-Mode Detailed Balance Formalism. The PI began to explore the feasibility of using the "weak turbulence," or detailed-balance theory of Tsytovich, Melrose, and others to encompass the relevant physics of the solar wind. This study did not go far, however, because if the "strong" MHD turbulence discussed above is a dominant player in the wind's acceleration region, this formalism is inherently not applicable to the corona. We will continue to study the various published approaches to the weak turbulence formalism, especially with an eye on ways to parameterize nonlinear wave reflection rates. 2. Building the Unified Model Code Architecture. We have begun developing the computational model of a time-steady open flux tube in the extended corona. The model will be "unified" in the sense that it will include (simultaneously for the first time) as many of the various proposed physical processes as possible, all on equal footing. To retain this generality, we have formulated the problem in two interconnected parts: a completely kinetic model for the particles, using the Monte Carlo approach, and a finite-difference approach for the self-consistent fluctuation spectra. The two codes are run sequentially and iteratively until complete consistency is achieved. The current version of the Monte Carlo code incorporates gravity, the zero-current electric field, magnetic mirroring, and collisions. The fluctuation code incorporates WKJ3 wave action conservation and the cascade/dissipation processes discussed above. The codes are being run for various test problems with known solutions. Planned additions to the codes include prescriptions for nonlinear wave steepening, kinetic velocity-space diffusion, and multi-mode coupling (including reflection and refraction).

Cranmer, Steven R.↗
DOE OSTI2020

Data-Driven Operator Theoretic Methods for Global Phase Space Learning

In this work, we propose to apply the recently developed Koopman operator techniques to explore the global phase space of a nonlinear system from time-series data. In particular, we address the problem of identifying various invariant subsets of a phase space from the spectral properties of the associated Koopman operator constructed from time-series data. Moreover, in the case when the system evolution is known locally in various invariant subspaces, then a phase space stitching result is proposed that can be applied to identify a global Koopman operator. A biological system, bistable toggle switch and a second-order nonlinear system example are considered to illustrate the proposed results. The construction of this global Koopman operator is very helpful in experimental works as multiple experiments can't be run at the same time starting from several initial conditions.

Nandanoori, Sai Pushpak↗
NASA NTRS1988

Optimizing advanced propeller designs by simultaneously updating flow variables and design parameters

A scheme is developed for solving constrained optimization problems in which the objective function and the constraint function are dependent on the solution of the nonlinear flow equations. The scheme updates the design parameter iterative solutions and the flow variable iterative solutions simultaneously. It is applied to an advanced propeller design problem with the Euler equations used as the flow governing equations. The scheme's accuracy, efficiency and sensitivity to the computational parameters are tested.

Rizk, Magdi H.↗
NASA NTRS1988

Aerodynamic optimization by simultaneously updating flow variables and design parameters with application to advanced propeller designs

A scheme is developed for solving constrained optimization problems in which the objective function and the constraint function are dependent on the solution of the nonlinear flow equations. The scheme updates the design parameter iterative solutions and the flow variable iterative solutions simultaneously. It is applied to an advanced propeller design problem with the Euler equations used as the flow governing equations. The scheme's accuracy, efficiency and sensitivity to the computational parameters are tested.

Rizk, Magdi H.↗
DOE OSTI2021

Data-driven optimization of mixed-integer bi-level multi-follower integrated planning and scheduling problems under demand uncertainty

The coordination of interconnected elements across the different layers of the supply chain is essential for all industrial processes and the key to optimal decision-making. Yet, the modeling and optimization of such interdependent systems are still burdensome. Here we address the simultaneous modeling and optimization of medium-term planning and short-term scheduling problems under demand uncertainty using mixed-integer bi-level multi-follower programming and data-driven optimization. Bi-level multi-follower programs model the natural hierarchy between different layers of supply chain management holistically, while scenario analysis and data-driven optimization allow us to retrieve the guaranteed feasible solutions of the integrated formulation under various demand considerations. We address the data-driven optimization of this challenging class of problems using the DOMINO framework, which was initially developed to solve single-leader single-follower bi-level optimization problems to guaranteed feasibility. This framework is extended to solve single-leader multi-follower stochastic formulations and its performance is characterized by well-known single and multi-product process scheduling case studies. Through our data-driven algorithmic approach, we present guaranteed feasible solutions to linear and nonlinear mixed-integer bi-level formulations of simultaneous planning and scheduling problems and further characterize the effects of the scheduling level complexity on the solution performance, which spans over several hundred continuous and binary variables, and thousands of constraints.

42 ENGINEERING↗
NASA NTRS1992

On the physics of waves in the solar atmosphere: Wave heating and wind acceleration

In the area of solar physics, new calculations of the acoustic wave energy fluxes generated in the solar convective zone was performed. The original theory developed was corrected by including a new frequency factor describing temporal variations of the turbulent energy spectrum. We have modified the original Stein code by including this new frequency factor, and tested the code extensively. Another possible source of the mechanical energy generated in the solar convective zone is the excitation of magnetic flux tube waves which can carry energy along the tubes far away from the region. The problem as to how efficiently those waves are generated in the Sun was recently solved. The propagation of nonlinear magnetic tube waves in the solar atmosphere was calculated, and mode coupling, shock formation, and heating of the local medium was studied. The wave trapping problems and evaluation of critical frequencies for wave reflection in the solar atmosphere was studied. It was shown that the role played by Alfven waves in the wind accelerations and the coronal hole heating is dominant. Presently, we are performing calculations of wave energy fluxes generated in late-type dwarf stars and studying physical processes responsible for the heating of stellar chromospheres and coronae. In the area of physics of waves, a new analytical approach for studying linear Alfven waves in smoothly nonuniform media was recently developed. This approach is presently being extended to study the propagation of linear and nonlinear magnetohydrodynamic (MHD) waves in stratified, nonisothermal and solar atmosphere. The Lighthill theory of sound generation to nonisothermal media (with a special temperature distribution) was extended. Energy cascade by nonlinear MHD waves and possible chaos driven by these waves are presently considered.

Musielak, Z. E.↗