NON DARCY Code V3 from George

[millerta]

This is new ND code from George and is on new branch gaz-nondarcy3
Files sent as is and fixed with good file names are on the server at
/project/eesdev/FEHM/FEHM_BRANCHES/FEHM_gaz-nondarcy3

From George:
NONDARCY_COMPLETE_FILES_071225/DOCS/src_changes_071225.filelist

Changes since 150625 (day.mo.yr)
The software changes since June 15, 2025(150625) are primarily due to adding non-Darcy flow to Air Water Heat simulations. The structure of the algorithm, mostly subroutine calls when generating the numerical balance equations, was simplified compared with the ND coding WH (2-phase water and heat) and isothermal 2-phase water. Use modules now provide more communication between the non-Darcy routines and subroutines that generate flowrates, assosiated parameter derivatives, and the numerical balance equations.
Several bugs were corrected. The routines listed below that the CHANGE was "none found" likely had some testing code that was used during development and hence have a newer date.

Files from SRC_051225:
air_cp.f CHANGE: none found
allocmem.f CHANGE: added and allocated small arrays needed for ND
co2ctr.f CHANGE: none found
com_nondarcy.f CHANGE: none found
com_nondarcy_AWH.f CHANGE: added variables for 2-phase and gas-only ND
comci.f CHANGE: added more non-darcy arrays
comji.f CHANGE: added more non-darcy arrays
commass_AWH.f CHANGE: added one non-darcy array
fehmn.f CHANGE: added call fluid_parameter_tests(:) before timestep
fluid_parameter_tests.f CHANGE: added subroutine fluid_parameter_tests to output fluid properties
geneq1_w_nondarcy.f CHANGE: added ND variables
geneqc_w_nondarcy.f CHANGE: This is a new subroutine that has the ND coding for the AWH fluid physics
gensl1_w_nondarcy.f CHANGE: added coding for ND
gensl2_w_nondarcy.f CHANGE: added coding for ND
gensl4_w_nondarcy.f CHANGE: added coding for ND
innondarcy_a.f added global variables to non-darcy input routine
input_w_nondarcy.f CHANGE: minor change to access variable n_awh_test
non_darcy_sub_AWH_liquid.f CHANGE: modified ND calculations to simplifly calling structure
non_darcy_sub_AWH_vapor.f CHANGE: new subroutine for ND gas flow
non_darcy_subroutines.f CHANGE: modification to local and global variables to simplifly calling structure
plotc1_w_test.f CHANGE: appended a test routine for particle tracking (not related to ND flow)
note: the calls to the test routine were commented out once particle tracing problem (actually input issues) was fixed
ptrac3.f CHANGE: coding added to check particle tracking
setparams.f CHANGE: added access to ND-related use modules
Sub_FEHM_tec_to_vtk.f CHANGE: corrected a bug; can now visualize triagles in 3D
tec_write_grid.f CHANGE: small change to identify case when no grid is available
thrmwc.f CHANGE: corrected mistake xrl =1 for gas only changed to xrl =0, corrected var_awh_param(:) values
varchk.f CHANGE: modified some phase change toterances
write_avs_node_mat_s.f CHANGE: increased length of variable grid_ghk to 200 characters
write_avs_node_v.f CHANGE: increased length of variable grid_ghk to 200 characters
wrtout.f CHANGE: corrected mispelled word in comment statement

[millerta]

In preparation of new code added to git, filenames were corrected and files compared to gaz-nondarcy2 latest ND code.
Note George's code runs on Windows and is missing some of the changes needed for Linux.

Fix file names:
mv gensl2_w_nondarcy_.f gensl2_w_nondarcy.f
mv geneq2_nondarcy_c3.f geneq2_nondarcy.f
mv thrair_c.f thrair.f

Files in gaz but missing from git (gaz-nondarcy2)
com_nondarcy_AWH.f - new file
fluid_parameter_tests.f - new file
geneqc_w_nondarcy.f - new file
gensl1_w_nondarcy.f - replaces gensl1.f
gensl2_w_nondarcy.f - replaces gensl2.f
gensl4_w_nondarcy.f - replaces gensl4.f
input_w_nondarcy.f - replaces input.f
non_darcy_sub_AWH_liquid.f - new file
non_darcy_sub_AWH_vapor.f - new file
plotc1_w_test.f - replaces plotc1.f
non_darcy_subroutines.f - replaces nd_props and nd_flow_vel in geneq2_nondarcy.f

71 files differ from previous version on git (see attached files)

diff_from_git_v3.txt

diffs_summary_v3.txt

[millerta]

Files from George:
/project/eesdev/FEHM/from_george/src_071225
es44.lanl.gov:% lt
total 39080
tamiller sft 4096 Dec 10 11:29 NONDARCY_COMPLETE_FILES_071225 (first iteration)
tamiller sft 4096 Dec 12 11:10 NONDARCY_COMPLETE_FILES_091225 (second with all files)
tamiller sft 4096 Dec 17 20:31 V+V_TESTS_NDWH
tamiller sft 4096 Dec 18 08:51 V+V_TESTS_NDAWH

ND Tests from George:
From: george zyvoloski gazyvoloski@gmail.com
Sent: Thursday, December 18, 2025 9:59 AM
To: Miller, Terry Ann tamiller@lanl.gov; Stauffer, Philip H stauffer@lanl.gov
Subject: [EXTERNAL] Test cases for ND and AWH and WH

Hi Terry and Phil. I am sending (via transfer.lanl.gov) two zip files containing V+V tests for non-Darcy flow with Air-Water-Heat (AWH) and Water-Heat (WH). They are 1D and designed to test all phases.
These should be designated as "truth". Each contains the result of two simulations: ND=OFF and ND=1.D12. Simply toggle the "off" in the input file:
ndar off
1 0 0 1.d12

end ndar

For future comparisons, I would use the history files like: "test_ND_AWH_LV_ND12_satr.his" for transient comparison, and "test_ND_AWH_LV_ND12.00002_sca_node.dat" for spatial comparison.

[millerta]

Dec 18 2025 George files renamed and pushed to gaz-nondarcy3
Windows build by Erica:
Build on windows using FEHM_gaz-nondarcy3
See /project/eesdev/elm/fehm-12-2025/src
And /project/eesdev/FEHM/FEHM_BRANCHES/FEHM_gaz-nondarcy3/src
Hi Phil, I tested Erica's executable on the historic V+V tests. Most tests passed. Some of the particle tests were off a little, much like the problem I had six months ago.

[millerta]

Make changes for Linux and check bug fixes from earlier.
Document listing differences between gaz-nondarcy2 and gaz-nondarcy3

diffs_summary_v3.docx