Showing papers in "Journal of Computational Physics in 2017"

TL;DR: Gaussian process priors are modified according to the particular form of such operators and are employed to infer parameters of the linear equations from scarce and possibly noisy observations, leading to model discovery from just a handful of noisy measurements.

TL;DR: A state-of-the-art platform in predictive image-based, multiscale modeling with co-designed simulations and experiments that executes on the world's largest supercomputers is discussed that can be the basis of Virtual Materials Testing standards and aids in the development of new material formulations.

TL;DR: This paper develops a first and second order time-stepping scheme based on the “Invariant Energy Quadratization” (IEQ) method, and proves that all proposed schemes are unconditionally energy stable.

TL;DR: An Optimized PS (OPS) scheme is presented for accurate and consistent implementation of particle shifting for free-surface flows, consistently resulting in perfect elimination of shifting normal to an interface and resolves the unphysical discontinuity beneath the interface.

TL;DR: It is found that combining polyharmonic splines (PHS) with multivariate polynomials offers an outstanding combination of simplicity, accuracy, and geometric flexibility when solving elliptic equations in irregular regions.

TL;DR: In numerical experiments carried out on a turbulent compressible-flow problem with over one million unknowns, it is shown that increasing the time step to an intermediate value decreases both the error and the simulation time of the LSPG reduced-order model by an order of magnitude.

TL;DR: This work develops data-driven algorithms for general linear equations using Gaussian process priors tailored to the corresponding integro-differential operators that circumvents the tyranny of numerical discretization as well as the consistency and stability issues of time-integration, and is scalable to high-dimensions.

TL;DR: A unified framework for designing high order DG methods which will satisfy entropy inequalities for any given single convex entropy, through suitable numerical quadrature which is specific to this given entropy.

TL;DR: This article proposes a new mechanical constraint functional, which mimics the layer by layer construction process featured by additive manufacturing technologies, and thereby appeals to the physical origin of the difficulties caused by overhangs.

TL;DR: A series of linear, unconditionally energy stable numerical schemes for solving the classical phase field crystal model based on the first order Euler method, second order backward differentiation formulas and the second order Crank-Nicolson method are developed.

TL;DR: Estimates of spectral resolution power for under-resolved turbulent Euler flows obtained with high-order discontinuous Galerkin (DG) methods are presented and are regarded as useful guidelines for no-model DG-based simulations of free turbulence at very high Reynolds numbers.

TL;DR: This work combines geometric methods on a finite-dimensional subspace with mesh-independent infinite-dimensional approaches to speed up MCMC mixing times, while retaining robust mixing times as the dimension grows by using pCN-like methods in the complementary subspace.

TL;DR: A finite-volume method is developed for simulating the mixing of turbulent flows at transcritical conditions and an entropy-stable formulation that combines high-order non-dissipative and low-order dissipative finite- volume schemes is proposed to preserve the physical realizability of numerical solutions across large density gradients.

TL;DR: A novel computationally efficient reaction-eikonal (R-E) model for modeling extracellular potential maps and electrograms and offers vast computational savings greater than three orders of magnitude.

TL;DR: A local Lax-Friedrichs type positivity-preserving flux for compressible Navier-Stokes equations is constructed, which can be easily extended to multiple dimensions for generic forms of equations of state, shear stress tensor and heat flux.

TL;DR: A new algorithm for specifying the three-phase contact angle on curved boundaries within the framework of structured Cartesian grids is proposed, which has superior computational accuracy compared with the common approach of approximating curved boundaries with stair cases.

TL;DR: It is aimed in this work to construct discrete formulations that conserve as many physical laws as possible without utilizing a strong enforcement of the divergence constraint, and doing so leads to a new formulation that conserves each of energy, momentum, angular momentum, enstrophy in 2D, helicity and vorticity.

TL;DR: This study systematically compares accuracy and cost of the high-order Flux Reconstruction solver PyFR running on GPUs and the industry-standard solver STAR-CCM+ running on CPUs when applied to a range of unsteady flow problems.

TL;DR: A Fourier pseudo-spectral method that conserves mass and energy is developed for a two-dimensional nonlinear Schrodinger equation and it is proved that the optimal rate of convergence is in the order of O in the discrete L 2 norm without any restrictions on the grid ratio.

TL;DR: A low-dissipation weakly-compressible SPH method for modeling free-surface flows exhibiting violent events such as impact and breaking that is compatible with the hydrostatic solution and exhibits considerably less numerical damping of the mechanical energy than previous methods.

TL;DR: A simple and efficient interface-fitted mesh generation algorithm which can produce a semi-structured interface- fitted mesh in two and three dimensions quickly is developed in this paper.

TL;DR: A simple approach is introduced to correct the implementation of this drag model for systems with two- way coupling, and it is shown that this corrected implementation provides accurate and grid-independent predictions of particle settling in two-way coupled flows at low particle Reynolds numbers.

TL;DR: The results show that the developed sparse regression technique is able to identify the most significant PC contributions describing the problem and the most important stochastic features are captured at a reduced computational cost compared to the LAR method.

TL;DR: A new arbitrary high order accurate semi-implicit spacetime discontinuous Galerkin (DG) method for the solution of the two and three dimensional compressible Euler and NavierStokes equations on staggered unstructured curved meshes is proposed and is able to deal with all Mach number flows.

TL;DR: The optimized cumulant lattice Boltzmann method with fourth order accurate diffusion is used to simulate the flow around a sphere up to Reynolds number 106 and it is demonstrated that the drag crisis corresponds to an almost discrete jump in the flow conditions.

TL;DR: In this paper, an arbitrary high-order accurate nodal discontinuous Galerkin spectral element approximation for the non-linear two-dimensional shallow water equations with non-constant, possibly discontinuous, bathymetry on unstructured, possibly curved, quadrilateral meshes was proposed.

TL;DR: This study presents an alternative way to perform large eddy simulation based on a targeted numerical dissipation introduced by the discretization of the viscous term, and it is shown that this regularisation technique is equivalent to the use of spectral vanishing viscosity.

TL;DR: SpECTRE as mentioned in this paper combines a discontinuous Galerkin method with a task-based parallelism model to achieve more accurate solutions for challenging relativistic astrophysics problems such as corecollapse supernovae and binary neutron star mergers.

TL;DR: This work is the first technique to include many-body hydrodynamic interactions (HIs), and the resulting fluid flows, in cellular assemblies of flexible fibers, and uses a pseudo-spectral representation of fiber positions and implicit time-stepping to resolve large fiber deformations.

TL;DR: This paper starts with a common no-model or implicit LES (iLES) DG approach with polynomial de-aliasing and Riemann solver dissipation, and finds that the strategy gives excellent results, but only when the resolution is such, that about 40% of the dissipation is resolved.