Skip to main content
33 events
when toggle format what by license comment
Jun 16, 2022 at 4:20 history edited user21
edited tags
Jun 16, 2020 at 9:23 history edited CommunityBot
Commonmark migration
Dec 13, 2019 at 23:25 review Close votes
Dec 15, 2019 at 14:56
S Dec 13, 2019 at 23:21 history bounty ended DMH16
S Dec 13, 2019 at 23:21 history notice removed DMH16
Dec 13, 2019 at 23:21 vote accept DMH16
Dec 13, 2019 at 9:06 answer added xzczd timeline score: 9
Dec 12, 2019 at 7:14 comment added DMH16 @xzczd I thought I implemented it correctly. Then, I am looking for a solution of the problem I presented above
Dec 12, 2019 at 6:52 comment added xzczd Well, what answer are you expecting? I've also checked your MATLAB code, the moving b.c. isn't considered therein, if I've understand it correctly. Translating the MATLAB code into Mathematica code won't help.
Dec 12, 2019 at 6:00 history tweeted twitter.com/StackMma/status/1205004528816406529
S Dec 12, 2019 at 5:43 history bounty started DMH16
S Dec 12, 2019 at 5:43 history notice added DMH16 Authoritative reference needed
Dec 11, 2019 at 11:28 comment added xzczd So, there's no moving boundary $h(t)$ in your MATLAB code? "I also looked at the links you provided but none handle the case with $x<0$" The $x=0$ boundary isn't moving, it's much easier to handle compared to the moving boundary, we also have several posts about the topic, for example: mathematica.stackexchange.com/a/121739/1871
Dec 11, 2019 at 10:39 comment added DMH16 @xzczd I simply created two functions for which when $x<s$ it returned a $\rho_1c_1$, if $x>s$ it returned $\rho_2c_2$ other wise it returned the average. I then did something similar to handle the other constants. Then I called these two functions in my iterative scheme that I presented. Although I am not sure how to do this in mathematica as I do no really know how to work with the implementations of matrices (for example a dynamic allocation of data on it), which is needed when implementing the iterative scheme. I also looked at the links you provided but none handle the case with $x<0$
Dec 11, 2019 at 9:58 comment added xzczd ……Then how did you handle it in MATLAB?
Dec 11, 2019 at 6:34 comment added DMH16 that's what I am having trouble with dong in mathematica
Dec 11, 2019 at 5:26 comment added xzczd Well, how is the moving boundary handled by this scheme?
Dec 10, 2019 at 18:23 comment added DMH16 @xzczd I included the scheme I followed
Dec 10, 2019 at 18:23 history edited DMH16 CC BY-SA 4.0
added 473 characters in body
Dec 10, 2019 at 15:55 review Close votes
Dec 12, 2019 at 5:45
Dec 10, 2019 at 7:05 comment added DMH16 @xzczd My goal here was to convert my matlab code (which works) to the mathematica one, mostly for comparison, so that is why I insisted on the difference scheme. The parameters, as I mentioned in the post, are picked arbitrarily. I will add the scheme as soon as I get back from work. Thanks anyways.
Dec 10, 2019 at 6:58 comment added xzczd The question is much better after the edits. (Still, you need to show us the parameters and i.c.s. ) I'd say the MATLAB code isn't quite related to this question in my view. Is the finite difference scheme used by the paper necessary? If so, it'll be better to add the difference scheme to the question directly, because it's easier for us to understand the difference formula in traditional math notation. If not, then I believe the combination of DChange, pdetoode and NDSolve is a better choice for solving the problem, and again, we don't need to know how the MATLAB code looks like.
Dec 10, 2019 at 6:43 comment added DMH16 @xzczd Unfortuntely, I had no other way of providing the code. Although I understand that a broken link may be an issue since such code would not be retrievable anymore. Do you think I should post it? I began translating the code myself, however being that I am somewhat inexperienced in Mathematica, I preferred posing the question to other people here, also because I did not get very far (and used many while and for loops which apparently are bad practice in Mathematica). If the question does not conform to the standards of the site, feel free to close it. I will look at the links u suggested
Dec 10, 2019 at 6:37 history edited b3m2a1 CC BY-SA 4.0
edited title
Dec 10, 2019 at 6:35 comment added xzczd Finally, we already have a few posts about moving b.c., have you read them?: 1. mathematica.stackexchange.com/q/58596/1871 2. mathematica.stackexchange.com/q/130962/1871 3. mathematica.stackexchange.com/q/95358/1871 4. mathematica.stackexchange.com/q/114967/1871 5. mathematica.stackexchange.com/q/184920/1871
Dec 10, 2019 at 6:35 comment added xzczd @b3m2a1 I've seen this post the moment it's posted, and have been fighting against the urge to vote to close since then. As mentioned by b3m2a1, this is a big ask, which is out of scope of this site in many cases. (There exist subtle cases of course. ) Also, it's a bad idea to provide code using external download link, because the link may be easily dead. (This aleady happened for several times in this site. )
Dec 10, 2019 at 6:35 history edited DMH16 CC BY-SA 4.0
Stated the model on the question
Dec 10, 2019 at 6:29 history edited DMH16 CC BY-SA 4.0
Stated the model on the question
Dec 10, 2019 at 6:13 comment added DMH16 @b3m2a1 will do. Thank you
Dec 10, 2019 at 6:09 comment added b3m2a1 I'd say most element-wise tasks can be done better as vectorized operations and once you do that there's really no difference between Mathematica and MATLAB. That's really neither here nor there, though. Looks like you just need to use the FDM, from what I glanced at in that image. Lots of questions here pertain to that but maybe xczd will come by and help you. Also how about as a compromise for not having some preliminary code, you post the relevant equations here in LaTeX form?
Dec 10, 2019 at 6:06 comment added DMH16 @b3m2a1 Yes what I meant was element-wise manipulation. And it appears to me it feels more natural in matlab when doing so, or maybe I’m not experienced enough in mathematica. Regarding the method, as I mentioned above, I used the finite difference method which there are a lot of snippets around the web. I will try to post some code, maybe the general idea. Anyways, any help is appreciated. It’s just a curiosity
Dec 10, 2019 at 6:03 comment added b3m2a1 I'd say Mathematica and MATLAB are about the same when it comes to 95% of the major matrix ops, but maybe you mean element-wise manipulation? Also I think it'd be nice to at least post some code detailing the method used. It's a big ask for us to read, translate, and do all the mental processing without any sort of starting point.
Dec 10, 2019 at 5:42 history asked DMH16 CC BY-SA 4.0