Quantcast
Channel: Solving, solving\ - Latest Discussions - COMSOL Forums
Viewing all 3988 articles
Browse latest View live

Heat transfer

$
0
0
Can anyone tell me what is the recomeded NLfail value in simulation Comsol solver.

Thank you inadvance for your help!

using the Amazon cloud

$
0
0
I was interested in hearing any opinions regarding the use of COMSOL in the Amazon cloud. We currently have limited computing resources in house and so we'll either purchase a new machine or switch over to using Amazon EC2. I have heard that Amazon can be really useful when the need is for LOTS of RAM or when performing parametric sweeps; I also heard that it is relatively slow (like a 10 minute computation performed locally might take hours when operating in the cloud). I'm interested in hearing anything about your experience using COMSOL in the Amazon cloud, including (but not limited to) how easy was it to setup/use COMSOL in the Amazon cloud? and what has been your experience with performance?

Thoughts/opinions regarding your experience with Amazon cloud is very much appreciated!

Best regards,

John

results depend on parameter sweep order?

$
0
0
Hi all,
Is it possible that in some cases the results of a structural mechanics simulation depends on the order in which a parametric sweep is done?

I am simulating the deformation of an object when it "hangs" from its support-points. The support can be tilted wrt. vertical about +/- 2 degrees.
So I am running a simulation where the body-force depends on the tilt-angle = range(-2,0.1,+2)
From this I would expect to get symmetric results, i.e. the -2 degrees tilt situation should look exactly symmetrical compared to the +2 degrees tilt situation.
However I seem to get "skewed" results and the -2 and +2 degrees tilt results do not look symmetrical.

What I think is happening is that comsol solves the FEA problem for one case, the first parameter value of -2 degrees, and then all other solutions are small (linear?) modifications to this original solution.
I have yet to try what happens if I run the same simulation with tilt-angle = range(+2,-0.1,-2), or with individual parameter values of -2, 0, +2 degrees of tilt - allways erasing the previous solution between runs.

My model consists of a piece of glass, and the deformations are small, 1e-12 to 1e-10 displacements, so there might be some numerics that should be tuned also..

thanks,
Anders

Ramp fuction or?

$
0
0
Hi,

I am simulating heat transfer with fluid flow, thus I need to slowly increase my heat source to make my simulation works. Until now I've been using step function but I feel liek ramp function fits it better. Id like to have my heat source value increasing from 0 to 100% in 20% of the time of simulation, then to stay at this level.
I'm not sure if I am doing correctly. For example, I have my heat source 1000 kw/m3 and 300 s time simulation. Should I set slope as: 1000/0.2*300 and cutoff: 0.2*300 ?
And then set my heat source as 1000*rm1(t[1/s])?

I appreciate any help cause I'm little confused now.

new to cluster computing, need help

$
0
0
Hi Elham,

My name is Ehsan a PhD student at University of Waterloo. I am trying to learn cluster computing with COMSOL and unfortunately nobody has done it here!
We have the right licence and we do have a server with windows server 2003 but apparently there are things that I don't get at all.
So I was wondering if I can ask you some questions regarding this issue.

Please let met me know if we can talk on phone.

Thanks,
Ehsan Asadi

PhD Candidate, Research Assistant
Mechanical Engineering Department
University of Waterloo
tel: +1 519 5025876
E-mail:
esn.asadi@gmail.com
easadi@uwaterloo.ca
www.linkedin.com/in/asadiehsan

Weak Form PDE: how to impose Neumann BCs?

$
0
0
Hi everyone,
I'm new in Comsol, so I apologise in advance if my question is trivial. I'm exploring Comsol capabilities in solving PDEs, and my current unidimensional equation is the follow:

A(u) u'' = f

where u = u(x). I'm trying to use the weak form, then by using a test function v, it becomes:

[v A ux] - (vx A ux + v A' ux^2)- v f = 0
-test(ux)*A(u)*ux - test(u)*dA(u)*ux^2 - test(u)*f(x) in Comsol

where [ ] denotes the boundary terms, and the apex ' is a derivative wrt the argument of the function (i.e. A' = dA/du = dA). Here (A ux) should be the flux term.
In order to check if what I do is right, I'm using two benchmark case:
- CASE (1): A(u) = u and u(x) = 2 sin(x) + (x/pi)^2 + 0.5
- CASE (2): A(u) = u and u(x) = x^2-x^3 then f = (2-6x)(x^2 - x^3)

By imposing Dirichlet BCs (by means of "Dirichlet Boundary Condition" modules and selecting "prescribed value of u") Comsol computes a solution that (1) does not match the benchmark function, even if the BCs are satisfied. In case (2) it finds the benchmark solution.

By imposing Neumann BCs, instead I have some trouble. A "flux/source" module is added, and then the equation is

-n flux = g - qu

I set q = 0. At the right side, for instance, I compute g as follows (R denotes right side, i.e x = xR):

flux(R) = A(uR) ux(R) = uR*duR = g

where uR denotes the value of u at x=xR whereas duR is the value of derivative at x = xR.
When I compute a solution using these settings (for instance a Dirichlet BC at the left side and a Neumann at right), in the case (1) Comsol reaches the maximum number of iterations allowed (even if I use as initial guess the exact function), whereas in case (2) the final solution breaks the BCs I wanted to impose.

Where is the mistake?
Thanks in advance for all your help.
Pietro

PS
In which direction is the normal vector n assumed to be positive? When does it point inside or outside the domain?

Ice Slurry flow in pipe

$
0
0
Dear Nagi Elabbasi

I am currently trying to build a model for ice slurry flow in a pipe using mixture model laminar flow, with my interest in velocity, dispersed phase fraction distribution. i initially tried building the my model based on the dense suspension benchmark in-built model in comsol but i was came up with inconsistent initial values problem
so i tried using a model based on sedimentation model which was applied in the modelling of contaminant- removal from wastewater in secondary clarifier, and i "Failed to find consistent initial values. Last time step is not converged" keeps showing

Pls i would be grateful for any help

Thank you

Finding the best step size

$
0
0
Hi!

I've been trying to simulate flow past a (hot) cylinder. I simulate Re-Numbers up to 190.

After running the simulation for various step sizes (0.5/1/10/20/50/100 s) I tried to visualize the relative errors of some specific values with respect to the value acquired with step size 0.5s.

The errors cocnerning temperature and v_max are insignificant - or at least not too bad (max is about 15%), but the errors concerning the velocity at some points are much worse (max is 95%).

I noticed different things: The errors for step sizes from 10s to 100s are almost the same, no big difference and the errors for step size 1s is sometimes even bigger than the ones for higher step sizes. Is it possible, that the step size 0.5s is too small? And is there a minimum step size?

The simulation time is 12400s and the water tube is about 1.20m in length. It may also be a problem that the vortexes that occur behind the cylinder develop at different times by using different step sizes and that I just cannot evaluate the same points at the same simulation times. But if that's the case, wouldn't the errors for the step sizes 10s to 100s be very different, too?

Best regards
Alexander Spranger

Joining Edges (results)

$
0
0
Hello All,

I have a situation where I have a parametric sweep, and I get a couple of solutions (all having the same geometry). I can use the join operator to subtract the solutions and I get something that makes sense. The issue is when I try to apply the same approach to one of the boundaries. I define it via "edge 2D" and I can plot the any variable over the arc length of the edge, however when I try to subtract two edges (via join), they don't even show up in the "data" drop down.

Does any one have any insight into why this isn't working, or how to go about subtracting 2 line plots (asides from exporting them)?

Thanks,
Aaron R. Shifman

How can I manage individual cell values of a mesh in a convergence process?

$
0
0
Using module Semiconductors for computing Electric Potential for a certain geometry, sometimes during convergence process one can observe how few cells close to both exterior limits of a domain and contacts (ohmic or Schottky), assume values clearly far from the expected. These sets of cells constitute a sort of splodges which prevent convergence uselessly. How to eliminate them? Can I manage the individual cell values in those cases when I am completely sure of a more appropriate ones?

the free vibration of the cantilever

$
0
0
Hi everyone:
The initial displacement D=1.05mm, is placed on the end of the cantilever beam, to solve the velocity response of the beam end. First, I select space dimension 2D, and add physics Beam, but the simulation result is wrong, when I select space dimension 2D, and add physics Solid Mechanics (solid) , the simulation result also is wrong, the velocity curve is increasing against time. Now I select space dimension 3D, and add physics Beam, but the simulation result is not the same to the experimental result and analytical result. I am confused. Hope some one give me a advice, thank you! Details see attachment.

Convergence problems simulating a thin film - structural mechanics coupling

$
0
0
Hey everybody!

My goal is to simulate the stresses during a load cycle in a bearing shell.

The problem:
I assembled a journal bearing in con-rod from which i cut off the major part so that the calculation time can be reduced and defined this as union. So I don't have to bother about contact problems. On the surface of the journal bearing I solve the hydrodynamic with help of the thin film module and a PDE which handles the equilibrium of an external force. On the top of the con-rod I fix the structure what seems to me as the best way. Assuming a rigid structure in the first case the simulation works fine until it reaches the physical borders when the time dependent load becomes to high and the fluid film breaks. Therefore I coupled the hydrodynamic with elastomechanic so that the bearing shell can deform and a higher load carrying capacity is possible.
Due to initial value problems I set up an extra study where I calculate an inital status, which acts as an inital solution for instationary time dependent main study.


I made first tests with small loads, checked how to link both physics best,... and it worked fine. Now I have the problem, that I can't reach convergence with the load I'd like to simulate.
The calculation takes very long and it stops at 4 %. Until this point the convergence behaviour looks very good but then it changes, see attached screenshots, and I can't figure out why. In my opinion its not a physical problem because at the time, when the calculation stops the maximum fluid pressure is about 20 bar, the minimum filmthickness at about 20µm whereby the stresses are quite low. I think the problem must be the way the solvers work. I dont have enough experience to cope with this kind of problem. Thats why I would like to ask help from u.

I add some screenshots of the model so u can imagine how it looks like and a screenshots of one of the two solvers. (both have the same specification)


With kind regards,
Philipp

Heat Transfer

$
0
0
Hello

I'm a student on the project, and I'm messing in COMSOL 4.0 a short time and have no notion of "heat transfer", I need to get the temperature to stabilize so you can compare with the real model!
We apply a certain power and determine how long it will work, but as time increases the temperature to stabilize her investment keeps rising gradually not getting any results!
Material I am using is to dispel "Al" and the power part'm using "Tin"!
If someone can help me be grateful, and what commands I can use to stabilize the material the right power exercised

Thank you!

Cyclic Voltammetry (CV)

$
0
0
1) I am using the electroanalysis physics module with the built-in CV study only. How does one define the two electrodes which the CV is being applied, in the case of a two-electrode setup?

2) Is there a way to conduct the CV using a three electrode setup (2 electrodes + reference point)?

initial Values of Variables Solved For

$
0
0
Hi, I have someting to know.

As you know, the initial condition set is as follows:

initial Values of Variables Solved For
Method: initial expression, Solution
Solution: Zero, Solver1, Solver2

in magnetic problem solving,
I made study1 for low frequency, and study2 for high frequency for the better convergence.
So I'd like to solve at first study1 and after that study2 by using initial condition from study1.
How do I set the initial values?

in study2
initial Values of Variables Solved For
Method:Solution
Solution: Solver1

is it right?

What's the difference between initial expression and Solution in Method set?
and when can i use Zero for Solution?

In addition, I changed Solution to Solver2(from Solver1), and I got the very diffenrent calculation result.
Could you explain the reason?

Thank you in advance.


Convergence problem in a simple geometry with Joule heating module

$
0
0
Hi guys,
I am modeling a simple geometry (some stacked layers) with um and nm dimensions in comsol 4.3b.
The convergence plot goes down rapidly at first and after some iterations it seems that the software changes the value of parameters. Then the convergence plot goes up and the cycle continues. I have many upward and downward trends in my convergence plot and can not get fully converged solution. The plot is attached as well. Could anyone suggest me what should I do? Should we play with different solution methods, which I did and couldn't get satisfactory results? I appreciate for any help.
Thank you.

Mhd model of a crucible furnace

$
0
0
Dear collegues,
i´m setting the procedure to evaluate the movement of the molten metal in a simple crucible furnace: reading the papers from comsol ( www.comsol.asia/paper/download...3709/quintana_presentation.pdf , www.comsol.com/paper/download/83615/roy_paper.pdf )and from internet (www.comsol.eu/paper/download/14619/Ernst_pres.pdf) i understood that the ac/dc, the moving mesh ale and the turbulent flow k-epsi module should be used. in few hours i hope to start the project by myself, so in the way to let you have a look at it. In any case, if you have knowledge about how to set the solve processing , are well accepted

Thanks in advance
Gianmarco

Weak form for 2-D elasticity

$
0
0
In an effort to learn COMSOL, I decided to model the two-dimensional elasticity equations using the weak form. I was surprised looking through the documentation that there are very few example problems, so I didn't have much to go on. I think I'm deriving these correctly, but I could use another pair of eyes:

-(lamb+2*mu)*test(ux)*ux-mu*test(uy)*uy-(lamb+mu)*test(ux)*vy

-(lamb+2*mu)*test(vy)*vy-mu*test(vx)*vx-(lamb+mu)*test(vy)*ux

lamb and mu are the lame coefficients. u=u(x,y) and v=v(x,y) are the displacement fields.

One question that I can't figure out is how to apply stress boundary conditions. Can anyone help with that?

Successive simulation in COMSOL 4.3

$
0
0
Hello everybody!
I am curently doing simulations of fluid flows in microfluidic channels. The total simulation domain is quite big and a mesh that could provide a solution with resonable precision would take too much memory.

Can I somehow partition the simulation domain so that I run computation only on a part of it and when it is done I use the outlet data as input to the next channel section and run the simulation on this section only, etc.? If yes, how? Has anyone experience with such tricks?

Thanks in advance,
Daniel

Optimization problem

$
0
0
Hi everybody!
Please give a hand.
Knowing the experimental value of the function (for example, the maximum temperature of the body at a specific point), I need to find the optimal value of parameter (for example, coefficient of thermal conductivity) at which calculated body temperature should be as close to experimental. This parameter can be in a certain range of values. How can I solve this task in Comsol?
Thank you.
Viewing all 3988 articles
Browse latest View live


Latest Images

<script src="https://jsc.adskeeper.com/r/s/rssing.com.1596347.js" async> </script>