Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonlinear approximation”

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 NTRS1981

On eigensolutions for discontinuous liners in a duct containing uniform mean flow

Sound attenuation in a rectangular acoustically lined duct containing uniform mean flow is analytically investigated using the generalized Wiener-Hopf technique. Uniqueness of the solution is enforced for lined sections of the finite axial extent by imposing edge conditions at the liner interface. Possible edge conditions are considered, including the Kutta condition, and the causal solution corresponding to edge conditions is considered the best choice. Solution methods such as the mode matching and singularity methods imply differing edge conditions, and results show that power attenuation is insensitive to the imposed edge conditions, although significant differences are observed for the reflection coefficient. The amplitude of the exponentially increasing instability mode in the lined section must be set to zero as a first approximation to the nonlinear situation, and results indicate that measurements of the reflection factor can be used to make a more definite decision about physically appropriate edge conditions.

Koch, W.↗
NASA NTRS1982

Numerical analysis of confined turbulent flow

The considered investigation is concerned with the development of an efficient computational method for obtaining a physical understanding of an internal turbulent field. The employed approach makes use of a 'two equation' type model for the turbulence to obtain the numerical solution of a two-dimensional confined turbulent flow. The mean flow governing equations are considered along with the governing equation of the mean temperature and concentrations, and the boundary conditions. The numerical procedure for solving the turbulent flow is discussed, taking into account an approximation to the nonlinear terms, and the inner and outer coupling. Attention is given to a stability convergence analysis, the stability characteristics, and computational examples.

Lin, A.↗
NASA NTRS1982

Unsymmetrically laminated composites

Classical lamination theory predicts that the room-temperature shapes of all elevated-temperature cure, unsymmetrically laminated composites are saddles. However, experimental observation indicates that the shapes are often cylindrical. In addition, a second cylindrical can sometimes be obtained from the first by a simple snap-through action. A geometrically nonlinear extension to classical lamination theory is used to explain this behavior. Approximate solutions to the nonlinear extension are obtained by using a Rayleigh-Ritz minimization of the laminate's total potential energy. A stability analysis explains the dual cylindrical shapes.

Hyer, M. W.↗
NASA NTRS1983

On the stellar rotation-activity connection

The relationship between rotation rates and surface activity in late-type dwarf stars is explored in a survey of recent theoretical and observational studies. Current theoretical models of stellar-magnetic-field production and coronal activity are examined, including linear kinematic dynamo theory, nonlinear dynamos using approximations, and full numerical simulations of the MHD equations; and some typical results are presented graphically. The limitations of the modeling procedures and the constraints imposed by the physics are indicated. The statistical techniques used in establishing correlations between various observational parameters are analyzed critically, and the methods developed for quasar luminosity functions by Avni et al. (1980) are used to evaluate the effects of upper detection bounds, incomplete samples, and missing data for the case of rotation and X-ray flux data.

Rosner, R.↗
NASA NTRS1986

Computation of turbulent pipe and duct flow using third order upwind scheme

The fully developed turbulence in a circular pipe and in a square duct is simulated directly without using turbulence models in the Navier-Stokes equations. The utilized method employs a third-order upwind scheme for the approximation to the nonlinear term and the second-order Adams-Bashforth method for the time derivative in the Navier-Stokes equation. The computational results appear to capture the large-scale turbulent structures at least qualitatively. The significance of the artificial viscosity inherent in the present scheme is discussed.

Kawamura, T.↗
NASA NTRS1989

Iterative procedures for space shuttle main engine performance models

Performance models of the Space Shuttle Main Engine (SSME) contain iterative strategies for determining approximate solutions to nonlinear equations reflecting fundamental mass, energy, and pressure balances within engine flow systems. Both univariate and multivariate Newton-Raphson algorithms are employed in the current version of the engine Test Information Program (TIP). Computational efficiency and reliability of these procedures is examined. A modified trust region form of the multivariate Newton-Raphson method is implemented and shown to be superior for off nominal engine performance predictions. A heuristic form of Broyden's Rank One method is also tested and favorable results based on this algorithm are presented.

Santi, L. Michael↗
NASA NTRS1992

Investigation of forced unsteady separated flows using velocity-vorticity form of Navier-Stokes equations

The phenomenon of forced unsteady separation and eruption of boundary-layer vorticity is a highly-complex, high-Reynolds number flow phenomenon, which abruptly leads to the formation of a dynamic stall vortex as demonstrated earlier by the authors for a NACA 0015 airfoil undergoing constant rate pitch-up motion. This, as well as the results of other researchers, have convincingly demonstrated a complex vortical structure within the state of unsteady separation prior to the evolution of dynamic stall. This phenomenon of vortex eruption, although observed in studying dynamic stall phenomena, is also associated with transition from laminar to turbulence flow and its generic nature has been stressed by many researchers including the present investigators. An unsteady Navier-Stokes (NS) analysis is developed for arbitrarily maneuvering bodies using velocity-vorticity variables; this formulation is nearly form-invariant under a generalized non-inertial coordinate transformation. A fully-implicit uniformly second-order accurate method is used, with the nonlinear convective terms approximated using a biased third-order upwind differencing scheme to be able to simulate higher-Re flows. No explicit artificial dissipation is added. The numerical method is fully vectorized and currently achieves a computational index of 7 micro-seconds per time step per mesh point, using a single processor on a CRAY Y-MP. The simulation results show that the energetic free shear from the leading edge is responsible for the wall viscous layer to abruptly erupt near the center of the counterclockwise rotating eddy in the unsteady boundary layer. Primary, secondary, tertiary and quaternary vortices have been observed before the dynamic stall vortex evolves and gathers its maximum strength. This study will discuss the simulation results of Reynolds number up to Re = 45,000 and will also discuss the efforts of initial acceleration in a specific maneuver, on the evolution of the stall vortex.

Ghia, K. N.↗
NASA NTRS1996

Analysis and Tests of Reinforced Carbon-Epoxy/Foam-Core Sandwich Panels with Cutouts

The results of a study of a low-cost structurally efficient minimum-gage shear-panel design that can be used in light helicopters are presented. The shear-panel design is based on an integrally stiffened syntactic-foam stabilized-skin with an all-bias-ply tape construction for stabilized-skin concept with an all-bias-ply tape construction for the skins. This sandwich concept is an economical way to increase the panel bending stiffness weight penalty. The panels considered in the study were designed to be buckling resistant up to 100 lbs/in. of shear load and to have an ultimate strength of 300 lbs/in. The panel concept uses unidirectional carbon-epoxy tape on a syntactic adhesive as a stiffener that is co-cured with the skin and is an effective concept for improving panel buckling strength. The panel concept also uses pultruded carbon-epoxy rods embedded in a syntactic adhesive and over-wrapped with a bias-ply carbon-epoxy tape to form a reinforcing beam which is an effective method for redistributing load around rectangular cutout. The buckling strength of the reinforced panels is 83 to 90 percent of the predicted buckling strength based on a linear buckling analysis. The maximum experimental deflection exceeds the maximum deflection predicted by a nonlinear analysis by approximately one panel thickness. The failure strength of the reinforced panels was two and a half to seven times of the buckling strength. This efficient shear-panel design concept exceeds the required ultimate strength requirement of 300 lbs/in by more than 100 percent.

Baker, Donald J.↗
NASA NTRS2007

NASA Langley's Approach to the Sandia's Structural Dynamics Challenge Problem

The objective of this challenge is to develop a data-based probabilistic model of uncertainty to predict the behavior of subsystems (payloads) by themselves and while coupled to a primary (target) system. Although this type of analysis is routinely performed and representative of issues faced in real-world system design and integration, there are still several key technical challenges that must be addressed when analyzing uncertain interconnected systems. For example, one key technical challenge is related to the fact that there is limited data on target configurations. Moreover, it is typical to have multiple data sets from experiments conducted at the subsystem level, but often samples sizes are not sufficient to compute high confidence statistics. In this challenge problem additional constraints are placed as ground rules for the participants. One such rule is that mathematical models of the subsystem are limited to linear approximations of the nonlinear physics of the problem at hand. Also, participants are constrained to use these models and the multiple data sets to make predictions about the target system response under completely different input conditions. Our approach involved initially the screening of several different methods. Three of the ones considered are presented herein. The first one is based on the transformation of the modal data to an orthogonal space where the mean and covariance of the data are matched by the model. The other two approaches worked solutions in physical space where the uncertain parameter set is made of masses, stiffnesses and damping coefficients; one matches confidence intervals of low order moments of the statistics via optimization while the second one uses a Kernel density estimation approach. The paper will touch on all the approaches, lessons learned, validation 1 metrics and their comparison, data quantity restriction, and assumptions/limitations of each approach. Keywords: Probabilistic modeling, model validation, uncertainty quantification, kernel density

Horta, Lucas G.↗
DOE OSTI2021

Design of double- and multi-bend achromat lattices with large dynamic aperture and approximate invariants

A numerical method to design nonlinear double- and multi-bend achromat (DBA and MBA) lattices with approximate invariants of motion is investigated. The search for such nonlinear lattices is motivated by Fermilab’s Integrable Optics Test Accelerator (IOTA), whose design is based on an integrable Hamiltonian system with two invariants of motion. While it may not be possible to design an achromatic lattice for a dedicated synchrotron light source storage ring with one or more exact invariants of motion, it is possible to tune the sextupoles and octupoles in existing DBA and MBA lattices to produce approximate invariants. In our procedure, the lattice is tuned while minimizing the turn-by-turn fluctuations of the Courant-Snyder actions J x and J y at several distinct amplitudes, while simultaneously minimizing diffusion of the on-energy betatron tunes. The resulting lattices share some important features with integrable ones, such as a large dynamic aperture, trajectories confined to invariant tori, robustness to resonances and errors, and a large amplitude-dependent tune-spread. Compared to the nominal NSLS-II lattice, the single- and multi-bunch instability thresholds are increased and the bunch-by-bunch feedback gain can be reduced.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Implicit shock tracking for unsteady flows by the method of lines

Here, a recently developed high-order implicit shock tracking (HOIST) framework for resolving discontinuous solutions of inviscid, steady conservation laws is extended to the unsteady case. Central to the framework is an optimization problem which simultaneously computes a discontinuity-aligned mesh and the corresponding high-order approximation to the flow, which provides nonlinear stabilization and a high-order approximation to the solution. This work extends the implicit shock tracking framework to the case of unsteady conservation laws using a method of lines discretization via a diagonally implicit Runge-Kutta method by “solving a steady problem at each timestep”. We formulate and solve an optimization problem that produces a feature-aligned mesh and solution at each Runge-Kutta stage of each timestep, and advance this solution in time by standard Runge-Kutta update formulas. A Rankine-Hugoniot based prediction of the shock location together with a high-order, untangling mesh smoothing procedure provides a high-quality initial guess for the optimization problem at each time, which results in rapid convergence of the sequential quadratic programing (SQP) optimization solver. This method is shown to deliver highly accurate solutions on coarse, high-order discretizations without nonlinear stabilization and recover the design accuracy of the Runge-Kutta scheme. We demonstrate this framework on a series of inviscid, unsteady conservation laws in both one- and two- dimensions. We also verify that our method is able to recover the design order of accuracy of our time integrator in the presence of a strong discontinuity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1978

An advanced structural analysis/synthesis capability - ACCESS 2

An advanced automated design procedure for minimum weight design of structures (ACCESS 2) is reported. Design variable linking, constraint deletion, and explicit constraint approximation are used to effectively combine finite element and nonlinear mathematical programming techniques. The approximation concepts approach to structural synthesis is extended to problems involving fiber composite structure, thermal effects and natural frequency constraints in addition to the usual static stress and displacement limitations. Sample results illustrating these new features are given.

Schmit, L. A.↗
NASA NTRS1986

Solving Nonlinear Coupled Differential Equations

Harmonic balance method developed to obtain approximate steady-state solutions for nonlinear coupled ordinary differential equations. Method usable with transfer matrices commonly used to analyze shaft systems. Solution to nonlinear equation, with periodic forcing function represented as sum of series similar to Fourier series but with form of terms suggested by equation itself.

Mitchell, L.↗
NASA NTRSDate not supplied

Bayesian Symbolic Regression: Addressing Challenges in Estimating Fractional Bayes Factors and Application to Fatigue Crack Growth Modeling

This research pioneers advancements in computational mechanics by integrating Bayesian-based uncertainty quantification into symbolic regression, specifically focusing on the critical task of accurately estimating the fractional Bayes factor for selecting arbitrary equations. In our exploration, we rigorously study two prominent methods—sequential Monte Carlo and the Laplace approximation—employed for computing the fractional Bayes factor. Our findings underscore the limitations of the Laplace approximation, revealing its diminished accuracy in nonlinear and multimodal scenarios. Specifically, the Laplace approximation is shown to underpredict fractional Bayes factor on a wide set of equations associated with a symbolic regression benchmark. This comparative analysis sheds light on the nuanced performance of these techniques, guiding researchers toward more informed choices in uncertainty quantification within symbolic regression. Furthermore, we showcase the practical utility of these enhanced symbolic regression tools through their application to a real-world problem in fatigue crack growth modeling, emphasizing their efficacy in capturing the complexities of mechanical systems.

Geoffrey Bomarito↗
DOE OSTI2023

Deep Koopman learning of nonlinear time-varying systems

Here this paper presents a data-driven approach to approximate the dynamics of a nonlinear time-varying system (NTVS) by a linear time-varying system (LTVS), which results from the Koopman operator and deep neural networks. Analysis of the approximation error between states of the NTVS and the resulting LTVS is presented. Simulations on a representative NTVS show that the proposed method achieves small approximation errors, even when the system changes rapidly. Furthermore, simulations in an example of quadcopters demonstrate the computational efficiency of the proposed approach.

97 MATHEMATICS AND COMPUTING↗