Michael

Hi everyone,

I'm running a large model with steady flow and transient transport. For the transport, I'm using a reverse flow field and I'm simulating transport in backward from 7 pumping wells in order to estimate the contributing area of each well.

Feflow allow an easy convergence of flow but the transport result is an horror. There's many log messages during the run indicating error with BICGSTABD (default) over than 1 e60. Results are incoherent with positive values and negatives (-1e60; 1e60).

I'm using a consistent mass and convective form of transport.

What can cause this problem ? Good idea is welcome because I'm trying different scenarios and it takes over than 12 h for each run.

Thanks, Michael

Hi,

Someone know how to make a backward advective-dispersive transport simulation.  My goal is to define the contributing area of a well using dispersive properties. Because of this, I don't want to use the backward particule tracking (advective properties only).

Thanks,

Hy,

I'm running a transient flow/transport model. So the borders of the contamianation plume are negative in mass. This situation is obviously impossible. Why this results and how to resolve this problem. If I have to put a constraint, please specify wich one.

Thanks ,

Michael

Peter,

Just another information. I modified my 8 slices to a constant elevation in all the model and I tried to run it. As I supposed, it worked well. It's just confirm the problem. Before to re-arrange my slices profile, I will wait for your advice ....

thanks, Michae;

Thanks Peter for your answer. Just another question.

With your answer, I supposed that the problem is caused by an important topographie on the roc body and by really thin layer superposition (8) in this section. Feflow indicate a mesh located in the center of a superelement. Beacause of this, I supposed the problem is not caused by superelement overlap. So, for all layers topography, Feflow must ''adapt '' the mesh in order elimate gap in the layer.

There is rule of tumbs that we should follow in order to eliminate this problem.....because I have to modify the model layers topography ?

Thanks, Michael

Hi everyone,

I have a really big model with 434 546 elements. When I'm trying to run the model, there is an a pop up window indicated :

Non positive determinant for element = 1947. Obviously, an extremely distorted element geometry as a result of false mesh design an generation.

I really dont know what to do with it . It fact, mesh generation is performed by FEFLOW and it autodetect a problem. How can I do to correct it considering that I have 434 000 elements....

Thanks,

Michael

Thanks Peter,

A little advice for other user. Mesh generation is frequently interrupted because ''add-ins lines''. If you have an important model with many rivers (50-100) it's difficult to pass all nodes and lines with superelements borders.

If you use MAPINFO, you can genrate points conresponding to lines at each nodes of the line. You can also generate points at a specific distance on the lines. After, you just have to use ''add-ins points'' with the nodes corresponding to river lines.

The result is not perfect but pretty closed to river lines. So it's more simple and it allow's us to save many time.

Michael,

Hy,

Each time I try to generate the mesh for a regional flow simulation model (30-60 km2) :-\, I have the same problem. I'm unable to genrate mesh with ''genrate automatically'' (obviously) and with ''generate areally''.

I tried different settings: with or without refinement, using deiffent nodes number....and nothing. Feflow cannot generate mesh.

I veriefied all nodes connections and there is not gap between superelement, lines and points. Someone can help me ?

Thank, Michael

Good,

It works properly.

Thanks

Hi,

I'm building a 3D model for water budget analysis. Now, I'm running in steady state in order to calibrate the model. After I wish to run it in transient  mode in order to simulate groundwater pumping effect. So, when I run the model with pumping well, I obviously observed a modification in mass balance, particularly for input comming from Head (1st Kind) Boundary situated in the border of the model.

Someone can explain how to constraint Head Boundary to a maximum influx. I'm planing to run the model in steady state and to check input and output from each border condition. After I supposed that I can constraint flux (m3/d) for each border section in order to fixed the maximum input from each border section.  Someone can tell me how to do that and if my planning is good ?

Thanks ,  Michael