Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Discontinuous solutions”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

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↗
NASA NTRSDate not supplied

Discontinuous Galerkin and Related Methods for ODE

Starting from the standard integral formulation, the DG method is derived here in differential form. The key ingredient is a polynomial called the correction function, which helps ‘correct’ the discontinuous solution by approximating the jump and yields a continuous one. Under the right Radau quadrature, this continuous solution is identical to the solutions by the right Radau collocation and the continuous Galerkin (CG) methods. Next, the correction function facilitates the construction of the associated implicit Runge-Kutta schemes (IRK-DG). Different quadratures for DG result in different IRK-DG methods: left Radau quadrature in Radau IA, right Radau quadrature in Radau IIA or right Radau collocation, and Gauss quadrature in a method called DG-Gauss. The construction of IRK-DG clarifies the meaning and facilitates the proofs of various 𝐵(𝑝), 𝐶(𝜂), and 𝐷(𝜁) conditions for accuracy. The two consequences of these conditions are that all 𝑠-stage IRK-DG methods are accurate to order 2𝑠− 1, and the IRK-DG methods of Radau type are unique. Numerical examples showing the behavior of the DG solutions are provided. In all, the correction function plays a key role and helps establish the relations among the DG, IRK DG, collocation, and CG methods.

Numerical Methods for Ordinary Differential Equati↗
NASA NTRS1991

On the discontinuity in kinetic solutions of the collisionless polar wind

The use of a Maxwellian ion velocity distribution in kinetic solutions of the collisionless polar wind regime results in a velocity distribution and associated moments that display sharp discontinuities at the boundary, which quantitatively affects the solution at high altitudes. This paper investigates the causes of this discontinuity, using both analytical and numerical approaches; it is shown that the choice of appropriate boundary conditions will eliminate the discontinuity.

Khoyloo, A.↗
NASA NTRS1991

Finite element computation of compressible flows with the SUPG formulation

Finite element computation of compressible Euler equations is presented in the context of the streamline-upwind/Petrov-Galerkin (SUPG) formulation. The SUPG formulation, which is based on adding stabilizing terms to the Galerkin formulation, is further supplemented with a shock capturing operator which addresses the difficulty in maintaining a satisfactory solution near discontinuities in the solution field. The shock capturing operator, which has been derived from work done in entropy variables for a similar operator, is shown to lead to an appropriate level of additional stabilization near shocks, without resulting in excessive numerical diffusion. An implicit treatment of the impermeable wall boundary condition is also presented. This treatment of the no-penetration condition offers increased stability for large Courant numbers, and accelerated convergence of the computations for both implicit and explicit applications. Several examples are presented to demonstrate the ability of this method to solve the equations governing compressible fluid flow.

Le Beau, G. J.↗
NASA NTRS1988

Nonoscillatory solution of the steady-state inviscid Burgers' equation by mathematical programming

In order to obtain the physically relevant discontinuous numerical solution, the steady-state inviscid Burgers' equation is singularly perturbed through the addition of a small amount of viscosity. A 'cell-centered' finite-difference scheme is proposed which employs two points for the inviscid part and four points for the viscous one. While difficulties are experienced in the capture of interior layers centered at node points, computational results for interior layers centered between node points, and for boundary layers, exhibit accurate nonoscillatory solutions whose discontinuities are captured in one cell on both coarse and fine grids.

Lavery, John E.↗
NASA NTRS1983

Design and implementation of a multigrid code for the Euler equations

The steady-state equations of inviscid fluid flow, the Euler equations, are a nonlinear nonelliptic system of equations admitting solutions with discontinuities (for example, shocks). The efficient numerical solution of these equations poses a strenuous challenge to multigrid methods. A multigrid code has been developed for the numerical solution of the Euler equations. In this paper some of the factors that had to be taken into account in the design and development of the code are reviewed. These factors include the importance of choosing an appropriate difference scheme, the usefulness of local mode analysis as a design tool, and the crucial question of how to treat the nonlinearity. Sample calculations of transonic flow about airfoils will be presented. No claim is made that the particular algorithm presented is optimal.

Jespersen, D. C.↗
NASA NTRS1994

Fourth order difference methods for hyperbolic IBVP's

Fourth order difference approximations of initial-boundary value problems for hyperbolic partial differential equations are considered. We use the method of lines approach with both explicit and compact implicit difference operators in space. The explicit operator satisfies an energy estimate leading to strict stability. For the implicit operator we develop boundary conditions and give a complete proof of strong stability using the Laplace transform technique. We also present numerical experiments for the linear advection equation and Burgers' equation with discontinuities in the solution or in its derivative. The first equation is used for modeling contact discontinuities in fluid dynamics, the second one for modeling shocks and rarefaction waves. The time discretization is done with a third order Runge-Kutta TVD method. For solutions with discontinuities in the solution itself we add a filter based on second order viscosity. In case of the non-linear Burger's equation we use a flux splitting technique that results in an energy estimate for certain different approximations, in which case also an entropy condition is fulfilled. In particular we shall demonstrate that the unsplit conservative form produces a non-physical shock instead of the physically correct rarefaction wave. In the numerical experiments we compare our fourth order methods with a standard second order one and with a third order TVD-method. The results show that the fourth order methods are the only ones that give good results for all the considered test problems.

Gustafsson, Bertil↗
NASA NTRS1990

Formal mathematical solutions of the force-free equations, spontaneous discontinuities, and dissipation in large-scale magnetic fields

Direct integration of the force-free field equation del x B = alpha B, in the simple case of the local deformation of a laminar field, produces field configurations containing tangential discontinuities (current sheets). Whereas continuous solutions allow only restricted field topologies, the discontinuities provide the necessary release from those restrictions in more general topologies. Magnetic fields in nature are strongly deformed by convection, so as to contain significant internal discontinuities. The bipolar magnetic fields containing the active X-ray corona of the sun are a case in point. It appears that the dissipation caused by the discontinuities may be the primary heat source producing the X-ray corona.

Parker, E. N.↗
NASA NTRSDate not supplied

Discontinuous Galerkin and Related Methods for ODE

A defining feature of the discontinuous Galerkin (DG) method for ODE is that the piecewise polynomial solution can have a jump discontinuity at the beginning of each step. Starting from the standard integral formulation, the DG method is derived here in differential form. The key ingredient is a polynomial called the correction function, which helps ‘correct’ the discontinuous solution by approximating the jump and yields a continuous one. Under the right Radau quadrature, this continuous solution is identical to the solutions by the right Radau collocation and the continuous Galerkin (CG) methods. Next, the correction function facilitates the construction of the associated implicit Runge-Kutta schemes (IRK-DG). Different quadratures for DG result in different IRK-DG methods: left Radau quadrature in Radau IA, right Radau quadrature in Radau IIA or right Radau collocation, and Gauss quadrature in a method called DG-Gauss. The construction of IRK-DG clarifies the meaning and facilitates the proofs of various 𝐵(𝑝), 𝐶(𝜂), and 𝐷(𝜁) conditions for accuracy. The two consequences of these conditions are that all 𝑠-stage IRK-DG methods are accurate to order 2𝑠 − 1, and the IRK-DG methods of Radau type are unique. Numerical examples showing the behavior of the DG solutions are provided. In all, the correction function plays a key role and helps establish the relations among the DG, IRK-DG, collocation, and CG schemes.

numerical methods↗
NASA NTRS1976

Computer modeling of electromagnetic problems using the geometrical theory of diffraction

Some applications of the geometrical theory of diffraction (GTD), a high frequency ray optical solution to electromagnetic problems, are presented. GTD extends geometric optics, which does not take into account the diffractions occurring at edges, vertices, and various other discontinuities. Diffraction solutions, analysis of basic structures, construction of more complex structures, and coupling using GTD are discussed.

Burnside, W. D.↗
NASA NTRS2016

High-Order Residual-Distribution Schemes for Discontinuous Problems on Irregular Triangular Grids

In this paper, we develop second- and third-order non-oscillatory shock-capturing hyperbolic residual distribution schemes for irregular triangular grids, extending our second- and third-order schemes to discontinuous problems. We present extended first-order N- and Rusanov-scheme formulations for hyperbolic advection-diffusion system, and demonstrate that the hyperbolic diffusion term does not affect the solution of inviscid problems for vanishingly small viscous coefficient. We then propose second- and third-order blended hyperbolic residual-distribution schemes with the extended first-order Rusanov-scheme. We show that these proposed schemes are extremely accurate in predicting non-oscillatory solutions for discontinuous problems. We also propose a characteristics-based nonlinear wave sensor for accurately detecting shocks, compression, and expansion regions. Using this proposed sensor, we demonstrate that the developed hyperbolic blended schemes do not produce entropy-violating solutions (unphysical stocks). We then verify the design order of accuracy of these blended schemes on irregular triangular grids.

Mazaheri, Alireza↗
NASA NTRS1997

Essentially Non-Oscillatory and Weighted Essentially Non-Oscillatory Schemes for Hyperbolic Conservation Laws

In these lecture notes we describe the construction, analysis, and application of ENO (Essentially Non-Oscillatory) and WENO (Weighted Essentially Non-Oscillatory) schemes for hyperbolic conservation laws and related Hamilton- Jacobi equations. ENO and WENO schemes are high order accurate finite difference schemes designed for problems with piecewise smooth solutions containing discontinuities. The key idea lies at the approximation level, where a nonlinear adaptive procedure is used to automatically choose the locally smoothest stencil, hence avoiding crossing discontinuities in the interpolation procedure as much as possible. ENO and WENO schemes have been quite successful in applications, especially for problems containing both shocks and complicated smooth solution structures, such as compressible turbulence simulations and aeroacoustics. These lecture notes are basically self-contained. It is our hope that with these notes and with the help of the quoted references, the reader can understand the algorithms and code them up for applications.

Shu, Chi-Wang↗
NASA NTRS1995

Structure of the magnetopause inferred from one-dimensional hybrid simulations

The structure of the magnetopause is investigated by studying the interaction between two plasmas with solar wind and magnetospheric properties. Both Hall-MHD and hybrid (particle ions and fluid electrons) simulations are performed to compare and contrast the nature of the solutions in the fluid and the kinetic limits. It is shown that, in accordance with previous studies, the fluid solutions consist of multiple discontinuities and waves such as slow shocks and waves, as well as contact and rotational discontinuities. In contrast, the kinetic solutions consist of fewer discontinuities and include non-MHD boundaries. The difference between the two types of solutions are attributed to the absence of contact surfaces in collisionless plasmas and the possible Landau damping of slow waves. The kinetic solutions are found to be in a much better agreement with magnetopause observations, which have shown little evidence for the presence of slow shocks or contact discontinuities. The results of kinetic calculations suggest that the presence of a small but finite normal component of the magnetic field allows for the mixing of magnetosheath and magnetospheric plasmas and may, in part, be responsible for the formation of the boundary layer even during the periods of northward interplanetary magnetic field. These results also show that the necessary changes in the fields and the plasma parameters across the magnetopause do not occur simultaneously (i.e., in the same spatial location). As such, the magnetopause cannot be viewed as a single boundary but instead is a region which includes both abrupt and gradual changes in the fields and plasma parameters.

Omidi, N.↗
NASA NTRS1991

Recent developments in shock-capturing schemes

The development of the shock capturing methodology is reviewed, paying special attention to the increasing nonlinearity in its design and its relation to interpolation. It is well-known that higher-order approximations to a discontinuous function generate spurious oscillations near the discontinuity (Gibbs phenomenon). Unlike standard finite-difference methods which use a fixed stencil, modern shock capturing schemes use an adaptive stencil which is selected according to the local smoothness of the solution. Near discontinuities this technique automatically switches to one-sided approximations, thus avoiding the use of discontinuous data which brings about spurious oscillations.

Harten, Ami↗
NASA NTRS1973

High speed flow past wings

The analytical solution to the transonic small perturbation equation which describes steady compressible flow past finite wings at subsonic speeds can be expressed as a nonlinear integral equation with the perturbation velocity potential as the unknown function. This known formulation is substituted by a system of nonlinear algebraic equations to which various methods are applicable for its solution. Due to the presence of mathematical discontinuities in the flow solutions, however, a main computational difficulty was to ensure uniqueness of the solutions when local velocities on the wing exceeded the speed of sound. For continuous solutions this was achieved by embedding the algebraic system in an one-parameter operator homotopy in order to apply the method of parametric differentiation. The solution to the initial system of equations appears then as a solution to a Cauchy problem where the initial condition is related to the accompanying incompressible flow solution. In using this technique, however, a continuous dependence of the solution development on the initial data is lost when the solution reaches the minimum bifurcation point. A steepest descent iteration technique was therefore, added to the computational scheme for the calculation of discontinuous flow solutions. Results for purely subsonic flows and supersonic flows with and without compression shocks are given and compared with other available theoretical solutions.

Norstrud, H.↗
NASA NTRS2004

Discontinuous Spectral Difference Method for Conservation Laws on Unstructured Grids

A new, high-order, conservative, and efficient discontinuous spectral finite difference (SD) method for conservation laws on unstructured grids is developed. The concept of discontinuous and high-order local representations to achieve conservation and high accuracy is utilized in a manner similar to the Discontinuous Galerkin (DG) and the Spectral Volume (SV) methods, but while these methods are based on the integrated forms of the equations, the new method is based on the differential form to attain a simpler formulation and higher efficiency. Conventional unstructured finite-difference and finite-volume methods require data reconstruction based on the least-squares formulation using neighboring point or cell data. Since each unknown employs a different stencil, one must repeat the least-squares inversion for every point or cell at each time step, or to store the inversion coefficients. In a high-order, three-dimensional computation, the former would involve impractically large CPU time, while for the latter the memory requirement becomes prohibitive. In addition, the finite-difference method does not satisfy the integral conservation in general. By contrast, the DG and SV methods employ a local, universal reconstruction of a given order of accuracy in each cell in terms of internally defined conservative unknowns. Since the solution is discontinuous across cell boundaries, a Riemann solver is necessary to evaluate boundary flux terms and maintain conservation. In the DG method, a Galerkin finite-element method is employed to update the nodal unknowns within each cell. This requires the inversion of a mass matrix, and the use of quadratures of twice the order of accuracy of the reconstruction to evaluate the surface integrals and additional volume integrals for nonlinear flux functions. In the SV method, the integral conservation law is used to update volume averages over subcells defined by a geometrically similar partition of each grid cell. As the order of accuracy increases, the partitioning for 3D requires the introduction of a large number of parameters, whose optimization to achieve convergence becomes increasingly more difficult. Also, the number of interior facets required to subdivide non-planar faces, and the additional increase in the number of quadrature points for each facet, increases the computational cost greatly.

Liu, Yen↗