SEARCH · Engineering Papers

Results for “conservative schemes”

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 379 records · Page 21

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Evolution of solitary density waves in stellar winds of early-type stars: A simple explanation of discrete absorption component behavior

We model the evolution of a density shell propagating through the stellar wind of an early-type star, in order to investigate the effects of such shells on UV P Cygni line profiles. Unlike previous treatments, we solve the mass, momentum, and energy conservation equations, using an explicit time-differencing scheme, and present a parametric study of the density, velocity, and temperature response. Under the assumed conditions, relatively large spatial scale, large-amplitude density shells propagate as stable waves through the supersonic portion of the wind. Their dynamical behavior appears to mimic propagating 'solitary waves,' and they are found to accelerate at the same rate as the underlying steady state stellar wind (i.e., the shell rides the wind). These hydrodynamically stable structures quantitatively reproduce the anomalous 'discrete absorption component' (DAC) behavior observed in the winds of luminous early-type stars, as illustrated by comparisons of model predictions to an extensive International Ultraviolet Explorer (IUE) time series of spectra of zeta Puppis (O4f). From these comparisons, we find no conclusive evidence indicative of DACs accelerating at a significantly slower rate than the underlying stellar wind, contrary to earlier reports. In addition, these density shells are found to be consistent within the constraints set by the IR observations. We conclude that the concept of propagating density shells should be seriously reconsidered as a possible explanation of the DAC phenomenon in early-type stars.

Waldron, Wayne L.↗

NASA NTRS1998

Annual Research Briefs, 1998

The topics contained in this progress report are direct numerical simulation of turbulent non-premixed combustion with realistic chemistry; LES of non-premixed turbulent reacting flows with conditional source term estimation; measurements of the three-dimensional scalar dissipation rate in gas-phase planar turbulent jets; direct simulation of a jet diffusion flame; on the use of interpolating wavelets in the direct numerical simulation of combustion; on the use of a dynamically adaptive wavelet collocation algorithm in DNS (direct numerical simulation) of non-premixed turbulent combustion; 2D simulations of Hall thrusters; computation of trailing-edge noise at low mach number using LES and acoustic analogy; weakly nonlinear modeling of the early stages of bypass transition; interactions between freestream turbulence and boundary layers; interfaces at the outer boundaries of turbulent motions; largest scales of turbulent wall flows; the instability of streaks in near-wall turbulence; an implementation of the v(sup 2) - f model with application to transonic flows; heat transfer predictions in cavities; a structure-based model with stropholysis effects; modeling a confined swirling coaxial jet; subgrid-scale models based on incremental unknowns for large eddy simulations; subgrid scale modeling taking the numerical error into consideration; towards a near-wall model for LES of a separated diffuser flow; on the feasibility of merging LES with RANS (Reynolds Averaging Numerical simulation) for the near-wall region of attached turbulent flows; large-eddy simulation of a separated boundary layer; numerical study of a channel flow with variable properties; on the construction of high order finite difference schemes on non-uniform meshes with good conservation properties; development of immersed boundary methods for complex geometries; and particle methods for micro and macroscale flow simulations.

Spinks, Debra↗

NASA NTRS1991

Cost-effective accurate coarse-grid method for highly convective multidimensional unsteady flows

A fundamentally multidimensional convection scheme is described based on vector transient interpolation modeling rewritten in conservative control-volume form. Vector third-order upwinding is used as the basis of the algorithm; this automatically introduces important cross-difference terms that are absent from schemes using component-wise one-dimensional formulas. Third-order phase accuracy is good; this is important for coarse-grid large-eddy or full simulation. Potential overshoots or undershoots are avoided by using a recently developed universal limiter. Higher order accuracy is obtained locally, where needed, by the cost-effective strategy of adaptive stencil expansion in a direction normal to each control-volume face; this is controlled by monitoring the absolute normal gradient and curvature across the face. Higher (than third) order cross-terms do not appear to be needed. Since the wider stencil is used only in isolated narrow regions (near discontinuities), extremely high (in this case, seventh) order accuracy can be achieved for little more than the cost of a globally third-order scheme.

Leonard, B. P.↗

NASA NTRS2000

Solution of Third Computational Workshop Internal Propagation Problems using Low Order Schemes

The problems solved are propagation of sound waves through a transonic nozzle and shock-sound interaction. The first problem is solved using a MacCormack scheme. The second problem is solved by perturbation of a conservative Euler equation solution to the steady state. The problems are solved on evenly spaced grids. While solutions were found, the methods selected and not using stretched grids lead to solutions which do not compare well with those found using more accurate schemes. However, the perturbation of the mean flow scheme used to solve the second problem shows promise.

Miles, Jeffrey Hilton↗

DOE OSTI2021

Multilevel Spectral Coarsening for Graph Laplacian Problems with Application to Reservoir Simulation

We extend previously developed two-level coarsening procedures for graph Laplacian problems written in a mixed saddle point form to the fully recursive multilevel case. The resulting hierarchy of discretizations gives rise to a hierarchy of upscaled models, in the sense that they provide approximation in the natural norms (in the mixed setting). This property enables us to utilize them in three applications: (i) as an accurate reduced model, (ii) as a tool in multilevel Monte Carlo simulations (in application to finite volume discretizations), and (iii) for providing a sequence of nonlinear operators in a full approximation scheme for solving nonlinear pressure equations discretized by the conservative two-point flux approximation. Finally, we illustrate the potential of the proposed multilevel technique in all three applications on a number of popular benchmark problems used in reservoir simulation.

multilevel Monte Carlo↗

NASA NTRS1982

Computations of transonic flow over an oscillating airfoil with shock-induced separation

A numerical study is presented on unsteady transonic flow over an airfoil in pitching oscillation at a moderately high incidence. The study was focused on the effects of unsteady shock-wave/boundary-layer interactions and subsequent flow separations. An implicit finite-difference scheme was used to solve the full time-dependent compressible Navier-Stokes equations in conservation-law form. The algorithm adopted for this study is the Beam-Warming spatially factored scheme. It is, however, constructed in a body-fitted curvilinear coordinate system for efficient unsteady flow computations. The computational results are compared with experimental data measured in the Ames 11- by 11-Foot Transonic Wind Tunnel.

Chyu, W. J.↗

DOE OSTI2022

Energy conserving and well-balanced discontinuous Galerkin methods for the Euler–Poisson equations in spherical symmetry

ABSTRACT This paper presents high-order Runge–Kutta (RK) discontinuous Galerkin methods for the Euler–Poisson equations in spherical symmetry. The scheme can preserve a general polytropic equilibrium state and achieve total energy conservation up to machine precision with carefully designed spatial and temporal discretizations. To achieve the well-balanced property, the numerical solutions are decomposed into equilibrium and fluctuation components that are treated differently in the source term approximation. One non-trivial challenge encountered in the procedure is the complexity of the equilibrium state, which is governed by the Lane–Emden equation. For total energy conservation, we present second- and third-order RK time discretization, where different source term approximations are introduced in each stage of the RK method to ensure the conservation of total energy. A carefully designed slope limiter for spherical symmetry is also introduced to eliminate oscillations near discontinuities while maintaining the well-balanced and total-energy-conserving properties. Extensive numerical examples – including a toy model of stellar core collapse with a phenomenological equation of state that results in core bounce and shock formation – are provided to demonstrate the desired properties of the proposed methods, including the well-balanced property, high-order accuracy, shock-capturing capability, and total energy conservation.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1984

A conservative treatment of zonal boundaries for Euler equation calculations

Finite-difference calculations require the generation of a grid for the region of interest. A zonal approach, wherein the given region is subdivided into zones and the grid for each zone is generated independently, makes the grid-generation process for complicated topologies and for regions requiring selective grid refinement a fairly simple task. This approach results in new boundaries within the given region, that is, zonal boundaries at the interfaces of the various zones. The zonal-boundary scheme (the integration scheme used to update the points on the zonal boundary) for the Euler equations must be conservative, accurate, stable, and applicable to general curvilinear coordinate systems. A zonal-boundary scheme with these desirable properties is developed in this study. The scheme is designed for explicit, first-order-accurate integration schemes but can be modified to accommodate second-order-accurate explicit and implicit integration schemes. Results for inviscid flow, including supersonic flow over a cylinder, blast-wave diffraction by a ramp, and one-dimensional shock-tube flow are obtained on zonal grids. The conservative nature of the zonal-boundary scheme permits the smooth transition of the discontinuities associated with these flows from one zone to another. The calculations also demonstrate the continuity of contour lines across zonal boundaries that can be achieved with the present zonal scheme.

Rai, M. M.↗

NASA NTRS1986

A Conservation Treatment of Zonal Boundaries for Euler Equation Calculations

Rai, Man Mohan↗

NASA NTRS1995

Nonlinear (time domain) and linearized (time and frequency domain) solutions to the compressible Euler equations in conservation law form

Two linearized solvers (time and frequency domain) based on a high resolution numerical scheme are presented. The basic approach is to linearize the flux vector by expressing it as a sum of a mean and a perturbation. This allows the governing equations to be maintained in conservation law form. A key difference between the time and frequency domain computations is that the frequency domain computations require only one grid block irrespective of the interblade phase angle for which the flow is being computed. As a result of this and due to the fact that the governing equations for this case are steady, frequency domain computations are substantially faster than the corresponding time domain computations. The linearized equations are used to compute flows in turbomachinery blade rows (cascades) arising due to blade vibrations. Numerical solutions are compared to linear theory (where available) and to numerical solutions of the nonlinear Euler equations.

Sreenivas, Kidambi↗

NASA NTRS1977

An implicit factored scheme for the compressible Navier-Stokes equations

An implicit finite difference scheme is developed for the numerical solution of the compressible Navier-Stokes equations in conservation-law form. The algorithm is second-order-time accurate, noniterative, and spatially factored. In order to obtain an efficient factored algorithm, the spatial cross-derivatives are evaluated explicitly. However, the algorithm is unconditionally stable and, although a three-time-level scheme, requires only two-time-levels of data storage. The algorithm is constructed in a 'delta' form (i.e., increments of the conserved variables and fluxes) that provides a direct derivation of the scheme and leads to an efficient computational algorithm. In addition, the delta form has the advantageous property of a steady-state (if one exists) independent of the size of the time step. Numerical results are presented for a two-dimensional shock boundary-layer interaction problem.

Beam, R. M.↗

NASA NTRS1985

A fast, time-accurate unsteady full potential scheme

The unsteady form of the full potential equation is solved in conservation form by an implicit method based on approximate factorization. At each time level, internal Newton iterations are performed to achieve time accuracy and computational efficiency. A local time linearization procedure is introduced to provide a good initial guess for the Newton iteration. A novel flux-biasing technique is applied to generate proper forms of the artificial viscosity to treat hyperbolic regions with shocks and sonic lines present. The wake is properly modeled by accounting not only for jumps in phi, but also for jumps in higher derivatives of phi, obtained by imposing the density to be continuous across the wake. The far field is modeled using the Riemann invariants to simulate nonreflecting boundary conditions. The resulting unsteady method performs well which, even at low reduced frequency levels of 0.1 or less, requires fewer than 100 time steps per cycle at transonic Mach numbers. The code is fully vectorized for the CRAY-XMP and the VPS-32 computers.

Shankar, V.↗

NASA NTRS1978

The artificial compression method for computation of shocks and contact discontinuities. III - Self-adjusting hybrid schemes

This paper presents a new computational method for the calculation of discontinuous solutions of hyperbolic systems of conservation laws, which deal effectively with both shock and contact discontinuities. The method consists of two stages: in the first stage a standard finite-difference scheme is hybridized with a nonoscillatory first order accurate method to provide for the monotonic variation of the solution near discontinuities, and in the second stage artificial compression is applied to sharpen transitions at discontinuities. This modification of a standard finite-difference method results in a scheme which preserves the order of truncation error of the original method and yet yields a sharp and oscillation free transition for both shocks and contact discontinuities. The modification can be easily implemented in existing computer codes.

Harten, A.↗

NASA NTRS1989

Comparison between pressure gradient method and MAC method on high Re calculation

A cavity flow driven by shear and buoyancy forces is used as a test problem in the application of a nonstaggered pressure gradient (PG) method in solving the two-dimensional incompressible Navier-Stokes equations. Twelve finite differencing schemes are used to solve the cavity flow problem. The schemes consist of various combinations of grid arrangements, upwinding treatments, and conservativeness of convection terms. An artificial source term is introduced, and the solutions are compared with those obtained by the conventional marker-and-cell (MAC) method. The comparisons favor the PG method. Numerical results obtained by the twelve schemes are compared with exact solutions in order to assess the stability and accuracy of each scheme.

Tan, C.-H.↗

NASA NTRS1993

Effects of spatial order of accuracy on the computation of vortical flowfields

The effect of the order-of-accuracy, used for the spatial discretization, on the resolution of the leading edge vortices over sharp-edged delta wings is investigated. The flowfield is computed using a viscous/inviscid zonal approach. The viscous flow in the vicinity of the wing is computed using the conservative formulation of the compressible, thin-layer Navier-Stokes equations. The leeward-side vortical flowfield and the other flow regions away from the surface are computed as inviscid. The time integration is performed with both an explicit fourth-order Runge-Kutta scheme and an implicit, factorized, iterative scheme. High-order-accurate inviscid fluxes are computed using both a conservative and a non-conservative (primitive variable) formulation. The nonlinear, inviscid terms of the primitive variable form of the governing equations are evaluated with a finite-difference numerical scheme based on the sign of the eigenvalues. High-order, upwind-biased, finite difference formulas are used to evaluate the derivatives of the nonlinear convective terms. Computed results are compared with available experimental data, and comparisons of the flowfield in the vicinity of the vortex cores are presented.

Ekaterinaris, J. A.↗

DOE OSTI2024

Asymptotic-preserving gyrokinetic implicit particle-orbit integrator for arbitrary electromagnetic fields

We extend the asymptotic preserving and energy conserving time integrator for charged-particle motion developed in Ricketson and Chacón (2020) to include finite Larmor-radius (FLR) effects in the presence of electric-field length-scales comparable to the particle gyro-radius (the gyro-kinetic limit). We introduce two modifications to the earlier scheme. The first is the explicit gyro-averaging of the electric field at the half time-step, along with an analogous modification to the current deposition, which we show preserves total energy conservation in implicit PIC schemes. The number of gyrophase samples is chosen adaptively, ensuring proper averaging for large timesteps and the recovery of full-orbit dynamics in the small time-step limit. The second modification is an alternating large and small time-step strategy that ensures the particle trajectory samples gyrophases evenly. We show that this strategy relaxes the time-step restrictions on the scheme, allowing even larger speed-ups than previously achievable. We demonstrate the new method with several single-particle motion tests in a variety of electromagnetic field configurations featuring gyro-scale variation in the electric field. Finally, the results demonstrate the advertised ability to capture FLR effects accurately even when significantly stepping over the gyration time-scale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1993

Positivity-preserving numerical schemes for multidimensional advection

This report describes the construction of an explicit, single time-step, conservative, finite-volume method for multidimensional advective flow, based on a uniformly third-order polynomial interpolation algorithm (UTOPIA). Particular attention is paid to the problem of flow-to-grid angle-dependent, anisotropic distortion typical of one-dimensional schemes used component-wise. The third-order multidimensional scheme automatically includes certain cross-difference terms that guarantee good isotropy (and stability). However, above first-order, polynomial-based advection schemes do not preserve positivity (the multidimensional analogue of monotonicity). For this reason, a multidimensional generalization of the first author's universal flux-limiter is sought. This is a very challenging problem. A simple flux-limiter can be found; but this introduces strong anisotropic distortion. A more sophisticated technique, limiting part of the flux and then restoring the isotropy-maintaining cross-terms afterwards, gives more satisfactory results. Test cases are confined to two dimensions; three-dimensional extensions are briefly discussed.

Leonard, B. P.↗

NASA NTRS1992

A hybrid multigrid technique for computing steady-state solutions to supersonic flows

Recently, Li and Sanders have introduced a class of finite difference schemes to approximate generally discontinuous solutions to hyperbolic systems of conservation laws. These equations have the form together with relevant boundary conditions. When modelling hypersonic spacecraft reentry, the differential equations above are frequently given by the compressible Euler equations coupled with a nonequilibrium chemistry model. For these applications, steady state solutions are often sought. Many tens (to hundreds) of super computer hours can be devoted to a single three space dimensional simulation. The primary difficulty is the inability to rapidly and reliably capture the steady state. In these notes, we demonstrate that a particular variant from the schemes presented can be combined with a particular multigrid approach to capture steady state solutions to the compressible Euler equations in one space dimension. We show that the rate of convergence to steady state coming from this multigrid implementation is vastly superior to the traditional approach of artificial time relaxation. Moreover, we demonstrate virtual grid independence. That is, the rate of convergence does not depend on the degree of spatial grid refinement.

Sanders, Richard↗