Linear stability of WENO schemes coupled with explicit Runge–Kutta schemes

This paper focuses on the results of the linear stability analysis of the finite‐difference weighted essentially non‐oscillatory (WENO) schemes with optimal weights. The standard WENO schemes between the third and 11th order, the order‐optimised WENO schemes of the sixth and eighth order and the bandwidth‐optimised WENO schemes of the third and fourth order are considered. Several explicit Runge–Kutta schemes including the recently published strong stability‐preserving explicit Runge–Kutta schemes are considered for time discretisation. The stability limits as well as dissipation and dispersion properties dependent on the Courant–Friedrichs–Lewy number are presented for a hyperbolic model equation. The different combinations of space and time discretisation schemes are compared in terms of their accuracy and efficiency. For a parabolic model equation, the viscous term is discretised with high‐order central differences. The stability limits for the parabolic problem are presented as well. Numerical results of linear test cases are shown. Copyright © 2011 John Wiley & Sons, Ltd.     

Uncertainty quantification in a chemical system using error estimate-based mesh adaption

This paper describes a rigorous a posteriori error analysis for the stochastic solution of non-linear uncertain chemical models. The dual-based a posteriori stochastic error analysis extends the methodology developed in the deterministic finite elements context to stochastic discretization frameworks. It requires the resolution of two additional (dual) problems to yield the local error estimate. The stochastic error estimate can then be used to adapt the stochastic discretization. Different anisotropic refinement strategies are proposed, leading to a cost-efficient tool suitable for multi-dimensional problems of moderate stochastic dimension. The adaptive strategies allow both for refinement and coarsening of the stochastic discretization, as needed to satisfy a prescribed error tolerance. The adaptive strategies were successfully tested on a model for the hydrogen oxidation in supercritical conditions having 8 random parameters. The proposed methodologies are however general enough to be also applicable for a wide class of models such as uncertain fluid flows.     

Bi-Sobolev solutions to the prescribed Jacobian inequality in the plane with Lp data and applications to nonlinear elasticity

We construct planar bi-Sobolev mappings whose local volume distortion is bounded from below by a given function $f\in L^p$ with $p>1$. More precisely, for any $1<q<(p+1)/2$ we construct $W^{1,q}$-bi-Sobolev maps with identity boundary conditions; for $f\in L^{\infty }$, we provide bi-Lipschitz maps. The basic building block of our construction are bi-Lipschitz maps which stretch a given compact subset of the unit square by a given factor while preserving the boundary. The construction of these stretching maps relies on a slight strengthening of the celebrated covering result of Alberti, Csörnyei, and Preiss for measurable planar sets in the case of compact sets. We apply our result to a model functional in nonlinear elasticity, the integrand of which features fast blowup as the Jacobian determinant of the deformation becomes small. For such functionals, the derivation of the equilibrium equations for minimizers requires an additional regularization of test functions, which our maps provide.     

Back to the future: A retrospective assessment of model‐based scenarios for the management of the shrimp fishery in French Guiana facing global change

While the number of models dedicated to predicting the consequences of alternative resource management strategies has increased, instances in which authors look back at past predictions to learn from discrepancies between these and observed developments are scarce. In the past decades, the French Guiana shrimp fishery has experienced shrimp market globalization and decreasing levels of shrimp recruitment due to environmental changes. In 2006, a bio‐economic model of this fishery was developed to simulate its possible responses to economic and environmental scenarios up to 2016. Here, we compare here these predictions to the observed trajectories. While the number of active vessels corresponds to that which was predicted, the estimated shrimp stock does not. Important driving factors had not been anticipated, including a general strike, natural disasters, and the end of the global financial crisis. These results show the importance of participative approaches involving stakeholders in the co‐construction and shared representation of scenarios. Recommendations for resource managers

• Effective fisheries resources management and a fortiori, the capacity of the fisheries to adapt to global change, requires understanding of both ecological and economics dynamics.
• The temporal trajectory of the trawling shrimp fisheries has been well monitored, and the decline of both stock and fleet is understood regarding ecological and economic changes: Changes in the environmental conditions of shrimp recruitment, and oil price increase and selling price decrease.
• However, our bio‐economic modeling work showed that, even with a good understanding of the dynamics explaining past trajectories, unpredictable events (strike, natural disasters…) have acted as other key driving factors altering the capacity of the model to represent possible futures.
• These results led us to recommend a better integration of the expertise of social and political scientists in developing models of bio‐economic systems to increase the quality of scenario predictions, and to argue for more participative approaches involving the stakeholders.

     

Conservative numerical methods for the Full von Kármán plate equations

This article is concerned with the numerical solution of the full dynamical von Kármán plate equations for geometrically nonlinear (large‐amplitude) vibration in the simple case of a rectangular plate under periodic boundary conditions. This system is composed of three equations describing the time evolution of the transverse displacement field, as well as the two longitudinal displacements. Particular emphasis is put on developing a family of numerical schemes which, when losses are absent, are exactly energy conserving. The methodology thus extends previous work on the simple von Kármán system, for which longitudinal inertia effects are neglected, resulting in a set of two equations for the transverse displacement and an Airy stress function. Both the semidiscrete (in time) and fully discrete schemes are developed. From the numerical energy conservation property, it is possible to arrive at sufficient conditions for numerical stability, under strongly nonlinear conditions. Simulation results are presented, illustrating various features of plate vibration at high amplitudes, as well as the numerical energy conservation property, using both simple finite difference as well as Fourier spectral discretizations. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1948–1970, 2015     

Numerical study of two-airfoil arrangements by a discrete vortex method

Theoretical and Computational Fluid Dynamics - The aerodynamic characteristics of two neighboring airfoils are greatly different from those of a single airfoil, for both attached and detached flow...     

Sample Entropy as a Tool to Assess Lumbo-Pelvic Movements in a Clinical Test for Low-Back-Pain Patients

Low back pain (LBP) obviously reduces the quality of life but is also the world’s leading cause of years lived with disability. Alterations in motor response and changes in movement patterns are expected in LBP patients when compared to healthy people. Such changes in dynamics may be assessed by the nonlinear analysis of kinematical time series recorded from one patient’s motion. Since sample entropy (SampEn) has emerged as a relevant index measuring the complexity of a given time series, we propose the development of a clinical test based on SampEn of a time series recorded by a wearable inertial measurement unit for repeated bending and returns (b and r) of the trunk. Twenty-three healthy participants were asked to perform, in random order, 50 repetitions of this movement by touching a stool and another 50 repetitions by touching a box on the floor. The angular amplitude of the b and r movement and the sample entropy of the three components of the angular velocity and acceleration were computed. We showed that the repetitive b and r “touch the stool” test could indeed be the basis of a clinical test for the evaluation of low-back-pain patients, with an optimal duration of 70 s, acceptable in daily clinical practice.     

The adsorption of a substituted benzene,the ethynyl-trifluoro-toluene on Si(100)-2 × 1

The adsorption of 3-ethynyl-trifluoro-toluene (ETFT) on Si(100)-2 × 1 surface in ultra high vacuum is studied in the low coverage regime, through a joint experimental and theoretical approach. The STM images of both filled and empty states revealed few distinct adsorption configurations. On the basis of Density Functional Theory (DFT) calculations the STM images were simulated and three main adsorption configurations were identified, with a predominance of di-sigma bonded species that leave the benzene ring unreacted. A discussion of the reactivity of the reconstructed silicon surface towards benzene derivatives is proposed by comparing the adsorption of ETFT close related molecules, like styrene and phenylacetylene.     

Improved and isotropic resolution in tomographic diffractive microscopy combining sample and illumination rotation

Tomographic Diffractive Microscopy is a technique, which permits to image transparent living specimens in three dimensions without staining. It is commonly implemented in two configurations, by either rotating the sample illumination keeping the specimen fixed, or by rotating the sample using a fixed illumination. Under the first-order Born approximation, the volume of the frequency domain that can be mapped with the rotating illumination method has the shape of a “doughnut”, which exhibits a so-called “missing cone” of non-captured frequencies, responsible for the strong resolution anisotropy characteristic of transmission microscopes. When rotating the sample, the resolution is almost isotropic, but the set of captured frequencies still exhibits a missing part, the shape of which resembles that of an apple core. Furthermore, its maximal extension is reduced compared to tomography with rotating illumination. We propose various configurations for tomographic diffractive microscopy, which combine both approaches, and aim at obtaining a high and isotropic resolution. We illustrate with simulations the expected imaging performances of these configurations.