 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1739 | 回復(fù): 7 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
markwu4銅蟲 (著名寫手)
|
[求助]
請有comsol license的人提供我Knowledgebase 103/952 已有2人參與
|
||
|
有l(wèi)icense的人請?zhí)峁┪襅nowledgebase 103, 952內(nèi)容 http://www.comsol.com/support/knowledgebase/103/ http://www.comsol.com/support/knowledgebase/952/ 我算出負(fù)的濃度 聽說里面有一些技巧可以解決這類問題 |
回復(fù)此樓» 本帖已獲得的紅花（最新10朵）
tang6920» 猜你喜歡
一志愿西電085401數(shù)一英一299求調(diào)劑 六級521
已經(jīng)有3人回復(fù)
-
生物考研337分求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑：085600材料與化工，考材科基，總分319
已經(jīng)有14人回復(fù)
-
285求調(diào)劑
已經(jīng)有10人回復(fù)
-
化學(xué)工程085602 305分求調(diào)劑
已經(jīng)有29人回復(fù)
-
理學(xué)07化學(xué) 303求調(diào)劑
已經(jīng)有5人回復(fù)
-
081200-11408-276學(xué)碩求調(diào)劑
已經(jīng)有3人回復(fù)
-
土木304求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿中國科學(xué)院大學(xué)265求調(diào)劑
已經(jīng)有7人回復(fù)
-
286求調(diào)劑
已經(jīng)有8人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- comsol里面AC/DC模塊-----Magnetic Fields(mf)仿真
已經(jīng)有3人回復(fù)
- 地質(zhì)地球所管理部門與儀器技術(shù)研發(fā)人員招聘啟事
已經(jīng)有0人回復(fù)
- 中國科學(xué)院地質(zhì)與地球物理研究所 招聘
已經(jīng)有0人回復(fù)
|
Solution Number: 952 Title: Avoiding negative concentrations [Requires COMSOL Access] Platform: All Platforms Applies to: All Products Versions: All versions Categories: Solver, Physics Keywords: convergence negative concentrations Problem Description Why do I get negative values of the concentration in my diffusion, convection and reaction model? Why do I see values that go below, or above, their lowest, or highest, possible values based upon the boundary conditions applied for problems involving diffusion and convection? It is clearly unphysical. Solution There are several possible reasons for, and ways to avoid the issue of over- and under-shoot of the solution: Numerical Errors A common reason is numerical noise: For example, when the variable for concentration approaches zero, the numerical noise may become significant in comparison to the value (which is close to zero). If you notice negative concentrations of very small magnitude, numerical noise is probably the cause and does not affect the problem much in diffusion/convection problems without reactions. Reacting Domains: Use Danckwerts inflow type If you have an Inflow condition at a domain containing a chemical reaction, or a domain with reactions on the walls, use Boundary Condition Type: Flux (Danckwerts), available in the Inflow settings. This inflow condition type is designed to avoid negative oscillations and also to speed up the solution of reacting flow problems. Danckwerts. For an example how to use this, see the Separation Through Dialysis example, and also the file provided for download below. Discontinuous Boundary and Initial Conditions: Smoothen the Settings Another cause for locally slightly negative concentrations is a discontinuity in space or time, for instance in the initial condition. As an example, consider the one dimensional time dependent Convection-Diffusion-Reaction equation The initial condition is zero concentration within the domain and the boundary conditions are set to 1 and 0 concentration at the left and right boundary, respectively. The physical interpretation of this is an initially sharp, gradually diffusing front moving in the positive x-direction. However, for finite element shape functions (second order Lagrange), only continuous functions are admissible as solutions, for which reason the discontinuous initial value is modified before the time-iterations can begin. This often results in a small dip in the solution for t=0. The dip traverses to the first time steps too. The concentration will locally be slightly negative at t=0, as shown in the figure. Negative concentrations Figure 1: Solutions to the time-dependent diffusion problem with a boundary condition at the left (c=1) that creates a discontinuity with the initial condition (c=0) at t=0. The finite element method gives rise to a swing below zero in the leftmost element. The dip close to t=0 sometimes results in oscillations and convergence problems. This problem can be avoided by using one of COMSOL Multiphysics's built in smoothed step functions smooths out the initial discontinuity in a controlled way. For instance, in the example above, you could, instead of the uniform zero initial condition, use a smoothed step transition as the initial condition. For more information please refer to Solution 905: Smoothed step transitions. Condition the Reaction Term to be Larger than Zero for Concentrations Approaching Zero A significant negative concentration often indicates that the underlying mathematical model does not correctly describe the physics. In this case the approach mentioned above does not take care of the problem. One potential cause is that you have a constant sink in your reaction term, which is an approximation that only works for large concentrations. When the concentration reaches zero, the reaction term continues to consume the species, finally resulting in negative concentration. This can be avoided by making sure that your reaction rate is such that when the concentration of the species approaches zero, then so does the species sink. This can be achieved for instance by writing max(eps^2,Q). eps is an internal COMSOL constant that is a very small number in the order or 10-15. This type of expression is also helpful if you want to avoid Q being 0, for example if you apply the logarithm to it somewhere within the geometry, or if the symbolic derivative runs a risk of containing division by zero. Make Sure the Mesh Resolves the Problem Insufficient mesh resolution can result in dips below zero. Convergence problems are often the underlying issue when negative concentrations are observed in high convection regimes (high Peclet number) and with large reaction terms or fast kinetics (high Damköhler number).It can also be useful to investigate whether the negative concentration gets better or worse with mesh refinement. If, better the problem is most likely mesh related. If worse, it is probably the mathematical model that causes the problem. In Knowledge Base solution 103 a few tricks are presented which may alleviate these types of problems without extensively refining the mesh. Note: The Danckwerts inflow condition requires any one of the add-on modules Chemical Reaction Engineering Module, Batteries & Fuel Cells Module, Corrosion MOdule, Electrodeposition Module, Electrochemistry Module, or the Subsurface Flow Module Module. |
8樓2015-11-10 21:41:56
2樓2014-10-14 09:48:38
markwu4
銅蟲 (著名寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 14530.5
- 紅花: 2
- 帖子: 1311
- 在線: 254.2小時(shí)
- 蟲號: 2301321
- 注冊: 2013-02-26
- 專業(yè): 電化學(xué)
3樓2014-10-14 16:05:33
albertwowwow
木蟲 (小有名氣)
- 應(yīng)助: 80 (初中生)
- 金幣: 3804.6
- 散金: 20
- 紅花: 12
- 帖子: 211
- 在線: 83.4小時(shí)
- 蟲號: 3292789
- 注冊: 2014-06-26
- 性別: GG
- 專業(yè): 煤炭地下開采
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
markwu4: 金幣+10, ★★★★★最佳答案, 謝謝 因?yàn)橐话胍呀?jīng)分給第一個(gè)應(yīng)助者 2014-10-16 06:48:48
markwu4: 金幣+10, ★★★★★最佳答案, 謝謝 因?yàn)橐话胍呀?jīng)分給第一個(gè)應(yīng)助者 2014-10-16 06:48:48
|
Description Stationary models with high reaction rates, Ri, or more generally, large source terms, may result in convergence problems. You can for instance get error messages such as "No convergence, even when using the minimum damping factor", "NaN repeatedly found in solution", or "Error: Failed to find a solution". Now what could you do to get around this problem? Try to solve the problem using reduced source terms, and then gradually increase them, using the previous solution when solving, until you reach the solution to your original problem. The following approach uses the parametric solver to gradually increase the source term. COMSOL Multiphysics 4 and later Multiply the reaction (source) term (Ri) by a multiplication factor variable, for example "k". Click the Stationary node and select the Auxiliary sweep (in version 4.4) or Continuation (in earlier versions) check box. Enter k as the Sweep parameter In Parameter values, enter a range of increasing values, starting at a low value and increasing to unity, for example "1e-3 1e-2 1e-1 1". Click Compute. COMSOL Multiphysics 3.5a and earlier Multiply the reaction (source) term (Ri) by a multiplication factor variable, for example "k". In the *'Solver Parameters** dialog box, select the Parametric Solver and specify Name of Parameter to that of the multiplication factor, k. Enter the List of Parameter Values as a range of increasing values, starting at a low value and increasing to unity, for example "1e-3 1e-2 1e-1 1". Click OK. Click Solve. COMSOL Multiphysics will now start to solve the problem using the first multiplication factor. Then, if that works, the solution will be used as the initial guess for the next parameter value. This process will continue until the last parameter value is reached. Because you specified the last parameter value to be unity, the last solution will correspond to your original problem formulation. Use the time dependent solver. By making the problem time dependent, you usually get a smoother convergence. COMSOL Multiphysics 4 and later Right-click the Model node and select Add Study. Select Time Dependent and click Finish. In the Times edit field, enter a range of times to solve for, for example "0 1e3". Click Compute. COMSOL Multiphysics 3.5a and earlier Open Solver Parameters. Select the Time dependent solver and specify a range of Times to solve for, for example 0 1e3". Click OK, then Solve. Make sure to solve the problem for large enough times, such that the solution does not vary any longer. This implies that the solution has reached steady-state. Thus, COMSOL Multiphysics has calculated the stationary solution! If the solution process never gets past t = 0, see Solution 964 for more information. If none of the previous methods would help, try using them in combination with refining the mesh in regions where the source term is large. 這個(gè)是103的 |
4樓2014-10-16 06:14:44
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +5 | 丟掉懶惰 2026-03-27 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +5 | hanamiko 2026-03-30 | 5/250 |
|
|
[考研] 本科211生物醫(yī)學(xué)工程085409求調(diào)劑339分 +5 | 里子木yy 2026-03-29 | 5/250 |
|
|
|
[考研] 環(huán)境工程 085701，267求調(diào)劑 +11 | minht 2026-03-29 | 11/550 |
|
|
|
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑，總分281 +4 | 橙子（勝意） 2026-03-30 | 4/200 |
|
|
|
[考研] 本科211總分289，08工學(xué)真心求調(diào)劑 +3 | utopiaE 2026-03-30 | 3/150 |
|
|
|
[有機(jī)交流]
 甲基亞磺磺酸鈉和甲基磺酸酯反應(yīng)機(jī)理
 10+3
|
kaobao456 2026-03-29 | 4/200 |
|
|
|
[考研] 085600 286分 材料求調(diào)劑 +11 | 麻辣魷魚 2026-03-27 | 12/600 |
|
|
|
[考研] 一志愿鄭州大學(xué)，080500學(xué)碩，總分317分求調(diào)劑 +10 | 舉個(gè)栗子oi 2026-03-24 | 11/550 |
|
|
|
[考研] 085601材料工程找調(diào)劑 +17 | oatmealR 2026-03-29 | 18/900 |
|
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-26 | 7/350 |
|
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-24 | 6/300 |
|
|
|
[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
|
[考研] 本科雙非材料，跨考一志愿華電085801電氣，283求調(diào)劑，任何專業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
|
[考研] 本科新能源科學(xué)與工程，一志愿華理能動(dòng)285求調(diào)劑 +7 | AZMK 2026-03-28 | 11/550 |
|
|
|
[考研] 339求調(diào)劑，想調(diào)回江蘇 +6 | 烤麥芽 2026-03-27 | 8/400 |
|
|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +6 | 阿姨復(fù)古不過 2026-03-27 | 6/300 |
|
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個(gè)好覺_16 2026-03-24 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|