SEARCH · Search NASA
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.
Boundary and interface conditions are derived for high order finite difference methods applied to multidimensional linear problems in curvilinear coordinates. The boundary and interface conditions lead to conservative schemes and strict and strong stability provided that certain metric conditions are met.
The upwind difference schemes of Godunov, Osher, Roe and van Leer are able to resolve one dimensional steady shocks for the Euler equations within one or two mesh intervals. Unfortunately, this resolution is lost in two dimensions when the shock crosses the computing grid at an oblique angle. To correct this problem, a numerical scheme was developed which automatically locates the angle at which a shock might be expected to cross the computing grid and then constructs separate finite difference formulas for the flux components normal and tangential to this direction. Numerical results which illustrate the ability of this method to resolve steady oblique shocks are presented.
The upwind difference schemes of Godunov, Osher, Roe and van Leer are able to resolve one dimensional steady shocks for the Euler equations within one or two mesh intervals. Unfortunately, this resolution is lost in two dimensions when the shock crosses the computing grid at an oblique angle. To correct this problem, a numerical scheme was developed which automatically locates the angle at which a shock might be expected to cross the computing grid and then constructs separate finite difference formulas for the flux components normal and tangential to this direction. Numerical results which illustrate the ability of this method to resolve steady oblique shocks are presented.
We outline a methodology for the simulation of particle-laden flows whereby the dispersed and fluid phases are two-way coupled. The drag force which couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. The method is general in that it can be extended to other partial differential equations which may be associated with particle-laden flows, such as the thermal energy equation or Maxwell's equations. In this work, we demonstrate the method of discrete Green's functions by obtaining these functions for the Navier-Stokes equations at low particle Reynolds number in a two-plane channel geometry. We perform verification at low and finite particle Reynolds number for the case of a point-particle settling under gravity parallel to a plane wall, for different wall normal separations. In comparing to other point-particle schemes the discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations and is the most accurate in the near wall region at finite Reynolds number. We discuss how the accuracy away from the wall at finite Reynolds number could be improved by appealing to Oseen-like discrete Green's functions. Finally we demonstrate that the discrete Green's function approach can have important implications on statistics of particle-laden turbulent channel flow.
An accurate numerical solution of diffusion problems containing large local gradients can be obtained with a significant reduction in computational time by using a multigrid computational scheme. The spatial domain is covered with sets of uniform square grids of different sizes. The finer grid patterns overlap the coarse grid patterns. The finite-difference expressions for each grid pattern are solved independently by iterative techniques. Two interpolation methods were used to establish the values of the potential function on the fine grid boundaries with information obtained from the coarse grid solution. The accuracy and computational requirements for solving a test problem by a simple multigrid and a multilevel-multigrid method were compared. The multilevel-multigrid method combined with a Taylor series interpolation scheme was found to be best.
The considerations and the resultant approach used to implement design sensitivity capability for grids into a large scale, general purpose finite element system (MSC/NASTRAN) are presented. The design variables are grid perturbations with a rather general linking capability. Moreover, shape and sizing variables may be linked together. The design is general enough to facilitate geometric modeling techniques for generating design variable linking schemes in an easy and straightforward manner. Test cases have been run and validated by comparison with the overall finite difference method. The linking of a design sensitivity capability for shape variables in MSC/NASTRAN with an optimizer would give designers a powerful, automated tool to carry out practical optimization design of real life, complicated structures.
Preliminary results are presented of a finite element/finite difference method (semidiscrete Galerkin method) used to calculate compressible boundary layer flow about airfoils, in which the group finite element scheme is applied to the Dorodnitsyn formulation of the boundary layer equations. The semidiscrete Galerkin (SDG) method promises to be fast, accurate and computationally efficient. The SDG method can also be applied to any smoothly connected airfoil shape without modification and possesses the potential capability of calculating boundary layer solutions beyond flow separation. Results are presented for low speed laminar flow past a circular cylinder and past a NACA 0012 airfoil at zero angle of attack at a Mach number of 0.5. Also shown are results for compressible flow past a flat plate for a Mach number range of 0 to 10 and results for incompressible turbulent flow past a flat plate. All numerical solutions assume an attached boundary layer.
Flow through a turbine annular cascade is calculated using a three-dimensional Euler method based on unstructured tetrahedral meshes. The equations are integrated in time using an explicit Runge-Kutta time-stepping scheme. The inviscid flux terms are discretized using a cell-centered finite-volume formulation with upwind flux-difference splitting. The tetrahedral meshes around the turbine blade are generated using an advancing-front technique with forced geometric periodicity between the blades. Good agreement is obtained between the present calculation and the experiment for both surface pressure distribution and flow behavior in the passage between the blades, demonstrating the capability of the present methodology for turbomachinery flow applications.
A rational interpolation method for approximating a frequency response is presented. The method is based on a product formulation of finite differences, thereby avoiding the numerical problems incurred by near-equal-valued subtraction. Also, resonant pole and zero cancellation schemes are developed that increase the accuracy and efficiency of the interpolation method. Selection techniques of interpolation points are also discussed.
We derive interpolation formulae for a third-order finite difference method in curvilinear, orthogonal coordinate systems. These formulae serve as a supplement to Colella and Woodward's PPM scheme for problems where the coordinate origin is included in the computational domain. Numerical examples of the improved accuracy of the advection scheme near coordinate singularities are shown.
A computational procedure is presented for evaluating the sensitivity coefficients of the dynamic axisymmetric response of viscoplastic shells of revolution. The analytical formulation is based on Reissner's large deformation shell theory with the effects of transverse shear deformation, rotatory inertia and moments turning around the normal to the middle surface included. The material model is chosen to be isothermal viscoplasticity, and an associated flow rule is used with a von Mises effective stress. A mixed formulation is used with the fundamental unknowns consisting of six stress resultants, three generalized displacements and three velocity components. Spatial discretization is performed using finite elements, with discontinuous stress resultants across element interfaces. The temporal integration is performed by using an explicit central difference scheme (leap-frog method) with an implicit constitutive update. The sensitivity coefficients are evaluated using a direct differentiation approach. Numerical results are presented for a spherical cap subjected to step loading, and a circular plate subjected to impulsive loading. The sensitivity coefficients are generated by evaluating the derivatives of the response quantities with respect to thickness, mass density, Young's modulus, and two of the material parameters characterizing the viscoplastic response. Time histories of the response and sensitivity coefficients are presented, along with spatial distributions of these quantities at selected times.
A numerical procedure was developed for the finite-volume solution of the Euler equations on unstructured triangular meshes based on a flux-difference split upwind method. Techniques for implementing Roe's (1985) approximate Reimann solver together with the preprocessing MUSCL differencing on unstructured grids are presented. Applications and comparisons with structured grid problems are carried out for a supersonic shock reflection problem, the supersonic flow over a blunt body, the transonic flow over NACA 0012 and RAE 2822 airfoils, and the flow about a double element Karman-Trefftz airfoil.
The method of pseudocompressibility is tested for its accuracy in solving the incompressible Navier-Stokes equations. An implicit, finite difference computer code is used to solve the equations in a three-dimensional, curvilinear coordinate system. The code employs artificial compressibiity for solving the pressure field, coupled with an implicit, approximate factorization scheme. This coupling is known as the pseudocompressibility method and introduces pressure waves of finite speed into the fluid medium that would otherwise have an infinite sound speed. The waves die out as the solution converges, and the steady state solution approaches a divergence-free condition. These waves, however, limit the time accuracy of the computations. The effects of these waves are analysed and criteria are set for choosing the pseudocompressibility parameter that governs the pressure wave speed; test cases are presented that verify these criteria. The code is tested by computing laminar flow over a two-dimensional, backward-facing step and over a two-dimensional, impulsively started circular cylinder.
A time accurate, three-dimensional, finite volume, high resolution scheme for solving the compressible full Navier-Stokes equations is presented. The present derivation is based on the upwind split formulas, specifically with the application of Roe's (1981) flux difference splitting. A high-order accurate (up to the third order) upwind interpolation formula for the inviscid terms is derived to account for nonuniform meshes. For the viscous terms, discretizations consistent with the finite volume concept are described. A variant of second-order time accurate method is proposed that utilizes identical procedures in both the predictor and corrector steps. Avoiding the definition of midpoint gives a consistent and easy procedure, in the framework of finite volume discretization, for treating viscous transport terms in the curvilinear coordinates. For the boundary cells, a new treatment is introduced that not only avoids the use of 'ghost cells' and the associated problems, but also satisfies the tangency conditions exactly and allows easy definition of viscous transport terms at the first interface next to the boundary cells. Numerical tests of steady and unsteady high speed flows show that the present scheme gives accurate solutions.
A filtering operation, based on B-spline discretizations, is introduced to target weakly growing mesh-scale oscillations that can arise in high-fidelity turbulence simulations. This is a spectral regularization that can be described using the singular values of a banded matrix operator, with the filtering strength set by a scalar- or vector-valued penalty parameter. The penalty parameter can be specified though it can also be advantageously selected to minimize the generalized cross validation (GCV) measure of distance between the pre- and post-filtered solutions. Efficient algorithms are developed to compute both the scalar and vector penalty parameters. The B-spline filter has a sharper localization to high-wavenumber than compact or explicit filters of the same stencil width and is demonstrated for solutions of the Burgers' equation, decaying Burgers' turbulence, and compressible Navier–Stokes turbulent channel flow. Furthermore, these simulations confirm the scheme's numerical stability and ability to narrowly target the high wavenumber components of numerical solutions. An advantage over finite-difference filters is that these B-spline filters are stable on bounded domains and even preserve formal order of accuracy.
This paper studies a power transmission system with both conventional generators (CGs) and distributed energy assets (DEAs) providing frequency control. We consider an operating condition with demand aggregating two dynamic components: one that switches between different values on a finite set, and one that varies smoothly over time. Such dynamic operating conditions may result from protection scheme activations, external cyber-attacks, or due to the integration of dynamic loads, such as data centers. Mathematically, the dynamics of the resulting system are captured by a system that switches between a finite number of vector fields - or modes -, with each mode having a distinct equilibrium point induced by the demand aggregation. To analyze the stability properties of the resulting switching system, we leverage tools from hybrid dynamic inclusions and the concept of ..omega.. -limit sets from sets. Specifically, we characterize a compact set that is semi-globally practically asymptotically stable under the assumption that the switching frequency and load variation rate are sufficiently slow. For arbitrarily fast variations of the load, we use a level-set argument with multiple Lyapunov functions to establish input-to-state stability of a larger set and with respect to the rate of change of the loads. The theoretical results are illustrated via numerical simulations on the IEEE 39-bus test system.
Vibration of a statically loaded, inherently compensated hydrostatic journal bearing due to oscillating exhaust pressure is investigated. Both angular and radial vibration modes are analyzed. The time-dependent Reynolds equation governing the pressure distribution between the oscillating journal and sleeve is solved together with the journal equation of motion to obtain the response characteristics of the bearing. The Reynolds equation and the equation of motion are simplified by applying regular perturbation theory for small displacements. The numerical solutions of the perturbation equations are obtained by discretizing the pressure field using finite-difference aproximations with a discrete, nonuniform line-source model which excludes effects due to feeding hole volume. An iterative scheme is used to simultaneously satisfy the equations of motion for the journal. The results presented include Bode plots of bearing-oscillation gain and phase for a particular bearing configuration for various combinations of parameters over a range of frequencies, including the resonant frequency.