From 9c3270d0350910861f4250a120f32bf03854f9c8 Mon Sep 17 00:00:00 2001 From: Yann Gaillard Date: Fri, 7 Feb 2025 05:35:12 +0100 Subject: [PATCH] Correct the centrifugal buoyancy term The centrifugal buoyancy in a rotating reference frame commonly defined as: beta*DeltaT*Omega^2*d^4*/nu^2*r*sin(theta)[-sin(theta),-cos(theta),0] is corrected according to the scaling [r^2;sin(theta)/r,sin(theta)/r]. --- src/get_nl.f90 | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/src/get_nl.f90 b/src/get_nl.f90 index 3c40ca2c..d343ed14 100644 --- a/src/get_nl.f90 +++ b/src/get_nl.f90 @@ -366,9 +366,11 @@ subroutine get_nl(this, time, nR, nBc, lRmsCalc) ! !-- neglect pressure contribution ! !& + polind*DissNb*oek*opressure0(nR)*this%pc(:,nPhi) ) !else - this%CAr(:,nPhi) = -dilution_fac*r(nR)*sinTheta(:)**4*ra*opr* & + ! Pr*Di/Ra*r^3*sin(theta)^2*T + this%CAr(:,nPhi) = -dilution_fac*r(nR)**3*sinTheta(:)**2*ra*opr* & & this%sc(:,nPhi) - this%CAt(:,nPhi) = -dilution_fac*r(nR)*sinTheta(:)**3*cosTheta(:)*ra*opr* & + ! Pr*Di/Ra*sin(theta)^2*cos(theta)*T + this%CAt(:,nPhi) = -dilution_fac*sinTheta(:)**2*cosTheta(:)*ra*opr* & & this%sc(:,nPhi) !end if end if ! centrifuge