How to calculate an atomic bomb's critical mass

Поделиться
HTML-код
  • Опубликовано: 27 июн 2024
  • Support the channel: ko-fi.com/jkzero
    Step-by-step guide on how to derive and solve the neutron diffusion equation to determine the critical mass of an atomic bomb. It involves solving a partial differential equation but every step is presented. In the end a general formula for a nuclear critical mass is found.
    Relevant links:
    • Critical Mass: when the atomic bomb got real • Critical Mass: when th...
    • Physics of a Nuclear Explosion • Physics of a nuclear e...
    • Nuclear Weapons Q&A #1 • Nuclear Weapons Q&A #1
    • Series on Nuclear Weapons Physics • Physics of Nuclear Wea...
    • "Critical Assembly Simulator" by Alex Wellerstein blog.nuclearsecrecy.com/misc/...
    • The Feynman Lectures on Physics www.feynmanlectures.caltech.edu/
    Errata:
    • at 10:10 the Laplacian of f should say: ∇²f= (1/r²)*d/dr(r² df/dr) (not 1/r² in the derivative bracket), this is a typo and not an error in the calculation; thanks to @fatassmonkey37 for pointing this out
    • at 15:58 the subscripts are reversed, it should say: λ_f = 16.89 cm, λ_t = 3.596 cm; thanks to @lucacastronuovo9789 for pointing this out
    Credits:
    "Critical Assembly Simulator" courtesy of Alex Wellerstein blog.nuclearsecrecy.com/misc/...
    Diffusion by JrPol under CC BY 3.0 Deed creativecommons.org/licenses/...
    Metal cube by JOGOS Public Assets under CC BY-SA 4.0 Deed creativecommons.org/licenses/...
    Chapters:
    00:00 Neutron diffusion equation
    05:41 Solution
    15:24 Critical mass
    20:13 Advanced solution
  • НаукаНаука

Комментарии • 368

  • @jakubstawarczyk
    @jakubstawarczyk 5 месяцев назад +61

    Amazing, as a theoretical physicist myself, I'm looking forward to see more mathematical videos about this topic 😁
    Edit: makes me wanna sit down with numerical methods and analyze more reallistic example of shrinking and expanding core (true R(t))

    • @pressaltf4forfreevbucks179
      @pressaltf4forfreevbucks179 5 месяцев назад +16

      Is that actually 2500 dollars?!?!? I should become a theoritical physicist too💀💀💀

    • @jakubstawarczyk
      @jakubstawarczyk 5 месяцев назад +1

      @@pressaltf4forfreevbucks179 Hong Kong $ ... so only 64 US$ I'm poor too 😭😉 (and given I'm now in Philippines it shows 2500PHP... not sure why as my visa was charged 500HKD that I've selected)

    • @jkzero
      @jkzero  5 месяцев назад +14

      WOW, thank you so much for your kind support, I am grateful for the appreciation. This drives me to keep making these videos. I enjoy making them (despite the consumption of my free time) so it is very gratifying to know that the content is appreciated. Thanks again!

    • @Keshav-il2hp
      @Keshav-il2hp 5 месяцев назад +1

      its pesos lol 2500 is $44
      still great@@pressaltf4forfreevbucks179

    • @multiarray2320
      @multiarray2320 5 месяцев назад +1

      ​@@pressaltf4forfreevbucks179 no, its about 40 bucks.

  • @hypercomms2001
    @hypercomms2001 5 месяцев назад +91

    Absolutely fascinating. I have always been interested how this was calculated ever since I did a school assignment on the Manhattan project when I was 14 in 1974.

    • @jkzero
      @jkzero  5 месяцев назад +17

      glad you liked it; similar story here: I read about the Manhattan Project when I was a kid and got obsessed and decided to be a physicist. That was a few decades back and here I am.

    • @-VinhKhang_yearsago
      @-VinhKhang_yearsago 5 месяцев назад

      @@jkzero congratulations

  • @kirkhamandy
    @kirkhamandy 5 месяцев назад +29

    17:15 _This will very quickly start blowing up_ LOL, I've often heard of exponentials "blowing up" but here, never a more true definition 😂

    • @jkzero
      @jkzero  5 месяцев назад +6

      glad to see that people appreciate the puns... there is a Taylor Swift joke hidden in a previous video than maybe was too subtle because only one person noticed it

  • @verdienthusiast3868
    @verdienthusiast3868 5 месяцев назад +11

    Man, I am a high schooler and I am trying my best to learn Physics, a subject that I love with all of my heart, I am self studing calculus to, one day, contribute a little to this wonderful field, the things that I find more interesting are the possibilities that the Island of stability can give us, and I want to dedicate my life to that.
    Take my money (earned by surveys so no problem for me) and buy a cup of something

    • @jkzero
      @jkzero  5 месяцев назад +1

      thanks so much for your support, the appreciation from strangers watching my videos really drives me to keep making them. Thank you!

    • @jacobharris5894
      @jacobharris5894 5 месяцев назад +1

      Sounds like your more proactive with your education than I was at your age. Good for you. I did take calculus in high school but not until the second half of senior year after taking pre calculus over the summer. I wish I did it sooner or tried testing into some more advanced classes because I could of gotten the entire calc chain including differential equations out of the way before college if I wasn’t afraid to take AP courses. I could of also gotten physics 1 and 2 as transfer credits I think if I took the calc based versions. It would have been intense but easier than doing it in college and put me 2 years ahead for my Bachelor’s. It may have lowered my GPA but I didn’t go to an elite University so it wouldn’t have mattered anyway. And if your aiming for an elite school, AP courses are a good way to stand out if you do well in your classes. In hindsight, I think it’s worthwhile to take some risks and not play it too safe in high school because if you have big goals your going to have to do that at some point.

    • @verdienthusiast3868
      @verdienthusiast3868 3 месяца назад

      ​@@jacobharris5894 hello, thank you for your words and your experiences.
      I'll surely try all I can, here we have many good universities but to attend the best one you need to pass a very difficult exam with at least 65-70% (only 32 people) else you are out.
      It is needed to start studing at least 1-2 years before.
      I still have more than 2 years, if my will is enough I’ll make it *fingers crossed*, if no there are many others very good.

  • @MrFib112358
    @MrFib112358 5 месяцев назад +48

    This is the first time that I am excited for a homework assignment and I can't believe it came from RUclips, I really want to try to do this! Brilliant video and thanks for educating us on this fascinating topic!

    • @jkzero
      @jkzero  5 месяцев назад +4

      happy to see that some viewers have accepted the challenge, let me know what you get

  • @Hexnano
    @Hexnano 5 месяцев назад +6

    Thank you for this video, I got stuck while working on building my atomic bomb and this was very helpful.

    • @jkzero
      @jkzero  5 месяцев назад +5

      I am glad you liked the video; however, I cannot wish you well in your endeavors :)

  • @KiwiExpressCream
    @KiwiExpressCream 5 месяцев назад +10

    Not going to lie, the maths is a bit beyond me, but I love how you worked out the critical masses. All of your videos have been amazing and I always look forward to the next one.

    • @jkzero
      @jkzero  5 месяцев назад +3

      thanks for the feedback, this video was math heavy; people seem to like these so I will make more but I will also continue with the standard video of some history sprinkled with a less-intense calculation on interesting topics, even beyond nuclear weapons; I don't want the FBI knocking on my door :D

  • @ion8264
    @ion8264 5 месяцев назад +8

    Shame how underrated this channel is😢

    • @jkzero
      @jkzero  5 месяцев назад +1

      thanks for that, you can support the channel by liking, commenting, and sharing the videos, that engagement drives the algorithm that gives the channel exposure. Thanks in advance.

  • @footbballa1784
    @footbballa1784 5 месяцев назад +11

    This channel deserves more views and subs, one of the best Physics YT channels out there keep it up!

    • @jkzero
      @jkzero  5 месяцев назад +3

      Thanks for the endorsement; I wish I could get to wider audiences but the RUclips algorithm is driven by engagement so you can actively help the channel by liking, subscribing, and sharing. This type of support is highly appreciated so I can continue making videos.

  • @fablearchitect7645
    @fablearchitect7645 5 месяцев назад +14

    Can you please do a video on deriving the transfer function for a nuclear reactor in terms of control rods vs output thermal power?

    • @terrycole472
      @terrycole472 Месяц назад

      You could always look at RUclips channel "Nuclear Engineering Lectures" NE560, here:
      ruclips.net/video/iSbLUBuaT78/видео.html
      - but remember transfer function analyses are perturbative, and the full range of control rod operation is not, so you have to hope things stay linear. Also remember transfer function theory often looks at the Zero-Power state (to avoid pesky thermal feedback considerations) of a reactor at criticality. It's a useful tool but extending it to more general cases is not trivial.

  • @Asterism_Desmos
    @Asterism_Desmos 5 месяцев назад +3

    You know it’s a good day when you get a notification for Dr. Jorge S. Diaz 🔥

  • @VaradMahashabde
    @VaradMahashabde 5 месяцев назад +5

    I believe that the elementary solution does not correspond to neutrons not escaping the core, but rather that neutrons that leave the core are permanently lost, so the function is 0 outside the sphere.
    I think an actual condition for neutrons not leaving the center would be the spatial derivative wrt to r being 0 at the center, as expected from Fick's law, instead of the actual value. Otherwise sub-criticality is not achievable because the number of neutrons would not only be conserved but also continuously multiplying.
    Other than that, loved the video, think your calm and collected teaching style is great.

    • @jkzero
      @jkzero  5 месяцев назад

      I am glad you liked the video, thanks for the feedback. Honestly, "elementary" boundary condition is quite artificial, in fact, it is only used because it allows an analytical solution for R_c in terms of the relevant parameters (
      u, cross-sections, density, etc.), it gives a glimpse on what is driving R_c and the critical mass. However, this is the only reason this condition is introduced, mostly for academic purposes. The "advanced" boundary condition is what has to be used in a practical calculation.

  • @omargaber3122
    @omargaber3122 5 месяцев назад +5

    Actually, I am amazed by your work. I don’t know what to say. I thought a lot about the effort you made to make this video. It must have been hard. I don’t think anyone can thank you as you deserve. In fact, I am traveling to Russia, where there are no means of support or Communication with the outside world, even RUclips ads are prohibited, so I turned on a VPN so that ads would appear on your videos and I could watch them. Indeed, I watched 14 ads and visited 11 advertiser’s websites. This took some time, but it is nothing compared to the time that you spent making this video, and perhaps this is worth a few cents to the channel, but I feel that it is more moral support than material support. Thank you from the bottom of my heart, you brilliant mind.❤

    • @jkzero
      @jkzero  5 месяцев назад +2

      Thanks for the endorsement; I wish I could get to wider audiences but the RUclips algorithm is driven by engagement so you can actively help the channel by liking, subscribing, and sharing. This type of support is highly appreciated so I can continue making videos. Thanks also for watching the ads, I do the same with small channels that I watch, if there is an ad I let it go until the end because I know I am contributing with a few cents to the creator.

    • @omargaber3122
      @omargaber3122 5 месяцев назад

      ​@@jkzeroDone❤

  • @Isaac-mt9hx
    @Isaac-mt9hx 5 месяцев назад +3

    This is an amazing video! Made me even more excited to learn DPEs next semester!

    • @jkzero
      @jkzero  5 месяцев назад +1

      If anything in my videos can make people excited about physics and math, I can say: mission accomplished! Thanks for sharing

  • @dylangabriel2703
    @dylangabriel2703 5 месяцев назад +2

    Thank you I was trying to find this for a long long time

  • @fatassmonkey37
    @fatassmonkey37 5 месяцев назад +1

    there's a mistake at 10:15, the part of the laplacian in spherical is equal to (1/r²)*d/dr(r² df/dr) (not 1/r² in the brackets)

    • @jkzero
      @jkzero  5 месяцев назад +2

      oh no... you are totally right, I messed up the r² term already at 10:10. Thanks so much for pointing this out, it is a good catch. Fortunately, it doesn't affect the calculation because it is correctly applied, this is just a typo on the video but not an error in the calculation; however, I will immediately add an erratum in the description. I do my best to avoid these typos but after watching everything many times some minor details slip through. Thanks again.

  • @davidrandell2224
    @davidrandell2224 5 месяцев назад +1

    A proton is a collection of 1836 expanding electrons and add a bouncing expanding electron makes a hydrogen atom. “G” calculated from first principles- the hydrogen atom- in 2002. No energy, charge, photons, waves, spin, fields, potential, quantum,space, time, space-time etc. All Standard Theory/Model was replaced by Expansion Theory in 2002. “The Final Theory: Rethinking Our Scientific Legacy “, Mark McCutcheon for proper physics.

  • @CarlosMartinez-ep5mi
    @CarlosMartinez-ep5mi 5 месяцев назад +1

    A little mistake, min 10:17. On radial term you have inside the brackets (1/r^2 df/dt) but it is ( r^2 df/dt ). After you use the function f=u/r all come back to be good.

    • @jkzero
      @jkzero  5 месяцев назад +1

      Thank for the notification, you are right and another viewer also let me know about this, I messed up the r² term already at 10:10 so I have included an erratum in the video description. Thanks so much for pointing this out and in such a courteous way (not very common on the internet), it is a good catch. Fortunately, as you proved, it doesn't affect the calculation because it is correctly applied, this is just a typo on the video but not an error in the calculation. I do my best to avoid these typos but after watching everything many times some minor details slip through. Thanks again.

  • @The_Fancy_Duck
    @The_Fancy_Duck 5 месяцев назад +3

    Very underrated RUclips channel.

    • @jkzero
      @jkzero  5 месяцев назад

      Thanks for the endorsement; I wish I could get to wider audiences but the RUclips algorithm is driven by engagement so you can actively help the channel by liking, subscribing, and sharing. This type of support is highly appreciated so I can continue making videos.

  • @paulm1241
    @paulm1241 5 месяцев назад +15

    Excellent! Could you make a video on the implosion physics ? For example, it is often said that although it was known that hollow cores would have a much better yield, the original bomb used a "conservative" design based on a solid core ("Christy core") because the maths of hollow core implosion was too difficult to solve at the time (but it was finally implemented in all designs a few years later). I could not find any more detailed info on that on the web.

    • @jkzero
      @jkzero  5 месяцев назад +8

      did you check my "Nuclear Weapons Q&A" video? I don't go into the details of the solid vs. hollow core but I do calculate how much the core is compressed because somebody asked about it.

    • @gaHuJIa_Macmep
      @gaHuJIa_Macmep 5 месяцев назад

      Read the NWFAQ by Carey Sublette.

  • @tomashernandez8711
    @tomashernandez8711 5 месяцев назад +1

    Man, this theory was taught to me by a professor the last semester, amazing the application and the smooth way this works to explain such a beauty phenomenon

    • @jkzero
      @jkzero  5 месяцев назад

      Awesome! I am curious to know what brings viewers to the channel, were you searching for something in particular or did the 'mighty algorithm' find you?

  • @ytashu33
    @ytashu33 5 месяцев назад +4

    My kind of YT video, just at the right intersection of math and physics, delivered at the right cadence, not too fast, not too slow, loved it. One question, if sin(kR) = 0, wouldn't that give us multiple solutions? ie, kR = n * pi? Any particular physical constraint that we must consider only the lowest numbered "harmonic", n = 1?

    • @1995error
      @1995error 5 месяцев назад +2

      If you used solutions with multiples of pi you would have areas with negative neutron density which I guess would not make much physical sense, but I am kinda guessing. The reason why you can use multiple solutions in QM is that you multiply the function with it’s complex conjugate before you get the probability, which means that it will always be positive:)

    • @jkzero
      @jkzero  5 месяцев назад +2

      nice catch, that is a very nice observation; you are right, the general solution is kappa*R=n*pi; however, the higher "harmonics" will only make the bomb core bigger, we are searching for the minimum radius. Also, a solution for n>1 can lead to a negative density, which is also a non-physical solution. Mathematically, your point is 100% correct; however, here is when physics guides the final solution. Just like setting B=0 to remove the non-physical cos(kappa*r) term.

    • @ytashu33
      @ytashu33 5 месяцев назад

      Thank you Prof, that makes sense. You did proactively point out all other, non-valid, solutions and explained why they are not applicable. So this stuck with me. This was the only point, everything else i was able to follow-up rather clearly. These are not easy thing to communicate, you struck absolutely the right balance between explaining the points in sufficient detail, without getting bogged down. That balance is hard to get right, you hit the ball out of the park. Thank you!

    • @jkzero
      @jkzero  5 месяцев назад +1

      @@ytashu33 in the video at 14:20 I refer to the solution used as "the simplest solution" this is precisely to address your question. I could have explained more details but the video was already 25 min long.
      Thanks for your appreciation, I love teaching and created this channel mostly because I wanted to share fun and education stuff, I just never thought that a video only about solving a hard equation would produce this level of interest, which I find great. Welcome to the channel and make sure to check the other videos.

  • @latefoolstalk676
    @latefoolstalk676 3 месяца назад

    Wow, so well explained. Such advanced physics explained step by step. I hope more people checks your wonderful channel

    • @jkzero
      @jkzero  2 месяца назад

      Thanks; I hope to get to wider audiences but the RUclips algorithm is driven by engagement so you can actively help the channel by liking, subscribing, and sharing. This type of support is highly appreciated so I can continue making videos.

  • @christopherwithers1131
    @christopherwithers1131 5 месяцев назад +2

    Very nice! I will use some of this as an exercise for my PDE course.!

    • @jkzero
      @jkzero  5 месяцев назад +2

      this is great, I have been told of others using my videos in their classes, this is such an honor; curious to know how it is received, please report back on how this interesting applied problem is received by your students

  • @stoneaa756
    @stoneaa756 Месяц назад

    very very clear explanation! thank you

    • @jkzero
      @jkzero  Месяц назад

      Thanks, I hope you also check the other videos.

  • @pawelbjanowski
    @pawelbjanowski 5 месяцев назад +5

    A coffee/beer from a nuclear colleague from Poland, with kind regards!⚛

    • @jkzero
      @jkzero  5 месяцев назад +1

      Thank you so much for your kind support, I am grateful for the appreciation, which drives me to keep making these videos. Thanks again!

  • @cosmicnomad8575
    @cosmicnomad8575 3 месяца назад

    Amazing video and amazing channel!

    • @jkzero
      @jkzero  3 месяца назад

      Thanks for the visit and the positive feedback

  • @FPSIreland2
    @FPSIreland2 5 месяцев назад

    Thanks. I'll make sure to employ this new information... wisely...

  • @cewkins721
    @cewkins721 5 месяцев назад +5

    What a great video! i truly want to thank you for going through the effort of laying out the math in the latex format you were talking about, it looks clean, the math is a bit higher than my level but your explanation makes total sense, i will see if i can try the questions you put in the end, your content is really enjoyable!

    • @jkzero
      @jkzero  5 месяцев назад +2

      I am glad you liked the video and thanks for appreciating the effort, writing that code was a lot of work. I hope you give the "assignments" a try and let me know if any help is needed; I am happy that others have also showed interest in attempting it

  • @strangercydonian
    @strangercydonian 5 месяцев назад +1

    i would love to see a complete playlist about the calculus for nuclear reactions

    • @jkzero
      @jkzero  5 месяцев назад

      thanks for the request; I cannot guarantee to be able to fulfill all the requests but I am collecting suggestions, thanks

  • @brahimhihi56
    @brahimhihi56 5 месяцев назад +1

    Thanks a lot!!

  • @billynomates920
    @billynomates920 5 месяцев назад +2

    clicked on your video just as the cider i'd had earlier informed me that actually, it has not yet fully diffused through my system and that perhaps i'd get more from your video if i delayed watching it until tomorrow. til then! thank you.

    • @jkzero
      @jkzero  5 месяцев назад +2

      stay safe and enjoy the video later

  • @pacoramirez2714
    @pacoramirez2714 5 месяцев назад +1

    grande maestro! buen video

  • @danieleambrosini1681
    @danieleambrosini1681 5 месяцев назад +12

    Awesome! Make a lot of videos on physics because you are really talented! :) Love from Italy

    • @jkzero
      @jkzero  5 месяцев назад +6

      thanks, I appreciate the feedback. I hope to keep making videos, I really enjoy this and the interaction with the community. I am curious to know what brings viewers to the channel, were you searching for something in particular or did the 'mighty algorithm' find you?

    • @ion8264
      @ion8264 5 месяцев назад

      ​@@jkzeroI searched relaxing physics video and found you😂

    • @jkzero
      @jkzero  5 месяцев назад +5

      @@ion8264 nothing more relaxing than taking time off and spending free time solving partial differential equations 😃

    • @danieleambrosini1681
      @danieleambrosini1681 5 месяцев назад +2

      ​@@jkzero i really love physic and math, i study engineering so i really like to Watch videos about what i like, i found you about 5 months ago when you did the video about the Enola Gay, today i opened RUclips and in the home i found your video, i really enjoyed It!😊

    • @jkzero
      @jkzero  5 месяцев назад +3

      @@danieleambrosini1681 the Enola Gay video was my official debut with the channel and having a community of nerds interested in physics stories and that don't shy away from some math is just great. Thanks for returning for more!

  • @spencerwenzel7381
    @spencerwenzel7381 4 месяца назад

    I really appreciate that in every formula you describe what each variable represents and how it relates to the physical interpretation. Even if I don't completely understand the math, the implications are not completely lost and I am still learning. Thanks for making these.

    • @jkzero
      @jkzero  4 месяца назад

      Thanks, I am glad you liked the video. I am curious to know what brings viewers to the channel, were you searching for something in particular or did the 'mighty algorithm' find you?

    • @spencerwenzel7381
      @spencerwenzel7381 4 месяца назад

      @@jkzero Your video about Planck caught my interest initially and then I started looking through your other videos. However I initially saw it on my feed.

    • @jkzero
      @jkzero  4 месяца назад

      @@spencerwenzel7381 Thanks for sharing and I am glad the algorithm has been showing my video to so many people, I hope you find the other videos of interest too and welcome to the channel.

  • @Sokol_
    @Sokol_ 5 месяцев назад

    Thank you, i need it.

  • @chuckgrigsby9664
    @chuckgrigsby9664 2 месяца назад

    Another interesting condition involves uranium dissolved in aqueous solution. In solution, the uranium is not nearly so dense as is the case with uranium metal, but the water slows the neutrons down. It turns out that the fission cross-section for U-235 is higher for slower neutrons. Now, we have to take neutron scattering off of water molecules into account, and how that scattering slows down the neutrons. There have been a number of serious nuclear criticality accidents involving inadvertently exceeding the maximum uranium concentration limits for a given solution geometry. That has led to the imposition of critically safe geometries for places where fissile materials are dissolved in aqueous solutions.

  • @jeremyvirin6532
    @jeremyvirin6532 5 месяцев назад +1

    Incredible this video😍

  • @giovannimariatanda9251
    @giovannimariatanda9251 5 месяцев назад +1

    Molto interessante e ben fatto. Complimenti vivissimi dottor Jorge.

    • @jkzero
      @jkzero  5 месяцев назад

      Grazie mille! Thanks for watching and make sure to check the other videos

    • @giovannimariatanda9251
      @giovannimariatanda9251 5 месяцев назад

      Farollo sicuramente.
      Grazie.

  • @jperez7893
    @jperez7893 5 месяцев назад

    this is a fantastic video. the math is beautiful. no doubt that this is also the same equations used for the generation of nuclear power for thorium reactors. so for nuclear reactors, nu prime is not allowed to reach zero but always kept negative.
    amazing work

    • @jkzero
      @jkzero  5 месяцев назад

      thanks, I am glad you liked the video. Honestly, I am not an expert on reactors, I just made one video showing that they cannot explode like a bomb but I have a good friend who is an expert on thorium reactors, if you have questions on this topic make sure to write them in this post, I will try to have him in the future ruclips.net/user/postUgkxkfK_xvCYFftK_Wsw6wh0L6ECtL87-Dcv

  • @Unique-Concepts
    @Unique-Concepts 5 месяцев назад +2

    Subscribed!Thank you for this video.Please do more videos on advanced mathematical modelling.

    • @jkzero
      @jkzero  5 месяцев назад +1

      thanks fro the sub, I am glad you found the content of interest. I am already working on the next video and more juicy calculations are coming. I want to balance the content with some historical context but always including some high-level calculation to complement the story, and more importantly, for viewers to follow along and reproduce themselves if desired. I am curious to know what brings viewers to the channel, were you searching for something in particular or did the 'mighty algorithm' find you?

    • @Unique-Concepts
      @Unique-Concepts 5 месяцев назад +2

      @@jkzero Thanks for the reply, recently I was searching on the topic detonation wave hydrodynamics. I found that topic fascinating because of welch lab channel. Now RUclips recommended this video and absolutely love it...thanks again 👍🙏🙏👌👌👏👏😍

    • @jkzero
      @jkzero  5 месяцев назад

      @@Unique-Concepts thanks for sharing and great that the searching algorithm brought you here and thanks again for your support, welcome to the channel!
      I might not know much about the topic that you were searching for but you might find my seminar on blast waves of interest (link below). It is about a paper that I published in 2022 extending the result from Sedov-Taylor-von Neumann but written for undergraduate-level students of physics and engineering. It required solving hydrodynamic equations for properly describing an expanding blast wave and I include all the steps in the seminar, check it out and let me know if this helps ruclips.net/video/JySY4bkW5wY/видео.html

  • @7177YT
    @7177YT 5 месяцев назад

    Very cool! Thank you!

    • @jkzero
      @jkzero  5 месяцев назад

      Thanks for the endorsement; I wish I could get to wider audiences but the RUclips algorithm is driven by engagement so you can actively help the channel by liking, subscribing, and sharing. This type of support is highly appreciated so I can continue making videos.

  • @pappaflammyboi5799
    @pappaflammyboi5799 5 месяцев назад +2

    What would be the equations for critical mass if you had a perfectly spherical neutron reflector and a neutron half-life of 11 minutes? Thanks again for the video. I loved the math derivation.

  • @krwada
    @krwada 4 месяца назад

    A very nice video.
    You explain very clearly the basic physics concepts of a chain reaction.
    Laplacian, Ficks diffusion law, Continuity, conservation of Mass. All these concepts are within reach of any undergraduate engineering student.
    Many of these concepts were taught to me while as an undergrad in Chemical Engineering.
    Very well done!

    • @jkzero
      @jkzero  4 месяца назад

      Thanks you, I am glad you liked the video. You have a good point, there are many concepts packed but I did my best make the calculation flow. I warned the general viewer about the math because this video is in fact a response to the request from viewers after I shared the historical aspects of the calculation of the critical mass in a video some months back (ruclips.net/video/LduH7613QXw/видео.html). Please check the other videos, in my latest I began a series on quantum mechanics that with your background you might enjoy and also be surprised (many of the physicists viewers were surprised), check it out here ruclips.net/video/gXeAp_lyj9s/видео.html

    • @krwada
      @krwada 4 месяца назад

      @@jkzero I will definitely check out the video on quantum mechanics.
      There are also other very good RUclips channels which very nicely explain complex topics such as what you do.
      For example, a couple of days ago, the Veritasium channel showed a very nice explanation of various concepts in solid state physics, including:
      1. band gap
      2. the various bands, (conduction and valence)
      3. A very nice animation showing the Poisson distribution and how it creates the PN junction discontinuity.
      Lots of good stuff.
      ruclips.net/video/AF8d72mA41M/видео.html
      I am now a subscriber to your very fine channel!

    • @jkzero
      @jkzero  4 месяца назад

      Yes, Veritasium offers excellent content. Thanks for watching and the positive feedback. Also thanks for supporting the channel with your subscription and welcome to the channel!

  • @natanfreire8692
    @natanfreire8692 5 месяцев назад +1

    I love this channel

    • @jkzero
      @jkzero  5 месяцев назад

      Thanks for watching; you can actively help the channel by liking, subscribing, and sharing. This type of support is highly appreciated so I can continue making videos.

  • @giovannimariatanda9251
    @giovannimariatanda9251 5 месяцев назад

    Chiarissimo dottor Jorge, ho notato, con grande sorpresa e ammirazione, che ha vergato i suoi appunti in data 25 dicembre. La stimo ancor più per questo: chi può lavorare il giorno di Natale? Solo un appassionato, un entusiasta (etimo: en + theos = un dio che ti parla dentro) può fare una cosa del genere, perciò, c'è ne fosse bisogno, rinnovo e attesto la mia perspicua ammirazione per la sua persona e il suo lavoro. Se lei permette, appena avrò un poco di tempo, riporterò su un testo la sua disanima corredandola di ulteriori semplificazioni e ampliamenti per renderla ancor più comprensibile per i meno competenti. E naturalmente citerò la fonte, cioè: lei.
    Grazie.
    Forte nello spirito, tenero nel cuore.
    G.

  • @andresfeliperamirezgaviria8976
    @andresfeliperamirezgaviria8976 5 месяцев назад

    Excellent.

  • @enzolamande5963
    @enzolamande5963 5 месяцев назад +4

    Great Vidéo! In cylindrical coordinates we are using the Bessel Function right?

    • @jkzero
      @jkzero  5 месяцев назад +4

      Yes, the cylindrical case is a bit more complicated than the spherical and cubic core, and yeah, if you encounter a Bessel equation then you are on the right track

  • @faisalsheikh7846
    @faisalsheikh7846 5 месяцев назад +1

    Incredible❤

  • @sturrum5250
    @sturrum5250 5 месяцев назад

    First of all, excellent video. It's refreshing to see a subject that is often *talked* about get the mathematical treatment. I was wondering if there's somewhere I could find a derivation of your expression for the diffusion constant D = 1/3*lambda_f*. Most explanations I've seen are very handwave-y and not particularly rigorous.

    • @jkzero
      @jkzero  5 месяцев назад +1

      Thanks, I am glad you found the content of interest. I agree with you, most derivation of the diffusion equation are quite fishy. Regarding the explicit form of the diffusion constant the only legit derivation that I have found is in the appendix of the book The Physics of the Manhattan Project by B. Cameron Reed. It is a long a quite tedious calculation involving several angular integrations.

    • @sturrum5250
      @sturrum5250 5 месяцев назад

      @@jkzero Alright, thanks a whole bunch for the reference.

  • @MisterTutor2010
    @MisterTutor2010 5 месяцев назад +3

    The Differential Equation of DOOM :)

    • @jkzero
      @jkzero  5 месяцев назад +1

      I read this as "The Differential Equation of BOOM" but yeah, both work, I guess

  • @thespiciestmeatball
    @thespiciestmeatball 5 месяцев назад

    Great video! I was surprised the derivation of the neutron diffusion equation didn’t require much (if any) quantum mechanics. I imagined physics at the atomic level required qm for an accurate description

    • @jkzero
      @jkzero  5 месяцев назад +3

      Thanks, I am glad you liked the video. You make an important observation, yes, the diffusion equation doesn't explicitly involve any quantum mechanics (QM) but it is the case even in QM: as you probably know, quantum probability is conserved in QM, if you replace neutron density by probability density you will find that it satisfies the continuity equation shown at the beginning (without the two terms on the right-hand side). In that case the QM effects are hidden in the fact that the probability density is given by the complex-squared of the wave function, which in turn satisfies the Schrödinger equation. So where is QM in the neutron diffusion equation? I would say that in the cross-sections (which have to be calculated using QM or experimentally measured, in which case all QM effects are already weighted in) and in the neutron speed (notice that only the average neutron speed enters here but a more refined calculation would include a speed distribution).

  • @gordonnash5894
    @gordonnash5894 3 месяца назад

    Thanks!

    • @jkzero
      @jkzero  3 месяца назад

      Thanks to you for supporting me to keep making these videos!

  • @JazzBerri
    @JazzBerri 5 месяцев назад

    Didn’t expect to find such an interesting video on youtube!

    • @jkzero
      @jkzero  5 месяцев назад +1

      I am glad you liked it, I am curious to know what brings viewers to the channel, were you searching for something in particular or did the 'mighty algorithm' find you?

    • @JazzBerri
      @JazzBerri 5 месяцев назад

      ​@@jkzeroyour video appeared on my feed randomly yeah, wasn't looking for anything related.
      I guess that being a new chanel plus the "oppenheimer" related theme pushed your videos to many people.
      I hope more people are able to watch your videos!

    • @jkzero
      @jkzero  5 месяцев назад +1

      @@JazzBerri thanks fro sharing, I am glad the algorithm brought you here. Make sure to check the other videos, there is one dedicated to the real meaning of the famous words of Oppenheimer0 Welcome to the channel, I hope to have back for more videos.

    • @JazzBerri
      @JazzBerri 5 месяцев назад

      @@jkzero yeah thanks! I liked your chanel so im staying hahah

  • @DensityMatrix1
    @DensityMatrix1 5 месяцев назад

    Can you cover any usage of Markov Chains? I believe they were used/developed during the Manhattan Project for neutron diffusion. Thanks. Great video bro.

  • @randym1954tx
    @randym1954tx 2 месяца назад +1

    🚀Thanks!

    • @jkzero
      @jkzero  2 месяца назад

      I cannot thank you enough for your continuous support. The only way, I guess, is by keeping the videos coming and the quality standards. Thank you so much!

  • @GeoffryGifari
    @GeoffryGifari 5 месяцев назад +1

    Hmm when we set the boundary condition N(r,t) =0 at the critical radius R, is it equivalent to the Uranium ball being placed inside a shell of perfect neutron mirror?
    I still haven't entirely grasp how allowing neutrons to leak out (advanced solution) end up reducing the critical mass compared to the case of neutrons being trapped inside (basic solution) while in the basic case fission reaction also outputs more neutron compared to the input....
    That being said, hopefully your channel gets enough support so us audience can get more cool stuff like this.

    • @jkzero
      @jkzero  5 месяцев назад +1

      the "elementary" boundary condition is quite artificial, in fact, it is only used because it allows an analytical solution for R_c in terms of the relevant parameters (
      u, cross-sections, density, etc.), it gives a glimpse on what is driving R_c and the critical mass. However, this is the only reason this condition is introduced, mostly for academic purposes. The "advanced" boundary condition is what has to be used in a practical calculation.
      Thanks for watching and the positive feedback. Welcome to the channel.

  • @CaixetaArthur
    @CaixetaArthur 5 месяцев назад

    I love the video @jkzero! But may I ask you a question? I know that the video is not supposed to extremely technical and formal, but as a matter of curiosity: when you say that sin(kR)=0 implies kR=\pi (because it is simpler), would the solution be more "complete" if you consider every possible solution for sin(kR)=0? I don't wanna get into technicalities (eigenfunctions, Fourier expansion, etc.), but it got me curious. Thanks in advance!
    Edit: Extremely underrated content. I will definitely subscribe!

    • @jkzero
      @jkzero  5 месяцев назад +2

      thanks for the feedback and the sub, I am glad the video has been of interest for many people. You have a good point: the general solution is kappa*R=n*pi; however, the higher "harmonics" will only make the bomb core bigger, we are searching for the minimum radius. Also, a solution for n>1 can lead to a negative density, which is also a non-physical solution. Mathematically, your point is 100% correct; however, here is when physics constrains the final solution, just like setting B=0 to remove the non-physical cos(kappa*r) term.

    • @CaixetaArthur
      @CaixetaArthur 5 месяцев назад

      ​@@jkzeroperfect. Thanks a lot and keep up with the amazing work!

    • @jkzero
      @jkzero  5 месяцев назад

      @@CaixetaArthur thanks to you for watching and following up with your question. Welcome to the channel!

  • @CaptainCalculus
    @CaptainCalculus 5 месяцев назад +6

    That's worth a coffee☕

    • @jkzero
      @jkzero  5 месяцев назад +1

      Thank you so much for your continuous support, I will make sure to enjoy a coffee while being grateful for the appreciation of strangers on the other side of the planet that find this content of interest. Thanks again!

  • @pawelbjanowski
    @pawelbjanowski 5 месяцев назад +4

    Very good job! I think I will have one question but maybe later ;)

    • @jkzero
      @jkzero  5 месяцев назад

      questions are more than welcome, I hope that I can help

  • @MrFib112358
    @MrFib112358 5 месяцев назад +2

    Hi @jkzero, here is my solution to the homework problem of the cube of uranium: I used the Laplacian operator in Cartesian coordinates, solved the equations for all three axes, and found the "critical side" a_c=19 cm, which gives me a critical mass of 538 kg. Is this correct? How can using a cube instead of a sphere increase the critical mass by a factor 4?

    • @jkzero
      @jkzero  5 месяцев назад

      first of all: congrats on accepting the challenge, you got the critical side of the uranium cube correct but the critical mass seems to be off. Could you share how you got the mass from the size of the "critical cube"?

  • @suncrafterspielt9479
    @suncrafterspielt9479 5 месяцев назад

    very nice video thanks

    • @jkzero
      @jkzero  5 месяцев назад

      Thanks for watching; you can actively help the channel by liking, subscribing, and sharing. This type of support is highly appreciated so I can continue making videos.

  • @berdigylychrejepbayev7503
    @berdigylychrejepbayev7503 3 месяца назад

    to be honest, i liked the way you present the video with formulas and graphs but it is very hard to grasp and track. at least i as a non-physicist had a hard time to understand and could not finish the video. but i liked your type and really love to see chemistry version of you (or not this much of (unheard) definition and formulas). but I know that some people really like this kind of videos so dont worry about my comment.

    • @jkzero
      @jkzero  3 месяца назад

      thanks for your honest comment, one thing missing at the beginning is that this video is the response to the request of the step-by-step calculation of the critical mass in a previous video. If you want a more detailed description of the story behind this video please check "Critical Mass: when the atomic bomb got real" (ruclips.net/video/LduH7613QXw/видео.html). The present video by its own is pure mathematical physics but the other video provides the context and historical relevance of this calculation. I hope you like that video better.

  • @BorisNVM
    @BorisNVM 5 месяцев назад

    really cool video

    • @jkzero
      @jkzero  5 месяцев назад +1

      thanks for watching, make sure to check the full list of more videos on this topic

  • @peppescala4113
    @peppescala4113 5 месяцев назад +1

    How is the boundary condition expression for the advanced solution derived (20:59)? I'd really be thrilled to read the derivation, can you provide something? Is it something like Neumann vs Dirichlet boundary condition? Thank you for all of your work!

    • @jkzero
      @jkzero  5 месяцев назад +1

      Thanks fro the feedback, I am so excited that his video got so many people interested. The elementary boundary condition (BC) is a Dirichlet BC but it is quite artificial, the correct way of doing the calculation is by using a Neumann BC. I am glad someone got hooked on that detail, it is a lengthy derivation presented for the first time by Robert Serber, one of Oppenheimer's protegés, on his annotated version of The Los Alamos Primer. He only did it in Cartesian coordinates but it can be generalized to spherical coordinates chooser.crossref.org/?doi=10.2307%2Fj.ctvw1d5pf.6

  • @stevebeal73
    @stevebeal73 5 месяцев назад +1

    I graduated in physics in 1971 in Bristol, UK. When we were shown the neutron diffusion equation, our lecturer said "you all have all of the mathematical tools that you need to solve this equation, but I'm not allowed to show you how to do it, for reasons of national security." How times have changed!

    • @jkzero
      @jkzero  5 месяцев назад

      maybe your lecturer was joking or trying to make it more interesting because the use of diffusion theory for weapons design was rendered obsolete pretty soon, it is of little use to someone attempting to build a nuclear weapon, and even the reports including these calculations were declassified in 1965

    • @stevebeal73
      @stevebeal73 5 месяцев назад

      @@jkzero many thanks for your reply. Klaus Fuchs gained his PhD at Bristol and worked on the theoretical calculations for nuclear weapons. Many years later was imprisoned in the UK for passing nuclear secrets to the Russians. This may have led to the staff at Bristol becoming over cautious in protecting "secrets."

    • @jkzero
      @jkzero  5 месяцев назад +1

      @@stevebeal73 oh... yeah, I can see that. Fuchs' case was a big thing. It was quite sad from the Peierls family side

  • @abrikos1100
    @abrikos1100 4 месяца назад

    I think N=0 at boundary of the ball is not "no neutrons escape" boundary condition, it is more like "all neutrons trying to escape are vanishing" because neutron current J is nonzero through boundary

    • @jkzero
      @jkzero  4 месяца назад +1

      I think you are totally right. Honestly, the "elementary" boundary condition is quite artificial, in fact, it is only used because it allows an analytical solution for R_c in terms of the relevant parameters (
      u, cross-sections, density, etc.), it gives a glimpse of what drives R_c and the critical mass. However, this is the only reason this condition is introduced, mostly for academic purposes. The "advanced" boundary condition is what has to be used in a practical calculation.

  • @cyrilguillouard1266
    @cyrilguillouard1266 5 месяцев назад

    Beautiful video on a nice subject❤

    • @jkzero
      @jkzero  5 месяцев назад

      Glad you think so! I am curious to know what brings viewers to the channel, were you searching for something in particular or did the 'mighty algorithm' find you?

    • @cyrilguillouard1266
      @cyrilguillouard1266 5 месяцев назад

      @@jkzero ah the mighty algorithm

    • @jkzero
      @jkzero  5 месяцев назад

      @@cyrilguillouard1266 Thanks for sharing and I am glad the algorithm is working, I hope you find the other videos of interest too and welcome to the channel

  • @tobiaspeelen4395
    @tobiaspeelen4395 4 месяца назад

    Hi, love the video so far, but could you explain why at 12:35 the r’s cancel

    • @jkzero
      @jkzero  4 месяца назад

      Thanks, I am glad you like the video. Regarding your question, we need to take the limit r→0 of sin(κr)/r. It is tempting to say that this is 0/0 but the limit has to be taken with care. First, you expand the sin(κr) in the numerator as a Taylor series sin(κr) = κr - (κr)³/6 + ..., then divide this by r and you get sin(κr)/r ≈ κ - κ³r²/6 + ... and now you take the limit r→0 and only the first term survives, all the others contain r so they vanish, and you get lim r→0 sin(κr)/r ≈ κ. I hope this helps.

  • @Revyy729
    @Revyy729 5 месяцев назад +1

    Pheonomenal content. Subscribed! as a side note, I am a highschool student aiming to pursue a career in theoretical physics. How do you suggest I proceed with that? All and any help would be appreciated

    • @okkam7078
      @okkam7078 5 месяцев назад

      Go to college for physics, pursue a Masters (i.e. get good grades in undergrad), maybe pursue a PhD (you’ll figure that out when you get closer to it)

    • @Revyy729
      @Revyy729 5 месяцев назад

      Thank you!​@@okkam7078

  • @lukesiemens5421
    @lukesiemens5421 4 месяца назад

    I think that there is a small error at 13:35. At this point in the video it is stated that the third boundary condition is that no neutrons escape the sphere and this is done by requiring that N(r=R, t)=0. Then there is a quick disclaimer that this condition is not completely true and the more advanced solution is briefly mentioned. I may just be misunderstanding it, but I understood this to mean that the condition N(r=R, t)=0 does imply that no neutrons will escape the core however this boundary condition is not practical or optimal. Using the solution that satisfies the condition N(r=R, t)=0, I calculated the total neutron flux at the boundary from the neutron density and found that it is not zero but proportional to R*N(r=0, t). I am not criticizing the use of the boundary condition N(r=R, t)=0, but I think it is incorrect to say that this condition is equivalent to requiring that no neutrons escape the core.

    • @jkzero
      @jkzero  4 месяца назад

      Thanks for your exposition. Honestly, the "elementary" boundary condition is quite artificial, in fact, it is only used because it allows an analytical solution for R_c in terms of the relevant parameters (
      u, cross-sections, density, etc.), it gives a glimpse on what is driving R_c and the critical mass. However, this is the only reason this condition is introduced, mostly for academic purposes. The "advanced" boundary condition is what has to be used in a practical calculation.

  • @NicholasMati
    @NicholasMati 4 месяца назад

    When doing problems like this, there are an infinite number of solutions that satisfy the boundary condition (just keep adding pi for the Neumann BC). The coefficients for these solutions depend on the starting conditions over the whole domain. However, the higher harmonics change the value of kappa which changes the exponential growth rate. For heat diffusion problems, the exponent gets more negative which causes fine structures to quickly damp out leaving the first few harmonics. Here, it looks like the negative in front of nu' also means that it is negative in the exponent (or less positive) for higher harmonics. However, it should still contribute somewhat to the overall reaction. For bombs that have a neutron generator in the pit, how important is the initial distribution to the critical mass?

    • @jkzero
      @jkzero  4 месяца назад

      You are right: the general solution is kappa*R=n*pi; however, the higher "harmonics" will only make the bomb core bigger, we are searching for the minimum radius. Also, a solution for n>1 can lead to a negative density, which is also a non-physical solution. Mathematically, your point is 100% correct; however, here is when physics constrains the final solution, just like setting B=0 to remove the non-physical cos(kappa*r) term.

  • @Oliveir51
    @Oliveir51 5 месяцев назад

    This is why neutron cross section value is so important. Was measured in Argonne or was it Fermi lab ?

  • @expomax01
    @expomax01 5 месяцев назад

    At 10:10, the expression for the Laplacian operator in spherical coordinates is wrong. Inside the parenthesis, there should be r squared, not 1 over r squared. Apart from this, very interesting video!

    • @jkzero
      @jkzero  5 месяцев назад

      Thanks, I am glad you liked the video. You are totally right, I messed up the r² term. Thanks so much for pointing this out, it is a good catch. Fortunately, it doesn't affect the calculation because it is correctly applied, this is just a typo on the video but not an error in the calculation. Others have also notified it and I included an erratum in the video description. I do my best to avoid these typos but after watching everything many times some minor details slip through. Thanks again.

  • @t850
    @t850 5 месяцев назад +1

    ...excellent video as always. One question immedietly pops up to my mind.
    If criticality condition equals exactly 0 does this theoretically mean that the fuel won't explode but "burn" as a self sustaining nuclear reactor until whole mass is "consumed"?
    On th side note by knowing critical mass we can now calculate minimal theoretical yield of this device. Since critical mass is about 50 kg and energy denistiy of uranium/plutonium is 17 kt/kg minimal theoretical yield would be around 850 kt of TNT.
    Realistic result is unfortunately much lower as the Trinity test (similar device design) produced yiel of only 25 kt. Mass/energy convertion factor was much lower since majority of the mass was blown up (unused).
    Btw, I always wondered what is "cost/benefit" analysis of uranium enrichment in terms of energy invested and energy gained. Is there any information as to how much energy was used to produce first atom bomb since enriching facility is generally massive. So massive in fact that even with this secret known, one needs a facility the size of a small town to produce a bomb. It would be interesting to know if one needs to invest more energy to enrich the uranium than the bomb would produce...:)

    • @kwzieleniewski
      @kwzieleniewski 5 месяцев назад

      IIRC fuel for nuclear power plants needed as much as 4% of energy released for the enrichment. There are equations on Wikipedia saying how much more energy it takes to enrich to 90% than to just 4%.
      However, this also depend on how you enrich. Zippe type centrifuges are far more efficient than gasous diffusion.

    • @t850
      @t850 5 месяцев назад

      @@kwzieleniewski ...you mean to say that 4% of (thermal) energy released from the nuclear fuel is needed to enrich that exact ammount of uranium ore to be graded as nuclear fuel? Thing is that 4% (3-5%) is basically average ammount of the enrichemnt that ore needs to be processed to to be classified as nuclear fuel so I don't want to be confused with this number similarity.
      Beside the obvious difference in enrichemnt methods there is probably some efficiency drop as the enrichment goes to the higher levels because extracting the remaining mass is that much harder since the saturation of U235 is pretty high and mass difference between U235 and U238 almost insignificant...:/

  • @Tim-Kaa
    @Tim-Kaa 5 месяцев назад

    Nice

  • @nathansudermann-merx4586
    @nathansudermann-merx4586 5 месяцев назад +1

    Great video, Jorge. Thank you! In the method of separating variables you assume that N(x,t) can be factorized. Is there a reason why this makes sense from a physics point of view?

    • @digxx
      @digxx 5 месяцев назад +2

      The reason is mathematical: Given certain conditions, you can prove uniqueness of the solution and thus the solution he found (obeying the boundary conditions) can be the only solution. On the other hand, if he would have chosen more complicated geometries and boundary conditions, the separation ansatz would probably not have worked.

    • @nathansudermann-merx4586
      @nathansudermann-merx4586 5 месяцев назад

      Makes sense, thank you!@@digxx

    • @jkzero
      @jkzero  5 месяцев назад

      Physicists are quite informal in general but specially when solving equations; instead of trying to prove that a solution exists and that it is unique, we just take an ansatz, work it out, and if a solution is found then that is the solution. Period. The factorization of partial solutions is just an ansatz that works. It is the first thing we always try. If it fails, then we try more formal methods.

    • @nathansudermann-merx4586
      @nathansudermann-merx4586 5 месяцев назад

      @@jkzero yeah, that is what you have mathematicians for who will prove that everything you have done was right ;). I was just wondering wether another ansatz would have yield another solution but that was answered by the reply of @digxx.

    • @jkzero
      @jkzero  5 месяцев назад +1

      @@nathansudermann-merx4586 we do many things in physics that would drive mathematicians mad, like "multiplying" or "dividing" by dx when integrating

  • @fightingforcatalonia
    @fightingforcatalonia 5 месяцев назад +2

    donde trabajeste con neutrones. Central nuclear? que gran trabajo

    • @jkzero
      @jkzero  5 месяцев назад

      soy físico teórico, todo esto lo aprendí por mi cuenta por gusto, simplemente me fascina el tema y creé este canal porque pensé que era hora de compartir la pasión, afortunadamente la recepción ha sido fenomenal lo que me anima a seguir creando este tipo de videos. Cuéntame que te trajo por acá, ¿te encontró el algoritmo?

    • @fightingforcatalonia
      @fightingforcatalonia 5 месяцев назад

      Algo parecido al algoritmo. eres el 3blue1brown de la fisica. Yo soy físico también. Un saludo@@jkzero

    • @jkzero
      @jkzero  5 месяцев назад

      "the mighty algorithm", me alegro que te haya traído por acá. ¿En tipo de física trabajas? ¿Eres estudiante, postdoc, académico o estás en la industria?

  • @peppescala4113
    @peppescala4113 5 месяцев назад

    Hi! Is there a way to predict theoretically the fission and capture cross-sections of our processes? You can get more or less the order of magnitude of this treating the U-235 as an hard sphere and using the Bohr-Wheeler theory that tells you that the higher Z²/A, the higher the fission probability. I tried to evaluate the U-235 fission cross section and I get something like 6.20 time 10^(-23) cm². But I'm not able to see how to improve this little model in order to study capture cross section. Also I have no idea how to distinguish slow and fast neutrons! By the way, your videos are beautiful.

    • @jkzero
      @jkzero  5 месяцев назад +1

      Thanks for the feedback and I am glad you liked the videos. Yes, you can motivate the order of magnitude for the cross-sections but at the scale of neutrons and nuclei, quantum-mechanical effects can become relevant and classical intuition falls short. Plus a nucleus with 235 nucleons is really hard to model and you can only use approximation methods, in the end the best way is just go a measure it, all you need is a few microgram of uranium-235, which is easy to "make" with a cyclotron

  • @alexkalug
    @alexkalug 5 месяцев назад

    Great video! However, I think that your explanation of he trivial b.c. is a bit misleading. We are not 'disallowing the neuterons to escape' the sphere, however we set the radius of the sphere where N=0. Hence a larger Rc, we artificially require the spere to end only where the neuterons would have stopped, assuming the same λ_t. But I'm nitpicking

    • @jkzero
      @jkzero  5 месяцев назад

      thanks, I am glad you liked the video. Honestly, "elementary" boundary condition is quite artificial, in fact, it is only used because it allows an analytical solution for R_c in terms of the relevant parameters (
      u, cross-sections, density, etc.), it gives a glimpse on what is driving R_c and the critical mass. However, this is the only reason this condition is introduced, mostly for academic purposes. The "advanced" boundary condition is what has to be used in a practical calculation.

  • @chromodyn7769
    @chromodyn7769 5 месяцев назад

    In honor of Barbenheimer.... calculate the crtitical mass for a doll shaped bomb LMAO

  • @dobotube
    @dobotube 3 месяца назад

    Amazing channel and videos. A question:
    I'm very puzzled by the various prompt critical excursions with the demon core and friends that occured in history.
    I understand nuclear reactors rely on delayed neutrons to keep a reaction in balance.
    I understand a prompt supercritical state will have a neutron population grow extremely fast.
    I understand an unconstrained core will tear itself apart before full combustion.
    What I don't understand is how can a core balance in prompt critical state (blue flash and all) and neither fizzle nor explode.
    I've been given explanations like Doppler broadening and thermal expansions which certainly could fit the bill but, without knowing how quickly temperature changes vs how quickly the neutron population grows, I find unconvincing.
    Is there a way to do back of the envelope calculations to think about these prompt critical balanced reactions?

    • @jkzero
      @jkzero  3 месяца назад

      my most sincere and honest answer is: I do not know. I am pretty sure a nuclear engineer would know the answer, they study these reactions in non-ideal scenarios in detail. As a physicist myself, we mostly do calculations in ideal conditions, prove the feasibility based on the physics, conservation laws, etc. and then leave the dirty (but crucial) dirty work to engineers, they really make things happen.

  • @parsecgilly1495
    @parsecgilly1495 5 месяцев назад

    for a more precise calculation, however, the transport equation should be used

    • @jkzero
      @jkzero  5 месяцев назад +1

      you are totally right, this is one of the things what I meant when I said that the calculation shown in the videos involves many approximations. But it could be argued that diffusion theory is not even applicable because the size of the material is barely bigger than the transport mean free path. Careful calculations show that many effects ignored here cancel each other out but this is just a lucky coincidence.

    • @parsecgilly1495
      @parsecgilly1495 5 месяцев назад

      in any case yours is an excellent video, I have a degree in nuclear engineering and in the nuclear reactor physics courses I had followed, the first approach to use for calculating the distribution of neutrons was to use the diffusion equation, we then had the most precise calculations with the transport equation
      @@jkzero

  • @henning4070
    @henning4070 4 месяца назад

    Thank you very much for your great work, your video explains the neutron diffusion in a clear way and I enjoyed watching it.
    I have a comment on the boundary condition N(r=R, t) = 0. This does not mean that no neutrons escape the core. Neutron flux is not induced by neutron density but by the gradient of the neutron density, just like heat flow isn't caused by temperature but by a temperature gradient. The gradient of N at R and therefore the flow of neutrons out of the core is not zero at v'=0, this can also be seen by realising, that the neutron density is constant in time for v'=0, but neutrons are produced by fission. Where do the produced neutrons go if their density remains constant in time? They must leave the core through its surface.
    Nevertheless, N(r=R, t) = 0 is a reasonable boundary condition because N(r=R, t) > 0 would cause a discontinuity of the neutron density at r=R, leading to an unphysical singularity in neutron flux.

    • @jkzero
      @jkzero  4 месяца назад +1

      Thanks for your extensive exposition, I appreciate that. Honestly, the "elementary" boundary condition is quite artificial, in fact, it is only used because it allows an analytical solution for R_c in terms of the relevant parameters (
      u, cross-sections, density, etc.), it gives a glimpse on what is driving R_c and the critical mass. However, this is the only reason this condition is introduced, mostly for academic purposes. The "advanced" boundary condition is what has to be used in a practical calculation.

    • @henning4070
      @henning4070 4 месяца назад

      ​@@jkzeroIt would be interesting to know where the term for the advanced boundary condition comes from.

    • @jkzero
      @jkzero  4 месяца назад +1

      @@henning4070 you can find a proof of the correct boundary condition in Cartesian coordinates in the appendix of the Los Alamos Primer, freely available online

    • @henning4070
      @henning4070 4 месяца назад

      @@jkzero Thank you very much, I will have a look.

  • @easondu9236
    @easondu9236 5 месяцев назад +1

    Is this made on manim? Very nice animation

    • @jkzero
      @jkzero  5 месяцев назад +3

      yes, I use Manim in all my videos. I know that Grant recommends to not use Manim for just showing equations but I can't help it, I like how LaTeX formulas are shown too much. Thanks for stopping by and make sure to check the other videos.

  • @hypercomms2001
    @hypercomms2001 28 дней назад

    I am interested in taking up the challenge and following through the procedure to attempt to calculate the critical mass for a cylinder. I do have a background in electrical engineering, and as part of that I studied quantum mechanics, and nuclear physics, and did extensive training in calculus.... However I am more interested in trying to calculate the critical mass that relates to a nuclear reactor instead of a nuclear weapon, would you be able to do a video on that or point to where one could find information how to make such a calculation? I am thinking something like calculating the critical mass for something like a Chicago pile, CP-1... minus the bottle of Chianti Bertolli wine... and Enrico fermi!

    • @jkzero
      @jkzero  27 дней назад

      I have the suspicion that the procedure followed here might serve as a guide for solving the neutron diffusion equation, but the controlled conditions in a reaction can be quite different. I am not a nuclear engineer but maybe the lectures on the MIT course Introduction to Nuclear Engineering and Ionizing Radiation might be more useful (lectures are on RUclips).

    • @hypercomms2001
      @hypercomms2001 26 дней назад

      @@jkzero Thank you Dr Diaz, I do appreciate your suggestion, and I have searched for the course, and bookmarked it, as it is something that I will attempt at the end of the year, as it will be an interesting and cool challenge, as I am currently completing a Masters of IT, and wont have the bandwidth until then.

  • @s5101953
    @s5101953 5 месяцев назад

    In the remarkable story two flaws are negating each others. First, Diffusion is over estimated for fast neutrons (of bomb core) and for thermal ones as well. Second, the source neutron term due to fission must be decreased for capturing neutrons. But at criticality threshold these mistakes happly are zeroesed. Good luck!

    • @jkzero
      @jkzero  5 месяцев назад

      it is interesting how the simple calculation in this video was enough despite the many limitations of the formalism. The most remarkable is the fact that two non-negligible effects have been ignored: 1. neutrons can also been captured by U235 without fissioning, this was completed ignored; and 2. the elastic-scattering cross section is not isotropic so care must be taken when calculating the corresponding cross-section. Nonetheless, when these two effects are incorporated they happen to cancel each other out. Just a lucky coincidence.

  • @erikziak1249
    @erikziak1249 Месяц назад

    I was somewhat capable to follow what this video is about. It is more than two decades since I studied math like this and I never needed it in my job. So all my knowledge, which was not particularly deep in the first place, is all but gone. A couple of years ago I opened my notes that I made during math courses at the university. I was surprised to find that I even learned about Hamiltonians. I did not forget just how to apply them and calculate them, I completely forgot that they exist! Math was not my cup of tea, I always treated it as a necessary evil.

    • @jkzero
      @jkzero  29 дней назад

      Thanks for sharing and great to learn that the content can be followed. This was a particularly hard video, it has no story behind, it is a purely mathematical-physics exercise. I forgot to point this out at the beginning but this video was a response to a popular request to show this calculation in the comments on my video with the history of the nuclear critical mass: ruclips.net/video/LduH7613QXw/видео.html

  • @tatytatytaty95
    @tatytatytaty95 5 месяцев назад

    Great Video!
    In the continuity equation dN/dt describes the change in neutron density due to neutron absorption and emission, while del J describes the change in neutron density due to neutrons diffusing away. So why is the entire expression (dN/dt+ del J) equal to the neutron sink and source term (S+ + S-)? Shouldn't there be an additional term on the right side, accounting for the neutrons being lost due to their motion out of the volume that is being considered?

    • @jkzero
      @jkzero  5 месяцев назад +1

      the term ∂N/∂t describes the change on neutron density in a infinitesimal volume with time, that is all, the reasons are given by other terms in the equation. The term ∇·J describes the differential number of neutrons moving in and out of the infinitesimal volume (they can move in and out by any reason, in our case by diffusion when we impose Fick's law). If the number of particles were fixed then the RHS of the equation would be zero, like in the case of electric currents or fluid motion, but here we need to tell the equation that inside the infinitesimal volume the neutron number is not conserved and in fact neutrons are created and destroyed, this is why the two terms on the RHS are needed.

    • @kwzieleniewski
      @kwzieleniewski 5 месяцев назад

      This equation is local, that is, it describes what happens in an infinitesimally small volume around a particular point (of coordinates x). Sphere boundary is far away.

    • @jkzero
      @jkzero  5 месяцев назад

      @@kwzieleniewski not sure if you are trying to make a particular point or asking for something in particular

  • @user-ek9vo2ub9b
    @user-ek9vo2ub9b 5 месяцев назад +1

    For the Pu239 density in g/cm^(3) - isn't it greater than the figure you provide?

    • @jkzero
      @jkzero  5 месяцев назад +2

      you are right that the standard density of plutonium is closer to 20 g/cm^3; however, depending on temperature plutonium exhibits different crystalline phases. The relevant phase used for weapons is the so-called “delta” phase and that is the density that I have provided.

    • @user-ek9vo2ub9b
      @user-ek9vo2ub9b 5 месяцев назад

      @jkzero Huh. Thanks for that - so basically the standard density is not in that phase at all. Interesting.

    • @jkzero
      @jkzero  5 месяцев назад

      @@user-ek9vo2ub9b exactly; I mean you can have plutonium metal in other phases too but for the bomb they used the delta phase due to stability and for machining purposes

  • @rustym.shackelford5546
    @rustym.shackelford5546 28 дней назад +1

    Okay, I am a little stumped at the laplacian part. So I already specifically got the form down pretty well: (nabla)^(2)*F + (kappa)^(2)*F = 0 but now how exactly do I plug in the "auxiliary function" to this figure in order to get the positional function. At least just give me a hint or point me in a general direction that I could learn how to do this (note: I am familiar with Calculus, however, I am not entirely versed in doing these equations on my own).

    • @jkzero
      @jkzero  27 дней назад +1

      Replace u(r) = r*f(f) in the Laplacian ∇²f= (1/r²)*d/dr(r² df/dr), expand all the terms, simplify what can be simplified, and you will get now get the harmonic-motion equation for u(r). Please note that there is a typo on the video: at 10:10 the Laplacian of f should say: ∇²f= (1/r²)*d/dr(r² df/dr) (not 1/r² in the derivative bracket), this is a typo and not an error in the calculation. I hope this helps. If not just let me know.

    • @rustym.shackelford5546
      @rustym.shackelford5546 27 дней назад

      @@jkzero Major thanks.

  • @user-ek9vo2ub9b
    @user-ek9vo2ub9b 5 месяцев назад

    Have a question: How would you input the N(r , t) function into a Graphing Calculator alongside the Criticality Factor equation?

    • @jkzero
      @jkzero  5 месяцев назад

      I have not touched a graphing calculator in decades (when I spent all my savings in a TI-92+ and barely used it). I would expect that you can define custom functions but I guess the easiest way would be explicitly replacing nu' (criticality factor) in the definition of g(t). Note that you will need to define a radius, say 30 cm, that is bigger than the critical value; otherwise, your solution will be constant (because nu'=0 for R=R_c)

    • @user-ek9vo2ub9b
      @user-ek9vo2ub9b 5 месяцев назад

      Can you simplify what you mean a bit? Still kind of new to this. ​@jkzero

    • @jkzero
      @jkzero  5 месяцев назад

      @@user-ek9vo2ub9b I honestly don't know how to code this in a graphing calculator but I have created a Python code to generate the plots shown in the video you can find in this link github.com/jsdiazpo/JK0/blob/main/11-Neutron-Diffusion-Eq.ipynb Please let me know if this helps

  • @Archiekunst
    @Archiekunst 4 месяца назад

    This is so interesting. I bet a similar equation can be written for a biological cell producing some biomolecule in its 'core' (cytoplasm or nucleus or some organelle) and then it is flowing out. What do you think Jorge?

    • @jkzero
      @jkzero  4 месяца назад

      This topic goes beyond my expertise so whatever I say should be taken with a grain of salt. Moreover, I do not like it when people comment too confidently on topics beyond their expertise. With that said, you can search for modeling of so-called "reaction-diffusion systems" that include many biological and chemical systems. There you will encounter Fick's law again and the diffusion equation for modeling, for example, patterns in ecological systems or the diffusion of ions in living systems. There is a whole branch that makes use of this type of equation called mathematical biology.

    • @Archiekunst
      @Archiekunst 4 месяца назад

      ​@jkzero The math should be equivalent. I do research on biochemical systems and can vouch for the biology side. Recently I've been teaching myself some modelling which is why I came across your videos. Perhaps I can pick your brain some time to validate my equations. I work at the max planck in dresden. I like the idea where you can set the boundary flux to be zero or nonzero. It's just like a cell's receptors not letting stuff out for zero flux or secreting something, for nonzero.

    • @jkzero
      @jkzero  4 месяца назад

      @@Archiekunst you can find my email on the "about" of the channel

  • @Tadesan
    @Tadesan 5 месяцев назад

    Separation of bariables

  • @chonchjohnch
    @chonchjohnch 5 месяцев назад

    How is the mass being destroyed or heavily damaged by it partially fissioning accounted for?

    • @jkzero
      @jkzero  5 месяцев назад

      in this calculation the critical mass is defined as the minimum amount of fissile material to trigger a self-sustained nuclear chain reaction. Of course as soon as the reaction starts, the material will evaporate and the whole assembly will self-destroy. This is not accounted for in the calculation shown in the video.

  • @001firebrand
    @001firebrand 5 месяцев назад

    "Boom! Science!" 💖

  • @lucacastronuovo9789
    @lucacastronuovo9789 5 месяцев назад +2

    Being sigma_t greater then sigma_f (min 5:31), how is it possible that lambda_t is greater than lambda_f (min 15:58)?

    • @jkzero
      @jkzero  5 месяцев назад

      oh no... you are totally right, I messed up the subscripts at 15:58, they should be reversed, it should say: lambda_f = 16.89 cm and lambda_t = 3.596 cm. I will add an errata in the description. Thanks so much for pointing this out, it is a good catch. Fortunately, it doesn't affect most of the elementary calculation because the two mean-free paths are multiplied but it does make as difference in the 'advanced' solution, where lambda_t appears isolated. Thanks again. I do my best to avoid these typos but after watching everything many times some minor details slip through. Thanks also for asking in such a nice way contrary to patronizing mode of writing comments favored on the internet.

    • @lucacastronuovo9789
      @lucacastronuovo9789 5 месяцев назад

      @@jkzero You're welcome! Anyway, thank you for the amazing video and the explanation!

  • @kotowaty3386
    @kotowaty3386 Месяц назад

    Hey. I like this explanation however I don't understand how to obtain tau term. I feel like one can get a formula for this however I'm stuck

    • @kotowaty3386
      @kotowaty3386 Месяц назад

      Should I just calculate an average speed from average neutron energy? What value is that average energy?

    • @jkzero
      @jkzero  Месяц назад

      tau can be determined in terms of the fission mean-free-path and the average neutron speed. I showed the detail of this in the video about the Physics of a Nuclear Explosion at the 4:04 mark ruclips.net/video/6VSrGDOrWXc/видео.htmlsi=5uq4tpjHFlhVPBah&t=244

    • @kotowaty3386
      @kotowaty3386 Месяц назад

      @@jkzero Everything clear now. Thank you for the answer

    • @jkzero
      @jkzero  Месяц назад

      @@kotowaty3386 excellent! I'm glad that helped. Thanks to you for watching and for following the calculations along. There are some typos later in the video, please check the errata in the video description.

  • @kondojunagaraju8899
    @kondojunagaraju8899 5 месяцев назад

    What is mean conceptual understanding? Please tell me example.

    • @jkzero
      @jkzero  5 месяцев назад

      this video is the technical response to an early video that I published titled "Critical Mass: when the atomic bomb got real." If you want the conceptual and historical account for what it is explained in the present video you check it out here ruclips.net/video/LduH7613QXw/видео.html

  • @claudiomenchinelli7465
    @claudiomenchinelli7465 10 часов назад

    DIFFICULTY LEVEL: TISE FOR THE HYDROGEN ATOM