Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical partial differential equations”

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 NTRS2018

Instability Wave–Streak Interactions in a High Mach Number Boundary Layer at Flight Conditions

The interaction of stationary streaks undergoing nonmodal growth with modally unstable instability waves in a hypersonic boundary-layer flow is studied using numerical computations. The geometry and flow conditions are selected to match a relevant trajectory location from the ascent phase of the HIFiRE-1 flight experiment; namely, a 7 degree half-angle, circular cone with 2.5 mm nose radius, freestream Mach number equal to 5.30, unit Reynolds number equal to 13.42 m -1 , and wall-to-adiabatic temperature ratio of approximately 0.35 over most of the vehicle. This paper investigates the nonlinear evolution of initially linear optimal disturbances that evolve into finite-amplitude streaks, followed by an analysis of the modal instability characteristics of the perturbed, streaky boundary-layer flow. The investigation is performed with stationary direct numerical simulations (DNS) and plane-marching parabolized stability equations (PSE), in conjunction with partial-differential-equation-based planar eigenvalue analysis. The overall effect of streaks is to reduce the peak amplification factors of instability waves, indicating a possible downstream shift in the onset of laminar-turbulent transition. The present study confirms previous findings that the mean flow distortion of the nonlinear streak perturbation reduces the amplification rates of the Mack-mode instability. More importantly, however, the present results demonstrate that the spanwise varying component of the streak can produce a larger effect on the Mack-mode amplification. The study with selected azimuthal wavenumbers for the stationary streaks reveals that a wavenumber of approximately 1.4 times larger than the optimal wavenumber is more effective in stabilizing the planar Mack-mode instabilities. In the absence of unstable first-mode waves for the present cold-wall condition, transition onset is expected to be delayed until the peak streak amplitude increases to nearly 35 percent of the freestream velocity, when intrinsic instabilities of the boundary-layer streaks begin to dominate the transition process. For streak amplitudes below that limit a significant net stabilization is achieved, yielding a potential transition delay that can exceed 100 percent of the length of the laminar region in the uncontrolled case.

Boundary layer stability↗
NASA NTRS1998

Flow Applications of the Least Squares Finite Element Method

The main thrust of the effort has been towards the development, analysis and implementation of the least-squares finite element method (LSFEM) for fluid dynamics and electromagnetics applications. In the past year, there were four major accomplishments: 1) special treatments in computational fluid dynamics and computational electromagnetics, such as upwinding, numerical dissipation, staggered grid, non-equal order elements, operator splitting and preconditioning, edge elements, and vector potential are unnecessary; 2) the analysis of the LSFEM for most partial differential equations can be based on the bounded inverse theorem; 3) the finite difference and finite volume algorithms solve only two Maxwell equations and ignore the divergence equations; and 4) the first numerical simulation of three-dimensional Marangoni-Benard convection was performed using the LSFEM.

Jiang, Bo-Nan↗
NASA NTRS1973

An improved numerical process for solution of solid mechanics problems

This paper gives an overview of the development and status of an improved numerical process for the solution of solid mechanics problems. The proposed process uses a mixed formulation with the fundamental unknowns consisting of both stress and displacement parameters. The problem is formulated either by means of first-order partial differential equations or in a variational form by using a Hellinger-Reissner-type mixed variational principle. For presentation purposes, the components of a numerical process are characterized and the criteria for an ideal process are outlined. Commonly used finite-difference and finite-element procedures are examined in the light of these criteria and it is shown that they fall short in a number of ways. The proposed numerical process, on the other hand, satisfies most of the optimality criteria and appears to be particularly suited for use with the forthcoming generation computers.

Noor, A. K.↗
NASA NTRS1991

Asynchronous and corrected-asynchronous numerical solutions of parabolic PDES on MIMD multiprocessors

A major problem in achieving significant speed-up on parallel machines is the overhead involved with synchronizing the concurrent process. Removing the synchronization constraint has the potential of speeding up the computation. The authors present asynchronous (AS) and corrected-asynchronous (CA) finite difference schemes for the multi-dimensional heat equation. Although the discussion concentrates on the Euler scheme for the solution of the heat equation, it has the potential for being extended to other schemes and other parabolic partial differential equations (PDEs). These schemes are analyzed and implemented on the shared memory multi-user Sequent Balance machine. Numerical results for one and two dimensional problems are presented. It is shown experimentally that the synchronization penalty can be about 50 percent of run time: in most cases, the asynchronous scheme runs twice as fast as the parallel synchronous scheme. In general, the efficiency of the parallel schemes increases with processor load, with the time level, and with the problem dimension. The efficiency of the AS may reach 90 percent and over, but it provides accurate results only for steady-state values. The CA, on the other hand, is less efficient, but provides more accurate results for intermediate (non steady-state) values.

Amitai, Dganit↗
NASA NTRS1977

Prediction of hydrodynamics and chemistry of confined turbulent methane-air flames with attention to formation of oxides of nitrogen

A formulation of the governing partial differential equations for fluid flow and reacting chemical species in a tubular combustor is presented. A numerical procedure for the solution of the governing differential equations is described, and models for chemical equilibrium and chemical kinetics calculations are presented. The chemical equilibrium model is used to characterize the hydrocarbon reactions. The chemical kinetics model is used to predict the concentrations of the oxides of nitrogen. The combustor consists of a cylindrical duct of varying cross sections with concentric streams of gaseous fuel and air entering the duct at one end. Four sample cases with specified inlet and boundary conditions are considered, and the results are discussed

Elghobashi, S.↗
NASA NTRS1998

Nonlinear Interaction of Detuned Instability Waves in Boundary-Layer Transition: Resonant-Triad Interaction

The non-equilibrium critical-layer analysis of a system of frequency-detuned resonant-triads is presented using the generalized scaling of Lee. It is shown that resonant-triads can interact nonlinearly within the common critical layer when their (fundamental) Strouhal numbers are different by a factor whose magnitude is of the order of the growth rate multiplied by the wavenumber of the instability wave. Since the growth rates of the instability modes become larger and the critical layers become thicker as the instability waves propagate downstream, the frequency-detuned resonant-triads that grow independently of each other in the upstream region can interact nonlinearly in the later downstream stage. In the final stage of the non-equilibrium critical-layer evolution, a wide range of instability waves with the scaled frequencies differing by almost an Order of (l) can nonlinearly interact. Low-frequency modes are also generated by the nonlinear interaction between oblique waves in the critical layer. The system of partial differential critical-layer equations along with the jump equations are presented here. The amplitude equations with their numerical solutions are given in Part 2. The nonlinearly generated low-frequency components are also investigated in Part 2.

Lee, Sang Soo↗
NASA NTRS1999

The Thin Oil Film Equation

A thin film of oil on a surface responds primarily to the wall shear stress generated on that surface by a three-dimensional flow. The oil film is also subject to wall pressure gradients, surface tension effects and gravity. The partial differential equation governing the oil film flow is shown to be related to Burgers' equation. Analytical and numerical methods for solving the thin oil film equation are presented. A direct numerical solver is developed where the wall shear stress variation on the surface is known and which solves for the oil film thickness spatial and time variation on the surface. An inverse numerical solver is also developed where the oil film thickness spatial variation over the surface at two discrete times is known and which solves for the wall shear stress variation over the test surface. A One-Time-Level inverse solver is also demonstrated. The inverse numerical solver provides a mathematically rigorous basis for an improved form of a wall shear stress instrument suitable for application to complex three-dimensional flows. To demonstrate the complexity of flows for which these oil film methods are now suitable, extensive examination is accomplished for these analytical and numerical methods as applied to a thin oil film in the vicinity of a three-dimensional saddle of separation.

Brown, James L.↗
NASA NTRS1992

On an origin of numerical diffusion: Violation of invariance under space-time inversion

The invariant properties of the convection equation du/dt + adu/dx = 0 (where d is the partial differential operator) with respect to spatial reflection, time reversal, and space-time inversion are studied. Generally, a finite-difference analog of this equation may possess some or none of these properties. It is shown that, under certain conditions, the von Neumann amplification factor of an analog satisfies a special relation for each invariant property this analog possesses. Particularly, an analog is neutrally stable and thus free of numerical diffusion if it possesses the invariant property related to space-time inversion. It is also explained why generally (1) an upwind scheme possesses neither the invariant property related to spatial reflection nor that related to space-time inversion, and (2) an explicit scheme possesses neither the invariant property related to time reversal nor that related to space-time inversion. Extension to the viscous case and a remarkable connection between the current work and a new numerical framework for solving conservation laws are also discussed.

Chang, Sin-Chung↗
NASA NTRS1980

Computation of three-dimensional viscous supersonic flow in inlets

A new approach has been developed for the computation of the three-dimensional viscous supersonic flow with embedded subsonic regions adjacent to solid boundaries and is applied to a mixed-compression supersonic inlet typical of current designs. The approach uses a reduced form of the three-dimensional Navier-Stokes equations so that the resultant equations can be treated as an initial boundary value problem and thus be solved by non-iterative forward marching in space. The numerical procedure utilizes an efficient consistently-split linearized block implicit technique to solve the finite difference analogues to the set of governing partial differential equations.

Buggeln, R. C.↗
NASA NTRS1983

Legendre-Tau approximations for functional differential equations

The numerical approximation of solutions to linear functional differential equations are considered using the so called Legendre tau method. The functional differential equation is first reformulated as a partial differential equation with a nonlocal boundary condition involving time differentiation. The approximate solution is then represented as a truncated Legendre series with time varying coefficients which satisfy a certain system of ordinary differential equations. The method is very easy to code and yields very accurate approximations. Convergence is established, various numerical examples are presented, and comparison between the latter and cubic spline approximations is made.

Ito, K.↗
NASA NTRS1995

Analysis of discretization errors in LES

All numerical simulations of turbulence (DNS or LES) involve some discretization errors. The integrity of such simulations therefore depend on our ability to quantify and control such errors. In the classical literature on analysis of errors in partial differential equations, one typically studies simple linear equations (such as the wave equation or Laplace's equation). The qualitative insight gained from studying such simple situations is then used to design numerical methods for more complex problems such as the Navier-Stokes equations. Though such an approach may seem reasonable as a first approximation, it should be recognized that strongly nonlinear problems, such as turbulence, have a feature that is absent in linear problems. This feature is the simultaneous presence of a continuum of space and time scales. Thus, in an analysis of errors in the one dimensional wave equation, one may, without loss of generality, rescale the equations so that the dependent variable is always of order unity. This is not possible in the turbulence problem since the amplitudes of the Fourier modes of the velocity field have a continuous distribution. The objective of the present research is to provide some quantitative measures of numerical errors in such situations. Though the focus of this work is LES, the methods introduced here can be just as easily applied to DNS. Errors due to discretization of the time-variable are neglected for the purpose of this analysis.

Ghosal, Sandip↗
NASA NTRS1963

Application of Dorodnitsyn's Integral Method to Nonequilibrium Flows over Pointed Bodies

Dorodnitsyn's integral method is used to obtain an approximate solution to the supersonic nonequilibrium flow over pointed bodies with attached shock waves. The partial differential equations governing the flow are converted to an approximate set of ordinary equations, which are solved by numerical integration starting at the body tip. Detailed analytical and numerical results for the first approximation are presented, considering the vibrational relaxation of a diatomic gas over a wedge or cone. It is shown that the first approximation yields: (1) The exact flow-variable gradients at the wedge tip (2) Expressions for the flow-variable gradients at the cone tip which are in agreement with extrapolations of characteristics calculations (3) A good approximate algebraic solution for frozen or equilibrium conical flow (4) An approximate expression for the nonequilibrium-flow stream function which affords a means of obtaining variations across the shock layer of the temperature and vibrational energy. Numerical results for both the wedge and cone compare favorably with identical cases computed by the method of characteristics.

NONEQUILIBRIUM FLOW↗
NASA NTRS1980

Errors in finite-difference computations on curvilinear coordinate systems

Curvilinear coordinate systems were used extensively to solve partial differential equations on arbitrary regions. An analysis of truncation error in the computation of derivatives revealed why numerical results may be erroneous. A more accurate method of computing derivatives is presented.

Mastin, C. W.↗
NASA NTRS1985

Bridgman growth of semiconductors

The purpose of this study was to improve the understanding of the transport phenomena which occurs in the directional solidification of alloy semiconductors. In particular, emphasis was placed on the strong role of convection in the melt. Analytical solutions were not deemed possible for such an involved problem. Accordingly, a numerical model of the process was developed which simulated the transport. This translates into solving the partial differential equations of energy, mass, species, and momentum transfer subject to various boundary and initial conditions. A finite element method with simple elements was initially chosen. This simulation tool will enable the crystal grower to systematically identify and modify the important design factors within her control to produce better crystals.

Carlson, F. M.↗
NASA NTRS1991

Calculations of separated 3-D flows with a pressure-staggered Navier-Stokes equations solver

A Navier-Stokes equations solver based on a pressure correction method with a pressure-staggered mesh and calculations of separated three-dimensional flows are presented. It is shown that the velocity pressure decoupling, which occurs when various pressure correction algorithms are used for pressure-staggered meshes, is caused by the ill-conditioned discrete pressure correction equation. The use of a partial differential equation for the incremental pressure eliminates the velocity pressure decoupling mechanism by itself and yields accurate numerical results. Example flows considered are a three-dimensional lid driven cavity flow and a laminar flow through a 90 degree bend square duct. For the lid driven cavity flow, the present numerical results compare more favorably with the measured data than those obtained using a formally third order accurate quadratic upwind interpolation scheme. For the curved duct flow, the present numerical method yields a grid independent solution with a very small number of grid points. The calculated velocity profiles are in good agreement with the measured data.

Kim, S.-W.↗