Search NASA⌕ Search

SEARCH · Search NASA

Results for “parabolic 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 NTRS1989

Goertler vortices in supersonic and hypersonic boundary layers

The problem of Goertler vortices in compressible boundary layers over concave walls is considered. At O(1) wavelengths, the instability is governed by parabolic partial differential equations that are solved numerically to determine the effect of various initial conditions on the development of Goertler vortex instability in compressible boundary layers. The results show that both the velocity and temperature fluctuations may lead to a Goertler vortex. The vortex growth rates determined from the present method are found to differ somewhat from those given by a normal mode solution. At both the supersonic and hypersonic Mach numbers, cooling has a small destabilizing effect. In addition, the most unstable disturbances shift toward lower wavelengths because of thinning of the boundary layer. The results also show that compressibility has a stabilizing effect on the Goertler instability, while the effect of an adverse pressure gradient is found to be destabilizing. The behavior of the Goertler vortex structure with Mach number is also examined. At hypersonic Mach numbers, vortices are located near the edge of the boundary layer for adiabatic wall conditions. However, the entire boundary layer is affected when the wall is cooled.

Spall, R. E.↗
NASA NTRS1979

Utility of a finite element solution algorithm for initial-value problems

The Galerkin criterion within a finite element Weighted Residuals formulation is employed to establish an implicit solution algorithm for an initial-value partial differential equation. Numerical solutions of a transient parabolic and a hyperbolic equation, obtained using linear, quadratic and two cubic finite element basis functions, are employed to quantize accuracy and confirm and refine theoretical convergence rate estimates. The linear basis algorithm for the hyperbolic equation displays excellent accuracy on a coarse computational grid and a high-order convergence rate with discretization refinement. Good accuracy and a strong convergence rate in surface flux are determined for a nonhomogeneous Neumann boundary constraint applied to a parabolic equation. The results amply demonstrate the impact of the nondiagonal finite element initial-value matrix structure on solution accuracy and/or convergence rate.

Baker, A. J.↗
NASA NTRS1984

Fully-coupled analysis of jet mixing problems. Part 1. Shock-capturing model, SCIPVIS

A computational model, SCIPVIS, is described which predicts the multiple cell shock structure in imperfectly expanded, turbulent, axisymmetric jets. The model spatially integrates the parabolized Navier-Stokes jet mixing equations using a shock-capturing approach in supersonic flow regions and a pressure-split approximation in subsonic flow regions. The regions are coupled using a viscous-characteristic procedure. Turbulence processes are represented via the solution of compressibility-corrected two-equation turbulence models. The formation of Mach discs in the jet and the interactive analysis of the wake-like mixing process occurring behind Mach discs is handled in a rigorous manner. Calculations are presented exhibiting the fundamental interactive processes occurring in supersonic jets and the model is assessed via comparisons with detailed laboratory data for a variety of under- and overexpanded jets.

Dash, S. M.↗
DOE OSTI2018

Solving high-dimensional partial differential equations using deep learning

Developing algorithms for solving high-dimensional partial differential equations (PDEs) has been an exceedingly difficult task for a long time, due to the notoriously difficult problem known as the “curse of dimensionality.” This paper introduces a deep learning-based approach that can handle general high-dimensional parabolic PDEs. To this end, the PDEs are reformulated using backward stochastic differential equations and the gradient of the unknown solution is approximated by neural networks, very much in the spirit of deep reinforcement learning with the gradient acting as the policy function. Numerical results on examples including the nonlinear Black–Scholes equation, the Hamilton–Jacobi–Bellman equation, and the Allen–Cahn equation suggest that the proposed algorithm is quite effective in high dimensions, in terms of both accuracy and cost. Furthermore, this opens up possibilities in economics, finance, operational research, and physics, by considering all participating agents, assets, resources, or particles together at the same time, instead of making ad hoc assumptions on their interrelationships.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

An approximation theory for the identification of linear thermoelastic systems

An abstract approximation framework and convergence theory for the identification of thermoelastic systems is developed. Starting from an abstract operator formulation consisting of a coupled second order hyperbolic equation of elasticity and first order parabolic equation for heat conduction, well-posedness is established using linear semigroup theory in Hilbert space, and a class of parameter estimation problems is then defined involving mild solutions. The approximation framework is based upon generic Galerkin approximation of the mild solutions, and convergence of solutions of the resulting sequence of approximating finite dimensional parameter identification problems to a solution of the original infinite dimensional inverse problem is established using approximation results for operator semigroups. An example involving the basic equations of one dimensional linear thermoelasticity and a linear spline based scheme are discussed. Numerical results indicate how the approach might be used in a study of damping mechanisms in flexible structures.

Rosen, I. G.↗
NASA NTRS1993

Stability and transition on swept wings

This paper describes the extension and application of the Parabolized Stability Equations (PSE) to the stability and transition of the supersonic three-dimensional laminar boundary layer on a swept wing. The problem formulation uses a general coordinate transformation for arbitrary curvilinear body-fitted computational grids. Some testing using these coordinates is briefly described to help validate the software used for the investigation. The disturbance amplitude ratios as a function of chord position for supersonic (Mach 1.5) boundary layers on untapered, untwisted wings of different sweep angles are then presented and compared with those obtained from local parallel analyses.

Stuckert, Greg↗
NASA NTRS1987

Approximation techniques for domain identification in two-dimensional parabolic systems under boundary observations

The identification of the geometrical structure of the system boundary for a two-dimensional diffusion system is considered. The domain identification treated is converted into an optimization problem based on a fit-to-data criterion and theoretical convergence results for approximate identification techniques are discussed. Results of numerical experiments to demonstrate the efficiency of the theoretical ideas are reported.

Banks, H. T.↗
NASA NTRS1978

Evaluation of the three-dimensional parabolic flow computer program SHIP

The three-dimensional parabolic flow program SHIP designed for predicting supersonic combustor flow fields is evaluated to determine its capabilities. The mathematical foundation and numerical procedure are reviewed; simplifications are pointed out and commented upon. The program is then evaluated numerically by applying it to several subsonic and supersonic, turbulent, reacting and nonreacting flow problems. Computational results are compared with available experimental or other analytical data. Good agreements are obtained when the simplifications on which the program is based are justified. Limitations of the program and the needs for improvement and extension are pointed out. The present three dimensional parabolic flow program appears to be potentially useful for the development of supersonic combustors.

Pan, Y. S.↗
DOE OSTI2026

Explicit Monotone Stable Super-Time-stepping Methods for Finite Time Singularities

We explore a novel way to numerically resolve the scaling behavior of finite-time singularities in solutions of nonlinear parabolic PDEs. The Runge–Kutta–Legendre (RKL) and Runge–Kutta–Gegenbauer (RKG) super-time-stepping methods were originally developed for nonlinear complex physics problems with diffusion. These are multistage single step second-order, forward-in-time methods with no implicit solves. The advantage is that the time-step size for stability scales with stage number 𝑠 as $\mathcal{O}$⁡(𝑠 2 ). Many interesting nonlinear PDEs have finite-time singularities, and the presence of diffusion often limits one to using implicit or semi-implicit time-step methods for stability constraints. Finite-time singularities are particularly challenging due to the large range of scales that one desires to resolve, often with adaptive spatial grids and adaptive time steps. Here, in this study, we show two examples of nonlinear PDEs for which the self-similar singularity structure has time and space scales that are resolvable using the RKL and RKG methods, without forcing even smaller time steps. Compared to commonly used implicit numerical methods, we achieve a significantly smaller run time while maintaining comparable accuracy. We also prove numerical monotonicity for both the RKL and RKG methods under their linear stability conditions for the constant coefficient heat equation, in the case of infinite domain and periodic boundary condition, leading to a theoretical guarantee of the superiority of the RKL and RKG methods over traditional super-time-stepping methods, such as the Runge-Kutta-Chebyshev and the orthogonal Runge-Kutta-Chebyshev methods. Code can be found at https://github.com/ZT220501/SRK-Singularity.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1982

The passage of a star by a massive black hole

The gridless, smoothed particle hydrodynamic method of Lucy (1977) and Gingold and Monaghan (1977, 1980, a, b, 1981) is used to solve the three-dimensional problem posed by the effects on a 1.0 solar mass star of a 10,000 solar mass black hole's passage, in an initially parabolic orbit with a variety of pericenter distances. It is found that the tidal forces induce rotation and radial and nonradial pulsations in the star, and that the loss of orbital energy to these internal motions leads to the star's capture by the black hole. While the outer layers of the star are disrupted at small pericenter distances, the entire star is destroyed at even smaller distances. These results are applicable to X-ray sources, active galactic nuclei, and quasars

Nolthenius, R. A.↗
NASA NTRS1988

Fully-coupled analysis of jet mixing problems. Three-dimensional PNS model, SCIP3D

Numerical procedures formulated for the analysis of 3D jet mixing problems, as incorporated in the computer model, SCIP3D, are described. The overall methodology closely parallels that developed in the earlier 2D axisymmetric jet mixing model, SCIPVIS. SCIP3D integrates the 3D parabolized Navier-Stokes (PNS) jet mixing equations, cast in mapped cartesian or cylindrical coordinates, employing the explicit MacCormack Algorithm. A pressure split variant of this algorithm is employed in subsonic regions with a sublayer approximation utilized for treating the streamwise pressure component. SCIP3D contains both the ks and kW turbulence models, and employs a two component mixture approach to treat jet exhausts of arbitrary composition. Specialized grid procedures are used to adjust the grid growth in accordance with the growth of the jet, including a hybrid cartesian/cylindrical grid procedure for rectangular jets which moves the hybrid coordinate origin towards the flow origin as the jet transitions from a rectangular to circular shape. Numerous calculations are presented for rectangular mixing problems, as well as for a variety of basic unit problems exhibiting overall capabilities of SCIP3D.

Wolf, D. E.↗
NASA NTRS1990

Explicit-implicit staggered procedure for multibody dynamics analysis

An explicit-implicit staggered time-integration procedure is presented for the solution of multibody dynamical equations involving large rotations and constraints. The algorithms adopts a two-stage modification of the central difference algorithm for integrating the translational coordinates and the angular velocity vector, and the midpoint implicit algorithm to solve the kinematical relation in terms of the Euler parameters for updating the angular orientations. The Lagrange multipliers to enforce the system constraints are obtained by implicitly integrating a parabolically regularized differential equation for the multipliers. The performance of the present procedure has been evaluated to applying the procedure to solve several sample problems. The results indicate that the procedure is robust in dealing with a variety of constraints and spatial kinematic motions, hence it is recommended for applications to general multibody dynamics analyses.

Park, K. C.↗
NASA NTRS1999

On the Stability of Three-Dimensional Boundary Layers: Linear and Nonlinear Stability - Part 1

The primary stability of incompressible three-dimensional boundary layers is investigated using the Parabolized Stability Equations (PSE). We compute the evolution of stationary and traveling disturbances in the linear and nonlinear region prior to transition. As model problems, we choose Swept Hiemenz Flow and the DLR Transition Experiment. The primary stability results for Swept Hiemenz Flow agree very well with computations by Malik et al. For the DLR Experiment, the mean flow profiles are obtained by solving the boundary layer equations for the measured pressure distribution. Both linear and nonlinear results show very good agreement with the experiment.

Janke, Erik↗
NASA NTRS1995

Wing-section optimization for supersonic viscous flow

To improve the shape of a supersonic wing, an automated method that also includes higher fidelity to the flow physics is desirable. With this impetus, an aerodynamic optimization methodology incorporating thin-layer Navier-Stokes equations and sensitivity analysis had been previously developed. Prior to embarking upon the wind design task, the present investigation concentrated on testing the feasibility of the methodology, and the identification of adequate problem formulations, by defining two-dimensional, cost-effective test cases. Starting with two distinctly different initial airfoils, two independent shape optimizations resulted in shapes with similar features: slightly cambered, parabolic profiles with sharp leading- and trailing-edges. Secondly, the normal section to the subsonic portion of the leading edge, which had a high normal angle-of-attack, was considered. The optimization resulted in a shape with twist and camber which eliminated the adverse pressure gradient, hence, exploiting the leading-edge thrust. The wing section shapes obtained in all the test cases had the features predicted by previous studies. Therefore, it was concluded that the flowfield analyses and sensitivity coefficients were computed and fed to the present gradient-based optimizer correctly. Also, as a result of the present two-dimensional study, suggestions were made for the problem formulations which should contribute to an effective wing shape optimization.

Item, Cem C.↗
NASA NTRS1975

Application of local linearization and the transonic equivalence rule to the flow about slender analytic bodies at Mach numbers near 1.0

The local linearization method for axisymmetric flow is combined with the transonic equivalence rule to calculate pressure distribution on slender bodies at free-stream Mach numbers from .8 to 1.2. This is an approximate solution to the transonic flow problem which yields results applicable during the preliminary design stages of a configuration development. The method can be used to determine the aerodynamic loads on parabolic arc bodies having either circular or elliptical cross sections. It is particularly useful in predicting pressure distributions and normal force distributions along the body at small angles of attack. The equations discussed may be extended to include wing-body combinations.

Tyson, R. W.↗
NASA NTRS1978

A high order accurate finite element algorithm for high Reynolds number flow prediction

A Galerkin-weighted residuals formulation is employed to establish an implicit finite element solution algorithm for generally nonlinear initial-boundary value problems. Solution accuracy, and convergence rate with discretization refinement, are quantized in several error norms, by a systematic study of numerical solutions to several nonlinear parabolic and a hyperbolic partial differential equation characteristic of the equations governing fluid flows. Solutions are generated using selective linear, quadratic and cubic basis functions. Richardson extrapolation is employed to generate a higher-order accurate solution to facilitate isolation of truncation error in all norms. Extension of the mathematical theory underlying accuracy and convergence concepts for linear elliptic equations is predicted for equations characteristic of laminar and turbulent fluid flows at nonmodest Reynolds number. The nondiagonal initial-value matrix structure introduced by the finite element theory is determined intrinsic to improved solution accuracy and convergence. A factored Jacobian iteration algorithm is derived and evaluated to yield a consequential reduction in both computer storage and execution CPU requirements while retaining solution accuracy.

Baker, A. J.↗
NASA NTRS1985

Two dimensional and axisymmetric heat transfer results with the CSCM upwind implicit algorithm

The single level conservative supra characteristic method (CSCM-S) is an upwind implicit method that is very robust and in a single level relaxation format makes efficient use of computer resources, and is very rapidly convergent. These factors render the unconditionally stable algorithm a useful compromise between Parabolized Navier-Stokes methods and two level linearized implicit time dependent methods and, thus, the method is very attractive for multidimensional Navier-Stokes problems. In this paper the heat transfer prediction capabilities of the CSCM schemes for flat plate, nozzle, blunt cone and bluff body flows are reported for the first time. Results are compared with available experimental heat transfer measurements.

Bardina, J.↗
NASA NTRS1989

A gravity assist primer

An account is given of the method by which the 'energy gain' accruing to a spacecraft as a result of its 'gravity-assist', parabolic-trajectory flyby of a massive body, such as a planet. The procedure begins with the solution of the two-body portion of the problem, and the results thus obtained are used to calculate changes with respect to the other massive body in the overall scenario, namely the sun. Attention is given to the 'vector diagram' often used to display the gravity-assist effect. The present procedure is noted to be reasonably accurate for flybys in which the plane of the spacecraft's trajectory is approximately the same as that of the planet's orbit around the sun, or the ecliptic plane; this reduces the problem to one in two dimensions.

Cesarone, R. J.↗