Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonlinear solution”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Computational mechanics analysis tools for parallel-vector supercomputers

Computational algorithms for structural analysis on parallel-vector supercomputers are reviewed. These parallel algorithms, developed by the authors, are for the assembly of structural equations, 'out-of-core' strategies for linear equation solution, massively distributed-memory equation solution, unsymmetric equation solution, general eigensolution, geometrically nonlinear finite element analysis, design sensitivity analysis for structural dynamics, optimization search analysis and domain decomposition. The source code for many of these algorithms is available.

Storaasli, Olaf O.↗
NASA NTRS1993

The barotropic normal modes in certain shear flows and the traveling waves in the atmosphere

It is shown analytically and numerically that in certain shear flows the linearized nondivergent barotropic vorticity equation has a limited number of neutral normal modes. The latitudinal structures of these shear flows can be expressed as polynomials of the sine of latitude. The first few such shear flows resemble the gross features of the zonal winds in the atmosphere of the earth at different times and altitudes. The spatial structures of the neutral normal modes in these shear flows are spherical harmonics, and, as a consequence, these modes are also the exact solutions of the fully nonlinear equation because the nonlinear interaction term vanishes identically. The spatial structures of the observed 5-, 4-, 2-, and 16-day free traveling waves in the atmosphere are often identified with the spherical harmonics with indices of (m, n) = ( 1, 2), (2, 3), (3, 3), and ( 1, 4), which are known previously as the neutral normal modes of the nondivergent barotropic vorticity equation in a motionless background state. Our results could explain why these free traveling waves can survive the shearing effects of zonal flows that are far different from rest because these spherical harmonics are also normal modes in certain shear flows that resemble the observations of the atmosphere.

Chen, Ping↗
NASA NTRS1988

Parameter identification of material constants in a composite shell structure

One of the basic requirements in engineering analysis is the development of a mathematical model describing the system. Frequently comparisons with test data are used as a measurement of the adequacy of the model. An attempt is typically made to update or improve the model to provide a test verified analysis tool. System identification provides a systematic procedure for accomplishing this task. The terms system identification, parameter estimation, and model correlation all refer to techniques that use test information to update or verify mathematical models. The goal of system identification is to improve the correlation of model predictions with measured test data, and produce accurate, predictive models. For nonmetallic structures the modeling task is often difficult due to uncertainties in the elastic constants. A finite element model of the shell was created, which included uncertain orthotropic elastic constants. A modal survey test was then performed on the shell. The resulting modal data, along with the finite element model of the shell, were used in a Bayes estimation algorithm. This permitted the use of covariance matrices to weight the confidence in the initial parameter values as well as confidence in the measured test data. The estimation procedure also employed the concept of successive linearization to obtain an approximate solution to the original nonlinear estimation problem.

Martinez, David R.↗
NASA NTRS1993

Numerical solution of flame sheet problems with and without multigrid methods

Flame sheet problems are on the natural route to the numerical solution of multidimensional flames, which, in turn, are important in many engineering applications. In order to model the structure of flames more accurately, we use the vorticity-velocity formulation of the fluid flow equations, as opposed to the streamfunction-vorticity approach. The numerical solution of the resulting nonlinear coupled elliptic partial differential equations involves a pseudo transient process and a steady state Newton iteration. Rather than working with dimensionless variables, we introduce scale factors that can yield significant savings in the execution time. In this context, we also investigate the applicability and performance of several multigrid methods, focusing on nonlinear damped Newton multigrid, using either one way or correction schemes.

Douglas, Craig C.↗
NASA NTRS1993

Direct time integration of Maxwell's equations in two-dimensional dielectric waveguides for propagation and scattering of femtosecond electromagnetic solitons

Solution of 2D vector nonlinear Maxwell's equations is described for material media with linear and nonlinear instantaneous and Lorentz dispersion effects in the electric polarization. Finite difference time domain method is used in the analysis. The method discussed here achieves robustness by enforcing vector-field boundary conditions at all interfaces of dissimilar media.

VECTOR FIELDS↗
NASA NTRS2002

Optimal Design of Functionally Graded Metallic Foam Insulations

The focus of our work has been on developing an insight into the physics that govern the optimum design of thermal insulation for use in thermal protection systems of launch vehicle. Of particular interest was to obtain optimality criteria for designing foam insulations that have density (or porosity) distributions through the thickness for optimum thermal performance. We investigate the optimum design of functionally graded thermal insulation for steady state heat transfer through the foam. We showed that the heat transfer in the foam has competing modes, of radiation and conduction. The problem assumed a fixed inside temperature of 400 K and varied the aerodynamic surface heating on the outside surface from 0.2 to 1.0 MW/sq m. The thermal insulation develops a high temperature gradient through the thickness. Investigation of the model developed for heat conduction in foams showed that at high temperatures (as on outside wall) intracellular radiation dominates the heat transfer in the foam. Minimizing radiation requires reducing the pore size, which increases the density of the foam. At low temperatures (as on the inside wall), intracellular conduction (of the metal and air) dominates the heat transfer. Minimizing conduction requires increasing the pore size. This indicated that for every temperature there was an optimum value of density that minimized the heat transfer coefficient. Two optimization studies were performed. One was to minimize the heat transmitted though a fixed thickness insulation by varying density profiles. The second was to obtain the minimum mass insulation for specified thickness. Analytical optimality criteria were derived for the cases considered. The optimality condition for minimum heat transfer required that at each temperature we find the density that minimizes the heat transfer coefficient. Once a relationship between the optimum heat transfer coefficient and the temperature was found, the design problem reduced to the solution of a simple nonlinear differential equation. Preliminary results of this work were presented at the American Society of Composites meeting, and the final version was submitted for publication in the AIAA Journal. In addition to minimizing the transmitted heat, we investigated the optimum design for minimum weight given an acceptable level of heat transmission through the insulation. The optimality criterion developed was different from that obtained for minimizing beat transfer coefficient. For minimum mass design, we had to find for a given temperature the optimum density, which minimized the logarithmic derivative of the insulation thermal conductivity with respect to its density. The logarithmic derivative is defined as the ratio of relative change in the dependent response (thermal conductivity) to the relative change in the independent variable (density). The results have been documented as a conference paper that will be presented at the upcoming AIAA.

Haftka, Raphael T.↗
NASA NTRS2006

Dynamic Theory of Relativistic Electrons Stochastic Heating by Whistler Mode Waves with Application to the Earth Magnetosphere

In the Hamiltonian approach an electron motion in a coherent packet of the whistler mode waves propagating along the direction of an ambient magnetic field is studied. The physical processes by which these particles are accelerated to high energy are established. Equations governing a particle motion by group symmetries of the problem were transformed in to a closed pair of nonlinear difference equations. The solutions of these equations have shown there exists the energetic threshold below that the electron motion is regular, and when the initial energy is above the threshold an electron moves stochastically. It is proved that the upper boundary of particle stochastic heating is conditioned by intrinsic property of the particle chaotic motion. Particle energy spectra and pitch angle electron scattering are described by the Fokker-Planck-Kolmogorov equations. It is shown that significant pitch angle diffusion occurs for the Earth radiation belt electrons with energies from a few keV up to a few MeV.

Khazanov, G. V.↗
NASA NTRS2007

Dynamic Theory of Relativistic Electrons Stochastic Heating by Whistler Mode Waves with Application to the Earth Magnetosphere

In the Hamiltonian approach an electron motion in a coherent packet of the whistler mode waves propagating along the direction of an ambient magnetic field is studied. The physical processes by which these particles are accelerated to high energy are established. Equations governing a particle motion were transformed in to a closed pair of nonlinear difference equations. The solutions of these equations have shown there exists the energetic threshold below that the electron motion is regular, and when the initial energy is above the threshold an electron moves stochastically. Particle energy spectra and pitch angle electron scattering are described by the Fokker-Planck-Kolmogorov equations. Calculating the stochastic diffusion of electrons due to a spectrum of whistler modes is presented. The parametric dependence of the diffusion coefficients on the plasma particle density, magnitude of wave field, and the strength of magnetic field is studies. It is shown that significant pitch angle diffusion occurs for the Earth radiation belt electrons with energies from a few keV up to a few MeV.

Khazanov, G. V.↗
NASA NTRS2015

Aeroelastic Analysis of SUGAR Truss-Braced Wing Wind-Tunnel Model Using FUN3D and a Nonlinear Structural Model

Considerable attention has been given in recent years to the design of highly flexible aircraft. The results of numerous studies demonstrate the significant performance benefits of strut-braced wing (SBW) and trussbraced wing (TBW) configurations. Critical aspects of the TBW configuration are its larger aspect ratio, wing span and thinner wings. These aspects increase the importance of considering fluid/structure and control system coupling. This paper presents high-fidelity Navier-Stokes simulations of the dynamic response of the flexible Boeing Subsonic Ultra Green Aircraft Research (SUGAR) truss-braced wing wind-tunnel model. The latest version of the SUGAR TBW finite element model (FEM), v.20, is used in the present simulations. Limit cycle oscillations (LCOs) of the TBW wing/strut/nacelle are simulated at angle-of-attack (AoA) values of -1, 0 and +1 degree. The modal data derived from nonlinear static aeroelastic MSC.Nastran solutions are used at AoAs of -1 and +1 degrees. The LCO amplitude is observed to be dependent on AoA. LCO amplitudes at -1 degree are larger than those at +1 degree. The LCO amplitude at zero degrees is larger than either -1 or +1 degrees. These results correlate well with both wind-tunnel data and the behavior observed in previous studies using linear aerodynamics. The LCO onset at zero degrees AoA has also been computed using unloaded v.20 FEM modes. While the v.20 model increases the dynamic pressure at which LCO onset is observed, it is found that the LCO onset at and above Mach 0.82 is much different than that produced by an earlier version of the FEM, v. 19.

Bartels, Robert E.↗
NASA NTRS2016

Transition Prediction in Hypersonic Boundary Layers Using Receptivity and Freestream Spectra

Boundary-layer transition in hypersonic flows over a straight cone can be predicted using measured freestream spectra, receptivity, and threshold values for the wall pressure fluctuations at the transition onset points. Simulations are performed for hypersonic boundary-layer flows over a 7-degree half-angle straight cone with varying bluntness at a freestream Mach number of 10. The steady and the unsteady flow fields are obtained by solving the two-dimensional Navier-Stokes equations in axisymmetric coordinates using a 5th-order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using a third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. The calculated N-factors at the transition onset location increase gradually with increasing unit Reynolds numbers for flow over a sharp cone and remain almost the same for flow over a blunt cone. The receptivity coefficient increases slightly with increasing unit Reynolds numbers. They are on the order of 4 for a sharp cone and are on the order of 1 for a blunt cone. The location of transition onset predicted from the simulation including the freestream spectrum, receptivity, and the linear and the weakly nonlinear evolutions yields a solution close to the measured onset location for the sharp cone. The simulations over-predict transition onset by about twenty percent for the blunt cone.

Balakumar, P.↗
NASA NTRS1995

Nonlinear interaction between a pair of oblique modes in a supersonic mixing layer: Long-wave limit

The nonlinear interaction between a pair of symmetric, oblique, and spatial instability modes is studied in the long-wave limit using asymptotic methods. The base flow is taken to be a supersonic mixing layer whose Mach number is such that the corresponding vortex sheet is marginally stable according to Miles' criterion. It is shown that the amplitude of the mode obeys a nonlinear integro-differential equation. Numerical solutions of this equation show that, when the obliqueness angle is less than pi/4, the effect of the nonlinearity is to enhance the growth rate of the instability. The solution terminates in a singularity at a finite streamwise location. This result is reminiscent of that obtained in the vicinity of the neutral point by other authors in several different types of flows. On the other hand, when the obliqueness angle is more than pi/4, the streamwise development of the amplitude is characterized by a series of modulations. This arises from the fact that the nonlinear term in the amplitude equation may be either stabilizing or destabilizing, depending on the value of the streamwise coordinate. However, even in this case the amplitude of the disturbance increases, though not as rapidly as in the case for which the angle is less than pi/4. Quite generally then, the nonlinear interaction between two oblique modes in a supersonic mixing layer enhances the growth of the disturbance.

Balsa, Thomas F.↗
DOE OSTI2025

Experimental observation of nonlinear relation between pressure and water flux is consistent with the solution-diffusion model

In several recent studies, it has been proposed that the fundamental understanding of penetrant transport in dense polymer membranes occurring via the solution-diffusion model, which has been the generally accepted theoretical framework for describing penetrant transport in such materials for the past several decades, is flawed. An alternate mechanistic framework based on the idea of two-phase flow in a porous medium (i.e., pore-flow) has been broadly advanced instead, with proponents of this approach claiming that the pore-flow theoretical framework provides the necessary mechanistic insight to design novel polymeric membrane materials for emerging applications. In this study, we show experimental results for hydraulic permeation of water that are entirely consistent with the solution-diffusion theory, without modification, for three dense polymeric membranes: crosslinked poly(ethylene glycol diacrylate) (XLPEGDA), Nafion 117 ionomer in the sodium counterion form (Nafion 117-Na), and cellulose acetate (CA). By measuring water flux at transmembrane pressures up to 240 bar, we observe a nonlinear relationship between the transmembrane pressure (TMP) and water flux, J w , for XLPEGDA and Nafion 117-Na, while this relationship is linear for CA. We demonstrate that the behavior of these three materials is described via the solution-diffusion model. According to the solution-diffusion model, flux is, to a good approximation, proportional to the transmembrane concentration difference induced by the pressure difference across the membrane, rather than to TMP itself. Water sorption isotherms are reported for all three materials. They further justify the nonlinear relationship between TMP and J w observed in XLPEGDA and Nafion 117-Na, emphasizing that the nonlinearity in the flux/TMP relationship stems from nonlinearities in the sorption isotherm with pressure. Additionally, the relationship between water flux and TMP can be predicted, a priori, with no adjustable parameters when a predictive model for the diffusion coefficient of water is employed in conjunction with the experimental water sorption isotherms in the solution-diffusion model. Furthermore, our results demonstrate the validity of the solution-diffusion model to describe transport of penetrants in dense polymer membranes, while highlighting the sensitivity of the solution-diffusion model to the many physical and mathematical simplifications commonly applied to the theory in literature.

materials↗
DOE OSTI2025

Uncertainty estimation of bifurcated solutions in the Rayleigh–Bénard problem for advanced nuclear reactors applications

Multiphysics models of nuclear reactors frequently comprise nonlinear systems of equations. The nonlinear nature of these models could lead to solution bifurcations, where a small change in a certain parameter, e.g., the thermophysical properties of the coolant, can lead to a sudden change in the system’s behavior. At the point in parameter space where this happens, called a critical point, the Jacobian matrix of the model’s nonlinear operator becomes singular potentially permitting multiple solutions to coexist. In this paper, we perform uncertainty estimation (UE) in a parameter range that includes bifurcated solutions within the context of Rayleigh–Bénard problem. We perform this analysis assuming uncertain temperature difference, and tilt angle for the iterative solution algorithm with a unit Prandtl number (Pr = 1). Also, we perform this analysis under uncertain thermophysical properties for both FLiBe molten salt and liquid sodium as working fluid. We deploy two approaches to compute statistical moments for the resulting distributions of selected flow-field variables. The first approach is the blind computation of the mean and the standard deviation without any consideration of solution bifurcation, while the second approach utilizes k-means clustering to cluster each branch’s solutions together and compute separate statistical moments for each branch. The statistical distributions are obtained by perturbing the selected parameters about nominal values that correspond to a solution on one of the valid branches, and that solution is used as initial guess for the iterative solution algorithm. We found that perturbation of any parameter when its nominal value is close to its critical point always leads to branch jumping, i.e., the iterations converge to a solution on a branch different from the branch of the initial guess. This produces a statistical ensemble comprised of fundamentally different solutions leading to wrong mean values and uncertainty estimates, whereas clustering provides an efficient way to deal with this type of computation. This work is important for developing Gen IV nuclear systems because many of these systems rely on natural convection for cooling especially in accident conditions.

97 - MATHEMATICS AND COMPUTING↗
NASA NTRS1987

A computational code for two-dimensional unsteady magnetohydrodynamics by the method of characteristics

A computational code for solving two-dimensional, time-dependent MHD equations by the method of characteristics is presented. Its capabilities are demonstrated by solving two very different problems for which analytical solutions exist: linearized, standing MHD wave motions in a magnetized cylindrical plasma, and nonlinear self-similar expansion of a magnetized plasma ball. The nonlinear development of standing MHD wave solutions in a cylindrical plasma is also studied. The method can be naturally embedded in the computational architecture of massively parallel processors.

Lou, Y. Q.↗
NASA NTRS1991

A vector unsymmetric eigenequation solver for nonlinear flutter analysis on high-performance computers

A finite element approach is presented for determining the nonlinear flutter characteristics of composite panels using unsteady, third-order piston theory aerodynamics. Both nonlinear structural (large-amplitude) and nonlinear aerodynamics terms are considered in the finite element formulation. Solution procedures are presented to solve the nonlinear panel flutter and the large-amplitude free vibration finite element equations. Nonlinear aerodynamic and linear structural finite element flutter results for composite panels are also presented. An efficient, vector-version generalized unsymmetric eigenequation solver is developed for large-amplitude vibration and nonlinear panel flutter analyses on high-performance computers.

Qin, Jiangning↗
NASA NTRSDate not supplied

Improvements in Iterative Convergence of FUN3D Solutions

This paper presents a hierarchical adaptive nonlinear iteration method (HANIM) implemented in the NASA computational fluid dynamics code, FUN3D, to improve robustness and computational efficiency. In contrast to the legacy FUN3D iterative solver that relies on an approximate Jacobian, a simple multicolor Gauss-Seidel point-implicit iteration scheme, and linear Courant-Friedrichs-Lewy number (CFL) ramping, HANIM is based upon a hierarchy of modules including preconditioner, generalized conjugate residual, realizability check, nonlinear control, and CFL adaption modules. HANIM performance is systematically compared with the performance of the legacy solver of FUN3D and a baseline solver based on a preconditioner alone. Iterative solutions are compared for three benchmark cases: a subsonic separated flow around a hemisphere cylinder, a supersonic flow through a long duct, and a subsonic flow over the NASA wing-fuselage juncture model. Two Reynolds-averaged Navier-Stokes turbulence models are used in these computations, namely, the negative variant of the linear one-equation Spalart-Allmaras model and its nonlinear extension based on quadratic constitutive relations.

CFD↗
NASA NTRS1993

Collisionless parallel shocks

Consideration is given to a collisionless parallel shock based on solitary-type solutions of the modified derivative nonlinear Schroedinger equation (MDNLS) for parallel Alfven waves. The standard derivative nonlinear Schroedinger equation is generalized in order to include the possible anisotropy of the plasma distribution and higher-order Korteweg-de Vies-type dispersion. Stationary solutions of MDNLS are discussed. The anisotropic nature of 'adiabatic' reflections leads to the asymmetric particle distribution in the upstream as well as in the downstream regions of the shock. As a result, nonzero heat flux appears near the front of the shock. It is shown that this causes the stochastic behavior of the nonlinear waves, which can significantly contribute to the shock thermalization.

Khabibrakhmanov, I. KH.↗
NASA NTRS1993

Efficient Minimum-Polynomial And Reduced-Rank Extrapolation

MPERRE computer program accelerates convergence of sequence of vectors by use of minimum-polynomial extrapolation (MPE) and reduced-rank extrapolation (RRE). Effective in accelerating convergences of such sequences of vectors as those obtained from iterative solution of systems of linear and nonlinear equations. In conjunction with various iterative techniques, successfully employed in finite-difference solution of large-scale elliptic boundary-value problems and problems in computational-fluid-dynamics. Only input required is sequence of vectors, convergence of which is accelerated. Program economical and easy to use. Written in FORTRAN 77.

Sidi, Avram↗