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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Hamilton's Equations with Euler Parameters for Rigid Body Dynamics Modeling

A combination of Euler parameter kinematics and Hamiltonian mechanics provides a rigid body dynamics model well suited for use in strongly nonlinear problems involving arbitrarily large rotations. The model is unconstrained, free of singularities, includes a general potential energy function and a minimum set of momentum variables, and takes an explicit state space form convenient for numerical implementation. The general formulation may be specialized to address particular applications, as illustrated in several three dimensional example problems.

Shivarama, Ravishankar↗
NASA NTRS2006

An Ensemble-Based Smoother with Retrospectively Updated Weights for Highly Nonlinear Systems

Monte Carlo computational methods have been introduced into data assimilation for nonlinear systems in order to alleviate the computational burden of updating and propagating the full probability distribution. By propagating an ensemble of representative states, algorithms like the ensemble Kalman filter (EnKF) and the resampled particle filter (RPF) rely on the existing modeling infrastructure to approximate the distribution based on the evolution of this ensemble. This work presents an ensemble-based smoother that is applicable to the Monte Carlo filtering schemes like EnKF and RPF. At the minor cost of retrospectively updating a set of weights for ensemble members, this smoother has demonstrated superior capabilities in state tracking for two highly nonlinear problems: the double-well potential and trivariate Lorenz systems. The algorithm does not require retrospective adaptation of the ensemble members themselves, and it is thus suited to a streaming operational mode. The accuracy of the proposed backward-update scheme in estimating non-Gaussian distributions is evaluated by comparison to the more accurate estimates provided by a Markov chain Monte Carlo algorithm.

Monte Carlo↗
NASA NTRS2007

Simulations of Boundary-Layer Transition

For incompressible benchmark flows, we have demonstrated the capability of the parabolized stability equations (PSE) to simulate the transition process in excellent agreement with microscopic experiments and direct Navier-Stokes simulations at modest computational cost. Encouraged by these results, we have developed the PSE methodology of three-dimensional boundary-layers in general curvilinear coordinates for the range from low to hypersonic speeds, and for both linear and nonlinear problems. For given initial and boundary conditions, the approach permits simulations from receptivity through linear and secondary instabilities into the late stages of transition where significant changes in skin friction and heat transfer coefficients occur. We have performed transition simulations for a variety of two- and three-dimensional similarity solutions and for realistic flows over swept wings at subsonic and supersonic speeds, the pressure ans suction side of turbine blades at low and medium turbulence levels, and over a blunt cone at Mach number Ma = 8. We present selected results for different transition mechanisms with emphasis on the late stage of transition and the evolution of wall-shear stress and heat transfer.

Herbert, Thorwald↗
NASA NTRS2010

Comparison of Sigma-Point and Extended Kalman Filters on a Realistic Orbit Determination Scenario

Sigma-point filters have received a lot of attention in recent years as a better alternative to extended Kalman filters for highly nonlinear problems. In this paper, we compare the performance of the additive divided difference sigma-point filter to the extended Kalman filter when applied to orbit determination of a realistic operational scenario based on the Interstellar Boundary Explorer mission. For the scenario studied, both filters provided equivalent results. The performance of each is discussed in detail.

Gaebler, John↗
NASA NTRS2015

Beamed-Energy Propulsion (BEP): Considerations for Beaming High Energy-Density Electromagnetic Waves Through the Atmosphere

A study to determine the feasibility of employing beamed electromagnetic energy for vehicle propulsion within and outside the Earth's atmosphere was co-funded by NASA and the Defense Advanced Research Projects Agency that began in June 2010 and culminated in a Summary Presentation in April 2011. A detailed report entitled "Beamed-Energy Propulsion (BEP) Study" appeared in February 2012 as NASA/TM-2012-217014. Of the very many nuances of this subject that were addressed in this report, the effects of transferring the required high energy-density electromagnetic fields through the atmosphere were discussed. However, due to the limitations of the length of the report, only a summary of the results of the detailed analyses were able to be included. It is the intent of the present work to make available the complete analytical modeling work that was done for the BEP project with regard to electromagnetic wave propagation issues. In particular, the present technical memorandum contains two documents that were prepared in 2011. The first one, entitled "Effects of Beaming Energy Through the Atmosphere" contains an overview of the analysis of the nonlinear problem inherent with the transfer of large amounts of energy through the atmosphere that gives rise to thermally-induced changes in the refractive index; application is then made to specific beamed propulsion scenarios. A brief portion of this report appeared as Appendix G of the 2012 Technical Memorandum. The second report, entitled "An Analytical Assessment of the Thermal Blooming Effects on the Propagation of Optical and Millimeter- Wave Focused Beam Waves For Power Beaming Applications" was written in October 2010 (not previously published), provides a more detailed treatment of the propagation problem and its effect on the overall characteristics of the beam such as its deflection as well as its radius. Comparisons are then made for power beaming using the disparate electromagnetic wavelengths of 1.06 microns and 2.0 millimeters..

Beamed Energy↗
NASA NTRS2015

De-Aliasing Through Over-Integration Applied to the Flux Reconstruction and Discontinuous Galerkin Methods

High-order methods are quickly becoming popular for turbulent flows as the amount of computer processing power increases. The flux reconstruction (FR) method presents a unifying framework for a wide class of high-order methods including discontinuous Galerkin (DG), Spectral Difference (SD), and Spectral Volume (SV). It offers a simple, efficient, and easy way to implement nodal-based methods that are derived via the differential form of the governing equations. Whereas high-order methods have enjoyed recent success, they have been known to introduce numerical instabilities due to polynomial aliasing when applied to under-resolved nonlinear problems. Aliasing errors have been extensively studied in reference to DG methods; however, their study regarding FR methods has mostly been limited to the selection of the nodal points used within each cell. Here, we extend some of the de-aliasing techniques used for DG methods, primarily over-integration, to the FR framework. Our results show that over-integration does remove aliasing errors but may not remove all instabilities caused by insufficient resolution (for FR as well as DG).

High-Order↗
NASA NTRS2019

eddy Users Manual

eddy is a collection of tools - nonlinear solvers, meshing, post-processing, visualization, optimization, etc. - for performing scale-resolving simulations of multi-physics applications. The framework is designed to enable advanced R&D on a variety of topics by leveraging a mature capability for scale resolving simulations, and simultaneously be an appropriate tool for application analysis and support. Currently, eddy is at a relatively low technical readiness level (TRL), and users and developers should maintain appropriate expectations. The technical details behind eddy are outlined in several publications which can be consulted for more information [1–10]. The solvers are built around an unstructured high-order capability, and heavily utilize the tensor product sum-factorization approach for efficiency. The unsteady formulation utilizes a fully implicit space-time approach with a matrix-free Newton- Krylov method. A primitive steady-state solver is available for testing purposes, but is not expected to converge for all but simple verification cases. The Navier-Stokes fluid solvers do not support either RANS or hybrid-RANS capability, only LES and wall-modeled LES approaches. All of the solvers within eddy support three modes of operation: a primal solve of the full nonlinear problem, and two linearization approaches of the primal solve - the ad joint and the tangent solution. Details on how to select and use these three modes are outlined in Sec. 3.

Murman, Scott M.↗
NASA NTRS2020

Multigrid Preconditioning for a Space-Time Spectral-Element Discontinuous-Galerkin Solver

In this work we examine a multigrid preconditioning approach in the context of a high- order tensor-product discontinuous-Galerkin spectral-element solver. We couple multigrid ideas together with memory lean and efficient tensor-product preconditioned matrix-free smoothers. Block ILU(0)-preconditioned GMRES smoothers are employed on the coarsest spaces. The performance is evaluated on nonlinear problems arising from unsteady scale- resolving solutions of the Navier-Stokes equations: separated low-Mach unsteady ow over an airfoil from laminar to turbulent ow. A reduction in the number of ne space iterations is observed, which proves the efficiency of the approach in terms of preconditioning the linear systems, however this gain was not reflected in the CPU time. Finally, the preconditioner is successfully applied to problems characterized by stiff source terms such as the set of RANS equations, where the simple tensor product preconditioner fails. Theoretical justification about the findings is reported and future work is outlined.

Franciolini, Matteo↗
NASA NTRSDate not supplied

Scalability of Cohesive Fatigue Analyses Using Explicit Solvers

A cohesive fatigue law has been integrated into a constitutive material model compatible with an explicit finite element solver. The cohesive fatigue model response is based on engineering approximations of the endurance limit and the Goodman diagram. This approach can predict stress-life diagrams for crack initiation, the Paris law regime, and transient effects of crack initiation and stable tearing. Simplified cyclic loading is utilized so that the applied load(or displacement) corresponds to the peak load of a fatigue cycle. Loads are held constant during fatigue while damage develops with increasing solution increments. An automatically-calculated ratio of fatigue cycles per solution increment controls the rate of damage growth, ensuring that damage growth is modeled with a sufficient minimum number of increments and damage growth advances to a minimum desired extent within the explicit analysis step time. The compatibility with an explicit finite element solver enables the analysis of structures that are computationally intractable for implicit finite element solvers. Scalability studies are conducted for geometrically nonlinear problems involving fiber-reinforced composite structures that exhibit fatigue damage growth of interacting matrix cracks and delaminations

Frank A Leone↗
NASA NTRS1984

Energy conservation in the transient response of nonlinear beam vibration problems subjected to pulse loading - A numerical approach

The nonlinear vibration response of a double cantilevered beam subjected to pulse loading over a central sector is studied. The initial response is generated in detail to ascertain the energetics of the response. The total energy is used as a gauge of the stability and accuracy of the solution. It is shown that to obtain accurate and stable initial solutions an extremely high spatial and time resolution is required. This requirement was only evident through an examination of the energy of the system. It is proposed, therefore, to use the total energy of the system as a necessary stability and accuracy criterion for the nonlinear response of conservative systems. The results also demonstrate that even for moderate nonlinearities, the effects of membrane forces have a significant influence on the system. It is also shown that while the fundamental response is contained in a first mode envelope, the fluctuations caused by the higher order modes must be resolved.

Moyer, E. T., Jr.↗
NASA NTRS1988

Galerkin approximation for inverse problems for nonautonomous nonlinear distributed systems

An abstract framework and convergence theory is developed for Galerkin approximation for inverse problems involving the identification of nonautonomous nonlinear distributed parameter systems. A set of relatively easily verified conditions is provided which are sufficient to guarantee the existence of optimal solutions and their approximation by a sequence of solutions to a sequence of approximating finite dimensional identification problems. The approach is based on the theory of monotone operators in Banach spaces and is applicable to a reasonably broad class of nonlinear distributed systems. Operator theoretic and variational techniques are used to establish a fundamental convergence result. An example involving evolution systems with dynamics described by nonstationary quasilinear elliptic operators along with some applications are presented and discussed.

Banks, H. T.↗
NASA NTRS1992

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

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

Hanson, R. J.↗
NASA NTRS1993

A nonlinear control design for energy sink simulation in the Euler-Poinsot problem

A nonlinear control design is presented for the purpose of quantitatively simulating the effects of internal damping mechanisms modeled as energy sinks on the attitude dynamics of rigid body spacecraft. Damping is important because it is often the driving mechanism behind passive attitude acquisition maneuvers. Introduction of the controller into the Euler attitude equations of motion allows for the explicit representation of damping without the introduction of additional degrees of freedom required for a physical damping mechanism. This result is significant because perturbation techniques which rely on the closed form solution of the unperturbed problem can then be used to analyze the effects of perturbations upon a damped system. The controller is designed to dissipate kinetic energy while maintaining the magnitude of the angular momentum vector. Control torques are nonlinear functions of the angular momentum components expressed in a body-fixed frame. A numerical simulation of an actual damping mechanism during a decay from minor axis spin into a flat spin is presented showing that the nonlinear controller gives a good qualitative representation and, in many instances, a good quantitative approximation of the attitude motion of a representative spacecraft containing a damping mechanism.

Kammer, Daniel C.↗