Skip to content

Reply To: Some doubts about codes in multiconponents example

OpenLB – Open Source Lattice Boltzmann Code Forums on OpenLB General Topics Some doubts about codes in multiconponents example Reply To: Some doubts about codes in multiconponents example

#4970
Zhangshi
Participant

Dear Mathias,

I’m sorry for asking too many questions.I’m glad to going to spring school for further learning.But it seems that I don’t have enough time to get there.
Some of the questions above have been solved by myself.But there are still two problems.Do you mind listenning for a while?

First is question 2.

`AnalyticalConst3D<T,T> one ( 1. );
SmoothIndicatorSphere3D<T,T> sphere( {nx/2., nx/2., nx/2.}, radius, 10.*alpha );
AnalyticalIdentity3D<T,T> rho( one );
AnalyticalIdentity3D<T,T> phi( one – sphere – sphere );
sLattice1.iniEquilibrium( superGeometry, 1, rho, zeroVelocity );
sLattice2.iniEquilibrium( superGeometry, 1, phi, zeroVelocity );`

AnalyticalIdentity3D<T,T> phi( one - sphere - sphere );
It’s used for constructing spherical erea.But “one” is not the same kind as the sphere.And why here subtract sphere for twice?

2.Something about question 1 is really confusing.How to get the value of these variable in lattice units?Would yuo give me some advice?

Thanks for you reply!
Best
Zhangshi

  • This reply was modified 4 years, 2 months ago by Zhangshi.