Skip to content

Re: Choosing a value for latticeU and charU


Hi Liliana,rnrnplease excuse the late answer, we’ve been busy publishing the new release.rnrnYou’re right, charU is the characteristic velocity in SI units (m/s), that you use for computing the Reynolds number. And is usually chosen to be a maximal velocity of your system. So you’ll be fine with the mean inlet velocity, if the velocity doesn’t increase much. Better would be to take the maximal velocity.rnrnLatticeU implicitly defines the time step size and is given in lattice units. It is connected to the Mach number, as 1/latticeU is approx the number of time steps a pressure wave needs to cross one lattice node. So for latticeU = 1 you get about the speed of sound (Mach 1), which would be a bad choice for LBMs don’t work well at high mach numbers. A better choice for a start is latticeU = 0.1 and slowly decreasing for more stability.rnrnCharU and latticeU are connected in a way, that charU/latticeU is a conversion factor between lattice units and SI units. Meaning that if you choose charU = 1. [m/s] and latticeU = 0.05 and during your simulation at some point you get a lattice velocity of 0.1 your physical velocity at that point is 0.1 * 1./0.05 = 2 [m/s].rnrnFurthermore, you can use CharU and latticeU to compute the time step size in SI units (seconds) as:rn# physTime = latticeU/charU * latticeLrnon the other hand, if you want a fixed time step, you can compute latticeU:rn# latticeU = physTime * charU / latticeLrnagain be careful as for large time steps your simulation will diverge.rnrnRegardsrnThomas