39 research outputs found
Incompressibility Enforcement for Multiple-fluid SPH Using Deformation Gradient
To maintain incompressibility in SPH fluid simulations is important for visual plausibility. However, it remains an outstanding challenge to enforce incompressibility in such recent multiple-fluid simulators as the mixture-model SPH framework. To tackle this problem, we propose a novel incompressible SPH solver, where the compressibility of fluid is directly measured by the deformation gradient. By disconnecting the incompressibility of fluid from the conditions of constant density and divergence-free velocity, the new incompressible SPH solver is applicable to both single- and multiple-fluid simulations. The proposed algorithm can be readily integrated into existing incompressible SPH frameworks developed for single-fluid, and is fully parallelizable on GPU. Applied to multiple-fluid simulations, the new incompressible SPH scheme significantly improves the visual effects of the mixture-model simulation, and it also allows exploitation for artistic controlling
Monolith: a monolithic pressure-viscosity-contact solver for strong two-way rigid-rigid rigid-fluid coupling
We propose Monolith, a monolithic pressure-viscosity-contact solver for more accurately, robustly, and efficiently simulating non-trivial two-way interactions of rigid bodies with inviscid, viscous, or non-Newtonian liquids. Our solver simultaneously handles incompressibility and (optionally) implicit viscosity integration for liquids, contact resolution for rigid bodies, and mutual interactions between liquids and rigid bodies by carefully formulating these as a single unified minimization problem. This monolithic approach reduces or eliminates an array of problematic artifacts, including liquid volume loss, solid interpenetrations, simulation instabilities, artificial "melting" of viscous liquid, and incorrect slip at liquid-solid interfaces. In the absence of solid-solid friction, our minimization problem is a Quadratic Program (QP) with a symmetric positive definite (SPD) matrix and can be treated with a single Linear Complementarity Problem (LCP) solve. When friction is present, we decouple the unified minimization problem into two subproblems so that it can be effectively handled via staggered projections with alternating LCP solves. We also propose a complementary approach for non-Newtonian fluids which can be seamlessly integrated and addressed during the staggered projections. We demonstrate the critical importance of a contact-aware, unified treatment of fluid-solid coupling and the effectiveness of our proposed Monolith solver in a wide range of practical scenarios.This work was supported in part by the Natural Sciences and Engineering Research Council of Canada (Grant RGPIN-04360-2014)
FrictionalMonolith: A Monolithic Optimization-based Approach for Granular Flow with Contact-Aware Rigid-Body Coupling
We propose FrictionalMonolith, a monolithic pressure-friction-contact solver for more accurately, robustly, and efficiently simulating two-way interactions of rigid bodies with continuum granular materials or inviscid liquids. By carefully formulating the components of such systems within a single unified minimization problem, our solver can simultaneously handle unilateral incompressibility and implicit integration of friction for the interior of the continuum, frictional contact resolution among the rigid bodies, and mutual force exchanges between the continuum and rigid bodies. Our monolithic approach eliminates various problematic artifacts in existing weakly coupled approaches, including loss of volume in the continuum material, artificial drift and slip of the continuum at solid boundaries, interpenetrations of rigid bodies, and simulation instabilities. To efficiently handle this challenging monolithic minimization problem, we present a customized solver for the resulting quadratically constrained quadratic program that combines elements of staggered projections, augmented Lagrangian methods, inexact projected Newton, and active-set methods. We demonstrate the critical importance of a unified treatment and the effectiveness of our proposed solver in a range of practical scenarios.Natural Sciences and Engineering Research Council of Canada, Grant RGPIN-2021-02524
Lagrangian Neural Style Transfer for Fluids
Artistically controlling the shape, motion and appearance of fluid
simulations pose major challenges in visual effects production. In this paper,
we present a neural style transfer approach from images to 3D fluids formulated
in a Lagrangian viewpoint. Using particles for style transfer has unique
benefits compared to grid-based techniques. Attributes are stored on the
particles and hence are trivially transported by the particle motion. This
intrinsically ensures temporal consistency of the optimized stylized structure
and notably improves the resulting quality. Simultaneously, the expensive,
recursive alignment of stylization velocity fields of grid approaches is
unnecessary, reducing the computation time to less than an hour and rendering
neural flow stylization practical in production settings. Moreover, the
Lagrangian representation improves artistic control as it allows for
multi-fluid stylization and consistent color transfer from images, and the
generality of the method enables stylization of smoke and liquids likewise.Comment: ACM Transaction on Graphics (SIGGRAPH 2020), additional materials:
http://www.byungsoo.me/project/lnst/index.htm
Development of GPU-based SPH Framework for Hydrodynamic Interactions With Non-spherical Solid Debris
μΌλ³Έμ νμΏ μλ§ μ¬κ³ μ΄ν μμλ‘ μ€λ μ¬κ³ μ λν μ°κ΅¬μ νμμ±κ³Ό λμ² λ₯λ ₯ ν보μ λν μ€μμ±μ΄ μ μ μ¦κ°νκ³ μλ€. μ¬κ³ μ λ°μν μ μλ λ
Έμ¬ μ©μ΅λ¬Ό κ±°λμ λν νκ°λ μ©μ΅λ¬Ό-μ½ν¬λ¦¬νΈ μνΈμμ©(MCCI, Molten Core Concrete Interaction)κ³Ό μ¦κΈ° νλ°λ‘λΆν°μ μμλ‘ λ
Έμ¬ λκ°μ± λ° κ±΄μ μ±μ λ°λ₯Έ μ¬μκ³ μΈ‘λ©΄μμ λ§€μ° μ€μνλ€. νΉν OPR 1000μ κ²½μ°, μ¬μ μΆ©μ 쑰건(Wet cavity condition)μ κΈ°λ³Έμ μΈ μμλ‘ μΈλ²½ λκ° λμ μ λ΅μΌλ‘ μ±νν¨μΌλ‘μ¨ ν΅μ°λ£-λκ°μ¬ μνΈμμ©(FCI, Fuel Coolant Interaction) λ°μμ΄ νμ°μ μΌλ‘ λ°μνλ κ²μΌλ‘ μλ €μ Έ μλ€. [Jin, 2014] FCI νμμ μμ ννμ ν΅μ°λ£ κ³ μ²΄ ννΈλ¬Όκ³Ό λκ°μ¬μ μνΈμμ©λΏλ§ μλλΌ, λκ°μ¬ λΉλ± νμ λ±λ ν¬ν¨νλ λ€μ 체, λ€μ νμμΌλ‘ κ·Έ νμμ΄ λ§€μ° λ³΅μ‘νλ€. μ΄ κ³Όμ μμ μμλ‘ κ±΄λ¬Ό νλΆμ κ³ μ²΄ ννΈλ¬Όμ΄ ν΄μ λμ΄ μν΄ μΈ΅μ΄ νμ±λκ³ , κ·Έ λκ°μ±μ λ°λΌ μ¬κ³ μ λ€μ μ§ν μν©μ μν₯μ μ€ μ μλ€. μ΄λ¬ν λΉκ΅¬ν κ³ μ²΄ ννΈλ¬Ό κ±°λμ λν μ΄ν΄λ₯Ό μν΄ κ°μ²΄ κ°λ
μ μ μ©ν κ³ μ²΄ ν΄μ 체κ³λ μ’μ μ κ·Όλ²μ΄ λ μ μλ€. λ°λΌμ λ³Έ μ°κ΅¬μμλ μ 체μ κ³ μ²΄ κ° μλ ₯νμ μνΈμμ© ν΄μμ μν΄ μ
μμ 체λμν(SPH, Smoothed Particle Hydrodynamics) κΈ°λ²κ³Ό κ°μ²΄μν(RBD, Rigid Body Dynamics) κΈ°λ²μ μ°κ³νμ¬ λΌκ·Έλμ§μ ν΄μ 체κ³λ₯Ό ꡬμΆνμλ€.
μνμ
μμ 체λμν κΈ°λ²μ ν΄μ μ 체λ₯Ό μ νκ°μ μ
μλ‘ ννν¨μΌλ‘μ¨ μ λμ ν΄μνλ λΌκ·Έλμ§μ ν΄μ κΈ°λ² μ€ νλμ΄λ€. κ°λ³ μ
μλ€μ μμ§μμΌλ‘ μ λμ ν΄μνλ―λ‘ λΉμ νμ λλ₯νμ λν κ³μ°μ΄ νμ μμΌλ©°, μ
μκ° μΆκ°λκ±°λ μ¬λΌμ§μ§ μλ ν ν΄μ κ³μ μ 체 μ§λμ μλμΌλ‘ 보쑴λλ€. μ΄λ¬ν λΌκ·Έλμ§μ κΈ°λ²μ νΉμ±μΌλ‘ SPH λ°©λ²μ μμ νλ©΄ μ λ, λ€μ 체 μ λ, λ€μ μ λ, νν λ³νκ° ν° μ λ λ±μ λν΄ ν΄μ μ₯μ μ κ°λλ€. λ³Έ μ°κ΅¬μμλ SPH κΈ°λ²μ μ μ©ν in-house SOPHIA μ½λλ₯Ό μ¬μ©νμ¬ λΉμμΆ λ€μ μ λ ν΄μμ μννμμΌλ©°, λ²€μΉλ§ν¬ λ°μ΄ν°μμ λΉκ΅μμ μ’μ κ²μ¦ ν΄μ κ²°κ³Όλ₯Ό 보μλ€.
κ°μ²΄μνμ μΈλ ₯μ μν΄ ννκ° λ³νμ§ μλ κ°μ²΄μ κ°λ
μ μ΄μ©νμ¬ κ³ μ²΄μ λ³μ§ μ΄λκ³Ό νμ μ΄λμ ν΄μνλ μ°κ΅¬ λΆμΌμ΄λ€. λ³Έ μ°κ΅¬μμλ μ΄μ°μμλ²(DEM, Discrete Element Method) λΆμΌμμ μ€λ μκ° λμ λ리 μ¬μ©λκ³ κ²μ¦λμλ Hertz-Mindlin μΆ©λ λͺ¨λΈμ μ μ©νμ¬ κ°μ²΄ κ° μΆ©λ ν΄μμ ꡬννμλ€. κ°μ²΄λ μ νκ°μ μ
μλ€λ‘ ννν μ μμΌλ©°, κ°μ²΄ κ° μΆ©λμ κ° κ°μ²΄λ₯Ό ꡬμ±νκ³ μλ μ
μμμ μμ μ€μ²©μ κΈ°λ°μΌλ‘ κ³μ°λλ€. λ³Έ μ°κ΅¬μμλ μ
μκΈ°λ°μ κ°μ²΄μν ν΄μ μ½λλ₯Ό μ΄μ©νμ¬ λ¨μΌ κ°μ²΄ λ° λ€μ€ κ°μ²΄ μΆ©λμ λν΄ κ²μ¦ ν΄μμ μννμμΌλ©°, ν΄μν΄ λ° λ²€μΉλ§ν¬ λ°μ΄ν° κ²°κ³Όμ μ μΌμΉν¨μ νμΈνμλ€.
μμλ ₯ λΆμΌμμ λ°μν μ μλ λΉκ΅¬ν κ³ μ²΄μ μ μ²΄κ° μνΈμμ© ν΄μμ μν΄ μμ μ€λͺ
ν SPH κΈ°λ²κ³Ό κ°μ²΄μν μ°κ³ ν΄μ μ½λλ₯Ό κ°λ°νμλ€. λ³Έ μ°κ΅¬μμ μ μ©ν μμ ν΄μ λ°©μ(Fully resolved approach)μ μ 체-κ³ μ²΄μ μμ΄ λΆλ¦¬λμ΄ μκ³ , μ 1 μ리λ₯Ό λ§μ‘±νλ―λ‘ κ³ μ²΄μ νμμ λ°λ₯Έ μκ΄μκ³Ό νλ©΄ μ λΆμ΄ νμνμ§ μλ€λ μ₯μ μ΄ μλ€. λν κ³ μ²΄ κ²½κ³λ©΄μμμ μ νν μλ ₯ κ³μ°μ μν΄ μ 체 μ
μ μ 보λ₯Ό κΈ°λ°μΌλ‘ λ
Έμ΄λ§ μλ ₯ κ²½κ³ μ‘°κ±΄μ μ μ©νμλ€. λ³Έ μ°κ΅¬μμλ μ΄λ¬ν ν΄μ λ°©μμ μ 체-κ°μ²΄ μ°κ³ ν΄μ μ½λλ₯Ό μ΄μ©νμ¬ λΉκ΅¬ν κ³ μ²΄μ μ 체μ μλ ₯νμ μνΈμμ©μ λν κ²μ¦ ν΄μμ μννμμΌλ©°, μ νλ μ€νκ³Όμ λΉκ΅μμ μ’μ κ²°κ³Όλ₯Ό 보μλ€.
μ λ ν΄μμ μν΄ λ³Έ μ°κ΅¬μ μ μ©ν SPH λ°©λ²μμλ μμλ€μ΄ λ§€μ° μ νμ μ΄κ³ μΈμ°μ (Explicit)μΌλ‘ κ³μ°μ μννκΈ° λλ¬Έμ κ° μ
μμ λν κ³μ°μ΄ λ
립μ μΌλ‘ μνλμ΄λ λ¬Έμ κ° μλ€. λ°λΌμ SPH λ°©λ²μ κ³μ° λ³λ ¬νμ μ΅μ νλ λ°©λ²μΌλ‘ μ μλ €μ Έ μμΌλ©°, λκ·λͺ¨ κ³ ν΄μλ ν΄μμ μν΄ μ΄λ νμμ μ΄λ€. λν μ
μ κΈ°λ°μ κ°μ²΄ κ³μ°μ μν΄μλ ν¨μ¨μ μΈ κ³μ° μκ³ λ¦¬μ¦μ΄ νμνλ€. λ°λΌμ λ³Έ μ°κ΅¬μμλ λκ·λͺ¨ κ³μ°κ³Ό λμ μ°μ° ν¨μ¨μ±μ μν΄ κ·Έλν½μ²λ¦¬μ₯μΉ(GPU, Graphic Processing Unit)λ₯Ό μ΄μ©νμ¬ κ³μ° λ³λ ¬νλ₯Ό μννμμΌλ©°, μ΄λ₯Ό μ΄μ©ν λ€μ€ κ³ μ²΄μ μ 체μ μνΈμμ© ν΄μμμ μ’μ κ³μ° μ±λ₯μ νμΈνμλ€.
λ³Έ μ°κ΅¬μμ μνν λΉκ΅¬ν κ³ μ²΄μ μ 체μ μλ ₯νμ μνΈμμ©μ μν GPU κΈ°λ°μ SPH ν΄μ μ½λ κ°λ° μ°κ΅¬λ₯Ό ν΅ν΄ μμλ‘ μ€λμ¬κ³ μ λ°μν μ μλ λκ°μ¬μ ν΅μ°λ£ κ³ μ²΄ ννΈλ¬Όμ μλ ₯νμ μνΈμμ© λΏλ§ μλλΌ, κ³ μ²΄ ννΈλ¬Ό κ° μνμ μνΈμμ©μ λν΄ ν¨μ¨μ μΈ ν΄μ 체κ³λ₯Ό κ°λ°νμλ€. μ΄λ₯Ό ν΅ν΄ μ΅μ 곡λ(wet cavity)μμ λ°μνλ ν΅μ°λ£ κ³ μ²΄ ννΈλ¬Όμ ν΄μ μμ©, μ°λλ―Έ μ¬κ³ λ‘ μΈν ν΄μ ꡬ쑰물μ μλ ₯νμ μνΈμμ©, κ·Έλ¦¬κ³ μΉ¨μ μ¬κ³ μ μμλ‘ κ±΄λ¬Ό λ΄ λΆμ λ¬Όμ κ±°λ λ± μμλ ₯ λΆμΌμμ λ°μν μ μλ λ€μν κ³ μ²΄-μ 체μ μλ ₯νμ μνΈμμ©μ λν ν΄μμ μ°κ΅¬μ μ μ©νκ³ κΈ°μ¬ν μ μμ κ²μΌλ‘ κΈ°λνλ€.Since the Fukushima accident, the necessity for researches on severe accidents and the importance of securing the ability to cope with the accidents have been increasing. The evaluation of the molten core behavior that may occur during the accident is very important in terms of re-criticality according to the coolability and integrity of the reactor core from the MCCI (Molten Core Concrete Interaction) and steam explosion. In the case of OPR 1000, especially, FCI (Fuel Coolant Interaction) is known to occur unconditionally because the wet cavity condition has been adopted as a basic strategy for ex-vessel cooling. [Jin, 2014] FCI is a highly complicated phenomenon, which includes multi-fluid, multi-phase interaction between the arbitrary shape of solid debris and coolant as well as coolant boiling. In this process, the debris bed is formed at the bottom of the containment, and its coolability influences the next phase of the accident. For the understanding on the solid debris behavior, a solid system with a rigid body can be a good approach for the non-spherical solid debris analysis. Therefore, in this study, Smoothed Particle Hydrodynamics (SPH) method and Rigid Body Dynamics (RBD) are coupled in a fully Lagrangian manner for the hydrodynamic interactions between fluid and solid.
Smoothed Particle Hydrodynamics (SPH) is one of the Lagrangian-based analysis methods which represents the fluid flow as a finite number of particles. Since the flow is analyzed by the motion of individual particles, there is no need to calculate the nonlinear convective term, and the total mass of the system is automatically conserved as long as particles are not added or removed. Through these Lagrangian nature, it is well known that the SPH method is effective for the free surface flow, multi-fluid and multi-phase flow, and highly deformable flow. In this study, the incompressible multi-phase flow analysis has been performed using the in-house SPH code, SOPHIA code, and V&V simulation results showed good agreement with the benchmark data.
Rigid Body Dynamics (RBD) is a research field that analyses the translation and rotation of a solid body by using the concept that a rigid body doesnβt change its shape by external forces. In this study, the collision calculation between rigid bodies is implemented by applying the Hertz-Mindlin contact force model commonly used and verified for a long time in the Discrete Element Method (DEM) field. A rigid body can be expressed as a group of finite particles, and the contact forces between solid bodies are calculated based on the small overlap of the particle pairs. Using the particle-based RBD analysis code implemented in this study, V&V simulations on single- and multi- rigid body collisions have been performed and showed good agreement with the analytical solution and the benchmark data.
To analyze the hydrodynamic interactions between non-spherical solids and fluids that can occur in the nuclear field, the integrated code has been developed by coupling RBD with SPH code. Since a fully resolved approach adopted in this study as a phase coupling method satisfies the 1st principle and the fluid-solid phase is entirely separated from each other, there is no need for the surface integral and empirical correlations depending on the solid geometry. In addition, the Neumann pressure boundary condition is implemented for accurate pressure estimation at the solid interface using the fluid particle properties. By applying the resolved SPH-RBD coupled code, V&V simulations were carried out on the hydrodynamic interactions of non-spherical solid-fluid and showed good agreement with the experimental data.
In the SPH method, since the numerical expression are highly linear and the calculations are performed explicitly, there is no problem even if the calculations for each particle are performed independently. Therefore, the SPH is well known as an optimized method for parallelization, and it is essential for large scale high-resolution simulations. In addition, an efficient computational algorithm is required for particle-based rigid body calculation. In this study, therefore, the parallelization was performed using a Graphical Processing Unit (GPU) for large-scale calculations and high computational efficiency, and it showed a good performance in analyzing the interactions of a large number of solids and fluids particles.
Through the researches on the development of a GPU-based SPH framework for the hydrodynamic interaction of non-spherical solids and fluids in this study, an efficient analysis system has been developed for not only the hydrodynamic interaction of solid corium debris with coolant but also the mechanical interaction between solid debris which can occur at the severe accidents in the nuclear reactor. By using this, it is expected that the integrated code will contribute to analytical researches on various accident scenarios that may occur in the nuclear field such as solid fuel debris sedimentation in the wet cavity, hydrodynamic interactions with coastal structures caused by the Tsunami, and the behavior of floating objects in the reactor building at the flooding accident, etc.Chapter 1 Introduction 1
1.1 Background and Motivation 1
1.2 Previous Studies 3
1.2.1 Numerical Studies on FCI Premixing Jet Breakup 3
1.2.2 Numerical Studies on Fluid-Solid Coupling with RBD 4
1.3 Objectives and Scope 5
Chapter 2 Smoothed Particle Hydrodynamics (SPH) 9
2.1 SPH Overview 9
2.1.1 Basic Concept of SPH 9
2.1.2 SPH Particle Approximation 10
2.1.3 SPH Kernel Function 12
2.1.4 SPH Governing Equations 13
2.2 SPH Multi-phase Models 16
2.2.1 Normalized Density Approach 16
2.2.2 Treatments for Multi-phase Flow 17
2.2.3 Surface Tension Force Model 18
2.3 SPH Code Implementation 20
2.3.1 Nearest Neighbor Particle Search (NNPS) 20
2.3.2 Algorithm of SPH Code 21
2.3.3 Time Integration 21
2.3.4 GPU Parallelization 22
Chapter 3 Rigid Body Dynamics (RBD) 30
3.1 RBD Overview 30
3.2 Collision Models of Rigid Body 31
3.2.1 Monaghan Boundary Force (MBF) Model 31
3.2.2 Ideal Plastic Collision Model 33
3.2.3 Impulse-based Boundary Force (IBF) Model 35
3.2.4 Penalty-based Contact Model 37
3.2.5 Determination of Collision Model 40
3.3 Algorithm of RBD 41
3.3.1 Calculation of Rigid Body Information 41
3.3.2 Contact Detection 42
3.3.3 Contact Normal Calculation 42
3.3.4 Contact Force Calculation 45
3.3.5 Summation of Rigid Body Particles 46
3.3.6 Time Integration 47
3.4 GPU Parallelization 48
3.4.1 Algorithm 1: Atomic Operation 49
3.4.2 Algorithm 2: Sorting 50
3.5 Code V&V Simulations 51
3.5.1 Conservation of Momentum & Angular Momentum 51
3.5.2 Conservation of Kinetic Energy in Elastic Collision 52
3.5.3 Bouncing Block 53
3.5.4 Sliding Block on a Slope 55
3.5.5 Collapse of Stacked Multi-body 57
Chapter 4 Two-way Coupling of SPH-RBD 75
4.1 Resolved Approach 75
4.2 Governing Equations 75
4.2.1 Solid Phase 75
4.2.2 Fluid Phase 78
4.3 Algorithm of SPH-RBD Code 78
4.4 Code V&V Simulations 81
4.4.1 Karman Vortex Problem 81
4.4.2 Water Entry 84
4.4.3 Sinking & Rotating Body 85
4.4.4 Floating & Falling Body 85
4.4.5 Collapse of Stacked Multi-body with Fluid 87
4.4.6 Code Application to Non-spherical Debris Sedimentation 89
Chapter 5 Conclusion 110
5.1 Summary 110
5.2 Recommendations 112
Nomenclature 114
Bibliography 117
κ΅λ¬Έ μ΄λ‘ 127λ°
A moving least square reproducing kernel particle method for unified multiphase continuum simulation
In physically based-based animation, pure particle methods are popular due to their simple data structure, easy implementation, and convenient parallelization. As a pure particle-based method and using Galerkin discretization, the Moving Least Square Reproducing Kernel Method (MLSRK) was developed in engineering computation as a general numerical tool for solving PDEs. The basic idea of Moving Least Square (MLS) has also been used in computer graphics to estimate deformation gradient for deformable solids. Based on these previous studies, we propose a multiphase MLSRK framework that animates complex and coupled fluids and solids in a unified manner. Specifically, we use the Cauchy momentum equation and phase field model to uniformly capture the momentum balance and phase evolution/interaction in a multiphase system, and systematically formulate the MLSRK discretization to support general multiphase constitutive models. A series of animation examples are presented to demonstrate the performance of our new multiphase MLSRK framework, including hyperelastic, elastoplastic, viscous, fracturing and multiphase coupling behaviours etc
Numerical modelling of aeration and hydroelasticity in slamming loads and responses of marine structures
Slamming plays a significant role in the ultimate and fatigue limit state design of marine structures. Despite a relatively long history of investigations, there are still gaps in knowledge and open questions in understanding the slamming phenomenon and the approach it needs in the design phase due to its complex nature and limitations of research tools. The so-called hydroelasticity effect, which is the coupled interaction of structural responses with the body of fluid on both global and local scales, is one of the main complex aspects of slamming. Variation of fluid compressibility due to the mixing of air bubbles with the fluid, called aeration, alters the slamming loads and could also affect the hydroelastic coupling. The possible interaction of the two mentioned processes affecting the slamming physics and how to approach it in the analysis of slamming is still not well understood and is the focus of investigation in this thesis.
The research methodology of this work is based on studying the details of pressure and the flow field around the slamming area and the evolution of slamming force and structural response employing numerical modeling. A numerical tool for this purpose was developed and validated against benchmark experimental data available in the literature.
In the study of hydroelasticity, local shell deformations, as well as global deformations of the structure, were studied. The interlinked effects of local hydroelasticity and aeration were investigated by performing two sets of numerical simulation campaigns on the water entry of elastic flat plates and cylindrical shell sections. Both studies revealed that local flexibility has a noticeable reducing effect on peak values of slamming pressure and forces. This reduction effect of flexibility disappears for plates in the presence of aeration, which shows a significant indication of interdependence in the roles of aeration and hydroelasticity in slamming dynamics. In plate entry and cylinder entry simulations, aeration shows a damping effect on the response strain oscillations, strengthening with increasing aeration.
Both water entry studies present new insights with valuable details on slamming load's major characteristics and local structural response of plates and cylinders. Noticeable differences between plate and cylinder entries were observed; for instance, aeration causes a substantial extension of slamming load duration in plate entries, but no meaningful change is observed in cylinder entries. Extensive parameter studies led to new functional relations to determine peak slamming pressure/force in pure and aerated water entries in terms of relatively simple power-law approximations, which have been derived for plates and cylinders. The study shows that hydroelasticity may not be an essential issue for locally stiff structures, but considering air entrainment and entrapment processes is important to determine local loading characteristics.
This thesis also presents a novel simplified model of wave slamming on an SDOF cylindrical structure. The model could reproduce the experimental slamming force and pressure time series of the large-scale wave slamming on a vertical monopile with a reasonable accuracy level. The validation study shows that the introduced simplified model could present valuable data on the physics of the interaction of a flexible cylindrical structure with impacting body of water. The simplified model was applied in a parameter study to investigate the effect of global structural vibration characteristics and aeration on wave slamming loads and structural response characteristics. The parameter study indicates that processes related to compressibility, such as aeration and air entrapment, are far more important than the structure's global flexibility. Since wave impact events in natural conditions may incorporate variable aeration levels in the water, which is shown to alter the structural response and duration of vibration, in both deterministic and stochastic studies of wave impact dynamics, the compressibility parameter is important and should be considered in the analysis
SURFSUP: Learning Fluid Simulation for Novel Surfaces
Modeling the mechanics of fluid in complex scenes is vital to applications in
design, graphics, and robotics. Learning-based methods provide fast and
differentiable fluid simulators, however most prior work is unable to
accurately model how fluids interact with genuinely novel surfaces not seen
during training. We introduce SURFSUP, a framework that represents objects
implicitly using signed distance functions (SDFs), rather than an explicit
representation of meshes or particles. This continuous representation of
geometry enables more accurate simulation of fluid-object interactions over
long time periods while simultaneously making computation more efficient.
Moreover, SURFSUP trained on simple shape primitives generalizes considerably
out-of-distribution, even to complex real-world scenes and objects. Finally, we
show we can invert our model to design simple objects to manipulate fluid flow.Comment: Website: https://surfsup.cs.columbia.edu
Implementation and Study of Optimal Parameters by Taichi of Particle-Based Simulation with Stable Fluid-Rigid Interaction
The author implements of the interlinked SPH method with the Taichi programming language and investigates of optimal parameters for simulations considering fluid-rigid-body interactions. The author implements of the interlinked SPH method proposed by Gissler et al. with the Taichi data-oriented programming language proposed by Hu et al. for particle-based simulations that simultaneously handle fluid and rigid interactions, to stabilize and speed up the simulation. Taichi is a programming language that uses hierarchical data structures to accelerate operations by utilizing the sparsity of particle data. The optimal values for the parameters of the hierarchical data structure differ from scene to scene. Therefore, the author investigates the optimal parameters for hierarchical data structures by measuring the processing speed of multiple scenes for each number of layers