Search NASA⌕ Search

SEARCH · Search NASA

Results for “Hyperbolic systems”

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 145 records · Page 8

A New Time-Space Accurate Scheme for Hyperbolic Problems: Quasi-Explicit Case - 1

This paper presents a new discretization scheme for hyperbolic systems of conservations laws. It satisfies the TVD property and relies on the new high-resolution mechanism which is compatible with the genuinely multidimensional approach proposed recently. This work can be regarded as a first step towards extending the genuinely multidimensional approach to unsteady problems. Discontinuity capturing capabilities and accuracy of the scheme are verified by a set of numerical tests.

Sidilkover, David↗

Gas-Kinetic Theory Based Flux Splitting Method for Ideal Magnetohydrodynamics

A gas-kinetic solver is developed for the ideal magnetohydrodynamics (MHD) equations. The new scheme is based on the direct splitting of the flux function of the MHD equations with the inclusion of "particle" collisions in the transport process. Consequently, the artificial dissipation in the new scheme is much reduced in comparison with the MHD Flux Vector Splitting Scheme. At the same time, the new scheme is compared with the well-developed Roe-type MHD solver. It is concluded that the kinetic MHD scheme is more robust and efficient than the Roe- type method, and the accuracy is competitive. In this paper the general principle of splitting the macroscopic flux function based on the gas-kinetic theory is presented. The flux construction strategy may shed some light on the possible modification of AUSM- and CUSP-type schemes for the compressible Euler equations, as well as to the development of new schemes for a non-strictly hyperbolic system.

Xu, Kun↗

High Order Numerical Simulation of Sound Generated by the Kirchhoff Vortex

An improved high order finite difference method for low Mach number computational aeroacoustics (CAA) is described. The improvements involve the conditioning of the Euler equations in perturbation form to minimize numerical cancellation error, and the use of a stable non-dissipative sixth-order central spatial differencing for the interior points and third-order at the boundary points. The spatial difference operator satisfies the summation-by-parts property to guarantee strict stability for linear hyperbolic systems. Spurious high frequency oscillations are damped by a third-order characteristic-based filter. The objective of this paper is to apply these improvements in the simulation of sound generated by the Kirchhoff vortex.

Mueller, Bernhard↗

Macroscopic Models of Radiative Transfer as Applied to Computation of the Radiation Field in the Solar Atmosphere

In this paper a new three dimensional half-moment model for radiative transfer is presented for a gray medium. It describes the evolution of the zeroth and first directional half moments of the radiative intensity. The closure is provided, similarly to Dubroca and Klar, by the maximum entropy concept. This work generalizes that model to three dimensions. The model presented here (the derivation being done in Ripoll and Wray, called the M(sup 1/2)(sub 1) model, is a hyperbolic system consisting of a total of eight equations in three dimensions, four equations for each direction. Each half model has the classical form of a macroscopic moment model in which the pressure tensor is constructed from the well-known Eddington tensor with a particular Eddington factor. Moreover, different source and border terms occur. The latter introduce couplings between the macroscopic and microscopic quantities and between the + and - streams, through the intensity in the plane perpendicular to the flux. The main theoretical application of the half moment model, treated in this paper, is its reduction to a full moment model, called M(sup +)(sub 1), for the particular but important case of a hot, opaque source radiating in a cold transparent (or semi-transparent) medium for very specific applications, such as stellar interiors or atmospheres, or combustion problems. The structure of the paper is as follows. In section 2, the model M(sup 1/2)(sub 1) is presented. In section 3, for the particular case of a hot, opaque source radiating into a cold medium, the half moment model is reduced to the M(sup +)(sub 1) model. In section 4, we first solve a simple and academic problem to validate the models, followed by a simplified solar atmosphere.

Ripoll, J.-F.↗

A Study of Multigrid Preconditioners Using Eigensystem Analysis

The convergence properties of numerical schemes for partial differential equations are studied by examining the eigensystem of the discrete operator. This method of analysis is very general, and allows the effects of boundary conditions and grid nonuniformities to be examined directly. Algorithms for the Laplace equation and a two equation model hyperbolic system are examined.

Roberts, Thomas W.↗

Evolution of Advection Upstream Splitting Method Schemes

This paper focuses on the evolution of advection upstream splitting method(AUSM) schemes. The main ingredients that have led to the development of modern computational fluid dynamics (CFD) methods have been reviewed, thus the ideas behind AUSM. First and foremost is the concept of upwinding. Second, the use of Riemann problem in constructing the numerical flux in the finite-volume setting. Third, the necessity of including all physical processes, as characterised by the linear (convection) and nonlinear (acoustic) fields. Fourth, the realisation of separating the flux into convection and pressure fluxes. The rest of this review briefly outlines the technical evolution of AUSM and more details can be found in the cited references. Keywords: Computational fluid dynamics methods, hyperbolic systems, advection upstream splitting method, conservation laws, upwinding, CFD

Liou, Meng-Sing↗

The High-Resolution Wave-Propagation Method Applied to Meso- and Micro-Scale Flows

The high-resolution wave-propagation method for computing the nonhydrostatic atmospheric flows on meso- and micro-scales is described. The design and implementation of the Riemann solver used for computing the Godunov fluxes is discussed in detail. The method uses a flux-based wave decomposition in which the flux differences are written directly as the linear combination of the right eigenvectors of the hyperbolic system. The two advantages of the technique are: 1) the need for an explicit definition of the Roe matrix is eliminated and, 2) the inclusion of source term due to gravity does not result in discretization errors. The resulting flow solver is conservative and able to resolve regions of large gradients without introducing dispersion errors. The methodology is validated against exact analytical solutions and benchmark cases for non-hydrostatic atmospheric flows.

Ahmad, Nashat N.↗

High-resolution Wave Propagation Method for Stratified Flows

The implementation of the multidimensional f-waves Riemann solver for the time-dependent, three-dimensional, nonhydrostatic, meso- and microscale atmospheric flows is described in detail. The Riemann solver employs flux-based wave decomposition (f-waves) for the calculation of Godunov fluxes in which the flux differences are written directly as the linear combination of the right eigenvectors of the hyperbolic system. The scheme incorporates the source term due to gravity without introducing discretization errors which is an important property in the context of atmospheric flows. The resulting flow solver is conservative, accurate, stable, and well-balanced. The implementation of the solver is evaluated using benchmark test cases for atmospheric dynamics.

Riemann problem↗

TPSAS-NF1676L-10658-DND

The systematic methodology for constructing fourth-order finite domain Energy Stable WENO schemes is developed. We prove that for hyperbolic systems, the finite domain ESWENO scheme is stable in the energy norm for both continuous and discontinuous solutions. The eigenvalues of the finite domain ESWENO dissipation operator are located in the left-half plane. Based on the rigorous truncation error analysis, the new weight functions are developed, which drastically improve the accuracy of the ESWENO scheme and provide excellent shock-capturing capabilities. Numerical experiments show that the new finite domain ESWENO scheme with the new weights outperform the conventional WENO schemes in terms of accuracy.

Travis Fisher↗

On the identification of continuous vibrating systems modelled by hyperbolic partial differential equations

This paper deals with the identification of spatially varying parameters in systems of finite spatial extent which can be described by second order hyperbolic differential equations. Two questions have been addressed. The first deals with 'partial identification' and inquires into the possibility of retrieving all the eigenvalues of the system from response data obtained at one location x-asterisk epsilon (0, 1). The second deals with the identification of the distributed coefficients rho(x), a(x) and b(x). Sufficient conditions for unique identification of all the eigenvalues of the system are obtained, and conditions under which the coefficients can be uniquely identified using suitable response data obtained at one point in the spatial domain are determined. Application of the results and their usefulness is demonstrated in the identification of the properties of tall building structural systems subjected to dynamic load environments.

Udwadia, F. E.↗

Guidance of Nonlinear Nonminimum-Phase Dynamic Systems

The research work has advanced the inversion-based guidance theory for: systems with non-hyperbolic internal dynamics; systems with parameter jumps; and systems where a redesign of the output trajectory is desired. A technique to achieve output tracking for nonminimum phase linear systems with non-hyperbolic and near non-hyperbolic internal dynamics was developed. This approach integrated stable inversion techniques, that achieve exact-tracking, with approximation techniques, that modify the internal dynamics to achieve desirable performance. Such modification of the internal dynamics was used (a) to remove non-hyperbolicity which is an obstruction to applying stable inversion techniques and (b) to reduce large preactuation times needed to apply stable inversion for near non-hyperbolic cases. The method was applied to an example helicopter hover control problem with near non-hyperbolic internal dynamics for illustrating the trade-off between exact tracking and reduction of preactuation time. Future work will extend these results to guidance of nonlinear non-hyperbolic systems. The exact output tracking problem for systems with parameter jumps was considered. Necessary and sufficient conditions were derived for the elimination of switching-introduced output transient. While previous works had studied this problem by developing a regulator that maintains exact tracking through parameter jumps (switches), such techniques are, however, only applicable to minimum-phase systems. In contrast, our approach is also applicable to nonminimum-phase systems and leads to bounded but possibly non-causal solutions. In addition, for the case when the reference trajectories are generated by an exosystem, we developed an exact-tracking controller which could be written in a feedback form. As in standard regulator theory, we also obtained a linear map from the states of the exosystem to the desired system state, which was defined via a matrix differential equation.

Devasia, Santosh↗

Output Tracking for Systems with Non-Hyperbolic and Near Non-Hyperbolic Internal Dynamics: Helicopter Hover Control

A technique to achieve output tracking for nonminimum phase linear systems with non-hyperbolic and near non-hyperbolic internal dynamics is presented. This approach integrates stable inversion techniques, that achieve exact-tracking, with approximation techniques, that modify the internal dynamics to achieve desirable performance. Such modification of the internal dynamics is used (1) to remove non-hyperbolicity which an obstruction to applying stable inversion techniques and (2) to reduce large pre-actuation time needed to apply stable inversion for near non-hyperbolic cases. The method is applied to an example helicopter hover control problem with near non-hyperbolic internal dynamic for illustrating the trade-off between exact tracking and reduction of pre-actuation time.

Devasia, Santosh↗

On the symmetric form of systems of conservation laws with entropy

The present investigation has the objective to review the general structure of system of conservation laws with entropy, giving particular attention to a presentation of symmetric formulations of the equations of gas dynamics. Systems of hyperbolic conservation laws are considered, taking into account the case in which a system of equations is symmetric hyperbolic, questions of symmetrization, and entropy functions. The Euler equations for a polytropic gas in conservation form are considered. The symmetric form retains the conservation properties of the equations. Weak solutions remain, therefore, unchanged. Of particular computational interest is the possibility of using the symmetric form.

Harten, A.↗

Finite difference numerical methods for boundary control problems governed by hyperbolic partial differential equations

This paper briefly reviews convergent finite difference schemes for hyperbolic initial boundary value problems and their applications to boundary control systems of hyperbolic type which arise in the modelling of vibrations. These difference schemes are combined with the primal and the dual approaches to compute the optimal control in the unconstrained case, as well as the case when the control is subject to inequality constraints. Some of the preliminary numerical results are also presented.

Chen, G.↗