Comsol-Pressure drop fluid flow in rectangular channel problem

AI Thread Summary
A user is attempting to model pressure drop in a rectangular channel using COMSOL but encounters a solver error when setting the inlet velocity to 0.5 m/s, while lower velocities work fine. Suggestions include ramping up the velocity gradually and adjusting boundary conditions to improve solver stability. The user successfully resolves the issue by switching to a time-dependent model, which allows for a smoother transition in flow conditions. The discussion highlights the importance of avoiding sudden changes in flow to prevent solver instability. Overall, proper modeling techniques can significantly enhance the performance of simulations in COMSOL.
baranii2000
Messages
4
Reaction score
0
Dear all,

I am trying to solve a simple model with Comsol to find the pressure drop in a rectangular channel W=5mm, h=0.6mm, L=10mm. When i try to set the inlet velocity to 0.5m/s the solver shows the error message "Maximum number of Newton iterations reached". The solver works fine if i set the input velocity to lower velocity(eg., 0.005m/s). when i increase the velocity to 0.05 m/s the solver ends up with the error message. I even tried to increase the height of the channel, still get the same error sequence.

I am using a Stationary, Direct (Pardiso)solver.

Any comments or help is much appreciated.

Thanks in advance

Arun
 
Engineering news on Phys.org
I'm not familiar with COMSOL, but here are a couple things you can try. First, try ramping the velocity up. Solvers hate suddenly seeing a flow.

Also, you can try changing the boundary condition to total pressure at the inlet, and static at the outlet for a first cut. These tend to be more well posed than others and can have better stability.
 
Thanks Minger,

Actually your comment on the solver crash on seeing sudden high velocity at inlet was very useful. I remodelled as a time dependent file and now it works fine...

Many thanks

Arun
 
Right, just to clarify a little more on what happens. Unless you normalize your mean flow at t=0 to some value, at that time 0, you suddenly impose a velocity at the inlet. One point inwards your trying to either take a derivative or calculate an integral with a discontinous flow.

You will calculate an artificially high flux at that point, which will then give you an artificially high time derivative. This will just continue until your solution blows up.
 
Thread 'Physics of Stretch: What pressure does a band apply on a cylinder?'
Scenario 1 (figure 1) A continuous loop of elastic material is stretched around two metal bars. The top bar is attached to a load cell that reads force. The lower bar can be moved downwards to stretch the elastic material. The lower bar is moved downwards until the two bars are 1190mm apart, stretching the elastic material. The bars are 5mm thick, so the total internal loop length is 1200mm (1190mm + 5mm + 5mm). At this level of stretch, the load cell reads 45N tensile force. Key numbers...
I'm trying to decide what size and type of galvanized steel I need for 2 cantilever extensions. The cantilever is 5 ft. The space between the two cantilever arms is a 17 ft Gap the center 7 ft of the 17 ft Gap we'll need to Bear approximately 17,000 lb spread evenly from the front of the cantilever to the back of the cantilever over 5 ft. I will put support beams across these cantilever arms to support the load evenly
Back
Top