...
(my Mod hat is off)
This will be my last public posting on this issue, but you've made public comments regarding my Phase I observations and so will I:
I find your assessment a bit perplexing since you had been away for so long and had no commentary during my experiment as it was being conducted...only well after the fact and it was very subjective/declarative. I accept specific quantifiable critique, not "generalized". This, of course, is a big problem on other forums and not usually here.
"That's why you resorted to a statistical analysis."
This is another problem...I did not seek/resort to Glenn's assistance. He volunteered it when other's like yourself fell silent. In fact, several others helped when I was struggling with the Laser-spot analysis. Other people here stepped up THEN, not weeks later with hindsight commentary.
So, to summarize, I do not agree with your belated assessment, mainly because you have no specific critique, but generalized statements such as "there are times at which the magnetron is turned on and there is no such reaction force". What about the majority of the time when something was noted?
Glenn has the statistical analysis which I don't believe you can dismiss off hand. Rather than take up bandwidth here, please PM him to sort things out. I only ask one thing:
Whether accidently or by intent, you have given people the impression that my Phase I observational tests were null. I do not agree with that, both from what I personally witnessed with my own eyes (ala Shell) and more importantly, what Glenn relayed to me thru statistical analysis. Check out my FT2B video where the "aha" moment was and compare that to FT1, FT2 and FT2A where there was none.
Look, if we don't have the computing power to do your fluid dynamics analysis...so be it, don't even bring it up. I happen to believe that lift can be predicted and accounted for statistically and there is no phantom downwards "turbulence" that must be accounted for. Did I have an ideal experiment? No. Did I do things simply and within a small budget? Yes.
Got that off my chest, now back to mod-mode.
I accept specific quantifiable critique, not "generalized". This, of course, is a big problem on other forums and not usually here.
I happen to believe that lift can be predicted and accounted for statistically and there is no phantom downwards "turbulence" that must be accounted for.
@rfmguy:
It should always be the case that any forces delivered by thermally lifted air will have a net upward component over the time of heating, although at times it could push in any other direction given a specific measurement (due to the chaos inherent in vortices etc.).
So one way to think of it would be as thus: if the magnetron-on effect you observed consistently pushed the weight of your frustrum pointed down (or resting position) above the resting weight, for greater than half the runtime, there is almost no way for it to be attributable to thermal activity.
I haven't seen your test numbers specifically, but from the sound of it, I'd say that you either observed thrust or some magnetic artifact, especially given the fact that hot objects at rest tend to create a barrier of static air around their surface that deflect incoming thermal convection. (This is why a convection oven runs a fan, because moving air actually maximizes thermal change around a radiative body.)
In other words you could solve this by putting an emdrive in a very chaotic air environment and testing then, which is why I have advocated in the past for heating the entire test assembly to several hundred degrees passively before running experiments with the magnetron/antenna on.
<snip>
So, a computational fluid dynamics model of rfmwguy's experiment natural convection interaction with his turning the magnetron on and off would be fun to see.
What CFD computer program do you have experience with, and are you willing to take your time to model the natural convection in these experiments?
Since this is read by a wide audience, this question is really addressed at large: it would be useful for somebody to model the low Reynolds number, natural convection in these experiments with a CFD model<snip>
Also, just on a side note, wouldn't this be rather hellish to model? It seems like having a fine mesh like rfmwguy's EMDrive would require an ultra fine grid to get near accurate results of the air interacting with the heated mesh, not to mention the difficulty of creating the model. I am sure there are tools that I am unaware of for modeling this kind of surface though.
What kind of models do you need? I made a simple IGES model of rfmwguy's magnetron housing located here:
http://forum.nasaspaceflight.com/index.php?topic=39004.msg1459224#msg1459224
I may have some free time in the next day or two to make more elaborate models, if given enough information.Thanks so much for doing the model. I hope someone here provides you more info. Dimensional is all I need. Perhaps @Rodal can help since he is wanting this type of data. Here is the video that shows the magnetron ignited:
...
See this message regarding Computational Fluid Dynamics (CD+FD) vis-a-vis rfmwguy's experiment:
https://forum.nasaspaceflight.com/index.php?topic=39004.msg1465399#msg1465399
There may be some confusion in rfmwguy's response (unless I am missing something):
1) The message deals only with a simple IGES model of rfmwguy's magnetron housing, not with a CFD analysis
2) The dimensions needed for rfmwguy's experiment must be given by rfmwguy since rfmwguy is the one that run the experiment. I don't understand on what basis rfmwguy writes << I hope someone here provides you more info. Dimensional is all I need.>>
3) The issues involved with doing a CFD analysis, as previously discussed, have to do with computational resources. Since even Meep models are presently run to only 0.01 microseconds because of lack of computational speed, and since as previously discussed a CFD transient 3D analysis with a full Navier-Stokes solution would make children's play of the Meep analysis. Compound that by the fact that a CFD analysis with the Navier Stokes equation would be nonlinear and hence much more demanding to formulate properly and to interpret properly.
Also, to cut down on storage requirements, some discussion on what exactly needs to be saved for analysis would help. Along the lines of "20 time slices every 100 time slices" or some such, (which meep supports if I remember correctly), and a discussion of single precision vs double precision.
.../...
3) During the extremely short runs of Meep up to now, the Poynting vector and stresses were increasing at an exponential rate. What does the equilibrium balance between Poynting vector field rate and stress gradient look like vs time?
.../...
By "exponential rate" do you mean exp(-t/tau) as a first order charge at constant power against leaking power proportional to stored energy, asymptotically reaching a plateau (would make sense), or exp(t/tau) as ever increasing (at increasing rates) values, i.e. diverging (would not make sense, but your phrasing leaves ambiguity). Sorry for the nitpicking.
At hollidays with limited connectivity, happy celebrations everyone.I showed this over multiple posts, using Mathematica to post-process the Meep data output by Aero (*). See for example http://forum.nasaspaceflight.com/index.php?topic=37642.msg1406309#msg1406309
that shows an excellent fit of the net force (on the flat ends) data (from which the force has been computed and shown during the last two cycles ending at 0.013 microseconds from the time at which the Microwave feed was turned on) using the following fit:
A t +B (exp(t/tau) Sin [C t + D] + E
where the circles represent the data points output by Meep and the solid line represents the fitted model (with excellent R^2 = 0.999353)
I also wrote:Quotethe present Finite Difference model (from which the force has been computed at the last two cycles ending at 0.013 microseconds from the time at which the Microwave feed was turned on), would have to be marched forward for 1,000 times longer, to a total of 10 microseconds, for the force to be magnified by the calculated exponential growth to a value of 10 microNewtons (for an inputPower of 43 watts). Given the fact that the present Meep model takes an hour to run on a good PC modern computer, 1,000 hours of computer time represents over 41 days of computing time. Thus running the Meep model to steady state is impractical. Rather than using a supercomputer to perform such a computation, I suggest to use an implicit (unconditionally stable) Finite Difference model in time (rather than the explicit time difference model presently used that is subject to stability problems that limit the maximum finite difference time step). Such implicit finite difference models are well known (I developed a version of them in my PhD thesis 35 years ago) and can be run much faster than explicit FD models. There are also numerous alternative numerical schemes that are more accurate than Finite Differences.
The need to march forward the analysis to the order of tens of microseconds (instead of the extremely small time of 0.013 microseconds in the present Meep analysis) to ascertain what is going on agrees with a recent estimate by TheTraveller (where he proposed 30something microseconds was needed).
__________
(*) I have not seen anyone else similarly analyzing the Meep data in these threads, as they usually just output the electric or magnetic field values at one point in (extremely early) time (usually around 0.01 microseconds
-many times the Meep electric or magnetic field output are shown without indicating at what time, or without showing the audience that the output is increasing with time -), without analysis of the stresses or analysis of the time variation of the Meep model. Also, I have not seen the Meep model validated by checking it against the experimental results and COMSOL FEA analysis for TM 212 of NASA.
.../...
3) During the extremely short runs of Meep up to now, the Poynting vector and stresses were increasing at an exponential rate. What does the equilibrium balance between Poynting vector field rate and stress gradient look like vs time?
.../...
By "exponential rate" do you mean exp(-t/tau) as a first order charge at constant power against leaking power proportional to stored energy, asymptotically reaching a plateau (would make sense), or exp(t/tau) as ever increasing (at increasing rates) values, i.e. diverging (would not make sense, but your phrasing leaves ambiguity). Sorry for the nitpicking.
At hollidays with limited connectivity, happy celebrations everyone.I showed this over multiple posts, using Mathematica to post-process the Meep data output by Aero (*). See for example http://forum.nasaspaceflight.com/index.php?topic=37642.msg1406309#msg1406309
that shows an excellent fit of the net force (on the flat ends) data (from which the force has been computed and shown during the last two cycles ending at 0.013 microseconds from the time at which the Microwave feed was turned on) using the following fit:
A t +B (exp(t/tau) Sin [C t + D] + E
where the circles represent the data points output by Meep and the solid line represents the fitted model (with excellent R^2 = 0.999353)
I also wrote:Quotethe present Finite Difference model (from which the force has been computed at the last two cycles ending at 0.013 microseconds from the time at which the Microwave feed was turned on), would have to be marched forward for 1,000 times longer, to a total of 10 microseconds, for the force to be magnified by the calculated exponential growth to a value of 10 microNewtons (for an inputPower of 43 watts). Given the fact that the present Meep model takes an hour to run on a good PC modern computer, 1,000 hours of computer time represents over 41 days of computing time. Thus running the Meep model to steady state is impractical. Rather than using a supercomputer to perform such a computation, I suggest to use an implicit (unconditionally stable) Finite Difference model in time (rather than the explicit time difference model presently used that is subject to stability problems that limit the maximum finite difference time step). Such implicit finite difference models are well known (I developed a version of them in my PhD thesis 35 years ago) and can be run much faster than explicit FD models. There are also numerous alternative numerical schemes that are more accurate than Finite Differences.
The need to march forward the analysis to the order of tens of microseconds (instead of the extremely small time of 0.013 microseconds in the present Meep analysis) to ascertain what is going on agrees with a recent estimate by TheTraveller (where he proposed 30something microseconds was needed).
__________
(*) I have not seen anyone else similarly analyzing the Meep data in these threads, as they usually just output the electric or magnetic field values at one point in (extremely early) time (usually around 0.01 microseconds
-many times the Meep electric or magnetic field output are shown without indicating at what time, or without showing the audience that the output is increasing with time -), without analysis of the stresses or analysis of the time variation of the Meep model. Also, I have not seen the Meep model validated by checking it against the experimental results and COMSOL FEA analysis for TM 212 of NASA.
Oh yes, I see. Was less than attentive when this analysis work was carried on...
So, sorry, your phrase don't leave any ambiguity : the fitted curve does show amplitude "increasing at an exponential rate" (didn't thought it was possible) in the form e+t/τ1.
From your attached plot, in pN (10-12 Newton), t in time slice (meep units) :
Fxx=6.17×10-6 t + 1.53×10-4 e0.0222 t sin(1.23 t+0.938)+1.85×10-4
Which leaves me in a state of confusion as to why ? Doesn't make physical sense to me.
Surely, such an exponentially increasing value can't go on forever, this is diverging (and quick !). I understand we are looking at the first .01 µs or so of a phenomenon that reaches stability rather above 10µs, so we are at a thousandth of a complete view of the initial dynamic (assuming an instantaneous power on of the microwave generator). So is it an initial knee at the very first periods, that is to be followed by an inflection point and an overall slow down in the form max_amplitude×(1-e-t/τ2) where τ2>>τ1 ?
The motive of my questioning is that I would expect, overall (notwithstanding a specific shape of increasing rate around the few very first periods) something like a first order "charge" (starting at t=0 with E=0) :
E(t) stored electromagnetic energy in the cavity
Power_in(t)=P (constant)
Power_loss(t)=(ω/Q)E(t) (on average around a cycle)
dE/dt=Power_in-Power_loss
=> dE/dt=P-(ω/Q)E
=> E(t)=τ2×P×(1-e-t/τ2) with τ2=Q/ω
Here for Q using the "other common definition for Q is the ratio of the energy stored in the oscillating resonator to the energy dissipated per cycle by damping processes" (from wikipedia Q factor)
The amplitude you are fitting is that of a (not averaged) radiation_pressure(t), isn't this amplitude proportional to stored energy E(t) ?
Am I at lost ? Are we too early in the time marching simulation to hope see (fit) such a τ2 ?
And what about the linear "offset" 6.17×10-6t, it is also diverging (albeit more gently) so it has also to reach a plateau at some point ? Wasn't the simulations on which you based those analysis done at the time with wrong parameters (now corrected) that gave unnaturally high Q ?


I am having trouble reconciling all the dimensions in SeeShell's blueprint. The waveguide is shown as being 9.97cm high and protruding from the frustrum wall by 7.6 cm at the top. The end of the waveguide is given as 17.79cm from the central axis. This means that the point where the top of the waveguide meets the frustrum is 17.79 - 7.6 = 10.19 cm from the axis. The large end diameter is 29.5cm, so bigR = 14.75cm. Thus the slope of the wall isatan((14.75 - 10.19) / 9.97) = atan( 4.56 / 9.97) = atan( .457 ) = 24.56 degrees.
But the frustrum itself is described as Height 24.79cm, smallD = 17cm, smallR = 8.5cm, for a slope ofatan( (14.75-8.5) / 24.79 ) = atan( 6.25 / 24.79 ) = atan( .252 ) = 14.15 degrees
Something is inconsistent. I notice the blueprint does not show the dimensions for height or small end diameter, so perhaps the numbers I am using for that, from an earlier post, are obsolete.
W.r.t. building the cavity... Cutting the template out of sheet copper is easy, but then making the actual frustum out of it - not so much. In fact I am still trying to find out a way... The next attempt is going to be with a "jig" made from 2 wood boards with pre-routed paths for both the big and the small sides...
Or is it all because my copper is C110 and not the O2 free? Would the latter be much more pliable and hence easier to make into a frustum?
...
(my Mod hat is off)
This will be my last public posting on this issue, but you've made public comments regarding my Phase I observations and so will I:
I find your assessment a bit perplexing since you had been away for so long and had no commentary during my experiment as it was being conducted...only well after the fact and it was very subjective/declarative. I accept specific quantifiable critique, not "generalized". This, of course, is a big problem on other forums and not usually here.
"That's why you resorted to a statistical analysis."
This is another problem...I did not seek/resort to Glenn's assistance. He volunteered it when other's like yourself fell silent. In fact, several others helped when I was struggling with the Laser-spot analysis. Other people here stepped up THEN, not weeks later with hindsight commentary.
So, to summarize, I do not agree with your belated assessment, mainly because you have no specific critique, but generalized statements such as "there are times at which the magnetron is turned on and there is no such reaction force". What about the majority of the time when something was noted?
Glenn has the statistical analysis which I don't believe you can dismiss off hand. Rather than take up bandwidth here, please PM him to sort things out. I only ask one thing:
Whether accidently or by intent, you have given people the impression that my Phase I observational tests were null. I do not agree with that, both from what I personally witnessed with my own eyes (ala Shell) and more importantly, what Glenn relayed to me thru statistical analysis. Check out my FT2B video where the "aha" moment was and compare that to FT1, FT2 and FT2A where there was none.
Look, if we don't have the computing power to do your fluid dynamics analysis...so be it, don't even bring it up. I happen to believe that lift can be predicted and accounted for statistically and there is no phantom downwards "turbulence" that must be accounted for. Did I have an ideal experiment? No. Did I do things simply and within a small budget? Yes.
Got that off my chest, now back to mod-mode.
W.r.t. building the cavity... Cutting the template out of sheet copper is easy, but then making the actual frustum out of it - not so much. In fact I am still trying to find out a way... The next attempt is going to be with a "jig" made from 2 wood boards with pre-routed paths for both the big and the small sides...
Or is it all because my copper is C110 and not the O2 free? Would the latter be much more pliable and hence easier to make into a frustum?
...
(my Mod hat is off)
This will be my last public posting on this issue, but you've made public comments regarding my Phase I observations and so will I:
I find your assessment a bit perplexing since you had been away for so long and had no commentary during my experiment as it was being conducted...only well after the fact and it was very subjective/declarative. I accept specific quantifiable critique, not "generalized". This, of course, is a big problem on other forums and not usually here.
"That's why you resorted to a statistical analysis."
This is another problem...I did not seek/resort to Glenn's assistance. He volunteered it when other's like yourself fell silent. In fact, several others helped when I was struggling with the Laser-spot analysis. Other people here stepped up THEN, not weeks later with hindsight commentary.
So, to summarize, I do not agree with your belated assessment, mainly because you have no specific critique, but generalized statements such as "there are times at which the magnetron is turned on and there is no such reaction force". What about the majority of the time when something was noted?
Glenn has the statistical analysis which I don't believe you can dismiss off hand. Rather than take up bandwidth here, please PM him to sort things out. I only ask one thing:
Whether accidently or by intent, you have given people the impression that my Phase I observational tests were null. I do not agree with that, both from what I personally witnessed with my own eyes (ala Shell) and more importantly, what Glenn relayed to me thru statistical analysis. Check out my FT2B video where the "aha" moment was and compare that to FT1, FT2 and FT2A where there was none.
Look, if we don't have the computing power to do your fluid dynamics analysis...so be it, don't even bring it up. I happen to believe that lift can be predicted and accounted for statistically and there is no phantom downwards "turbulence" that must be accounted for. Did I have an ideal experiment? No. Did I do things simply and within a small budget? Yes.
Got that off my chest, now back to mod-mode.
I have some knowledge of the difficulty of characterizing complex fluid mechanics behavior so can't offer any more guidance than Dr. Rodal has already offered. What I would suggest is to perform a series of experiments that will reveal how much of what you have observed is due to thermal effects. For example what happens when the fustrum is removed and a box made of Copper mesh is put in its place? The microwave energy will still be safely contained but there should be no em-drive thrust effect. The box can be made so it has the same weight as the fustrum. Another experiment would be to replace the magnetron with a weighted box (to equal the weight of the magnetron) covered with Aluminum heat sink fins and with a 200 Watt halogen bulb inside for heating. If we recall the experiments done to prove the existence of the ether; they instead showed the wave nature of light.
My overly simple analysis of rfmwguy's em-drive experiment is that the intense heating of the magnetron results in a lowering of the air pressure above it. This causes lift, just laike a plane's wing generates lift when air flows over it. However that can change because of air currents so it's possible to see a downward thrust as well. I have done my share of odd experiments. I have never been convinced the effect I am hoping to see is there if it is not always observable each and every time I throw the switch. If there is a natural pendulum motion can I stop that movement by throwing the switch at the right time? Can I increase the oscillation by throwing the switch at the opposite side of the cycle? And most of the time even after this has been observed there are other reasons for doubting a positive conclusion. Quoting Richard Feynman: "Science is a way of trying not to fool yourself. The first principle is that you must not fool yourself, and you are the easiest person to fool."
W.r.t. building the cavity... Cutting the template out of sheet copper is easy, but then making the actual frustum out of it - not so much. In fact I am still trying to find out a way... The next attempt is going to be with a "jig" made from 2 wood boards with pre-routed paths for both the big and the small sides...
Or is it all because my copper is C110 and not the O2 free? Would the latter be much more pliable and hence easier to make into a frustum?
You'll find there are several different dimensions due to the confusion of the tuning chamber. But it's as simple as my small top plate is 165mm my large bottom plate is 295mm, I'm seeking TE modes and that determines the plate spacing along with the angle. That said my frustum isn't the exact solution from the calculations because I have a tuning chamber on top of the small section of the frustum.
So are you just trying to build what I've done, or model it, knowing will help me help you.
...Take your time and make sure your surfaces are kept clean, I wore gloves to keep the oils and acids from my hands off the copper.
...Take your time and make sure your surfaces are kept clean, I wore gloves to keep the oils and acids from my hands off the copper.
Shell,
Any chance to see your plans and photos of your magnetron to coax adapter build?
Phil
You'll find there are several different dimensions due to the confusion of the tuning chamber. But it's as simple as my small top plate is 165mm my large bottom plate is 295mm, I'm seeking TE modes and that determines the plate spacing along with the angle. That said my frustum isn't the exact solution from the calculations because I have a tuning chamber on top of the small section of the frustum.
So are you just trying to build what I've done, or model it, knowing will help me help you.
Now I see what you are getting at. It is not a pure frustrum, but a frustrum with a cylinder on top. So the slope of the frustrum walls matches the ratio given by the waveguide dimensions in your blueprint, but the actual overall height is adjustable by moving the circular plate inside the cylinder. (See picture below, with exagerated cylinder height)
That probably does interesting things to the impendance. Watch out for heating where the two shapes meet.
I am trying to replicate your design in MEEP.