A compatibly differenced total energy conserving form of SPH

We describe a modified form of smoothed particle hydrodynamics (SPH) in which the specific thermal energy equation is based on a compatibly differenced formalism, guaranteeing exact conservation of the total energy. We compare the errors and convergence rates of the standard and compatible SPH formalisms on a variety of shock test problems with analytic answers. We find that the new compatible formalism reliably achieves the expected first‐order convergence for these analytic shock tests and, in all cases, improves the accuracy of the numerical solution over the standard formalism. We also examine the performance of our new formalism on a more complicated applied problem: the diversion of an asteroid by a kinetic impactor. We find the compatible discretization demonstrates measurable improvement in the convergence of properties such as the deflection velocity in this kind of applied problem as well. Copyright © 2014 John Wiley & Sons, Ltd.     

Monotone nonlinear finite‐volume method for nonisothermal two‐phase two‐component flow in porous media

This article presents a new nonlinear finite‐volume scheme for the nonisothermal two‐phase two‐component flow equations in porous media. The face fluxes are approximated by a nonlinear two‐point flux approximation, where transmissibilities nonlinearly depend on primary variables. Thereby, we mainly follow the ideas proposed by Le Potier combined with a harmonic averaging point interpolation strategy for the approximation of arbitrary heterogeneous permeability fields on polygonal grids. The behavior of this interpolation strategy is analyzed, and its limitation for highly anisotropic permeability tensors is demonstrated. Moreover, the condition numbers of occurring matrices are compared with linear finite‐volume schemes. Additionally, the convergence behavior of iterative solvers is investigated. Finally, it is shown that the nonlinear scheme is more efficient than its linear counterpart. Copyright © 2016 John Wiley & Sons, Ltd.     

A refined r‐factor algorithm for TVD schemes on arbitrary unstructured meshes

A refined r‐factor algorithm for implementing total variation diminishing (TVD) schemes on arbitrary unstructured meshes, referred to henceforth as a face‐perpendicular far‐upwind interpolation scheme for arbitrary meshes (FFISAM), is proposed based on an extensive review of the existing r‐factor algorithms available in the literature. The design principles, as well as the respective advantages and disadvantages, of the existing algorithms are first systematically analyzed before presenting the FFISAM. The FFISAM is designed to combine the merits of various existing r‐factor algorithms. The performance of the FFISAM, implemented in 10 classical TVD schemes, is evaluated against four two‐dimensional pure‐advection benchmark test cases where analytical solutions are available. The numerical results clearly show that the FFISAM leads to a better overall performance than the existing algorithms in terms of accuracy and convergence on arbitrary unstructured meshes for the 10 classical TVD schemes. Copyright © 2015 John Wiley & Sons, Ltd.     

Periodic excitation of a buckled beam using a higher order semianalytic approach

This paper considers the steady-state behavior of a transversally excited, buckled pinned–pinned beam, which is free to move axially on one side. This research focuses on higher order single-mode as well as multimode Galerkin discretizations of the beam’s partial differential equation. The convergence of the static load-paths and eigenfrequencies (of the linearized system) of the various higher-order Taylor approximations is investigated. In the steady-state analyses of the semianalytic models, amplitude–frequency plots are presented based on 7th order approximations for the strains. These plots are obtained by solving two-point boundary value problems and by applying a path-following technique. Local stability and bifurcation analysis is carried out using Floquet theory. Dynamically interesting areas (bifurcation points, routes to chaos, snapthrough regions) are analyzed using phase space plots and Poincaré plots. In addition, parameter variation studies are carried out. The accuracy of some semianalytic results is verified by Finite Element analyses. It is shown that the described semianalytic higher order approach is very useful for fast and accurate evaluation of the nonlinear dynamics of the buckled beam system.     

悬索在外激励作用下的1∶3内共振分析Ⅰ∶离散法

研究了悬索在受到外激励作用和考虑1∶3内共振情况下的两模态非线性响应。对于一定范围内的悬索弹性-几何参数而言,悬索第三阶面内对称模态的固有频率接近于第一阶面内对称模态的固有频率的3倍,从而导致1∶3内共振的存在。首先利用Galerkin方法把悬索的面内运动方程进行离散,然后利用多尺度法对离散的运动方程进行摄动,可得到两组不同主共振情况下的平均方程。     

A Robust High-Resolution Split-Type Compact FD Scheme for Spatial Direct Numerical Simulation of Boundary-Layer Transition

In this paper an efficient split-type Finite-Difference (FD) scheme with high modal resolution – most important for the streamwise convection terms that cause wave transport and interaction – is derived for a mixed Fourier-spectral/FD method that is designed for the spatial direct numerical simulation (DNS) of boundary-layer transition and turbulence. Using a relatively simple but thorough and instructive modal analysis we discuss some principal trouble sources of the related FD discretization. The new scheme is based on a 6th-order compact FD discretization in streamwise and wall-normal direction and the classical 4th-order Runge–Kutta time-integration scheme with symmetrical final corrector step. Exemplary results of a fundamental-(K-) type breakdown simulation of a strongly decelerated Falkner–Skan boundary layer (Hartree parameter _H = – 0.18) using 70 mega grid points in space are presented up to the early turbulent regime (Re_,turb 820). The adverse pressure gradient gives rise to local separation zones during the breakdown stage and intensifies final breakdown by strong amplification of (background) disturbances thus enabling rapid transition at moderate Reynolds number. The appearance and dynamics of small-scale vortical structures in early turbulence basically similar to the large-scale structures at transition can be observed corroborating Kachanov's hypothesis on the importance of the K-regime of breakdown for coherent structures in turbulence.     

Explicit one-step time discretizations for discontinuous Galerkin and finite volume schemes based on local predictors

We consider a family of explicit one-step time discretizations for finite volume and discontinuous Galerkin schemes, which is based on a predictor-corrector formulation. The predictor remains local taking into account the time evolution of the data only within the grid cell. Based on a space–time Taylor expansion, this idea is already inherent in the MUSCL finite volume scheme to get second order accuracy in time and was generalized in the context of higher order ENO finite volume schemes. We interpret the space–time Taylor expansion used in this approach as a local predictor and conclude that other space–time approximate solutions of the local Cauchy problem in the grid cell may be applied. Three possibilities are considered in this paper: (1) the classical space–time Taylor expansion, in which time derivatives are obtained from known space-derivatives by the Cauchy–Kovalewsky procedure; (2) a local continuous extension Runge–Kutta scheme and (3) a local space–time Galerkin predictor with a version suitable for stiff source terms. The advantage of the predictor–corrector formulation is that the time evolution is done in one step which establishes optimal locality during the whole time step. This time discretization scheme can be used within all schemes which are based on a piecewise continuous approximation as finite volume schemes, discontinuous Galerkin schemes or the recently proposed reconstructed discontinuous Galerkin or P_NP_M schemes. The implementation of these approaches is described, advantages and disadvantages of different predictors are discussed and numerical results are shown.     

Efficient solution strategy for the semi-implicit discontinuous Galerkin discretization of the Navier-Stokes equations

We deal with the numerical solution of the system of the compressible Navier–Stokes equations with the aid of the interior penalty Galerkin method. We employ a semi-implicit time discretization which leads to the solution of a sequence of linear algebraic systems. We develop an efficient strategy for the solution of these systems. It is based on a simple adaptive technique for the choice of the time step and a relatively weak stopping criterion for iterative linear algebraic solvers. The presented numerical experiments show that the proposed strategy is efficient for steady-state problems using various grids, polynomial degrees of approximations and flow regimes. Finally, we apply this strategy with a minor modification to an unsteady flow.     

Geometric, variational discretization of continuum theories

This study derives geometric, variational discretization of continuum theories arising in fluid dynamics, magnetohydrodynamics (MHD), and the dynamics of complex fluids. A central role in these discretizations is played by the geometric formulation of fluid dynamics, which views solutions to the governing equations for perfect fluid flow as geodesics on the group of volume-preserving diffeomorphisms of the fluid domain. Inspired by this framework, we construct a finite-dimensional approximation to the diffeomorphism group and its Lie algebra, thereby permitting a variational temporal discretization of geodesics on the spatially discretized diffeomorphism group. The extension to MHD and complex fluid flow is then made through an appeal to the theory of Euler-Poincaré systems with advection, which provides a generalization of the variational formulation of ideal fluid flow to fluids with one or more advected parameters. Upon deriving a family of structured integrators for these systems, we test their performance via a numerical implementation of the update schemes on a cartesian grid. Among the hallmarks of these new numerical methods are exact preservation of momenta arising from symmetries, automatic satisfaction of solenoidal constraints on vector fields, good long-term energy behavior, robustness with respect to the spatial and temporal resolution of the discretization, and applicability to irregular meshes.     

Hyperfinite construction of G-expectation

The hyperfinite G-expectation is a nonstandard discrete analogue of G-expectation (in the sense of Robinsonian nonstandard analysis). A lifting of a continuous-time G-expectation operator is defined as a hyperfinite G-expectation which is infinitely close, in the sense of nonstandard topology, to the continuous-time G-expectation. We develop the basic theory for hyperfinite G-expectations and prove an existence theorem for liftings of (continuous-time) G-expectation. For the proof of the lifting theorem, we use a new discretization theorem for the G-expectation (also established in this paper, based on the work of Dolinsky et al. [Weak approximation of G-expectations, Stoch. Process. Appl. 122(2) (2012), pp. 664–675]).