SEARCH · Search NASA

Results for “Nonlinear Solvers”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS1983

Postbuckling of long orthotropic plates in combined shear and compression

The nonlinear large-deflection partial differential equations of von karman for orthotropic plates loaded in combined shear and compression are converted into a set of first-order nonlinear ordinary differential equations by assuming trigonometric functions in one direction. These equations are solved numerically using a two point boundary problem solver which makes use of Newton's method. Results are obtained which determine the postbuckling behavior of rectangular plates with loading up to about three times the buckling load. Both isotropic and orthotropic composite plates are considered. Results show that orthotropic plates may behave quite differently than isotropic plates and that in-plane boundary conditions are important for plates loaded in shear.

Stein, M.↗

NASA NTRS1984

Crashworthiness of light aircraft fuselage structures: A numerical and experimental investigation

The dynamic behavior of aircraft fuselage structures subject to various impact conditions was investigated. An analytical model was developed based on a self-consistent finite element (CFE) formulation utilizing shell, curved beam, and stringer type elements. Equations of motion were formulated and linearized (i.e., for small displacements), although material nonlinearity was retained to treat local plastic deformation. The equations were solved using the implicit Newmark-Beta method with a frontal solver routine. Stiffened aluminum fuselage models were also tested in free flight using the UTIAS pendulum crash test facility. Data were obtained on dynamic strains, g-loads, and transient deformations (using high speed photography in the latter case) during the impact process. Correlations between tests and predicted results are presented, together with computer graphics, based on the CFE model. These results include level and oblique angle impacts as well as the free-flight crash test. Comparisons with a hybrid, lumped mass finite element computer model demonstrate that the CFE formulation provides the test overall agreement with impact test data for comparable computing costs.

Nanyaro, A. P.↗

NASA NTRS2019

Accuracy, Scalability, and Efficiency of Mixed-Element USM3D for Benchmark Three-Dimensional Flows

The unstructured, mixed-element, cell-centered, finite-volume flow solver USM3D is enhanced with new capabilities including parallelization, line generation for general unstructured grids, improved discretization scheme, and optimized iterative solver. The paper reports on the new developments to the flow solver and assesses the accuracy, scalability, and efficiency. The USM3D assessments are conducted using a baseline method and the recent hierarchical adaptive nonlinear iteration method framework. Two benchmark turbulent flows, namely, a subsonic separated flow around a three-dimensional hemisphere-cylinder configuration and a transonic flow around the ONERA M6 wing are considered.

Pandya, Mohagna J.↗

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 NTRS1993

An iterative multidisciplinary analysis for rotor blade shape determination

A CFD solver called ADPAC-APES is coupled with a NASTRAN structural analysis and a MARC thermal/heat transfer analysis to determine rotor blade shape. Nonlinear blade displacements due to centrifugal loads, aerodynamic pressures, and nonuniform temperature distribution are determined simultaneously. The effect of blade displacements on aerodynamic pressures and temperatures is then analyzed. These calculations are iterated till a steady state is reached across all the disciplines. This iterative procedure is applied to a ducted fan rotor blade and the manufactured shape is determined from a given operating shape. Effect of a part-span shroud on blade deflections is also analyzed.

Mahajan, Aparajit J.↗

NASA NTRS2019

Aero-Structural Modeling of the Truss-Braced Wing Aircraft Using Potential Method with Correction Methods for Transonic Viscous Flow and Wing-Strut Interference Aerodynamics

This paper describes an aero-structural modeling method for the Transonic Truss-Braced Wing (TTBW) aircraft using VSPAERO. A vortex-lattice model of the TTBW aircraft is developed, and a transonic and viscous flow correction method is implemented in the VSPAERO models to account for transonic and viscous flow effects. A correction method for the wing-strut interference aerodynamics is developed and applied to the VSPAERO solver. Also, a structural dynamic finite-element model of the TTBW aircraft is developed. This finite-element model includes the geometric nonlinear effect due to the tension in the struts which cause a deflection dependent nonlinear stiffness. The VSPAERO models are coupled to the finite-element model to provide a rapid capability for aero-structural modeling and flutter analysis. A flight-optimized jig twist model is being developed and will be applied for the purpose of generating a full flight dynamic model of the TTBW aircraft.

Fugate, Jason↗

NASA NTRS1989

Assessment of numerical methods for the solution of fluid dynamics equations for nonlinear resonance systems

The capability of accurate nonlinear flow analysis of resonance systems is essential in many problems, including combustion instability. Classical numerical schemes are either too diffusive or too dispersive especially for transient problems. In the last few years, significant progress has been made in the numerical methods for flows with shocks. The objective was to assess advanced shock capturing schemes on transient flows. Several numerical schemes were tested including TVD, MUSCL, ENO, FCT, and Riemann Solver Godunov type schemes. A systematic assessment was performed on scalar transport, Burgers' and gas dynamic problems. Several shock capturing schemes are compared on fast transient resonant pipe flow problems. A system of 1-D nonlinear hyperbolic gas dynamics equations is solved to predict propagation of finite amplitude waves, the wave steepening, formation, propagation, and reflection of shocks for several hundred wave cycles. It is shown that high accuracy schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time.

Przekwas, A. J.↗

NASA NTRS1999

Flutter and Forced Response Analyses of Cascades using a Two-Dimensional Linearized Euler Solver

Flutter and forced response analyses for a cascade of blades in subsonic and transonic flow is presented. The structural model for each blade is a typical section with bending and torsion degrees of freedom. The unsteady aerodynamic forces due to bending and torsion motions. and due to a vortical gust disturbance are obtained by solving unsteady linearized Euler equations. The unsteady linearized equations are obtained by linearizing the unsteady nonlinear equations about the steady flow. The predicted unsteady aerodynamic forces include the effect of steady aerodynamic loading due to airfoil shape, thickness and angle of attack. The aeroelastic equations are solved in the frequency domain by coupling the un- steady aerodynamic forces to the aeroelastic solver MISER. The present unsteady aerodynamic solver showed good correlation with published results for both flutter and forced response predictions. Further improvements are required to use the unsteady aerodynamic solver in a design cycle.

Reddy, T. S. R.↗

NASA NTRS1991

Parallelized solvers for heat conduction formulations

Based on multilevel partitioning, this paper develops a structural parallelizable solution methodology that enables a significant reduction in computational effort and memory requirements for very large scale linear and nonlinear steady and transient thermal (heat conduction) models. Due to the generality of the formulation of the scheme, both finite element and finite difference simulations can be treated. Diverse model topologies can thus be handled, including both simply and multiply connected (branched/perforated) geometries. To verify the methodology, analytical and numerical benchmark trends are verified in both sequential and parallel computer environments.

Padovan, Joe↗

NASA NTRS1996

Practical Aspects of Krylov Subspace Iterative Methods in CFD

Implementation issues associated with the application of Krylov subspace iterative methods, such as Newton-GMRES, are presented within the framework of practical computational fluid dynamic (CFD) applications. This paper categorizes, evaluates, and contrasts the major ingredients (function evaluations, matrix-vector products, and preconditioners) of Newton-GMRES Krylov subspace methods in terms of their effect on the local linear and global nonlinear convergence, memory requirements, and accuracy. The discussion focuses on Newton-GMRES in both a structured multi-zone incompressible Navier-Stokes solver and an unstructured mesh finite-volume Navier-Stokes solver. Approximate versus exact matrix-vector products, effective preconditioners, and other pertinent issues are addressed.

Pulliam, Thomas H.↗

NASA NTRS2014

Static Aeroelastic Analysis with an Inviscid Cartesian Method

An embedded-boundary, Cartesian-mesh flow solver is coupled with a three degree-of-freedom structural model to perform static, aeroelastic analysis of complex aircraft geometries. The approach solves a nonlinear, aerostructural system of equations using a loosely-coupled strategy. An open-source, 3-D discrete-geometry engine is utilized to deform a triangulated surface geometry according to the shape predicted by the structural model under the computed aerodynamic loads. The deformation scheme is capable of modeling large deflections and is applicable to the design of modern, very-flexible transport wings. The coupling interface is modular so that aerodynamic or structural analysis methods can be easily swapped or enhanced. After verifying the structural model with comparisons to Euler beam theory, two applications of the analysis method are presented as validation. The first is a relatively stiff, transport wing model which was a subject of a recent workshop on aeroelasticity. The second is a very flexible model recently tested in a low speed wind tunnel. Both cases show that the aeroelastic analysis method produces results in excellent agreement with experimental data.

Aeroelastic↗

NASA NTRS2015

Toward a Nonlinear Acoustic Analogy: Turbulence as a Source of Sound and Nonlinear Propagation

An acoustic analogy is proposed that directly includes nonlinear propagation effects. We examine the Lighthill acoustic analogy and replace the Green's function of the wave equation with numerical solutions of the generalized Burgers' equation. This is justified mathematically by using similar arguments that are the basis of the solution of the Lighthill acoustic analogy. This approach is superior to alternatives because propagation is accounted for directly from the source to the far-field observer instead of from an arbitrary intermediate point. Validation of a numerical solver for the generalized Burgers' equation is performed by comparing solutions with the Blackstock bridging function and measurement data. Most importantly, the mathematical relationship between the Navier- Stokes equations, the acoustic analogy that describes the source, and canonical nonlinear propagation equations is shown. Example predictions are presented for nonlinear propagation of jet mixing noise at the sideline angle

Miller, Steven A. E.↗

NASA NTRS2015

Toward a Nonlinear Acoustic Analogy: Turbulence as a Source of Sound and Nonlinear Propagation

An acoustic analogy is proposed that directly includes nonlinear propagation effects. We examine the Lighthill acoustic analogy and replace the Green's function of the wave equation with numerical solutions of the generalized Burgers' equation. This is justified mathematically by using similar arguments that are the basis of the solution of the Lighthill acoustic analogy. This approach is superior to alternatives because propagation is accounted for directly from the source to the far-field observer instead of from an arbitrary intermediate point. Validation of a numerical solver for the generalized Burgers' equation is performed by comparing solutions with the Blackstock bridging function and measurement data. Most importantly, the mathematical relationship between the Navier-Stokes equations, the acoustic analogy that describes the source, and canonical nonlinear propagation equations is shown. Example predictions are presented for nonlinear propagation of jet mixing noise at the sideline angle.

Miller, Steven A. E.↗

NASA NTRS1986

A full potential rotor analysis with wake influence using an inner-outer domain technique

A three-dimensional, quasi-steady, full potential flow solver was adapted to include realistic rotor wake influence for the aerodynamic analysis of helicopter rotors. The method uses an inner/outer domain technique to accommodate wake effects. Nonlinear flow is computed in the inner domain using a finite difference solution method. The wake is modeled using prescribed wake techniques to allow for the inclusion of realistic wake geometries. Portions of the wake passing inside the inner domain are treated using an embedded vortex technique. The procedure couples the wake influence with the inner domain solution in a consistent and efficient solution process. Correlation with measured lifting transonic data in hover and forward flight is shown which demonstrates the merits of the approach.

Egolf, T. A.↗

NASA NTRS2021

Efficiency of Mixed-Element USM3D for Benchmark Three-Dimensional Flows

The unstructured, mixed-element, cell-centered, finite-volume flow solver USM3D is enhanced with new capabilities such as line generation for general unstructured grids and optimization of the iterative solver. The efficiency of the new developments is assessed in a massively parallel environment. The assessments are conducted using a baseline method and the recent hierarchical adaptive nonlinear iteration method framework. Three benchmark turbulent flows, namely, a subsonic separated flow around a three-dimensional hemisphere cylinder configuration, a transonic flow around the ONERA M6 wing, and a supersonic flow through a square duct are considered.

Mohagna J Pandya↗

NASA NTRS1986

Predicting the Dynamic Behavior of High-Speed Roller Bearings

Computer program developed for solving nonlinear differential equations of motion of high-speed, lightly loaded, cylindrical roller bearing. Program organized into three major parts: Main program, controls input and output; differential equation solver; and subroutines, used to evaluate vector function. Program has capability of performing either twodimensional or three-dimensional simulation.

Conry, T. F.↗

NASA NTRS2019

Predicting Transonic Flutter Using Nonlinear Computational Simulations

This paper builds on the computational aeroelastic results published previously and generated in support of the second Aeroelastic Prediction Workshop for the NASA Benchmark Supercritical Wing configuration. The computational results are obtained using FUN3D, an unstructured flow solver developed at the NASA Langley Research Center. Flutter analyses are performed using various aerodynamic models including linear doublet lattice, Euler solutions, Reynolds averaged Navier-Stokes and Delayed Detached Eddy Simulations. The analyses are performed across the transonic Mach range for a range in angle of attack from 0° to 5°, with focus on identifying the transonic dips in the flutter onset boundaries.

Heeg, Jennifer↗

NASA NTRS2000

Applications of the Space-Time Conservation Element and Solution Element (CE/SE) Method to Computational Aeroacoustic Benchmark Problems

The Internal Propagation problems, Fan Noise problem, and Turbomachinery Noise problems are solved using the space-time conservation element and solution element (CE/SE) method. The problems in internal propagation problems address the propagation of sound waves through a nozzle. Both the nonlinear and linear quasi 1D Euler equations are solved. Numerical solutions are presented and compared with the analytical solution. The fan noise problem concerns the effect of the sweep angle on the acoustic field generated by the interaction of a convected gust with a cascade of 3D flat plates. A parallel version of the 3D CE/SE Euler solver is developed and employed to obtain numerical solutions for a family of swept flat plates. Numerical solutions for sweep angles of 0, 5, 10, and 15 deg are presented. The turbomachinery problems describe the interaction of a 2D vortical gust with a cascade of flat-plate airfoils with/without a downstream moving grid. The 2D nonlinear Euler Equations are solved and the converged numerical solutions are presented and compared with the corresponding analytical solution. All the comparisons demonstrate that the CE/SE method is capable of solving aeroacoustic problems with/without shock waves in a simple and efficient manner. Furthermore, the simple non-reflecting boundary condition used in the CE/SE method which is not based on the characteristic theory works very well in 1D, 2D and 3D problems.

Wang, Xiao-Yen↗