nv4dll
Forum Replies Created
-
AuthorPosts
-
nv4dllParticipant
I had resend the email via gmail , please check for it.
February 13, 2020 at 3:18 am in reply to: Undesired velocity at corner when changing the contact angle in free energy #4768nv4dllParticipantDear Sam,
Thanks for your help and time! It helped me a lot.
Best wishes,
nv4dll
February 12, 2020 at 2:11 pm in reply to: Undesired velocity at corner when changing the contact angle in free energy #4765nv4dllParticipantIf 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
February 12, 2020 at 2:07 pm in reply to: Undesired velocity at corner when changing the contact angle in free energy #4764nv4dllParticipantDear 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,
nv4dllnv4dllParticipantI had send you a email, please check your inbox.
You are welcome:)February 12, 2020 at 6:13 am in reply to: Undesired velocity at corner when changing the contact angle in free energy #4759nv4dllParticipantAnd 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
February 12, 2020 at 6:05 am in reply to: Undesired velocity at corner when changing the contact angle in free energy #4758nv4dllParticipantthe velocity filed: https://imgur.com/weEpqk4
nv4dllParticipantDear 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
February 3, 2020 at 11:00 am in reply to: Why am i keep getting inf velocity when using shan-chen model? #4729nv4dllParticipantDear 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. -
AuthorPosts