CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

Using value at certain location inside equation

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   March 23, 2024, 00:00
Default Using value at certain location inside equation
Song Young Ik
Join Date: Apr 2022
Location: South Korea
Posts: 58
Rep Power: 4
songyi719 is on a distinguished road
I am trying to include some turbulent flux term of scalar by modifying scalartransport, but having some difficulty. There are two parts that is difficult to implement, one is differentiating for only single direction, and other is normalizing by center value.

I am trying to neglect turbulent flux in other direction, but leave only x direction u's' term, so I have to calculate only ds/dx, not fvm::div ds/dx+ds/dy. Thankfully, ChatGPT tells me that I can use fvm::grad().component(DIM_Y), so I think I will try this one.

Actually, real issue is that closure equation of this scalar use z direction velocity value at certain point, such as (0,0,z) when equation is solved at location (x,y,z). I first thought of using findCell, but it wouldn't work for parallel computation. How can I implement this to work at parallel condition? Luckily, z dimension of cell is equal along the x and y axis. But I can't just divide coordinate by dz to find index because height ratio exists along z axis.
songyi719 is offline   Reply With Quote


findcell, fvm, parallel

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
sliding mesh problem in CFX Saima CFX 46 September 11, 2021 07:38
decomposePar problem: Cell 0contains face labels out of range vaina74 OpenFOAM Pre-Processing 37 July 20, 2020 05:38
[snappyHexMesh] SHM error: inside mesh not possible Naresh yathuru OpenFOAM Meshing & Mesh Conversion 1 January 10, 2017 23:58
Some problem of "Qcriterion.mcr& yuyuxuan Tecplot 9 February 12, 2016 03:27
Question about heat transfer coefficient setting for CFX Anna Tian CFX 1 June 16, 2013 06:28

All times are GMT -4. The time now is 22:26.