
[Sponsors] 
November 12, 2012, 13:30 
Transient convective heat transport  temperature diverges

#1 
Senior Member
Håkon Strandenes
Join Date: Dec 2011
Location: Norway
Posts: 111
Rep Power: 12 
Hi,
I am trying to solve a temperature transport problem in a steady flow. The case in question is groundwater flow, where hot water (303 K) is injected through some wells, and water is extracted at some other wells (cold wells). The flow is governed by the laplace equation (laplacian(p) = 0), and then the volume fluxes q is q =  K & grad(p) (volume fluxes = average velocity). We assume that the water and soil phases is in thermal equilibrium at all times. For details regarding this, please see the USGS Hydrotherm manual: ftp://brrftp.cr.usgs.gov/pub/klkipp/...HT_Doc3.1.pdf pages 22 to 24. My intention was to solve this problem in OpenFOAM, and I created a solver for this purpose with basis in the "scalarTransportFoam" solver. The flow/pressure equation is solved once at the beginning, and then this vector field is used to solve the transient thermal equation: Code:
( (porosity*rhow*Cpw + (1porosity)*rhos*Cps)*fvm::ddt(T)  fvm::laplacian(Gamma, T) + rhow*Cpw*fvm::div(phi, T) );
The mesh is 3D, and I have tried a lot of different resolutions, adding of prism layers around the wells, making the wells square etc. The attached screenshot is a slice through the middle of the domain. The wells are created by setting a fixed pressure on the well surfaces, and then adjusting this to the desired volume flow. As previously mentioned, this works great. The BC for the temperature is a fixed temperature of 303 kelvin at the inlet wells (left) and zeroGradient at the outlet wells (right). What I find most strange is the extremely lowtemperature regions around the hot wells, almost like a flower... I have figured out that this must be a problem with the convective heat transport (the diffusion is not very important here). Can any of you help me understand why this fails so miserably? BTW: I have also tried to set the wells as sources and sinks directly in the equation with SuSpterms instead of resolving the mesh around the wells, and again, the flow was nice and smooth, while the temperature went bananas... 

November 8, 2013, 04:30 
Source terms to Laplacian Equations

#2  
Member
Kumudu
Join Date: Oct 2013
Posts: 63
Rep Power: 6 
Quote:
Hi, I am interested in your problem. Please tell me how you add these source terms to the laplacian equation. I am really stuck in my thesis , and can you send me the all the files you changed in the LaplacinFoam. It is really, helpful if you help me in this regards Kumudu 

November 8, 2013, 04:35 

#3 
Senior Member
Håkon Strandenes
Join Date: Dec 2011
Location: Norway
Posts: 111
Rep Power: 12 
Unfortunately I cannot help you. I never found out this problem, and I hadn't much time to work on it either, so I gave it up. I might have the files laying around in my archives somewhere, but I do not think that I will share them at the moment as I never finished this project.
When it comes to the modification of solvers and equations, I know that there is a lot of useful resources around on the Internet, and I seriously doubt that I can write anything better here. Good luck with your thesis. 

Tags 
groundwater, openfoam, temperature, thermal, transient 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Water subcooled boiling  Attesz  CFX  7  January 5, 2013 04:32 
Heat Flux , Convective + Radiant Heat flux  Mandeep  CFX  3  October 25, 2012 16:57 
Inverse and Transient Heat Transfer Problem on commercial software: is it possible?  rogbrito  Main CFD Forum  0  January 28, 2012 06:54 
anybody know how to define latent heat depend on variable temperature?  alloveyou  CFX  0  July 21, 2011 23:19 
Convective / Conductive Heat Transfer in Hypersonic flows  enigma  Main CFD Forum  2  November 1, 2009 23:53 