Source author record

Jochen Fröhlich

Jochen Fröhlich appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.flu-dyn cond-mat.soft math.NA Numerical Analysis physics.class-ph

Catalog footprint

What is connected

4works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

A collision model for DNS with ellipsoidal particles in viscous fluid

The article proposes an algorithm to model the collision between arbitrary ellipsoids in viscous fluid. It is composed of several steps, each improving upon the standard procedure employed in the current literature. First, an efficient contact detection algorithm is presented. Then, the semi-implicit Immersed Boundary Method of Tschisgale et al. (2018) is enhanced by the collision forces using the hard-sphere approach, so that the resulting model accounts for fluid forces throughout the entire collision process. Additionally, a new lubrication model is proposed that applies a constant lubrication force in the region where hydrodynamic forces are not resolved by the spatial grid. The new collision model is validated against benchmark test cases of particle wall collisions with excellent agreement. Furthermore, the collision of an ellipsoidal particle with a wall is investigated. The normal restitution coefficient in the case of ellipsoidal particles does not solely depend upon the Stokes number as in the case of spherical particles but is also a function of its shape, and the orientation before the collision. Using the new model, the effect of these parameters on the rebound trajectory is studied. It is found that the maximum normal restitution coefficient decreases significantly as the flatness of the particle increases. Also, the coefficient of restitution depends on the particle orientation, a tendency increasing with particle flatness. >>>>>>>>>>>>>>> Note, that the distance algorithm described in the paper is subject of a Corrigendum: https://www.sciencedirect.com/science/article/pii/S0301932222000313/pdfft?md5=1c1b8e790efdc6785c47297d89906fe6&pid=1-s2.0-S0301932222000313-main.pdf <<<<<<<<<<<<<<<

preprint2020arXiv

Linearizing the hybridizable discontinuous Galerkin method: A linearly scaling operator

This paper proposes a matrix-free residual evaluation technique for the hybridizable discontinuous Galerkin method requiring a number of operations scaling only linearly with the number of degrees of freedom. The method results from application of tensor-product bases on cuboidal Cartesian elements, a specific choice for the penalty parameter, and the fast diagonalization technique. In combination with a linearly scaling, face-wise preconditioner, a linearly scaling iteration time for a conjugate gradient method is attained. This allows for solutions in 1 $μs$ per unknown on one CPU core - a number typically associated with low-order methods.

preprint2015arXiv

Efficient modelling of particle collisions using a non-linear viscoelastic contact force

In this paper the normal collision of spherical particles is investigated. The particle interaction is modelled in a macroscopic way using the Hertzian contact force with additional linear damping. The goal of the work is to develop an efficient approximate solution of sufficient accuracy for this problem which can be used in soft-sphere collision models for Discrete Element Methods and for particle transport in viscous fluids. First, by the choice of appropriate units, the number of governing parameters of the collision process is reduced to one, which is a simple combination of known material parameters as well as initial conditions. It provides a dimensionless parameter that characterizes all such collisions up to dynamic similitude. Next, a rigorous calculation of the collision time and restitution coefficient from the governing equations, in the form of a series expansion in this parameter is provided. Such a calculation based on first principles is particularly interesting from a theoretical perspective. Since the governing equations present some technical difficulties, the methods employed are also of interest from the point of view of the analytical technique. Using further approximations, compact expressions for the restitution coefficient and the collision time are then provided. These are used to implement an approximate algebraic rule for computing the desired stiffness and damping in the framework of the adaptive collision model (Kempe & Fröhlich, J. Fluid Mech., 709: 445-489, 2012). Numerical tests with binary as well as multiple particle collisions are reported to illustrate the accuracy of the proposed method and its superiority in terms of numerical efficiency.

preprint2014arXiv

An analytic solution to the equations governing the motion of a point mass with quadratic resistance and generalizations

The paper is devoted to the motion of a body in a fluid under the influence of gravity and drag. Depending on the regime considered, the drag force can exhibit a linear, quadratic or even more general dependence on the velocity of the body relative to the fluid. The case of quadratic drag is substantially more complex than the linear case, as it nonlinearly couples both components of the momentum equation. Careful screening of the literature on this classical topic showed that, unexpectedly, the solutions reported do not directly provide the particle velocity as a function of time but use auxiliary quantities or apply to special cases only. No explicit solution using elementary operations on analytical expressions is known for a general trajectory. After a detailed account of the literature, the paper provides such a solution in form of a ratio of two series expansions. This result is discussed in detail and related to previous approaches. In particular, it is shown to yield, as limiting cases, certain approximate solutions proposed in the literature. The solution technique employs a strategy to reduce systems of ordinary differential equations with a triangular dependence of the right-hand side on the vector of unknowns to a single equation in an auxiliary variable. For the particular case of quadratic drag, the auxiliary variable allows an interpretation in terms of canonical coordinates of motion within the framework of Hamiltonian mechanics. Another result of the paper is the extension of the solution technique to more general drag laws, such as general power laws or power laws with an additional linear contribution. Furthermore, generalization to variable velocity of the surrounding fluid is addressed by considering a linear velocity profile for which a solution is provided as well. Throughout, the results obtained are illustrated by numerical examples.