[liu]

thank you

[liu]

thank you

[liu]

thank you

[liu]

This line of code is on line 411 of the channeld version 1.7

[liu]

The version I am using is 1.7. thank you

[liu]

How should I pass 𝑦 plus from the NS lattice to the scalar lattice?(temperature)? thank you

[liu]

ok,thank you

[liu]

I want to add a temperature wall function in the slattice where temperature is applied. How should I proceed,My idea is to set a power law boundary function. “Implement through this code.
“ setWallFunctionBoundary<T, DESCRIPTOR>(sLattice, superGeometry, 3, converter, wallFunctionParam);”But this is a velocity wall function; how do I modify it. thank you

[liu]

thank you

[liu]

thank you， I solved the problem using your method.

[liu]

thank you, I’ve tried to include that code in the channel3d，“#define DISABLE_CSE”. But I found that I still can’t output taueff.

[liu]

namespace detail {

template <typename COLLISION, typename DESCRIPTOR, typename MOMENTA, typename EQUILIBRIUM>
struct SmagorinskyEffectiveOmega {
using MomentaF = typename MOMENTA::template type<DESCRIPTOR>;
using CollisionO = typename COLLISION::template type<DESCRIPTOR, MOMENTA, EQUILIBRIUM>;

template <typename CELL, typename PARAMETERS, typename V=typename CELL::value_t>
V computeEffectiveOmega(CELL& cell, PARAMETERS& parameters) any_platform {
V piNeqNormSqr { };
MomentaF().computePiNeqNormSqr(cell, piNeqNormSqr);
const V rho = MomentaF().computeRho(cell);
const V omega = parameters.template get<descriptors::OMEGA>();
std::cout << “omega1111111: ” << parameters.template get<descriptors::OMEGA>() << std::endl;
const V smagorinsky = parameters.template get<collision::LES::Smagorinsky>();
V piNeqNorm = util::sqrt(piNeqNormSqr);
V preFactor = smagorinsky*smagorinsky
* descriptors::invCs2<V,DESCRIPTOR>()*descriptors::invCs2<V,DESCRIPTOR>()
* 2 * util::sqrt(2);
/// Molecular realaxation time
V tauMol = V{1} / omega;

std::cout << “tauMol: ” << tauMol << std::endl << std::flush;
/// Turbulent realaxation time
V tauTurb = V{0.5} * (util::sqrt(tauMol*tauMol + preFactor / rho * piNeqNorm) – tauMol);
std::cout << “tauTurb: ” << tauTurb << std::endl << std::flush;
/// Effective realaxation time
V tauEff = tauMol + tauTurb;
return V{1} / tauEff;
} Why can’t I output taueff here, thanks

[liu]

thanks

[liu]

superGeometry.rename(2,1,{1,1} );I’d also like to ask what the {1, 1} in the sentence means. thank you very much!

[liu]

Dear Community,
I have another simple question how to define its boundary conditions for an edge in a 2D model, thanks!

[Author]

Posts