From cfe76264d11037f3a6dad02315a1eb1b56b3f4c1 Mon Sep 17 00:00:00 2001 From: Sven Berger Date: Fri, 13 Dec 2024 12:56:45 +0100 Subject: [PATCH] reset to original version --- .../fluid/sphere_surface_tension_wall_2d.jl | 37 +++---------------- 1 file changed, 6 insertions(+), 31 deletions(-) diff --git a/examples/fluid/sphere_surface_tension_wall_2d.jl b/examples/fluid/sphere_surface_tension_wall_2d.jl index 25c2d9b8d..4a4eef531 100644 --- a/examples/fluid/sphere_surface_tension_wall_2d.jl +++ b/examples/fluid/sphere_surface_tension_wall_2d.jl @@ -1,5 +1,3 @@ -# In this example we try to approach the static shape of a water droplet on a horizontal plane. -# The shape of a static droplet can be calculated from the Young-Laplace equation. using TrixiParticles using OrdinaryDiffEq @@ -21,7 +19,7 @@ sound_speed = 120.0 state_equation = StateEquationCole(; sound_speed, reference_density=fluid_density, exponent=1) -sphere_radius = 0.04 +sphere_radius = 0.05 sphere1_center = (0.25, sphere_radius) sphere1 = SphereShape(fluid_particle_spacing, sphere_radius, sphere1_center, @@ -29,11 +27,8 @@ sphere1 = SphereShape(fluid_particle_spacing, sphere_radius, sphere1_center, # ========================================================================================== # ==== Fluid -# fluid_smoothing_length = 1.0 * fluid_particle_spacing - eps() -# fluid_smoothing_kernel = SchoenbergCubicSplineKernel{2}() - -fluid_smoothing_length = 3.25 * fluid_particle_spacing -fluid_smoothing_kernel = WendlandC2Kernel{2}() +fluid_smoothing_length = 1.0 * fluid_particle_spacing - eps() +fluid_smoothing_kernel = SchoenbergCubicSplineKernel{2}() # For perfect wetting # nu = 0.0005 @@ -45,34 +40,14 @@ alpha = 8 * nu / (fluid_smoothing_length * sound_speed) # `adhesion_coefficient = 0.001` and `surface_tension_coefficient = 2.0` for no wetting viscosity = ArtificialViscosityMonaghan(alpha=alpha, beta=0.0) -# sphere_surface_tension = WeaklyCompressibleSPHSystem(sphere1, ContinuityDensity(), -# state_equation, fluid_smoothing_kernel, -# fluid_smoothing_length, -# viscosity=viscosity, -# acceleration=(0.0, -gravity), -# surface_tension=SurfaceTensionMorris(surface_tension_coefficient=10.0), -# reference_particle_spacing=fluid_particle_spacing, -# surface_normal_method=ColorfieldSurfaceNormal(boundary_contact_threshold=0.05, -# ideal_density_threshold=0.75)) - sphere_surface_tension = WeaklyCompressibleSPHSystem(sphere1, ContinuityDensity(), state_equation, fluid_smoothing_kernel, fluid_smoothing_length, viscosity=viscosity, acceleration=(0.0, -gravity), - surface_tension=SurfaceTensionMomentumMorris(surface_tension_coefficient=1.0), - reference_particle_spacing=fluid_particle_spacing, - surface_normal_method=ColorfieldSurfaceNormal(boundary_contact_threshold=0.05, - ideal_density_threshold=0.75)) - -# sphere_surface_tension = WeaklyCompressibleSPHSystem(sphere1, ContinuityDensity(), -# state_equation, fluid_smoothing_kernel, -# fluid_smoothing_length, -# viscosity=viscosity, -# acceleration=(0.0, -gravity), -# surface_tension=SurfaceTensionAkinci(surface_tension_coefficient=2.0), -# correction=AkinciFreeSurfaceCorrection(fluid_density), -# reference_particle_spacing=fluid_particle_spacing) + surface_tension=SurfaceTensionAkinci(surface_tension_coefficient=2.0), + correction=AkinciFreeSurfaceCorrection(fluid_density), + reference_particle_spacing=fluid_particle_spacing) trixi_include(@__MODULE__, joinpath(examples_dir(), "fluid", "falling_water_spheres_2d.jl"),