 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 8 | 1/1 | 返回列表 |
| 查看: 1716 | 回復(fù): 7 | ||
markwu4銅蟲 (著名寫手)
|
[求助]
請(qǐng)有comsol license的人提供我Knowledgebase 103/952 已有2人參與
|
|
有l(wèi)icense的人請(qǐng)?zhí)峁┪襅nowledgebase 103, 952內(nèi)容 http://www.comsol.com/support/knowledgebase/103/ http://www.comsol.com/support/knowledgebase/952/ 我算出負(fù)的濃度 聽(tīng)說(shuō)里面有一些技巧可以解決這類問(wèn)題 |
回復(fù)此樓» 本帖已獲得的紅花（最新10朵）
tang6920» 猜你喜歡
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有4人回復(fù)
-
生物學(xué)調(diào)劑招人�。�！
已經(jīng)有4人回復(fù)
-
0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有17人回復(fù)
-
一志愿武漢理工材料工程專碩調(diào)劑
已經(jīng)有5人回復(fù)
-
0856調(diào)劑，是學(xué)校就去
已經(jīng)有7人回復(fù)
-
申博26年
已經(jīng)有3人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有10人回復(fù)
-
0703化學(xué)調(diào)劑 ，六級(jí)已過(guò)，有科研經(jīng)歷
已經(jīng)有12人回復(fù)
-
294求調(diào)劑材料與化工專碩
已經(jīng)有13人回復(fù)
-
288求調(diào)劑
已經(jīng)有14人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- comsol里面AC/DC模塊-----Magnetic Fields(mf)仿真
已經(jīng)有3人回復(fù)
- 地質(zhì)地球所管理部門與儀器技術(shù)研發(fā)人員招聘啟事
已經(jīng)有0人回復(fù)
- 中國(guó)科學(xué)院地質(zhì)與地球物理研究所 招聘
已經(jīng)有0人回復(fù)
2樓2014-10-14 09:48:38
markwu4
銅蟲 (著名寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 14530.5
- 紅花: 2
- 帖子: 1311
- 在線: 254.2小時(shí)
- 蟲號(hào): 2301321
- 注冊(cè): 2013-02-26
- 專業(yè): 電化學(xué)
3樓2014-10-14 16:05:33
albertwowwow
木蟲 (小有名氣)
- 應(yīng)助: 80 (初中生)
- 金幣: 3804.6
- 散金: 20
- 紅花: 12
- 帖子: 211
- 在線: 83.4小時(shí)
- 蟲號(hào): 3292789
- 注冊(cè): 2014-06-26
- 性別: GG
- 專業(yè): 煤炭地下開(kāi)采
【答案】應(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
albertwowwow
木蟲 (小有名氣)
- 應(yīng)助: 80 (初中生)
- 金幣: 3804.6
- 散金: 20
- 紅花: 12
- 帖子: 211
- 在線: 83.4小時(shí)
- 蟲號(hào): 3292789
- 注冊(cè): 2014-06-26
- 性別: GG
- 專業(yè): 煤炭地下開(kāi)采
【答案】應(yīng)助回帖
|
Numerical errors The most common reason for negative concentrations is numerical noise: when the species concentration approaches zero, the numerical noise becomes significant in comparison to the concentration. If you are seeing negative concentrations of very small magnitude, numerical noise is probably the cause. In a pure convection-diffusion problem, the scale of the concentration does not matter, so in order to avoid the problem, you can add an arbitrary baseline concentration to keep the resulting overall concentration above zero. However, with a reaction term which depends on the concentration, the scale and origin do matter, which means that you need to think of other ways to keep the concentration strictly positive. Discontinuous concentrations Another common cause for 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 and Diffusion equation where the convection is in the positive x-direction (i.e. the direction vector u=1), with a uniformly zero initial condition, and boundary conditions which set the concentration at the end nodes to one and zero respectively. The physical interpretation of this PDE is an initially sharp, gradually diffusing front moving in the positive x-direction. However, for the default shape function (second order Lagrange), only continuous functions are admissible as FEM 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, and in the above example, the concentration will locally be slightly negative at t=0, as shown in the figure below. Solutions to the time-dependent Convection and diffusion equation at times t=0, 0.01, 0.1, 0.2, 0.3. This type of behavior can also result in wildly oscillating solutions and convergence problems. This problem can be avoided by using one of COMSOL Multiphysics's built in smoothed step functions in order to smooth out the initial discontinuity in a controlled manner. For example, in the Convection and Diffusion problem described above, you could, instead of the uniformly zero initial condition, use a smoothed step transition as the initial condition in order to avoid negative values for concentration at t=0. For more information please refer to Knowledge Base solution 905. Incorrect reaction term Usually a significantly negative concentration (i.e. not noise around zero) indicates that the underlying mathematical model does not correctly describe the physics. In this case "fixing" the numerics 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. In order to avoid this problem you need to make 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 good if you want to avoid Q being 0, for example if you apply the logarithm to it somewhere, or if the symbolic derivative runs a risk of containing division by zero. Remember that COMSOL symbolically differentiates all expressions that contribute to the Jacobian. Mesh resolution One other thing that can be indicated by significantly negative concentrations is the lack of mesh resolution. The resulting convergence problems are often the underlying issue when negative concentrations are observed in high convection regimes (high Peclet number) and in those with large reaction terms or fast kinetics (high Damkohler number). In Knowledge Base entry 103 a few tricks are presented which may alleviate these types of problems without extensively refining the mesh. It can also be useful to investigate whether the negative concentration problem gets better or worse with mesh refinement. If better, then you know in which direction to go. If worse, then the physics of the model probably need to be checked. Formulate logarithmic concentrations A nice way of eliminating mesh resolution problems and negative dips is to use the logarithm of the concentration and not the concentration itself as the dependent variable. The reason for this is that a linearly varying mesh can sometimes not capture the exponential behavior of the concentration changes. In addition, modeling the logarithm of the concentration ensures that the real concentration can never become negative during the solution process. 這個(gè)是192的 |
5樓2014-10-16 06:15:31
6樓2015-03-20 22:06:42
小帆遠(yuǎn)航
木蟲 (正式寫手)
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 3953.7
- 散金: 20
- 紅花: 4
- 帖子: 420
- 在線: 124.5小時(shí)
- 蟲號(hào): 3303719
- 注冊(cè): 2014-07-02
- 性別: GG
- 專業(yè): 機(jī)械結(jié)構(gòu)強(qiáng)度學(xué)
7樓2015-07-08 20:10:53
|
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
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 markwu4 的主題更新
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
