SEARCH · Search NASA

Results for “nonlinear conservation laws”

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 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS1983

Application of second-order-accurate Total Variation Diminishing (TVD) schemes to the Euler equations in general geometries

A one-parameter family of explicit and implicit second-order-accurate, entropy satisfying, total variation diminishing (TVD) schemes was developed by Harten. These TVD schemes were the property of not generating spurious oscillations for one-dimensional nonlinear scalar hyperbolic conservation laws and constant coefficient hyperbolic systems. Application of these methods to one- and two-dimensional fluid flows containing shocks (in Cartesian coordinates) yields highly accurate nonoscillatory numerical solutions. The goal of this work is to expand these methods to the multidimensional Euler equations in generalized coordinate systems. Some numerical results of shock waves impinging on cylindrical bodies are compared with MacCormack's method.

Yee, H. C.↗

NASA NTRS2001

Rheology of Foam Near the Order-Disorder Transition

The first part of our research results are summarized in the recent journal publication: J. Gregory McDaniel and R. Glynn Holt, 'Measurement of aqueous foam rheology by acoustic levitation', Phys. Rev. E 61, 2204 (2000). This aspect of the work was a combination of experiment and analysis. We built a levitation system capable of acoustically levitating small samples of aqueous foam of arbitrary gas and liquid volume fractions. We then modulated the acoustic field to induce normal mode oscillations of the foam samples. The observables from the experiment were frequency and mode number. For dry (roughly > 70% gas by volume) foams and small deformations, we developed an effective medium, normal-modes analysis which took the frequency and mode number from experiment, and gave us the shear elastic modulus of the foam as a function of Poisson's ratio. The second part of our results may be found in a soon-to-be submitted manuscript 'Dynamics of aqueous foam drops', I.Sh. Akhatov, J.G. McDaniel and R.G. Holt, describing our modeling in the wet foam limit by considering the acoustic problem. This aspect of the research is purely theoretical. Beginning from a mass-conserving mixture law, the fully nonlinear equations of motion for a wet (roughly < 10% gas by volume) foam drop of initially spherical shape were derived. The frequencies for normal mode oscillations were derived in the linear inviscid limit. The nonlinear equations were numerically solved to elicit the motion of a foam drop under acoustic excitation. The role of the time-varying void fraction in breathing-mode oscillations is of particular interest. As of the end of the current (NAG#3-2121) grant, this work was not yet concluded. We continue to work on this aspect in order to extend the analysis to cover the transition regime of gas volume fractions, as well as to compare to experiments in the wet regime.

Holt, R. Glynn↗

NASA NTRS1999

Entropy Splitting and Numerical Dissipation

A rigorous stability estimate for arbitrary order of accuracy of spatial central difference schemes for initial-boundary value problems of nonlinear symmetrizable systems of hyperbolic conservation laws was established recently by Olsson and Oliger (1994) and Olsson (1995) and was applied to the two-dimensional compressible Euler equations for a perfect gas by Gerritsen and Olsson (1996) and Gerritsen (1996). The basic building block in developing the stability estimate is a generalized energy approach based on a special splitting of the flux derivative via a convex entropy function and certain homogeneous properties. Due to some of the unique properties of the compressible Euler equations for a perfect gas, the splitting resulted in the sum of a conservative portion and a non-conservative portion of the flux derivative. hereafter referred to as the "Entropy Splitting." There are several potential desirable attributes and side benefits of the entropy splitting for the compressible Euler equations that were not fully explored in Gerritsen and Olsson. The paper has several objectives. The first is to investigate the choice of the arbitrary parameter that determines the amount of splitting and its dependence on the type of physics of current interest to computational fluid dynamics. The second is to investigate in what manner the splitting affects the nonlinear stability of the central schemes for long time integrations of unsteady flows such as in nonlinear aeroacoustics and turbulence dynamics. If numerical dissipation indeed is needed to stabilize the central scheme, can the splitting help minimize the numerical dissipation compared to its un-split cousin? Extensive numerical study on the vortex preservation capability of the splitting in conjunction with central schemes for long time integrations will be presented. The third is to study the effect of the non-conservative proportion of splitting in obtaining the correct shock location for high speed complex shock-turbulence interactions. The fourth is to determine if this method can be extended to other physical equations of state and other evolutionary equation sets. If numerical dissipation is needed, the Yee, Sandham, and Djomehri (1999) numerical dissipation is employed. The Yee et al. schemes fit in the Olsson and Oliger framework.

Yee, H. C.↗

NASA NTRS1982

High-resolution schemes for hyperbolic conservation laws

A class of new explicit second order accurate finite difference schemes for the computation of weak solutions of hyperbolic conservation laws is presented. These highly nonlinear schemes are obtained by applying a nonoscillatory first order accurae scheme to an appropriately modified flux function. The so derived second order accurate schemes achieve high resolution while preserving the robustness of the original nonoscillatory first order accurate scheme.

Harten, A.↗

NASA NTRS1983

High resolution schemes for hyperbolic conservation laws

A class of new explicit second order accurate finite difference schemes for the computation of weak solutions of hyperbolic conservation laws is presented. These highly nonlinear schemes are obtained by applying a nonoscillatory first order accurate scheme to an appropriately modified flux function. The so-derived second order accurate schemes achieve high resolution while preserving the robustness of the original nonoscillatory first order accurate scheme. Numerical experiments are presented to demonstrate the performance of these new schemes.

Harten, A.↗

NASA NTRS1988

Total Variation Diminishing (TVD) schemes of uniform accuracy

Explicit second-order accurate finite-difference schemes for the approximation of hyperbolic conservation laws are presented. These schemes are nonlinear even for the constant coefficient case. They are based on first-order upwind schemes. Their accuracy is enhanced by locally replacing the first-order one-sided differences with either second-order one-sided differences or central differences or a blend thereof. The appropriate local difference stencils are selected such that they give TVD schemes of uniform second-order accuracy in the scalar, or linear systems, case. Like conventional TVD schemes, the new schemes avoid a Gibbs phenomenon at discontinuities of the solution, but they do not switch back to first-order accuracy, in the sense of truncation error, at extrema of the solution. The performance of the new schemes is demonstrated in several numerical tests.

Hartwich, PETER-M.↗

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 NTRS1989

Nonlinear filters for efficient shock computation

A new type of methods for the numerical approximation of hyperbolic conservation laws with discontinuous solution is introduced. The methods are based on standard finite difference schemes. The difference solution is processed with a nonlinear conservation form filter at every time level to eliminate spurious oscillations near shocks. It is proved that the filter can control the total variation of the solution and also produce sharp discrete shocks. The method is simpler and faster than many other high resolution schemes for shock calculations. Numerical examples in one and two space dimensions are presented.

Engquist, Bjorn↗

NASA NTRS2000

Detection of Extraterrestrial Ecology (Exoecology)

Researchers in the Astrobiology Technology Branch at Ames Research Center have begun investigating alternate concepts for the detection of extraterrestrial life. We suggest searching for extraterrestrial ecology, exoecology, as well as for extraterrestrial biology, exobiology. Ecology describes the interactions of living things with their environment. All ecosystems are highly constrained by their environment and are constrained by well-known system design principles. Ecology could exist wherever there is an energy source and living I things have discovered some means to capture, store, and use the available energy. Terrestrial ecosystems use as energy sources, light, organic molecules, and in thermal vents and elsewhere, simple inorganic molecules. Ecosystem behavior is controlled by matter and energy conservation laws and can be described by linear and nonlinear dynamic systems theory. Typically in an ecosystem different molecules are not in chemical equilibrium and scarce material is conserved, stored, or recycled. Temporal cycles and spatial variations are often observed. These and other -eneral principles of exoecology can help guide the search for extraterrestrial life. The chemical structure observed in terrestrial biology may be highly contingent on evolutionary accidents. Oxygen was not always abundant on Earth. Primitive sulfur bacteria use hydrogen sulfide and sulfur to perform photosynthesis instead of water and oxygen. Astrobiologists have assumed, for the sake of narrowing and focusing our life detection strategies, that extraterrestrial life will have detailed chemical similarities with terrestrial life. Such assumptions appear very reasonable and they allow us to design specific and highly sensitive life detection experiments. But the fewer assumptions we make, the less chance we have of being entirely wrong The best strategy for the detection of extraterrestrial life could be a mixed strategy. We should use detailed assumptions based on terrestrial biology to guide some but not all future searches for alien life. The systems principles of exoecology seem much more fundamental and inescapable than the terrestrial biology analogies of exobiology. We should search for exoecology as well as exobiology.

Jones, Harry↗

NASA NTRS1987

A generalized procedure for constructing an upwind-based TVD scheme

A generalized formulation for constructing second- and higher-order accurate TVD (total variation diminishing) schemes is presented. A given scheme is made TVD by limiting antidiffusive flux differences with some nonlinear functions, so-called limiters. The general idea of the formulation and its mathematical proof of Harten's TVD conditions is shown by applying the Lax-Wendroff method to a scalar nonlinear equation and constant-coefficient system of conservation laws. For the system of equations, several definitions are derived for the argument used in the limiter function and present their performance to numerical experiments. Then the formulation is formally extended to the nonlinear system of equations. It is demonstrated that use of the present procedure allows easy conversion of existing central or upwind, and second- or higher-order differencing schemes so as to preserve monotonicity and to yield physically admissible solutions. The formulation is simple mathematically as well as numerically; neither matrix-vector multiplication nor Riemann solver is required. Roughly twice as much computational effort is needed as compared to conventional scheme. Although the notion of TVD is based on the initial value problem, application to the steady Euler equations of the formulation is also made. Numerical examples including various ranges of problems show both time- and spatial-accuracy in comparison with exact solutions.

Liou, Meng-Sing↗

NASA NTRS1995

Symmetries of the TDNLS equations for weakly nonlinear dispersive MHD waves

In this paper we consider the symmetries and conservation laws for the TDNLS equations derived by Hada (1993) and Brio, Hunter and Johnson, to describe the propagation of weakly nonlinear dispersive MHD waves in beta approximately 1 plasmas. The equations describe the interaction of the Alfven and magnetoacoustic modes near the triple umbilic, where the fast magnetosonic, slow magnetosonic and Alfven speeds coincide and a(g)(exp 2) = V(A)(exp 2) where a(g) is the gas sound speed and V(A) is the Alfven speed. We discuss Lagrangian and Hamiltonian formulations, and similarity solutions for the equations.

Webb, G. M.↗

NASA NTRS1989

An Entropy-Based Approach to Nonlinear Stability

Many numerical methods used in computational fluid dynamics (CFD) incorporate an artificial dissipation term to suppress spurious oscillations and control nonlinear instabilities. The same effect can be accomplished by using upwind techniques, sometimes augmented with limiters to form Total Variation Diminishing (TVD) schemes. An analysis based on numerical satisfaction of the second law of thermodynamics allows many such methods to be compared and improved upon. A nonlinear stability proof is given for discrete scalar equations arising from a conservation law. Solutions to such equations are bounded in the L sub 2 norm if the second law of thermodynamics is satisfied in a global sense over a periodic domain. It is conjectured that an analogous statement is true for discrete equations arising from systems of conservation laws. Analysis and numerical experiments suggest that a more restrictive condition, a positive entropy production rate in each cell, is sufficient to exclude unphysical phenomena such as oscillations and expansion shocks. Construction of schemes which satisfy this condition is demonstrated for linear and nonlinear wave equations and for the one-dimensional Euler equations.

Merriam, Marshal L.↗

NASA NTRS1987

A generalized procedure for constructing an upwind based TVD scheme

A generalized formulation for constructing second- and higher-order accurate TVD (total variation diminishing) schemes is presented. A given scheme is made TVD by limiting antidiffusive flux differences with some linear functions, so-called limiters. The general idea of the formulation and its mathematical proof of Harten's TVD conditions is shown by applying the Lax-Wendroff method to scalar nonlinear equations and a constant-coefficient system of conservation laws. For the system of equations, several definitions are derived for the argument used in the limiter function and present their performance in numerical experiments. The formulation is extended to the nonlinear system. It is demonstrated that the present procedure can easily convert existing central or upwind, and second- or higher-order differencing schemes to preserve monotonicity and yield physically admissible solutions. The formulation is simple mathematically as well as numerically; both matrix-vector multiplication and Riemann solver are avoided. Although the notion of TVD is based on the initial value problem, application to the steady Euler equations of the formulation is also made.

Liou, Meng-Sing↗

NASA NTRS1995

Evolution of Vapor Bubbles Nucleation Sites in Low Gravity

When liquid is expelled by a vapor bubble growing at a nucleation site on a superheated surface, a thin microlayer underneath the bubble is left behind. It is evaporated from the free microlayer surface that provides for bubble growth. The average thickness of the microlayer determining the evaporation rate increases with time if the latter does not exceed a threshold value associated with the burn-out crisis. The bubble is described as a spherical segment with its flattened part adjoining the microlayer. This introduces two independent variables - the radius of the spherical part of the bubble surface and the polar angle that defines the relative area of the flattened part. They are to be found out from a set of two strongly nonlinear equations resulting from mass and momentum conservation laws. The first one depends on both microlayer thickness and nonmonotonously changing bubble base area. The second involves two major factors favoring bubble detachment - the buoyancy and a force due to the initial momentum of vapor input into the bubble. The former force depends on gravity whereas the latter one does not. It is why the limiting regimes of bubble evolution that correspond to normal or moderately reduced gravity and to microgravity feature drastically different properties. In the first case, the buoyancy dominates and the bubble evolves in such a manner as to become a full sphere at a moment that can be viewed as that of detachment. The detachment volume grows as gravity decreases. In the second case, the buoyancy is negligible and the bubble stays near the surface, while its volume continues to increase for a sufficiently long time. The findings are discussed in connection with experimental data obtained under different gravity conditions, some unpublished experiments being included. They help to understand why the pool boiling heat transfer coefficient frequently increases as gravity falls down and eventually vanishes.

Buyevich, Yu A.↗

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↗

NASA NTRS2003

New Type of the Interface Evolution in the Richtmyer-Meshkov Instability

We performed systematic theoretical and numerical studies of the nonlinear large-scale coherent dynamics in the Richtmyer-Meshkov instability for fluids with contrast densities. Our simulations modeled the interface dynamics for compressible and viscous uids. For a two-fluid system we observed that in the nonlinear regime of the instability the bubble velocity decays and its surface attens, and the attening is accompanied by slight oscillations. We found the theoretical solution for the system of conservation laws, describing the principal influence of the density ratio on the motion of the nonlinear bubble. The solution has no adjustable parameters, and shows that the attening of the bubble front is a distinct property universal for all values of the density ratio. This property follows from the fact that the RM bubbles decelerate. The theoretical and numerical results validate each other, describe the new type of the bubble front evolution in RMI, and identify the bubble curvature as important and sensitive diagnostic parameter.

Abarzhi, S. I.↗

NASA NTRS2023

Chemical Thermodynamics and the Mathematical Integration of Reaction Kinetics

Key advances in the development of numerical methods for non-reacting compressible flows have been enabled by translating physical requirements into concrete numerical guidelines, such as the satisfaction of entropy inequalities for shock-capturing techniques [Lax, Contributions to Nonlinear Functional Analysis (1971) 603-634]. In the present work, we present nonlinear numerical analysis tools that draw from Chemical Thermodynamics , the branch of Nonequilibrium Thermodynamics that deals with chemical reactions. Through Gibbs formalism, chemical thermodynamics provides a well-known theoretical expression for the chemical equilibrium constant of a reaction in terms of reduced chemical potentials. A less-known, yet extremely valuable result, due to [Krambeck, Arch. Ration. Mech. Anal. , 38 (1970) 317], states that when this expression is implemented, mass-action kinetic models are consistent with the dynamical prescriptions of the 2nd law of thermodynamics. For fixed-temperature ordinary differential equations modeling constant-volume reacting gas mixtures, this leads to a decreasing Helmholtz free energy. If the temperature is allowed to vary in accordance with conservation of energy (1st law), this leads to the statement of increasing entropy. These nonlinear prescriptions can, and should be, used to further develop temporal integration techniques for reaction kinetics. We demonstrate that Krambeck's result holds even when the equilibrium constants are approximated from data. We prove this result by constructing the implicit free energy and the implicit entropy inherent to a given approximation. This is first done for a 5-species, 17-reaction model problem for air. With this structure established, elements of discrete entropy-stability theory [Tadmor, Acta Numer. , 12 (2003) 451] are leveraged to examine the consistency of time-integration schemes with these prescriptions. Using chemical potentials, one can compute the respective contributions of the kinetics model and of the temporal scheme to free energy/entropy variations. We introduce a nonlinear-stable version of the Discontinuous-Galerkin (DG) scheme in time which shows robustness improvements over the standard linearly-stable version. Most notably, the maximum timestep that can be resolved with the nonlinearly-stable variant tends to grow with polynomial order, in contrast to the linearly-stable variant. We generalize our constructions to arbitrary systems of reversible chemical reactions, ultimately showing that the compressible reacting Euler system admits the opposite of the implicitly constructed thermodynamic entropy as a mathematical entropy . This lays important theoretical foundations towards robust scheme development [Harten, J. Comput. Phys. 49 (1983) 151-164].

STMD↗

NASA NTRS1983

Implicit Total Variation Diminishing (TVD) schemes for steady-state calculations

The application of a new implicit unconditionally stable high resolution total variation diminishing (TVD) scheme to steady state calculations. It is a member of a one parameter family of explicit and implicit second order accurate schemes developed by Harten for the computation of weak solutions of hyperbolic conservation laws. This scheme is guaranteed not to generate spurious oscillations for a nonlinear scalar equation and a constant coefficient system. Numerical experiments show that this scheme not only has a rapid convergence rate, but also generates a highly resolved approximation to the steady state solution. A detailed implementation of the implicit scheme for the one and two dimensional compressible inviscid equations of gas dynamics is presented. Some numerical computations of one and two dimensional fluid flows containing shocks demonstrate the efficiency and accuracy of this new scheme.

Yee, H. C.↗