Reputation
1,599
Top tag
Next privilege 2,000 Rep.
Edit questions and answers
Badges
12 31
Impact
~46k people reached

Jul
29
comment Potential flow over cylinder with FEM: Pressure field is flipped
Is there some way I could choose the range of data represented by the PlotLegend? Thank you for your answer; I'll try this very soon.
Jul
29
comment Specifying NeumannValue for Axial load instead of transverse load (Plane Stress situation)
@user21 You make a good point. I'll try this with MeshOrder
Jul
28
comment Potential flow over cylinder with FEM: Pressure field is flipped
@marcob yikes my mistake. The code I commented out may be treated as "old version". Sorry about that.
Jul
10
comment Specifying NeumannValue for Axial load instead of transverse load (Plane Stress situation)
Yes, I just realized that. Now I am on the hunt for a generalized Hookes law that describes a cantilever with transverse load (original mma example). Would you mind putting the final statement in your answer "The results are unchanged from those produced by the code in the question but differ from that of the test case, because the differential equations are different." in bold? This is vital. I was under the impression that the plane stress operator in full would still solve the linear deflection situation while it solves only the transverse load condition.
Jul
10
comment Specifying NeumannValue for Axial load instead of transverse load (Plane Stress situation)
I suppose I am unable to understand what the x and y directions are. Is x along the thickness of this bar element and y along its span? If so, should the load be an +xload? (of course the difference between - and + loads would be the sign of the deflection alone)
Jul
9
comment How to solve heat equation with Robin type conditions with NDSolve?
Are you sure of your boundary conditions? This example has Robin boundary conditions. (it is a link to a pdf file)
Jul
9
comment Displaying a matrix with integral entries
What I have noticed is m.n would "work" to produce a meaningful result for matrices m and n but m//MatrixForm . n//MatrixForm does not compute. This would seem to be not explicit in the documentation?
Jul
9
comment Specifying NeumannValue for Axial load instead of transverse load (Plane Stress situation)
@bbgodfrey yes. But in a purely axial load setting, there will be deflection only in axial x direction (u here)
Jul
9
comment Specifying NeumannValue for Axial load instead of transverse load (Plane Stress situation)
@bbgodfrey But with inactive, the operators do execute. I still do get an interpolating function result. The real question for me is to understand why the example problem in the mma documentation uses inactive if it is of no use for this situation. I'll clean up the code like you say.
Jul
9
comment Specifying NeumannValue for Axial load instead of transverse load (Plane Stress situation)
@bbgodfrey I have not used the Dirichlet condition gamma although I could have. I'm not sure about the "inactive", I forget what it does. Yes a rectangular mesh works better but I want to compare rect to tri :)
Jul
9
comment Specifying NeumannValue for Axial load instead of transverse load (Plane Stress situation)
@bbgodfrey The transverse direction, although included should not affect the results. The object here is that the plane stress operator could be used for this example. If you look at modern FEM software, this axial problem includes transverse direction but only an axial load. The question here is how to model the axial load in terms of Neumann or other values mma understands.
May
18
comment Choosing FEM element type and mesh refinement
@MichaelE2 Yes, I just found the QuadElement. However, am I to understand that I cannot force Mathematica to use a QuadElement on an irregular shape? It always defauls to a TriangleElement.
May
18
comment Choosing FEM element type and mesh refinement
@MichaelE2 Yes. But how can I ensure that "MeshRefinementFunction" allows me to specify regions of refinement? Also, the "MeshElementType" -> "TriangleElement": Are there other "MeshElementType" that aren't mentioned in the Help menu?
May
18
comment Choosing FEM element type and mesh refinement
@AlexeiBoulbitch I meant tri along CURVED boundaries and Quad everywhere else.
May
18
comment Choosing FEM element type and mesh refinement
@AlexeiBoulbitch yes, I forgot to mention that. Tri around the boundary and quads everywhere else. This refinement might not have physical significance but will serve only as an example.
May
18
comment Choosing FEM element type and mesh refinement
@Nasser I did that. However it pastes weird symbols. I'll try that again just to make sure. Thank you!
Feb
26
comment Changing FrameTicks of ListPlot
Thanks for your answer though! :) I still keep track of questions I started a few years ago ;)
Aug
25
comment Proving (self) similarity with Mathematica - Reccurrence Plots, Similarity Plots etc
@s.s.o That is an interesting point you make. I considered this later in 2012 but dropped the idea entirely as I had to defend and finish my PhD! Perhaps I will do this on an other day to include in my lectures! :) Thank you for your answer/comment!
Apr
28
comment NDSolve - sampling for result during the computation
@xzczd Didn't know one could do that. However, how would one control the time step size? I imagine it would be thro' MaxStepFraction or MaxStepSize?
Apr
25
comment NDSolve - sampling for result during the computation
I performed simulations on MD for my Master's degree and I am curious, if you don't mind. Are you simulating Argon? With periodic BCs? Is the eqn you are solving the convection-diffusion equation? In 2012 I had a similar question on SE and I am wondering if it would help you. Good luck!