Hi,
I'm modeling fluid flow through a 2D porous media using Darcy's Law. I intend to estimate fluid flux values across the top boundary of model, however,my desired outputs will be one in which the results obtained at this boundary are dependent on the domain solution.
For my simulations so far, I have observed that the domain solution depends on the prescribed values of the boundary condition, i.e I have to specify Dirichlet boundary conditions of stress, pressure or velocity along the top and bottom boundaries. The Side boundaries are specified as No- flow.
The Comsol user guide suggested that the use of Open Boundaries ( with No viscous stress) might be useful in obtaining solutions in cases where pressure cannot be predicated ahead of time, as such I am trying to run the model using the Brinkman's Interface in order to take advantage of the open boundary physics for the top boundary, along with a specified condition along the bottom of the model. However, I have been unable to get my model to converge with this specification. (The mph file and model diagram file is attached)
I do not understand why this is happening and would appreciate insights/tips into this problem, either under the original Darcy's Law formulation or using the Brinkman's equations.
Thanks.
I'm modeling fluid flow through a 2D porous media using Darcy's Law. I intend to estimate fluid flux values across the top boundary of model, however,my desired outputs will be one in which the results obtained at this boundary are dependent on the domain solution.
For my simulations so far, I have observed that the domain solution depends on the prescribed values of the boundary condition, i.e I have to specify Dirichlet boundary conditions of stress, pressure or velocity along the top and bottom boundaries. The Side boundaries are specified as No- flow.
The Comsol user guide suggested that the use of Open Boundaries ( with No viscous stress) might be useful in obtaining solutions in cases where pressure cannot be predicated ahead of time, as such I am trying to run the model using the Brinkman's Interface in order to take advantage of the open boundary physics for the top boundary, along with a specified condition along the bottom of the model. However, I have been unable to get my model to converge with this specification. (The mph file and model diagram file is attached)
I do not understand why this is happening and would appreciate insights/tips into this problem, either under the original Darcy's Law formulation or using the Brinkman's equations.
Thanks.