Skip to content

nv4dll

Forum Replies Created

Viewing 9 posts - 1 through 9 (of 9 total)
  • Author
    Posts
  • in reply to: OpenLB including #4783
    nv4dll
    Participant

    I had resend the email via gmail , please check for it.

    nv4dll
    Participant

    Dear Sam,

    Thanks for your help and time! It helped me a lot.

    Best wishes,

    nv4dll

    nv4dll
    Participant

    If I want to convert the surface tension converted into lattice units which fator shold I use? Since The units of surface tension is N/M.
    Should I use 0.133333/0.0144?

    
    Voxel length(m):                  physDeltaX=     0.133333
    Time step(s):                     physDeltaT=     0.148148
    Velocity factor(m/s):             physVelocity=   0.9
    Density factor(kg/m^3):           physDensity=    1
    Mass factor(kg):                  physMass=       0.00237037
    Viscosity factor(m^2/s):          physViscosity=  0.12
    Force factor(N):                  physForce=      0.0144
    Pressure factor(N/m^2):           physPressure=   0.81
    
    nv4dll
    Participant

    Dear Sam,
    Thank you very much! When I set kappa1=kappa2 it do solve the problem!
    But I got another question, let’s take contactAngle2d for example. If I use the formula you give me then the kappa1&2 should be 30, but in the code it is 0.005. Why it that? Should I convert it again?

    Best wishes,
    nv4dll

    in reply to: OpenLB including #4763
    nv4dll
    Participant

    I had send you a email, please check your inbox.
    You are welcome:)

    nv4dll
    Participant

    And I use only 2 lattice in the \microFluidics2d example and deleted all the velocity boundary except inlet1 ( dx, yl1, {0., ny/2.} );

    And I used \contactAngle2d to find some h1&h2 for different angles:

    
    h1 = 0.0002948;                  // Contact angle 120 degrees   
    h2 = -0.0002948;                  // Contact angle 120 degrees   
    
    h1 = 0.0002048;                  // Contact angle 110 degrees   
    h2 = -0.0002048;                  // Contact angle 110 degrees   
    
    h1 = 0.0001448;                  // Contact angle 100 degrees   
    h2 = -0.0001448;                  // Contact angle 100 degrees   
    
    nv4dll
    Participant

    the velocity filed: https://imgur.com/weEpqk4

    in reply to: OpenLB including #4756
    nv4dll
    Participant

    Dear iJokerify,
    I’m using Clion for debug, here is my cmakelist:

    cmake_minimum_required(VERSION 3.10)
    project(LBM)
    SET (CMAKE_USE_RELATIVE_PATHS ON)
    #OpenLB库文件
    set(LINK_DIR ${PROJECT_SOURCE_DIR}/lib/olb/build/precompiled/lib)
    link_directories(${LINK_DIR})
    set(CMAKE_CXX_STANDARD 14)
    #OpenLB库文件
    include_directories(${PROJECT_SOURCE_DIR}/lib/olb/src)
    include_directories(${PROJECT_SOURCE_DIR}/lib/olb/src/external/zlib)
    include_directories(${PROJECT_SOURCE_DIR}/lib/olb/src/external/tinyxml)
    
    include_directories(SYSTEM ${MPI_INCLUDE_PATH})
    set(CMAKE_CXX_COMPILE_FLAGS ${CMAKE_CXX_COMPILE_FLAGS} 
    #your source code PATH
    add_executable(LBM ${PROJECT_SOURCE_DIR}/src/main3d.cpp )
    
    target_link_libraries(LBM olb)
    target_link_libraries(LBM z)

    you should precompiled opeblb first and put everthing in opeblb libary under the /lib/olb/ (make new dir named lib&olb)

    or you can just give me your email I can send you files

    in reply to: Why am i keep getting inf velocity when using shan-chen model? #4729
    nv4dll
    Participant

    Dear Mathias,
    I’m very appreciate for your help. Sorry for the late reply ,I definit want to attend the school I was plan for it ,but I’m in China right now. I will see if there will be any chance depending on whether the epidemic situation has improved.

Viewing 9 posts - 1 through 9 (of 9 total)