如何人为设定变量的最小值

xiaoyaoyou

比如PDE方程中的一个变量是rho，计算过程中可能会出现rho小于零的情况，那么此时我想人为将此时的rho强制赋值为e-5,请问如何操作。多谢。

Vivi_Song

可以用判断性的语句来实现：rho*(rho>=0)+1e-5*(rho<0).它表示的意思是当满足括号里的不等式时，括号项为1，否则为0.这样就可以实现密度值的切换。

xiaoyaoyou

首先谢谢Vivi_Song朋友的回答。
我的目的是限制rho本身的值（使其不能为负数），比如我的计算周期为100s，计算周期为1s，当我计算到50s时，空间中有些单元的rho由正变负，那么我想要对这些为负数的rho重新赋值为e-5，以便进行第51秒的计算。
我现在的初步想法是给PDE方程组另外联立一个global ODEs and DAEs, 写入rho-if(rho-10^-5>=0,rho,10^-5),但是计算下来并无效果。
请问有何建议？

Vivi_Song

xiaoyaoyou 发表于 2012-6-5 16:58
首先谢谢Vivi_Song朋友的回答。
我的目的是限制rho本身的值（使其不能为负数），比如我的计算周期为100s， ...

看您是用rho来表示的，所以下意识的把它想象成密度了。好吧，您是在用自定义PDE来求解方程组，当rho由正变负以后，计算要继续，是不是因为rho还被代入到其他方程中计算其他未知量，而rho这个方程就是为了求其他未知量？如果是这样的话，那您为什么不把用到rho的方程中直接把rho用rho*(rho>=0)+1e-5*(rho<0)替换掉？

xiaoyaoyou

Vivi_Song 发表于 2012-6-5 17:26
看您是用rho来表示的，所以下意识的把它想象成密度了。好吧，您是在用自定义PDE来求解方程组，当rho由正 ...

rho确实如您所说是质量密度，只是我的物理模型中质量密度的变化范围会很大，会有rho=0(即真空)的情况，于是计算PDE方程组的过程中难免有rho<=0的情况出现（这会直接引起计算结果不收敛），因此此时我想对这些单元赋予e-5的值来代表物理上的真空。
如您所说rho是为了求其他未知量，但是下一步的计算中仍要求解质量守恒方程，因此替换掉rho为负值的那一部分也是必须的。