7 results
Now showing1 - 7 of 7
- Some of the metrics are blocked by yourconsent settings
Item type:Publication, Access status: without files , An incremental singular value decomposition approach for large-scale spatially parallel & distributed but temporally serial data – applied to technical flowsThe article presents a strategy and its algorithm to compile a simulation-accompanying, incremental Singular Value Decomposition (SVD) for time-evolving, spatially parallel discrete data sets. The framework addresses state-of-the-art PDE solvers for computational science and engineering applications. An important characteristic of such applications is that the spatial size of the data is often time-invariant and significantly exceeds the temporal size due to the large computational grid in 3D applications. Typical examples, which are also considered in this article, relate to results extracted from unsteady flow simulations. Herein, the flow data, which progresses over time, is frequently calculated spatially parallel based on domain decomposition strategies, which allow to parallelize the simulation on distributed memory machines following a Single Instruction Multiple Data (SIMD) concept. With a view to the memory-efficient reuse of (compressed) simulation results and their CPU time-saving, sufficiently accurate generation, the paper scrutinizes the efficiency of incremental/parallel SVD approaches for such simulation examples. To improve the computational efficiency, the introduction of a bunch matrix is proposed, which enables the aggregation of multiple time steps and SVD updates, and significantly increases the efficiency. The suggested strategy is verified and validated by simple 2D laminar single-phase flows and subsequently applied to more complex 2D and 3D turbulent two-phase flows. Emphasis is given to (a) the accuracy of SVD-based reconstruction, (b) the physical realizability of the reconstructed quantities, (c) the independence of domain partitioning, (d) an efficient snapshot bunching, and (e) related implementation aspects. In addition, the influence of lower and (adaptive) upper rank thresholds on the effort and accuracy is evaluated. A final application renders the practical benefits of the approach and refers to a merchant ship in head waves at Re = 1.4×107 and Fn = 0.26. The simulation involves 2880 processor cores and the related full-rank snapshot matrix has (108×104) entries. With a numerical overhead of O(10%), this snapshot matrix can be incrementally generated and compressed by O(95%). The compression is accompanied by only small errors in the integral force and local wave elevation of O(10−2%). This qualifies the method for an efficient subsequent data processing.Publicationtype: Journal Article31 - Some of the metrics are blocked by yourconsent settings
Item type:Publication, Access status: with files , Two-phase flow simulations of surface waves in wind-forced conditions(American Institute of Physics (AIP), 2023-07-01); ; ; Scopus© Citations 2 97 138 - Some of the metrics are blocked by yourconsent settings
Item type:Publication, Access status: with files , Cahn-Hilliard Navier-Stokes simulations for marine free-surface flows(2022-09); ; The paper is devoted to the simulation of maritime two-phase flows of air and water. Emphasis is put on an extension of the classical Volume-of-Fluid (VoF) method by a diffusive contribution derived from a Cahn-Hilliard (CH) model and its benefits for simulating immiscible, incompressible two-phase flows. Such flows are predominantly simulated with implicit VoF schemes, which mostly employ heuristic downwind-biased approximations for the concentration transport to mimic a sharp interface. This strategy introduces a severe time step restriction and requires pseudo time-stepping of steady flows. Our overall goal is a sound description of the free-surface region that alleviates artificial time-step restrictions, facilitates an efficient and robust numerical framework and inherently includes surface tension effects when needed. The approach is verified for an analytical Couette-flow example and the bubble formation under the influence of surface tension forces. 2D Validation examples are concerned with laminar standing waves reaching from gravity to capillary scale as well as a submerged hydrofoil flow. The final application refers to the 3D flow around an experimentally investigated container vessel at fixed floatation for Re=1.4E+07 and Fn=0.26. Results are compared with data obtained from VoF approaches, supplemented by analytical solutions and measurements. The study indicates the superior efficiency, resharpening capability and wider predictive realm of the CH-based extension for free surface flows with a confined spatial range of interface Courant numbers.Publicationtype: Journal ArticleScopus© Citations 4 125 99 - Some of the metrics are blocked by yourconsent settings
Item type:Publication, Access status: without files , Adjoint Complement to the Volume-of-Fluid Method for Immiscible Flows(2021-09-01); ; ; ; The paper is concerned with an adjoint complement to the Volume-of-Fluid (VoF) method for immiscible two-phase flows, e.g. air and water, which is widely used in marine engineering due to its computational efficiency. The particular challenge of the primal and the corresponding adjoint VoF-approach refers to the sharp interface treatment featuring discontinuous physical properties. Both the continuous adjoint two-phase system (integration-by-parts) and the corresponding dual compressive convection schemes (summation-by-parts) are derived for two prominent compressive convection schemes, namely the High Resolution Interface Capturing Scheme (HRIC) and Compressive Interface Capturing Scheme for Arbitrary Meshes (CICSAM). The dual scheme rigorously mirrors the primal Normalized-Variable-Diagram (NVD) stencils. Attention is restricted to steady state applications. Thus both the primal and the dual procedures are performed in pseudo time and the backward integration of the dual approach is performed around the (pseudo-temporal) converged primal field. The paper analyses the primal and adjoint equations for an engineering model problem. An analytical solution to the model problem is initially presented, which displays that the adjoint part does not offer a unique, non-trivial solution. As a remedy, an additional diffusive concentration term is introduced to the adjoint concentration equation. Numerical results obtained from the modified approach are benchmarked against the analytical solution for the model problem. Supplementary, the influence of the modification on the sensitivities obtained from simulations for the two-dimensional flow around a submerged hydrofoil are discussed. The final application refers to a shape-optimization of a generic 3D underwater vehicle and underlines a negligible influence of the free mobility parameter.Publicationtype: Journal Article86 - Some of the metrics are blocked by yourconsent settings
Item type:Publication, Access status: with files , A novel p-harmonic descent approach applied to fluid dynamic shape optimization(Springer, 2021-08-26); ; ; ; We introduce a novel method for the implementation of shape optimization for non-parameterized shapes in fluid dynamics applications, where we propose to use the shape derivative to determine deformation fields with the help of the p- Laplacian for p> 2. This approach is closely related to the computation of steepest descent directions of the shape functional in the W1,∞- topology and refers to the recent publication Deckelnick et al. (A novel W1,∞ approach to shape optimisation with Lipschitz domains, 2021), where this idea is proposed. Our approach is demonstrated for shape optimization related to drag-minimal free floating bodies. The method is validated against existing approaches with respect to convergence of the optimization algorithm, the obtained shape, and regarding the quality of the computational grid after large deformations. Our numerical results strongly indicate that shape optimization related to the W1,∞-topology—though numerically more demanding—seems to be superior over the classical approaches invoking Hilbert space methods, concerning the convergence, the obtained shapes and the mesh quality after large deformations, in particular when the optimal shape features sharp corners.Publicationtype: Journal ArticleScopus© Citations 17 120 193 - Some of the metrics are blocked by yourconsent settings
Item type:Publication, Access status: without files , Decoupling of control and force objective in adjoint-based fluid dynamic shape optimizationPublicationtype: Journal Article303 - Some of the metrics are blocked by yourconsent settings
Item type:Publication, Access status: without files , Adjoint shape optimization for fluid–structure interaction of ducted flows(2018-03-01); ; ; Based on the coupled problem of time-dependent fluid–structure interaction, equations for an appropriate adjoint problem are derived by the consequent use of the formal Lagrange calculus. Solutions of both primal and adjoint equations are computed in a partitioned fashion and enable the formulation of a surface sensitivity. This sensitivity is used in the context of a steepest descent algorithm for the computation of the required gradient of an appropriate cost functional. The efficiency of the developed optimization approach is demonstrated by minimization of the pressure drop in a simple two-dimensional channel flow and in a three-dimensional ducted flow surrounded by a thin-walled structure.Publicationtype: Review Article217