Heat Transport Extreme Values

A discussion forum for HYDRUS users. HYDRUS is a software package for simulating water, heat and solute movement in one-, two- and three-dimensional variably saturated media. Happy Posting!
Post Reply
austin1025
Posts: 7
Joined: Wed Oct 30, 2024 7:12 pm
Location: USA

Heat Transport Extreme Values

Post by austin1025 » Fri Feb 14, 2025 2:39 am

Hello community,

I am having trouble getting my heat transport simulation to provide reasonable results. When my model runs, it initially calculates unreasonably high temperature values (10^20 degrees Celcius) for the first few minutes near the bottom of my domain and then returns nan(ind) for the rest of the simulation.

I am attempting to model an infiltration basin that receives runoff from a nearby parking lot. I currently have the parking lot inflow defined as a flux BC on a 1m section at the top of my domain. Since I am running my model for short periods of time, I have ignored evapotranspiration and define the rest of my domain as another flux BC to model direct precipitation onto the basin. Associated with each of these water flow BC are my heat transport BCs. The parking lot inflow and precipitation have different temperatures so I am currently modeling them both as separate type 3 BC with each assigned their respective vector to boundary conditions. I also have two draintiles that I am modeling as seepage faces near the bottom of my domain.

I am assuming I have set up something incorrectly for heat transport causing the model to fail, but am not sure what to fix. I have attached my model below. Any help/advice would be greatly appreciated.

Thank you in advance.

Austin
Attachments
Heat Transport Forum Post.zip
(198.83 KiB) Downloaded 29 times

Jirka
Posts: 6193
Joined: Sat Mar 16, 2002 3:47 pm
Location: USA
Location: Riverside, CA

Re: Heat Transport Extreme Values

Post by Jirka » Fri Feb 28, 2025 12:39 am

You use unrealistic (ridiculous) parameters. For example, you had the saturated water content for the third material equal to 6.8 (while it obviously must be smaller than 1). Also, you had a negative solid fraction for heat transport parameters. When I fixed these errors, the program seemed to run OK. You also did not have any heat transport BC on the seepage faces. J.

Post Reply