-
Boundary conditions through shapefiles and constraints
I have a large model with a large river network. The top of my model is defined by the DTM. I have a shapefile of the river network that I want to import as a Fluid-transfer into my model. How I do this in a way that the BCs will be assigned to the valleys of the DTM?<br> I tried the method I use for the hydraulic properties, selecting the whole slice and it applied the BCs to all the nodes of the slice which is not what I want.<br> <br> I also want to apply constraints to those BCs, once FEFLOW floods large areas of the model, which is not correct. However, and once the river network is big, is defined by a large number of nodes with this BC type. So it is too time-consuming define the constraints node-by-node. How can I define constraints for the whole river network in a fast way?<br> <br> And just to be sure, the values of the Fluid-transfer BCs are the water height of the river, right? So if my river has 1 m of water, that is the value I should define in the BCs, correct?<br> <br> By the way, so far I am simulating this river network as head-constant BCs and when I check the imbalances by adding the BCs values and recharge (entered through In/outflow on top/bottom) in an Excel sheet and compare with the values shown in FEFLOW Rate Budget panel I get this: I f the 1st slices is defined as Phreatic, then the Imbalance in Excel match that of the Rate Budget panel. However if I define the 1st  slice as Free, the Imbalance in Excel and in the Rate Budget panel do not match and have a mismatch from 2 to 3 orders magnitude. Why is this so? -
constraints not working
<span style="font-size:11pt"><span style="line-height:107%"><span style="font-family:"Calibri",sans-serif">I have a regional model that with the assumed K of the formations and the calculated recharge (even if divided by half) <b>floods with heads of 8 km above the model top</b> (= topographic surface), for most part of the model area.</span></span></span><br> <span style="font-size:11pt"><span style="line-height:107%"><span style="font-family:"Calibri",sans-serif">So <b><u>I removed recharge</u></b> and put the model with the FEFLOW default K = 10 m/d in all the slices of the model (31 slices). But in spite of recharge (input through <u>In/Ourtflow on top/bottom</u> applied to slice 1) = 0 m3/d and Kxx = Kyy = Kzz = 10 m/d and the other flow parameters also as FEFLOW default, <b>the model still floods</b>, altough now head above topography values range from circa 1 to 10 m. <b>It is not in the whole model area as before but it floods around the boundary conditions</b>.</span></span></span><br> <span style="font-size:11pt"><span style="line-height:107%"><span style="font-family:"Calibri",sans-serif">My boundary conditions are <u>hydraulic head = elevation</u>, are set in the 1st slice and represent the river network, once I have no data on river stages.</span></span></span><br> <span style="font-size:11pt"><span style="line-height:107%"><span style="font-family:"Calibri",sans-serif">My initial head = 0 m in the whole model and the lowest Elevation value in the 1st slice is = 0.6 m. So my initial head is always below the topographic elevation of the model top.</span></span></span><br> <span style="font-size:11pt"><span style="line-height:107%"><span style="font-family:"Calibri",sans-serif">The error tolerance is = 0.00005727</span></span></span><br> <span style="font-size:11pt"><span style="line-height:107%"><span style="font-family:"Calibri",sans-serif">The model is defined as <u>Unconfined aquifer</u>, Slice 1 as Phreatic and the remaining ones as Dependent, <u>Storage change in phreatic top layer where water table exceeds the surface</u> as Extend storage of unconfined layer to the water table (once I assume there can be areas where flloding might occour, however not with the extension and location of the present floods in the model) and I set <u>Head limits for uncofined conditions</u> as Unconstrained on Top and Bottom of the model, once below my model area there are limestones which also are aquifers.</span></span></span><br> <span style="font-size:11pt"><span style="line-height:107%"><span style="font-family:"Calibri",sans-serif">The rate budget has an imbalance of around 2550 m3/d and the BCs have outflows of circa 6.182x10<sup>6</sup> and inflows of circa 6.1824x10<sup>6</sup>.</span></span></span><br> <span style="font-size:11pt"><span style="line-height:107%"><span style="font-family:"Calibri",sans-serif">I also made an attempt with <u>Head limits for uncofined conditions</u> as Constrained head on Top of model domain (and recharge = 0) but <b>the model sitll flooded around the boundary conditions</b>.</span></span></span><br> <span style="font-size:11pt"><span style="line-height:107%"><span style="font-family:"Calibri",sans-serif">By the way, I also tried to change the K values, setting the whole model at Kx = Ky = Kz = 100 m/d and in another simulation at 1 m/d. The <b>flood occured in the same areas as when Kx = Ky = Kz = 10 m/d</b>, what changed where the values in the rate budget. Awhen there was K = 100 m/d these where 10 times larger that when K = 10 m/d and when K = 1 m/d where 10 times smaller as when K = 10 m/d.</span></span></span><br> <span style="font-size:11pt"><span style="line-height:107%"><span style="font-family:"Calibri",sans-serif">I also tried to use Seepage face but the flloding problem remained.</span></span></span><br> <span style="font-size:11pt"><span style="line-height:107%"><span style="font-family:"Calibri",sans-serif">All these simulations – except the initial one with the “real” K and recharge – were performed with recharge = 0.</span></span></span><br> <span style="font-size:11pt"><span style="line-height:107%"><span style="font-family:"Calibri",sans-serif">From all this I suspect that there is too much water entering the model through the boundary conditions. So I tried to set up constraints. However when I select the nodes with the BCs I want to constrain, then click on Hydraulic head BC in the Data panel, and click on Add paramenter and choose for instance max. flow-rate constraint in the window</span></span></span><br> <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAYwAAAAjCAYAAACHBD3FAAAgAElEQVR4Xu2dB3hU1daG35lkUiAJCSX0HjpIR6kivYOgqFx67016kd5CC5CACKEKgoIICDaKSAABCyiK9DRSCenJpEzmX3tCIJAghov/czVnP85jEmb2Ofubc9a3V/nW0fn4+Jjr1atHamoqKSkp6HQ6y0sNs9lseWX8bvmjNnIVAhnfv16v166DXPXNa4vVEMiKgO5JwrC1tcXa2pqkpCQLiWhkoV02GZuHtLQ0DQwNAQ2BXIyA7uzZs+Y6deo89DDy5MmD0Wjk1q1b2NnZUbJkSUwmUy6GKHcv3WAwEBkZyZEjR4iNjcXKyip3A6KtXkMgFyPwGGEoj0IZCEUWN2/epEyZMlSsWNFCJtrInQjY2NgQHBzMhx9+SOPGjcmfP///JBCZw6j/kyeonZSGwL8AgSwehrrx7t+/T0JCAi4uLjg6OqKFIv4F3/RzLiGDMA4cOMCQIUNwcnJ6zpm0j2kIaAj80xHIQhgqf6HyGCrJqUJTGYnwf/pCtfN/PgQyCOPTTz+lT58+FC5c+Pkm+ps/dfToUaKjo3nnnXf+5iNp02sI5F4Ess1hKDhCQ0MtqDg7O2seRu69Pvi7CCMpMZGE+xEkRkaAFOXZuhQkX6FCWBtsngtt5f3cuXOHY8eOPdfntQ9pCGgIPBuBLIShEt0hISH8/vvvFClShMqVK1u8jL9zmM06rKxtsHMwkBafQFKqCSnmfXGHNItNspL5HW0wi6FKTk7FJPM/zxHUuRoEIxtbHalxCST/nYVD6rz1BjlvWzAmkvRfnPfzgvmiCUOqtEmMiyfC7zZR134hPuAWKWlm8pavQvnGzclX0PW5TnX8+PH4+/ujPCFtaAhoCPw9CGRJeqtQlNqp3bt3z1IhVbRo0adUSVlhk9cWXUoyKckpFgOMGFNrWxtJnOswJUg4S4zDs4cOvZUJY+w9gq9HYVOuDIXy2WEwpwlpvIghxKBPI9UYSdDVcHQlSuFawAnHPNbo0kR7YkxOP/e/NITYrJKJDgomLNSEU+WyFBDi0Cu9yl/6fE7elH7epuQogn4Lg2IlcHV1wk7OVuyrZViIVnbktnlsMWRMLeeSkhQv5GJ6SIpm9BgkzGhna5CfMj6bJmuX96UonJ++/hdJGGly4sqzCAv0JeT3X0j1u4Z1VBgJ8UbMhUpQqctbFK1UBf0DHVBO0MpKGGbSUiLwu3QV30jZhJjScC5XhyrlC5PPOgczJ4Vx+49rXAuKx8omL/lcy1OlRjG0TE4OMNTe+q9BIIuHkVmglSHcy7JaZV/MiUQGRGLK64xTfgdsldnRp5IgJZgx0WnYF3PF8S/dmGoHHYPv2Y9Z3etrinmupF/LsrimiSfwQnbv1tjYJ3L/zlFWtN+F3cwFDOhVBYfrt4m0zkf+EgVx0D0ywn/+zdrgkC+Ik6vWsWNzLE33e/BGSfE40sQ4v3DGsBJPJomYoK9Z2WYrjJ9N//6NcDPHYbRUOQtR6U0kxUUSHhhOpNLNKLysbSlStiKuDjZYp6WQqrfB3sZEXGggt/3DSbRgKp+V97mWq4RrHiuszLL+pyz8RRFGhmcRec+f8JA/uB8YgCnkHrqQcOLvRZPm5Eq1N3tRunYdbA3PunBMxPiHkuzgIteePSqIlZUwUogLPsB79SZzQCr/DHbW1B/qyfQhrajumIP7N/xbNrmvYOXhaxhjrHAuOZA156bwagbz5mAq7a0aAv90BLINSanEd6LsBFUoKmvdvV52tWIgk77Hq+P7RHQYRK/xnaiBEX3eYE57vs/u7XE02b2C7qXFWKcmSYjpgbegEyNorcMsBjZVLKzOylpEgnbksYvhzvldLHv9c4qs96B/azeKWYvnkqpiMrKL1utIEy2IKU3tqK0e7pIV+GmpyZa51C7ZSgyNTjyTNAlpKQOo04sxlFCXwTqG8Fufs6z1h9jPXcLQ/nm5MGwtVwq9SqdpPamljyMFOZfMc5tSSJZdqTJ0j0Y6YZxY4cHWD+JofngDvcvJmtR6UlPkfDO8ImWQZT4rUUdbXAHxcARLC6no9HIcA1aZNvVmU6qs4RFpqfO2Nthiqwjj7lGWNNmIaeIcBgxsTIWHhGFDHqf7XPlsI4sGbOBsPkfsrAQHwaPWiJ1M79OUmgXEyzOJpxJ5lr2eXqzafYUEg1Luy8vOgeZjdjGle0VKOqSSJt+NTojDpLDLtOYXQRgWsoiXMFSoL7H3r2JK8bNoOoJuxXLvj/sQnYRLWTeqdexEueovkcfeLvv7SjzCxGjx7iLOsPGtveiHTGPw8Jcpky1hGIn0/ZyFzY5RbssoOtUuhou9o8wtRPpXHUrLBZZMolQMJiSFcGnPPj7eYOatn2fS3KB/7Dr8pxsC7fw1BP4KAtkmvePl5r527RoODg6ULVv2CR1GBmGcZnmjVYS/PoZBM7pTi0SsHO5y3H01WzfG8urhdfQoasJGDFFGxEOoQkhI8hMS9rJRu0hllI1pQjQx+J//FM+B31LK4mGUoYCEVVJlF2cSDUiyhE2sxMha682kJqnw14MhxtfG3h6JgFmMssXYyW1sZZ1uqBUxmUxW4sGkEOWb7mHYiocxsFdVnG75EmXliHOxAuSVec2W8FSSeDUWq47eRgy2zKM2ko/sZ4aH4cnOzTE0+didDq4mCZ+BlYR8DJZQivpFjzlZiFKF6uRXvbVBch62WMnZqXNKluMoLkwfKoxna8FDp5TUgg1iuFMSkjHbCGGEHGdtl11YT3yPAf0yexg25M0XxoUtXqyZdolyG5fyRpkkUk57MdDDli4rxzGhWw0KXT/A4rnT2eBbg6GTFtGtfKKcTxzhQT/y2XY7ukzuQP1qLqRFGzEJnmqDIPz8cDyLMMzmFCHzpMdQeuzCk7BZiuReIkKDiAz/FZ3pjhjs+4SERvDT5WgC70Ahp6K41axNlfr1KFWuHI5y3WU3TNEBHF/Zhel7Qrh5w4U3PtjB3CH1KPlUwjjK0lbnqXNgMj1qFJQtgRoSqrK4YvJ/wVnlhXRWBiFnG/G0UiRVJCFKs5C65KlsM7M6MVzZu4OtS6LpdHE6r2qE8Vfsi/aefxkCWXIYyrtQoj2Vx1Bk4ebm9hTCOMPq19YS0XUU/ad2oSZJWOWRcM1qT3bviqXagnak7dvIzVfm8W7XlyjlKJH0Pz5itPddnBu+zoL/FJUbUHa98z7hJwmtpCTEEBf7EkP3rGfIK/f5auh2ol8qRnjwJY58fY+XBw2kU80QDo3dwjmxaMqm6W0daD1lO6PbVKC00yU+GrONPxwb0HbaOzR0Cubyof3s3xRGLc+xNDAdZ127HdjMWsDgvk5cencTVws0pM3YbtR1CufqFzvxmLyVM8qQWzti32EqW0e2oHphG+IlH5Buawyyq7/LybWbObD2Zxz6lufwnm/F4NhTffAq5vSsSwUXMfh2ARxbtAKvXT7clo85V25Gv5nr6V8zlti7J9jQdxmfRcQJYumEUen1yYwZ2ofGBRIx28fj5/Mxa4av56QQjCk1ThLEhXnz/fWMf7MupTN5GOmE8QGblgTQ+OgOBlUQT+WSF93bfU8VdyGMvjUJ8ejNjP0h5O/jhWf/WjgblLEUEk6Mwu96DA7F8+Ns9MHj3eX8UrQvgyYOpHmhOBIe1Dk8izCSEvxIir2WThiZcw/KrZD/lIYnNiaJqPv35JdQIeFIIsLCuXwlmDvBOvI4l6dMWSG2IqUoXboMpUqVsmh/shvmVCMRvj/yh995tg45ivO7Cxk96hXK/ilhnKH6ngl0q1eC9FkDOb58F1d/uUt4pST27zmLzrUyXcYsYrzTAQYv2MOtCBMl31zA/GFdqf+wijiMH3fuZufKeLpqhPEvM4Pacv4qAtkK91S4IDk52eJhqKoplct4NDI8DEUYnkS9MY4hs9pRzfKGYM54rWbblkTquXcidds4FiWMYtvy3jQv7YjvpvYMPVGG5kPnMdi8ieXrv+Z20f70fU1PyLUzHF0fRuNtEpJq4sf2VwbhGVqR1mN70LpyUYqWKEPx/EYCLtwiUrwKO108MTeOsfezsnR370vPJmFsaT6Tc85t6bt1Eu1dfDn9gSees3xpctKTzrYn8Wq7HcOshQwZYOZQk2mcL/E6w3aNptKZlcxbsp2TVt2Y/Z/ykocJ5NR5WwZN703TykXQi5ejzKxJsvg29oH4rF/C7MkniRsygXE1xCCe/5i5PvWZsWk0gxsX4PqW0bjvi8BQtzfdavhy8dBlriZ0YeaHHXAzhnLt+GUxW+K9SEjOPvEce72jKdC6F5OmNcV8YSPeazZz2qoXvdoVgJifOTD/Z8qrUFrvlymXLWH40/DIdiEM8ZBOufPq9Ei6zBjGlDbReL31Lp+ntWTkhvn0KBCN2G7ZWatCA9Hb2IvPI6EyfcBuhrYaxunyM1m0aSFvlYgmNp3NnllWGysx/rh732GwL26p6FKenuzdLV6SKcVIXEw4MVFR4jhZk1fyJRHhYngvBeEXYqZYqbLUqFESB8f8xKeUpFAhNyGNpxPGo2vwLO61lnC//1QGj2+CW44I4w4HRw5hypZ7lJ4zmyFuMYRe3MfGIwkUaNGVfo2L45Bwis2eOppNHsGod6ribDmwRhh/1aho7/v3IpCFMFRrECXcU0M1IMwq3FOEIXmClDOsbT2NjyNtcJLKEydlJKyMklwOJDlvJybteZcqV9czY1wgHbbM461X0/iqWx++KdWbIe+1JXzuDLx/L06ffat4O38K4b/uYGGnQ7hucKfPq0HseW0qnxfoxmivSfQsZ0Nakkk8CogPPc3uqd6cSjZiDLvJqeN5eHvXOqb3MnCk8yzOObbkba+xtHHx4+zWD/hgYQANP3enrSETYfTX82W7OfxQojsj1rclYsFo5n5spNnWj1nWWMIhYpRvnrnC/YBDfHBEjHt4qhhBR6p1GMbwcS8Rvm4OC5f64ea9U7yWfBivbKJ3q8OUXj6fCV3T2Nt9LhervMW0df0lVAcBhxeyYMEJyq45wPB6+XAx+HNsrgf7boQSnxTAxWPBFH59Bms2tCdq7XxWfhhDuwPbGFXeGlP058x72RPj6Bn0G9A4U9JbhaQiuCRe2sKBO7hZtwbF7YXYwm+R2FH+NrQlzR0PM73rx9ws25HpW/pQOS72ico11ZlYPDVjAD9+9wP38lalet1qlLCXZPmDLPizPIyYkCPER1/CsWAruS4cZHMhoTRJtpskdxUVFU5YwFVJA9wlv7O1eJDxnD57m19vpFDarQrNm7lR1CWMqMgwIlOa41rsFUqXKY2TY/YhqfTbUPJBicdY3GAN8YNnMHRcY8rnkDAOjXkXj+8K0+3ABsaWE07+3p3+bXZjXr2PvYPcJIl+gaV1lxDYdiBjFnemkkYY/14LqK0sRwhkm8NQYYTAwEBLPLtgwYJPCPceEcaalu/xtW0pyjZww1XCHDpDHP4/XObu/ToMP7KWVknf8dHgMVzp6cW0xgEsGnmMPD1GsKTPPVa8s5Vj9p14/+BA8U5Cue2zl1W9jlFy43L6NAtid7MVhHcbSZ8pHamWFEtqHhPhl/eza+kOjgaWpXYjZ3RJvlzYHUq9ZWuYOMSe411mcdYhB4RRuidjZ5fk1KSlbAxsytqvptJAcjFS9yQ7axP3Luxlx5e/EByrPKw8lG3YjR5vl+DmqpXs3JpAi4PreL1MAnGBX7KmxRZMk+cxouk1lgxYj0+SC9UbFMVeLG+U73mpUEqkyerLzG54kx88FrL2qBnXOuUpUTCSXz/7CatXxvDegupcXbCKZT/Vx+vMGF42GIn1+4Ll2VZJKcK4x08frcd95EGiurammmrzFHGVwKTmdBjWhzdq/MSKN7dxpUQnIYyBVE+MybbU2SyVVA6O9lJVlSwlrqItUYnyB3mMZxJG6BcY42/hXLS75I6chDCkRFbmSZKS3YCgMAkp+uNsG0ya8Q4Xf7jGD1fiyZu/DC2a16R2dSfRslwkIPA2cbrulCjbkgoVK2Br87BIOJuL+b8ljFvsG76IE3er8s7BSTTVx+N/bger/nOaavuF+GvnxZ7zrG26gBsvD2DQih4W0tc8jBzZFe3N/1IEsq2SCggI4OrVqxQvXpwqVao8IdzLFJJq7kXsO5MYPq0FFSwARXJxozvve4TS+HNv3sh7g192DGLkjf8wqeRBNt+uQ7veE5hR+1tm9trDiTzt8TowgJeIkCqpfXi++QWFhTB6K8J4dRkBbYbQd3oPakvFlCH1MnvnzmX2Div6fb2XWbXtMRu/ZlEtD+LHzGTQKAdOdJzL906t6LVhFK2c73JxtzebZt+k7tM8DAthFOfkhA3sjGzJ5uODcIszyiZWKMPeIF6FlKRaMuoPhtp2W9/mm1Uq6a3Kaj3pVTGJOKnGcW+7C8PUeYyu6sOo4Vs4b1WcZnVc0UtsR58nH4Uq1qNLj46U8n2fQT3Ek1q2EY/RDSmBL5/0nsnxpEb0X1WPG4tWs+znOqw5OZpGstOPvfsNa1tsImXcLPpIWW3mKimVwzgvOQzvpXdp+tU2BqhgfsQ+xtaeTOBbq5kzvSo/j+3HjsgmDFwnuEqoKSYxIyRlwD6vXnJH4kWmqdyMFA88L2HE3SBfkc5CGI5CFklCGknEx8dx7VYId4PuY5UciP+Nc/xxIwyXwuWp36AKL1Vxxdkpmdh7F2Rz4k+cvpsQRgvZpOhFMFrY0mEg+/HfE8b+EQv44k5Feh6YQRvJGfmf2cqy/5yiwkebGNLQmbycY3WTJdxpNIDB7q/L9akRxr/U/mnLyiECWZLeqvmg8i5US2sl2nN1dX1CuJc5h7GWe11G0G9KN0l6qyqpRxqFZp9JyWmJCKIvetN55A9Ym36g5NCljOgp+YwiPnj1WS6eQk0mfDufVlZR+H++QkokL9Fw0yqGNLvLh5kJw1aqrcK/wnPqWtx/bcyOy7NpLrVSxste9O56ANepi5kyohCn3hwsJNSZvh+8S0vDVY4unMOyHfa8c9id9oYTeGbkMB4LSbWU8NgcVnwBHfeuZWIlSY2mRnHjolTzlC5GPqc8cu4Z5bJP6jAyE8aHFsIY0SIUj94r+KXBEDzW9JFyYzVSSAy6I34UxB11p/PkRIYdXcXEhoWxST7Ooi5L+c6xG/N2diRq2Vw8PpT8ybmNvFlQqnjUGrt9SsHJCxnXryFln8hhZBBGky+2MVDFZgK20bfhZELe8GL+4p4UEUIes+488e1WsW1aO0pZKlaVaC+Qi98EU7BmGYo4RPLLKRWSqkzVWlUltJWDkFTolyTFK8LoKiEpRRgSkhIvIz4ujqvXA7l6LYBAv+sEB1ynQIH8NGlci2qVC0nYSarn0qKJCr9AkBBGNJ3BUFFae3xDx44deO211/7kUv6eFbWWEjd6AWMG10AyPYwfNw5/2eg8UnqrslpVJfVk0vsWGmHk0Epob9cQeIBAFg9DhaEspZVS3qnammd9FkYGYZxiaf3lhPcYz5DZb1DnQVntN0tXsHl9DC2+9KZPJSmDDTsleo3ezLndhFkfuTOufQUc04I47zEVzx23ML09hbfLBnPj9EHWbk2m975NjGnmh3fd+fh3HMWguW9R30pKQbmLz+a1eHhfwWXAMNoUlqT3ddEoLL1Jk5XrmTu8NjFrerL8Ox2OLfrT0fkXvtyxh4O/1WTasbV0NnzFysabsVm4glFD9RxuMvVR0ttnOXMk6f2tdQ/m9nGzJL2//d5A/0lvUL9iYXRSapke0rfBMd9dMtb42hfe9KuSROztAyxo5I31e8sYOaoeyduGsPiT+5L0fovudR2kBcbv3I1Opfzg4TQMP8Lq0R9wp2F32jSWGH7CadbNOUJU3ZGs2T+MIgfn4TH3ELfaz6FfjXjibx5l8dIbNFm+npkDFGHEPhDuKfJSHoaQ6JhvcZo7kTbFRLtyYQeLTxWh45RpTOxRnZL3fPBeOZNVZ2x5ueto2pYxopOy2vvhv3POpwYDZ3WgkfMXDG80kG/LzWSxJL175iDpHS05jESVw3BtLeW6knsQD0NqiqVxZYK0mAnHLzCE4LsBUkSRQOkShahUQfJdjnIN6ZReJobYCMkRBYVwI7QhFy/HiMHfz+LFixgwYECWm9ScHC+q92P8cF2qpEZ8irF9Lzr2akv7V+rguWwaATfvsP9ha5CnEcYN9vSfyaFblen91Xw65InD99RG5r9+nEoHdzGqqQsOnGZprbncajqMket6UlvzMDSDqSFgQSALYdhLBZIijJiYGAtZqAT4k1VSemvROyRfYuegXUQ270nXAS2oaCmrDefCjt0c/jSeuqun0bakHbqoO3yzqAVTIkeyeMIwer5kS3SyLU73z7Bn01rm7/1N2kKI2rpyHVpeN1F58WTaSPns0QFbCGnUg86DW1EZEdbJ7lUfJJVUW+bz3r4AzPKgJ+tXmtDiciyVJ4yldceauN39jCXu69n5XSCOdWpStUR1Sv2WSI0143nF+gJ7hh3CeugYenTVcWbEen4r1IxOU3vSwCmU345uZ/WkbZxV6Qo5ln37aWwbpcpqDcRJwj19SOjG4dEa66yeQYfSEvcP+o6tQw9jNWAEXbvUwi3Pdb5YsIJ1UlYrMgMJY5WgQpcxrJzZnhKpkdw5OJsJXqfxj5aDNWnIS75m6tVrwcszpSor6Vd89ngwzPOclNhKXqFh0/Q1jh9Dq/Y1KGWWNhcW9pKQkkM014/tZt1Eb05L3skiIZGeWS0mbWds55qUt48lQSLy+cx/8NUmTxZu/I7wjJU4FKT99A8Z27I4ReOOsGScO78U68fgSYN5LQdltVHBh4kM+kyU1CUtVVJmqY5SmpY0SXynpIgWRTQnKclSLiweoZ0II21tlbJFmrGIeFASGyTESvjzppRRnyvC1yf9pB/UHby9vRk4cGCWW9QUE8gJjx7MOxBPVKrkWaR6L61gO6aunsNPBxYReOU6+54kjNYXqHtwKm9UdXkgtJOS54Ub8QksQ5s1g2lkm0jwz5+xZfoPlBbvrkdN+e65zM5+m7lbsxvdJ7aUa1uNeK5+sh3vRZF01MpqNfOZSxHINukdJWWQKoeh4shZdRgZSEksWURPsq18KJSziKGUIltUaRbhm50OY+gFdo4cyuWWXgzv1UY0DzFSQillnWIwUsWgJEovpzQRjOlFka1qs3QiYNOLkM6UZW6ViFX9gcTIGEUsppoHSh8lFWHRiTBOL7FvtWtNlsouJZgzixjLIH+3VmI4SaLqlaZBzSl/Uw+NS5Of1XGVctxKieXEwCUnJMruXSmh5Xxs7SR/ka4qf7zrR+Y1PlBsi5FMVVoNy9xK/Ca9tJKkiitDoCfiRWsbmc9W1NQKPqUeThTxnlKoyxpsJAxoLeputXalQlZCuATJLVgaJD62xqwNEy2CwyTxLDKU4jKPXZ682Cj1nSpxVWtUy5Oqsng5ZnoLE9WnSi89qPIi+rP0nlop6WpvpcPJiXAvLuI894MOCO5C6kqUKcJEVTSR3l1LsLQcL/3n9J+USDLdX1OEEhDgx4nTNzh6MplbfpJgke9py5Yt2XoYFsW8CDqNcv3oRFQpBxKSFBISZfiUSRPx9/N7IiR1mAWNj1Bk3SDavlSUgi7qpUKMslbVh8tG1qrOSdafKuJQnUUcqr4hdRypjFOdAnRGoiMiiIgN5tf9hzi43ZY+l2fxmibcy6UmM3cvO0sOQ3kXSrinEt9KuJdV6Z0BmKhhbZTlTW8nYWlgp+r7rVVIS/1ZbkhzJMHfb2XC0Gs08ppC73aVcYlNkIi+eqt8XrQABnXjq+Z96ua32G0Jg6kWIDK3LvPcFmGYavUhRKAUuA8/I21DpC1HmuQZ0nRqPtWSQ7UfUcZEGUz5WbQGFgW4EICcrHhOwnNKWS1/tRzvgS7BWgyI5QGklrlF3KbW9ThbPL5Gy2eV/VXtPjLmVsdVx1JtTx60MVE7YcuOOz0XopTFBvVvyjYpQ69anyhDq+azKNWlpYr8e9Y1PtHk0AKJHCfjvB+cu8nSLuVRU0FL2xSZUx0zI4WvvEbVViXF0lZFmhMqUjWL4VSaE8EjYzyrSkoJC1OMIdI+I17I3yiGV3kXqqWKOrn0WR4eM2NS+TejEPstuc6++vprjhw9zh3/uIfH3Lx5M4MGDcrRnTlm9Bj8/P04dOjQg8+lEh9ykAUSevxYiEwn30/94RuYOaw1NXLSSyrsBBuXubP84HVSEu1xrTCYVScn0DQn7UVytBLtzRoC/7sIZCvcUw9OyghHKV3G4yGpv7YYvXT2TAu/wldr3uHdWwNYP78fHao6SvWM7P7/2hTau/4HEHgWYWScokkISr0sxKcINwvTpr9TL+SqiNzP1499+z7ho4/2cO36NYtnkjGeFpL6MzjGqaS3tDdXTwZMH4oQYwi7HUBIrPTUknidQ1E3ShV1IW9OHkueEiUJe0ncR0i1m/Q9y+tclFJlC0iRtTY0BHIfAlkIQxkI9VK7RCXce97Hs6oGemnGWELvXMLX1o0axaWMUlqBqxtXG/8cBDITRt++fS1Vcy9iJEj4LygoiLCwMEvZtj5THEw9Rz6nT/a7ffu25XpVZeDa0BDQEPh7EMg2h6HCEr6+vvJcC8OfPA/jGSekdpiS37Cxd8Be4sAqfi6RgcfaDf09S9JmfZEIZCYM9VS7fPnyvcjptbk0BDQE/kEIZMlhKAOhGg+qbrWlS5f+f3ni3j8Ir1x3qup6UE9g3Lp1q+V6eFpjwFwHjLZgDYFciIDuzJkz5rp161qSneqlhHvqed6qrLaQPGPZxUVaXz+oasmF+OT6JasiCNXu/scffyQ6OtpyfWhDQ0BDIHcioPPx8THXq1fPQhYZjQYzhHsq8a2RRe68MDKvWpGE0uNoZKFdCxoCuRuBbAkjd0OirV5DQENAQ0BDIDsELCGpBg0aPOZhaFBpCGgIaAhoCGgIPImA7vvvvzfXqFHjsRyGBpOGgIaAhoCGgIZAFsLYv3+/WbWDeHU3vkUAAABVSURBVB5xnganhoCGgIaAhkDuQUDXvn176aKhOjNpQ0NAQ0BDQENAQ+DpCOh+/uknTXqtXSEaAhoCGgIaAs9EQHfm2j2NMJ4Jk/YGDQENAQ0BDYH/A4brPZJEut3JAAAAAElFTkSuQmCC"><br> <span style="font-size:11pt"><span style="line-height:107%"><span style="font-family:"Calibri",sans-serif"> there is no indication of a flow rate but only a head value that I can change. So I cannot define the flow rate. I’ve checked with the Inspection after entering a value in this above window and there is no constraint defined.</span></span></span><br> <span style="font-size:11pt"><span style="line-height:107%"><span style="font-family:"Calibri",sans-serif">So, how can I define this flow rate?</span></span></span><br> <span style="font-size:11pt"><span style="line-height:107%"><span style="font-family:"Calibri",sans-serif">And is it realy that my problem stems from the fact that there are no constraints in the BCs or is there something else?</span></span></span><br> <span style="font-size:11pt"><span style="line-height:107%"><span style="font-family:"Calibri",sans-serif">One last question: when FEFLOW finds that the error is below the error tolerance and ends the simulation, this error below the errro tolerance is the average of errors in all the nodes between one simulation and the next one or is the highest error between simulations in at least one node of the model?</span></span></span><br> -
My panel Data has no option for inactive cells
My panel Data hás no option to define inactive cells. There shoukd be an option Activew/inactive cells but there isn't. How can I define inactive cells in this situation? -
How to import boundary conditions from a previous model
I have a regional model that is an expansion of a earlier one. There are slight changes in the node distribution of the mesh, once I had to add the new areas. I want to define boundary conditions in this new model using the boundary conditions of the former one. So I exported the BC in question (head constant along rivers) from the earlier model using the export Data and Current slice (where this BC was defined). I opened the exported file in a SIG and, as expected, it was a point shapefile with all the points along the rivers, as it should be. Aftewards I imported this file into the panel Maps and linked Z (the height of the terrain) to head constant BC, as in the former model. I mantained the default snap distance (2091,61 m). Then I went into the panel Data, selected head constant BC and when did the procedure to apply the BC along the nodes, I found that it did not defined the BC as it was in the shapefile – that is along the direction of the rivers – but applied it to all the nodes of the area I selected. The way I found to overcome this was to select node by node, for it to apply the BC in the correct nodes. Once this BC has to be applied to the whole river network, so theere is a lot of nodes, I would like to know how can I import the BC shapefile in a way that just the rightful nodes are defined, instead of all the nodes of the model. I have the same problem with the multilayer wells shapefile -
delauney violations
I would like to know what is more important for an efficient performance of a model: eliminate all Delauney violations or equilateralise the triangles? Once I have to do this manualy I have no time to do both, so which should I choose?
By the way, is there a way to autmatically eliminate all Delauney violations? Smoothing the mesh helps but is far from elimminate all violations.
Hoping to ear from you, Emília Novo -
error message PCG has a matrix with 13127th diagonal element <=0
I am running my model, and I got the message "PCG has a matrix with 13127th diagonal element <=0". The model seems run and seems to converge but does not calcularte the piezometric heads. I defined the conductivity, head-constant BCs, recharge on inflow/outflow on top/bottom as 300 mm/year and the remaining material properties I used the default values.
I send a file with some more info.
How can I overcome this and get the model calculating the heads? -
RE: disapering mesh
I did hat you suggest and the problem remains.
-
How to import and link a shapefile of boundary conditions?
<span style="font-size:11pt"><span style="line-height:107%"><span style="font-family:"Calibri",sans-serif">I have a regional model that is an expansion of a earlier one. There are slight changes in the node distribution of the mesh, once I had to add the new areas.</span></span></span><br> <span style="font-size:11pt"><span style="line-height:107%"><span style="font-family:"Calibri",sans-serif">I want to define boundary conditions in this new model using the boundary conditions of the former one. So I exported the BC in question (head constant along rivers) from the earlier model using the export Data > Current slice (where this BC was defined). Just to check, I opened the exported file in a SIG and, as expected, it was a point shapefile with all the points along the rivers, as it should be.</span></span></span><br> <span style="font-size:11pt"><span style="line-height:107%"><span style="font-family:"Calibri",sans-serif">Aftewards I imported this file into the panel Maps and linked Z (the height of the terrain) to head constant BC, as in the former model. I mantained the default snap distance (2091,61 m).</span></span></span><br> <span style="font-size:11pt"><span style="line-height:107%"><span style="font-family:"Calibri",sans-serif">Then I went into the panel Data, selected head constant BC and when did the procedure to apply the BC along the nodes, I found that it did not defined the BC as it was in the shapefile – that is along the direction of the rivers – but applied it to all the nodes of the area I selected. The way I found to overcome this was to select node by node, for it to apply the BC in the correct nodes.</span></span></span><br> <span style="font-size:11pt"><span style="line-height:107%"><span style="font-family:"Calibri",sans-serif">Once this BC has to be applied to the whole river network, so theere is a lot of nodes, I would like to know how can I import the BC shapefile in a way that just the rightful nodes are defined, instead of all the nodes of the model.</span></span></span><br> <span style="font-size:11pt"><span style="line-height:107%"><span style="font-family:"Calibri",sans-serif">I have the same exact problem with the Multi-layer wells BC.</span></span></span><br> -
How to add new areas to a model
I have a structured mesh model already developed and I need to add a neighboring area to that model. Is that possible? If so, How it is done, can you guide me to documentation that explains how it is done?
Hopping you can answer, Emília Novo -
how to correct Delauney criterion violations in n automatic fashion?
i have a very large mesh with lots of Delauney criterion violations and editing them manually is an highly time-consuming task. So, I would like to know if there is a way to automatically correct Delauney violations and if so, how it is done.
And after the Delauney violation corrections, how to automatically perform the equilateralization of the triangles?