Uniformly high-order accurate non-oscillatory schemes, 1

The construction and the analysis of nonoscillatory shock capturing methods for the approximation of hyperbolic conservation laws was begun. These schemes share many desirable properties with total variation diminishing schemes (TVD), but TVD schemes have at most first order accuracy, in the sense of truncation error, at extreme of the solution. A uniformly second order approximation was constucted, which is nonoscillatory in the sense that the number of extrema of the discrete solution is not increasing in time. This is achieved via a nonoscillatory piecewise linear reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell

Some physical implications of recent solar wind measurements

The physical implications of the existence at about 1 AU of a quiet solar wind particle flux about 90 percent larger than that suggested in the past is investigated within the framework of the two-fluid solar wind model equations. During the spherically symmetric radial expansion of the quiet solar wind, the particle flux is conserved quantity. It is found that a pure collisional two-fluid model provides good particle density and streaming velocity at 1 AU, but predicts too large an electron temperature and too small a proton temperature. When noncollisional contributions to the transport coefficients are incorporated in the model equations, a complete satisfactory agreement with the available observations is obtained. Upper limits to the effective coupling between electrons and protons, as well as to the effective proton thermal conductivity, and both upper and lower limits to the effective electron thermal conductivity in the quiet solar wind, required to provide agreement with observations, are given

Uniformly high order accurate essentially non-oscillatory schemes 3

In this paper (a third in a series) the construction and the analysis of essentially non-oscillatory shock capturing methods for the approximation of hyperbolic conservation laws are presented. Also presented is a hierarchy of high order accurate schemes which generalizes Godunov's scheme and its second order accurate MUSCL extension to arbitrary order of accuracy. The design involves an essentially non-oscillatory piecewise polynomial reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell. The reconstruction algorithm is derived from a new interpolation technique that when applied to piecewise smooth data gives high-order accuracy whenever the function is smooth but avoids a Gibbs phenomenon at discontinuities. Unlike standard finite difference methods this procedure uses an adaptive stencil of grid points and consequently the resulting schemes are highly nonlinear

On the Nonlinearity of Modern Shock-Capturing Schemes

The development is reviewed of shock capturing methods, paying special attention to the increasing nonlinearity in the design of numerical schemes. The nature is studies of this nonlinearity and its relation to upwind differencing is examined. This nonlinearity of the modern shock capturing methods is essential, in the sense that linear analysis is not justified and may lead to wrong conclusions. Examples to demonstrate this point are given

A high order compact scheme for hypersonic aerothermodynamics

A novel high order compact scheme for solving the compressible Navier-Stokes equations has been developed. The scheme is an extension of a method originally proposed for solving the Euler equations, and combines several techniques for the solution of compressible flowfields, such as upwinding, limiting and flux vector splitting, with the excellent properties of high order compact schemes. Extending the method to the Navier-Stokes equations is achieved via a Kinetic Flux Vector Splitting technique, which represents an unusual and attractive way to include viscous effects. This approach offers a more accurate and less computationally expensive technique than discretizations based on more conventional operator splitting. The Euler solver has been validated against several inviscid test cases, and results for several viscous test cases are also presented. The results confirm that the method is stable, accurate and has excellent shock-capturing capabilities for both viscous and inviscid flows

Myofibrillar Protein Status of the Gastrocnemius in Male Rats: Effect of Mild Undernutrition

The aim of this work was the determination of the myofibrillar protein profiles in the fed and the mildly underfed rat. Sixteen male rats were divided into 2 groups: CR (control) fed ad libitum and MR (mildly undernourished) fed 75% of energetic maintenance needs. The animals were sacrificed at day 23 and the gastrocnemius muscle was taken for myofibrillar protein characterisation. The myofibrillar protein profiles were found to be very similar in the two groups revealing the lack of preferred catabolism of myofibrillar proteins and consequently that the muscle structure is maintained even in situations of mild undernutrition

Dispersive wave runup on non-uniform shores

Historically the finite volume methods have been developed for the numerical integration of conservation laws. In this study we present some recent results on the application of such schemes to dispersive PDEs. Namely, we solve numerically a representative of Boussinesq type equations in view of important applications to the coastal hydrodynamics. Numerical results of the runup of a moderate wave onto a non-uniform beach are presented along with great lines of the employed numerical method (see D. Dutykh et al. (2011) for more details).Comment: 8 pages, 6 figures, 18 references. This preprint is submitted to FVCA6 conference proceedings. Other author papers can be downloaded at http://www.lama.univ-savoie.fr/~dutykh