Jump to content
GMS, SMS, and WMS User Forum
Sign in to follow this  
Paul

FESWMS Sediment Transport

Recommended Posts

Hi:

Been using FESWMS and the hydrodynamics run great! However, sediment transport is not good at all. I get the RowCount > MaxFrontWidth error. I have a simple mesh, mostly rectangles, yet it stops when assembling elements on semi-coupled calculations. I have used all the functions and receive the same error. I know it has to do with the .upp matrix file, burt the program deletes these after the erroneous run, so I cannot check into tyhe problem here. There is nothing special about the element that causes the error except that it is next to a dry node. Please help me with this issue. I looked in info, and my max front width is 133, however, the sediment calculations in the .prt file say it is 103. This is less than the info SMS is giving. Is there a way to change the max front width to accomidate for this.

Any help or suggestions would be greatly appreciated. Any other users have this error and found out how to fix it or be able to provide more answers. I am stumped.

Thanks

Share this post


Link to post
Share on other sites
Hi:

Been using FESWMS and the hydrodynamics run great! However, sediment transport is not good at all. I get the RowCount > MaxFrontWidth error. I have a simple mesh, mostly rectangles, yet it stops when assembling elements on semi-coupled calculations. I have used all the functions and receive the same error. I know it has to do with the .upp matrix file, burt the program deletes these after the erroneous run, so I cannot check into tyhe problem here. There is nothing special about the element that causes the error except that it is next to a dry node. Please help me with this issue. I looked in info, and my max front width is 133, however, the sediment calculations in the .prt file say it is 103. This is less than the info SMS is giving. Is there a way to change the max front width to accomidate for this.

Any help or suggestions would be greatly appreciated. Any other users have this error and found out how to fix it or be able to provide more answers. I am stumped.

Thanks

Hi all:

After lots of research and a great deal of testing, I finally have a stable sediment transport run. I will post and faq with the detailed procedures on how I was able to get it to run for others who may have had similar problems. I should be able to get something posted later today.

I hope anyone else with problems will find this helpful and if you have questions for me please reply to this topic and I should be able to get your model going in the right direction.

Best,

Paul

Share this post


Link to post
Share on other sites
Hi all:

After lots of research and a great deal of testing, I finally have a stable sediment transport run. I will post and faq with the detailed procedures on how I was able to get it to run for others who may have had similar problems. I should be able to get something posted later today.

I hope anyone else with problems will find this helpful and if you have questions for me please reply to this topic and I should be able to get your model going in the right direction.

Best,

Paul

As promised, here is the procedure I used. Use with caution as this probably will not work for all scenarios. It is a good step since the documentation on this is quite lacking.

Procedures for getting sediment transport to work:

1. Make sure the mesh is very stable when running hydrodynamics. i.e. its important to have as small of an elevation change as possible across elements. In addition, elements near the location of boundary conditions, especially the upstream one, are prone to cause instabilities. Here, it’s a good idea to use as small elements as possible with a maximum of 1-foot in elevation change across each element. Also, element area changes greater than 2 or 0.5 can cause instabilities elsewhere. This is not as vital as boundary elements, but it helps stabilize the model. If the steering module reaches a solution in 5 or less runs, then the mesh is ready for testing the application of sediment transport. In addition, if the velocities are greater than 15 fps anywhere, your model is not stable enough. Mine had a maximum velocity of about 7 fps, but I imagine larger velocities up to supercritical flow will also work.

2. With a strong mesh, enter your bed control parameters in the model control sediment control tab for feswms. Here enter the particle distributions with the sum of each percentage equaling one. Use 1-8 particle sizes. More particle sizes does not create instabilities, it only increases computational time. Also, active bed layer, deposition depth, and original depth are very important to stable sediment runs. The active bed layer should be approximately 2 times the median particle size (50-percent finer) of the study reach. Deposition layer of 1-3 feet was stable for all of my runs. The original bed layer is only stable if the depth is equal to or greater than the bed elevation of your upstream boundary condition. For example, if your upstream bed elevation is 600 feet, use 601 feet for the bed thickness. Values lower than this will cause your model to fail.

3. Next, the sediment transport function (under parameters) should be chosen based upon hydraulic conditions of the reach. I have used Engelund-Hansen, Laursen and Yang’s sand and gravel with reasonable results. Since Yang’s equation was the best method for my hydraulic geometry (fine sand to coarse gravel bed), I used this.

4. Still under sediment control, I used a single iteration for each sediment calculation. Each time I used more than one, the model became unstable more quickly. Balance diffusion was also check, with bed convergence parameters set to 0.1 and 0.1. More stringent convergence criteria may be possible, but for my purpose, this was adequate.

5. Under the general tab in model control, semi coupled simulation was used with the initial conditions file being the latest hydrodynamic solution produced by the steering module. This is important, because for each successive run under semi coupled, the initial conditions will be updated and allowed to achieve final convergence. Also, set the run to dynamic, not steady state.

6. Under parameter times, I used 15 iterations for the hydrodymamics, with a 6.0 hour run at a time step on 0.1 hours. For sediment transport larger time steps are necessary to show bed evolution. Time steps are not important since the solution method is implicit meaning the Courant time step is not applicable. The time integration factor was left at the default.

7. Exit the model control parameters and check the material properties under the FESWMS menu. Manning’s values are arbitrary and work for all simulations. However, under turbulence parameters, artificially high values are necessary to ensure the hydrodynamics are as stable as possible.

8. For the upstream boundary condition use a constant flow with sediment options set to equilibrium conditions applied and initial bed elevation fixed. I tried clear water and it worked, but more realistic results occurred with equilibrium conditions. For the downstream boundary condition, equilibrium rates are also applied.

9.The model is now ready to run. Go to FESWMS menu and select run simulation. I always received the bad aspect triangles, but I ignored it and ran the model. Wait a bit and the model will hit an error (for me it was around 0.5 hours). This is not a problem, since the *.FLO file as the initial conditions will now have 5 times steps in it. Just rerun the model and it should progress further. Mine reached approximately 2.0 hours the second time, 3.5 hours the third, etc, until it finally ran to completion. With a complete *.SDI file, use this as the initial sediment condition if the sediment transport function is changed.

10. It is also recommended for maximum stability that you have zero dry nodes at all time steps. This is very important especially at the boundaries. Every time I had even one dry node, the sediment transport run became unstable.

I hope this helps for anyone else who was as frustrated as I was.

Good luck,

Paul

Share this post


Link to post
Share on other sites

Hi Paul

I'm a begginer in SMS and I have interest in purchasing the 10.0 SMS program with the Finite Element Riverin Package. The software with the modulo map, mesh, scatter, RMA2, RMA4 and FESWMS.

The application of the model will be in hydraulic and transport of sediments in reservoirs. This package would have the necessary modules for this application or will be necessary another type?

Best regards,

Denis

Share this post


Link to post
Share on other sites

Aquaveo employees have successfully run Sediment Transport test cases using FESWMS and the following options:

* Equilibrium inflow

* Clear water inflow

We have been unsuccessful in our attempts to create test cases using the other options.

See:

FESWMS known issues

Ipson, Mark K. (2006). Analysis of the Sediment Transport Capabilities of FESWMS FST2DH. Thesis, Brigham Young University.

Share this post


Link to post
Share on other sites

Join the conversation

You can post now and register later. If you have an account, sign in now to post with your account.

Guest
Reply to this topic...

×   Pasted as rich text.   Paste as plain text instead

  Only 75 emoji are allowed.

×   Your link has been automatically embedded.   Display as a link instead

×   Your previous content has been restored.   Clear editor

×   You cannot paste images directly. Upload or insert images from URL.

Loading...
Sign in to follow this  

×
×
  • Create New...