CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

COMSOL: optimization applied to a laminar flow

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   April 30, 2012, 08:57
Default COMSOL: optimization applied to a laminar flow
New Member
Join Date: Apr 2012
Posts: 2
Rep Power: 0
tdh89 is on a distinguished road
Hi all,i'm new here and i really hope that you could help me with this problem that is blocking the progression of my graduation thesis.

I'm working on a 3D stationary simulation of brachiocefalic artery; i'm using laminar flow physics with a known velocity in inlet and a normal stress in outlet. The model geometry is a simple cylinder.

For the normal stress i'm using this formula:

S=Q*C+po (where S= normal stress; Q= flow rate; po= outlet pressure (in this case venous pressure; C= resistance coefficient)

All the data are known, except the coeffcient C that is the parameter which i want to find: so i'm using the OPTIMIZATION module.

My objective is to find the C value for which the average pressure (p) in the domain is near to a reference pressure (called Pref), so i used these settings:

INTEGRAL OBJECTIVE (applied to domain):
Objective expression: abs(p-Pref)

CONTROL VARIABLE FIELD (applied to outlet boundary) :
Control variable name: C
Initial value: 0

All other are default settings.

The problem is that the solver returns always the initial value as optimal value, even if the resultant pressure is very far from the reference one.

I tried to change mesh, initial value for variable C, i used also probe objective instead of integral objective, but the problem persists: optimization solver terminate WITHOUT errors, returning back whatever initial value used for C as the optimal value.

Pref=15000Pa (so p should be near to this value)
C initial value = 1000

Log of simulation:
Optimization Solver 1 in Solver 1 started at 28-apr-2012 19:02:53.
Optimization solver (SNOPT)
Itns Major Minor Step nPDE Error Objective
0 0 0 - 1 2.48e-011 0.3085
Number of optimization variables: 75.
Number of objective function evaluations: 3.
Number of Jacobian evaluations: 2.
Final objective function value: 0.3085111759.
Optimality conditions satisfied.
Optimization Solver 1 in Solver 1: Solution time: 110 s. (1 minute, 50 seconds)

p average on the domain= 1460 Pa ( instead of 15000!)

What's the problem?

P.S.Tell me if you need more informations or data.

tdh89 is offline   Reply With Quote


cfd, comsol 4, laminar flow, optimization

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
boundary layer specifications for laminar flow ranjith CFX 7 October 16, 2017 08:09
Boundary Layer of Laminar Flow over a Flat Plate Blasius_Pohlhausen_Crocco Main CFD Forum 12 September 30, 2013 17:35
pressure distribution in water flow, differences in icoFoam and COMSOL deniggo OpenFOAM Running, Solving & CFD 14 September 30, 2010 03:48
Laminar flow mech FLUENT 4 February 1, 2007 15:18
Laminar flow or Turbulent flow mech FLUENT 0 January 27, 2007 19:51

All times are GMT -4. The time now is 04:50.